{"id":"5622fb96-5868-4249-a310-ba1809611047","arxiv_id":"2506.23239","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"For 136Xe, the predicted 0νββ half-life is (1.3-3.0)×10^31 y at ⟨mν⟩=1 meV, roughly ten times longer than the standard compilation, based on unifying the effective axial coupling of the 2ν and 0ν modes.","lead":"The authors predict that the neutrinoless double-beta decay half-life of xenon-136 is (1.3-3.0) times 10^31 years for a neutrino mass of 1 meV, about ten times longer than most earlier estimates. If correct, this changes the experimental exposure needed to discover the decay and supports a mild axial-vector quenching in heavy nuclei.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"First-order ratio R≈1.14 is interaction-dependent (SGII gives 0.54); without a second-order check the transferred g_A^eff and the predicted (1.3–3.0)×10^31 y half-life are unsecured.","rationale":"The reader's weakest_assumption identifies the load-bearing premise as the reliability of first-order perturbation theory for the transition operator. My independent reading agrees: the entire half-life prediction hinges on the coefficient R being close to the all-order ratio, and the paper provides no estimate of higher-order corrections. I also note that the interaction dependence undermines the 'universality' claim: SGII gives R = 0.54, so the cancellation of the neutrino potential is not robust. The proposed concrete test—computing second-order corrections—would directly settle whether the first-order ratio is stable. Because the reader's verdict is already CONDITIONAL and my concern reinforces it rather than moving it, the verdict should remain unchanged.","tokens_in":10707,"tokens_out":2061,"duration_ms":24306,"concrete_test":"Compute the second-order perturbative corrections to the 0νββ transition operator (diagrams with two NN-potential insertions) for SkM* using the same QRPA framework of Ref. [31], and recalculate g_A^eff,0ν(ld;pt), g_A^eff,2ν(ld;pt), and R via Eqs. (13)–(15). If R moves by more than 15% from 1.14, the transfer in Eqs. (20)–(21) is not stable and the predicted half-life should be revised or accompanied by a quantitative uncertainty.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on Eq. (15), the near-equality g_A^eff,0ν(ld;pt) ≃ g_A^eff,2ν(ld;pt) for SkM*, which is then used in Eqs. (20)–(21) to transfer the experimentally calibrated 2ν effective coupling to the 0ν mode. This transfer assumes that the first-order perturbative ratio R ≡ g_A^eff,0ν(ld;pt)/g_A^eff,2ν(ld;pt) reproduces the ratio of the true all-order effective couplings. The paper explicitly states that higher-order corrections are ignored, so no evidence is given that the first-order result survives beyond the leading correction. The analytical argument that the neutrino potential cancels in the ratio is approximate and interaction-dependent: the same formalism with SGII gives R = 0.54 (Table III), far from 1.14 for SkM*. Thus the 'closeness' is not a universal feature of the operator but a property of the selected interaction, and the selection of SkM* as reliable is justified by binding-energy systematics rather than by convergence of the ratio. The convergence claim based on Eq. (16) is also weak, since ld;ld is trivially 1 for both modes and provides no information about all-order behavior. If higher-order terms alter R for SkM* significantly, the estimated g_A^eff,0ν values (0.482 and 0.921) and the resulting half-life range (1.3–3.0)×10^31 y would not follow.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper calculates the neutrinoless double-beta-decay (0νββ) half-life of 136Xe using effective axial-vector coupling constants obtained from the two-neutrino double-beta-decay (2νββ) mode. The NMEs include leading-order, vertex-correction, and two-body-current contributions computed in lowest-order perturbation theory with Skyrme interactions SkM* and SGII. The key observation is the near equality g_A,0ν^eff(ld;pt) ≃ g_A,2ν^eff(ld;pt) for SkM* (Eq. 15), which is used in Eqs. (20)–(21) to transfer the experimentally calibrated 2νββ effective coupling to the 0νββ mode. The authors predict T_{1/2}^{0ν} = (1.3–3.0)×10^31 y for ⟨mν⟩ = 1 meV, about an order of magnitude longer than the compilation of Ref. [27].","tokens_in":11148,"tokens_out":5662,"duration_ms":59678,"significance":"If the central premise is valid, the paper provides a falsifiable prediction for the 0νββ half-life of 136Xe and supports the idea of a common effective axial coupling for 2νββ and 0νββ modes, consistent with recent claims of g_A^eff ≈ 1. The inclusion of vertex corrections and two-body currents is a step beyond leading-order treatments, and the explicit use of the measured 2νββ half-life as a non-perturbative calibration is a useful idea. The predicted half-life range is directly testable by upcoming 136Xe experiments. However, the significance is conditional on the validity of the first-order perturbative transfer, which is not yet established.","major_comments":[{"comment":"The load-bearing step of the paper is the transfer of the experimentally calibrated 2νββ effective coupling to the 0νββ mode. This transfer depends on the ratio R = g_A,0ν(ld;pt)/g_A,2ν(ld;pt), which is 1.14 for SkM* but 0.54 for SGII. The paper argues that SkM* is more reliable because of binding-energy systematics, but that argument concerns the overall interaction strength, not the convergence or universality of the ratio R. Since the closeness in Eq. (15) is not a robust property across the two interactions, the paper needs a second-order perturbative calculation or at least an explicit estimate of the neglected higher-order terms before the half-life band (1.3–3.0)×10^31 y can be considered reliable.","section":"Eqs. (20)–(21), Table III"},{"comment":"The claimed convergence of g_A,0ν/g_A,2ν toward unity is not demonstrated. Eq. (16) states that g_A,0ν(ld;ld) = g_A,2ν(ld;ld) = g_A^bare, which is true by definition of the leading-order calculation; it contains no dynamical information about convergence. The first-order values being close for SkM* is a single data point. Table I shows that the perturbative corrections are large: for SkM*, the two-body-current GT term for 0νββ is -2.731 compared with the leading term 3.095, and the 2νββ GT matrix element changes sign under perturbation (0.102 to -0.035). A perturbative series with corrections of this size cannot be assumed to have converged at first order.","section":"Table I and Eq. (16)"},{"comment":"The 'stability' of Comparisons 3–5 is presented as evidence of convergence, but the spread among these three SkM* half-lives is a factor of about 2.4 (127, 140, and 304 in units of 10^29 y), and the selection of these comparisons as most reliable is made after excluding Comparisons 1 and 2. The analogous SGII results do not cluster at all (337, 3990, and 4210 in units of 10^29 y). A quantitative convergence criterion is needed, and the paper should explain why the SGII failure does not also undermine the SkM* result beyond the binding-energy argument.","section":"Table IV and Fig. 2"},{"comment":"The final 0νββ half-life is, through Eq. (20), essentially proportional to the experimental 2νββ half-life, because g_A,2ν(ld;exp) is defined to reproduce that measured value. This is a legitimate calibration strategy, but it means the predictive content of the method rests entirely on the correction factor R and the NME ratios. The paper should state this dependence more explicitly and should propagate the uncertainty in R—including its strong interaction dependence—into an error estimate for the final half-life, rather than selecting one interaction and one cluster of comparisons.","section":"Eq. (20) and surrounding text"}],"minor_comments":[{"comment":"The phrase 'close to the convergence at the first-order perturbation' is unclear; the authors likely mean that the series appears to be converging, but the sentence should be rewritten to say what is actually being compared.","section":"Eq. (16) and following sentence"},{"comment":"The caption contains the sentence 'these results were taken from Ref. [31]' with a lowercase initial letter; this should be 'These results were taken from Ref. [31].'","section":"Table II caption"},{"comment":"The notation 'ld' and 'pt' in g_A^eff is explained in the text, but a one-line definition in the table caption would improve readability for the reader who jumps directly to the table.","section":"Table II"},{"comment":"The figure would benefit from an explicit statement that the filled and open symbols correspond to SkM* and SGII, respectively, and that the y-axis is logarithmic.","section":"Fig. 2 caption"},{"comment":"The sentence 'If ⟨mν⟩ = 10 meV, the half-lives are two orders of magnitude shorter' is correct because T_{1/2} ∝ 1/⟨mν⟩^2, but it could be clarified that this applies to all entries in Table IV.","section":"after Eq. (17)"}],"recommendation":"major_revision","confidential_remarks":"The paper is essentially an application of the previous result in Ref. [31], and the referee should verify that the overlap with that article is not excessive for a letter. The comparison with the works of Rho (Refs. [29,30]) is presented as supportive evidence, but it is more of an analogy than a rigorous connection; this should not be used to bolster the perturbative-convergence claim. The main technical weakness—absence of a second-order check of the ratio R—is, in principle, addressable within the manuscript's scope, which is why I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the half-life range (1.3–3.0)×10^31 y at ⟨mν⟩=1 meV for 136Xe, about a factor of ten longer than the Engel–Menéndez compilation value. The correction-factor scheme in Eqs. (20)–(21) is a transparent way to transfer the experimentally calibrated 2ν coupling to the 0ν mode, and the paper compares its result with recent independent calculations, where it overlaps. It also states its own limitations plainly: the NMEs come from an earlier paper, and the authors admit higher-order corrections are ignored. That is useful and honest.\n\nThe soft spot is load-bearing. The transfer assumes that the first-order perturbative ratio R = g_A,0ν(ld;pt)/g_A,2ν(ld;pt) ≈ 1.14 for SkM* equals the ratio of the true all-order couplings. The paper gives no estimate of the higher-order terms, and the perturbative corrections are not small: for SkM* the 2bc GT term (−2.731) nearly cancels the leading term (3.095), and the perturbed 2ν GT NME flips sign. The interaction dependence is also a problem. The same formalism with SGII gives R = 0.54, so the closeness is not a universal property of the operator; it is a property of the chosen interaction. The convergence argument based on Eq. (16) is weak, because the leading-order ratio is trivially 1 and says nothing about all-order behavior. Since the half-life scales as g_A^−4, a 20% error in g_A moves the prediction by roughly a factor of two, and the unquantified R is what makes the half-life range unsecured.\n\nThere is also a circularity smell, though not a fatal one. Up to NME ratios and the correction factor, the final half-life reduces to the experimental 2νββ half-life. That is not automatically wrong, but it means the prediction inherits both the measured value and the assumption that R is universal. The paper's use of binding-energy systematics to prefer SkM* is a reasonable choice, but it does not establish convergence of R.\n\nWho is this for? Experimentalists planning 136Xe exposures (KamLAND-Zen, nEXO, NEXT) and theorists working on g_A quenching. It is a short paper with a clear, falsifiable prediction, and it engages with the literature. I would send it to peer review rather than desk reject it, with the request to either supply a second-order estimate of the corrections or otherwise show that R is robust, and to put an uncertainty band on the final half-life range. Without that, the central claim is a plausible hypothesis, not a demonstrated result.","headline":"A new half-life prediction for 136Xe 0νββ that deserves a careful referee but rests on an unquantified first-order transfer of g_A between the two decay modes.","tokens_in":11703,"tokens_out":2447,"would_cite":true,"duration_ms":25597,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["23.40.-s","21.60.Jz","14.60.Pq"],"model":"deepseek-v4-flash","headline":"This paper predicts that the neutrinoless double-beta decay half-life of 136Xe is (1.3–3.0)×10^31 years for a Majorana neutrino mass of 1 meV, about an order of magnitude longer than most earlier calculations.","keywords":["neutrinoless double-beta decay","two-neutrino double-beta decay","effective axial-vector coupling","136Xe","nuclear matrix elements","quasiparticle random-phase approximation","Skyrme energy density functional","Majorana neutrino mass"],"falsifier":"A direct measurement or stronger experimental limit on the $0\\nu\\beta\\beta$ half-life of $^{136}$Xe, combined with an independent determination of the neutrino mass scale, would settle the $(1.3{-}3.0)\\times10^{31}$ y window, since a true half-life far outside it at $\\langle m_\\nu\\rangle = 1$ meV would falsify the ratio transfer. Theoretically, recomputing the matrix elements with the perturbation summed beyond lowest order would test whether $R$ stays near 1.14; the SGII value 0.54 already shows the ratio is not interaction-independent, so an all-order SkM* calculation is a direct arbiter.","tokens_in":10428,"feed_emoji":"⚛️","tokens_out":13328,"duration_ms":118142,"temperature":0.7,"pith_summary":"This paper aims to pin down the half-life of the neutrinoless double-$\\beta$ decay ($0\\nu\\beta\\beta$) of $^{136}$Xe, the decay whose observation would prove that the neutrino is its own antiparticle. Its strategy is to take the effective axial-vector coupling constant ($g_A^{\\rm eff}$) that reproduces the measured two-neutrino double-$\\beta$ half-life and transfer it to the neutrinoless mode, using the ratio of perturbed to leading-order matrix elements computed for each mode. The central claim is that $g_A^{\\rm eff}$ is nearly the same for the two modes — about 1.14 for the SkM* interaction — because the neutrino potential enters both the numerator and the denominator of the ratio and approximately cancels. With this unified coupling the paper predicts $T_{1/2}^{0\\nu} = (1.3{-}3.0)\\times10^{31}$ y for an assumed Majorana mass of $\\langle m_\\nu\\rangle = 1$ meV, about an order of magnitude longer than the older compilation range it is measured against. If the prediction is right, the half-life that experiments must reach is longer and the neutrino-mass bounds drawn from any given experimental limit are correspondingly weaker.","feed_headline":"A 1.14 ratio sets the 136Xe 0νββ half-life near 10^31 years","feed_subtitle":"The unified g_A coupling lengthens the predicted half-life tenfold, weakening the mass limits from experiments.","key_machinery":"The load-bearing object is the ratio $R = g_{A,0\\nu}^{\\rm eff}({\\rm ld;pt})/g_{A,2\\nu}^{\\rm eff}({\\rm ld;pt})$, built from the effective couplings that make each mode's leading-order matrix element reproduce its perturbed matrix element; the perturbed matrix elements include the leading Gamow-Teller and Fermi terms, vertex corrections, and the two-body current correction shown in Fig. 1. For the SkM* Skyrme interaction the ratio is 1.14, while the less realistic SGII interaction gives 0.54. The transfer equations (20) and (21) multiply the experimentally calibrated two-neutrino couplings $g_{A,2\\nu}^{\\rm eff}({\\rm ld;exp})$ and $g_{A,2\\nu}^{\\rm eff}({\\rm pt;exp})$ by $R$ to obtain the estimated neutrinoless-mode couplings, a non-perturbative step because the experimental half-life carries the higher-order physics that the one-step perturbed matrix elements miss.","core_discovery":"The central discovery is that the effective axial-vector couplings of the two double-$\\beta$ decay modes nearly coincide: $g_{A,0\\nu}^{\\rm eff}({\\rm ld;pt}) \\simeq g_{A,2\\nu}^{\\rm eff}({\\rm ld;pt})$ for the SkM* interaction. The authors compute the $0\\nu\\beta\\beta$ and $2\\nu\\beta\\beta$ matrix elements of $^{136}$Xe in the quasiparticle random-phase approximation, perturbing the decay operator once by the nucleon-nucleon potential through vertex corrections and a two-body current term. For each mode an effective coupling is defined as the value that makes the leading-order matrix element reproduce the perturbed one, and they find the two effective couplings are close even though the neutrinoless mode's virtual neutrino can carry arbitrarily high momentum. They trace this to an approximate cancellation: the neutrino potential sits in both the numerator and the denominator of the ratio, so it factorizes and drops out. That near-equality licenses the paper's main move — using the experimentally calibrated two-neutrino coupling, rescaled by the ratio 1.14, as the coupling for the neutrinoless mode. Five such estimation methods yield half-lives that cluster, and the paper concludes that the reliable $0\\nu\\beta\\beta$ half-life of $^{136}$Xe is $(1.3{-}3.0)\\times10^{31}$ y at $\\langle m_\\nu\\rangle = 1$ meV.","pith_inferences":["The ratio-transference scheme is portable in principle to other double-beta emitters such as $^{76}$Ge or $^{130}$Te, but because the paper shows the ratio is interaction-dependent (1.14 for SkM* versus 0.54 for SGII), each nucleus and each energy functional would need its own validation before its half-life is trusted.","A sharp cross-framework test would be to compute the same ratio $R$ in shell-model or ab initio frameworks that include two-body currents; if $R\\approx 1$ survives there, the universality of the effective axial-vector coupling is a genuine nuclear-physics fact rather than a feature of this QRPA-plus-Skyrme setup.","The calibration rides on the measured two-neutrino half-life, so more precise $2\\nu\\beta\\beta$ measurements, including spectral shapes, would either tighten or destabilize the predicted $0\\nu\\beta\\beta$ half-life.","An independent determination of the in-medium axial coupling, for instance from super-allowed Gamow-Teller decays, would arbitrate between the SkM* branch (which supports $g_A^{\\rm eff}\\approx 1$) and the SGII branch of the argument."],"forward_implications":["The $0\\nu\\beta\\beta$ half-life of $^{136}$Xe is $(1.3{-}3.0)\\times10^{31}$ y at $\\langle m_\\nu\\rangle = 1$ meV, an order of magnitude longer than the $(3{-}30)\\times10^{29}$ y compilation of older calculations it is compared with.","The same effective axial-vector coupling, near $g_A^{\\rm eff}\\approx 1$, works for both decay modes, so the neutrinoless mode inherits the calibration of the measured two-neutrino mode instead of an ad hoc quenching.","The predicted half-life is stable across the five extraction methods (Comparisons 3–5 for SkM*), so the spread seen among earlier predictions is attributed to using the bare coupling with under-corrected matrix elements.","Because the predicted half-life is longer, any given experimental limit on $0\\nu\\beta\\beta$ decay converts to a weaker upper bound on the Majorana neutrino mass.","The predicted range partially overlaps recent quenched shell-model and chiral effective-field-theory results at the $10^{31}$ y scale, so the lengthening is not isolated to this calculation."],"supporting_citations":[{"why":"Supplies the perturbed 0ν and 2ν matrix elements, the effective couplings of Table II, and the analytical argument that the neutrino potential cancels in the ratio.","marker":"[31]"},{"why":"Provides the experimental 2νββ half-life (2.18×10^21 y) that fixes g_A,2ν^eff and calibrates the 0ν transfer.","marker":"[35]"},{"why":"The compilation of 0νββ half-lives in the (3–30)×10^29 y range that sets the baseline the paper's prediction is compared against.","marker":"[27]"},{"why":"Defines the SkM* Skyrme interaction used for the realistic nuclear-structure calculation.","marker":"[32]"},{"why":"Defines the SGII interaction used to show the method's dependence on the force, yielding the smaller ratio 0.54.","marker":"[33]"},{"why":"The QRPA framework used to compute the nuclear wave functions of 136Xe.","marker":"[34]"},{"why":"Provides the Skyrme-QRPA formulation with pairing interactions on which the NME calculation is built.","marker":"[17]"},{"why":"Shell-model NMEs with and without quenching (half-lives 0.1×10^31 y and 3.8×10^31 y) that overlap the paper's range.","marker":"[9]"},{"why":"Independent estimate of g_A^eff ≈ 1 from nuclear-matter arguments, cited as corroboration of the unquenched value.","marker":"[30]"}],"fun_headline_variants":["136Xe 0νββ half-life pinned at ~2×10^31 years","Unified g_A for both ββ modes sets 136Xe half-life near 10^31 y","Two ββ modes share g_A; 136Xe 0νββ half-life converges to ~10^31 y","Ratio of g_A for 2ν and 0ν ββ fixes 136Xe half-life at ~10^31 y","Closeness of g_A links two ββ modes, setting 136Xe T1/2 near 10^31 y"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that one round of perturbation on the decay operator is a faithful guide to the full physics, so that the SkM* ratio 1.14 equals the ratio of the true non-perturbative effective couplings; the paper states that higher-order corrections are ignored, and the SGII result of 0.54 shows the ratio is sensitive to the interaction.","fun_headline_variants_meta":{"raw":{"variants":["136Xe 0νββ half-life pinned at ~2×10^31 years","Unified g_A for both ββ modes sets 136Xe half-life near 10^31 y","Two ββ modes share g_A; 136Xe 0νββ half-life converges to ~10^31 y","Ratio of g_A for 2ν and 0ν ββ fixes 136Xe half-life at ~10^31 y","Closeness of g_A links two ββ modes, setting 136Xe T1/2 near 10^31 y"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000977,"raw_usage":{"total_tokens":4267,"prompt_tokens":1180,"completion_tokens":3087,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":796,"completion_tokens_details":{"reasoning_tokens":2940}},"tokens_in":796,"tokens_out":3087,"duration_ms":23114,"temperature":1.0,"reasoning_tokens":2940,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:50:03.833149+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement or stronger experimental limit on the $0\\nu\\beta\\beta$ half-life of $^{136}$Xe, combined with an independent determination of the neutrino mass scale, would settle the $(1.3{-}3.0)\\times10^{31}$ y window, since a true half-life far outside it at $\\langle m_\\nu\\rangle = 1$ meV would falsify the ratio transfer. Theoretically, recomputing the matrix elements with the perturbation summed beyond lowest order would test whether $R$ stays near 1.14; the SGII value 0.54 already shows the ratio is not interaction-independent, so an all-order SkM* calculation is a direct arbiter.","supporting_citations":[{"cited_title":"Terasaki and O","cited_arxiv_id":null,"evidence_quote":"Supplies the perturbed 0ν and 2ν matrix elements, the effective couplings of Table II, and the analytical argument that the neutrino potential cancels in the ratio."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental 2νββ half-life (2.18×10^21 y) that fixes g_A,2ν^eff and calibrates the 0ν transfer."},{"cited_title":"Bartel, P","cited_arxiv_id":null,"evidence_quote":"Defines the SkM* Skyrme interaction used for the realistic nuclear-structure calculation."},{"cited_title":"Van Giai and H","cited_arxiv_id":null,"evidence_quote":"Defines the SGII interaction used to show the method's dependence on the force, yielding the smaller ratio 0.54."},{"cited_title":"Terasaki and J","cited_arxiv_id":null,"evidence_quote":"The QRPA framework used to compute the nuclear wave functions of 136Xe."},{"cited_title":"Terasaki, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the Skyrme-QRPA formulation with pairing interactions on which the NME calculation is built."},{"cited_title":"Coraggio, A","cited_arxiv_id":null,"evidence_quote":"Shell-model NMEs with and without quenching (half-lives 0.1×10^31 y and 3.8×10^31 y) that overlap the paper's range."},{"cited_title":"From Nuclear Matter with Quenched $g_A$ to Compact-Star Matter with a Signal for Emergent Hidden Scale Symmetry","cited_arxiv_id":"2507.04939","evidence_quote":"Independent estimate of g_A^eff ≈ 1 from nuclear-matter arguments, cited as corroboration of the unquenched value."}],"review_version":1}