First-principles rovibrational analysis of the H3+-molecule

Investigation of Nonadiabatic Effects for the Vibrational Spectrum of a Triatomic Molecule: The Use of a Single Potential Energy Surface with Distance-Dependent Masses for H3. J Phys Chem A 2017;121:7016-7030. [PMID: 28820589 DOI: 10.1021/acs.jpca.7b04703]

On the basis of first-principles, the influence of nonadiabatic effects on the vibrational bound states of H₃⁺ has been investigated using distance-dependent reduced masses and only one single potential energy surface. For these new vibrational calculations, potentials based on explicitly correlated wave functions are used where, in addition, adiabatic corrections and relativistic contributions are taken into account. For the first time, several different fully distance-dependent reduced mass surfaces in three dimensions have been incorporated in the vibrational calculations.

A systematic investigation of the ground state potential energy surface of H3+

Based on different ab initio electronic structure calculations (CI-R12 and Gaussian Geminals) of the Born-Oppenheimer electronic energy E(BO) of H(3)(+) from high to highest quality, we build up a potential energy surface which represents a highly reliable form of the topology of the whole potential region, locally and globally. We use the CI-R12 method in order to get within reasonable CPU-time a relatively dense grid of energy points. We demonstrate that CI-R12 is good enough to give an accurate surface, i.e., Gaussian Geminals are not absolutely necessary. For different types of potential energy surface fits, we performed variational calculations of all bound vibrational states, including resonances above the dissociation limit, for total angular momentum J = 0. We clarify the differences between different fits of the energy to various functional forms of the potential surface. Small rms-values (<1 cm(-1)) of the fit do not provide precise information about the interpolatory behaviour of the fit functions.

Rovibrational energy levels of H3(+) with energies above the barrier to linearity

The H(3)(+) potential energy surface is sampled at 5900 geometries with the emphasis on the nonequilibrium and asymptotic points. Apart from the Born-Oppenheimer energy converged to the accuracy better than 0.02 cm(-1), the adiabatic and the leading relativistic corrections are computed at each geometry. To represent analytically the potential energy surface, the parameters of a power series are determined by fitting to the computed energy points. Possible choice of nuclear masses simulating the nonadiabatic effects in solving the nuclear Schrodinger equation is analyzed. A set of theoretically predicted rovibrational transitions is confronted with the experimental data in the 10,700-13,700 cm(-1) window of the spectra.

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.

Explicitly correlated potential energy surface of H3+, including relativistic and adiabatic corrections

After a short historical account of the theory of the H3+ ion, two ab initio methods are reviewed that allow the computation of the ground-state potential energy surface (PES) of H3+ in the Born-Oppenheimer (BO) approximation, with microhartree or even sub-microhartree accuracy, namely the R12 method and the method of explicitly correlated Gaussians. The BO-PES is improved by the inclusion of relativistic effects and adiabatic corrections. It is discussed how non-adiabatic effects on rotation and vibration can be simulated by corrections to the moving nuclear masses. The importance of the appropriate analytic fit to the computed points of the PES for the subsequent computation of the rovibronic spectrum is addressed. Some recent extensions of the computed PES in the energy region above the barrier to linearity are reviewed. This involves a large set of input geometries and the correct treatment of the dissociation asymptotics, including the coupling with the first excited singlet state. Some comments on this state as well as on the lowest triplet state of H3+ are made. The paper ends with a few remarks on the ion H5+.

Near-infrared spectroscopy of H3+ above the barrier to linearity

Since the Royal Society Discussion Meeting on H3+ in 2000, the laboratory spectroscopy of H3+ has entered a new regime. For the first time, transitions of H3+ above the barrier to linearity have been observed. A highly sensitive near-infrared spectrometer based on a titanium:sapphire laser and incorporating a dual-beam, double-modulation technique with bidirectional optical multi-passing has been developed in order to detect these transitions, which are more than 4600 times weaker than the fundamental band. We discuss our recent work on the 2v1 + 2v2(2) <-- 0, 3v1 + v2(1) <-- 0, v1 + 4v2(2) <-- 0, v1 + 4V2(4) <-- 0 and 2v1 + 3v2(1) <-- 0 combination bands and the 5v2(1) <-- 0, 5v2(3) <-- 0, 52(5) <-- 0 and 6v2(2) <-- 0 overtone bands. Experimentally determined energy levels provide a critical test of ab initio calculations in this challenging energy regime (greater than 10,000 cm(-1)). By comparing the experimental energy levels and theoretical energy levels from ab initio calculations in which the adiabatic and relativistic corrections are incorporated, the extent of higher-order effects such as non-adiabatic and radiative corrections is revealed.

Time-Dependent Reactive Scattering for the System H- + D2 ↔ HD + D- and Comparison with H- + H2 ↔ H2 + H-

This work presents results of quantum mechanical calculations of reaction probabilities for the ion-neutral molecule collisions H- + D2 <--> HD + D-. Time-dependent wave packet propagations for total angular momentum J not equal to 0, including the full Coriolis coupling, are performed. The calculated state-to-state reaction probabilities using product Jacobi coordinates are compared with energy-resolved reaction probabilities calculated with the flux-operator using reactant Jacobi coordinates and with time-independent calculations. Differences between nearly converged integral cross sections and those using the J-shifting method and centrifugal sudden approximation and comparison with experimental results will be presented.

Accurate double many-body expansion potential energy surface for triplet H3+. II. The upper adiabatic sheet (2 3A′)

We report on a global potential energy hypersurface for the upper sheet of the lowest triplet state of H3+. The analytic representation is based on the double many-body expansion theory. The ab initio data points, calculated with a large cc-pV5Z basis, are represented with a root mean square deviation of only 5.54 cm(-1) in the energy region below the H(+)+2H(2S) dissociation threshold. The quasi-bound vibronic states supported by this surface have also been calculated.