{"id":"ed1f3bed-0ca3-4a3b-bd59-99ca3dcd9643","arxiv_id":"2607.12168","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Improved all-order-in-αZ Wichmann-Kroll correction values for the interelectronic interaction in He- and Li-like ions, including a new finite-nuclear-size term and a flagged discrepancy with Ref. [31].","lead":"This paper sharpens the calculation of a subtle quantum-electrodynamics correction to how electrons interact inside two- and three-electron heavy ions, cutting numerical error bars by orders of magnitude. It adds a new finite-nuclear-size term and reports a dispute with a 1999 published calculation that the field will need to settle.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Multipole tail extrapolation at |κ|=16 is unbenchmarked and the Li-like discrepancy rests on an uninspectable private communication; both need independent checks before the claimed precision can be trusted.","rationale":"The reader's verdict correctly identifies the numerical extrapolation and the reliance on a private communication as the weakest assumptions. My analysis agrees with that assessment. The reader recommends CONDITIONAL, and I see no reason to change that verdict. The concern is not that the paper is wrong, but that its central precision claim is not independently verified: the multipole tail is the dominant numerical uncertainty and its estimation procedure is only self-consistent, not benchmarked. The Li-like discrepancy is an unresolved issue that the paper itself acknowledges, and the only supporting evidence is uninspectable. These are sufficient to keep the verdict CONDITIONAL, requiring the authors to either ship the numerical machinery/data or provide an independent benchmark for the tail and a public version of Malyshev's calculation. No red flags of misconduct or circularity are present; the paper is transparent about its limitations and the He-like agreement with Ref. [30] provides a sanity check. Thus the conditional verdict stands.","tokens_in":12427,"tokens_out":4932,"duration_ms":54038,"concrete_test":"For a representative high-Z ion (e.g., He-like uranium, Z=92), recompute the partial-wave contributions to E^(2+)_eVPe up to κ,κ′=±24 or ±32, and compare the sum of the added terms (|κ|>16) with the tail predicted by the paper's polynomial extrapolation. If the difference exceeds the paper's quoted uncertainty, the extrapolation is unreliable. Additionally, request from A. Malyshev a public, inspectable account (preprint or data table) of the independent Li-like calculation; if the values cannot be independently reproduced, the Li-like discrepancy remains unresolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central deliverable is a set of numerical values for E^(2+)_eVPe with uncertainties reduced by 'several orders of magnitude' relative to Refs. [30,31]. This requires that the truncation of the Dirac–Coulomb propagator partial-wave expansion at κ,κ′=±16 and the subsequent extrapolation (Sec. III) are reliable. The extrapolation is performed by fitting polynomials in inverse powers of |κ| along diagonals, with the error estimated by varying the truncation point by 25% and the polynomial degree. This is an internal consistency check, not a validation: it assumes the tail has a polynomial form and that the observed variation under truncation changes captures the true error. If the asymptotic behavior contains terms not represented by the chosen polynomials (e.g., logarithmic or oscillatory contributions), the tail could be systematically biased. The He-like agreement with Ref. [30] (Table I) provides a benchmark only at the coarser precision of that earlier calculation; it does not test the tail at the claimed higher precision. A second unresolved issue is the Li-like discrepancy: the paper explicitly states it could not identify the origin of the disagreement with Ref. [31], and the only independent support is a private communication from A. Malyshev (Ref. [47]), which the reader cannot inspect. If the numerical extrapolation is off, all the new central values are suspect; if the Li-like agreement with Malyshev is incorrect or misinterpreted, the Li-like half of the payload, including the implicit claim that Ref. [31] is superseded, is unsupported. Both are load-bearing because they are required for the central claim of improved, reliable numbers.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an all-order-in-αZ calculation of the Wichmann–Kroll (higher-order vacuum-polarization) correction to the one-photon exchange interelectronic interaction in He- and Li-like ions, denoted E^(2+)_eVPe. The method uses the Dirac–Coulomb Green's function, an angular multipole expansion, and a subtraction of spurious non-gauge-invariant terms (Sec. II.C) to improve convergence of the radial and loop-energy integrations. Numerical results are reported for Z=20–100 for point nuclei and for finite nuclear size (Tables I and II). He-like results are in agreement with Artemyev et al. [30] with substantially smaller quoted uncertainties; Li-like results disagree with Artemyev et al. [31] but agree with an independent calculation communicated privately by A. Malyshev [47]. The finite-nuclear-size correction is found to exceed 5% for some heavy ions. The authors argue the method transfers to the two-loop self-energy–vacuum-polarization diagram.","tokens_in":12679,"tokens_out":8489,"duration_ms":83395,"significance":"If correct, this work materially improves the precision of a hard QED correction in few-electron high-Z ions, with direct relevance to ongoing Lamb-shift and g-factor experiments and to the all-order treatment of the two-loop SVPE diagram. The derivation is self-contained, the Uehling part is handled analytically, and the error budget is explicit: quadrature convergence, multipole truncation with polynomial tail extrapolation, and nuclear-model dependence. The He-like benchmark against Artemyev et al. [30] gives confidence in the low-multipole part. However, the two load-bearing pillars — the extrapolation beyond κ,κ′ = ±16 and the Li-like discrepancy — are not independently verified in the manuscript: the former is tested only by an internal consistency check, and the latter rests on a private communication. Both need to be strengthened before the claimed accuracy can be accepted.","major_comments":[{"comment":"The claimed improvement of 'several orders of magnitude' over Refs. [30,31] rests entirely on the reliability of the tail estimate for the partial-wave expansion of the two Dirac-Coulomb propagators, truncated at κ,κ′ = ±16. The paper estimates the tail by fitting polynomials in inverse powers of |κ| and assigns an uncertainty by varying the cutoff by 25% and the polynomial degree. That is an internal consistency check, not a validation: it assumes the tail is a polynomial in 1/|κ| with no logarithmic or oscillatory contributions, and it cannot detect a systematic bias of that ansatz. The He-like agreement in Table I validates the low multipoles at the older, coarser precision; it does not test the tail at the new claimed precision. I ask the authors to provide an independent check for at least one or two high-Z cases — for example, recomputation with κmax = 24 and 32, or a comparison wi","section":"Sec. III, multipole truncation"},{"comment":"The Li-like values are a central deliverable of the paper, but the disagreement with Ref. [31] for uranium is left unresolved, and the only supporting evidence is a private communication from A. Malyshev [47] that the reader cannot inspect. This is load-bearing because the manuscript offers no other way to decide between Ref. [31] and the present calculation. I strongly recommend that the authors either (i) include the numerical data from [47] in a table or appendix, or arrange a citable, inspectable publication of that independent calculation; (ii) identify and resolve the source of the discrepancy with Ref. [31] (for example, the sign of ω, the angular reduction, or the summation over core magnetic substates); or (iii) clearly mark the Li-like results as provisional and reduce the stated accuracy. As written, the Li-like half of the paper cannot be independently verified.","section":"Sec. IV, Li-like results"}],"minor_comments":[{"comment":"The Ref. [31] values are quoted without their original uncertainties and without specifying whether they include finite nuclear size. Since the Li-like discrepancy is a key issue, the complete comparison with errors should be shown.","section":"Table II, Z=92a row"},{"comment":"The ω-dependent Uehling expression is not the standard static form; please add a brief derivation or an explicit reference for the branch and prefactor conventions in the exponential.","section":"Eq. (21)"},{"comment":"FLINT is described as 'Fast Library for Number Theory'; the Whittaker-function evaluation likely uses Arb [44] or another special-functions package. Please verify the citations and state which software actually evaluates the special functions.","section":"Refs. [43,44]"},{"comment":"The convergence criterion (50% more integration points changes the result by <1e-6) should be stated as applying after the spurious-term subtraction and for all Z and states; also clarify how the radial infinity cutoff is chosen.","section":"Sec. III / Fig. 3"},{"comment":"The asserted Z^6 scaling in the low-Z region is not quantified. A fit or an explicit leading-Z^6 formula would make the statement testable and more informative.","section":"Sec. IV, low-Z scaling"},{"comment":"The notation E^(2+)_eVPe is used before it is defined; please define it explicitly at first use and state that it is the Wichmann–Kroll part after Uehling subtraction.","section":"Eq. (18)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a plausible strong-field QED calculation, and the He-like agreement with Ref. [30] is reassuring. The main obstacles are verification rather than formalism: the multipole-tail extrapolation must be validated by an independent method at the claimed precision, and the Li-like comparison currently rests on an uninspectable private communication. I would encourage the editor to ask that the Malyshev data be made available as supplementary material or in a citable form, and that at least one high-Z case be recomputed with a larger κ cutoff or an independent method before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth your attention if you care about bound-state QED in heavy ions, but don't treat the Li-like numbers as settled. The He-like half is solid: it reproduces Artemyev's old results, adds a new finite-nuclear-size correction, and the explicit error budget is believable. The Li-like half is the problem: the paper finds a discrepancy with the 1999 table and supports its values only with a private communication from Malyshev. That's transparent, but it's not checkable by a referee.\n\nWhat's new: improved numerical values for the Wichmann-Kroll screening correction to one-photon exchange in He- and Li-like ions, with uncertainties reduced by orders of magnitude; a new finite-nuclear-size term that can shift the effect by more than 5% for high Z; and a subtraction scheme for spurious non-gauge-invariant terms that speeds up the radial and energy integrals. The subtraction method is clearly explained and demonstrated (Fig. 3). The paper is honest about its limitations—it explicitly says it couldn't locate the source of the Li-like discrepancy.\n\nThe soft spots are two. First, the claimed accuracy gain rests on a multipole tail extrapolation beyond |κ|=16 using polynomial fits, and the uncertainty is estimated by internal consistency checks (varying truncation point and polynomial degree). That's standard practice in this field, but it isn't an independent validation. The He-like agreement with Artemyev only benchmarks the low multipoles at the old precision, not the tail at the new claimed precision. So the 'several orders of magnitude' improvement is plausible, not proven. Second, the Li-like discrepancy is resolved only by an uninspectable private communication. If that comparison is wrong, the Li-like table is wrong. The authors deserve credit for flagging it, but the burden is real.\n\nNo p-hacking or circular reasoning: the calculation is self-contained, and the results are not fitted to any target. Citation pattern is normal for the group.\n\nBottom line: I'd send this to peer review rather than desk-reject. The method and He-like values are publishable now; the Li-like claims need either more detail from Malyshev or an independent calculation. If you need the He-like numbers or the finite-nuclear-size correction, cite this. For the Li-like values, wait until the discrepancy is resolved.","headline":"Solid He-like values and a useful method; Li-like numbers rest on a private communication and an unbenchmarked tail extrapolation—needs referee scrutiny, not desk rejection.","tokens_in":13315,"tokens_out":3222,"would_cite":true,"duration_ms":35150,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Improved all-order Wichmann-Kroll values for interelectronic interaction in He- and Li-like ions.","keywords":["Wichmann-Kroll","vacuum polarization","interelectronic interaction","highly charged ions","Lamb shift","bound-state QED","finite nuclear size","Dirac-Coulomb Green function"],"falsifier":"Recompute the correction with the multipole truncation raised to κ,κ′ = ±32 (or with a different subtraction/regularization) and compare the extrapolated tail contribution to the quoted uncertainty; a deviation larger than the stated error would falsify the precision claim. For the Li-like discrepancy, an independently published recalculation that resolves the two-body values would settle which set of numbers is correct.","tokens_in":12203,"feed_emoji":"⚛️","tokens_out":5546,"duration_ms":51062,"temperature":0.7,"pith_summary":"This paper tries to pin down the higher-order vacuum-polarization (Wichmann-Kroll) contribution to the interelectronic interaction in helium-like and lithium-like highly charged ions—the diagram in which a virtual electron-positron loop is inserted into the exchanged photon. It develops a subtraction of spurious, non-gauge-invariant pieces inside the integrand, which accelerates numerical convergence by roughly two orders of magnitude, and then computes the correction for Z = 20–100 with uncertainties several orders smaller than previous values. The He-like results agree with the earlier calculation, while the Li-like results disagree with it and are instead checked against an unpublished independent calculation. The paper also computes the finite-nuclear-size correction, which reduces the effect by more than 5% for heavy ions. If correct, these numbers replace the earlier references for the interelectronic Wichmann-Kroll contribution and tighten Lamb-shift predictions in strong-field QED tests.","feed_headline":"Interelectronic QED correction recalculated to new precision","feed_subtitle":"Wichmann-Kroll loop values for Z=20–100, with finite-nuclear-size shifts up to 5%, tighten Lamb-shift predictions.","key_machinery":"The central object is the Wichmann-Kroll loop in the photon line, expressed through two Dirac-Coulomb Green's functions in the nuclear field, expanded in multipoles κ,κ′ and photon partial waves. The load-bearing mechanism is the subtraction scheme: replacing the two bound propagators by a second-derivative expression with free propagators and the Coulomb potential removes non-gauge-invariant spurious terms that vanish after integration, shortening the integrand's tail and accelerating energy/radial integration by about two orders of magnitude. The multipole tail is handled by polynomial extrapolation in inverse powers of |κ| beyond κ,κ′ = ±16; a Wick rotation to the imaginary axis and Gauß","core_discovery":"Central claim: all-order-in-αZ numerical values for E^(2+)_eVPe, the Wichmann-Kroll correction to the one-photon-exchange interelectronic interaction, for the He-like ground state and the 2s1/2, 2p1/2, 2p3/2 states of Li-like ions. Using the Dirac-Coulomb Green's function, the authors compute point-nucleus values plus a finite-nuclear-size difference, find a sign change near Z≈74 in the He-like ground state, and report a finite-size correction exceeding 5% for heavy ions. He-like values agree with the earlier calculation; Li-like values disagree and are supported by a private independent calculation. The subtraction method transfers, the authors argue, to the two-loop self-energy-vacuum-pola","pith_inferences":["A natural test of the tail extrapolation would be to push the multipole series to κ,κ′ = ±32 or use a mixed-gauge subtraction; if the extrapolated tail shifts by more than the quoted uncertainty, the 'several orders of magnitude' claim would need revision.","If the Li-like discrepancy is due to an error in the old calculation, existing experimental data on high-Z lithium-like ions may already prefer the new values; reanalyzing those data with the new correction would be a check.","The same subtraction idea may accelerate other multi-loop bound-state QED diagrams where spurious gauge-dependent terms slow convergence, not just the SVPE diagram mentioned in the paper.","The sign change near Z≈74 in the He-like ground state is a nontrivial prediction of the all-order treatment that could serve as a sensitive fingerprint for the Wichmann-Kroll contribution at moderate Z."],"forward_implications":["The new He-like values supersede the earlier numbers for the interelectronic Wichmann-Kroll contribution, replacing their error bars with uncertainties several orders of magnitude smaller.","The Li-like values, if adopted, shift predicted transition energies in lithium-like ions by a few meV relative to the earlier calculation, which is relevant for x-ray and dielectronic-recombination measurements.","The finite-nuclear-size correction, up to and beyond 5% for heavy ions, must be included in any comparison of high-Z Lamb-shift theory with experiment.","The subtraction scheme provides a template for all-order calculations of the two-loop self-energy-vacuum-polarization diagram, currently a dominant theory uncertainty."],"fun_headline_variants":["He-like QED agrees; Li-like discrepancy found","Wichmann-Kroll: heavy-ion finite-size shifts >5%","New Lamb-shift numbers for He- and Li-like ions","All-order QED loop: sign change near Z=74","Improved interelectronic QED: tighter Lamb shift"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The several-orders-of-magnitude uncertainty reduction rests on a polynomial extrapolation of the multipole tail beyond κ,κ′=±16 that no independent calculation has verified.","fun_headline_variants_meta":{"raw":{"variants":["He-like QED agrees; Li-like discrepancy found","Wichmann-Kroll: heavy-ion finite-size shifts >5%","New Lamb-shift numbers for He- and Li-like ions","All-order QED loop: sign change near Z=74","Improved interelectronic QED: tighter Lamb shift"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000184,"raw_usage":{"total_tokens":1200,"prompt_tokens":837,"completion_tokens":363,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":581,"completion_tokens_details":{"reasoning_tokens":294}},"tokens_in":581,"tokens_out":363,"duration_ms":3815,"temperature":1.0,"reasoning_tokens":294,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T06:39:57.600289+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the correction with the multipole truncation raised to κ,κ′ = ±32 (or with a different subtraction/regularization) and compare the extrapolated tail contribution to the quoted uncertainty; a deviation larger than the stated error would falsify the precision claim. For the Li-like discrepancy, an independently published recalculation that resolves the two-body values would settle which set of numbers is correct.","supporting_citations":[],"review_version":2}