Theoretical Determination of the Dissociation Energy of Molecular Hydrogen

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.

Bound-State Relativistic Quantum Electrodynamics: A Perspective for Precision Physics with Atoms and Molecules

Precision physics aims to use atoms and molecules to test and develop the fundamental theory of matter, possibly beyond the Standard Model. Most of the atomic and molecular phenomena are described by the quantum electrodynamics (QED) sector of the Standard Model. Do we have the computational tools, algorithms, and practical equations for the most possible complete computation of atoms and molecules within the QED sector? What is the fundamental equation to start with? Is it still Schrödinger's wave equation for molecular matter, or is there anything beyond that? This paper provides a concise overview of the relativistic QED framework and recent numerical developments targeting precision physics and spectroscopy applications with common features of the robust and successful relativistic quantum chemistry methodology.

Precision spectroscopy of molecular hydrogen

Precision measurements on the hydrogen molecule are of fundamental importance in understanding molecular theory. Comparison of accurate experimental data and theoretical results are used to test the quantum electrodynamics theory and determine physical constants used in the calculation. We review recent advances and perspectives in the precision spectroscopy of molecular hydrogen, representing state-of-the-art molecular spectroscopy methods and cutting-edge high-precision calculations.

Evaluation of the Bethe Logarithm: From Atom to Chemical Reaction

A general computational scheme for the (nonrelativistic) Bethe logarithm is developed, opening the route to "routine" evaluation of the leading-order quantum electrodynamics correction (QED) relevant for spectroscopic applications for small polyatomic and polyelectronic molecular systems. The implementation relies on the Schwartz method and minimization of a Hylleraas functional. In relation with electronically excited states, a projection technique is considered, which ensures positive definiteness of the functional over the entire parameter (photon momentum) range. Using this implementation, the Bethe logarithm is converged to a relative precision better than 1:10³ for selected electronic states of the two-electron H₂ and H₃⁺, and the three-electron He₂⁺ and H+H₂ molecular systems. The present work focuses on nuclear configurations near the local minimum of the potential energy surface, but the computations can be repeated also for other structures.

Dipolar physics: a review of experiments with magnetic quantum gases

Since the achievement of quantum degeneracy in gases of chromium atoms in 2004, the experimental investigation of ultracold gases made of highly magnetic atoms has blossomed. The field has yielded the observation of many unprecedented phenomena, in particular those in which long-range and anisotropic dipole-dipole interactions (DDIs) play a crucial role. In this review, we aim to present the aspects of the magnetic quantum-gas platform that make it unique for exploring ultracold and quantum physics as well as to give a thorough overview of experimental achievements. Highly magnetic atoms distinguish themselves by the fact that their electronic ground-state configuration possesses a large electronic total angular momentum. This results in a large magnetic moment and a rich electronic transition spectrum. Such transitions are useful for cooling, trapping, and manipulating these atoms. The complex atomic structure and large dipolar moments of these atoms also lead to a dense spectrum of resonances in their two-body scattering behaviour. These resonances can be used to control the interatomic interactions and, in particular, the relative importance of contact over dipolar interactions. These features provide exquisite control knobs for exploring the few- and many-body physics of dipolar quantum gases. The study of dipolar effects in magnetic quantum gases has covered various few-body phenomena that are based on elastic and inelastic anisotropic scattering. Various many-body effects have also been demonstrated. These affect both the shape, stability, dynamics, and excitations of fully polarised repulsive Bose or Fermi gases. Beyond the mean-field instability, strong dipolar interactions competing with slightly weaker contact interactions between magnetic bosons yield new quantum-stabilised states, among which are self-bound droplets, droplet assemblies, and supersolids. Dipolar interactions also deeply affect the physics of atomic gases with an internal degree of freedom as these interactions intrinsically couple spin and atomic motion. Finally, long-range dipolar interactions can stabilise strongly correlated excited states of 1D gases and also impact the physics of lattice-confined systems, both at the spin-polarised level (Hubbard models with off-site interactions) and at the spinful level (XYZ models). In the present manuscript, we aim to provide an extensive overview of the various related experimental achievements up to the present.

Ab Initio Partition Functions and Thermodynamic Quantities for the Molecular Hydrogen Isotopologues

In this work, we calculate the partition functions and thermodynamic quantities of molecular hydrogen isotopologues using the rovibrational energy levels provided by the highly accurate ab initio adiabatic potential energy functions recently determined by Pachucki and Komasa (Pachucki, K.; Komasa, J. J. Chem. Phys. 2014, 141, 224103). The partition functions are calculated by including all bound energy levels of the isotopologues, up to their dissociation limits, plus the quasi-bound levels lying below the centrifugal potential barriers. For the homonuclear isotopologues, H2, D2, and T2, we also determine the partition functions and thermodynamic quantities of the normal mixtures using the statistical treatment recently proposed by Colonna et al. (Colonna, G.; D'Angola, A.; Capitelli, M. Int. J. Hydrogen Energy 2012, 37, 9656) based on the definition of the partition function of the mixture, which avoids inconsistencies in the values of the thermodynamic quantities depending directly on the internal partition function, in the high-temperature limit.

Nonadiabatic, Relativistic, and Leading-Order QED Corrections for Rovibrational Intervals of ^{4}He_{2}^{+} (X ^{2}Σ_{u}^{+})

The rovibrational intervals of the ^{4}He_{2}^{+} molecular ion in its X ^{2}Σ_{u}^{+} ground electronic state are computed by including the nonadiabatic, relativistic, and leading-order quantum-electrodynamics corrections. Good agreement of theory and experiment is observed for the rotational excitation series of the vibrational ground state and the fundamental vibration. The lowest-energy rotational interval is computed to be 70.937 69(10)  cm^{-1} in agreement with the most recently reported experimental value, 70.937 589(23)(60)_{sys}  cm^{-1} [L. Semeria et al., Phys. Rev. Lett. 124, 213001 (2020)PRLTAO0031-900710.1103/PhysRevLett.124.213001].

Analysis of QED and non-adiabaticity effects on the rovibrational spectrum of H3 + using geometry-dependent effective nuclear masses

The influence of QED effects (including one- and two-electron Lamb-shift, Araki-Sucher term, one-loop self-energy, and finite nuclear size correction) together with non-adiabatic effects on the rovibrational bound states of H₃ ⁺ has been investigated. Non-adiabaticity is modeled by using geometry-dependent effective nuclear masses together with only one single potential energy surface. In conclusion, for rovibrational states below 20 000 cm^-1, QED and relativistic effects do nearly compensate, and a potential energy surface based on Born-Oppenheimer energies and diagonal adiabatic corrections has nearly the same quality as the one including relativity with QED; the deviations between the two approaches for individual rovibrational states are mostly below 0.02 cm^-1. The inclusion of non-adiabatic effects is important, and it reduces deviations from experiments mostly below 0.1 cm^-1.

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.

Nonadiabatic QED Correction to the Dissociation Energy of the Hydrogen Molecule

The quantum electrodynamic correction to the energy of the hydrogen molecule has been evaluated without expansion in the electron-proton mass ratio. The obtained results significantly improve the accuracy of theoretical predictions reaching the level of 1 MHz for the dissociation energy, in very good agreement with the parallel measurement [Hölsch et al., Phys. Rev. Lett. 122, 103002 (2019)PRLTAO0031-900710.1103/PhysRevLett.122.103002]. Molecular hydrogen has thus become a cornerstone of ultraprecise quantum chemistry, which opens perspectives for determination of fundamental physical constants from its spectra.

Benchmarking Theory with an Improved Measurement of the Ionization and Dissociation Energies of H_{2}

The dissociation energy of H_{2} represents a benchmark quantity to test the accuracy of first-principles calculations. We present a new measurement of the energy interval between the EF ^{1}Σ_{g}^{+}(v=0,N=1) state and the 54p1_{1} Rydberg state of H_{2}. When combined with previously determined intervals, this new measurement leads to an improved value of the dissociation energy D_{0}^{N=1} of ortho-H_{2} that has, for the first time, reached a level of uncertainty that is 3 times smaller than the contribution of about 1 MHz resulting from the finite size of the proton. The new result of 35 999.582 834(11)  cm^{-1} is in remarkable agreement with the theoretical result of 35 999.582 820(26)  cm^{-1} obtained in calculations including high-order relativistic and quantum-electrodynamics corrections, as reported in the following Letter [M. Puchalski, J. Komasa, P. Czachorowski, and K. Pachucki, Phys. Rev. Lett. 122, 103003 (2019)PRLTAO0031-900710.1103/PhysRevLett.122.103003]. This agreement resolves a recent discrepancy between experiment and theory that had hindered a possible use of the dissociation energy of H_{2} in the context of the current controversy on the charge radius of the proton.

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.

Fundamental vibration frequency and rotational structure of the first excited vibrational level of the molecular helium ion ( He 2 + )

The term values of the rotational levels of the first excited vibrational state of the electronic ground state of He 2 + with a rotational quantum number N ⁺ ≤ 13 have been determined with an accuracy of 1.2 × 10^-3 cm^-1 (∼35 MHz) by multichannel-quantum-defect-theory-assisted Rydberg spectroscopy of metastable He₂. Comparison of the experimental term values with the most accurate ab initio results for He 2 + available in the literature [W.-C. Tung, M. Pavanello, and L. Adamowicz, J. Chem. Phys. 136, 104309 (2012)] reveals inconsistencies between the theoretical and experimental results that increase with increasing rotational quantum numbers. The fundamental vibrational wavenumber of He 2 + was determined to be 1628.3832(12) cm^-1 by fitting effective molecular constants to the obtained term values.

Nonadiabatic Relativistic Correction to the Dissociation Energy of H_{2}, D_{2}, and HD

The relativistic correction to the dissociation energy of H_{2}, D_{2}, and HD molecules has been accurately calculated without expansion in the small electron-nucleus mass ratio. The obtained results indicate the significance of nonadiabatic effects and resolve the discrepancy of theoretical predictions with recent experimental values for H_{2} and D_{2}. While the theoretical accuracy is now significantly improved and is higher than the experimental one, we observe about 3σ discrepancy for the dissociation energy of HD, which requires further investigation.

Dissociation Energy of the Hydrogen Molecule at 10^{-9} Accuracy

The ionization energy of ortho-H_{2} has been determined to be E_{I}^{o}(H_{2})/(hc)=124 357.238 062(25)  cm^{-1} from measurements of the GK(1,1)-X(0,1) interval by Doppler-free, two-photon spectroscopy using a narrow band 179-nm laser source and the ionization energy of the GK(1,1) state by continuous-wave, near-infrared laser spectroscopy. E_{I}^{o}(H_{2}) was used to derive the dissociation energy of H_{2}, D_{0}^{N=1}(H_{2}), at 35 999.582 894(25)  cm^{-1} with a precision that is more than one order of magnitude better than all previous results. The new result challenges calculations of this quantity and represents a benchmark value for future relativistic and QED calculations of molecular energies.

Generalized elimination of the global translation from explicitly correlated Gaussian functions

This paper presents the multi-channel generalization of the center-of-mass kinetic energy elimination approach [B. Simmen et al., Mol. Phys. 111, 2086 (2013)] when the Schrödinger equation is solved variationally with explicitly correlated Gaussian functions. The approach has immediate relevance in many-particle systems which are handled without the Born-Oppenheimer approximation and can be employed also for Dirac-type Hamiltonians. The practical realization and numerical properties of solving the Schrödinger equation in laboratory-frame Cartesian coordinates are demonstrated for the ground rovibronic state of the H₂⁺={p⁺,p⁺,e^-} ion and the H₂ = {p⁺, p⁺, e^-, e^-} molecule.

Relativistic and QED Effects in the Fundamental Vibration of T_{2}

The hydrogen molecule has become a test ground for quantum electrodynamical calculations in molecules. Expanding beyond studies on stable hydrogenic species to the heavier radioactive tritium-bearing molecules, we report on a measurement of the fundamental T_{2} vibrational splitting (v=0→1) for J=0-5 rotational levels. Precision frequency metrology is performed with high-resolution coherent anti-Stokes Raman spectroscopy at an experimental uncertainty of 10-12 MHz, where sub-Doppler saturation features are exploited for the strongest transition. The achieved accuracy corresponds to a 50-fold improvement over a previous measurement, and it allows for the extraction of relativistic and QED contributions to T_{2} transition energies.

Toward a Determination of the Proton-Electron Mass Ratio from the Lamb-Dip Measurement of HD

Precision spectroscopy of the hydrogen molecule is a test ground of quantum electrodynamics (QED), and it may serve for the determination of fundamental constants. Using a comb-locked cavity ring-down spectrometer, for the first time, we observed the Lamb-dip spectrum of the R(1) line in the overtone of hydrogen deuteride (HD). The line position was determined to be 217 105 182.79±0.03_{stat}±0.08_{syst}  MHz (δν/ν=4×10^{-10}), which is the most accurate rovibrational transition ever measured in the ground electronic state of molecular hydrogen. Moreover, from calculations including QED effects up to the order m_{e}α^{6}, we obtained predictions for this R(1) line as well as for the HD dissociation energy, which are less accurate but signaling the importance of the complete treatment of nonadiabatic effects. Provided that the theoretical calculation reaches the same accuracy, the present measurement will lead to a determination of the proton-to-electron mass ratio with a precision of 1.3 parts per billion.

Sub-Doppler Frequency Metrology in HD for Tests of Fundamental Physics

Weak transitions in the (2,0) overtone band of the hydrogen deuteride molecule at λ=1.38  μm were measured in saturated absorption using the technique of noise-immune cavity-enhanced optical heterodyne molecular spectroscopy. Narrow Doppler-free lines were interrogated with a spectroscopy laser locked to a frequency comb laser referenced to an atomic clock to yield transition frequencies [R(1)=217105181895(20)  kHz; R(2)=219042856621(28)  kHz; R(3)=220704304951(28)  kHz] at three orders of magnitude improved accuracy. These benchmark values provide a test of QED in the smallest neutral molecule, and they open up an avenue to resolve the proton radius puzzle, as well as constrain putative fifth forces and extra dimensions.

Quantum Forces from Dark Matter and Where to Find Them

We observe that sub-GeV dark matter (DM) induces Casimir-Polder forces between nucleons that can be accessed by experiments from nuclear to molecular scales. We calculate the nucleon-nucleon potentials arising in the DM effective theory and note that their main features are fixed by dimensional analysis and the optical theorem. Molecular spectroscopy and neutron scattering turn out be DM search experiments and are found to be complementary to nucleon-based DM direct detection. Existing data set limits on DM with mass up to ∼3-50  MeV and with effective interaction up to the O(10-100)  MeV scale, constraining a region typically difficult to reach for other experiments.

Deep-Ultraviolet Frequency Metrology of H_{2} for Tests of Molecular Quantum Theory

Molecular hydrogen and its isotopic and ionic species are benchmark systems for testing quantum chemical theory. Advances in molecular energy structure calculations enable the experimental verification of quantum electrodynamics and potentially a determination of the proton charge radius from H_{2} spectroscopy. We measure the ground state energy in ortho-H_{2} relative to the first electronically excited state by Ramsey-comb laser spectroscopy on the EF^{1}Σ_{g}^{+}-X^{1}Σ_{g}^{+}(0,0) Q1 transition. The resulting transition frequency of 2 971 234 992 965(73) kHz is 2 orders of magnitude more accurate than previous measurements. This paves the way for a considerably improved determination of the dissociation energy (D_{0}) for fundamental tests with molecular hydrogen.

Determination of the Spin-Rotation Fine Structure of He_{2}^{+}

Measuring spin-rotation intervals in molecular cations is challenging, particularly so when the ions do not have electric-dipole-allowed rovibrational transitions. We present a method, based on an angular-momentum basis transformation, to determine the spin-rotational fine structure of molecular ions from the fine structure of high Rydberg states. The method is illustrated by the determination of the so far unknown spin-rotation fine structure of the fundamentally important He_{2}^{+} ion in the X ^{2}Σ_{u}^{+} state. The fine-structure splittings of the v^{+}=0, N^{+}=1, 3, and 5 levels of He_{2}^{+} are 7.96(14), 17.91(32), and 28.0(6) MHz, respectively. The experiment relies on the use of single-mode cw radiation to record spectra of high Rydberg states of He_{2} from the a ^{3}Σ_{u}^{+} metastable state.

Complete α^{6} m Corrections to the Ground State of H_{2}

We perform the calculation of all relativistic and quantum electrodynamic corrections of the order of α^{6} m to the ground electronic state of a hydrogen molecule and present improved results for the dissociation and the fundamental transition energies. These results open the window for the high-precision spectroscopy of H_{2} and related low-energy tests of fundamental interactions.

Precision measurement of the rotational energy-level structure of the three-electron molecule He2. J Chem Phys 2016;145:204301. [PMID: 27908106 DOI: 10.1063/1.4967256]

The term values of all rotational levels of the ⁴He₂⁺ X^{+  2}Σ_u⁺ (ν⁺=0) ground vibronic state with rotational quantum number N⁺ ≤ 19 have been determined with an accuracy of 8 × 10^-4 cm^-1 (∼25 MHz) by multichannel-quantum-defect-theory-assisted Rydberg spectroscopy of metastable He₂^∗. Comparison of these term values with term values recently calculated ab initio by Tung et al. [J. Chem. Phys. 136, 104309 (2012)] reveals discrepancies that rapidly increase with increasing rotational quantum number and reach values of 0.07 cm^-1 (∼2.1 GHz) at N⁺ = 19.

An Evaluation of Gas Phase Enthalpies of Formation for Hydrogen-Oxygen (HxOy) Species

We have compiled gas phase enthalpies of formation for nine hydrogen-oxygen species (HxOy) and selected recommended values for H, O, OH, H2O, HO2, H2O2, O3, HO3, and H2O3. The compilation consists of values derived from experimental measurements, quantum chemical calculations, and prior evaluations. This work updates the recommended values in the NIST-JANAF (1985) and Gurvich et al. (1989) thermochemical tables for seven species. For two species, HO3 and H2O3 (important in atmospheric chemistry) and not found in prior thermochemical evaluations, we also provide supplementary data consisting of molecular geometries, vibrational frequencies, and torsional potentials which can be used to compute thermochemical functions. For all species, we also provide supplementary data consisting of zero point energies, vibrational frequencies, and ion reaction energetics.

Observation and Calculation of the Quasibound Rovibrational Levels of the Electronic Ground State of H2+

Although the existence of quasibound rotational levels of the X^{+} ^{2}Σ_{g}^{+}  ground state of H_{2}^{+} was predicted a long time ago, these states have never been observed. Calculated positions and widths of quasibound rotational levels located close to the top of the centrifugal barriers have not been reported either. Given the role that such states play in the recombination of H(1s) and H^{+} to form H_{2}^{+}, this lack of data may be regarded as one of the largest unknown aspects of this otherwise accurately known fundamental molecular cation. We present measurements of the positions and widths of the lowest-lying quasibound rotational levels of H_{2}^{+} and compare the experimental results with the positions and widths we calculate using a potential model for the X^{+} state of H_{2}^{+} which includes adiabatic, nonadiabatic, relativistic, and radiative corrections to the Born-Oppenheimer approximation.

Precision Spectroscopy in Cold Molecules: The Lowest Rotational Interval of He_{2}^{+} and Metastable He_{2}

Multistage Zeeman deceleration was used to generate a slow, dense beam of translationally cold He_{2} molecules in the metastable a ^{3}Σ_{u}^{+} state. Precision measurements of the Rydberg spectrum of these molecules at high values of the principal quantum number n have been carried out. The spin-rotational state selectivity of the Zeeman-deceleration process was exploited to reduce the spectral congestion, minimize residual Doppler shifts, resolve the Rydberg series around n=200 and assign their fine structure. The ionization energy of metastable He_{2} and the lowest rotational interval of the X^{+} ^{2}Σ_{u}^{+} (ν^{+}=0) ground state of ^{4}He_{2}^{+} have been determined with unprecedented precision and accuracy by Rydberg-series extrapolation. Comparison with ab initio predictions of the rotational energy level structure of ^{4}He_{2}^{+} [W.-C. Tung, M. Pavanello, and L. Adamowicz, J. Chem. Phys. 136, 104309 (2012)] enabled us to quantify the magnitude of relativistic and quantum-electrodynamics contributions to the fundamental rotational interval of He_{2}^{+}.

Communication: Test of quantum chemistry in vibrationally hot hydrogen molecules

Precision measurements are performed on highly excited vibrational quantum states of molecular hydrogen. The v = 12, J = 0 - 3 rovibrational levels of H2 (X(1)Σg (+)), lying only 2000 cm(-1) below the first dissociation limit, were populated by photodissociation of H2S and their level energies were accurately determined by two-photon Doppler-free spectroscopy. A comparison between the experimental results on v = 12 level energies with the best ab initio calculations shows a good agreement, where the present experimental accuracy of 3.5 × 10(-3) cm(-1) is more precise than theory, hence providing a gateway to further test theoretical advances in this benchmark quantum system.

Modelling non-adiabatic effects in H₃⁺: solution of the rovibrational Schrödinger equation with motion-dependent masses and mass surfaces

Introducing different rotational and vibrational masses in the nuclear-motion Hamiltonian is a simple phenomenological way to model rovibrational non-adiabaticity. It is shown on the example of the molecular ion H3(+), for which a global adiabatic potential energy surface accurate to better than 0.1 cm(-1) exists [M. Pavanello, L. Adamowicz, A. Alijah, N. F. Zobov, I. I. Mizus, O. L. Polyansky, J. Tennyson, T. Szidarovszky, A. G. Császár, M. Berg et al., Phys. Rev. Lett. 108, 023002 (2012)], that the motion-dependent mass concept yields much more accurate rovibrational energy levels but, unusually, the results are dependent upon the choice of the embedding of the molecule-fixed frame. Correct degeneracies and an improved agreement with experimental data are obtained if an Eckart embedding corresponding to a reference structure of D(3h) point-group symmetry is employed. The vibrational mass of the proton in H3(+) is optimized by minimizing the root-mean-square (rms) deviation between the computed and recent high-accuracy experimental transitions. The best vibrational mass obtained is larger than the nuclear mass of the proton by approximately one third of an electron mass, m(opt,p)((v))=m(nuc,p)+0.31224m(e). This optimized vibrational mass, along with a nuclear rotational mass, reduces the rms deviation of the experimental and computed rovibrational transitions by an order of magnitude. Finally, it is shown that an extension of the algorithm allowing the use of motion-dependent masses can deal with coordinate-dependent mass surfaces in the rovibrational Hamiltonian, as well.

Non-Born-Oppenheimer calculations of the pure vibrational spectrum of T₂ including relativistic corrections

We report very accurate calculations of the complete pure vibrational spectrum of the T2 molecule with an approach where the Born-Oppenheimer (BO) approximation is not assumed. As the considered states correspond to the zero total angular momentum, their non-BO wave functions are spherically symmetric and are expanded in terms of all-particle, one-center, spherically symmetric explicitly correlated Gaussian functions multiplied by even nonnegative powers of the internuclear distance. The nonrelativistic energies of the states obtained in the non-BO calculations are corrected for the relativistic effects of the order of α(2) (where α is the fine structure constant) calculated as expectation values of the operators representing these effects.