Rovibrational levels of helium hydride ion

A nearly complete treatment of the effect of non-adiabaticity on rovibrational energies of H3+ (Part III)

In this article, significant contributions of non-adiabaticity for the rovibrational bound states up to 25 000 cm-1 and total angular momentum J = 0-20 of H3+ are investigated. A coupled-perturbed full configuration interaction (CP-FCI) treatment is applied to calculate all couplings between electronic states caused by the nuclear motion. These derivative couplings were evaluated up to the second order by means of a perturbation treatment and include all nuclear Cartesian first and second derivatives of the electronic wavefunctions. In particular, the coupling of special derivatives with respect to r and R in the Jacobi coordinate representation is more significant than thought. The perturbation approach is especially optimal for the treatment of weak non-adiabaticity in case of rovibrational energies in H3+ and had not been available before for H3+ or other triatomics. Using exclusively Gaussian basis functions for CP-FCI appears to be sufficient, because explicit correlated wavefunctions are already used for all other potential energy contributions. Our work is an extension of earlier non-adiabatic investigations based on first derivative couplings of electronic states that led to the concept of geometry-dependent effective nuclear masses and which needs only a single potential energy surface for the dynamics. The implementation allows us to include all non-adiabatic effects up to the order of O(μ-2), μ being the reduced nuclear mass. Our treatment works for any isotopologue and for the whole potential energy curve or surface. By this treatment, a further reduction in deviations to experimental data for most rovibrational levels to less than 0.1 cm-1 is possible. For the related transition frequencies, 1366 of 1720 known rovibrational transitions in H3+ have deviations less than 0.1 cm-1 without using any empirically adjustable parameters or optimizing the nuclear mass for a specific transition. For many questionable assignments (deviations >0.3 cm-1) of observed transitions in H3+, a new labeling is proposed.

Effective non-adiabatic Hamiltonians for the quantum nuclear motion over coupled electronic states

The quantum mechanical motion of the atomic nuclei is considered over a single- or a multidimensional subspace of electronic states which is separated by a gap from the rest of the electronic spectrum over the relevant range of nuclear configurations. The electron-nucleus Hamiltonian is block-diagonalized up to O(εⁿ⁺¹) through a unitary transformation of the electronic subspace, and the corresponding nth-order effective Hamiltonian is derived for the quantum nuclear motion. Explicit but general formulas are given for the second- and the third-order corrections. As a special case, the second-order Hamiltonian corresponding to an isolated electronic state is recovered which contains the coordinate-dependent mass-correction terms in the nuclear kinetic energy operator. For a multidimensional, explicitly coupled electronic band, the second-order Hamiltonian contains the usual Born-Oppenheimer terms and nonadiabatic corrections, but generalized mass-correction terms appear as well. These, earlier neglected terms, perturbatively account for the outlying (discrete and continuous) electronic states not included in the explicitly coupled electronic subspace.

Non-adiabatic mass-correction functions and rovibrational states of 4 He 2 + ( X 2 Σ u + )

The mass-correction functions in the second-order non-adiabatic Hamiltonian are computed for the ⁴ He 2 + molecular ion using the variational method, floating explicitly correlated Gaussian functions, and a general coordinate-transformation formalism. When non-adiabatic rovibrational energy levels are computed using these (coordinate-dependent) mass-correction functions and a highly accurate potential energy and diagonal Born-Oppenheimer correction curve, significantly improved theoretical results are obtained for the nine rotational and two rovibrational intervals known from high-resolution spectroscopy experiments.

Non-adiabatic mass correction to the rovibrational states of molecules: Numerical application for the H 2 + molecular ion

General transformation expressions of the second-order non-adiabatic Hamiltonian of the atomic nuclei, including the kinetic-energy correction terms, are derived upon the change from laboratory-fixed Cartesian coordinates to general curvilinear coordinate systems commonly used in rovibrational computations. The kinetic-energy or so-called "mass-correction" tensor elements are computed with the stochastic variational method and floating explicitly correlated Gaussian functions for the H 2 + molecular ion in its ground electronic state. {Further numerical applications for the ⁴ He 2 + molecular ion are presented in the forthcoming paper, Paper II [E. Mátyus, J. Chem. Phys. 149, 194112 (2018)]}. The general, curvilinear non-adiabatic kinetic energy operator expressions are used in the examples, and non-adiabatic rovibrational energies and corrections are determined by solving the rovibrational Schrödinger equation including the diagonal Born-Oppenheimer as well as the mass-tensor corrections.