• Free to Read
  • Editors Choice
Quantum Electronic Structure

From First-Principles to Quantum Electrodynamics: Pushing the Limits of Theory with the Hydrogen Molecule
Click to copy article linkArticle link copied!

Open PDFSupporting Information (1)

Journal of Chemical Theory and Computation

Cite this: J. Chem. Theory Comput. 2025, 21, 24, 12664–12673
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jctc.5c01702
Published December 5, 2025

Copyright © 2025 American Chemical Society. This publication is licensed under these Terms of Use.

Abstract

Click to copy section linkSection link copied!

Modern spectroscopic techniques enable the determination of the spacing between rovibrational levels of H2 with a relative accuracy of approximately 10–11. At this extreme level of precision, subtle quantum electrodynamic (QED) effects, such as electron self-interaction and vacuum polarization, are probed. A theoretical model aiming to achieve similar accuracy must precisely describe not only these relatively small QED effects but also the more significant contributions related to electron correlation, coupling between electronic and nuclear motions, and relativistic effects. Although the hydrogen molecule exhibits most of the phenomena found in larger molecules, it is simple enough to meet the requirements mentioned above. In this article, we report on enhancements to the current capabilities of quantum mechanical calculations for the hydrogen molecule. We present a method based on exponential functions that fully captures electron correlation or, more broadly, interparticle correlation, enabling a comprehensive description of effects related to nuclear motion. Specifically, we solve the four-particle Schrödinger equation without invoking commonly used approximations such as the one-electron or the Born–Oppenheimer approximation. The only source of nonrelativistic energy error comes from the finite size of the basis set. The explicitly correlated nonadiabatic wave function used here is then employed to determine the relativistic and QED effects. As a result, the dissociation energy for the lowest rovibrational levels in the electronic ground state of H2 has been obtained with a relative accuracy of 7 × 10–10, while the frequencies of intervals between these levels have been determined with sub-MHz accuracy, corresponding to a relative accuracy of 3 × 10–9. In consequence, the discrepancies between the highest precision measurements and earlier theoretical predictions have been resolved.

This publication is licensed for personal use by The American Chemical Society.

Copyright © 2025 American Chemical Society

Special Issue

Published as part of Journal of Chemical Theory and Computation special issue “100 Years of Quantum Mechanics-All About Molecules”.

1. Introduction

Click to copy section linkSection link copied!

The hydrogen molecule has served as a fundamental benchmark for theoretical calculations since the advent of quantum mechanics more than a hundred years ago. In 1927, Heitler and London (1) provided the first qualitative explanation of molecular binding based on quantum mechanical principles. A few years later, James and Coolidge (2) reached an agreement with the then-available experimental data using an explicitly correlated wave function. Early calculations, including those by Kołos and Roothaan, (3) were based on the clamped nucleus approximation. Efforts to incorporate nuclear motion began with the evaluation of adiabatic and nonadiabatic corrections, initiated by Kołos and Wolniewicz in the early 1960s. (4−6) The most accurate calculations of these corrections to date have been performed within the framework of nonadiabatic perturbation theory (NAPT), (7) employing explicitly correlated exponential wave functions and expanding the perturbational series up to the second order in the electron-to-nucleus mass ratio. (8) The first direct nonadiabatic calculation of molecular energy, bypassing the intermediate step of evaluating internuclear potentials, is also attributed to Kołos and Wolniewicz, who used a product of their electronic wave function and a harmonic oscillator nuclear function. (9) A breakthrough in accurate direct nonadiabatic calculations of the nonrelativistic component of the ground-state energy was achieved by Kinghorn and Adamowicz, (10) who employed explicitly correlated Gaussian (ECG) functions, reaching the relative accuracy of 7 × 10–9. The accuracy was later improved to 3 × 10–11 by Puchalski et al., (11) also using ECGs. Currently, the most precise direct nonadiabatic results are obtained with exponential-type wave functions, yielding the nonrelativistic dissociation energy with an accuracy of 10–8 cm–1, corresponding to the relative precision of 10–13. (12)
In parallel with studies on the coupling of nuclear and electronic motion, attempts were made to estimate the contribution of relativistic effects to the energy of the hydrogen molecule. In 1961, Kołos and Wolniewicz published a pioneering work (4) that provided the first estimation of these relativistic effects within the framework of the Born–Oppenheimer (BO) approximation. Relativistic calculations were not pursued until 1993, when Wolniewicz presented ground-state relativistic energies augmented by approximate QED corrections. (13) A significant advancement occurred in 2009 with the introduction of explicitly correlated Gaussian (ECG) functions for evaluating relativistic matrix elements, which enabled accurate determination of the dissociation energy, including the complete leading-order quantum electrodynamic (QED) correction. (14,15) The development of the so-called rECG functions, incorporating the cusp 1 + r12/2 prefactor, facilitated efficient quadrature of nonstandard two-center integrals (16) and enabled evaluation of the higher-order QED correction potential. Subsequent studies addressed finite nuclear mass relativistic corrections using NAPT, leading to the derivation of the relativistic recoil potential, (17) followed by improvements in the numerical accuracy of the QED correction components. (18)
A fully nonadiabatic treatment of the recoil effect in the relativistic correction for H2, incorporating finite nuclear mass terms directly into the relativistic operators, was first rigorously implemented in 2018. (11,19) The complete expression for the leading nonadiabatic QED correction was subsequently introduced in 2019. (11) Nonetheless, the numerical evaluation of both types of corrections remained confined to the ground rovibrational state.
A more detailed account of the centennial history behind the effort to determine the precise total energy of the hydrogen molecule can be found in the review. (20) A key aspect of this theoretical research has been the ongoing interaction with the experimental results. Over the years, the fruitful interplay between theory and experiment has led to an increased precision in both fields. Of particular interest are measurable separations between energy levels including the dissociation energy. Thanks to the long lifetimes of rovibrational states, counted in days, (21) H2 is an attractive subject for highly precise spectroscopic measurements. The highly accurate rovibrational ladder of energy levels and the infrared (IR) spectrum of the hydrogen molecule are widely utilized across various disciplines, including astronomy, (22−26) plasma physics, (27) physical chemistry, (28,29) and metrology. (30,31)
Throughout a century of research on the hydrogen molecule, the agreement between theory and experiment has undergone continuous, albeit diminishing, fluctuations. For instance, just a decade ago, the theoretical predictions and measurements of the energy of rovibrational transitions in H2 isotopologues agreed to within 10 MHz (≈0.00033 cm–1). (32−52) This precision opened the possibility for testing QED effects on rotational transitions, (53) setting bounds on a hypothetical fifth force, (54) and on the number of extra spatial dimensions. (55,56) Recent advancements in spectroscopic techniques have enabled measurements to achieve an accuracy of the order of 10 kHz. (57−73) As a result, discrepancies between theory and experiment at the level of several megahertz have been revealed. It was suspected that the incomplete treatment of relativistic and QED recoil effects could be the main reason for these discrepancies.
In this paper, we report the results of our calculations, which were performed using a nonadiabatic exponential wave function. Such functions possess correct asymptotic properties but are extremely rarely used in molecular calculations due to difficulties in evaluating corresponding integrals. We have overcome these difficulties (74,75) and developed a method applicable to arbitrary rovibrational excitation levels of diatomic two-electron molecules, like H2, HeH+, and their isotopologues. (12,76−78) We also present results from using the direct nonadiabatic method to estimate the relativistic and QED corrections. With this approach, we were able to bring the theoretical predictions back into agreement with the measurements at a new sub-MHz level of accuracy.

2. Theoretical Background

Click to copy section linkSection link copied!

Our calculations are based on the nonrelativistic quantum electrodynamic (NRQED) theory, which describes the energy of a bound rovibrational level of a light molecule as a function of fine-structure constant α. This energy can be expanded in powers of α and its logarithm (atomic units are used for E(n) coefficients)
E(α)=mα2E(2)+mα4E(4)+mα4Efs(4)+mα5E(5)+mα6E(6)+mα7E(7)+...
(1)
Subsequent terms of this expansion can be identified as the nonrelativistic energy, relativistic correction, finite nuclear size correction, QED correction, and higher-order corrections. An important feature of NRQED is that the coefficients E(n) of this expansion can be expressed as expectation values of some effective Hamiltonians H(n) with the nonrelativistic wave function Ψ. These expectation values can be evaluated using a fully four-body wave function or, if this is infeasible, using approximations that involve the separation of nuclear and electronic variables, as in, e.g., NAPT. Obviously, the former approach, referred to as the direct nonadiabatic method, is preferable because it fully accounts for the effects of the finite nuclear mass.

2.1. The Wave Function

The wave function must reflect the rotational and electronic degrees of freedom present in the ground electronic state (Σg+1) of the hydrogen molecule. Ψ should also account for the coupling between the rotational angular momentum of nuclei and orbital angular momentum of electrons L. This coupling results in the total spatial angular momentum, represented by J. Considering these factors, we can construct the wave function Ψ as follows
ΨJ,M=ΛΨΛJ,M
(2)
where J and M are the total rotational and magnetic quantum numbers, respectively. The latter is the projection of J onto the laboratory frame Z axis. The summation index Λ is an eigenvalue of the projection of L on the internuclear axis and runs over subsequent electronic states Σg, Πg, Δg, ..., which correspond to |Λ| = 0, 1, 2, ..., respectively.
The functions ΨΛJ,M are represented as the linear combination
ΨΛJ,M={k}c{k}(1+P12)ψ{k},ΛJ,M
(3)
of four-particle basis functions
ψ{k},ΛJ,M=QΛeαΛRβΛ(ζ1+ζ2)Rk0r12k1η1k2η2k3ζ1k4ζ2k5
(4)
with elliptic-like coordinates ζ1 = r13 + r14, η1 = r13r14, ζ2 = r23 + r24, η2 = r23r24, and R=r34. From this point onward, indices 1 and 2 assigned to the distance between particles r, and later also to the particle momentum p and mass m, will refer to the electrons, while indices 3 and 4 will correspond to protons. The exponents ki are non-negative integers collectively denoted as {k}. The basis function is symmetric with respect to the inversion when k2 + k3 is even and antisymmetric when this sum is odd. The operator P12 permutes electrons, while the factor QΛ determines the electronic angular momentum associated with a given basis function. Explicit expressions for QΛ corresponding to the lowest angular momentum values are provided below
QΣ=YJMforJ0
(5)
QΠ=2J(J+1)ρiR(RiYJM)forJ1
(6)
QΔ=4(J1)J(J+1)(J+2)×12(ρiρj+ρjρiδijρ·ρ)×R2(RiRjYJM)forJ2
(7)
In addition to a J-dependent normalization constant, the QΛ consists of the electronic and the nuclear parts. The electronic part is constructed from combinations of the vectors ρρa (a = 1 or 2), which correspond to the components of the electron-nucleus vectors ra3 that are perpendicular to the molecular axis. These components are defined as ρai = (δijninj)ra3j, where n=R/R denotes the normalized internuclear vector. In this context, we utilize the Einstein summation convention. The nuclear part, in turn, involves derivatives with respect to nuclear coordinates of the solid harmonic YJM(n).
The basis function described in eq 4 is termed the nonadiabatic James–Coolidge (naJC) function. This name reflects its similarity to the classic two-electron function introduced initially by James and Coolidge. (2) Their inspiration came from the Hylleraas function, which was designed for the helium atom in 1929. (79)
The nonlinear (αΛ, βΛ) and linear c{k} parameters are determined variationally by solving the four-body Schrödinger equation H(2)ΨJ,M = E(2)ΨJ,M with the nonrelativistic Hamiltonian
H(2)=T+V
(8)
T=a=14m2mapa2
(9)
V=1r121r131r231r141r24+1r34
(10)
According to the notation introduced earlier, for H2, m1 = m2 = m is the electron mass, and m3 = m4 = mp denotes the mass of the proton. A more detailed description of the Ansatz and the results of nonrelativistic calculations performed using this wave function can be found in refs (12and74). At this point, we merely note that the nonrelativistic energy can be determined with an accuracy limited only by the uncertainties in physical constants such as the proton-to-electron mass ratio and the Rydberg constant.
The nonrelativistic wave function described above is utilized to calculate the expectation values of operators present in relativistic and QED Hamiltonians H(n). Relevant formulas are provided in the following subsections.

2.2. Relativistic Correction, E(4)

The relativistic correction is evaluated as an expectation value of the spin-independent Breit–Pauli Hamiltonian H(4). It consists of the mass velocity term, the Breit orbit–orbit terms, and the Dirac delta terms
E(4)=Ψ|a=14m38ma3pa4|Ψ12Ψ|a=13b=a+14m2sa,bmambpai(δijrab+rabirabjrab3)pbj|ΨπΨ|δ3(r12)|Ψ+π2(1+m2mp2)Ψ|a=12b=34δ3(rab)|Ψ
(11)
In the above, sa,b = 2[δa,1·δb,2 + δa,3·δb,4] – 1 provides a proper sign. To achieve the target of sub-MHz accuracy, it is essential to consider another small relativistic effect arising from the spatial distribution of the proton charge. The leading correction due to the finite nuclear size, which is formally of the order of α4, is given by (80)
Efs(4)=2π3Ψ|a,Xδ3(raX)|Ψrp2ƛ2
(12)
where rp is the root-mean-square charge radius of the proton and ƛ is the reduced electron Compton wavelength. All of the expectation values mentioned above were calculated using the naJC wave function. Detailed information regarding the newly developed integration methods, which are suitable for evaluating these expectation values, along with their regularization and numerical convergence, can be found in refs (75and81). Preliminary results for the lowest rotational levels of H2 are also shown therein.
In eq 11, the nucleus–nucleus Dirac delta interaction is omitted due to its negligible value. Contributions from the coupling between the nuclear spin and rotation as well as between nuclear spins are small because they depend quadratically on the electron-to-nucleus mass ratio. These effects are usually undetectable in infrared spectroscopy. Nonetheless, they were measured by Ramsey in 1971, (82) which allowed us, using the four-body Gaussian wave function, to determine the quadrupole moment of the deuteron to an accuracy better than for any other element of the periodic table. (30) In this work, the energy levels are averaged over the nuclear spin orientations; therefore, no hyperfine coupling is considered.

2.3. QED Correction, E(5)

The next term of the α-expansion (1) is the leading spin-independent quantum electrodynamic correction represented by the following expression, including terms up to the first order in the electron-to-nucleus mass ratio (11)
E(5)=23πm2μ2Dlnk076πΨ|1r123|Ψ76πa=12b=34mmbΨ|1rab3|Ψ+43a=12b=34{19302lnα+mmp(31612lnα)}Ψ|δ3(rab)|Ψ+(16415+143lnα)Ψ|δ3(r12)|Ψ
(13)
The atomic reduced mass, defined as μ = mpm/(mp + m), is used in this context. Equation 13 consists of five terms. The first term contains the so-called Bethe logarithm
lnk0=1DΨ|(a=12pa)(HE)ln[2(HE)](a=12pa)|Ψ
(14)
where
D=2πa=12b=34Ψ|δ3(rab)|Ψ
(15)
The next two expressions represent the electron–electron and electron–nucleus Araki–Sucher terms, respectively. These expectation values are defined as
Ψ|1r3|Ψ=limϵ0[Ψ|θ(rϵ)r3|Ψ+4π(γ+lnϵ)Ψ|δ3(r)|Ψ]
(16)
where the symbol γ denotes the Euler–Mascheroni constant, and θ is the Heaviside function. Finally, there are two terms in eq 13 that involve electron–electron and electron–nucleus Dirac deltas. Currently, only these two numerically dominant terms have been evaluated by using the four-body naJC wave function. The remaining three terms of eq 13 were estimated within the framework of the BO approximation, which significantly influenced the overall error budget. The Bethe logarithm and the electron–electron Araki–Sucher BO potentials were taken from ref (18), while the electron–nucleus Araki–Sucher potential was from ref (83).

2.4. Higher-Order QED Correction, E(6)

The higher-order QED contribution, calculated within the BO approximation using the wave function Ψ(r,R)=ϕ(r)χ(R), is denoted as E(6,0). The energy E(6,0) is expressed in terms of two internuclear potentials: the Breit–Pauli potential E(4,0)(R) and the higher-order relativistic potential E(6,0)(R), which are defined as follows (16)
E(4,0)(R)=ϕ|H(4,0)|ϕel
(17)
E(6,0)(R)=ϕ|H(6,0)|ϕel+ϕ|H(4,0)1(EelHel)H(4,0)|ϕel
(18)
The Hamiltonians H(4,0) and H(6,0) are the leading and higher-order relativistic Hamiltonians, respectively, both in the nonrecoil limit indicated by 0 in the superscript. This limit entails neglecting all terms that depend on the nuclear mass. The expectation value symbol ⟨···⟩el means integration over electronic variables only. HelH(2,0) and Hn represent the electronic and nuclear components of nonrelativistic Hamiltonian H(2), respectively. The explicit formulas for E(6,0)(R) are too extensive to be presented here─they can be found in ref (16). The total correction is obtained by integration over the nuclear variables
E(6,0)=χvJ|E(6,0)(R)|χvJ+χvJ|E(4,0)(R)1(EvJ(2,0)Hn)E(4,0)(R)|χvJ
(19)
where χ is the nuclear wave function obtained by solving the radial Schrödinger equation. More details on how to evaluate this correction can be found in the original article. (16)

2.5. Estimation of E(7)

The E(7) term is a particularly problematic correction, because its complete form is still unknown. It is therefore estimated by a dominating term inferred from the atomic hydrogen theory (80)
E(7)πΨ|a=12b=34δ3(rab)|Ψ{1π[A60+A61lnα2+A62ln2α2]+1π2B50+1π3C40}
(20)
This approximation adds another significant contribution to the error budget.

2.6. Atomic Limits

This section provides a comprehensive list of atomic thresholds relevant to the subsequent terms of eq 1. These atomic limits can be used in the calculation of dissociation energy Dv,J = 2 E(H) – E(H2) for a specific rovibrational level (v, J).
E(2)(H)=μ2m
(21)
E(4)(H)=μ38m3[1+3mmp+m2mp2]
(22)
Efs(4)(H)=2μ33m3rp2ƛ2
(23)
E(5)(H)=43πμm(lnk0H+lnμm)+76πmmp4ln2+43πμ3m3[19302lnα+mmp(31612lnα)]
(24)
E(6)(H)=116+(427962ln2)+(9ζ(3)4π22179648π2+3ln221027)
(25)
E(7)(H)=μ3m3[1π(A60+A61lnα2+A62ln2α2)+1π2B50+1π3C40]
(26)
ζ is the Riemann zeta function. The atomic Bethe logarithm is defined as follows
lnk0H=12p(HE)ln[2(HE)]p
(27)
with the atomic nonrelativistic Hamiltonian H = p2/2 – 1/r. The numerical value of the hydrogenic Bethe logarithm is known with high accuracy, (84) lnk0H=2.984128555765498. In the present work, eqs 2123 were employed in molecular calculations performed using the nonadiabatic wave function. The Bethe logarithm, in turn, was evaluated in the BO approximation, which means that the first term in eqs 13 and 24 is modified by taking the infinite nuclear mass limit. The expression E(6)(H) represents the atomic limit in the BO approximation, too. Three distinct components of eq 25 arise from the Dirac theory and one-loop and two-loop QED corrections, respectively. Finally, the factor μ3/m3 in E(7)(H) reflects the usage of the Dirac delta calculated with the fully nonadiabatic molecular wave function.

3. Results

Click to copy section linkSection link copied!

Numerical calculations were performed using the CODATA 2022-recommended physical constants. (80) Specifically, the following constants were used: the proton-to-electron mass ratio, mp/m = 1836.152 673 426(32), the inverse fine structure constant, 1/α = 137.035999177(21), the Rydberg constant, R=109,737.31568157(12)cm1, the proton charge radius, rp = 0.84075(64) fm, the reduced electron Compton wavelength, ƛ = 386.15926744(12) fm, and the speed of light in vacuum, c = 299,792,458 m s–1. To convert our numerical results from atomic units to wavenumbers, we multiply the E(n) by 2Rαn2, while conversion to MHz requires multiplication by 2×104Rcαn2.
In accordance with eq 1, the total energy of a rovibrational level is calculated as the sum of the nonrelativistic energy and successive corrections of increasing order in α, as outlined in the previous section. The convergence of individual relativistic expectation values, evaluated within the framework of the direct nonadiabatic approach, is demonstrated in the Supporting Information tables. Expectation values involving the electron–nucleus Dirac delta, which appear in corrections Efs(4), E(5), and E(7), can also be inferred from this data. Analyzing this convergence provides a foundation for estimating the uncertainties associated with individual corrections and the total energy. Furthermore, the results presented in these tables may serve as benchmarks for other calculation methods.

3.1. Dissociation Energy

Table 1 provides the final results of the dissociation energy Dv,J for the lowest vibrational (v = 0, 1, 2) and rotational (J = 0, ..., 4) quantum numbers. In addition to the total dissociation energy, the table includes all components of Dv,J that pertain to subsequent terms E(n) of the α-expansion. These components were calculated as 2 E(n)(H) – E(n)(H2). The numerical values of E(n)(H), discussed in Sec. 2.6, are considered exact, meaning they do not affect the uncertainty of Dv,J.
Table 1. Theoretically Predicted Dissociation Energy Dv,J for the Five Lowest Rotational Levels in the Vibrational States v = 0, 1, 2 of H2 [cm–1]b
componentD0,0D1,0D2,0
E(2)36,118.797744716(4)a31,957.63367420(2)28,031.79525303(2)
E(4)–0.531217263(61)–0.554783333(62)–0.572886764(67)
Efs(4)–0.000030900(33)–0.000027745(30)–0.000024797(27)
E(5)–0.19491021(15)–0.173629(10)–0.153814(20)
E(6)–0.0020577(66)–0.0018664(60)–0.0016868(54)
E(7)0.000101(25)0.000090(23)0.000081(20)
total36,118.069630(26)31,956.903458(26)28,031.066922(29)
componentD0,1D1,1D2,1
E(2)36,000.31248423(1)31,845.06061442(1)27,925.00480794(2)
E(4)–0.53380004(31)–0.55714556(38)–0.57502713(41)
Efs(4)–0.000030749(33)–0.000027603(30)–0.000024664(27)
E(5)–0.19388906(52)–0.172669(11)–0.152912(20)
E(6)–0.0020488(66)–0.0018581(60)–0.0016789(54)
E(7)0.000100(25)0.000090(22)0.000080(20)
total35,999.582816(26)31,844.329004(25)27,924.275245(29)
componentD0,2D1,2D2,2
E(2)35,764.42923080(1)31,620.96551340(1)27,712.44065580(2)
E(4)–0.53890776(38)–0.56181401(38)–0.57925386(41)
Efs(4)–0.000030451(33)–0.000027322(29)–0.000024399(26)
E(5)–0.1918636(15)–0.170765(12)–0.151123(21)
E(6)–0.0020312(65)–0.0018415(59)–0.0016633(53)
E(7)0.000099(25)0.000089(22)0.000079(20)
total35,763.696497(26)31,620.231155(26)27,711.708670(29)
componentD0,3D1,3D2,3
E(2)35,413.28801899(1)31,287.41651194(2)27,396.10305521(2)
E(4)–0.54642754(38)–0.56868046(38)–0.58546271(41)
Efs(4)–0.000030009(32)–0.000026906(29)–0.000024007(26)
E(5)–0.1888668(30)–0.167948(13)–0.148477(22)
E(6)–0.0020052(64)–0.0018169(58)–0.0016401(53)
E(7)0.000098(24)0.000088(22)0.000078(20)
total35,412.550787(25)31,286.678128(26)27,395.367529(30)
componentD0,4D1,4D2,4
E(2)34,950.01523845(2)30,847.43371054(2)26,978.91200753(2)
E(4)–0.55619744(35)–0.57758917(37)–0.59350370(40)
Efs(4)–0.000029431(32)–0.000026362(28)–0.000023495(25)
E(5)–0.1849460(49)–0.164264(15)–0.145018(24)
E(6)–0.0019711(63)–0.0017848(57)–0.0016098(52)
E(7)0.000096(24)0.000086(21)0.000076(19)
total34,949.272190(25)30,846.690132(26)26,978.171929(31)
a

This energy differs from that in our previous studies (11,12,74) due to differences in the physical constants used.

b

Components of Dv,J related to subsequent terms E(n) of the α-expansion (1) are also included.

As mentioned in the previous section, there are three components of the leading QED correction, E(5), that were determined within the BO approximation. To enhance the precision of this correction’s final value, we utilized the precise values of these components obtained for the ground level (0,0) from the direct nonadiabatic calculations with the naECG wave function. (11) For all levels (v, J), we applied a uniform shift to these three components based on the difference between the nonadiabatic value and the value obtained from the BO approximation. Thus, the shifted contribution to the dissociation energy from the component Q can be expressed as Qv,JBOQ0,0BO + Q0,0naECG. The total shift from all three components is −1.1 × 10–5 cm–1. Because this is a uniform shift applied to all energy levels, it does not alter the spacing between them.
The uncertainties related to E(2) and E(4) are purely numerical; they arise from the finite size of the naJC basis set and have a negligible influence on the overall error budget. The most considerable contributions to the uncertainty of Dv,J come from the E(5) and E(7) terms. In the first case, the primary source of error is the omission of finite mass effects in the three aforementioned components. Their uncertainties were estimated as (Qv,JBOQ0,0BO)/μn, where μn = mp/2 is the nuclear reduced mass. In the latter case, the error comes from the incompleteness of the expression for E(7). In summary, the uncertainties associated with Dv,J are less than 3.3 × 10–5 cm–1(≈1 MHz).

3.2. Rovibrational Transition Energy

Establishing reliable uncertainty for a transition frequency ν is challenging due to the partial error cancellation between two rovibrational levels. The overall uncertainty u is calculated as the root-mean-square of the individual uncertainties u(n), which arise from the terms ν(n) = ΔE(n) in the α-expansion. To determine the u(n) values, we applied distinct approaches depending on how each ν(n) term was evaluated. For ν(n) obtained with the direct nonadiabatic approach (as ν(2) and ν(4)), the uncertainty was estimated as the maximum uncertainty of both levels. In contrast, for ν(n) calculated using the BO approximation (as ν(5) and ν(6)), the uncertainty was estimated by the formula ν(n)n. This scaling accounts for the missing finite nuclear mass effects. Since the complete ν(7) is currently unknown, and its dominant term from the atomic hydrogen theory only provides an estimate, we assumed a conservative error margin of 25% for this contribution. The uncertainty of νfs(4) is determined by that of the proton charge radius and is relatively small compared to the main error components.
Having determined the uncertainties, we are now in a position to compare our results with experimental values─a crucial step in validating the accuracy and reliability of the theoretical model. Among the numerous rovibrational intervals measured for H2, only nine are known with sub-MHz precision. These measurements, taken within the last two years, (69−71,73,85) encompass both purely rotational transitions and those belonging to the fundamental vibrational band as well as the first overtone band. Table 2 presents a compilation of these experimental frequencies alongside their theoretical counterparts. Each panel in the table is dedicated to a specific transition line and includes all components that contribute to the final theoretical frequency along with their individual error bars. The bottom row of each panel displays the difference between the theoretical and experimental frequencies, accompanied by the combined uncertainty, σ. We observe agreement for all nine lines with the largest difference being 1.4 σ for the Q1(1) line. Figure 1 illustrates these differences relative to the ±1 MHz band, clearly demonstrating that the theoretical predictions and spectroscopic measurements are in sub-MHz agreement. In some cases, the experimental results are an order of magnitude more accurate than the theoretical predictions, which highlights the need for further development of computational methods, particularly for ν(7).
Table 2. Comparison of Theoretically Predicted Transition Frequencies (in MHz) for H2 with Nine Experimental Results That Were Measured with Sub-MHz Accuracyf
componentS0(0): (0, 2) → (0, 0)S0(1): (0, 3) → (0, 1)Q1(1): (1, 1) → (0, 1)
ν(2)10,623,700.7825(3)17,598,550.7339(4)124,571,317.1657(4)
ν(4)230.555(11)378.562(11)699.881(11)
νfs(4)–0.01346(2)–0.02219(3)–0.0943(1)
ν(5)–91.335(47)–150.565(77)–636.16(33)
ν(6)–0.793(3)–1.308(4)–5.719(18)
ν(7)0.040(10)0.070(17)0.310(78)
theory10,623,839.229(49)17,598,777.469(80)124,571,375.38(34)
experiment10,623,839.09(39)a17,598,777.46(75)a124,571,374.73(31)b
difference0.14(39)0.01(75)0.65(46)
componentS1(0): (1, 2) → (0, 0)Q2(1): (2, 1) → (0, 1)Q2(2): (2, 2) → (0, 2)
ν(2)134,841,618.0297(4)242,091,633.7378(5)241,392,544.6686(5)
ν(4)917.267(11)1235.957(12)1209.546(12)
νfs(4)–0.1072(2)–0.1824(3)–0.1814(3)
ν(5)–723.86(37)–1228.46(63)–1221.38(63)
ν(6)–6.482(21)–11.089(36)–11.031(35)
ν(7)0.340(85)0.59(15)0.59(15)
theory134,841,805.19(38)242,091,630.55(65)241,392,522.22(65)
experiment134,841,805.102(15)c242,091,630.140(9)d241,392,522.00(34)a
difference0.09(38)0.41(65)0.22(73)
componentQ2(3): (2, 3) → (0, 3)S2(0): (2, 2) → (0, 0)S2(1): (2, 3) → (0, 1)
ν(2)240,349,158.6533(5)252,016,245.4511(5)257,947,709.3872(5)
ν(4)1170.245(12)1440.101(12)1548.807(12)
νfs(4)–0.1799(3)–0.1949(3)–0.2021(3)
ν(5)–1210.85(62)–1312.71(68)–1361.42(70)
ν(6)–10.945(35)–11.824(38)–12.253(39)
ν(7)0.59(15)0.63(16)0.66(17)
theory240,349,107.52(64)25,2016,361.45(69)25,7947,884.98(72)
experiment240,349,107.15(71)a25,2016,361.164(8)e25,7947,884.597(30)a
difference0.36(96)0.28(69)0.39(72)
a

Fleurbaey et al., 2023. (69)

b

Lamperti et al., 2023. (70)

c

Stankiewicz et al., 2025. (73)

d

Diouf et al., 2024. (85)

e

Cozijn et al., 2023. (71)

f

The components ν(n) = ΔE(n) of the theoretical frequency, corresponding to subsequent terms of the α-expansion (1), are also included. CODATA 2022-recommended physical constants were used. (80)

Figure 1

Figure 1. Differences between theoretically predicted and experimental line positions shown against the ±1 MHz error band (in blue) and the theoretical error band (in green). The error bars around the dots come from the experiment.

Unlike earlier calculations (86) that estimated the relativistic and quantum electrodynamic corrections using the Born–Oppenheimer approximation, the current direct nonadiabatic results include the finite nuclear mass effects in both the relativistic and QED corrections. For frequencies in the first overtone, the difference between current and previous calculations reaches 3 MHz; in the fundamental band, it is approximately 1.5 MHz, and for purely rotational transitions, it is about 0.1 MHz. The small third value demonstrates the effective cancellation of recoil effects within a single oscillation band. The change in the total recoil effect is primarily due to relativistic and leading QED components of the frequency.
In addition to the nine transition frequencies listed in Table 2, other energy gaps of interest for future measurements can be evaluated by calculating the difference in dissociation energies from Table 1. An upper bound for the uncertainty associated with these transitions can be inferred from Table 2: 0.1 MHz for rotational transitions, 0.4 MHz for the fundamental band, and 0.8 MHz for the first overtone. The energies for both the levels and transitions will be included in the new release of the publicly accessible H2Spectre program. (86)

4. Conclusion

Click to copy section linkSection link copied!

We have shown that the relative uncertainty of the theoretically predicted rovibrational intervals in H2 is about three parts per billion. This indicates that the calculated vibrational transition frequencies have an uncertainty at the tenth significant figure. The results obtained demonstrate that finite nuclear mass effects significantly impact not only nonrelativistic energy but also relativistic and QED corrections. By incorporation of these effects, the theoretical predictions matched the experimental results with sub-MHz accuracy.
Attaining the current level of experimental accuracy is essential for testing QED within molecular systems as well as for interpreting any discrepancies in terms of hypothetical unknown forces. There are only a few instances of physical systems where such a level of agreement between molecular or atomic spectroscopy and first-principles calculations has been achieved. Examples include the transitions in hydrogen (87,88) and helium atoms, (89) as well as those in the hydrogen molecular ion. (90)
Further advancements will require a fully nonadiabatic calculation of the Bethe logarithm and an Araki–Sucher correction for excited rovibrational levels. We believe this task is feasible since it has already been accomplished at the ground rovibrational level. (11)

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.5c01702.

  • Convergence of relativistic expectation values for individual rovibrational levels of H2 (PDF)

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Authors
    • Notes
      The authors declare no competing financial interest.

    Acknowledgments

    Click to copy section linkSection link copied!

    This research was supported by the National Science Center (Poland) Grant No. 2021/41/B/ST4/00089. A computer grant from the Poznań Supercomputing and Networking Center was used to carry out the numerical calculations.

    References

    Click to copy section linkSection link copied!

    This article references 90 other publications.

    1. 1
      Heitler, W.; London, F. Wechselwirkung neutraler Atome und homeopolare Bindung nach der Quantenmechanik. Z. Phys. 1927, 44, 455,  DOI: 10.1007/BF01397394
    2. 2
      James, H. M.; Coolidge, A. S. The Ground State of the Hydrogen Molecule. J. Chem. Phys. 1933, 1, 825,  DOI: 10.1063/1.1749252
    3. 3
      Kolos, W.; Roothaan, C. C. J. Accurate Electronic Wave Functions for the H2 Molecule. Rev. Mod. Phys. 1960, 32, 219232,  DOI: 10.1103/RevModPhys.32.219
    4. 4
      Kołos, W.; Wolniewicz, L. The coupling between electronic and nuclear motion and the relativistic effects in the ground state of the H2 molecule. Acta Phys. Polonica 1961, 20, 129140
    5. 5
      Kołos, W.; Wolniewicz, L. Accurate Adiabatic Treatment of the Ground State of the Hydrogen Molecule. J. Chem. Phys. 1964, 41, 36633673,  DOI: 10.1063/1.1725796
    6. 6
      Kołos, W.; Wolniewicz, L. Accurate Computation of Vibronic Energies and of Some Expectation Values for H2, D2, and T2. J. Chem. Phys. 1964, 41, 36743678,  DOI: 10.1063/1.1725797
    7. 7
      Pachucki, K.; Komasa, J. Nonadiabatic corrections to rovibrational levels of H2. J. Chem. Phys. 2009, 130, 164113,  DOI: 10.1063/1.3114680
    8. 8
      Pachucki, K.; Komasa, J. Leading order nonadiabatic corrections to rovibrational levels of H2, D2, and T2. J. Chem. Phys. 2015, 143, 034111,  DOI: 10.1063/1.4927079
    9. 9
      Kołos, W.; Wolniewicz, L. Nonadiabatic Theory for Diatomic Molecules and Its Application to the Hydrogen Molecule. Rev. Mod. Phys. 1963, 35, 473,  DOI: 10.1103/RevModPhys.35.473
    10. 10
      Kinghorn, D. B.; Adamowicz, L. Improved Nonadiabatic Ground-State Energy Upper Bound for Dihydrogen. Phys. Rev. Lett. 1999, 83, 25412543,  DOI: 10.1103/PhysRevLett.83.2541
    11. 11
      Puchalski, M.; Komasa, J.; Czachorowski, P.; Pachucki, K. Nonadiabatic QED Correction to the Dissociation Energy of the Hydrogen Molecule. Phys. Rev. Lett. 2019, 122, 103003,  DOI: 10.1103/PhysRevLett.122.103003
    12. 12
      Pachucki, K.; Komasa, J. Nonadiabatic rotational states of the hydrogen molecule. Phys. Chem. Chem. Phys. 2018, 20, 247255,  DOI: 10.1039/C7CP06516G
    13. 13
      Wolniewicz, L. Relativistic energies of the ground state of the hydrogen molecule. J. Chem. Phys. 1993, 99, 1851,  DOI: 10.1063/1.465303
    14. 14
      Piszczatowski, K.; Łach, G.; Przybytek, M.; Komasa, J.; Pachucki, K.; Jeziorski, B. Theoretical Determination of the Dissociation Energy of Molecular Hydrogen. J. Chem. Theory Comput. 2009, 5, 3039,  DOI: 10.1021/ct900391p
    15. 15
      Komasa, J.; Piszczatowski, K.; Łach, G.; Przybytek, M.; Jeziorski, B.; Pachucki, K. Quantum Electrodynamics Effects in Rovibrational Spectra of Molecular Hydrogen. J. Chem. Theory Comput. 2011, 7, 31053115,  DOI: 10.1021/ct200438t
    16. 16
      Puchalski, M.; Komasa, J.; Czachorowski, P.; Pachucki, K. Complete α6 m Corrections to the Ground State of H2. Phys. Rev. Lett. 2016, 117, 263002,  DOI: 10.1103/PhysRevLett.117.263002
    17. 17
      Czachorowski, P.; Puchalski, M.; Komasa, J.; Pachucki, K. Nonadiabatic relativistic correction in H2, D2, and HD. Phys. Rev. A 2018, 98, 052506,  DOI: 10.1103/PhysRevA.98.052506
    18. 18
      Siłkowski, M.; Pachucki, K.; Komasa, J.; Puchalski, M. Leading-order QED effects in the ground electronic state of molecular hydrogen. Phys. Rev. A 2023, 107, 032807,  DOI: 10.1103/PhysRevA.107.032807
    19. 19
      Wang, L. M.; Yan, Z.-C. Relativistic corrections to the ground state of H2 calculated without using the Born-Oppenheimer approximation. Phys. Rev. A 2018, 97, 060501,  DOI: 10.1103/PhysRevA.97.060501
    20. 20
      Puchalski, M.; Moszynski, R.; Komasa, J. In Polish Quantum Chemistry from Kołos to Now. In Advances in Quantum Chemistry; Musiał, M., Grabowski, I., Eds.; Academic Press, 2023; Vol. 87, pp 136.
    21. 21
      Pachucki, K.; Komasa, J. Magnetic dipole transitions in the hydrogen molecule. Phys. Rev. A 2011, 83, 032501,  DOI: 10.1103/PhysRevA.83.032501
    22. 22
      Margolis, J.; Hunt, G. On the level of H2 quadrupole absorption in the Jovian atmosphere. Icarus 1973, 18, 593598,  DOI: 10.1016/0019-1035(73)90061-4
    23. 23
      Roueff, E.; Abgrall, H.; Czachorowski, P.; Pachucki, K.; Puchalski, M.; Komasa, J. The full infrared spectrum of molecular hydrogen. Astron. Astrophys. 2019, 630, A58,  DOI: 10.1051/0004-6361/201936249
    24. 24
      Tchernyshyov, K. A Detection of H2 in a High-velocity Cloud toward the Large Magellanic Cloud. Astrophys. J. 2022, 931, 78,  DOI: 10.3847/1538-4357/ac68e0
    25. 25
      Jacob, A. Small molecules, big impact: a tale of hydrides past, present, and future. Astrophys. Space Sci. 2023, 368, 76,  DOI: 10.1007/s10509-023-04229-8
    26. 26
      Tennyson, J. The 2024 release of the ExoMol database: Molecular line lists for exoplanet and other hot atmospheres. J. Quant. Spectrosc. Radiat. Transfer 2024, 326, 109083,  DOI: 10.1016/j.jqsrt.2024.109083
    27. 27
      Bruno, D. Uncertainty quantification of a collisional-radiative model for hydrogen plasmas. Phys. Scr. 2025, 100, 035611,  DOI: 10.1088/1402-4896/adb227
    28. 28
      Zúñiga, J.; Bastida, A.; Requena, A.; Cerezo, J. Ab Initio Partition Functions and Thermodynamic Quantities for the Molecular Hydrogen Isotopologues. J. Phys. Chem. A 2021, 125, 92269241,  DOI: 10.1021/acs.jpca.1c06468
    29. 29
      Leachman, J. W.; Wilhelmsen, Ø.; Matveev, K. I. Cool Fuel: The Science and Engineering of Cryogenic Hydrogen; Oxford University Press, 2025.
    30. 30
      Puchalski, M.; Komasa, J.; Pachucki, K. Hyperfine Structure of the First Rotational Level in H2, D2 and HD Molecules and the Deuteron Quadrupole Moment. Phys. Rev. Lett. 2020, 125, 253001,  DOI: 10.1103/PhysRevLett.125.253001
    31. 31
      Puchalski, M.; Komasa, J.; Spyszkiewicz, A.; Pachucki, K. Nuclear magnetic shielding in HD and HT. Phys. Rev. A 2022, 105, 042802,  DOI: 10.1103/PhysRevA.105.042802
    32. 32
      Kassi, S.; Campargue, A. Electric quadrupole and dipole transitions of the first overtone band of HD by CRDS between 1.45 and 1.33 μm. J. Mol. Spectrosc. 2011, 267, 36,  DOI: 10.1016/j.jms.2011.02.001
    33. 33
      Campargue, A.; Kassi, S.; Pachucki, K.; Komasa, J. The absorption spectrum of H2: CRDS measurements of the (2–0) band, review of the literature data and accurate ab initio line list up to 35000 cm–1. Phys. Chem. Chem. Phys. 2012, 14, 802,  DOI: 10.1039/C1CP22912E
    34. 34
      Kassi, S.; Campargue, A.; Pachucki, K.; Komasa, J. The absorption spectrum of D2: Ultrasensitive cavity ring down spectroscopy of the (2–0) band near 1.7 μm and accurate ab initio line list up to 24 000 cm–1. J. Chem. Phys. 2012, 136, 184309,  DOI: 10.1063/1.4707708
    35. 35
      Cheng, C.-F.; Sun, Y. R.; Pan, H.; Wang, J.; Liu, A.-W.; Campargue, A.; Hu, S.-M. Electric-quadrupole transition of H2 determined to 10–9 precision. Phys. Rev. A 2012, 85, 024501,  DOI: 10.1103/PhysRevA.85.024501
    36. 36
      Dickenson, G. D.; Niu, M. L.; Salumbides, E. J.; Komasa, J.; Eikema, K. S. E.; Pachucki, K.; Ubachs, W. Fundamental Vibration of Molecular Hydrogen. Phys. Rev. Lett. 2013, 110, 193601,  DOI: 10.1103/PhysRevLett.110.193601
    37. 37
      Niu, M.; Salumbides, E.; Dickenson, G.; Eikema, K.; Ubachs, W. Precision spectroscopy of the rovibrational splittings in H2, HD and D2. J. Mol. Spectrosc. 2014, 300, 4454,  DOI: 10.1016/j.jms.2014.03.011
    38. 38
      Tan, Y.; Wang, J.; Cheng, C.-F.; Zhao, X.-Q.; Liu, A.-W.; Hu, S.-M. Cavity ring-down spectroscopy of the electric quadrupole transitions of in the 784–852 nm region. J. Mol. Spectrosc. 2014, 300, 6064,  DOI: 10.1016/j.jms.2014.03.010
    39. 39
      Kassi, S.; Campargue, A. Electric quadrupole transitions and collision-induced absorption in the region of the first overtone band of H2 near 1.25μm. J. Mol. Spectrosc. 2014, 300, 5559,  DOI: 10.1016/j.jms.2014.03.022
    40. 40
      Mondelain, D.; Kassi, S.; Sala, T.; Romanini, D.; Gatti, D.; Campargue, A. Sub-MHz accuracy measurement of the S(2) 2–0 transition frequency of D2 by Comb-Assisted Cavity Ring Down spectroscopy. J. Mol. Spectrosc. 2016, 326, 58,  DOI: 10.1016/j.jms.2016.02.008
    41. 41
      Wcisło, P.; Thibault, F.; Zaborowski, M.; Wójtewicz, S.; Cygan, A.; Kowzan, G.; Masłowski, P.; Komasa, J.; Puchalski, M.; Pachucki, K.; Ciuryło, R.; Lisak, D. Accurate deuterium spectroscopy for fundamental studies. J. Quant. Spectrosc. Radiat. Transfer 2018, 213, 4151,  DOI: 10.1016/j.jqsrt.2018.04.011
    42. 42
      Trivikram, T. M.; Schlösser, M.; Ubachs, W.; Salumbides, E. J. Relativistic and QED Effects in the Fundamental Vibration of T2. Phys. Rev. Lett. 2018, 120, 163002,  DOI: 10.1103/PhysRevLett.120.163002
    43. 43
      Tao, L.-G.; Liu, A.-W.; Pachucki, K.; Komasa, J.; Sun, Y. R.; Wang, J.; Hu, S.-M. Toward a Determination of the Proton-Electron Mass Ratio from the Lamb-Dip Measurement of HD. Phys. Rev. Lett. 2018, 120, 153001,  DOI: 10.1103/PhysRevLett.120.153001
    44. 44
      Fasci, E.; Castrillo, A.; Dinesan, H.; Gravina, S.; Moretti, L.; Gianfrani, L. Precision spectroscopy of HD at 1.38μm. Phys. Rev. A 2018, 98, 022516,  DOI: 10.1103/PhysRevA.98.022516
    45. 45
      Martinez, R. Z.; Bermejo, D.; Wcisło, P.; Thibault, F. Accurate wavenumber measurements for the S0(0), S0(1), and S0(2) pure rotational Raman lines of D2. J. Raman Spectrosc. 2019, 50, 127129,  DOI: 10.1002/jrs.5499
    46. 46
      Hölsch, N.; Beyer, M.; Salumbides, E. J.; Eikema, K. S. E.; Ubachs, W.; Jungen, C.; Merkt, F. Benchmarking Theory with an Improved Measurement of the Ionization and Dissociation Energies of H2. Phys. Rev. Lett. 2019, 122, 103002,  DOI: 10.1103/PhysRevLett.122.103002
    47. 47
      Beyer, M.; Hölsch, N.; Hussels, J.; Cheng, C.-F.; Salumbides, E. J.; Eikema, K. S. E.; Ubachs, W.; Jungen, C.; Merkt, F. Determination of the Interval between the Ground States of Para- and Ortho-H2. Phys. Rev. Lett. 2019, 123, 163002,  DOI: 10.1103/PhysRevLett.123.163002
    48. 48
      Lai, K.-F.; Czachorowski, P.; Schlösser, M.; Puchalski, M.; Komasa, J.; Pachucki, K.; Ubachs, W.; Salumbides, E. J. Precision tests of nonadiabatic perturbation theory with measurements on the DT molecule. Phys. Rev. Res. 2019, 1, 033124,  DOI: 10.1103/PhysRevResearch.1.033124
    49. 49
      Lai, K.-F.; Hermann, V.; Trivikram, M.; Diouf, M.; Schlösser, M.; Ubachs, W.; Salumbides, E. Precision measurement of the fundamental vibrational frequencies of tritium-bearing hydrogen molecules: T2, DT, HT. Phys. Chem. Chem. Phys. 2020, 22, 89738987,  DOI: 10.1039/D0CP00596G
    50. 50
      Mondelain, D.; Kassi, S.; Campargue, A. Transition frequencies in the (2–0) band of D2 with MHz accuracy. J. Quant. Spectrosc. Radiat. Transfer 2020, 253, 107020,  DOI: 10.1016/j.jqsrt.2020.107020
    51. 51
      Jóźwiak, H.; Cybulski, H.; Wcisło, P. Hyperfine structure of quadrupole rovibrational transitions in tritium-bearing hydrogen isotopologues. J. Quant. Spectrosc. Radiat. Transfer 2020, 256, 107255,  DOI: 10.1016/j.jqsrt.2020.107255
    52. 52
      Kassi, S.; Fleurbaey, H.; Campargue, A. First detection and absolute transition frequencies in the (3–0) band of D2. J. Chem. Phys. 2024, 160, 094306,  DOI: 10.1063/5.0196903
    53. 53
      Salumbides, E. J.; Dickenson, G. D.; Ivanov, T. I.; Ubachs, W. QED Effects in Molecules: Test on Rotational Quantum States of H2. Phys. Rev. Lett. 2011, 107, 043005,  DOI: 10.1103/PhysRevLett.107.043005
    54. 54
      Salumbides, E. J.; Koelemeij, J. C. J.; Komasa, J.; Pachucki, K.; Eikema, K. S. E.; Ubachs, W. Bounds on fifth forces from precision measurements on molecules. Phys. Rev. D 2013, 87, 112008,  DOI: 10.1103/PhysRevD.87.112008
    55. 55
      Salumbides, E. J.; Schellekens, A. N.; Gato-Rivera, B.; Ubachs, W. Constraints on extra dimensions from precision molecular spectroscopy. New J. Phys. 2015, 17, 033015,  DOI: 10.1088/1367-2630/17/3/033015
    56. 56
      Ubachs, W.; Koelemeij, J.; Eikema, K.; Salumbides, E. Physics beyond the Standard Model from hydrogen spectroscopy. J. Mol. Spectrosc. 2016, 320, 112,  DOI: 10.1016/j.jms.2015.12.003
    57. 57
      Drouin, B. J.; Yu, S.; Pearson, J. C.; Gupta, H. Terahertz spectroscopy for space applications: 2.5–2.7 THz spectra of HD, H2O and NH3. J. Mol. Spectrosc. 2011, 1006, 212,  DOI: 10.1016/j.molstruc.2011.05.062
    58. 58
      Cozijn, F. M. J.; Dupré, P.; Salumbides, E. J.; Eikema, K. S. E.; Ubachs, W. Sub-Doppler Frequency Metrology in HD for Tests of Fundamental Physics. Phys. Rev. Lett. 2018, 120, 153002,  DOI: 10.1103/PhysRevLett.120.153002
    59. 59
      Diouf, M. L.; Cozijn, F. M. J.; Darquié, B.; Salumbides, E. J.; Ubachs, W. Lamb-dips and Lamb-peaks in the saturation spectrum of HD. Opt. Lett. 2019, 44, 47334736,  DOI: 10.1364/OL.44.004733
    60. 60
      Diouf, M. L.; Cozijn, F. M. J.; Lai, K.-F.; Salumbides, E. J.; Ubachs, W. Lamb-peak spectrum of the HD (2–0) P(1) line. Phys. Rev. Res. 2020, 2, 023209,  DOI: 10.1103/PhysRevResearch.2.023209
    61. 61
      Hua, T.-P.; Sun, Y. R.; Hu, S.-M. Dispersion-like lineshape observed in cavity-enhanced saturation spectroscopy of HD at 1.4 μm. Opt. Lett. 2020, 45, 48634866,  DOI: 10.1364/OL.401879
    62. 62
      Wójtewicz, S.; Gotti, R.; Gatti, D.; Lamperti, M.; Laporta, P.; Jóźwiak, H.; Thibault, F.; Wcisło, P.; Marangoni, M. Accurate deuterium spectroscopy and comparison with ab initio calculations. Phys. Rev. A 2020, 101, 052504,  DOI: 10.1103/PhysRevA.101.052504
    63. 63
      Fast, A.; Meek, S. A. Sub-ppb Measurement of a Fundamental Band Rovibrational Transition in HD. Phys. Rev. Lett. 2020, 125, 023001,  DOI: 10.1103/PhysRevLett.125.023001
    64. 64
      Castrillo, A.; Fasci, E.; Gianfrani, L. Doppler-limited precision spectroscopy of HD at 1.4 μm: An improved determination of the R(1) center frequency. Phys. Rev. A 2021, 103, 022828,  DOI: 10.1103/PhysRevA.103.022828
    65. 65
      Fast, A.; Meek, S. A. Precise measurement of the D2 S1(0) vibrational transition frequency. Mol. Phys. 2022, 120, e1999520  DOI: 10.1080/00268976.2021.1999520
    66. 66
      Kassi, S.; Lauzin, C.; Chaillot, J.; Campargue, A. The (2–0) R(0) and R(1) transition frequencies of HD determined to a 10–10 relative accuracy by Doppler spectroscopy at 80 K. Phys. Chem. Chem. Phys. 2022, 24, 2316423172,  DOI: 10.1039/D2CP02151J
    67. 67
      Cozijn, F. M. J.; Diouf, M. L.; Hermann, V.; Salumbides, E. J.; Schlösser, M.; Ubachs, W. Rotational level spacings in HD from vibrational saturation spectroscopy. Phys. Rev. A 2022, 105, 062823,  DOI: 10.1103/PhysRevA.105.062823
    68. 68
      Cozijn, F. M. J.; Diouf, M. L.; Ubachs, W. Saturation spectroscopy of R(0), R(2) and P(2) lines in the (2–0) band of HD. Eur. Phys. J. D 2022, 76, 220,  DOI: 10.1140/epjd/s10053-022-00552-x
    69. 69
      Fleurbaey, H.; Koroleva, A. O.; Kassi, S.; Campargue, A. The high-accuracy spectroscopy of H2 rovibrational transitions in the (2–0) band near 1.2 μm. Phys. Chem. Chem. Phys. 2023, 25, 1474914756,  DOI: 10.1039/D3CP01136D
    70. 70
      Lamperti, M.; Rutkowski, L.; Ronchetti, D.; Gatti, D.; Gotti, R.; Cerullo, G.; Thibault, F.; Jóźwiak, H.; Wójtewicz, S.; Masłowski, P.; Wcisło, P.; Polli, D.; Marangoni, M. Stimulated Raman scattering metrology of molecular hydrogen. Commun. Phys. 2023, 6, 67,  DOI: 10.1038/s42005-023-01187-z
    71. 71
      Cozijn, F. M. J.; Diouf, M. L.; Ubachs, W. Lamb Dip of a Quadrupole Transition in H2. Phys. Rev. Lett. 2023, 131, 073001,  DOI: 10.1103/PhysRevLett.131.073001
    72. 72
      Cozijn, F. M. J.; Diouf, M. L.; Ubachs, W.; Hermann, V.; Schlösser, M. Precision Measurement of Vibrational Quanta in Tritium Hydride. Phys. Rev. Lett. 2024, 132, 113002,  DOI: 10.1103/PhysRevLett.132.113002
    73. 73
      Stankiewicz, K.; Makowski, M.; Słowinski, M.; Sołtys, K. L.; Bednarski, B.; Jóźwiak, H.; Stolarczyk, N.; Narożnik, M.; Kierski, D.; Wójtewicz, S.; Cygan, A.; Kowzan, G.; Masłowski, P.; Piwinski, M.; Lisak, D.; Wcisło, P. Cavity-enhanced spectroscopy in the deep cryogenic regime – new hydrogen technologies for quantum sensing. arXiv 2025, arXiv:2502.12703,  DOI: 10.48550/arXiv.2502.12703
    74. 74
      Pachucki, K.; Komasa, J. Schrödinger equation solved for the hydrogen molecule with unprecedented accuracy. J. Chem. Phys. 2016, 144, 164306,  DOI: 10.1063/1.4948309
    75. 75
      Pachucki, K.; Komasa, J. Integrals for relativistic nonadiabatic energies of H2 in an exponential basis. Phys. Rev. A 2024, 109, 032822,  DOI: 10.1103/PhysRevA.109.032822
    76. 76
      Pachucki, K.; Komasa, J. Nonrelativistic energy levels of HD. Phys. Chem. Chem. Phys. 2018, 20, 2629726302,  DOI: 10.1039/C8CP05493B
    77. 77
      Pachucki, K.; Komasa, J. Nonrelativistic energy levels of D2. Phys. Chem. Chem. Phys. 2019, 21, 1027210276,  DOI: 10.1039/C9CP01308C
    78. 78
      Pachucki, K.; Komasa, J. Nonrelativistic energy of tritium-containing hydrogen molecule isotopologues. Mol. Phys. 2022, 120, e2040627  DOI: 10.1080/00268976.2022.2040627
    79. 79
      Hylleraas, E. Neue Berechnung der Energie des Heliums im Grundzustande, sowie des tiefsten Terms von Ortho-Helium. Z. Physik 1929, 54, 347366,  DOI: 10.1007/BF01375457
    80. 80
      Mohr, P. J.; Newell, D. B.; Taylor, B. N.; Tiesinga, E. CODATA recommended values of the fundamental physical constants: 2022. Rev. Mod. Phys. 2025, 97, 025002,  DOI: 10.1103/RevModPhys.97.025002
    81. 81
      Pachucki, K.; Komasa, J. Relativistic Correction from the Four-Body Nonadiabatic Exponential Wave Function. J. Chem. Theory Comput. 2024, 20, 86448651,  DOI: 10.1021/acs.jctc.4c00861
    82. 82
      Code, R. F.; Ramsey, N. F. Molecular-Beam Magnetic Resonance Studies of HD and D2. Phys. Rev. A 1971, 4, 19451959,  DOI: 10.1103/PhysRevA.4.1945
    83. 83
      Komasa, J. Electron-nucleus Araki-Sucher QED correction to rovibrational levels of H2. arXiv 2025, arXiv:2511.07053v2.
    84. 84
      Drake, G. W. F.; Swainson, R. A. Bethe logarithms for hydrogen up to n = 20, and approximations for two-electron atoms. Phys. Rev. A 1990, 41, 12431246,  DOI: 10.1103/PhysRevA.41.1243
    85. 85
      Diouf, M. L.; Cozijn, F. M. J.; Ubachs, W. Hyperfine structure in a vibrational quadrupole transition of ortho-H2. Mol. Phys. 2024, 122, e2304101  DOI: 10.1080/00268976.2024.2304101
    86. 86
      H2Spectre ver. 7.4 Fortran source code. 2022, https://www.fuw.edu.pl/~krp/, https://qcg.home.amu.edu.pl/H2Spectre.html (accessed Oct 10, 2025).
    87. 87
      Parthey, C. G.; Matveev, A.; Alnis, J.; Bernhardt, B.; Beyer, A.; Holzwarth, R.; Maistrou, A.; Pohl, R.; Predehl, K.; Udem, T.; Wilken, T.; Kolachevsky, N.; Abgrall, M.; Rovera, D.; Salomon, C.; Laurent, P.; Hänsch, T. W. Improved Measurement of the Hydrogen 1S–2S Transition Frequency. Phys. Rev. Lett. 2011, 107, 203001,  DOI: 10.1103/PhysRevLett.107.203001
    88. 88
      Pachucki, K.; Jentschura, U. D. Two-Loop Bethe-Logarithm Correction in Hydrogenlike Atoms. Phys. Rev. Lett. 2003, 91, 113005,  DOI: 10.1103/PhysRevLett.91.113005
    89. 89
      Clausen, G.; Scheidegger, S.; Agner, J. A.; Schmutz, H.; Merkt, F. Imaging-Assisted Single-Photon Doppler-Free Laser Spectroscopy and the Ionization Energy of Metastable Triplet Helium. Phys. Rev. Lett. 2023, 131, 103001,  DOI: 10.1103/PhysRevLett.131.103001
    90. 90
      Alighanbari, S.; Schenkel, M.; Korobov, V.; Schiller, S. High-accuracy laser spectroscopy of and the proton–electron mass ratio. Nature 2025, 644, 6975,  DOI: 10.1038/s41586-025-09306-2

    Cited By

    Click to copy section linkSection link copied!

    This article is cited by 3 publications.

    1. K. Stankiewicz, M. Makowski, M. Słowiński, K. L. Sołtys, B. Bednarski, H. J. Jóźwiak, N. Stolarczyk, M. Narożnik, D. Kierski, S. Wójtewicz, A. Cygan, G. Kowzan, P. Masłowski, M. Piwiński, D. Lisak, P. Wcisło. Cavity-enhanced spectroscopy in the deep cryogenic regime for quantum sensing and metrology. Nature Physics 2026, 47 https://doi.org/10.1038/s41567-026-03204-8
    2. Michał Siłkowski, Paweł Czachorowski, Jacek Komasa, Mariusz Puchalski. Progress toward spectroscopic accuracy: Relativistic and QED corrections for HeH + . Physical Review A 2026, 113 (3) https://doi.org/10.1103/k4qf-thhm
    3. Jacek Komasa. Electron-nucleus Araki-Sucher correction in the hydrogen molecule isotopologues. Physical Review A 2026, 113 (1) https://doi.org/10.1103/qyz4-n1tc

    Journal of Chemical Theory and Computation

    Cite this: J. Chem. Theory Comput. 2025, 21, 24, 12664–12673
    Click to copy citationCitation copied!
    https://doi.org/10.1021/acs.jctc.5c01702
    Published December 5, 2025

    Copyright © 2025 American Chemical Society. This publication is licensed under these Terms of Use.

    Article Views

    1985

    Altmetric

    -

    Citations

    Learn about these metrics

    Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

    Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

    The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

    • Abstract

      Figure 1

      Figure 1. Differences between theoretically predicted and experimental line positions shown against the ±1 MHz error band (in blue) and the theoretical error band (in green). The error bars around the dots come from the experiment.

    • References


      This article references 90 other publications.

      1. 1
        Heitler, W.; London, F. Wechselwirkung neutraler Atome und homeopolare Bindung nach der Quantenmechanik. Z. Phys. 1927, 44, 455,  DOI: 10.1007/BF01397394
      2. 2
        James, H. M.; Coolidge, A. S. The Ground State of the Hydrogen Molecule. J. Chem. Phys. 1933, 1, 825,  DOI: 10.1063/1.1749252
      3. 3
        Kolos, W.; Roothaan, C. C. J. Accurate Electronic Wave Functions for the H2 Molecule. Rev. Mod. Phys. 1960, 32, 219232,  DOI: 10.1103/RevModPhys.32.219
      4. 4
        Kołos, W.; Wolniewicz, L. The coupling between electronic and nuclear motion and the relativistic effects in the ground state of the H2 molecule. Acta Phys. Polonica 1961, 20, 129140
      5. 5
        Kołos, W.; Wolniewicz, L. Accurate Adiabatic Treatment of the Ground State of the Hydrogen Molecule. J. Chem. Phys. 1964, 41, 36633673,  DOI: 10.1063/1.1725796
      6. 6
        Kołos, W.; Wolniewicz, L. Accurate Computation of Vibronic Energies and of Some Expectation Values for H2, D2, and T2. J. Chem. Phys. 1964, 41, 36743678,  DOI: 10.1063/1.1725797
      7. 7
        Pachucki, K.; Komasa, J. Nonadiabatic corrections to rovibrational levels of H2. J. Chem. Phys. 2009, 130, 164113,  DOI: 10.1063/1.3114680
      8. 8
        Pachucki, K.; Komasa, J. Leading order nonadiabatic corrections to rovibrational levels of H2, D2, and T2. J. Chem. Phys. 2015, 143, 034111,  DOI: 10.1063/1.4927079
      9. 9
        Kołos, W.; Wolniewicz, L. Nonadiabatic Theory for Diatomic Molecules and Its Application to the Hydrogen Molecule. Rev. Mod. Phys. 1963, 35, 473,  DOI: 10.1103/RevModPhys.35.473
      10. 10
        Kinghorn, D. B.; Adamowicz, L. Improved Nonadiabatic Ground-State Energy Upper Bound for Dihydrogen. Phys. Rev. Lett. 1999, 83, 25412543,  DOI: 10.1103/PhysRevLett.83.2541
      11. 11
        Puchalski, M.; Komasa, J.; Czachorowski, P.; Pachucki, K. Nonadiabatic QED Correction to the Dissociation Energy of the Hydrogen Molecule. Phys. Rev. Lett. 2019, 122, 103003,  DOI: 10.1103/PhysRevLett.122.103003
      12. 12
        Pachucki, K.; Komasa, J. Nonadiabatic rotational states of the hydrogen molecule. Phys. Chem. Chem. Phys. 2018, 20, 247255,  DOI: 10.1039/C7CP06516G
      13. 13
        Wolniewicz, L. Relativistic energies of the ground state of the hydrogen molecule. J. Chem. Phys. 1993, 99, 1851,  DOI: 10.1063/1.465303
      14. 14
        Piszczatowski, K.; Łach, G.; Przybytek, M.; Komasa, J.; Pachucki, K.; Jeziorski, B. Theoretical Determination of the Dissociation Energy of Molecular Hydrogen. J. Chem. Theory Comput. 2009, 5, 3039,  DOI: 10.1021/ct900391p
      15. 15
        Komasa, J.; Piszczatowski, K.; Łach, G.; Przybytek, M.; Jeziorski, B.; Pachucki, K. Quantum Electrodynamics Effects in Rovibrational Spectra of Molecular Hydrogen. J. Chem. Theory Comput. 2011, 7, 31053115,  DOI: 10.1021/ct200438t
      16. 16
        Puchalski, M.; Komasa, J.; Czachorowski, P.; Pachucki, K. Complete α6 m Corrections to the Ground State of H2. Phys. Rev. Lett. 2016, 117, 263002,  DOI: 10.1103/PhysRevLett.117.263002
      17. 17
        Czachorowski, P.; Puchalski, M.; Komasa, J.; Pachucki, K. Nonadiabatic relativistic correction in H2, D2, and HD. Phys. Rev. A 2018, 98, 052506,  DOI: 10.1103/PhysRevA.98.052506
      18. 18
        Siłkowski, M.; Pachucki, K.; Komasa, J.; Puchalski, M. Leading-order QED effects in the ground electronic state of molecular hydrogen. Phys. Rev. A 2023, 107, 032807,  DOI: 10.1103/PhysRevA.107.032807
      19. 19
        Wang, L. M.; Yan, Z.-C. Relativistic corrections to the ground state of H2 calculated without using the Born-Oppenheimer approximation. Phys. Rev. A 2018, 97, 060501,  DOI: 10.1103/PhysRevA.97.060501
      20. 20
        Puchalski, M.; Moszynski, R.; Komasa, J. In Polish Quantum Chemistry from Kołos to Now. In Advances in Quantum Chemistry; Musiał, M., Grabowski, I., Eds.; Academic Press, 2023; Vol. 87, pp 136.
      21. 21
        Pachucki, K.; Komasa, J. Magnetic dipole transitions in the hydrogen molecule. Phys. Rev. A 2011, 83, 032501,  DOI: 10.1103/PhysRevA.83.032501
      22. 22
        Margolis, J.; Hunt, G. On the level of H2 quadrupole absorption in the Jovian atmosphere. Icarus 1973, 18, 593598,  DOI: 10.1016/0019-1035(73)90061-4
      23. 23
        Roueff, E.; Abgrall, H.; Czachorowski, P.; Pachucki, K.; Puchalski, M.; Komasa, J. The full infrared spectrum of molecular hydrogen. Astron. Astrophys. 2019, 630, A58,  DOI: 10.1051/0004-6361/201936249
      24. 24
        Tchernyshyov, K. A Detection of H2 in a High-velocity Cloud toward the Large Magellanic Cloud. Astrophys. J. 2022, 931, 78,  DOI: 10.3847/1538-4357/ac68e0
      25. 25
        Jacob, A. Small molecules, big impact: a tale of hydrides past, present, and future. Astrophys. Space Sci. 2023, 368, 76,  DOI: 10.1007/s10509-023-04229-8
      26. 26
        Tennyson, J. The 2024 release of the ExoMol database: Molecular line lists for exoplanet and other hot atmospheres. J. Quant. Spectrosc. Radiat. Transfer 2024, 326, 109083,  DOI: 10.1016/j.jqsrt.2024.109083
      27. 27
        Bruno, D. Uncertainty quantification of a collisional-radiative model for hydrogen plasmas. Phys. Scr. 2025, 100, 035611,  DOI: 10.1088/1402-4896/adb227
      28. 28
        Zúñiga, J.; Bastida, A.; Requena, A.; Cerezo, J. Ab Initio Partition Functions and Thermodynamic Quantities for the Molecular Hydrogen Isotopologues. J. Phys. Chem. A 2021, 125, 92269241,  DOI: 10.1021/acs.jpca.1c06468
      29. 29
        Leachman, J. W.; Wilhelmsen, Ø.; Matveev, K. I. Cool Fuel: The Science and Engineering of Cryogenic Hydrogen; Oxford University Press, 2025.
      30. 30
        Puchalski, M.; Komasa, J.; Pachucki, K. Hyperfine Structure of the First Rotational Level in H2, D2 and HD Molecules and the Deuteron Quadrupole Moment. Phys. Rev. Lett. 2020, 125, 253001,  DOI: 10.1103/PhysRevLett.125.253001
      31. 31
        Puchalski, M.; Komasa, J.; Spyszkiewicz, A.; Pachucki, K. Nuclear magnetic shielding in HD and HT. Phys. Rev. A 2022, 105, 042802,  DOI: 10.1103/PhysRevA.105.042802
      32. 32
        Kassi, S.; Campargue, A. Electric quadrupole and dipole transitions of the first overtone band of HD by CRDS between 1.45 and 1.33 μm. J. Mol. Spectrosc. 2011, 267, 36,  DOI: 10.1016/j.jms.2011.02.001
      33. 33
        Campargue, A.; Kassi, S.; Pachucki, K.; Komasa, J. The absorption spectrum of H2: CRDS measurements of the (2–0) band, review of the literature data and accurate ab initio line list up to 35000 cm–1. Phys. Chem. Chem. Phys. 2012, 14, 802,  DOI: 10.1039/C1CP22912E
      34. 34
        Kassi, S.; Campargue, A.; Pachucki, K.; Komasa, J. The absorption spectrum of D2: Ultrasensitive cavity ring down spectroscopy of the (2–0) band near 1.7 μm and accurate ab initio line list up to 24 000 cm–1. J. Chem. Phys. 2012, 136, 184309,  DOI: 10.1063/1.4707708
      35. 35
        Cheng, C.-F.; Sun, Y. R.; Pan, H.; Wang, J.; Liu, A.-W.; Campargue, A.; Hu, S.-M. Electric-quadrupole transition of H2 determined to 10–9 precision. Phys. Rev. A 2012, 85, 024501,  DOI: 10.1103/PhysRevA.85.024501
      36. 36
        Dickenson, G. D.; Niu, M. L.; Salumbides, E. J.; Komasa, J.; Eikema, K. S. E.; Pachucki, K.; Ubachs, W. Fundamental Vibration of Molecular Hydrogen. Phys. Rev. Lett. 2013, 110, 193601,  DOI: 10.1103/PhysRevLett.110.193601
      37. 37
        Niu, M.; Salumbides, E.; Dickenson, G.; Eikema, K.; Ubachs, W. Precision spectroscopy of the rovibrational splittings in H2, HD and D2. J. Mol. Spectrosc. 2014, 300, 4454,  DOI: 10.1016/j.jms.2014.03.011
      38. 38
        Tan, Y.; Wang, J.; Cheng, C.-F.; Zhao, X.-Q.; Liu, A.-W.; Hu, S.-M. Cavity ring-down spectroscopy of the electric quadrupole transitions of in the 784–852 nm region. J. Mol. Spectrosc. 2014, 300, 6064,  DOI: 10.1016/j.jms.2014.03.010
      39. 39
        Kassi, S.; Campargue, A. Electric quadrupole transitions and collision-induced absorption in the region of the first overtone band of H2 near 1.25μm. J. Mol. Spectrosc. 2014, 300, 5559,  DOI: 10.1016/j.jms.2014.03.022
      40. 40
        Mondelain, D.; Kassi, S.; Sala, T.; Romanini, D.; Gatti, D.; Campargue, A. Sub-MHz accuracy measurement of the S(2) 2–0 transition frequency of D2 by Comb-Assisted Cavity Ring Down spectroscopy. J. Mol. Spectrosc. 2016, 326, 58,  DOI: 10.1016/j.jms.2016.02.008
      41. 41
        Wcisło, P.; Thibault, F.; Zaborowski, M.; Wójtewicz, S.; Cygan, A.; Kowzan, G.; Masłowski, P.; Komasa, J.; Puchalski, M.; Pachucki, K.; Ciuryło, R.; Lisak, D. Accurate deuterium spectroscopy for fundamental studies. J. Quant. Spectrosc. Radiat. Transfer 2018, 213, 4151,  DOI: 10.1016/j.jqsrt.2018.04.011
      42. 42
        Trivikram, T. M.; Schlösser, M.; Ubachs, W.; Salumbides, E. J. Relativistic and QED Effects in the Fundamental Vibration of T2. Phys. Rev. Lett. 2018, 120, 163002,  DOI: 10.1103/PhysRevLett.120.163002
      43. 43
        Tao, L.-G.; Liu, A.-W.; Pachucki, K.; Komasa, J.; Sun, Y. R.; Wang, J.; Hu, S.-M. Toward a Determination of the Proton-Electron Mass Ratio from the Lamb-Dip Measurement of HD. Phys. Rev. Lett. 2018, 120, 153001,  DOI: 10.1103/PhysRevLett.120.153001
      44. 44
        Fasci, E.; Castrillo, A.; Dinesan, H.; Gravina, S.; Moretti, L.; Gianfrani, L. Precision spectroscopy of HD at 1.38μm. Phys. Rev. A 2018, 98, 022516,  DOI: 10.1103/PhysRevA.98.022516
      45. 45
        Martinez, R. Z.; Bermejo, D.; Wcisło, P.; Thibault, F. Accurate wavenumber measurements for the S0(0), S0(1), and S0(2) pure rotational Raman lines of D2. J. Raman Spectrosc. 2019, 50, 127129,  DOI: 10.1002/jrs.5499
      46. 46
        Hölsch, N.; Beyer, M.; Salumbides, E. J.; Eikema, K. S. E.; Ubachs, W.; Jungen, C.; Merkt, F. Benchmarking Theory with an Improved Measurement of the Ionization and Dissociation Energies of H2. Phys. Rev. Lett. 2019, 122, 103002,  DOI: 10.1103/PhysRevLett.122.103002
      47. 47
        Beyer, M.; Hölsch, N.; Hussels, J.; Cheng, C.-F.; Salumbides, E. J.; Eikema, K. S. E.; Ubachs, W.; Jungen, C.; Merkt, F. Determination of the Interval between the Ground States of Para- and Ortho-H2. Phys. Rev. Lett. 2019, 123, 163002,  DOI: 10.1103/PhysRevLett.123.163002
      48. 48
        Lai, K.-F.; Czachorowski, P.; Schlösser, M.; Puchalski, M.; Komasa, J.; Pachucki, K.; Ubachs, W.; Salumbides, E. J. Precision tests of nonadiabatic perturbation theory with measurements on the DT molecule. Phys. Rev. Res. 2019, 1, 033124,  DOI: 10.1103/PhysRevResearch.1.033124
      49. 49
        Lai, K.-F.; Hermann, V.; Trivikram, M.; Diouf, M.; Schlösser, M.; Ubachs, W.; Salumbides, E. Precision measurement of the fundamental vibrational frequencies of tritium-bearing hydrogen molecules: T2, DT, HT. Phys. Chem. Chem. Phys. 2020, 22, 89738987,  DOI: 10.1039/D0CP00596G
      50. 50
        Mondelain, D.; Kassi, S.; Campargue, A. Transition frequencies in the (2–0) band of D2 with MHz accuracy. J. Quant. Spectrosc. Radiat. Transfer 2020, 253, 107020,  DOI: 10.1016/j.jqsrt.2020.107020
      51. 51
        Jóźwiak, H.; Cybulski, H.; Wcisło, P. Hyperfine structure of quadrupole rovibrational transitions in tritium-bearing hydrogen isotopologues. J. Quant. Spectrosc. Radiat. Transfer 2020, 256, 107255,  DOI: 10.1016/j.jqsrt.2020.107255
      52. 52
        Kassi, S.; Fleurbaey, H.; Campargue, A. First detection and absolute transition frequencies in the (3–0) band of D2. J. Chem. Phys. 2024, 160, 094306,  DOI: 10.1063/5.0196903
      53. 53
        Salumbides, E. J.; Dickenson, G. D.; Ivanov, T. I.; Ubachs, W. QED Effects in Molecules: Test on Rotational Quantum States of H2. Phys. Rev. Lett. 2011, 107, 043005,  DOI: 10.1103/PhysRevLett.107.043005
      54. 54
        Salumbides, E. J.; Koelemeij, J. C. J.; Komasa, J.; Pachucki, K.; Eikema, K. S. E.; Ubachs, W. Bounds on fifth forces from precision measurements on molecules. Phys. Rev. D 2013, 87, 112008,  DOI: 10.1103/PhysRevD.87.112008
      55. 55
        Salumbides, E. J.; Schellekens, A. N.; Gato-Rivera, B.; Ubachs, W. Constraints on extra dimensions from precision molecular spectroscopy. New J. Phys. 2015, 17, 033015,  DOI: 10.1088/1367-2630/17/3/033015
      56. 56
        Ubachs, W.; Koelemeij, J.; Eikema, K.; Salumbides, E. Physics beyond the Standard Model from hydrogen spectroscopy. J. Mol. Spectrosc. 2016, 320, 112,  DOI: 10.1016/j.jms.2015.12.003
      57. 57
        Drouin, B. J.; Yu, S.; Pearson, J. C.; Gupta, H. Terahertz spectroscopy for space applications: 2.5–2.7 THz spectra of HD, H2O and NH3. J. Mol. Spectrosc. 2011, 1006, 212,  DOI: 10.1016/j.molstruc.2011.05.062
      58. 58
        Cozijn, F. M. J.; Dupré, P.; Salumbides, E. J.; Eikema, K. S. E.; Ubachs, W. Sub-Doppler Frequency Metrology in HD for Tests of Fundamental Physics. Phys. Rev. Lett. 2018, 120, 153002,  DOI: 10.1103/PhysRevLett.120.153002
      59. 59
        Diouf, M. L.; Cozijn, F. M. J.; Darquié, B.; Salumbides, E. J.; Ubachs, W. Lamb-dips and Lamb-peaks in the saturation spectrum of HD. Opt. Lett. 2019, 44, 47334736,  DOI: 10.1364/OL.44.004733
      60. 60
        Diouf, M. L.; Cozijn, F. M. J.; Lai, K.-F.; Salumbides, E. J.; Ubachs, W. Lamb-peak spectrum of the HD (2–0) P(1) line. Phys. Rev. Res. 2020, 2, 023209,  DOI: 10.1103/PhysRevResearch.2.023209
      61. 61
        Hua, T.-P.; Sun, Y. R.; Hu, S.-M. Dispersion-like lineshape observed in cavity-enhanced saturation spectroscopy of HD at 1.4 μm. Opt. Lett. 2020, 45, 48634866,  DOI: 10.1364/OL.401879
      62. 62
        Wójtewicz, S.; Gotti, R.; Gatti, D.; Lamperti, M.; Laporta, P.; Jóźwiak, H.; Thibault, F.; Wcisło, P.; Marangoni, M. Accurate deuterium spectroscopy and comparison with ab initio calculations. Phys. Rev. A 2020, 101, 052504,  DOI: 10.1103/PhysRevA.101.052504
      63. 63
        Fast, A.; Meek, S. A. Sub-ppb Measurement of a Fundamental Band Rovibrational Transition in HD. Phys. Rev. Lett. 2020, 125, 023001,  DOI: 10.1103/PhysRevLett.125.023001
      64. 64
        Castrillo, A.; Fasci, E.; Gianfrani, L. Doppler-limited precision spectroscopy of HD at 1.4 μm: An improved determination of the R(1) center frequency. Phys. Rev. A 2021, 103, 022828,  DOI: 10.1103/PhysRevA.103.022828
      65. 65
        Fast, A.; Meek, S. A. Precise measurement of the D2 S1(0) vibrational transition frequency. Mol. Phys. 2022, 120, e1999520  DOI: 10.1080/00268976.2021.1999520
      66. 66
        Kassi, S.; Lauzin, C.; Chaillot, J.; Campargue, A. The (2–0) R(0) and R(1) transition frequencies of HD determined to a 10–10 relative accuracy by Doppler spectroscopy at 80 K. Phys. Chem. Chem. Phys. 2022, 24, 2316423172,  DOI: 10.1039/D2CP02151J
      67. 67
        Cozijn, F. M. J.; Diouf, M. L.; Hermann, V.; Salumbides, E. J.; Schlösser, M.; Ubachs, W. Rotational level spacings in HD from vibrational saturation spectroscopy. Phys. Rev. A 2022, 105, 062823,  DOI: 10.1103/PhysRevA.105.062823
      68. 68
        Cozijn, F. M. J.; Diouf, M. L.; Ubachs, W. Saturation spectroscopy of R(0), R(2) and P(2) lines in the (2–0) band of HD. Eur. Phys. J. D 2022, 76, 220,  DOI: 10.1140/epjd/s10053-022-00552-x
      69. 69
        Fleurbaey, H.; Koroleva, A. O.; Kassi, S.; Campargue, A. The high-accuracy spectroscopy of H2 rovibrational transitions in the (2–0) band near 1.2 μm. Phys. Chem. Chem. Phys. 2023, 25, 1474914756,  DOI: 10.1039/D3CP01136D
      70. 70
        Lamperti, M.; Rutkowski, L.; Ronchetti, D.; Gatti, D.; Gotti, R.; Cerullo, G.; Thibault, F.; Jóźwiak, H.; Wójtewicz, S.; Masłowski, P.; Wcisło, P.; Polli, D.; Marangoni, M. Stimulated Raman scattering metrology of molecular hydrogen. Commun. Phys. 2023, 6, 67,  DOI: 10.1038/s42005-023-01187-z
      71. 71
        Cozijn, F. M. J.; Diouf, M. L.; Ubachs, W. Lamb Dip of a Quadrupole Transition in H2. Phys. Rev. Lett. 2023, 131, 073001,  DOI: 10.1103/PhysRevLett.131.073001
      72. 72
        Cozijn, F. M. J.; Diouf, M. L.; Ubachs, W.; Hermann, V.; Schlösser, M. Precision Measurement of Vibrational Quanta in Tritium Hydride. Phys. Rev. Lett. 2024, 132, 113002,  DOI: 10.1103/PhysRevLett.132.113002
      73. 73
        Stankiewicz, K.; Makowski, M.; Słowinski, M.; Sołtys, K. L.; Bednarski, B.; Jóźwiak, H.; Stolarczyk, N.; Narożnik, M.; Kierski, D.; Wójtewicz, S.; Cygan, A.; Kowzan, G.; Masłowski, P.; Piwinski, M.; Lisak, D.; Wcisło, P. Cavity-enhanced spectroscopy in the deep cryogenic regime – new hydrogen technologies for quantum sensing. arXiv 2025, arXiv:2502.12703,  DOI: 10.48550/arXiv.2502.12703
      74. 74
        Pachucki, K.; Komasa, J. Schrödinger equation solved for the hydrogen molecule with unprecedented accuracy. J. Chem. Phys. 2016, 144, 164306,  DOI: 10.1063/1.4948309
      75. 75
        Pachucki, K.; Komasa, J. Integrals for relativistic nonadiabatic energies of H2 in an exponential basis. Phys. Rev. A 2024, 109, 032822,  DOI: 10.1103/PhysRevA.109.032822
      76. 76
        Pachucki, K.; Komasa, J. Nonrelativistic energy levels of HD. Phys. Chem. Chem. Phys. 2018, 20, 2629726302,  DOI: 10.1039/C8CP05493B
      77. 77
        Pachucki, K.; Komasa, J. Nonrelativistic energy levels of D2. Phys. Chem. Chem. Phys. 2019, 21, 1027210276,  DOI: 10.1039/C9CP01308C
      78. 78
        Pachucki, K.; Komasa, J. Nonrelativistic energy of tritium-containing hydrogen molecule isotopologues. Mol. Phys. 2022, 120, e2040627  DOI: 10.1080/00268976.2022.2040627
      79. 79
        Hylleraas, E. Neue Berechnung der Energie des Heliums im Grundzustande, sowie des tiefsten Terms von Ortho-Helium. Z. Physik 1929, 54, 347366,  DOI: 10.1007/BF01375457
      80. 80
        Mohr, P. J.; Newell, D. B.; Taylor, B. N.; Tiesinga, E. CODATA recommended values of the fundamental physical constants: 2022. Rev. Mod. Phys. 2025, 97, 025002,  DOI: 10.1103/RevModPhys.97.025002
      81. 81
        Pachucki, K.; Komasa, J. Relativistic Correction from the Four-Body Nonadiabatic Exponential Wave Function. J. Chem. Theory Comput. 2024, 20, 86448651,  DOI: 10.1021/acs.jctc.4c00861
      82. 82
        Code, R. F.; Ramsey, N. F. Molecular-Beam Magnetic Resonance Studies of HD and D2. Phys. Rev. A 1971, 4, 19451959,  DOI: 10.1103/PhysRevA.4.1945
      83. 83
        Komasa, J. Electron-nucleus Araki-Sucher QED correction to rovibrational levels of H2. arXiv 2025, arXiv:2511.07053v2.
      84. 84
        Drake, G. W. F.; Swainson, R. A. Bethe logarithms for hydrogen up to n = 20, and approximations for two-electron atoms. Phys. Rev. A 1990, 41, 12431246,  DOI: 10.1103/PhysRevA.41.1243
      85. 85
        Diouf, M. L.; Cozijn, F. M. J.; Ubachs, W. Hyperfine structure in a vibrational quadrupole transition of ortho-H2. Mol. Phys. 2024, 122, e2304101  DOI: 10.1080/00268976.2024.2304101
      86. 86
        H2Spectre ver. 7.4 Fortran source code. 2022, https://www.fuw.edu.pl/~krp/, https://qcg.home.amu.edu.pl/H2Spectre.html (accessed Oct 10, 2025).
      87. 87
        Parthey, C. G.; Matveev, A.; Alnis, J.; Bernhardt, B.; Beyer, A.; Holzwarth, R.; Maistrou, A.; Pohl, R.; Predehl, K.; Udem, T.; Wilken, T.; Kolachevsky, N.; Abgrall, M.; Rovera, D.; Salomon, C.; Laurent, P.; Hänsch, T. W. Improved Measurement of the Hydrogen 1S–2S Transition Frequency. Phys. Rev. Lett. 2011, 107, 203001,  DOI: 10.1103/PhysRevLett.107.203001
      88. 88
        Pachucki, K.; Jentschura, U. D. Two-Loop Bethe-Logarithm Correction in Hydrogenlike Atoms. Phys. Rev. Lett. 2003, 91, 113005,  DOI: 10.1103/PhysRevLett.91.113005
      89. 89
        Clausen, G.; Scheidegger, S.; Agner, J. A.; Schmutz, H.; Merkt, F. Imaging-Assisted Single-Photon Doppler-Free Laser Spectroscopy and the Ionization Energy of Metastable Triplet Helium. Phys. Rev. Lett. 2023, 131, 103001,  DOI: 10.1103/PhysRevLett.131.103001
      90. 90
        Alighanbari, S.; Schenkel, M.; Korobov, V.; Schiller, S. High-accuracy laser spectroscopy of and the proton–electron mass ratio. Nature 2025, 644, 6975,  DOI: 10.1038/s41586-025-09306-2
    • Supporting Information

      Supporting Information


      The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.5c01702.

      • Convergence of relativistic expectation values for individual rovibrational levels of H2 (PDF)


      Terms & Conditions

      Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.