Binding Energy of Muonic Beryllium: Perturbative versus All--Order Calculations

This paper demonstrates that both perturbative and all-order relativistic approaches yield consistent ground-state binding energies for muonic beryllium to within one part per million, thereby providing a precise parametrization for extracting the $^9$Be charge radius and bridging theoretical methodologies between light and heavy muonic systems.

The Gist

The Big Picture: A Tiny Planet and a Heavy Star

Imagine the atom as a solar system. Usually, we have a heavy sun (the nucleus) and a tiny, fast planet (the electron) orbiting it. But in muonic atoms, we swap the electron for a muon.

A muon is like a "heavy electron." It's about 200 times heavier. Because it's so heavy, it doesn't orbit far away; it spirals in tight, hugging the nucleus like a satellite in a very low orbit. In fact, it gets so close that it can actually "feel" the texture of the nucleus itself, not just its electric charge.

This paper is about Beryllium-9, a specific type of atom. The scientists wanted to calculate exactly how much energy it takes to keep this heavy muon stuck to the Beryllium nucleus. Why? Because if we know the energy perfectly, we can work backward to measure the size of the nucleus (its "charge radius") with incredible precision. This helps us understand the building blocks of the universe and even test if there are new, hidden forces in nature.

The Problem: Two Different Maps for the Same Territory

The scientists faced a dilemma. To calculate the energy of this system, there are two main ways to do the math:

1. The "Light System" Approach (Perturbative): This is like calculating the distance between two cities by adding up small, straight-line segments. It works great for light atoms (like Hydrogen) where the nucleus is tiny and the muon is far away. You treat the nucleus as a point and add small corrections for its size.
2. The "Heavy System" Approach (All-Order): This is like using a high-resolution GPS that accounts for every bump, hill, and curve in the road immediately. This is usually used for heavy atoms (like Gold or Lead) where the nucleus is huge and the muon is squished right up against it. You treat the nucleus as a big, fuzzy cloud and solve the equations all at once.

The Conflict: For a medium-sized atom like Beryllium, scientists weren't sure which map was better. The "Light" method might miss some big bumps, and the "Heavy" method might get confused by the recoil (the wobble) of the nucleus.

The Experiment: A "Tug-of-War" of Math

The authors decided to run a massive comparison. They calculated the energy of the muonic Beryllium atom using both methods independently.

• Method A (The Light Approach): They started with a simple formula and added correction after correction, like stacking blocks, to account for the nucleus's size and the muon's wobble.
• Method B (The Heavy Approach): They used a super-complex, "all-in-one" equation that treats the nucleus as a real, physical object with a specific shape (a Gaussian cloud) and solves the physics without breaking it into small steps.

The Result: It was a tie! Both methods gave the exact same answer, down to a tiny fraction of a percent (better than one part per million).

The Analogy: Imagine two chefs trying to bake the exact same cake.

• Chef A measures every ingredient with a teaspoon, one by one (Perturbative).
• Chef B dumps everything into a high-tech mixer that calculates the perfect blend instantly (All-Order).
• They taste the cakes, and they are identical. This proves that both methods are reliable, and we can trust the "All-Order" method to work even for lighter atoms, not just heavy ones.

The "Wobble" Factor (Recoil)

One tricky part of the math is recoil. When the muon orbits, it pulls on the nucleus. Since the nucleus isn't infinitely heavy, it wobbles a little bit (like a mother holding a spinning child).

• In the "Heavy" approach, this wobble is often ignored or simplified.
• In the "Light" approach, the wobble is a major part of the calculation.

The paper shows that even though the "Heavy" approach usually ignores the wobble, when you add the specific corrections for Beryllium, it matches the "Light" approach perfectly. This bridges the gap between the two communities of physicists who usually work on different types of atoms.

The Payoff: A New Ruler for the Universe

The main goal wasn't just to do the math; it was to create a tool for experimentalists.

The authors created a simple formula (a "parametrization"). Think of it like a conversion chart.

• Before: If an experiment measured an energy, scientists had to run complex simulations to guess the size of the nucleus.
• Now: Scientists can plug the measured energy into this simple equation, and it instantly tells them the size of the Beryllium nucleus.

This is crucial because recent experiments are measuring these atoms with lasers and X-rays. If we know the energy perfectly, we can measure the size of the nucleus with extreme precision. This helps answer big questions:

• Is the Standard Model of physics correct?
• Are there new, invisible particles (like a "dark boson") messing with the muon?

Summary in a Nutshell

1. The Goal: Calculate the energy of a heavy muon orbiting a Beryllium nucleus to measure the nucleus's size.
2. The Method: They compared two different mathematical "languages" (one for light atoms, one for heavy atoms).
3. The Discovery: Both languages say the exact same thing. They agree perfectly.
4. The Impact: This gives scientists confidence to use the powerful "heavy atom" math for lighter atoms, and provides a simple formula for experimentalists to extract precise nuclear sizes from their data.

It's a victory for precision, proving that whether you look at the universe through a microscope or a telescope, the laws of physics remain consistent.

Technical Summary

1. Problem Statement

Muonic atoms (atoms where an electron is replaced by a muon) serve as precision probes for nuclear structure and tests of Quantum Electrodynamics (QED). However, theoretical approaches for calculating their energy levels have historically diverged based on the nuclear charge number ($Z$):

• Light Systems ($Z=1, 2$): Calculations predominantly use Nonrelativistic QED (NRQED), treating the nucleus as a point charge initially and adding finite-size effects perturbatively. This approach handles recoil corrections (finite nuclear mass) rigorously.
• Heavy Systems ($Z \gtrsim 30$): Calculations use fully relativistic approaches (solving the Dirac equation) where the finite nuclear size is treated non-perturbatively (all-order) from the start, often assuming an infinite nuclear mass.
• The Gap: For intermediate systems like muonic Beryllium ($Z=4$), neither approach is clearly superior without cross-validation. There is a need to bridge the theoretical communities and ensure that the "all-order" methods used for heavy ions are robust enough for lighter systems, while verifying that perturbative expansions remain valid.

The specific objective of this work is to compute the ground-state binding energy of muonic $^9$Be using both methods to compare them term-by-term, validate their agreement, and provide a high-precision parametrization for extracting the nuclear charge radius from experimental data.

2. Methodology

The authors performed two independent calculations for the ground-state ($1S$) binding energy of muonic $^9$Be:

A. The "All-Order" (Hybrid/Relativistic) Approach

• Framework: Solves the Dirac-Coulomb equation self-consistently for a muon bound to a nucleus with a finite charge distribution (modeled as a Gaussian).
• Finite Size: Treated non-perturbatively (to all orders in the finite-size parameter $\rho$).
• Vacuum Polarization: The leading electronic vacuum polarization (eVP) is included via the Uehling potential added to the nuclear potential, solved self-consistently.
• Recoil: Nonrelativistic recoil is handled by numerical calculations using the reduced mass ($m_r$) in the limit of infinite light speed.
• Tools: Utilized the Multiconfiguration Dirac Fock and General Matrix Elements (MDFGME) code.

B. The "Fully Perturbative" (NRQED) Approach

• Framework: Starts with the nonrelativistic Bohr energy for a point nucleus.
• Finite Size: Expanded as a perturbation series in the finite-size parameter $\rho$ (one-photon, two-photon, and three-photon exchange terms).
• Vacuum Polarization: Calculated using perturbation theory (First-Order and Second-Order) on top of the point-nucleus wavefunctions.
• Recoil: Treated by scaling terms with the reduced mass ratio ($m_r/m_\mu$) and calculating the difference between finite-mass and infinite-mass results.

C. Comparison Strategy

The authors calculated every contribution (Dirac-Coulomb, Uehling, Self-Energy, Källén-Sabry, Recoil, etc.) using both methods. They compared the results term-by-term to identify discrepancies and trace them to specific physical approximations (e.g., missing higher-order photon exchanges or model dependencies).

3. Key Contributions

A. Validation of Theoretical Frameworks

The study demonstrates that the fully relativistic, all-order approach (typically reserved for heavy ions) can be successfully adapted for light systems ($Z=4$) and agrees with the perturbative NRQED approach to better than one part-per-million of the total energy.

• Total Energy Agreement: The difference between the two methods is approximately 41 meV out of a total binding energy of $\approx -44,513$ eV.
• Source of Discrepancy: The small difference is attributed to:
  1. Missing elastic four-photon (and higher) exchange corrections in the perturbative expansion.
  2. Model dependence (e.g., using a hard-sphere model vs. Gaussian model for subleading mixed finite-size/eVP terms in the perturbative approach).

B. Comprehensive Error Budget

The paper provides a rigorous breakdown of all energy contributions, including:

• Dirac-Coulomb-Uehling: $\approx -45,072$ eV.
• Nonrelativistic Recoil: $\approx +556.5$ eV.
• Radiative Corrections: Includes Källén-Sabry (two-loop eVP), Muon Self-Energy (SE), Muonic Uehling, and Hadronic Uehling.
• Recoil-Radiative Corrections: Calculated for eVP and SE, showing that while small, they are necessary for sub-eV precision.

C. Parametrization for Experimental Extraction

The authors derived a simple, high-precision parametrization equation to relate the $2P-1S$ transition energy to the nuclear charge radius ($r_C$).

• Equation: $E_{2P-1S}(r_C) = 33,392.55(2) - 14.32(r_C^2 - 2.519^2) + 0.389(r_C^3 - 2.519^3) + \Delta E_{NS}$.
• Validity: Accurate within 10 meV for $r_C$ variations of $\pm 5\%$ around the reference value of 2.519 fm.
• Uncertainty: The pure QED uncertainty is quoted as 20 meV, which is significantly smaller than the current uncertainties in nuclear polarization corrections ($\approx 0.1-0.2$ eV).

4. Key Results

• Recommended Binding Energy: The final recommended ground-state binding energy for muonic $^9$Be is:
  $$E_{1S} = -44,513.07(2) \text{ eV}$$
  The uncertainty (20 meV) arises from numerical convergence and ambiguities in relativistic quadratic recoil and nuclear self-energy definitions.
• Transition Energy: The $2P-1S$ transition energy is calculated as $33,392.55(2)$ eV (excluding nuclear structure corrections $\Delta E_{NS}$).
• Recoil Effects: The nonrelativistic recoil correction is substantial ($+556.5$ eV), while relativistic recoil corrections are negligible ($< 0.1$ eV) for this system, confirming that the nonrelativistic recoil approximation is highly accurate for $Z=4$.

5. Significance

• Bridging Communities: This work acts as a "bridge" between theorists specializing in light muonic atoms (using NRQED) and those working on heavy ions (using all-order Dirac methods). It proves that the all-order relativistic framework is robust and can be extended to cover all charge numbers $Z$ with high precision.
• Experimental Utility: The provided parametrization allows experimentalists at facilities like PSI (Paul Scherrer Institute) to extract the $^9$Be charge radius from recent and future X-ray spectroscopy measurements with high precision, limited primarily by nuclear structure uncertainties rather than QED theory.
• Methodological Benchmark: By systematically comparing term-by-term, the paper identifies exactly where perturbative expansions break down (e.g., the need for all-order treatment of finite size for higher $Z$) and validates the numerical stability of all-order calculations.
• New Physics Potential: High-precision binding energies in muonic atoms are sensitive to physics beyond the Standard Model (e.g., new bosons). Establishing a reliable QED baseline for $Z=4$ is crucial for isolating potential new physics signals.

In conclusion, the paper establishes that for intermediate $Z$ systems like muonic Beryllium, the choice between perturbative and all-order methods is less critical than previously thought, provided that recoil and higher-order radiative corrections are handled with care. The resulting theoretical predictions are precise enough to support next-generation nuclear structure experiments.