Accurate ab initio determination of the adiabatic potential energy function and the Born–Oppenheimer breakdown corrections for the electronic ground state of LiH isotopologues

Deformed explicitly correlated Gaussians

Deformed explicitly correlated Gaussian (DECG) basis functions are introduced, and their matrix elements are calculated. All matrix elements can be calculated analytically in a closed form, except the Coulomb one, which has to be approximated by a Gaussian expansion. The DECG basis functions can be used to solve problems with nonspherical potentials. One example of such potential is the dipole self-interaction term in the Pauli-Fierz Hamiltonian. Examples are presented showing the accuracy and necessity of deformed Gaussian basis functions to accurately solve light-matter coupled systems in cavity QED.

Diagonal Born-Oppenheimer corrections to the ground electronic state potential energy surfaces of ozone: improvement of ab initio vibrational band centers for the 16O3, 17O3 and 18O3 isotopologues

Mass-dependent diagonal Born-Oppenheimer corrections (DBOCs) to the ab initio electronic ground state potential energy surface for the main ¹⁶O₃ isotopologue and for homogeneous isotopic substitutions ¹⁷O₃ and ¹⁸O₃ of the ozone molecule are reported for the first time. The system being of strongly multiconfigurational character, multireference configuration interaction wave function ansatz with different complete active spaces was used. The reliable DBOC calculations with the targeted accuracy were possible to carry out up to about half of the dissociation threshold D₀. The comparison with the experimental band centers shows a significant improvement of the accuracy with respect to the best Born-Oppenheimer (BO) ab initio calculations reducing the total root-mean-squares (calculated-observed) deviations by about a factor of two. For the set of ¹⁶O₃ vibrations up to five bending and four stretching quanta, the mean (calculated-observed) deviations drop down from 0.7 cm^-1 (BO) to about 0.1 cm^-1, with the most pronounced improvement seen for bending states and for mixed bending-stretching polyads. In the case of bending band centers directly observed under high spectral resolutions, the errors are reduced by more than an order of magnitude down to 0.02 cm^-1 from the observed levels, approaching nearly experimental accuracy. A similar improvement for heavy isotopologues shows that the reported DBOC corrections almost remove the systematic BO errors in vibrational levels below D₀/2, though the scatter increases towards higher energies. The possible reasons for this finding, as well as remaining issues are discussed in detail. The reported results provide an encouraging accuracy validation for the multireference methods of the ab initio theory. New sets of ab initio vibrational states can be used for improving effective spectroscopic models for analyses of the observed high-resolution spectra, particularly in the cases of accidental resonances with "dark" states requiring accurate theoretical predictions.

Non-adiabatic effects in the H3+ spectrum

The effect of non-adiabatic coupling on the computed rovibrational energy levels amounts to about 2 cm-1 for H3+ and must be included in high-accuracy calculations. Different strategies to obtain the corresponding energy shifts are reviewed in the article. A promising way is to introduce effective vibrational reduced masses that depend on the nuclear configuration. A new empirical method that uses the stockholder atoms-in-molecules approach to this effect is presented and applied to H3+. Furthermore, a highly accurate potential energy surface for the D3+ isotopologue, which includes relativistic and leading quantum electrodynamic terms, is constructed and used to analyse the observed rovibrational frequencies for this molecule. Accurate band origins are obtained that improve existing data. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H3+, H5+ and beyond'.

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.

Electronic structure and time-dependent description of rotational predissociation of LiH

The adiabatic potential energy curves of the ¹Σ⁺ and ¹Π states of the LiH molecule were calculated. They correlate asymptotically to atomic states, such as 2s + 1s, 2p + 1s, 3s + 1s, 3p + 1s, 3d + 1s, 4s + 1s, 4p + 1s and 4d + 1s. A very good agreement was found between our calculated spectroscopic parameters and the experimental ones. The dynamics of the rotational predissociation process of the 1¹Π state were studied by solving the time-dependent Schrödinger equation. The classical experiment of Velasco [Can. J. Phys., 1957, 35, 1204] on dissociation in the 1¹Π state is explained for the first time in detail.

Matrix elements of explicitly correlated Gaussian basis functions with arbitrary angular momentum

A new algorithm for calculating the Hamiltonian matrix elements with all-electron explicitly correlated Gaussian functions for quantum-mechanical calculations of atoms with arbitrary angular momentum is presented. The calculations are checked on several excited states of three and four electron systems. The presented formalism can be used as unified framework for high accuracy calculations of properties of small atoms and molecules.

Algorithms for calculating mass-velocity and Darwin relativistic corrections with n-electron explicitly correlated Gaussians with shifted centers

Algorithms for calculating the leading mass-velocity (MV) and Darwin (D) relativistic corrections are derived for electronic wave functions expanded in terms of n-electron explicitly correlated Gaussian functions with shifted centers and without pre-exponential angular factors. The algorithms are implemented and tested in calculations of MV and D corrections for several points on the ground-state potential energy curves of the H2 and LiH molecules. The algorithms are general and can be applied in calculations of systems with an arbitrary number of electrons.

Polarization of core orbitals and computation of nuclear quadrupole coupling constants using Gaussian basis sets

Most standard Gaussian basis sets for first row atoms, even large sets designed to converge on a 'complete basis set' limit, systematically overestimate the electric field gradient at nuclear sites for first row atoms, resulting in errors of up to 15% in the computation of nuclear quadrupole coupling constants. This error results from a failure to include tight d functions, which permit the core 1s orbitals to distort under the influence of the field of the nuclear quadrupole. Augmentation of standard basis sets with a single set of single-exponent d functions, matched to the reciprocal square of the nominal 1s radius, reduces these errors by up to 90%.

Non-adiabatic effects in thermochemistry, spectroscopy and kinetics: the general importance of all three Born–Oppenheimer breakdown corrections

Analytical and numerical solutions describing Born–Oppenheimer breakdown in a simple, widely applicable, model depict shortcomings in modern computational methods.

Studies of singlet Rydberg series of LiH derived from Li(nl) + H(1s), with n ≤ 6 and l ≤ 4

The 50 singlet states of LiH composed of 49 Rydberg states and one non-Rydberg ionic state derivable from Li(nl) + H(1s), with n ≤ 6 and l ≤ 4, are studied using the multi-reference configuration interaction method combined with the Stuttgart/Köln group's effective core potential/core polarization potential method. Basis functions that can yield energy levels up to the 6g orbital of Li have been developed, and they are used with a huge number of universal Kaufmann basis functions for Rydberg states. The systematics and regularities of the physical properties such as potential energies, quantum defects, permanent dipole moments, transition dipole moments, and nonadiabatic coupling matrix elements of the Rydberg series are studied. The behaviors of potential energy curves and quantum defect curves are explained using the Fermi approximation. The permanent dipole moments of the Rydberg series reveal that they are determined by the sizes of the Rydberg orbitals, which are proportional to n(2). Interesting mirror relationships of the dipole moments are observed between l-mixed Rydberg series, with the rule Δl = ±1, except for s-d mixing, which is also accompanied by n-mixing. The members of the l-mixed Rydberg series have dipole moments with opposite directions. The first derivatives of the dipole moment curves, which show the charge-transfer component, clearly show not only mirror relationships in terms of direction but also oscillations. The transition dipole moment matrix elements of the Rydberg series are determined by the small-r region, with two consequences. One is that the transition dipole moment matrix elements show n(-3/2) dependence. The other is that the magnitudes of the transition dipole moment matrix elements decrease rapidly as l increases.