REVIEW 3 major objections 5 minor 1 cited by
Two-loop electron self-energy for low nuclear charges
T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read All-order two-loop self-energy calculations for hydrogen-like ions, extrapolated to hydrogen, give a value that differs from the accepted one by 2.8σ and shifts the Rydberg constant by one standard deviation.
desk verdict Solid all-order numerics, but the hydrogen extrapolation is hostage to unverified analytic coefficients, so the 2.8σ discrepancy is not yet credible. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The enabling technique is an accelerated subtraction scheme for the two-loop self-energy diagrams, generalizing the one-loop method of Ref. [18] to the overlapping and nested diagrams. By subtracting and re-adding diagrams with two explicit interactions with the nuclear binding field, and computing the subtracted contributions in closed form in momentum space without partial-wave expansion, the authors improve the partial-wave convergence by one to two orders of magnitude. This permits calculations down to Z = 5 and reduces the numerical uncertainty at, for example, Z = 10 by a factor of 30.
What would settle it
A direct calculation of the missing part of the B60 coefficient (or an analytic evaluation of B71) that gives a value near the Zα-expansion prediction would shift the extrapolated G60(1) back toward the accepted value; alternatively, a nonperturbative two-loop self-energy calculation for Z = 4 or lower, performed by an independent group, would settle whether the deviation persists.
Extended reading notes
Core claim
The central discovery is that the higher-order remainder function G60(Zα) of the two-loop self-energy, extracted from all-order numerical calculations for nuclear charges Z = 5 to 50, lies outside the band predicted by the Zα expansion and extrapolates to −104.1(3.8) at Z = 1. This disagrees with the previously recommended value of −75.8(9.4), obtained by combining Zα-expansion and older all-order results, by 2.8σ. The authors argue that the improved numerical precision challenges the assumed consistency between the nonperturbative and Zα-expansion approaches, so the theoretical hydrogen Lamb shift must be reexamined, with the updated SESE contribution shifting the 1S–2S prediction by −2.5 kHz and decreasing the Rydberg constant by one standard deviation.
Load-bearing premise
The extrapolation from the calculated nuclear charges Z ≥ 8 to hydrogen assumes the analytic constraints on the logarithmic coefficients B71, B83, B82, B81 of the Zα expansion; the numerical data alone cannot determine these logarithms, so the hydrogen value rests on the very Zα-expansion input whose consistency the paper questions.
Editorial extensions
If this is right
- The theoretical prediction for the hydrogen 1S–2S transition frequency shifts by −2.5 kHz, changing the value of the Rydberg constant by about one standard deviation.
- The previously recommended SESE value for hydrogen, G60(1) = −75.8(9.4), is replaced by −104.1(3.8), indicating that the consistency between all-order and Zα-expansion approaches no longer holds.
- Improved accuracy for the 1S state also refines other nS energy predictions through the relation n^3 E(nS) = ν_n + E(1S).
- The method opens the way to improved calculations of other two-loop QED corrections to the Lamb shift and to the bound-electron g factor.
Reading between the lines
- If the disagreement is real, one of the logarithmic coefficients constrained analytically (e.g., B71) or the partially computed coefficient B60 must be substantially modified; the paper estimates that a missing contribution δB60 ≈ −40 or a 10% adjustment in B61 would reconcile the two approaches.
- A fully independent nonperturbative calculation at Z = 4–7, where the present data are sparse and extrapolation carries the largest weight, would provide a sharp test of the claimed hydrogen value.
- The same accelerated subtraction scheme may reduce the numerical cost of other two-loop bound-state QED calculations, including the g-factor corrections needed for ongoing experiments.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper reports new all-order-in-Zalpha numerical calculations of the two-loop electron self-energy for the 1S state of hydrogen-like ions, extending the accessible range down to Z=5 and improving numerical accuracy by more than an order of magnitude through an accelerated subtraction scheme. The authors extract the higher-order remainder function G60(Zalpha) defined in Eq. (4), extrapolate it to Z=1 using a polynomial fit supplemented by analytic constraints on logarithmic coefficients, and obtain G60(1)=-104.1(3.8). This value differs by 2.8 standard deviations from the previously accepted value of -75.8(9.4), and the paper argues that the all-order and Zalpha-expansion approaches are no longer consistent, with implications for the 1S-2S hydrogen transition and the Rydberg constant.
Significance. If the hydrogen extrapolation is robust, this is an important result: it would improve the two-loop self-energy contribution to the hydrogen Lamb shift, change the theoretical prediction for the 1S-2S transition, and shift the Rydberg constant by about one standard deviation. The numerical all-order data themselves are a substantial advance, with clear improvements in accuracy and internal cross-checks between different computational schemes. However, as detailed below, the headline hydrogen value and the claimed inconsistency depend sensitively on unverified Zalpha-expansion coefficients, so the central claim is not yet established.
major comments (3)
- [Main text, Eq. (4) and Eq. (7)] The extracted G60 values in Table I and the extrapolated hydrogen value G60(1)=-104.1(3.8) are obtained after subtracting the logarithmic terms B63 L^3 + B62 L^2 + B61 L, where B62 and B61 are taken from Pachucki (2001) and are described in the text as "yet to be independently verified." Since L(Z=1)=ln((Zalpha)^-2) is about 9.8 and L^2 about 96.5, a change of about 3 in B61 or about 0.3 in B62 shifts G60(1) by roughly 28, which is comparable to the entire discrepancy with the accepted value -75.8(9.4). The uncertainty quoted in Eq. (7) accounts only for the constraints in Eq. (6) and polynomial extrapolation scatter, not for this dependence. The authors' own statement that "a 10% adjustment in the B61 coefficient ... could restore the consistency" confirms that the central 2.8-sigma discrepancy is not robust to plausible variation of unverified input coefficients. The paper should propagate the uncertainty of B61 and B62 into G60(1), or provide an explicit sensitivity analysis, before claiming a discrepancy.
- [Main text, Eq. (5)-(6); Supplemental Material, 'Details of the extrapolation procedure'] The extrapolation to Z=1 is not a direct all-order result. It uses only data at Z=8,...,50, with Z=5,6,7 excluded, and fits a polynomial to G'_60 after subtracting the B72, B71, and B84 terms using analytic values from the Zalpha expansion. The constraints for B71, B83, B82, and B81 in Eq. (6) come from Ref. [16], the same Zalpha-expansion framework whose consistency with the all-order results is the central claim. The text acknowledges that "the current numerical data cannot constrain the logarithmic coefficients." Thus the claimed inconsistency is partly a joint test of the all-order calculation and these specific analytic inputs, not a pure all-order result. The extrapolation uncertainty should include the uncertainty of these inputs without presupposing their validity, or the conclusion should be reframed accordingly.
- [Supplemental Material, Eqs. (22)-(23)] The exclusion of the three lowest-Z data points (Z=5,6,7) and the choice of weights w_i^2=delta_i^2+(a Z_i^{n+1})^2, with the parameter a selected by requiring that variations from removing data points "are small and approximately equal," are heuristic. These excluded points are the ones closest to hydrogen and the only direct all-order anchor near the extrapolation region. The reported polynomial-extrapolation uncertainty is estimated as twice the maximal difference among P2, P3, and P4 fits, which does not account for model-selection uncertainty beyond these three polynomials. A more systematic assessment, such as exploring dependence on the excluded points, the weight parameter a, and higher-degree polynomials, would strengthen the reliability of the extrapolated G60(1).
minor comments (5)
- [Summary and Conclusion] The phrase "break the previously assumed consistency" is stronger than the evidence supports; consider "challenge" or "question" in light of the caveats discussed in the major comments.
- [Supplemental Material, first paragraph] The text says formulas and figures of the main text will be referred to as "Eq. (?M)" and "Fig. ?M"; these placeholders should be resolved before publication.
- [Main text, Eq. (11)] The ellipsis in Eq. (11) hides the right-hand side of Eq. (9); the notation is understandable from context but could be made explicit.
- [Main text, Fig. 4 caption] The caption does not identify the shaded region or the meaning of the triangle; the text should define these elements clearly.
- [Main text, after Eq. (4)] The numerical values of B62 and B61 are not given; adding them would help readers assess the sensitivity discussed in the major comments.
Circularity Check
Extrapolated hydrogen value G60(1) = -104.1(3.8) is defined by subtracting unverified Zα-expansion coefficients (B61, B62) from the same expansion the paper claims to falsify; a 10% change in B61 would erase the 2.8σ discrepancy.
-
other
[Eq. (4) and following paragraph; extrapolation discussion near Eqs. (5)-(7); Supplemental Material 'Details of the extrapolation procedure'.]
"The coefficients B62 and B61 derived by Pachucki [8] have yet to be independently verified. We now use the analytical results to extract the higher-order remainder function G60(Zα) from our calculations ... Specifically, a 10% adjustment in the B61 coefficient, δB61 ≈ −5, could restore the consistency with the nonperturbative calculations."
G60(Zα) is not a raw all-order observable: Eq. (4) defines it as F(Zα) minus Zα-expansion terms, including B61L and B62L^2 taken from the very approach whose consistency with all-order results is the central claim. Because L ≈ 9.8 at Z = 1, an unverified shift δB61 ≈ −5 changes every extracted G60(Zi) by δB61·L(Zi) and moves G60(1) by roughly −49, which the authors state would restore consistency and erase the headline 2.8σ disagreement. The quoted ±3.8 uncertainty covers polynomial-extrapolation scatter and the B71/B82/B83/B81 constraints of Eq. (6), but not this B61/B62 sensitivity.
full rationale
The all-order numerical F(Zα) data in Table I are genuinely independent of the Zα expansion: they come from a nonperturbative Dirac-Coulomb calculation with an improved subtraction scheme, and they agree with previous all-order results while improving precision by an order of magnitude. There is no self-citation chain: the extrapolation form (Eq. 5) and constraints (Eq. 6) come from Karshenboim et al. (Ref. [16]), not from the present authors, and the comparison with B60 (Ref. [15]) is not used in the fit. The circularity is therefore partial and located at the extraction step: the quantity G60 that is extrapolated to hydrogen is defined by subtracting the unverified B61 and B62 coefficients of the Zα expansion, and the paper itself concedes that a 10% B61 adjustment would restore consistency. The discrepancy and derived Rydberg shift are consequently not robust to the unverified input, but they are not forced by construction; a genuine, independent all-order dataset underlies the claim. Hence score 4 rather than 6 or 8.
Assumptions & free parameters
free parameters (3)
- Polynomial coefficients P_n(Z) =
Not given; P_3 yields G60(0) = -101.1370
- Weight coefficient a =
Not specified
- B71 central value =
-12
assumptions (4)
- domain assumption The Zα-expansion of G60(Zα) has the functional form of Eq. (5) with the specified logarithmic terms.
- domain assumption The constraints on the logarithmic coefficients in Eq. (6) from Ref. [16] are correct and sufficient.
- domain assumption The renormalization/subtraction scheme and the partial-wave expansion converge to the correct two-loop self-energy.
- ad hoc to paper The data at Z=8-50, with weights chosen by Eq. (23), are sufficient to determine the polynomial fit; data at Z=5-7 can be excluded.
Cite this review
Pith. "Pith review of Two-loop electron self-energy for low nuclear charges." pith.science (2026). https://pith.science/paper/N6FHU46U
@misc{pith2026241112459,
author = {Pith},
title = {Pith review of: Two-loop electron self-energy for low nuclear charges},
year = {2026},
howpublished = {\url{https://pith.science/paper/N6FHU46U}},
note = {Machine review of arXiv:2411.12459}
}
abstract
Calculations of the two-loop electron self-energy for the $1S$ Lamb shift are reported, performed to all orders in the nuclear binding strength parameter $Z\alpha$ (where $Z$ is the nuclear charge number and $\alpha$ is the fine structure constant). Our approach allows calculations to be extended to nuclear charges lower than previously possible and improves the numerical accuracy by more than an order of magnitude. Extrapolation of our all-order results to hydrogen yields a result twice as precise as the previously accepted value [E. Tiesinga et al. Rev. Mod. Phys. 93, 025010 (2021)], differing from it by 2.8 standard deviations. The resulting shift in the theoretical prediction for the $1S$-$2S$ transition frequency in hydrogen decreases the value of the Rydberg constant by one standard deviation.
Figures
Forward citations
Cited by 1 Pith paper
-
One-loop electron self-energy with accelerated partial-wave expansion in Coulomb gauge
Coulomb-gauge accelerated partial-wave methods now compute one-loop electron self-energies for hydrogen-like ions from Z=1 to Z=100, including high excited states.
Reference graph
Works this paper leans on
-
[16]
S. G. Karshenboim, A. Ozawa, and V. G. Ivanov, Phys. Rev. A 100, 032515 (2019)
work page 2019
-
[1]
A. Matveev, C. G. Parthey, K. Predehl, J. Alnis, A. Beyer, R. Holzwarth, T. Udem, T. Wilken, N. Ko- lachevsky, M. Abgrall, D. Rovera, C. Salomon, P. Lau- 5 rent, G. Grosche, O. Terra, T. Legero, H. Schnatz, S. Weyers, B. Altschul, and T. W. H¨ ansch, Phys. Rev. Lett. 110, 230801 (2013)
work page 2013
-
[2]
A. Antognini, F. Nez, K. Schuhmann, F. D. Amaro, F. Biraben, J. M. R. Cardoso, D. S. Covita, A. Dax, S. Dhawan, M. Diepold, L. M. P. Fernandes, A. Giesen, A. L. Gouvea, T. Graf, T. W. H¨ ansch, P. Indelicato, L. Julien, C.-Y. Kao, P. Knowles, F. Kottmann, E.- O. L. Bigot, Y.-W. Liu, J. A. M. Lopes, L. Ludhova, C. M. B. Monteiro, F. Mulhauser, T. Nebel, P....
work page 2013
-
[3]
Pachucki, V
K. Pachucki, V. Lensky, F. Hagelstein, S. S. Li Muli, S. Bacca, and R. Pohl, Rev. Mod. Phys. 96, 015001 (2024)
2024
-
[4]
Tiesinga, P
E. Tiesinga, P. J. Mohr, D. B. Newell, and B. N. Taylor, Rev. Mod. Phys. 93, 025010 (2021)
2021
- [5]
- [6]
-
[7]
Weinberg, The Quantum Theory of Fields, Cambridge University Press, 4th edition, 2010
S. Weinberg, The Quantum Theory of Fields, Cambridge University Press, 4th edition, 2010
work page 2010
Show all 40 references
-
[8]
Pachucki, Phys
K. Pachucki, Phys. Rev. A 63, 042503 (2001)
2001
-
[9]
M. I. Eides, H. Grotch, and V. A. Shelyuto, Phys. Rep. 342, 63 (2001)
2001
-
[10]
V. A. Yerokhin, P. Indelicato, and V. M. Shabaev, Phys. Rev. A 71, 040101(R) (2005)
2005
-
[11]
V. A. Yerokhin, Phys. Rev. A 80, 040501(R) (2009)
2009
-
[12]
Beiersdorfer, A
P. Beiersdorfer, A. L. Osterheld, J. H. Scofield, J. R. Crespo L´ opez-Urrutia, and K. Widmann, Phys. Rev. Lett. 80, 3022 (1998)
1998
-
[13]
Beiersdorfer, H
P. Beiersdorfer, H. Chen, D. B. Thorn, and E. Tr¨ abert, Phys. Rev. Lett. 95, 233003 (2005)
2005
-
[14]
V. A. Yerokhin, P. Indelicato, and V. M. Shabaev, Phys. Rev. Lett. 97, 253004 (2006)
2006
-
[15]
Pachucki and U
K. Pachucki and U. D. Jentschura, Phys. Rev. Lett. 91, 113005 (2003)
2003
-
[17]
V. A. Yerokhin, K. Pachucki, and V. Patk´ oˇ s, Ann. Phys. (Leipzig) 531, 1800324 (2019)
2019
-
[18]
Sapirstein and K
J. Sapirstein and K. T. Cheng, Phys. Rev. A 108, 042804 (2023)
2023
-
[19]
A. V. Malyshev, E. A. Prokhorchuk, and V. M. Shabaev, Phys. Rev. A 109, 062802 (2024)
2024
-
[20]
Mallampalli and J
S. Mallampalli and J. Sapirstein, Phys. Rev. Lett. 80, 5297 (1998)
1998
-
[21]
V. A. Yerokhin, Phys. Rev. A 62, 012508 (2000)
2000
-
[22]
Mallampalli and J
S. Mallampalli and J. Sapirstein, Phys. Rev. A 57, 1548 (1998)
1998
-
[23]
V. A. Yerokhin and V. M. Shabaev, Phys. Rev. A 64, 062507 (2001)
2001
-
[24]
V. A. Yerokhin, P. Indelicato, and V. M. Shabaev, Phys. Rev. Lett. 91, 073001 (2003)
2003
-
[25]
V. A. Yerokhin, P. Indelicato, and V. M. Shabaev, Eur. Phys. J. D 25, 203 (2003)
2003
-
[26]
V. A. Yerokhin, Eur. Phys. J. D 58, 57 (2010)
2010
-
[27]
V. A. Yerokhin, Phys. Rev. A 97, 052509 (2018)
2018
-
[28]
See Supplemental Material at [URL will be inserted by publisher] for details of the extrapolation procedure
-
[29]
S. G. Karshenboim and V. G. Ivanov, Phys. Rev. A 98, 022522 (2018)
2018
-
[30]
S. G. Karshenboim and V. A. Shelyuto, Phys. Rev. A 100, 032513 (2019)
2019
-
[31]
2022 CODATA recommended values, https://pml.nist.gov/cuu/Constants/, 2024
2022
-
[32]
V. A. Yerokhin, P. Indelicato, and V. M. Shabaev, Phys. Rev. A 77, 062510 (2008)
2008
-
[33]
Sikora, V
B. Sikora, V. A. Yerokhin, N. S. Oreshkina, H. Cakir, C. H. Keitel, and Z. Harman, Phys. Rev. Research 2, 012002 (2020)
2020
-
[34]
standard
J. Morgner, B. Tu, C. M. K¨ onig, T. Sailer, F. Heiße, H. Bekker, B. Sikora, C. Lyu, V. A. Yerokhin, Z. Har- man, J. R. Crespo L´ opez-Urrutia, C. H. Keitel, S. Sturm, K. Blaum, Nature 622, 53 (2023). SUPPLEMENT AL MA TERIAL: DET AILS OF THE SUBTRACTION SCHEME Here we describe...
2023
-
[35]
(1M) by applying the following subtractions (as shown in the uppermost line of Fig
M term In the standard scheme, the overlapping M term is obtained from Eq. (1M) by applying the following subtractions (as shown in the uppermost line of Fig. 2M): (O, M) : G G G→ G G G− G G(0) G(0) − G(0) G(0) G + G(0) G(0) G(0) − G(0) G(1) G(0) . (9) Here, the three G’s on t...
-
[36]
P term In the standard scheme, the overlapping P term takes care of the two (equivalent) subtraction terms G G(0) G(0) and G(0) G(0) G in Eq. (9). It is defined by applying the following subtractions (see the third line from the top in Fig. 2M): (O, P) : G G(0) G(0) → G G(0) G...
-
[37]
The additional terms subtracted in the M and P terms are (of course) added back and evaluated in the momentum space without any partial-wave expansion
F term The F term takes care of all subtracted terms which now contain only G(0) (and its derivatives) and G(1), but not G’s. The additional terms subtracted in the M and P terms are (of course) added back and evaluated in the momentum space without any partial-wave expansion....
-
[38]
(2M) by applying the following subtractions (see the uppermost line of Fig
M term In the standard scheme, the nested M term is obtained from Eq. (2M) by applying the following subtractions (see the uppermost line of Fig. 3M): (N, M) : G G G→ G G G− G G(0) G − G G(1) G . (14) 7 Here, we ignore the reference-state infrared (IR) divergencies present in ...
-
[39]
(14); they are referred to as the first and as the second, respectively
P term In the standard scheme, there are two nested P terms that take care of the two subtraction terms G G(0) G and G G(1) G in Eq. (14); they are referred to as the first and as the second, respectively. In the standard scheme, the first nested P term is defined by the subtr...
-
[40]
The corresponding diagrams are shown in the two bottom lines of Fig
F term Similarly to the overlapping F term, the nested F term takes care of all subtracted term in the nested M and P terms. The corresponding diagrams are shown in the two bottom lines of Fig. 3M. Their calculation is performed by the technique developed for the F term in our...
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.