{"id":"2233d37b-abac-4a88-8bc2-0fbd81f8441d","arxiv_id":"2607.16026","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"A finite-basis numerical Green function reproduces one-loop self-energy results for hydrogenlike atoms, reaching 10^-4 (Feynman gauge) and 10^-5 (Coulomb gauge) relative precision for hydrogen.","lead":"Using a numerical Green function built from an exponential basis set, the authors compute the one-loop self-energy of hydrogenlike atoms, matching known values to 10^-5 for uranium and 10^-4 for hydrogen. The method is a proof-of-principle step toward self-energy calculations in molecular hydrogen ions, where analytical Green functions do not exist.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Partial-wave extrapolation may be biased by basis-set errors that grow with |κ|; the 20% range-variation test does not bound this systematic error.","rationale":"The reader's weakest assumption correctly identifies the partial-wave extrapolation as the leading source of uncertainty. My stress test sharpens this: the uncertainty estimate itself may be too optimistic because it does not account for the systematic growth of basis-set error with |κ|, which the paper's own data show (e.g., the |κ|=16–35 discrepancy for Z=92 and the increasing Diff. values in Table IV). This does not overturn the paper's core feasibility demonstration, but it directly affects the strength of the central precision claim. The proposed test would distinguish between a benign extrapolation-modeling issue and a real bias, and would either validate or refute the claimed uncertainties. Since the paper is already CONDITIONAL, this concern reinforces the need for the stated reproducibility checks (code, parameters, benchmark data) and adds a specific numerical validation. I therefore recommend keeping the verdict unchanged rather than moving to REJECT or ACCEPT, pending the test.","tokens_in":15861,"tokens_out":6405,"duration_ms":66460,"concrete_test":"Apply the same polynomial-in-1/|κ| extrapolation to the exact analytical Green function partial-wave contributions (from Ref. [28] for Z=1, Ref. [23] for Z=92) over the same fitting windows used in the paper (|κ|=11–15 for Z=1 Coulomb, |κ|=16–20 for Z=92). Compare the extrapolated tail to the exactly computed tail from the next |κ| to infinity. If the exact-data extrapolation error is comparable to or larger than the paper's quoted uncertainty (about 1e-5 relative), the numerical extrapolation is unreliable. If it is small, then repeat the numerical fitting at two basis sizes (e.g., n_G=600 and n_G=800) and check that the resulting tail and its uncertainty are stable; instability would indicate contamination by basis-set error.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The advertised precision for Z=92 (about 10^-5, Table II) and for Z=1 Coulomb (10^-5, Table IV) depends on the polynomial-in-1/|κ| extrapolation of the partial-wave tail. In Table IV, the fitting window |κ|=11–15 is exactly where the differences with the analytical Green function grow from about 3×10^-6 to 1.2×10^-5, and the paper notes that the error is 'dominated by the basis set limitation.' For Z=92, Sec. V A reports that the sum |κ|=16–35 differs from the analytical result by about 2×10^-5, 'indicating that this trend continues at higher |κ|.' The uncertainty is estimated by varying the maximum |κ| by 20%, which probes only the sensitivity to the number of fitted points, not the bias from a systematically growing basis-set error that could be absorbed into the fit. The extrapolated tail contributes the dominant part of the total uncertainty (e.g., the 0.0002(1) tail versus a total 10.3169(1) in Coulomb gauge), so a biased tail would directly undermine the claimed precision. The paper honestly flags the extrapolation as the leading uncertainty, but the quoted error bar may be optimistic.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a finite-basis numerical Green function approach to the one-loop self-energy in hydrogenlike atoms, using exponential basis functions and performing all radial integrations numerically. The many-potential term is computed in both Feynman and Coulomb gauges, while the zero- and one-potential terms are taken from analytical calculations. The method is benchmarked for hydrogenlike uranium (Z=92) and hydrogen (Z=1). For Z=92 in Feynman gauge the total F(Zα) is 1.49090(2), to be compared with the analytical Green function result 1.490916(3); for Z=1 the paper reports F=10.318(1) in Feynman gauge and F=10.3169(1) in Coulomb gauge, compared with 10.31679365(1). Additional contributions include the demonstration that DKB exponential basis sets improve precision substantially over NKB, the implementation of a convergence acceleration scheme, and a complexity reduction for ΔE^1_κ.","tokens_in":16185,"tokens_out":8928,"duration_ms":88162,"significance":"If the precision claims hold, this is a meaningful step toward extending numerical-Green-function self-energy calculations to non-hydrogenic and molecular systems, where analytical Dirac-Coulomb Green functions are unavailable. The benchmarks against independent analytical Green function calculations are a genuine strength, and no parameter is fitted to reproduce the final self-energy. The DKB-basis improvement and the O(max(n_G,n_q)^4) scaling for ΔE^1_κ are concrete, useful results. However, the advertised 10^-5-level uncertainties rest on a partial-wave tail extrapolation whose systematic uncertainty is not fully controlled; this is the main barrier to accepting the precision claims as stated.","major_comments":[{"comment":"The Z=92 total F=1.49090(2) depends on a polynomial-in-1/|κ| extrapolation of the |κ|=16-20 partial-wave contributions. The paper itself notes that the sum |κ|=16-35 differs from the analytical Green function result by about -2.1e-5 and attributes this to basis-set limitations that continue at higher |κ|. The uncertainty is estimated by varying the maximum fitted |κ| by 20%, but this procedure does not bound a systematic trend that grows with |κ| and could be absorbed into the fitted tail. Since the extrapolated tail is the dominant source of uncertainty, please provide a direct cross-check of the tail against analytical high-κ results or otherwise quantify how the observed P35 discrepancy is reflected in the quoted uncertainty.","section":"Sec. V A, Table II"},{"comment":"For Z=1 in Coulomb gauge, the final F=10.3169(1) relies on the extrapolated tail 0.0002(1) for |κ|≥11. The fit window |κ|=11-15 is exactly where deviations from the analytical Green function grow from about 3e-6 to 1.2e-5 and are described as dominated by basis-set limitations. Varying the upper fitted |κ| by 20% tests the sensitivity to the number of fitted points, not the systematic bias from basis-set errors that increase with |κ|. Because the tail contributes the dominant uncertainty, a separate estimate of this systematic bias should be given before the 10^-5 relative-uncertainty claim can be considered established.","section":"Sec. V B, Table IV"}],"minor_comments":[{"comment":"The Diff. entries for the extrapolated rows with 10≤|κ|≤15 appear to be off by a factor of ten. For example, for κ=10, 0.0000426 - 0.0000472 = -4.6e-6, not -4.6e-5; similarly for κ=11 the extrapolated row gives -3.8e-6, not -3.8e-5. Please correct the notation or clarify the convention.","section":"Table IV"},{"comment":"The expression for G^1_κ is written under the assumption x2 < x1, but this is stated rather briefly. It would help to explicitly specify the domain of the (r,y) variables after the transformation (26) so that the quadrature ranges are unambiguous.","section":"Sec. IV B, Eq. (28)"},{"comment":"The text states that the total error on the sum of the first 20 terms is smaller than 4e-6, but the comparison shown is for the sum up to |κ|=35. It would be useful to state explicitly how the error on the first 20 terms is estimated, especially because the P35 discrepancy is attributed to a growing basis-set trend.","section":"Sec. V A, Table II"},{"comment":"The reference values labeled 'Anal. [43]' in Table I come from a private communication. Providing a public source or an ancillary file with these numerical values would improve verifiability.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest and technically credible, and the benchmarks against analytical Green function results are convincing. My main reservation is the partial-wave tail extrapolation: the 20% range-variation test does not by itself bound a systematically growing basis-set error, and the tail is the dominant uncertainty in the final precision claims. If the authors add a direct cross-check of the high-κ tail or demonstrate that the observed discrepancies are already included in the error budget, I would be happy to support acceptance. I see no reason to doubt the basic validity of the method."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a genuine advance: the first calculation of the hydrogen 1s one-loop self-energy using a numerical (basis-set) Green function, reaching 1e-4 in Feynman gauge and 1e-5 in Coulomb gauge, plus a high-Z benchmark at Z=92 at about 2e-5. The methodological pieces are credible — exponential DKB basis, convergence acceleration, Coulomb gauge, and an efficient G1 algorithm that cuts the cost from O(n_q^3 n_G^3) to O(max(n_q,n_G)^4). The paper is honest about its limits and does not oversell.\n\nWhat is solid: the benchmarks against analytical Green function results are the right test, and the final numbers agree within the quoted uncertainties. The partial-wave tail extrapolation is the leading uncertainty, and the paper says so. The 20% range-variation test is a heuristic, and the worry that basis-set errors growing with |kappa| could bias the fit is legitimate. But the actual deviations from the analytical benchmarks sit within the quoted error bars (e.g., Z=92 difference 1.6e-5 vs quoted 2e-5; Z=1 Coulomb difference 1.1e-4 vs quoted 1e-4). So the error bars look roughly right, not obviously optimistic. Still, this is the softest part, and a sharper estimate of the extrapolation bias would strengthen the central claim.\n\nThe main reproducibility gap: no code and no explicit basis-set exponent parameters. Since the method is a proof of feasibility for molecular ions, the community needs to reproduce it. The paper cites a private communication for one intermediate benchmark, but the final comparisons are to published results, so that is not a substantive flaw.\n\nThis is worth a serious referee. It is a methods paper with a clear target audience — precision QED and molecular ion spectroscopy people. I would bring it to a reading group and would cite it if I worked on numerical Green functions for bound-state QED.","headline":"A solid, honest methods paper that delivers the first Z=1 self-energy calculation with a numerical Green function; the error bars match benchmarks, though the partial-wave extrapolation remains the weakest link.","tokens_in":16663,"tokens_out":2829,"would_cite":true,"duration_ms":27835,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A numerical Green function built from exponential basis functions computes the one-loop self-energy of hydrogenlike atoms to about 10^-5 relative precision, extending to hydrogen with a convergence acceleration and Coulomb gauge, opening a","keywords":["self-energy","Lamb shift","numerical Green function","Dirac equation","exponential basis set","hydrogenlike atoms","QED","convergence acceleration"],"falsifier":"Compute the total one-loop self-energy for the hydrogen ground state with an independent high-precision method (e.g., analytic Green function with partial waves carried to millions) and compare to the authors' extrapolated value; a discrepancy larger than their quoted 10^-5 or 10^-4 uncertainty would falsify the extrapolation procedure.","tokens_in":15730,"feed_emoji":"⚛️","tokens_out":3286,"duration_ms":33553,"temperature":0.7,"pith_summary":"The paper demonstrates that a numerical Green function, constructed by solving the radial Dirac equation in an exponential basis set, can evaluate the one-loop self-energy of hydrogenlike atoms with precision comparable to analytic Green function methods. For hydrogenlike uranium (Z=92) in the Feynman gauge, the total self-energy is obtained with about 10^-5 relative uncertainty; for hydrogen (Z=1), the Feynman gauge yields 10^-4 and the Coulomb gauge reaches 10^-5. The central trick is a subtraction and convergence-acceleration scheme that isolates the slowly converging partial-wave tail and handles it via extrapolation. If the method holds up, it opens a practical path to self-energy calculations in molecular ions, where no analytic Green function exists.","feed_headline":"Numerical Green function hits 10^-5 self-energy precision","feed_subtitle":"A finite-basis Dirac Green function with convergence acceleration reaches analytic-level accuracy, opening the door to molecular-ion QED","key_machinery":"The central object is the numerical Dirac radial Green function G_kappa(x1,x2,z) = sum_n phi_{kappa,n}(x1) phi_{kappa,n}(x2)/(z - E_{kappa,n}), where phi and E come from a finite exponential basis set that enforces dual kinetic balance (DKB). The many-potential term of the self-energy is extracted by the subtraction G^{2+} = G - G0 - G1, and the slow partial-wave convergence is tamed by subtracting an analytically known two-potential approximation G^{2+}_a and adding it back in momentum space. This combination makes the partial-wave sum short enough to extrapolate reliably, and the error cancellation among the subtracted terms is the reason the final result is far more accurate than any sing","core_discovery":"The authors show that a finite-basis numerical Green function, built from exponential functions satisfying the dual kinetic balance condition, can reproduce the many-potential term of the one-loop self-energy to high precision. By subtracting the known zero- and one-potential terms and an approximate two-potential term, they convert the partial-wave expansion into a rapidly converging form, then extrapolate the residual tail to infinity using a polynomial in 1/|kappa|. This yields, for the hydrogen ground state, a total self-energy with 10^-4 relative uncertainty in the Feynman gauge and 10^-5 in the Coulomb gauge, and for hydrogenlike uranium about 10^-5 in the Feynman gauge. The achievemen","pith_inferences":["If the exponential basis set can be extended to two-center problems (as the authors imply), the same numerical Green function approach could bring all-order self-energy calculations within reach for molecular ions, where currently only nonrelativistic Z-alpha expansions are used.","The extrapolation method might be sharpened by combining information from both gauges, or by using a more physical model for the tail of the partial-wave series instead of a generic polynomial fit.","The success of this approach hints that many bound-state QED corrections, not just the self-energy, could be computed with numerical Green functions in a basis-set framework, provided the extrapolation tail is handled carefully.","A systematic comparison of DKB exponential basis sets versus Gaussian and B-spline bases, with optimized parameters, could quantify which family is most efficient for high-precision QED calculations in non-hydrogenic systems."],"forward_implications":["The method provides a feasible way to compute one-loop self-energies in systems with no analytic Green function, notably the molecular hydrogen ions H2+ and HD+.","The DKB exponential basis set performs on par with Gaussian basis sets for high-Z atoms, and despite slower convergence it can be extended to low Z with appropriate acceleration.","The Coulomb gauge offers a significant precision advantage over the Feynman gauge for low-Z self-energy calculations, even with a numerical Green function.","The observed strong cancellation of errors between the zero-, one-, and many-potential terms suggests that the final uncertainty is dominated by the partial-wave extrapolation, not by the basis set itself.","Optimizing the nonlinear basis parameters for each partial wave, and pushing the extrapolation to higher |kappa|, could further reduce the uncertainty below 10^-5."],"fun_headline_variants":["Numerical Green function achieves 10^-5 self-energy","Finite-basis Dirac Green function hits 10^-5 accuracy","Self-energy via numerical Green function to 10^-5","Convergence-accelerated Green function reaches 10^-5"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The precision claim rests on the assumption that the partial-wave tail follows the fitted polynomial in 1/|kappa| beyond the computed partial waves; if the tail behaves differently, the extrapolated total self-energy and its quoted uncertainty are unreliable.","fun_headline_variants_meta":{"raw":{"variants":["Numerical Green function achieves 10^-5 self-energy","Finite-basis Dirac Green function hits 10^-5 accuracy","Self-energy via numerical Green function to 10^-5","Convergence-accelerated Green function reaches 10^-5"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000221,"raw_usage":{"total_tokens":1228,"prompt_tokens":630,"completion_tokens":598,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":374,"completion_tokens_details":{"reasoning_tokens":527}},"tokens_in":374,"tokens_out":598,"duration_ms":5876,"temperature":1.0,"reasoning_tokens":527,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T21:33:37.936525+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the total one-loop self-energy for the hydrogen ground state with an independent high-precision method (e.g., analytic Green function with partial waves carried to millions) and compare to the authors' extrapolated value; a discrepancy larger than their quoted 10^-5 or 10^-4 uncertainty would falsify the extrapolation procedure.","supporting_citations":[],"review_version":1}