🔬 atomic physics

Accurate Wichmann-Kroll vacuum polarization calculations with Gaussian basis set

This paper demonstrates that large Gaussian basis sets can accurately calculate all-order Wichmann-Kroll vacuum polarization corrections in highly-charged ions, achieving precision comparable to traditional methods while enabling faster convergence for difference-based quantities and facilitating the calculation of screening corrections in lithium-like systems.

Original authors: V. K. Ivanov, D. A. Glazov, A. V. Volotka

Published 2026-10-22

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: V. K. Ivanov, D. A. Glazov, A. V. Volotka

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the realm of the very small, where atoms are stripped of their outer electrons to become highly charged ions, the laws of physics behave in ways that seem almost alien. To predict the exact energy levels of these ions, scientists must account for the subtle, jittery nature of the vacuum itself. Even empty space is not truly empty; it is a seething foam of virtual particles that constantly pop in and out of existence. When a powerful electric field from a heavy atomic nucleus passes through this foam, it distorts the vacuum, creating a cloud of virtual charges that slightly alters the energy of the real electrons orbiting the nucleus. This phenomenon, known as vacuum polarization, is a crucial piece of the puzzle for understanding the fundamental constants of nature. For decades, calculating the most complex part of this effect required methods that were either incredibly difficult to apply or limited in their precision, leaving a gap between theory and the ultra-precise measurements now possible in modern laboratories.

A team of researchers has now bridged that gap by developing a new way to calculate these vacuum effects with unprecedented accuracy. They focused on a specific, difficult-to-calculate contribution called the Wichmann-Kroll correction, which represents the non-linear, higher-order effects of the vacuum distortion. Instead of using the traditional mathematical tools that had been the standard for years, the team employed a massive collection of mathematical building blocks known as a Gaussian basis set. Imagine trying to paint a perfect circle using only square tiles; the more tiles you use, the smoother the curve becomes. The researchers used up to 500 of these mathematical "tiles" to construct the electron's wave function, allowing them to model the vacuum distortion with a level of detail that had previously been out of reach for this specific type of calculation.

The results of this effort are striking. By pushing the size of their mathematical basis set to 500, the team achieved a level of precision that matches the most trusted, traditional methods used in the field. They found that while the total energy correction takes a long time to settle into a final, stable number as they added more mathematical tiles, the differences between various models of the atomic nucleus settled down much faster. This discovery is significant because it means scientists can now determine how the specific shape and size of a nucleus affects these quantum corrections without needing to run the most computationally expensive simulations for every single variation. It is like realizing that while you need a high-resolution camera to see the entire landscape clearly, you only need a standard lens to see how the lighting changes when you move a single tree.

The team also applied their new method to lithium-like ions, which are atoms with three electrons, to study how the presence of other electrons screens or shields the nuclear charge. They calculated the vacuum polarization screening correction for the ground state of a uranium ion and found their results to be in excellent agreement with previous independent calculations. This confirms that their approach is not only accurate for simple, single-electron systems but also robust enough to handle the added complexity of multi-electron atoms. The researchers demonstrated that by using their large basis set, they could estimate the errors in their calculations and confirm that their results were converging toward the true physical value, providing a reliable new tool for the community.

Ultimately, this work provides a powerful, independent way to verify the most precise predictions in atomic physics. By confirming that their results align with those obtained through different, established techniques, the researchers have validated a method that can be used to explore even more complex quantum diagrams in the future. Their findings ensure that as experiments continue to push the boundaries of precision, the theoretical tools needed to interpret them are equally sharp, allowing for a clearer understanding of the fundamental forces that govern our universe.

Technical Summary: Accurate Wichmann-Kroll Vacuum Polarization Calculations with Gaussian Basis Set

Problem Statement
High-precision predictions of spectra for highly-charged ions require rigorous treatment of Quantum Electrodynamics (QED) corrections, specifically electron self-energy and vacuum polarization. While standard perturbative QED methods suffice for light systems, heavy ions necessitate non-perturbative techniques due to the strong interaction with the nuclear electric field. A critical component of these calculations is the all-order Wichmann-Kroll (WK) vacuum polarization correction, which represents the non-linear contribution beyond the Uehling potential. Historically, accurate numerical results for WK corrections have been obtained via Green's function integration. Although finite-basis-set (Rayleigh-Ritz) methods using B-splines or Gaussian functions have been successfully applied to self-energy and other QED effects, their application to vacuum polarization has faced challenges. Previous attempts using Gaussian bases were limited to sizes of $N \approx 120$ due to linear dependence issues requiring high arithmetic precision, and the resulting numerical values often deviated from standard methods.

Methodology
The authors developed a new C++ program utilizing the Eigen3, Boost::math, and Boost::multiprecision libraries to calculate the WK vacuum polarization correction using an extended Gaussian basis set.

Basis Set Construction: The wave function components (large and small) are approximated by linear combinations of Gaussian functions with parameters $\zeta_i$ forming a geometric progression. The study employs basis sizes up to $N = 500$ .
Precision Handling: To overcome the linear dependence problem inherent in large Gaussian bases, the authors utilized multiple-precision arithmetic (up to 140 significant decimal digits) for the calculation of the vacuum polarization density, converting to double precision only for the final potential and energy shift.
Theoretical Framework: The calculation follows the Furry picture, solving the Dirac equation for heavy hydrogen-like ions. The WK contribution is isolated by subtracting the linear (Uehling) contribution. The authors employ a C-symmetry preserving density formulation to ensure renormalization is handled correctly.
Nuclear Models: Calculations were performed using three nuclear charge distribution models: shell, sphere, and two-parameter Fermi (2pF).
Extrapolation and Error Estimation: The authors analyzed the convergence of results with respect to the basis size $N$ and the geometric progression parameter $\beta$ . They tested inverse polynomial extrapolations and power-law fits to estimate errors arising from incomplete convergence.
Screening Application: The method was extended to lithium-like ions to calculate vacuum polarization screening corrections, utilizing screened potentials derived from the electron core density.

Key Results

Convergence and Accuracy: With a basis size of $N=500$ , the calculated WK correction for the main partial component ( $|\kappa|=1$ ) in hydrogen-like uranium ( $Z=92$ ) converged to four correct decimal digits ($4.4733$ eV), matching the precision of Green's function integration methods.
Partial Wave Convergence: While the absolute value of the correction converges slowly for higher $|\kappa|$ , the authors demonstrated that the difference between results obtained with different nuclear models converges significantly faster. This allows for the calculation of nuclear model corrections using smaller basis sets without performing full, costly calculations for every model.
Error Estimation: The study evaluated various extrapolation techniques. Simple linear or polynomial fits were found to be insufficient for high $|\kappa|$ values. However, using the monotonicity of the energy shift with respect to $\beta$ , the authors established conservative bounds for the error. They found that results up to $|\kappa| \le 10$ are in good agreement with independent Green's function results, with the total WK correction for uranium agreeing with Ref. [2] within the estimated uncertainties.
Screening Effects: For lithium-like uranium, the calculated screening correction for the $2s_{1/2}$ state was found to be $\Delta E_{scr} \approx -0.036$ eV (using core-Hartree potentials) and $\approx -0.031$ eV (using Kohn-Sham potentials), consistent with previous literature values.

Significance and Claims
The paper claims that the application of a large Gaussian basis set ( $N \le 500$ ) with high-precision arithmetic provides an independent computational technique for calculating the all-order Wichmann-Kroll correction that is competitive with the established Green's function integration method.

Validation: The results validate the finite-basis-set approach for vacuum polarization, showing that with sufficient basis size and precision, it can reproduce standard results.
Efficiency: The finding that nuclear model corrections converge faster than absolute values offers a practical advantage, enabling efficient sensitivity studies of nuclear charge distributions.
Future Utility: The authors suggest that this method is promising for molecular QED calculations and for evaluating complex QED diagrams involving vacuum polarization loops, where traditional Green's function integration may be less convenient.
Limitations: The authors modestly note that while the method is accurate for low $|\kappa|$ , convergence slows for higher partial waves, and extrapolation techniques must be applied with caution to avoid subjectivity or overestimation of precision.

More like this