Reactive Scattering with Row-Orthonormal Hyperspherical Coordinates. 1. Transformation Properties and Hamiltonian for Triatomic Systems

Computation and analysis of bound vibrational spectra of the neon tetramer using row orthonormal hyperspherical coordinates

This paper presents the first implementation of the row-orthonormal hyperspherical coordinate formalism for the computation of the vibrational spectrum of a tetratomic system. The wavefunction of Ne₄ is expanded on a large basis set of hyperspherical harmonics generated numerically. This method not only provides spectra with reasonable accuracy, but also gives physical insight into the vibrational dynamics of the system. The characteristics of the spectra are related to the symmetry and localization of the wavefunction in configuration space.

Roles of dynamical symmetry breaking in driving oblate-prolate transitions of atomic clusters

This paper explores the driving mechanisms for structural transitions of atomic clusters between oblate and prolate isomers. We employ the hyperspherical coordinates to investigate structural dynamics of a seven-atom cluster at a coarse-grained level in terms of the dynamics of three gyration radii and three principal axes, which characterize overall mass distributions of the cluster. Dynamics of gyration radii is governed by two kinds of forces. One is the potential force originating from the interactions between atoms. The other is the dynamical forces called the internal centrifugal forces, which originate from twisting and shearing motions of the system. The internal centrifugal force arising from twisting motions has an effect of breaking the symmetry between two gyration radii. As a result, in an oblate isomer, activation of the internal centrifugal force that has the effect of breaking the symmetry between the two largest gyration radii is crucial in triggering structural transitions into prolate isomers. In a prolate isomer, on the other hand, activation of the internal centrifugal force that has the effect of breaking the symmetry between the two smallest gyration radii is crucial in triggering structural transitions into oblate isomers. Activation of a twisting motion that switches the movement patterns of three principal axes is also important for the onset of structural transitions between oblate and prolate isomers. Based on these trigger mechanisms, we finally show that selective activations of specific gyration radii and twisting motions, depending on the isomer of the cluster, can effectively induce structural transitions of the cluster. The results presented here could provide further insights into the control of molecular reactions.

Photoabsorption spectrum of helium trimer cation--theoretical modeling

The photoabsorption spectrum of He3(+) is calculated for two semiempirical models of intracluster interactions and compared with available experimental data reported in the middle UV range [H. Haberland and B. von Issendorff, J. Chem. Phys. 102, 8773 (1995)]. Nuclear delocalization effects are investigated via several approaches comprising quantum samplings using either exact or approximate (harmonic) nuclear wavefunctions, as well as classical samplings based on the Monte Carlo methodology. Good agreement with the experiment is achieved for the model by Knowles et al., [Mol. Phys. 85, 243 (1995); Mol. Phys. 87, 827 (1996)] whereas the model by Calvo et al., [J. Chem. Phys. 135, 124308 (2011)] exhibits non-negligible deviations from the experiment. Predictions of far UV absorption spectrum of He3(+), for which no experimental data are presently available, are reported for both models and compared to each other as well as to the photoabsorption spectrum of He2(+). A simple semiempirical point-charge approximation for calculating transition probabilities is shown to perform well for He3(+).

Theoretical methods for ultrafast spectroscopy

Time-resolved spectroscopy in the femtosecond and attosecond time domain is a tool to unravel the dynamics of nuclear and electronic motion in molecular systems. Theoretical insight into the underlying physical processes is ideally gained by solving the time-dependent Schrödinger equation. In this work, methods currently used to solve this equation are reviewed in a compact presentation. These methods involve numerical representations of wavefunctions and operators, the calculation of time evolution operators, the setting up of the Hamiltonian operators and the types of coordinates to be used hereto. The advantages and disadvantages of some methods are discussed.

Vibrational spectrum of Ar3(+) and relative importance of linear and perpendicular isomers in its photodissociation

The photodissociation dynamics of the argon ionized trimer Ar(3)(+) is revisited in the light of recent experimental results of Lepère et al. [J. Chem. Phys. 134, 194301 (2009)], which show that the fragment with little kinetic energy is always a neutral one, thus the available energy is shared by a neutral and ionic fragments as in Ar(2)(+). We show that these results can be interpreted as the photodissociation of the linear isomer of the system. We perform a 3D quantum computation of the vibrational spectrum of the system and study the relative populations of the linear (trimer-core) and perpendicular (dimer-core) isomers. We then show that the charge initially located on the central atom in the ground electronic state of the linear isomer migrates toward the extreme ones in the photoexcitation process such that photodissociation of the linear isomer produces a neutral central atom at rest in agreement with measured product state distributions.

The H(3) (+) rovibrational spectrum revisited with a global electronic potential energy surface

In this paper, we have computed the rovibrational spectrum of the H(3) (+) molecule using a new global potential energy surface, invariant under all permutations of the nuclei, that includes the long range electrostatic interactions analytically. The energy levels are obtained by a variational calculation using hyperspherical coordinates. From the comparison with available experimental results for low lying levels, we conclude that our accuracy is of the order of 0.1 cm(-1) for states localized in the vicinity of equilateral triangular configurations of the nuclei, and changes to the order of 1 cm(-1) when the system is distorted away from equilateral configurations. Full rovibrational spectra up to the H(+)+H(2) dissociation energy limit have been computed. The statistical properties of this spectrum (nearest neighbor distribution and spectral rigidity) show the quantum signature of classical chaos and are consistent with random matrix theory. On the other hand, the correlation function, even when convoluted with a smoothing function, exhibits oscillations which are not described by random matrix theory. We discuss a possible similarity between these oscillations and the ones observed experimentally.

Modeling of HeN+ clusters. II. Calculation of He3+ vibrational spectrum

We have computed the vibrational spectrum of the helium ionized trimer He(3)(+) using three different potential energy surfaces [D. T. Chang and G. L. Gellene, J. Chem. Phys. 119, 4694 (2003); E. Scifoni et al., ibid. 125, 164304 (2006); I. Paidarova et al., Chem. Phys. 342, 64 (2007)]. Differences in the details of these potential energy surfaces induce discrepancies between bound state energies of the order of 0.01 eV. The effects of the geometric phase induced by the conical intersection between the ground electronic potential energy surface and the first excited one are studied by computing vibrational spectra with and without this phase. The six lowest vibrational bound states are negligibly affected by the geometric phase. Indeed, they correspond to wavefunctions localized in the vicinity of the linear symmetric configurations and can be assigned well defined vibrational quantum numbers. On the other hand, higher excited states are delocalized, cannot be assigned definite vibrational quantum numbers, and the geometric phase shifts their energies by approximately 0.005 eV.

Calculation of argon trimer rovibrational spectrum

Rovibrational spectra of Ar3 are computed for total angular momenta up to J=6 using row-orthonormal hyperspherical coordinates and an expansion of the wave function on hyperspherical harmonics. The sensitivity of the spectra to the two-body potential and to the three-body corrections is analyzed. First, the best available semiempirical pair potential (HFDID1) is compared with our recent ab initio two-body potential. The ab initio vibrational energies are typically 1-2 cm-1 higher than the semiempirical ones, which is related to the slightly larger dissociation energy of the semiempirical potential. Then, the Axilrod-Teller asymptotic expansion of the three-body correction is compared with our newly developed ab initio three-body potential. The difference is found smaller than 0.3 cm-1. In addition, we define approximate quantum numbers to describe the vibration and rotation of the system. The vibration is represented by a hyper-radial mode and a two-degree-of-freedom hyperangular mode, including a vibrational angular momentum defined in an Eckart frame. The rotation is described by the total angular momentum quantum number, its projection on the axis perpendicular to the molecular plane, and a hyperangular internal momentum quantum number, related to the vibrational angular momentum by a transformation between Eckart and principal-axes-of-inertia frames. These quantum numbers provide a qualitative understanding of the spectra and, in particular, of the impact of the nuclear permutational symmetry of the system (bosonic with zero nuclear spin). Rotational constants are extracted from the spectra and are shown to be accurate only for the ground hyperangular mode.

Gyration-radius dynamics in structural transitions of atomic clusters

This paper is concerned with the structural transition dynamics of the six-atom Morse cluster with zero total angular momentum, which serves as an illustrative example of the general reaction dynamics of isolated polyatomic molecules. It develops a methodology that highlights the interplay between the effects of the potential energy topography and those of the intrinsic geometry of the molecular internal space. The method focuses on the dynamics of three coarse variables, the molecular gyration radii. By using the framework of geometric mechanics and hyperspherical coordinates, the internal motions of a molecule are described in terms of these three gyration radii and hyperangular modes. The gyration radii serve as slow collective variables, while the remaining hyperangular modes serve as rapidly oscillating "bath" modes. Internal equations of motion reveal that the gyration radii are subject to two different kinds of forces: One is the ordinary force that originates from the potential energy function of the system, while the other is an internal centrifugal force. The latter originates from the dynamical coupling of the gyration radii with the hyperangular modes. The effects of these two forces often counteract each other: The potential force generally works to keep the internal mass distribution of the system compact and symmetric, while the internal centrifugal force works to inflate and elongate it. Averaged fields of these two forces are calculated numerically along a reaction path for the structural transition of the molecule in the three-dimensional space of gyration radii. By integrating the sum of these two force fields along the reaction path, an effective energy curve is deduced, which quantifies the gross work necessary for the system to change its mass distribution along the reaction path. This effective energy curve elucidates the energy-dependent switching of the structural preference between symmetric and asymmetric conformations. The present methodology should be of wide use for the systematic reduction of dimensionality as well as for the identification of kinematic barriers associated with the rearrangement of mass distribution in a variety of molecular reaction dynamics in vacuum.

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.

Dynamical and statistical effects of the intrinsic curvature of internal space of molecules

The Hamilton dynamics of a molecule in a translationally and/or rotationally symmetric field is kept rigorously constrained in its phase space. The relevant dynamical laws should therefore be extracted from these constrained motions. An internal space that is induced by a projection of such a limited phase space onto configuration space is an intrinsically curved space even for a system of zero total angular momentum. In this paper we discuss the general effects of this curvedness on dynamics and structures of molecules in such a manner that is invariant with respect to the selection of coordinates. It is shown that the regular coordinate originally defined by Riemann is particularly useful to expose the curvature correction to the dynamics and statistical properties of molecules. These effects are significant both qualitatively and quantitatively and are studied in two aspects. One is the direct effect on dynamics: A trajectory receives a Lorentz-like force from the curved space as though it was placed in a magnetic field. The well-known problem of the trapping phenomenon at the transition state is analyzed from this point of view. By showing that the trapping force is explicitly described in terms of the curvature of the internal space, we clarify that the physical origin of the trapped motion is indeed originated from the curvature of the internal space and hence is not dependent of the selection of coordinate system. The other aspect is the effect of phase space volume arising from the curvedness: We formulate a general expression of the curvature correction of the classical density of states and extract its physical significance in the molecular geometry along with reaction rate in terms of the scalar curvature and volume loss (gain) due to the curvature. The transition state theory is reformulated from this point of view and it is applied to the structural transition of linear chain molecules in the so-called dihedral angle model. It is shown that the curvature effect becomes large roughly linearly with the size of molecule.