{"id":"f19b78ce-87e9-4dee-a033-9163f215d781","arxiv_id":"2509.07310","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The first quantum Monte Carlo evaluation of the nuclear-structure-dependent radiative correction in carbon-10 confirms the NCSM dispersion result, with the residual uncertainty set by two undetermined low-energy constants.","lead":"Physicists simulated the nucleus carbon-10 with quantum Monte Carlo methods to compute a tiny electromagnetic correction in its beta decay, and found the correction agrees with two earlier independent methods. The correction feeds the most precise test of the Standard Model's CKM unitarity rule, which currently shows a possible hint of new physics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed 'good agreement' with prior δ_NS evaluations rests on the arbitrary g̃=±1 LEC band: LEC-free GFMC central values (−3.3 to −3.7×10⁻³) sit ~2–3σ above the NCSM result −4.22(32)×10⁻³, and only the ±(0.5–0.8)×10⁻³ LEC prior converts this into agreement.","rationale":"The reader's weakest_assumption identifies the same load-bearing point I find: the headline 'good agreement' is obtained only after adding the ±1 LEC band. My read sharpens it in two ways. First, the quantitative deficit of the LEC-free part: combining Eq. (34) with the O(α²) term (Eq. 37) and δ_E (Eq. 39) gives no-LEC δ_NS central values −3.29 to −3.67×10⁻³, which are 1.7–2.9σ above the NCSM central value −4.22(32)×10⁻³ (Eq. 36), and 0.3–0.7×10⁻³ from HT's −4.0(5)×10⁻³. The LEC band ±(0.48–0.77)×10⁻³ (Eq. A14) is what renders these 'compatible within error.' Because the LEC signs are unknown, the band only widens the error; it never shifts the central values toward the comparison targets. The abstract's 'good agreement' is therefore a consistency check with a dimensional-analysis prior, not a parameter-free outcome. Second, an internal tension: Sec. V.A's two-nucleon-amplitude check (<2% agreement between AV18 and NV2-Ia) argues for small contact terms, while Table V's ±1 naturalness assignment allows contact matrix elements ~10% of the magnetic ones. The two statements can coexist, but they point in opposite directions on whether the LEC band is honest or conservative.\n\nI checked the other candidate concerns and none is more load-bearing. The four-Hamiltonian bracket, the M_F = √2(Z−1) Fermi check (Table II), the explicit corrections of the LS sign and L_CM term relative to Ref. [48], the stated omission of the three-body O(α²) term (Eq. A18), and the acknowledged regulator inconsistency for AV18 contacts are all transparent limitations, not hidden defects. The GFMC mixed-estimate error is propagated in τ and is second-order in the trial-state correction; it is secondary to the LEC effect. The V_ud part of the abstract is weak but not false: the ¹⁰C-only extraction has a 66×10⁻⁵ experimental error that dominates, so compatibility there is nearly guaranteed. The regulator dependence of the magnetic matrix elements (noted in Sec. V.A) actually reinforces the LEC concern: the long-range/contact split is scheme-dependent, so even the paper's 3σ comparison of the bare long-range part to δ_NS,B is meaningful only after the LEC scheme is fixed.\n\nThe verdict stays CONDITIONAL (no change from the reader): the computation is careful, reproducible in structure, and honestly labeled, but the central claim of agreement is not yet a prediction. The concrete test — a two-isotope fit of the LECs using ¹⁰C and ¹⁴O, or a lattice/matching determination from two-nucleon amplitudes — would settle whether the EFT+QMC machinery genuinely agrees with Refs. [45,50] or sits ~3σ lower.","tokens_in":23434,"tokens_out":22066,"duration_ms":237235,"concrete_test":"Two-isotope LEC determination: g_V1 and g_V2 in Eq. (A10) are process-independent, and the same EFT formalism was already applied to ¹⁴O in Ref. [48]. Treat V_ud as known (neutron decay, or the global superallowed fit), solve for the two LECs using the measured ft values of ¹⁰C and ¹⁴O, recompute the LEC column of Table V with the fitted values, and re-test agreement with Eqs. (35)-(36). If the fitted |g̃| deviate strongly from 1, or if δ_NS shifts away from the HT/NCSM central values, the headline agreement fails; if they land near ±1 and agreement persists, it survives. Cross-check: match the two-nucleon weak scattering amplitudes pp(¹S₀)→np(¹S₀)e+ν and pn(¹S₀)→nn(¹S₀)e+ν to lattice-QCD two-nucleon matrix elements (Refs. [89,90]).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's central claim — 'good agreement' with both the Hardy-Towner extraction and the NCSM/dispersive evaluation (Ref. [50]) — is carried almost entirely by the arbitrary LEC band, not by the calculation.\n\nThe LEC-independent content is fixed by the paper's own numbers. The mag+LS part of δ_NS^(0) is −[4.06,4.43]×10⁻³ (Eq. 34), which the paper notes is ~40% larger than, and ~3σ away from, HT's δ_NS,B = −3.06(35)×10⁻³ (Eq. 32); the deviation is attributed to GFMC correlations. Adding the O(α²) piece (−[0.21,0.40]×10⁻³, Eq. 37) and δ_E (+[0.97,1.17]×10⁻³, Eq. 39) gives no-LEC δ_NS central values of −3.29 to −3.67×10⁻³ across the four interactions. These sit 0.55–0.93×10⁻³ less negative than the NCSM value −4.22(32)×10⁻³ (Eq. 36) — 1.7–2.9σ of the NCSM error alone, ~2σ even if the model spread is counted — and 0.3–0.7×10⁻³ from HT's −4.0(5)×10⁻³ (Eq. 35). The ±(0.48–0.77)×10⁻³ LEC band (Eq. 33, Eq. A14 with g̃=±1) is what converts this into 'compatible within error.' Since the LEC signs are unknown, that band never moves the central values toward the comparison targets; it only inflates the error until overlap occurs. The agreement is therefore an overlap between a dimensional-analysis prior and an external result, not a prediction. That is underdetermination rather than inconsistency, and the paper is honest about it — but the abstract's 'good agreement' goes beyond what the evidence supports.\n\nA second, internal tension strengthens the point: Sec. V.A reports that the two-nucleon weak amplitudes computed with AV18 and NV2-Ia agree to better than 2%, arguing that contact contributions are unlikely to explain the model spread (small LECs), while Table V's ±1 naturalness assignment allows contact matrix elements ~10% of the magnetic ones. Which regime is realized decides whether this approach agrees with Ref. [50] at the ~1σ or ~3σ level, and nothing in the paper determines it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents the first QMC calculation of the nuclear-structure-dependent radiative correction \\bar{\\delta}_{NS} for the superallowed 10C→10B* decay, using the EFT formalism of Refs. [48,49]. The authors evaluate the relevant two-body operator matrix elements (magnetic, tensor, spin-orbit, contact, and O(α^2) Fermi) with VMC and GFMC wave functions generated from four Hamiltonians: AV18+UX, AV18+IL7, NV2+3-Ia, and NV2+3-Ia*. They benchmark the Fermi matrix element (M_F = 7.07), analyze the transition densities, and study the role of OPE correlations. The GFMC long-range magnetic-plus-spin-orbit contribution is -[4.06,4.43]×10^{-3} (Eq. (34)); the O(α^2) two-body contribution is -[0.21,0.40]×10^{-3} (Eq. (37)); and the energy-dependent term is +[0.97,1.17]×10^{-3} (Eq. (39)). Adding an arbitrary dimensional-analysis range for the unknown contact LECs (Eq. (A14) with \\tilde{g}=±1) makes the total \\bar{\\delta}_{NS} compatible with the Hardy-Towner and NCSM values, and leads to the V_ud extractions in Eqs. (46)–(49).","tokens_in":23998,"tokens_out":8869,"duration_ms":104582,"significance":"If the underlying EFT derivation is accepted, this is a useful independent many-body evaluation of \\bar{\\delta}_{NS}. The paper is careful in its operator definitions, provides a transparent decomposition of uncertainties, corrects a sign in the spin-orbit operator, restores the L_CM term, and includes a simple Fermi-matrix-element sanity check. The radial-density analysis and the correlation study in Sec. V are valuable diagnostics. The main limitation is structural rather than mathematical: the two contact LECs are undetermined, and the quoted agreement with previous evaluations is obtained by adding an arbitrary LEC band to the calculated central values. The authors are transparent about this, but the abstract's \"good agreement\" language goes beyond what the LEC-free calculation supports. The work is nonetheless a meaningful first step in benchmarking QMC against existing shell-model and dispersive results, provided the claims are recalibrated.","major_comments":[{"comment":"The claim of \"good agreement\" is not supported by the LEC-independent part of the calculation. Summing the GFMC entries in Table V (mag+LS, O(α²), and δ_E^NS) without the contact-LEC term gives δ_NS central values of approximately −3.3 to −3.7×10^{-3} across the four interactions. These lie 0.55–0.93×10^{-3} above the NCSM value −4.22(32)×10^{-3} (Eq. (36)), i.e., about 2σ, and 0.3–0.7×10^{-3} from the Hardy–Towner value −4.0(5)×10^{-3} (Eq. (35)). The ±(0.48–0.77)×10^{-3} band from Eq. (A14) with \\tilde{g}=±1 is an arbitrary dimensional-analysis prior; it only widens the error bars and does not move the central values toward the comparison targets. The abstract and the conclusion (\"very good agreement\" with Ref. [50]) should be revised to state that the results are compatible only after including this prior, and the LEC-free comparison should be displayed separately.","section":"§VI, Table V, Abstract"},{"comment":"The O(α^2) three-body transition operator is not evaluated. Because the O(α^2) two-body contribution is already −[0.21,0.40]×10^{-3} (Eq. (37)), an uncomputed three-body term of comparable size would enter at the quoted precision. The manuscript should either compute/estimate this contribution or explicitly identify it as a missing piece in the uncertainty budget. As written, \"which we have not evaluated\" leaves a known omission that is not reflected in the error bars.","section":"Appendix A, Eq. (A18)"}],"minor_comments":[{"comment":"The tensor operator definition in Eq. (7) contains a typo: S^(jk)(r̂) = 3 r̂·σ^(j) r̂·σ^(k) − σ^(i)·σ^(j) should read − σ^(j)·σ^(k), as in Eq. (A4).","section":"Eq. (7) and Eq. (A4)"},{"comment":"The text says the L_CM term is \"actually non-negligible\" in relation to Ref. [48], but later states that the LS matrix element is negligible because of a cancellation between L and L_CM. Please clarify that the individual contributions are non-negligible while their sum is small.","section":"Sec. V.A"},{"comment":"The regulator is written as R_S = 0.8 fm^{-1}, but the Gaussian in Eq. (A13) requires R_S to have dimensions of length. This is presumably a typo for R_S = 0.8 fm; since the contact matrix elements M_CT depend on R_S, this should be corrected.","section":"Eq. (A13)"},{"comment":"The error labeled σ_{g_NN_V} in Eqs. (46)–(49) is the width of an assumed LEC range, not a determined uncertainty. The manuscript should explicitly state this is a prior and that the central values would shift if the LECs were found to lie outside the assumed ±1 range.","section":"Eqs. (33), (46)–(49)"},{"comment":"In the sentence following Eq. (44), \"the first error\" and \"the second error\" are clear, but the notation (56)gV(87)µ could be made more explicit by saying \"the first error is from g_V and the second from the scale µ.\"","section":"Sec. VI, Eq. (44)"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically sound in its many-body machinery and transparent about its main limitation. My recommendation of major_revision is driven by the gap between the abstract's \"good agreement\" claim and the LEC-free numerical comparison, plus the unquantified three-body O(α^2) operator. These are fixable with revised wording and an explicit uncertainty statement; the calculation itself does not need to be redone. The authors should also double-check the R_S units in Eq. (A13)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first QMC calculation of the nuclear-structure correction δ_NS in 10C, and it is a careful, transparent piece of work. The catch is that the abstract's 'good agreement' with previous evaluations is bought with an arbitrary band for two low-energy constants, not with the calculation itself.\n\nThe genuinely new content is real: VMC and GFMC matrix elements for δ_NS with four Hamiltonians, including the corrected spin-orbit operator (sign and center-of-mass term) and an explicit check that the Fermi matrix element comes out at 7.07. The operator definitions in Appendix A are complete enough to reproduce, and the paper honestly reports that the dominant uncertainty is the undetermined LECs g_V1^NN and g_V2^NN. The density plots and the VMC/GFMC comparison are useful diagnostics. This is a solid methodological step, not a precision extraction.\n\nThe soft spots are in how the result is framed. The LEC-free content says one thing; after adding the ±1 naturalness range it says another. The magnetic-plus-spin-orbit piece alone is about 3σ above the Hardy-Towner δ_NS,B, and the full LEC-free central value sits ~0.5–0.9×10^-3 (roughly 2σ) above the NCSM result. The ±(0.5–0.8)×10^-3 band from g̃=±1 converts this into 'agreement within error.' That is underdetermination, not prediction, and the abstract oversells it. There is also a tension inside the paper: the two-nucleon weak amplitudes from AV18 and NV2-Ia agree to better than 2%, which suggests small contact LECs, while the ±1 assignment permits contact terms ~10% of the magnetic ones. The paper never resolves which regime is physical.\n\nSo for whom? This is for people working on V_ud, on chiral EFT currents, and on ab initio nuclear structure. It is a useful benchmark and a start toward fitting the LECs from multiple isotopes or two-nucleon data. I would send it to a serious referee. The referee should ask for a less loaded abstract and a more direct treatment of the LEC range — for example, a figure showing δ_NS as a function of the LEC values, or constraints from two-nucleon amplitudes.","headline":"First QMC calculation of δ_NS in 10C, technically solid and transparent, but the 'good agreement' with prior results is carried by an arbitrary LEC band.","tokens_in":24653,"tokens_out":4274,"would_cite":true,"duration_ms":50242,"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":"This paper computes the nuclear-structure correction to 10C superallowed beta decay with quantum Monte Carlo wave functions and chiral EFT operators, and finds it agrees with both the traditional survey value and the recent no-core shell mo","keywords":["superallowed beta decay","nuclear-structure radiative corrections","V_ud extraction","chiral effective field theory","quantum Monte Carlo","Green's function Monte Carlo","carbon-10 decay","CKM unitarity"],"falsifier":"Measure the charge radius of 10B*(0+;1) and check the paper's claimed linear correlations with M_F^E and M_F^+, and/or compute the two contact couplings g_V1^NN and g_V2^NN on the lattice: if the true LECs land outside the +/- (1/m_N)(2F_pi)^-2 band, the central values and error bars of delta_bar_NS and V_ud shift by more than the quoted ranges.","tokens_in":23303,"feed_emoji":"⚛️","tokens_out":8750,"duration_ms":99622,"temperature":0.7,"pith_summary":"This paper tries to establish that the nuclear-structure-dependent radiative correction to the superallowed beta decay of carbon-10 can be computed from first principles with quantum Monte Carlo wave functions and chiral effective field theory operators, and that the result agrees with the two existing determinations. If correct, it provides an independent, ab initio cross-check of the correction that feeds the most precise determination of the CKM matrix element V_ud, whose first-row unitarity is currently in tension with the Standard Model. The calculation splits the correction into an energy-independent piece carried by magnetic, tensor, spin-orbit, contact, and O(alpha^2) two-body operators, plus a smaller energy-dependent piece. Across four nuclear Hamiltonians, the QMC results give a long-range contribution of -(4.06 to 4.43) x 10^-3 and, once the unknown contact couplings are assigned an arbitrary range, a full delta_bar_NS compatible with the survey value -4.0(5) x 10^-3 and the NCSM value -4.22(32) x 10^-3. The largest remaining uncertainty is not the many-body method but two undetermined low-energy constants.","feed_headline":"New calculation confirms the nuclear correction in 10C beta decay","feed_subtitle":"Quantum Monte Carlo result agrees with both existing extractions; two unknown constants dominate the error.","key_machinery":"The load-bearing object is the EFT decomposition of the radiative correction into an energy-independent piece delta_NS^(0) and an energy-dependent piece delta_NS^E, with the isospin structure separated into spectator-proton and spectator-neutron contributions. Each piece is a sum over nuclear matrix elements of two-body transition operators labelled Fermi, Gamow-Teller, tensor, and spin-orbit, defined through radial functions h(r) such as the ~1/r magnetic and spin-orbit pieces, the delta-function contact terms carrying the two unknown LECs g_V1^NN and g_V2^NN, and the logarithmic O(alpha^2) Fermi term. The many-body matrix elements and their radial densities C(r) are computed with variation","core_discovery":"The paper's central claim is that the nuclear-structure-dependent radiative correction delta_bar_NS for 10C -> 10B* can be evaluated using chiral EFT two-body current operators and quantum Monte Carlo (VMC and GFMC) wave functions, and that the resulting correction is consistent with the standard survey value and the more recent dispersive no-core shell model calculation. Concretely, the GFMC magnetic-plus-spin-orbit part of delta_NS^(0) sits in the range -[4.06,4.43] x 10^-3, roughly 40 percent larger in magnitude than the delta_NS,B = -3.06(35) x 10^-3 used in the traditional survey, with the difference attributed to many-body correlations in the GFMC wave functions. The energy-dependent p","pith_inferences":["If the two contact LECs were determined from lattice QCD or from modeled two-nucleon weak amplitudes, the dominant uncertainty in this approach would drop, likely making EFT+QMC competitive with or more precise than the dispersive NCSM evaluation.","The strong empirical correlation between radius and the spin-independent matrix elements suggests a cheap experimental route: measuring the charge radius of the short-lived 10B* 0+ state would pin down the model dependence of the largest energy-dependent term without waiting for a full QCD calculation of the LECs.","The sensitivity of the spin-dependent matrix elements to OPE correlations, and the nodal densities that drive GFMC changes, warn that calculations in heavier superallowed emitters need wave functions that reproduce spin-isospin correlations, not just energies and radii.","Clarifying how pion-range contributions depend on the definition of the isospin limit, as the paper notes, is a prerequisite for consistently separating delta_C from delta_bar_NS at higher chiral orders; until that is settled, the current delta_C input is taken from the survey analysis rather than computed in the same framework."],"forward_implications":["The long-range magnetic-plus-spin-orbit part of delta_NS^(0) comes out at -[4.06,4.43] x 10^-3 in GFMC, about 40 percent larger in magnitude than the survey's delta_NS,B, so the difference is attributed to many-body correlations rather than to physics missing from the EFT.","With the arbitrary contact-LEC range included, EFT+QMC, the traditional survey analysis, and the NCSM dispersion calculation all agree on delta_bar_NS within errors, removing a reason to suspect a large nuclear-structure error in V_ud extraction.","The 10C-only extraction of V_ud from EFT+QMC falls between 0.97336 and 0.97355, compatible with the survey-based (0.97318) and NCSM-based (0.97317) values within the experimental error; the theory error from the unknown LECs is larger than the spread from the four Hamiltonians.","Because the spin-independent matrix elements M_F^E and M_F^+ correlate strongly with the charge radius of 10B*, a measured radius of the 0+ daughter would directly reduce the model dependence of the energy-dependent correction.","The corrected spin-orbit operator and the non-negligible L^CM term change the operator set of Ref. [48]; these changes matter most in systems where the LS and L^CM contributions do not cancel."],"supporting_citations":[{"why":"Supplies the chiral EFT operator formalism and master formulas for delta_bar_NS that the paper applies to 10C, correcting its spin-orbit sign and adding the L^CM term.","marker":"[48]"},{"why":"The no-core shell model plus dispersion calculation of delta_NS in 10C that the paper benchmarks against.","marker":"[50]"},{"why":"The survey providing the 10C half-life, Q_EC, delta_C, and the standard value delta_NS = -4.0(5) x 10^-3.","marker":"[45]"},{"why":"Provides g_V at the pion scale and the nonperturbative W-gamma box uncertainty used in the V_ud master formula.","marker":"[9]"},{"why":"The standard review of VMC and GFMC methods, including the mixed-estimate and off-diagonal propagation formulas used here.","marker":"[53]"},{"why":"The phenomenological AV18 two-nucleon potential used as one of the four Hamiltonians.","marker":"[58]"},{"why":"The NV2-Ia chiral two-body interaction underlying the two chiral Hamiltonians.","marker":"[61]"},{"why":"The NV2+3 three-nucleon force fit used for the NV2+3-Ia Hamiltonian.","marker":"[62]"},{"why":"The alternative three-nucleon fit (Ia*) used to test sensitivity to the three-body force.","marker":"[63]"},{"why":"The dispersive framework that Ref. [50] combines with NCSM; supplies the comparison method.","marker":"[47]"}],"fun_headline_variants":["QMC-EFT calculation matches both extractions for 10C decay","10C beta decay: QMC-EFT agrees with both extractions","Two unknown constants dominate error in 10C correction","New 10C correction from QMC agrees with both V_ud routes"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"Two short-range coupling constants in the nuclear operator are not known from experiment or QCD; the paper arbitrarily sets them to plus or minus one in dimensionless form, and this choice drives the largest uncertainty in the final correction to V_ud.","fun_headline_variants_meta":{"raw":{"variants":["QMC-EFT calculation matches both extractions for 10C decay","10C beta decay: QMC-EFT agrees with both extractions","Two unknown constants dominate error in 10C correction","New 10C correction from QMC agrees with both V_ud routes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001371,"raw_usage":{"total_tokens":5375,"prompt_tokens":707,"completion_tokens":4668,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":451,"completion_tokens_details":{"reasoning_tokens":4592}},"tokens_in":451,"tokens_out":4668,"duration_ms":39932,"temperature":1.0,"reasoning_tokens":4592,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T22:26:16.398453+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the charge radius of 10B*(0+;1) and check the paper's claimed linear correlations with M_F^E and M_F^+, and/or compute the two contact couplings g_V1^NN and g_V2^NN on the lattice: if the true LECs land outside the +/- (1/m_N)(2F_pi)^-2 band, the central values and error bars of delta_bar_NS and V_ud shift by more than the quoted ranges.","supporting_citations":[{"cited_title":"A model with vectorlike fermions and $U(1)_X$ symmetry: CKM unitarity, $b \\rightarrow s$ transitions, and prospect at Belle II","cited_arxiv_id":"2303.14913","evidence_quote":"The survey providing the 10C half-life, Q_EC, delta_C, and the standard value delta_NS = -4.0(5) x 10^-3."},{"cited_title":"Dispersive formalism for the nuclear structure correction $\\delta_\\mathrm{NS}$ to the $\\beta$ decay rate","cited_arxiv_id":"2211.10214","evidence_quote":"The standard review of VMC and GFMC methods, including the mixed-estimate and off-diagonal propagation formulas used here."}],"review_version":1}