Separation of variables and search algorithm for internal rotation coordinates and kinematic coefficients in the theory of nuclear motions in polyatomic molecules

The construction algorithm is proposed for the internal rotation coordinates in polyatomic molecules. It is based on the properties of the matrix of kinematic coefficients when an excessive system of natural coordinates is introduced. The approximations providing the separation of variables are considered. The exact form of the kinetic energy operator is given.     

Large amplitude vibration in polyatomic molecules. I. A polar representation of orthogonal relative coordinates

A general expression for the rovib kinetic energy operator in a system of orthogonal internal coordinates is derived for arbitrary amplitudes of motion. Plots of the molecular potential in such coordinates serve to determine coordinate subsets in which the potential approaches separability. A revision of previous theory for the three-body problem is presented. A rovib hamiltonian for non-branched four-atom molecules (such as C₂H₂) is derived. Certain features of application to polyatomic (five-atom) molecules are considered.     

Expressing kinetic‐energy for vibration–rotation motions in general rotating systems of axes and introducing quasirectilinear vibrational coordinates to simplify Hamiltonian forms

To understand how the internal and rotational motions of a polyatomic system depend on which rotating system of axes is selected, we derived the explicit form of the atomic velocities determined by an observer stationed on the general rotating system of axes. Using the derived velocities, we formulated the kinetic energy expression for vibration–rotation motions with respect to the rotating system of axes. From this expression, we clarified covariant metric tensors under zero angular momentum, which have been confused with an erroneous expression even in the professional literature, and the relationship between the kinetic energy expression and the rotating system of axes. Furthermore, to simplify the Hamiltonian form, we introduced quasirectilinear vibrational coordinates to describe the Hamiltonian. The resulting Hamiltonian form is superior to those of the previous studies in that the kinetic and potential energy expressions are simple and the vibrational frequencies are independent of the original internal coordinates used. In fact, we show that its application for three examples is useful. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 83: 22–29, 2001     

Natural coordinates for polyatomic reactions

The curvilinear natural coordinates (NRC) are considered to describe the polyatomic reaction dynamics. An invariant equation is found for the reaction path curve (RP), allowing its calculation in an arbitrary initial coordinate system. Apart from common properties, well known for triatomic reactions, the polyatomic NRC take explicit account of a smooth change of normal vibrational modes of a chemical system during its evolution along RP. This generates special terms in the kinetic energy operator that induce nonadiabatic vibrational transitions and do not vanish in a classical limit, no matter how small the curvature of the RP is. For the case when the amplitudes of vibrations perpendicular to the RP are small it is shown that, as a first approximation, the vibrational transitions can be treated independently of the rotational degrees of freedom. This conclusion is obtained as a generalization of Eckart theory stating the possibility of vibration—rotational separation for stable molecules.     

Rovibrational energy levels of hydrogen peroxide, studied by MULTIMODE with a reaction path Hamiltonian

Recently, Carter and Handy [J. Chem. Phys. 113 (2000) 987] have introduced the theory of the reaction path Hamiltonian (RPH) [J. Chem. Phys. 72 (1980) 99] into the variational scheme MULTIMODE, for the calculation of the J=0 vibrational levels of polyatomic molecules, which have a single large-amplitude motion. In this theory the reaction path coordinate s becomes the fourth dimension of the moment-of-inertia tensor, and must be treated separately from the remaining 3N-7 normal coordinates in the MULTIMODE program. In the modified program, complete integration is performed over s, and the M-mode MULTIMODE coupling approximation for the evaluation of the matrix elements applies only to the 3N-7 normal coordinates. In this paper the new algorithm is extended to the calculation of rotational-vibration energy levels (i.e. J>0) with the RPH, following from our analogous implementation for rigid molecules [Theoret. Chem. Acc. 100 (1998) 191]. The full theory is given, and all extra terms have been included to give the exact kinetic energy operator. In order to validate the new code, we report studies on hydrogen peroxide (H2O2), where the reaction path is equivalent to torsional motion. H2O2 has previously been studied variationally using a valence coordinate Hamiltonian; complete agreement for calculated rovibrational levels is obtained between the previous results and those from the new code, using the identical potential surface. MULTIMODE is now able to calculate rovibrational levels for polyatomic molecules which have one large-amplitude motion.     

A new method for the derivation of exact vibration-rotational kinetic energy operator in internal coordinates

An efficient angular momentum method is presented and used to derive analytic expressions for the vibration-rotational kinetic energy operator of polyatomic molecules.The vibration-rotational kinetic energy operator is expressed in terms of the total angular momentum operator J,the angular momentum operator J and the momentum operator p conjugate to Z in the molecule-fixed frame Not only the method of derivation is simpler than that in the previous work,but also the expressions ot the kinetic energy operators arc more compact.Particularly,the operator is easily applied to different vibrational or rovibrational problems of the polyatomic molecules by variations of matrix elements Gn of a mass-dependent constant symmetric matrix     

A new exact quantum mechanical rovibrational kinetic energy operator for penta-atomic systems in internal coordinates

The concrete molecule-fixed (MF) kinetic energy operator for penta-atomic molecules is expressed in terms of the parameterδ, the matrix element G_(?), and angular momentum operator (?). The applications of the operator are also discussed. Finally, a general compact form of kinetic energy operator suitable for calculating the rovibrational spectra of polyatomie molecules is presented.     

Algebraic methods in molecular structure

The use of algebraic methods in molecular structure is briefly reviewed. The fundamental algebra, ?? ≡ U(4), of rotation–vibration spectra is introduced and its implications on spectra of di-, tri-, and polyatomic molecules are discussed.     

Franck–Condon factors in curvilinear coordinates: the photoelectron spectrum of ammonia

An approach to the calculation of Franck–Condon factors in curvilinear coordinates is outlined. The approach is based on curvilinear normal coordinates, which allows for an easy extension of Duschinsky’s transformation to the case of curvilinear coordinates, and on the power series expansion of the kinetic energy operator. Its usefulness in the case of molecules undergoing large displacements of their equilibrium nuclear configurations upon excitation is then demonstrated by an application to the vibrational structure of the photoelectron spectrum of ammonia, using an anharmonic potential only for the symmetric stretching and bending coordinates of the radical cation.     

Vibrational relaxation induced population inversions in laser pumped polyatomic molecules

Conditions for population inversion in laser pumped polyatomic molecules are described. For systems which exhibit metastable vibrational population distributions [slow vibration—translation/rotation (V—T/R) relaxation], large, long lived inversions are possible even when the vibrational modes are strongly coupled by rapid collisional vibration—vibration (V—V) energy transfer. Overtone states of a hot mode are found to invert with respect to fundamental levels of a cold mode even at V—V steady state. Inversion persists for a V—T/R relaxation time. A gain of 4 m^?1 for the 2v₃ → v₂ transition in CH₃F (λ ≈ 15.9 μ) was found assuming a spontaneous emission lifetime of 10 s for this transition. General equations are derived which can be used to determine the magnitude of population inversion in any laser pumped, vibrationally metastable, polyatomic molecule. A discussion of factors controlling the population maxima of different vibrational states in optically pumped, V—V equilibrated metastable polyatomics is also given.     

Automatic computer procedure for generating exact and analytical kinetic energy operators based on the polyspherical approach

We develop a new general code to automatically derive exact analytical kinetic energy operators in terms of polyspherical coordinates. Computer procedures based on symbolic calculations are implemented. Sets of orthogonal or non-orthogonal vectors are used to parametrize the molecular systems in space. For each set of vectors, and whatever the size of the system, the exact analytical kinetic energy operator (including the overall rotation and the Coriolis coupling) can be derived by the program. The correctness of the implementation is tested for different sets of vectors and for several systems of various sizes.     

Global analysis of composite particles

The theory of vibrations of a composite particle when vibrational amplitudes are not constrained to be small according to the Eckart conditions is developed using the methods of differential topology. A global classical Hamiltonian appropriate for this system is given, and for the case of the molecular vibration–rotation problem, it is transformed into a global quantum Hamiltonian operator. It is shown that the zeroth-order term in the global Hamiltonian operator is identical to the Wilson–Howard Hamiltonian; higher-order terms are shown to give successively better approximations to the large amplitude problem. Generalized Eckart conditions are derived for the global classical Hamiltonian; the quantum equivalent of these conditions along with the quantum equivalent of the Eckart conditions are given. The spectrum of the global Hamiltonian operator is discussed and it is shown that the calculation of the vibration–rotation energy states of the system reduces to the same straight-forward procedure, the solution of a secular determinant, as was carried out for the Wilson–Howard Hamiltonian at a later time by Nielsen.     

Numerical generation of hyperspherical harmonics for tetra-atomic systems

A numerical generation method of hyperspherical harmonics for tetra-atomic systems, in terms of row-orthonormal hyperspherical coordinates-a hyper-radius and eight angles-is presented. The nine-dimensional coordinate space is split into three three-dimensional spaces, the physical rotation, kinematic rotation, and kinematic invariant spaces. The eight-angle principal-axes-of-inertia hyperspherical harmonics are expanded in Wigner rotation matrices for the physical and kinematic rotation angles. The remaining two-angle harmonics defined in kinematic invariant space are expanded in a basis of trigonometric functions, and the diagonalization of the kinetic energy operator in this basis provides highly accurate harmonics. This trigonometric basis is chosen to provide a mathematically exact and finite expansion for the harmonics. Individually, each basis function does not satisfy appropriate boundary conditions at the poles of the kinetic energy operator; however, the numerically generated linear combination of these functions which constitutes the harmonic does. The size of this basis is minimized using the symmetries of the system, in particular, internal symmetries, involving different sets of coordinates in nine-dimensional space corresponding to the same physical configuration.     

位移谐振子模型对三(2-苯基吡啶)合铱磷光光谱影响的理论研究

采用谐振子模型理论探讨了振动模式对Ir(ppy)₃配合物的磷光光谱的影响.多原子分子发射光谱的一般形式可以从两个绝热电子态之间的热振动关联函数推导出,相应地势能面之间的位移和Duschinsky转动的影响也被包含在多维谐振子模型的表达式中,所得关系式模拟出了Ir(ppy)₃较为精细的磷光发射光谱.计算结果表明T₁态到S₀态之间的0→1振动跃迁对发射光谱贡献较大,尤其振动频率小于1600 cm^-1的振动模贡献更多,配体中苯和吡啶环上C=C和C=N的呼吸振动,是Ir(ppy)₃出现肩峰的主要原因.玻耳兹曼分布使得主峰和肩峰的强度下降,并且两峰相互接近.该谐振子模型与密度泛函理论(DFT)结合,可以较好地定量描述多原子分子光物理过程的发射光谱以及详细了解光谱谱图的细节.     

A simple method for derivation of the rovibrational Hamiltonian: application to orthogonal radau coordinate systems as special cases

The molecular Hamiltonian of polyatomic molecules has been obtained. A general choice of internal coordinates depending on external parameters was considered. The rovibrational Hamiltonian for this set of coordinate system was derived in general terms as a function of the external parameters a and b. This procedure is also applicable to various kinds of internal coordinates in a straightforward way. The rovibrational Hamiltonian of triatomic molecules is considered as an application of this general formulation. In addition, orthogonal Radau coordinates are considered as cases of this new approach     

Use of a single trajectory to study product energy partitioning in unimolecular dissociation: mass effects for halogenated alkanes

A single trajectory (ST) direct dynamics approach is compared with quasiclassical trajectory (QCT) direct dynamics calculations for determining product energy partitioning in unimolecular dissociation. Three comparisons are made by simulating C(2)H(5)F-->HF + C(2)H(4) product energy partitioning for the MP26-31G(*) and MP26-311 + + G(**) potential energy surfaces (PESs) and using the MP26-31G(*) PES for C(2)H(5)F dissociation as a model to simulate CHCl(2)CCl(3)-->HCl + C(2)Cl(4) dissociation and its product energy partitioning. The trajectories are initiated at the transition state with fixed energy in reaction-coordinate translation E(t) (double dagger). The QCT simulations have zero-point energy (ZPE) in the vibrational modes orthogonal to the reaction coordinate, while there is no ZPE for the STs. A semiquantitative agreement is obtained between the ST and QCT average percent product energy partitionings. The ST approach is used to study mass effects for product energy partitioning in HX(X = F or Cl) elimination from halogenated alkanes by using the MP26-31G(*) PES for C(2)H(5)F dissociation and varying the masses of the C, H, and F atoms. There is, at most, only a small mass effect for partitioning of energy to HX vibration and rotation. In contrast, there are substantial mass effects for partitioning to relative translation and the polyatomic product's vibration and rotation. If the center of mass of the polyatomic product is located away from the C atom from which HX recoils, the polyatomic has substantial rotation energy. Polyatomic products, with heavy atoms such as Cl atoms replacing the H atoms, receive substantial vibration energy that is primarily transferred to the wag-bend motions. For E(t) (double dagger) of 1.0 kcalmol, the ST calculations give average percent partitionings to relative translation, polyatomic vibration, polyatomic rotation, HX vibration, and HX rotation of 74.9%, 6.8%, 1.5%, 14.4%, and 2.4% for C(2)H(5)F dissociation and 39.7%, 38.1%, 0.2%, 16.1%, and 5.9% for a model of CHCl(2)CCl(3) dissociation.     

Hamiltonian for rovibrational spectra of “hot” molecules

A variational method based on the determination of rotation as a state with a constant angular momentum has been proposed for calculating rovibrational energy levels of a polyatomic molecule. By using this method, energy level values can be determined for any vibrational state with arbitrary values of vibrational quantum numbers. This makes it possible to calculate the rovibrational energy levels of the “hot” molecules.     

Quantum dynamics of dissociative chemisorption of CH(4) on Ni(111): Influence of the bending vibration

Two-dimensional, three-dimensional, and four-dimensional quantum dynamic calculations are performed on the dissociative chemisorption of CH(4) on Ni(111) using the multiconfiguration time-dependent Hartree (MCTDH) method. The potential energy surface used for these calculations is 15-dimensional (15D) and was obtained with density functional theory for points which are concentrated in the region that is dynamically relevant to reaction. Many reduced dimensionality calculations were already performed on this system, but the molecule was generally treated as pseudodiatomic. The main improvement of our model is that we try to describe CH(4) as a polyatomic molecule by including a degree of freedom describing a bending vibration in our three-dimensional and four-dimensional models. Using a polyspherical coordinate system, a general expression for the 15D kinetic energy operator is derived, which discards all the singularities in the operator and includes rotational and Coriolis coupling. We use seven rigid constraints to fix the CH(3) umbrella of the molecule to its gas phase equilibrium geometry and to derive two-dimensional, three-dimensional, and four-dimensional Hamiltonians, which were used in the MCTDH method. Only four degrees of freedom evolve strongly along the 15D minimum energy path: the distance of the center of mass of the molecule to the surface, the dissociative C[Single Bond]H bond distance, the polar orientation of the molecule, and the bending angle between the dissociative C[Single Bond]H bond and the umbrella. A selection of these coordinates is included in each of our models. The polar rotation is found to be important in determining the mode selective behavior of the reaction. Furthermore, our calculations are in good agreement with the finding of Xiang et al. [J. Chem. Phys. 117, 7698 (2002)] in their reduced dimensional calculation that the helicopter motion of the umbrella symmetry axis is less efficient than its cartwheel motion for promoting the reaction. The effect of pre-exciting the bend modes is qualitatively incorrect at higher energies, suggesting the necessity of including additional rotational and vibrational degrees of freedom in the model.     

Rotating full- and reduced-dimensional quantum chemical models of molecules

A flexible protocol, applicable to semirigid as well as floppy polyatomic systems, is developed for the variational solution of the rotational-vibrational Schro?dinger equation. The kinetic energy operator is expressed in terms of curvilinear coordinates, describing the internal motion, and rotational coordinates, characterizing the orientation of the frame fixed to the nonrigid body. Although the analytic form of the kinetic energy operator might be very complex, it does not need to be known a priori within this scheme as it is constructed automatically and numerically whenever needed. The internal coordinates can be chosen to best represent the system of interest and the body-fixed frame is not restricted to an embedding defined with respect to a single reference geometry. The features of the technique mentioned make it especially well suited to treat large-amplitude nuclear motions. Reduced-dimensional rovibrational models can be defined straightforwardly by introducing constraints on the generalized coordinates. In order to demonstrate the flexibility of the protocol and the associated computer code, the inversion-tunneling of the ammonia ((14)NH(3)) molecule is studied using one, two, three, four, and six active vibrational degrees of freedom, within both vibrational and rovibrational variational computations. For example, the one-dimensional inversion-tunneling model of ammonia is considered also for nonzero rotational angular momenta. It turns out to be difficult to significantly improve upon this simple model. Rotational-vibrational energy levels are presented for rotational angular momentum quantum numbers J = 0, 1, 2, 3, and 4.     

Hierarchical expansion of the kinetic energy operator in curvilinear coordinates for the vibrational self-consistent field method

A new hierarchical expansion of the kinetic energy operator in curvilinear coordinates is presented and modified vibrational self-consistent field (VSCF) equations are derived including all kinematic effects within the mean field approximation. The new concept for the kinetic energy operator is based on many-body expansions for all G matrix elements and its determinant. As a test application VSCF computations were performed on the H(2)O(2) molecule using an analytic potential (PCPSDE) and different hierarchical approximations for the kinetic energy operator. The results indicate that coordinate-dependent reduced masses account for the largest part of the kinetic energy. Neither kinematic couplings nor derivatives of the G matrix nor its determinant had significant effects on the VSCF energies. Only the zero-point value of the pseudopotential yields an offset to absolute energies which, however, is irrelevant for spectroscopic problems.