{"id":"b268b182-e6cc-40c4-8e9a-499eb32ac2ea","arxiv_id":"2509.01303","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Nuclear polarizability, enhanced in odd-odd nuclei by spin-isospin symmetry, explains the Zemach radius discrepancy between effective and elastic values in 6Li, 7Li, 2H, and 3He.","lead":"Using neural-network quantum Monte Carlo with chiral nuclear forces, this paper computes the magnetic-size radii of lithium-6 and lithium-7 and attributes the observed discrepancy in lithium-6 to nuclear polarizability from virtual excitations. The result offers a common symmetry-based explanation for isotope-dependent hyperfine corrections in light atoms, relevant for extracting nuclear properties from precision spectroscopy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Closure-energy fit to the same HFS data makes the absolute polarizability-dominance claim for 6Li non-predictive; an independent ωbar is required.","rationale":"The reader's weakest_assumption correctly identifies the closure approximation with per-nucleus fitted ωbar as the load-bearing element. My independent reading of the manuscript reaches the same conclusion: the central quantitative claim—that polarizability dominates in 6Li—rests on the closure approximation, and the only evidence for that approximation in 6Li is that adjusting ωbar makes the experimental reff_Z match. This is circular in the sense that the parameter is chosen to force agreement, while no uncertainty quantifies the closure error. The elastic Zemach radii rZ are a genuinely independent result and agree with experiment, and the SU(4) ratio pattern in Fig. 3 is a useful parametric trend that is robust to the choice of ωbar; these provide real partial support for the paper's message. However, the abstraction to 'ab initio' effective Zemach radii is weakened by the fit. The 3He residual is also an important internal indicator: the same closure framework cannot reproduce 3He within 6% even when ωbar is varied freely, showing that missing physics (two-body currents, relativistic corrections) is present at the same order as the claimed polarizability effect. A concrete, feasible test—computing ωbar from the exact excitation spectrum for 6Li—would settle whether the closure approximation is quantitatively reliable. In the absence of that test, the reader's CONDITIONAL verdict is appropriate, so I recommend no change to the verdict.","tokens_in":12172,"tokens_out":3256,"duration_ms":44682,"concrete_test":"Compute the closure energy for 6Li from the actual excitation spectrum rather than from HFS: with the same N2LO(1.0) Hamiltonian and the VMC-NN ground state, evaluate the energy-weighted average ωbar = <0|[∇·j, ((r−r')×j)_z] (H−E0)|0> / <0|[∇·j, ((r−r')×j)_z]|0> using a complete set of excited states generated by, e.g., the Lorentz Integral Transform or stochastic Lanczos. Then recompute reff_Z = rZ + recoil + single-nucleon TPE + δrpol(ωbar_ab initio) without fitting. If the predicted reff_Z deviates from the HFS-extracted value by more than the combined uncertainties (roughly 0.05 fm), the closure approximation is not validated and the polarizability-dominance claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that nuclear polarizability dominates in 6Li and explains the discrepancy between effective and elastic Zemach radii. The absolute size of δrpol is computed in the closure approximation (Eqs. 9–11) using a single state-independent excitation energy ωbar. In Fig. 4, the experimental reff_Z are reproduced by adjusting ωbar to 23 MeV for 6Li and 29 MeV for 7Li. That adjustment is a two-parameter fit to exactly the data the claim is supposed to explain: the agreement is built in, and no closure-error estimate is attached. The only independent anchor is for 2H, where ωbar=14.4(4) MeV from pionless EFT matches the HFS-determined 16 MeV; no such anchor exists for 6Li/7Li. The SU(4) ratios in Fig. 3 are much less sensitive to ωbar and support the relative trend, but they do not establish the absolute statement that polarizability—rather than density-shape or omitted two-body currents/relativistic corrections—dominates in 6Li. The acknowledged 3He residual (about 6%, not removable by varying ωbar because δrpol is always negative) is direct evidence of systematic omissions in the framework; their size for 6Li is unquantified. Without an independent determination of ωbar, the headline quantitative agreement is not a prediction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents ab initio variational Monte Carlo calculations of elastic and effective Zemach radii for 6,7Li using a novel neural-network wave function, together with a closure-approximation treatment of nuclear polarizability contributions to hyperfine splitting. The elastic radii from realistic densities are found to closely match the values obtained from Gaussian-density assumptions, ruling out density-shape effects as the source of the reffZ/rZ discrepancy in 6Li. The central claim is that nuclear polarizability, computed via a state-independent closure energy ωbar, dominates in 6Li and accounts for the observed difference between effective and elastic Zemach radii, while being negligible in 7Li; the same mechanism is argued to explain trends in 2H and 3He, with an SU(4) enhancement factor for odd-odd nuclei. The manuscript includes validation of elastic radii against experiments and a cross-check of the closure energy for 2H against pionless EFT.","tokens_in":12533,"tokens_out":4466,"duration_ms":58012,"significance":"If the central claim holds, this is a notable step: it would establish nuclear polarizability as the dominant source of the long-standing 6Li effective-vs-elastic Zemach radius discrepancy, and it provides a concrete SU(4)-based pattern for light nuclei that is testable in other systems. The elastic Zemach radii themselves are a solid result: the VMC-NN densities are shown to reproduce experimental radii, and the direct comparison with Gaussian shapes is a useful resolution of a debated point. The closure-approximation expression relates polarizability to ground-state observables, and the 2H cross-check against pionless EFT gives some independent support for the method. The SU(4) ratio analysis is largely insensitive to the chosen ωbar and provides a robust qualitative explanation of the isotope trend. However, the absolute polarizability contribution in 6Li/7Li is not predicted, because the closure energies are fitted to the same HFS data the claim is supposed to explain; this undermines the headline statement and needs to be addressed before the paper can be accepted.","major_comments":[{"comment":"The absolute polarizability contribution in 6Li and 7Li is not predicted: the closure energies ωbar=23 MeV and 29 MeV are adjusted to reproduce the experimental reffZ values. The agreement in Fig. 4 is therefore partly built in and cannot by itself establish that nuclear polarizability dominates the reffZ−rZ discrepancy in 6Li. The 2H cross-check (ωbar=14.4(4) MeV from pionless EFT versus 16 MeV from HFS) is valuable but does not validate the closure energy for A=6,7. Please provide an independent determination of ωbar (e.g., from the computed excitation spectrum or a sum-rule/energy-weighting calculation) or a quantified closure-error band, and state which conclusions are independent of the fitting.","section":"Eqs. (9)–(11) and Fig. 4"},{"comment":"The acknowledged ~6% deviation for 3He, which cannot be removed by varying ωbar because δrpol is always negative while experimental reffZ>rZ, signals missing contributions (the authors suggest two-body currents and relativistic corrections). Since the same formalism is used for lithium, the potential size of these omitted effects in 6Li is unquantified. This does not invalidate the SU(4) ratio trend, but it undercuts the absolute claim that polarizability, rather than omitted currents/relativistic effects, is responsible for the 6Li discrepancy. Please quantify or bound the analogous residual for 6Li, or temper the absolute claim accordingly.","section":"Fig. 4, lower panel (3He)"},{"comment":"The SU(4)-limit formula predicts δrpol=0 for even-Z odd-N nuclei, which would make the 3He polarizability contribution zero, yet the manuscript treats 3He with a nonzero fitted ωbar=25 MeV and describes the observed reffZ−rZ discrepancy for 3He as smaller than for 6Li/2H but not zero. Please clarify how the full VMC calculation produces a nonzero 3He polarizability and how Fig. 3 represents 3He relative to the SU(4) prediction. As written, Eq. (12) appears inconsistent with the 3He discussion.","section":"Eq. (12) and Fig. 3"}],"minor_comments":[{"comment":"The text 'using (H−E0)/ωN ≡ 1' is informal: this is an identity only when acting on excited states with excitation energy ωN. Please reformulate to avoid implying an operator identity in the ground-state subspace.","section":"Eq. (9)"},{"comment":"The derivation of the coordinate-space kernel f(ω,r) and the impulse-approximation current operators is relegated to the Supplemental Material. For a Letter, at least the definitions of f and the current operator conventions should be stated in the main text or in an appendix.","section":"Eqs. (10)–(11)"},{"comment":"The 'experimental' δrpol values are obtained by subtracting recoil and single-nucleon TPE contributions from reffZ−rZ, but the magnitudes and uncertainties of those subtraction terms are not shown. Please display them or list them so the error bars on the experimental δrpol can be assessed.","section":"Fig. 3, right panel"},{"comment":"The statement that the chosen ωbar values 'lie close to known resonances' is too vague. Provide a quantitative comparison (e.g., excitation-energy centroids) and explain why proximity to a resonance validates a state-independent closure approximation.","section":"p. 4, closure-energy justification"}],"recommendation":"major_revision","confidential_remarks":"The paper has strong components: the elastic Zemach radii and the Gaussian-shape test are solid, and the SU(4) ratio pattern is a useful contribution. The blocking issue is the fitted closure energy for 6Li/7Li; unless the authors supply an independent anchor or reframe the claim as a consistency check within the closure model, the abstract's 'dominate' assertion is not supported. The 3He residual reinforces the need for a closure-error or missing-currents estimate. I would be willing to review a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth a serious referee. The genuinely new pieces are the first ab initio Zemach radii for 6,7Li and a compact SU(4) ratio rule that organizes the isotope pattern across 2H, 3He, 6Li, and 7Li. The elastic radii look solid, the neural-network VMC benchmarks well against GFMC and AFDMC, and the direct check that Gaussian densities are not the source of the 6Li discrepancy is a real step forward.\n\nThe SU(4) expression in Eq. (12) is elegant and the ratios in Fig. 3 are insensitive to the closure energy, which gives me confidence in the relative trends. The 2H anchor from pionless EFT (omega_bar=14.4 MeV close to the HFS-derived 16 MeV) is a meaningful independent check that the closure approximation is not crazy. For that part, the authors earn their claims.\n\nThe soft spot is the absolute scale of the polarizability correction. For 6Li and 7Li, omega_bar is adjusted to reproduce the experimental effective Zemach radii—23 and 29 MeV, respectively. So the agreement in Fig. 4 is a two-parameter fit to the very data the headline claim is supposed to explain. The 2H anchor is good but it is one nucleus, not a substitute for an independent omega_bar in lithium. The acknowledged ~6% residual in 3He, which cannot be removed by varying omega_bar because delta_rpol is always negative, shows that omitted two-body currents and relativistic effects are not negligible; their size in 6Li is unquantified. The authors are transparent about the fitting, which I appreciate, but the abstract's phrasing overstates what has actually been demonstrated.\n\nWho is this for? People at the atomic/nuclear interface working on hyperfine structure, Zemach moments, and ab initio methods for light nuclei. It deserves peer review, not desk rejection. A good referee should ask for an explicit excited-state calculation or an independent estimate of omega_bar for 6,7Li, or at least a clear statement that the absolute dominance claim is contingent on the closure parametrization. With that caveat addressed, it is publishable; even without, the elastic radii and SU(4) systematics justify publication after revision.\n\nSend it to peer review, with a request to either provide an independent omega_bar or soften the absolute claim.","headline":"First ab initio 6,7Li Zemach radii and a neat SU(4) ratio rule, but the headline 6Li polarizability dominance rests on a closure energy fitted to the data it explains.","tokens_in":13008,"tokens_out":2924,"would_cite":true,"duration_ms":35925,"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":"Nuclear polarizability, not nuclear density shape, explains the 6Li Zemach-radius discrepancy.","keywords":["Zemach radius","hyperfine splitting","nuclear polarizability","lithium isotopes","chiral effective field theory","neural-network variational Monte Carlo","SU(4) symmetry","closure approximation"],"falsifier":"Perform an explicit sum over excited nuclear states, rather than the closure average, for 6Li and 7Li using the same chiral Hamiltonians and currents, and compare the resulting polarizability correction with the fitted-closure values (ωbar ≈ 23 and 29 MeV). If the explicit correction is much smaller, the claim that polarizability dominates in 6Li would fail.","tokens_in":12073,"feed_emoji":"⚛️","tokens_out":8229,"duration_ms":87311,"temperature":0.7,"pith_summary":"The paper sets out to explain why the effective Zemach radius of 6Li extracted from hyperfine splitting is about 40% smaller than the elastic Zemach radius built from measured charge and magnetic radii, even though the two agree in 7Li. It argues that virtual nuclear excitations — nuclear polarizability — are negligible in 7Li but dominate in 6Li, and that this same mechanism accounts for similar discrepancies in 2H and 3He. To make the case, the authors compute both radii from first principles using a new neural-network variational Monte Carlo method with chiral nuclear forces, and they derive a compact closure-approximation expression that turns the polarizability correction into a ground-state expectation value. If the argument holds, isotope-dependent nuclear corrections to hyperfine splitting in light atoms are no longer mysterious: polarizability is the dominant and largely predictable source.","feed_headline":"Nuclear polarizability explains 6Li's Zemach radius anomaly","feed_subtitle":"Ab initio calculation shows virtual nuclear excitations dominate the 40% gap between effective and elastic radii in 6Li.","key_machinery":"The load-bearing object is the closure-approximation polarizability operator δrpol(ωbar), which rewrites the two-photon-exchange current using current conservation and replaces every virtual excitation energy ωN with a single state-independent average ωbar, turning δrpol into a ground-state expectation value built from single-nucleon charge and magnetic densities. On top of this, approximate SU(4) symmetry reduces the correction to a counting rule: only unpaired nucleons contribute, so δrpol vanishes for even-Z, odd-N nuclei, scales as 1/Z for odd-Z nuclei, and is enhanced by µp/(µp+µn)≈3 in odd-odd nuclei. The neural-network variational Monte Carlo wave function provides the ground states w","core_discovery":"The central claim is that the 6Li mismatch between reffZ (from hyperfine splitting) and rZ (from elastic nuclear densities) is caused by nuclear polarizability, not by the Gaussian-shape assumption used to convert charge and magnetic rms radii into rZ. Using ab initio wave functions, the paper computes rZ directly from realistic densities and shows it agrees with the Gaussian-based elastic radius, ruling out the density-shape explanation. It then evaluates the polarizability correction δrpol in the closure approximation and demonstrates that it is large in 6Li and small in 7Li; adjusting a single closure energy reproduces the experimental effective Zemach radii. The paper further claims that","pith_inferences":["Because the fitted closure energies (23 and 29 MeV for 6Li and 7Li) sit near known resonances, a testable extension is that ωbar tracks the energy-weighted centroid of the magnetic-dipole excitation spectrum, allowing it to be fixed by photoabsorption or inelastic-scattering data instead of by the hyperfine measurements themselves.","The residual 6% deviation in 3He, where the calculated reffZ exceeds the measured value while δrpol is always negative in closure, suggests a positive missing term of comparable size; computing two-body currents and relativistic corrections would provide a sharp test of the closure approximation at the few-percent level.","If the SU(4) ratio rule extends to beryllium, the effective Zemach radius of 9Be (odd-odd) should deviate from its elastic value in the opposite sense relative to 7Be (odd-even), a signature that precision Be+ hyperfine spectroscopy could look for."],"forward_implications":["The discrepancy between effective and elastic Zemach radii in 6Li is resolved as a nuclear-polarizability effect, so the Gaussian density assumption in earlier extractions is not the source of the anomaly.","Hyperfine-structure measurements in light atoms can be used as probes of nuclear excitation spectra, with polarizability corrections estimated from ground-state wave functions alone through the closure formula.","The SU(4) ratio rule gives a predictive grid: among light isotopes, odd-odd nuclei should show roughly three times larger polarizability corrections than odd-even ones, with the correction decreasing as proton number grows.","The same unified framework reproduces the effective Zemach radii of 2H and 3He along with 6Li and 7Li, connecting previously separate anomalies in deuterium and helium to the lithium puzzle."],"supporting_citations":[{"why":"Defines the Zemach radius rZ as the density convolution that serves as the elastic baseline.","marker":"[8]"},{"why":"Provides the 6,7Li atomic hyperfine measurements from which the effective Zemach radius reffZ is extracted.","marker":"[9]"},{"why":"Precision theory of 6,7Li+ hyperfine splitting used for an independent extraction of reffZ from ion data.","marker":"[10]"},{"why":"Recent 6Li+ optical Ramsey measurement that fixes the experimental effective Zemach radius compared in Fig. 4.","marker":"[11]"},{"why":"Gives the Gaussian-density estimate REM from charge and magnetic rms radii, the elastic value that disagrees with reffZ in 6Li.","marker":"[12]"},{"why":"Previous ab initio nuclear-polarizability calculation for deuterium, supplying the independent closure-energy benchmark around 14.4 MeV.","marker":"[21]"},{"why":"Earlier Low-term treatment of polarizability for A≤3 that the closure formalism is designed to improve upon.","marker":"[23]"},{"why":"Supplies the local chiral N2LO two- and three-nucleon interactions used in the VMC-NN Hamiltonians.","marker":"[40]"},{"why":"Provides the single-nucleon charge and magnetic form factors entering the impulse-approximation densities and currents in Eq. (10).","marker":"[48]"}],"fun_headline_variants":["6Li hyperfine anomaly traced to nuclear polarizability","Virtual nuclear excitations explain 6Li Zemach radius gap","Polarizability, not shape, explains 6Li Zemach gap","Nuclear polarizability dominates 6Li's Zemach radius deviation","Why 6Li's Zemach radius deviates: polarizability"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The calculation replaces all the many possible excitation energies of the nucleus with one single average number, and that number is chosen to reproduce the very hyperfine measurements being explained.","fun_headline_variants_meta":{"raw":{"variants":["6Li hyperfine anomaly traced to nuclear polarizability","Virtual nuclear excitations explain 6Li Zemach radius gap","Polarizability, not shape, explains 6Li Zemach gap","Nuclear polarizability dominates 6Li's Zemach radius deviation","Why 6Li's Zemach radius deviates: polarizability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000545,"raw_usage":{"total_tokens":2457,"prompt_tokens":768,"completion_tokens":1689,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":512,"completion_tokens_details":{"reasoning_tokens":1602}},"tokens_in":512,"tokens_out":1689,"duration_ms":15083,"temperature":1.0,"reasoning_tokens":1602,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T12:39:57.181591+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform an explicit sum over excited nuclear states, rather than the closure average, for 6Li and 7Li using the same chiral Hamiltonians and currents, and compare the resulting polarizability correction with the fitted-closure values (ωbar ≈ 23 and 29 MeV). If the explicit correction is much smaller, the claim that polarizability dominates in 6Li would fail.","supporting_citations":[{"cited_title":"To go beyond 2H, cer- tain approximations have to be adopted to make numer- ical calculations tractable","cited_arxiv_id":null,"evidence_quote":"Defines the Zemach radius rZ as the density convolution that serves as the elastic baseline."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the 6,7Li atomic hyperfine measurements from which the effective Zemach radius reffZ is extracted."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Recent 6Li+ optical Ramsey measurement that fixes the experimental effective Zemach radius compared in Fig. 4."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Gaussian-density estimate REM from charge and magnetic rms radii, the elastic value that disagrees with reffZ in 6Li."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Previous ab initio nuclear-polarizability calculation for deuterium, supplying the independent closure-energy benchmark around 14.4 MeV."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier Low-term treatment of polarizability for A≤3 that the closure formalism is designed to improve upon."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the single-nucleon charge and magnetic form factors entering the impulse-approximation densities and currents in Eq. (10)."}],"review_version":1}