QED vacuum polarization in the Coulomb field of a nucleus: a method of high-order calculation
Pith reviewed 2026-05-18 11:07 UTC · model grok-4.3
The pith
Corrections to QED vacuum polarization in a nuclear Coulomb field up to order α²(Zα)⁷ are obtained by reducing two-loop bound-state diagrams to free QED graphs with up to eight loops.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The calculation of the QED vacuum polarization potential corrections up to order α² (Zα)⁷ is achieved through a reduction of the two-loop Feynman graphs with proper propagators of fermions in the external field to a set of free QED Feynman graphs with up to eight independent loops.
What carries the argument
Reduction of bound-state vacuum polarization diagrams to free QED Feynman graphs with up to eight loops, which carries the evaluation of the α²(Zα)⁷ corrections.
If this is right
- The method makes it possible to compute vacuum-polarization shifts in atoms with large nuclear charge Z at the required precision.
- Treating Zα as a separate parameter allows controlled expansions for heavy nuclei where the field is strong.
- The same reduction technique supplies a template for other two-loop QED corrections that involve external Coulomb fields.
- Results obtained this way can be inserted directly into calculations of atomic energy levels and transition frequencies.
Where Pith is reading between the lines
- The reduction may simplify calculations of vacuum polarization at still higher orders in α and Zα.
- Similar diagram reductions could be tested on other QED observables in strong static fields.
- Precision measurements on highly charged ions offer a direct experimental check of the extracted corrections.
Load-bearing premise
The reduction from diagrams that contain external-field propagators to a set of free QED graphs fully accounts for every contribution at the stated order.
What would settle it
An independent numerical evaluation of one two-loop vacuum-polarization diagram performed both with the external-field propagators and through the reduced free QED graphs, followed by a mismatch at order α²(Zα)⁷, would refute the equivalence.
Figures
read the original abstract
A calculation of the QED vacuum polarization potential in the Coulomb field of a pointlike nucleus was presented in an earlier publication by the author and his collaborators. Corrections up to order $\alpha^2 (Z\alpha)^7$ were evaluated, where $Z$ is the nuclear charge number and $Z\alpha$ is treated as an independent variable. These corrections correspond to two-loop Feynman graphs with proper propagators of fermions in the external field. The calculation employed a reduction to free QED, leading to free QED Feynman graphs with up to eight independent loops. The method of calculation is described here in detail.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript details a method for calculating the QED vacuum polarization potential in the Coulomb field of a pointlike nucleus. It describes the evaluation of corrections up to order α²(Zα)⁷, obtained by reducing two-loop bound-state vacuum polarization diagrams (with external-field fermion propagators) to free QED Feynman graphs containing up to eight independent loops, while treating Zα as an independent expansion parameter.
Significance. If the reduction procedure is complete, the work provides a systematic framework for high-order QED calculations in strong external fields. This is relevant for precision atomic spectroscopy and tests of QED in high-Z systems. The mapping to free QED diagrams is a strength, as it permits cross-checks against independent free-QED results and supports reproducibility through explicit construction of the diagram set.
major comments (1)
- [Reduction to free QED] The reduction mapping (detailed in the section describing the procedure) must explicitly demonstrate that all contributions through (Zα)⁷ are captured by the free-QED graphs; an omitted term from the expansion of the external-field propagators would undermine the order-by-order claim.
minor comments (2)
- The abstract refers to an earlier publication but does not include a citation; adding the reference would allow readers to connect the present method description to the prior numerical results.
- [Introduction] Clarify in the introduction how the independent treatment of Zα affects the power counting and ensures no double-counting occurs when combining the two-loop structure with the external-field expansion.
Simulated Author's Rebuttal
We thank the referee for the positive overall assessment and the recommendation of minor revision. We address the single major comment below.
read point-by-point responses
-
Referee: [Reduction to free QED] The reduction mapping (detailed in the section describing the procedure) must explicitly demonstrate that all contributions through (Zα)⁷ are captured by the free-QED graphs; an omitted term from the expansion of the external-field propagators would undermine the order-by-order claim.
Authors: We agree that an explicit demonstration of completeness strengthens the presentation. The reduction is performed by expanding each external-field fermion propagator in powers of Zα and converting the resulting insertions into additional photon lines attached to free propagators. In the revised manuscript we add a short subsection that makes the power counting explicit: the expansion of each propagator is truncated after terms that would produce a total power of Zα greater than 7, and the complete set of free-QED diagrams (up to eight loops) that survive this truncation is enumerated. This construction guarantees that every contribution through α²(Zα)⁷ is retained and that no lower-order term is omitted. revision: yes
Circularity Check
No significant circularity; derivation maps to independent free QED graphs
full rationale
The paper presents a systematic reduction of two-loop bound-state vacuum polarization diagrams (with external-field propagators) to a collection of free QED diagrams containing up to eight independent loops, expanded in powers of Zα treated as an independent parameter. This mapping is constructed explicitly to account for all contributions through order α²(Zα)⁷, with the resulting free QED graphs being standard objects whose evaluation is independent of the original bound-state problem and can be verified against separate calculations. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations that collapse the central claim are present in the described procedure. The derivation remains self-contained against external benchmarks in free QED.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
QED vacuum polarization in the Coulomb field of a nucleus: a method of high-order calculation
useful formulas concerning ˜V13(r). An explicit expres- sion for ˜V13(r) was reported in [4]. The Wichmann-Kroll ∗ sergey.volkov.1811@gmail.com, volkoff˙sergey@mail.ru, svolkov@mpi-hd.mpg.de 1 Different corrections of this type may have different factors f (m, M ), where m and M are the fermion and nuclear masses. 2 Strictly speaking, we use the renormalized...
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[2]
and other scientists [21–26]; there exist also applica- tions to other high-order quantum field theory problems [27]. In contrast to the case of free-electron g − 2, the free- QED Feynman graphs of VVP do not have infrared (IR) divergences. Therefore, a special subtraction procedure for removing IR divergences is not required: we can ap- ply the BPHZ 3 ren...
-
[3]
The added vertex, marked by a boxed dot, embodies momentum conservation for Coulomb lines; it does not contribute any factor and carries no tensor index. We refer to it as the fictitious vertex . Dashed lines denote Coulomb propagators. The Feynman gauge is used; the 3 FIG. 2. All Feynman graphs contributing to V25. FIG. 3. All Feynman graphs contributing ...
-
[4]
001 0 . 062303(73) 0 . 32 0 . 5 0 . 0039 − 5. 6 × 10− 8 1 0 . 044156(40) 0 . 013 0 . 00069 − 0. 00001 0 1000 0 . 0000013361(10) 0 . 0000011 2 . 6 × 10− 7 1. 1 × 10− 9 0 100000 2 . 3622(18) × 10− 10 1. 1 × 10− 11 9. 2 × 10− 11 4. 5 × 10− 13 − 9. 5 × 10− 22 Evaluation of V25
-
[5]
001 0 . 014832(67) 0 . 044 0 . 023 0 . 0021 − 1. 6 × 10− 15 1 0 . 012235(52) 0 . 013 0 . 0012 6 × 10− 6 0 1000 6 . 649(25) × 10− 7 8. 7 × 10− 7 − 7. 2 × 10− 8 − 3. 5 × 10− 10 − 6. 3 × 10− 69 100000 1 . 1819(44) × 10− 10 3. 5 × 10− 10 − 2. 8 × 10− 11 − 1. 2 × 10− 13 − 1. 4 × 10− 65 Evaluation of V27
-
[6]
001 0 . 00824(14) 0 . 0043 − 0. 00097 − 7. 3 × 10− 6 − 6. 1 × 10− 22 1 0 . 00708(12) 0 . 0032 − 0. 00089 − 6. 4 × 10− 6 0 1000 4 . 596(70) × 10− 7 3. 8 × 10− 7 − 6. 3 × 10− 8 − 4. 9 × 10− 10 2. 5 × 10− 67 100000 8 . 10(12) × 10− 11 7. 2 × 10− 11 − 1. 8 × 10− 11 − 1. 3 × 10− 13 − 2. 2 × 10− 67 TABLE II. Statistics of round-off errors in the integrand eva lu...
-
[7]
001 6 × 1011 71 × 109 47 × 109 29 × 108 44 × 106 34 × 105 1 21 × 109 12 × 108 23 × 107 5 × 105 0 0 1000 15 × 1010 13 × 109 27 × 108 72 × 105 0 0 100000 5 × 1011 82 × 109 75 × 108 21 × 106 2 0 Evaluation of V25
-
[8]
001 7 × 1010 12 × 109 31 × 108 2 × 108 67 × 105 45 × 104 1 11 × 109 17 × 108 12 × 107 22 × 104 0 0 1000 33 × 109 7 × 109 49 × 107 14 × 105 569 1 100000 87 × 109 23 × 109 18 × 108 11 × 106 15 × 104 3232 Evaluation of V27
-
[9]
The column “Object” refers to the graph number from Fig
001 42 × 108 15 × 108 46 × 106 16 × 105 33870 1456 1 4 × 109 14 × 108 27 × 106 49427 0 0 1000 82 × 108 36 × 108 84 × 106 288675 200 0 100000 19 × 109 87 × 108 4 × 108 29 × 105 44360 941 TABLE III: Contributions of individual Feynman graphs to ˜V23(p) for several values of p and ( m0)2 from (4) and (5). The column “Object” refers to the graph number from F...
work page 1984
-
[10]
in tabular form. Their application to high-precison calculations of energy levels of atomic systems was also discussed there. These values agree well with the result of order α 2(Zα )3 at p = 0 from [13]. The values of orders α 2(Zα )5 and α 2(Zα )7, as well as those of order α 2(Zα )3 at p ̸= 0 were presented for the first time. The methodology of the cal...
-
[11]
Indelicato, QED tests with highly charged ions, J
P. Indelicato, QED tests with highly charged ions, J. Phys. B 52, 232001 (2019)
work page 2019
-
[12]
M. I. Eides, H. Grotch, and V. A. Shelyuto, Theory of Light Hydrogenic Bound States (Springer Berlin, Heidelberg, 2007)
work page 2007
-
[13]
V. A. Yerokhin and V. M. Shabaev, Lamb shift of n=1 and n=2 states of hydrogen-like atoms, 1 ≤ Z ≤ 110, J. Phys. Chem. Ref. Data 44, 033103 (2015)
work page 2015
-
[14]
Blomqvist, Vacuum polarization in exotic atoms, Nuclear Physics B 48, 95 (1972)
J. Blomqvist, Vacuum polarization in exotic atoms, Nuclear Physics B 48, 95 (1972)
work page 1972
-
[15]
E. A. Uehling, Polarization Effects in the Positron The- ory, Phys. Rev. 48, 55 (1935)
work page 1935
-
[16]
E. H. Wichmann and N. M. Kroll, Vac- uum Polarization in a Strong Coulomb Field, Phys. Rev. A 101, 843 (1956)
work page 1956
- [17]
-
[18]
N. L. Manakov, A. A. Nekipelov, and A. G. Fainshtein, Vacuum polarization by a strong coulomb field and its contribution to the spectra of multiply-charged ions, Zh. Eksp. Teor. Fiz. 95, 1167 (1989), [Sov. Phys. JETP 68, 673 (1989)]
work page 1989
-
[19]
A. G. Fainshtein, N. L. Manakov, and A. A. Nekipelov, Vacuum polarization by a Coulomb field. Analytical approximation of the polarization potential, J. Phys. B 24, 559 (1991)
work page 1991
-
[20]
N. L. Manakov and A. A. Nekipelov, Proceedings of Voronezh State University 2, 53 (2012), [in Russian]
work page 2012
-
[21]
G. K¨ allen and A. Sabry, Fourth Order Vacuum Polarization, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 29 (1955)
work page 1955
-
[22]
L. W. Fullerton and J. G.A. Rinker, Accurate and efficient methods for the evaluation of vacuum- polarization potentials of order Zα and Zα 2, Phys. Rev. A 13, 1283 (1976)
work page 1976
-
[23]
P. A. Krachkov and R. N. Lee, O(mα 2(Zα )6) contribution to Lamb shift from radiative cor- rections to the Wichmann-Kroll potential, Journal of High Energy Physics 2023, 1 (2023)
work page 2023
-
[24]
S. Volkov and V. A. Yerokhin and Z. Harman and C. H. Keitel, Two-loop vacuum polarization in a Coulomb field, 2509.03284 (2025)
-
[25]
V. A. Yerokhin, P. Indelicato, and V. M. Shabaev, Evaluation of the two-loop self-energy correction to the ground state energy of H-like ions to all orders in Zα , European Physical Journal D 25, 203 (2003)
work page 2003
-
[26]
S. Mallampalli and J. Sapirstein, Fourth-order vacuum-polarization contribution to the Lamb shift, Phys. Rev. A 54, 2714 (1996)
work page 1996
-
[27]
I. Goidenko, L. Labzowsky, A. Nefioodv, G. Plunien, and G. Soff, Second-Order Electron Self-Energy in Hydrogen- like Ions, Physical Review Letters 83, 2312 (1999)
work page 1999
-
[28]
L. N. Labzowsky and A. O. Mitrushenkov, Renormal- ization of the second-order electron self-energy for a tightly bound atomic electron: A detailed derivation, Phys. Rev. A 53, 3029 (1996)
work page 1996
- [29]
-
[30]
Volkov, Calculation of the total 10th order QED contribution to the electron magnetic moment, Phys
S. Volkov, Calculation of the total 10th order QED contribution to the electron magnetic moment, Phys. Rev. D 110, 036001 (2024), 2404.00649
-
[31]
M. J. Levine and J. Wright, Anomalous magnetic mo- ment of the electron, Phys. Rev. D 8, 3171 (1973)
work page 1973
-
[32]
R. Carroll and Y. Yao, α 3 contributions to the anomalous magnetic moment of an electron in the mass-operator formalism, Phys. Lett. B 48, 125 (1974)
work page 1974
-
[33]
P. Cvitanovi´ c and T. Kinoshita, Sixth Order Magnetic Moment of the electron, Phys. Rev. D 10, 4007 (1974)
work page 1974
-
[34]
T. Kinoshita and W. B. Lindquist, Eighth-Order Anomalous Magnetic Moment of the Electron, Phys. Rev. Lett. 47, 1573 (1981)
work page 1981
-
[35]
T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Tenth-Order QED Contribution to the Electron g − 2 and an Improved Value of the Fine Structure Constant, Phys. Rev. Lett. 109, 111807 (2012), 1205.5368
work page internal anchor Pith review Pith/arXiv arXiv 2012
- [36]
-
[37]
L. T. Adzhemyan and M. V. Kompaniets, Five-loop nu- merical evaluation of critical exponents of the ϕ 4 theory, in Journal of Physics: Conference Series , 15th Interna- tional Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT2013), Vol. 523 (IOP Publishing, 2014) p. 012049, 1309.5621
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[38]
N. N. Bogoliubov and O. S. Parasiuk, On the Multipli- cation of the causal function in the quantum theory of fields, Acta Math. 97, 227 (1957), [in German]
work page 1957
-
[39]
Hepp, Proof of the Bogolyubov-Parasiuk theorem on renormalization, Commun
K. Hepp, Proof of the Bogolyubov-Parasiuk theorem on renormalization, Commun. Math. Phys. 2, 301 (1966)
work page 1966
-
[40]
Zimmermann, Convergence of Bogoliubov’s method of renormalization in momentum space, Commun
W. Zimmermann, Convergence of Bogoliubov’s method of renormalization in momentum space, Commun. Math. Phys. 15, 208 (1969)
work page 1969
-
[41]
S. Zschocke, G. Plunien, and G. Soff, Renormalization of the two-photon vacuum polarization and the self en- ergy vacuum polarization for a tightly bound electron, Eur. Phys. J. D 19, 147 (2002)
work page 2002
-
[42]
O. I. Zavialov and B. M. Stepanov, Yad. Fys. 1, 922 (1965), [in Russian; English trans. : Soviet J. Nucl. Phys. 1, 658 (1965)]
work page 1965
-
[43]
P. Cvitanovi´ c and T. Kinoshita, New Approach to the Separation of Ultraviolet and Infrared Divergences of Feynman - Parametric Integrals, Phys. Rev. D 10, 3991 (1974)
work page 1974
-
[44]
S. Volkov, Subtractive procedure for calculating the anomalous electron magnetic moment in QED and its ap- plication for numerical calculation at the three-loop leve l, Zh. Eksp. Teor. Fiz. 149, 1164 (2016), [JETP 122, 1008 (2016)]
work page 2016
-
[45]
S. Volkov, Calculation of lepton magnetic mo- ments in quantum electrodynamics: A justifica- tion of the flexible divergence elimination method, 18 Phys. Rev. D 109, 036012 (2024), 2308.11560
-
[46]
N. N. Bogoliubov and D. V. Shirkov, Introduction to the Theory of Quantized Fields (John Wiley & Sons, Inc., 1980)
work page 1980
-
[47]
O. I. Zavyalov, Renormalized quantum field theory (Dor- drecht, Netherlands: Kluwer, 1990)
work page 1990
-
[48]
V. A. Smirnov, Renormalization and Asymptotic Expansions (Birkh¨ auser, 2000)
work page 2000
-
[49]
J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields (McGraw-Hill College, New York, 1965)
work page 1965
-
[50]
Speer, Analytic Renormalization, J
E. Speer, Analytic Renormalization, J. Math. Phys. 9, 1404 (1968)
work page 1968
-
[51]
S. Volkov, New method of computing the contri- butions of graphs without lepton loops to the electron anomalous magnetic moment in QED, Phys. Rev. D 96, 096018 (2017), 1705.05800
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[52]
Borinsky,Tropical Monte Carlo quadrature for Feynman integrals, Ann
M. Borinsky, Tropical Monte Carlo quadrature for Feynman integrals, Annales de l’Institut Henri Poincar´ e D 10, 635 (2023), 2008.12310
-
[53]
M. Borinsky, H. J. Munch, and F. Tellander, Trop- ical Feynman integration in the Minkowski regime, Computer physics communications 292, 108874 (2023), 2302.08955
-
[54]
Volkov, Infrared and ultraviolet power count- ing on the mass shell in quantum electrodynamics, Nucl
S. Volkov, Infrared and ultraviolet power count- ing on the mass shell in quantum electrodynamics, Nucl. Phys. B 961, 115232 (2020), 1912.04885
-
[55]
Numerical calculation of high-order QED contributions to the electron anomalous magnetic moment
S. Volkov, Numerical calculation of high-order QED con - tributions to the electron anomalous magnetic moment, Phys. Rev. D 98, 076018 (2018), 1807.05281
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[56]
S. Volkov, Calculating the five-loop QED contribution to the electron anomalous magnetic moment: Graphs without lepton loops, Phys. Rev. D 100, 096004 (2019), 1909.08015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.