Measuring the sterile neutrino CP phase at DUNE and T2HK

The CP phases associated with the sterile neutrino cannot be measured in the dedicated short-baseline experiments being built to test the sterile neutrino hypothesis. On the other hand, these phases can be measured in long-baseline experiments, even though the main goal of these experiments is not to test or measure sterile neutrino parameters. In particular, the sterile neutrino phase δ24\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta _{24}$$\end{document} affects the charged-current electron appearance data in long-baseline experiment. In this paper we show how well the sterile neutrino phase δ24\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta _{24}$$\end{document} can be measured by the next-generation long-baseline experiments DUNE, T2HK (and T2HKK). We also show the expected precision with which this sterile phase can be measured by combining the DUNE data with data from T2HK or T2HKK. The T2HK experiment is seen to be able to measure the sterile phase δ24\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta _{24}$$\end{document} to a reasonable precision. We also present the sensitivity of these experiments to the sterile mixing angles, both by themselves, as well as when DUNE is combined with T2HK or T2HKK.


Introduction
Neutrino oscillation physics has reached precision era. Observation of neutrino oscillations by the solar [1] and atmospheric [2] neutrinos have been confirmed independently at accelerator [3,4] and reactor [5] experiments. The two mass squared differences and all three mixing angles of the three-generation neutrino oscillation theory have been well determined and the three-generation paradigm well established [6]. The only remaining questions which are still to be answered are CP violation in the leptonic sector, neutrino mass hierarchy (that is, whether ν 3 is the lightest or the heaviest) and whether θ 23 lies in the lower or in the higher octant. Future experiments like DUNE [7][8][9][10] and T2HK a e-mail: sandhya@hri.res.in b e-mail: debajyotidutta@hri.res.in c e-mail: dipyamanpramanik@hri.res.in [11,12] are going to explore these questions and are expected to come with definite answers.
Even though neutrino oscillation with three-generations is well established, there are some hints of neutrino oscillations at a higher frequency corresponding to a mass-squared difference m 2 ∼ 1 eV 2 [13]. LSND experiment [14,15] in Los Alamos, USA, first showed evidence for such oscillations, where aν μ beam was sent to a detector and the observations showed a 3.8σ excess in the positrons, which could be explained in terms ofν μ →ν e oscillations. For the L and E applicable for the LSND experiment, this oscillation corresponds to a mass-squared difference of m 2 ∼ 1 eV 2 . Experiments like KARMEN [16] and MiniBooNE [17][18][19] tested the claim. While KARMEN data did not show any evidence forν μ →ν e oscillations, it could not rule out the entire allowed region from LSND. More recently, the Mini-BooNE experiment ran in both neutrino as well as antineutrino mode. MiniBooNE did not find any significant excess in their neutrino mode, however they reported some excess in the antineutrino mode consistent with the LSND result. Apart from these, MiniBooNE also reported some excess in the low energy bins for both neutrino and antineutrino appearance channels, but these cannot be explained in terms of neutrino flavour oscillations. The one additional sterile neutrino can be fitted along with the three active neutrinos in the so-called 3 + 1 [20] type neutrino mass spectrum. The global fit of all the relevant short-baseline data shows severe tension between the appearance and disappearance data sets and the overall goodness of fit is only 31 %. However, if one consider appearance and disappearance channels separately, the goodness of fit improves slightly and it becomes 50 % and 35 %, respectively [21].
In the 3 + 1 scenario we have 3 mass squared differences 1 m 2 41 , m 2 31 and m 2 21 , 6 mixing angles and 3 phases. ln addition to the mixing angles θ 12 , θ 23 and θ 13 which appear also in the three-generation sector we have 3 additional angles θ 14 , θ 24 and θ 34 involving the fourth generation. Also, there are two new CP phases δ 24 and δ 34 in addition to the standard CP phase δ 13 . At short-baseline experiments, the oscillation probabilities P ee , P μμ and P μe in the 3 + 1 scenario can be written in an effective two-generation framework which depends only on m 2 41 and an effective mixing term given as a combination of the sterile mixing angles θ 14 , θ 24 and θ 34 . As a result, the short-baseline experiments are completely insensitive to the sterile CP phases δ 24 and δ 34 . On the other hand, it is now well known that even though the m 2 41 -driven oscillations get averaged out, these phases show up in the oscillation probabilities at the long-baseline experiments. A lot of effort in the last couple of years has gone into estimating the impact of the sterile neutrino mixing angles and phases on the measurement of standard oscillation parameters at longbaseline experiments [22][23][24][25][26][27][28][29][30][31]. 2 These papers showed that in presence of sterile neutrino mixing the sensitivity to the measurement of CP violation, mass hierarchy, as well as octant of θ 23 becomes a band, where the width of the band comes from the uncertainty on both the values of the sterile mixing angles as well as the sterile phases. While the sterile neutrino mixing angles are constrained by the global short-baseline data, there are no constraints on the sterile phases. In the future, bounds on the sterile neutrino mixing angles are expected to improve by the data from forthcoming short-baseline experiments [54][55][56]. Studies have shown that the long-baseline experiments could also give constraints on the sterile neutrino mixing at their near [57] and far [24,58] detectors. The sterile phases on the other hand, can be constrained only in the longbaseline experiments. A short discussion on the study of sterile phases were done at T2K+reactor [29] and T2K+NOvA [28] and the sensitivity was shown to be poor. In this paper, we study how well the next-generation experiments DUNE and T2HK will be able to measure the sterile phases. We give the expected sensitivity of DUNE and T2HK alone as well as by combining data from the two experiments. To the best of our knowledge, this is the first time such a complete study is being performed. While the authors in [27] did attempt to present the expected precision on the sterile phase (which in their parametrisation was δ 14 ) in the DUNE experiment, the analysis they performed has several short-comings. In their analysis, the authors of [27] keep the sterile mixing angles θ 14 , θ 24 and θ 34 fixed in the fit at their assumed true values. We allow these angles to vary freely in the fit we perform in this paper. This allows the uncertainty due to both the mixing angles as well as the phases to impact out final results. We also keep our sterile neutrino mixing angles within the currently allowed limits, which have been updated following the results from the NEOS, MINOS and MINOS+ experiments 2 Some recent studies on other new physics scenarios in the context of DUNE and T2HK can be found here . [21]. Expected precision on the sterile phase from T2HK has never been studied before and we present them for the first time. We will also show the combined expected sensitivity of DUNE and T2HK to constrain the sterile mixing angles θ 24 and θ 14 , both when the 3 + 1 scenario is true as well as when there is no positive evidence for sterile neutrino oscillations.
The paper is organised as follows: In Sect. 2 we briefly describe the sterile neutrino hypothesis and the simulation procedure. In the same section we provide the details of the DUNE, T2HK and T2HKK experiments. In Sect. 3 we give our main results on how well the future long-baseline experiments will constrain the sterile phase δ 24 if the 3 + 1 scenario is indeed true. In Sect. 4 we present our results on the expected constraints from future long-baseline experiments if the experiments do not see any positive signal for sterile neutrino oscillations. Finally, we conclude in Sect. 5.

Sterile neutrino mixing and simulation
As discussed above, extra light neutrino state(s) with mass ∼ O(1eV 2 ) have been proposed to explain the LSND results. In this article we have considered one additional sterile neutrino within the so-called 3 + 1 scenario, where the three active neutrinos are separated from the sterile neutrinos by a mass gap of ∼ 1 eV. In this scenario the mixing matrix is 4 × 4 and hence is defined in terms of six mixing angles and three phases. The mixing matrix can be parametrised in the following way: Here O(θ i j , δ i j ) are 4×4 orthogonal matrices with associated phase δ i j in the i j sector, and R(θ i j ) are the rotation matrices in the i j sector. There are three mass-squared differences in the 3+1 scenario -the solar mass-squared difference m 2 21 7.5×10 −5 eV 2 , the atmospheric mass-squared difference m 2 31 2.5 × 10 −3 eV 2 , and the LSND masssquared difference m 2 41 1 eV 2 . One can of course write the most general oscillation probabilities in terms of these three mass-squared differences, six mixing angles and three phases. However, since oscillation driven by a given masssquared difference depends on the L/E of the experiment concerned, the expression for the oscillation probabilities simplify accordingly. In particular, the short-baseline experiments have a very small L/E such that sin 2 ( m 2 i j L/4E) ∼ 0 for m 2 21 and m 2 31 and only the terms for m 2 41 survive. This is the one-mass-scale-dominance case, where only one oscillation frequency due to one mass scale survives. As a result the short-baseline experiments depend on only "effective" sterile mixing angles, which are combinations of the mixing angles θ i j in Eq. (1). More importantly, since they have only one oscillation frequency, they do not depend on any CP violation phase. Hence, short-baseline experiments are completely insensitive to the sterile phases for the 3 + 1 scenario. 3 In long-baseline experiments such as T2HK and DUNE, the oscillations driven by m 2 31 dominate while those driven by m 2 21 are sub-dominant, while the very fast oscillations driven by m 2 41 ∼ O(1eV 2 ) get averaged out. The transition probability P μe in the limit sin 2 ( m 2 41 L/4E) ∼ 1/2 and neglecting earth matter effect is [22]: Here P 1 is the term independent of any phase, P 2 (δ 13 ) depends only on δ 13 , P 3 (δ 24 ) depends only on δ 24 and P 4 (δ 13 + δ 24 ) depends on the combination (δ 13 + δ 24 ). The full expression of the different terms in Eq. (2) are as follows: where, sin 2θ 3ν μe = sin 2θ 13 sin θ 23 b = cos θ 13 cos θ 23 sin 2θ 12 sin 2θ 4ν μe = sin 2θ 14 sin θ 24 a = cos θ 14 cos θ 24 .
We can see from Eq. (2) that even though the m 2 41 -driven oscillations are averaged out, the CP phases associated with the sterile sector still appear in the neutrino oscillation probability P μe . This dependence comes in term P 3 (δ 24 ) that depend only on the sterile phase δ 24 as well as in term P 4 (δ 13 + δ 24 ) which depends on combination of δ 13 and δ 24 . Hence, we can expect the next-generation long-baseline experiments to be sensitive to the sterile phases. We will see the anti-correlation between δ 13 and δ 24 manifest in our results on measurement of these phases in the long-baseline experiments. In fact, as has been pointed above, the sterile CP phases cannot be measured in the short-baseline experiments which are dedicated to measuring the sterile neutrino mixing. Hence, experiments like DUNE and T2HK are the only place where δ 24 can be measured in the 3 + 1 scenario. Note that in Eq. (2) the probability P μe does not depend on the mixing angle θ 34 , hence the corresponding phase associated with this angle δ 34 also does not appear. Once earth matter effects are taken into account the probability P μe picks up a θ 34 dependence and hence depends on δ 34 as well. However, for DUNE and T2HK experiments earth matter effects are rather weak and hence their corresponding sensitivity to δ 34 cannot be expected to be strong. Therefore, as we will see in the Results section, these experiments are mainly able to put constraints on δ 24 .
Prior to proceeding, we briefly discuss our simulation procedure as well as the present statues of the neutrino oscillation parameters. For our analysis we have used GLoBES (Global Long Baseline Experiment Simulator) [59,60] along with the additional codes [61,62] for calculating probabilities in the 3 + 1 scenario. We have used constant matter density for all the cases. Throughout the analysis we choose the true values 4 of the standard oscillation parameters as: 5×10 −5 eV 2 and δ 13 = − 90 • unless specified otherwise. This choice of parameters are consistent with the current limits [63]. Although one should marginalize over all the free parameters whenever one introduces some new physics, but new physics scenarios often give large number of parameters. Marginalisation over these large number of parameters is computationally challenging, so one has to do some approximation. For our analysis we have checked that the effect of marginalisation over the standard three neutrino parameters other than δ 13 have no significant effect. So to save computational time we did not marginalise over these not so relevant set of parameters.
For the sterile neutrino mixing, we have considered two scenarios. We first start by assuming that active-sterile neutrino oscillations indeed exist and find the expected constraints on the sterile neutrino mixing angles and phase δ 24 assuming non-zero sterile neutrino mixing angles in the "data". For this case we generate the "data" at true values m 2 41 ≈ 1.7 eV 2 , θ 14 ≈ 8.13 • , θ 24 ≈ 7.14 • , θ 34 = 0 • , which are the current best-fit values taken from [21]. The true values of the sterile phases will be specified. We marginalise our χ 2 over all sterile mixing parameters except m 2 41 in the fit. The χ 2 is marginalised by varying θ 14 , θ 24  Purpose of DUNE is to address all the three unknowns of the neutrino oscillation sector -the leptonic CP violation, mass-hierarchy and the octant of the θ 23 . DUNE will consist of a source facility at Fermilab and a far detector at Sanford Underground Research facility in South Dakota at a distance of 1300 km from the source. Hence, the baseline of the experiment is 1300 km. The accelerator facility at Fermilab will give a proton beam of energy 80-120 GeV at 1.2-2.4 MW which will eventually give a wide-band neutrino beam of energy range 0.5-8.0 GeV. The far site will consist of 4 identical detector of 10 kt each which will give fiducial mass of 34 kt. All the detectors will be LArTPC (Liquid Argon Time Projection Chamber).
In this work we have considered a Liquid Argon detector of fiducial mass 34 kt at a baseline of 1300 km. The neutrino flux is given by the 120 GeV, 1.2 MW proton beam. Here we have considered 5 years of neutrino and 5 years of antineu-trino. Appearance and disappearance channels are combined for the analysis. The energy resolutions for the μ and e are taken to be 20%/ √ E and 15%/ √ E, respectively. The signal efficiency is taken to be 85%. The backgrounds are taken from [8]. In the neutrino (antineutrino) mode, the signal normalization error is 2% (5%), the background normalization error is 10% (10%) and the energy calibration error is 5% (5%). This choice of systematics is conservative compared to the projected systematics in [8]. 5

T2HK
The Hyper-Kamiokande (HK) [11,12] is the upgradation of the Super-kamiokande (SK) [65] program in Japan, where the detector mass is projected to be increased by about twenty times the fiducial mass of SK. HK will consist of two 187 kt water Cherenkov detector modules, to be place near the current SK site about 295 km away from source. The detector will be 2.5 • off-axis from the J-PARC beam which is currently being used by the T2K experiment [66]. T2HK has similar physics goals as DUNE, but since it will employ a narrow-band beam, it can be complimentary to the DUNE experiment.
For our analysis we take a beam power of 1.3 MW and the 2.5 • off-axis flux. We consider a baseline of 295 km and the total fiducial mass of 374 kt (two tank each of which is 187 kt). We consider 2.5 years of neutrino and 7.5 years of anti-neutrino runs in both appearance and disappearance channels. The energy resolution is taken to be 15%/ √ E. The number of events are matched with the Tables 2 and 3 of [12]. The signal normalization error in ν e (ν e ) appearance and ν μ (ν μ ) disappearance channel are 3.2% (3.6%) and 3.9% (3.6%), respectively. The background and energy calibration errors in all channels are 10 and 5%, respectively.

T2HKK
In [12], the collaboration has also discussed the possibility of shifting one of the water tanks to a different location in Korea at a distance of about 1100 km from the source. This proposed configuration will consist of the same neutrino source but with one detector of fiducial mass 187 kt at a baseline of 295 km and another similar detector of fiducial mass 187 kt at a baseline of 1100 km. The second oscillation maximum takes place near E ν = 0.6 GeV at the second detector. Both detectors are taken 2.5 • off-axis in our study. Table 1 The 2σ allowed ranges of δ 24 (test) for the three experiments in the 3 + 1 scenario. The assumed true value of δ 13 is − 90 • . We give the allowed ranges of δ 24  Since the flux peaks at the same energy for both detector locations, the Japan detector sees the flux at the first oscillation maximum while the Korea detector sees it at the second oscillation maximum. This whole setup is called T2HKK. In our analysis, we have considered signal normalization error of 3.2% (3.6%) in ν e (ν e ) appearance channel and 3.9% (3.6%) in ν μ (ν μ ) disappearance channel, respectively. The background and energy calibration errors are 10% and 5% in all the channels, respectively.

Measurement of the sterile phases
In this section we discuss the ability of the long baseline experiments to constrain the sterile phases. The "data" is generated for the 3 +1 scenario for the values of mixing parameters discussed in Sect. 2. In particular, for the sterile neutrino parameters we take the following values: m 2 41 (true) = 1.7 eV 2 , θ 14 (true) = 8.13 • , θ 24 (true) = 7.14 • , θ 34 (true) = 0 • . The true values of the phases δ 24 will be taken at some benchmark values and will be mentioned whenever needed. The true values of standard oscillation parameters are taken at their current best-fit values, mentioned in Sect. 2. The χ 2 is marginalised over the relevant oscillation parameters in the 3 + 1 scenario, as discussed in Sect. 2, where the parameters are allowed to vary within their current 3σ ranges. Although there are three phases in the 3 + 1 scenario, the role of the phase δ 34 is weak. As was discussed in the previous section, the mixing angle θ 34 affects the oscillation probability P μe only when matter effects become important. For DUNE and T2HKK earth matter effects are not very strong while for T2HK the effect of earth matter is even weaker. Since the impact of the phase δ 34 on P μe is proportional to the mixing angle θ 34 , the phase δ 34 is also less important for P μe for the same reason. Moreover, the current global best-fit for the angle θ 34 turns out to be zero [21]. Therefore, in this work we set θ 34 (true) = 0 • in the data. As a result the phase δ 34 is not expected to be very crucial in our analysis and hence we take and δ 34 (true) = 0 • in the data and show our results only in the δ 13 -δ 24 plane. We reiterate that the χ 2 is marginalised over the mixing angle θ 34 and phase δ 34 in the fit, where the mixing angle is allowed to vary between [0 • , 12 • ] [21].
The Fig. 1 shows the capability of DUNE, T2HK and T2HKK to measure the phase δ 24 . We show the plots of  Table 1 shows that T2HK can better constrain the phase δ 24 as compared to DUNE, while T2HKK is expected to perform better than DUNE but worse than T2HK. Note that the sensitivity of DUNE and T2HKK is marginally better for δ 24 (true) = 90 • than for δ 24 (true) = − 90 • while the reverse is true in case of T2HK (see Table 1).
In order to understand why the measurement of δ 24 is expected to be better at T2HK than DUNE, we show in Fig . The upper panels of the figure reveal that the two experiments behave in almost the same way as far as the dependence of the probability P μe to δ 24 is concerned. However, there is a clear difference between the two when it comes to the overall statistics. T2HK expects to see nearly 14 times more events than DUNE due to its bigger detector size. Hence the corresponding χ 2 for T2HK is also expected to be higher. Of course the systematic uncertainty for DUNE is considerably less than for T2HK and that compensates the effect of the lower statistics, however, the effect of statistic shows up in a non-trivial way for δ 24 measurement at the long baseline experiments and T2HK with its bigger detector emerges as a better option in this regard.
The lower panels of Fig. 2 show the event rate at the oscillation maximum for DUNE (lower left panel) and T2HK (lower right panel) as a function of δ 24 . The black solid curves show the expected number of events whereas the green and yellow bands show the 1σ and 3σ statistical deviation. The black short-dashed straight lines show the event rate at oscillation maxima for the four benchmark values of δ 24 = 0 • , 180 • , 90 • and − 90 • . The Fig. 1 had revealed that the χ 2 corresponding to δ 24 (true) = 0 • and δ 24 (true) = 180 • are much lower compared to that for δ 24 (true) = ± 90 • . This can be understood from the lower panels of Fig. 2 as follows. Figure 2 shows that the predicted number of events at oscillation maximum for δ 24 (true) = 0 • and 180 • lie between the predicted events for δ 24 (true) = ± 90 • . Therefore, for the cases where data is generated for δ 24 (true) = 0 • and 180 • , it is easier for other δ 24 values to fit the data and give a smaller χ 2 . However, data corresponding to δ 24 (true) = ± 90 • takes a more extreme value and the difference between the data and fit for other values of δ 24 for these cases becomes larger, giving larger χ 2 .
Another interesting feature visible in the event plots in Fig. 2 is that the maxima and minima of the events are not at δ 24 = 90 • or −90 • . Rather they are slightly shifted towards the right. For the same reason the χ 2 plots in Fig. 1 are also asymmetric about the true value of δ 24 . One can explain this using Eq. (2). By inspecting the probability one can see that the correlation between m 2 21 and δ 24 is negligible. Also, we have taken δ 13 = − 90 • everywhere. Hence, for m 2 21 = 0 and δ 13 = − 90 • , the Eq. (2) can be rearranged as, where A and B are independent of δ 24 . In the absence of the last term, the probability would be a sine function shifted by the constant A. However, the presence of the cosine term shifts the curve and the shift is towards right because of the minus sign in front of the cosine term. In particular, the extrema of the probability in Eq. (8) is given by the condition, which corresponds to minimum at δ 24 = −63.4 • and maximum at δ 24 = 116.6 • . This agrees very well with the event plots in the lower panels of Fig. 2 which is obtained using the exact numerical probability. We next show in Fig. 3 the expected 95% CL allowed areas in the δ 13 (test)−δ 24 (test) plane, expected to be measured by the next generation long-baseline experiments T2HK (or T2HKK) and DUNE, and by combining them. The four panels of Fig. 3  In all the four cases we have considered 3+1 scenario both in the 'data' and in the 'fit' or 'theory'. The χ 2 thus generated is then marginalised over the sterile mixing angles θ 14 , θ 24 , θ 34 and δ 34 , as discussed before. The black dotted, red dash-dotted and blue dashed contours are for T2HK, T2HKK and DUNE, respectively, while the grey and magenta solid contours are for DUNE+T2HK and DUNE+T2HKK. As in Fig. 1 we note that T2HK can constrain the phase δ 24 much better than DUNE, while T2HKK performs better than DUNE but worse than T2HK. We also see, as before, that for DUNE the precision . The black dotted curve is for T2HK, the red dash-dotted curve is for T2HKK, the blue dashed curve is for DUNE, the grey solid curve is for DUNE+T2HK and the magenta solid curve is for DUNE+T2HKK on δ 24 is expected to be better for δ 24 (true)= ±90 • compared to when δ 24 (true)= 0 • or 180 • . For T2HK this dependence of precision on δ 24 measurement on δ 24 (true) is less pronounced. The effect of θ 34 on the measurement of δ 24 is also minimal. Finally, note that there is an anti-correlation between δ 13 and δ 24 . This comes from the term P 4 (δ 13 + δ 24 ) of Eq. (2).
The Fig. 3 also shows how the measurements of δ 24 and δ 13 improve as we combine DUNE with either T2HK or T2HKK. We see that combining DUNE with T2HKK improves the precision considerably, with the combined precision of DUNE and T2HKK becoming slightly better than the precision expected from T2HK alone. Combining DUNE with T2HK improves the precision even further, albeit only marginally, since T2HK alone can measure the phases rather precisely.
The question on how the measurement of the standard CP phase δ 13 gets affected by the sterile mixing angle phases in the 3 + 1 scenario is another pertinent question that one can ask. The Fig. 4 shows how the expected precision on δ 13 changes in presence of sterile neutrinos. The left panel is for T2HK, middle panel is for T2HKK and right panel is for DUNE. The blue dashed curves are for the standard 3+0 case with no sterile neutrinos while the red dash-dotted curves The other standard and sterile neutrino oscillation parameters are taken in data as described above and the fit performed as before. The Fig. 4 shows that the expected precision on δ 13 worsens when the sterile neutrino is present. From Eq. (6) one can see that there is an anti-correlation between δ 13 and δ 24 which makes the δ 13 precision worse. For DUNE the effect is more compared to T2HK and T2HKK. For T2HK the δ 13 measurement is seen to be nearly unaffected. DUNE measures δ 24 worse than T2HK and T2HKK and hence the corresponding measurement of δ 13 worsens due to the anti-correlation mentioned above. Table 2 summarises the expected precision on δ 13 for the 3 + 0 and 3 + 1 scenario for four benchmark values of δ 24 (true). We see that DUNE's measurement of δ 13 gets affected for all δ 24 while effect on T2HK's measurement of δ 13 is negligible.

Measurement of the mixing angles
Prospects of measuring the sterile neutrino mixing angles at long-baseline experiments DUNE [24] and T2HK [58] has been studied before. Here we study how well the sterile mixing can be constrained by combining data from these experiments. We also present the sensitivity of the individual experiment DUNE, T2HK and T2HKK. Here we consider two complementary approaches. We first assume that sterile neutrino mixing does exist (as in the last section) and see how precisely the data from long-baseline experiments can measure and constrain the angles θ 14 and θ 24 . 6 We next consider the alternate situation where the active-sterile oscillations do not really exist and then we see how well DUNE, T2HK and T2HKK, as well as their combination, could put upper bounds on the sterile neutrino mixing angles θ 14 and θ 24 .

Measuring the sterile mixing angles when 3 + 1 is true
In this section we assume that the 3 + 1 scenario is indeed true in nature and the mixing angles θ 14 and θ 24 are indeed non-zero. We perform a χ 2 analysis with prospective data generated in the 3 + 1 scenario and fitted within the 3 + 1 scenario and give expected allowed CL regions in the sterile neutrino parameter spaces. As before, we take the true sterile oscillation parameters at the following benchmark values: m 2 41 (true) = 1.7 eV 2 , θ 14 (true) = 8.13 • , θ 24 (true) = 7.14 • , θ 34 (true)= 0 • which are consistent with [21]. The standard oscillation parameters are taken and treated as discussed before. The χ 2 is marginalised over all relevant parameters in the fit and no Gaussian priors are included.
In Fig. 5, we show the contours in sin 2 θ 14 (test)m 2 41 (test) plane (left panel) and sin 2 θ 24 (test)m 2 41 (test) plane (right panel). The colour code is same as Fig. 3. The point where the data is generated is shown by the black star in the two panels. The results show that in both panels, T2HK gives better results than both DUNE and T2HKK. Again, T2HKK is better than DUNE. Combining T2HK/T2HKK and DUNE experiments improves the results and the expected allowed ranges for the sterile neutrino mixing parameters shrink. The precision expected from the combined data from T2HK and DUNE is nearly the same as that from T2HKK and DUNE, which for the former is only marginally better. Figure 6 shows the contours in the sin 2 θ 14 (test)-sin 2 θ 24 (test). The colour code for the different data sets considered is the Table 2 The 2σ allowed δ 13 (test) ranges for the three experiments both in 3 + 0 and 3 + 1 scenario. In both the scenarios, assumed true value of δ 13 is −90 • while in 3 + 1 case, we give the allowed ranges of δ 13

Fig. 5
The left panel shows the expected 95% C.L contours in the sin 2 θ 14 (test)m 2 41 (test) plane, while the right panel shows the 95% CL contours in the sin 2 θ 24 (test)m 2 41 (test) plane. The colour code is same as Fig. 3 Fig . 6 The expected 95% C.L contours in sin 2 θ 24 (test)-sin 2 θ 14 (test) plane. The colour code is same as Fig. 3 same as Fig. 3. We note that the expected upper limit on sin 2 θ 24 is the same for all the three individual experiments for 0.001 < sin 2 θ 14 < 0.03. Also, the expected upper bound on sin 2 θ 14 is seen to be better for T2HK than DUNE. Combining the experiments can improve the measurement of the sterile neutrino mixing angles as seen from the solid contours in Fig. 6. In particular, we now see a lower bound on sin 2 θ 24 .

Excluding the sterile hypothesis when 3 + 1 is not true
If the sterile neutrino hypothesis was wrong and there was no mixing between the active and sterile neutrinos the nextgeneration experiments would falsify it. There are a series of new short-baseline experiments planned which will be testing this hypothesis [54,55,57]. Even the near detector of planned long-baseline experiments are well-suited to check the sterile neutrino mixing as their baseline and energy match well to correspond to the maximum of m 2 41 -driven oscillations [57]. In the same vein it is pertinent to ask how well the next-generation long-baseline experiments could constrain this hypothesis, since the oscillation probabilities for long-baseline experiments also depend on the sterile neu- correspond to standard three-generation oscillation scenario with no sterile mixing while the fit is done in the 3 + 1 framework to obtain the exclusion contours. The colour code is same as Fig. 3 trino mixing and phases even though the m 2 41 -driven oscillations themselves average out. While some work in this direction has already been done in the literature [24,58], we will present here, for the first time, the sensitivity of T2HKK set-up to the sterile neutrino mixing angles θ 24 and θ 14 . We will also present the expected sensitivity from the combined prospective data-sets of T2HK (or T2HKK) and DUNE, which has not been studied before.
In Fig. 7 we show the exclusion curves for the 3 + 1 hypothesis expected from the next-generation long-baseline experiments. The left panel of Fig. 7 shows the 95% CL exclusion plots in the m 2 41 (test)-sin 2 θ 14 (test) plane while the right panel shows the results in m 2 41 (test)-sin 2 θ 24 (test) plane. Here we generate the data assuming the standard three-generation neutrino scenario and then fit it with the 3 + 1 scenario. The blue dashed, black dotted and red dashdotted curves show the exclusion plots for DUNE, T2HK and T2HKK, while the magenta and grey solid curves show the expected exclusion sensitivity for DUNE+T2HKK and DUNE+T2HK, respectively. If only three neutrinos exist in the nature, then the parameter region in the top-right of the plots are excluded at 95% CL. Again, as in the previous results, T2HK constrains θ 14  is small and lies in the range 10 −3 eV 2 < m 2 41 < 0.01 eV 2 , T2HK gives better constraint on sin 2 θ 24 than both DUNE and T2HKK. But for higher values of m 2 41 , in the range of 0.014 eV 2 < m 2 41 < 0.1 eV 2 , the performance of DUNE is better than both T2HK and T2HKK. For 0.1 eV 2 < m 2 41 < 10.0 eV 2 , performance of T2HKK is almost similar to that of DUNE. Similar behaviour can be seen in the combined case. In the lower m 2 41 region, DUNE+T2HK could constrain sin 2 θ 24 slightly better than DUNE+T2HKK. But in the higher m 2 41 region, DUNE+T2HKK is expected to perform better than DUNE+T2HK. The expected exclusion sensitivity for DUNE, T2HKK and T2HK in the range 0.1 eV 2 < m 2 41 < 10.0 eV 2 are given as sin 2 θ 24 ∼ > 0.026, sin 2 θ 24 ∼ > 0.026 and sin 2 θ 24 ∼ > 0.03, at 95% CL Similarly, the expected exclusion limit for DUNE+T2HKK and DUNE+T2HK at 95% CL are sin 2 θ 24 ∼ > 0.017 and sin 2 θ 24 ∼ > 0.019, respectively, for 0.1 eV 2 < m 2 41 < 10.0 eV 2 .
In Fig. 8 we show the expected exclusion contour in sin 2 θ 14 (test)-sin 2 θ 24 (test) plane. Here, the region outside the contour is excluded at 95% CL. The figure represents the slice at m 2 41 = 1.7 eV 2 of the contour in the sin 2 14 , Fig. 8 The expected 95% C.L exclusion curves in the sin 2 θ 14 (test)sin 2 θ 24 (test) plane for a fixed m 2 41 = 1.0 eV 2 . The colour code is same as Fig. 3 sin 2 θ 24 , m 2 41 space. The colour code is the same as in Fig. 7. Here also, we observe better capability of T2HK to constrain sin 2 θ 14 -sin 2 θ 24 parameter space than DUNE and T2HKK in most regions of the parameter space. The plot also shows that constraint on sin 2 θ 24 is complicated. We see that T2HK is better than DUNE and T2HKK in constraining sin 2 θ 24 for sin 2 θ 14 ∼ > 10 −2 . However, for sin 2 θ 14 ∼ < 10 −2 DUNE and T2HKK perform better than T2HK in constraining sin 2 θ 24 . Combining the data-sets improves the expected sensitivity on both the sterile mixing angles.

Conclusions
There are a number of observational hints that support the existence of neutrino oscillations at short baselines. Since the m 2 needed for these frequencies is inconsistent with the m 2 needed to explain the solar and atmospheric neutrino anomalies -both of which have been confirmed by earth-based experiments -it has been postulated that there are additional light neutrino states which are mixed with the three standard neutrino states. Since the Z -decay width restricts the number of light neutrino states coupled to the Z boson to 3, the additional neutrino states should be "sterile". In this paper we considered one extra such sterile neutrino in the so-called 3 + 1 mass spectrum. In the 3 + 1 scenario the neutrino oscillation parameter space is extended by one new mass squared difference m 2 41 , three new active-sterile mixing angles θ 14 , θ 24 and θ 34 and two new CP phases δ 24 and δ 34 . We work within a parametrisation of the mixing matrix such that the phase δ 24 is associated with the mixing angle θ 24 and δ 34 is associated with θ 34 . It is now well known that even though the m 2 41 -driven oscillations are averaged out in the longbaseline experiments, the active-sterile mixing angles and the additional phases appear in the oscillation probabilities and modify it. The sensitivity of the long-baseline experiments to the active-sterile mixing angles has been studied before. The impact of the sterile neutrino parameters on the physics reach of these experiments for standard parameter measurement such as CP violation, mass hierarchy measurement and octant of θ 23 measurement has been investigated in details before. In this work, for the first time, we looked at the prospects of measuring the sterile CP phase δ 24 in the long-baseline experiment T2HK (and T2HKK) and DUNE as well as when data from them is combined.
Dedicated short-baseline experiments are being built to test the active-sterile neutrino oscillation hypothesis. However, these experiments are sensitive to oscillation probabilities that have one-mass-scale-dominance. In other words, these experiments mainly work within effective twogeneration scenarios and are sensitive to m 2 41 and effective two-generation mixing angles, that can be written as combination of the mixing angles θ 14 and θ 24 . Therefore, they are completely insensitive to the CP phases δ 24 and δ 34 . On the other hand, these phases do affect the oscillation probabilities of the long-baseline experiments even though they are not sensitive to m 2 41 since oscillations corresponding to this frequency averages out. In particular, the probability P μe that affects the electron appearance data in long-baseline experiments depends on δ 24 . We exploited this dependence to show that the phase δ 24 can be measured at DUNE, T2HK and T2HKK and estimated the expected precision on this parameter. Since P μe depends on θ 34 only through earth matter effects, the effect of δ 34 on the long-baseline data is expected to be very small and we cannot constrain it easily. Hence, in this paper we concentrated only on δ 24 . We performed a χ 2 analysis of the prospective data at these experiments and presented the expected precision on δ 24 expected at these experiments. We showed that with (7.5 + 2.5) years of running in (neutrino, antineutrino) mode, T2HK could constrain δ 24  are allowed in T2HK. The corresponding constraints from T2HKK and DUNE were seen to be weaker, with DUNE measurement on δ 24 expected to be the weakest. We explained why T2HK is expected to perform better than DUNE, with the main reason being the higher statistics in T2HK. We also showed the expected allowed areas in the δ 13 (test)-δ 24 (test) plane from prospective data from DUNE, T2HK and T2HKK. We again reiterated that the expected constraints on δ 24 were seen to be strongest from the T2HK experiment, while DUNE was seen to be weakest. We also presented the allowed areas in this plane expected from combined prospective data of DUNE and T2HK (or T2HKK). We showed that the combined data set could constrain the CP parameters better. Expected constraints from DUNE+T2HK was seen to be better than constraints from DUNE+T2HKK. We also discussed the impact of the sterile neutrino mixing angles and phases on the measurement of the standard CP phase δ 13 . We showed that for DUNE the expected δ 13 precision worsens more than for T2HK.
We also presented constraints on the mixing angles θ 14 and θ 24 . Again, we did not consider θ 34 since the appearance and disappearance data in long-baseline experiments depend only mildly on this angle, with the dependence coming solely from matter effects. We took two complementary approaches in this study. First we assumed the 3 + 1 scenario to be correct and generated data assuming non-zero values of the sterile neutrino oscillation parameters. This was used to present allowed areas in the sterile neutrino parameter space expected from full run of DUNE, T2HK, T2HKK and combinations of DUNE+T2HK and DUNE+T2HKK. We made a comparison between the expected precision reach of the different experiments. We next took the complementary approach where we assumed that the 3 + 1 scenario was not true. The data in this case was generated for no sterile mixing and fitted with the 3 + 1 hypothesis to yield expected exclusion limits on the sterile neutrino mixing parameters. Again, we did this analysis for DUNE, T2HK, T2HKK and combinations of DUNE+T2HK and DUNE+T2HKK and made a comparative analysis of the different data-sets.
In conclusion, the sterile neutrino phases can be measured only in experiments that are sensitive to more than one oscillation frequency other than m 2 41 . The long-baseline experiments therefore are the best place to measure δ 24 . We showed that the sterile phase can be measured to reasonable precision in the next-generation long-baseline experiments. The T2HK set-up is better suited to measure δ 24 compared to DUNE due to larger statistics.