{"id":"c414758a-d90e-49bf-8436-96e90651a2fb","arxiv_id":"2412.14160","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Chiral EFT corrections to the 2νββ electron spectrum from weak magnetism and pion exchange appear at next-to-leading order and must be included in searches for new physics.","lead":"This paper computes quantum chromodynamics-derived chiral corrections to the electron energy spectrum of two-neutrino double beta decay. The corrections are large enough that experiments using 2νββ spectra to search for new physics should include them, and they connect 2νββ to the matrix elements of neutrinoless double beta decay.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The pionic spectral distortion is not a controlled prediction until the short-range LECs are included: natural O(1) values in Eq. (11) can dominate or cancel the computed pion-range effect, so the quantitative per-mille claim rests on an unverified assumption.","rationale":"The reader identified the same weakest point: the unknown short-range LECs gNN_2ν,F and gNN_2ν,GT. My stress test confirms this is the single most load-bearing assumption. It is not a manufactured concern: using the paper's own NMEs, natural O(1) values of these LECs produce contact-term contributions to ϵF and ϵGT that are comparable to or larger than the computed pion-range terms in Eqs. (14)-(15), so the pionic spectral distortions in Section 4 are currently an estimate with an unquantified error, not a controlled prediction. The weak-magnetism piece is robust because it depends only on known couplings (gM, gA, mN). The 0νββ-NME extraction is explicitly conditional on the same LECs. The paper is transparent about this limitation, which supports a CONDITIONAL rather than REJECT verdict. The concrete scan over gNN = ±1 would quantify whether natural-size short-range effects alter the plotted pionic shape; if they do, the abstract's per-mille pionic phrasing should be softened. No change to the reader's verdict is required.","tokens_in":21880,"tokens_out":16433,"duration_ms":152786,"concrete_test":"Using Eqs. (14), (15), and Table 1, recompute δSπ(ϵ) of Figs. 4/7 with the contact terms included for the four combinations gNN_2ν,F, gNN_2ν,GT = ±1, keeping all other inputs fixed. If the resulting normalized pionic distortion changes by more than the plotted pion-only band or changes sign, then the Section 4 pionic estimates are not robust under the stated naturalness assumption, and the paper's quantitative pionic claims should be reported as LEC-dependent rather than as per-mille predictions. A secondary lattice check would be to compute the two LECs along the lines of Refs. [54,55] to fix their actual values.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4's numerical pionic corrections and the derived 0νββ-NME connection depend on the assumption in Eqs. (10)-(11) that the short-range LECs gNN_2ν,F and gNN_2ν,GT are naturally O(1) and subdominant. The paper explicitly omits them (\"we do not control... only provide an incomplete estimate\"), but the omission is not a small effect. In Eq. (14), the contact term enters as [M_F(mπ) + (mπ^2/Fπ^2) gNN_2ν,F M_F,sd]. Using Table 1 for 76Ge QRPA, this is approximately [-1.16 + 2.26*gNN_2ν,F*(-3.46)], so even gNN_2ν,F=+1 changes ϵF by roughly a factor of seven, while gNN_2ν,F=-1 changes its sign. For ϵGT in Eq. (15), the long-range combination M_GT^AA - 2M_GT^AP + 4M_GT^PP ≈ 8.5, while the contact term is -3(mπ^2/(gA^2Fπ^2))gNN_2ν,GT M_F,sd ≈ +14.6*gNN_2ν,GT, again of order or larger than the long-range piece for |gNN|~1. Consequently the pionic spectral distortions in Figs. 4-7 and the lifetime shifts in Table 2 carry an unquantified, potentially dominant error. The weak-magnetism contribution is controlled; the pionic part is not yet a prediction. The 0νββ-NME extraction is explicitly gated on the same LECs, so that component of the central claim is prospective.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a chiral effective field theory description of two-neutrino double beta decay (2νββ), focusing on the electron energy spectrum. It derives next-to-leading-order corrections from weak magnetism and from double-weak pion-exchange diagrams, together with the associated short-range counterterms. Using the lepton-energy expansion of Šimkovic et al. and nuclear matrix elements from QRPA and shell-model calculations for 76Ge and 136Xe, it computes the impact of these corrections on the normalized spectrum and on total rates. It compares the chiral distortions with those from sterile neutrinos and non-standard charged-current interactions, and argues that the pion-exchange contributions involve matrix elements related to 0νββ, potentially allowing 0νββ NMEs to be constrained by precision 2νββ spectral measurements.","tokens_in":22194,"tokens_out":8833,"duration_ms":77084,"significance":"If the central claims hold, precision 2νββ spectral analyses will need to incorporate the weak-magnetism correction, which is a genuine parameter-free prediction at the per-mille level and is robust within the paper's framework. The explicit propagation of nuclear-structure uncertainties in the lepton-energy expansion and the validation of that expansion for the shell model are strengths. The pion-exchange part and the proposed link to 0νββ NMEs are interesting and potentially important, but they are currently conditional on unknown short-range low-energy constants and on the poorly known NME ratios ξ31 and ξ51; as presented, they do not yet constitute a controlled numerical prediction of the spectral distortion.","major_comments":[{"comment":"The numerical pion-exchange results are not controlled by the naturalness assumption in Eq. (11). For 76Ge QRPA, the long-range Fermi matrix element in Eq. (14) is M_F(mπ) ≈ -1.16, while the contact term is (mπ^2/Fπ^2) gNN_2ν,F M_F,sd ≈ 2.3 × (-3.46) gNN_2ν,F ≈ -7.9 gNN_2ν,F; hence gNN_2ν,F = +1 changes ϵF by roughly a factor 8 and gNN_2ν,F = -1 reverses its sign. In Eq. (15), the long-range GT combination M_GT^AA - 2 M_GT^AP + 4 M_GT^PP ≈ 8.5 is confronted by a contact contribution ≈ -14.6 gNN_2ν,GT. Consequently the pionic spectral distortions in Figs. 4-7 and the lifetime shifts in Table 2 carry an unquantified, potentially dominant error. The manuscript acknowledges this in Section 2 (after Table 1) and in the Fig. 4 discussion, but the abstract and Section 6 present the pion-exchange corrections as a quantitative result. The revision should either provide estimates of the short-range LECs or explicitly restrict the quantitative 'should be included' claim to the weak-magnetism contribution.","section":"Section 4, Eqs. (14)-(15), Table 1"},{"comment":"The claim that 0νββ NMEs can be isolated from 2νββ spectra is conditional on the same unknown short-range LECs and on the nuclear-structure ratios ξ31 and ξ51. The paper states that with current LECs this program cannot be carried out, yet the concluding bullet '0νββ NMEs can, in principle, be isolated in 2νββ measurements' goes beyond what is demonstrated. For 136Xe, the QRPA and NSM values of ξ31 differ by roughly a factor of 2 (Table 1 and Eqs. (38)-(39)), and Figs. 6-7 show the pionic component is washed out by this uncertainty. Please soften this conclusion to reflect the gating, or specify the LEC and ξ precision required to make the extraction viable.","section":"Section 6 and abstract"}],"minor_comments":[{"comment":"The known limitations of the Fermi-function treatment (finite nuclear size, electron screening, radiative corrections) are acknowledged, but their ~0.1% effect on the normalized spectrum is not included in the uncertainty bands of Figs. 8-10; a brief statement that these corrections constitute a common systematic for all curves would make the comparison fairer.","section":"Section 2, Eq. (24)"},{"comment":"The NSM uncertainty estimate is based on multiplying individual matrix elements by random factors and imposing the half-life range of Eq. (36); please clarify whether this procedure is meant to cover systematic method spread, since for 136Xe it gives ξ31 ≈ 0.16, a factor of 2 below the QRPA value ξ31 ≈ 0.32.","section":"Section 3, Eq. (39)"},{"comment":"The lifetime shifts labeled T_π^(2) and T_χ^(2) include only the long-range pion piece; the table caption should state explicitly that the short-range LEC contribution is omitted, so these entries are not complete chiral predictions.","section":"Table 2 and Section 4"},{"comment":"There are several minor typographical issues, including 'T able 1' in the Table 1 caption and 'incomplete estimate the size' in Section 2 after Eq. (20); these should be corrected.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"This is a solid theory paper from an experienced group, and the weak-magnetism result should be credited. The central issue is the overstatement of the pion-exchange corrections, which the authors themselves partly acknowledge. A revision that either estimates the short-range LECs or clearly limits the quantitative claims to the weak-magnetism contribution would make the paper publishable. The QRPA/NSM discrepancy for 136Xe ξ31 is worth highlighting to the authors as a point that needs more discussion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the weak-magnetism part of this paper is solid and should be taken seriously by anyone doing precision 2νββ spectral work. The pion-exchange part is genuinely new formalism, but its numerical spectral distortions are not yet a controlled prediction. The authors say this themselves, which is to their credit.\n\nWhat's new: the double-weak pion-exchange two-neutrino potential, and the observation that its matrix elements are related to 0νββ NMEs at the pion scale. That link is concrete and I don't think it appears in the earlier literature. The paper also validates the lepton-energy expansion of Šimkovic et al. for the shell model, and shows that weak-magnetism effects survive nuclear uncertainties near the nodes. The comparison with sterile-neutrino and tensor BSM distortions is useful and well presented.\n\nSoft spots: the short-range LECs gNN_2ν,F and gNN_2ν,GT enter at the same chiral order as the pion-range pieces, and the stress-test arithmetic checks out. For 76Ge QRPA, Eq. (14) gives a long-range MF of about -1.16 while the contact contribution is roughly -7.8 × gNN_2ν,F, so natural O(1) values change epsilon_F by a factor of several or flip its sign. The GT channel behaves similarly. That means the pionic spectral distortions in Figs. 4-7 and the lifetime shifts in Table 2 carry an unquantified, potentially dominant uncertainty. The paper explicitly omits these LECs and calls the estimate incomplete, so the flaw is not hidden, but the abstract's claim that the corrections \"should be included in analyses\" overstates what is controlled. The weak-magnetism correction is free of this problem and is the paper's real quantitative result.\n\nAlso missing, and acknowledged, are finite-size, screening, and radiative electromagnetic corrections at the same per-mille level as the chiral effects, so the error budget is not closed. No code or data is shipped. The nuclear ratios ξ31 and ξ51 agree between QRPA and shell model for Ge but differ by about a factor of two for Xe, which weakens the pionic extraction further. Still, the central EFT derivation and multipole organization are sound.\n\nThis paper is for people analyzing 2νββ spectra for BSM physics and for those working on 0νββ NME correlations. It deserves a serious referee. I would ask the authors to either estimate the LECs from lattice QCD or naturalness, or clearly label the pionic numbers as a range, and to add a systematic estimate of missing electromagnetic corrections. With those caveats stated up front, the weak-magnetism result is publishable and the pionic part is worth publishing as a first, incomplete step.","headline":"A technically solid chiral-EFT derivation whose weak-magnetism prediction is robust, while the pion-exchange numbers are an honest but uncontrolled estimate until short-range LECs are determined.","tokens_in":22861,"tokens_out":2501,"would_cite":true,"duration_ms":25776,"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":"The electron spectrum of two-neutrino double beta decay carries per-mille-level chiral corrections that precision searches for new physics must include.","keywords":["two-neutrino double beta decay","chiral effective field theory","weak magnetism","pion-exchange currents","two-neutrino potential","nuclear matrix elements","neutrinoless double beta decay","electron spectrum"],"falsifier":"Compute the two short-range couplings on the lattice (the Gamow-Teller one is related to an isotensor axial polarizability already studied in lattice QCD): if either is far from order one, the predicted pionic distortion changes size or sign. Separately, a normalized 136Xe spectrum measured to about 0.1% precision should show the predicted pattern—mid-spectrum enhancement, attenuations near two-fifths and four-fifths of the energy range—after weak magnetism and the lepton-energy expansion are subtracted; its absence would rule out the claimed pionic correction.","tokens_in":21627,"feed_emoji":"⚛️","tokens_out":8039,"duration_ms":71798,"temperature":0.7,"pith_summary":"This paper argues that the electron energy spectrum of two-neutrino double beta decay, the slowest nuclear process ever measured, carries small but systematic corrections that any precision analysis must include. Working in chiral effective field theory, the authors identify the next-to-leading-order contributions: weak magnetism from the nucleon current, and double-weak pion exchange between nucleons together with the short-range contact terms that renormalize it. These corrections distort the normalized spectrum at the per-mille level, comparable to the distortions expected from sterile neutrinos, right-handed currents, or tensor interactions, so omitting them can bias or fake beyond-Standard-Model signals. The pion-exchange part is governed by nuclear matrix elements of the same two-body operators that appear in neutrinoless double beta decay, opening a possible path to constrain those matrix elements from ordinary 2νββ spectra once two unknown short-range couplings are determined.","feed_headline":"Chiral corrections reshape double-beta spectra at the per-mille level","feed_subtitle":"Searches for new physics in two-neutrino double beta decay need these next-to-leading-order corrections.","key_machinery":"The load-bearing object is the two-neutrino potential $\\mathcal{V}_{2\\nu}(q)$, built from double-weak pion-exchange diagrams plus a short-range contact potential regulated by two unknown low-energy constants $g^{NN}_{2\\nu,F}$ and $g^{NN}_{2\\nu,GT}$. This potential generates the pionic spectral corrections through Fermi and Gamow-Teller nuclear matrix elements at pion-mass momentum transfers, the same operator structure that appears in 0νββ. The argument is carried also by the lepton-energy expansion, which factorizes lepton phase space from nuclear ratios $\\xi_{31}$ and $\\xi_{51}$; the weak-magnetism correction enters through an energy-dependent amplitude that requires no new nuclear matrix elements. The paper's master result is the decay-rate kernel $C_{2\\nu}$ of Eq. (23), which combines the leading double Gamow-Teller transition with weak magnetism and the pion-exchange/contact corrections.","core_discovery":"Using chiral effective field theory, the paper derives the 2νββ decay amplitude and differential rate through next-to-leading order and isolates three new ingredients: subleading one-nucleon currents (weak magnetism), a two-nucleon two-neutrino potential from pion exchange, and short-range double-weak contact operators. It shows that weak magnetism shifts the spectral peak to higher energies by a few per-mille and is theoretically well controlled, while the pion-exchange terms produce a spectral modulation that is small for 76Ge and at the few per-mille level for 136Xe. The pion-exchange and contact contributions enter through nuclear matrix elements evaluated at momentum transfers of order the pion mass, which are related by stable ratios to the matrix elements of 0νββ mediated by light Majorana neutrino exchange; the paper therefore asks whether detailed 2νββ spectra can constrain 0νββ matrix elements. It concludes that the short-range low-energy constants are needed first, and that with present nuclear uncertainties the pionic terms are difficult to isolate, whereas weak magnetism is distinguishable, especially at spectral nodes where lepton-energy-expansion uncertainties vanish.","pith_inferences":["Beyond the paper's own analysis, the node structure of Eq. (40) suggests an experimental strategy: measuring the spectrum near those energies, where nuclear-structure uncertainties vanish, gives the cleanest window for isolating weak magnetism and any beyond-Standard-Model distortion.","If lattice QCD provides the unknown short-range couplings, the operator-level connection found here implies that 2νββ spectral fits could yield a data-driven, model-independent extraction of 0νββ Gamow-Teller matrix elements; the paper stops short of asserting this is achievable, but its operator-level mapping points in that direction.","The near-degeneracy between the tensor-coupling and weak-magnetism spectral shapes implies that a combined analysis across several isotopes with different Q-values could break the degeneracy, since the relative weight of the two contributions varies with isotope; the paper does not perform this multi-isotope fit."],"forward_implications":["A 0.1%-precision measurement of the 2νββ spectrum can constrain leptonic right-handed charged currents more tightly than current global fits, but only if the chiral corrections are included in the fit.","Weak magnetism mimics a tensor interaction with $|\\epsilon_T| \\sim 0.0014$; a fit that omits weak magnetism will misinterpret such a signal or cancel it.","The pion-exchange part of the spectrum shares its nuclear matrix elements with 0νββ, so 2νββ spectra could become an indirect handle on 0νββ matrix elements once the unknown short-range couplings are known.","Sterile-neutrino distortions and chiral distortions pull the spectral peak in opposite directions, so searches for sterile neutrinos in 2νββ data must subtract the chiral corrections to set reliable limits.","Total decay rates receive corrections of about 3% for 76Ge and 10% for 136Xe, though these are hidden by nuclear uncertainties in the absolute rate; the normalized shape is where the corrections become observable."],"supporting_citations":[{"why":"Supplies the lepton-energy expansion in terms of the nuclear ratios ξ31 and ξ51 and the QRPA matrix elements on which the spectral-shape analysis is built.","marker":"[11]"},{"why":"Provides the earlier 2νββ amplitude including weak magnetism that the paper extends within chiral power counting.","marker":"[29]"},{"why":"Provides the phase-space factors and Fermi-function prescription the paper reproduces and uses for decay-rate integrals.","marker":"[30]"},{"why":"Establishes the chiral-EFT treatment of 0νββ and the short-range operator needed to renormalize the pion-range neutrino potential, which the paper adapts to 2νββ.","marker":"[33]"},{"why":"Lattice calculation of the isotensor axial polarizability that is related to the short-range Gamow-Teller low-energy constant.","marker":"[54]"},{"why":"QRPA 0νββ nuclear matrix elements used to evaluate the pion-exchange corrections.","marker":"[58]"},{"why":"Shell-model 0νββ nuclear matrix elements used for the same evaluation in the nuclear shell model.","marker":"[59]"},{"why":"Shell-model 2νββ matrix elements used to validate the lepton-energy expansion and to compute the ξ ratios.","marker":"[64]"}],"fun_headline_variants":["Double-beta spectra need chiral corrections","Chiral EFT adds weak-magnetism and pion terms to 2νββ","Per-mille chiral shifts in two-neutrino double-beta decay","Chiral corrections key for double-beta new physics searches"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative size of the pion-exchange spectral distortion assumes that two unknown constants, which describe how two nucleons emit electrons and neutrinos at short distances, are naturally of order one; the paper leaves them out of its numerical estimates and states they could change the overall size of the effect.","fun_headline_variants_meta":{"raw":{"variants":["Double-beta spectra need chiral corrections","Chiral EFT adds weak-magnetism and pion terms to 2νββ","Per-mille chiral shifts in two-neutrino double-beta decay","Chiral corrections key for double-beta new physics searches"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00047,"raw_usage":{"total_tokens":2329,"prompt_tokens":922,"completion_tokens":1407,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":1347}},"tokens_in":538,"tokens_out":1407,"duration_ms":12340,"temperature":1.0,"reasoning_tokens":1347,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:25:08.748710+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the two short-range couplings on the lattice (the Gamow-Teller one is related to an isotensor axial polarizability already studied in lattice QCD): if either is far from order one, the predicted pionic distortion changes size or sign. Separately, a normalized 136Xe spectrum measured to about 0.1% precision should show the predicted pattern—mid-spectrum enhancement, attenuations near two-fifths and four-fifths of the energy range—after weak magnetism and the lepton-energy expansion are subtracted; its absence would rule out the claimed pionic correction.","supporting_citations":[{"cited_title":"The isotensor axial polarisability and lattice QCD input for nuclear double-$\\beta$ decay phenomenology","cited_arxiv_id":"1701.03456","evidence_quote":"Lattice calculation of the isotensor axial polarizability that is related to the short-range Gamow-Teller low-energy constant."}],"review_version":1}