Quasi-Classical Trajectory Study of Atom-Diatomic Molecule Collisions in Symmetric Hyperspherical Coordinates: The F + HCl Reaction as a Test Case

We investigate the reactive dynamics of the triatomic system F + HCl → HF + Cl for total angular momentum equal zero and for different low-lying rovibrational states of the diatomic molecule. For each of the initial vibrational quantum numbers, the time evolution of the atom-diatom collision process is investigated for a wide range of impact angles and collision energies. To this purpose, the Quasi-Classical Trajectories (QCT) method was implemented in a hyperspherical configuration space. The Hamilton equations of motion are solved numerically in an intermediate effective Cartesian space to exploit the relative simplicity of this intermediate representation. Interatomic interactions are described by a London-Eyring-Polanyi-Sato potential energy surface, specifically developed for the title reaction, and the results of the QCT simulations are discussed in terms of the time-evolution of the hyperangles. The analysis of the collision dynamics using symmetric hyperspherical coordinates provides, in addition to the description in terms of a natural reaction coordinate (the hyperradius), a more striking representation of the exchange dynamics, in terms of the time-dependent probability distribution along the kinematic rotation hyperangle, and a precise distinction between direct and indirect mechanisms of the reaction.

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.

Incorporating the Geometric Phase Effect in Triatomic and Tetraatomic Hyperspherical Harmonics

Hyperspherical harmonics in the democratic row-orthonormal hyperspherical coordinates are very appropriate basis sets for performing reactive scattering calculations for triatomic and tetraatomic systems. The mathematical conditions for incorporating the geometric phase effect in these harmonics are given. These conditions are implemented for triatomic systems, and their explicit analytical expressions in terms of Jacobi polynomials, in both the absence and presence of the geometric phase effect, are given.