Self-energy corrections to the ionization energies in sodium-like ions: comparison of the \textit{ab initio} QED and model-QED-operator approaches

This paper presents a comparative study of self-energy corrections to ionization energies in sodium-like ions using rigorous bound-state QED and model-QED-operator approaches, demonstrating their good agreement and validating the efficiency of the latter method for many-electron systems.

The Gist

Imagine you are trying to weigh a single grain of sand (an electron) sitting on a giant, heavy boulder (an atomic nucleus). In the world of physics, this isn't just about gravity; it's about the intense electromagnetic forces that bind them together.

This paper is a report from a team of physicists who are trying to measure the "weight" (specifically, the ionization energy) of that grain of sand when it's part of a specific type of atom called a sodium-like ion. These are atoms that have been stripped of most of their electrons, leaving just one "valence" electron hanging out on the outside, surrounded by a tight core of inner electrons.

Here is the breakdown of their work, explained through simple analogies:

1. The Problem: The "Self-Reflection" Tax

In the quantum world, an electron isn't just a tiny ball; it's a fuzzy cloud of probability. Because it's charged, it constantly interacts with the vacuum of space around it, creating and destroying tiny "ghost" particles. This interaction changes the electron's energy. Physicists call this the Self-Energy.

Think of it like this: Imagine you are walking through a crowded room. As you move, you bump into people, they bump into you, and the crowd shifts around you. This interaction makes it harder for you to move than if you were walking alone in an empty hallway. That extra effort is the "self-energy."

For heavy atoms (like the ones in this study, with nuclear charges $Z=30, 50, 70, 92$), this "crowd" effect is so intense that standard math breaks down. The forces are so strong that you can't just add up small corrections; you have to calculate the whole interaction from scratch.

2. The Two Methods: The "Master Chef" vs. The "Smart Shortcut"

The authors wanted to calculate this energy shift for sodium-like ions. They used two different approaches to see if they got the same answer.

Method A: The "Rigorous QED" (The Master Chef)

This is the ab initio (from the beginning) approach.

• The Analogy: Imagine a Master Chef trying to bake a perfect cake. They don't use a recipe; they start with raw flour, eggs, and sugar, measuring every single molecule and accounting for the humidity, the oven's heat distribution, and the chemical reaction of every ingredient.
• In Physics: This method uses the full, complex equations of Quantum Electrodynamics (QED). It calculates every possible interaction between the electron, the nucleus, and the electromagnetic field, order by order. It is incredibly accurate but also incredibly slow and computationally expensive. It's like trying to simulate every single air molecule in a room to predict the wind.

Method B: The "Model-QED Operator" (The Smart Shortcut)

This is the approximate approach.

• The Analogy: Imagine a Smart Shortcut Chef. They know the physics of baking, but instead of measuring every molecule, they use a "magic spice blend" (a mathematical operator) that simulates the effect of the complex interactions. They add this spice to a standard recipe, and it gives a result that is 99% as good as the Master Chef's, but in a fraction of the time.
• In Physics: This method replaces the complex, messy calculations with a pre-calculated "operator" (a mathematical tool) that approximates the self-energy. It's much faster and easier to use for complex atoms with many electrons.

3. The Experiment: Testing the Shortcut

The team calculated the energy for sodium-like ions with different nuclear charges ($Z=30$ to $Z=92$). They used various "screening potentials" (different ways of modeling how the inner electrons shield the outer one) to see how the results changed.

The Results:

• For Heavy Atoms ($Z=92$): Both the Master Chef and the Shortcut Chef agreed almost perfectly. The "magic spice" worked great.
• For Lighter Atoms ($Z=30$): The results started to drift apart slightly. Why? Because in lighter atoms, the electrons interact with each other in complex, chaotic ways (correlation effects) that are hard to predict with a simple shortcut. The "Master Chef" method struggled a bit too, showing that the math is tricky here.

4. The "Secret Sauce": Mixing the Methods

The authors found that the best way to get a reliable answer for all atoms was to combine the Shortcut with a powerful technique called Configuration Interaction (CI).

• The Analogy: Instead of just adding the "magic spice" to a simple recipe, they added it to a complex, multi-layered cake that already accounted for how the ingredients interacted with each other.
• The Result: When they did this, the "Shortcut" method became incredibly robust. It stopped caring about which specific "screening potential" they used and gave consistent, accurate results that matched the "Master Chef" method.

The Big Takeaway

This paper is a "quality control" report. The authors wanted to prove that the Model-QED Operator (the shortcut) is a reliable tool for scientists.

They demonstrated that:

1. It works: The shortcut gives results very close to the super-accurate, super-slow method.
2. It's efficient: It allows scientists to study complex, multi-electron systems without needing a supercomputer to run for a year.
3. It needs a partner: To get the best results, especially for lighter atoms, you should combine this shortcut with advanced correlation methods (like CI).

In short: They proved that you don't always need to be a Master Chef to bake a perfect cake. If you have the right "magic spice" and a good understanding of how your ingredients mix, you can get a delicious result much faster. This is great news for physicists who need to calculate energy levels for many different types of atoms quickly and accurately.

Technical Summary

1. Problem Statement

The paper addresses the challenge of accurately calculating self-energy (SE) corrections to the ionization energies of sodium-like ions (electronic configuration $1s^2 2s^2 2p^6 v$, where $v$ is a valence electron in the $3s$, $3p_{1/2}$, or $3p_{3/2}$ state).

• Context: For highly charged ions (HCIs), the expansion in the electron-nucleus coupling parameter $\alpha Z$ breaks down, necessitating non-perturbative, all-order calculations in $\alpha Z$.
• Difficulty: While rigorous ab initio Quantum Electrodynamics (QED) methods exist, they are computationally complex and become extremely difficult for systems with multiple valence electrons or neutral systems where perturbative treatments of interelectronic interactions may fail.
• Gap: Approximate methods, such as the model-QED-operator approach, are widely used for many-electron systems but lack sufficient validation against rigorous ab initio calculations for sodium-like configurations. The authors aim to fill this gap by comparing these two independent methodologies.

2. Methodology

The study employs two distinct theoretical frameworks to calculate SE corrections for ions with nuclear charge numbers $Z = 30, 50, 70,$ and $92$.

A. Rigorous Ab Initio QED Approach

• Formalism: Based on the Furry picture of QED, where the unperturbed system is described by the Dirac equation with an effective potential ($V_{eff} = V_{nucl} + V_{scr}$).
• Perturbation Theory: Evaluates SE corrections up to the second order in the QED perturbation series.
  • First Order: Includes the one-electron self-energy diagram (Fig. 1) and mass counterterms.
  • Second Order: Includes screened self-energy corrections (Fig. 2) to account for interelectronic interactions.
• Screening Potentials: To accelerate convergence, various screening potentials are used in the zeroth-order approximation:
  • $x_a$-potentials: Based on core electron charge density with parameters $x_a = 0, 1/3, 2/3, 1$ (representing Core-Hartree, Kohn-Sham, and Dirac-Slater potentials).
  • SC Potential: A specific Kohn-Sham potential using total charge density (core + valence), adopted for comparison with previous literature (Sapirstein and Cheng).
• Technical Implementation: Uses potential-expansion methods to handle ultraviolet divergences and improved partial-wave expansions (Ref. [21, 22]) to ensure numerical convergence.

B. Model-QED-Operator Approach

• Framework: Utilizes the QEDMOD package based on the Two-Time Green's Function (TTGF) method.
• Operator Structure: The SE operator ($h_{SE}$) is decomposed into:
  • Semilocal part ($h_{SE}^{sl}$): Depends on angular momentum quantum number $\kappa$, modeled as an exponential function.
  • Nonlocal part ($h_{SE}^{nl}$): A separable operator constructed to reproduce exact hydrogenic matrix elements.
• Implementation Variants:
  1. QEDMOD(av): Calculates the expectation value of $h_{SE}$ using valence wavefunctions derived from the Dirac equation with screening potentials.
  2. CI-QEDMOD: Incorporates $h_{SE}$ directly into the Dirac-Coulomb-Breit (DCB) Hamiltonian within a Configuration-Interaction (CI) framework. The SE contribution is the difference in ionization energies with and without the operator.
  3. CI-QEDMOD(scf): A self-consistent variant where the model-QED operator is included in the one-electron Hamiltonian ($h_\Lambda$) used to construct the CI basis, accounting for single-particle excitations into the negative-energy continuum.

3. Key Contributions

• Systematic Comparison: Provides the first detailed, state-by-state comparison between rigorous ab initio QED (up to 2nd order) and the model-QED-operator method for sodium-like ions across a wide range of nuclear charges ($Z=30$ to $92$).
• Validation of Approximations: Demonstrates that the model-QED-operator approach, when combined with the CI method, yields results highly consistent with rigorous QED, validating its use for complex many-electron systems.
• Analysis of Correlation Effects: Investigates how the choice of the initial screening potential affects the results. It highlights that while first-order results vary significantly with the potential choice, the inclusion of second-order corrections (in rigorous QED) or the CI treatment (in the model approach) stabilizes the results.

4. Key Results

• Agreement: There is excellent agreement between the rigorous ab initio results and the CI-QEDMOD results for all studied states ($3s, 3p_{1/2}, 3p_{3/2}$) and nuclear charges.
• Dependence on Initial Approximation:
  • For low Z (e.g., $Z=30$), the first-order SE corrections calculated via the rigorous method show significant scatter depending on the chosen screening potential. This indicates that correlation effects are strong and not fully captured by first-order perturbation theory.
  • The CI-QEDMOD results remain stable and consistent regardless of the screening potential used, as the CI method treats correlation effects to all orders in $1/Z$.
• Second-Order Effects: The inclusion of second-order screened self-energy corrections (SE+ScrSE) significantly alters the total values, particularly for the $3p$ states and lower $Z$ values, emphasizing the necessity of accounting for electron-electron correlations in QED calculations.
• Specific Discrepancy: A discrepancy was noted between the authors' rigorous results and those of Sapirstein and Cheng (Ref. [45]) for the SC potential in the second-order corrections, though the origin remains unclear. However, the authors' results are cross-validated by the independent CI-QEDMOD method.
• Self-Consistency: The CI-QEDMOD(scf) approach yields results nearly identical to the standard CI-QEDMOD, suggesting that for sodium-like ions, the self-consistent modification of the projector operators does not drastically change the outcome, though this was not predictable a priori.

5. Significance

• Methodological Validation: The study confirms that the model-QED-operator approach is a reliable, accurate, and efficient tool for evaluating QED effects in many-electron systems where full ab initio QED calculations are prohibitively difficult.
• Handling Correlations: It demonstrates that combining the model-QED operator with the Configuration-Interaction (CI) method effectively mitigates the dependence on the choice of the zeroth-order approximation, providing robust estimates even when $Z$ and the number of electrons $N$ are comparable (i.e., in less highly charged ions).
• Benchmarking: The provided numerical data for $Z=30, 50, 70, 92$ serves as a critical benchmark for future theoretical developments and experimental spectroscopy of sodium-like ions.
• Practical Application: The findings support the use of the QEDMOD package for high-precision atomic structure calculations in complex ions, facilitating research in atomic physics, plasma physics, and tests of fundamental physics.