Neutrino Mass Ordering (NMO)

• Neutrino mass ordering (NMO) is the arrangement of three neutrino mass states—normal (m₁ < m₂ < m₃) or inverted (m₃ < m₁ < m₂)—that underpins key insights in particle physics and cosmology.
• Experimental strategies using atmospheric, long-baseline, and reactor neutrinos exploit oscillation signatures and matter effects, with upgrades aiming for 3σ to 5σ sensitivity in the coming years.
• Cosmological observations and constraints on non-standard interactions complement oscillation data, providing additional pathways via CMB measurements, neutrinoless double beta decay, and supernova neutrino detection.

The neutrino mass ordering (NMO), also called the neutrino mass hierarchy, refers to the global pattern of the three neutrino mass eigenstates $m_1, m_2, m_3$. Determining whether $m_3$ is heavier (Normal Ordering, NO: $m_1 < m_2 < m_3$) or lighter (Inverted Ordering, IO: $m_3 < m_1 < m_2$) than the other two is a fundamental open problem in neutrino physics. Resolution of the NMO has broad implications for the flavor structure of the Standard Model, lepton-number violating processes, the interpretation of neutrinoless double beta decay results, cosmological constraints on absolute neutrino masses, and the modeling of core-collapse supernovae.

1. Theoretical Framework and Oscillation Signatures

Neutrino oscillation phenomena are governed by the propagation Hamiltonian in the flavor basis,

$$H = U\, \text{diag}(0,\, \Delta m^2_{21}/2E,\, \Delta m^2_{31}/2E)\, U^\dagger + \text{diag}(V_{CC},\,0,\,0),$$

where $U$ is the PMNS mixing matrix, $\Delta m^2_{ij} = m_i^2 - m_j^2$, and $V_{CC} = \sqrt{2} G_F n_e$ parametrizes coherent forward scattering on electrons. The ordering appears through the sign of $\Delta m^2_{31}$, which controls the resonance behavior of $\nu_e$ in matter (the MSW effect). For atmospheric and long-baseline accelerator neutrinos in the $\gtrsim$ few-GeV range, matter-induced alterations of the oscillation pattern depend sensitively on whether the mass ordering is normal ($\Delta m^2_{31} > 0$) or inverted ($\Delta m^2_{31} < 0$) (Aartsen et al., 2019).

Under constant-density, the two-flavor matter-modified mixing angle is given by

$$\sin^2 2\theta_{13}^m = \frac{\sin^2 2\theta_{13}}{[\cos 2\theta_{13} - A/\Delta m^2_{31}]^2 + \sin^2 2\theta_{13}},$$

with $A=2EV_{CC}$. Resonant enhancement of $\nu_\mu \to \nu_e$ occurs for $A \Delta m^2_{31} > 0$ (i.e., for $\nu$ in NO and for $\bar\nu$ in IO), a critical observable enabling experimental discrimination.

2. Experimental Probes of the Mass Ordering

Multiple, largely independent observables constrain the NMO.

2.1 Atmospheric Neutrino Experiments

Detectors such as IceCube DeepCore measure the energy- and zenith-dependent rates of atmospheric neutrinos traversing the Earth. The O(10^4) events in the 5–100 GeV range contain a subtle mass-ordering signature: MSW and parametric matter effects modify flavor oscillations, producing slight but distinctive distortions in the reconstructed energy and direction distributions. Statistical extraction is performed via a binned likelihood fit over $(\log_{10} E_\nu, \cos \theta_\nu, \text{PID})$, with extensive marginalization over systematic uncertainties (total fluxes, $\nu/\bar{\nu}$ ratio, energy scale, cross sections, etc.). In three years of DeepCore data, a mild preference for NO was observed, but the dataset remained fully compatible with both orderings, with $p_\mathrm{IO}=15.3\%$ and $CL_\mathrm{s}=53.3\%$ for IO (Aartsen et al., 2019).

The limiting sensitivity of DeepCore (with current data and resolutions) is below $1\sigma$, but upgrades to lower thresholds (IceCube Upgrade, PINGU) are projected to reach $3\sigma$ discrimination of the ordering within $\sim 5$ years by targeting the resonance region with enhanced statistics and reduced event reconstruction uncertainties.

2.2 Long-Baseline Accelerator Measurements

Long-baseline experiments, such as T2K and NO$\nu$A, utilize $\nu_\mu \to \nu_e$ appearance (and $\bar\nu_\mu \to \bar\nu_e$) to exploit the matter effect at $E_\nu \sim 0.6$–$2$ GeV over hundreds of km. Oscillation probabilities depend on the relative sign of $\Delta m^2_{31}$ and matter potentials, inducing $\delta_{CP}$–ordering degeneracies: the observed appearance event rates can be matched with different combinations of ordering and CP phase. Nonetheless, global fits to T2K and NO$\nu$A data show a $\sim2.4\sigma$ hint in favor of normal ordering under the standard 3-flavor paradigm (Capozzi et al., 2020, Capozzi et al., 2019). Discrimination power is fundamentally bounded by the parameter degeneracy.

Crucially, the sensitivity in $P_{\mu e}$ (and in other channels) can, in principle, be mimicked by certain non-standard neutrino–matter interactions (NSI), so robust ordering determination requires either stringent external bounds on NSI or independent confirmation from probes not primarily relying on the MSW mechanism (Capozzi et al., 2020, Capozzi et al., 2019).

3. Complementarity of Reactor, Accelerator, and Synergistic Approaches

Medium-baseline reactor experiments such as JUNO are designed to resolve the NMO by a different observable: the interference between slow ($\Delta m^2_{21}$) and fast ($|\Delta m^2_{31}|, |\Delta m^2_{32}|$) oscillation modes in $\bar\nu_e$ disappearance over $L \sim 52.5$ km. Matter effects here are negligible; the ordering is imprinted as a phase offset in the oscillation spectrum, and reconstructed $\Delta m^2_{ee}$ is subtly shifted depending on the true hierarchy (Parke et al., 2024, Cabrera et al., 2020).

At sub-percent energy resolution, correlation (“sum rule”) between JUNO's $\Delta m^2_{ee}$ measurement and precision $\nu_\mu$ disappearance measurements from T2K/NO$\nu$A ($\Delta m^2_{\mu\mu}$) breaks the NO/IO degeneracy at $3\sigma$ within one year of JUNO operation; $5\sigma$ resolution requires additional precision, achievable in the DUNE/LBL era (Parke et al., 2024, Cabrera et al., 2020).

Combining reactor and accelerator results is particularly effective: any tension in best-fit $|\Delta m^2_{31}|$ values for IO (or NO) between the spectrally-resolved reactor disappearance and the matter-enhanced appearance/disappearance LBL channels amplifies the total test statistic significantly, providing an avenue for $\gtrsim5\sigma$ discrimination even under less-than-optimal experimental conditions (Choubey et al., 2022, Parke et al., 2024).

4. Systematics, Parameter Degeneracies, and Non-Standard Physics

The robustness of NMO extraction is closely linked to the control of experimental and theoretical systematics and to the possible presence of physics beyond the Standard Model.

• Parameter degeneracies: Leading ordering sensitivity in LBL experiments is reduced in regions where the CP phase $\delta_{CP}$ can compensate for the sign change in $\Delta m^2_{31}$. For atmospheric neutrinos and reactor experiments, the ordering dependence is independent of $\delta_{CP}$, providing crucial complementarity (Aartsen et al., 2019, Choubey et al., 2022).
• Non-Standard Interactions (NSI): The oscillation signal encoding the ordering can be erased by flavor-changing NSI (especially $\epsilon_{e\tau}$) unless strongly constrained. Explicit calculations demonstrate that the $\sim2.4\sigma$ indication for NO in T2K+NO$\nu$A is fully obfuscated if moderate $\epsilon_{e\tau}$ is allowed. Only by constraining $|\epsilon_{e\tau}|\lesssim 0.02-0.05$ (i.e., much below current bounds) can robust NMO sensitivity be guaranteed in oscillation experiments (Capozzi et al., 2020, Capozzi et al., 2019).
• Scalar NSI: Non-standard scalar-mediated interactions can generate a resonant enhancement of the mixing angle $\theta_{12}^{\rm eff}$, producing a new type of degeneracy where the IO spectrum with finite scalar NSI exactly mimics the NO spectrum without NSI, nullifying NMO sensitivity in reactor experiments such as JUNO for critical values $|\eta_{ee}| \gtrsim 5.7\times 10^{-3}$ (for $m_l=0.01$ eV). Resolving the ordering thus requires tight constraints on scalar NSI from other observables and/or cross-channel synergy (Choubey et al., 5 Feb 2026, Devi et al., 2024).

5. Cosmology, Absolute Masses, and Beyond-Standard Probes

Cosmological observables are sensitive to the sum of the neutrino masses $\Sigma = m_1 + m_2 + m_3$. The minimum possible value differs between orderings: $\Sigma\geq 0.059$ eV (NO) and $\Sigma\geq 0.101$ eV (IO). Current analyses yield $2:1$–$3:2$ posterior odds for NO over IO from CMB and large-scale structure data; achieving $>95\%$ exclusion of IO requires errors on $\Sigma$ below $0.02$ eV, within reach of next-generation surveys (Hannestad et al., 2016, Ge et al., 2024).

Recently, cosmic gravitational focusing (CGF) effects, whose amplitude scales as $m_i^4$, have been proposed as a highly sensitive cosmological probe. CGF, when combined with galaxy clustering, enhances the discrimination between normal and inverted ordering, with projected $P_\mathrm{NO} \approx 98.2\%$ versus $P_\mathrm{IO} \approx 1.8\%$, improving upon the clustering-only case (Ge et al., 2024).

Other beyond-oscillation probes under active study include:

• Neutrinoless double beta decay ($0\nu\beta\beta$): Effective Majorana mass $m_{\beta\beta}$ is hierarchy-dependent, and precise $|m_{\beta\beta}|$ limits, correlated with oscillation and cosmological data, can in principle exclude IO or provide unique access to Majorana phases in the NO case (Ge et al., 2019, Ge et al., 2024).
• Direct beta-decay endpoint measurements: Experiments like KATRIN or Project 8 can probe $m_\beta^2 = \sum |U_{ei}|^2 m_i^2$, allowing comparison with NMO-driven cosmological and oscillation expectations (Ge et al., 2024).
• Supernova neutrino detection: Both energy/angle spectral analysis in water Cherenkov detectors and early-time event timing across multiple channels or detectors enable statistically robust mass ordering determination, especially when exploiting differences in MSW flavor conversion during core collapse (Jesús-Valls, 2022, Brdar et al., 2022).

6. Global Statistical Status and Bayesian Considerations

Combined likelihood analyses from oscillation, cosmological, and beta-decay data consistently yield a moderate preference for normal ordering. However, the statistical significance is highly sensitive to prior choices and parameterization. Ordering-agnostic priors (e.g., flat in $m_\mathrm{lightest}$ and $\Delta m^2$) yield a final NO preference of $2.7-3.7\sigma$, with truly agnostic setups converging to the lower end of this range (Gariazzo et al., 2022, Gariazzo et al., 2018, Salas et al., 2018).

The best fit is currently driven by oscillation data (not cosmology or $0\nu\beta\beta$), implying that conclusive $>5\sigma$ resolution must await synergistic results from the full suite of next-generation oscillation facilities (JUNO, DUNE, Hyper-Kamiokande, IceCube Upgrade, ORCA) and continued tightening of external constraints on NSI and absolute masses.

7. Outlook: Future Facilities, Synergy, and Open Challenges

Prospects for definitive NMO resolution rest on:

• Exquisite sub-percent measurements of $\Delta m^2_{ee}$ (JUNO) and $\Delta m^2_{\mu\mu}$ (DUNE, T2HK), with synergy between vacuum and matter-effect–driven probes providing cross-validation (Parke et al., 2024, Choubey et al., 2022, Cabrera et al., 2020).
• High-statistics atmospheric neutrino arrays (PINGU, ORCA) achieving $>3\sigma$ sensitivity within $3$–$5$ years (Wren, 2016, Olavarrieta et al., 2024).
• Disentangling non-standard physics, notably flavor-changing or scalar NSI, with dedicated constraints and combining information-overlapping but orthogonally sensitive experiments (Capozzi et al., 2020, Choubey et al., 5 Feb 2026, Devi et al., 2024).
• Cosmological datasets (DESI, Euclid, CMB-S4) closing the IO window via either improved clustering/sum-of-mass or CGF observations (Hannestad et al., 2016, Ge et al., 2024).
• Complementary channels: supernova neutrino bursts, collider searches for heavy neutral lepton flavor structure, and possible discovery of $0\nu\beta\beta$ at the tens-of-meV level.

The unambiguous determination of the neutrino mass ordering constitutes a critical milestone for neutrino physics, with broad phenomenological, cosmological, and model-building ramifications. The current global program, combining oscillation vacuum and matter channels, precision experiment, and cosmology, is poised to deliver $5\sigma$ resolution in the coming decade, provided that systematic and non-standard effects are comprehensively controlled and cross-checked across all relevant observables (Aartsen et al., 2019, Parke et al., 2024, Capozzi et al., 2020, Cabrera et al., 2020, Ge et al., 2024, Choubey et al., 5 Feb 2026, Devi et al., 2024).