Recognition: unknown
Visible Neutrino Decay As An Open Quantum System
Pith reviewed 2026-05-10 16:35 UTC · model grok-4.3
The pith
Neutrino oscillations combined with decays are fully described using open quantum system methods like Lindblad equations and Kraus operators.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Decays of heavier neutrino mass eigenstates into lighter ones, while very slow in the Standard Model, can be significantly enhanced in scenarios with more than three neutrino flavours, or in models with new ultra-light particles such as Majorons. A full theoretical description is challenging due to the intricate interplay between oscillations and decay, interference between different decay channels, and the possibility of multi-step decay cascades. We develop a fully general description of arbitrarily complex systems of oscillating and decaying neutrinos using methods from the theory of open quantum systems. Notably, we demonstrate how such systems can be implemented using the Lindblad maste
What carries the argument
Application of open quantum system tools (Lindblad master equation, Liouvillian superoperator, and Kraus operators) to the density-matrix evolution of multi-flavor neutrino states, capturing both unitary oscillations and non-unitary decay effects in one formalism.
If this is right
- Allows consistent modeling of interference between multiple decay channels and multi-step cascades.
- Provides numerically efficient implementations that avoid solving differential equations at each time step.
- Extends directly to models with more than three neutrino flavors or new ultra-light mediators.
- Preserves information about coherence and off-diagonal density matrix elements during evolution.
Where Pith is reading between the lines
- The same formalism could incorporate matter potentials or other environmental couplings by adding appropriate terms to the Lindblad operators.
- It suggests a route to treat neutrino flavor evolution as a quantum information process, with decay acting like a controlled decoherence channel.
- Benchmarks against known analytic limits (pure oscillations or pure exponential decay) would quickly validate the implementation for more complex cases.
Load-bearing premise
The interplay of oscillations, decays, interference, and cascades in multi-flavor or Majoron-extended neutrino models can be fully and accurately captured by standard open quantum system formalisms without additional unstated approximations or loss of unitarity.
What would settle it
A direct numerical comparison for a two-flavor neutrino system with decay where the probabilities computed from the Lindblad or Kraus evolution differ from those obtained by solving the exact time-dependent non-Hermitian Hamiltonian.
Figures
read the original abstract
Decays of heavier neutrino mass eigenstates into lighter ones, while very slow in the Standard Model, can be significantly enhanced in scenarios with more than three neutrino flavours, or in models with new ultra-light particles such as Majorons. A full theoretical description is challenging due to the intricate interplay between oscillations and decay, interference between different decay channels, and the possibility of multi-step decay cascades. In this paper, we develop a fully general description of arbitrarily complex systems of oscillating and decaying neutrinos using methods from the theory of open quantum systems. Notably, we demonstrate how such systems can be implemented using the Lindblad master equation, the Liouvillian superoperator, as well as Kraus operators. The last two methods eschew the need for solving a differential equation, thereby showing superior numerical performance.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a fully general description of arbitrarily complex systems of oscillating and decaying neutrinos (including visible decays, interference, and multi-step cascades) by mapping them onto open quantum system formalisms. It demonstrates implementations via the Lindblad master equation, the Liouvillian superoperator, and Kraus operators, noting that the latter two avoid solving differential equations and offer superior numerical performance.
Significance. If the explicit constructions correctly encode visible decay channels (retaining daughter neutrinos in the reduced density matrix) while preserving off-diagonal coherences and positivity, the framework would provide a systematic, verifiable tool for BSM neutrino models with enhanced decays. The use of three equivalent formalisms for cross-checks and the focus on numerical efficiency are clear strengths.
major comments (1)
- [Abstract and implementation sections] Abstract and the section describing the dissipator/Kraus maps: the central claim that standard Lindblad/Liouvillian/Kraus constructions fully capture visible multi-step cascades without extra damping or loss of interference requires explicit jump operators (or Kraus operators) that re-inject daughter neutrinos into the tracked system while maintaining the off-diagonal terms responsible for oscillations. The provided description does not specify these operators, leaving open whether the implementation treats visible decays as irreversible loss (inappropriate here) or applies an unstated secular approximation.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for acknowledging the potential utility of the open quantum system framework for BSM neutrino models with enhanced visible decays. We address the major comment below.
read point-by-point responses
-
Referee: [Abstract and implementation sections] Abstract and the section describing the dissipator/Kraus maps: the central claim that standard Lindblad/Liouvillian/Kraus constructions fully capture visible multi-step cascades without extra damping or loss of interference requires explicit jump operators (or Kraus operators) that re-inject daughter neutrinos into the tracked system while maintaining the off-diagonal terms responsible for oscillations. The provided description does not specify these operators, leaving open whether the implementation treats visible decays as irreversible loss (inappropriate here) or applies an unstated secular approximation.
Authors: We thank the referee for this observation. The manuscript constructs the open quantum system such that all neutrino mass eigenstates participating in the decays (including daughters) are retained in the Hilbert space; visible decays are therefore not modeled as irreversible loss. The Lindblad jump operators are explicitly of the form L_{j←i} = √Γ_{i→j} |ν_j⟩⟨ν_i| (and their Hermitian conjugates), which map parent to daughter states inside the same space. When inserted into the dissipator, these operators generate both population transfer and off-diagonal contributions that preserve the relative phases responsible for oscillations. The full (non-secular) Liouvillian is retained, with no approximation invoked. Multi-step cascades are handled by composing the corresponding Kraus operators, which are products of the individual jump maps. The abstract summarizes the general construction; the implementation sections derive the Lindblad, Liouvillian-superoperator, and Kraus representations from these operators. To make the explicit forms and their action on coherences more transparent, we will add a dedicated subsection containing the concrete jump/Kraus operators for a two-step visible cascade together with a short numerical check confirming coherence preservation. revision: yes
Circularity Check
No circularity: standard open-quantum-system formalisms applied to neutrino decay without self-referential reduction
full rationale
The paper's central contribution is a mapping of oscillating and decaying neutrino systems onto independently established open-quantum-system tools (Lindblad master equation, Liouvillian superoperator, Kraus operators). These formalisms pre-exist the present work and are not derived from neutrino-specific data or prior self-citations within the paper. No parameters are fitted to the target observables, no predictions are generated by construction from the same inputs, and the derivation chain consists of explicit implementation steps rather than tautological re-labeling. The abstract and described framework remain self-contained against external benchmarks in quantum mechanics.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Neutrino systems with decay and oscillation can be modeled as open quantum systems using Lindblad, Liouvillian, or Kraus formalisms
Reference graph
Works this paper leans on
-
[1]
Visible Neutrino Decay As An Open Quantum System
INTRODUCTION Among all the unusual properties that neutrinos pos- sess, their decay is perhaps the least studied one. This is not surprising, given that in the Standard Model, the rate for loop-suppressed radiative decays of the form νi →ν j +γis≲10 −42 yr−1 [1–6]. (Here,ν i,ν j denote neutrino mass eigenstates.) νi νj γ W ∓ ℓ± ℓ± However, neutrino decays...
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[2]
regeneration
EXISTING RESUL TS ON NEUTRINO DECA Y 2.1. Decay Rates While neutrino decay occurs even in the Standard Model (ν j →ν k +γ), the corresponding rates are too small to be phenomenologically relevant. This changes in the presence of heavier, sterile, neutrinos or in the presence of additional light particles that open up new decay modes. In the following, we ...
-
[3]
Lindbladian
OPEN QUANTUM SYSTEMS We now introduce the concepts from the theory of open quantum systems which we will use to reformulate the neutrino oscillation+decay problem. 3.1. T oy Example To set the stage and introduce the basic idea, we con- sider as a toy example theamplitude dampingmodel of an excited state of an atom (|1⟩) decaying to the ground state (|0⟩)...
-
[4]
Lindblad Approach In this section, we pose the visible neutrino decay sys- tem as an open quantum system, deriving the main re- sults of the present work
NEUTRINO DECA Y AS AN OPEN QUANTUM SYSTEMS 4.1. Lindblad Approach In this section, we pose the visible neutrino decay sys- tem as an open quantum system, deriving the main re- sults of the present work. Some steps in this direction have been taken in the existing literature by invoking a modified von Neumann equation of the form [58, 59] dρ(E) dt =−i[H(E)...
-
[5]
We have one unstable neutrino, and we allow for the possibility of helicity-violating decays (see relevant formulae in Section 2.1
or KamLAND [66]. We have one unstable neutrino, and we allow for the possibility of helicity-violating decays (see relevant formulae in Section 2.1. violating decays are allowed. This requires us to track six neutrino species (three neutrino species, three anti- neutrino species) with four possible decay modes. The density matrix in flavour space therefor...
-
[6]
(48), that is, to decays from a particular initial energy bin,n, to a particular final energy bin,m, withn≥m
Single Decay For a single decay,ν 3 →ν 1, there areN E(NE + 1)/2 Lindblad operators, each corresponding to one block of Eq. (48), that is, to decays from a particular initial energy bin,n, to a particular final energy bin,m, withn≥m. The (m, n) block of each Lindblad operator will have the structure 0 0■ 0 0 0 0 0 0 .(51) (with all other blocks be...
-
[7]
Dual Decay Consider now multiple decays from a single unstable neutrinoν 3 →ν 2, ν 3 →ν 1. Now, each of theN E(NE + 1)/2 Lindblad operators will have the form 0 0■ 0 0■ 0 0 0 .(56) By similar arguments as in Section 4 4.2 1, the form of the Kraus operators in the absence of oscillations can again be inferred: M(mn) osc 0 (L) =δ nm 1 0 0 0 1 0 ...
-
[8]
0■0 0 0 0 0 0 0 .(62) The first matrix drivesν 3 →ν 2 transitions and the sec- ond one drivesν 2 →ν 1 decays
Cascade Decay Consider now a decay chainν 3 →ν 2 →ν 1, for which theN E(NE + 1)/2 pairs of Lindblad operators take the form 0 0 0 0 0■ 0 0 0 . 0■0 0 0 0 0 0 0 .(62) The first matrix drivesν 3 →ν 2 transitions and the sec- ond one drivesν 2 →ν 1 decays. The blocks of the dy- namical mapE(L) now have the structure ■0 0 0■0 0 0...
-
[9]
DISCUSSION 50 100 200 300 500 700 Number of energy bins, NE 100 101 102 103 104 105 Wall time [s] O(N3 E) O(N2 E) g = 0.1 L = 50 km Nt = 100 Nν = 3 Lindbladian ODE O(N2.71 E ) Dynamical Map/Kraus O(N2.03 E ) FIG. 6. Numerical complexity of different algorithms for com- puting combined neutrino oscillation+decay probabilities as a function of the number of...
-
[10]
CONCLUSIONS In summary, we have investigated the theory of neu- trino oscillation and decay. By applying methods from the theory of open quantum systems, we have found that arbitrarily complicated oscillation+decay problems can be tackled with high numerical efficiency, either by solv- ing the Lindblad master equation or by computing a set of Kraus operat...
-
[11]
Shrock, DecayL 0 →ν lγin gauge theories of weak and electromagnetic interactions, Phys
R. Shrock, DecayL 0 →ν lγin gauge theories of weak and electromagnetic interactions, Phys. Rev. D9, 743 (1974)
1974
-
[12]
B. W. Lee and R. E. Shrock, Natural Suppression of Sym- metry Violation in Gauge Theories: Muon - Lepton and Electron Lepton Number Nonconservation, Phys. Rev. D 16, 1444 (1977)
1977
-
[13]
W. J. Marciano and A. I. Sanda, Exotic Decays of the Muon and Heavy Leptons in Gauge Theories, Phys. Lett. B67, 303 (1977)
1977
-
[14]
P. B. Pal and L. Wolfenstein, Radiative Decays of Mas- sive Neutrinos, Phys. Rev. D25, 766 (1982)
1982
-
[15]
R. E. Shrock, Electromagnetic Properties and Decays of Dirac and Majorana Neutrinos in a General Class of Gauge Theories, Nucl. Phys. B206, 359 (1982)
1982
-
[16]
Xing and S
Z.-z. Xing and S. Zhou,Neutrinos in particle physics, astronomy and cosmology, Advanced Topics in Science and Technology in China (Springer Berlin, Heidelberg, 2011)
2011
-
[17]
Chikashige, R
Y. Chikashige, R. N. Mohapatra, and R. D. Peccei, Are There Real Goldstone Bosons Associated with Broken Lepton Number?, Phys. Lett. B98, 265 (1981)
1981
-
[18]
G. B. Gelmini and M. Roncadelli, Left-Handed Neutrino Mass Scale and Spontaneously Broken Lepton Number, Phys. Lett. B99, 411 (1981)
1981
-
[19]
J. W. F. Valle, Fast Neutrino Decay in Horizontal Ma- joron Models, Phys. Lett. B131, 87 (1983)
1983
- [20]
- [21]
-
[22]
K. Abazajian, G. M. Fuller, and W. H. Tucker, Direct detection of warm dark matter in the X-ray, Astrophys. J.562, 593 (2001), arXiv:astro-ph/0106002
-
[23]
A White Paper on keV Sterile Neutrino Dark Matter
M. Dreweset al., A White Paper on keV Sterile Neutrino Dark Matter, JCAP01, 025, arXiv:1602.04816 [hep-ph]
-
[24]
A. Boyarsky, M. Drewes, T. Lasserre, S. Mertens, and O. Ruchayskiy, Sterile neutrino Dark Matter, Prog. Part. Nucl. Phys.104, 1 (2019), arXiv:1807.07938 [hep-ph]
-
[25]
B. Dasgupta and J. Kopp, Sterile Neutrinos, Phys. Rept. 928, 1 (2021), arXiv:2106.05913 [hep-ph]
- [26]
-
[27]
T. Abrah˜ ao, H. Minakata, H. Nunokawa, and A. A. Quiroga, Constraint on Neutrino Decay with Medium- Baseline Reactor Neutrino Oscillation Experiments, JHEP11, 001, arXiv:1506.02314 [hep-ph]
-
[28]
P. Coloma and O. L. G. Peres, Visible neutrino decay at DUNE, (2017), arXiv:1705.03599 [hep-ph]
- [29]
-
[30]
S. Choubey, D. Dutta, and D. Pramanik, Invisible neu- trino decay in the light of NOvA and T2K data, JHEP 08, 141, arXiv:1805.01848 [hep-ph]
- [31]
-
[32]
M. Maltoni and W. Winter, Testing neutrino oscilla- tions plus decay with neutrino telescopes, JHEP07, 064, arXiv:0803.2050 [hep-ph]
-
[33]
P. Baerwald, M. Bustamante, and W. Winter, Neutrino Decays over Cosmological Distances and the Implications for Neutrino Telescopes, JCAP10, 020, arXiv:1208.4600 [astro-ph.CO]
- [34]
-
[35]
M. Archidiacono and S. Hannestad, Updated constraints on non-standard neutrino interactions from Planck, JCAP07, 046, arXiv:1311.3873 [astro-ph.CO]
-
[36]
M. Escudero and M. Fairbairn, Cosmological Constraints on Invisible Neutrino Decays Revisited, Phys. Rev. D 100, 103531 (2019), arXiv:1907.05425 [hep-ph]
-
[37]
G. Barenboim, J. Z. Chen, S. Hannestad, I. M. Old- engott, T. Tram, and Y. Y. Y. Wong, Invisible neu- trino decay in precision cosmology, JCAP03, 087, arXiv:2011.01502 [astro-ph.CO]
-
[38]
M. Lindner, T. Ohlsson, and W. Winter, A Combined treatment of neutrino decay and neutrino oscillations, Nucl. Phys. B607, 326 (2001), arXiv:hep-ph/0103170
-
[39]
A. Abdullahi and P. B. Denton, Visible Decay of Astro- physical Neutrinos at IceCube, Phys. Rev. D102, 023018 (2020), arXiv:2005.07200 [hep-ph]
-
[40]
K. Stankevich, A. Studenikin, and M. Vyalkov, Gener- alized Lindblad master equation for neutrino evolution, Phys. Rev. D111, 036014 (2025), arXiv:2411.19303 [hep- ph]
-
[41]
K. Stankevich and A. Studenikin, Neutrino quantum de- coherence engendered by neutrino radiative decay, Phys. Rev. D101, 056004 (2020), arXiv:2002.02621 [hep-ph]
-
[42]
J. Bernabeu, N. E. Mavromatos, and P. Villanueva-Perez, Consistent Probabilistic Description of the Neutral Kaon System: Novel Observable Effects, Phys. Lett. B724, 269 (2013), arXiv:1208.3572 [hep-ph]
- [43]
- [44]
-
[45]
F. Benatti and R. Floreanini, Open system approach to neutrino oscillations, JHEP02, 032, arXiv:hep- ph/0002221
- [46]
- [47]
-
[48]
Ohlsson, Equivalence between neutrino oscillations and neutrino decoherence, Phys
T. Ohlsson, Equivalence between neutrino oscillations and neutrino decoherence, Phys. Lett. B502, 159 (2001), arXiv:hep-ph/0012272
- [49]
- [50]
- [51]
- [52]
-
[53]
V. De Romeri, C. Giunti, T. Stuttard, and C. A. Ternes, Neutrino oscillation bounds on quantum decoherence, JHEP09, 097, arXiv:2306.14699 [hep-ph]
-
[54]
E. Akhmedov and A. Y. Smirnov, Damping of neutrino oscillations, decoherence and the lengths of neutrino wave packets, JHEP11, 082, arXiv:2208.03736 [hep-ph]
-
[55]
J. Alexandre, K. Farakos, N. E. Mavromatos, and P. Pasipoularides, Neutrino oscillations in a stochastic model for space-time foam, Phys. Rev. D77, 105001 (2008), arXiv:0712.1779 [hep-ph]
-
[56]
Neutrino decoherence from quantum gravitational stochastic perturbations,
T. Stuttard and M. Jensen, Neutrino decoherence from quantum gravitational stochastic perturbations, Phys. Rev. D102, 115003 (2020), arXiv:2007.00068 [hep-ph]
-
[57]
Kopp and G
J. Kopp and G. A. Parker,nuDICEPython package (2026),https://github.com/george-parker/nuDICE
2026
-
[58]
C. W. Kim and W. P. Lam, Some remarks on neutrino decay via a Nambu-Goldstone boson, Mod. Phys. Lett. A5, 297 (1990)
1990
- [59]
- [60]
- [61]
-
[62]
Gorini, A
V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, Completely Positive Dynamical Semigroups of N Level Systems, J. Math. Phys.17, 821 (1976)
1976
-
[63]
Lindblad, On the Generators of Quantum Dynamical Semigroups, Commun
G. Lindblad, On the Generators of Quantum Dynamical Semigroups, Commun. Math. Phys.48, 119 (1976)
1976
-
[64]
T. F. Havel, Robust procedures for converting among Lindblad, Kraus and matrix representations of quantum dynamical semigroups, J. Math. Phys.44, 534 (2003)
2003
- [65]
-
[66]
Choi, Completely positive linear maps on complex matrices, Linear Algebra Appl.10, 285 (1975)
M.-D. Choi, Completely positive linear maps on complex matrices, Linear Algebra Appl.10, 285 (1975)
1975
-
[67]
Jamio lkowski, Linear transformations which preserve trace and positive semidefiniteness of operators, Rept
A. Jamio lkowski, Linear transformations which preserve trace and positive semidefiniteness of operators, Rept. Math. Phys.3, 275 (1972)
1972
- [68]
-
[69]
Decaying Sterile Neutrinos and the Short Baseline Oscillation Anomalies,
M. Dentler, I. Esteban, J. Kopp, and P. Machado, De- caying Sterile Neutrinos and the Short Baseline Oscil- lation Anomalies, Phys. Rev. D101, 115013 (2020), arXiv:1911.01427 [hep-ph]
- [70]
- [71]
-
[72]
Weckesser,odeintwPython package (2025),https: //github.com/WarrenWeckesser/odeintw
W. Weckesser,odeintwPython package (2025),https: //github.com/WarrenWeckesser/odeintw
2025
-
[73]
A. C. Hindmarsh, ODEPACK, a systematized collection of ODE solvers, inScientific Computing, edited by R. S. Stepleman (North-Holland, Amsterdam, 1983) pp. 55–64, https://computing.llnl.gov/projects/odepack
1983
-
[74]
Virtanenet al., SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python, Nature Methods17, 261 (2020)
P. Virtanenet al., SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python, Nature Methods17, 261 (2020)
2020
-
[75]
A. Abuslemeet al.(JUNO), JUNO physics and de- tector, Prog. Part. Nucl. Phys.123, 103927 (2022), arXiv:2104.02565 [hep-ex]
-
[76]
Abeet al.(KamLAND), Precision Measurement of Neutrino Oscillation Parameters with KamLAND, Phys
S. Abeet al.(KamLAND), Precision Measurement of Neutrino Oscillation Parameters with KamLAND, Phys. Rev. Lett.100, 221803 (2008), arXiv:0801.4589 [hep-ex]
-
[77]
J. Wanget al.(JUNO), Damping signatures at JUNO, a medium-baseline reactor neutrino oscillation experiment, JHEP06, 062, arXiv:2112.14450 [hep-ex]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.