Exploring non standard interactions effects in T2HK and DUNE

Neutrino oscillations in matter offer a novel path to investigate new physics. One of the main goals of neutrino experiments is to determine the CP phase, and the presence of new physics can alter the scenario. We assume that the observed difference, if any, in the CP phase is due to the possible non-standard interactions. We derive the relevant coupling strengths using the results of NOν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu $$\end{document}A and T2K and study their effects in the next generation of long-baseline experiments: T2HK and DUNE. Our analysis reveals a significant impact on the sensitivity of atmospheric mixing angle θ23\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _{23}$$\end{document} in the normal and inverted orderings. Furthermore, we observe discernible differences in probabilities for both experiments when non-standard interaction from e-μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e-\mu $$\end{document} and e-τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e-\tau $$\end{document} sectors are included.


Introduction
Accelerator-based neutrino experiments offer exciting avenues to study neutrino physics. They travel long distances (a few hundred kilometers for the long baseline experiments currently underway) and are detected far from the source. Neutrinos change their flavor while going from one place to another and mix among the various mass eigenstates. Interestingly, neutrino oscillations [1,2] provide us with indirect signature of physics beyond the standard model. On entering the Earth's atmosphere and travelling through the Earth's crust, the neutrino gets influenced by a matter potential known as the Wolfenstein matter effect. Wolfenstein in addition to the neutrino mass matrix [3], introduced non-standard interaction (NSI) to investigate new physics. There have been extensive studies of neutrino phenomenology in the literature [4][5][6][7][8][9][10][11][12]. NSI describes new physics models at low energies in which neutrino interactions with ordinary matter are parama e-mail: ph19resch11001@iith.ac.in (corresponding author) b e-mail: giria@phy.iith.ac.in eterized in terms of effective coupling αβ (defined later in the text) [13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29]. For the latest review on NSI and related work, one can see Ref. [30] and the references therein.
NSI effects can exhibit interesting signatures that can be differentiated from the standard model predictions and thus provide golden avenues to decipher new physics. NSI arises naturally from the high-energy physics of new heavy states [30][31][32] or light mediators [33][34][35]. It is quite crucial to comprehend how NSI affect standard neutrino oscillation in matter. However, in this article, we will not resort to any specific scenario but assume that new physics arises only from the NSI and is responsible for any deviation. In general, NSI affects neutrino propagation in matter not only through neutral-current interactions but also through charged-current interactions, which affect neutrino generation and detection. Neutral current NSI interacting with matter fields can yield observable effects. As model-independent bounds on the production and detection of NSI are typically orders of magnitude higher than those of matter NSI, we ignore the production and detection of NSI in this study and focus solely on NSI due to propagation.
The standard model (SM) CP phase promises to help us understand the baryon asymmetry of the universe and is the most sought-after observable in the currently running and future neutrino experiments. The most recent results from the two long-baseline accelerator experiments, NOνA and T2K, show some kind of tension in the standard 3-flavor scenario. NOνA detects neutrinos at the far detector, which is 810 km away from the source, with an energy of approximately 2 GeV. T2K detects neutrinos with an energy of 0.6 GeV at far detector, which is 295 km away from the source. Both NOνA and T2K are off-axis experiments and they detect a stream of neutrinos with a very narrow energy distribution.
According to the recent results, NOνA prefer the CP phase to be close to δ C P ≈ 0.8π [36] whereas T2K hints a value of δ C P around 1.5π [37] in the case of normal ordering.
There appears to be no disagreement in the case of inverted ordering. Once the NSI from the e − μ sector is taken into account the tension concerning the δ C P parameter for NOνA and T2K becomes placid, but one can see a difference for θ 23 [38,39]. NOνA prefers lower octant, whereas T2K prefers higher octant. We extracted datasets of NOνA [40] and T2K [41] from the recent data release in order to find the constraints on NSI contributions. Thereafter, we use the same coefficients to see if we can get any discernible result in future long-baseline (LBL) neutrino experiments, such as the DUNE and T2HK. In particular, T2HK, with its short baseline and low energy, will not be as sensitive to the matter effect as will be DUNE, with its relatively larger baseline and high energy. The objective here is to determine whether or not the degeneracy for the standard model parameter θ 23 persists in the presence of NSI arising from both e − μ and e − τ sectors for DUNE and T2HK. In addition to that, we also explore the question of mass ordering through oscillation probability plots and sensitivity to the CP violating parameter δ C P .

Formalism
The NSI can be characterised by six-dimensional fourfermion ( f f ) operators of the form [3]: where α, β = e, μ, τ indicate the neutrino flavor, superscript C = L , R refers to the chirality of f f current, f = u, d, e denotes the matter fermions and f C αβ are dimensionless parameters that measure the new interaction's strength in relation to the SM. The neutrino propagation Hamiltonian in the presence of matter, NSI, can be expressed as 1 , m 2 and m 3 are the different mass eigenstates. V is written as: ⎤ ⎦ N e is the number density of electrons and for neutrino propagation in the Earth, G F is Fermi coupling constant, N e , N f being the number density of f fermion. The αβ are real and φ αβ = 0 for α = β. We concentrate on flavour nondiagonal NSI ( αβ 's with α = β). Here, we consider single NSI parameter eμ or eτ (one at a time) to examine the conversion probability of ν μ → ν e for the LBL studies which can be stated as the sum of three (plus higher order; cubic and beyond) terms in the presence of NSI [42][43][44]: the above Eq. (2), similar to [45] takes the following form: and, . Furthermore, here we used:

Analysis details and results
In our analysis, we used the software GLoBES [46,47] and its additional public tool [48]. The best fit values of the standard model parameters along with their corresponding uncertainties are taken from nuFIT v5. 1 [49] and PDG [50]. For example, the parameter values taken (for normal ordering) are: sin 2 θ 12 = 0.304 +0.013 −0.012 ; sin 2 θ 23 = 0.573 +0.018 −0.023 ; sin 2 θ 13 = 0.02220 +0.00068 −0.00062 ; δ C P = 194 +52 −25 ; 10 −3 eV 2 = +2.517 +0.028 −0.028 . We utilised GLoBES to combine the extracted datasets of T2K and NOνA. Using the obtained NSI constraints we discuss the sensitivity as well as the oscillation probabilities for the two next generation LBL experiments: DUNE and T2HK. We used the AEDL (a comprehensive abstract experiment definition language) files available for simulating experiments like T2HK and DUNE [51]. For our analysis purpose, we used DUNE and T2HK running for 3.5 years and 3 years in ν mode and similarly 3.5 years and 4 years inν mode, respectively.
In the case of DUNE, it will have a 40 kiloton liquid argon detector that will use a 1.2 MW proton beam to generate neutrino and antineutrino beams from in-flight pion decays. The proton beam will originate 1300 km upstream at Fermilab. The neutrino energy ranges will be between 0.5 and 20 GeV In Fig. 1 (top panel), the results of the analysis for the combination of T2K and NOνA are displayed. The left panel shows the allowed region in the plane spanned by eμ and the CP-phase δ C P , whereas the right panel displays the allowed region for eμ and the NSI phase φ eμ . For the left panel plot the non-standard CP-phase φ eμ , θ 13 , and θ 23 are marginalized away whereas for the right panel plot θ 13 , θ 23 , and δ C P are marginalized. The similar plots for IO case are displayed in Fig. 1 (bottom panel) Similarly, in Fig. 2, the left panel shows the allowed region in the plane spanned by eτ and the CP-phase δ C P , whereas the right panel displays the allowed region for eτ and the NSI phase φ eτ . For the left panel plot, the non-standard CPphase φ eτ , θ 13 , and θ 23 are marginalized away whereas for the right panel plot θ 13 , θ 23 , and δ C P are marginalized.  From the right panel of Figs. 1 and 2 we can visualize that both in NO as well as in IO cases there is a preference for a non-zero value of the coupling | eμ | and | eτ | and their corresponding phases φ eμ and φ eτ , whose values are listed out in Table 1. These values are consistent with the global constraints on neutral current NSI parameters [52]. We found δ C P value for e − μ sector around 1.12π (for NO case) as evident from the top left panel of Fig. 1. Interestingly, for e − τ sector we obtained similar value of δ C P (Fig. 2).
In Fig. 3 (top panel), we display the allowed regions in the plane spanned by the standard CP-phase δ C P and the atmospheric mixing angle θ 23 in the NO case for DUNE. The left panel refers to the SM case, while the middle and right panels concern the SM+NSI scenario with NSI arising from the e − μ and e − τ sectors, respectively. The mixing angle θ 13 and m 2 31 are marginalized away in the SM case whereas along with θ 13 and m 2 31 relevant NSI coupling ( eμ / eτ ) and non-standard CP-phase (φ eμ /φ eτ ) are marginalized in SM+NSI case. In the middle and right panels we have taken the NSI parameters with their best fit values from the combined analysis of NOνA and T2K. More specifically, | eμ | = 0.1, φ eμ = 0.2π (middle panel) and | eτ |= 0.1, φ eτ = 1.47π (right panel).
In Fig. 3 (bottom panel), similarly, we display the allowed regions in the plane spanned by the standard CP-phase δ C P and the atmospheric mixing angle θ 23 in the NO case but now for T2HK. The left panel refers to the SM case, while the middle and right panels concern the SM+NSI scenario with NSI from the e − μ and e − τ sectors, respectively. Comparing the SM scenario with that of SM+NSI arising from e − μ sector, we found distinct parameter space in the determination of θ 23 for both DUNE and T2HK. When NSI is included with SM, the allowed region corresponding to the higher octant disappears and we are left only with the allowed region from the lower octant. Whereas in SM+NSI scenario from e − τ sector, we find that both the lower as well as the higher octants are allowed for T2HK and DUNE with increased parameter space.
Concerning the θ 23 octant, we note that in the SM and SM+NSI case arising from e−μ sector there is a clear preference for lower octant for DUNE ( χ 2 = 3.69) and similarly for T2HK ( χ 2 = 0.81), where χ 2 = χ 2 SM − χ 2 SM+N SI . Corresponding one-dimensional projection plots are given in

Effect of NSI parameters on oscillation probability
In order to understand clearly the effect of NSI on LBL experiments, DUNE and T2HK, we discuss next the corresponding probability plots for both neutrino and anti-neutrino modes.
In Fig. 6 (top panel), the oscillation probability plots for DUNE in neutrino mode in the SM (left panel), SM+NSI from the e − μ sector (middle panel), and SM+NSI from the e − τ sector (right panel) are shown. We see a good separation between NO-IO for both δ C P = 90 • and δ C P = −90 • in the SM scenario. For SM+NSI scenario from the e − μ sector, we still have some separation between NO-IO for δ C P = 90 • in mid energy region, and they gradually merges around 4 GeV. Whereas δ C P = −90 • has good NO-IO separation. For SM+NSI scenario from the e−τ sector, we see a reasonable separation between NO-IO for δ C P = 90 • . In the case of δ C P = −90 • , there is some NO-IO separation in mid energy region, which gradually decreases as energy increases.
The oscillation probability plots for T2HK in neutrino mode in the SM (left panel), SM+NSI from the e − μ sector (middle panel) and SM+NSI from the e − τ sector (right panel) are shown in Fig. 7. We see a perceptible separation between NO-IO for both δ C P = 90 • and δ C P = −90 • until 1 GeV energy in the SM scenario. For the SM+NSI case from e − μ sector, we see a better separation between NO-IO for δ C P = −90 • . The NO-IO separation continuously decreases for δ C P = 90 • crossing each other around 0.7 GeV. For the SM+NSI case, from e − τ sector, we see a separation between NO-IO for δ C P = 90 • until 1.5 GeV, whereas there is no NO-IO separation for δ C P = −90 • after 0.7 GeV energy. We have repeated the exercise for anti-neutrino case in DUNE displayed in Fig. 6 (bottom panel) and find similar striking differences for NO-IO.

CP violation sensitivity
As mentioned before, one of important objectives of the current and future LBL neutrino experiments is to determine the CP phase δ C P , as precisely as possible. In the standard framework of three neutrino oscillation, we discuss here about CP violation sensitivity. The signal indicating CP violation in the lepton sector will be seen if the true values of δ C P dif-fers from the CP conserving values by a considerable amount [53]. Here, We found that for both DUNE and T2HK in Fig. 8 there is appreciable difference in the sensitivities for SM+NSI case  in comparison to SM prediction. In the case of the DUNE there appears to be better sensitivity to NSI than T2HK.

Conclusions
In this article, we assumed that new physics occurs in the form of NSI. Following that, we obtained the constraints on NSI parameters by combining the NOνA and T2K datasets.
We used the derived constraints (we have considered here mostly the case for normal ordering but checked that similar results also follow in the case of inverted ordering) and have shown that for θ 23 when we use NSI arising from the e − μ sector, both DUNE and T2HK prefer the lower octant, whereas inclusion of NSI arising from the e − τ sector brings back the degeneracy of both the lower and higher octants. Moreover, using the same set of constraints, we see striking differences in oscillation probabilities for both neutrino and anti-neutrino channels in DUNE and T2HK, which can help us understand the neutrino mass ordering problem. Furthermore, we have shown the CP discovery potential for both SM and SM+NSI scenarios and observed that the effect of NSI reduces the sensitivity, which is prominent in DUNE. Future data from NOνA and T2K will determine the fate of the existing tension in δ C P and clear the picture. If the tension persists, as we have shown in this analysis, it could probably signal the existence of new physics. Nonetheless, future studies may enable us to disentangle the NSI effects for cleaner extraction of the neutrino parameters.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.