arxiv: 2604.16652 · v2 · submitted 2026-04-17 · ✦ hep-ph · astro-ph.HE· hep-ex

Recognition: unknown

Investigating the Neutrino Mass Ordering Problem via Ternary Plots

Alexander Migala , Kate Scholberg

Authors on Pith no claims yet

Pith reviewed 2026-05-10 07:40 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.HEhep-ex

keywords neutrino mass orderingsupernova neutrinosternary plotsflavor compositioncore-collapse supernovaeneutrino detectorsmass hierarchy

0 comments

The pith

Ternary plots of supernova neutrino flavors separate normal from inverted mass ordering across models.

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

The paper examines whether neutrinos from core-collapse supernovae can reveal the unknown neutrino mass ordering. It applies ternary plots to show how the mix of neutrino flavors evolves during the burst for multiple models. A simplified unfolding step accounts for detector responses before checking the plotted paths. The analysis finds that normal and inverted orderings consistently land in different parts of the ternary space, even when the supernova model changes.

Core claim

We explore what may be deduced about the neutrino mass ordering problem from the observation of core-collapse supernova burst neutrinos in modern terrestrial detectors. We employ ternary plots in a novel way to visualize the time evolution of the flavor composition of various supernova neutrino flux models. Through our analysis of several models using a simplified unfolding process, we have explored potential robust discriminants between the normal and inverted mass orderings. We find that the normal and inverted mass orderings tend to occupy different regions in ternary space across different models.

What carries the argument

Ternary plots that display the time-dependent fractions of electron, muon, and tau neutrinos in the supernova flux.

If this is right

Supernova neutrino data could distinguish the two mass orderings without depending on details of any single model.
The regional separation in ternary space persists after the simplified unfolding step that models detector response.
Multiple independent supernova models all exhibit the same pattern of separation between orderings.
Future terrestrial detectors could use this visual method as one route to ordering information from a burst event.

Where Pith is reading between the lines

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

The separation points to a difference in how the two orderings shape the flavor ratios that reach Earth.
Combining ternary-plot results with data from long-baseline oscillation experiments could cross-check the ordering.
The same plotting approach might be tested on other high-energy neutrino transients once suitable models exist.

Load-bearing premise

The simplified unfolding process accurately captures detector effects without introducing biases that would erase the separation between orderings.

What would settle it

A real supernova detection whose unfolded data places normal and inverted ordering trajectories in overlapping regions of the same ternary plot would show that the separation does not survive actual observations.

Figures

Figures reproduced from arXiv: 2604.16652 by Alexander Migala, Kate Scholberg.

**Figure 2.** Figure 2: Ternary diagram for the truth flux of Nakazato index 0 and 2 (see App. B) under NMO/IMO [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

**Figure 3.** Figure 3: For CC interactions on nuclei, the final-state leptons tend to retain most of the energy of the neutrino, which is observable with reasonable resolution in many detectors. For scattering interactions, energy is lost; however, the recoil energy distribution is well understood so that information on the neutrino energy spectrum can be determined statistically. In general, practical detector time resoluti… view at source ↗

**Figure 3.** Figure 3: The computed supernova signal as seen in different detection technologies for a Nakazato flux [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗

**Figure 4.** Figure 4: Flux profiles for the Nakazato s0 model under different oscillation prescriptions. The time [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗

**Figure 5.** Figure 5: Nakazato submodel 0 channel event rate evolution for both IMO (right) and NMO (left); axes [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗

**Figure 6.** Figure 6: Unfolded trajectories for Nakazato 2013. By row: model index 0, 2 (see Appendix B). For [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗

**Figure 7.** Figure 7: Unfolded trajectories for Nakazato 2013, plotted in the same way as Fig. 6. By row: model index [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗

**Figure 8.** Figure 8: Unfolded trajectories for Warren 2020 [41] (see Appendix B), plotted in the same way as Fig. 6. [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗

**Figure 10.** Figure 10: Superimposed cumulatively summed unfolded flux for five Warren [41] models, plotted as in [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗

**Figure 11.** Figure 11: Superimposed cumulatively summed unfolded flux for five Zha [42] models, plotted as in [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗

read the original abstract

We explore what may be deduced about the neutrino mass ordering problem from the observation of core-collapse supernova burst neutrinos in modern terrestrial detectors. We employ ternary plots in a novel way to visualize the time evolution of the flavor composition of various supernova neutrino flux models from the SNEWPY software package. Through our analysis of several models using a simplified unfolding process, we have explored potential robust discriminants between the normal and inverted mass orderings. We find that the normal and inverted mass orderings tend to occupy different regions in ternary space across different models.

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

Ternary plots of unfolded supernova neutrino fluxes suggest a separation between normal and inverted mass orderings, though the simplified unfolding step requires closer scrutiny for robustness.

read the letter

The main point is that this paper applies ternary plots to time-evolving flavor fractions from SNEWPY supernova models and reports that normal and inverted mass orderings land in different regions after a simplified unfolding step meant to approximate detector response. The plots make the time dependence visible as paths or clusters in the triangle, and the separation appears across the models they examined. This is a straightforward visual approach to an important open question in neutrino physics. The work is honest in its use of public models and does not hide behind complicated fitting. It gives a clear picture of how the flavor composition evolves and where the two orderings sit relative to each other. That visual clarity is the useful part. The limitation sits in the unfolding. Because the process is described as simplified, it likely skips full energy resolution, flavor misidentification, and background effects. Those omissions could shift the points enough to reduce or remove the apparent separation once a more complete detector simulation is used. The abstract supplies no overlap metrics, sensitivity tests, or error bands, so the separation looks suggestive but not yet demonstrated to be stable. The stress-test concern about bias is therefore on target until the full details are checked. This is for people in the supernova neutrino and mass-ordering community who want a new way to look at burst data. A reader interested in quick visual comparisons across models will find the figures helpful even if the claim needs more support. It is not a finished analysis, but the idea is clear enough and the topic important enough that it should go to peer review so referees can examine the unfolding procedure and request the missing quantitative checks.

Referee Report

2 major / 2 minor

Summary. The manuscript explores the use of ternary plots to visualize the time evolution of neutrino flavor composition from core-collapse supernova models generated via the SNEWPY package. After applying a simplified unfolding procedure to incorporate detector effects, the authors report that normal and inverted neutrino mass orderings tend to occupy distinct regions in ternary space across the models examined, suggesting a potential visual discriminant for the mass ordering problem.

Significance. If the separation proves robust under detailed detector modeling, the ternary-plot approach would supply an intuitive, model-independent visualization tool for supernova neutrino data analysis that complements traditional spectral methods. The explicit use of public SNEWPY models is a clear strength, supporting reproducibility and allowing direct comparison with future work. The current significance remains provisional because the central separation result rests on an incompletely specified unfolding step and lacks quantitative metrics.

major comments (2)

[Methods (unfolding procedure)] The claim that normal and inverted orderings occupy different regions rests on the simplified unfolding process (described in the methods section following the SNEWPY model introduction). No equations, efficiency functions, or validation against full Monte Carlo are supplied, so it is impossible to determine whether omitted energy-dependent flavor misidentification or background terms would shift the loci into overlap, directly undermining the robustness asserted in the abstract.
[Results] The results section states that the orderings 'tend to occupy different regions' but reports neither quantitative separation measures (e.g., centroid distances, overlap integrals, or Kolmogorov-Smirnov statistics between the two sets of points) nor error bands arising from model variations or unfolding assumptions. Without these, the visual impression cannot be translated into a falsifiable statement about discrimination power.

minor comments (2)

Figure captions should explicitly list the SNEWPY model names and the precise unfolding assumptions used for each panel to allow readers to reproduce the plotted points.
The abstract and introduction would benefit from a brief statement of the specific supernova detectors (e.g., Hyper-Kamiokande, DUNE) whose response is approximated by the unfolding, together with a reference to the relevant detector technical design reports.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and valuable feedback on our manuscript. We address each of the major comments below and outline the revisions we plan to make.

read point-by-point responses

Referee: [Methods (unfolding procedure)] The claim that normal and inverted orderings occupy different regions rests on the simplified unfolding process (described in the methods section following the SNEWPY model introduction). No equations, efficiency functions, or validation against full Monte Carlo are supplied, so it is impossible to determine whether omitted energy-dependent flavor misidentification or background terms would shift the loci into overlap, directly undermining the robustness asserted in the abstract.

Authors: We agree that the unfolding procedure requires more detailed specification to allow proper assessment of its robustness. In the revised version, we will provide the explicit equations for the simplified unfolding, describe the efficiency functions employed, and elaborate on the assumptions regarding energy-dependent effects and backgrounds. We will also discuss potential limitations and how they might affect the observed separation, while noting that a comprehensive Monte Carlo simulation is outside the current scope but could be addressed in future work. revision: yes
Referee: [Results] The results section states that the orderings 'tend to occupy different regions' but reports neither quantitative separation measures (e.g., centroid distances, overlap integrals, or Kolmogorov-Smirnov statistics between the two sets of points) nor error bands arising from model variations or unfolding assumptions. Without these, the visual impression cannot be translated into a falsifiable statement about discrimination power.

Authors: We acknowledge the need for quantitative support to the visual observations. We will revise the results section to include quantitative separation metrics, such as the distances between centroids of the normal and inverted ordering point distributions in ternary space, as well as overlap integrals. We will also incorporate error bands or shaded regions representing the spread due to different SNEWPY models and variations in the unfolding assumptions. This will allow for a more rigorous evaluation of the discrimination power. revision: yes

Circularity Check

0 steps flagged

No circularity: external models and visualization yield independent empirical observation

full rationale

The paper processes external public SNEWPY supernova neutrino flux models through a simplified unfolding step and visualizes flavor composition evolution in ternary plots. No load-bearing equations, fitted parameters, or self-citations reduce the claimed separation between normal and inverted orderings to a definitional or constructed result. The separation is reported as an observed pattern across the models rather than a prediction forced by the analysis inputs themselves. The derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based solely on abstract; no explicit free parameters, new entities, or non-standard axioms are stated. The work assumes standard three-flavor neutrino oscillations and the validity of the SNEWPY supernova flux models.

axioms (1)

domain assumption Standard three-flavor neutrino oscillation physics and core-collapse supernova neutrino emission models from SNEWPY are sufficiently accurate for the purpose of this visualization study.
The analysis relies on these models without deriving or testing them.

pith-pipeline@v0.9.0 · 5378 in / 1197 out tokens · 51840 ms · 2026-05-10T07:40:13.168569+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

44 extracted references · 29 canonical work pages

[1]

pinched-thermal

Neutrino Energy Spectra At any given moment in time, the neutrino spectrum for each flavor is well described by a three-parameter expression, known as the “pinched-thermal” or “Garching” parameteri- zation [31], ϕ(Eν) =N0 (α + 1)α+1 ⟨Eν⟩Γ(α + 1) Eν ⟨Eν⟩ α × exp − (α + 1) Eν ⟨Eν⟩ (3) 4 Figure 1. Supernova neutrino flux evolution in the Nakazato index 0 mod...

2013
[2]

Navas et al

S. Navas et al. (Particle Data Group), Phys. Rev. D 110, 030001 (2024)

2024
[3]

and Michael Turner, The Early Universe (Perseus Books Publishing, LLC, 1990)

Kolb, Edward W. and Michael Turner, The Early Universe (Perseus Books Publishing, LLC, 1990)

1990
[4]

DESI DR2 Results II: Measurements of Baryon Acoustic Oscillations and Cosmological Constraints

M. Abdul Karim et al. (DESI), Phys. Rev. D 112, 083515 (2025), arXiv:2503.14738 [astro- ph.CO]

work page Pith review arXiv 2025
[5]

Elbers et al

W. Elbers et al. , Phys. Rev. D 112, 083513 (2025), arXiv:2503.14744 [astro-ph.CO]

work page arXiv 2025
[6]

Abe et al

K. Abe et al. (T2K, Super-Kamiokande), Phys. Rev. Lett. 134, 011801 (2025), arXiv:2405.12488 [hep-ex]

work page arXiv 2025
[7]

Abubakar et al

S. Abubakar et al. (T2K, NOvA), Nature 646, 818 (2025), arXiv:2510.19888 [hep-ex]

work page arXiv 2025
[8]

Abi et al

B. Abi et al. (DUNE), Eur. Phys. J. C 81, 322 (2021), arXiv:2008.12769 [hep-ex]

work page arXiv 2021
[9]

Abe et al

K. Abe et al. (Hyper-Kamiokande), Eur. Phys. J. C 86, 170 (2026), arXiv:2505.15019 [hep-ex]

work page arXiv 2026
[10]

Abbasi et al

R. Abbasi et al. (IceCube), (2025), arXiv:2509.13066 [hep-ex]

work page arXiv 2025
[11]

Aiello et al

S. Aiello et al. (KM3NeT), Eur. Phys. J. C 82, 26 (2022), arXiv:2103.09885 [hep-ex]

work page arXiv 2022
[12]

Potential to identify the neutrino mass ordering with reactor antineutrinos in juno

A. Abusleme et al. (JUNO), Chin. Phys. C 49, 033104 (2025), arXiv:2405.18008 [hep-ex]

work page arXiv 2025
[13]

Scholberg, Journal of Physics G: Nuclear and Particle Physics 45, 014002 (2017)

K. Scholberg, Journal of Physics G: Nuclear and Particle Physics 45, 014002 (2017)

2017
[14]

Ikeda et al

M. Ikeda et al. (Super-Kamiokande), Astro- phys. J. 669, 519 (2007), arXiv:0706.2283 [astro-ph]

work page arXiv 2007
[15]

Simpson et al

C. Simpson et al. (Super-Kamiokande), The Astrophysical Journal 885, 133 (2019)

2019
[16]

Abe et al

K. Abe et al. (Super-Kamiokande), Astropart. Phys. 81, 39 (2016), arXiv:1601.04778 [astro- ph.HE]

work page arXiv 2016
[17]

M. G. Aartsen et al. , arXiv e-prints , arXiv:1510.05227 (2015), arXiv:1510.05227 [astro-ph.HE]

work page arXiv 2015
[18]

Wikipedia contributors, Ternary plot — Wikipedia, The Free Encyclopedia (2026)

2026
[19]

G. L. Fogli, E. Lisi, and G. Scioscia, Phys. Rev. D 52, 5334 (1995), arXiv:hep-ph/9506350

work page arXiv 1995
[20]

Abbasi et al

R. Abbasi et al. (IceCube), Eur. Phys. J. C 82, 1031 (2022), arXiv:2011.03561 [hep-ex]

work page arXiv 2022
[21]

Quigg, in 35th SLAC Summer Institute on Particle Physics: Dark matter: From the cos- mos to the Laboratory (2008) arXiv:0802.0013 [hep-ph]

C. Quigg, in 35th SLAC Summer Institute on Particle Physics: Dark matter: From the cos- mos to the Laboratory (2008) arXiv:0802.0013 [hep-ph]

work page arXiv 2008
[22]

G. L. Fogli, E. Lisi, and D. Montanino, Phys. Rev. D 54, 2048 (1996), arXiv:hep- ph/9605273

work page arXiv 2048
[23]

Fogli, E

G. Fogli, E. Lisi, A. Marrone, and I. Tamborra, JCAP 10, 002, arXiv:0907.5115 [hep-ph]

work page arXiv
[24]

Capanema, Y

A. Capanema, Y. Porto, and M. M. Saez, Phys. Rev. D 111, 103021 (2025), arXiv:2403.14762 [hep-ph]

work page arXiv 2025
[25]

Capanema, Y

A. Capanema, Y. Porto, and M. M. Saez (2025) arXiv:2512.17447 [hep-ph]

work page arXiv 2025
[26]

Hirata et al

K. Hirata et al. (Kamiokande-II), Phys. Rev. Lett. 58, 1490 (1987)

1987
[27]

K. S. Hirata et al., Phys. Rev. D38, 448 (1988)

1988
[28]

R. M. Bionta et al. , Phys. Rev. Lett. 58, 1494 (1987)

1987
[29]

E. N. Alexeyev, L. N. Alexeyeva, I. V. Krivosheina, and V. I. Volchenko, Physics Let- 15 ters B 205, 209 (1988)

1988
[30]

2013, The Astrophysical Journal Supplement Series, 205, 2, doi: 10.1088/0067-0049/205/1/2

K. Nakazato et al. , The Astrophysical Journal 205, 10.1088/0067-0049/205/1/2 (2013)

work page doi:10.1088/0067-0049/205/1/2 2013
[31]

Supernova Neutrinos: Production, Oscillations and Detection

A. Mirizzi, I. Tamborra, H.-T. Janka, N. Sa- viano, K. Scholberg, R. Bollig, L. Hudepohl, and S. Chakraborty, Riv. Nuovo Cim. 39, 1 (2016), arXiv:1508.00785 [astro-ph.HE]

work page Pith review arXiv 2016
[32]

High-resolution supernova neutrino spectra represented by a simple fit

I. Tamborra, B. Muller, L. Hudepohl, H.-T. Janka, and G. Raffelt, Phys. Rev. D86, 125031 (2012), arXiv:1211.3920 [astro-ph.SR]

work page Pith review arXiv 2012
[33]

Johns, S

L. Johns, S. Richers, and M.-R. Wu, Ann. Rev. Nucl. Part. Sci. 75, 399 (2025), arXiv:2503.05959 [astro-ph.HE]

work page arXiv 2025
[34]

A. L. Baxter et al. , Astrophys. J. 925, 107 (2022)

2022
[35]

Albert et al

J. Albert et al. , ”snowglobes 1.2”, https:// github.com/SNOwGLoBES/snowglobes (2021)

2021
[36]

H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi, Progress of Theoretical Physics 100, 1013 (1998), arXiv:nucl-th/9806095 [nucl- th]

work page arXiv 1998
[37]

H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi, Nuclear Physics A 637, 435 (1998), arXiv:nucl-th/9805035 [nucl-th]

work page arXiv 1998
[38]

Scholberg, Annual Review of Nuclear and Particle Science 62, 81 (2012)

K. Scholberg, Annual Review of Nuclear and Particle Science 62, 81 (2012)

2012
[39]

Shinokiet al.(Super-Kamiokande), Phys

M. Shinoki et al. (Super-Kamiokande), Phys. Rev. D 107, 092009 (2023), arXiv:2212.10801 [hep-ex]

work page arXiv 2023
[40]

Abe et al

K. Abe et al. (Super-Kamiokande), Nucl. Instrum. Meth. A 1027, 166248 (2022), arXiv:2109.00360 [physics.ins-det]

work page arXiv 2022
[41]

Abe et al

K. Abe et al. (Super-Kamiokande), Nucl. Instrum. Meth. A 1065, 169480 (2024), arXiv:2403.07796 [physics.ins-det]

work page arXiv 2024
[42]

M. L. Warren, S. M. Couch, E. P. O’Connor, and V. Morozova, The Astrophysical Journal 898, 139 (2020)

2020
[43]

S. Zha, E. P. O’Connor, and A. d. S. Schneider, The Astrophysical Journal 911, 74 (2021)

2021
[44]

Das and M

S. Das and M. H. Reno, JCAP 02, 076, arXiv:2508.16510 [hep-ph]. Appendix A: Estimated Uncertainties in Ternary Space We assume only uncertainty due to statistics of observed counts at 10 kpc and therefore that the fractional uncertainty of the inferred flux scales as √ N. In order to paint an estimated er- ror region in ternary space for ϕest, the error f...

work page arXiv 2013