Relativistic nuclear recoil effects in hyperfine splitting of hydrogenic systems

This paper calculates finite nuclear mass corrections to the hyperfine splitting in hydrogenic systems using a combined relativistic and nonrelativistic QED approach, finding results that disagree with previous calculations and reveal a $2\sigma$ discrepancy with hydrogen measurements that may point to issues with proton structure corrections.

The Gist

Imagine you are trying to tune a radio to a very specific, perfect frequency. In the world of atoms, this "frequency" is the hyperfine splitting—a tiny energy difference in the hydrogen atom that acts like a super-precise clock. Scientists have measured this clock with incredible accuracy, but when they try to calculate what the clock should read using the laws of physics, the numbers don't quite match. There's a tiny gap, a "static" in the signal.

This paper by Jakub Hevler, Andrzej Maroń, and Krzysztof Pachucki is about fixing a specific part of that calculation to see if we can clear up the static.

Here is the story of their work, explained without the heavy math.

The Problem: The "Heavy" Dancer

In a hydrogen atom, you have a tiny electron dancing around a much heavier proton (the nucleus).

• The Old Way of Thinking: For a long time, physicists treated the proton like a giant, immovable anchor. They assumed the electron danced around a fixed point, and the proton just sat there.
• The Reality: The proton isn't an anchor; it's a heavy dancer. When the electron spins and jumps, the proton wobbles back and forth to keep the balance. This "wobble" is called nuclear recoil.

Calculating this wobble is hard because it involves two things fighting each other:

1. Relativity: The electron moves so fast it needs Einstein's rules.
2. Quantum Mechanics: The proton is heavy but still follows quantum rules.

The Conflict: The "Bodwin and Yennie" Map

For decades, the best map for calculating this wobble was drawn by two physicists, Bodwin and Yennie, back in 1988. Everyone used their map. But recently, as our atomic clocks got more precise, the map started to look wrong. The calculated "wobble" didn't match the real-world measurements.

The authors of this paper decided to redraw the map. They asked: "Did we calculate the recoil effect correctly, or did we miss a step in the dance?"

The New Calculation: Two Different Lenses

To get the answer right, the authors used two different "lenses" (mathematical frameworks) to look at the same problem, hoping they would agree with each other.

1. Lens 1: The Heavy Particle View (HPQED)
   Imagine looking at the proton as a massive boulder and the electron as a pebble. This method treats the proton as a "heavy particle" and calculates how the electron's energy changes when the boulder moves slightly.
2. Lens 2: The Effective Field View (NRQED)
   This is like looking at the dance floor from a distance. Instead of tracking every single step of the electron and proton, this method creates a simplified "effective" rulebook that summarizes how they interact at low energies.

The Result: Both lenses gave the same answer, but it was different from the old 1988 map.

The Discovery: A New Twist in the Dance

The authors found that the old calculation missed a subtle interaction. It's like realizing that when the heavy dancer (proton) wobbles, it doesn't just move back and forth; it also changes the shape of the dance floor slightly in a way the old map didn't account for.

Their new calculation changes the predicted energy of the hydrogen atom.

• Old Prediction: The map said the energy should be $X$.
• New Prediction: The map says the energy should be $X + \text{a tiny bit}$.

The "2 Sigma" Mystery

Here is the twist: Even with this new, more accurate map, the numbers still don't perfectly match the experimental measurements.

• The gap has shrunk, but it's still there. It's a "2-sigma" discrepancy. In the world of science, this is like hearing a faint hum in a quiet room. It's not loud enough to be a scream (which would be a 5-sigma discovery), but it's definitely there and annoying.

What does this mean?
The authors suggest that the remaining error isn't in their math about the "wobble" (recoil). Instead, the error is likely in our understanding of the proton itself.

Think of the proton not as a solid marble, but as a fuzzy cloud of quarks and gluons. We don't know the exact shape and "squishiness" of this cloud perfectly. The authors suspect that the "fuzziness" of the proton (its internal structure) is the missing piece of the puzzle.

Why Should You Care?

1. Testing the Universe: Hydrogen is the simplest atom in the universe. If our math doesn't work for the simplest thing, it means our fundamental understanding of how the universe works (Quantum Electrodynamics) might have a crack in it.
2. The Muon Clue: The paper suggests a way to solve the mystery: look at Muonic Hydrogen. This is an atom where the electron is replaced by a muon (a heavier cousin of the electron). Because the muon is heavier, it orbits closer to the proton, making the "proton fuzziness" much more obvious. If we measure Muonic Hydrogen with the new math, we might finally figure out the true shape of the proton.

The Bottom Line

The authors fixed a specific error in how we calculate the "wobble" of the proton in a hydrogen atom. They proved the old map was wrong. However, the mystery isn't fully solved yet. The remaining mismatch suggests that the proton is more complex and "fuzzy" than we thought. It's a reminder that even in the simplest atom, there are still secrets waiting to be discovered.

Technical Summary

1. Problem Statement

The paper addresses a persistent discrepancy between theoretical predictions and experimental measurements of the ground-state hyperfine splitting (HFS) in hydrogen.

• The Discrepancy: Previous theoretical calculations of relativistic nuclear recoil effects (finite nuclear mass corrections) of order $(Z\alpha)^2 (m/M) E_F$ disagreed with the experimental value by several standard deviations ($\sigma$).
• The Conflict: The authors identify that their new calculations contradict the seminal work by Bodwin and Yennie (1988), which has been the standard reference for these corrections.
• The Implication: The existing discrepancy was often attributed to uncertainties in proton structure (specifically the Zemach radius and polarizability). However, if the recoil corrections themselves are incorrect, the attribution of the error to proton structure is flawed. Resolving this is crucial for testing the Standard Model and determining fundamental constants.

2. Methodology

The authors employ a dual-approach strategy, combining Nonrelativistic QED (NRQED) and Heavy Particle QED (HPQED), to calculate the recoil corrections to order $(Z\alpha)^2 (m/M) E_F$.

• Framework:

  • HPQED (High-Energy Part): They utilize exact non-perturbative formulas derived by Shabaev and extended by Pachucki for relativistic recoil. This involves calculating the two-photon exchange forward scattering amplitude. The calculation is split into high-energy ($k \sim m$) and low-energy ($k \sim m\alpha$) parts.
  • NRQED (Low-Energy Part): For the low-energy sector, they construct an effective Hamiltonian ($H^{(6)}$) using NRQED formalism. This allows for a systematic expansion in $Z\alpha$ and $m/M$, handling the bound-state dynamics more efficiently than direct QED integration.

• Key Technical Steps:

  • Dimensional Regularization: The calculations are performed in $d = 3 - 2\epsilon$ dimensions to handle ultraviolet and infrared divergences.
  • Separation of Scales: The recoil correction is separated into:
    • High-Energy ($E_H$): Calculated via HPQED using hard three-photon exchange approximations.
    • Low-Energy ($E_L$): Calculated via NRQED using effective operators (spin-spin, spin-orbit, and kinetic terms).
  • Form Factors: The nuclear structure is treated via electromagnetic form factors ($G_E, G_M$). The authors explicitly address the definition of the Zemach radius correction, arguing for the use of the lepton mass $m$ rather than the reduced mass $\mu$ in specific terms to avoid double-counting or missing terms in the expansion.

3. Key Contributions

• Correction of Bodwin and Yennie (1988): The primary contribution is the recalculation of the $(Z\alpha)^2 (m/M)$ recoil term. The authors find that the result in Bodwin and Yennie is incorrect.
  • They derived a new analytical formula for the $1S$ state recoil correction (Eq. 89) which differs significantly in the coefficients of the $\ln 2$ and $\ln(Z\alpha)$ terms compared to the 1988 result (Eq. 90).
• Unified Derivation: The paper provides a rigorous derivation that unifies the HPQED and NRQED approaches, verifying that the sum of the high-energy and low-energy parts yields a consistent, finite result.
• Re-evaluation of Proton Structure: By correcting the recoil term, the authors re-evaluate the "proton structure" contribution. They demonstrate that the remaining discrepancy between theory and experiment is smaller ($\approx 2\sigma$) but still present, suggesting the issue likely lies in the poorly known proton structure (polarizability) rather than the QED recoil calculation itself.

4. Results

• New Recoil Formula:
  The corrected relativistic recoil correction for the $1S$ state is given by:
  $$E_{rec}^{(6)}(1S) = \left[ \frac{65}{18} + \frac{13}{18}\kappa + \frac{31}{36}\kappa^2 - \left(8 + 2\kappa - \frac{1}{4}\kappa^2\right)\ln 2 - \left(2 + 2\kappa + \frac{7}{4}\kappa^2\right)\ln(Z\alpha) \right] (Z\alpha)^2 \frac{m}{M} \frac{E_F}{1+\kappa}$$
  This differs from the Bodwin-Yennie result, particularly in the coefficients of the logarithmic terms.

• Impact on Hydrogen HFS:

  • The new theoretical prediction for the hydrogen ground state HFS reduces the discrepancy with the experimental value ($1420405.751768$ kHz) from a larger gap to approximately $2\sigma$.
  • The authors attribute this remaining $2\sigma$ difference to uncertainties in the proton structure corrections (specifically the Zemach radius and polarizability), rather than QED recoil errors.

• Specific Difference $D_{21}$:
  The paper also analyzes the specific difference $D_{21} = 8 E_{hfs}(2S) - E_{hfs}(1S)$.

  • For Hydrogen ($^1$H): The new calculation worsens the agreement with current experimental data for $D_{21}$, suggesting a need for more precise measurements of the $2S$ HFS.
  • For Helium Ion ($^3$He$^+$): The new result significantly improves the agreement between theory and experiment, validating the calculation for systems with different nuclear charges.

• Muonic Hydrogen ($\mu$H):
  The authors note that their recoil corrections are sufficiently accurate for muonic hydrogen, which is critical for future tests of fundamental interactions where proton structure effects are amplified.

5. Significance

• Resolution of a Decades-Old Discrepancy: The paper resolves a long-standing conflict in atomic physics by correcting a fundamental error in the 1988 Bodwin-Yennie calculation.
• Refining Fundamental Constants: By isolating the recoil error, the authors clarify that the remaining tension in hydrogen HFS is likely due to proton structure. This directs future experimental efforts (e.g., in muonic hydrogen) toward better constraining the proton's electromagnetic form factors and polarizability.
• Methodological Benchmark: The successful combination of HPQED and NRQED provides a robust framework for calculating higher-order recoil effects in other hydrogenic and muonic systems, including heavy ions.
• Future Directions: The paper highlights that the $2\sigma$ discrepancy in hydrogen can be definitively tested by comparing hydrogen and muonic hydrogen ($\mu$H) measurements, where the recoil and structure effects scale differently, potentially canceling out large uncertainties.

In conclusion, this work represents a critical update to the theoretical foundation of atomic spectroscopy, ensuring that tests of the Standard Model using hydrogen-like systems are based on accurate QED recoil corrections.