Snowmass white paper: beyond the standard model effects on neutrino flavor

Neutrinos are one of the most promising messengers for signals of new physics Beyond the Standard Model (BSM). On the theoretical side, their elusive nature, combined with their unknown mass mechanism, seems to indicate that the neutrino sector is indeed opening a window to new physics. On the experimental side, several long-standing anomalies have been reported in the past decades, providing a strong motivation to thoroughly test the standard three-neutrino oscillation paradigm. In this Snowmass21 white paper, we explore the potential of current and future neutrino experiments to explore BSM effects on neutrino flavor during the next decade.


Contents
Executive Summary Neutrinos are one of the most promising messengers for signals of new physics Beyond the Standard Model (BSM). On the theoretical side, their elusive nature, combined with their unknown mass mechanism, seems to indicate that the neutrino sector is indeed opening a window to new physics. On the experimental side, several long-standing anomalies have been reported in the past decades, providing a strong motivation to thoroughly test the standard three-neutrino oscillation paradigm. This can be done in three main ways. First, neutrino oscillation experiments are very precise interferometers, sensitive to subleading effects from new physics affecting neutrino flavor transitions. On a separate front, neutrino telescopes provide a unique avenue to probe BSM effects, given the very long distances traveled by the detected neutrinos (which range from the Earth radius to several gigaparsecs), as well as their ultra-high energies. Finally, astrophysical observations (such as a nearby core-collapse supernovae) in the upcoming decade will provide invaluable information and allow us to test neutrino propagation in extremely dense environments.
In the past, neutrino physics has been driven by data, from the postulation of the neutrino by Pauli to the discovery of neutrino oscillations, which were awarded the Nobel Prize in 2015. While the upcoming generation of oscillation experiments aims to measure the leptonic CP phase and the neutrino mass ordering, it will also test the standard picture with an unprecedented level of precision. For the community to succeed in this goal, a rich experimental program is required, and should be complemented with a similar effort by the community working in phenomenology and theory. We highlight the following main issues, which can largely impact our ability to constrain BSM scenarios in the next two decades: Neutrino oscillations should be tested in different environments. This is extremely useful to break degeneracies between standard and BSM parameters, and to make sure that the three-neutrino paradigm is robust. This entails contrasting the oscillation pattern in matter and vacuum, as well as on different oscillations channels, and/or for experiments relying on different detection mechanisms, among others. In fact, since the last Snowmass process new data from T2K and NOvA have become available. This is already enabling us to test the consistency of the standard neutrino oscillation picture. In the future, DUNE, JUNO, and T2HK will continue this endeavor, with a much higher level of precision. This also extends to non-neutrino experiments: as an example, dark matter direct detection experiments are sensitive to BSM effects on the neutrino floor.
Global fits between different data sets should be encouraged. Global fits to neutrino data are very powerful to unveil subleading effects in the oscillation probabilities, and to test the robustness of the three-neutrino framework. In fact, these provided evidence for a non-zero θ 13 before it was determined experimentally. Similarly, global fits will be critical to unveil new physics effects on neutrino flavor. Experimental collaborations are already working on joint fits for the determination of the neutrino mixing parameters (between T2K and NOvA, ongoing), and analogous efforts have been carried out by reactor experiments searching for eV-scale sterile neutrinos. These should be extended to constrain other new physics models as well.
Accurate simulations, and control of systematic uncertainties, is required. New-generation neutrino oscillation facilities feature beams of unprecedented luminosity which, when observed by their near detectors, will provide a very powerful tool to explore new physics effects (as long as these do not require long baselines to develop). Given the very high statistics, the bottleneck for the sensitivity to these searches is always the level of understanding of the signal sample. That is, systematic uncertainties are critical and must be evaluated and modelled thoroughly in order to derive reliable constraints. Similar restrictions apply to BSM searches using neutrino telescopes, which are subject to large uncertainties stemming from our poor knowledge of neutrino fluxes and cross sections at such high energies. At neutrino telescopes, improving neutrino flavor identification will be key in order to improve their capabilities to test for BSM effects.
The use of general theoretical frameworks and parametrizations is desirable. The use of general parametrizations allows to easily recast the obtained bounds into specific BSM models. An example is the use of effective operators to parametrize possible effects from new interactions at low energies (a similar approach is also employed in other areas of particle physics). Other examples include the parametrizations used to describe deviations from unitarity, or SME coefficients in scenarios where Lorentz symmetry is violated. This has to be complemented with an effort on the theory side, to make sure that the bounds on effective parameters are correctly interpreted, and matched onto viable models.
Experiments should try to unlock their full potential, should try everything within reach. So far, the new physics seems to evade our searches in the energy frontier, and thus it could yield unexpected signals. Neutrino experiments covered in this white paper span over 20 orders of magnitude in energy, from meV to EeV, and all of them should be encouraged to look for new signals from BSM. Because of this, experiments should be flexible and not single-purpose, and those experiments still in their designing stages should take into account all possibilities. For example, at long-baseline oscillation experiments, strong emphasis is placed on the charged-current measurements for the standard oscillation program; however, for BSM searches neutral-current measurements are equally important.
Data should be handled with great care not to miss any signals. Experimental data may hide the next discovery in unknown ways. For example, at long-baseline oscillation experiments the neutrino flux measurements at the near detectors are assumed to be unoscillated so they can be used to extract the oscillation parameters at the far detector. However, near detectors themselves are sensitive to certain BSM effects; thus, unexpected anomalies (e.g., on the event rate normalization at near detectors) should not be completely dismissed. Collaborations should also facilitate information to make it possible for the community to analyze their data to look for BSM effects.
Two notable examples in this context are the Icecube and COHERENT collaborations: their data release efforts have empowered the theoretical community to reinterpret their results and to derive new bounds on BSM scenarios, or to propose new experimental searches to extract their full potential.

Introduction
The discovery of neutrino masses stands as one of the strongest evidence pointing to the existence of physics beyond the Standard Model (BSM). Our progress in the field of neutrino physics has been mainly driven by data, which lead to the discovery of neutrino oscillations and the award of the 2015 Nobel Prize. This trend seems to continue to the future, and the momentum of both theoretical and experimental studies of neutrino is unstoppable [1].
In the near future, we aim to significantly improve our understanding of neutrino flavor transitions, since these are uniquely positioned to discover BSM effects. The basis for this assessment is three-fold: Neutrino physics is one of the most promising places to discover BSM effects. The violation of individual lepton flavor symmetries already supports the prevailing view that global symmetries are accidental. On top of this, the smallness of neutrino masses seems anomalous compared with other Standard Model (SM) fermions, and a number of BSM mechanisms have been proposed to explain it including models allowing for the violation of total lepton-number conservation. The observed mixing pattern in the neutrino sector (with large mixing angles) is also unexpected a priori, as it is drastically different from the small mixing observed in the quark sector. These peculiar features suggest neutrinos may be the key to discovering BSM physics which can explain this so-called flavor puzzle, and boost our understanding of particle physics. Persistent data anomalies in neutrino physics also imply a big discovery is just around the corner.
Neutrino oscillations are extremely precise interferometers. The tiny squared mass-splittings in the neutrino sector allows us to observe quantum interference at macroscopic scales. Similarly, any new additional couplings of neutrinos to new degrees of freedom, albeit tiny, may lead to measurable consequences and an anomalous pattern in flavor transitions. This can be seen as a correction on the dispersion relation due to new BSM effects on neutrino masses. Or, quantum states of neutrinos themselves may exhibit non-standard decoherence due to the BSM effects. New physics may also modify neutrino flavor at production, propagation, or detection processes.
Many corresponding test facilities exist or are being commissioned. The main goal of current and next-generation long-baseline neutrino oscillation experiments is to provide an unambiguous test for CP violation in the leptonic sector, and to measure the neutrino mass ordering. The current long-baseline experiments (T2K and NOvA) will continue taking data in the upcoming years, and new long-baseline facilities (T2HK, DUNE, JUNO) will start running in the next decade. These experiments will also serve as solar and atmospheric neutrino observatories, and will be prepared to detect neutrinos from core-collapse supernovae. In parallel, the observed short-baseline anomalies have motivated a strong experimental effort with many new projects running, being built, or conceived. These will not only study the properties of neutrinos, but also use neutrinos as unique probes for BSM effects. New proposals for neutrino telescopes, such as IceCube-Gen2, have the potential to test for BSM physics over a wide range of observables, from neutrino properties to multi-messenger astronomy. Overall, measurements performed in neutrino oscillations will be complemented by new, precise tests of the impact of neutrinos in cosmology and upper bounds on the effective electron neutrino mass. A strong effort dedicated to neutrinoless double-beta decay searches will also provide complementary information.
Stimulated with these motivations, this white paper explores the phenomenological implications on neutrino flavor observables which stem from four main categories of BSM effects.
The existence of additional neutrino states. Additional neutrino states, usually referred to as sterile neutrinos, have attracted great attention in the community in the past decades. Two dedicated white papers will be submitted to the Snowmass process (on eV-scale sterile neutrinos [2] and Heavy Neutral Lepton searches [3]). Here, we will address the indirect phenomenological consequences derived from these additional neutrino states, through deviations from the unitarity of the leptonic mixing matrix.
New neutrino interactions. The impact of new neutrino interactions, or non-standard interactions (NSI) is covered in this white paper in some detail since in the past few years it is a topic that has gathered a lot of attention in the community. This is motivated by a renewed interest in NSI from the theoretical point of view: on one side, on viable NSI models induced by light mediators (well below the EW scale), which could lead to potentially large effects in neutrino oscillations; and on the other, on the use of the SMEFT formalism to study low-energy NSI induced by heavy mediators.
New neutrino properties such as neutrino decay. As we are unsure about the mechanism that is behind neutrino masses, it is important to test the neutrino lifetime experimentally. Neutrino decay may also cause anomalous neutrino flavor transitions, affecting the flavor pattern for ultra-high energy neutrinos and neutrino oscillations.
Small violations of fundamental laws of physics. Finally, the exquisite sensitivity of neutrino flavor transitions allows to search for the violation of fundamental space-time symmetries, the very first principles of quantum field theory and general relativity: Lorentz and CPT symmetries, and quantum decoherence.
The phenomenological consequences and experimental sensitivities to these four types of effects depend largely on the neutrino energy. Thus, experimental technologies and facilities covered in this white paper have been categorized according to three main energy regions, outlined below.
Low-energy ( 100 MeV) neutrino experiments include solar, reactor, and supernova neutrinos, as well as new efforts to measure coherent elastic neutrino-nucleus scattering (CEνNS). We will comment on the interplay and complementarity between different data sets, and highlight potential new avenues for discovery. Large scale direct dark matter detectors may also be sensitive to new effects in the neutrino sector and provide additional sensitivity to BSM effects.
Medium-energy (100 MeV∼ 100 GeV) neutrino experiments include both long-baseline experiments and atmospheric neutrino observatories. The highest sensitivity on neutrino oscillation parameter measurements makes them a natural place to look for effects coming from new interactions or nonunitarity of the leptonic mixing matrix. Detectors range from few meter scale near detectors of long-baseline neutrino experiments, to km-scale low-energy arrays for neutrino telescopes.
High-energy ( 100 GeV) neutrino experiments are neutrino telescopes mainly designed for multimessenger astronomy, but with great potential to discover BSM effects in the neutrino sector. The downside is that flavor identification is extremely challenging. In the next decade, we may see the first cosmogenic neutrinos with unseen information of particle physics and astrophysics.
This white paper is structured in two separate sections, covering the main phenomenological (Sec. 2) and experimental (Sec. 3) aspects of neutrino flavor in current and future facilities. Each section has been further divided into the main categories outlined above. Finally, Sec. 4 presents a summary of the white paper and outlook for the future.

Theoretical aspects
A possible way to classify the plethora of new physics models can be as follows. New physics effects can stem from the inclusion of: (i) additional neutrinos, (ii) new neutrino interactions, (iii) direct consequences from the fact that neutrinos are massive, such as neutrino decay, and (iv) small violations of fundamental symmetries, such as CPT or Lorentz violation. In this section we review the main phenomenological implications for these four categories separately. Among all possible extensions of the Standard Model (SM) of particle physics to accommodate neutrino masses and mixings to address the evidence for neutrino flavor change from the observed oscillation phenomenon, the addition of extra, heavy neutrino states is arguably the simplest [4][5][6][7]. Indeed, right-handed neutrinos, in analogy to all other fermions, allow for neutrino masses but their gauge singlet sterile nature implies their mass scale could be very different and of a Majorana, leptonnumber-violating, nature. Depending on the value of this new mass scale, the phenomenology of these new states can be very different.
Given their contribution to the active neutrino masses, typically through regular Yukawa couplings, a very general consequence of these models is that the new states will mix with the active neutrino flavors so that the full neutrino mixing matrix is larger than 3 × 3. In particular the flavor states ν α = U αi ν i with α = e, µ, τ and i running not only over the known 3 light neutrino mass eigenstates but also through the extra heavy ones. As a consequence, the PMNS matrix that describes the W interactions with the charged leptons and the 3 light neutrinos, that is the 3 × 3 upper-left submatrix of U, is not unitary. We will dub this matrix N to stress its non-unitary nature. One of the possible general ways to parameterize these unitarity deviations in N is through a triangular matrix [8] with U a unitary matrix that tends to the usual PMNS matrix when the non-unitary parameters α βγ → 0 2 . The triangular matrix in this equation accounts for the non-unitarity of the 3×3 matrix for any number of extra neutrino species. This parameterization has been shown to be particularly wellsuited for oscillation searches [8,11] since, compared to other alternatives, it minimizes the departures of its unitary component U from the mixing angles that are directly measured in neutrino oscillation experiments when unitarity is assumed.

Neutrino oscillations in presence of heavy sterile neutrinos
If the new states are too heavy to be produced in the neutrino beam, the flavor states in which the neutrinos are produced and detected are given by |ν α = N * αi |ν i . That is, the truncated sum over the light accessible eigenstates of the full matrix U. Hence, the flavor basis is no longer orthonormal given that N is non-unitary. It is therefore convenient to study the evolution of the states and the oscillation phenomenon on the mass basis so that the Hamiltonian is simply given by: where V CC = √ 2G F n e and V NC = −G F n n / √ 2 are the charged-current (CC) and neutral-current (NC) matter potentials, respectively. The amplitude for a neutrino in the mass eigenstate j to interact as a neutrino of flavor β is given by the mixing matrix element N βj , which means that the oscillation probability will be given by We can remark two important differences with the standard case when N is non-unitary. Firstly, in the limit L → 0, P αβ = δ αβ . This zero distance effect is a direct consequence of the flavor eigenstates not being orthonormal. Another interesting difference is that the NC matter potential now plays a relevant role, since N † N is no longer the identity in eq. (2). Non-standard matter effects as well as zero distance effects can thus be a way to probe for unitarity deviations in neutrino oscillations.

Neutrino oscillations in presence of light sterile neutrinos
Conversely, if the extra neutrinos are light enough to be produced in the beam, these new mass eigenstates will introduce new oscillation frequencies and will participate in the neutrino oscillation phenomenon. Furthermore, oscillations will be sensitive to the whole mixing matrix U and not just to the 3 × 3 truncated matrix N . This scenario has been studied in great detail in connection to the standing short-baseline neutrino oscillation anomalies and is covered in a dedicated white paper.
Here we will instead consider the simplified scenario in which all the new oscillation frequencies introduced by the additional states are large such that ∆m 2 iJ /2E L −1 , V CC , V NC . In this limit, the new oscillation frequencies are averaged out and the oscillation probability becomes the same as eq. (3) at leading order in the mixing between the new mass eigenstates and the active flavors [11]. Thus, "high-energy non-unitarity" and "averaged-out light steriles" have exactly the same impact in neutrino flavor oscillations, up to quartic terms in the heavy-active mixing 3 or, equivalently, O(α 2 ). When the propagation is through matter, differences quadratic in the heavy-active mixing are also present, although further suppressed by the matter potential over the new mass splittings. See Ref. [13] for a thorough discussion.

Normalization of the oscillation probability
Notice that, even though the evolution operator in Eq. (3) is unitary, N is not and as such the "probabilities" P αβ do not add up to 1 and at L = 0 P αα = 1. Therefore, this amplitude cannot be directly interpreted as an oscillation probability of ν α transitioning to a ν β , at least not in the usual sense. In order to make the connection with a given experimental result, the relevant question is how the expected number of events is estimated. If the measured number of events associated with a charged lepton of flavor β from a neutrino source produced with a charged lepton of flavor α is compared with the SM expectations, then Eq. (3) does indeed correspond to the "oscillation probability" inferred from this data [14].
However, in most situations this is not the case and the neutrino flux and/or the detection cross section are rather estimated in a data-driven way. For instance, a common scenario is to have a near detector to calibrate the neutrino flux and cross section for the far detector measurement. In this situation, the unitarity deviations might also affect the measurement at the near detector. This will always be the case when the sterile neutrinos are too heavy to be produced or when they are heavy enough to be in the averaged-out regime also for the shorter baseline of the near detector L N D such that ∆m 2 iJ /2E L −1 N D . Thus, the oscillation probability inferred in this scenario would rather correspond to [11,15]: And some of the sensitivity to the non-unitarity parameters would cancel in the ratio. Thus, it is very important to consider the impact of the unitarity deviations not only at the detector probing the neutrino oscillations, but also at the measurements used to calibrate that expectation, since they might also be affected depending on the mass of the new sterile states as well as the corresponding baseline.
Finally, below the MeV scale, the sterile neutrinos may be produced in neutrino sources and participate, together with their SM active neutrino counterparts in the neutrino oscillation phenomenon. Depending on their mass, the beam energy and baseline, they may introduce new oscillation frequencies that can be searched for (or decohere, or average out their oscillations), leading to the effects discussed in Eq. (3).
In the intermediate mass regime, between the MeV 4 and the electroweak scales, the constraints from direct production are very strong and preclude their observation in the neutrino oscillation phenomenon. Above the electroweak scale, the constraints from the PMNS unitarity are stringent but may still allow to probe for them at very precise future facilities. These constraints have been summarized in the middle column of Table 1. For much lighter sterile neutrino masses, bounds can be derived through their impact in the neutrino oscillation phenomenon [11,24,27,[33][34][35][57][58][59][60][61][62][63]. We list in the last column of Tab. 1 the most stringent present constraints we found on the different elements for the new mass splittings ∆m 2 0.1 eV 2 would be in the averaged-out regime for the corresponding observable. This way the bounds become independent from the new mass scale as long as ∆m 2 0.1 eV 2 . Nevertheless, for particular values of the new mass scale, slightly stronger constraints might apply if they correspond to an oscillation maximum of some facility. Thus, the bounds in Tab. 1 should be regarded as conservative but applicable over a long range of masses. The 4 Or down to the keV scale for the electron flavor from β decays "flavor+electroweak" "Averaged-out oscillations" m > EW (2σ limit)   2 100 eV 2 and are taken to a global fit to flavor and electroweak precision observables from Ref. [36]. The bounds in the right column (shown at 90% CL, for 1 d.o.f.) apply for ∆m 2 0.1 eV 2 and come from neutrino oscillation searches, when the sterile neutrinos would be in the averaged-out regime for such a mass scale. The second number quoted in parenthesis for the α µe element includes the µ → eγ observable, which can in principle be evaded [28]. The bounds without a reference are indirectly obtained from constraints on the diagonal parameters via α αβ ≤ 2 √ α αα α ββ , which follows from the unitarity of the full U mixing matrix.
bound on α ee comes from solar neutrino data from Ref. [55]. The bound quoted is, conservatively, for the GS98 solar model while the AGSS09 model would lead to constraints about a factor 2 stronger. The new estimation of fluxes for reactor antineutrinos point to a solution of the reactor anomaly and the possibility to use reactor data also as constraints for α ee . These constraints would be similar to those from solar data but more dependent on the actual flux assumed [64]. The most stringent bound on α µµ comes from MINOS/MINOS+ data [65] and we quote a result from a global fit from Ref. [15]. For α τ τ the bound comes from a combination of different atmospheric neutrino data from Ref. [56]. Finally, the off-diagonal bounds are obtained indirectly from constraints on the diagonal parameters via α αβ ≤ 2 √ α αα α ββ , which is a consequence of the unitarity of the full U mixing matrix (these bounds turn out to be stronger than direct constraints of the α → β). Apart from constraining the unitarity-violating parameters α, an interesting question is what are the present bounds on the different elements of N when unitarity is not assumed. This has been explored in Refs. [33,62,63,66]. Notice however that the present bounds on the possible unitarity deviations encoded by α from Tab. 1 are of the same order or stronger than the present accuracy on the mixing angles of the unitary parametrization (see for instance Ref. [67]). Thus, when these are taken into account, the constraints on the individual matrix entries are not significantly affected by not assuming unitarity.

Neutrino oscillations in the presence of Large Extra Dimensions (LED)
The main motivation for introducing extra space-time dimensions was to alleviate the so-called hierarchy problem, i.e. the large difference between the electroweak and the GUT [68,69] or the Planck energy scales [70][71][72]. More interestingly, models with large extra dimensions (LED) can also accommodate non-zero neutrino masses [73][74][75][76][77]. Since right-handed neutrinos are singlets under the Standard Model (SM) gauge group, it is often assumed that they can propagate in the bulk while the SM particles are confined in the four-dimensional brane. In this context, neutrinos acquire a Dirac mass that is naturally small, since the neutrino four-dimensional neutrino Yukawa couplings are suppressed relative to charged-fermion Yukawa couplings by a factor proportional to the volume of the extra dimensions. This provides an alternative to the usual seesaw mechanism [73] for the generation of the neutrino masses, while avoiding the need for large energy scales.
In the four-dimensional brane, the Kaluza-Klein (KK) modes of the right-handed neutrino behave as an infinite tower of sterile neutrinos. If these are light enough, they will induce modifications of the neutrino oscillation probabilities. The LED-induced oscillation probability is shown in Fig. 1 left, for the case of muon neutrino disappearance.  [78]. Right: Future sensitivities to LED in neutrino oscillation experiments. DUNE sensitivity line taken from Ref. [79] and JUNO+TAO sensitivity line taken from Ref. [80], in comparison to the MINOS sensitivity from Ref. [81]. All lines correspond to 90% C.L. contours for normal mass ordering. The excluded region is towards the right of the lines. Finally, note that the actual bound derived from MINOS data is stronger than the expected sensitivity shown in this figure, see text for details.
From the probability plots in Fig. 1 left, the main features present in LED-induced oscillations are: • A global reduction of survival probabilities respect to the three-flavor case, which are typically noticeable at high energies. • Appearance of modulations and fast oscillations to Kaluza-Klein states on top of the main oscillation. So far, the only experimental collaboration that has constrained the LED parameter space is MI-NOS [81]. Since evidence for deviations from the standard three-flavor oscillations was not found, a limit on the compactification radius R < 0.45 µm at 90% of C.L. (for the conservative case m 0 → 0) was reported. An analysis including MINOS/MINOS+ data [82] could set even stronger limits. In a recent analysis of the KATRIN, Daya Bay and MINOS data [83], the bound R < 0.20 µm (R < 0.10 µm) at 90% C.L. for normal (inverted) neutrino mass ordering was reported. Before the MINOS analysis, the same LED-induced oscillations were probed at accelerator [84], reactor [84][85][86][87], and atmospheric [88] neutrino experiments, producing a similar bound.
As shown in Fig. 1 left, LED-induced oscillations imprint particular features in the event energy spectrum which can be exploited at the analysis level. With the expected energy resolution, future neutrino experiments can be sensitive to the main modulation of the fast oscillations to Kaluza-Klein states (green line in Fig. 1 left). Sensitivity to LED-induced oscillations is expected to be limited by systematical uncertainties that impact the shape of the spectrum, especially at the near detectors. This has been implemented in forecast studies, through simulations of data at different experiments like DUNE [79,89], SBND [90], and JUNO [80]. Since LED oscillations occur at short and long baselines simultaneously, a two-detector analysis has better prospects for testing this hypothesis.
In Fig. 1 (right panel), a comparison of the DUNE [79], JUNO [80], and MINOS [81] sensitivity to LED-induced oscillations is shown. For the DUNE analysis only the far detector spectral information was used to test the LED hypothesis, and 5% systematic errors were assumed in the simulation. As can be seen from the figure, DUNE will be sensitive to R 0.55µm at 90% of C.L. The same figure shows the sensitivity for a reactor neutrino experiment which uses a combination of near (TAO) and far (JUNO) detectors. Under the assumption of 1% systematic errors, the JUNO+TAO combination would be sensitive to R 0.65µm at 90% of C.L. This figure shows a synergy of future reactor and accelerator neutrino experiments, which will use different detector technologies and complementary oscillation channels. Based on these results, there are good prospects for probing LED-induced oscillations at DUNE and JUNO+TAO.
Besides neutrino oscillations [91], models with LED can be tested in beta-decay experiments [92] and tabletop experiments [93][94][95]. However, neutrino oscillation experiments provide a stronger constraint on R, which is at the sub-micrometer scale.
In the forthcoming precision Era, new physics signals might emerge as subleading effects of the three-neutrino paradigm or as a new oscillation phase(s). Since the LED model predicts a large number of sterile states, it is worth testing this model in the same footing as it is being done for the 3 + 1 scenario (for a discussion on the equivalence between a 3 + N sterile neutrino framework and neutrino oscillations in presence of LED see, e.g. Ref. [88]). This section discusses the impact of new interactions in neutrino production, detection and propagation. We start reviewing the introduction of such effects within the Effective Field Theory (EFT) formalism, for operators up to dimension 6. Section 2.2.1 discusses the inclusion of operators leading to charged-current (CC) interactions, while Sec. 2.2.2 focuses on the formalism of neutral-current (NC) Non-Standard Interactions (NSI) of neutrinos with matter fields. As we will see in Sec. 2.2.3 it is theoretically challenging to build a model that can lead to large NSI effects in the neutrino sector, while respecting bounds from charged lepton experiments. A common way out of this is to invoke the existence of light weakly coupled mediators. This possibility is entertained in Sec. 2.2.4, which includes a thorough discussion of current constraints on such scenarios. Finally, the last part of this section focuses on operators linking the active neutrinos to the dark sector: Sec. 2.2.5 discusses the phenomenological implications of new operators involving heavy neutral leptons, while Sec. 2.2.6 presents the main implications of neutrino interactions with dark matter.

General effective interactions between neutrinos: the charged-current case
From an EFT point of view, the most general Lagrangian which describes the new physics above the electroweak scale is the Standard Model Effective Field Theory (SMEFT). It has the exact same particle content and local symmetry as the SM, the difference is that higher dimensional non-renormalizable operators are added to the SM Lagrangian. In this way we can provide an effective description of the physical effects caused by heavy BSM particles, by e.g. investigating how the vertices of the SM gauge bosons and fermions are modified [96]. At lower energies, e.g. below the Z boson mass m Z , the relevant Lagrangian is the one of the weak effective field theory (WEFT), for which the heavy SM particles (Z and W bosons, Higgs and top quark) are integrated out from the theory. The parameters of WEFT can be matched to the ones of SMEFT at a scale m Z [97].
At the SM level, the interaction between neutrinos and charged leptons or quarks is described by the Fermi interaction. Focusing only on the left handed neutrinos, the Charged-Current (CC) part of the WEFT Lagrangian is: where G F is the Fermi constant, V is the CKM matrix, P L,R = 1 2 (1 ∓ γ 5 ) are the chirality projection operators, σ µν = i 2 [γ µ , γ ν ], and α are the charged leptons, while u j (d k ) stand for the up-type (downtype) quarks (j and k are quark mass eigenstates). The new interactions between neutrinos, charged leptons and quarks are parameterized in terms of the dimensionless Wilson coefficients [ jk X ] αβ , where X = L, R, S, P, T stands for the left-and right-handed, scalar, pseudo-scalar and tensor interactions, respectively.
We can systematically study these new interactions at the neutrino oscillation experiments using the approach introduced in Refs. [97,98]. From the reactor neutrino experiments Daya Bay and RENO and taking into account the EFT effect on nuclear beta and inverse beta decay interactions we can get constraints on the right handed, scalar and tensor interactions using theν e →ν e transition: [ ud R ] eµ,eτ ∼ (0.05−0.1), where X = R, S, T [97]. Using the same approach the SMEFT coefficients were studied at the FASERν detector at LHC, where new physics effects on the neutrino production through meson decays and neutrino detection through Deep Inelastic Scattering (DIS) were studied [99]. It was shown that FASERν detector has the capability to constraint new physics at the multi-TeV scale, for some 4-fermion interactions giving constraints that are much better than the existing constraints in the literature (see e.g. Fig. 2).
The resulting bounds can then be connected to the parameters within a UV-complete scenario. For example, left-right symmetric models can generate the effective operators proportional to R [100,101]. On the other hand, the scalar and pseudo-scalar interactions are usually generated in models with leptoquarks [102][103][104]. More generally, these Wilson coefficients can also be matched to the SMEFT parameters and, in this way, the low-energy neutrino experiments can be sensitive to heavy states in a specific BSM theory.

Neutral-Current Non-Standard neutrino Interactions
As discussed in the previous section, CC operators involving neutrino fields may be tested using a variety of experiments. However, effective operators including two neutrino fields (NC-like) are harder to constrain experimentally. In the case of purely vector operators, the best bounds come from oscillation experiments, which are sensitive to modifications in the matter potential that neutrinos feel as they propagate in a medium. The effect of NC-like vector interactions involving two neutrino fields is conveniently introduced using four-fermion operators of the form where f ∈ {u, d, e} and P ≡ P L , P R . Similarly to the standard weak interaction between neutrinos and matter fields, these Non-Standard Interaction (NSI) operators will induce a position-dependent    [99], see [99] for more details.
effective potential in matter: where n f (x) stands for the f fermion density in the medium at a given point x along the neutrino path, and the "1" in the first entry of the matrix refers to the SM contribution to the potential. At this point it should be stressed that, since neutrino oscillations are only determined by the differences of the eigenvalues of the full Hamiltonian, they are only sensitive to two linear combinations of the diagonal NSI parameters. Global fits to oscillation data have been performed for NSI with quarks [105][106][107] and set relatively strong bounds on the size of the off-diagonal NSI coefficients. For NSI with electrons, NSI limits from single experiments are available (e.g. for Borexino [108] or IceCube [109]) as well as from the combination of ν − e scattering experiments and collider measurements [110]. The combination of solar and reactor neutrino data has also been studied in e.g. Ref. [111], but so far no global oscillation analysis (including also atmospheric and long-baseline oscillation data) is available in the literature.
In the presence of NSI new degeneracies typically arise between the standard and non-standard parameters. In many cases these are only approximate and can be lifted (at least partially) through combinations of different data sets. The problem of new degeneracies arising in presence of NSI has been extensively studied in the literature (for an incomplete list, see Refs. ). Strategies to solve them typically involve a combination of data taken at low energies or short baselines (which are less sensitive to NSI in propagation), or from different oscillation channels (which are sensitive to a different combination of parameters). As an explicit example, it has been demonstrated that even if data are compatible with the SM, the sensitivity of future long-baseline experiments to leptonic CP violation and the neutrino mass ordering can be compromised by the existence of NSI within present bounds [114,115,134]. It is thus pressing to understand if complementary channels can probe NSI and resolve them. Fortunately, these degeneracies are not very relevant to present data: the current determination of leptonic CP violation and the mass ordering is dominated by T2K and reactor experiments, less sensitive to matter effects; and atmospheric constraints on NSI are strong enough to avoid large distortions [67,[135][136][137].
A more profound degeneracy involves the neutrino mass ordering, the solar mixing angle, and the NSI parameters. This was first observed in Refs. [138,139] where it was found that, for large enough NSI, solar data can be explained by a solar mixing angle in the second octant: the so-called "LMAdark" solution. Such degeneracy stems from an invariance at the level of the Hamiltonian describing neutrino evolution in the presence of NSI [106,107,140,141]. Although a partial breaking of this socalled "generalized mass ordering degeneracy" is a priori feasible by combining data with sufficiently different matter composition [141] (i.e., n n (x)/n e (x)), global fits to oscillation data have shown that this only leads to a mild lifting of the degenerate solution [105,142]. Consequently, additional input is required besides oscillation data.
Neutrino scattering data is also sensitive to effective operators in Eq. (6). Most importantly, while the matter potential affecting oscillations only depends on two linear combinations of the diagonal parameters, scattering experiments are sensitive to the diagonal NSI parameters individually and thus provide complementary information. Neutrino scattering experiments such as NuTeV, CHARM, and COHERENT have measured NC neutrino scattering on nuclei and are therefore sensitive to NC-NSI with up-and down-type quarks [143][144][145][146]. Of particular relevance in this context are measurements of coherent elastic neutrino-nucleus scattering (CEνNS) at very low momentum transfers (for the derived bounds see e.g. Refs. [142,[147][148][149][150][151][152]) since these are sensitive to a wider set of NSI models, as will be explained in more detail in Sec. 2.2.3. Regarding NSI with electrons, in this case the strongest bounds come from ν − e elastic scattering measurements (see e.g. Ref. [110,153] and references therein). In the future, the IsoDAR proposal [154] may be able to improve over current constraints and test NSI with electrons at the few percent level.
Besides improving the overall sensitivity to NSI parameters, the combination of oscillation and scattering data also disfavors the LMA-Dark solution. In fact, the analysis in Ref. [147] shows that the global combination of oscillation and COHERENT CsI data rejects it at more than 3σ, for NSI involving only up-or down-quarks. However, for more general NSI models allowing NSI with upand down-quarks simultaneously the dependence of the CEνNS cross section allows for a cancellation of the effective couplings to neutrons and protons, leading to a vanishing effect for CEνNS on a given nucleus (see e.g. Refs. [105,151] for a discussion). This can be partially broken with the addition of oscillation data [105,142], or from the combination of CEνNS data on different nuclei, see e.g. Refs. [146,[155][156][157]. For example, the combination of oscillations, present COHERENT data and future data obtained using a silicon target at European Spallation Source [158] will be able to discriminate between the standard LMA and LMA-Dark solutions at more than 4 σ C.L., regardless of the relative strength for the NSI couplings to up-and down-type quarks [159].
NSI effects can also be induced by new scalar mediator, leading to the effective four-Fermion interaction [160] √ where η αβ are dimensionless parameters. In a medium such as the Earth or the Sun where its constituents are non-relativistic, this term will induce an effective mass δm αβ = √ 2G F η αβ (n f + nf ), with the same Lorentz structure as a Dirac neutrino mass. 5 Thus, the correction to the dispersion relation E p + m 2 /(2E) will be through a shift as m 2 /(2E) → (m + δm) 2 /(2E) which is suppressed by the energy. This contrasts with the vector-mediated NSI, where the contribution to the dispersion relation is of form (6)) and therefore not suppressed with the energy. Given the current constraints, the effects induced by the operators in Eq. (8) will be negligible for solar, atmospheric and terrestrial neutrinos [161,162]. However, since the effect is proportional to n f + nf , in very dense environments such as the supernova core it can be important [162]. Also, unlike the case of NSI induced by a vector mediator, the scalar-mediated contributions from matter and antimatter add up in the effective mass, leading to sizable effects in the early Universe [162].

Model building aspects: heavy vs light mediators
It is worth pointing out that the interactions included in Eq. (6) are not gauge invariant. If these operators are obtained from a new theory at high energies, gauge invariance generically implies the simultaneous generation of similar operators involving charged leptons, for which tight experimental constraints exist [25,163]. Unless fined-tuned cancellations are invoked [163], this makes it hard to build a model that leads to sizable NSI effects in neutrino oscillation experiments.
A possibility to generate large NSI effects arises in radiative neutrino mass models, which explain the smallness of neutrino masses by allowing the Majorana mass terms to be induced only at the loop level. In this case, if the new mediators inside the loop have masses below the TeV scale they can induce NSI effects as well, which can have observable effects at the production, detection and oscillation of neutrinos. Reference [164] focuses on ultraviolet (UV) completion of a series of models leading to large NSI in the neutrino sector. In certain leptoquark scenarios, for example, the NSI with SM fermions are generated at tree level which leads to values as large as O(0.1). At the same time these values are consistent with constraints from Charged Lepton Flavor Violation (CLFV), direct collider searches as well as electroweak precision measurements, see Ref. [164] for details.
A second possibility to avoid the tight bounds from CLFV observables is to consider that the new physics is weakly coupled to the SM, via new mediators with masses well below the EW scale. In this regard it is relevant to note that, as a neutrino propagates in a medium, if the matter constituents are charged under the new U (1) symmetry 6 the amplitude of the new interaction will be proportional to the propagator g 2 /(q 2 + M 2 Z ) → g 2 /M 2 Z since in forward coherent scattering the momentum transfer is very small. Thus, for light enough Z masses, the new interaction may still generate large NSI-like effects even for small enough couplings to evade the bounds from charged lepton experiments. At the same time, for weakly-coupled Z models, bounds from scattering experiments with momentum transfer above the mediator mass (q 2 M 2 Z ) will be suppressed since the cross section in this case is proportional to ∼ g 2 /q 2 . For example, references [165][166][167][168][169][170] show that it is possible to build viable NSI models with αβ as large as 0.1−1 (i.e., NSI with a strength comparable to that of weak interactions) by invoking a new light Z with a mass of O(MeV) or below.
Oscillation bounds apply to NSI induced by light mediators as long as the range of the new interaction is shorter than the scale over which the matter density extends. This condition is typically fulfilled for Z masses above M Z 10 −12 eV. Below this mass range the contact interaction approximation is no longer valid, as the new force becomes long-range. For such ultra-light mediators, the impact of the new interaction can still be described in terms of a modified matter potential, using the NSI formalism. However, in this case the matter potential should be computed taking the average of the matter density within a radius r ∼ 1/M Z around the neutrino position as it propagates through a medium [161,[171][172][173].

New neutrino interactions with light mediators
Once we abandon the EFT approach, in order to write down the renormalizable interactions between neutrino fields and new light mediators, a choice has to be made regarding the nature of the new particle introduced. A possibility is to include a new light scalar φ, which is coupled to the neutrino field via a Yukawa-like interaction. In the literature, φ is sometimes referred to as a Majoron [174][175][176]. Its coupling to neutrinos can be either lepton-number conserving (φνν) or lepton-number violating (φν c ν). A second possibility is to couple the neutrinos to a light gauge boson in the form Z µν γ µ ν, which arises naturally in scenarios where the SM gauge group is enlarged with a local U (1) symmetry. A usual choice in this case is to gauge a linear combination of lepton-flavor and baryon numbers since it would be anomaly-free. In this case, the coupling will be flavor-diagonal (although it is also possible to generate also off-diagonal interactions, as shown for example in Ref. [169]) but not necessarily flavor-universal, which would lead to an effect in neutrino oscillations.
Besides neutrino oscillation experiments (see Secs. 2.2.2 and 2.2.3), the existence of new light mediators coupled to neutrinos can be probed in multiple ways.
Low-energy experiments: CEνNS and 0νββ -From a purely phenomenological perspective, an interaction of the form φν c e ν e with φ with a mass of O(MeV) or lighter can be tested by neutrinoless double beta-decay experiments [177,178]. If the mediator also couples to quarks, it can be tested by Coherent Elastic ν Nucleus Scattering (CEνNS) measurements at the COHERENT experiment [148,152,157,[179][180][181][182][183][184]. Relevant limits on light mediators have also been obtained from CEνNS searches using reactor neutrinos at CONUS [185], CONNIE [186], and the Dresden-II reactor experiment [148,149,187]. Proposals for similar searches have been put forward at the European Spallation Source [158,188] and at Los Alamos National Laboratory [189].
Fixed-target experiments -The new mediator may be produced in meson decays (e.g., , a scenario tightly constrained by KLOE, NA62 and E494 [190][191][192]. Such processes can also take place in the source of neutrino beam experiments such as DUNE, resulting in a minor but detectable distortion of the flavor composition and energy spectrum of the neutrino flux [193]. If the new mediator is produced on-shell, it may also propagate to the detector and decay back into SM particles. For a Z coupled to both quarks and leptons, this is strongly constrained by CHARM [194,195], ν-CAL [195,196], NOMAD [197,198] and PS191 [197,199]. If the Z is coupled to L e very tight constraints are obtained from the electron beam-dump experiments E137, E141 and E774 [200][201][202][203]. Additional bounds come from searches for π 0 → γZ decays at NA62 [204], and eN → eN Z at NA64 [205,206]. For compilations of bounds from fixed-target experiments on Z in the MeV range, see e.g. Refs. [195,207,208]. Neutrino trident searches -These are searches for the production of a µ + µ − pair from the scattering of a ν µ off the Coulomb field of a nucleus, known as neutrino trident production [209]. While the strongest bounds are currently obtained from CCFR [210] and CHARM-II [211], the DUNE Near Detector (ND) is expected to collect a large sample of trident events and to improve significantly over these [212,213].
Precision measurements of ν − e scattering -A light Z coupled to L e could be probed using ν − e scattering measurements, for which the cross section can be computed with a very good level of precision. Very competitive limits are obtained from low-energy solar neutrinos in Borexino [214], or from reactor neutrinos at the GEMMA [215] and TEXONO [216] experiments (for recent analyses see e.g. Refs. [217,218]). Finally, LSND and CHARM-II data on ν − e scattering would also be sensitive to light mediators, but the constraints obtained in this case are less stringent (see e.g. Ref. [218]). Future measurements of ν − e scattering at DUNE [219,220] may also be sensitive to this scenario.
Colliders -For light Z masses, the most relevant constraints come mainly from: (1) Z production in e + e− collisions in Babar, leading to both visible [221] and invisible [222] final states; (2) LHCb searches for Z bosons decaying into di-muon final states (both prompt and displaced vertex decay searches have been performed [223]).
Cosmology -Due to the limited lifetime of the neutron, Big Bang Nucleosynthesis (BBN) is greatly sensitive to the Hubble expansion rate and, therefore, to the number of relativistic degrees of freedom N ef f [224][225][226]. This severely constrains the possibility of having new light particles in thermal equilibrium with neutrinos at temperatures of a few MeV. On the other hand, Large Scale Structure (LSS) data sets bounds on the couplings to light particles with masses at or below the eVscale since these may hinder the free streaming of the neutrinos at matter-radiation equality Era, such as ν i → ν j φ [227][228][229]. Cosmological bounds can be evaded in two very different range of mediator masses. First, in the range above > 5 MeV the contribution to the neutrino energy density from the Z decay into neutrinos is sufficiently suppressed by the Boltzmann factor. Second, for M Z O(eV) and a very weak coupling (g 10 −10 ), the production rate of Z in the early Universe will be sufficiently small, efficiently suppressing the contribution to N ef f .
Supernovae -New light particles with mass below < few × 10 MeV can be abundantly produced in core-collapse supernovae, which would lead to an additional cooling mechanism, see e.g. Refs. [230][231][232]. Figure 3, shows the derived constraints on the coupling of a scalar mediator from supernova cooling considerations. Note that, when the coupling to the light scalar is too large, the new particles become trapped inside the supernova core and cannot freely transfer energy outside, so the bound vanishes. The bounds are valid as long as the new particles do not decay back to neutrinos outside the collapsing star. Additional effects caused by new light particles on supernovae include: (1) If the mean free path of produced light particles inside the supernova core is comparable to the core size (i.e., ∼ 1 km−15 km), it can alter supernova core evolution. (2) If the decay length of new particles is much larger than the core size but smaller than the size of the progenitor star, it can warm up the outer layers, affecting shock revival and leading to observable effects. (3) The new couplings provide new scattering channels for neutrinos, prolonging the neutrino diffuse time. (4) New couplings can alter the flavor composition as well as the spectrum of emitted neutrinos which can be tested by future supernova observations. In the future, the observation of the neutrino flux from a nearby supernova by Hyper-Kamiokande [233], DUNE [234] and JUNO [235] would greatly improve these constraints.
Star cooling -By requiring that the energy losses induced by the new interaction in stars do not exceed those from standard neutrino emission, strong bounds can be set on its coupling strength (see e.g. Ref. [236] for an analysis using white dwarfs). This can be used to set rather strong constraints Figure 3: The parameter range excluded by supernova cooling consideration for a scalar, S, coupled to neutrinos. In the region above the dashed lines, the new particle gets trapped inside the supernova and the obtained bounds do not apply. Figure taken from Ref. [232].
on light Z models, see e.g. Refs. [180,214,237,238]. Note, however, that these bounds typically rely on the assumption that the new mediator interacts with electrons. Fifth-force searches and equivalence principle tests -For ultra-light mediators (< 10 −13 eV) the force induced by the new interaction becomes long-range. Fifth-force searches try to observe a deviation from the standard Newton (Coulomb) potential ∝ 1/r for gravitational (electromagnetic) interactions [239][240][241][242][243]. Equivalence principle tests, on the other hand, try to search for differences in the potential felt by different materials [242]. Bounds from both types of constraints have been recasted to flavored U (1) models in Refs. [180,214,237]).
From the discussion above it becomes clear that models which are not coupled directly to L e will be constrained by fewer laboratory experiments, as in this case BaBar, and most bounds from ν − e scattering and fixed-target experiments do not apply. However, strong constraints are still obtained from COHERENT data on CEνNS as well as from neutrino oscillation data, as long as the interaction takes place with quarks. A comparison between the bounds derived from oscillations [105,180] to those from other experiments is shown in Fig. 4, using U (1) B−3Lµ as an example (see Refs. [180,214,219,237,276] for similar comparisons using other U (1) charges). The left panel shows the case for ultra-light mediators in the long-range domain, while the right panel shows the comparison in the NSI limit (for m Z > 5 MeV, where cosmological bounds on N ef f do not apply). As can be seen, oscillation data set very powerful constraints on the interaction in the both cases and considerably improves over other bounds (sometimes by several orders of magnitude). Figure taken from Ref. [180]. Results are shown for U (1) B−3Lµ as an example (see Refs. [180,214,219,237] for similar comparisons with other U (1) charges).

Effective operators involving extra neutrino states
So far before this section we have discussed new interactions involving SM neutrinos. However, we know that the SM needs to be extended with new particles to give neutrinos a mass. In this context, it is worth considering the addition of heavy sterile neutrino states to the SM particle content, since these are one of the simplest SM extensions that can serve for this purpose. In this case new operators may be generated, involving the sterile neutrino and the SM neutrino fields. In this section we will briefly discuss the phenomenological consequences of operators of dimension d = 5, 6 in this context. At d = 5, the introduction of a right-handed heavy neutral lepton (HNL) allows to write a dipole moment interaction of the form where µ ν is the transition magnetic moment, F µν is the electromagnetic field strength tensor, and σ µν ≡ i 2 [γ µ , γ ν ]. A dipole interaction of this form leads to a multitude of phenomenological consequences in direct dark matter searches, neutrino experiments and neutrino scattering data, as well as in astrophysics and cosmology (see e.g. Ref. [277] for a recent compilation of bounds). Although such operator may be probed at neutrino oscillation experiments [278][279][280][281] it does not have an impact on neutrino flavor, which is the scope of the present document. For details on the phenomenological consequences of this operator, see the Snowmass White Paper on Heavy Neutral Leptons, Ref. [3]. At d = 6, the same operators as in Sec. 2.2.2 may be generated, replacing a light neutrino by a heavy neutrino field. For example, in the case of a vector mediator. The first operator in Eq. (10) is generated, for example, from the mixing between the light and heavy neutrinos once the Z boson is integrated out. In this case, the HNL may be produced in meson decays at fixed-target or at neutrino-beam experiments. Once the HNL has been produced, it may enter a near detector and decay through the same effective interaction back into SM particles, leaving a visible signal (see e.g. Refs. [45,46,48,54] for recent studies on the production and decay of the HNL at DUNE, or Refs. [282][283][284][285] for FASERν). The second operator in Eq. (10), on the other hand, leads to interactions between the HNL and the SM fermions. As an example, operators of this form were obtained in Ref. [286] in the context of a seesaw model with an extra U (1) B symmetry. In this case, the new operators generate a new MSW potential that the new (not-so) sterile neutrinos experience while propagating in matter [287]. This scenario could also lead to novel signals in direct dark matter and solar neutrino experiments [214,286,288], thanks to their non-vanishing scattering rate with nuclei.

Neutrino interactions with dark matter
Small but nonzero neutrino mass and the nature of Dark Matter (DM) are two most pressing questions before particle physicists. A coupling between neutrinos and DM would link these two mysteries to each other. In this section, we discuss the observable consequences from such couplings at neutrino oscillation experiments, which strongly depends on the mass of the DM: PeV mass range -Decay of PeV dark matter into neutrino pairs has been suggested as the source of PeV neutrino events observed by IceCube. There are however alternative sources such as AGN [289], Tidal Disruption Events [290][291][292][293] and galactic PeVatrons [294]. Further data in future will help to reconstruct the energy spectrum and flavor ratio of events more accurately, providing a discriminant between the sources. Moreover, multimessenger methods (seeking temporal and directional correlation with the photons of different wavelength from the same source) will be very helpful to test DM as the source of ultra high energy neutrinos.
GeV-TeV mass range -DM in the GeV-TeV mass range are usually WIMP candidates. The SK collaboration has performed a dedicated analysis for the WIMP induced DM annihilation in the galactic center in the mass 1 − 10 4 GeV [295]. Also, Ref. [296] has studied the neutrino portal to DM in this mass range with the possibility that DM couples to the SM particles with either scalar or vector couplings.
MeV mass range -DM in the mass range of few MeV to tens of MeV coupled to neutrinos is well motivated DM candidate as it can be related to the neutrino mass generation mechanism [297]. For such light DM candidate the number density of DM (n DM = ρ DM /m DM ) will be sizable. The annihilation of galactic dark matter into neutrino pairs can lead to a monochromatic flux of neutrinos with energy equal to the DM mass to be detected by future experiments such as Hyper-Kamiokande [298], as well as the proposed liquid scintilator detector LENA [299]. As shown in [264], the scattering of relic supernova neutrinos off the background DM can lead to a dip in the energy spectrum of the diffuse supernova neutrino backgrounds which can be discerned by hyper-Kamiokande. Finally, the interaction of neutrinos with light dark matter in the early universe can impact the structure formation [300]. Also, by using the Borexino data one can put an upper bound on the monochromatic antineutrino flux for DM mass in the range of 2-17 MeV [300].
keV mass range -Sterile neutrino (ν s ) of a few keV mass has been widely advocated as DM candidate [301]. Its characteristic observational signature is the monochromatic γ line from ν s → γν a which can be detected by X-ray observatories focused on the galactic center(s) and galaxy clusters. See e.g. Ref. [302] for a recent review.
Ultralight Dark Matter -Thermally produced dark matter cannot be lighter than a few keV. However, if the production is not thermal the lower bound on the DM mass can be dramatically relaxed. In particular, bosonic DM can be as light as 10 −21 − 10 −20 eV. For such light dark matter mass the de Broglie wavelength exceeds the average distance between two close-by DM particles (ρ DM /m DM ) −1/3 , and the DM should be described by a classic field which oscillates in time. For example, a real non-relativistic scalar DM can be described as where ρ DM is the DM density at spacetime coordinate x = (t, x), and v is the virialized DM velocity. Coupling of neutrinos to the ultralight DM background opens up the possibility of very exciting observable effects. As shown in Refs. [303][304][305][306][307], a Yukawa interaction between real scalar DM with neutrinos induces a time varying effective mass in the dark matter background which can lead to time modulation of the parameters measured in neutrino experiments, averaged distortions of the neutrino oscillation probabilities, or even time-dependent modifications to the matter potential felt as neutrinos propagate. If φ is a complex scalar, we can define a current for it. The current-current interaction between φ and neutrinos also induces an effective mass for neutrinos but without time modulation. Such a current-current interaction can come from gauging lepton flavor symmetries, L α − L β and assigning lepton flavor to φ [308]. The Lorenz structure of the effective mass will be of form ν † ν so like the MSW effects, it can be dominant only for very high energy neutrinos. The flavor ratio extracted from the IceCube data already sets a bound on the current-current coupling but the predicted non-trivial flavor ratio can be tested by combining information from various extremelyhigh-energy neutrino telescopes (sec. 3.3.3).
Finally, neutrino interaction with DM is also advocated as a solution to cosmological tensions but recent studies show that it can only ease the so-called σ 8 tension but it cannot solve the H 0 tension [309]. As neutrinos have mass, they have radiative decay channels, e.g., ν j → ν i + γ [310,311], but their expected lifetimes are far too long to be tested [312,313]. However, if neutrinos couple to a new light or massless particle, they may decay into it, and the neutrino decay rate would be enhanced and could be probed in a variety of environments. A relatively minimal model for neutrino decay involves a light or massless Majoron, which is a spin-0 gauge singlet with non-zero lepton number, possibly related to neutrino mass generation [174][175][176]314]. Additional models include mirror models [315], SUSY models [316][317][318], left-right symmetric models [319], neutrino masses generated by a topological formulation of gravitational anomaly [320], unparticles [321,322], and others.

Neutrino decay
The decay rate of ν i , with mass m i and lifetime τ i depends on the dampening factor were t is the time elapsed since the production of ν i , L t is the distance traveled, γ i ≡ E/m i is the Lorentz boost, and E is the neutrino energy. Neutrino decay constraints are often expressed in terms of the combination τ i /m i , since the absolute neutrino mass scale is unknown, but the neutrino energy is; thus, the m i factor accounts for the Lorentz boost. Phenomenologically, neutrino decay is often classified into two main categories: invisible and visible decay 7 . Invisible decay is where the decay products are undetected because either they are sterile neutrinos or they are too low of energy to be detected in a given experiment. Visible decay involves the detection of the regenerated lower-energy neutrinos.
Neutrino decay, invisible and visible, has been probed in a wide range of experiments. In addition, the atmospheric mass ordering often has a significant impact on neutrino decay phenomenology. Simply put, in the normal ordering neutrino decay tends to lead to a deficit in muon and tau neutrinos and a possible increase in electron neutrinos for visible decay while in the inverted ordering the opposite is expected. Invisible neutrino decay can potentially enhance Earth matter effects of core-collapse supernova neutrinos besides the similarity of average energies and fluxes between electron and other flavors of neutrinos [323]. Finally, we note that if the decay basis is different from the flavor or mass bases then additional care is needed and neutrino decay may appear as unitarity violation [324].

Theoretical Formalism -Invisible and Visible Decay
The phenomenology of neutrino decay is related to the underlying physics model governing the decay and depends on whether the field is a scalar or pseudoscalar, whether it carries lepton number, and whether the neutrinos are Dirac or Majorana, making the phenomenology very rich.
Following from [325], we first discuss the Dirac case. If φ n is a scalar field with lepton number n, then the Lagrangian is modified by the terms, where Λ is the effective energy scale of the dimension-five and -six operators, v is the Higgs vacuum expectation value, y ij = y ji and h ij = h ji . One can additionally have pseudoscalar couplings with an iγ 5 term included. The decay width for mass state i is [325][326][327][328], where Additional expressions exist for the pseudoscalar case. The outgoing neutrinos are emitted with a particular spectrum depending on the nature of the decay and the decay product (helicity-flipping or helicity-conserving). For the lepton number-conserving scalar with g ij , some of the daughter neutrinos will be right-handed and thus invisible, but none will be neutrinos. For the lepton number-violating scalar, the y ij decays will all be invisible and the h ij decays will all be visible, but anti-neutrinos. If neutrinos are Majorana a similar procedure follows, with where f ij =f ij v 2 /Λ 2 and f ij = f ji and a neutrino will decay to both neutrinos and anti-neutrinos, both visible. There are other ways to generate neutrino decay from more exotic scenarios, however. E.g., in Ref. [320], neutrino masses are generated by a gravitational condensate which leads to an enhanced decay over the existing expected decay.

Current Bounds -Invisible and Visible Decay
The bounds for invisible decay that can be imposed by the neutrinos observed from SN1987A [335] are rather loose, given that only ν e were detected, and that there are significant uncertainties in the emission of neutrinos. In Section 2.3.3, we discuss the expected sensitivity to the neutrino lifetimemass ratio from the observation of a future Galactic core-collapse supernova.
Using supernova relic neutrino observations, and considering visible decay plus normal hierarchy, an upper bound for the neutrino lifetime-mass ratio of τ 2 /m 2 , τ 3 /m 3 10 10 s/eV can be extracted [336]. However, as noticed in Ref. [336], the robustness of these upper bounds should be treated with caution, given the large uncertainties embodied in the predictions of supernova relic neutrino flux. Further discussion on the theoretical treatment of the supernova relic neutrino flux, including the three-neutrino flavor transitions and decays can be found in Ref. [337].
Additionally, measurements of the CMB constrain neutrino decay, in a somewhat model-dependent fashion, to be τ i /m i 10 11 s/eV [228,[338][339][340], although these bounds may be relaxed by an additional 3-4 orders of magnitude in a separate analysis [229].
It is known that neutrino decay can produce rather different astrophysical neutrino flavor ratios compared to those obtained when only the standard oscillation mechanism is considered [341,342]. Therefore, using the IceCube flavor ratios a lower bound (τ 2 /m 2 , τ 3 /m 3 ) 10 s/eV is imposed, considering normal hierarchy [331,343]. These bounds are tied to the uncertainties in the prediction and measurement of the high-energy neutrino flavor composition. Figure 5 summarizes the invisible neutrino decay picture and table 2 shows the current experimental lower bounds on the neutrino lifetime-to-mass ratio τ 3 /m 3 using data from atmospheric experiments, long-baseline experiments, and constraints on the solar anti-neutrino flux. These bounds correspond to normal ordering, i.e., ν 3 is the decaying neutrino. Given that the experimental power for constraining τ /m is proportional to L/E, it is not surprising to see that the best limit in Table 2, for invisible decay, involves atmospheric neutrino data: τ 3 /m 3 ≥ 2.9×10 −10 s/eV. It is also clear from this table that, for the same data set (e.g., MINOS+T2K), the limits improve by one order of magnitude from invisible to visible decay. This is due to the fact that visible decay produces lower-energy active neutrinos and allows chirality-changing transitions (ν i →ν j ). Within this context, it is remarkable the very stringent limit of τ 3 /m 3 10 −5 s/eV obtained from anti-neutrino data, which represents a rather novel approach for getting it. This bound was first estimated in Ref. [344]. However, in Tab. 2 we quote the values obtained from the dedicated analysis in Ref. [345], which provide a more accurate result.
On the other hand, when ν 2 decays into invisible daughters, we find the lower bound τ 2 /m 2 ≥ 7.2×10 −4 s/eV at 95.45 % C.L. in [351], a similar value appears in [352], but given at 99 % C.L., see also the similar constraint from SNO [330]. Both limits are extracted using solar and/or reactor neutrino data. In [351] it is also found the lower bound τ 1 /m 1 ≥ 4.2 × 10 −3 s/eV, assuming ν 1 as the decaying mass neutrino eigenstate.
A number of low-to modest-significance anomalies exist which have been interpreted in the context of neutrino decay. Starting with the hints at the shortest lifetimes, they are: 1. Comparing different long-baseline accelerator experiments with sensitivity leads to a few verylow-significance hints around τ 3 /m 3 ∼ 10 −12 s/eV: 1.0(7.0) × 10 −5 † [345] † These two bounds corresponds to ν 3 →ν 2 + X and ν 3 →ν 1 + X Table 2: Current experimental bounds on τ 3 /m 3 at 90 % C.L.. Values are provided for visible/invisible decay modes and obtained from atmospheric, long-baseline neutrino data, and from KamLAND and Borexino anti-neutrino data, as indicated.
3. Comparing Planck's measurements of the CMB in TT, TE, and EE correlations, lensing measurements, and BAO measurements leads to low significance hints for invisible decay around τ /m ∼ 10 10 s/eV. Cosmological measurements prefer an invisible decay scenario with τ (2 − 16) × 10 9 s ×(m/0.05 eV) 3 for all three states at ∆χ 2 = 1.3 − 2.0 where the ranges depend on which data sets are included; any combination of data lead to some hint for neutrino decay [228]. We note that these hints are largely incompatible with each other. That is, if, e.g., the cosmological hint is verified then the accelerator and astrophysical hints would be ruled out due to the large hierarchy in lifetimes considered. This does not apply, however, if the model leading to neutrino decay does not apply in certain environments (e.g., early Universe, extragalactic space, etc.), in which case the shorter lifetime scenarios would become viable again. In addition, the data here has only been examined in the invisible case except for the astrophysical data; the scenario may be different for the visible case. Table 3 shows the expected sensitivities for τ 3 /m 3 , at 90 % C.L. and for visible and invisible decay, that will be achieved in in-development or projected atmospheric, reactor, and long-baseline neutrino experiments. In this table, the most restrictive sensitivity limit for invisible decay is τ /m ∼ 10 −10 s/eV that would be reached in future atmospheric neutrino experiments. This expected bound is similar to the one obtained with atmospheric neutrino data (see Tab. 2), while the sensitivity for invisible decay placed by DUNE will be improved in one order of magnitude with respect to the best one placed using long-baseline neutrino data. The latter observation will also be valid in the case of visible decay. These sensitivities corresponds to: * CC events and CC+ NC events, * * two different baselines, and ‡ slightly different theoretical treatments. As follows, we discuss the expected sensitivity potential for neutrino decay of astrophysical neutrinos from different sources, and that will be detected at future facilities:

Future Perspectives
• Galactic supernova neutrinos: Within the context of visible decay it has been shown in Ref. [325] that DUNE will be able to test τ i /m i 10 6 s/eV for neutrinos emitted from a corecollapse supernova located at a distance of 10 kpc. Meanwhile, Hyper-Kamiokande will probe τ i /m i 10 7 s/eV. Both sensitivities are valid for normal (ν 3 → ν 1 ) and inverted (ν 2 → ν 3 ) hierarchy (i.e. i = 3, 2).
• High-energy astrophysical neutrinos: IceCube, combining the current flavor ratios and the high-energy shower rate, will be able to set a limit of (τ 1 /m 1 , τ 2 /m 2 ) 10 s/eV in IceCube [331], for inverted hierarchy. On the other hand, the hint of neutrino decay that we already mentioned, which is supported by the fact that the cascade spectrum is softer than the track spectrum, will be probed in IceCube through the neutrino energy spectrum per flavor plus information about tau neutrinos [327,334]. It is interesting to note that IceCube-Gen2 will have the ability to discriminate between neutrino decay and standard neutrino mixing/other non-standard ones. [363].
• Cosmic neutrino background: Through the direct detection of relic neutrinos it will be possible to reach sensitivities for neutrino lifetime-to-mass ratio such as: τ i /m i ∼ 4.36 × 10 18 s/eV for neutrino masses of O(0.1 eV), invisible decays and normal hierarchy [333].

Tests of fundamental physics principles
Main authors: G. Barenboim, T. Katori, R. Lehnert Poincaré invariance and quantum mechanics belong to the most foundational principles of physics. Their synthesis is intimately associated with local Lorentz-invariant quantum field theory [364], which provides the framework for our best description of nongravitational physics, the Standard Model (SM). Moreover, versions of both of these principles form the springboard for most theoretical explorations of physics beyond the Standard Model (BSM). To identify examples of this fact in the neutrino context, one has to look no further than this document: the above sections discuss additional neutrino states, nonstandard interactions, and neutrino decay, which can all be formulated within local Lorentzinvariant quantum field theory.
The experimental nature of physics combined with the singular significance of spacetime symmetries and quantum mechanics provide ample impetus for the continued scrutiny of these principles. At the same time, certain BSM physics ideas based on these principles, and often also involving aspects of gravity, open exciting prospects for small, effective modifications to these two cornerstones at currently attainable energies. Examples include spontaneous CPT and Lorentz breaking in string theory, through noncommutative field theory, and through cosmologically varying scalars [365][366][367][368][369][370][371][372][373][374]. Likewise, small, effective departures from quantum-mechanical coherence have been predicted in the context of black-hole physics [375,376], string theory [377], and D-brane foam [378,379], etc. Various formalisms for the theoretical description of these effects in presently accessible physical systems are now available and widely used [380][381][382][383][384][385]. In light of these latter developments, investigations of CPT and Lorentz symmetry as well as quantum mechanics are not only necessary, but also acquire an element of urgency.
The remainder of this subsection provides some key ideas and possibilities for testing such fundamental physics principles showcasing the bright future for such endeavors.

CPT violation
CPT symmetry, the combination of Charge Conjugation, Parity and Time reversal, and it guarantees particle and antiparticle have the same properties such as mass, lifetime, and mixing angles. CPT symmetry is the cornerstone of our model building strategy. We describe the elementary world in terms of local relativistic quantum field theories, and CPT theorem [386] is based on only several assumptions from them, including Lorentz invariance, hermiticity of the Hamiltonian and local commutativity. Therefore the repercussions of potential violation of CPT will severely threaten the most extended tool we currently use to describe physics, i.e. local relativistic quantum fields. This is precisely the reason why testing the CPT symmetry is so important. If not found, one of the three ingredients named before must be false and all our model building strategy would need to be revisited.
A summary of the results of the search of CPT violation can be found in Figure 6. The most stringent limit is quoted in the Particle Data Group [95] in terms of relative precision, and it comes from the neutral kaon system and seems so robust that leaves little room for imagining a CPT violating world. |m(K 0 ) − m(K 0 )|/m K < 0.6 × 10 −18 .
This idea, however, may be deceptive. First, the strength of this limit arises from the choice of scale in the denominator, however, we do not have any model of CPT violation so far and therefore the choice of scale is arbitrary. We could have obtained a equally meaningful bound by choosing the Planck scale. Therefore, before having a full theory of CPT violation to address the question of which is the appropriate scale of the problem, a more meaningful bound would be; |m(K 0 ) − m(K 0 )| < 0.6 × 10 −18 m K 10 −9 eV. Second, the kaon is not an elementary particle and its mass is determined by QCD and therefore the most we could say with this bound is that QCD is CPT invariant. To test CPT violation in elementary particles, lepton probes are mandatory. Such a test was carried out with electrons; |m(e + )−m(e − )| < 8×10 −9 m e 4×10 −3 eV. Third, on top of that being the kaon a boson, the parameter which enters the Lagrangian is its mass squared and the bound should be re-written as H-H HFS H-H 1S 2S K 0 -K 0 mass p-p mass p-p q/m p-p magnetic moment e -e + mass e -e + charge ± magnetic moment e ± magnetic moment d-d mass 3 He-3 He mass there is a potential to improve it further using cosmological data [389]. Fig. 7 shows experimental results of neutrino and anti-neutrino oscillation parameter fits. In the long-baseline neutrino oscillation experiments, beams are produced either neutrino dominant or antineutrino dominant mode, and this and simulation are used to extrapolate neutrino and anti-neutrino oscillation parameters separately. T2K [390] relies on this to perform neutrino and anti-neutrino parameter extractions. Since they are consistent, the limit of CPT is usually given by the error dominated by statistics. A similar result would be expected from NOvA [391]. A mild tension between T2K and NOvA oscillation measurements can be eased by nonzero CPT-odd Lorentz violation [392], although such possibility is severely constrained by other results [393]. MINOS [394] had a magnetized far detector and unlikely other experiments it allowed a direct comparison of neutrinos and anti-neutrinos. In future oscillation experiments, statistical errors are expected to be reduced and oscillation parameters will be systematic error dominant. We expect, for example, DUNE can push the CPT test down to |∆m 2 31 − ∆m 2 31 | < 8.1 × 10 −5 eV 2 if the claimed statistics and systematic error are achieved [388]. It should also be noted, that no meaningful CP bound can be obtained before CPT symmetry is tested to the same level [395].
On top of that there are plenty of reasons to believe neutrinos might be an ideal probe for CPT violation: quantum gravity is assumed to be non-local, opening the door to a potential CPT violation. Its effects however are expected to be Planck suppressed, i.e. v 2 /M P , exactly in the right ballpark for neutrino experiments to see them. Besides, neutrinos enjoy a unique mass generation mechanism, the see-saw, and therefore their masses are sensitive to new physics and new scales. Scales where non-locality can be expected to show up.
In summary, if there is a sector where CPT violation can show up, this is in neutrino physics. Amazingly enough, this sector also offers the best bounds. CPT violation in the neutrino sector may also have important implications for the generation of the baryon asymmetry in the Universe: if CPT is violated the baryon asymmetry can be generated in equilibrium (the Sakharov conditions do not  [390,394]. Note that the MINOS far detector is magnetized to perform charge separations of events, while T2K relies on simulations of neutrino and anti-neutrino dominant beams. One cannot ignore CPT symmetry to measure the CP-violation phase in neutrino oscillation experiments. need to be fulfilled).
Note as well that the CPT violating observable chosen is also motivated by the fact that the mass squared is the parameter entering the dispersion relation E 2 = p 2 + m 2 and the natural parameter in relativistic kinematics. Therefore, deviations from this standard dispersion relation are a way to explore CPT violations triggered by Lorentz invariance breaking.

Lorentz-Invariance Violation
Lorentz symmetry guarantees that the orientation and state of uniform motion of physical systems do not affect the physical properties of the system. Small violations of this principle are accommodated in various approaches to BSM physics. For example, a nonzero vacuum expectation value of a Lorentz tensor field selects a preferred spacetime direction incompatible with Lorentz invariance. In the absence of a unique and realistic underlying theory, observable signals of such Lorentz breakdown may be described in a model-independent and general way via realistic Lagrangian effective field theory (EFT). We remark in passing that Lorentz symmetry is a key ingredient of the CPT theorem [386], and it is therefore unsurprising that such an EFT also contains corrections to CPT invariance [396]. This type of CPT violation maintains m =m, but typically generates different Lorentz-violating dispersion relations for particle and antiparticle and other CPT-violating effects.
The resulting EFT framework is known as the Standard-Model Extension (SME) [380,381,383]. The SME contains both the SM and General Relativity along with all operators for Lorentz violation. This provides a powerful, realistic, and calculable framework for the identification of Lorentz tests and the analysis of experimental data. Each Lorentz-violating term in the SME Lagrange density is composed of a Lorentz-breaking operator contracted with a coefficient for Lorentz violation that represents, e.g., the nonzero vacuum expectation value of a putative Lorentz tensor field; these coefficients control the size of the violation, and each one governs a physically distinct effect.
All SME coefficients governing the propagation and flavor oscillations of neutrinos have been systematically classified and enumerated for operators of arbitrary mass dimension [412,413] Table 4: Limits on SME coefficients obtained from neutrino experiments. Four types of tests are performed, including sidereal-variation searches in the oscillation signals (sidereal), spectrum distortions due to Lorentz-violating neutrino oscillations or mixings (spectrum), seasonal variations of Lorentz-violating solar mixings (seasonal), and searches for anomalous astrophysical neutrino flavor ratios (flavor ratio). Quoted limits are the characteristic orders of d = 3 CPT-odd SME best limits, and details are described in the references. Note that only MINOS near detector (ND) performed direct CPT tests of SME coefficients. Note also that IceCube astrophysical neutrino flavor-ratio limits depend on the production model.
Lorentz-violating SME coefficients. The upper (lower) 3×3 block governs the Lorentz-violating mixing of (anti)neutrinos. As a result of possible CPT violation, these two blocks are not identical. Nonzero entries in the two off-diagonal 3×3 pieces of H AB LIV are also possible and result from the presence of Majorana-like contributions to the SME Lagrangian.
Several unique effects originate from this EFT expression for H AB LIV . It adds distinct CPT-and Lorentz-breaking contributions to the ν a ↔ ν b and to theν a ↔ν b mixing probabilities, where a, b = e, µ, τ . A qualitatively new signal is also predicted: the mixing ν a ↔ν b of neutrinos with antineutrinos, described by the off-diagonal 3×3 pieces of H AB LIV . CPT-and Lorentz-violating contributions to oscillation-free neutrino propagation are also present. Each of the three kinds of mixing effects contains both CPT-even and CPT-odd parts. An additional feature is associated with novel energy dependences in the mixing and propagation with EFT operators of different mass dimension leading to different energy dependencies.
This large variety of interesting LV and CPTV effects implies that a correspondingly broad range of experimental parameters, such as neutrino-beam flavor composition, length, direction, and energy, as well as detector set-up and capabilities, all provide different SME sensitivities. A comprehensive experimental search for CPT and Lorentz violation in neutrinos must therefore involve an array of different, complementary neutrino experiments, each one providing access to a particular region of SME parameter space.
This assessment is borne out in past searches for CPT-and Lorentz-violating SME effects including at Daya Bay [414], Double Chooz [402,415], EXO-200 [416], IceCube [408,409], LSND [397], MiniBooNE [398,417], MINOS [401,405,418,419], SuperKamiokande [420,421], and T2K [422]. A few sample measurements and their overall sensitivity reach are listed in Table 4. For a subset of SME coefficients, sensitivities have been attained at or surpassing a level at which Planck-suppressed effects might be anticipated to manifest. Nonetheless, much of the available SME parameter space in the neutrino sector remains unexplored.
This situation, coupled with the expected abundance of incoming neutrino-flavor data in the Studies of DUNE atmospheric neutrinos can obtain first measurements (red) on some nonminimal isotropic coefficients and improved results (green) on others. Gray bars display existing limits from the IceCube atmospheric neutrino data [409]. Figure adopted from the DUNE Technical Design Report [234].
coming decade, foreshadows a surge in future CPT and Lorentz tests based on neutrino-oscillation measurements. For example, DUNE is poised to conduct various first-ever searches for SME coefficients and to improve existing SME measurements by up to six orders of magnitude utilizing beam, atmospheric, and supernova neutrinos. A representative set of projected DUNE sensitivities for atmospheric neutrino-flavor investigations can be found in the DUNE Technical Design Report [234] and is reproduced here in Fig. 8. High-energy astrophysical neutrinos observed by IceCube is similarly positioned for key advances in CPT and Lorentz tests [423][424][425][426][427][428][429][430][431][432][433]. Its detection capabilities of highest-energy atmospheric and astrophysical neutrinos translates directly to unmatched sensitivities to CPT-and Lorentz-violating operators of higher mass dimensions d because their relative effects grow with correspondingly higher powers of the neutrino energy. Although dependent on the neutrino production model, recent IceCube limits on d = 6 SME coefficients [411], for instance, are of order 10 −42 GeV −2 . This already surpasses naïve Planck sensitivity of E −2 Pl ∼ 10 −38 GeV −2 and is expected to be improved at future high-statistics neutrino telescopes, such as IceCube-Gen2 [434]. In general, the multitude of such future studies will sharpen our understanding of these spacetime symmetries, and with their projected experimental reach, discovery potential for CPT and Lorentz breakdown exists.

Quantum Decoherence
Quantum mechanics states that isolated pure states can never evolve into mixed states. An examination of a black hole's ultimate fate, on the other hand, may lead us to consider the potential of such development. Quantum fluctuations in the gravitational field due to micro black holes may lead to the loss of quantum coherence on a microscopic level [435]. Such quantum foam is inspired by quantum gravity theories including string theory [376][377][378][379]. Quantum theory should therefore be changed in some way if pure states may indeed evolve into incoherent mixtures. Fortunately it was shown that decoherence effects may be effectively implemented by introducing a phenomenological modification in the Liouville equation for the density operator of a quantum system [384]. Quantum dynamics semigroups and the master equation formalism may both be used to explain the temporal evolution of an open system [385,436] . The former is a highly generic treatment for systems with non-reversible time evolution, relying on only a few assumptions: probability conservation, rising entropy with time, and full positivity.
The decoherence term causes the damping of coherence between two mass states i and j (∼ e −γ ij ·L ). Since the effect grows with the baseline L, long-baseline neutrino oscillation experiments offer unique searches of quantum decoherence in macroscopic coherent systems. Super-K oscillatory shape measurement rejected quantum decoherence as the primary source of neutrino flavor conversion [437,438], thus, decoherence has been studied as a potential sub-dominant phenomenon. Limits on quantum decoherence parameters are studied from available neutrino data [437,[439][440][441][442]. A class of models with nonzero decoherence is especially interesting because that could be used to explain neutrino data anomalies [443][444][445][446], or be a source of CPT violation [447], or for certain models, it would give us a hint of the existence of non-zero CP-Majorana phase [448]. Neutrino coherent oscillation data are also used to test alternative quantum mechanics models [449,450].
The decoherence parameter γ ij can be a function of energy E, γ ij ≡ γ 0 ij (E/GeV ) n where n is normally treated as an integer. An experiment is sensitive to the decoherence when γ ij · L ∼ 1 so the sensitivity of an experiment with the decoherence parameter γ 0 ij can be written in the following way [451].
Although n = −1, −2, . . . is studied [452], n > 0 is predicted by some string-theory inspired models [376,378,379] and n = 0, 1, 2, . . . are more studied. This makes high-energy neutrinos with longbaseline experiments have advantages to explore the decoherence effect. This is part of a motivation for the DUNE to use a high-energy mode flux to explore new physics [132,442,453]. Atmospheric neutrinos [451,454,455] and astrophysical neutrinos [456][457][458][459][460][461] have even higher energy and longer baseline which can explore further parameter spaces. Fig. 9 shows oscillation probability differences between the standard oscillation and decoherence cases. For n = 0, the effect is accessible from neutrinos around < 20 GeV. For n = 2, interesting region are above TeV, and neutrino telescopes such IceCube, KM3NeT, and GVD have more chance to explore this case. Figure 9: Atmospheric ν µ +ν µ disappearance probability difference between the standard case and decoherence cases. Differences are shown for n = 0 (left) and n = 2 (right). Here, the x-axis is energy, and the y-axis is cosine of zenith angle which is related to the propagation distance. Here, the coherent length is set ∼ 2R ⊕ and the natural scale is chosen to be 1 TeV, corresponding to γ 0 sensitivity of 1.3 × 10 −23 GeV (n = 0) and 1.3 × 10 −29 GeV (n = 2). Figures are taken from [455].

Experimental overview and prospects
A multitude of experiments are planned to detect neutrinos and set constraints on BSM physics. Here we divide them according to the energies of the neutrinos observed in each case, since this determines the experimental techniques and challenges that have to be met.

Low energies: neutrino experiments below the GeV
Main authors: I. Esteban and A. M. Suliga Historically, sub-GeV neutrinos have been a key to precisely explore neutrino properties. And, as we demonstrate in this section, they will keep being vital. There are two major reasons for this. On the one hand, low-energy neutrino interactions are relatively well understood [462][463][464][465][466][467]. On the other hand, many natural and artificial nuclear processes produce a large flux of low-energy neutrinos. Figure 10 illustrates

Solar and reactor neutrinos
The sub-GeV neutrino flux at Earth is dominated by solar and reactor neutrinos (c.f. Fig. 10). They also come "for free", i.e., they are generated by sources not specifically designed for neutrino detection, but we can take advantage of them to carry out powerful neutrino experiments. Furthermore, the fluxes are quite well-understood [471][472][473][474][475], which allows for precise measurements of neutrino properties.
But, most importantly, solar and reactor neutrino experiments are complementary. They probe similar physics, despite having very different production mechanisms, traversing different baselines with different matter densities, and being detected with different techniques. This complementarity extends to standard 3-neutrino oscillations [67,78,476,477], but also to new physics searches: many scenarios introduce different distortions in reactor and solar neutrino oscillations [56,105,107].  Figure 11: Complementarity between solar and reactor neutrinos in measuring neutrino matter effects. We show ∆χ 2 for solar (red) and KamLAND reactor (blue) data as a function of an overall rescaling of the neutrino matter potential by a factor a MSW (in the SM, a MSW = 1). The red dashed line is obtained neglecting Earth matter effects. Figure adapted from Ref. [478]. Figure 11 illustrates the complementarity between solar and reactor neutrino experiments to determine new physics, in this case as an overall rescaling of the neutrino matter potential by a factor a MSW . As we see, solar neutrino data allows for several quasidegenerate values of a MSW , out of which reactor neutrino data picks out a MSW 1.
This research field has been extremely successful in the last decades, and it is still ongoing. As an example, the Borexino collaboration has recently detected neutrinos from the CNO reaction chain in the Sun [479]. The future is even brighter: the JUNO and JUNO-TAO experiment aims to detect reactor neutrinos with unprecedented accuracy [235,480]. Meantime, other multi-baseline reactor neutrino experiments, such as PROSPECT [481], STEREO [482], and DANSS [483] can help to disentangle different facets of possibly complex BSM effects. And, in the solar sector, next generation experiments, such as DUNE [476], Hyper-Kamiokande [484], JUNO [235], Super-Kamiokande-Gd [485], SNO+ [486], and THEIA [487], will collect more events with higher precision. Together, they will inaugurate the precision era of solar and reactor neutrino detection. From the theory and phenomenology side, it is essential that we scrutinize the opportunities that this program will bring about.

Exploring new avenues
Within the last decade, solar and reactor neutrinos detected with dedicated kiloton experiments have been the paradigm of sub-GeV neutrino physics. But this does not mean that they are the only option. There exist very promising experiments that could test any possible signal in a complementary way, or provide hints for the larger-scale experiments to scrutinize [214,278,488]. These are mostly experiments with a different main purpose, which illustrates the trend in particle physics of employing the same experiment to look for different physics.

Coherent Neutrino-Nucleus Elastic Scattering
First theorized in 1974 [489], in 2017 the COHERENT collaboration detected for the first time Coherent Neutrino Nucleus Elastic Scattering (CEνNS) using neutrinos from a neutron spallation source [490]. In this process, a neutrino interacts coherently with an entire atomic nucleus, boosting the cross-section by about 2 orders of magnitude. Although this would allow to detect neutrinos with smaller detectors (and, furthermore, the cross-section is quite well understood), the process is very challenging to detect as the only signature is a nucleus with a ∼ keV recoil energy. However, recent developments in detector technologies motivated by dark matter searches have finally made detecting this process feasible.

Dark matter detectors
Direct detection dark matter experiments are approaching the multi-ton scale. At this level, they should start being sensitive to solar neutrinos (DARWIN [527,528], DarkSide-20k [529,530]). It is therefore very appealing to explore the capabilities of these experiments that, although not mainly built to explore neutrino physics, will provide useful data. Furthermore, they will be uniquely sensitive to low-energy pp neutrinos. Figure 12 illustrates the complementarity between reactor neutrino experiments and solar neutrino data from DUNE and the DARWIN dark matter experiment (labelled as DWN). The figure shows the combined sensitivity on θ 12 and the mixing angle θ 14 parametrizing the mixing with putative sterile neutrinos.
Exploring sterile neutrinos in ν e detection has become more pressing now that the Gallium anomaly has reached the > 5σ level [531,532]. Solar neutrino data exclude an sterile neutrino interpretation of the anomaly with ∼ 3σ, and provide the only constraint in a large sterile neutrino mass range. It is therefore very appealing that a combination of different experimental approaches (a priori with As we see, different experiments probe different directions in parameter space. The best results are thus obtained combining them. Figure from Ref. [55]. different systematic uncertainties) will significantly increase the current sensitivity; and it could be interesting to explore if the complementarity also extends to alternative models that may explain the Gallium anomaly.
The findings of numerous studies indicate that BSM physics may considerably alter the evolution of core-collapse supernovae. In turn, these supernovae should serve as compelling testing facilities to discover or constrain BSM physics. This includes but is not limited to sterile neutrinos [547,[562][563][564][565][566][567][568][569][570][571][572][573][574][575], non-standard mediators coupling to neutrinos [232,262,263,576], or also charged fermions [231,577,578] and quarks [179,[579][580][581], where any flavor-changing scenarios play a role as well. But core-collapse supernovae are complex phenomena that incorporate symmetries and interplays between particle, nuclear physics as well as hydrodynamics [582]. Due to this, to make sure we get reliable limits, modeling of the source evolution has to be performed meticulously, taking into account the arising feedback from invoked BSM physics, as that might drastically modify the obtained limits, e.g., [575].
Luckily, for the next Milky Way core collapse supernova we will have much better detectors [583] that will also complement advances in astrophysics. Furthermore, we should be able to detect all neutrino flavors: Super-Kamiokande [584] (or its successor, Hyper-Kamiokande [484,585]) will lead the statistics by detectingν e , DUNE will detect ν e [234], and JUNO [235] together with large-scale direct dark matter and CEνNS detectors will observe other flavors. By looking at the excess over large background rates, IceCube and KM3NeT should also be able to detect supernova neutrinos [586,587]. It is thus interesting to explore the complementarities that detecting different flavors will bring about, and whether uncertainties related to collective neutrino flavor effects can be overcome (see, e.g., Refs. [325,[588][589][590]). In addition, this multi-channel observation will serve well in probing and constraining non-standard physics in the neutrino sector, described in the previous paragraph.
Finally, even if we are unlucky and a galactic supernova does not take place within the near future, supernovae have been exploding throughout the entire history of the Universe, generating a Diffuse Supernova Neutrino Background (DSNB) [591][592][593]. Its detection is also valuable by itself regardless of the occurrence of the next galactic supernova, as it is essential to observe neutrinos from multiple supernova events, which is by definition the DSNB, to understand the whole core-collapse supernova population [594][595][596][597][598]. It also opens up a window to probe non-standard interactions occurring between DSNB and cosmic neutrino background [599] or dark matter [264]. Although backgrounds affecting the detection of this flux are large, the enrichment of Super-Kamiokande with gadolinium should dramatically reduce them and make detection ofν e from the DSNB an actual possibility [485,[600][601][602][603][604][605][606]. In the future, Hyper-Kamiokande [484], DUNE [234] and JUNO [235] could detect bothν e and ν e ; and together with dark matter detectors put an upper bound on contributions from other flavors [607,608]. Both the generic DSNB detection as well as its flavor content would bring about interesting physics studies [336,337,607,609].

Medium energies: neutrinos with energies between the GeV and TeV scales
Main authors: D. V. Forero and M. Ross-Lonergan As we increase in neutrino energy, we enter the realm of man-made neutrino sources from accelerators, and atmospheric neutrino measurements. In this section we focus on the energy range between the GeV and TeV scales, where oscillation effects are observable, leaving for Sec. 3.3 the discussion of ultra-high-energy neutrinos and neutrino telescopes.

Past and Contemporary Neutrino Beam Experiments
At neutrino beams, high-energy proton beams are collided with a stationary target to produce charged mesons (predominantly pions and kaons). These are then focused inside a magnetic horn where their subsequent decays produce a relatively collimated and intense beam of neutrinos. By choosing the polarity of the magnetic focusing horns either positive (π + ,K + ) or negative (π − , K − ) mesons can be selectively chosen, leading to a predominately neutrino (ν µ ) or anti-neutrino (ν µ ) beam respectively. The fact that the initial proton beam, the magnetic horn and target are all carefully designed and controlled means that the resulting spectrum of neutrino species and kinematics can be well understood in comparison to many natural sources of neutrinos. This is further bolstered by the use of a neardetector, located close to the neutrino beam source, that measures the expected spectrum of neutrino well before any oscillations are expected to take place under the assumption of three neutrinos. Precise measurement of the neutrino interactions at the near detector allows the spectrum at the far detector be estimated to a much higher precision, leading to higher sensitivity to oscillatory effects that take place in the interim. This section summarizes past and contemporary neutrino beam experiments which could be sensitive to BSM effects on neutrino flavor, while a review of future proposals is provided in Sec. 3.2.2.

MINOS/MINOS+
MINOS [610] was an on-axis, long-baseline experiment with a near detector located at Fermilab, and a far detector at the Soudan Underground Laboratory in northern Minnesota. It ran originally from 2005 to 2012, and for a further 3 years as MINOS+. The neutrino energy of the NuMI beam is configurable, with the majority of the original MINOS data being taken in "low-energy" mode with a peak neutrino energy of ≈ 3 GeV, and MINOS+ taking data in the "medium-energy" configuration whose peak neutrino energy is ≈ 6 GeV. Both the far and near detectors were functionally identical magnetized steel plastic scintilator detectors [611]. The near detector had a mass of 980 ton, located 1km from the target, with the 5.4 kt far detector located 735km from the target. The primary goal of MINOS/MINOS+ was the measurement of the atmospheric oscillation parameters though both ν µ disappearance and ν e neutrino appearance channels [612,613].
Alongside the standard three-neutrino paradigm, MINOS/MINOS+ had a significant program of BSM searches that highlighted the power of long baseline experiments in probing exotic models. These results include searches for Lorentz invariance and CPT violation using both the far detector [404] and near detector [401], neutrino decoherence [612], non-standard neutrino interactions [614,615], light sterile neutrinos searches [65] and large-extra dimensions [81]. Alongside their own light sterile neutrino searches, MINOS+ also published a combined study with the Daya-Bay collaboration and Bugey-3 experiment [616].

NOvA
The NuMI Off-Axis ν e Appearance (NOvA) Experiment is a long-baseline neutrino experiment looking for oscillations in neutrinos originating in the "Neutrinos at the Main Injector" (NuMI) beam [617] at Fermilab in Illinois. The experiment consists of two detectors, one 300 ton near detector located at Fermilab and a second 14 kt far detector located 810 km away on the surface in Ash River, Minnesota. The location of the far detector is 0.8 • off the primary NuMI neutrino beam axis, leading to a predicted neutrino spectrum peaked strongly at 2 GeV energies. There are several planned searches for BSM physics at NOνA, with a mature analysis on NSI focusing on the µτ element by looking for deviations in the CC ν µ spectrum, when compared to three ν oscillations. The combination of this with a CC ν e disappearance study can be sensitive to µ and τ individually. NOνA has also published searches for 3+1 sterile neutrino oscillations in neutral current interactions both for neutrino [618] and antineutrino [619] running mode, where no evidence for oscillations was observed in either channel.

T2K
The Tokai-to-Kamioka (T2K) experiment is a 275km long-baseline experiment, using neutrinos from the J-PARC in Tokai, operating since 2009. The primary near detector is located 280m downstream, ND280, and samples the neutrino flux before oscillations. Alongside ND280 there are additional detectors such as INGRID, WAGASCI+BabyMIND used to constrain the flux and better understand interaction cross-sections. The Far detector consists of the 50kt water Cerenkov Super-Kamiokande detector, located 2.5 • off the neutrino beam axis. This results in the neutrino beam forming a narrow peak at 650MeV, at the expected maximum oscillation point assuming three neutrinos.
In the past T2K has probed Lorentz and CPT violation [399] by searching for sidereal modulations in the observed neutrino interaction rate, with plans to improve with larger data sets in the future. T2K has also published bounds on light sterile neutrinos [620] using a combination of five charged current and three neutral current samples, observing no evidence for new oscillation frequencies. Future searches will be further improved by a series of planned improvements to both the beam and near detector complex. By studying the charged pion yields in a T2K-replica target measured by the NA61/SHINE experiment [621] the beam flux prediction can be constrained and has can reduce the flux uncertainty from about 10% to around 5% near the flux peak. This is expected to be further improved as more data is obtained by NA61/SHINE experiment on the replica target. In addition, a magnet power supply upgrade as well as RF upgrades mean the neutrino beam itself is expected to increase in power from 500 kW to 1MW by Japanese Fiscal Year 2025.

Future and Potential Neutrino Beam Experiments
The overall energy and baseline ranges for both contemporary and future planned long-baseline experiments is highlighted in Fig. 13. All contemporary experiments are situated at or near the first oscillation maxima assuming the standard three neutrino paradigm. DUNE and T2HK also are sitting directly on the first oscillation maxima, but T2HKK (detector situated in Korea) and ESSνSB target the second oscillation maxima. By sitting at the second oscillation maximum in terms of L/E ν , the CP violating term in ν µ to ν e oscillations is approximately three times larger than at the first oscillation maxima.

T2HK/T2HKK
The T2HK [622] (Tokai to Hyper-K) experiment continues the usage of the of the T2K beam of neutrinos from the J-PARC in Tokai, upgraded to 2.6 times the power (1.3 MW), but with the far detector being the next generation underground water Cherenkov detector Hyper-Kamiokande located near Kamioka town. Hyper-Kamiokande build on the extensive experience and knowledge of water Cherekov reconstruction developed during the successful running of the Super-Kamiokande experiment. For more details on the Hyper-Kamiokande detector, see the following section on atmospheric neutrinos. The T2HKK ("Tokai to Hyper-K detector in Korea") [233] is a proposed expansion featuring both a detector in the Kamioka mine and a second Hyper-Kamiokande detector located in Korea. When considering the Korean candidate site, the longer baseline and thus larger matter effect significantly enhances sensitivities to non-standard interactions [134,623].

DUNE
The Deep Underground Neutrino Experiment (DUNE) is a next-generation, long-baseline neutrino oscillation experiment, designed to be sensitive to ν µ to ν e oscillation. The experiment consists of a 1.2 MW, broadband neutrino beam around 1 − 5 GeV, with a near detector complex located at Fermilab, and a > 20 kt liquid argon time-projection chamber far detector located at the 4850 ft level of Sanford Underground Research Facility in Lead, South Dakota. The combination of high powered beam and comprehensive suite of both near and far detectors provides a huge opportunity for a rich and diverse set of BSM physics searches, which is explored extensively in [79] and [234].
Due to the particularly long 1300km baseline DUNE is sensitive to non-standard interactions that affect neutrinos propagation thought the earth and can substantially improve bounds on NC matter NSI parameters [114,116,121]. New CC NSI interactions which can produce BSM effects at the source and detector have been shown to not be discoverable at DUNE [121,624]. DUNE will also be able to place world leading bounds on Lorentz and CPT violation [388,625] in the neutrino sector, improving bounds by several orders of magnitude. Searches for non-unitary at DUNE [11,35] have been shown to produce bounds comparable with other constraints from present oscillation experiments. Further studies on sensitivities to new physics at DUNE include neutrino decay, see Tab. 3.

ESSνSB
The ESS neutrino Super Beam (ESSνSB) is a proposed neutrino beam utilizing the 5 MW proton linear accelerator of the European Spallation Source (ESS) in Lund, Sweden, that is currently under construction. Although the final site for the associated far detector is not set, current proposed sites are the Garpenberg mine (baseline of 540 km) or the Zinkgruvan mine (baseline of 360km) both located in Sweden. In combination with relatively low energy expected flux peak of ≈0.25 GeV [626], this allows for the probing of the second oscillation maxima of ν µ → ν e oscillations.
Studies of the potential of the ESSνSB have been performed to bound non-unitarity [627], with the best sensitivities to the non-unitary parameters studied (α 21 ) are achieved for the shortest 200 km baseline that was studied, with 360km performing somewhat better that 540km although remaining comparable to current bounds. Studies have shown that ESSνSB would also be sensitive to CPT violation [628] and to neutrino decay (see Tab. 3). Due to the shorter baselines, the ESSνSB is less sensitive to matter effects and therefore to NSI parameters.

Precision beams: nuSTORM and Iso-DAR
While conventional superbeams (neutrino beams from meson decay-at-rest) from accelerators provide well understood fluxes in comparison to many natural sources, there remains still an overall O(5 → 10%) uncertainty, growing larger at lower energies. Alternative approaches, for example neutrino factory (neutrino beams from muon decays in a storage ring) and beta beam (neutrino beams from boosted radioactive decays) have known spectrum and flavor structures. These beams can contribute BSM physics search from neutrino flavor beyond what superbeams can.
NuSTORM [629] -The Neutrinos from Stored Muons (nuSTORM) facility forms the neutrino beam from the decay of stored muons and pions in a low-energy muon decay ring. By placing appropriate instrumentation in the decay ring, the total integrated neutrino flux can be estimated with a precision of ≈ 1%. The initial pion injection will lead to an intense flux of neutrinos from primarily pion decay, with subsequent rotations of the stored muons providing a neutrino beam that can be calculated precisely using the known three-body Michel decay parameters alongside the instrumentation of the muon decay ring.
IsoDAR [630] -An accelerator-driven single-isotope antineutrino source (IsoDAR) uses a highpower cyclotron to deliver protons on a beryllium target, and ejected neutrons are captured by a surrounding 7 Li sleeve. Then subsequent decays ( 8 Li → 8 Be + e − +ν e ) produce a pureν e beam with known spectrum. Currently, there is a ongoing effort to build such facility at the Yemilab underground facility, South Korea. Here, 60 MeV cyclotron is proposed to operate from 2027, with 2.3 kt of liquid scintillator counter.
MOMENT [631] -The MuOn-decay MEdium-baseline NeuTrino beam facility is proposed to generate a very intense beam of neutrinos and antineutrinos (from a proton beam expected to reach 15 MW power at 1.5 GeV energy) via decay of positive and negative muons to match the second oscillation maximum with the neutrino beam energy peaked at around 200-300 MeV and a baseline of 150 km, although different baselines can be considered [632]. BSM physics potential studies for nonstandard neutrino interactions [633] and invisible neutrino decays [357] has been performed assuming a large gadolinium-doped water Cherenkov detector of 500 kt.

Near detectors
New generation neutrino oscillation facilities feature beams of unprecedented luminosity which, when observed by their near detectors (ND), will provide a very powerful tool to explore extremely small new physics effects as long as they do not require long baselines to develop. These searches can be divided in three main categories: • Searches for the decay products of long-lived new particles produced in the beam. This is the case of new light and weakly interacting gauge bosons, dark fermions, or even heavy neutrinos in the MeV-GeV range as outlined in Section 2.1. These scenarios are discussed in more detail in the dedicated Snowmass White Papers in Refs. [3] and [634].
• Precision measurements of neutrino scattering processes. Due to the large cross section uncertainties involved in neutrino-nucleus scattering, better prospects are expected for ν − e elastic scattering. As discussed in Sec.2.2.4, these measurements can be used to constrain models with light mediators, as long as they are coupled to electrons. For example, in flavored U (1) models such as L e − L µ , a new mediator opens up a new t-channel diagram for ν − e scattering through the Z , which leads to interference with the existing SM processes. This has been studied for the DUNE ND, both in the context of a total rate analysis [219,635] and for a binned energy spectrum [220] which shows that in a bin-by-bin analysis the effect of destructive interference is reduced compared to the total rate approach. The DUNE ND could provide competitive bounds for mass ranges (≈ 0.02 − 0.2 GeV) of the L µ − L e Z .
• Additionally, new physics effects that may lead to neutrino flavour change at very short or even zero baselines are best probed at near detectors. These include the zero distance effects discussed for non-unitarity and LED, but also light sterile neutrinos with large enough ∆m 2 so as to be in the averaged-out regime as discussed in Sec. 2.1 as well as CC NSI affecting neutrino production and detection processes, described in Sec. 2.2. Indeed, the most stringent constraints in the right column of Tab. 1 stem from the non-observation of this zero-distance effect and in the case of α µµ the MINOS/MINOS+ ND plays a crucial role, as discussed in [15]. Future prospects for the DUNE ND to some of these BSM scenarios have been studied e.g. in Refs. [636][637][638].
Given the very high statistics that can be collected with intense beams at near detectors, the bottleneck for the sensitivity to these searches is always the level of understanding of the signal sample. That is, systematic uncertainties are critical and must be evaluated and modelled thoroughly in order to derive reliable constraints.

Neutrinos from colliders: FASERνnd SND@LHC
Two experiments at CERN aim to measure for the first time the intense flux of neutrinos created in pp collisons at the LHC. Constructed using 1000 emulsion layers interleaved with 1-mm tungsten plates, FASERν [639,640] sits 480m from the ATLAS interaction point and was installed during the LHC Long Shutdown 2. FASErν aims to measure primarily ν µ interactions but with sensitivity to ν e and ν τ interactions also. The energy of the expected ν µ neutrinos peaks ≈ 0.5 TeV, with ν τ and ν e components being slightly higher O(TeV) energies. Its primary goal is the detection of these collider neutrinos across all three flavors and measure their cross sections at O(TeV) energies, higher than any man-made neutrino beam thus far. Beyond this, it has been shown that FASERν will provide strong bounds on neutrino non-standard interactions [99] as described in Sec. 2.2.1. Complementary to FASERν, the Scattering and Neutrino Detector at the LHC (SND@LHC) [641] will search for collider neutrinos off-axis angles relative to the beamline larger than those covered by FASERν, with the the corresponding neutrinos mostly originating from charm decays with energies between 350 GeV and a few TeV, with sensitivity to all neutrino flavours. The detector will be a small-scale prototype of the scattering and neutrino detector of the SHiP experiment [642].

Atmospheric neutrinos
Neutrinos produced by the interaction of incoming cosmic rays with the upper atmosphere of the Earth provide us with an abundant source of neutrinos. These neutrinos, referred to as atmospheric neutrinos, are the byproduct of the decay of mesons that are produced by energetic protons colliding with the elements that composes the Earth atmosphere. Not only does Nature provide us with a natural neutrino source, but its flux spans several orders of magnitude in energy, ranging from few hundred MeV to hundreds of TeV. These atmospheric neutrinos were the first natural neutrinos to ever be observed, at experiments in the East Rand Proprietary Mines in South Africa [643] and the Kolar Gold Fields mine in India [644]. Atmospheric neutrinos has since been observed in many experiments including MACRO [645], Soudan2 [646], Super-Kamiokande [647] and IceCube-DeepCore [648]. Historically, the first observation of the neutrino flavor conversion, consistent with the neutrino oscillation phenomenon, was performed in Super-Kamiokande. This observation was made possible thanks to both a wide atmospheric neutrino flux and that neutrinos can travel different distances from the production point to the detector, which ranges from few tens of kilometers to the Earth diameter. In the following we describe the main experimental features of current atmospheric neutrino experiments and the BSM studies that have been performed so far to search for any modifications to the expected neutrino flavor profile. This will help to introduce the future experiments and their experimental prospects.

The Super-Kamiokande (SK) experiment
SK has been observing both solar and atmospheric neutrinos for more than two decades providing valuable information that has converged in our current understanding of the three-flavor neutrino oscillation mechanism.
The SK detector, located underground inside a zinc mine in Kamioka-Japan, is a cylindrical 50 kt water Cerenkov neutrino detector with nominal fiducial volume of 22.5 kt. SK is able to observe atmospheric neutrino events spanning four orders of magnitude. Fully-contained events have energies ranging from a few hundred MeV to ∼ 10 GeV. Partially contained events (with vertices inside the fiducial volume, but leptons that leave the inner detector) have energies ranging from a few GeV to tens of GeV. Finally, up-going muon events (those which start in the surrounding rock and pass through both the outer and inner detector volumes) have energies between a few GeVs and TeVs. The data collected since SK started operating in 1996 (divided into four different running periods) have been used to search for a wide set of BSM effects. In particular, tests of non-standard neutrino interactions (NSI) [649], tests of Lorentz invariance [650] and light sterile neutrino searches [651] have been performed.
In 2018, the Super-Kamiokande experiment entered in a new phase, after the water purification system was upgraded to allow operation while doped with gadolinium sulfate (SK-Gd). This improves the differentiation between neutrinos and antineutrinos, thanks to neutron capture by gadolinium (neutron tagging), and is expected to impact the physics reach of SK program. Operation with 0.01% loading started in 2020. In addition, the atmospheric neutrino analysis is expected to be improved by better modeling of neutrino interactions from measurements at the T2K near detector [652] and MINERvA [653].

The Hyper-Kamiokande (HK) experiment
Hyper-Kamiokande is the successor to the SK experiment with a 8.4 times larger effective volume. Hyper-Kamiokande is expecting to start operation in 2027. The detector tank will have a total volume of 258 kt. As SK it is separated into an inner and an outer veto detector. The inner region containing 217 kt of water and will have an array of 40,000 high QE Box & Line (B&L) PMTs with 50 cm diameter. The new PMT type improves the photon detection efficiency with a better timing resolution, and charge resolution of SK PMTs [654].
With a larger detector mass than SK and running for 20 years (high statistics) plus the mentioned detection improvements plus better modeling of neutrino interactions (reduction of systematical uncertainties), precise measurements with atmospheric neutrinos are expected in HK. This can be used to test for BSM scenarios that affect neutrino flavor, similar to what has been done at SK, such as Lorentz invariance violation [650], non-standard neutrino interactions [649], and light sterile neutrino oscillations [651].

IceCube DeepCore
The IceCube neutrino detector is a 1 km 3 of ice instrumented with 86 strings with 60 digital optical modules (DOMs) each, located at the South Pole, Antarctica. The basic mechanism is similar to that of Water Cerenkov detectors: the DOMs contain the PMTs that detect the Cerenkov radiation emitted by charged particles (produced in neutrino interactions) traveling faster than the speed of light in the ice. The strings are arranged in a hexagonal grid, buried below the ice surface at depths ranging from 1450 m to 2450 m, and with a typical separation between strings of 125 m (or, 50 m between DOMs).
At the center of the IceCube detector, DeepCore is comprised of eight additional strings with more efficient PMTs (increasing light collection), and with a typical separation of 17 m between strings (7 m between DOMs). With an increased DOM density, DeepCore is sensitive to lower neutrino energies than Icecube. Moreover, its geometry was optimized to maximize the detection efficiency for neutrino events in the range of 10 GeV to 100 GeV. In this energy range, the high cosmic-ray background is mitigated by using IceCube surrounding strings as an active muon veto [648]. At DeepCore, muon tracks are reconstructed with an angular resolution of 12 • at 10 GeV. The muon energy can be reconstructed from the muon track length and adding the energy of the hadronic shower (estimated from the total amount of light collected by the detector) provides a proxy for the neutrino energy [655]. The median energy resolution is about 30% at 8 GeV [656].
Using three years of atmospheric neutrino data from the DeepCore detector [657] a search for a light sterile neutrino was performed in Ref. [656], resulting in limits on the sterile-active mixing matrix elements |U µ4 | 2 < 0.11 and |U τ 4 | 2 < 0.15 (at 90% C.L), for a value of the neutrino mass splitting ∆m 2 41 = 1 eV 2 . Using the same data sample, a limit on the NC NSI parameter ε µτ was found to be in the range −0.0067 < ε µτ < 0.081 (at 90% C.L.) [655], performing a similar analysis as Super-Kamiokande [649]. Recently, a more general analysis considering all NSI couplings with quarks (including the effect of the NSI phases) was performed in Ref. [109], using three years of Icecube DeepCore data.

IceCube Upgrade
An additional array of 7 strings, with an separation between DOMs of 3 m, will be added to DeepCore within the following year. This upgrade will impact important experimental features as the energy and directional resolution, lowering the detector energy threshold while increasing the efficiency below 10 GeV. The detector and ice calibration will also be improved by the deployment of new calibration devises. A significant reduction of systematic uncertainties is expected which, together with the high IceCube statistics, will improve precision on the determination of standard oscillation parameters and at the same time extend the BSM physics sensitivity.

High energies: neutrinos above the TeV scale
Main authors: C. Argüelles and T. Katori The quest to observe neutrinos from high-energy astrophysical sources, where hadronic activity has been inferred to be present, has led to the development of a category of detectors known as neutrino telescopes. In order to achieve the large effective masses needed to observe these very small fluxes, these detectors use naturally occurring materialse.g. lakes, oceans, glaciers, mountains, etc -as targets for the neutrino interactions.
Current and proposed neutrino telescopes are listed in Table 5, where they have been organized by the detector technology and their acceptance to different neutrino flavors. These experiments can be roughly grouped into two categories: those that aim to observe and characterize high-energy neutrinos (1 TeV < E ν < 10 PeV) and those that aim to discover extremely-high-energy neutrinos (E ν > 10 PeV). The high-energy neutrino telescopes (HENT) can measure the high-energy part of the atmospheric neutrino flux and astrophysical neutrino sources. HENT will be able to perform new physics searches using both of these sources, which we discuss in Sec. 3.3.1 and Sec. 3.3.2 respectively. Extremely high-energy neutrino telescopes (EHENT) aim to observe for the first time neutrinos produced in interactions of cosmic rays with the cosmic microwave background. These so-called cosmogenic neutrinos have not yet been observed; however, as we will discuss in Sec. 3.3.3 observation of anomalous events on detectors hunting for cosmogenic neutrinos has produced significant discussion.
Neutrino flavor conversions are quantum effects happening on a macroscopic scale. Neutrino flavor conversion probabilities are the solution of the effective neutrino Hamiltonian, which is given by Here, the first term is the neutrino mass term which is proportional to E −1 , V is the matter potential, and the last term represents new physics one can explore from the flavor conversion. Since new physics terms are in general zero or positive power of energy at larger energies, the mass term is suppressed, potentially making new physics terms become the dominant term. On top of this, longer baseline oscillation can explore smaller couplings. The quantum effect, long-baseline, and high-energy combination make HENT and EHENT extremely sensitive to many new physics models.
The neutrino telescopes shown in Table 5 can also be organized by their capacity to identify different neutrino flavors. In this respect, neutrino telescopes can be roughly organized as all-flavor detectors or tau-neutrino-specific detectors. The latter are in a special category since tau-neutrinos interactions generate secondary neutrinos that carry a significant fraction of the primary energy [658][659][660][661][662][663][664][665][666][667][668][669][670][671][672] The capacity to identify flavor in these experiments will be crucial for the target models of this white paper. At present multi-flavor neutrino detectors can differentiate well between muonneutrino charged-current interactions and all other neutrino interactions but have limited capacity to separate between tau neutrinos and electron neutrinos. The latter point motivates tau-neutrinospecific experiments that aim to single out that component.  Table 5: Neutrino detectors and neutrino telescopes sorted from smaller to larger energies. These are grouped into three categories separated by the double lines: large neutrino detectors (top), high-energy neutrino telescopes (HENT, middle), and extremely high-energy neutrino telescopes (EHENT, bottom). The detector technology and the neutrino interaction that they are sensitive to are shown in the right most columns with references. The label 'All Flavors' implies that they can also detect both chargedand neutral-current interactions. Table adapted from [691]. Figure 14 shows the energy ranges and the relevant neutrino fluxes. The three fluxes discussed in this section are shown there: atmospheric (orange), high-energy astrophysical (blue), and ultra-highenergy (red). Since neutrino flavor changing depends on the distance traversed by the neutrino the distance from the neutrino source to us is shown in the vertical axis. These spans distances from the Earth radius to several gigaparsecs in length. The region probed by experiments today is signaled by a pink background, while the region that is expected to be probed by next-generation detectors is shaded light blue. Figure 14: Energy and distance scales relevant for neutrino telescopes. Three high-energy neutrino fluxes are labeled as atmospheric (orange), high-energy astrophysical (blue), and ultra-high-energy (red). The region explored by the current experiment is shown in pink, while next-generation is in light blue. Adapted from [692].
Due to the fact that neutrino telescopes need to instrument large volumes, they are sparsely instrumented. This implies that the information we can obtain on each individual event is limited compared to what is available in hodoscopic detectors such as DUNE. Four observables can be extracted out of the neutrino events observed by these detectors. These are: direction, morphology, energy, and time of arrival of the event. Of these four observables, the one more accurately measured is the time of arrival.

High-energy atmospheric neutrinos
Atmospheric neutrinos are produced in collisions of primary cosmic rays and the Earth's atmosphere. These neutrinos have varying propagation distances from their production point, on average above 20 km above the Earth's surface, to their detection. The incident angle of detection of these neutrinos is then a proxy for the propagation distance, where the maximum propagation distance is an up-going neutrino travelling the Earth's diameter. The atmospheric neutrino flux has a steeply falling powerlaw energy distribution, which at low energies follows the cosmic-ray spectrum and at higher energies becomes on spectral index softer due to the interactions of mesons with the atmosphere. These neutrinos, predominantly produced by kaon decay, have been observed to have energies as large as hundreds TeVs; making them the highest energy neutrinos from terrestrial origin. Thus atmospheric neutrinos are the only system to study neutrino flavor transitions at baselines up to ∼ 12, 800 km and energies up to ∼ 100 TeV. A regime far beyond what accelerator-based neutrino experiments can achieve. Since neutrino oscillations act as natural interferometers, therefore, atmospheric neutrino interferometry is unique to test various exotic flavor mixing due to new physics.
The energy region of 1 TeV or lower can be measurable by neutrino detectors with tanks. Super-Kamiokande reported 318 up-going showering muon events from 1645 days of data [771]. These highenergy muons lose energy stochastically most likely originated from muon neutrino over 1 TeV. Thus, we expect similar or a few times larger sample sizes of such events in neutrino tank detectors such as JUNO [235], DUNE [673], THEIA [487], Super-and Hyper-Kamiokande [484,647]. The advantages of these detectors are their resolution compared with neutrino telescopes without detector tanks [772]. In this energy region, however, low-energy arrays of neutrino telescopes, such as ANTARES [674], DeepCore [648], IceCube-Upgrade [434], and ORCA [676] have overwhelmingly larger sample sizes. For example, DeepCore has an effective 10 Mton volume recorded over 300, 000 neutrino events in the eight years run. Neutrino telescopes can use large sample sizes to study the ν µ disappearance and ν τ appearance channels to look for new physics through neutrino oscillations [409,[773][774][775][776], and tank detectors use all channels however much smaller sample sizes. Thus, detector tanks designed for accelerator-based neutrino experiments and neutrino telescopes designed for extraterrestrial neutrino detection overlap in this energy region but differences in technologies make them complimentary.
Currently, the TeV or larger energy region is only accessible by neutrino telescopes such as IceCube/IceCube-Gen2 [434,675], KM3NeT [676], and Baikal-GVD [677]. Over the next decade, these experiments will be joined by next-generation detectors such as P-ONE [678] and TAMBO [679]. This is a unique region to look for exotic effects which disturb standard neutrino oscillation. Figure 16 shows examples from IceCube. The left figure shows an oscillogram with eV-scale sterile neutrino, where ν µ disappearance peak shows up in the TeV range [774,775]. The right figure shows the double ratio of ν µ -disappearance oscillation probability [409]. The numerator is the vertical events (long propagation) and the denominator is horizontal events (short propagation and no exotic oscillation). If the new physics, for example Lorentz violation or quantum decoherence [455], are small effects, we only expect these at high energy where neutrino mass effect is suppressed, then neutrino telescopes with atmospheric neutrinos can be the only place to explore this.

High-energy astrophysical neutrinos
High-energy astrophysical neutrinos were discovered by the IceCube Neutrino Observatory in the South Pole [777]. To date, no other neutrino telescope has measured astrophysical neutrinos, though hints exist in the ANTARES data [778]. IceCube has established and characterized the astrophysical neutrino flux using all available morphologies [354,777,[779][780][781][782][783][784][785]. First, they were measured using high-energy starting events [777,779], which are predominantly cascades, and, in their latest analysis, reported that the flux energy distribution is compatible with an unbroken single power-law with a spectral index of 2.87 [784]. Second, evidence was found using northern sky muon-neutrino events, which are morphologically tracks, and, in their latest analysis, this channel reported a spectral index of 2.37 [785]. Finally, the astrophysical flux was measured by using all-sky cascades and obtained a spectral index of 2.53 [783]. Other notable measurements of the diffuse astrophysical neutrino flux include measurement of it by selecting starting muon-neutrino events [354], by selecting medium energy starting events [780], by searching for double bang signatures [695], and finally by searching for signatures for resonant W -production in electron antineutrino scattering [699].
Studies of the arrival direction of IceCube's neutrinos indicate that the diffuse astrophysical component is compatible with an isotropic distribution [786]. Searches for contributions from the galactic center have so far yielded only constraints and the galactic component has been shown to be no larger than approximately 10% of the diffuse astrophysical flux [787]. Searches for neutrino sources by either looking for clustering, comparing with catalogues of known sources, and looking for real-time and spatial correlation with gamma-rays have also been performed. Evidence has been found for a neutrino source in the direction of the TXS0506+056 blazar by correlating the arrival time of a neutrino alert Earth-core-through going muon-antineutrinos (cos θ = −1) experience a significant matter-resonant-enhanced depletion at TeV energies for an eV sterile neutrino with non-zero muon-neutrino mixing. Figure is taken from Ref. [774]. Right panel shows the ratio of the muon-neutrino survival probability between a scenario with Lorentz violation and the standard case as a function of the neutrino energy. Different color lines indicate different values of the dimensionful coupling of the effective operator that governs the neutrino interaction with the Lorentz violating field. Figure is taken from [409]. and a gamma-ray flare [788]. To date, the sources of bulk astrophysical neutrino flux are unknown and they are the targets of current and next-generation HENTs.
Measurements of the flavor composition of astrophysical neutrinos using different samples are shown in Fig. 17, left. Current measurements are limited by the small sample size and difficulty in separating between different flavors. The expected flavor composition from standard production scenarios and known neutrino oscillation parameters yields a democratic flavor composition [308,343,363,692,747,769,[789][790][791][792][793]. The current constraints are compatible with this expectation and only atearth flavor configurations that significantly deviate from the democratic behaviour are ruled out, e.g. 100% electron-neutrinos or 100% muon neutrinos [411,784]. Despite this limitation, astrophysical neutrino flavor composition has already been used to put some of the most stringent constraints in dark matter-neutrino interactions [270,308,[753][754][755][756][757][758][759][760], and low-energy manifestations of quantum gravity [411], among others [327,331,334,363,692,692,[741][742][743][744][745]. Figure 17 right shows the expected progression of high-energy neutrino telescopes effective volume as a function of time. Notably, within this decade we will five-fold increase the neutrino effective volume and in the next decade the aggregated neutrino telescopes effective volume will be more than ten-fold larger. Not only we expect an increase in the sample size, but also improvements in reconstruction using machine learning are expected to yield improved flavor identification. This is complemented by the development of tau-neutrino-flavor-specific detectors, e.g. TAMBO or Trinity, whose independent measurement would break the degeneracy between electron and tau neutrino flavors. Figure 17 left shows the expected flavor composition measurements at Earth when combining the experiments listed in Fig. 17 right (IceCube [675], IceCube-Gen2 [434], Baikal GVD [677], KM3NeT [676], P-ONE [678], TAMBO [679]). In this figure, we can also see the expected regions from standard production scenarios, which are also expected to become smaller as both the sources are identified and the neutrino oscillation parameters uncertainties are reduced. The effect of the latter is shown in this figure by comparing the size of the same color region with the two different hues, the dark hue is with current uncertainties and the lighter one with expected uncertainties. As illustrated in this figure by 2040 the combined neutrino telescopes, with morphological mis-identification, is expected to be able to disentangle between the different production scenarios and also constraint new physics that produces anomalous flavor compositions. Meantime, as shown in Fig. 17 right, nextgeneration long-baseline oscillation experiments (JUNO [235], IceCube-Upgrade [434], DUNE [673], Hyper-Kamiokande [484]) will improve the oscillation parameter measurements. This is relevant because BSM effects are often sub-dominant processes of the standard flavor conversions, and the ability to access new physics through astrophysical neutrino flavor is limited by oscillation parameter errors. This is shown in various shaded regions in Fig. 17 left, these predicted flavor ratios on the Earth with different astrophysical neutrino flavor production models. In what follows we briefly described the HENT that are expected to contributed to the astrophysical neutrino flavor measurements in the 100 TeV to 10 PeV energy range.
IceCube [675] (build): The IceCube Neutrino Observatory is a gigaton scale ice-Cherenkov detector build in the antartic continent. IceCube is comprised of 86 strings which are arranged on a hexagonal grid instrumenting approximately one cubic kilometer of ice. Each string hosts sixty digital optical modules (DOMs) each of which comprises of one photomultiplier tube with its on data acquisition system.
KM3NeT [676] (underconstruction): The cubic-kilometer scale water Cherenkov detector is under construction in the Mediterranean sea. KM3NeT has two separate sites optimized for different energy ranges. The low-energy site is called KM3NeT/ORCA and is been deployed south of Toulon, France. The high-energy site is called KM3NeT/ARCA and will be deployed close to the sicilian shore. In its final configuration KM3NeT/ARCA is expected to see 11 muon-neutrinos, 41 electron-neutrinos, and 26 tau-neutrinos when assuming an E −2 flux compatible with IceCube's measurements.
Baikal-GVD [677] (underconstruction): The Baikal Gigaton Volume Detector is the next generation detector to be deployed in lake Baikal following the NT-200 detector. The detector has recently reached an effective volume of approximately 0.4km 3 . The full detector will contain 10,386 opticalmodules arranged in 27 clusters of strings and is expected to have an instrumented volume of 1.5km 3 . Due to the relatively large spacing between clusters Baikal-GVD is expected to have a larger energy threshold than IceCube. P-ONE [678] (proposed): The Pacific-Ocean Neutrino Experiment is a proposed water Cherenkov detector to be deployed off the coast of Vancouver island in Canada. The full detector is expected to have 70 strings and be completed by 2030.
TAMBO [679] (proposed): The Tau Air-shower Mountain-Based Observatory is a proposed water-Cherenkov detector to be deployed in the Colca Valley in Peru. The Colca Valley geometry is ideal for the detection of Earth-skimming PeV tau neutrinos. In its final configuration the TAMBO tauneutrino effective area is expected to be ten times larger than IceCube's at a PeV and thirty times larger at ten PeV.
IceCube-Gen2 (ice) [434] (proposed): The next-generation ice-Cherenkov detector to be deployed in Antartica is expected to have a detection volume a factor of ten times greater than IceCube.

Extremely-high-energy astrophysical neutrinos
Extremely-high-energy neutrino telescopes (EHENTs) [689,[795][796][797] aim to discovery neutrinos produced in interactions of cosmic rays with the cosmic microwave background [798,799]. Radio emission is expected from the neutrino interaction at these extremely high energies. Radio has the advantage of traveling longer distances than optical light and allows for the detector to be more spaced, increasing the effective area. Additionally, the polarization of the radio emission can be used to reject downward-going cosmic-ray showers. Several experiments have used this technique to search for cosmogenic neutrinos, yielding so far constraints [685,687,800,801] and anomalous results observed by ANITA [802,803]. These events cannot be explained by standard neutrino interactions [804], since the required flux overshoots the present constraints and due to the lack of observation of the correlated low-energy neutrinos produced in the neutrino transport through the Earth [670,805]. Several explanations have been put forward for this events, these include dark matter [806][807][808][809][810], supersymmetric partners [811,812], leptoquarks [813], sterile neutrinos [814,815], secret interactions [816], axions [817], as well as more mundane explanations [818]. Figure 18 shows the projected sensitivity of next-generation EHENTs; see Ref. [819] for a recent discussion on new physics searches for very high energy tau neutrinos. Using different detection techniques that vary from earth-skimming neutrinos, Cherenkov detectors, and radio. Combining up-going and down-going measurements allow disentangling particle physics (cross-sections) from astrophysics (flux) [820,821], and constrain both the properties of astrophysical sources and the interactions of neutrinos far above the weak scale [822]. In what follows, we briefly describe the detectors that are currently probing this energy range and those proposed to follow up: Auger [687,688] (build): The Pierre Auger Observatory is a hybrid detector consisting of both an array of water Cherenkov surface detectors and atmospheric fluorescence detectors. The main surface array is comprised of 1660 water Cherenkov tanks spread over 3000 squared kilometers and overlooked by 24 air-fluorescence telescopes. The observatory has been in operation since 2008 and so far collected an exposure that exceeds 40,000 km 2 sryr. The Global Cosmic Ray Observatory (GCOS) is a conceptual design of similar observatory but covering 40,000 km 2 , roughly 13 times bigger than Auger.
ANITA/PUEO [685,686] (build/under construction): The ANtarctic Impulsive Transient Antenna (ANITA) is an array of radio antennas installed on a helium balloon surveys Antarctica for ∼ 30 days at a time. ANITA has completed four flights and set the strongest constraints on the astrophysical neutrino flux above 10 10 GeV, as well as observed anomalous events that are yet unexplained. The follow-up of ANITA is called PUEO, and it has recently been approved for construction. PUEO is expected to have a sensitivity that is an order of magnitude larger than its predecessor below 30 EeV and with lower threshold for tau neutrinos [686].
ARIANNA [682] (proposed): Antartic Ross Ice-Shelf ANtenna Neutrino Array is the larger form of ongoing ARIANNA effort. It also aims to detect radio wave from neutrino interactions, however, it utilizes the radio wave reflection by the Ross Ice Shelf, allowing to locate autonomous independent detectors on the ice surface.
RNO-G [684] (under construction): The Radio Neutrino Observatory in Greenland [823] is an array of antennas to be deployed in ice designed to measure the neutrino flux above 10 16 eV. The array consists of surface and deep antenna arrays organized in 35 stations. The deep array will survey a larger effective volume, while the surface array will be used to reject incident cosmic rays.
Trinity [680] (proposed): Trinity is comprised of three stations, each of which deploys an imaging atmospheric Cherenkov telescope on a mountaintop. These are expected to observe the Cherenkov radiation produced by the decay of Earth-skimming tau-neutrinos.
BEACON [689] (proposed): The Beamforming Elevated Array for COsmic Neutrinos (BEACON) uses a sparse array of clustered antennas to search for the radio signal from upgoing tau neutrinos. The concept uses interferometric triggering on top of a high-elevation mountains to maximize each station's sensitivity. Stations of about ten antennas are separated by 1 kilometer and view a large area of the horizon.
RET-N [681] (proposed): Radio Echo Telescope-Neutrino observes radio wave reflection by the clouds of ionized media produced by high-energy neutrinos. This is the second stage of RET-CR (Cosmic Ray).
POEMMA [683] (proposed): POEMMA consists of two identical satellites orbiting the Earth at an altitude of 525 km. These satellites can operate in two different observation modes. First, they can operate in POEMMA-stereo mode, which aims to detect cosmic rays or neutrinos up to 20 EeV by observing the fluorescence emission in up-ward going atmospheric air showers. Second, they can operate in POEMMA-limb mode, where they aim to detect the Cherenkov light produced by an upward going tau neutrino decay. These two modes are complementary as they provide all-flavor and single-flavor detection capabilities.
GRAND [690] (proposed): GRAND is a radio antenna array expected to be deployed in China over a mountain slope. It aims to detect the radio emission produced by high-energy air-shower showers produced by Earth-skimming tau neutrino interactions. The distance between the antennas is expected to be one kilometer; these are then arranged into twenty GRAND10k clusters. Each GRAND10k cluster is comprised of ten thousand antennas that cover approximately 100 kilometerssquared patch of the mountain slope.
IceCube-Gen2 (radio) [434] (proposed): Expected to be located in the antarctic glacier near the current IceCube array, the radio component of IceCube-Gen2 is expected to be composed of stations that combine shallow, sub-surface antennas with antennas deployed at depths down to 100 m. The shallower antennas are expected to have improved signal sensitivity and capacity to reject cosmic-ray airshower. These are complemented by the deeper antennas that provide a larger effective volume and sky coverage. Two hundred stations are expected to be deployed in the South Pole over an area of 500 square kilometers.

Summary and outlook
In this White Paper we have reviewed the theoretical aspects of Beyond the Standard Model (BSM) scenarios that can impact neutrino flavor measurements, leading to deviations with respect to the three-neutrino standard paradigm in neutrino oscillations, neutrino telescopes and core-collapse supernovae. The list of models affecting neutrino flavor patterns is wide and difficult to cover exhaustively in a single document. Therefore, we have defined four general categories accounting for different phenomenological consequences expected in general extensions of the SM. The addition of extra neutrino states (Sec. 2.1) is well-motivated from the theoretical point of view: right-handed neutrinos are singlets of the SM gauge group and pose one of the simplest extensions of the SM able to generate neutrino masses and mixings. These may lead to deviations from non-unitarity, to which neutrino oscillation experiments could be sensitive, or could impact the neutrino oscillation pattern through the addition of new oscillation frequencies. This white paper provides an updated summary of current constraints on non-unitarity parameters and LED-induced oscillations. We have also reviewed the phenomenological consequences derived from the plausible existence of new neutrino interactions in Sec. 2.2. Theory challenges associated to model building which could lead to large effects in the neutrino sector have been reviewed, pointing out possible ways out and model building strategies. We also include a detailed review of constraints applicable to flavored models with light and ultralight mediators, to which neutrino oscillations in matter are sensitive, as well as a discussion of the phenomenological consequences from dark matter-neutrino interactions. Next, we have reviewed the main consequences of neutrino decay in a variety of setups in Sec. 2.3.3. Neutrino decay searches offer an indirect way to test the neutrino mass generation mechanism, which is a priori unknown. The phenomenological consequences are diverse, and here we review the main consequences from both visible and invisible neutrino decay, which impact neutrino flavor measurements in different ways. Tables 2 and 3 summarize the current bounds and expected sensitivities for a variety of experiments. Finally, given the singular quantum-mechanical neutrino properties, neutrinos also provide a unique avenue for testing fundamental symmetries as CPT or Lorentz invariance, as well as non-standard quantum decoherence. The phenomenological consequences derived from these scenarios are reviewed in Sec. 2.4, which includes a review of current bounds (see, e.g. Tab. 4) as well as a summary of future prospects to improve over these using neutrino oscillation experiments and neutrino telescopes.
Neutrinos from different sources spans several orders of magnitude in neutrino energy (ranging from keVs to EeVs), and therefore to be detected, facilities with different experimental requirements have been (and will be) built. At the same time, BSM effects might be present at neutrino production, propagation, and detection affecting each one of this stages in different ways. In particular, BSM scenarios that impact neutrino flavor leave signals at more than one experiment. Fortunately, several future neutrino experiments have been approved or are already being built, and will allow for exciting BSM searches. In the future, identifying all neutrino flavors is important, which calls for multi-channel observation.
Understanding complementarity between different techniques and facilities is also relevant. For instance, developments in detector technologies motivated by direct dark matter searches have finally made detecting coherent neutrino nucleus elastic scattering process feasible. Additionally, direct dark matter searches are starting to be sensitive to solar and atmospheric neutrino scattering and therefore can provide additional constraints on new physics affecting the neutrino floor, as discussed in Sec. 3.1. Section 3.2 discussed medium-energy experiments, including neutrino beam experiments and atmospheric neutrino detectors. Neutrino beam experiments can provide beams with relatively well-known spectra, flavor content, and timing information, which offer multiple advantages with respect to neutrino fluxes from natural sources. At these, the use of near detectors is crucial to reduce the level of uncertainties and allow to search for BSM effects using their far detector data. Past and current experiments have performed BSM searches using this two-detector setup, and future beambased experiments will build upon past experience and efforts in this regard, leading to state-of-the-art facilities, with higher power, and a reduced level of systematic errors. Atmospheric neutrino experiments are complementary to low-energy and beam-based experiments. The atmospheric neutrino flux spans several orders of magnitude in energy (tens of MeV to hundreds of TeV) and neutrinos can travel different distances from the production point to the detector (tens of km to the Earth diameter). BSM searches benefit from both higher energies and baselines in comparison to long-baseline experiments, and also from the experience gained with multi-purpose detectors used to look for different physics. Upgrades of current atmospheric experiments include increasing their volume (to gather more statistics) or improving the detector technology in order to gain energy and angular resolution as well as particle identification, which aids in BSM searches. Finally, neutrino telescopes (described in Sec. 3.3) provide a unique window for high-and extremely high-energy neutrinos. These facilities aim to measure the high-energy part of the atmospheric neutrino flux and astrophysical neutrino sources. Extremely high-energy neutrino telescopes aim to detect particles with the highest energies on the Earth surface. Given the energy ranges and the very long distances traveled by the neutrinos detected (from the Earth radius to several gigaparsecs), neutrino telescopes provide a unique avenue to probe BSM scenarios.
Clearly, the best way to confront BSM models with data is by exploiting the synergies between different experiments and probing BSM effects across a wide range of energies, and for multiple detector technologies. Thus, next generations experiments will be crucial for discovering or ruling out the new physics scenarios reviewed here, that are of interest of the particle physics community.