{"id":"43b64ea0-80d2-4fd1-8075-d4a5f03b14b8","arxiv_id":"2411.19135","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"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.","lead":"Electrons in atoms feel a tiny correction from their own electromagnetic field, and this paper computes that correction for hydrogen and hydrogen-like ions more efficiently than before. The new scheme should improve comparisons between quantum electrodynamics theory and precision atomic spectroscopy, including measurements used to set fundamental constants.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unresolved 50-sigma discrepancy for 2S1/2 at Z=10 is the load-bearing weak spot; the paper calls it 'minor' but the difference is 2.8e-5 in F(Zα), many times the quoted uncertainty.","rationale":"The reader's conditional verdict is appropriate, but the decisive concern is more specific than the general extrapolation model. Table III contains a direct, quantitative disagreement with a benchmark cited as validation: 2S1/2 at Z=10, 4.8944161(2) vs 4.8944444(6). This is not an extrapolation artifact in the abstract; it is a concrete inconsistency in the paper's own result table, acknowledged only as 'minor.' The 1S Z=10 entry also differs by about 4-5 sigma, reinforcing the pattern. The hydrogen results agreeing with Jentschura and Mohr and many other entries matching literature are real independent support, so a full REJECT is not warranted. However, until the Z=10 2S discrepancy is either resolved or traced to a specific numerical step, the paper's broad claim of validated accuracy for arbitrary excited states remains conditional. The proposed test isolates whether the tail-extrapolation/uncertainty procedure is at fault, which is the load-bearing numerical premise.","tokens_in":31460,"tokens_out":6800,"duration_ms":58846,"concrete_test":"Recompute the 2S1/2 Z=10 entry with the same YPS Coulomb-gauge code but extend explicit partial waves from |κmax|=35-40 to |κmax|=80-100, and repeat the tail extrapolation with a wider family of (l0,l1) fits and a larger |κmax| variation than 20%. If the total moves from 4.894416 toward 4.894444, the error budget is systematically underestimated and the central claim needs revision; if it stays at 4.8944161(2), the discrepancy with Indelicato and Mohr must be resolved by an independent code or by a documented Z-alpha expansion cross-check before 'arbitrary excited states' can be claimed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the Coulomb-gauge YPS/SC scheme is accurate for arbitrary excited states down to Z=1. The most direct threat is a benchmark inside the paper: Table III gives F(2S1/2,Z=10)=4.8944161(2) from this work, versus 4.8944444(6) from Indelicato and Mohr [19]. The difference is 2.8e-5, roughly 50 times the stated combined uncertainty. Section IV dismisses this as a 'minor discrepancy,' but no explanation or error analysis is provided. If the present value is correct, a high-precision literature value is wrong; if the literature value is correct, the paper's uncertainty estimate for this state is wrong by orders of magnitude. Either way, the claim of validated accuracy across excited states is not established at this point. The likely source is the partial-wave tail extrapolation of Section III (fits to 1/|κ|^l, uncertainty from 20% variation of |κmax|), since 2S at Z=10 is precisely where numerical cancellations and state-dependent tail behavior are nontrivial. That the discrepancy is not reconciled, and is labeled 'minor' without quantitative support, is a missing piece of evidence for the headline claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript generalizes the potential-expansion approach for the one-loop electron self-energy to the general covariant gauge and the Coulomb gauge, and implements the YPS and Sapirstein–Cheng accelerated partial-wave convergence schemes in the Coulomb gauge. It presents numerical results for hydrogen-like ions in the range Z=1–100 for states up to n=5, including D5/2, F, and G states, and reports comparisons with Jentschura–Mohr, Mohr–Kim, and Indelicato–Mohr values. The central claims are that the Coulomb gauge avoids low-Z numerical cancellations, that the accelerated schemes give high accuracy down to Z=1, and that the approach is applicable to arbitrary excited states. Most comparisons agree with the literature at the quoted level, but the 2S1/2 Z=10 entry of Table III disagrees with the high-precision result of Ref. [19] by far more than the combined uncertainties, and the paper does not analyze this discrepancy.","tokens_in":31762,"tokens_out":6948,"duration_ms":76876,"significance":"If validated, this would be a significant methodological and numerical contribution: it provides an independent all-order cross-check of the Jentschura–Mohr hydrogen values, extends direct all-order self-energy calculations to D5/2, F, and G states for Z<60, and shows a practical accuracy gain from combining the Coulomb gauge with YPS/SC acceleration. The paper does not ship code, but the extended tables in the supplementary material are useful for benchmarks. However, the unresolved Table III discrepancy and the heuristic nature of the partial-wave extrapolation uncertainty prevent full confidence in the headline claim until those points are addressed.","major_comments":[{"comment":"The value F(2S1/2, Z=10)=4.8944161(2) differs from the Indelicato–Mohr result 4.8944444(6) by 2.8×10^-5, roughly 45 times the combined quoted uncertainty. The text in §IV dismisses this as a 'minor discrepancy' without quantitative analysis. This is not a minor issue because the entry lies in the benchmark region (moderate Z, an excited state with nontrivial numerical cancellations) that is supposed to validate the method, and it is far outside the error bars of both calculations. Please provide a concrete explanation: either correct the present value, show with a controlled numerical experiment that the quoted uncertainty is underestimated, or demonstrate a specific reason to question the literature value.","section":"§IV, Table III (2S1/2 row, Z=10)"},{"comment":"The uncertainty of the final values is set by comparing extrapolations obtained when varying |κmax| by 20% and by choosing among fits of the form δE_|κ|=Σ c_l/|κ|^l. This procedure assumes that the asymptotic model Eq. (55) is an accurate representation of the tail; if the true tail contains state-dependent logarithmic or oscillatory terms, the extrapolated contribution can be biased by more than the estimated uncertainty. In view of the 2S1/2 Z=10 discrepancy, the paper should demonstrate with a known-benchmark state (for example, the 2P states with accurate literature values) that the extrapolated tail is stable under changes of (l0,l1), the number of fitted terms, and the |κmax| window, and that the quoted errors cover the resulting spread.","section":"§III, partial-wave extrapolation paragraph (Eq. (55))"},{"comment":"The estimate A70=0±8A60 is calibrated by comparing the Zα expansion with the all-order results for the hydrogen 2P states. This tuned uncertainty is then used to assess whether the new all-order results for other states at Z=1 and Z=5 agree with the Zα-expansion predictions. The cross-check is therefore partially circular: agreement cannot be stronger than the calibration chosen for A70. Please either use an independent estimate of A70 (for example, from known higher-order coefficients or from a different data set) or formulate the low-Z comparisons without relying on this tuned uncertainty.","section":"§IV, discussion of A70 (around Eq. (57))"}],"minor_comments":[{"comment":"The statement that the Coulomb gauge is 'the optimal choice' is stronger than the evidence in Table I, which compares gauges for a single state (1s, Z=10) plus the hydrogen 1s case; a more precise phrasing would be 'advantageous for low-Z states tested here'.","section":"§IV"},{"comment":"The phrase 'arbitrary excited reference states' overstates the demonstrated scope; the results cover n≤5 and l≤4 for the listed states. Qualify this claim or provide results for higher-n/higher-l states.","section":"§IV"},{"comment":"The sentence 'with only minor discrepancies observed for the 1S and 2S states at Z=10' should be revised to reflect the quantitative analysis requested in the major comments; at present the statement is not supported by the numbers in Table III.","section":"§IV"},{"comment":"The contour parameters δy, δx1, δx2 are chosen empirically; a sentence reporting the sensitivity of the final values to their variation would improve reproducibility.","section":"§III"},{"comment":"The reference to the supplementary material [54] still contains a placeholder URL; please insert the final link.","section":"Supplementary material"}],"recommendation":"major_revision","confidential_remarks":"This is a strong technical paper from a group with a long track record in this area, and most benchmarks agree at the stated level. The unresolved 2S1/2 Z=10 discrepancy is the single decisive obstacle; in my view it should be addressed with a quantitative error analysis and, if possible, an independent recalculation of that entry, before the paper is accepted. I would not recommend rejection unless the authors are unable to resolve the discrepancy."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a substantial numerical paper with real new content, but the unresolved 2S1/2 discrepancy at Z=10 is a genuine weak spot that the paper shrugs off far too quickly.\n\nWhat's new and good: implementing the YPS and SC convergence-acceleration schemes in Coulomb gauge works, and the paper convincingly shows that Coulomb gauge avoids the spurious α(Zα)^2 cancellations that plague Feynman gauge at low Z, with better partial-wave convergence. The hydrogen cross-check against Jentschura and Mohr is a legitimate independent verification at the 1e-7 level. The first direct all-order self-energy values for D5/2, F, and G states below Z=60 fill a real gap, and the extended tables for Z=1-100, n≤5 are a usable resource. The presentation is clear, with enough detail to reproduce the derivations, and the citation pattern is appropriate.\n\nSoft spots: Table III lists 2S1/2 at Z=10 as 4.8944161(2) against 4.8944444(6) from Indelicato and Mohr. That's a difference of 2.8e-5, about 50 times the combined uncertainty. The paper calls it a 'minor discrepancy' and moves on. That's not defensible. Either the uncertainty estimate for this state is wrong by orders of magnitude, or the literature value is wrong, and the reader has no way to tell which. The most plausible source is the partial-wave tail extrapolation — fits to 1/|κ|^l with uncertainty from varying |κmax| by 20% — but the paper doesn't give enough detail to audit it. There are also smaller but still uncomfortable differences at 1S and 2P3/2 at Z=10 with the same reference. The rest of the comparisons, including hydrogen and the higher-Z data, look fine, so the method is probably sound, but the Z=10 anomaly is load-bearing for the claim of validated accuracy across arbitrary excited states.\n\nRecommendation: this deserves peer review. It's a serious, careful piece of work with valuable results, but the authors should be required to either fix the Z=10 2S value, properly inflate the uncertainty, or provide a quantitative explanation for why the Indelicato-Mohr number is off. I wouldn't cite that particular value until it's sorted out, but the high-l results and the gauge comparison are worth engaging with.","headline":"A careful, valuable extension of accelerated partial-wave methods to Coulomb gauge, but the 2S1/2 Z=10 discrepancy with Indelicato and Mohr is too large to dismiss as 'minor' and needs a real explanation before the headline accuracy claim is credible.","tokens_in":32286,"tokens_out":3724,"would_cite":false,"duration_ms":36954,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.30.jr","31.15.-p"],"model":"deepseek-v4-flash","headline":"The Coulomb gauge makes one-loop electron self-energy calculations work all the way down to hydrogen, without expansions in the nuclear binding strength.","keywords":["electron self-energy","Coulomb gauge","partial-wave expansion","hydrogen-like ions","QED corrections","Lamb shift","accelerated convergence","bound-state QED"],"falsifier":"Recompute the hydrogen $1s$ self-energy without the tail extrapolation—for instance, by summing partial waves explicitly to $|\\kappa|\\sim 10^4$ in multiprecision arithmetic, as was done in earlier benchmark calculations—and compare with the paper's value $F_{\\rm SE}=10.3167935(7)$; a difference at the $10^{-7}$ level would invalidate the extrapolation-based uncertainty estimate.","tokens_in":31275,"feed_emoji":"⚛️","tokens_out":8218,"duration_ms":64545,"temperature":0.7,"pith_summary":"The paper's central claim is that the Coulomb gauge, not the customary Feynman gauge, is the right setting for all-order (in $Z\\alpha$) numerical computations of the one-loop electron self-energy in low-nuclear-charge ions. In this gauge the spurious $\\alpha(Z\\alpha)^2$-order terms that cause severe numerical cancellations for $Z<5$ are absent, and the partial-wave sum converges faster. By combining the Coulomb gauge with two accelerated-convergence schemes (the YPS and SC schemes), the authors obtain accurate self-energy values for hydrogen-like ions across the whole range of $Z$, including hydrogen itself and excited states up to $n=5$. The method's reach is concrete: it extends direct computations to $D_{5/2}$, $F$, and $G$ states that were previously out of reach for $Z<60$, and it provides an independent cross-check of the most precise hydrogen self-energy values used in the determination of the Rydberg constant.","feed_headline":"Coulomb gauge tames hydrogen's electron self-energy","feed_subtitle":"One-loop self-energy now computed to $10^{-7}$ for hydrogen and all excited states up to $n=5$.","key_machinery":"The central object is the partial-wave expansion of the many-potential Dirac-Coulomb Green function $G^{(2+)}$ inside the self-energy loop; its slow convergence and the tail extrapolation of its partial-wave terms limit numerical accuracy. The machinery is the potential-expansion decomposition of the self-energy into zero-, one-, and many-potential terms, together with a subtraction-and-readdition trick: an explicitly computable approximate Green function $G^{(2+)}_a$ that captures the slowest-converging part is subtracted from $G^{(2+)}$ in the partial-wave sum and re-added in closed form. The YPS scheme uses $G^{(2+)}_a=G^{(0)}(\\varepsilon+\\Omega)-G^{(0)}(\\varepsilon)-\\Omega\\,\\partial_\\varepsilon G^{(0)}$ with $\\Omega=2Z\\alpha/(x_1+x_2)$; the SC scheme uses $G^{(2+)}_a=\\frac12 V(x_1)\\partial^2_\\varepsilon G^{(0)} V(x_2)$. In the Coulomb gauge these manipulations are carried through with the transverse photon propagator split into $D_1$ and $D_2$ parts, and the resulting radial integrals are evaluated with Whittaker-function representations of the Coulomb-Dirac Green function.","core_discovery":"The paper establishes that previous low-$Z$ difficulties in self-energy calculations were in large part a gauge artifact. In the Feynman gauge, individual zero-, one-, and many-potential contributions contain spurious terms of order $\\alpha(Z\\alpha)^2$ that cancel only after summation, producing catastrophic cancellations as $Z\\alpha$ shrinks; the physical result starts at order $\\alpha(Z\\alpha)^4$. Working in the Coulomb gauge removes these spurious contributions and simultaneously improves the convergence of the partial-wave expansion. Implementing the accelerated-convergence schemes of YPS (summing multiple commutators with an effective potential $\\Omega=2Z\\alpha/(x_1+x_2)$) and of SC (commuting two Coulomb potentials out of the Green function) in the Coulomb gauge, the authors compute the dimensionless self-energy function $F_{\\rm SE}(Z\\alpha)$ for $Z=1$ through 100 and principal quantum numbers through $n=5$. For hydrogen's $1s$ state they obtain $F_{\\rm SE}=10.3167935(7)$, consistent with the earlier benchmark value $10.316793650(1)$, and they present the first direct all-order values for $D_{5/2}$, $F$, and $G$ states for $Z<60$.","pith_inferences":["The same Coulomb-gauge acceleration approach should transfer to self-energy calculations with a finite nuclear size or with a screening potential, since the acceleration acts on the free-electron part of the Green function and does not depend on the specific bound-state potential; this would extend the method beyond the point-nucleus hydrogen-like ions treated here.","For high-$\\ell$ states the $Z$-dependence of $F_{\\rm SE}(Z\\alpha)$ is smooth, so interpolating the tabulated values would give accurate self-energy corrections at arbitrary $Z$ in the range 1–100 without recomputing each point.","The improved low-$Z$ values can tighten tests of the $Z\\alpha$-expansion, in particular pinning down the unknown $A_{70}$ coefficient, which the paper currently bounds only by a scale estimate of $0\\pm 8A_{60}$.","A comparison of YPS and SC results at fixed $Z$ for states where the two schemes differ by more than the quoted uncertainty would give a direct empirical check of the extrapolation error, since the two rely on different approximate Green functions."],"forward_implications":["For hydrogen-like ions with $Z=1$ through roughly 10, the method delivers self-energy values with 2–7 more significant digits than previous all-order approaches for states beyond the low-$n$, low-$\\ell$ set, enabling sharper comparisons with spectroscopic data.","First direct all-order results for $D_{5/2}$, $F$, and $G$ states at $Z<60$ become available, providing benchmarks for the $Z\\alpha$-expansion coefficients (such as $A_{60}$ and $A_{70}$) that enter these high-angular-momentum levels.","The agreement with the benchmark hydrogen values provides an independent check of the theoretical input to the Rydberg constant determination.","Because the Coulomb-gauge formulation is free of the low-$Z$ cancellations, the same subtraction machinery can be applied to the two-loop self-energy and vertex corrections in the low-$Z$ region, a direction the paper explicitly identifies as the next target."],"supporting_citations":[{"why":"Supplies the potential-expansion method that splits the self-energy into zero-, one-, and many-potential terms, and the analysis of spurious $\\alpha(Z\\alpha)^2$ contributions.","marker":"[21]"},{"why":"Provides the detailed Feynman-gauge scheme and integration-contour prescriptions that the present work generalizes to Coulomb gauge.","marker":"[25]"},{"why":"Introduces the YPS accelerated-convergence approximation for the many-potential Green function that is implemented here in Coulomb gauge.","marker":"[30]"},{"why":"Introduces the SC accelerated-convergence approach used here as a second, independent acceleration scheme.","marker":"[31]"},{"why":"Supplies the high-precision hydrogen self-energy benchmark values against which the new results are checked.","marker":"[13]"},{"why":"Provides the comparison values at $Z=10$ and other low-$Z$ ions used to validate the new calculations.","marker":"[19]"},{"why":"Provides the Coulomb-gauge radial integrals and regularized Bessel functions on which the many-potential term is built.","marker":"[44]"},{"why":"Supplies the Coulomb-gauge renormalized free-electron vertex operator used in the one-potential term.","marker":"[38]"}],"fun_headline_variants":["Coulomb gauge cuts self-energy cancellations","Coulomb gauge unlocks hydrogen self-energy precision","Self-energy tamed by Coulomb gauge choice","Coulomb gauge fixes low-Z self-energy calculations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quoted numerical results and their error bars rely on the assumption that the partial-wave tail is well described by a polynomial in $1/|\\kappa|$, so that fitting the last few partial-wave terms and extrapolating to infinity yields the true remainder.","fun_headline_variants_meta":{"raw":{"variants":["Coulomb gauge cuts self-energy cancellations","Coulomb gauge unlocks hydrogen self-energy precision","Self-energy tamed by Coulomb gauge choice","Coulomb gauge fixes low-Z self-energy calculations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000665,"raw_usage":{"total_tokens":3035,"prompt_tokens":946,"completion_tokens":2089,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":2033}},"tokens_in":562,"tokens_out":2089,"duration_ms":16510,"temperature":1.0,"reasoning_tokens":2033,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:29:17.308481+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the hydrogen $1s$ self-energy without the tail extrapolation—for instance, by summing partial waves explicitly to $|\\kappa|\\sim 10^4$ in multiprecision arithmetic, as was done in earlier benchmark calculations—and compare with the paper's value $F_{\\rm SE}=10.3167935(7)$; a difference at the $10^{-7}$ level would invalidate the extrapolation-based uncertainty estimate.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the potential-expansion method that splits the self-energy into zero-, one-, and many-potential terms, and the analysis of spurious $\\alpha(Z\\alpha)^2$ contributions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the YPS accelerated-convergence approximation for the many-potential Green function that is implemented here in Coulomb gauge."},{"cited_title":"Sapirstein and K","cited_arxiv_id":null,"evidence_quote":"Introduces the SC accelerated-convergence approach used here as a second, independent acceleration scheme."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the high-precision hydrogen self-energy benchmark values against which the new results are checked."},{"cited_title":"Indelicato and P","cited_arxiv_id":null,"evidence_quote":"Provides the comparison values at $Z=10$ and other low-$Z$ ions used to validate the new calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Coulomb-gauge renormalized free-electron vertex operator used in the one-potential term."}],"review_version":1}