1
\section{Geometry Optimization and Vibrational Frequency Analysis} \label{opt}
1
\section{Geometry Optimization} \label{opt}
3
3
\PSIthree\ is capable of carrying out geometry optimizations (minimization
4
4
only, at present) for a variety of molecular structures using either analytic
5
5
and numerical energy gradients.
7
When analytic gradients are available (see Table \ref{table:methods}),
8
\PSIthree\ will automatically generate and use redundant, simple internal
7
When present, internal coordinates provided in the INTCO: section of the
8
input will be read and used by \PSIthree. If these are missing, \PSIthree\
9
will automatically generate and use redundant, simple internal
9
10
coordinates for carrying out the optimization. These simple stretch, bend,
10
torsion, and linear bend coordinates may be read from the intco.dat file
11
or determined completely by distance criteria using the input geometry.
11
torsion, and linear bend coordinates are determined by distance
12
criteria using the input geometry.
12
14
By default, optimization is performed in redundant internal coordinates
13
15
regardless of how the geometry was provided in the input. Alternatively,
14
16
the user may specify zmat\_simples=true, in which case the simple internal
15
17
coordinates will be taken from the ZMAT given in the input file. Also,
16
the user may specify optimization in delocalized internal coordinates
18
the user may specify optimization in non-redundant, delocalized internal coordinates
19
with delocalize=true. In this case, the automatically generated simple
20
coordinates are delocalized and redandancies are removed. Advanced users
21
may wish to specify the simple internal coordinates in the intco.dat file, and
22
then allow \PSIthree\ to delocalize them.
19
Geometrical constraints may be imposed by the addition of a section
20
with nearly the same format as intco.dat. For example, to fix the distance
24
Only those coordinates or combinations of coordinates that are specified
25
by the "symm =" keyword in the INTCO: section are optimized. Coordinates can
26
be approximately frozen by commenting them out within the "symm =" section.
27
Geometrical constraints may be precisely imposed by the addition of a section
28
with nearly the same format as in INTCO:. For example, to fix the distance
21
29
between atoms 1 and 2, as well as the angle between atoms 2, 1 and 3
22
30
in an optimization, add the following to your input file.
31
42
The constrained simple internals must be ones present (either manually or
32
automatically) among the simple internals in intco.dat. Alternatively,
43
automatically) among the simple internals in the INTCO: section. Alternatively,
33
44
the z-matrix input format may be used to specify constrained optimizations.
34
45
If zmat\_simples=true, then variables in the z-matrix which end in
35
46
a dollar sign will be taken as simple internals to be optimized, and
36
47
all other variables will be taken as simple internals to keep frozen.
38
For methods for which only energies are available, \PSIthree\ will use
49
To aid optimizations, force constants may be computed using "jobtype = symm\_fc".
50
The determined force constants will be saved in a binary file PSIF\_OPTKING
51
(currently file 1). Subsequent optimizations will read and use these force
52
constants. In general, \PSIthree\ looks for force constants in the following
53
order: in this binary file, in the FCONST: section of the input, and in the fconst.dat
54
file. If no force constants are found in any of these, then an empirical
55
diagonal force constant matrix is generated.
57
For methods for which only energies are available, \PSIthree\ will use non-redundant,
39
58
symmetry-adapted delocalized internal coordinates to generate geometrical
40
59
displacements for computing finite-difference gradients. The simple
41
60
coordinates can be linearly combined by hand or automatically. The goal
42
61
is to form 3N-6(5) symmetry-adapted internal coordinates. The automated
43
62
delocalized coordinates may work for low-symmetry molecules without
44
linear angles, but has not been extensively tested. For both analytic-
63
linear angles, but have not been extensively tested. For both analytic-
45
64
and finite-difference-gradient optimization methods, Hessian updates are
46
65
performed using the BFGS method.
67
The list below shows which coordinates are used by default for different types of jobs. \\
68
jobtype=freq dertype=first symmetry-adapted cartesians \\
69
jobtype=freq dertype=none symmetry-adapted cartesians \\
70
jobtype=fc dertype=first delocalized internals (or user-defined SALCs) \\
71
jobtype=symm\_fc dertype=first delocalized internals (or user-defined SALCs) \\
72
jobtype=opt dertype=first redundant internals \\
73
jobtype=opt dertype=none delocalized internals (or user-defined SALCS) \\
76
The following keywords are pertinent for geometry optimizations.
78
\item[JOBTYPE = string]\mbox{}\\
79
This keyword must be set to {\tt OPT} for geometry optimizations and
80
{\tt SYMM\_FC} to compute force constants.
81
\item[DERTYPE = string]\mbox{}\\
82
This keyword must be set to {\tt NONE} if only energies are available
83
for the chosen method and {\tt FIRST} if analytic gradients are available.
84
\item[CONV = integer]\mbox{}\\
85
The maximum force criteria for optimization is $10^{-conv}$.
86
\item[BFGS = boolean]\mbox{}\\
87
If true (the default), a BFGS Hessian update is performed.
88
\item[BFGS\_USE\_LAST = integer]\mbox{}\\
89
This keyword is used to specify the number of gradient step for the BFGS
90
update of the Hessian. The default is six.
91
\item[SCALE\_CONNECTIVITY = float]\mbox{}\\
92
Determines how close atoms must be to be considered bonded in the automatic
93
generation of the bonded list. The default is 1.3.
94
\item[DELOCALIZE = integer]\mbox{}\\
95
Whether to delocalize simple internal coordinates to attempt to produce
96
a symmetry-adapted, non-redundant set.
97
\item[MIX\_TYPES = boolean]\mbox{}\\
98
If set to false, different types of internal coordinates are not allowed
99
to mix in the formation of the delocalized coordinates. Although this
100
produces cleaner coordinates, often the resulting delocalized coordinates
101
form a redundant set.
102
\item[ZMAT\_SIMPLES = boolean]\mbox{}\\
103
If set to true, the simple internal coordinates are taken from the zmat
104
entry in the input file. The default is false.
105
\item[POINTS = 3 or 5]\mbox{}\\
106
Specifies a 3-point or a 5-point formula for optimization by energy points.
107
\item[EDISP = float]\mbox{}\\
108
The default displacment size (in au) for finite-difference computations. The
110
\item[FRAGMENT\_DISTANCE\_INVERSE = boolean]\mbox{}\\
111
For interfragment coordinates. If true, then 1/R(AB) is used, if false,
112
then R(AB) is used. The default is true.
113
\item[FIX\_INTRAFRAGMENT = boolean]\mbox{}\\
114
If true, all intrafragment coordinates are constrained.
115
\item[FIX\_INTERFRAGMENT = boolean]\mbox{}\\
116
If true, all interfragment coordinates are constrained.
117
\item[DUMMY\_AXIS\_1 = 1 or 2 or 3]\mbox{}\\
118
Specifies the axis for the location of a dummy atom for the definition
119
of a linear bending coordinate. The default is 2.
120
\item[DUMMY\_AXIS\_2 = 1 or 2 or 3]\mbox{}\\
121
Specifies the axis for the location of a dummy atom for the definition
122
of a linear bending coordinate. The default is 3.
123
\item[TEST\_B = boolean]\mbox{}\\
124
If set to true, a numerical test of the B-matrix is performed.
125
\item[PRINT\_FCONST = boolean]\mbox{}\\
126
If set to true and jobtype=symm\_fc, then the force constants will
127
be written to the fconst.dat file. This allows force constants to be
128
reused even if the binary PSIF\_OPTKING file is no longer present.
129
\item[Print options]\mbox{}\\
130
The following when set to true, print additional information to the
131
output file: PRINT\_SIMPLES, PRINT\_PARAMS, PRINT\_DELOCALIZE,
132
PRINT\_SYMMETRY, PRINT\_HESSIAN, PRINT\_CARTESIANS.
133
\item[DISPLACEMENTS = ( (integer float ...) ...)]\mbox{}\\
134
A user may specify displacments along internal coordinates using this
135
keyword. For example, displacements = ( (2 0.01 3 0.01) ) will compute
136
a new cartesian geometry with the second and third internal coordinates
141
\section{Vibrational Frequency Computations} \label{freq}
48
142
\PSIthree\ is also capable of computing harmonic vibrational frequencies
49
143
for a number of different methods using energy points or analytic energy first or
50
144
second derivatives. (At present, only RHF-SCF analytic second derivatives
51
145
are available.) If analytic energy second derivatives are not available,
52
146
\PSIthree\ will generate displaced geometries along symmetry adapted cartesian
53
coordinates, compute the appropriate energies or first derivatives, and use finite-difference
54
methods to compute the Hessian.
147
coordinates, compute the appropriate energies or first derivatives, and use
148
finite-difference methods to compute the Hessian.
56
The following keywords are pertinent for geometry optimizations and
57
vibrational frequency analyses:
150
The following keywords are pertinent for vibrational frequency analyses:
58
151
\begin{description}
59
152
\item[JOBTYPE = string]\mbox{}\\
60
This keyword (described earlier in this manual) must be set to
61
{\tt OPT} for geometry optimizations and {\tt FREQ} for frequency analyses.
153
This keyword must be set to {\tt FREQ} for frequency analyses.
62
154
\item[DERTYPE = string]\mbox{}\\
63
This keyword (also described earlier) must be set to {\tt NONE} is only
64
energies are available for the chosen method and {\tt FIRST} if analytic
65
gradients are available.
66
\item[BFGS\_USE\_LAST = integer]\mbox{}\\ This keyword is used to specify
67
the number of gradient step for the BFGS update of the Hessian. The default
155
This keyword may be set to {\tt NONE} if only energies are available
156
for the chosen method, or {\tt FIRST} if analytic gradients are available.
157
\item[POINTS = 3 or 5]\mbox{}\\
158
Specifies whether frequencies are determined by a 3-point or a 5-point
159
formula of gradient differences. If only energy points are used, more
160
displacements are required, but the effect of this keyword in terms of
161
accuracy is the same.
165
Note: In some situations, vibrational frequency analysis via finite
166
differences may fail if the full point group symmetry is specified via
167
the {\tt symmetry} keyword. This happens because the user-given
168
{\tt symmetry} value can become incompatible with the actual symmetry
169
of the molecule when energies or gradients are evaluated for
170
symmetry-lowering displacements. In such situations, the user is
171
advised to let the program determine the symmetry automatically, rather
172
than specifying {\tt symmetry} manually. Otherwise, an error such as the
173
following may result:
177
error: problem assigning number of operations per class
178
*** stopping execution ***
181
The manual pages for the \PSInormco\ and \PSIintder\ modules contain
182
information on additional tools useful in vibrational frequency analysis
183
and coordinate transformation.