Two-body $P$-state energies at $α^6$ order

This paper presents an analytical calculation of the complete $\alpha^6$ correction to the energies of $nP$-levels in two-body systems with arbitrary masses and magnetic moments, revealing a previously overlooked correction for positronium and providing essential theoretical data for extracting nuclear charge radii in light muonic atoms.

The Gist

Imagine the universe as a giant, cosmic dance floor. On this floor, particles like electrons, protons, and muons are constantly pairing up, spinning, and orbiting each other. These pairs are called two-body systems. The most famous dancer is the hydrogen atom (one proton, one electron), but there are exotic couples too, like positronium (an electron and its anti-matter twin, a positron) or muonic helium (a heavy muon orbiting a helium nucleus).

For decades, physicists have been trying to predict exactly how much energy these dancing pairs have. Why? Because if our predictions are perfect, we confirm our understanding of the universe's rules (Quantum Electrodynamics, or QED). If there's a tiny mismatch between our math and reality, it might mean there's a new, undiscovered force or particle hiding in the shadows.

This paper by Patkóš, Yerokhin, and Pachucki is like a master chef refining a recipe for the most complex dish imaginable: the energy levels of "P-states" (a specific type of orbit where the dancers are spinning in a figure-eight pattern) with extreme precision.

Here is the breakdown of what they did, using simple analogies:

1. The Problem: The "Fine-Tuning" Knob

Think of the energy of these atomic pairs as a radio station.

• The Main Signal ($E^{(0)}$): This is the basic tune you hear. It's the energy of the two particles just sitting there.
• The Static ($E^{(2)}$ to $E^{(4)}$): As you get closer to the station, you hear static. In physics, this is the "relativistic correction." It accounts for the fact that particles move fast and have spin. We've known how to tune this out for a long time.
• The Whisper ($E^{(6)}$): This is the paper's focus. It's a whisper so quiet it's buried under the static. It represents the sixth-order correction (denoted as $\alpha^6$). To hear this whisper, you need a microphone so sensitive it can detect the vibration of a single atom.

For a long time, scientists had a recipe for this whisper, but it had a flaw. It worked perfectly for dancers spinning in wide circles ($l > 1$), but when the dancers got into a tight, specific figure-eight spin ($l = 1$, or P-states), the recipe broke down. It was like having a map that works for driving on a highway but fails when you try to park in a tight spot.

2. The Solution: The "Microscope"

The authors built a new, ultra-precise mathematical microscope. They used a framework called NRQED (Non-Relativistic Quantum Electrodynamics).

Imagine you are trying to calculate the exact cost of a road trip.

• Old Method: You calculated the cost of gas, food, and hotels.
• The New Method: They realized they missed the cost of the toll booths and the tiny friction of the tires on the road.
• The "Contact Terms": In the old math, when the two particles got very close to each other (almost touching), the math got messy. It was like trying to count the grains of sand on a beach by looking at the whole beach at once. The authors realized that for P-states, you have to look at the "grains of sand" (the contact terms) individually. They found that previous calculations had accidentally canceled out two mistakes, which made the answer look right for simple cases (like positronium) but was actually wrong in the details. They fixed the math to get the true answer.

3. The Ingredients: Mass, Spin, and Size

The paper provides a "universal recipe" that works for almost any two-particle couple, regardless of:

• How heavy they are: Whether it's a light electron and a heavy proton, or two equal-mass particles.
• How they spin: Whether they are spin-0 (like a bowling ball) or spin-1/2 (like a spinning top).
• How "fuzzy" they are: Real particles aren't perfect points; they have a size (like a fuzzy tennis ball vs. a sharp pin). The authors included corrections for this "fuzziness" (finite size), which is crucial for heavy atoms like muonic helium.

4. Why It Matters: The "Ruler" of the Universe

Why do we care about this tiny whisper of energy?

• Measuring the Nucleus: In muonic atoms (where a muon orbits a nucleus), the muon is 200 times heavier than an electron. It orbits much closer to the nucleus, acting like a super-sensitive probe. By measuring the energy of these orbits, scientists can measure the size of the nucleus (the charge radius) with incredible precision.
• The "Proton Radius Puzzle": There was a famous mystery where measurements of the proton's size using electrons didn't match measurements using muons. This new, precise math helps resolve that mystery.
• Testing the Rules: If the experimental measurements of these energy levels match the authors' new formulas, it proves our current laws of physics are solid. If they don't match, it could be the first clue to "New Physics" beyond our current understanding.

5. The "Recoil" Effect

The paper also looked at nuclear recoil. Imagine a heavy boxer (the nucleus) and a lightweight boxer (the muon). When the lightweight boxer punches, the heavy boxer wobbles slightly. This wobble changes the energy of the system.
The authors compared their new "tight-parking" formula against a "heavy-duty" numerical simulation. They found that for light atoms, their simple formula works great. But for heavier atoms, the "wobble" gets complicated, and the simple formula starts to drift. Their work helps scientists know exactly when to use the simple recipe and when to switch to the heavy-duty computer simulation.

Summary

In short, this paper is a masterclass in precision. The authors took a complex, messy problem in quantum physics—calculating the energy of spinning particle pairs—and cleaned up the math. They fixed a hidden error in previous recipes, added the missing ingredients (like particle size and spin interactions), and provided a new, universal tool.

This tool allows scientists to measure the size of atomic nuclei with unprecedented accuracy and test the fundamental laws of the universe, ensuring that our "cosmic dance floor" is understood down to the very last step.

Technical Summary

1. Problem Statement

The paper addresses the need for high-precision theoretical predictions of energy levels in two-body atomic systems (such as hydrogen, positronium, muonium, and muonic atoms) to test Quantum Electrodynamics (QED) and determine fundamental constants (e.g., nuclear charge radii).

Specifically, the authors focus on calculating the complete $\alpha^6$ correction (where $\alpha$ is the fine-structure constant) to the binding energies of $nP$-states ($l=1$).

• Context: Previous calculations for $l > 1$ states existed, but $l=1$ states present unique challenges due to the presence of contact terms (delta-function interactions) that vanish for higher angular momenta.
• Gap: There were discrepancies and potential errors in previous literature regarding the treatment of these contact terms for $l=1$ states, particularly in positronium calculations. Additionally, a complete analytical formula for arbitrary mass ratios and extended-size particles (finite nuclear size) at this order was lacking.

2. Methodology

The authors employ Non-Relativistic QED (NRQED) formalism to derive the effective Hamiltonian for the $\alpha^6$ correction.

• Framework: The energy expansion is expressed as $E(\alpha) = E^{(0)} + E^{(2)} + E^{(4)} + E^{(5)} + E^{(6)} + \dots$. The paper focuses on deriving $E^{(6)}$.
• Hamiltonian Derivation:
  • They start with the NRQED Hamiltonian for particles with arbitrary spin ($s=0$ or $1/2$), mass, and magnetic moments.
  • They derive the effective operator $H^{(6)}$ by evaluating one-photon exchange diagrams between two particles.
  • The calculation is performed in the Coulomb gauge using the non-retardation approximation ($k_0 = 0$) for the majority of terms, with retardation corrections ($\delta E_7, \delta E_8, \delta E_9$) treated separately.
• Operator Classification: The effective Hamiltonian $H^{(6)}$ is decomposed into a sum of 10 distinct contributions ($\delta H_0$ to $\delta H_9$), including:
  • Kinetic energy corrections ($p^6$ terms).
  • Corrections to Coulomb interaction involving finite-size parameters (charge radius $r_E$, magnetic radius $r_M$, polarizability $\alpha_E$).
  • Spin-orbit and spin-spin interactions.
  • Retardation effects.
• Matrix Elements: The authors calculate the expectation values of these operators for $nP$-states. A critical part of the derivation involves handling contact terms (terms proportional to $\delta^3(\vec{r})$), which are non-zero only for $l=1$.
• Verification:
  • The results are compared against previous general formulas for $l>1$ (Ref. [10]).
  • The authors identify and correct a "double error" cancellation in previous positronium calculations (Ref. [13–15]), confirming that while the previous results were numerically correct for point particles, the derivation was flawed.
  • Numerical validation is performed by comparing the derived analytical formulas against all-order in $Z\alpha$ numerical calculations of the nuclear recoil correction for muonic atoms.

3. Key Contributions

1. Complete Analytical Formulas: The paper provides the first complete analytical expressions for the $\alpha^6$ energy correction for $nP$-states in two-body systems with:
   • Arbitrary mass ratios ($m_1/m_2$).
   • Arbitrary spins ($s=0, 1/2$).
   • Finite particle sizes (extended charge and magnetic radii).
2. Resolution of Contact Term Discrepancies: The authors rigorously derive the contact terms for $l=1$, clarifying that they differ from the $l>1$ case. They demonstrate that previous literature results for positronium were correct by accident due to the cancellation of two distinct errors (neglect of a local term and improper inclusion of intermediate states).
3. Finite-Size Recoil Corrections: New terms describing the interplay between finite nuclear size and recoil effects are derived, which are crucial for light muonic atoms.
4. Validation via Numerical Comparison: The analytical results are benchmarked against all-order numerical calculations for muonic atoms, showing rapid convergence as the nuclear charge $Z$ decreases, thereby validating the perturbative approach.

4. Key Results

The final energy correction is expressed as:
$$E^{(6)} = (Z\alpha)^6 \left[ E_{NS} + \vec{s}_1 \cdot \vec{s}_2 E_{SS} + \vec{L} \cdot \vec{s}_1 E_{L1} + \vec{L} \cdot \vec{s}_2 E_{L2} + (L_i L_j)_{(2)} s_{1i} s_{2j} E_{LL} \right]$$

• Specific Cases:
  • Spin-0/Spin-0: The result differs from the general $l>1$ formula only by finite-size terms.
  • Spin-0/Spin-1/2: Formulas derived for systems like muonic hydrogen (where the nucleus is spin-0 or spin-1/2).
  • Spin-1/2/Spin-1/2: The most complex case (e.g., positronium, muonium), including hyperfine splitting and mixing terms.
• Muonic Helium ($\mu$He):
  • The authors applied their formulas to calculate the $2P$ fine structure of $\mu^3$He$^+$ and $\mu^4$He$^+$.
  • Table I shows the theoretical prediction for the fine structure splitting: $144.785(3)$ meV for $\mu^3$He and $146.182(3)$ meV for $\mu^4$He.
  • These values agree with previous calculations (Ref. [22, 23]) and experimental data, supporting the extraction of nuclear charge radii.
• Nuclear Recoil:
  • Table II and Figure 1 present the nuclear recoil corrections for muonic atoms up to $Z=40$.
  • The study reveals that for $Z > 10$ (for $2p_{1/2}$) and $Z > 20$ (for $2p_{3/2}$), the finite-nuclear-size (fns) recoil correction dominates over the point-nucleus recoil correction.
  • The lowest-order analytical approximation overestimates the fns correction by a factor of ~2 for $Z=40$, highlighting the necessity of all-order numerical methods for heavy muonic ions.

5. Significance

• Precision Physics: The derived formulas are essential for reducing theoretical uncertainties in the determination of nuclear charge radii from muonic atom spectroscopy. This is particularly relevant for resolving the "proton radius puzzle" and discrepancies between electronic and muonic spectroscopy.
• QED Tests: The results provide a rigorous test of QED in the strong-field regime (high $Z$) and for systems with comparable masses (positronium, muonium).
• Future Applications: The framework is applicable to exotic systems like protonium and other hadronic atoms, offering a tool to probe long-range interactions between hadronic particles.
• Methodological Correction: By identifying and correcting the accidental error cancellation in previous positronium calculations, the paper ensures the theoretical foundation for future high-precision tests in these systems is sound.

In summary, this work completes the theoretical description of $nP$-state energies at order $\alpha^6$ for a broad class of two-body systems, resolving previous ambiguities and providing the necessary tools for next-generation precision spectroscopy experiments.