pith. sign in

arxiv: 2604.23660 · v1 · submitted 2026-04-26 · ⚛️ physics.atom-ph · quant-ph

Wichmann-Kroll vacuum polarization correction to lithium-like systems in a Gaussian basis set

Haisum Hayat , Harry M. Quiney This is my paper

Pith reviewed 2026-05-08 04:56 UTC · model grok-4.3

classification ⚛️ physics.atom-ph quant-ph
keywords Wichmann-Krollvacuum polarizationGaussian basis setsHartree-Fock potentialslithium-like systemsQED correctionsatomic structure
0
0 comments X

The pith

Finite Gaussian basis sets with self-consistent Hartree-Fock potentials compute Wichmann-Kroll vacuum polarization corrections for lithium-like systems.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends finite Gaussian basis sets to the multi-electron case for the alpha(Z alpha) to the n greater than or equal to 3 vacuum polarization. Self-consistent Hartree-Fock potentials are inserted into the calculation for lithium-like systems, producing energy shifts that agree with published results. The work shows the basis-set approach succeeds for atomic potentials whose Green's functions lack simple analytic or numerical forms. A sympathetic reader would see this as a practical route to QED corrections in atoms beyond the one-electron limit.

Core claim

Hartree-Fock potentials obtained self-consistently are used to treat the vacuum polarization for lithium-like systems and are found to be in good agreement with comparable results in the literature. The results presented demonstrate the use of Gaussian basis sets for atomic potentials whose Green's functions expressions cannot be simply obtained via analytic or numerical methods.

What carries the argument

Finite Gaussian basis sets applied to self-consistent Hartree-Fock potentials to evaluate the Wichmann-Kroll term of vacuum polarization.

If this is right

  • The same Gaussian-basis procedure applies directly to other multi-electron ions where analytic Green's functions are unavailable.
  • Energy shifts for s and p states can be obtained reliably once the Hartree-Fock potential is fixed self-consistently.
  • The method supplies a numerical alternative when both analytic expressions and direct numerical integration of the Green's function are impractical.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach could be inserted into larger many-body calculations to add these QED corrections to heavier atoms without rebuilding the entire potential from scratch.
  • Direct comparison of the same Wichmann-Kroll shifts against B-spline or finite-element bases would test whether Gaussian sets offer any systematic advantage in convergence speed or accuracy.
  • Extending the same framework to higher-order vacuum-polarization terms or to the Lamb shift in the same lithium-like ions would reveal how far the basis-set route scales before truncation errors become visible.

Load-bearing premise

The finite Gaussian basis set converges adequately for the multi-electron Hartree-Fock potentials so that truncation errors do not affect the reported agreement with literature values.

What would settle it

A recomputation of the same energy shifts using a substantially larger Gaussian basis set or an independent numerical method that produces values differing by more than the stated agreement would show the basis is insufficient.

Figures

Figures reproduced from arXiv: 2604.23660 by Haisum Hayat, Harry M. Quiney.

Figure 1
Figure 1. Figure 1: FIG. 1. Leading order contributions to the total unrenor view at source ↗
read the original abstract

Recent developments have seen the application of finite Gaussian basis sets to the $\alpha(Z\alpha)^{n\geq3}$ vacuum polarization. The energy shift for $s$ and $p$ electron states have been tabulated and their convergence investigated. In this work, we extend this problem to the multi-electron case. Hartee-Fock potentials obtained self-consistently are used to treat the vacuum polarization for lithium-like systems and are found to be in good agreement with comparable results in the literature. The results presented in this work demonstrate the use of Gaussian basis sets for atomic potentials whose Green's functions expressions cannot be simply obtained via analytic or numerical methods.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript extends the use of finite Gaussian basis sets for computing the Wichmann-Kroll vacuum polarization correction from single-electron cases to lithium-like multi-electron systems. Self-consistent Hartree-Fock potentials are generated numerically and inserted into the VP integral, with the resulting energy shifts reported to be in good agreement with literature values. The work emphasizes that this approach enables treatment of atomic potentials whose Green's functions lack simple analytic forms.

Significance. If the basis-set convergence for the multi-electron HF potentials is shown to be under control, the method offers a practical numerical route to higher-order QED corrections in systems where analytic Green's functions are unavailable. This is a natural and useful extension of prior single-electron Gaussian-basis VP work. The absence of quantitative convergence data, however, leaves the reliability of the reported agreement unverified at present.

major comments (2)
  1. [Abstract] Abstract: the statement that the lithium-like results 'are found to be in good agreement with comparable results in the literature' is unsupported by any numerical values, tables, or error estimates. Without these data the central claim cannot be assessed.
  2. [Results and method sections] Results and method sections: while the single-electron case is stated to have had its convergence investigated, no basis-size enlargement tests, extrapolation, or error bounds are provided for the self-consistent HF potentials used in the multi-electron Wichmann-Kroll integral. Because the potential is generated numerically, any residual incompleteness enters the VP correction directly; this is load-bearing for the validity of the reported agreement.
minor comments (2)
  1. [Abstract] Abstract contains the typo 'Hartee-Fock' (should be 'Hartree-Fock').
  2. [Method] The manuscript would benefit from an explicit statement of the Gaussian basis parameters (exponents, number of functions) employed for the HF potentials and for the subsequent VP integration.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment below and will revise the manuscript to incorporate the suggested improvements.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the statement that the lithium-like results 'are found to be in good agreement with comparable results in the literature' is unsupported by any numerical values, tables, or error estimates. Without these data the central claim cannot be assessed.

    Authors: We agree that the abstract statement would be stronger with explicit support. The results section contains tables with direct numerical comparisons to literature values for the lithium-like systems, including the specific energy shifts obtained. To address the concern, we will revise the abstract to include representative numerical values and a brief mention of the associated uncertainties, allowing the agreement to be assessed directly from the abstract. revision: yes

  2. Referee: [Results and method sections] Results and method sections: while the single-electron case is stated to have had its convergence investigated, no basis-size enlargement tests, extrapolation, or error bounds are provided for the self-consistent HF potentials used in the multi-electron Wichmann-Kroll integral. Because the potential is generated numerically, any residual incompleteness enters the VP correction directly; this is load-bearing for the validity of the reported agreement.

    Authors: We concur that explicit convergence data for the numerically generated self-consistent Hartree-Fock potentials is essential. Although the single-electron convergence was documented in our prior work, we will add to the revised manuscript a dedicated subsection with basis-size enlargement tests for the lithium-like cases. This will include results for successively larger Gaussian basis sets, any extrapolations performed, and quantitative error bounds on the Wichmann-Kroll shifts, confirming that incompleteness effects are under control. revision: yes

Circularity Check

0 steps flagged

No circularity; numerical VP corrections validated against independent external literature

full rationale

The paper extends prior Gaussian-basis treatments of Wichmann-Kroll vacuum polarization to multi-electron lithium-like ions by replacing the nuclear potential with self-consistent Hartree-Fock potentials and evaluating the resulting energy shifts numerically. The reported values are obtained directly from the basis-set integrals and then compared to separate literature results for validation; no quantity is defined in terms of itself, no fitted parameter is relabeled as a prediction, and no load-bearing step reduces to a self-citation chain or an ansatz imported from the authors' own prior work. The derivation therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Assessment based on abstract only; no explicit free parameters, axioms, or invented entities are stated. The method inherits standard Hartree-Fock and Gaussian-basis assumptions from prior literature without introducing new postulates.

pith-pipeline@v0.9.0 · 5401 in / 1208 out tokens · 61890 ms · 2026-05-08T04:56:37.802290+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages

  1. [1]

    and 1/rij is the Coulomb repulsion between electrons iandj. The energy minimization is achieved by solving the set of coupled equations for the Fock operator fi =h i + X j (Jj −K j),(7) in which the electron-electron interaction is expressed through the direct operator Ji(x1)ϕ(x1) =⟨ψ i(x2)| 1 r12 |ψi(x2)⟩ϕ(x 1),(8) and the exchange operator Ki(x1)ϕ(x1) =...

    work page

  2. [2]

    W. E. Lamb and R. C. Retherford, Fine structure of the hydrogen atom by a microwave method, Phys. Rev.72, 241 (1947)

    work page 1947

  3. [3]

    P. J. Mohr, D. B. Newell, B. N. Taylor, and E. Tiesinga, CODATA recommended values of the fundamental phys- ical constants: 2022*, Journal of Physical and Chemical Reference Data54, 033105 (2025)

    work page 2022

  4. [4]

    D. P. Aguillard, T. Albahri, D. Allspach, J. Annala, K. Badgley, S. Baeßler, I. Bailey, L. Bailey, E. Barlas- Yucel, T. Barrett, E. Barzi, F. Bedeschi, M. Berz, M. Bhattacharya, H. P. Binney, P. Bloom, J. Bono, E. Bottalico, T. Bowcock, S. Braun, M. Bressler,et al. 6 (Muong−2 Collaboration), Measurement of the positive muon anomalous magnetic moment to 12...

    work page 2025

  5. [5]

    Shabaev, Two-time Green’s function method in quantum electrodynamics of high-Zfew-electron atoms, Physics Reports356, 119 (2002)

    V. Shabaev, Two-time Green’s function method in quantum electrodynamics of high-Zfew-electron atoms, Physics Reports356, 119 (2002)

    work page 2002

  6. [6]

    Hangele and M

    T. Hangele and M. Dolg, Coupled-cluster and DFT stud- ies of the Copernicium dimer including QED effects, Chem. Phys. Lett.616-617, 222 (2014)

    work page 2014

  7. [7]

    Nguyen, J

    T. Nguyen, J. Lowe, T. Pham, I. Grant, and C. Chantler, Electron self-energy corrections using the Welton concept for atomic structure calculations, Radiat. Phys. Chem. 204, 110644 (2023)

    work page 2023

  8. [8]

    E. A. Uehling, Polarization effects in the positron theory, Phys. Rev.48, 55 (1935)

    work page 1935

  9. [9]

    Serber, Linear modifications in the Maxwell field equa- tions, Phys

    R. Serber, Linear modifications in the Maxwell field equa- tions, Phys. Rev.48, 49 (1935)

    work page 1935

  10. [10]

    E. H. Wichmann and N. M. Kroll, Vacuum polarization in a strong Coulomb field, Phys. Rev.101, 843 (1956)

    work page 1956

  11. [11]

    W. H. Furry, A symmetry theorem in the positron theory, Phys. Rev.51, 125 (1937)

    work page 1937

  12. [12]

    G. A. Rinker and L. Wilets, Vacuum polarization in high- Z, finite-size nuclei, Phys. Rev. Lett.31, 1559 (1973)

    work page 1973

  13. [13]

    Gyulassy, Higher order vacuum polarization for finite radius nuclei, Nucl

    M. Gyulassy, Higher order vacuum polarization for finite radius nuclei, Nucl. Phys. A244, 497 (1975)

    work page 1975

  14. [14]

    Soff and P

    G. Soff and P. J. Mohr, Vacuum polarization in a strong external field, Phys. Rev. A38, 5066 (1988)

    work page 1988

  15. [15]

    Persson, I

    H. Persson, I. Lindgren, S. Salomonson, and P. Sunner- gren, Accurate vacuum-polarization calculations, Phys. Rev. A48, 2772 (1993)

    work page 1993

  16. [16]

    A. N. Artemyev, T. Beier, G. Plunien, V. M. Shabaev, G. Soff, and V. A. Yerokhin, Vacuum-polarization screen- ing corrections to the energy levels of heliumlike ions, Phys. Rev. A62, 022116 (2000)

    work page 2000

  17. [17]

    Beier, A

    T. Beier, A. N. Artemyev, G. Plunien, V. M. Shabaev, G. Soff, and V. A. Yerokhin, Vacuum-Polarization Screening Corrections to the Low-Lying Energy Levels of Heliumlike Ions., Hyperfine Interact.132, 369 (2001)

    work page 2001

  18. [18]

    A. N. Artemyev, V. M. Shabaev, and V. A. Yerokhin, Vacuum polarization screening corrections to the ground- state energy of two-electron ions, Phys. Rev. A56, 3529 (1997)

    work page 1997

  19. [19]

    Persson, S

    H. Persson, S. Salomonson, P. Sunnergren, and I. Lind- gren, Two-electron Lamb-shift calculations on heliumlike ions, Phys. Rev. Lett.76, 204 (1996)

    work page 1996

  20. [20]

    V. A. Yerokhin, Z. Harman, and C. H. Keitel, QED cal- culations of the 2p−2stransition energies in Li-like ions, Phys. Rev. A112, 042801 (2025)

    work page 2025

  21. [21]

    A. N. Artemyev, T. Beier, G. Plunien, V. M. Shabaev, G. Soff, and V. A. Yerokhin, Vacuum-polarization screen- ing corrections to the energy levels of lithiumlike ions, Phys. Rev. A60, 45 (1999)

    work page 1999

  22. [22]

    I. P. Grant, The Dirac Equation, inRelativistic Quan- tum Theory of Atoms and Molecules: Theory and Com- putation(Springer New York, New York, NY, 2007) pp. 121–179

    work page 2007

  23. [23]

    V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii,Quantum Electrodynamics: Volume 4, Vol. 4 (Butterworth-Heinemann, 1982)

    work page 1982

  24. [24]

    V. G. Bagrov and D. Gitman,The Dirac Equation and its Solutions(De Gruyter, 2014)

    work page 2014

  25. [25]

    Szabo and N

    A. Szabo and N. S. Ostlund,Modern quantum chem- istry: introduction to advanced electronic structure theory (Courier Corporation, 2012)

    work page 2012

  26. [26]

    Shavitt and R

    I. Shavitt and R. J. Bartlett,Many-body methods in chemistry and physics: MBPT and coupled-cluster the- ory(Cambridge university press, 2009)

    work page 2009

  27. [27]

    S. P. Goldman, Variational representation of the Dirac- Coulomb Hamiltonian with no spurious roots, Phys. Rev. A31, 3541 (1985)

    work page 1985

  28. [28]

    P. A. M. Dirac, Hamiltonian methods and quantum me- chanics, Math. Proc. R. Ir. Acad.63, 49 (1963)

    work page 1963

  29. [29]

    Salman and T

    M. Salman and T. Saue, Calculating the many-potential vacuum polarization density of the Dirac equation in the finite-basis approximation, Phys. Rev. A108, 012808 (2023)

    work page 2023

  30. [30]

    F. Jiyu, C. Shouwan, and H. Taihua, Solution to the Dirac equation using the finite difference method, Nucl. Sci. Tech.31, 10.1007/s41365-020-0728-6 (2025)

    work page doi:10.1007/s41365-020-0728-6 2025

  31. [31]

    Zhang, Y

    Y. Zhang, Y. Bao, H. Shen, and J. Hu, Resolving the spurious-state problem in the dirac equation with the finite-difference method, Phys. Rev. C106, L051303 (2022)

    work page 2022

  32. [32]

    Wallmeier and W

    H. Wallmeier and W. Kutzelnigg, Use of the squared dirac operator in variational relativistic calculations, Chem. Phys. Lett.78, 341 (1981)

    work page 1981

  33. [33]

    Kutzelnigg, Basis set expansion of the dirac operator without variational collapse, Int

    W. Kutzelnigg, Basis set expansion of the dirac operator without variational collapse, Int. J. Quantum Chem.25, 107 (1984)

    work page 1984

  34. [34]

    V. M. Shabaev, I. I. Tupitsyn, V. A. Yerokhin, G. Plu- nien, and G. Soff, Dual Kinetic Balance Approach to Basis-Set Expansions for the Dirac Equation, Phys. Rev. Lett.93, 130405 (2004)

    work page 2004

  35. [35]

    V. K. Ivanov, S. S. Baturin, D. A. Glazov, and A. V. Volotka, Vacuum-polarization Wichmann-Kroll correc- tion in the finite-basis-set approach, Phys. Rev. A110, 032815 (2024)

    work page 2024

  36. [36]

    I. P. Grant and H. M. Quiney, GRASP: The Future?, Atoms10, 10.3390/atoms10040108 (2022)

    work page doi:10.3390/atoms10040108 2022

  37. [37]

    K. G. Dyall and K. Faegri,Introduction to Relativis- tic Quantum Chemistry(Oxford University Press, 2007) Chap. One-Electron Atoms

    work page 2007

  38. [38]

    Kim, Relativistic self-consistent-field theory for closed-shell atoms, Phys

    Y.-K. Kim, Relativistic self-consistent-field theory for closed-shell atoms, Phys. Rev.154, 17 (1967)

    work page 1967

  39. [39]

    R. M. Pitzer, Relativistic Electronic Structure Theory: Part 1, Fundamentals, J. Am. Chem. Soc.125, 13304 (2003)

    work page 2003

  40. [40]

    I. P. Grant, ed., Computation of atomic structures, in Relativistic Quantum Theory of Atoms and Molecules: Theory and Computation(Springer New York, New York, NY, 2007) pp. 393–432

    work page 2007

  41. [41]

    Benazzouk, M

    R. Benazzouk, M. Salman, and T. Saue, Wichmann-Kroll vacuum polarization density in a finite Gaussian basis set, Phys. Rev. A113, 032813 (2026)

    work page 2026

  42. [42]

    Sapirstein and K

    J. Sapirstein and K. T. Cheng, Vacuum polarization calculations for hydrogenlike and alkali-metal-like ions, Phys. Rev. A68, 042111 (2003)

    work page 2003

  43. [43]

    M. W. Schmidt and K. Ruedenberg, Effective conver- gence to complete orbital bases and to the atomic Hartree–Fock limit through systematic sequences of Gaussian primitives, J. Chem. Phys.71, 3951 (1979)

    work page 1979

  44. [44]

    E. U. Condon and G. H. Shortley,The Theory of Atomic Spectra(Cambridge university press, 1935)

    work page 1935

  45. [45]

    Sapirstein and K

    J. Sapirstein and K. T. Cheng, Calculation of the Lamb shift in neutral alkali metals, Phys. Rev. A66, 042501 (2002)

    work page 2002