New sources of leptonic CP violation at the DUNE neutrino experiment

We check the capability of the DUNE neutrino experiment to detect new sources of leptonic CP violation beside the single phase expected in the Standard Model. We illustrate our strategy based on the measurement of CP asymmetries in the case New Physics will show up as Non-Standard neutrino Interactions and sterile neutrino states and show that the most promising one, once the experimental errors are taken into account in both scenarios, is the one related to the $\nu_\mu \to \nu_e$ transition.


Introduction
Neutrino flavor oscillations have been one of the most important particle physics discoveries of the last decades [1]. Experiments using natural and artificial neutrino sources were able to measure mixing parameters with a good precision [2,3]. In particular, the three mixing angles, as well as the absolute values of the mass splittings, have been determined with a few percent uncertainty. However, some degeneracies still exist and future experiments aim to solve them. In particular, the atmospheric mass splitting sign (hierarchy problem) and the octant in which the atmospheric angle θ 23 lies are still unknown. Moreover, the phase of the PMNS matrix, which is the only source of CP violation in neutrino oscillation 1 , has not been measured with a satisfactory precision [4,5].
Even though most of the experimental results are in agreement with the three neutrino paradigm, the long-standing short baseline [6]- [8] and reactor [9]- [12] anomalies gave some hint of New Physics (NP) phenomena such as the presence of additional sterile neutrino states [13]- [24]. Moreover, uncertainties on the oscillation parameters leave rooms for the presence of additional effects in neutrino interaction [25]- [27] and propagation in matter not contemplated in the Standard Model (SM) and that can be described in a model independent way by four fermions effective operators, namely the Non Standard Interactions (NSI) operators [28]- [39].
Both types of NP models introduce new sources of CP violation in neutrino oscillations. In particular, when one sterile state is contemplated, the 4 × 4 PMNS matrix contains two more phases than in the SM case. On the other hand, if we take into account the possibility that during their travel through matter neutrinos can undergo NSI, three new phases in the state evolution equations emerge.
The NP phases inevitably affect CP violation in neutrino oscillations, whose presence can be highlighted from CP-odd observables; among them, asymmetries of the type A αβ ∼ P (ν α → ν β ) − P (ν α →ν β ) can be directly measured at future neutrino experiments capable of distinguishing neutrino from antineutrino events. Such quantities are generally dependent on the CP phases, and become a smoking gun for new CP violation as soon as their values deviate from the SM predictions, including matter effects [40].
The goal of the present paper is to test whether neutrino experiments will perform sufficiently well as to establish that other phases beside the single one expected in the SM are responsible for CP violation in the lepton sector. To achieve our aim, we adopt the following strategy: • compute the leptonic asymmetries in the SM; • evaluate the experimental uncertainty on them; • re-compute the asymmetries including the effects of NP; • check whether the new results are sufficiently away from the SM predictions.
The outcome of our procedure is the identification, if any, of the most promising asymmetry capable to produce an experimental signature well beyond the SM expectations.
We validate our strategy studying the asymmetries as they can be measured at the Deep Underground Neutrino Experiment (DUNE) [41]- [44]. An unprecedented feature of DUNE will be the possibility to study three different transition channels (namely ν µ → ν µ , ν µ → ν e , ν µ → ν τ ) and the Neutral Current (NC) neutrino interactions using both neutrino and antineutrino fluxes; thus, this experiment will be able to measure four independent asymmetries.
The manuscript is organized as follows: in Sect.2, we discuss the analytic structure of the relevant asymmetries in the SM, by means of the perturbation theory in the small parameters (to be defined later in the manuscript) and in the regime of small matter effects (a perfectly viable approximation for DUNE); we then repeat the same calculation in the NSI scenario (Sect.3) and in the 3 + 1 sterile neutrino model (Sect.4). Sect.5 is devoted to a description of the DUNE experiment and the impact of NP of the neutrino energy spectra while in Sect. 6 we show the foreseen experimental values on the asymmetries built in terms of number of events (integrated asymmetries) and their related uncertainties, both in the SM and in the NP scenarios. Eventually, in Sect.7 we draw our conclusions. Expressions of the probabilities in the NP models analyzed here are provided in the Appendices A and B.

CP asymmetries in the Standard Model
The neutrino state evolution equations when they travel through matter can be written as: where U is the usual neutrino mixing matrix and A CC ≡ √ 2G F n e , with n e being the electron density in the Earth crust. Defining P (ν α → ν β ) as the transition probability from a flavor α to a flavor β, one can construct the CP-odd asymmetries as: It is well known that matter effects modify the behaviour of the asymmetries as a function of the Standard Model CP phase δ (see, e.g., [45]): in fact, the passage of neutrinos through matter introduces fake CP-violating effects which allows A αβ = 0 even when sin δ = 0. In principle, to extract genuine CP violating effects, one could defines the subtracted asymmetries as A sub αβ (δ) = A αβ (δ) − A αβ (δ = 0). However, we prefer to deal with more directly measurable quantities and we will use eq.(2) which, for non negligible matter effects, are non vanishing when δ = 0, ±π.
Being the numerator and the denominator of eq.(4) doubly suppressed by small quantities, we expect A SM 0 µe ∼ O(1). For the ν µ → ν τ channel, on the other hand, we find that the leading contribution to the asymmetry is given by a simpler expression: which is clearly smaller than A µe . Notice also that, differently from A µe , this asymmetry becomes negative if δ > 180 • , as emerging from fits to neutrino oscillation data [2,3]. A third possible asymmetry, namely A µµ , is obviously vanishing in vacuum because of CPT conservation but can assume a relevant role when matter effects are taken into account (as we will discuss later on).
As it is well known, the inclusion of matter effects complicates the analytic expressions of the transition probabilities and, more importantly, that of the asymmetries. In order to deal with readable formulae, we can work in the regime of week matter potential V CC 1 which, as outline before, is a good approximation in the case of DUNE. Thus, we can organize our perturbative expansion as follows: where A SM 1 αβ represents the first order correction to the vacuum case V CC = 0. Thus, the asymmetries considered in this study acquire the following corrections: It is evident that A µe increases because a term proportional to r 2 /f 1 appears, which is of O(1). Since at the atmospheric peak sin ∆ 31 cos ∆ 31 , the r 2 /f 1 correction is positive and adds an O(V CC ) contribution to the total A µe , that at the DUNE peak energy becomes roughly 1/2. On the other hand, A SM 1 µτ contains only terms proportional to V CC r 2 and V CC α 2 which are not balanced by any small denominator. Thus, both contributions set a correction to the vacuum asymmetry. A similar situation arises for A µµ , where only terms proportional to V CC rα and V CC α 2 appear.

NSI and CP asymmetries
As mentioned in the Introduction, the uncertainties on the mixing parameters leave room for the possibility of the presence of Non Standard Interactions between neutrinos and the particles they meet travelling through the Earth. The strength of such new interactions can be parameterized in terms of the complex couplings ε αβ = | αβ |e iφ αβ , which modify the matter potential of eq.(1) to: Since the Hamiltionian has to be Hermitian, the three diagonal couplings ε αα must be real. Moreover, we can always subtract a matrix proportional to the identity without changing the transition probabilities. If we choose to subtract ε µµ I, only two independent diagonal parameters (ε ee = ε ee − ε µµ and ε τ τ = ε τ τ − ε µµ ) will appear in the NSI matrix 2 . Thus, beside the standard oscillation angles and phases, the parameter space is enriched by five more moduli |ε αβ | and three more phases φ αβ , which could provide new sources of CP violation in the lepton sector.

Asymmetries in the NSI framework
Since NSI effects are strongly intertwined with standard matter effects driven by V CC , the asymmetries can be cast in a form which generalizes eq. (7): where A SM 0,1 αβ refers to the pure Standard Model results and all the effects of the NSI are included in the A N SI αβ term. Bounds on the magnitude of the NSI couplings have been widely discussed [34]; even though some of them could in principle be of O(1) and give rise, for example, to degeneracies leading to the so-called LMA-Dark solution [48], we decided nonetheless to consider all ε αβ 's on the same footing and of the same order of magnitude as the other small standard parameters a, s, r, α and V CC . In this way, we are able to catch the leading dependence on NP carried on by the CP asymmetries.
For the ν µ → ν e channel, the leading order NSI contributions can be arranged as follows: where the a s functions are given by: The only NP parameters appearing at the considered perturbative level are ε eµ and ε eτ which, in turn, carry the dependence on the CP phases δ eµ , δ eτ . All in all, the NSI contributions set an O(V CC ) correction to A SM 0 µe . We also notice that the largest of the considered terms, namely the ones linear in r in the numerator, have similar expressions in both a εeµ µe and in a εeτ µe , apart from the sign in front of cos ∆ 31 . This means that, around the atmospheric peak, the phases δ eµ and δ eτ are equally important even though the magnitude of their impact strongly depends on the value of the standard CP phase δ.
For the asymmetry in the µτ -channel, we found the following structure: In this case, four different NSI parameters enter the leading order corrections, namely ε µτ , ε eτ , ε eµ (together with their phases) and ε τ τ . Contrary to the µe case, the largest correction to the vacuum expression is given by the first order term ε µτ in the first line of eq.(15), which is not suppressed by any of the standard small parameters a, r, s and α. Considering that A SM 0 µτ ∼ O(rα), this makes the µτ -channel very promising for searching for NP, at least at the probability level where possible complications due to small τ statistics do not enter.
Finally, for the µµ channel, matter effects generate a substantial difference in the propagation of neutrinos versus antineutrinos, which results in the following NSI contributions: As expected from unitarity relations, we get an opposite linear dependence on ε µτ but with a coefficient proportional to tan ∆ 31 which, close to the atmospheric peak, gives an important correction to A SM 1 µµ .

Sterile neutrinos and CP asymmetries
The next NP scenario under discussion is the so-called 3+1 model, in which a sterile neutrino state supplements the three standard active neutrinos. Even though the new state cannot interact with the ordinary matter, it can have a role in neutrino oscillations thanks to the mixing with the active partners. The long-standing reactor, gallium and short-baseline anomalies [21] suggested that, if present, the fourth mass eigenstate m 4 should have a mass such that ∆m 2 that is orders of magnitude larger than the solar and the atmospheric mass splittings, thus capable to drive very fast oscillations visible at accordingly small L/E. In addition to the new mass splitting ∆m 2 41 , the PMNS matrix becomes a 4 × 4 matrix which can be parametrized in terms of 6 angles and 3 phases. In this manuscript, we adopt the following multiplication order of the rotation matrices R(θ ij ) [49]- [51]: Apart from δ 2 , which becomes the standard CP phase for m 4 → 0, we have two potential new sources of CP violation, encoded in two phases δ 1 and δ 3 . In the description of neutrino propagation in matter, we cannot disregard the role of the NC interactions because the sterile state does not feel at all the presence of matter; this results in the following evolution equations: where ν s is the new sterile state, A CC is the usual matter charged current potential and A N C is the matter NC potential, A N C ≡ 1/ √ 2G F n n , with n n being the neutron density in the Earth crust.

Asymmetries in 3+1 framework
The parameter space of the 3+1 model is enlarged compared to the SM case by three new mixing angles θ i4 , two more CP phases δ 1,3 and the mass-squared difference ∆m 2 41 . Thus, in addition to the expansion parameters used in the previous sections (r, s, a), we also expand in the small s 14 , s 24 and s 34 (where s i4 = sin θ i4 ) that we can still assume of O(0.1). To further simplify the analytic expressions of the asymmetries, we also introduce V N C = A N C L/2∆ 31 . It is useful to present the results in a form similar to eq.(11): where A SM αβ are the SM asymmetries and the symbol λ represents a common order of magnitude of all small quantities used in our perturbation theory, including V CC (but not V N C , whose dependence in A αβ is exact). The exponent amount to n = 3 for A µτ , A µµ and n = 2 for A µe . Notice that, due to the parametrization adopted in this manuscript, the SM phase δ of eqs.(4)-(10) must be replaced by the combination δ 2 − δ 1 − δ 3 . Averaging out all the fast oscillations driven by ∆m 2 41 , the various A 3+1 αβ have the following expressions: To avoid large expressions, for A µe we only quote the corrections due to the new mixing angles.
First of all, we notice that the corrections to the µe asymmetry are only linearly suppressed compared to the leading order results; thus, we expect such an asymmetry to be quite sensitive to new sources of CP violation. Then, both corrections to the µτ and the µµ asymmetries are linear in the combination s 24 s 34 . Since the angle θ 34 has weak constraints (values of 20−30 • are still allowed), these corrections can be relatively large. Notice also that, since V N C is roughly of the same order of magnitude as V CC , A 3+1 µτ is expected to provide a large correction to the standard model asymmetries, making the ν τ appearance channel, at least in principle, very sensitive to NP effects. As for the PMNS phases, all leading order corrections depend only on the new phase δ 3 . This means that a long baseline experiment is mostly sensitive only to the combination δ 2 − δ 1 − δ 3 and to the single phase δ 3 .
Beside the results of eq.(20), it is worth considering a new asymmetry corresponding to the ν µ → ν s transition. Even though sterile neutrinos cannot be directly detected, the probability P (ν µ → ν s ) is a measure of the NC events in the detector. Indeed, being the NC interactions flavor independent, the number of events is proportional to the sum of the transition probabilities from the starting flavor (ν µ ) to the three active final flavors (ν e,τ,µ ) because of the unitarity relation P (ν µ → ν s ) = 1 − P (ν µ → ν e,µ,τ ). The new asymmetry has vanishing matter corrections and, at the leading non-vanishing order, reads: This is clearly an O(1) result since both numerator and denominator are of O(λ 2 ). In Tab.
(1) we summarize the outcome of our analytic considerations on the magnitude of the NP corrections to the asymmetries discussed in this paper. Table 1: Order of magnitude estimates of the various contributions to the asymmetries discussed in this paper. λ is a common order parameter such that: r, s, a, ∆ 21 , V CC , ε αβ , θ i4 ∼ O(λ).

The DUNE experiment
The DUNE (Deep Underground Neutrino Experiment) experiment is a proposed long-baseline experiment based in the USA [41]- [44]. The accelerator facility and the Near Detector will be located at Fermilab, while the Far Detector is going to be built at the SURF (Sanford Underground Research Facility) laboratories in South Dakota, 1300 km away from the neutrino source. Both far and near detectors will be LAr-TPCs, namely detectors with very good imaging capabilities, which are expected to collect a huge number of different neutrino interaction events.
The neutrino beam will be a ν µ beam with a small ν e contamination. Focusing horns, which are able to select particles with a given electric charge before the decay tunnel, will produce particle and antiparticles beams, allowing the experiment to run in two different modes, namely the ν-mode and theν-mode. The neutrino flux energy spectra should be peaked at E peak = 2.5 GeV, however different proposal have been promoted for higher energy fluxes. Indeed, even though E peak is the energy of the atmospheric peak of the oscillation probabilities at 1300 km baseline, we are below the τ production threshold (E thr = 3.1 GeV). Thus, with such a flux, CC interactions of the huge number of ν τ -s arriving at the far detector are forbidden. A broader and more energetic flux would overcome this problem, allowing the tau neutrinos to be energetic enough to produce τ leptons. This τ -optimized flux [52, 53] would be less performing in constraining oscillation parameters from ν e appearance channel (due to the increased number of background events such as misidentified ν τ -s) but at the same time very useful for NP searches thanks to the increased number of ν τ events. In order to simulate the DUNE experiment, we used the GLoBES package [54,55]. Initial fluxes, far detector efficiencies and energy resolutions, backgrounds and systematic uncertainties have been provided by the DUNE collaboration for the ν e appearance and ν µ disappearance channels [56,57]. For the former, backgrounds are misidentified ν µ , ν τ and NC events, as well as ν e -s from the flux contamination; for the latter, the main background source are NC events. Systematic normalization uncertainties at the far detector are 2% for the ν e appearance signal and 5% for the ν µ disappearance signal, numbers proposed by the DUNE collaboration from the foreseen performance of the Near Detector.
In the last years, the possibility to study also the ν τ appearance and the NC channels as signals have been taken into account. For the first one, the use of the hadronic and electronic decays of τ leptons to identify the event topology have been discussed in [58,59] Table 3: Best fits for oscillation parameters obtained by the global analysis in [2].
the ν µ events as a background have been proposed in [60] and adopted here.

Effects of NSI and sterile neutrinos on DUNE spectra
Before discussing the sensitivity on the CP asymmetries, it is useful to have a look at the effects of New Physics on the neutrino spectra, which will help in the interpretation of our numerical results. For the NSI case, different global analyses on oscillation experiments have been done [30,33,34]. The 2σ current limits on the various ε αβ from ref. [30] have been summarized in Tab.2. Scanning the NSI parameters in the allowed ranges and taking the new phases δ eµ ,δ eτ and δ µτ in the range [0, 2π] (both sets of parameters extracted randomly flat), we get the neutrino and antineutrino spectra at the far detector, as shown in Fig.1. The number of events is normalized by the bin width and the exposure for both ν-mode andν-mode is 3.5 years. In all panels, the blue dashed lines refer to the spectra obtained in the standard oscillation framework with normal hierarchy (NH), for which the best fit values are summarized in Tab.3. The figures clearly show that the ν e andν e spectra are the ones affected the most by the NSI parameters. This is because at the DUNE energies and baseline, the SM ν e appearance probability is suppressed and the effect of the NSI parameters, in particular that of ε eτ which has weaker bounds compared to the others, results more evident. Conversely, in the ν µ spectra the NSI have a very small impact: in fact, in the disappearance probability the first term is of O(1), and the largest NSI corrections is driven by ε µτ which has very strong bounds ∼ O(10 −3 ). Finally, in the ν τ appearance channel the modifications in the spectra due the NSI's are evident but, due to the small number of expected events, the changes with respect to the SM case are difficult to observe. If we repeat the same study using the high energy τ -optimized flux, (see Fig.2 3 ), we get very similar features as before but for the ν τ events, which are obviously much larger.
In the 3+1 framework, bounds on the parameters strongly depend on the adopted parametrization of the 4 × 4 mixing matrix U . In addition, transition probabilities are affected by many degeneracies among standard and non-standard mixing angles that makes the extraction of the allowed/excluded ranges more complicated. However, in the parametrization of eq.(17), if we allow the standard parameters to vary only in their allowed ranges (Tab.3) and we fix the new mass splitting to be ∆m 2 41 ∼ 1 eV 2 , different studies [21,22] suggest that θ 14 and θ 24 can be taken in the range [0 − 10] • while θ 34 in the range [0 − 30] • . The variations of the neutrino spectra for the 3+1 model are shown in Figs.3 and 4. As for the NSI case, the largest deviation from the SM results are found in the ν e spectra, since, as showed in eqs. (20), the 3+1 corrections to the standard asymmetry are at the first order of our expansion.
As for the other flavors, the ν µ spectra are not really affected by NP, while the deviation of the ν τ spectra from the SM predictions can be mostly ascribed to θ 34 and its relatively large allowed range. We also see that the contributions to NSI is mailny negative, as it should be because the corrections to the probabilities driven by θ 34 are negative indeed. Computing the same spectra using the τ -optimized flux (Fig.4), the changes in the ν τ ap-

Numerical evaluation of asymmetries
The relevant question now is related to the experimental capability to measure the asymmetries we are considering: in fact, if the CP violating quantities will not be measured with a sufficient precision, then we cannot distinguish the deviation from the SM results due to NP. Instead of considering the asymmetries at the probability level, we deal with the experimentally relevant integrated asymmetries built from the number of expected events N β and N β : where the event rates for the ν α → ν β and the CP conjugateν α →ν β transitions are computed from: in which σ β(β) is the cross section for producing the lepton β(β), ε β(β) the detector efficiency to reveal that lepton and φ α(ᾱ) the initial neutrino flux at the source. Since in the SM the only dependence on the CP phase is carried on by δ, the correlations between the pair of asymmetries, for instance (A µτ , A µe ) and (A µτ , A µµ ), is maximal and a close curve appears in the related physical planes. If, in addition, we also take into account the experimental errors on angles and mass differences, the curves are scattered as reported in Fig.5, for the DUNE standard flux and in Fig.6 for the optimized flux. The blue dots are obtained using parameters in the normal hierarchy, while the orange ones are obtained using the inverse hierarchy hypothesis. The first striking features of the integrated asymmetries is related to the fact that their sign is always positive; in fact, being integrated quantities, they are influences not only by the relative differences among ν andν probabilities, but also by the differences among ν andν fluxes and cross sections. As we can observe in Fig.1, the SM spectra (blue lines) ofν is always lower, and this helps in justifying the observed signs. The other important observation is that, as discussed above, A µe is the asymmetry that changes the most with a change of the CP phase. On the other hand, the other two asymmetries A µτ and A µµ change at a much slower rate.
Eventually, it is worth mentioning that, for each pair of asymmetries, the closed curves corresponding to NH and IH never overlap. This means that, at least in principle, one could be able to solve the neutrino hierarchy problem simply looking at the CP asymmetries. However, in DUNE as well as in other future experiments, the foreseen experimental errors on such asymmetries will probably be too large to allow for such a discrimination, as we will discuss later on in the manuscript.

Numerical evaluation of the asymmetries in presence of NSI
Now we are ready to apply our strategy to check whether other sources of CP violation carried on by NP can be sufficiently distinguished from the SM phase. In order to do that, we first need to evaluate the experimental errors on the SM asymmetries and then recompute them as predicted by the NSI and the 3+1 sterile models. From Fig.5 we see that the uncertainties on the standard angles and mass splittings are not playing an important role. A simple but accurate estimate from error propagation gives: where δN is the uncertainty related to the number of expected events which receives contributions from the systematic error (normalization errors cited in Sect.5) and the statistical error. For the ν µ disappearance channel, the first source of uncertainty is always dominating, since the number of events is very large and the statistical error is reduced. On the other hand, in the other two channels both terms are important. In particular, in the ν τ appearance, systematic errors are quite large (due to the poorly known cross section and to the systematics related to the complicated event reconstruction) and the number of events is small. Thus we expect δA µτ to be particularly large.
In Fig.7 we show the values of the asymmetries where the effects of NSI are taken into account, computed by using the standard neutrino flux. The blue stars represent the asymmetries in the standard case (fixing all the standard parameters to their best fits 4 but varying the values of δ), while the orange dots are the results obtained in presence of NSI, computed from the number of events corresponding to random flat extraction of the couplings in the ranges shown in Tab.2. The sides of the grey rectangles represent the maximum 1σ error bars on the standard asymmetries at different chosen values of δ as computed from eq. (25). For the sake of illustration, we do not show here the error bars associated to the NSI points because the number of events is not much different from the standard case, thus the error bars in the NSI framework are of the same order of magnitude as the displayed ones. It is clear that A µτ is very sensitive to New Physics. Indeed, the SM asymmetry has almost a fixed value A µτ ∼ 0.245, as showed in Fig.5, while the NSI contributions can turn A µτ into the range [0.21,0.27]. However, the error bars are much larger than the produced variation, making this asymmetry at the DUNE conditions not useful for discerning new CP phases. Even though the A µe asymmetry gets very different values in the standard case (in the range [0.28,0.55]), the inclusion of the NSI is able to even extend the foreseen asymmetry beyond such a range, enough to reach values outside the error bars of the standard asymmetries. The problem in this case is that, as discussed before, we should also take into account the error bars on the orange dots so that, when we include them, also A µe cannot give hints of NP at DUNE. Finally, for the A µµ the very same analysis done for A µτ applies.
With a higher energy flux, the results partially differ from what illustrated above (see Fig.8). Even though the larger number of events reduces the error bars, the A µτ and A µµ with NSI do not change enough in such a way to be clearly distinguished at an acceptable confidence level from the SM case. On the other hand, A µe can assume values very different For sake of simplicity, we do not report here the error bars on the orange dots. Standard neutrino flux has been employed to compute the number of events.
from the SM ones, in particular, a sets of NSI parameters can push it toward negative values. Indeed, as it is clear from eqs. (12)(13)(14), NSI corrections to the asymmetries can be comparable to the SM case when ε eµ and ε eτ are of O(0.1). With higher energy fluxes, the appearance transition probabilities are mainly evaluated off peak, making the cosine of ∆ 31 in eqs. (13,14) no longer negligible. Thus NSI corrections become more and more important, causing an opposite sign of the asymmetry with respect to the SM case when cos(δ −δ eµ,τ ) terms become negative.
6.2 Numerical evaluation of the asymmetries in the 3 + 1 sterile neutrino model In Fig.9, we report our numerical results for the 3+1 case, obtained for fixed ∆m 2 41 = 1 eV 2 and all mixing angles and phases extracted randomly flat in the ranges discussed in Sect.5.1. Standard neutrino fluxes have been employed. As previously mentioned, we have four independent asymmetries. Three of them (A µe , A µµ and A µτ ) are accessible through the corresponding oscillation channels. The other one, namely A µs , can be measured looking at the NC events. Indeed, since the NC interactions are flavor independent, the number of events in this channel depends on the sum: which, from the unitarity relation, corresponds to 1 − P (ν µ → ν s ). Thus, the integrated asymmetry is closely related to the µs asymmetry. We present our results in the (A µτ , A µe ) and (A µµ , A N C ) planes, see Fig.9, for the standard flux. The situation is quite clear: even though the analytic corrections to A µe ∼ O(λ) and to A µµ,µτ ∼ O(λ 2 ), the relatively large uncertainties do no allow the 3+1 points to spread outside the error bars. As before, the use of the higher energy flux reduces the error bands and increases the number of points outside the SM uncertainties, see Fig.10. Furthermore, as for the NSI case, the asymmetry which vary the most when NP enters into the game is A µe since, as shown in (20), the correction to SM asymmetry is at first order in our perturbative expansion. It is clear from the left panel of Fig.10 that there are some points at more than two sigmas away from the standard values but, differently from the previous case, A µe never becomes negative. On the other hand, in the (A µµ , A N C ) plane, no orange point lies outside the grey rectangle.

Conclusions
The search of Physics beyond the Standard Model has become an attractive research field in the neutrino sector thanks to the huge experimental efforts in the measurement of the standard oscillation parameters. Beside few exceptions, they are known with a very good precision and the relevant question is now to establish whether New Physics is hidden within the experimental uncertainties. In this respect, the quest for new sources of CP violation beyond the single phase already present in the Standard Model is a pressing one. Contrary to many similar studies in the literature, we avoided to focus on the sensitivity of a given experiment to a particular phase, as many assumptions are usually done in the fit procedure (marginalization over a subset of parameters while other are kept fixed, and so on) which obscure the true sensitivity to CP violation. Instead, we designed a more objective strategy based on the evaluation within the Standard Model of the integrated CP-odd quantities, the asymmetries A αβ , and compare them with the same observables evaluated in the New Physics scenarios under consideration. In this paper, we performed a perturbative analytic evaluation of all asymmetries accessible at the DUNE experiment under the assumptions that Non-Standard Interactions and one sterile state affect the standard neutrino oscillation framework. We then apply our procedure to the more realistic asymmetries built from the expected number of event in DUNE, reaching the conclusion that, for both New Physics scenarios, the A µe asymmetry can reach values well beyond the Standard Model expectation, including the foreseen statistics and systematic uncertainties, when an high energy flux is employed. A special mention should be devoted to A µτ : while analytic considerations indicates that New Physics sets large corrections compared to the Standard Model results, the uncertainties involved in the evaluation of the number of expected events obscure this important feature. An experimental effort should be carried out to reduce the uncertainties in τ detection.

Appendix A: perturbative expressions of probabilities
We provide here the perturbative expressions of the probabilities expanded as discussed in Sects.2, 3 and 4, in both NSI and the sterile neutrino models. The standard model probabilities can be obtained putting all new physics parameters to zero. For the NSI case we have: For the 3+1 model, we found the following expressions (δ = δ 2 − δ 1 − δ 3 ): In this appendix, we will provide the perturbative expressions for the asymmetries in the 3+1 model in the case the oscillations driven by ∆m 2 41 cannot be averaged out. The leading order of A µe is unchanged, thus there are no corrections to the SM at the chosen perturbative order.
For the other asymmetries, we give here a first order expansion in the matter potentials V CC and V N C . Thus, we put the corrections in the following form: For the µτ asymmetry we have: For the µµ asymmetry: