Bound-electron self-energy calculations in Feynman and Coulomb gauges: detailed analysis

This paper presents a comprehensive comparative analysis of partial-wave expansion convergence for bound-electron self-energy calculations in Feynman and Coulomb gauges to address accuracy limitations in Lamb shift predictions for highly charged ions, while also discussing methods to improve convergence.

The Gist

Imagine you are trying to measure the weight of a single, tiny feather floating in a hurricane. That's essentially what physicists do when they calculate the energy of an electron orbiting a heavy atomic nucleus. The electron isn't just sitting there; it's constantly interacting with the "wind" of quantum fields around it. This interaction changes the electron's energy slightly, a phenomenon known as the Lamb shift.

This paper is a detailed "tuning guide" for the mathematical tools used to measure that tiny energy change. The authors, a team of physicists, are asking a very specific question: Which mathematical "lens" gives us the clearest, most accurate picture of this interaction?

Here is a breakdown of their work using everyday analogies:

1. The Problem: The "Static" in the Signal

In the quantum world, an electron is constantly emitting and re-absorbing tiny packets of light (photons). This creates a "self-energy" effect. To calculate this, physicists use a method called Partial-Wave Expansion (PWE).

Think of the electron's energy as a complex musical chord. To understand the chord, you have to break it down into individual notes (the "partial waves").

• The Issue: In many cases, the lower notes are loud and easy to hear, but the higher notes are faint and get lost in the noise. If you stop listening after the 20th note, your calculation of the chord is slightly wrong. To get the perfect pitch, you need to listen to thousands of notes.
• The Bottleneck: Calculating thousands of notes takes a massive amount of computer power. If the notes don't "settle down" (converge) quickly, the calculation becomes imprecise.

2. The Two Lenses: Feynman vs. Coulomb

Physicists can look at this problem through two different mathematical "lenses" (or gauges): the Feynman gauge and the Coulomb gauge.

• The Feynman Gauge (The "All-Encompassing" Lens): This is the standard, go-to lens used for decades. It treats all forces symmetrically. However, the authors found that when using this lens, the "higher notes" (the partial waves) are very loud and chaotic. You have to listen to many of them to get a clear signal. It's like trying to hear a whisper in a crowded stadium; the background noise is huge.
• The Coulomb Gauge (The "Focused" Lens): This lens focuses specifically on the electric force first. The authors discovered that in this view, the "background noise" is much quieter. The higher notes are naturally smaller and fade away faster. It's like moving to a quiet room to hear that whisper.

The Surprise: While the Coulomb gauge produces much smaller numbers (which sounds good), the authors found that the rate at which the numbers settle down isn't always faster than the Feynman gauge. It's a trade-off. Sometimes the quiet room is better, but sometimes you still need a good microphone.

3. The Magic Tricks: Acceleration Schemes

Since listening to thousands of notes is slow, the authors tested two "magic tricks" (acceleration schemes) to speed up the process. These tricks are like noise-canceling headphones for the calculation.

• The "Two-Potential" Trick: This method tries to predict the "loud" part of the noise mathematically and subtracts it out before you start counting the notes. You only have to calculate the tiny, quiet remainder.
• The "Sapirstein-Cheng" (SC) Trick: This is a clever approximation. Instead of calculating the exact "loud" part (which is hard), it calculates a very close guess that is easy to do. It turns out this guess is so good that it works almost as well as the exact calculation but is much faster.

The Verdict: The authors found that the SC trick combined with the Coulomb gauge is the "Golden Combo." It's like putting on noise-canceling headphones in a quiet room. You get the cleanest signal with the least amount of effort.

4. The Tools: Two Ways to Build the Map

To do these calculations, the team used two different construction methods to build their mathematical map of the electron:

1. The Green's Function (GF) Method: This is like using a high-resolution satellite map. It's incredibly precise but can be computationally heavy and tricky to handle because the map has "jumps" or discontinuities.
2. The Finite-Basis Set (FBS) Method: This is like building a map out of Lego bricks. It's flexible and easy to build, but you have to keep adding more bricks to get the edges smooth. If you stop too early, your map is a bit blocky.

The authors showed that while the "Lego" method is great for flexibility, the "Satellite" method (Green's Function) generally gives the most precise results for this specific problem.

The Big Picture Conclusion

Why does this matter?

• Precision is King: In modern physics, we are testing the fundamental laws of the universe. If our calculations of the electron's energy are off by a tiny fraction, we might think we've discovered a new particle or a new force when we've actually just made a math error.
• The Recommendation: The paper concludes that for the most accurate results, especially for lighter atoms, physicists should use the Coulomb gauge combined with the Sapirstein-Cheng acceleration scheme.

In a nutshell: The authors didn't just find a new way to calculate an electron's energy; they optimized the entire process. They figured out which mathematical lens to use and which noise-canceling trick to apply to get the most precise measurement possible, ensuring that our understanding of the quantum world remains rock-solid.

Technical Summary

1. Problem Statement

The accurate theoretical description of atomic systems, particularly highly charged ions (HCIs), requires precise calculations of Quantum Electrodynamics (QED) effects. The one-electron self-energy (SE) correction is the leading contribution to the Lamb shift and is critical for determining atomic energy levels, $g$-factors, hyperfine structure, and transition amplitudes.

The primary challenge in calculating the SE correction lies in the many-potential term (the part of the SE involving two or more interactions with the nuclear binding potential).

• Convergence Issues: Conventional approaches rely on Partial-Wave (PW) expansions to handle the bound-electron propagator in coordinate space. However, these series often converge slowly, limiting the achievable numerical accuracy.
• Gauge Dependence: While the total SE correction is gauge-invariant, its individual components (zero-, one-, and many-potential terms) are gauge-dependent. Previous studies suggested the many-potential term is smaller in the Coulomb gauge than in the Feynman gauge, but it was unclear if this translated to better convergence rates for the PW expansion.
• Acceleration Schemes: Various schemes exist to accelerate PW convergence (e.g., subtracting the two-potential term), but their comparative effectiveness across different gauges and nuclear models (point vs. extended nuclei) required rigorous re-evaluation.

2. Methodology

The authors performed a comprehensive comparative analysis using two distinct numerical approaches and two photon gauges (Feynman and Coulomb).

A. Theoretical Framework

• Formalism: Calculations were performed within the Furry picture of QED, treating the electron non-perturbatively in the nuclear potential $V(r)$.
• Renormalization: UV divergences in the zero- and one-potential terms were handled via potential expansion and dimensional regularization in momentum space. The remaining many-potential term was calculated in coordinate space.
• Gauges: The study compared the Feynman gauge (predominantly used historically) and the Coulomb gauge (known for smaller many-potential contributions).

B. Numerical Approaches

Two state-of-the-art methods were employed to evaluate the coordinate-space integrals:

1. Green's Function (GF) Method: Uses analytical solutions (Whittaker functions) for point nuclei or numerical solutions for extended nuclei. It constructs the propagator using solutions bounded at the origin and infinity. This method avoids basis-set extrapolation but requires careful handling of discontinuities in the Green's function.
2. Finite-Basis-Set (FBS) Method: Uses B-splines within the Dual-Kinetic-Balance (DKB) approach to construct the basis. This allows for flexible handling of arbitrary potentials and non-spherical systems but requires extrapolation to an infinite basis size ($N \to \infty$).

C. Convergence Acceleration Schemes

The paper analyzed three strategies to improve PW convergence:

1. Direct (DIR): Standard calculation without subtraction.
2. Two-Potential (TP) Scheme: Subtracts the exact two-potential term (calculated separately) from the many-potential term.
3. Sapirstein-Cheng (SC) Scheme: Subtracts a "quasi-two-potential" term. This approximation is calculated efficiently in momentum space and serves as a proxy for the exact two-potential term, leaving a remainder with faster convergence.

D. Extrapolation Procedures

• Basis Extrapolation (FBS): Partial sums were extrapolated to infinite basis size ($N \to \infty$) using a non-linear least-squares ansatz.
• Partial-Wave Extrapolation: The remainder of the PW series (truncated at $k_{max}$) was estimated using statistical extrapolation based on the asymptotic behavior $\Delta \epsilon(k) \sim 1/k^p$.

3. Key Contributions

• Comprehensive Gauge Comparison: The authors provided a detailed, independent analysis of PW convergence in both Feynman and Coulomb gauges for both point and extended nuclei, challenging some previous assumptions about the superiority of the Coulomb gauge.
• Validation of Acceleration Schemes: They rigorously tested the TP and SC schemes, demonstrating that the SC scheme provides an excellent approximation to the TP scheme while being computationally more tractable.
• Asymptotic Analysis: The study confirmed that the leading asymptotic behavior of the PW expansion terms follows $\Delta \epsilon(k) \sim 1/k^3$ (i.e., $p=3$) for both gauges, justifying specific extrapolation ansatzes.
• Data Compilation: The paper collects and provides a self-consistent set of formulas and numerical results for SE corrections across a wide range of nuclear charges ($Z=10$ to $Z=100$) and states ($1S, 2S, 2P$).

4. Results

• Convergence Rates:
  • In the Coulomb gauge, the absolute magnitude of the many-potential contribution is significantly smaller (often by orders of magnitude for low $Z$) compared to the Feynman gauge.
  • However, the rate of convergence of the PW series is not inherently faster in the Coulomb gauge for large $Z$. The numerical coefficients for higher-order terms are larger in the Coulomb gauge, which partially offsets the benefit of the smaller absolute value.
  • For low-$Z$ systems, the Coulomb gauge shows slightly better convergence behavior regarding basis-set size.
• Effectiveness of Acceleration Schemes:
  • Both TP and SC schemes significantly reduce the magnitude of the remainder to be extrapolated.
  • The SC scheme is found to be highly effective, providing results nearly identical to the TP scheme but with easier implementation (momentum-space calculation of the subtracted term).
  • Using the SC scheme allows for the extraction of two additional significant digits in the final result.
• Method Comparison:
  • The GF method generally yields higher accuracy for the many-potential term because it does not suffer from basis-set extrapolation errors.
  • The FBS (DKB) method is more flexible for complex systems (non-spherical) but requires careful extrapolation.
• Optimal Strategy: The most accurate results are obtained by combining the Coulomb gauge with the SC acceleration scheme and the GF method (or highly converged FBS). This combination minimizes numerical cancellations and maximizes convergence speed.
• Benchmarking: The calculated SE corrections for $1S_{1/2}$, $2S_{1/2}$, $2P_{1/2}$, and $2P_{3/2}$ states in ions from Argon ($Z=18$) to Uranium ($Z=92$) and beyond show excellent agreement with previous benchmark calculations (e.g., Refs. [26, 69]), with minor deviations attributed to different nuclear models or constants.

5. Significance

• Precision Metrology: The refined methods and data presented are essential for interpreting high-precision spectroscopy experiments in HCIs, which are used to test QED in strong fields and determine fundamental constants (like the fine-structure constant $\alpha$).
• Methodological Guidance: The paper resolves ambiguities regarding gauge choice and convergence acceleration. It establishes that while the Coulomb gauge offers smaller absolute errors in the many-potential term, the SC acceleration scheme is the critical factor for achieving high precision, regardless of the gauge.
• Future Applications: The authors note that the SC scheme and the derived formulas are generalizable to more complex diagrams, such as two-electron SE and two-loop SE corrections, facilitating future high-order QED calculations.
• Resource Efficiency: By identifying the most efficient combination of gauge and acceleration scheme, the work reduces the computational cost required to achieve specific accuracy targets in atomic physics simulations.

In conclusion, this work provides a definitive guide for calculating bound-electron self-energy corrections, demonstrating that the combination of the Coulomb gauge and the Sapirstein-Cheng convergence-acceleration scheme offers the most robust path to high-precision results for hydrogen-like ions.