Hyperspherical Harmonics for Triatomic Systems

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(+).

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.

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.

Phase-space invariants for aggregates of particles: Hyperangular momenta and partitions of the classical kinetic energy

Rigorous definitions are presented for the kinematic angular momentum K of a system of classical particles (a concept dual to the conventional angular momentum J), the angular momentum L(xi) associated with the moments of inertia, and the contributions to the total kinetic energy of the system from various modes of the motion of the particles. Some key properties of these quantities are described-in particular, their invariance under any orthogonal coordinate transformation and the inequalities they are subject to. The main mathematical tool exploited is the singular value decomposition of rectangular matrices and its differentiation with respect to a parameter. The quantities introduced employ as ingredients particle coordinates and momenta, commonly available in classical trajectory studies of chemical reactions and in molecular dynamics simulations, and thus are of prospective use as sensitive and immediately calculated indicators of phase transitions, isomerizations, onsets of chaotic behavior, and other dynamical critical phenomena in classical microaggregates, such as nanoscale clusters.

Three-body problem in quantum mechanics: Hyperspherical elliptic coordinates and harmonic basis sets

Elliptic coordinates within the hyperspherical formalism for three-body problems were proposed some time ago [V. Aquilanti, S. Cavalli, and G. Grossi, J. Chem. Phys. 85, 1362 (1986)] and recently have also found application, for example, in chemical reaction theory [see O. I. Tolstikhin and H. Nakamura, J. Chem. Phys. 108, 8899 (1998)]. Here we consider their role in providing a smooth transition between the known "symmetric" and "asymmetric" parametrizations, and focus on the corresponding hyperspherical harmonics. These harmonics, which will be called hyperspherical elliptic, involve products of two associated Lame polynomials. We will provide an expansion of these new sets in a finite series of standard hyperspherical harmonics, producing a powerful tool for future applications in the field of scattering and bound-state quantum-mechanical three-body problems.