Rydberg states of muonic helium in quantum electrodynamics

This paper employs the variational method with Gaussian wave functions to analytically calculate the energy levels and relativistic corrections of Rydberg states in muonic helium, providing a theoretical basis for future experimental studies.

The Gist

Imagine a tiny, exotic solar system. In our normal world, a helium atom has a heavy nucleus in the center with two light electrons buzzing around it. But in this paper, the authors are studying a strange, temporary version of this atom called muonic helium.

Here, one of the electrons has been swapped out for a muon. A muon is like a "heavy electron"—it has the same charge but is about 200 times heavier. Because it's so heavy, it doesn't just orbit the nucleus; it crashes deep into the inner layers, usually getting very close to the center.

The "Rydberg" Twist: A High-Flying Dancer

Usually, when a muon gets captured by a helium atom, it falls down to the lowest, most stable orbit (the ground state) very quickly. However, the authors are interested in a special, rare scenario where the muon gets stuck in a Rydberg state.

Think of a Rydberg state like a dancer spinning on a very high stage, far away from the center. In this specific study, the muon is in a high-energy orbit (around level 14) where it is roughly the same distance from the nucleus as the remaining electron. It's as if the heavy muon and the light electron are holding hands and dancing in a wide circle around the nucleus, keeping an equal distance from the center.

The Problem: Calculating the Dance

Calculating the energy of this three-part dance (Nucleus + Muon + Electron) is incredibly difficult. It's like trying to predict the exact path of three people holding hands while running on a trampoline, where everyone is pulling on everyone else.

The authors used a mathematical tool called the Variational Method. Imagine you are trying to guess the shape of a complex, wobbly jelly. Instead of trying to solve the exact physics of every molecule in the jelly, you build a model out of smooth, simple shapes (in this case, Gaussian curves, which look like perfect bell curves or soft hills). You stack these smooth hills together to approximate the wobbly jelly.

By adjusting the size and shape of these "hills," they found the best mathematical fit for the energy of this exotic atom.

Adding the "Fine Print"

Once they had the basic shape of the energy levels, they had to add the "fine print" corrections. In the quantum world, things aren't perfectly simple. They added three specific corrections to their calculation:

1. Relativism: Because the particles are moving fast, they have to account for Einstein's theory of relativity (like a speedometer that changes as you get closer to the speed of light).
2. Vacuum Polarization: In quantum physics, empty space isn't truly empty; it's filled with popping in and out of "virtual" particles. The authors calculated how this "quantum foam" slightly pushes or pulls on the muon and electron.
3. Contact Interactions: This accounts for what happens when the particles get extremely close to each other, almost touching.

The Results: A New Map

The paper provides a detailed map of energy levels for these high-flying muonic helium atoms. They calculated exactly how much energy is needed to move the muon between these high orbits.

Why does this matter?

• Precision: These calculations are so precise that they can be used by experimentalists to check their measurements. If scientists shine a laser at these atoms and see a specific color of light (energy), they can compare it to this paper's map to see if their math matches reality.
• Solving Mysteries: The introduction mentions a "proton radius puzzle" (a disagreement between how big we think the proton is based on different experiments). While this paper focuses on helium, the methods used here help refine our understanding of fundamental constants, which helps solve those bigger puzzles.
• Measuring Mass: The authors note that measuring the transition frequencies (the "notes" the atom sings) in these Rydberg states could help scientists determine the mass of the muon with extreme accuracy.

The Bottom Line

This paper is a theoretical blueprint. The authors didn't build the atom; they built the mathematical model for it. They showed us exactly what the energy levels should look like for a muon and an electron dancing in a wide circle around a helium nucleus. This blueprint is now ready for experimental physicists to use as a reference to test their own high-precision laser experiments.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the theoretical investigation of muonic helium atoms ($\mu^- - e^- - \text{He}$), specifically focusing on Rydberg states where the muon is in a highly excited state with large principal ($n$) and orbital ($l$) quantum numbers ($n \sim l + 1 \sim 14$).

• Scientific Context: Rydberg states are crucial for refining fundamental constants (e.g., the Rydberg constant) and resolving discrepancies like the "proton radius puzzle." In these states, the electron and muon are approximately equidistant from the nucleus, minimizing the influence of nuclear structure effects on the energy spectrum.
• The Gap: While low-lying energy levels of muonic helium have been studied extensively, there is a lack of precise theoretical data for high-lying Rydberg states where the muon and electron coexist at similar distances from the nucleus. These states are relevant for interpreting X-ray spectroscopy data from muon capture experiments and for potential future laser spectroscopy to determine the muon mass with higher precision.
• Challenge: The three-body nature of the system (nucleus, muon, electron) combined with the large orbital angular momentum of the muon makes standard analytical perturbation theory difficult to apply.

2. Methodology

The authors employ a rigorous variational method within the framework of Quantum Electrodynamics (QED) to calculate energy levels and structural properties.

• Coordinate System: The system is described using Jacobi coordinates ($\rho, \lambda$), which separate the relative motion of the particles:
  • $\rho$: Vector between the muon and the nucleus.
  • $\lambda$: Vector between the electron and the center of mass of the nucleus-muon pair.
• Wave Function Ansatz:
  • The trial wave functions are constructed as a superposition of Gaussian exponentials multiplied by bipolar spherical harmonics to account for the angular dependence.
  • The radial part is modeled as $e^{-\frac{1}{2}(A_{11}\rho^2 + 2A_{12}\rho\lambda + A_{22}\lambda^2)}$, where $A_{ij}$ are nonlinear variational parameters optimized for the ground and excited states.
  • The angular part uses spherical harmonics $Y_{lm}(\theta_\rho, \phi_\rho)$ to describe the muon's orbital motion, while the electron is assumed to be in the ground state ($1S$).
• Hamiltonian and Corrections:
  • Non-relativistic Hamiltonian ($\hat{H}_0$): Includes kinetic energy terms and pairwise Coulomb interactions. Matrix elements are calculated analytically.
  • Relativistic Corrections: Included via the Breit potential, specifically the $p^4$ term (mass-velocity correction) for the electron and muon.
  • Vacuum Polarization: Modeled using a $\delta$-function approximation (Uehling potential) suitable for the short Compton wavelengths of the particles in these states.
  • Contact Interactions: Included via $\delta$-function terms representing short-range interactions.
• Computational Approach: The Schrödinger equation is reduced to a generalized matrix eigenvalue problem ($H \cdot C = E B \cdot C$). All matrix elements for kinetic, potential, and correction terms are derived analytically and computed using a custom program (previously validated for pion and kaonic helium).

3. Key Contributions

• Extension to High-$n$ States: The study extends variational calculations to muonic helium states with principal quantum numbers up to $n=14$ and orbital angular momenta $l=10, 11, 12, 13$.
• Analytical Precision: The authors provide closed-form analytical expressions for all matrix elements, including complex relativistic and vacuum polarization corrections, ensuring high numerical precision.
• Structural Analysis: The paper calculates the mean-square distances between particles and radial distribution densities ($W(\rho)$, $W(\lambda)$, $W(\rho, \lambda)$), offering a detailed picture of the spatial configuration of the three-body system in Rydberg states.
• Benchmarking: The results serve as a necessary theoretical benchmark for experimentalists studying muon cascades and X-ray emissions in helium targets.

4. Results

• Energy Levels: Table I presents the calculated energy levels for both $^3\text{He}$ and $^4\text{He}$ isotopes.
  • The non-relativistic energies ($E_{nr}$) are provided in electron atomic units (e.a.u.).
  • Specific corrections are quantified:
    • Relativistic correction ($-\frac{\alpha^2}{8}p_e^4$): Ranges from $\approx -0.0002$ to $-0.0003$ e.a.u.
    • Vacuum polarization ($\Delta U_{vp}$): Small negative contributions ($\approx -3 \times 10^{-7}$ e.a.u.).
    • Contact interaction ($\Delta U_{cont}$): Small positive contributions ($\approx 3 \times 10^{-6}$ e.a.u.).
• Spatial Structure:
  • For the $(14, 13)$ state in $^4\text{He}$, the root-mean-square distances are calculated as:
    • Nucleus-Muon ($\delta_{12}$): $1.34$a.u. ($0.71 \times 10^{-8}$ cm)
    • Nucleus-Electron ($\delta_{13}$): $0.47$a.u. ($0.25 \times 10^{-8}$ cm)
    • Muon-Electron ($\delta_{23}$): $1.08$a.u. ($0.57 \times 10^{-8}$ cm)
  • Key Finding: The radial distribution densities confirm that in these Rydberg states, the muon and electron are located at comparable distances from the nucleus, validating the "equidistant" assumption crucial for minimizing nuclear structure effects.
• Isotope Shift: The calculations show distinct energy differences between $^3\text{He}$ and $^4\text{He}$ isotopes, which are sensitive to nuclear mass effects.

5. Significance and Implications

• Muon Mass Determination: The calculated transition frequencies between Rydberg states are proposed as a new avenue for determining the muon mass with high accuracy, potentially improving upon current measurements by orders of magnitude (similar to recent successes with pionic helium).
• Resolving Fundamental Puzzles: By studying states where nuclear structure effects are suppressed, these calculations contribute to the broader effort of resolving the "proton radius puzzle" and refining the Rydberg constant.
• Experimental Guidance: The results provide essential data for interpreting X-ray spectra from muon capture experiments. The predicted energies help identify specific transitions in the cascade process, aiding in the analysis of initial muon capture probabilities and state populations.
• Methodological Advancement: The successful application of the Gaussian variational method with Jacobi coordinates to high-$l$ Rydberg states in three-body systems demonstrates the robustness of this approach for exotic atoms (pionic, kaonic, and muonic helium).

In summary, this paper provides a high-precision theoretical framework for the Rydberg states of muonic helium, bridging the gap between non-relativistic quantum mechanics and QED corrections, and offering critical data for future high-precision spectroscopy experiments aimed at testing the Standard Model.