Recoil corrections to $\mu$H hyperfine splitting

This paper presents a complete theoretical calculation of the muonic hydrogen ($\mu$H) ground-state hyperfine splitting, incorporating direct quantum electrodynamic and recoil corrections alongside proton structure effects derived from hydrogen data, to yield a prediction of $182\,626(5)$ $\mu$eV with precision exceeding 1 ppm.

The Gist

Imagine the hydrogen atom as a tiny solar system. Usually, it has a heavy sun (the proton) and a light, speedy planet (the electron) orbiting it. But in Muonic Hydrogen (µH), we swap the electron for a muon.

A muon is like a "heavy electron." It's about 200 times heavier than an electron. Because it's so heavy, it doesn't orbit far away; it dives deep, hugging the proton tightly, almost touching it. This makes the atom a perfect, ultra-sensitive laboratory for testing the laws of physics.

The Big Mystery: The "Hyperfine Splitting"

Inside this tiny atom, the proton and the muon both spin like little tops. When they spin in the same direction, the atom has a tiny bit more energy than when they spin in opposite directions. The difference between these two energy states is called Hyperfine Splitting.

Think of it like a piano key. If you press it just right, it makes a specific note. Physicists want to know the exact frequency of that note for Muonic Hydrogen. If our theoretical prediction matches the real-world measurement perfectly, our understanding of the universe (the Standard Model) is correct. If they don't match, it means there's a hidden piece of the puzzle we're missing.

The Problem: The "Fuzzy" Proton

For regular hydrogen (with an electron), we've measured this note incredibly precisely. But when we try to calculate it for Muonic Hydrogen, things get messy.

Why? Because the muon is so close to the proton that it can "feel" the proton's internal structure. The proton isn't a perfect, solid marble; it's a fuzzy cloud of quarks and gluons.

• The Analogy: Imagine trying to measure the weight of a bowling ball by bouncing a ping-pong ball off it. The ping-pong ball (electron) bounces off the surface and doesn't care what's inside. But if you use a tiny, heavy cannonball (muon) that crashes right into the bowling ball, the result depends entirely on how squishy or hard the inside of the bowling ball is.

The paper by Maroń, Pańtak, and Pachucki is essentially a massive accounting project. They are trying to calculate every single tiny force that could shift that "piano note" by more than 1 part per million (1 ppm).

What Did They Do?

The authors went through the "receipts" of physics, calculating corrections from three main categories:

1. The Vacuum is Not Empty (Vacuum Polarization):
   In quantum physics, empty space is actually bubbling with virtual particles popping in and out of existence. These particles act like a fog around the proton.

   • Analogy: Imagine the proton is a lighthouse. The "fog" of virtual particles bends the light. The authors calculated exactly how much this fog shifts the muon's orbit. They found that the "electron fog" and the "muon fog" both play a role, and they had to add them up precisely.

2. The "Kick" (Recoil Corrections):
   When the muon orbits, the proton isn't perfectly still; it wobbles a bit, like a mother holding a spinning child. Because the muon is heavy, this wobble is much more significant than in regular hydrogen.

   • Analogy: If you push a shopping cart (light electron), the cart moves, but you barely feel a kickback. If you push a heavy boulder (muon), the boulder pushes back hard. The authors calculated this "kickback" effect with extreme precision.

3. The Proton's Shape (Finite Nuclear Size):
   Since the muon gets so close, it sees the proton's actual size and shape, not just a mathematical point.

   • Analogy: If you are far away, a person looks like a dot. If you get right up in their face, you see their nose, eyes, and ears. The authors calculated how the muon "seeing" the proton's nose and ears changes the energy of the atom.

The "Cheat Code": Using Regular Hydrogen

One of the cleverest parts of this paper is how they handled the biggest unknown: The Proton's Fuzziness.

Calculating the exact "fuzziness" (Zemach radius) of a proton from scratch is incredibly hard, like trying to predict the weather on a planet you've never visited.

• The Solution: They used the fact that we already know the answer for regular hydrogen very well. They took the difference between the "fuzzy" proton effect in regular hydrogen and the "fuzzy" effect in muonic hydrogen.
• The Analogy: Imagine you are trying to guess the weight of a mystery box. You know the weight of a similar box with a known item inside. By comparing the two, you can cancel out the "box weight" and isolate the "mystery item." They used the known data from regular hydrogen to cancel out the biggest uncertainties in the muonic calculation.

The Result

After adding up thousands of tiny corrections (some positive, some negative, some canceling each other out), they arrived at a final prediction:
The energy difference is 182,626 micro-electron-volts.

They claim this prediction is accurate to within 5 parts per million.

Why Does This Matter?

This paper provides the "theoretical target" for experimentalists.

• If scientists measure the muonic hydrogen note and it matches this number, it confirms our understanding of the universe is solid.
• If it doesn't match, it could mean there is New Physics waiting to be discovered—perhaps a new particle or a new force that we haven't seen yet.

In short, these authors built the most precise "map" possible for a tiny, exotic atom, helping us navigate the frontier between what we know and what we are about to discover.

Technical Summary

1. Problem Statement

The ground-state hyperfine splitting (HFS) in regular hydrogen ($H$) serves as a precision test of the Standard Model. However, a persistent $2\sigma$ discrepancy exists between the most recent theoretical predictions and extremely accurate experimental measurements. The primary source of this discrepancy is believed to be the uncertainty in proton structure corrections (specifically the Zemach radius and proton polarizability), which are difficult to calculate from first principles (e.g., Lattice QCD is not yet precise enough).

To resolve this, the authors propose studying muonic hydrogen ($\mu$H), where the electron is replaced by a muon. Due to the muon's mass being $\approx 200$ times larger than the electron's, the muon orbits much closer to the proton, making the system significantly more sensitive to nuclear structure effects. However, this also amplifies nuclear recoil corrections (due to the muon-proton mass ratio $m_\mu/M \approx 0.1$, compared to $m_e/M \approx 0.0005$ in regular hydrogen).

The paper aims to provide a complete theoretical calculation of the $\mu$H ground-state HFS, including all contributions larger than 1 ppm, with a specific focus on direct calculations of Quantum Electrodynamics (QED) and recoil corrections.

2. Methodology

The authors employ a rigorous perturbative approach, expanding the HFS energy in terms of the fine-structure constant $\alpha$ and the mass ratio $m/M$. The total energy is expressed as:
$$E_{hfs} = E_F (1 + \delta)$$
where $E_F$ is the Fermi energy and $\delta$ represents the sum of all corrections.

Key methodological components include:

• Vacuum Polarization (VP):
  • Electron VP (EVP): Calculated directly for a point-like nucleus using the Uehling potential (Yukawa-like) and numerical integration. Both one-loop and two-loop EVP corrections are computed.
  • Muon VP ($\mu$VP): Calculated similarly, including combined effects with EVP.
  • Relativistic EVP: Calculated using Dirac wave functions to account for relativistic recoil effects.
• Self-Energy (SE):
  • Calculated using the two-photon exchange forward scattering amplitude formalism.
  • The authors derive formulas for SE combined with recoil and Finite Nuclear Size (FNS) effects without expanding in the mass ratio $m/M$.
• Two-Photon Exchange (Recoil & FNS):
  • The recoil correction is derived using the temporal gauge and Wick rotation, separating the amplitude into point-nucleus, finite-size, and recoil components.
  • Proton Structure: Instead of calculating proton structure from scratch (which is uncertain), the authors utilize the difference between $\mu$H and $H$. They extract the proton structure contribution from the precise experimental HFS of regular hydrogen and apply a scaling factor ($m_\mu/m_e$) to $\mu$H. This technique cancels out much of the uncertainty associated with the Zemach radius and polarizability.
• Form Factors: The calculations utilize a dipole parametrization for proton form factors, adjusted to match the known Zemach radius ($r_Z = 1.054$ fm).

3. Key Contributions

• Comprehensive Recoil Calculation: This is likely the first comprehensive calculation of $\mu$H HFS that explicitly identifies and calculates all contributions above 1 ppm, with a heavy emphasis on nuclear recoil corrections which are negligible in regular hydrogen but dominant here.
• Direct QED Calculation: Unlike previous works that relied heavily on extrapolations from muonium or hydrogen, this work calculates QED corrections (EVP, SE, VP+SE) directly for the $\mu$H system.
• Hybrid Theory-Experiment Approach: The authors introduce a method to improve the theoretical prediction for $\mu$H by leveraging the known discrepancy in regular hydrogen. By calculating the difference $\Delta = \delta(\mu H) - (m_\mu/m_e)\delta(H)$, they effectively isolate the nuclear structure uncertainties, reducing the error budget significantly.
• New Formulas: The paper derives new analytical and numerical formulas for:
  • VP combined with recoil and FNS.
  • SE combined with recoil and FNS.
  • Relativistic recoil corrections to order $(Z\alpha)^2$.

4. Key Results

The authors present a table of contributions (Table I) and a final theoretical prediction:

• Total Correction Factor: The sum of all corrections yields $\delta \approx 0.000870(60)$.
• Corrected Prediction: After applying the correction derived from the $H$ vs. $\mu$H difference (Eq. 110), the final theoretical prediction for the ground-state HFS is:
  $$E_{hfs}(\mu H) = 182,626(5) \, \mu\text{eV}$$
  (or $0.182626(5)$ eV).
• Specific Corrections:
  • Leading Order EVP: $\delta^{(1)}_{evp} \approx 6075$ ppm.
  • Finite Nuclear Size (FNS): A massive correction of $\approx -149.8$ ppm (EVP) and $+16.2$ ppm (SE).
  • Recoil: The leading recoil correction is $\approx 1672$ ppm.
  • Proton Polarizability: Estimated at $\approx 200$ ppm, but the difference between $\mu$H and $H$ reduces this uncertainty to a negligible level in the final corrected value.
• Comparison: The result is close to the non-relativistic formula with the muon magnetic moment anomaly ($0.182656$ eV), suggesting a cancellation of higher-order corrections. The authors note discrepancies with previous work by Faustov and Martynenko regarding term-by-term classification.

5. Significance

• Standard Model Test: The predicted value of $182,626(5) \, \mu\text{eV}$ provides a precise target for future experiments (e.g., the FAMU experiment). A measurement with 1 ppm accuracy would constitute a stringent test of fundamental interactions and QED in the presence of a heavy lepton.
• Proton Structure Determination: If the experimental measurement matches the theory, it could be used to determine the proton Zemach radius with unprecedented precision, potentially resolving the "proton radius puzzle" and the discrepancies in regular hydrogen.
• Methodological Benchmark: The paper establishes a rigorous framework for calculating recoil and finite-size effects in muonic atoms, which is essential for interpreting future high-precision spectroscopy data.
• Future Work: The authors identify that the remaining uncertainties lie in the $(Z\alpha)^2$ corrections (evaluated without mass ratio expansion) and the weighted difference of nuclear structure terms, which require further independent verification.

In summary, this work provides the most complete theoretical baseline for muonic hydrogen hyperfine splitting to date, bridging the gap between QED theory and nuclear structure physics to enable high-precision tests of the Standard Model.