{"id":"41a275b8-5a0e-4f59-be42-69bdcc1ba991","arxiv_id":"2509.07738","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Energy levels for l > 1 states of antiprotonic atoms are computed to order (Z alpha)^6 with nonperturbative vacuum polarization, giving the most accurate theoretical transition energies to date.","lead":"Physicists calculated the energy levels of atoms in which an antiproton orbits a nucleus, including vacuum polarization effects directly in the quantum wave equation. Their predictions for excited rotational states are the most precise so far and could eventually let nuclear sizes be measured from antiproton X-ray transitions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The antiproton g-factor is never specified, yet Eq. (19) and Eqs. (28)-(32) make Table I fine-structure energies critically dependent on it.","rationale":"The reader's weakest assumption correctly identified the undisclosed antiproton g-factor as the most fragile premise. My independent reading confirms that Eq. (19), Eqs. (28)-(32), and the E(4) entries of Table I all depend on g1, and the manuscript gives no numerical value. This is a genuine load-bearing concern because the fine-structure splitting between j states is directly proportional to (g1-1) in the spin-orbit term, and the antiproton's anomalous magnetic moment is large. The issue is not an internal inconsistency or a disagreement with consensus; it is a missing input parameter that prevents verification of the 'most accurate to date' claim. The reader's conditional verdict is therefore appropriate: the paper should be accepted only after this parameter is disclosed and checked. I recommend no change to the verdict because the concern, while serious, is potentially addressable by inspecting the code or asking the authors; it does not yet prove the results wrong.","tokens_in":11816,"tokens_out":10467,"duration_ms":124737,"concrete_test":"Open the Supplemental Material PbarSpectr Mathematica code and locate the definition of the antiproton g-factor (likely named g1, gPbar, or similar). If it is set to 2, recompute the Table I entries with the physical value g1≈5.585 (or −5.585 with the sign convention of Eq. 19) and verify that the transition energies shift. If it is already set to the physical value, rerun the code after temporarily changing g1 to 2 and compare: any shift larger than the quoted uncertainty means the reported values depend critically on this undisclosed parameter, and the paper must state the value used.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that Table I contains the most accurate transition energies for l>1 antiprotonic atoms depends on the value of the antiproton gyromagnetic ratio g1 used in the Breit-Pauli Hamiltonian (Eq. 19) and in the (Zα)^6 formula (Eqs. 28-32). The paper never states this value. For an antiproton, g1 is not 2: the measured magnetic moment anomaly gives g1 ≈ 5.585 (with sign convention depending on charge). The spin-orbit terms in Eq. 19 scale as (g1-1)/(2m^2), so replacing g1=5.585 with g1=2 changes the leading relativistic fine-structure contribution by a factor of about 4.6. For the 184W transition in Table I, the E(4) column is 37.26 eV, and a change of this magnitude would shift the transition energy by tens to hundreds of eV, completely swamping the stated 0.02 eV uncertainty. The text explicitly uses g=2 for E(7) (after Eq. 33) but is silent for E(4) and E(6), raising the real possibility that the default g=2 was used there. Without a stated value, the calculation is not reproducible, and the 'most accurate to date' claim cannot be checked.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents an NRQED treatment of antiprotonic atoms with spinless nuclei, aimed at circular (l>1) states where the effective expansion parameter is Zalpha/n. The authors solve the radial Schrödinger equation numerically with the Coulomb potential plus the Uehling, two-loop vacuum-polarization, and Wichmann-Kroll potentials included nonperturbatively. Using the resulting wave functions they evaluate the Breit-Pauli relativistic correction, combine this with analytic E5 and E6 results with full mass dependence, and estimate E7 and E8. Table I gives transition energies for 20Ne, 40Ar, 132Xe, and 184W and claims these are the most accurate predictions to date for l>1 antiprotonic atoms. A Mathematica code, PbarSpectr, is provided in the supplemental material.","tokens_in":12199,"tokens_out":10954,"duration_ms":136981,"significance":"The approach is significant: it extends practical NRQED calculations to high-Z two-body systems by exploiting the small parameter Zalpha/n for rotational states, and it incorporates vacuum polarization nonperturbatively while retaining exact finite-mass corrections through order (Zalpha)^6. The reported agreement with Ref. [8] at lower accuracy and the release of the code are concrete strengths. However, the central numerical claim cannot be verified from the manuscript as written because the antiproton g-factor entering the spin-dependent parts of the Hamiltonian is never specified; the same is true for the antiproton charge radius in the finite-size correction. These are not cosmetic omissions, because the tabulated transition energies are between states of definite j and therefore depend on the spin-orbit and fine-structure terms.","major_comments":[{"comment":"The numerical value of the antiproton g-factor g1 is never stated. For l>1 states with definite j, Eq. (19) contains spin-orbit terms proportional to [(g1-1)/(2m^2)+g1/(2mM)] L·s V'/r, and the E6 coefficients in Eqs. (28)-(32) also depend on g1. Table I lists transitions between different j states, so E(4) is fine-structure sensitive. The physical antiproton has |g|≈5.585, not 2; the text only says 'with g=2' in the context of the approximate E(7) estimate in Eq. (34). If g1=2 was used in E(4) and E(6), the Table I entries, e.g. E(4)=37.26 eV for 184W, would be shifted by an amount that swamps the quoted 0.02 eV uncertainty. Please state g1 and its sign convention explicitly, and confirm that Table I was computed with the physical value.","section":"II, Eq. (19) and IV, Table I"},{"comment":"E_fns in Table I includes the finite charge radii of both the nucleus and the antiproton, and the rows E_fns(fs N) are used to demonstrate nuclear charge radius determination. However, no numerical value for the antiproton charge radius r_C1 is given in the paper. This makes the 'Total (point N + fs pbar)' entries and the separation between nuclear and antiproton finite-size effects non-reproducible. The value of r_C1 and its uncertainty should be stated explicitly, preferably in a table of input constants.","section":"IV, Table I"},{"comment":"The uncertainty estimate for the omitted three-loop vacuum polarization, δE3loop ≈ (α/π)^2 [E + (Zα)^2/(2n^2)], is not defined in a way that reproduces the quoted digits in Table I. If E is the tabulated E(2) in eV, (α/π)^2 E is about 0.16 eV for 20Ne and about 5 eV for 184W, whereas the table lists uncertainties of 0.001 eV and 0.02 eV. If a different convention is intended (e.g. Hartree units without the reduced mass), it must be stated. This matters because the three-loop VP uncertainty is claimed to dominate and sets the final precision.","section":"II, Eq. (17)"}],"minor_comments":[{"comment":"Please define all symbols in the uncertainty formula; in particular, specify the units of E and the origin of the (Zα)^2/(2n^2) term.","section":"II, Eq. (17)"},{"comment":"The sums over i=−1 use notation that is easy to misread; please clarify the lower limits and the meaning of the A_r and B_r terms.","section":"II, Eq. (37)"},{"comment":"The number of decimal places is inconsistent across rows (e.g. E(5) is given to 4 decimals for Xe/W but E(7) to 3 decimals). A consistent convention would improve readability.","section":"IV, Table I"},{"comment":"The paper states that electric dipole polarizabilities are neglected, but the E6 formula in Eq. (28) includes them. A sentence quantifying the expected size of this neglected contribution would be useful.","section":"IV, Table I"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is potentially publishable, but the missing antiproton g-factor is a load-bearing omission: if the authors used g1=2 rather than the physical value, the Table I fine-structure energies would be substantially wrong and the 'most accurate to date' claim would collapse. The revision should include a complete input-constants table (g1, r_C1, nuclear radii, grid parameters) and a sensitivity check showing that the physical g1 was used. If the code already contains these values, they must be brought into the paper text; the reader should not have to open supplemental code to know the central physical inputs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper makes a genuine contribution: it extends NRQED predictions for antiprotonic atoms to l>1 rotational states, includes one- and two-loop vacuum polarization nonperturbatively in the Schrödinger equation, and treats finite nuclear mass exactly through (Zα)^6. The transition energies in Table I are new, the code is included, and the numbers agree with the less detailed results of Ref. [8] at the expected lower accuracy. The expansion parameter Zα/n for circular states is a sensible way to keep the perturbation series under control for high Z. That part deserves credit.\n\nThe soft spots are real. The most serious is that the antiproton g-factor is never stated. The Breit Hamiltonian in Eq. (19) and the (Zα)^6 coefficients in Eqs. (28)-(32) depend on g1, and for l>1 the spin-orbit terms are dominated by (g1-1)/(2m^2). The antiproton's g-factor is about 5.58, not 2. If the authors used g1=2 anywhere in E(4) or E(6), the fine-structure splittings in Table I would be off by large multiples of the quoted 10^-6 accuracy. The text explicitly uses g=2 only for the approximate E(7) estimate, which is fine. But for the main calculation no value is given. This makes the calculation irreproducible and the central 'most accurate to date' claim uncheckable. This is not a minor omission; it affects every term in the table.\n\nA second, lesser issue is that nuclear electric dipole polarizabilities are neglected without an uncertainty estimate. The paper itself notes these enter at (Zα)^6 and can be accounted for, so a simple bound or a numerical estimate should be there. A third issue is the unresolved dispute with Adkins and Jentschura over the completeness of the Breit-Pauli Hamiltonian from Ref. [6]. The authors defend their formula, and the historical verification is decent, but a short note addressing the specific 'gauge' criticism would help an independent referee.\n\nThe numerics and the structure of the expansion look sound. The code exists and the agreement with Ref. [8] is external. If the g-factor is clarified—and the table recomputed if needed—the results will be a solid reference for the PAX program. As written, the paper deserves peer review, but with a request for a major revision on the g-factor before it can be trusted.\n\nI'd bring it to a reading group focused on precision calculations, but I'd flag the g-factor at the start.\n\nRecommendation: send to peer review, but the referee should demand the g-factor value and a demonstration that the fine structure is insensitive or corrected.","headline":"Solid NRQED calculation with a genuine new dataset, but the unspecified antiproton g-factor is a load-bearing omission that must be fixed before the results can be taken at face value.","tokens_in":12650,"tokens_out":6129,"would_cite":false,"duration_ms":58271,"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":"This paper computes antiprotonic-atom transition energies to meV accuracy by solving the Schrödinger equation with vacuum polarization built in.","keywords":["antiprotonic atoms","vacuum polarization","NRQED","circular Rydberg states","finite nuclear mass","nuclear charge radius","Breit-Pauli Hamiltonian","Uehling potential"],"falsifier":"Recompute the E(4) and E(6) contributions in Table I with the antiproton g-factor set to its measured value g ≈ 5.58 rather than whatever value was silently adopted, and compare the resulting 12o to 11n transition in 184W: if the shift exceeds the quoted 0.02 eV uncertainty, the paper's accuracy claim fails. A direct measurement of the same transition at meV precision would settle the issue empirically.","tokens_in":11724,"feed_emoji":"⚛️","tokens_out":4876,"duration_ms":51135,"temperature":0.7,"pith_summary":"The paper claims that for highly excited rotational states (l>1) of an antiproton circling a spinless nucleus, nonrelativistic QED still converges even at large nuclear charge because the effective expansion parameter is Z alpha / n rather than Z alpha. The authors place the one- and two-loop electron vacuum-polarization potentials directly into the radial Schrödinger equation, solve it numerically, and use the resulting wave function to evaluate relativistic and higher-order QED corrections with exact finite-nuclear-mass dependence through order (Z alpha)^6. Their tabulated transition energies, such as 180553.71(2) eV for the 12o to 11n transition in antiprotonic 184W, are claimed to be the most accurate theoretical predictions to date for these systems. This matters because precise antiprotonic X-ray spectroscopy could then determine nuclear charge radii more accurately than electronic or muonic atoms, and because the systems test QED in a high-field, two-body setting.","feed_headline":"Antiprotonic atom energies computed to meV with vacuum polarization built in","feed_subtitle":"New NRQED calculation is the most accurate to date for l>1 states and can extract nuclear charge radii.","key_machinery":"The central object is the radial Schrödinger equation with the potential V(r) = -Z alpha/r + V_VP(r), where V_VP contains the one-loop Uehling potential, the two-loop Källén-Sabry potential, and the Wichmann-Kroll correction. Solving this equation numerically, with small-r logarithmic singularities handled by a power-log series ansatz and large-r by an asymptotic expansion, gives nonrelativistic energies and wave functions that absorb vacuum polarization nonperturbatively. Those wave functions then evaluate the Breit-Pauli Hamiltonian with vacuum polarization and the analytic (Z alpha)^5 and (Z alpha)^6 corrections. The expansion parameter Z alpha / n keeps the series convergent for circular","core_discovery":"Using NRQED, the paper demonstrates that circular Rydberg states of antiprotonic atoms can be treated nonperturbatively with respect to vacuum polarization: the Uehling, Källén-Sabry, and Wichmann-Kroll potentials are included in the Schrödinger equation, and the Breit-Pauli Hamiltonian modified by those potentials supplies the leading relativistic correction. Finite nuclear mass is included exactly to order (Z alpha)^6 through analytic formulas valid for arbitrary mass ratio. The result is a set of theoretical transition energies, computed by the accompanying PbarSpectr code, that are claimed to be the most accurate to date for l>1 states of antiprotonic atoms with a spinless nucleus, with","pith_inferences":["The PbarSpectr code could be adapted to expose the antiproton g-factor as a tunable input; scanning it would reveal how strongly the quoted fine-structure energies depend on this parameter.","Because the dominant uncertainty is estimated rather than computed, a direct numerical inclusion of the recently derived three-loop vacuum-polarization density would settle whether the stated meV accuracy holds.","The analytic (Z alpha)^6 formulas used here are valid for arbitrary constituent masses, so the same machinery may apply to other exotic two-body atoms with heavy orbiting particles, not only antiprotonic ones.","The l-dependence of transition energies gives a cross-check on the assumed suppression of strong-interaction effects: if radius extraction from different l values disagrees, hadronic corrections would need to be reintroduced."],"forward_implications":["Table I gives meV-level predictions for antiprotonic transitions in 20Ne, 40Ar, 132Xe, and 184W that upcoming X-ray experiments can test directly.","The finite-size contribution E_fns = c r_C^2 to each transition means circular-state energies directly probe mean-square nuclear charge radii.","If experimental precision reaches the meV level, extracted nuclear radii could compete with values from muonic atoms and electron scattering.","Including the three-loop vacuum polarization potential should improve the theoretical accuracy by about two orders of magnitude.","The same method extends naturally to rotational states of muonic atoms, and later to l=0,1 states once additional QED contributions are added."],"supporting_citations":[{"why":"Source of the nonperturbative strategy of putting vacuum polarization contributions directly into the Schrödinger equation.","marker":"[5]"},{"why":"Supplies the Breit-Pauli Hamiltonian including vacuum-polarization potentials, used for the leading relativistic correction.","marker":"[6]"},{"why":"Provides the analytic (Z alpha)^6 energy formulas for two-body systems with arbitrary masses and g-factors.","marker":"[7]"},{"why":"Lists the antiprotonic transitions targeted by upcoming precision measurements, the comparison set for Table I.","marker":"[8]"},{"why":"Defines the two-loop Källén-Sabry vacuum-polarization potential included in the Hamiltonian.","marker":"[12]"},{"why":"Supplies the Wichmann-Kroll correction to the vacuum-polarization potential.","marker":"[13]"},{"why":"Gives the recently derived three-loop spectral density used to estimate the dominant omitted uncertainty.","marker":"[14]"},{"why":"Provides Bethe logarithm values for n up to 20, needed for the (Z alpha)^5 correction.","marker":"[17]"},{"why":"Gives the nuclear charge radii used to compute finite-size contributions in Table I.","marker":"[24]"}],"fun_headline_variants":["Most accurate antiprotonic atom energies for high-l states","Vacuum polarization included exactly in antiprotonic atoms","Antiprotonic atoms promise better nuclear charge radii","PbarSpectr: open code for antiprotonic atom transitions","Rotational states of antiprotonic atoms probe nuclear sizes"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The fine-structure energies depend on the antiproton's magnetic moment through spin-orbit terms, but the paper never states which g-factor value was used; if it is not the real antiproton value of about 5.58, the tabulated 'most accurate' transition energies are wrong.","fun_headline_variants_meta":{"raw":{"variants":["Most accurate antiprotonic atom energies for high-l states","Vacuum polarization included exactly in antiprotonic atoms","Antiprotonic atoms promise better nuclear charge radii","PbarSpectr: open code for antiprotonic atom transitions","Rotational states of antiprotonic atoms probe nuclear sizes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00033,"raw_usage":{"total_tokens":1692,"prompt_tokens":776,"completion_tokens":916,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":831}},"tokens_in":520,"tokens_out":916,"duration_ms":10899,"temperature":1.0,"reasoning_tokens":831,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T21:41:41.931721+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the E(4) and E(6) contributions in Table I with the antiproton g-factor set to its measured value g ≈ 5.58 rather than whatever value was silently adopted, and compare the resulting 12o to 11n transition in 184W: if the shift exceeds the quoted 0.02 eV uncertainty, the paper's accuracy claim fails. A direct measurement of the same transition at meV precision would settle the issue empirically.","supporting_citations":[{"cited_title":"Borie, Vacuum polarization corrections and spin-orbit splitting in antiprotonic atoms, Phys","cited_arxiv_id":null,"evidence_quote":"Source of the nonperturbative strategy of putting vacuum polarization contributions directly into the Schrödinger equation."},{"cited_title":"Veitia and K","cited_arxiv_id":null,"evidence_quote":"Supplies the Breit-Pauli Hamiltonian including vacuum-polarization potentials, used for the leading relativistic correction."},{"cited_title":"Zatorski, V","cited_arxiv_id":null,"evidence_quote":"Provides the analytic (Z alpha)^6 energy formulas for two-body systems with arbitrary masses and g-factors."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the two-loop Källén-Sabry vacuum-polarization potential included in the Hamiltonian."},{"cited_title":"Wichmann, N.M","cited_arxiv_id":null,"evidence_quote":"Supplies the Wichmann-Kroll correction to the vacuum-polarization potential."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the recently derived three-loop spectral density used to estimate the dominant omitted uncertainty."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides Bethe logarithm values for n up to 20, needed for the (Z alpha)^5 correction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the nuclear charge radii used to compute finite-size contributions in Table I."}],"review_version":1}