Rovibrational computations for the He$_2$ a $^3Σ_\mathrm{u}^+$ state including non-adiabatic, relativistic, and QED corrections

This paper presents a highly accurate theoretical study of the He$_2$ a $^3\Sigma_\mathrm{u}^+$ state, incorporating relativistic and QED corrections to compute potential energy curves and rovibrational levels that show remarkable agreement with high-resolution spectroscopic data.

The Gist

Imagine you are trying to build the most perfect, microscopic model of a helium molecule ($He_2$) that has ever existed. But this isn't just any helium molecule; it's a "triplet" state, meaning its electrons are spinning in a specific, slightly unstable way that makes it a rare and fascinating ghost in the world of chemistry.

This paper is the story of a team of scientists who built a digital twin of this molecule so accurate that it matches real-world experiments down to the billionth of a unit. Here is how they did it, explained without the heavy math.

1. The Challenge: The "Fuzzy" Map

To understand how a molecule moves, vibrates, and spins, you first need a map of its energy. In physics, this is called a Potential Energy Curve (PEC). Think of this map like a topographical map of a mountain valley.

• The Valley: Where the molecule likes to sit (the stable spot).
• The Walls: How hard it is to pull the atoms apart.

For decades, scientists had a "rough sketch" of this map for helium. It was good enough to see the mountain, but not good enough to navigate the tiny pebbles on the path. Previous maps had errors so big that if you tried to predict the molecule's behavior, your prediction would be off by a mile.

2. The Solution: The "Super-Microscope"

The authors used a technique called Explicitly Correlated Gaussians.

• The Analogy: Imagine trying to describe a complex shape made of clay. Old methods used a few large, blocky Lego bricks to approximate the shape. No matter how many bricks you added, the shape always looked a bit jagged.
• The New Method: This team used "liquid clay" that could flow and mold perfectly around every curve. They didn't just guess the shape; they calculated the exact interaction between every single electron and every single nucleus simultaneously. This allowed them to draw the energy map with a precision of 1 part in a million (sub-ppm).

3. Adding the "Invisible" Details

Once they had the perfect map, they realized the map wasn't enough. To match reality, they had to add corrections for things that usually don't matter, but for a precision experiment, they are everything. They added three layers of "fine-tuning":

• Relativistic Corrections (The "Fast" Effect): Electrons move incredibly fast. According to Einstein, fast-moving things gain mass and behave differently. The team added these "speed bumps" to their model.
• QED Corrections (The "Quantum Foam" Effect): In the quantum world, empty space isn't empty. It's a bubbling foam of virtual particles popping in and out of existence. These particles nudge the electrons slightly. The team calculated these tiny nudges, which are like the wind blowing on a feather.
• Non-Adiabatic Corrections (The "Dancing" Effect): Usually, we assume the heavy nuclei (the core of the atom) stand still while the light electrons dance around them. But in reality, the nuclei wobble a little bit too, and the electrons react to that wobble instantly. The team modeled this "dance" perfectly.

4. The Result: A Perfect Match

After building this incredibly complex model, they asked: "Does our digital twin behave like the real molecule?"

They compared their calculations to the most precise laser experiments ever done on helium.

• The Outcome: It was a home run. Their predictions for how the molecule vibrates, rotates, and splits its energy levels (fine structure) matched the experimental data almost perfectly.
• The Scale: The energy levels they calculated span a massive range. They got the "big picture" (ionization energy) right, the "medium picture" (vibrations) right, and even the "tiny picture" (fine-structure splittings) right.
  • Analogy: It's like predicting the height of a skyscraper, the length of a hallway inside it, and the thickness of a single brick, all with the same ruler, and hitting the mark every time.

5. Why Does This Matter?

You might ask, "Who cares about a helium molecule?"

• Testing the Laws of Physics: Helium is the simplest molecule with more than one electron. If our theories of physics (Quantum Mechanics and Relativity) can't explain helium perfectly, they are broken. This paper proves our current theories are rock solid.
• The "Ruler" of Science: By knowing the exact energy of this molecule, scientists can use it as a standard ruler to measure other things, potentially helping to refine the values of fundamental constants of the universe.
• Future Tech: Understanding these "triplet" states helps scientists figure out how to control and cool these molecules with lasers, which could lead to new types of quantum computers or ultra-precise sensors.

Summary

In short, this paper is a masterclass in precision. The authors took a simple helium molecule, built a digital model so detailed it accounted for the jitter of virtual particles and the wobble of nuclei, and proved that our understanding of the quantum world is accurate enough to predict reality with microscopic perfection. They didn't just see the forest; they counted the leaves, measured the veins, and predicted the color of the next leaf to fall.

Technical Summary

1. Problem Statement

The helium dimer (He$_2$) in its lowest triplet excited state ($a$ $^3\Sigma^+_u$) is a critical benchmark system for testing fundamental physical theories, including Quantum Electrodynamics (QED) and the precision of physical constants. While experimental spectroscopy for this few-electron molecule has reached extreme precision (uncertainties down to sub-kHz or $10^{-9}$ cm$^{-1}$), theoretical predictions have historically lagged behind.

Previous theoretical attempts suffered from:

• Insufficient Basis Sets: Standard uncorrelated atomic basis sets failed to reduce non-relativistic electronic energy errors below $\sim 1$ mE$_h$, which is insufficient for high-precision spectroscopy.
• Incomplete Physics: Many studies neglected or approximated non-adiabatic, relativistic, and QED effects, which are crucial for matching experimental fine-structure splittings and ionization energies.
• Lack of Explicit Correlation: There was a lack of fully variational, explicitly correlated potential energy curve (PEC) computations for this specific state.

The goal of this work is to compute the rovibrational energy levels of He$_2$ ($a$ $^3\Sigma^+_u$) with an accuracy of a fraction of 1 ppm ($10^{-6}$), explicitly accounting for all relevant physical corrections to achieve agreement with high-resolution experimental data.

2. Methodology

The authors employed a multi-step, high-precision computational framework:

A. Electronic Structure and Potential Energy Curve (PEC)

• Wavefunction: The electronic wavefunction was parametrized using a truncated expansion of floating explicitly correlated Gaussians (fECGs). This approach explicitly includes inter-electronic distances ($r_{ij}$) in the basis functions, allowing for rapid convergence toward the complete basis set limit.
• Variational Optimization: The PEC was generated over the internuclear distance range $\rho/a_0 \in [1, 100]$. The optimization utilized a stochastic variational approach with the Powell method.
• Basis Size: Calculations used up to 2,500 basis functions. The convergence error for the Born-Oppenheimer (BO) energy at the equilibrium distance ($\rho \approx 2$ a$_0$) was estimated at $\sim 1$ $\mu$E$_h$ (micro-Hartree).

B. Corrections to the Electronic Energy
The corrected PEC ($W$) was constructed by adding several perturbative corrections to the BO energy ($U$):

1. Diagonal Born-Oppenheimer Correction (DBOC): Included electronic, vibrational, and angular momentum parts to account for finite nuclear mass effects within the adiabatic approximation.
2. Non-Adiabatic Mass Corrections: Vibrational ($\delta m_{vib}$) and rotational ($\delta m_{rot}$) corrections were calculated to account for coupling with distant electronic states, effectively modifying the reduced mass in the nuclear kinetic energy term.
3. Relativistic Corrections: Calculated as the expectation value of the spin-independent Breit-Pauli Hamiltonian ($\hat{H}_{BP,sn}$), including mass-velocity, Darwin, and orbit-orbit terms.
4. QED Corrections:
   • Leading Order ($U_{lQED}$): Included self-energy, vertex corrections, and vacuum polarization. The Bethe logarithm ($\ln k_0$) was approximated using an ion-core approximation (using values from He$_2^{3+}$) to avoid prohibitively expensive direct calculations, justified by the weak dependence of $\ln k_0$ on electron number.
   • Higher Order ($U_{hQED}$): Estimated using hydrogenic theory.
   • Finite Nuclear Size: Included but found to be negligible for this state.
5. Spin-Dependent Corrections: The magnetic dipole-dipole interaction ($\hat{H}_{SS,dp}$) was treated separately to calculate fine-structure splittings. The electron $g$-factor was corrected for QED effects (anomalous magnetic moment).

C. Rovibrational Solution

• The nuclear Schrödinger equation was solved using the corrected PEC and non-adiabatic mass corrections.
• Method: A Discrete Variable Representation (DVR) using associated Laguerre polynomials was employed on a grid of 800 points.
• Symmetry: Only states symmetric under the exchange of the two $^4$He nuclei (bosons) were retained, restricting rotational quantum numbers to odd values.
• Fine Structure: The spin-rotation coupling was neglected in the direct computation (consistent with the effective Hamiltonian fitting used for comparison), focusing on the magnetic spin-spin splitting.

3. Key Contributions

• First High-Precision PEC: This work provides the first potential energy curve for the He$_2$ $a$ $^3\Sigma^+_u$ state computed with explicitly correlated Gaussians, achieving sub-ppm relative precision.
• Comprehensive Physics: It is the first study to simultaneously include non-adiabatic, relativistic, and QED corrections (including the anomalous magnetic moment) for the full rovibrational spectrum of this state.
• Ion-Core Approximation Validation: The authors validated the use of the ion-core approximation for the Bethe logarithm in molecular systems, demonstrating its sufficiency for high-precision rovibrational intervals.
• Regularization Techniques: The use of "numerical Drachmanization" (numDr) allowed for the stable computation of expectation values for highly singular operators (e.g., $\delta$-functions) which typically converge slowly in Gaussian bases.

4. Results

The computed results show remarkable agreement with experimental data across various energy scales:

• Ionization Energy:
  • Computed: $34,301.07(25)$ cm$^{-1}$.
  • Experimental: $34,301.20700(4)$ cm$^{-1}$.
  • The deviation ($\sim 0.14$ cm$^{-1}$) is dominated by the remaining convergence error in the BO energy, not by missing physics.
• Vibrational Intervals:
  • Agreement is excellent for low vibrational levels ($v=1, 2$), with deviations of $\sim 0.03-0.04$ cm$^{-1}$.
  • Higher levels show slightly larger discrepancies, attributed to non-uniform PEC errors or experimental uncertainties in older data.
• Rotational Intervals:
  • For the $v=0$ band, the agreement with high-precision data is exceptional, with deviations as small as $0.00004$ cm$^{-1}$ ($\sim 1$ MHz).
  • The theoretical uncertainty is estimated at $\sim 10^{-3}$ cm$^{-1}$, largely due to PEC convergence errors canceling out in interval differences.
• Fine-Structure Splittings:
  • The computed splittings agree with experimental values within 250–500 kHz ($\sim 0.000008$ cm$^{-1}$).
  • Crucial Finding: Including the QED correction to the electron $g$-factor (anomalous magnetic moment) was essential. Omitting this would introduce a uniform error of $\sim 3$ MHz.
  • The remaining small discrepancy is attributed to neglected non-adiabatic-relativistic coupling and higher-order relativistic terms.

5. Significance

• Benchmark for Theory: This work establishes a new benchmark for ab initio molecular quantum chemistry, demonstrating that explicitly correlated methods can achieve spectroscopic accuracy comparable to the most precise experiments.
• Validation of QED in Molecules: The successful reproduction of fine-structure splittings confirms the necessity and accuracy of including QED corrections (specifically the anomalous magnetic moment) in molecular calculations.
• Pathway to New Physics: The high accuracy of these computations allows for the extraction of fundamental constants (like the proton-to-electron mass ratio or nuclear charge radii) from molecular spectra with unprecedented precision.
• Future Applications: The methodology paves the way for studying more complex excited states (e.g., $b$ $^3\Pi_g$, $c$ $^3\Sigma^+_g$) and potentially aids in the development of laser cooling techniques for helium dimers, as suggested by the authors.

In summary, this paper represents a milestone in theoretical molecular physics, bridging the gap between high-precision theory and experiment for the helium dimer by rigorously treating all relevant physical interactions from the Born-Oppenheimer approximation up to QED corrections.