
2
design report [31]. In this article, the physics potential
of such an analysis is discussed.
C. Neutrino flux predictions
The explosion mechanism of CCSN is still poorly un-
derstood [32]. However, the gradual increase in available
computational power has allowed the set of simplifying
assumptions about said mechanism to be reduced over
time, and in recent years CCSN models have begun to
achieve realistic self-triggered explosions [32]. Further-
more, a general concordance has emerged among the neu-
trino flux predictions from different research teams [33–
37]. New data from a future CCSN would be game-
changing to better understand the dynamics of the ex-
plosion. Hyper-Kamiokande will have great sensitivity to
discriminate between different explosion models [5, 38]
and might be complemented by studies from JUNO [39]
and DUNE [30].
D. Flavor transformations
During the neutronization burst or shock period, i.e.
the first ≲50 ms [12], the matter potential is anticipated
to be dominant over the neutrino-neutrino potential.
Consequently, flavor transformations can be described by
the standard Mikheyev-Smirnov-Wolfenstein (MSW) ef-
fect [5, 12, 30, 40, 41]. Since sin2θ13 >10−3[42], a non-
oscillatory adiabatic flavor conversion is expected [43],
sensitive to the neutrino MO. At later times, 50 ≲t≲
200 ms [12], during the so-called accretion phase and
along the cooling phases, that describe the remainder of
the burst, non-trivial effects such as SASI (standing ac-
cretion shock instability), turbulence, and neutrino self-
interactions might change significantly the flavor compo-
sition of the flux [32, 44–47].
E. Relevant interaction cross-sections
Supernova neutrino energies are typically of the order
of a few tens of MeV [20]. At such energies, the main in-
teractions with detector targets consist of neutrino- and
antineutrino-electron elastic scattering (eES) [48], elec-
tron antineutrino inverse beta decay (IBD) [49] with un-
bound protons such as Hydrogen in water, and neutrino-
nucleon charged- and neutral-current interactions with
bound nucleons (e.g. νe-CC 16O [50, 51]). Due to nu-
clear effects, neutrino interactions with bound nucleons
are poorly understood [52, 53], instead, eES and IBD
predictions are well known.
II. METHODOLOGY
A. Flux models
To study different flux models, the open-sourced
software package SNEWPY [54, 55] is used. The
models under consideration are, following SNEWPY’s
nomenclature, Bollig 2016 (27 M⊙) [56], Fornax
2021 (20 M⊙) [32], Kuroda 2020 (9.6M⊙) [57],
Nakazato 2013 (20 M⊙) [58], OConnor
2015 (40 M⊙) [59], Sukhbold 2015 (9.6M⊙) [60],
Warren 2020 (13 M⊙) [61] and Zha 2021 (16 M⊙) [62].
This set of models aims to reflect the variability in the
existing neutrino predictions among different models,
computational approaches and progenitor masses. To
account for flavor transformations, the neutrino flux
predictions are modified according to AdiabaticMSW
transformationsiusing SNEWPY. These transformations
depend on θ13 and θ23 [11], with values chosen from
the Particle Data Group [63]. Since the uncertainty
on these parameters is small the variations inflicted to
the expected flavor predictions are minor, especially
when compared to CCSN model-to-model variations.
Henceforth, the uncertainty on these parameters is
neglected, such as in Ref. [5].
Neutrino flux models, as seen in Fig. 1, have the neutrino
luminosity divided into four flavor categories: Lνe,L¯νe,
Lνxand L¯νx, where x≡µ+τ. The time evolution of
the supernova explosion is markedly different between
models such that L(t) is not a model-robust observable.
This is also true for the total neutrino luminosity
integrated over time (RLν(t)dt) due to, among other
effects, scale differences, such as the progenitor mass.
However, as presented in Fig. 1, the time-integrated
fraction for each flavor is consistently different as a
function of the true neutrino MO across flux models.
Since each model uses a different time reference defini-
tion, small time offsets are applied by setting the t0of
each model as the time at which the neutrino luminosity
reaches its maximum. Notably, this time frame could
also be set for data by analyzing the time spectrum of
the events. Although this calculation would contribute
to the detector systematic uncertainty its role is later
neglected as it is arguably a sub-leading correctionii.
B. Observable definition
To be sensitive to the MO, it is necessary to define an
observable that changes with the flux flavor content such
iImplementation details are available in Appendix A of Ref. [54].
ii If a low number of interactions is recorded, statistical errors dom-
inate. Else, σt0could be determined precisely, and a small time
offset would translate into a minor variation of the integrated
flavor composition, due to the small flavor gradients in the cu-
mulative distributions observed at around 50 ms in Fig. 1.