{"id":"51a2edda-4f18-4a2a-b653-cf33678480b4","arxiv_id":"2608.07657","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Renormalization-group mixing lets heavy-quark dimension-seven operators feed neutrinoless double-beta decay, giving some of the strongest current bounds on these new-physics operators.","lead":"The authors calculate how quantum corrections mix hidden new-physics operators into neutrinoless double-beta decay, turning a rare nuclear process into a sharper probe of TeV-scale physics. The result is a set of improved bounds on many dimension-seven operators, in several cases stronger than limits from meson decays.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Neutrino-mass cancellation invoked against Ref. [58] is not applied or demonstrated for the paper's own operators that mix into the Weinberg operator, so the quoted 0νββ limits for those operators may be over-stated.","rationale":"The central claim is that RG mixing makes 0νββ the most stringent probe of many d=7 operators. The cleanest way this could fail is if part of the signal is an artifact of an inconsistent treatment of neutrino masses. The authors themselves identify this failure mode for Ref. [58] in Sec. V, so the concern is internal to the paper's logic, not an outside disagreement. Nuclear-matrix-element uncertainties are real, but they affect all 0νββ bounds roughly multiplicatively and would not change the qualitative mechanism; the neutrino-mass cancellation can remove entire contributions for specific operators, changing the ordering with meson-decay limits. The paper does not state which operators are affected, nor whether the quoted constraints are obtained after cancellation, nor how the cancellation is implemented. This is an addressable consistency check rather than a fatal flaw, so the verdict stays CONDITIONAL; since the reader already reached CONDITIONAL, no verdict change is needed. Independent support: the use of the published one-loop ADM and the public nuDoBe framework is appropriate, and no algebraic error was found in the illustrative examples of Sec. III.","tokens_in":17300,"tokens_out":7190,"duration_ms":74368,"concrete_test":"Take the operator C_dlqlH1^(7) with (i,j)=(3,3), the example named in footnote 8. (1) Compute the one-loop running from Lambda = 30 TeV to m_W using the ADM of Ref. [23] to obtain the induced d=5 Weinberg coefficient and the resulting m_nu; compare with m_nu < 0.1 eV. (2) Recompute the 0νββ limit shown in Fig. 4 under two treatments: (a) no neutrino-mass constraint, and (b) a counterterm tuned to cancel the induced m_nu, dropping the light-neutrino-exchange amplitude while retaining short-range d=7 contributions. If the Lambda limit in treatment (b) differs from the published value by more than ~30%, the Fig. 4 constraint for this operator depends on an unspecified cancellation and the 'most stringent' claim needs qualification. Repeat for all operators flagged as mixing into O_LH^(5/7).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's numerical constraints for operators that mix into the dimension-five Weinberg operator (and into O_LH^(7)) are not derived consistently with the neutrino-mass bound the authors invoke against Ref. [58]. In the final paragraph of Sec. V the authors state that several d=7 operators induce sizable contributions to neutrino masses and that these 'must be canceled out by other EFT contributions' to keep m_nu <~ 0.1 eV. But no such cancellation is implemented or demonstrated for the authors' own Fig. 4 bounds on exactly those operators, e.g., C_dlqlH1 with i=j=3, which footnote 8 names as generating a contribution to the Weinberg operator proportional to y_b. If the dominant 0νββ signal from such an operator is the one-loop RG-generated Weinberg contribution, imposing the neutrino-mass bound removes or reweights that signal, and the quoted Lambda limits and the 'most stringent' comparison to meson decays can change by orders of magnitude for the affected operators. The manuscript thus uses a fine-tuning assumption against a competing tool without showing that its own constraints are invariant under the same assumption.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies neutrinoless double-beta decay (0νββ) constraints on the twelve lepton-number-violating dimension-seven SMEFT operators. The authors use the one-loop anomalous dimension matrix of Ref. [23] to evolve the full set of d=7 Wilson coefficients from the new-physics scale to the electroweak scale, including Yukawa/CKM-induced flavor mixing and keeping first- and second-logarithmic terms in the RG expansion. They feed the resulting low-energy coefficients into the νDoBe 'Master Formula' chain (Refs. [12,13,43]) to compute 0νββ half-lives, and compare the resulting limits on Λ/∛C with limits from pseudoscalar meson decays (P→e+inv, P→P'e^-e^-, P→P'inv). The main results are that RG-induced mixing makes 0νββ the most stringent probe for many ψ4H and ψ4D operators, especially those with heavy-quark flavor indices, with scales up to about 30 TeV, and that for several operators the leading contribution arises at second logarithmic order through two-step mixing. The paper also compares its results with the numerical tool of Ref. [58] and argues that the more stringent limits obtained there for operators that mix into the Weinberg operator rely on unrealistic induced neutrino masses.","tokens_in":17541,"tokens_out":20181,"duration_ms":213002,"significance":"If the results are correct, this is a useful systematic contribution: it extends earlier single-operator studies to the complete d=7 LNV operator set with general quark flavor, identifies operators for which second-logarithmic mixing is the leading effect, and gives concrete predictions for the sensitivity of next-generation 0νββ experiments. Strengths include the use of an independent published anomalous-dimension matrix, the explicit classification of tree-level, first-log, and second-log dominance in Fig. 4, and the use of the public νDoBe pipeline, which makes the numerical analysis reproducible in principle. The comparison with meson decays is a valuable benchmark. The main caveats are the model-dependent treatment of d=7→d=5 Weinberg mixing and the absence of quantitative uncertainty estimates for the nuclear input; these affect the strength of the 'most stringent' claim for a subset of the operators.","major_comments":[{"comment":"The paper excludes, for a class of operators, the contribution obtained by mixing into the dimension-five Weinberg operator, invoking a cancellation of the induced neutrino masses. However, the text does not state which entries in Figs. 3 and 4 are affected by this exclusion, nor does it show that the cancellation required to satisfy mν≲0.1 eV leaves the remaining non-Weinberg contributions unchanged. The concrete case named in footnote 8, C^{(7)}_{dlqlH1} with i=j=3, appears in Fig. 4; if its displayed 0νββ bound had included the y_b-proportional Weinberg contribution, imposing the neutrino-mass bound would remove or strongly reweight that limit. The comparison with Ref. [58] suggests that this contribution was dropped, but the paper should say so explicitly for each affected operator and quantify the resulting change. As written, the 'most stringent' claim for the affected ψ4H and ψ4D entries is not fully supported.","section":"Sec. V, final paragraph and footnote 8; Figs. 3–4"},{"comment":"The abstract and introduction describe the analysis as a 'complete one-loop RG analysis' of 0νββ, but Eq. (3) drops the 'contributions from insertions of lower-dimensional operators' and the numerical analysis does not include the d=7→d=5 mixing that is precisely the subject of the final comparison paragraph. This is a physical choice: one can marginalize over a tuned cancellation of induced neutrino masses, but then the quoted bounds are conditional on that tuning. The paper should either include the Weinberg-mixing terms in the numerical pipeline or qualify the 'complete' claim and state the tuning assumption at the point where νDoBe is adopted (Sec. IV B), rather than only in the closing discussion.","section":"Sec. III, Eq. (3); Sec. IV B; Sec. V"},{"comment":"The quoted limits are central values with no uncertainty bands. The 0νββ half-lives inherit nuclear matrix elements and chiral EFT coefficients from νDoBe; order-one changes in these inputs could reorder the comparison with meson-decay limits for operators where the two probes are close. To support the 'most stringent' claim, the paper should provide an estimate of the dominant nuclear/EFT uncertainties, or identify how large a shift would be needed to change the ordering between 0νββ and meson-decay constraints.","section":"Sec. IV B; Figs. 3–4"}],"minor_comments":[{"comment":"Refs. [24] and [58] are the same paper (Y. Liao et al., JHEP 08 (2025) 138); the duplicate citation should be merged.","section":"References"},{"comment":"There are two typos: 'automized' should be 'automated', and 'out SMEFT analysis' should be 'our SMEFT analysis'.","section":"Sec. IV B and Sec. IV A"},{"comment":"The coefficient C^{(7)}_{dlqlH1} in footnote 8 uses a different notation from the operator O^{(7)}_{\\bar{d}LQLH1} in Table I; please unify the naming convention for this operator.","section":"Table I and footnote 8"},{"comment":"The notation P^-→P'^+ + inv is confusing because both mesons are charged; please define the charge assignments and state explicitly that 'inv' denotes missing neutrinos in the process.","section":"Sec. IV C and Fig. 2 caption"},{"comment":"The Belle-II B^+→K^++inv result is described as a 'mild excess' and then translated into Λ≃3 TeV as though it were an upper bound; please specify whether this is treated as a limit, a two-sided constraint, or a signal, and what confidence level is used.","section":"Sec. V, paragraph after Fig. 4"},{"comment":"The term 'second leading-logarithm' is used for the log² term built from two one-loop anomalous-dimension entries; this should be distinguished explicitly from a genuine two-loop anomalous-dimension contribution, which the conclusion mentions as future work.","section":"Sec. III, Eq. (5); Sec. VI"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is squarely within the scope of a hep-ph journal and the central RG machinery is sound and based on published anomalous dimensions. The main issue is the consistency and transparency of the treatment of operators that mix into the Weinberg operator: the authors invoke a cancellation of induced neutrino masses to explain disagreement with Ref. [58], but they do not identify which of their own Fig. 4 bounds are affected or demonstrate that the cancellation leaves those bounds unchanged. This appears fixable by clarifying the assumption and adding a comparison of limits with and without the Weinberg contribution for the affected operators. The lack of nuclear-uncertainty estimates is a second, more standard caveat. I recommend major revision rather than rejection; the novelty and technical content are adequate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read this one carefully. The new content is real: a complete one-loop RG analysis of all d=7 LNV operators in 0νββ, with general quark flavor, including the subleading-logarithm terms. Earlier work did diagonal running or subsets; this paper systematically maps which operators get constrained through mixing, and identifies flavor combinations where the leading effect appears only at second logarithmic order—e.g., C_duLLD with (3,1). That is a useful result for model builders. The machinery is standard, the ADM is taken from Ref [23], and the numerical rates go through the public νDoBe chain. The representative analytic expressions check out internally.\n\nThe soft spot is the neutrino-mass consistency. In the last paragraph of Sec. V, the authors explain that Ref [58] gets overly stringent bounds on operators that mix into the Weinberg operator because the induced neutrino masses are unrealistic and 'must be canceled out by other EFT contributions.' Fine. But they do not apply that same reasoning to their own Fig. 4 limits on exactly those operators. Footnote 8 even names C_dlqlH1 with i=j=3 as generating a Weinberg contribution proportional to y_b. If the 0νββ signal for that coefficient is dominated by the one-loop RG-generated Weinberg term, then imposing m_nu < 0.1 eV would remove or drastically reweight that signal. The quoted Λ limits and the comparison to meson decays could change by orders of magnitude for those operators. The paper needs to either demonstrate that its constraints are invariant under the same cancellation, or re-derive the affected limits with the neutrino-mass bound imposed. This doesn't sink the whole paper—many operators don't feed the Weinberg operator—but it is a load-bearing issue for a subset of the headline claims.\n\nA separate, minor issue is that the limits are quoted without uncertainty bands despite inheriting nuclear matrix elements and chiral coefficients from νDoBe. That's standard practice in this field, but it means the 'most stringent' ordering against meson decays could shift for some operators if the nuclear input moves by order one.\n\nWho's this for? EFT practitioners and model builders working on TeV-scale LNV. The core methodology is sound and the flavor atlas in Figs. 3–4 is valuable. With the Weinberg-mixing consistency fixed, this is a solid paper. I'd send it to a serious referee, but the revision should be asked to address that point before acceptance.","headline":"Solid d=7 SMEFT RG atlas for 0νββ, but the Weinberg-mixing limits need to be redone under the same neutrino-mass bound they use against a competing tool.","tokens_in":18040,"tokens_out":4274,"would_cite":false,"duration_ms":39381,"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":"Through one-loop renormalization-group mixing, lepton-number-violating operators that contribute nothing at tree level to neutrinoless double-beta decay still acquire constraints reaching scales near 30 TeV.","keywords":["neutrinoless double-beta decay","SMEFT","dimension-seven operators","renormalization group evolution","operator mixing","lepton number violation","leading-logarithm expansion","quark flavor mixing"],"falsifier":"Compute the two-loop anomalous dimensions for the dimension-seven lepton-number-violating operators and re-evaluate the quoted limits. For the flavor combinations where the paper places the leading term at second-log order, such as $O^{(7)}_{\\bar d u L L D}$ with quark indices $(3,1)$, a two-loop result that is not small compared with the iterated one-loop term would change the bounds; separately, a meson-decay measurement that bounds the same operator more tightly than $0\\nu\\beta\\beta$ would break the claimed ranking.","tokens_in":17130,"feed_emoji":"☢️","tokens_out":15358,"duration_ms":129549,"temperature":0.7,"pith_summary":"The paper asks whether neutrinoless double-$\\beta$ decay ($0\\nu\\beta\\beta$), a lepton-number-violating process in which a nucleus emits two electrons and no neutrinos, can constrain the full set of dimension-seven operators of the Standard Model Effective Field Theory, not just the few that contribute at tree level. The answer it argues for is yes: through the complete one-loop renormalization-group running, operators that look inert at the high scale, especially those involving heavy quarks, mix into operators that drive the decay. With the current KamLAND-Zen limit on the $^{136}$Xe half-life, $T^{0\\nu}_{1/2} > 3.8\\times 10^{26}~\\mathrm{yr}$, this makes $0\\nu\\beta\\beta$ the most restrictive probe of many such operators, reaching new-physics scales around $\\Lambda \\sim 30~\\mathrm{TeV}$ and improving on limits from meson decays. The paper also shows that for some operators and flavor combinations the leading contribution is a second logarithm, so a first-log approximation would miss the effect entirely.","feed_headline":"Mixing turns neutrinoless double-beta decay into a 30 TeV probe","feed_subtitle":"Operators inert at tree level still drive the decay after quantum mixing, beating meson-decay limits.","key_machinery":"The load-bearing object is the one-loop anomalous dimension matrix $\\gamma^{(7)}$ — the matrix controlling how operator coefficients change as the energy scale runs — for the twelve dimension-seven lepton-number-violating SMEFT operators. Its off-diagonal entries, which are proportional to electroweak gauge couplings and to quark Yukawa couplings with CKM factors, convert 'inert' operators into contributors to $0\\nu\\beta\\beta$. The analysis organizes the running through the leading-logarithm expansion of the evolution matrix $U(\\mu_{\\mathrm{ew}}, \\Lambda)$, and it keeps track of cases where a single one-loop mixing step is suppressed, so that the dominant effect comes from a two-step, second-logarithm mixing.","core_discovery":"The paper's central claim is that a complete one-loop renormalization-group analysis changes which lepton-number-violating operators neutrinoless double-$\\beta$ decay can see. Even when a dimension-seven SMEFT operator has no tree-level contribution to $0\\nu\\beta\\beta$, in particular operators built from heavy-quark fields, the one-loop anomalous dimension matrix mixes it into operators that do contribute, through the quark Yukawa couplings and the flavor misalignment encoded in the Cabibbo–Kobayashi–Maskawa (CKM) matrix. These RG-induced contributions make $0\\nu\\beta\\beta$ the strongest existing probe of numerous dimension-seven operators, with constraints reaching $\\Lambda \\sim 30~\\mathrm{TeV}$ for order-one Wilson coefficients, and for several flavor combinations the leading effect arises only at second logarithmic order rather than first.","pith_inferences":["The same Yukawa-induced mixing should also feed other lepton-number-violating searches, such as same-sign dilepton events at colliders and rare kaon decays, so the operators singled out here should appear with correlated coefficients in those channels.","The second-logarithm cases give a concrete target for future two-loop anomalous-dimension calculations: those operators are where the one-loop resummation is least protected and where the quoted limits are most likely to shift.","The paper's comparison with an independent automated implementation shows that some operators also feed neutrino masses; a realistic ultraviolet model must cancel those contributions to keep light-neutrino masses below about 0.1 eV, which can weaken constraints that naively look strongest.","Measuring $0\\nu\\beta\\beta$ in several isotopes could help identify which operator dominates, because RG-induced operators enter through different chiral and nuclear suppressions and would produce distinct isotope ratios relative to standard light-neutrino exchange."],"forward_implications":["Neutrinoless double-beta decay becomes the most stringent published constraint on a broad class of dimension-seven lepton-number-violating operators, surpassing limits from kaon, $D$- and $B$-meson decays for most flavor combinations.","Operators with third-generation quark flavors, which have no tree-level contribution to the decay, can be probed up to effective scales of about $\\Lambda \\sim 30~\\mathrm{TeV}$ through RG-induced mixing.","For specific operators and flavor combinations, such as $O^{(7)}_{\\bar d u L L D}$ with quark indices $(2,1)$ and $(3,1)$, the leading contribution is a second logarithm, so analyses that stop at the first leading logarithm would miss the dominant effect.","The operator $O^{(7)}_{\\bar d L u e H}$ with third-family indices receives an RG-induced constraint comparable to the tree-level valence-quark constraint, because the loop suppression is compensated by a milder chiral suppression.","Next-generation $0\\nu\\beta\\beta$ experiments are expected to nearly double the reach in $\\Lambda$ for these operators."],"supporting_citations":[{"why":"Supplies the complete one-loop anomalous dimension matrix for dimension-seven SMEFT operators that generates all the mixing effects.","marker":"[23]"},{"why":"Defines the operator basis in which the twelve dimension-seven lepton-number-violating operators are classified.","marker":"[19]"},{"why":"Provides the Master Formula that maps SMEFT Wilson coefficients onto 0νββ half-lives.","marker":"[12]"},{"why":"Supplies the companion Master Formula inputs for chiral and nuclear matrix elements used in the numerical rates.","marker":"[13]"},{"why":"Sets the KamLAND-Zen 136Xe half-life limit used as the experimental constraint.","marker":"[4]"},{"why":"Automates the conversion of SMEFT and LEFT coefficients into 0νββ decay rates for the numerical analysis.","marker":"[43]"},{"why":"Provides an independent automated implementation of LNV running used as a comparison target.","marker":"[58]"}],"fun_headline_variants":["Operator mixing gives neutrinoless double-beta decay a 30 TeV reach","RG mixing makes dimension-seven operators visible to 0νββ up to 30 TeV","Dimension-seven operators: 0νββ beats meson decays at 30 TeV","Second-order logs drive some 0νββ decays after RG mixing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical limits inherit the nuclear-physics calculations that turn each operator's strength into a predicted double-beta decay rate; if those calculations are off by a factor of a few, the claim that double-beta decay gives the strongest limits on some operators could fail.","fun_headline_variants_meta":{"raw":{"variants":["Operator mixing gives neutrinoless double-beta decay a 30 TeV reach","RG mixing makes dimension-seven operators visible to 0νββ up to 30 TeV","Dimension-seven operators: 0νββ beats meson decays at 30 TeV","Second-order logs drive some 0νββ decays after RG mixing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001596,"raw_usage":{"total_tokens":6292,"prompt_tokens":809,"completion_tokens":5483,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":425,"completion_tokens_details":{"reasoning_tokens":5393}},"tokens_in":425,"tokens_out":5483,"duration_ms":34621,"temperature":1.0,"reasoning_tokens":5393,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:26:28.076408+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the two-loop anomalous dimensions for the dimension-seven lepton-number-violating operators and re-evaluate the quoted limits. For the flavor combinations where the paper places the leading term at second-log order, such as $O^{(7)}_{\\bar d u L L D}$ with quark indices $(3,1)$, a two-loop result that is not small compared with the iterated one-loop term would change the bounds; separately, a meson-decay measurement that bounds the same operator more tightly than $0\\nu\\beta\\beta$ would break the claimed ranking.","supporting_citations":[{"cited_title":"RGE solver for the complete dim-7 SMEFT interactions and its application to $0\\nu\\beta\\beta$ decay","cited_arxiv_id":"2505.06499","evidence_quote":"Provides an independent automated implementation of LNV running used as a comparison target."}],"review_version":1}