{"id":"7f8f8ab6-0acb-4ecc-aea3-2d1c16adcad1","arxiv_id":"2411.13634","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Using time-of-flight delays of the three neutrino mass eigenstates across sharp supernova features, a galactic supernova at 10 kpc could individually constrain or measure neutrino masses at JUNO, with precision that depends strongly on the assumed mass scale.","lead":"A nearby supernova would let a large neutrino detector like JUNO clock how fast each of the three neutrino mass states travels, potentially measuring their individual masses. The paper works out the general equations and forecasts JUNO's sensitivity using sharp bursts in the supernova signal as timing markers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"High-mass benchmark contradicts 'each mass measured': Table II gives only upper bounds for m1, so the conclusion that each mass state can be uniquely identified is not supported by the paper's own results.","rationale":"The reader's verdict is CONDITIONAL with the weakest assumption identified as the fixed-mixing ansatz. I instead find a more concrete, internal problem: the paper's headline conclusion is contradicted by its own Table II for the high-mass benchmark, where m1 is an upper limit. This is not a matter of wording; it affects the central claim of individually measuring all three masses. The mixing ansatz is an explicitly stated assumption and a reasonable first step, whereas the m1 upper bound is a direct inconsistency between the numerical results and the conclusions. The proposed test quantifies at what confidence m1 is excluded; if it is not, the conclusion must be softened. This reinforces the reader's CONDITIONAL verdict but for a different reason.","tokens_in":22469,"tokens_out":25870,"duration_ms":1062929,"concrete_test":"Recompute the 1D profile likelihood for m1 in the 'High mass' benchmark (Table I) using a Feldman-Cousins or full Monte Carlo interval at 90% and 99% confidence; if m1=0 remains allowed at these levels, then the Section VIII conclusion that 'each separate mass state can be uniquely identified' is not supported, and the conclusions should be revised to state that only m2 and m3 are individually measured in this scenario.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central prose claim (Section VIII) is that 'if new physics indicates that neutrino masses are different elsewhere in the galaxy, as may be the case due to DM interactions, each separate mass state can be uniquely identified.' Table II, however, shows that in the 'High mass' benchmark (m1=0.2, m2=1.0, m3=1.8 eV) the 1σ constraint on m1 is an upper bound in every detection channel: m1 < 0.40 eV (QCD phase transition), m1 < 0.54 eV (black hole formation), and m1 < 0.80 eV (neutronization burst). These intervals do not exclude m1=0, so the lightest mass state is not measured even in the benchmark specifically constructed to display individual mass sensitivity. The likely reason is that for m1=0.2 eV the time delay at 10 MeV is ~0.2 ms, which is degenerate with the unknown emission time of the sharp feature; the fit absorbs this into the start-time nuisance parameter. Thus the conclusion overclaims: only m2 and m3 are individually resolved in the new-physics benchmark, and the 'measure every mass' statement is at best valid for the degenerate KATRIN-scale masses, and only for the speculative QCD phase-transition feature.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops a framework for using a galactic supernova observed by JUNO to constrain the three individual neutrino masses. The key theoretical step is rewriting the standard MSW-flavor-transformed supernova fluxes as sums over the three mass eigenstates (Eqs. 11-14) and applying an energy-dependent time delay to each state separately (Eqs. 27-28), with the mass-state labels assigned by the electron-row PMNS elements rather than by mass. The authors check that jump probabilities are negligible for the relevant parameter space, then simulate event rates for three sharp timing features (neutronization burst, QCD phase transition, black-hole formation) using a single 27 solar-mass fiducial model and compute 1-sigma constraints for several benchmark mass scenarios (Table II), including a speculative new-physics scenario with large masses. They conclude that a galactic supernova can measure the mass of each neutrino state individually.","tokens_in":22739,"tokens_out":11292,"duration_ms":125915,"significance":"If the standard-case sensitivity holds, the paper offers a new, orthogonal probe of the absolute neutrino mass scale, and the generalized treatment of separate mass states in Eqs. (11)-(14) is a useful conceptual contribution for future supernova analyses. Strengths include the transparent derivation of the flux formulas, a careful numerical check of jump probabilities, explicit treatment of arbitrary mass orderings through the H/M/L notation, and the release of customized SNEWPY code. The falsifiable prediction—that JUNO would observe energy-dependent delays of individual mass states after a sharp feature—is clearly stated. However, the headline claim that each mass state can be measured is not supported by the paper's own Table II for most benchmarks, and the new-physics scenario rests on an unexamined ansatz about the PMNS matrix. The central physics is sound and the shortcomings are addressable in revision.","major_comments":[{"comment":"The conclusion that 'a detection of a galactic SN can measure the mass of each neutrino state' is stronger than the numerical results in Table II. In the Oscs:NO and Planck:NO rows every entry is a one-sided upper limit; in the High mass row m1 is only bounded from above (m1 < 0.40 eV for QCD, m1 < 0.54 eV for black-hole formation, m1 < 0.80 eV for the neutronization burst), and none of these intervals excludes m1 = 0. Two-sided intervals for all three masses appear only in the KATRIN benchmark, and only for the QCD phase-transition feature. The abstract and conclusions should be reworded to state which scenarios yield measurements and which yield upper bounds; otherwise the central advertised capability is overstated.","section":"Section VIII, Table II"},{"comment":"The new-physics benchmark in Table I and the associated conclusion that 'each separate mass state can be uniquely identified' rest on the ansatz stated after Eq. (9) that the PMNS matrix is unchanged when neutrino masses vary. This ansatz is not derived or tested. If neutrino masses are modified by a coupling to dark matter, as Section II motivates, the mass matrix—and hence the mixing angles in the galactic center—can also be modified; the paper's own statement in Section II that spatially varying masses during propagation are not modeled adds a further unquantified effect. Without either a concrete model in which the PMNS matrix is exactly environment-independent or a demonstration that the Table II intervals are robust to O(1) changes in |U_ei|^2, the claimed new-physics sensitivities are not established.","section":"Section III.A, Section II"},{"comment":"The time-delay formula uses the ultra-relativistic approximation D (m_i/E)^2 / 2. For the High mass benchmark (Table I: m2 = 1.0 eV, m3 = 1.8 eV) and for the 5 eV illustrative cases in Figs. 5 and 10, the relevant low-energy events have m_i/E not much smaller than unity; at E = 5 MeV the exact delay D (1/sqrt(1 - (m_i/E)^2) - 1) differs from the approximation by about 10% for m_i = 1.8 eV and much more at lower energies. Since the sensitivity to masses is driven by low-energy events, the numerical implementation via Eq. (28) should either use the exact time delay or impose a validity cut. Otherwise the high-mass rows of Table II and the corresponding corner plots are biased.","section":"Eq. (27), Section VI.B"},{"comment":"The quantitative 1-sigma intervals in Table II are computed for one fiducial 27 solar-mass model from Ref. [21] and for fixed parameterized shapes of the QCD and black-hole features; the only systematic effect marginalized is the feature start time. The paper acknowledges in Sections III.C and IV that supernova models and feature details are uncertain, but it does not propagate these uncertainties. A robustness study varying the supernova model and the feature width, amplitude, and turn-off times is needed before the Table II intervals can be read as realistic predictions. As written, the numbers are conditional on the assumed feature and model.","section":"Section VI.B, Section III.C"}],"minor_comments":[{"comment":"The statements that the relevant angle 'is probably θ13' and that 'the other factor is most likely' the given expression are vague; since the numerical conclusion depends on jump probabilities being negligible, these three-flavor mapping choices should be stated more definitively or explicitly verified.","section":"Section III.B.2"},{"comment":"The symbol s is used for sin θ in Eq. (18), while s2ij is used for sin(2θij) elsewhere in the paper; this overloaded notation could confuse readers. Using sin θ in Eq. (18) would remove the ambiguity.","section":"Eq. (18)"},{"comment":"The 'Heaviest by Planck' benchmark values are rounded to three decimals, so the implied squared-mass splittings are not exactly consistent with the quoted oscillation parameters; specifying the exact masses used in the simulation would improve reproducibility.","section":"Table I"},{"comment":"The sentence 'some of the low-energy events are somewhat delayed due to m3 large, but the effect of the remaining events is smaller' is vague; it would be clearer to specify which mass state produces which visible delay feature in the panel.","section":"Fig. 5 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid phenomenology study with a transparent derivation and useful code release. The main issue is that the headline claim needs to be aligned with Table II, and the new-physics scenario needs either a concrete model or a robustness test of the fixed-PMNS ansatz. I would be comfortable with publication after a major revision addressing these points."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first treatment I've seen of a galactic SN time-delay that keeps the three mass eigenstates separate instead of folding them into a single effective mass. Equations (11)-(14) are the real work, and they're correct—the MSW/jump-probability check is careful and the authors even provide a customized SNEWPY fork, which makes the numerics reproducible. For anyone planning a SN neutrino analysis, those equations and the associated time-delay implementation (eq. 28) are worth taking seriously.\n\nThe JUNO sensitivity study is a standard forecast: event rates are generated with injected benchmark masses and then fit, so the intervals in Table II are expected closures, not measurements. That's fine for a forecast, but it means the numbers inherit the chosen 27 M☉ model and the assumed sharp-feature shapes. The authors are upfront about marginalizing over the feature start time, and they don't pretend the QCD phase transition is guaranteed—but they also don't propagate the SN model uncertainty into the sensitivities. That's a real limitation, not a fatal one.\n\nWhere I part company with the paper's conclusions is the title-level claim. Table II shows that in the 'High mass' new-physics benchmark, m1 is constrained to an upper bound in every channel (m1 < 0.40 eV for QCD, <0.54 for BH, <0.80 for NB). An upper bound is not a measurement of m1. The only benchmark where all three masses get two-sided intervals is the KATRIN-scale case with a QCD phase transition. So the abstract's 'masses of each of the three mass eigenstates individually' and the conclusion's 'each separate mass state can be uniquely identified' are too strong. The stress-test note on this is right, and it lands on reading the paper.\n\nThe other soft spot is the ansatz in Section III.A: the mixing matrix stays at the terrestrial values while the masses vary arbitrarily. For a DM-coupled new-physics scenario in the galactic center that's an assumption, not a consequence. The authors flag it, but it means the high-mass benchmark program—the part that would actually measure individual masses—rests on unvalidated physics.\n\nNet: the derivation is solid and the framework is a real step forward. The paper deserves a serious referee. It needs the conclusions and abstract toned down to match Table II, a discussion of model-systematic propagation, and ideally a clearer statement that for the standard light-mass case this is an upper-bound probe, not a measurement. I'd send it to review and would cite the three-state time-delay equations.","headline":"A genuine three-mass-state framework for SN time-delay, with a clean derivation and honest caveats, but the conclusion overclaims: Table II only measures all three masses in the KATRIN/QCD corner, and the new-physics benchmark leaves m1 as an upper bound.","tokens_in":23340,"tokens_out":2664,"would_cite":true,"duration_ms":28553,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.60.Pq","95.85.Ry","97.60.Bw"],"model":"deepseek-v4-flash","headline":"A galactic supernova could measure each of the three neutrino mass states individually.","keywords":["neutrino mass eigenstates","supernova neutrinos","time-of-flight delay","MSW effect","JUNO detector","neutronization burst","QCD phase transition","black hole formation"],"falsifier":"Measure the energy-dependent arrival delays of a sharp supernova feature (for example, the black-hole turnoff) in a galactic supernova at a known distance with JUNO. If the low-energy tail of events is not delayed relative to the sharp cutoff by the $\\propto 1/E^2$ pattern with the per-state weights 0.68/0.30/0.02, or if the delays imply mass-squared splittings that disagree with the oscillation-measured $\\Delta m^2_{21}=7.4\\times10^{-5}$ eV$^2$ and $\\Delta m^2_{31}=2.5\\times10^{-3}$ eV$^2$, the per-state time-delay picture would be falsified.","tokens_in":22196,"feed_emoji":"🔭","tokens_out":11999,"duration_ms":116129,"temperature":0.7,"pith_summary":"This paper claims that one nearby supernova, seen by the JUNO detector, can constrain the mass of each of the three neutrino mass eigenstates individually, rather than only the sum of masses or the squared mass differences that oscillation experiments measure. The idea is that each mass state travels at a slightly different speed, so a sharp feature in the supernova's neutrino signal — the neutronization burst, a QCD phase transition, or black hole formation — arrives as three time-staggered copies whose delays grow at low energy. The paper writes the supernova flavor flux with the three mass states kept separate (equations 11–14) and folds in the energy-dependent time delay of equation (27) for each state. For KATRIN-scale or heavier masses, all three masses can be measured; at the light masses allowed by oscillations or cosmology, the same data give upper bounds rather than measurements. This matters because the absolute scale and individual ordering of neutrino masses are among the last unknown parameters of the standard model and are inaccessible to oscillation experiments alone.","feed_headline":"A supernova could weigh each neutrino mass on its own","feed_subtitle":"Timing sharp supernova features lets JUNO separate the three mass states instead of measuring only their sum.","key_machinery":"The load-bearing object is the set of per-mass-state flux equations (11)–(14): the supernova flux at Earth is written as a sum over the heavy, medium, and light mass eigenstates, ordered by $m_H > m_M > m_L$, with each state's mixing weight shifted by its own time delay $\\Delta t_i(E) \\simeq (D/2)(m_i/E)^2$ (equation 28). What makes it work is that the mixing weights along the electron row (0.68, 0.30, 0.02) are never collapsed through unitarity into $1-|U_{e3}|^2$, so each mass state remains separately addressable in time. The accompanying checks — production above all MSW resonances (equation 15) and negligible jump probabilities (equations 16–22, figures 2–4) — keep the mapping from source emission to mass state to detected flavor clean enough for the time-of-flight effect to be read off.","core_discovery":"On its own terms, the paper shows that the standard supernova-neutrino flux formulas change character once the three mass eigenstates are followed separately. The familiar equations for $\\Phi_{\\nu_e}$, $\\Phi_{\\bar\\nu_e}$, and the non-electron flavors are rewritten using the heavy/medium/light mass labels $H, M, L$ and the projections $|U_{eH}|^2, |U_{eM}|^2, |U_{eL}|^2$ (equations 11–14), so that each mass state carries its own mixing weight and its own arrival time. Each weight is then shifted in time by $\\Delta t_i(E) \\simeq (D/2)(m_i/E)^2$ (equation 27), and the paper shows with resonance-density and jump-probability estimates that neutrinos are produced above all level crossings, so this single-flavor mapping is valid for any mass ordering. Using the 27 $M_\\odot$ supernova model of [21], JUNO's detector response, and a Poisson log-likelihood over the three sharp timing features, the paper finds (Table II) that at 10 kpc the KATRIN-scale benchmark yields 1$\\sigma$ ranges such as $m_1 = 0.25$–$0.54$ eV, $m_2 = 0.32$–$0.84$ eV, and $m_3 = 0.34$–$1.52$ eV for a QCD phase transition, while the light oscillation/Planck benchmarks lead to upper bounds. The conclusion is stated plainly: a detection of a galactic supernova can measure the mass of each neutrino state, subject to statistical uncertainties.","pith_inferences":["If the fixed-mixing ansatz survives scrutiny, the same per-state time-delay logic should extend to any sharp neutrino transient with a known emission time — neutron-star merger signals, for example — making the technique a general probe of absolute neutrino mass rather than a supernova-specific one.","A multidetector analysis combining JUNO with the other next-generation water and liquid-argon detectors would split the difference between low thresholds and high statistics, and the paper's per-state flux formulas give those collaborations a ready-made framework; the paper itself only projects one detector.","The most interesting outcome would be a supernova whose time-delay pattern and MSW flavor pattern demand different mass orderings; under the paper's assumptions that cannot happen, so such an event would be a clean test of spatially varying neutrino masses and would falsify the fixed-mixing ansatz.","Because the delay scales as $m_i^2/D$, an independent distance measurement to the supernova (for example, from gravitational waves or the expanding shock front) is as important as event statistics for turning a detected delay into an absolute mass; the paper fixes $D = 10$ kpc rather than marginalizing over it."],"forward_implications":["For a supernova at 10 kpc observed by JUNO, each of the three mass states is constrained independently: KATRIN-scale masses are measured in all three states under a QCD phase transition, while oscillation- or Planck-scale masses give upper bounds such as $m_1 < 0.28$ eV in the normal-ordering benchmark.","Closer supernovae sharpen the measurement: at 1 kpc or less the ordering $m_3 > m_2$ becomes distinguishable even in the KATRIN scenario where it cannot be resolved at 10 kpc.","Because the three timing features come from unrelated physics, their likelihoods can be added at each point in parameter space, so a supernova that shows both a neutronization burst and a later sharp feature yields combined constraints stronger than any single feature.","The paper's per-state description implies that previous supernova analyses that assumed degenerate masses were missing the main handle: the distinct arrival-time patterns of each eigenstate, not just a single effective mass.","Part of the mass-ordering sensitivity comes from the MSW flavor conversion inside the supernova, so the same event can distinguish normal from inverted ordering even when the absolute masses are too small to time-resolve directly."],"supporting_citations":[{"why":"Supplies the standard MSW flavor-transformation fluxes (equations 1–8) that the paper rewrites so each mass eigenstate keeps its own weight.","marker":"[45]"},{"why":"Provides the fiducial 27 $M_\\odot$ supernova model and the neutronization burst used for all event-rate and sensitivity numbers.","marker":"[21]"},{"why":"Sets the earlier supernova time-delay mass analysis whose equal-mass assumption the paper replaces with separate mass states.","marker":"[13]"},{"why":"Specifies the JUNO detector configuration (volume, thresholds, timing, energy resolution) assumed throughout the sensitivity projections.","marker":"[20]"},{"why":"Gives the improved jump-probability formula with the $f$ correction used to show the resonances stay adiabatic.","marker":"[56]"},{"why":"Provides the three-flavor corrections to the two-flavor jump-probability integrands used in the numerical adiabaticity check.","marker":"[57]"},{"why":"Provides the black-hole-formation neutrino flux turn-off model used as one of the three sharp timing features.","marker":"[31]"},{"why":"Supplies the detector channels, cross sections, and energy smearing used to convert supernova fluxes into JUNO event rates.","marker":"[72]"},{"why":"Supplies the dark-matter-mediated long-range force model that motivates the high-mass, spatially varying neutrino mass benchmarks.","marker":"[35]"}],"fun_headline_variants":["Supernova timing weighs each neutrino mass separately","JUNO to measure all three neutrino masses from a supernova","One supernova, three individual mass measurements","Galactic supernova reveals individual neutrino masses"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's strongest loaded assumption is that the neutrino mixing matrix in the galactic center is identical to the one measured on Earth — the electron-row weights 0.68, 0.30, and 0.02 — even while the three masses themselves are allowed to change freely; if a dark-matter coupling alters the mixing angles together with the masses, the mapping from detected flavors to time-delayed mass states breaks down.","fun_headline_variants_meta":{"raw":{"variants":["Supernova timing weighs each neutrino mass separately","JUNO to measure all three neutrino masses from a supernova","One supernova, three individual mass measurements","Galactic supernova reveals individual neutrino masses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000264,"raw_usage":{"total_tokens":1665,"prompt_tokens":1069,"completion_tokens":596,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":685,"completion_tokens_details":{"reasoning_tokens":536}},"tokens_in":685,"tokens_out":596,"duration_ms":6865,"temperature":1.0,"reasoning_tokens":536,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:02:15.262256+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the energy-dependent arrival delays of a sharp supernova feature (for example, the black-hole turnoff) in a galactic supernova at a known distance with JUNO. If the low-energy tail of events is not delayed relative to the sharp cutoff by the $\\propto 1/E^2$ pattern with the per-state weights 0.68/0.30/0.02, or if the delays imply mass-squared splittings that disagree with the oscillation-measured $\\Delta m^2_{21}=7.4\\times10^{-5}$ eV$^2$ and $\\Delta m^2_{31}=2.5\\times10^{-3}$ eV$^2$, the per-state time-delay picture would be falsified.","supporting_citations":[{"cited_title":"A Return To Neutrino Normalcy","cited_arxiv_id":"2003.04319","evidence_quote":"Gives the improved jump-probability formula with the $f$ correction used to show the resonances stay adiabatic."},{"cited_title":"Parameter symmetries of neutrino oscillations in vacuum, matter, and approximation schemes","cited_arxiv_id":"2106.12436","evidence_quote":"Provides the three-flavor corrections to the two-flavor jump-probability integrands used in the numerical adiabaticity check."}],"review_version":1}