REVIEW 4 major objections 5 minor 39 references
Search for neutron decay into an antineutrino and a neutral kaon in 0.401 megaton-years exposure of Super-Kamiokande
T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A search with 0.401 megaton-years of Super-Kamiokande data finds no evidence for the baryon-number-violating decay of a bound neutron into an antineutrino and a neutral kaon, and sets a 90% C.L.
desk verdict Solid SK search with a factor-of-six improvement, but the limit's absolute scale leans on an unvalidated kaon-interaction model. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the neutral kaon: in the two-body decay $n\to\bar\nu+K^0$ it carries roughly 300 MeV/c of momentum, and its $K^0_S$ component decays promptly in the water into $\pi^+\pi^-$ or $2\pi^0$, producing Cherenkov rings whose total invariant mass and momentum can be reconstructed. The new element is a binned spectrum fit to the $K^0$ invariant mass distributions, rather than a fixed two-dimensional cut, using a Poisson $\chi^2$ with 84 pull terms for systematic uncertainties; this handles the large atmospheric-neutrino background that would otherwise dominate the signal region. The simulation chain combines Fermi motion and binding energy of the bound neutron, nuclear and water kaon-nucleon scattering, and coherent $K^0_L\to K^0_S$ regeneration, together with an improved charged-pion momentum reconstruction that sharpens the reconstructed kaon mass peak.
What would settle it
A beam measurement of neutral-kaon scattering and regeneration on oxygen at 300–500 MeV/c would settle the central systematic assumption; the measured cross sections would either validate the carbon-derived model used here or force a rescaling of the detection efficiency and hence the lifetime limit by the ratio of efficiencies, and any candidate event with an invariant mass peak at the kaon mass inside the fiducial volume would directly falsify the no-signal conclusion.
Extended reading notes
Core claim
The paper's central result is a lower limit on the partial lifetime of a bound neutron decaying via $n\to\bar\nu+K^0$. The analysis selects $K^0_S\to\pi^+\pi^-$ and $K^0_S\to2\pi^0$ candidates in five detector periods, reconstructs the $K^0$ invariant mass and momentum, and fits the mass spectra with a background-plus-signal model; the best fit favors 1.7 signal events, consistent with the zero-signal expectation, with no significant excess. The 90% C.L. upper limit of 27.7 signal events, divided by the simulated detection efficiencies and the $2.68\times10^{32}$ neutrons per kiloton of water, yields $\tau/B > 7.8\times10^{32}$ years. This is six times stronger than the previous published limit and stands as the most stringent constraint on $n\to\bar\nu K^0$ decay.
Load-bearing premise
The limit is proportional to the simulated signal detection efficiency, and that efficiency rests on a Monte Carlo kaon-nucleon interaction model whose neutral-kaon regeneration probability was measured with a carbon target and then applied to oxygen and water inside Super-Kamiokande; if the real kaon interaction rates in water differ from the model, the efficiency and the lifetime limit shift by the same factor.
Editorial extensions
If this is right
- If the limit holds, any supersymmetric GUT in which $n\to\bar\nu K^0$ is a dominant neutron decay channel must have a partial lifetime above $7.8\times10^{32}$ years.
- The spectrum-fit method demonstrated here can be carried over to other nucleon decay searches with irreducible backgrounds, such as $n\to\bar\nu\pi^0$ and $p\to\bar\nu\pi^+$.
- The updated kaon interaction model and pion momentum reconstruction feed into future Super-Kamiokande kaon searches, including the related $p\to\bar\nu K^+$ channel.
- Combining this limit with the proton lifetime limit on $p\to\bar\nu K^+$ tightens the allowed parameter space of SUSY GUT models that predict both modes.
Reading between the lines
- A separate and formally identical bound could be quoted for a $B+L$-conserving mode $n\to\nu K^0$, since the paper notes the search is equally sensitive to a neutrino in the final state; the paper does not derive that separate limit.
- With an order of magnitude more exposure, the same method would reach partial lifetimes around $10^{34}$ years, the range in which some supersymmetric GUT models place the dominant neutron decay; future detectors could test this channel at that scale.
- The largest unvalidated ingredient is the neutral-kaon interaction model; a dedicated measurement of $K^0$ cross sections on oxygen around 300–500 MeV/c would directly reduce the dominant systematic and make the lifetime limit less model-dependent.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper reports a search for bound-neutron decay n→νbar+K0 using 0.401 Mton-years of Super-Kamiokande pure-water data (SK-I through SK-V), 4.4 times the exposure of the previous SK-I search. The analysis selects K0S→π+π− and K0S→2π0 candidates and performs a binned Poisson likelihood fit with 84 pull-term systematic parameters to the reconstructed K0 invariant mass distributions, rather than applying the narrower mass window used previously. No significant excess is found: the best-fit signal yield is 1.7 events, and the 90% C.L. upper limit is 27.7 signal events. Using per-phase efficiencies in Eq. (3), this translates to τ/B > 7.8×10^32 years at 90% C.L., a factor of six improvement over the previous limit of 1.3×10^32 years.
Significance. If correct, this is the world's most stringent constraint on n→νbar+K0 and is the first Super-K search for this mode to use a spectrum fit over the full pure-water dataset. The statistical treatment is standard for this field: a Poisson chi-square with pull terms for correlated systematics, and the limit follows from the fitted signal normalization rather than from a circular definition. The collaboration gives a legible account of the event selections, background composition, and systematic pulls, and the analysis is a natural extension of its nucleon-decay program. The principal caveat, discussed below, is that the absolute lifetime scale is set by the MC signal efficiency, whose kaon-nucleon interaction and regeneration model is not validated inside the detector; however, the quoted 25–50% systematic uncertainties at least acknowledge the size of this dependence.
major comments (4)
- [Section III A and Eq. (3)] The central limit inherits its absolute scale from the signal efficiency η_i, which depends on the kaon-nucleon interaction model and on the coherent K0L→K0S regeneration probability derived from carbon-target data [24] and applied to oxygen and water. The fit cannot constrain this model: the best-fit pulls for kaon FSI and SI are only 0.030 and 0.031 (Table II), and neutral-kaon production in the atmospheric-neutrino background is below 1% (Section III B). Since Eq. (3) divides by η_i, any unmodeled bias in the kaon interaction or regeneration rescales the headline limit linearly. The paper should either validate the model with an in-situ control sample, quantify how the limit changes under alternative kaon-interaction prescriptions, or at minimum state explicitly that the 7.8×10^32 year value scales with these model assumptions.
- [Section VI B and Eq. (3)] The numerical result 7.8×10^32 years is obtained by summing N^i_90CL/η_i over SK-I through SK-V, but the paper does not provide a table of the per-phase efficiencies η_i or the per-phase upper limits N^i_90CL. Without these inputs the central number cannot be reproduced, and the reader cannot see which phase dominates the sum. Please include a table with these values, and state whether η_i are the nominal MC efficiencies and how the 25% and 50% kaon uncertainties enter the quoted N^i_90CL.
- [Section V C and Section VI B] The text states that the roughly 30% deficit of the atmospheric-neutrino MC relative to data after selection B-6 is 'resolved' by systematic uncertainties, yet the best-fit pulls in Tables II and III are all well below 1σ. It is not demonstrated which combination of pull terms absorbs the excess, nor how stable the 27.7-event 90% C.L. upper limit is if the background normalization is treated as a free parameter or constrained differently. This matters because the limit is extracted from the same fit. Please provide a robustness check, such as the fitted β and N90 under a free background normalization, and identify the pulls that contribute most to absorbing the excess.
- [Section VI B] The derivation of the 90% C.L. upper limit from the fit is not described: the paper reports a best-fit β and an upper limit of 27.7 events, but does not state whether the limit is obtained from Δχ2=2.71 with the systematics profiled, from a different χ2-based construction, or from a Feldman-Cousins calculation. Since the central claim is a limit, the statistical procedure used to set it should be specified explicitly.
minor comments (5)
- [Section V A and V B] The statement that the K0 invariant mass cut is not applied in this analysis is confusing because A-6 and B-6 include the broad cut 300 < W_K0 < 700 MeV/c^2; please clarify that the narrow previous-analysis cuts (400–600 and 450–550 MeV/c^2) are replaced by the broad 300–700 MeV/c^2 fit range.
- [Section V] The sentence beginning 'Since the invariant mass distribution has a narrower peak...' is repeated in consecutive sentences; one of the two occurrences should be removed.
- [Section III A] The text says '98% of K0 from nucleon decay exit from oxygen nucleus as K0' and then '98% of K0S that exit the nucleus promptly decay...'; it would be clearer to state explicitly that the 98% refers to K0 survival without hadronic interaction, and to define whether K0L regeneration is included in the 0.1% of MC events mentioned in the same paragraph.
- [Equations (1) and (2)] The normalization convention of N_sig_i and the units of β should be stated explicitly; from the text β is a signal yield in events rather than a dimensionless multiplier of a normalized template, so the convention should be defined.
- [Figure 4] The caption notes that the error bars are too small to be visible; please provide the numerical signal efficiencies after the final selection in a table or in the caption so that the reader can connect them to Eq. (3).
Circularity Check
No significant circularity: the lifetime limit is a direct experimental upper bound obtained from a fit to data divided by a simulated efficiency, with no step reducing to its own inputs.
full rationale
The paper derives its headline limit using Eq. (3), which converts the fitted 90% C.L. upper limit on signal events into a lifetime by dividing by the simulated signal efficiency. The efficiency is an MC input, not a fitted parameter; the fit determines only the signal normalization beta. The kaon-interaction and regeneration models affect the efficiency, but the paper reports these as explicit systematic uncertainties with fit pulls. This is model dependence, not circularity: no step defines the target quantity in terms of itself, and the central claim is a measured limit against the simulated expectation.
Assumptions & free parameters
assumptions (6)
- standard math Poisson likelihood chi-square with pull terms is used for the spectrum fit.
- domain assumption Fermi momentum distribution of bound neutrons in oxygen is taken from electron scattering off 12C (Nakamura et al., Ref [19]).
- domain assumption The binding energy is modeled as a Gaussian with mean 39.0 (15.5) MeV and sigma 10.2 (3.82) MeV for the s-state (p-state), applied equally to p3/2 and p1/2 states.
- domain assumption Neutron population ratio in oxygen is 1:3 for s-state to p-state, per nuclear shell model, and 10% of decays are correlated with a neighboring nucleon.
- domain assumption Neutral kaon interactions in the nucleus and water, including elastic scattering, charge exchange, and coherent regeneration, are modeled using carbon-target regeneration data and the updated kaon interaction model from the p -> mu + K0 search.
- domain assumption Atmospheric neutrino flux and neutrino-nucleus interactions are simulated with the Honda flux and NEUT 5.4.0.1, with three-flavor oscillation parameters from PDG 2022.
Cite this review
Pith. "Pith review of Search for neutron decay into an antineutrino and a neutral kaon in 0.401 megaton-years exposure of Super-Kamiokande." pith.science (2026). https://pith.science/paper/4E545TVW
@misc{pith2026250614406,
author = {Pith},
title = {Pith review of: Search for neutron decay into an antineutrino and a neutral kaon in 0.401 megaton-years exposure of Super-Kamiokande},
year = {2026},
howpublished = {\url{https://pith.science/paper/4E545TVW}},
note = {Machine review of arXiv:2506.14406}
}
abstract
We searched for bound neutron decay via $n\to\bar{\nu}+K^0$ predicted by the Grand Unified Theories in 0.401 Mton$\cdot$years exposure of all pure water phases in the Super-Kamiokande detector. About 4.4 times more data than in the previous search have been analyzed by a new method including a spectrum fit to kaon invariant mass distributions. No significant data excess has been observed in the signal regions. As a result of this analysis, we set a lower limit of $7.8\times10^{32}$ years on the neutron lifetime at a 90% confidence level.
Figures
Reference graph
Works this paper leans on
-
[24]
P. H. Eberhard and F. Uchiyama, Neutral kaon regen- eration probabilities at asymmetric phi factory energies, Nucl. Instrum. Meth. A350, 144 (1994)
work page 1994
-
[1]
H. Georgi and S. L. Glashow, Unity of All Elementary- Particle Forces, Phys. Rev. Lett.32, 438 (1974)
work page 1974
-
[2]
Fritzsch and P
H. Fritzsch and P. Minkowski, Unified interactions of lep- tons and hadrons, Annals of Physics93, 193 (1975)
1975
- [3]
-
[4]
H. Goh, R. Mohapatra, S. Nasri, and S.-P. Ng, Proton decay in a minimal SUSY SO(10) model for neutrino mix- ings, Physics Letters B587, 105 (2004)
work page 2004
-
[5]
S. Antusch, C. Hohl, and V. Susiˇ c, Employing nucleon decay as a fingerprint of SUSY GUT models using SusyTCProton, J. High Energy Phys.06, 022 (2021)
work page 2021
-
[6]
K. Abeet al.(Super-Kamiokande Collaboration), Search for proton decay viap→νK + using 260 kiloton·year data of Super-Kamiokande, Phys. Rev. D90, 072005 (2014)
work page 2014
-
[7]
Vissani, (B+L)-conserving nucleon decays in super- symmetric models, Phys
F. Vissani, (B+L)-conserving nucleon decays in super- symmetric models, Phys. Rev. D52, 4245 (1995)
work page 1995
Show all 39 references
-
[8]
K. S. Babu and R. N. Mohapatra,B−Lviolating nucleon decay and GUT scale baryogenesis inSO(10), Phys. Rev. D86, 035018 (2012)
2012
-
[9]
Kobayashiet al.(Super-Kamiokande Collaboration), Search for nucleon decay via modes favored by supersym- metric grand unification models in Super-Kamiokande-I, Phys
K. Kobayashiet al.(Super-Kamiokande Collaboration), Search for nucleon decay via modes favored by supersym- metric grand unification models in Super-Kamiokande-I, Phys. Rev. D72, 052007 (2005)
2005
-
[10]
Fukudaet al.(Super-Kamiokande Collaboration), The Super-Kamiokande detector, Nucl
Y. Fukudaet al.(Super-Kamiokande Collaboration), The Super-Kamiokande detector, Nucl. Instrum. Meth. A 501, 418 (2003)
2003
-
[11]
Abeet al., Calibration of the Super-Kamiokande De- tector, Nucl
K. Abeet al., Calibration of the Super-Kamiokande De- tector, Nucl. Instrum. Meth. A737, 253 (2014)
2014
-
[12]
Mine, The Super-Kamiokande and Other Detectors: A Case Study of Large Volume Cherenkov Neutrino Detec- tors, inInstrumentation and Techniques in High Energy Physics, edited by D
S. Mine, The Super-Kamiokande and Other Detectors: A Case Study of Large Volume Cherenkov Neutrino Detec- tors, inInstrumentation and Techniques in High Energy Physics, edited by D. Lincoln (2024) pp. 251–290
2024
-
[13]
R. Brun, F. Bruyant, F. Carminati, S. Giani, M. Maire, A. McPherson, G. Patrick, and L. Ur- ban, GEANT Detector Description and Simulation Tool 10.17181/CERN.MUHF.DMJ1 (1994)
1994 doi
-
[14]
Hayato and L
Y. Hayato and L. Pickering, The NEUT neutrino inter- action simulation program library, Eur. Phys. J. ST230, 4469 (2021)
2021
-
[15]
Zeitnitz and T
C. Zeitnitz and T. A. Gabriel, The GEANT - CALOR interface and benchmark calculations of ZEUS test calorimeters, Nucl. Instrum. Meth. A349, 106 (1994)
1994
-
[16]
HETC: Monte Carlo High-Energy Nucleon-Meson Trans- port Code, Oak Ridge National Laboratory RSIC Com- puter Code Collection CCC-178
-
[17]
Fass´ o, G
A. Fass´ o, G. Stevenson, J. Zazula, A. Ferrari, P. Sala, and J. Ranft, Fluka92, inProc. 1st Workshop on Simulating Accelerator Radiation Environments(1993)
1993
-
[18]
J. O. Johnson and T. A. Gabriel,A user’s guide to MI- CAP: A Monte Carlo Ionization Chamber Analysis Pack- age, Tech. Rep. (1988)
1988
-
[19]
Nakamuraet al., The Reaction C-12 (e, e’ p) at 700- MeV and DWIA Analysis, Nucl
K. Nakamuraet al., The Reaction C-12 (e, e’ p) at 700- MeV and DWIA Analysis, Nucl. Phys. A268, 381 (1976)
1976
-
[20]
Hiramatsuet al., inProc
S. Hiramatsuet al., inProc. Int. Conf. on nuclear struc- ture studies using electron scattering and photoreaction (1972) p. 429
1972
-
[21]
Yamazaki and Y
T. Yamazaki and Y. Akaishi, Nuclear medium effects on invariant mass spectra of hadrons decaying in nuclei, Phys. Lett. B453, 1 (2000)
2000
-
[22]
Mayer and J
M. Mayer and J. Jensen,Elementary Theory of Nuclear Shell Structure, Structure of matter series (Wiley, 1955)
1955
-
[23]
Matsumotoet al.(Super-Kamiokande Collaboration), Search for proton decay viap→µ +K 0 in 0.37 megaton- years exposure of Super-Kamiokande, Phys
R. Matsumotoet al.(Super-Kamiokande Collaboration), Search for proton decay viap→µ +K 0 in 0.37 megaton- years exposure of Super-Kamiokande, Phys. Rev. D106, 072003 (2022)
2022
-
[25]
Honda, T
M. Honda, T. Kajita, K. Kasahara, and S. Midorikawa, Improvement of low energy atmospheric neutrino flux cal- culation using the JAM nuclear interaction model, Phys. 11 TABLE III. Systematic uncertainties related to neutrino flux, interaction, and oscillation for atmospheric n...
2011
-
[26]
Rein and L
D. Rein and L. M. Sehgal, Neutrino-excitation of baryon resonances and single pion production, Annals of Physics 133, 79 (1981)
1981
-
[27]
Berger and L
C. Berger and L. M. Sehgal, Lepton mass effects in single pion production by neutrinos, Phys. Rev. D76, 113004 (2007), Erratum: Phys. Rev. D77, 059901 (2008)
2007
-
[28]
→2𝜋" Total momentum [MeV/c]01002003004005006007008009001000 Number of events 0 20 40 60 80 100 120 140 160 180 200 220 Data MC0Kν→n MCνAtm- 𝐾!
assuming normal mass ordering. The updates of the kaon-nucleon interaction in nuclei and water are also included in the atmospheric neutrino MC used in this analysis. However, the rate of neutral kaon production is less than 1% in the remaining background events due to its sma...
-
[29]
R. L. Workmanet al.(Particle Data Group), Review of Particle Physics, PTEP2022, 083C01 (2022)
2022
-
[30]
Ashieet al.(Super-Kamiokande Collaboration), Mea- surement of atmospheric neutrino oscillation parameters by Super-Kamiokande I, Phys
Y. Ashieet al.(Super-Kamiokande Collaboration), Mea- surement of atmospheric neutrino oscillation parameters by Super-Kamiokande I, Phys. Rev. D71, 112005 (2005)
2005
-
[31]
Shiozawa, Reconstruction algorithms in the Super- Kamiokande large water Cherenkov detector, Nucl
M. Shiozawa, Reconstruction algorithms in the Super- Kamiokande large water Cherenkov detector, Nucl. In- strum. Meth. A433, 240 (1999)
1999
-
[32]
Abeet al.(Super-Kamiokande Collaboration), Search for Nucleon Decay vian→¯νπ 0 andp→¯νπ + in Super- Kamiokande, Phys
K. Abeet al.(Super-Kamiokande Collaboration), Search for Nucleon Decay vian→¯νπ 0 andp→¯νπ + in Super- Kamiokande, Phys. Rev. Lett.113, 121802 (2014). 12
2014
-
[33]
V. Takhistovet al.(Super-Kamiokande Collaboration), Search for Nucleon and Dinucleon Decays with an Invis- ible Particle and a Charged Lepton in the Final State at the Super-Kamiokande Experiment, Phys. Rev. Lett. 115, 121803 (2015)
2015
-
[34]
G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, and A. Palazzo, Getting the most from the statistical analy- sis of solar neutrino oscillations, Phys. Rev. D66, 053010 (2002)
2002
-
[35]
L. L. Salcedo, E. Oset, M. J. Vicente-Vacas, and C. Garcia-Recio, Computer Simulation of Inclusive Pion Nuclear Reactions, Nucl. Phys. A484, 557 (1988)
1988
-
[36]
de Perio, NEUT Pion FSI, AIP Conference Proceed- ings1405, 223 (2011)
P. de Perio, NEUT Pion FSI, AIP Conference Proceed- ings1405, 223 (2011)
2011
-
[37]
de Perio,Joint Three-Flavour Oscillation Analysis of νµ Disappearance andν e Appearance in the T2K Neu- trino Beam, Ph.D
P. de Perio,Joint Three-Flavour Oscillation Analysis of νµ Disappearance andν e Appearance in the T2K Neu- trino Beam, Ph.D. thesis, Toronto U. (2014)
2014
-
[38]
Westeret al.(Super-Kamiokande Collaboration), At- mospheric neutrino oscillation analysis with neutron tagging and an expanded fiducial volume in Super- Kamiokande I–V, Phys
T. Westeret al.(Super-Kamiokande Collaboration), At- mospheric neutrino oscillation analysis with neutron tagging and an expanded fiducial volume in Super- Kamiokande I–V, Phys. Rev. D109, 072014 (2024)
2024
-
[39]
Wester,Discerning the Neutrino Mass Ordering using Atmospheric Neutrinos in Super-Kamiokande I-V, Ph.D
T. Wester,Discerning the Neutrino Mass Ordering using Atmospheric Neutrinos in Super-Kamiokande I-V, Ph.D. thesis, Boston University, Boston U. (2023)
2023
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.