{"id":"79bc4abc-cb85-43ad-9b4d-209caaed473c","arxiv_id":"2607.04815","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"NRQED with nonperturbative vacuum polarization supplies controlled, mass-exact energy levels for L>1 states of muonic, kaonic and antiprotonic atoms, with clear paths to two-order accuracy gains for nuclear-property extraction.","lead":"The paper unifies NRQED formulas for two-body rotational levels (L>1) of spin-0 or spin-1/2 particles with arbitrary masses and magnetic moments, and ships an extended PbarSpectr code that yields state-of-the-art energies for muonic, kaonic and antiprotonic atoms. These predictions, improvable by roughly two orders of magnitude once three-loop vacuum polarization is included, open precision extractions of nuclear charge radii and polarizabilities and tests of long-range hadro","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The manuscript is a careful, openly documented extension of the authors’ prior NRQED program. All formulas needed for the L>1, spin-0/½, arbitrary-mass case are collected, the code is supplied, and residual uncertainties (three-loop EVP, polarizability) are quantified and separated. The polarizability estimate is indeed crude, yet it is treated as an external nuclear-physics input rather than a theoretical pillar; removing or refining it does not alter the QED hierarchy or the claim that three-loop EVP would improve theory by ~two orders. No load-bearing technical flaw is present, so the reader’s ACCEPT / high-confidence verdict stands.","tokens_in":11043,"tokens_out":456,"duration_ms":3964,"concrete_test":"Recompute the 6h–5g centroid of ¯p 28Si (Table IV) after replacing the phenomenological α_E2 by a modern ab-initio or experimental nuclear polarizability (or by zero); if the shift remains inside the quoted (13) eV polarizability band and the first (QED) uncertainty is unchanged, the claim that residual theory error is dominated by the deferred three-loop EVP is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper’s central claim—that NRQED yields state-of-the-art, mass-exact predictions for L>1 rotational levels of muonic/kaonic/antiprotonic atoms, with residual theory error reducible by ~two orders via three-loop EVP and EVP-on-E^(5)—is supported by the explicit formulas (Eqs. 7–18), the shipped PbarSpectr code, and the tabulated breakdowns. The reader correctly flags the phenomenological polarizability (Eq. 19 + 50 % ad-hoc) as the dominant second uncertainty for heavier systems, but the manuscript already isolates it, quotes it separately, and does not rest the “state-of-the-art” or “improvable-by-two-orders” claims on its accuracy. No derivation gap, circularity, or internal inconsistency appears that would undermine the strongest claim.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript presents a unified NRQED treatment of two-body hydrogenic systems with spin-0 or spin-1/2 constituents, arbitrary masses and magnetic moments, restricted to rotational states with L>1. It collects and implements the complete set of operators through O(α^6) (with nonperturbative vacuum polarization in E^(2) and E^(4)), extends the PbarSpectr code accordingly, and supplies state-of-the-art numerical predictions for selected levels and transitions in μ^20Ne, p-bar^28,29Si, K^19F and K^20Ne. Explicit formulas for the Breit Hamiltonian, the O(α^5) recoil and Bethe-logarithm terms, the spin-dependent O(α^6) recoil corrections, and the electric-dipole polarizability contribution are given; residual uncertainties are isolated into three-loop vacuum polarization and nuclear polarizability. The authors argue that inclusion of three-loop EVP and the EVP correction to E^(5) would improve theory by roughly two orders of magnitude, enabling high-accuracy extractions of nuclear charge radii, polarizabilities and possible long-range hadronic forces.","tokens_in":11267,"tokens_out":1134,"duration_ms":8818,"significance":"The work supplies mass-exact, spin-generalized NRQED predictions for exotic atoms that are currently or soon to be measured (J-PARC muonic neon, PAX antiprotonic silicon, SIDDHARTA-type kaonic systems). The explicit operators (Eqs. 7–18), the tabulated breakdowns with two distinct uncertainty sources, and the publicly released PbarSpectr code constitute a concrete, reproducible advance over reduced-mass Dirac/Klein-Gordon calculations. If the residual theory error can indeed be reduced by the two orders claimed, precision spectroscopy of these systems becomes a competitive route to nuclear radii, polarizabilities and tests of long-range hadronic interactions. The paper therefore has clear and timely impact for both atomic and nuclear physics.","major_comments":[{"comment":"Section II, paragraph after Eq. (18) and the uncertainty discussion: the claim that three-loop EVP plus the EVP correction to E^(5) would improve accuracy by “approximately two additional orders of magnitude” is asserted but not quantified for any of the concrete systems. A short estimate (or a reference to an existing three-loop calculation for a related system) for at least one transition (e.g., 5g–4f in μ^20Ne or 6h–5g in p-bar Si) would make the central “improvable-by-two-orders” statement falsifiable rather than qualitative.","section":null},{"comment":"Eq. (19) and the accompanying 50 % uncertainty: the phenomenological polarizability formula is used for every tabulated energy and dominates the second uncertainty for the heavier systems. While the manuscript correctly isolates this contribution, the text should state more clearly that the formula is taken from Ref. [7] without re-derivation and that the 50 % figure is an ad-hoc assignment; otherwise readers may misinterpret the second error bar as a controlled theoretical uncertainty.","section":null}],"minor_comments":[{"comment":"Table I caption: the kaon mass is quoted from the 2008 PDG; a more recent value (or an explicit statement that the 2008 value is retained for consistency with earlier work) would avoid confusion.","section":null},{"comment":"Eqs. (11)–(12): the lengthy spin-dependent O(α^6) expressions are given only for the NS component; a brief pointer that the remaining spin-orbit, spin-spin and tensor pieces are taken unchanged from Ref. [6] would help readers who do not have that paper open.","section":null},{"comment":"Section III.C (K^19F): the strong-interaction shift of −2 eV for the 4f–3d transition is mentioned only in the text and not in Table V; adding a footnote or a separate column would make the comparison with experiment transparent.","section":null},{"comment":"Typographical: “EXP ANSION” in the section heading II; “RESUL TS” in heading III; and the inconsistent use of “Wichman-Kroll” versus the more common “Wichmann-Kroll”.","section":null},{"comment":"The Supplemental Material is cited as [4] but the arXiv version does not yet contain a permanent link or DOI; a stable repository identifier would improve long-term reproducibility.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid, incremental but useful extension of the authors’ earlier NRQED work. The two major points I raise are easily addressable and do not threaten the central claim. I see no reason for a second full review cycle once the polarizability uncertainty is clarified and a quantitative estimate of the three-loop improvement is supplied."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a clean, practical extension of the authors’ NRQED program. They take the existing mass-exact operators through α^6 (their earlier papers plus the Breit Hamiltonian with finite-size and g-factors), put them into an updated PbarSpectr code that now handles spin-0 or 1/2 on either particle with arbitrary magnetic moments, and produce concrete numbers for the systems people are actually measuring: μ20Ne, p-bar Si isotopes, K19F, K20Ne.\n\nWhat is new is the unified implementation and the tables. The Schrödinger solution includes Uehling + two-loop EVP non-perturbatively in both E(2) and E(4); every contribution is broken out with two separate uncertainty sources (three-loop EVP estimated conservatively as α^{2}V(1), and nuclear polarizability). The isotope shift for p-bar 29Si–28Si is a nice, clean test case. Comparisons with earlier Dirac/Klein-Gordon-plus-reduced-mass work show the expected small but real mass-ratio differences once relativistic corrections matter. The code is shipped as supplemental material, so the results are reproducible.\n\nThe soft spot is exactly what the reader flagged: nuclear static polarizability is taken from a phenomenological formula with an ad-hoc 50 % error bar. That term dominates the second uncertainty for the heavier systems. They quote it separately and do not hide it, so it does not undermine the “state-of-the-art and improvable by two orders with three-loop EVP” claim. E(7) for spin-0 is borrowed from the spin-1/2 result; that is a minor approximation they acknowledge. Strong-interaction shifts for the lowest kaonic levels are left out, again stated.\n\nCitation pattern is normal for this group; the self-cites are to the operator papers that are independently checkable. No free parameters are being sold as predictions.\n\nThis is for people doing precision exotic-atom spectroscopy or extracting nuclear radii/polarizabilities. It deserves a serious referee. I would accept it, use the tables, and watch for the three-loop EVP update.","headline":"Solid, usable NRQED tables and open code for L>1 exotic atoms; the polarizability estimate is the only soft spot and they already isolate it.","tokens_in":11861,"tokens_out":539,"would_cite":true,"duration_ms":4654,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A unified NRQED treatment of rotational levels for spin-0 and spin-1/2 exotic atoms yields state-of-the-art energies and a clear path to two-order theory gains.","keywords":["NRQED","exotic atoms","muonic atoms","kaonic atoms","antiprotonic atoms","vacuum polarization","nuclear polarizability","rotational levels"],"falsifier":"Measure any of the listed transitions (for example the 6h–5g line in antiprotonic silicon or the 5g–4f line in muonic neon) at the few-meV level and compare with the tabulated NRQED prediction after the three-loop vacuum-polarization term has been added; a persistent discrepancy larger than the remaining nuclear-polarizability uncertainty would falsify the claimed accuracy.","tokens_in":11959,"feed_emoji":"⚛️","tokens_out":746,"duration_ms":5734,"temperature":0.7,"pith_summary":"Exotic atoms that bind a muon, kaon or antiproton to a light nucleus sit outside the reach of the ordinary Dirac or Klein–Gordon equations once the two masses become comparable. This paper collects and implements the complete set of nonrelativistic-QED formulas that work for any mass ratio and for either spin-0 or spin-1/2 constituents, restricted only to states with orbital angular momentum L greater than 1. The resulting code supplies the most accurate theoretical energies currently available for the transitions that experiments are already measuring or planning. The authors further show that two still-missing pieces—three-loop vacuum polarization and the vacuum-polarization correction to the leading recoil term—would shrink the residual theory error by roughly two orders of magnitude. With that improvement, precision spectroscopy of these systems could extract nuclear charge radii and electric-dipole polarizabilities at high accuracy and could place tight limits on long-range hadronic forces beyond the Standard Model.","feed_headline":"NRQED predicts exotic-atom levels and shows a path to 100\times tighter theory","feed_subtitle":"Muonic, kaonic and antiprotonic spectroscopy can then extract nuclear radii and polarizabilities cleanly","key_machinery":"The NRQED power series for the energy, E = E² + E⁴ + E⁵ + E⁶ + …, evaluated with non-perturbative inclusion of the Uehling and two-loop vacuum-polarization potentials inside the Schrödinger equation and the Breit Hamiltonian, and with exact mass-ratio dependence retained through order α⁶.","core_discovery":"Nonrelativistic quantum electrodynamics, applied uniformly to two-body systems of spin-0 or spin-1/2 particles with arbitrary masses and magnetic moments, already furnishes state-of-the-art predictions for the rotational (L>1) levels of muonic, kaonic and antiprotonic atoms; the same framework can be tightened by two further orders of magnitude once three-loop vacuum polarization and the vacuum-polarization correction to the O(α⁵) recoil term are included.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["NRQED unifies spin-0/1/2 exotic atoms for L>1 rotational level predictions","State-of-the-art NRQED spectra for muonic kaonic and antiprotonic atoms","NRQED framework yields precision L>1 levels in two-body exotic systems","Path to tighter theory for nuclear radii from exotic-atom NRQED levels","Unified NRQED of arbitrary-mass spin-0 or 1/2 particles in rotational states"],"cache_read_input_tokens":128,"weakest_assumption_plain":"Nuclear electric-dipole polarizabilities are taken from a simple phenomenological formula that carries an arbitrary 50 percent uncertainty; that estimate currently dominates the second error bar on every tabulated energy for the heavier systems.","fun_headline_variants_meta":{"raw":{"variants":["NRQED unifies spin-0/1/2 exotic atoms for L>1 rotational level predictions","State-of-the-art NRQED spectra for muonic kaonic and antiprotonic atoms","NRQED framework yields precision L>1 levels in two-body exotic systems","Path to tighter theory for nuclear radii from exotic-atom NRQED levels","Unified NRQED of arbitrary-mass spin-0 or 1/2 particles in rotational states"]},"model":"grok-4.5","effort":"low","cost_usd":0.006436,"raw_usage":{"total_tokens":1568,"prompt_tokens":695,"num_sources_used":0,"completion_tokens":121,"cost_in_usd_ticks":64360000,"prompt_tokens_details":{"text_tokens":695,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":752,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":695,"tokens_out":121,"duration_ms":5834,"temperature":1.0,"reasoning_tokens":752,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T13:05:54.088830+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Measure any of the listed transitions (for example the 6h–5g line in antiprotonic silicon or the 5g–4f line in muonic neon) at the few-meV level and compare with the tabulated NRQED prediction after the three-loop vacuum-polarization term has been added; a persistent discrepancy larger than the remaining nuclear-polarizability uncertainty would falsify the claimed accuracy.","supporting_citations":[],"review_version":1}