{"id":"dce9b7e8-8e5f-4115-8dce-20313d979158","arxiv_id":"2411.19303","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Derives a generalized Lindblad master equation for neutrino evolution with momentum-changing decay and absorption of massless particles, and translates a reactor decoherence bound into a neutrino lifetime limit.","lead":"Neutrinos can lose quantum coherence by decaying or absorbing massless particles, and this paper derives a master equation that includes such processes while allowing neutrinos to change momentum during the transition. The authors use an experimental bound on decoherence from the KamLAND reactor experiment to set a new, though weak, limit on how long a neutrino can live before this decay.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The scalar-decay specialization (33)/(34) has the gain projector reversed: Eq. (34) transfers population from ν1 to ν2 (absorption), not from ν2 to ν1 (decay), so the KamLAND lifetime bound (41) does not follow from Eq. (27).","rationale":"The reader's weakest assumption focused on the rotating-wave Kronecker deltas and the dropped principal-value term in the derivation of Eq. (27). Those are legitimate concerns about unquantified approximations, but the more decisive and concrete problem is in the application section: Eq. (34) has the gain projector reversed, so the two-flavour equation does not describe the decay needed for the lifetime bound. This is independent of the RWA debate and can be checked by a simple trace-conservation computation. The general master equation (27) may still be salvageable, and the index error in Eq. (33)/(34) is plausibly a correcting typo rather than a conceptual failure. However, as written, the advertised constraint τ2/m2 > 1.83×10^-10 s/eV does not follow from the derived equation, and the dissipative matrix used to map KamLAND constraints is based on the wrong direction of population transfer. I therefore retain the reader's CONDITIONAL verdict: the paper should be accepted only after the projector ordering in Eq. (33)/(34) is fixed, the trace-conservation identity is verified, and the lifetime bound is re-derived. The RWA/principal-value issue remains a secondary concern that should also be addressed in revision.","tokens_in":98,"tokens_out":18866,"duration_ms":358086,"concrete_test":"Trace Eq. (34) over flavor and momentum. With the printed projectors, ∂_t ∫dp Trρ = -∫dp Γ^S_21(p) ρ22(p) + ∫dp ∫dq K(q,p) ρ11(q), which is not identically zero. If the gain projector is corrected to Π12ρΠ21 and the kernel satisfies ∫dp K(q,p)=Γ^S_21(q), the trace is conserved. This single analytic check decides whether the gain term describes decay or absorption and whether the KamLAND bound (41) is supported.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most load-bearing flaw is in the specialization of the general master equation to scalar decay, not in the principal-value truncation. In Section III, Eq. (33) and its two-flavour version Eq. (34) are meant to implement the decay of a heavier state i into a lighter state f. With the paper's convention Π_ab=|a⟩⟨b| (stated after Eq. (27)), the loss term in Eq. (34) is -1/2 Γ^S_21 {Π_22, ρ(|p|)}, which correctly depletes the heavier state ν2. But the gain term is written as g^2_12 ∫ ... Π_21 ρ(|q|) Π_12, i.e., |2⟩⟨1| ρ |1⟩⟨2|, which populates the heavier state ν2 from the lighter ν1: that is the inverse absorption process, not the decay ν2→ν1. The correct decay gain would be Π_12 ρ Π_21, as also required by the paper's own Eq. (37), where L1=Π12 with rate Γ^S_21. The same reversed ordering appears in Eq. (33) for every pair i>f. As written, Eq. (34) is not a valid Lindblad equation for spontaneous decay: ρ11 has no loss term, while ρ22 gains from ρ11, so the total trace ∫dp Trρ is not conserved and the claimed equality ∫ d|p| Trρ = const is violated. Therefore the dissipative matrices (39)/(40) and the constraint (41) are not consequences of Eq. (27). This is an internal inconsistency, not merely an issue of matching experimental conventions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a QFT-based open-systems framework to derive a generalized Lindblad master equation for neutrinos interacting with a bath of massless particles. The central equation (27) includes dissipative loss and gain terms that describe neutrino decay into a lighter state plus a massless particle, and the inverse absorption process, with transitions between different neutrino momenta. The authors then specialize this equation to the case of Dirac neutrino decay into a lighter neutrino and a scalar particle in the degenerate mass hierarchy, obtaining dissipative matrices (39) and (40). Using the KamLAND constraint on the decoherence parameter Gamma21 from Ref. [9], they derive a constraint on the neutrino lifetime ratio tau_2/m_2 > 1.83e-10 s/eV for visible Dirac scalar decay.","tokens_in":11329,"tokens_out":10720,"duration_ms":84163,"significance":"If correct, the derivation of Eq. (27) would provide a novel and physically useful bridge between neutrino decay and quantum decoherence, going beyond the common Lindblad treatments that ignore momentum-changing transitions. The paper also explicitly derives the dissipative operators and parameters from a microphysical Hamiltonian, which is a strength. However, the phenomenological application is compromised by an internal inconsistency in the scalar-decay reduction: the gain term in Eqs. (33) and (34) is written with the projector ordering reversed relative to the decay direction stated in Eq. (27). This error propagates into the dissipative matrices and the lifetime bound (41), so the headline constraint does not follow from the derived master equation as it stands.","major_comments":[{"comment":"The gain term for the decay of a heavier state i into a lighter state f is written with the projector ordering reversed. With the convention Pi_ab = |a><b| stated after Eq. (27), the decay i->f should have the gain term Pi_{fi} rho Pi_{if} = |f><i| rho |i><f|, which transfers population from the heavier state i to the lighter state f. As written, Eq. (33) contains Pi_{if} rho Pi_{fi}, which transfers population from the lighter state f to the heavier state i, i.e., the absorption process. This is confirmed by Eq. (34), where the gain term g^2_12 ... Pi_21 rho Pi_12 populates state 2 from state 1, while the loss term depletes only state 2; consequently rho_11 has no loss term and the trace of the dissipative part is not zero. This directly contradicts the paper's own statement that the total number of neutrinos of all momenta is conserved (paragraph after Eq. (35)). The error is internal, because Eq. (37) lists L1 = Pi_12 with rate Gamma^S_21, which is the correct decay operator and would produce the gain term Pi_12 rho Pi_21, not Pi_21 rho Pi_12. As a result, the dissipative matrices (39) and (40) and the lifetime constraint (41) are not consequences of Eq. (27).","section":"Section III, Eqs. (33) and (34)"},{"comment":"The rotating-wave approximation is implemented by inserting Kronecker delta symbols 'by hand' into Eq. (22) to obtain Eq. (24). No derivation or validity criterion is given for these deltas, which enforce diagonality conditions such as delta_{nn'} delta_{ij} (2pi)^3 delta^(4)(q-q'). Since the final master equation (27) is intended to describe transitions between different momenta, the secular approximation for the momentum-dependent phase factors needs a quantitative justification, especially because the decaying/absorbing processes involve energy differences that may not be small. Without such a derivation, the precise form of the gain terms in Eq. (27) is not established.","section":"Section II, Eq. (24)"},{"comment":"The principal-value part of the time integral in Eq. (26) is dropped with the justification that it contributes only to the coherent Hamiltonian and is 'not the focus of this paper.' This term is a systematic, energy-dependent correction to the oscillation Hamiltonian (a Lamb-shift-type term), and its magnitude is never estimated. Since Eq. (27) is presented as the complete evolution equation, the omission should be justified either by an explicit estimate of its effect on the observables of interest (such as the survival probabilities used for the KamLAND constraint) or by absorbing it into a redefined Hamiltonian with a clear statement. As it stands, the derivation of the dissipative part is incomplete.","section":"Section II, paragraph after Eq. (26)"}],"minor_comments":[{"comment":"There is a typo in the sentence 'experimentrs wigth neutrino fluxes' that should read 'experiments with neutrino fluxes.'","section":"Section III, first paragraph"},{"comment":"The statement that 'there is no superposition between states with different momentum' (Eq. (6)) is in tension with the later development of transitions between different momenta; the authors should clarify that Eq. (6) is an initial or approximate condition for the reduced density matrix, not a restriction on the dynamics.","section":"Section II, paragraph before Eq. (6)"},{"comment":"The notation for the time limits of the integrals is inconsistent: Eq. (22) writes the second time integral as ∫_{t1}^{t0} d4x2, while Eq. (24) writes ∫_{t}^{t0} d4x2 in some places and ∫_{t1}^{t0} in others. This makes it difficult to follow the Markovian-limit argument.","section":"Section II, Eqs. (22) and (24)"},{"comment":"The replacement of the integral over |q| by the decay width Gamma^S_if uses the degenerate limit, but the integration limits depend on (m_i/m_f)^2, which is not exactly 1; the accuracy of approximating rho(|q|,t) by rho(|p|,t) over this interval should be stated with a quantitative estimate.","section":"Section III, Eq. (36)"}],"recommendation":"major_revision","confidential_remarks":"The central framework in Section II is plausible and worth developing, but the error in the scalar-decay reduction is load-bearing for the main phenomenological claim. The authors should be asked to correct the projector ordering in Eqs. (33) and (34), re-derive the dissipative matrices, and recompute the lifetime constraint; they should also address the principal-value truncation and the RWA justification. The final claim of obtaining the 'first' direct constraint on visible Dirac neutrino decay is currently unsupported. I recommend major revision rather than rejection because the general formalism may be salvageable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The general framework is the reason to read this. Eq. (27), with the decay/absorption widths (28)-(31), is a genuine extension of the usual Lindblad treatment: it allows transitions between neutrino states of different momenta, and it derives the dissipative operators from a QFT open-systems setup rather than postulating them. The reduction to the standard Lindblad form in the degenerate limit is sensible, and the comparison with Moss et al. is useful. That part deserves a serious referee.\n\nThe problems are concentrated in Section III, and they are load-bearing. In Eq. (34), the gain term is written as Π21 ρ Π12. With the paper's own convention Π_ab=|a⟩⟨b|, this takes population from ν1 and puts it into ν2. That is the absorption process, not the decay ν2→ν1. Together with the loss term −(Γ21/2){Π22,ρ}, the trace of ρ is not conserved, contradicting the paper's claim that ∫d|p| Trρ is constant. Eq. (37) later uses L1=Π12, which would give the correct gain Π12 ρ Π21; so the paper contradicts itself. As written, the dissipative matrices (39)/(40) and the bound (41) do not follow from Eq. (27). This looks fixable by swapping the projector order in (33)/(34), but it has to be fixed before any phenomenological claim can stand. Separately, the D matrix in (39) does not conserve trace even for the L1=Π12 choice, so the Bloch-sphere expansion needs a careful re-check.\n\nThe derivation of the general equation also has soft spots. The rotating-wave approximation in Eq. (24) inserts Kronecker deltas by hand, without a derivation for momentum-changing transitions. The Markovian limit is stated rather than justified, and the principal-value term in (26) is dropped with no estimate of its effect on the coherent Hamiltonian. These are not fatal to the framework, but they need to be quantified before Eq. (27) can be called complete.\n\nThe lifetime bound, even if repaired, is a weak reinterpretation of an existing KamLAND constraint on Γ21: 1.83×10⁻¹⁰ s/eV is orders of magnitude weaker than the visible-decay bounds they quote themselves (6.7×10⁻⁴ s/eV for Majorana). Calling it a 'first direct constraint' overstates matters.\n\nBottom line: the general equation is worth refereeing; Section III needs major revision. I would not cite the lifetime bound, but I would cite Eq. (27) once it is cleaned up.","headline":"The generalized master equation (27) is a real addition to the neutrino-decoherence toolbox, but a reversed projector in the scalar-decay section sinks the lifetime bound as stated.","tokens_in":11851,"tokens_out":8711,"would_cite":false,"duration_ms":74526,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A generalized Lindblad master equation derived from open-systems QFT ties neutrino decoherence to decay and absorption and yields a direct decay-width limit.","keywords":["neutrino quantum decoherence","Lindblad master equation","neutrino decay","massless particle absorption","open quantum systems","neutrino lifetime bound","KamLAND"],"falsifier":"Compute the discarded principal-value term in Eq. (26) for the two-flavour scalar-decay case and compare its magnitude with the oscillation frequency $\\Delta m^2/(2E)$; if the correction is comparable for reactor energies, Eq. (27) is not the complete evolution. Alternatively, a reactor measurement of the decoherence parameter $\\Gamma_{21}(E)$ whose energy dependence deviates from the assumed $\\Gamma_{21}(E)\\propto E^{-1}$ scaling would invalidate the derived lifetime bound.","tokens_in":10756,"feed_emoji":"⚛️","tokens_out":8995,"duration_ms":67785,"temperature":0.7,"pith_summary":"The paper develops a quantum field theory of open systems for neutrinos and derives a generalized Lindblad master equation that governs the density matrix of a neutrino propagating through a bath of massless particles. The new element is that the evolution explicitly includes transitions between neutrino stationary states with different momenta, caused by the neutrino decaying into a lighter state plus a massless particle and by the reverse absorption process. The dissipative operators and rates are not put in by hand; they are obtained from the interaction Hamiltonian, the bath correlation function, and the Markovian and rotating-wave approximations. As an application, the authors specialize to visible Dirac neutrino decay into a scalar particle in the degenerate mass hierarchy, use the KamLAND decoherence bound to constrain the decay width, and obtain $\\tau_2/m_2 > 1.83 \\times 10^{-10}~\\mathrm{s/eV}$. If the derivation is right, reactor and accelerator decoherence searches become a direct probe of neutrino decay dynamics.","feed_headline":"Neutrino decay sets a 1.83e-10 s/eV lifetime bound","feed_subtitle":"KamLAND decoherence data now give a direct neutrino decay-width limit.","key_machinery":"The machinery is the open-quantum-systems expansion of the evolution operator to second order in the coupling, combined with the rotating-wave (secular) approximation and the Markovian limit. The central technical identity is (26), $(2\\pi)^3\\delta^{(3)}(\\mathbf p)\\int_0^\\infty d\\tau\\, e^{-iE_p\\tau}=\\tfrac12(2\\pi)^4\\delta^{(4)}(p)-i(2\\pi)^3\\delta^{(3)}(\\mathbf p)\\,\\mathcal P(1/E_p)$, whose real part produces the energy-conserving delta functions in the decay and absorption widths and whose imaginary part is discarded as a coherent Hamiltonian correction. The rotating-wave approximation is implemented by inserting Kronecker deltas $\\delta_{nn'}$ and momentum delta functions by hand in Eq. (24). The derived widths $\\Gamma^{d}_{ip}$, $\\Gamma^{a}_{ip}$ and transition widths $\\Gamma^{d}_{jq\\to ip}$, $\\Gamma^{a}_{jq\\to ip}$ then supply the explicit dissipative operators and parameters in Eq. (27). In the scalar-decay application, the degenerate-limit decay width $\\Gamma^{S}_{if} = (g^2_{if}/\\pi)(m^2_i-m^2_f)/|\\mathbf p|$ parametrizes the Lindblad operators and the dissipative matrix.","core_discovery":"The central claim is that neutrinos undergoing decay into a lighter mass state and a massless particle, and the inverse absorption process, evolve according to the generalized Lindblad master equation (27), which is derived from first principles rather than assumed. The equation has the Lindblad structure with an anti-commutator damping term built from the decay width $\\Gamma^d_{ip}$ and absorption width $\\Gamma^a_{ip}$, plus a feeding term that transfers population from state $|jq\\rangle$ to state $|ip\\rangle$ with different momenta, integrated over the phase space of the massless particle. The dissipative operators are the projectors $\\Pi_{ij}=|i\\rangle\\langle j|$, and the rates are explicitly given by the on-shell matrix elements of the neutrino current contracted with the bath polarization tensor and the Bose-Einstein factors $1+N(\\omega)$ and $N(\\omega)$. For the scalar-decay case in the degenerate limit, the equation reduces to the standard Lindblad form with three dissipative operators and a dissipative matrix $D^{(3)}_{kl}$ that is not diagonal, and the element $\\Gamma_{21}$ is related to the decay width $\\Gamma^S_{21}$, yielding the lifetime constraint (41).","pith_inferences":["If the same derivation is repeated for Majorana neutrinos, the absorption terms and spin structure would differ, likely yielding a different lifetime bound; the paper's method provides the template for that calculation.","The discarded imaginary principal-value term may produce a momentum-dependent correction to the neutrino Hamiltonian; testing its size could either confirm the Lindblad truncation or reveal a new coherent decoherence effect.","The framework suggests a general recipe: any neutrino two-body process with a massless final or initial particle in a thermal bath generates a Lindblad-type dissipative sector whose rates are computable from on-shell amplitudes, connecting decoherence searches to decay searches across channels."],"forward_implications":["Neutrino quantum decoherence no longer needs to be parametrized ad hoc: decay and absorption widths fix the Lindblad dissipative operators and rates.","Because transitions between different momenta are included, the master equation applies to decays in neutrino fluxes from reactors, accelerators, and supernovae where momentum changes matter.","The KamLAND constraint on $\\Gamma_{21}$ translates into $\\tau_2/m_2 > 1.83\\times10^{-10}$ s/eV for visible Dirac scalar decay in the degenerate hierarchy, a direct oscillation-data limit on the decay width.","The explicitly derived dissipative matrix $D^{(3)}_{kl}$ contains off-diagonal structure, so future analyses should use it rather than the commonly assumed diagonal decoherence parameters."],"supporting_citations":[{"why":"Defines the Lindblad form (1) that the derived master equation generalizes.","marker":"[2]"},{"why":"Provides the complete-positivity framework that makes the Lindblad structure a valid open-system evolution.","marker":"[3]"},{"why":"Supplies the experimental KamLAND bound $\\Gamma_{21}(E_0)<1.8\\times10^{-24}$ GeV used to derive the lifetime constraint.","marker":"[9]"},{"why":"Earlier open-systems treatment of neutrino radiative decay that this paper extends to general massless-particle emission and absorption.","marker":"[24]"},{"why":"Standard open quantum systems theory underlying the interaction-picture expansion and trace over the environment.","marker":"[25]"},{"why":"Gives the simple Lindblad derivation and the integral identity (26) with the principal-value term the paper discards.","marker":"[26]"},{"why":"Provides the degenerate-limit scalar decay width $\\Gamma^S_{if}$ and the comparison Majorana lifetime bound.","marker":"[30]"},{"why":"Derives a related evolution equation from a non-Hermitian Hamiltonian; the paper contrasts Eq. (27) with it to highlight the absorption term and the extra feeding term.","marker":"[34]"},{"why":"Specifies the allowed energy range for the daughter neutrino used in the integration limits of Eq. (33).","marker":"[35]"}],"fun_headline_variants":["New Lindblad equation tightens neutrino lifetime bound","Neutrino decay width limited via KamLAND decoherence","First-principles master equation sets neutrino lifetime bound","Generalized Lindblad equation yields neutrino decay limit","Neutrino lifetime constraint from quantum decoherence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes that the imaginary principal-value part of the time integral, which would shift the coherent oscillation Hamiltonian, is negligible, and that the rotating-wave approximation holds for transitions between states of different momenta.","fun_headline_variants_meta":{"raw":{"variants":["New Lindblad equation tightens neutrino lifetime bound","Neutrino decay width limited via KamLAND decoherence","First-principles master equation sets neutrino lifetime bound","Generalized Lindblad equation yields neutrino decay limit","Neutrino lifetime constraint from quantum decoherence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000168,"raw_usage":{"total_tokens":1273,"prompt_tokens":968,"completion_tokens":305,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":584,"completion_tokens_details":{"reasoning_tokens":230}},"tokens_in":584,"tokens_out":305,"duration_ms":4018,"temperature":1.0,"reasoning_tokens":230,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:19:22.892203+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the discarded principal-value term in Eq. (26) for the two-flavour scalar-decay case and compare its magnitude with the oscillation frequency $\\Delta m^2/(2E)$; if the correction is comparable for reactor energies, Eq. (27) is not the complete evolution. Alternatively, a reactor measurement of the decoherence parameter $\\Gamma_{21}(E)$ whose energy dependence deviates from the assumed $\\Gamma_{21}(E)\\propto E^{-1}$ scaling would invalidate the derived lifetime bound.","supporting_citations":[{"cited_title":"Lindblad, On the Generators of Quantum Dynamical Semigroups, Commun","cited_arxiv_id":null,"evidence_quote":"Defines the Lindblad form (1) that the derived master equation generalizes."},{"cited_title":"Stankevich and A","cited_arxiv_id":null,"evidence_quote":"Earlier open-systems treatment of neutrino radiative decay that this paper extends to general massless-particle emission and absorption."},{"cited_title":"Breuer and F","cited_arxiv_id":null,"evidence_quote":"Standard open quantum systems theory underlying the interaction-picture expansion and trace over the environment."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the simple Lindblad derivation and the integral identity (26) with the principal-value term the paper discards."}],"review_version":1}