{"id":"b3eac66c-91cf-4875-9a27-440215cec771","arxiv_id":"2607.24684","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"pnQRPA, shell model, IBM-2, and heavy-nucleus EFT give mutually consistent 2ν half-life ranges for 126Xe ECEC and 134Xe ββ, with 134Xe lower bounds below ~2×10^24 y.","lead":"Four nuclear many-body methods predict half-lives for two rare xenon double-weak decays that have not yet been seen. The lower end of the 134Xe range sits near next-generation experimental reach, giving a concrete target for ongoing xenon detectors.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The headline \"within reach of next-generation experiments\" claim rests on the unquenched/minimal-quenching endpoints of each band — values the paper's own quenching systematics disfavor — and for the NSM the lower edge (2.04×10^24 y) actually exceeds the quoted 2×10^24 y threshold.","rationale":"The reader identified phenomenological quenching/LEC normalization as the weakest assumption and set CONDITIONAL with medium correctness risk. I agree with that diagnosis but locate the load-bearing point one step further downstream: not merely that quenching dominates the uncertainty bands, but that the paper's most quotable claim — next-generation experimental reachability of 134Xe — is carried specifically by the unquenched endpoints (pnQRPA g_A^eff=1.27, IBM-2 q=1/0.788) and, for the NSM, is numerically not quite true as stated (2.04 vs 2×10^24 y). This is a concern about the framing of the strongest claim rather than about the correctness of the NME calculations themselves: the tables and figures are internally consistent, the uncertainty treatment is unusually careful for this literature (explicit EFT truncation errors, truncation-convergence tables for the NSM, separate SSD/HSD handling for IBM-2), and the core result — overlapping multi-method bands around 10^24–10^25 y for 134Xe and an order of magnitude longer for 126Xe — is directly supported by Table 1 and survives any reasonable quenching choice. The overlap claim is somewhat weakened by shared calibration data across methods, but the paper is transparent about every normalization. Because the science is accept-shaped and the fix is a wording/qualification change plus the proposed B(GT) measurements (which the authors themselves request), I recommend the verdict remain CONDITIONAL, unchanged from the reader.","tokens_in":23011,"tokens_out":2654,"duration_ms":112174,"concrete_test":"Recompute the Table 1 / Fig. 1 lower edges for 134Xe using only each method's preferred quenched coupling: pnQRPA with g_A^eff = 0.8, IBM-2 with q = A^-0.18 and q = 0.788, NSM with the central quenching of each interaction (GCN5082 q≈0.50, QX q≈0.75), and EFT at central log ft values. Count how many methods retain T_1/2^lower below the projected LZ sensitivity of 1.7×10^24 y. If two or fewer survive, the abstract's \"for all calculations ... within the reach of next-generation experiments\" should be qualified to name the specific quenching choices that produce reachable half-lives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's punchline is that \"for all calculations the lower range of the predicted 134Xe half-life is shorter than T≈2×10^24 y.\" Inspection of Tables 1, B.3, and D.7 shows where those lower edges come from. pnQRPA: the 1.22×10^24 y edge corresponds to g_A^eff = 1.27, i.e., q = 1, no quenching at all; with the quenched value g_A^eff = 0.8 the half-life is 3.8×10^24 y (Table B.3). IBM-2 HSD: the short edges (0.081–0.21×10^24 y) come from q = 1 (bare g_A) and the mild q = 0.788; the maximally quenched case q = A^-0.18 gives 2.74×10^24 y, above threshold. NSM: the lower edge is 2.04×10^24 y — marginally *longer* than the 2×10^24 y the abstract claims all methods fall below. EFT: the 0.424×10^24 y edge is the extreme of a band built from hypothetical log ft values and an empirical β−/EC ratio r, which the authors themselves identify as the dominant and least-controlled uncertainty (§3.4, Fig. 2). So the \"experimental reachability\" statement is systematically driven by the least-quenched, least-favored corner of each method's parameter space; the authors argue elsewhere (§3.1–3.3) that quenching is required to reproduce neighboring measured decays, meaning the physically preferred points mostly sit above the projected 1.7×10^24 y LZ sensitivity. This does not make the claim false — the bands do include short half-lives, and the cross-method overlap around 2×10^24 y is real — but the abstract overstates how generic the reachable region is. A secondary, related point: the four methods are not independent confirmations of the low edges, since pnQRPA g_pp, NSM q, IBM-2 q, and EFT LECs are all normalized to the same small pool of neighboring data (124Xe, 136Xe, 130Te, GT systematics), so their agreement at the low end partly reflects shared calibration rather than independent structure physics.","agreement_with_reader":"partial"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The paper predicts nuclear matrix elements and half-lives for the two-neutrino double-electron capture of 126Xe and the two-neutrino double-beta decay of 134Xe using four many-body approaches: pnQRPA (with isospin-restoration g_pp fitting to measured 124Xe/136Xe decays), the nuclear shell model (GCN5082 and QX interactions, with a jump=3 truncation for 126Xe), the microscopic IBM-2 (under both SSD and HSD closure assumptions), and an EFT for heavy nuclei (LECs fitted to hypothetical log ft values with an EFT truncation uncertainty). Each method carries an estimated theoretical uncertainty, dominated by quenching (phenomenological methods) or LEC/ratio uncertainties (EFT). The central findings are: (i) all predictions are mutually consistent within uncertainties except the IBM-2 SSD case; (ii) 126Xe half-lives are typically an order of magnitude longer than 134Xe ones; (iii) the lower edges of the 134Xe bands approach or fall below ~2×10^24 y, near the projected 1.7×10^24 y LZ sensitivity.","tokens_in":23543,"tokens_out":3501,"duration_ms":121864,"significance":"If the bands hold, this is a useful, timely benchmark: 134Xe is directly accessible to running and planned xenon experiments (PandaX, LZ, XLZD), and a measurement would discriminate among the four methods in the same way the 124Xe ECEC measurement did. The paper's main strengths are methodological transparency and uncertainty honesty: each method's fitting procedure is documented in dedicated appendices (pnQRPA g_pp fits in App. B/Fig. B.6, NSM truncation systematics in App. C/Table C.4, IBM-2 Hamiltonian parameters in App. D, and the EFT hypothetical-log ft construction in App. E including a LogFT-vs-BetaShape comparison and an independent cross-check against the 136Xe(3He,t) B(GT) datum). The running-sum comparison of pnQRPA and NSM (Figs. 3–4) is a genuinely informative diagnostic of how two methods arrive at similar NMEs through very different strength distributions. The work is explicitly not first-principles — g_pp, quenching factors, and LECs are all anchored to neighboring measured decays or constructed log ft systematics — but the authors state this clearly, and the resulting predictions are falsifiable on a realistic experimental timescale, which is the standard by which ","major_comments":[{"comment":"The abstract states that 'for all calculations the lower range of the predicted 134Xe half-life is shorter than T≈2×10^24 y.' This is contradicted by the paper's own Table 1 and Table C.5: the NSM lower edge is 2.04×10^24 y (GCN5082 with maximal quenching), i.e., marginally *longer* than 2×10^24 y. The body text (§4) handles this correctly, saying the LZ projected sensitivity of 1.7×10^24 y is 'within the range of the pnQRPA, IBM-2 (HSD) and EFT predictions, and very close to the NSM one' — but the abstract does not. Since this is the paper's headline phenomenological claim, the abstract should be brought into agreement with Table 1 (e.g., 'shorter than or comparable to ~2×10^24 y').","section":"Abstract vs. Table 1 / Table C.5"},{"comment":"Relatedly, the reachability statement in the abstract should be qualified by where the short half-lives come from. Inspecting Tables B.3 and D.7: the pnQRPA lower edge (1.22×10^24 y) corresponds to g_A^eff = 1.27, i.e., zero quenching, while with the quenched g_A^eff = 0.8 the half-life is 3.8×10^24 y; the short IBM-2 HSD edges come from q = 1 and q = 0.788, while the maximally quenched q = A^{-0.18} gives 2.74×10^24 y; and the EFT low edge (0.424×10^24 y) is the extreme of the LEC band, which the authors themselves identify (§3.4, Fig. 2) as dominated by the least-controlled ingredient, the empirical β−/EC log ft ratio r. Since the paper argues elsewhere that quenching is *required* to reproduce neighboring measured decays, the physically preferred points within each band sit systematically above the quoted lower edges. The claim is not false — the bands do extend below threshold and th","section":"Abstract / §5, Tables B.3, D.7"},{"comment":"For 126Xe the NSM uses a jump=3 truncation and states explicitly that no truncation uncertainty is included. Table C.4 shows the NME is still dropping steeply with truncation: for GCN5082, M_eff goes 0.061–0.113 (jump=0) → 0.049–0.091 (jump=2) → 0.024–0.044 (jump=3), roughly a factor-of-two reduction in the last step, with no demonstration that jump=4 is converged or computationally inaccessible. Fig. C.7 shows the change is concentrated in the lowest 1+ state, so this is not a diffuse many-state effect that can be argued away. Given that the cross-method consistency claim for 126Xe (Table 1, Fig. 1 upper panel) relies on the NSM band sitting at 11.2–38.1×10^24 y, an unquantified factor-~2 systematic at the last truncation step is load-bearing. At minimum the authors should estimate the residual truncation error (e.g., from the jump=2→3 trend or the 124Xe experience in Ref. [93]) and eit","section":"Appendix C / Table C.4"}],"minor_comments":[{"comment":"Table B.3 lists single values (not ranges) of M^2ν for 134Xe (0.053 and 0.037) while 126Xe carries ranges from the unknown 124I 1+_1 energy. A brief note on why the 136Xe-anchored g_pp fit for 134Xe yields no analogous range (the 134Cs 1+_1 energy is known, Table A.2) would help the reader.","section":"Table B.3"},{"comment":"pnQRPA and NSM compute non-closure NMEs with explicit energy denominators but are paired with the HSD (average-energy) PSF; the SSD/HSD PSFs differ by only ~2% (Table A.2), so this is numerically harmless, but one sentence explaining the choice would preempt confusion about double-counting the closure energy.","section":"Table 1 / Appendix A"},{"comment":"Fig. 1 would be more informative if the edges of each band were annotated with the quenching factor or g_A^eff value that produces them, since the paper's own discussion shows the band positions are driven primarily by q. This would also make the point of major comment 2 visible at a glance.","section":"Fig. 1"},{"comment":"In §3.3, the HSD closure energy systematics '1.12 A^{1/2} MeV' appears without a reference; please add the source (or state it is from prior IBM-2 work). Also, Eq. (5) has a stray comma after '⟨E_k⟩, is simply replaced'.","section":"§3.3, Eq. (5)"},{"comment":"The semi-empirical-formula comparison (§4) is a useful sanity check; it would strengthen the paper to note explicitly that the SEF values (4.61×10^25 y and 4.03×10^24 y) lie inside all non-SSD bands, reinforcing the consistency claim.","section":"§4"},{"comment":"Typographical: 'we use four different widely-used many-body methods' (§1, 'different'/'used' repetition); 'half-live values' in §5 should be 'half-life values'; affiliation f lists 'Department of Physics and Astronomy at UNC' which reads oddly ('University of North Carolina at Chapel Hill'?).","section":"§1, §5, affiliations"}],"recommendation":"minor_revision","confidential_remarks":"The stress-test concern about the reachability claim resting on unquenched band edges does land, but only at the level of abstract framing: the body text (§4) already hedges correctly for the NSM, and the tables disclose the q-dependence of every edge. I therefore regard it as a wording correction rather than a substantive flaw. The genuinely new technical gap I found is the unquantified NSM jump=3 truncation for 126Xe, which the authors acknowledge but do not propagate. None of these require new calculations beyond a truncation-error estimate, hence minor revision. The author team overlaps heavily with the groups whose methods are being compared, so the 'independent methods' framing is somewhat softer than it appears, though the methods themselves are genuinely distinct."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful part of this paper is straightforward: first NSM numbers for 126Xe ECEC and 134Xe ββ, plus updated pnQRPA (larger bases, partial isospin restoration), IBM-2 SSD/HSD splits, and an EFT band with explicit LEC+truncation errors. Table 1 and Fig. 1 give overlapping half-life ranges once each method’s (non-statistical) uncertainty is folded in. That cross-check is what the xenon and 2ν communities will actually use.\n\nThey do the bookkeeping carefully. Running sums (Figs. 3–4) show why pnQRPA and NSM end up similar despite very different intermediate-state structure. Appendices spell out the g_pp fits, jump=3 truncation, IBM-2 parameters (including new 126Te), and the hypothetical log-ft construction for the EFT. Citations to prior pnQRPA/IBM-2/PHFB work are fair; the new pieces are clearly marked.\n\nSoft spots are real but proportionate. Quenching, g_pp, and EFT LECs are all normalized to the same small pool of neighboring 2ν/GT data, so the four methods are not fully independent confirmations. The abstract’s claim that every method’s lower edge sits below ~2×10^24 y is technically true for the bands as drawn, but those short edges come from the unquenched or lightly quenched corners (pnQRPA g_A^eff=1.27, IBM-2 q=1 or 0.788, EFT extreme LEC). The NSM lower edge is already 2.04×10^24 y, and the authors themselves argue elsewhere that quenching is required to match measured neighbors. Preferred central values mostly sit above the projected 1.7×10^24 y LZ sensitivity. IBM-2 SSD is an outlier the data already disfavor. None of this breaks the paper; it just means the “may be within reach” sentence should be read as “the short end of the uncertainty band,” not “the expected half-life.”\n\nNo code or input decks are released, which is normal for this literature but limits reproducibility.\n\nThis is for people who calculate or measure Xe double-weak rates. It deserves a serious referee. I would cite the Table 1 ranges and would bring it to reading group if we are discussing 2ν systematics or next-gen Xe sensitivity.","headline":"Solid multi-method half-life bands for two unobserved Xe 2ν modes; the “within reach” abstract line is real but driven by the least-quenched corners of each band.","tokens_in":24366,"tokens_out":591,"would_cite":true,"duration_ms":10492,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["23.40.-s","23.40.Hc","21.60.-n","27.60.+j"],"model":"grok-4.5","headline":"Four nuclear methods agree that 134Xe two-neutrino double-beta decay could be shorter than 2×10^24 years and within next-generation reach.","keywords":["two-neutrino double-beta decay","double-electron capture","nuclear matrix elements","xenon isotopes","pnQRPA","nuclear shell model","IBM-2","effective field theory"],"falsifier":"A measured half-life (or a firm lower limit tighter than ~10^24 y) for the two-neutrino double-beta decay of 134Xe that falls outside the overlapping theoretical window reported by the four methods.","tokens_in":23840,"feed_emoji":"⚛️","tokens_out":897,"duration_ms":19161,"temperature":0.7,"pith_summary":"This paper predicts the nuclear matrix elements and half-lives for two unobserved two-neutrino double-weak decays in xenon: double-electron capture of 126Xe and double-beta decay of 134Xe. It applies four independent many-body approaches—pnQRPA, the nuclear shell model, IBM-2, and an effective field theory for heavy nuclei—and reports theoretical uncertainty bands for each. The half-life ranges overlap once those uncertainties are included. For 134Xe the lower edge of every band lies below about 2×10^24 years, a sensitivity that next-generation xenon experiments project they can reach. 126Xe is typically predicted an order of magnitude slower. A positive detection, or a tighter experimental limit, would therefore benchmark the same nuclear methods used for neutrinoless double-beta decay searches.","feed_headline":"134Xe double-beta half-life may sit under 2×10^24 years","feed_subtitle":"Four nuclear methods agree the decay could be in reach of next-generation xenon experiments","key_machinery":"The two-neutrino nuclear matrix element M^{2ν} (and its quenched effective form M_eff^{2ν}), evaluated with four distinct many-body frameworks whose uncertainty bands are generated from quenching factors, pairing parameters, closure assumptions, or EFT truncation and low-energy constants.","core_discovery":"When theoretical uncertainties are assigned consistently inside each method, the four many-body calculations of the two-neutrino half-lives of 126Xe and 134Xe are mutually consistent; every calculation places the lower edge of the 134Xe half-life below ≈2×10^24 y, a window that projected next-generation xenon experiments can probe.","pith_inferences":["Because 2ν and 0ν matrix elements are known to correlate inside these frameworks, a measured 134Xe 2ν rate would tighten the nuclear uncertainty that currently limits 136Xe neutrinoless searches.","Charge-exchange measurements of the Gamow-Teller strength connecting the intermediate 1+ states to the initial and final nuclei would collapse the dominant EFT low-energy-constant uncertainty and shrink that band dramatically.","If future data favor the short half-life edge, the single-state-dominance closure used by IBM-2 would be preferred over higher-state dominance for this mass region."],"forward_implications":["Next-generation xenon experiments that reach ~1.7×10^24 y sensitivity can test part of every theory band for 134Xe.","A positive 134Xe detection would supply a new calibration point for the same nuclear methods used to predict neutrinoless double-beta matrix elements.","126Xe is predicted roughly ten times slower, so experimental priority naturally falls on 134Xe.","Consistency across pnQRPA, shell model, IBM-2 and EFT strengthens confidence that the shared nuclear structure input is not method-specific."],"fun_headline_variants":["Four methods place 134Xe 2ν half-life under 2×10^24 y","134Xe double-beta lower edge falls below 2×10^24 years","All four calculations put 134Xe within next-gen reach","Consistent 134Xe 2νββ half-lives sit under 2×10^24 y","134Xe two-neutrino decay may be next-gen accessible"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The half-life bands rest on phenomenological quenching factors and pairing or low-energy constants that are fitted to neighboring measured decays or to hypothetical log ft values; if those effective couplings misrepresent the target nuclei the predicted ranges shift by more than an order of magnitude.","fun_headline_variants_meta":{"raw":{"variants":["Four methods place 134Xe 2ν half-life under 2×10^24 y","134Xe double-beta lower edge falls below 2×10^24 years","All four calculations put 134Xe within next-gen reach","Consistent 134Xe 2νββ half-lives sit under 2×10^24 y","134Xe two-neutrino decay may be next-gen accessible"]},"model":"grok-4.5","effort":"low","cost_usd":0.004789,"raw_usage":{"total_tokens":1352,"prompt_tokens":722,"num_sources_used":0,"completion_tokens":91,"cost_in_usd_ticks":47888000,"prompt_tokens_details":{"text_tokens":722,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":539,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":722,"tokens_out":91,"duration_ms":10072,"temperature":1.0,"reasoning_tokens":539,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T07:55:54.452150+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A measured half-life (or a firm lower limit tighter than ~10^24 y) for the two-neutrino double-beta decay of 134Xe that falls outside the overlapping theoretical window reported by the four methods.","supporting_citations":[],"review_version":1}