Rci-Q: an improved QED correction model for the GRASP2018 package

The Rci-Q package extends the GRASP2018 suite by implementing an improved QED correction model that incorporates the Flambaum-Ginges radiative potential with new fitting prefactors, finite nucleus size corrections to self-energy, and the Wichmann-Kroll vacuum polarization potential.

The Gist

Imagine you are trying to predict the exact path of a tiny, fast-moving marble (an electron) as it zooms around a heavy bowling ball (an atomic nucleus). In the world of heavy atoms, this isn't just a simple game of pool; it's a high-stakes dance where the rules of physics get a little weird.

This paper introduces a new tool called Rci-Q, which is an upgrade to a famous computer program called Grasp2018. Think of Grasp2018 as a very sophisticated GPS for atoms. It helps scientists calculate where electrons are and how much energy they have. But for heavy atoms (like gold, uranium, or lead), the old GPS had a few blind spots. It was missing some subtle, invisible forces that only show up when things get really heavy and fast.

Here is the breakdown of what this paper does, using some everyday analogies:

1. The Problem: The "Ghost" Forces

In the world of Quantum Electrodynamics (QED), electrons don't just sit still. They are constantly interacting with "ghosts."

• The Self-Energy (The Echo): Imagine an electron shouting in a canyon. The sound bounces back, and the echo pushes the electron slightly. This is the electron interacting with its own radiation field. Calculating this "echo" is incredibly hard. The old Grasp2018 program used a rough guess for this, which worked okay for light atoms but got messy for heavy ones.
• The Vacuum Polarization (The Fog): Imagine the space around the nucleus isn't empty, but filled with a fog of virtual particles popping in and out of existence. This fog changes how the electron feels the pull of the nucleus. The old program handled the thick part of the fog well, but missed the thinner, more complex layers.

2. The Solution: A Better Map (Rci-Q)

The author, Karol Kozioł, has built a new module called Rci-Q that acts like a high-definition lens for the GPS. It doesn't just guess; it uses a specific, highly accurate mathematical method (the Flambaum–Ginges radiative potential) to calculate these "ghost" forces on the fly.

Think of it like upgrading from a paper map to a real-time satellite navigation system.

• The "Fitting" Magic: The author didn't just invent new math; they took the best existing math and "tuned" it. Imagine you have a recipe for a cake that works great for small cakes but fails for giant wedding cakes. The author ran thousands of simulations to create a new set of "adjustment knobs" (fitting coefficients) that make the recipe perfect for any size cake, from tiny hydrogen atoms to massive uranium atoms. These new knobs are listed in the tables in the paper.

3. New Features Added to the Toolkit

The Rci-Q upgrade adds three specific improvements:

• The "Finite Nucleus" Correction: In the old model, the nucleus was treated like a perfect, tiny point. In reality, the nucleus has a size (it's a fuzzy ball). For heavy atoms, this fuzziness matters. Rci-Q accounts for the fact that the nucleus takes up space, refining the calculation.
• The "Wichmann-Kroll" Layer: This is a deeper, more complex layer of the "fog" (vacuum polarization). The old program ignored this deep layer. Rci-Q digs it up and includes it, making the map much more accurate for heavy elements.
• On-the-Fly Calculation: Previously, if you wanted these corrections, you had to run the program, stop, do some extra math manually, and run it again. Rci-Q does all this automatically while it's running, saving time and reducing human error.

4. Does It Work? (The Test Drive)

The author tested this new GPS against real-world data and other super-computers:

• Hydrogen-like Atoms: It matched the known "gold standard" numbers almost perfectly (within a tiny fraction of a percent).
• Helium-like Atoms: When looking at transitions in heavy helium atoms, the new model (Rci-Q) got much closer to the experimental measurements than the old model did. For the heaviest atoms tested (Uranium), the difference was huge—about 10 electron-volts, which is a massive gap in atomic physics.
• Fluorine-like Ions: It successfully predicted the "fine splitting" (the tiny gap between energy levels) in heavy ions, matching experimental results better than previous attempts.

5. The Cost: A Little Slower, Much Better

There is a small trade-off. Because the computer is doing more complex math, the calculations take about 20% longer.

• Analogy: It's like driving a car. The old way was driving a sports car on a straight road (fast, but you might miss a turn). The new way is driving a heavy-duty truck with a full GPS and terrain mapping (slightly slower, but you are guaranteed to get exactly where you need to go without crashing).

The Bottom Line

This paper presents a vital upgrade for scientists studying heavy atoms. By refining how we calculate the subtle interactions between electrons and the vacuum of space, Rci-Q allows us to see the atomic world with much sharper clarity. It ensures that when we study heavy elements (which are crucial for nuclear physics, medical imaging, and understanding the universe), our theoretical models actually match what we see in the lab.

Technical Summary

1. Problem Statement

Quantum Electrodynamics (QED) corrections are critical for accurately calculating energy levels in high-$Z$ (high atomic number) highly charged ions. While the leading-order vacuum polarization (Uehling potential) is standard in atomic structure codes like Grasp2018, the self-energy (SE) correction is significantly more complex to implement.

• Limitations of Current Models: The default SE model in the Grasp2018 rci program is simple but lacks accuracy, particularly for inner shells and high-$Z$ atoms.
• Workflow Disruption: Existing alternative methods (e.g., rci4, Qedmod) provide better accuracy but require separate, non-integrated steps, making them cumbersome for routine calculations.
• Need: There is a requirement for an improved, on-the-fly SE correction model integrated directly into the rci workflow that maintains computational efficiency while enhancing precision for heavy elements.

2. Methodology

The paper introduces rci-q, a modified version of the Grasp2018 rci program that implements the Flambaum–Ginges radiative potential method. This approach approximates the leading self-energy correction as a local potential, $V_{SE}(r)$, which is computationally efficient.

A. Theoretical Framework

The self-energy potential is decomposed into three components:
$$V_{SE}(r) = V_{el}(r) + V_{low}(r) + V_{mag}(r)$$

1. Electric Part ($V_{el}$): The high-frequency contribution, dominant for $l=0$ orbitals.
2. Low-Frequency Part ($V_{low}$): Dominant for $l > 0$ orbitals.
3. Magnetic Part ($V_{mag}$): Arising from the electron's anomalous magnetic moment interaction.

The potential is parameterized using fitting coefficients ($A$ and $B$) derived from reference data for hydrogen-like ions.

B. Key Improvements and Implementations

1. Refitted Prefactors ($A$ and $B$):
   • The original Flambaum–Ginges coefficients were optimized for $n=5$ subshells and performed poorly for inner shells.
   • New Fits: The author presents new fitting coefficients for $A(Z, n, l)$ and $B(Z, n, \kappa)$ derived from high-precision reference data (Mohr, Yerokhin, etc.) for $Z$ ranges from 5 to 110.
   • Polynomial Fits:
     • $A(Z)$ for $l=0$ is fitted using 4th-order polynomials separately for $Z \ge 20$ and $Z < 20$.
     • $B(Z)$ for $l=1, 2, 3$ is fitted using 3rd and 4th-order polynomials in terms of $Z\alpha$.
2. Finite Nucleus Size (FNS) Correction:
   • A new correction term for the self-energy arising from the finite size of the nucleus is implemented.
   • It is parameterized using exponential fitting functions based on data from Refs. [22, 23], covering subshells $s, p, d$ for various $Z$ ranges.
3. Wichmann-Kroll (WK) Vacuum Polarization:
   • The implementation includes the Wichmann-Kroll part of the vacuum polarization potential (higher-order terms in $\alpha(Z\alpha)$).
   • This is based on the algorithm by Fainshtein et al., replacing the older Källén–Sabry implementation in the original code.
   • Note: The implementation assumes a point-like nucleus for the WK term, leading to a systematic overestimation (3–10%) compared to literature using uniform charge distributions, attributed to the lack of FNS correction in the WK algorithm itself.

3. Key Contributions

• Integration: Seamless integration of advanced QED corrections (SE and WK) into the standard Grasp2018 rci workflow without requiring external post-processing steps.
• Parameterization: Provision of a comprehensive set of new fitting coefficients (Tables I–IV) for $A(Z)$ and $B(Z)$ prefactors, valid for a wide range of atomic numbers ($Z$) and quantum numbers ($n, l, \kappa$).
• Code Availability: The rci-q package is provided as a drop-in replacement for the standard rci executable, with a computational overhead of only ~20%.
• FNS Implementation: First inclusion of the finite nucleus size correction specifically for the self-energy potential within this radiative potential framework.

4. Results and Benchmarking

The paper validates the rci-q model through extensive benchmarking against reference data and experimental results:

• Hydrogen-like Systems:
  • Comparison of calculated self-energy shifts ($F(Z\alpha)$) against reference values shows excellent agreement.
  • Differences are within 0.03% for $s$-shells, 0.5% for $p$-shells, and up to 3% for $f$-shells at $Z=50$.
• Helium-like Systems:
  • Calculations for the $1s2p \ ^1P_1 \to 1s^2 \ ^1S_0$ transition (K$\alpha_1$ line) were performed for ions from $Z=26$ to $Z=92$.
  • For high-$Z$ (e.g., Uranium, $Z=92$), the difference between the original Grasp2018 model and rci-q reaches ~10 eV.
  • The rci-q results align significantly better with the weighted mean of experimental data than the original model.
• Fluorine-like Systems:
  • The $2p_{3/2} - 2p_{1/2}$ fine splitting was tested for $Z=42$ (Mo) and $Z=92$ (U).
  • For $Z=92$, the rci-q result (2.47 eV) matches closely with high-precision theoretical models (Volotka et al., Shabaev et al.) and experimental deductions, outperforming the original Grasp2k model (1.33 eV).

5. Significance

• High-$Z$ Accuracy: The rci-q package significantly improves the reliability of atomic structure calculations for heavy elements where QED effects dominate. This is crucial for testing fundamental physics in strong fields.
• Efficiency: By embedding the Flambaum–Ginges potential directly into the radial grid integration, the method avoids the computational cost of separate perturbative steps, making high-precision QED corrections accessible for large-scale Multi-Configuration Dirac–Hartree–Fock (MCDHF) and Configuration Interaction (CI) calculations.
• Standardization: The new fitting coefficients and the rci-q implementation provide a standardized, robust tool for the atomic physics community, reducing the need for ad-hoc corrections or external packages.

In conclusion, rci-q represents a major step forward in the practical application of QED corrections in atomic structure codes, offering a balance of high accuracy and computational feasibility for heavy atoms.