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

This paper presents highly accurate rovibrational computations for the ground electronic state of the helium dimer cation (He$_2^+$) by calculating comprehensive potential energy corrections—including non-adiabatic, relativistic, and QED effects—over a broad configuration range, thereby determining all bound states with an estimated accuracy of 0.005 cm$^{-1}$.

The Gist

Imagine you are trying to build the most perfect, miniature model of a tiny molecule called Helium-2 Plus ($He_2^+$). This isn't just any molecule; it's a "three-electron" system (two helium nuclei sharing three electrons), and it's the smallest molecule with a real chemical bond.

For decades, scientists have been trying to predict exactly how this molecule vibrates and rotates, hoping their math matches reality perfectly. Think of it like trying to predict the exact pitch of a tiny guitar string. If your prediction is off by even a tiny fraction, it means your understanding of the universe's laws is slightly incomplete.

This paper is the story of a team of scientists who decided to build the ultimate, ultra-precise model of this molecule. Here is how they did it, explained simply:

1. The Blueprint: The "Energy Map"

First, the team needed a map of the molecule's energy. Imagine the two helium nuclei are like two magnets. If you pull them apart, they want to snap back together; if you push them too close, they repel.

• The Old Map: Previous models were like a low-resolution GPS. They knew the general shape of the terrain but missed the tiny bumps and dips.
• The New Map: The authors created a high-definition, 4K-resolution map. They used a super-complex mathematical technique (involving thousands of "Gaussian clouds" to represent where the electrons hang out) to calculate the energy at every single point between the nuclei. This is the Potential Energy Curve.

2. The Fine-Tuning: Adding the "Invisible Forces"

In the old days, scientists would stop at the main map. But in the quantum world, there are tiny, invisible forces that nudge the molecule just a little bit. To get the perfect prediction, you have to add these "corrections" to your map. The authors added four layers of fine-tuning:

• The "Wiggle" Correction (Non-Adiabatic): Usually, we pretend the heavy nuclei are stationary while the light electrons zip around. But in reality, the nuclei wiggle a bit too, and the electrons react instantly. It's like a dance where the partners are constantly adjusting their steps. The authors calculated exactly how this dance changes the energy.
• The "Heavy" Correction (Relativistic): Electrons move incredibly fast. When things move that fast, Einstein's theory of relativity kicks in, making them act slightly heavier. The authors added a "speed tax" to their calculation to account for this.
• The "Quantum Jitter" Correction (QED): This is the weirdest one. According to Quantum Electrodynamics (QED), empty space isn't empty; it's bubbling with virtual particles popping in and out of existence. These particles "jitter" against the electrons, changing their energy slightly. It's like the molecule is sitting in a gentle, invisible rainstorm that pushes it around. The authors calculated the force of this rain.
• The "Fuzzy Nucleus" Correction: We often think of the nucleus as a tiny, hard marble. But it's actually a fuzzy cloud of protons. The authors adjusted their math to account for this fuzziness.

3. The Result: A Perfect Match

Once they built this super-detailed map with all the corrections, they solved the equations to find out exactly how the molecule vibrates and rotates.

• The Achievement: They predicted the energy levels of every possible vibration and rotation of the molecule with an accuracy of 0.005 cm⁻¹.
• The Comparison: To put that in perspective, if the molecule's energy were the height of a skyscraper, their prediction would be accurate to within the width of a human hair.
• The Validation: When they compared their numbers to the most recent, ultra-precise experiments, the numbers matched almost perfectly. In fact, their new calculations were so good that they helped explain tiny discrepancies in older data.

Why Does This Matter?

You might ask, "Why bother calculating a tiny helium molecule so precisely?"

1. Testing the Laws of Physics: This molecule is a "laboratory" for testing the Standard Model of physics. If our calculations (which use all our current laws of physics) don't match the experiment, it means we are missing a piece of the puzzle. Since they matched, it confirms our current understanding of quantum mechanics is rock solid.
2. Refining Constants: These calculations help scientists pin down the exact values of fundamental constants (like the mass of an electron or the strength of the electromagnetic force) with greater precision.
3. The "Spin" Mystery: The authors noticed a tiny remaining difference in the most excited states of the molecule. They suspect this is due to a magnetic interaction between the electron's spin and the molecule's rotation—a "magnetic handshake" they haven't fully calculated yet. This gives them a new target for future research.

The Bottom Line

Think of this paper as the team that built the ultimate simulation of a helium molecule. They didn't just guess; they accounted for every tiny wobble, speed effect, and quantum jitter. The result is a theoretical model so accurate that it rivals the best experiments in the world, proving that with enough math and computing power, we can understand the universe down to its tiniest, most fundamental details.

Technical Summary

1. Problem Statement

The paper addresses the need for ultra-high-precision theoretical descriptions of the helium dimer cation (He$_2^+$), specifically its ground electronic state ($X$ $^2\Sigma_u^+$). While He$_2^+$ is a simple three-electron system, its rovibrational spectrum is rich and suitable for precision spectroscopy.

• Context: Experimental techniques have advanced to measure rovibrational intervals with uncertainties as low as $\sim 10^{-4}$ cm$^{-1}$.
• Challenge: To match this experimental precision, theoretical models must account for subtle physical effects often neglected in standard quantum chemistry: non-adiabatic couplings, relativistic effects, Quantum Electrodynamics (QED) corrections, and finite nuclear size effects.
• Goal: To compute the potential energy curve (PEC) and all rovibrational bound states of He$_2^+$ with an estimated accuracy of 0.005 cm$^{-1}$, improving upon previous theoretical works (specifically Refs. [18–20]) by refining the error control and expanding the internuclear distance range.

2. Methodology

The authors employ a rigorous variational approach based on the Born-Oppenheimer (BO) framework, corrected via perturbation theory.

A. Electronic Structure Calculation

• Wave Function: The non-relativistic electronic wave function is expanded as a linear combination of Explicitly Correlated Gaussian (ECG) functions with floating centers (fECGs).
  • Basis size: Up to 2,250 basis functions were used.
  • Optimization: Parameters (non-linear shift vectors and matrix exponents) were optimized using the Stochastic Variational Method (SVM) followed by the Powell minimization approach.
  • Symmetry: The wave function is projected onto the $\Sigma_u^+$ symmetry of the $D_{\infty h}$ point group.
• Software: Computations were performed using the in-house code QUANTEN.

B. Correction Terms

The total energy ($W$) is constructed by adding various corrections to the BO electronic energy ($U$):
$$W = U + U_{\text{DBOC}} + U_{\text{f.nuc.}} + U_{\text{rel+QED}}$$

1. Diagonal Born-Oppenheimer Correction (DBOC): Computed as expectation values of electronic kinetic energy and angular momentum operators.
2. Non-Adiabatic Mass Corrections: Effective rotational ($\mu_{\text{rot}}$) and vibrational ($\mu_{\text{vib}}$) reduced masses are calculated to account for coupling with excited electronic states.
3. Relativistic Corrections ($U_{\text{rel}}$): Calculated as the expectation value of the spin-independent Breit-Pauli Hamiltonian (mass-velocity, Darwin, orbit-orbit terms).
4. QED Corrections:
   • Leading-order ($\alpha^3 E_h$): Includes the Bethe logarithm ($\ln k_0$) and the Araki-Sucher term ($P(1/r^3)$).
   • Higher-order ($\alpha^4 E_h$): Approximated based on hydrogen-like extensions.
   • Finite Nuclear Size: Correction for the non-point-like nature of the $^4$He nucleus.
5. Regularization Techniques: To handle the slow convergence of singular operators (like $\delta$-functions and $1/r^3$) in Gaussian basis sets, the authors utilized two regularization schemes:
   • Integral Transformation (IT): Splits integrals into low/high-range regions.
   • Numerical Drachmanization (numDr): Rewrites singular operators into less singular forms using numerical integration. The numDr approach was selected for the full PEC generation due to its robustness.

C. Rovibrational Solution

• The corrected PEC and coordinate-dependent reduced masses were used to solve the radial Schrödinger equation for rotational quantum numbers $N$.
• Method: Discrete Variable Representation (DVR) using Laguerre polynomials scaled to the interval $[0.3, 100]$ bohr.
• Range: Calculations covered internuclear distances from 0.3 to 100.5 bohr, a significant expansion over previous works.

3. Key Contributions

1. Improved Error Control: The study achieves a convergence of the electronic energy within 10 nEh of the exact value and controls numerical errors in corrections to ensure a total rovibrational uncertainty of 0.005 cm$^{-1}$.
2. Advanced Regularization: The successful application and comparison of Numerical Drachmanization for singular QED and relativistic operators in molecular systems with floating ECGs.
3. Comprehensive Dataset: Generation of a complete PEC and correction curves over a broad range ($0.3 - 100$ bohr), including high-lying vibrational states near the dissociation limit.
4. Bethe Logarithm Approximation: Validated the approximation of using the He$_2^{3+}$ (one-electron) Bethe logarithm for the three-electron He$_2^+$ system, confirmed by explicit calculation at specific points.

4. Results

• Potential Energy Curve: The new PEC shows superior convergence compared to previous studies, particularly near the potential minimum.
• Rovibrational Energies:
  • Low-lying states ($v=0, 1$): The computed intervals show excellent agreement with experimental data (Refs. [4, 7, 10]), with deviations often within the experimental uncertainty ($\sim 0.001$ cm$^{-1}$).
  • High-lying states ($v=22, 23$): While previous computations were already close to experiment, this work refines the values. The remaining small deviations ($\sim 0.002$ cm$^{-1}$) are attributed to the electron spin-rotation coupling, which was not included in the current calculation but is present in the experimental data.
• Comparison with Experiment:
  • Rotational intervals ($v=0$) match experiment with differences $< 0.003$ cm$^{-1}$.
  • Vibrational fundamentals ($v=1$) show significant improvement over Ref. [20], reducing discrepancies from $\sim 0.003$ cm$^{-1}$ to $< 0.001$ cm$^{-1}$.
  • The total improvement allows the theoretical results to align with the comprehensive experimental dataset (0.1 cm$^{-1}$ uncertainty) published recently (Ref. [10]).

5. Significance

• Precision Physics Benchmark: This work provides one of the most accurate theoretical descriptions of a few-body molecular system to date. It serves as a critical benchmark for testing fundamental physical constants and the validity of QED in molecular environments.
• Methodological Advancement: The successful integration of numerical Drachmanization with variational fECG methods sets a new standard for computing relativistic and QED corrections in molecules, overcoming previous convergence hurdles.
• Future Directions: The authors identify the non-adiabatic relativistic coupling and electron spin-rotation coupling as the next frontiers. The current results highlight that these effects are now the dominant sources of discrepancy between theory and experiment for He$_2^+$, guiding future high-precision spectroscopic investigations.

In summary, this paper represents a milestone in molecular quantum electrodynamics, delivering a theoretical model for He$_2^+$ that is accurate enough to resolve the finest details of its spectrum, thereby closing the gap between theory and modern high-precision experiments.