Support and discussions for Molcas and OpenMolcas users and developers
You are not logged in.
Please note: The forum's URL has changed. The new URL is: https://molcasforum.univie.ac.at. Please update your bookmarks!
You can choose an avatar and change the default style by going to "Profile" → "Personality" or "Display".Dear Users,
How do I pass module-specific keywords to invoked modules?
In particular, I have in mind making MCKINLEY call SLAPAF with a CART keyword.
&SLAPAF; CART does the job, but &MCKINLEY; CART does not - cuz CART is not a MCKINLEY keyword.
&SLAPAF; CART ; &MCKINLEY; CART results in Cartesian gradients (the first part), but then vibrations are still done in internals.
So, what do I do?
Offline
You can't. What you can do is call SLAPAF directly, instead letting MCKINLEY do it for you:
> do while
&SLAPAF
Numerical
Iterations = 0
Thermochemistry
>end do
But if what you want is to obtain the full 3N set of frequencies (including rotations and translations) instead of the 3N-6 pure vibrations, I don't that will help. For that you can just add "VarT ; VarR" in GATEWAY.
Offline
Thanks.
Is this a new functionality?
My version (24.02, 455-gff0b5d1cb-dirty) does not like it:
>>>>> Input file for module SLAPAF <<<<<
Numerical
Iterations
0
Thermochemistry
****** Error *******
End of input
Edit: today's 334-gcea4da31a also does not like it .
Last edited by andrewshyichuk (2024-09-30 15:20:34)
Offline
Sorry, I just copied sloppily. It should be a block ended by "End of PT". The syntax is in the documentation: "The keyword must be followed by different entries containing the Rotational Symmetry Number, the Pressure (in atm), and one entry per Temperature (in K) for which the thermochemistry will be calculated." You can skip it if just want the frequencies and no thermochemistry.
Offline
It worked wonders! No imaginary frequencies at all.
Offline