{"id":"6c60be75-dba4-4f97-9ebe-7be62c4185fc","arxiv_id":"2505.09679","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A comprehensive scan of the 3+3 Type-I seesaw model finds that current and next-generation neutrinoless double beta decay experiments have broad discovery potential in both normal and inverted neutrino mass ordering.","lead":"This paper calculates the chances that upcoming neutrinoless double beta decay experiments will detect the Type-I seesaw model, a minimal theory of neutrino masses, and finds broad discovery potential in both neutrino mass orderings. The result matters because it tells experimentalists which parts of this well-motivated theory their detectors can probe or rule out, even in regions where the usual light-neutrino analysis sees nothing.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Uncertainty in the nuclear matrix elements and short-range LEC is not propagated into the half-life profile, so the quantitative boundary of the claimed 0νββ discovery potential is not yet established.","rationale":"The paper's central qualitative result, that the 3+3 Type-I seesaw can produce 0νββ rates beyond the high-scale seesaw expectation, is well motivated and consistent with earlier restricted studies, and the EFT framework of Ref. [26] provides a credible basis. The soft spot is not internal inconsistency but the unpropagated theory error in the nuclear and short-distance inputs. Because the discovery-potential statement is quantitative, referring to specific half-life windows and profile-likelihood boundaries, those uncertainties are load-bearing: they can move the T1/2 boundary by the square of the input error. The reader's CONDITIONAL verdict already captures this; no adjustment is needed, but a sensitivity study should be a condition of acceptance.","tokens_in":11177,"tokens_out":11480,"duration_ms":121901,"concrete_test":"Recompute the profile likelihood in Fig. 5 for both isotopes while varying M_ν^(3) over the published spread (e.g., 1.5-4.5 for 136Xe) and varying gNNν(mi) by ±50% and with the opposite sign, keeping all other scan inputs fixed. If the 95% C.L. interval in log10 T1/2 still overlaps the 1e26-1e28 yr window and still extends above the light-exchange band in every case, the concern is not load-bearing; if the overlap or the extension disappears in any case, the paper must quote the resulting range of reach or qualify the claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eqs. (9)-(11) convert the computed 0νββ amplitude into a half-life and then into |meff|. The normalization is fixed by the light-neutrino NME M_ν^(3)=2.7 (136Xe) and 3.4 (76Ge) and by the short-distance LEC gNNν(mi) from Ref. [29]. These inputs enter the rate quadratically: T1/2 ∝ |Σ U^2 m_i Aν(m_i)|^-2, so a multiplicative uncertainty in Aν shifts the T1/2 profile by the inverse square of that uncertainty. Published NME values for these isotopes differ by factors of roughly 2-3 depending on the many-body method, and no uncertainty is assigned to gNNν; the paper uses central values only. The 95% C.L. boundary in Fig. 5, and thus the portion of the (M1, |meff|) plane that overlaps the 1e26-1e28 yr experimental window, is conditional on these choices. If the true NME is smaller, or if the short-range sterile contribution is weaker or interferes destructively at a different level, the 'broad' tail of the profile can move by orders of magnitude, and the precise meaning of 'current and next-generation discovery potential' changes. This does not invalidate the model setup or the scan, but it means the quantitative reach claim is not yet robust.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the 3+3 Type-I seesaw extension of the Standard Model ('nuSM') with three right-handed neutrinos whose masses range from keV to 10 TeV. It performs an 18-dimensional parameter scan using the GAMBIT framework, incorporates updated constraints from neutrino oscillation data, electroweak precision observables, CKM unitarity, rare decays, cosmology (DESI 2024), direct sterile-neutrino searches, and BBN (with and without), and computes the neutrinoless double beta decay (0νββ) half-life using the EFT formalism of Ref. [26]. Results are presented as profile-likelihood plots for |meff_ββ| and T1/2 for 76Ge and 136Xe, and the central claim is that current and next-generation 0νββ experiments have broad discovery potential in both normal and inverted neutrino mass orderings.","tokens_in":11401,"tokens_out":5024,"duration_ms":58285,"significance":"If the central claim holds, this would be a useful and fairly comprehensive phenomenological map of the seesaw parameter space relevant to 0νββ, updating earlier N=2 and restricted N=3 studies. The paper's strengths include the use of a public, well-tested scanning framework (GAMBIT), a validated parametrization of the seesaw mixing matrix, an up-to-date set of experimental constraints, and a state-of-the-art EFT treatment of 0νββ from Ref. [26]. The predicted half-life profiles are falsifiable and directly relevant to current and next-generation experiments. The main weakness is that the quantitative reach claim is conditional on fixed nuclear matrix elements and a short-range low-energy constant whose uncertainties are not propagated.","major_comments":[{"comment":"The half-life and |meff_ββ| profiles are computed with central values only for the nuclear matrix elements M_ν^(3) = 2.7 (136Xe) and 3.4 (76Ge) and for the short-range LEC gNN_ν(m_i) from Ref. [29]. Because T1/2 scales as the inverse square of the total 0νββ amplitude, a factor-of-2 to 3 uncertainty in the NME translates into a factor-of-4 to 9 (roughly an order of magnitude) shift in the T1/2 profile. The location of the 95% C.L. boundary in Fig. 5 directly controls the overlap with the 10^26–10^28 yr experimental window, so the claim of broad current and next-generation discovery potential is conditional on these unpropagated inputs. The authors should propagate these uncertainties or, at minimum, show how the boundaries in Figs. 5 and 6 move under a conservative variation of M_ν^(3) and gNN_ν.","section":"Sec. 4, Eqs. (9)–(11); Figs. 5–6"},{"comment":"The 1D profile likelihood for T1/2 appears flat at Λ = 1 over many orders of magnitude (e.g., from about 10^20 to 10^40 yr in the no-BBN cases), and this flat region is the basis for describing the discovery potential as 'broad'. Because the scan is a random differential-evolution run in an 18-dimensional parameter space, the apparent plateau could in principle reflect incomplete coverage rather than a genuine likelihood maximum. Please provide convergence diagnostics (for example, reproducibility of L_max across independent runs or comparison with a second scanning method) or otherwise justify that the flat profile is physical. If the plateau is robust, the central claim is strengthened; if it is a sampling artifact, the discovery-potential statement would need to be qualified.","section":"Sec. 5, Fig. 5"}],"minor_comments":[{"comment":"The sign convention for Δm²_3ℓ is only specified in the table caption ('with '+' and ℓ=1 ... and '−' and ℓ=2'). It would be clearer to state the NH/IH assignment explicitly in the table itself or in a dedicated footnote.","section":"Table 1"},{"comment":"Eq. (11) defines |meff_ββ| using the light-neutrino NME M_ν^(3), while the actual half-life in Eq. (9) receives sterile-neutrino contributions through A_ν(m_i). The text should state more explicitly that |meff_ββ| is a convenient rescaling of T1/2 for presentation purposes and not the physical amplitude, to avoid confusion for readers who identify |meff| with the standard light-neutrino effective mass.","section":"Sec. 4, Eq. (11)"},{"comment":"The phrase 'bounds predicted by the light neutrino exchange mechanism' is slightly misleading, since the dashed lines are envelopes of the light-neutrino-only region rather than sharp bounds. Consider replacing 'bounds' with 'allowed region' or 'envelope'.","section":"Captions of Figs. 1–4"},{"comment":"The paper should state the GAMBIT version used and provide or link the scanner configuration files (including likelihood definitions and parameter ranges) so that the 18-dimensional profile likelihood can be reproduced by independent groups.","section":"Reproducibility"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for a hep-ph journal and is a solid phenomenological analysis. The main issue is the unpropagated nuclear-theory uncertainty, which affects the quantitative boundary of the central claim but is addressable with additional robustness checks. I would be willing to reconsider after a revision that addresses the two major points."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take on Cirigliano, Dekens, and Urrutia Quiroga, arXiv:2505.09679. It is a genuinely comprehensive GAMBIT-driven scan of the 3+3 Type-I seesaw over sterile masses from keV to 10 TeV, with updated constraints (DESI 2024, LHCb RK/RK*, CMS W decays, direct beta-decay limits) and the state-of-the-art EFT treatment of 0νββ from the same group. The main qualitative result—that current and next-generation 0νββ experiments can probe broad regions of seesaw parameter space well beyond the light-neutrino-exchange-only scenario, in both normal and inverted ordering—holds up. It is consistent with earlier, narrower studies, and the paper does a real service by making that case on a broader, updated footing.\n\nWhat the paper does well: exact parametrization of the seesaw without a small-mixing expansion, profile likelihoods with Diver, careful BBN-on/BBN-off separation, and a clear presentation of the statistical framework. The 18-dimensional scan is argued with appropriate caution about convergence.\n\nSoft spots, in order of weight. First, the quantitative boundary of the discovery potential is controlled by nuclear matrix elements and the short-range LEC, and the paper uses central values only. The NME uncertainty of a factor 2–3 shifts T1/2 by the inverse square, i.e., up to an order of magnitude, which does not kill the qualitative claim but makes the 95% C.L. contours in Fig. 5 and their overlap with the 1e26–1e28 yr window conditional on those inputs. Second, no likelihood tables, code, or scan configurations are released, so the profile likelihoods are not independently reproducible without reconstructing GAMBIT modules. Third, the ad hoc scan ranges (e.g., R-matrix parameters in [-10,10]) and the EWPO tension treatment are not stress-tested; the authors do not explore sensitivity to these choices. None of this undermines the main conclusion, but it does mean the paper's strongest statement—broad discovery potential—is best read as a qualitative result with uncertain quantitative boundaries.\n\nWho this is for: anyone working on seesaw phenomenology, 0νββ experimental strategy, or global fits. It deserves a serious referee. The referee should ask the authors to propagate or at least bound the NME/LEC uncertainty, release the scan inputs and outputs, and add a short sensitivity statement for the scan-range choices. I would trust the central qualitative claim in the meantime, but not the precise reach numbers.","headline":"Comprehensive νSM global scan with a credible central claim about 0νββ reach, though the numeric boundaries depend on unpropagated nuclear uncertainties and the results are not yet independently reproducible.","tokens_in":12116,"tokens_out":2757,"would_cite":true,"duration_ms":29257,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.60.Pq","23.40.-s","12.60.-i"],"model":"deepseek-v4-flash","headline":"In the minimal 3+3 Type-I seesaw model, current and next-generation neutrinoless double beta decay searches can discover the model in either mass ordering.","keywords":["neutrinoless double beta decay","Type-I seesaw","sterile neutrinos","effective Majorana mass","normal and inverted mass ordering","profile likelihood","nuclear matrix elements","lepton number violation"],"falsifier":"Compute the $0\\nu\\beta\\beta$ nuclear matrix elements for $^{136}$Xe and $^{76}$Ge and the short-range four-nucleon constant $g_{\\rm NN}^{\\nu}$ with controlled uncertainties; if, for example, the $^{136}$Xe matrix element turned out to be 1.0 rather than 2.7, the predicted half-lives in the sterile-dominated region would lengthen by roughly a factor of seven, pushing the 95% C.L. discovery boundary beyond next-generation ton-scale sensitivity. On the experimental side, a $0\\nu\\beta\\beta$ signal with $|m_{\\rm eff}^{\\beta\\beta}|$ above the maximum allowed by light-neutrino exchange in the relevant ordering would confirm that sterile neutrinos contribute.","tokens_in":10839,"feed_emoji":"⚛️","tokens_out":10719,"duration_ms":102518,"temperature":0.7,"pith_summary":"This paper tries to establish that the simplest viable seesaw extension of the Standard Model—three active neutrinos plus three sterile neutrinos, with sterile masses spanning keV to 10 TeV—predicts neutrinoless double $\\beta$ decay half-lives far outside the band set by light-neutrino exchange alone. If true, current experiments probing half-lives near $10^{26}$ years and next-generation ton-scale searches near $10^{28}$ years have a real chance of discovering the model, in both normal and inverted neutrino mass ordering. This matters because neutrinoless double $\\beta$ decay is the most direct practical probe of lepton-number violation and of the Majorana nature of neutrinos, and a discovery or exclusion would discriminate between the high-scale seesaw picture and the richer low-scale variant studied here.","feed_headline":"Sterile neutrinos widen double-beta decay discovery window","feed_subtitle":"In the minimal 3+3 seesaw model, half-lives within reach of current and next-generation searches are allowed in both mass orderings.","key_machinery":"The central object is the mass-dependent amplitude $A_\\nu(m_i)$ in Eq. (10), which splits the contribution of each neutrino mass eigenstate into three momentum regimes: for $m_i<100$ MeV, potential, hard, and ultrasoft pieces; for $100$ MeV $\\le m_i < 2$ GeV, potential and hard pieces; and for $m_i\\ge 2$ GeV, a dimension-nine operator controlled by the short-range low-energy constant $g_{\\rm NN}^{\\nu}$. The inverse half-life in Eq. (9) is proportional to the square of the weighted sum of these six amplitudes, so cancellations between active and sterile contributions can suppress the rate while individual sterile contributions can greatly enhance it. The paper also defines the effective Majorana mass $|m_{\\rm eff}^{\\beta\\beta}|$ in Eq. (11), which reduces to $m_{\\beta\\beta}$ in the high-scale seesaw limit; the deviation of this quantity from the light-exchange-only band is the signature that sterile neutrinos participate in the decay. The statistical machinery is a frequentist profile likelihood over the full parameter space, which converts the scan into confidence regions for $|m_{\\rm eff}^{\\beta\\beta}|$ and $T_{1/2}^{0\\nu}$.","core_discovery":"The paper's central claim is that in the 3+3 Type-I seesaw model (the $\\nu$SM), the neutrinoless double $\\beta$ decay rate is not governed by the light-neutrino-only effective mass $|m_{\\beta\\beta}|$ once sterile neutrinos are light enough to participate. The full amplitude sums contributions from all six Majorana neutrino mass eigenstates, with the sterile contribution entering through mass-dependent momentum regions: potential, hard, and ultrasoft. As a result, the allowed values of the effective Majorana mass $|m_{\\rm eff}^{\\beta\\beta}|$ and the half-life $T_{1/2}^{0\\nu}$ spread over many orders of magnitude, extending well beyond the high-scale seesaw band. The profile likelihoods, built from an 18-dimensional scan subject to a broad set of low- and high-energy constraints, show 68% and 95% confidence regions reaching half-lives that current and next-generation experiments can test, in both normal and inverted ordering. This is what the authors call the model's broad discovery potential.","pith_inferences":["If the model is right, a future neutrinoless double beta decay signal cannot be unambiguously interpreted as pure light-Majorana-neutrino exchange: any value of $|m_{\\rm eff}^{\\beta\\beta}|$ above the light-exchange ceiling would be direct evidence for sterile neutrinos, while a value inside the light band would still allow a sterile component hidden by destructive interference.","The quantitative reach boundaries are conditional on the assumed nuclear matrix elements and the short-range low-energy constant; because the half-life scales as the inverse square of these inputs, improved nuclear-structure calculations could move the 95% C.L. edges by a factor of several, so the reach statements should be revisited as those uncertainties shrink.","If future cosmology closes the entropy-production or axion-like-particle escape routes that avoid big-bang nucleosynthesis bounds, the no-BBN scenarios that open up the light-sterile and short-half-life region would be disfavored, concentrating the discovery potential in the heavier sterile mass range."],"forward_implications":["In the $\\nu$SM, the 95% C.L. region for $|m_{\\rm eff}^{\\beta\\beta}|$ extends well beyond the band allowed by light-neutrino exchange alone, in both mass orderings and with or without big-bang nucleosynthesis constraints.","Current experiments sensitive to half-lives near $10^{26}$ years and next-generation ton-scale searches near $10^{28}$ years already overlap a substantial part of the high-likelihood region, so a discovery is possible in either ordering.","A signal with $|m_{\\rm eff}^{\\beta\\beta}|$ above the maximum allowed by light-neutrino exchange would indicate sterile-neutrino contributions, and the correlation with the lightest sterile mass $M_1$ could help identify the mass scale responsible.","Imposing big-bang nucleosynthesis constraints narrows the allowed parameter space, especially for light sterile neutrinos, and shifts the profile likelihood toward longer half-lives, making the discovery potential more conservative but still present.","The same framework predicts that long half-lives above $10^{28}$ years remain allowed, so a null result at ton-scale sensitivity would not rule the model out but would erode the discovery-potential claim."],"supporting_citations":[{"why":"Provides the effective field theory decomposition of the 0$\\nu\\beta\\beta$ amplitude into potential, hard, and ultrasoft neutrino contributions across the sterile mass range.","marker":"[26]"},{"why":"Fixes the short-distance four-nucleon low-energy constant $g_{\\rm NN}^{\\nu}$ that controls the hard-neutrino contribution for sterile masses above about 100 MeV.","marker":"[29]"},{"why":"Introduces the effective Majorana mass definition and supplies the light-neutrino nuclear matrix elements used in Eq. (11).","marker":"[30]"},{"why":"Gives the earlier global fit of the same model whose constraints and sampling strategy are extended and updated here.","marker":"[12]"},{"why":"Provides the sampling and likelihood machinery used to scan the 18-dimensional parameter space and compute profile likelihoods.","marker":"[13]"},{"why":"Supplies the cosmological upper bound on the sum of active neutrino masses that shapes the allowed $m_{\\nu\\min}$ and $\\sum_i m_{\\nu i}$ regions.","marker":"[20]"},{"why":"Supplies the direct nuclear beta-decay limits on $|U_{ei}|^2$ for sterile neutrinos below about 1 MeV, newly added in this analysis.","marker":"[15]"},{"why":"Define the current experimental half-life exclusions near $10^{26}$ years that mark the present reach in the profile likelihood plots.","marker":"[42, 43]"},{"why":"Define the projected ton-scale sensitivity near $10^{28}$ years that sets the next-generation discovery reach.","marker":"[44–49]"}],"fun_headline_variants":["Sterile neutrinos expand double-beta half-life range","Seesaw model broadens double-beta discovery potential","Double-beta decay sensitivity beyond light neutrinos","New double-beta window from sterile neutrinos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative boundaries of the claimed discovery reach assume the nuclear matrix elements $M_\\nu^{(3)}=2.7$ for $^{136}$Xe and $3.4$ for $^{76}$Ge, together with the short-range constant $g_{\\rm NN}^{\\nu}$, and these inputs are used without propagating their uncertainties; if they differ substantially from the true values, the predicted half-lives shift by the square of that difference and the claimed reach moves with them.","fun_headline_variants_meta":{"raw":{"variants":["Sterile neutrinos expand double-beta half-life range","Seesaw model broadens double-beta discovery potential","Double-beta decay sensitivity beyond light neutrinos","New double-beta window from sterile neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000353,"raw_usage":{"total_tokens":1901,"prompt_tokens":907,"completion_tokens":994,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":523,"completion_tokens_details":{"reasoning_tokens":944}},"tokens_in":523,"tokens_out":994,"duration_ms":9675,"temperature":1.0,"reasoning_tokens":944,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:27:53.162643+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the $0\\nu\\beta\\beta$ nuclear matrix elements for $^{136}$Xe and $^{76}$Ge and the short-range four-nucleon constant $g_{\\rm NN}^{\\nu}$ with controlled uncertainties; if, for example, the $^{136}$Xe matrix element turned out to be 1.0 rather than 2.7, the predicted half-lives in the sterile-dominated region would lengthen by roughly a factor of seven, pushing the 95% C.L. discovery boundary beyond next-generation ton-scale sensitivity. On the experimental side, a $0\\nu\\beta\\beta$ signal with $|m_{\\rm eff}^{\\beta\\beta}|$ above the maximum allowed by light-neutrino exchange in the relevant ordering would confirm that sterile neutrinos contribute.","supporting_citations":[],"review_version":1}