Frames and vibrational coordinates¶
Here we introduce different ingredients available for triatomic molecules, including
Molecular frames ;
() vibrational coordinates ;
For the linearised coordinates, the default frame is Eckart. The equilibrium structures, required for the definition of the linearised KEO and PEF, are chosen as the principal axis system (PAS). For the curvilinear KEOs, the frames are defined by the construction of the KEOs in their analytic representations.
Triatomics¶
XY2 type molecules¶
A molecule type is defined by the keyword MolType
. For the XY2 example it is
MolType XY2
in the curvilinear KEO, it is common in TROVE to use the bisector frame for the XY2 molecules, with the axis bisecting the bond angle and the in the plane of the molecule, but other embeddings are possible. The PAS frame coincides with the bisector frame at the equilibrium or non-rigid reference configuration (i.e. symmetric). In TROVE, the definition of the frame is combined with the definition of the internal coordinates via the keyword transform
. In the following, these are described.
There are currently at least two exact, curvilinear KEO forms are provided for a quasi-linear XY2 molecules, MLkinetic_xy2_bisect_EKE
, MLkinetic_xy2_bisect_EKE_sinrho
, see below.
R-RHO-Z
¶
R-RHO-Z
is used for (quasi-)linear molecules of the XY2 type. it defined the curvilinear vibrational coordinates as the two bond angles and with the bending mode described by the angle , where is the interbond angle (). For the rigid reference frame (REFER-CONF RIGID
), the actual internal coordinates are the displacements of , and from the corresponding equilibrium:
where is the equilibrium bond length. If the non-rigid reference frame is used (REFER-CONF NON-RIGID
), the bending mode is given on an equidistant grid, typically of 1000-2000 points, while the stretching modes are the displacements from the given point along the non-rigid reference frame, the latter is usually defined as the principal axes system with the bond length fixed to the equilibrium:
TROVE uses Z-matrix coordinates to build any user-defined coordinates. In this case, the Z=matrix is given by
ZMAT
S 0 0 0 0 31.97207070
H 1 0 0 0 1.00782505
H 1 2 0 0 1.00782505
end
Alternatively, the reference value of the bond length can also vary with as e.g. in the minimum energy path (MEP) definition with being the optimised value at the given value of corresponding to the local energy minimum. In this case, the non-rigid frame must be defined using the MEP
block (see the corresponding section).
For the linearised coordinates type (COORDS Linear
), the actual internal coordinates are the linearised versions of above. More specifically, for the non-rigid
reference configuration, the bending coordinate :math`rho` is kept curvilinear on a grid of :math`rho_k` points as before, while the stretching coordinates are defined by linearly expanding :math`r_1` and :math`r_2` in terms of the Cartesian displacement around the corresponding reference values . In the Rigid
case, the bending coordinate is also linearised.
The advantage of the linearised coordinates is that the corresponding KEO can be constructed on the fly as part of the TROVE generalised procedure as a Taylor type expansion. The main disadvantage however is that the approximate linearised KEO operator is less accurate than the (exact) curvilinear EKO. Besides, the convergence of the variational solution is also poorer for the linearised case (see [15YaYu]).
R-RHO-Z-ECKART
¶
This Transform
type is very similar to R-RHO-Z
, but with the molecular frame define using the Eckart conditions.
R-ALPHA-Z
¶
R-ALPHA-Z
is very similar toR-RHO-Z
with the difference in the bending coordinate, which in the interbond angle in this case. In the ``Rigid reference configuration, it is a displacement from the equilibrium value :
In the Non-rigid
reference configuration, is given on a grid of points ranging from to and including the equilibrium value. In the linearised Rigid
case, the bending coordinated is defined as a linear expansion of at in terms of the Cartesian displacements.
TROVE input example:
COORDS local (curvilinear coordinates) TRANSFORM r-rho-z (r1, r2, rho with the x parallel to the bisector) MOLTYPE XY2 REFER-CONF non-RIGID (Reference configuration)
Note
The text in brackets is used for comments.
R-RHO-Z-M2-M3¶
A ‘bisecting’ XY2 frame used for isotopologies with slightly different masses of Y1 and Y2, for example O16CO17. Although this is an XYZ molecule, in this case it is formally treated as XY:sub:`2 but with non-symmetric masses and the Cs symmetry, e.g.:
TRANSFORM R-RHO-Z-M2-M3
MOLTYPE XY2
MOLECULE CO2
REFER-CONF non-RIGID
SYMGROUP Cs(M)
ZMAT
C 0 0 0 0 11.996709
O 1 0 0 0 16.995245
O 1 2 0 0 15.9905256
end
XYZ type molecules¶
The main embedding here is the ‘bond’-embedding, with the axis placed parallel to the bond Y-Z with a heavier atom Z comparing to X (second bond).
For molecules XYZ with comparable masses X and Z (e.g. in similar isotopologues), the bisector frames and associated TRANSFORM
can be used.
R1-Z-R2-RHO
¶
This is a ‘bond’-embedding with the same vibrational coordinates as in R-RHO-Z
and along the axis. The coordinates are givem as above:
Here is an example of the Z-matrix for NNO.
ZMAT
N 0 0 0 0 14.00307401
N 1 0 0 0 14.00307401
O 1 2 0 0 15.994915
end
R1-Z-R2-ALPHA
¶
This is another ‘bond’-embedding with the same vibrational coordinates as in R-ALPHA-Z
.
Tetratomics¶
XY3 rigid molecules (PH3 type)¶
Linearized KEOs use the Eckart frame with the PAS at the equilibrium configuration. The latter has the axis along the axis of symmetry with the axis chosen in plane containing the X-Y1 bond and passing through .
R-ALPHA
¶
For the rigid XY3, like PH3, the logical coordinate choice of the valence coordinates consists of three bond lengths , , , , and . For the linearised KEO, these valence are used to form the linearised coordinates in the same way as before (1st order expansion in terms of the Cartesian displacement). For the curvilinear KEO (local
), the vibrational coordinates are then defined as displacement from the corresponding equilibrium (or non-rigid reference) values:
The underlying Z-matrix coordinates are defined using the following Z-matrix:
ZMAT
N 0 0 0 0 14.00307401
H 1 0 0 0 1.00782505
H 1 2 0 0 1.00782505
H 1 2 3 1 1.00782505
end
This representation has been used for PH3 [15SoAlTe], SbH3 [10YuCaYa], AsH3 [19CoYuKo], PF3 [19MaChYa].
XY3 non-rigid with umbrella motion (NH3 type)¶
MolType XY3
Consider the Ammonia molecule NH33 with a relatively small barrier to the planarity. The three bending angles are not suitable in this case as they cannot distinguish the two opposite inversion configurations above and below the planarity. Instead, an umbrella mode has to be introduced as one of the bending modes. An example of an umbrella coordinate is an angle between the symmetry axis and the bond X-Y, see Figure. It is natural to use the non-rigid reference configuration along the umbrella, inversion motion and build the KEO as an expansion around it. For two other bending modes, in principle one can use two inter-bond angles, e.g. and , two dihedral angles and . However, for symmetry reasons, TROVE employs the symmetry-adapted bending pair and , defined as follows:
or
The umbrella mode for any instantaneous configuration of the nuclei is defined in TROVE as the angle between a trisector
Linearized KEOs use the Eckart frame with the PAS at the equilibrium configuration. The latter has the axis along the axis of symmetry with the axis chosen in plane containing the X-Y1 bond and passing through .
R-S-DELTA
¶
For this TRANSFORM
case, the following valence-based coordinates are used:
The umbrella mode :math:\delta
is defined as an angle between the trisector and any of the bonds X-Y. The other 5 coordinates are then used to construct the corresponding linearised vibrational coordinates (see above) for the linearised (linear
) representation.
ZXY2 (Formaldehyde type)¶
MolType ZXY2
The common valence coordinate choice for ZXY2 includes three bond lengths , two bond angles and a dihedral angle . The latter can be treated as the reference for a non-rigid reference configuration in TROVE on a grid of ranging from :math`[-tau_{0}ldots tau_{0}]`, while other 5 modes are treated as displacement from their equilibrium values at each grid point . The reference configuration is always in the principle axis sysetm, i.e. for each value of the book angle , TROVE solve the PAS conditions to reorient the molecule.
R-THETA-TAU
¶
Isotopologues of XY3 as ZXY2 type¶
The Z type can be used to define signle or double deturated isotopologues of an XY3 molecule such as a rigid PH3 or non-rigid NH3. For PDH2, we use R-THETA-TAU
in combination with the Z-matrix given as follows:
ZMAT
P 0 0 0 0 14.00307401
D 1 0 0 0 2.01410178
H 1 2 0 0 1.007825032
H 1 2 3 2 1.007825032
end
Here, the equilibrium frame coinsides with the principle axis system with the axis in the plane conteining PD and bisetcing the angle betwen two PH bonds.
For a PH:sub:`2`D type isotopologue, the Z-matrix is given by
- ZMAT
P 0 0 0 0 14.00307401 H 1 0 0 0 1.007825032 D 1 2 0 0 2.01410178 D 1 2 3 2 2.01410178
end
ZXY3 (Methyl Chloride type)¶
MolType ZXY3
Similarilly, for the ZXY3 type molecule we use valence coordinates consisting of four bond lengths , (), three bond angles and two symmetry adapted dihedral coordinates constructed from three dihedral angles , where . This is a rigid
type, where all coordinates are treated as displacements from the corresponding equilibrium values. Currently, only the standard linearised KEO is available in TROVE.
R-BETA-SYM
¶
The Z-matrix coordinates (underlying basic TROVE coordinates) are as given by the Z-matrix rules:
ZMAT
C 0 0 0 0 12.000000000
Cl 1 0 0 0 34.968852721
H 1 2 0 0 1.007825035
H 1 2 3 0 1.007825035
H 1 2 3 4 1.007825035
end
are as follows:
,
, ,
, ,
where alpha_{12}` and are interbond angles between the bonds X-Yi. The Z-matrix coordinates are transformed to :math:`tau_{12}, tau_{23}, tau_{13} ` via the following trigonometric rules:
A chain ABCD type molecule (hydrogen peroxide type)¶
MolType ABCD
R-ALPHA-TAU
¶
The six internal coordinates for the Transform R-ALPHA-TAU
type consist of three stretching, two bending and one dihedral coordinates as given by
The non-rigid reference frame such that the axis bisects the dihedral angle.
For this embedding, in order to be able to separate the vibrational and rotational bases into a product form, it is important to the an extended range for the dihedral angle . Otherwise the eigenfunction is obtained double valued due to the axis appearing in the opposite direction to the two bonds after one revolution.
A minimum energy path (MEP) as a non-rigid reference configuration¶
In MEP, the 5 internal coordinate displacements are defined around -dependent reference values. The latter are obtained as oprmised geometries by minimised molecule’s energy:
where :math: the MEP values are given by a parameterised expansion, for example
where .
Fuive-atomic molecules¶
The XY4 molecule (Td) and the XY4
type¶
MolType XY4
The frame for the tetrahedral molecule XY4 spanning the Td(M) symmetry group is chosen with the axes orthogonal to the faces of the box containing the molecule with the four atoms at its vertices, as shown in the figure, with the Cartesian coordinates at equilibrium given by
R-ALPHA
¶
The tetrahedral five-atomic molecule XY4 has 9 vibrational degrees of freedom. For a semi-rigid molecule (i.e. ignoring any isomerisation that can occur at higher energies), they can be characterised by four bond lengths and six inter-bond angles . For the equilibrium value of the tetrahedral angle , = which explains the factor in the definition of the Cartesian coordinates. There should, however, be only 9 independent vibrational degrees of freedom in a 5 atomic molecule. One of the inter-bond angles is redundant as there should be only five independent bending vibrations, with the following redundancy condition:
(1)¶
XY4 belongs to the Td(M) molecular symmetry group, which consists of five irreducible representations, , , , and . One way to define independent bending modes is to reduce the six inter-bond angles to five symmetry-adapted irreducible combinations, which, together with four bond lengths form nine independent vibrational modes as follows: four stretches
(2)¶
two -symmetry bends
(3)¶
and three -symmetry bends
(4)¶
where the corresponding symmetries of the bending modes are indicated.
The stretching modes can also be in principle combined into symmetry-adapted coordinates in Td(M):
(5)¶
Six-atomic molecules¶
The C2H:sub:4 molecule and the C2H4
type¶
MolType C2H4
C2H4_2BETA_1TAU
¶
The internal coordinates are defined using the following 12 valence coordinates: 5 stretching (molecular bond) coordinates, 4 bending (inter-bond angles) and 3 dihedral coordinates, with the last mode as a book angle describing the relative motion of two moieties:
where
This type can be used both for rigid and non-rigid molecule types. The non-rigid coordinate is in the latter case.