arxiv: 2604.24613 · v1 · submitted 2026-04-27 · ✦ hep-ph

Recognition: unknown

Nucleon decays into one lepton plus two non-strange mesons

Wei-Qi Fan , Yi Liao , Xiao-Dong Ma

Authors on Pith no claims yet

Pith reviewed 2026-05-08 02:44 UTC · model grok-4.3

classification ✦ hep-ph

keywords nucleon decaybaryon number violationeffective field theorythree-body decaysproton decayneutron decaylepton final states

0 comments

The pith

Correlations between two- and three-body nucleon decays produce much tighter lifetime bounds.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that two-body and three-body nucleon decays can be related through a global analysis in low-energy effective field theory without assuming any single operator dominates. This correlation yields improved lower limits on fifteen three-body processes that produce one lepton plus two non-strange mesons. A sympathetic reader would care because these modes test baryon number violation and the new bounds use well-measured two-body rates to constrain processes that are harder to observe directly. The results raise limits on charged-lepton modes by more than three orders of magnitude and set the first bounds above 10^34 years for several neutrino modes while also sharpening three two-body constraints by a factor of two.

Core claim

By correlating two- and three-body nucleon decay processes model-independently in low-energy effective field theory and performing a global analysis, the authors derive significantly improved lower limits on partial lifetimes for fifteen three-body modes with a lepton and two non-strange mesons. For charged-lepton modes these limits surpass current Particle Data Group values by more than three orders of magnitude, while for five neutrino modes they establish partial lifetimes greater than or equal to 10^34 years for the first time. The same analysis also improves constraints on three two-body processes by about a factor of two.

What carries the argument

Model-independent correlations of decay rates through the same low-energy effective operators that govern both two-body and three-body nucleon decays.

If this is right

Lower limits on partial lifetimes for charged-lepton three-body modes exceed current values by more than three orders of magnitude.
First lower limits of order 10^34 years are placed on five neutrino three-body modes.
Constraints on the two-body decays n to e+ pi-, n to mu+ pi-, and p to neutrino pi+ improve by roughly a factor of two.
Well-measured two-body modes can be used to set stronger indirect bounds on three-body modes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Future proton-decay searches could prioritize two-body channels to indirectly tighten limits on three-body channels.
The same correlation technique could be applied to decays involving strange mesons or other baryon-number-violating processes.
A mismatch between predicted and observed rates would point to additional operators or higher-scale physics not captured in the low-energy fit.

Load-bearing premise

The effective operators that mediate two-body nucleon decays also fully control the three-body modes without significant extra contributions from higher-scale physics or operator mixing that would break the rate correlations.

What would settle it

A measured three-body decay rate that exceeds the upper limit predicted by the global fit from two-body data.

Figures

Figures reproduced from arXiv: 2604.24613 by Wei-Qi Fan, Xiao-Dong Ma, Yi Liao.

**Figure 1.** Figure 1: FIG. 1: Correlation diagram of different processes induced by view at source ↗

**Figure 2.** Figure 2: FIG. 2: Leading-order Feynman diagrams contributing to view at source ↗

**Figure 4.** Figure 4: FIG. 4: The improved bounds on 18 two- and three-body view at source ↗

read the original abstract

Nucleon decays into a lepton and two pseudoscalar mesons represent key channels for probing baryon number violation, complementing conventional two-body modes. In this Letter, we model-independently correlate two- and three-body processes within the framework of low-energy effective field theory, performing a global analysis that avoids single-operator-dominance assumption. We derive significantly improved bounds on 15 three-body modes with a lepton ($e^+,\,\mu^+,\hat\nu=\nu/\bar\nu$) and two non-strange mesons ($\pi,\eta$). For charged-lepton modes, our lower limits on partial lifetimes ($\Gamma^{-1}$) surpass current Particle Data Group (PDG) values by more than three orders of magnitude. For 5 (anti)neutrino modes, we establish for the first time $\Gamma^{-1}\gtrsim 10^{34}\,\rm yr$. Additionally, our analysis improves constraints on two-body processes $n\to e^+\pi^-$, $n\to \mu^+\pi^-$, and $p\to \hat\nu \pi^+$ by approximately a factor of 2 compared to the PDG limits. These results highlight the importance of leveraging correlations among different processes to better probe new physics, enabling more stringent constraints on experimentally challenging processes from well-measured ones.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper correlates two- and three-body nucleon decays through a global low-energy EFT fit on baryon-number-violating operators and reports substantially tighter lifetime limits, but the gains rest on the accuracy of the hadronic matrix elements.

read the letter

The central result is a set of improved lower limits on nucleon partial lifetimes for 15 three-body channels with a lepton and two non-strange mesons. By fitting Wilson coefficients of baryon-number-violating operators to two-body data and then applying the same coefficients to three-body processes, they report bounds more than three orders of magnitude stronger than current PDG values for charged-lepton modes and new 10^34 year limits for neutrino modes. They also get modest improvements on a few two-body channels. This is a sensible way to use the correlations that the EFT provides. The global fit without single-operator dominance is a step beyond the usual approach, and it produces concrete numbers that could guide future experiments. The paper appears to have done the necessary work to list the relevant operators and compute or reference the matrix elements. The main question is how robust those matrix elements are. Three-body decays involve more complicated hadronic transitions, and if the uncertainties there are larger than assumed, the propagated bounds will not be as strong as stated. It would help to see explicit tables of the matrix elements used and how errors were handled in the fit. The claim of model independence holds only within the chosen operator basis; any missing higher-scale effects could change the picture. This kind of analysis is useful for the BNV phenomenology community. Experimental groups looking for these rare decays will find the updated limits practical. I would send it to peer review because the framework is standard and the results are falsifiable with better data, even though the hadronic inputs deserve close checking.

Referee Report

3 major / 2 minor

Summary. The manuscript performs a global fit in low-energy effective field theory that correlates two-body nucleon decays (e.g., n→e⁺π⁻) with three-body decays into one lepton plus two non-strange mesons (π, η) via shared baryon-number-violating Wilson coefficients. It reports lower limits on 15 three-body partial lifetimes that exceed current PDG values by more than three orders of magnitude for charged-lepton modes, first-ever bounds Γ⁻¹ ≳ 10³⁴ yr for five neutrino modes, and roughly factor-of-two improvements on the two-body channels n→e⁺π⁻, n→μ⁺π⁻ and p→ν̂π⁺.

Significance. If the operator correlations, hadronic matrix-element ratios, and error propagation are robust, the results would substantially strengthen experimental constraints on baryon-number violation, especially for channels lacking direct searches, and illustrate the utility of EFT-based global analyses for transferring limits between kinematically related processes.

major comments (3)

[§2] §2 (EFT operator basis): the assertion that a single set of dimension-6 operators fully governs both two- and three-body rates requires an explicit listing of the complete basis and a demonstration that no additional three-body operators appear at the same chiral order; without this, the linear relations used to propagate bounds are not guaranteed.
[§4] §4 (global fit and matrix elements): the manuscript must specify the experimental inputs, the treatment of hadronic matrix-element uncertainties (including their ratios between two- and three-body channels), and the full covariance matrix of the fit; the claimed orders-of-magnitude improvements cannot be assessed without these details and a sensitivity study to input variations.
[§5] §5 (neutrino-mode results): the derivation of Γ⁻¹ ≳ 10³⁴ yr for the five neutrino three-body modes relies on the absence of direct limits and on the two-body data constraining the same Wilson coefficients; the paper should quantify the impact of possible higher-dimensional operators or mixing that could break this correlation, as this is the load-bearing assumption for the new bounds.

minor comments (2)

[Abstract] Notation for the neutrino symbol (ν̂) should be defined once at first use and used consistently.
[Figure captions] Figure captions should explicitly state which experimental limits are used as inputs versus which are derived outputs.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The points raised will help improve the clarity and robustness of the presentation. We respond to each major comment below and indicate the planned revisions.

read point-by-point responses

Referee: [§2] §2 (EFT operator basis): the assertion that a single set of dimension-6 operators fully governs both two- and three-body rates requires an explicit listing of the complete basis and a demonstration that no additional three-body operators appear at the same chiral order; without this, the linear relations used to propagate bounds are not guaranteed.

Authors: We agree that an explicit listing strengthens the manuscript. In the revision we will add a table enumerating the complete set of dimension-6 baryon-number-violating operators in the low-energy EFT relevant to non-strange channels, following the standard basis in the literature. We will also include a short chiral-power-counting argument showing that no independent three-body operators enter at the leading order used for our matrix-element relations, thereby justifying the linear correlations between two- and three-body amplitudes. revision: yes
Referee: [§4] §4 (global fit and matrix elements): the manuscript must specify the experimental inputs, the treatment of hadronic matrix-element uncertainties (including their ratios between two- and three-body channels), and the full covariance matrix of the fit; the claimed orders-of-magnitude improvements cannot be assessed without these details and a sensitivity study to input variations.

Authors: We acknowledge the need for greater transparency. The revised manuscript will explicitly list all experimental inputs (PDG two-body limits and any additional constraints), describe the estimation and propagation of hadronic matrix-element uncertainties (including how ratios between two- and three-body channels are obtained from lattice QCD and phenomenological models), provide the full covariance matrix of the Wilson-coefficient fit, and add a sensitivity study varying the dominant inputs within their quoted uncertainties. These additions will allow quantitative assessment of the reported improvements. revision: yes
Referee: [§5] §5 (neutrino-mode results): the derivation of Γ⁻¹ ≳ 10³⁴ yr for the five neutrino three-body modes relies on the absence of direct limits and on the two-body data constraining the same Wilson coefficients; the paper should quantify the impact of possible higher-dimensional operators or mixing that could break this correlation, as this is the load-bearing assumption for the new bounds.

Authors: The neutrino-mode bounds rest on the dimension-6 correlations. In the revision we will expand the discussion to estimate the suppression of dimension-7 and higher operators by additional powers of m_N/Λ_NP and argue that such effects remain subdominant for the lifetime scales considered. We will make the assumption explicit and note its implications. A full global fit including higher-dimensional mixing lies beyond the present scope but could be addressed in future work; we therefore mark this revision as partial. revision: partial

Circularity Check

0 steps flagged

No significant circularity; standard EFT correlation of decay channels

full rationale

The paper performs a global fit of Wilson coefficients in low-energy EFT using experimental limits on two-body nucleon decays, then applies the same coefficients and computed hadronic matrix elements to derive bounds on three-body modes. This is a conventional model-independent analysis that transfers constraints via theory rather than redefining or fitting the output to itself. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations appear in the abstract or described derivation. The claimed improvements follow from the EFT operator basis and matrix-element ratios, which are independent of the target three-body limits. The analysis remains self-contained against external experimental inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the low-energy EFT description of baryon-number-violating operators and on the assumption that two- and three-body rates are linked solely through those operators.

free parameters (1)

Wilson coefficients of baryon-number-violating operators
These coefficients are bounded or fitted in the global analysis and directly determine the derived lifetime limits.

axioms (1)

domain assumption Low-energy effective field theory captures all relevant baryon-number-violating physics at nucleon decay scales
Invoked throughout the abstract as the framework for correlating decay modes.

pith-pipeline@v0.9.0 · 5535 in / 1446 out tokens · 51085 ms · 2026-05-08T02:44:47.876677+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

13 extracted references · 7 canonical work pages · 1 internal anchor

[1]

Nucleon decays into one lepton plus two non-strange mesons

BNV interactions in low-energy EFT (LEFT) [9, 10] as our starting point. Using chiral perturbation theory, we calculate the decay widths for all 24 modes as functions of the LEFT Wilson coefficients (WCs). In contrast to the traditional analyses that assume a single-operator or single-parameter dominance, i.e., turn- ing on one operator at a time, we perf...

work page internal anchor Pith review Pith/arXiv arXiv 2026
[2]

Three-body Decays of the Proton,

M. B. Wise, R. Blankenbecler, and L. F. Abbott, “Three-body Decays of the Proton,”Phys. Rev. D23 (1981) 1591

1981
[3]

MESON EXCHANGE CURRENTS IN THE P DECAY IN NUCLEI,

E. Oset, “MESON EXCHANGE CURRENTS IN THE P DECAY IN NUCLEI,”Nucl. Phys. B304(1988) 820–831

1988
[4]

Search for nucleon decay using the IMB-3 detector,

C. McGrewet al., “Search for nucleon decay using the IMB-3 detector,”Phys. Rev. D59(1999) 052004

1999
[5]

Low-Energy Effective Field Theory below the Electroweak Scale: Operators and Matching,

E. E. Jenkins, A. V. Manohar, and P. Stoffer, “Low-Energy Effective Field Theory below the Electroweak Scale: Operators and Matching,”JHEP03 (2018) 016, [arXiv:1709.04486[hep-ph]]. [Erratum: JHEP 12, 043 (2023)]

work page arXiv 2018
[6]

Extending low energy effective field theory with a complete set of dimension-7 operators,

Y. Liao, X.-D. Ma, and Q.-Y. Wang, “Extending low energy effective field theory with a complete set of dimension-7 operators,”JHEP08(2020) 162, [arXiv:2005.08013[hep-ph]]

work page arXiv 2020
[7]

Comprehensive investigation of nucleon decays into one lepton plus two mesons

W.-Q. Fan, Y. Liao, and X.-D. Ma, “Comprehensive investigation of nucleon decays into one lepton plus two mesons.”. (in preparation)
[8]

Chiral Lagrangian for Deep Mine Physics,

M. Claudson, M. B. Wise, and L. J. Hall, “Chiral Lagrangian for Deep Mine Physics,”Nucl. Phys. B195 (1982) 297–307

1982
[9]

New Chiral Structures for Baryon Number Violating Nucleon Decays,

Y. Liao, X.-D. Ma, and H.-L. Wang, “New Chiral Structures for Baryon Number Violating Nucleon Decays,”Phys. Rev. Lett.135no. 16, (2025) 161801, [arXiv:2504.14855[hep-ph]]

work page arXiv 2025
[10]

Baryon number violating hydrogen decay,

W.-Q. Fan, Y. Liao, X.-D. Ma, and H.-L. Wang, “Baryon number violating hydrogen decay,”Phys. Lett. B862(2025) 139335, [arXiv:2412.20774[hep-ph]]

work page arXiv 2025
[11]

Chiral Lagrangians for Massive Spin 1 Fields,

G. Ecker, J. Gasser, H. Leutwyler, A. Pich, and E. de Rafael, “Chiral Lagrangians for Massive Spin 1 Fields,”Phys. Lett. B223(1989) 425–432

1989
[12]

Interaction of the vector-meson octet with the baryon octet in effective field theory,

Y. Unal, A. Kucukarslan, and S. Scherer, “Interaction of the vector-meson octet with the baryon octet in effective field theory,”Phys. Rev. C92no. 5, (2015) 055208, [arXiv:1510.02205[nucl-th]]

work page arXiv 2015
[13]

Chiral perturbation theory for baryon-number-violating nucleon decay into a vector meson,

Y. Liao, X.-D. Ma, and H.-L. Wang, “Chiral perturbation theory for baryon-number-violating nucleon decay into a vector meson,”Phys. Rev. D112no. 3, (2025) L031704, [arXiv:2506.05052[hep-ph]]. [18]Super-KamiokandeCollaboration, M. Tanakaet al., “Search for proton decay into three charged leptons in 6 0.37 megaton-years exposure of the Super-Kamiokande,” Ph...

work page arXiv 2025