Exploring fake solutions in the sterile neutrino sector at long-baseline experiments

Active-sterile neutrino mixing is known to affect the neutrino oscillation probabilities at both short as well as long-baselines. In particular, constraints on active-sterile neutrino oscillation parameters can be obtained from long-baseline experiments such as T2HK and DUNE. We present here existence of fake solution in the appearance channel for the 3+1 scenario at long-baseline experiments. We show that the appearance probability is same for values of $\Delta m_{41}^2$ for which the fast oscillations are averaged out and for $\Delta m_{41}^2=(1/2)\Delta m_{31}^2$. The fake solution does not appear for the disappearance channel.


I. INTRODUCTION
Neutrino oscillation physics has entered the precision era. Now we know beyond all doubts that there are at least two massive neutrino mass eigenstates and three non-zero mixing angles, θ 12 , θ 13 and θ 23 . We also have achieved reasonable precision in the measurement of the two mass squared difference and the mixing angles θ 12 and θ 13 . While we have measured θ 23 , some work remains to be done in this sector. In particular, although there are hints from recent experiments like NOvA [1] and T2K [2] that θ 23 is close to maximal, it is still ambiguous. Also, if it was to be non-maximal, we do not know the octant of θ 23 , .i.e., whether it is less than 45 • or greater than 45 • . While there are hints of δ CP = −90 • from the current experiments, we still need to confirm that there is CP-violation in the leptonic sector. Finally, we still need to determine the sign of ∆m 2 31 , aka, the neutrino mass hierarchy. Dedicated future experiments such as INO [3], JUNO [4,5], DUNE [6][7][8][9] and T2HK [10,11] are being planned to resolve these issues.
Another unsettled issue is the question regarding the existence of a light sterile neutrino mixed with the active ones [12]. Hints of this came when LSND [13,14] in USA reported 3.8σ excess in their positron events, pointing atν µ →ν e oscillations. This excess could not be explained within the standard three-generation paradigm because the LSND L/E demanded a ∆m 2 ∼ eV 2 while, the solar and atmospheric neutrino data demand ∆m 2 ∼ 10 −5 eV 2 and 10 −3 eV 2 , respectively. Therefore, one has to postulate the existence of a fourth light neutrino such that it corresponds to a mass-squared difference ∆m 2 ∼ 1 eV 2 . Since the bound on number of light neutrinos coupled to the Z-boson is 2.9840 ± 0.0082 from LEP [15], the fourth generation must be sterile. Existence of these oscillations were tested at the short-baseline experiments KARMEN [16] and MiniBooNE [17][18][19] which looked forν µ →ν e appearance, as well as at CDHS [20], MINOS [21], MiniBooNE (disappearance search) [22] and SuperKamiokande atmospheric experiment [23] which looked at ν µ (ν µ ) disappearance.
Data from neutrinos with longer baseline such as MINOS [24], IceCube [25] and NOvA [26] have also reported null signal for ν µ → ν e oscillations at the LSND scale. While KARMEN did not find any evidence for neutrino oscillations consistent with LSND, it was unable to eliminate the entire parameter space favored by LSND. None of the disappearance searches observed any evidence for neutrino oscillations consistent with LSND, making the LSND signal anomalous, which reflects in a very low goodness-of-fit reported from global analyses of all experiments [27]. On the other hand, the MiniBooNE experiment designed specifically to test the LSND anomaly, ran in both neutrino and antineutrino mode and the collaboration have recently presented the complete analysis of their entire data set [28]. They reported an observed 4.8σ excess of electron-like events at MiniBooNE. When combined with LSND this constitutes a 6.1σ excess, putting the "LSND anomaly" at a statistically compelling level.
Hints for active-sterile neutrino oscillations have also come from the ν e andν e sectors.
While the Gallium experiments [29][30][31][32][33] require oscillations of ν e with ∆m 2 ∼ 1 eV 2 , there appear to be inconsistency between the reactor flux calculations and reactor data [34] which favors existence of a fourth neutrino. Confirmation for active-sterile mixing drivenν e disappearance searches have been performed by NEOS [35], DANSS [36] and STEREO [37] and so-far there is no convincing evidence for active-sterile mixing.
The impact of active-sterile mixing in long-baseline experiments have been studied extensively in the literature. The effect of active-sterile mixing on measurement of standard three-generation neutrino oscillation parameters at long-baseline experiments NOvA, T2HK and DUNE was studied in detail in [38][39][40][41][42][43][44][45]. On the other hand, possibility of measuring the active-sterile mixing parameters was explored extensively for DUNE near detector in [46] and NOvA, T2HK and DUNE far detectors in [47][48][49][50][51]. In this paper we point out a new fake solution that exists in the appearance channel at long-baseline experiments. We will show both analytically as well numerically that the appearance probability is the same for values of ∆m 2 41 for which the fast oscillations are averaged out and for ∆m 2 41 = (1/2)∆m 2 31 . The fake solution does not appear for the disappearance channel.
The paper is organised as follows: In section II we will show the existence of the fake solution analytically and numerically. In section III we will present the main χ 2 results of the paper. Finally, we will conclude in section IV.

II. FAKE SOLUTION IN THE APPEARANCE PROBABILITY
We will consider one additional sterile neutrino, a scenario referred to 3+1 [52] in the literature. Here we have a 4×4 mixing matrix described by 6 mixing angles: θ 12 , θ 13 , θ 23 , θ 14 , θ 24 and θ 34 and 3 phases: δ 13 , δ 24 and δ 34 . We use the parametrisation convention where the mixing matrix is given as: where O(θ ij , δ ij ) are 4×4 orthogonal matrices with phase δ ij associated with the ij sector, and R(θ ij ) are the rotation matrix with the ij sector. In the 3+1 scenario, in addition to ∆m 2 21 and ∆m 2 31 , there is a third mass-squared difference, ∆m 2 41 . For the long-baseline experiments, fast oscillations due this third mass squared difference would be averaged out for values of ∆m 2 41 ∼ O(1eV 2 ) due to the finite resolution of the detectors. Hence, for these values of ∆m 2 41 ∼1 eV 2 the ν µ → ν e oscillation probability can be written as [49], where 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 ). Note that while we will numerically calculate the exact probability in presence of earth matter effect in all our results presented in the next section, here for simplicity we give the analytic expressions for vacuum only. We stress that the same expressions will be valid for neutrinos propagating in constant density matter, just by replacing the oscillation parameters θ 13 and ∆m 2 31 in vacuum with those in constant density matter. We will discuss towards the end of this section the role of matter in the fake solution. The full expression of the different terms in Eq. (2) are as follows: P 4 (δ 13 + δ 24 ) = X sin 2θ 3ν µe sin 2θ 4ν µe cos(δ 13 + δ 24 ) − 1 2 sin 2 2θ 12 cos 2 θ 13 sin 2 ∆ 21 + cos 2θ 13 (cos 2 θ 12 sin 2 ∆ 31 + sin 2 θ 12 sin 2 ∆ 32 ) + 1 2 sin(δ 13 + δ 24 ) cos 2 θ 12 sin 2∆ 31 + sin 2 θ 12 sin 2∆ 32 , where, sin 2θ 3ν µe = sin 2θ 13 sin θ 23 (7) sin 2θ 4ν µe = sin 2θ 14 sin θ 24 (8) and, If we put the approximation ∆m 2 21 = 0 and the condition δ 13 + δ 24 = 0 in Eq. (3-6), we get, On the other hand, for values of ∆m 2 41 <<1 eV 2 , the oscillations due to this mass scale will survive the detector resolutions and show-up at the long-baseline detector. If we continue using the approximation ∆m 2 21 = 0 while allowing the ∆m 2 41 -driven oscillatory terms, the expression for the probability becomes, If we use the condition δ 24 + δ 13 = 0 in Eq. (13) we get, Now, for the particular choice ∆m 2 41 = 1 2 ∆m 2 31 the oscillation probability reduces to, We can see from Eq. (12) and Eq. (15) that the oscillation probability for these two cases are equal, except for the first term. But at the oscillation maximum, ∆ 31 = π 2 , the first term of the Eq. (15) becomes 1 2 sin 2 2θ 4ν µe . So at the oscillation maximum the oscillation probability for the ∆m 2 41 >1 eV 2 case given by Eq. (12) becomes exactly equal to the oscillation probability for the ∆m 2 41 = 1 2 ∆m 2 31 case given by Eq. (15). This produces a fake solution in ∆m 2 41 at the energy corresponding to ∆m 2 31 oscillation maximum for a given baseline L. Note that the above condition for the fake solution was obtained for oscillations in vacuum. As mentioned before, using the corresponding expression in constant density matter we would get the fake solution in matter at the energy that gives the oscillation maximum for ∆m 2 31 m for a given baseline L. One can check that for earth matter effect in experiments like T2HK and DUNE the energy at which we get oscillation maximum for ∆m 2 31 m is very similar to that for ∆m 2 31 . We will quantify the effect of matter on the fake solution shortly in Fig. 1.
For energies away from the ∆m 2 31 -driven first oscillation maximum, the oscillatory term sin 2 ( ∆ 31 2 ) = 1/2 exactly and hence, the above mentioned fake solution is not exact. However, we see that the fake solution is achieved even outside the ∆m 2 31 -driven first oscillation maximum, albeit approximately. To explore this further, let P ∆m 2 31 /2 be the oscillation probability for ∆m 2 41 = ∆m 2 31 /2 and P ave for ∆m 2 41 ∼ 1eV 2 , and let ∆P be the difference between P ∆m 2 31 /2 and P ave : The term sin 2 ∆m 2

L 8E
takes values only in the range [0,1]. Hence, ∆P is also bounded, and lies between −0.5 sin 2 (2θ 4ν µe ) ≤ ∆P ≤ 0.5 sin 2 (2θ 4ν µe ). Therefore, for sin 2 (2θ 4ν µe ) small, ∆P max is also expected to be small. For example, if sin 2 (2θ 4ν µe ) ∼ 10 −3 , then ∆P max ∼ O(10 −3 ), whereas the oscillation probability P µe ∼ O(10 −2 ), so we can say that there is an approximate degeneracy for all energies for small values of sin 2 (2θ 4ν µe ). In our simulations we have used the value of sin 2 (2θ 4ν µe ) that lies within the currently allowed range obtained from global analysis of all short baseline data including LSND and MiniBooNE [27] and corresponds to 10 −3 in our case. can be easily seen. The probabilities for the two cases mentioned above are almost same. Very small differences are seen at energies away from the oscillation maximum, for reasons discussed above. Due to the finite resolution of the detector, this tiny difference is not expected to be differentiated by the detector and hence we expect to see two degenerate solutions for two values of ∆m 2 41 . A final comment is in order. Fake solutions in the active-sterile sector have been studied in the context of Daya Bay [53], RENO and Double Chooz [54] experiments. The Daya Bay paper [53] discusses the degeneracy between sin 2 θ 13 and sin 2 θ 14 for the case when ∆m 2 31 ≈ ∆m 2 41 . In [54] the degeneracy between sin 2 θ 13 and sin 2 θ 14 appears in the Pēē channel in the 3+2 scenario when ∆m 2 41 and ∆m 2 51 is taken to be in the 1 eV 2 regime. In this paper we will look into the degeneracy in the P eµ channel in the context of acceleratorbased long-baseline experiments. The muon neutrino disappearance is seen to not have any active-sterile degeneracy.

III. NUMERICAL RESULTS
We have used GLoBES [55,56] and the necessary code [57,58] 2 31 in the data. In all panels the data have been generated for ∆m 2 41 = 1.6 eV 2 . From the left panels we can see that when δ 13 = −90 • and δ 24 = 90 • , we get two (nearly degenerate) solutions in ∆m 2 41 . One solution corresponds to all values ∆m 2 41 ∼ > 0.01 eV 2 for which ∆m 2 41 -driven oscillations are averaged out at the far detector of the experiment. This is the true solution since for all these values of ∆m 2 41 the fast oscillations are averaged out as for ∆m 2 41 = 1.6 eV 2 in the data. In addition, there is a fake solution at ∆m 2 41 ≃ (1/2)∆m 2 31 at which the value of χ 2 is nearly zero. This degeneracy corresponding to (δ 24 + δ 13 ) = 0 was explained in the previous section II in terms of the probabilities, both analytically as well as numerically. The fact that the χ 2 at the fake solution is slightly different from zero comes mainly due to two reasons -(i) in this plot we have taken ∆m 2 21 non-zero, and (ii) as discussed before, the χ 2 s have been obtained using the full energy range of T2HK and DUNE, while the exact degeneracy exists only at the oscillation maximum.
Note from the figure that for all other combinations of δ 13 + δ 24 , the χ 2 at ∆m 2 41 ≃ (1/2)∆m 2 31 is far from zero. However, for both the other combinations of δ 13 + δ 24 plotted in the figure, there are dips in the value of χ 2 at values of ∆m 2 41 other than the true values corresponding to the data. Therefore, one can expect to see islands of allowed ∆m 2 41 at values close to (1/2)∆m 2 31 for all combinations of δ 13 + δ 24 . While for δ 13 + δ 24 = 0, this fake island is expected to appear as a nearly degenerate solution, for all other combinations it would appear only at a higher C.L. Note also that though for T2HK there is dip at only one fake value of ∆m 2 41 , for DUNE there are multiple dips. This is because DUNE is a wide-band beam experiment while T2HK is a narrow-beam experiment.
The fake solution is further explored in the right panels of Fig. 2. Here the data is generated for only the case δ 13 +δ 24 = 0 with δ 13 = −90 • and δ 24 = 90 • but for different values of ∆m 2 31 . The red solid, blue dashed and the black dotted curves are for ∆m 2 31 = 2.5 × 10 −3 eV 2 , ∆m 2 31 = 1.0 × 10 −3 eV 2 and ∆m 2 31 = 5.0 × 10 −3 eV 2 , respectively. As the δ 13 + δ 24 = 0, there are degeneracies present in all of the three curves. And it is interesting to note that the fake solution is always at half of the value ∆m 2 31 . In Fig. 2 we saw that the fake solution of ∆m 2 41 always appears at half of ∆m 2 31 . In fits the red points very well, so indeed the fake solution is at half of the value ∆m 2 31 . The Fig. 4 shows the effect of combining the data from a narrow band beam experiment (T2HK) with that from a wide band beam experiment (DUNE). The black dotted curve shows the χ 2 for the narrow band experiment, the blue dashed curve shows the χ 2 for the wide band experiment, and the grey solid curve shows the combination of them. We can see from Fig. 4 that even though the grey curve has a higher χ 2 than the blue and black curves for some values of ∆m 2 41 , but the fake solution is still present. This is expected because the fake solution is seen to be present in both the narrow as well as the wide band beam experiments at the same value of ∆m 2 41 . Therefore, the combination of these two types of experiments can not remove the fake solution. Note however that while the fake solution survives, its range is reduced if we combine the narrow and wide band beam experiments because of increased total statistics. The Fig. 5 shows the effect of the near detector. The red solid curve gives the χ 2 for the DUNE far detector, the magenta dashed dotted curve is for the near detector and the blue dashed curve is for both near and far detector together. This figure shows the importance of the near detector. We can see from the figure that the fake solution completely goes away as soon as we add the near detector data in the analysis. This is expected since we do not expect the degeneracy in the detector. As a result, due to the much higher statistics in the near detector, the fake solution coming from the far detector data completely disappears.
In all figures shown so-far, we have only considered the appearance channel in the longbaseline experiments. This is because the fake solution presented in this paper comes only in the appearance channel. The disappearance channel does not have this fake solution and hence should water down the emergence of the fake ∆m 2 41 island. In Fig. 6 we show the impact of adding the disappearance data into the analysis. The blue curve shows the χ 2 as a function of ∆m 2 41 for the appearance channel only, while the red curve shows the results when data from both appearance and disappearance channels are combined to obtain the χ 2 . We can see from the figure that introduction of disappearance channels removes the fake solution. The left panel is for T2HK and the right panel is for DUNE. The impact of the disappearance channel in lifting the degeneracy is visible in both experiments.

IV. CONCLUSIONS
Although the three neutrino oscillation paradigm is well established, there are tantalising hints for the possible existence of sterile neutrino state(s). While the question on whether sterile neutrinos driving such short-baseline neutrino flavor oscillations exist or not is far from settled, the presence of active-sterile mixing has been shown to impact neutrino oscillation at long-baseline experiments. Detailed analyses have been performed in the literature looking at the potential of long-baseline experiments in measuring the active-sterile oscillation parameters. In this paper we showed, for the first time, the existence of fake solution in the P µe appearance probability, coming from ∆m 2 41 . We showed that if the condition δ 13 + δ 24 = 0 is satisfied, fake solution is obtained in P µe at the oscillation maximum for ∆m 2 41 = (1/2)∆m 2 31 in addition to the true solution. We showed this both numerically as well as analytically. It was shown from a comparison of the numerical and analytical solutions that ∆m 2 21 does not affect the position of the ∆m 2 41 fake solution. We discussed the breaking of the degeneracy for values of energy away from the oscillation maximum and showed that this breaking is very mild. The impact of earth matter effects on the degeneracy was also discussed and shown to be a very mild shift in the energy of oscillation maximum.
We next showed the appearance of the degenerate ∆m 2 41 solution from a χ 2 analysis of the prospective data at both T2HK and DUNE. The degeneracy was shown to appear at ∆m 2 41 = (1/2)∆m 2 31 , while other degenerate solutions in ∆m 2 41 appeared at higher C.L. It was also seen that DUNE had multiple fake solutions at different ∆m 2 41 , in addition to the degenerate solution. Finally, we discussed that the degeneracy between large ∆m 2 41 and ∆m 2 41 = (1/2)∆m 2 31 does not exist for the disappearance channel. As a result combining of the appearance and disappearance data resulted in erasing of the fake solution at ∆m 2 41 = (1/2)∆m 2 31 . In conclusion, an fake solution in ∆m 2 41 exists in the appearance probability at longbaseline experiments. This would result in one or more fake islands appearing in the ∆m 2 41 allowed C.L. contours coming from the appearance channel of long-baseline experiments.
The disappearance channel does not have any such fake solution.