Abstract

The future long baseline experiments such as DUNE and T2HKK have promising prospects to determine the neutrino mass hierarchy and measuring standard phase . However, presence of possible nonstandard interactions of neutrinos with matter may intricate this picture and is the subject matter of the present work. We have studied the standard parameter degeneracies in presence of nonstandard interactions (NSI) with DUNE and T2HKK experiments. We examine the mass hierarchy degeneracy assuming (i) all NSI parameters to be nonzero and (ii) one NSI parameter () and its corresponding phase () to be nonzero. We find that the latter case is more appropriate to resolve mass hierarchy degeneracy with DUNE and T2HKK experiments due to relatively small uncertainties emanating from the NSI sector. We have, also, investigated the octant degeneracy with neutrino () and antineutrino () mode separately. We find that to resolve this degeneracy the long baseline experiment with combination of neutrino and antineutrino mode is essential. Furthermore, we have considered DUNE in conjunction with T2HKK experiment to study phase degeneracy due to standard () and nonstandard () phases. We find that DUNE and T2HKK, in conjunction, have more sensitivity for violation effects (10 for true NH and 8.2 for true IH).

1. Introduction

The discovery of nonzero neutrino masses and lepton flavor mixing by the reactor [1], accelerator [2], atmospheric [3], and solar [4] neutrino oscillation experiments have revealed the values of oscillation parameters such as mass squared differences and mixing angles [5], to an unprecedented accuracy. At present, there are some unknown quantities in standard three-neutrino framework, namely, sign of , the octant of , and the -phase , the determination of which is the prime objective of current and future neutrino oscillation experiments. The difficulty in the determination of these unknowns is the existence of degeneracies in neutrino oscillation parameters. To overcome these degeneracies, one of the methods is to combine data from different neutrino oscillation experiments. Recently, this procedure has been adopted by various studies [68], where the synergy between current and future experiments has been considered. In principle, future neutrino oscillation experiments have sensitivity reach to perform precision test of standard neutrino oscillation paradigm and to probe new physics beyond standard model (SM). In neutrino oscillation experiments, one model-independent way to study new physics (NP) is given by the framework of nonstandard interactions (NSI) [9, 10].

An alternative phenomenon to explain neutrino flavor transitions, on the basis of NSI, was first proposed by Wolfenstein [11]. Although we know that they will show their effect in neutrino oscillation experiments at subleading level, they are important, with the emergence of next generation experiments like Tokai-to-Hyper-Kamiokande (T2HK) [12], Deep Underground Neutrino Experiment (DUNE) [13], Tokai-to-Hyper-Kamiokande-and-Korea (T2HKK) [14], etc., where such type of interactions can be probed. In general, the NSIs may manifest itself in propagation of neutrino through matter and the processes involved in its creation and detection. The possible manifestations of NSIs have been widely studied in the literature and bounds on NSI parameters have been derived from various experiments [9, 10, 15, 16]. Furthermore, the model-independent bounds on NSI in production and detection regions are an order of magnitude stronger than the matter NSI [15]. In this work, we focus on matter NSI which can be defined by dimension-six four-fermion operators given by [11, 17]where , , , and are dimensionless parameters indicating the strength of the new interaction having units of . To probe matter NSI long baseline neutrino experiments (LBNE) are ideal and the neutral current (NC) interactions which affect the neutrino propagation coherently can also be studied at far detectors. The next generation LBNE such as DUNE, T2HK, and T2HKK may reach the sensitivity to reveal NSI in neutrino sector.

In the leptonic sector, the violation can render leptogenesis mechanism which in turn may shed light on baryogenesis [19, 20]. It is very difficult to measure leptonic -violation in presence of NSI as it will get bewildered by the existence of possible -violation generated by NSI itself. Undoubtedly, the existence of NSI has opened an entirely new window to explore NP beyond standard model.

Previously, the authors of [6] have explored the ability to disentangle the violating effects due to standard and nonstandard contributions under the assumption that only one NSI parameter or is present. In [7], the parameter degeneracies in LBNE originating from nonstandard interactions have been studied and [21] has focused on NSI at DUNE, T2HK, and T2HKK and has concluded that overall DUNE has the best sensitivity to the magnitude of the NSI parameters, while T2HKK has the best violation sensitivity with or without NSI. Furthermore, in [22], the authors have studied the impact of nonzero NSI on the precision of DUNE. The authors of [2325] have explored the effects of NSI on violation sensitivity and hierarchy sensitivity at DUNE, respectively. In [8] the authors have studied the sensitivity to mass hierarchy, the octant of , and phase in the future long baseline experiments T2HK and DUNE assuming standard interactions (SI) only. In general, earlier studies on standard parameters degeneracies with SI or matter NSI have mostly focused on DUNE [2326]. Motivated by the long baseline of T2HKK experiment, it is imperative to study physics potential of T2HKK and DUNE+T2HKK, in resolving standard parameter degeneracies in presence of NSI. In the present work, we have investigated prospects for lifting mass hierarchy degeneracy (sign degeneracy), -octant degeneracy, and -phase degeneracy in DUNE, T2HKK, and DUNE+T2HKK with matter NSIs.

T2HKK is a long baseline experiment proposed to enhance the hierarchy sensitivity of T2HK by setting one of the two tanks of HK detector at a site in Korea. This multidetector setup is advantageous as it gives access to a longer baseline of 1100 km and simultaneously boosts the data at the T2HK with baseline of 295 km [27]. The neutrino oscillation probabilities are strongly affected by the matter effects in long baseline experiments. These matter effects can be beneficial in lifting up the standard parameter degeneracies. Therefore, we have considered the T2HKK setup with larger baseline of 1100 km in the analysis. In present work, we have studied the standard parameter degeneracies, i.e., mass hierarchy degeneracy and octant degeneracy in presence of matter NSI with DUNE and T2HKK experiment. We have assumed all NSI parameters to be nonzero in one case and only one off-diagonal NSI parameter to be nonzero, in another case. We find that the latter case is better at resolving standard parameter degeneracies in case of both DUNE and T2HKK experiments. Due to the larger baseline, T2HKK is found to have similar sensitivity as DUNE experiment to resolve standard parameter degeneracies including NSI. Furthermore, we have investigated the phase degeneracy occurring due to the contribution from standard and nonstandard phases. We observe that it is difficult to disentangle the effects due to SI phase from NSI phase at DUNE+T2HKK experiment as this conjunction is more sensitive to study violation effects [27].

We organize the paper as follows: in Section 2, we present the formalism to write oscillation probability in presence of matter NSIs. We discuss about the long baseline experiments DUNE, T2HKK, and corresponding simulation details in Section 3. In Section 4, we discuss the prospects to resolve standard parameter degeneracies in these LBNEs. We have presented our results and, subsequent, discussion in Section 5. Finally, we conclude in Section 6.

2. Formalism: Oscillation Probabilities

The Hamiltonian for the neutrino propagation in presence of matter NSI can be written aswhere is the PMNS mixing matrix containing three mixing angles and one phase , . is the matter potential due to interaction of neutrino with matter, viz.,where . The unit contribution in the first element of the matrix is due to the matter term contribution from standard charged-current interactions. The diagonal element of is real, i.e., (where ) for and . The oscillation probability for channel can be written as [7]where , . Similar expression can be obtained for inverted hierarchy (IH) by replacing (i.e., and ). The expression for antineutrino oscillation probability can be written by replacing and in (4).

3. Experimental Setups

Considering the sensitivity reach of the present and future long baseline neutrino oscillation experiments (for example, DUNE and T2HKK), it is very important to study the individual and collective effects of NSI parameters on parameter degeneracies. We have used GLoBES package [28, 29] with best-fit values and ranges of the standard neutrino mixing parameters, as given in [30] to simulate the DUNE and T2HKK. The current bounds on NSI parameters used in present analysis are [15]. The phases of the off-diagonal NSI parameters are still unconstrained and can lie in the range .

The experimental configurations, energy resolutions, and systematic uncertainties considered in the present work are as follows.

3.1. DUNE

The DUNE experiment [13], situated in the USA, is a globally synchronized endeavor of neutrino physicists around the world. Out of many others, the neutrino physics goals of the experiment are to unravel the sign of neutrino mass hierarchy and to measure the phase(s). The experiment is planned to direct neutrino beam from Fermilab to Homestake mine in South Dakota providing an optimum baseline of 1300 km for manifestation of matter effects in neutrino oscillations. Unlike Hyper-K, DUNE is an on-axis experiment. We have used DUNE CDR [13, 31] with 35 kt LAr far detector. The optimized beam design that employs 80 GeV beam of protons having 1.0 MW power has been used to simulate the experiment. We have considered 5(+5) years of run in neutrino (antineutrino) mode resulting in an exposure of 350 kt.MW.years. The appearance efficiency () and energy resolutions () taken in the present analysis are 80% (, ), respectively. The normalization and energy calibration uncertainty for signal ()/background () is taken to be , %, %, and %. For signal ()/background () the values are %, %, %, and %.

3.2. T2HKK

The T2HKK experiment [14], an extension T2HK [12], is proposed to be stretched over Japan and Korea. The neutrino beam will be directed from J-PARC facility in Japan to two water-Cherenkov detectors: (i) first detector at Kamioka mine in Japan with a baseline of 295 km; (ii) second detector to be built in Korea providing a baseline of 1100 km. In the present work, we have considered 1100 km baseline (also, referred to as T2HKK), with detector at off-axis with respect to the neutrino beam, where matter effects will be large. We choose 13 MW.years beam power which is similar to that of T2HK. The running time, in ratio 1:3 for neutrino and antineutrino mode, is 10 years amounting to total exposure of protons on target (POT). The appearance efficiency () and energy resolutions () taken in the present analysis are 50% (0.085, 0.085), respectively. The normalization and energy calibration uncertainty for signal)/background are taken to be , , , and . For signal)/background) the values are , , , and .

4. Parameter Degeneracies

In general, different set of oscillation parameters may give identical predictions for oscillation probability resulting in parameter degeneracies and making it difficult to uniquely determine these parameters. The biprobability plots are well known constructions to study parameter degeneracies in presence of NSI parameters and phases. In presence of off-diagonal NSI the neutrino/antineutrino oscillation probability exhibits degeneracy between SI and NSI phase when and , where is Dirac phase in a model with NSI. If we presume that mixing angles and mass squared differences are known from some other experiment, then for every value of SI phase there will be three unknowns (, , and the phase ) which generate an off-diagonal NSI degeneracy. Accordingly, a measurement of and in an experiment provides two constraints; for each value of , a solution for and will exist for any value of resulting in parameter degeneracy. As a representative plot to depict standard parameter degeneracies we have shown, in Figure 1, mass hierarchy and octant degeneracy assuming NSI parameters and varying from to for DUNE experiment. All other NSI parameters are assumed to be zero. The solid and dashed ellipses correspond to higher octant (HO) and lower octant (LO) of , respectively, for normal hierarchy (NH). Similarly, the dotted and dash-dotted ellipses correspond to higher octant (HO) and lower octant (LO) of , respectively, for inverted hierarchy (IH). The significant overlapping of the ellipses shows fourfold octant and mass hierarchy degeneracy. For example, points of intersection of solid (dotted) and dashed (dash-dotted) ellipses exhibit octant degeneracy as neutrino and antineutrino oscillation probabilities are the same in both cases for normal (inverted) hierarchy.

4.1. Mass Hierarchy Degeneracy

In the determination of unknown neutrino mixing parameters an ambiguity exists in correlated way between and sign of . The sign of can be determined by measuring interference between the vacuum and the matter effects. The simultaneous determination of and sign of can be done in long baseline experiments in both standard [8] and nonstandard cases [6]. Thus, we study the prospects to resolve this degeneracy with two future long baseline experiments DUNE and T2HKK involving NSIs. Throughout this work we have considered the channel to study various parameter degeneracies. We have obtained plots for DUNE and T2HKK for two cases assuming (i) all NSI parameters to be nonzero, i.e., and NSI phases ; (ii) only to be nonzero () (Figure 2).

We have, also, shown sensitivity plots for mass hierarchy for cases (i) and (ii) in Figures 3(a) and 3(b) and Figures 3(c) and 3(d), respectively. We plot significance as a function of to study hierarchy sensitivity of DUNE and T2HKK experiments for true NH(true IH) in left panel(right panel) of Figure 3. The true NH (true IH) sensitivity plot is obtained by considering NH (IH) in the true spectrum and IH (NH) in test spectrum. We have marginalized over and by considering them in test spectrum. The statistical definition of for understanding the aspects of the mass hierarchy sensitivity plot in case of true NH is as follows: where and denote the number of true events and test events for NH and IH in the bin, respectively. The index runs over 1 to ; is the number of bins for particular experiment. For DUNE, , of 250 MeV width in 0.5-10 GeV range and for T2HKK , of 40 MeV width in 0.4-1.2 GeV range. The index describes type of mode, i.e., neutrino or antineutrino. for neutrino (antineutrino) mode. The index describes type of channel considered, i.e., appearance or disappearance. for appearance (disappearance) channel. The NSI parameter . We have minimized over the marginalization ranges of NSI parameters as given in Table 1.

4.2. Octant Degeneracy

The oscillation probabilities for disappearance channels ) (neutrino) and (antineutrino) show main contribution from . In case, is not maximal, we have two possibilities: either or . This ambiguity creates two solutions and . The resolution of uncertainty in due to octant degeneracy is important for precise measurement of . We study this degeneracy in both neutrino and antineutrino mode with DUNE and T2HKK assuming NSI parameters and (Figure 4).

To study the octant degeneracy, we have obtained octant sensitivity plot for true NH(true IH) in Figure 5. We plot significance () as a function of . One can define octant sensitivity by considering LO in true spectrum (HO in true spectrum) and HO in test spectrum(LO in test spectrum). The octant sensitivity tells us about the ability of an experiment to distinguish the lower octant from its higher octant. In order to obtain sensitivity plots for octant degeneracy, we have marginalized over sign (), , and NSI parameter and the in true NH case is

4.3. Phase Degeneracy

In presence of new physics, there may appear additional sources of violation other than due to Dirac-type phase . The effects will, in general, include contributions from both standard and nonstandard phases. As the value of phases and is not known and possible existence of violation in nature, we look for all possible values of which are distinct from the conserving values of and . For this purpose, we define fraction [32] as the fraction of total permitted range of , i.e., [] where violation effects corresponding to standard(nonstandard) phases can be explored. Also, we have excluded the conserving values by marginalizing over for which implies that the value of get fixed to its conserving values or (Figure 6(b)). So, we can write fraction at / which will be discussed in detail, in Section 5.

To resolve phase degeneracy, we need to find out the sensitivity with which an experiment can distinguish between conserving cases and violating cases. In standard oscillations there is only one degree of freedom in -function corresponding to standard phase . However, in case of neutrino oscillations with NSI there are two degrees of freedom due to standard and nonstandard phases . We define the -function as where are the number of true events for in the range and are the number of test events in bin for with .

To study violation discovery, we obtain plot for as a function of for true NH and true IH for DUNE, T2HKK, and DUNE+T2HKK experiments (Figure 7).

5. Results and Discussion

The degeneracy in sign of and can be resolved with the experiments which involve matter effects such as DUNE and T2HKK. In left (right) panel of Figure 2, we have shown the mass hierarchy degeneracy assuming all () NSI parameters, along with corresponding phases, to be nonzero for DUNE and T2HKK experiments. The meshed region between solid and dashed lines corresponds to the normal hierarchy (NH) whereas the forward-diagonal region between dash-dotted and dotted lines corresponds to the inverted hierarchy (IH) of neutrino masses. The normal and inverted hierarchy regions overlap for DUNE and T2HKK experiments due to contributions from all the nonzero NSI parameters (left panel) making it difficult to resolve the mass hierarchy degeneracy in these experiments. However, there is no overlap in the NH and IH regions for DUNE and T2HKK for energy range 1 to 4 GeV assuming only and corresponding phase to be nonzero (right panel). Furthermore, the future long baseline experiments such as DUNE and T2HKK with higher statistics and better energy resolution, focusing on neutrino beam energy between 1 and 4 GeV energy bracket, may provide better opportunity to resolve mass hierarchy degeneracy (Figures 2(b) and 2(d)). However, for the same energy range DUNE has brighter prospects to resolve the mass degeneracy than T2HKK.

In Figure 3, we have obtained the sensitivity plots for mass hierarchy in presence of all() NSI parameters for DUNE and T2HKK in first (second) row with true NH and true IH. Here, all NSI parameters referred to with corresponding phases and their true and test values have been taken from Table 1. For ready reference, we have, also, shown the sensitivity plots for the SI case in Figures 3(e) and 3(f). We find from Figures 3(a) and 3(b) that there is no distinction between true NH and true IH sensitivities for any value of at nearly () C.L. for DUNE (T2HKK) experiments. The hierarchy sensitivity get enhanced when only one NSI parameter and its corresponding phase is present as shown in Figures 3(c) and 3(d). The DUNE experiment in case of true NH shows stronger hierarchy sensitivity for as compared to true IH and is maximum at ( C.L.), whereas for the hierarchy sensitivity in true IH case is stronger and is maximum at ( C.L.). For T2HKK experiment the hierarchy sensitivity of true NH case is stronger than true IH case for the region and lies nearly at C.L. It can be noted from Figures 3(a) and 3(b) that decreases as changes from lower to higher octant which is in contradistinction to the case with one off-diagonal NSI (Figures 3(c) and 3(d)) and SI case (Figures 3(e) and 3(f)). Consequent to the presence of multiple NSI parameters and long baselines of DUNE and T2HKK, measurement will be severely affected by the NSI-modified matter effects due to degeneracies between SI and NSI parameters and between NSI parameters. Also, the mass hierarchy sensitivity will be decreased as compared to SI case due to cancellation effects induced by the off-diagonal NSI parameters. In particular, for all NSI parameters, at lower value of the cancellation effects will be relatively small as compared to SI matter effects which are reflected as increase in the value. However, cancellation effects will be appreciable for larger value of resulting in low .

In Figure 4, we have shown octant degeneracy in DUNE and T2HKK with NSI parameters and . It is evident from Figure 4 that octant degeneracy can be resolved with DUNE and T2HKK experiments using combination of neutrino and antineutrino oscillation modes. In Figure 4, the meshed region between solid and dashed lines represents the lower octant (LO) and the forward-diagonal region between dash-dotted and dotted lines represents the higher octant (HO). The left (right) panel in Figure 4 represents neutrino (antineutrino) mode of DUNE and T2HKK experiments. For LO(), with and for HO() with , both DUNE and T2HKK can resolve the octant degeneracies with neutrino mode only. Moreover, for LO with and HO with , octant degeneracies can be resolved with both DUNE and T2HKK with antineutrino mode only. Thus, in general, the neutrino and antineutrino modes are exigent to resolve octant degeneracy in DUNE and T2HKK experiments with matter NSI. Also, we find that the neutrino (antineutrino) beam energy bracket of 1 to 4 GeV can, simultaneously, resolve the mass hierarchy and octant degeneracies in DUNE and T2HKK.

In Figure 5, we have obtained the sensitivity plots of octant with true NH and true IH for DUNE and T2HKK. It can be seen from Figure 5 that the octant degeneracy can be resolved for both DUNE and T2HKK experiments for true NH and true IH. The DUNE experiment shows strong sensitivity for LO in true NH ( C.L.) and for HO in true IH ( C.L.). The T2HKK experiment has weak sensitivity to resolve octant degeneracy in both true NH and true IH case. We have, also, checked that there is not much improvement in octant sensitivity (over DUNE case) if we take DUNE and T2HKK conjunctively. For comparison, we have, also, shown the SI sensitivity curves of octant with true NH and true IH for DUNE and T2HKK.

The current and future neutrino oscillation experiments are diligently aiming at measuring neutrino mass hierarchy and violating phase . In presence of NSI (for example, assuming and nonzero) the situation becomes more complicated due to presence of additional sources of violation. The nature may intromit violation for wide range of ().

In Figure 6(a), we have shown asymmetry with neutrino beam energy for both and for DUNE+T2HKK. In Figure 6(a), the solid (dashed) line represents case whereas in Figure 6(b), the solid (dashed) line represents the fraction corresponding to the . In Figure 6(b), instead of focusing on measurement of phases () we have obtained all possible and values which are different from conserving values of and at 3 C.L. for DUNE+T2HKK. We have shown the effect of real NSI parameter () on the discovery reach of violation due to in Figure 6(b).

For real NSI case (), we calculate fraction , i.e., fraction of which is distinguishable from its conserving values, at 3 C.L. The parameter space for which significance is less than 3 C.L. has been excluded; thus, has been fixed to its conserving values . The has been minimized for value from to and the parameter space for which significance is considered. It is evident from Figure 6(b) that real NSI parameter () has no effect on the discovery reach of violation due to (), within the available bound on , at DUNE+T2HKK.

For complex NSI (), with , i.e., violation is only due to NSI phase , we have calculated fraction on similar lines as that for in real NSI case (replace by ). It can be seen from Figure 6(b) that for . However, for longer baseline may be larger even for small value(s) of . In this case, if violation is not observed at shorter baselines then the larger value of implies violation due to matter NSI. Also, it is observed from Figure 6(b) that fraction implies that there exist certain range(s) of for which violation is undetectable because significance is less than 3. For the case when violation is due to both and and lies below 3 significance, violation may be observed for certain range(s) of NSI parameters and . However, it will be difficult to disentangle the source of violation i.e., whether it is due to SI phase or NSI phase .

Figure 7 depicts the violation sensitivity plots for DUNE, T2HKK, and DUNE+T2HKK experiments with true NH and true IH. The violation sensitivity of T2HKK is stronger than DUNE in true NH irrespective of the value of (Figure 7(a)). For true IH, violation sensitivity of T2HKK is stronger than DUNE except for range of in region and region for which violation sensitivity of DUNE is stronger than T2HKK. The violation sensitivity of DUNE experiment is stronger in true NH than in true IH and is maximum at (), in true NH (IH). In combined analysis DUNE+T2HKK, the violation sensitivity increases to and for true NH and true IH, respectively. In Figure 7(b), the dip in the T2HKK sensitivity is due to the hierarchy- degeneracy. The values at the dip correspond to wrong hierarchy solutions. SI violation sensitivity curves for DUNE, T2HKK, and DUNE + T2HKK experiments with true NH and true IH have, also, been shown.

6. Conclusions

In conclusion, we have investigated the sensitivities of DUNE and T2HKK experiments to resolve mass hierarchy and octant degeneracies in presence of matter NSI. We have, also, studied the phase degeneracy due to standard and nonstandard phases for DUNE+T2HKK. The results are in consonance with the earlier studies on DUNE. We have analyzed the standard parameter degeneracies in presence of matter NSI for T2HKK experiment. We find that the mass hierarchy degeneracy cannot be resolved in presence of all NSI parameters due to their large experimental uncertainties (Figures 2(a) and 2(c)). However, it can be resolved for neutrino beam energy range 1 to 4 GeV in case of one nonzero NSI parameter and corresponding NSI phase for DUNE and T2HKK experiments (Figures 2(b) and 2(d)).

DUNE and T2HKK show poor sensitivity to resolve mass hierarchy degeneracy in presence of all NSI parameters (Figures 3(a) and 3(b)). However, the sensitivity to mass hierarchy gets enhanced when only one NSI parameter and its corresponding phase are present (Figures 3(c) and 3(d)). DUNE shows stronger hierarchy sensitivity for in true NH than true IH case with maximum sensitivity at ( C.L.), whereas for the hierarchy sensitivity in true IH case is stronger and is maximum at ( C.L.). T2HKK shows stronger hierarchy sensitivity in true NH case than true IH case for the region and lies nearly at C.L.

Furthermore, for LO (), with and for HO () with , both DUNE and T2HKK can resolve the octant degeneracies with neutrino mode only. Moreover, for LO with and HO with , octant degeneracies can be resolved with both DUNE and T2HKK with antineutrino mode only. Thus, combination of neutrino and antineutrino mode of DUNE and T2HKK can resolve the octant degeneracy.

The octant degeneracy can be resolved for both DUNE and T2HKK experiments for true NH and true IH (Figure 5). The DUNE experiment shows stronger sensitivity for LO in true NH ( C.L.) and for HO in true IH ( C.L.). The T2HKK experiment has weak sensitivity to resolve octant degeneracy in both true NH and true IH case. We have, also, checked that there is not much improvement in octant sensitivity (over DUNE case) if we take DUNE and T2HKK conjunctively.

The asymmetry from nonstandard phase is more than the standard phase for neutrino beam energy 1.5 to 7 GeV (Figure 6(a)). From Figure 6(b) we observe that real NSI parameter () has no effect on the discovery reach of violation due to (), within the available bound on , at DUNE+T2HKK. Also, for complex NSI () with , for . It increases with increase in and is 0.27 when . Also, it is observed that implying that there exist certain range(s) of for which violation is undetectable because significance is less than 3. For and lies below 3 significance, then violation may be observed for certain range(s) of NSI parameters and . However, it will be difficult to disentangle the source of violation, i.e., whether it is due to SI phase or NSI phase .

The violation sensitivity of T2HKK is stronger than DUNE in true NH (Figure 7(a)). However, there exists a range of in case of true IH for which violation sensitivity of DUNE is stronger than T2HKK. In combined analysis DUNE+T2HKK, the violation sensitivity increases to and for true NH and true IH, respectively.

Data Availability

The data used to support the findings of the present work are taken from [29, 30] and are openly accessible.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank Jogesh Rout for valuable discussions during this work. Surender Verma acknowledges the financial support provided by University Grants Commission (UGC), Basic Science Research (BSR), Government of India vide Grant no. F.20-2(03)/2013(BSR). Shankita Bhardwaj acknowledges the financial support provided by the Central University of Himachal Pradesh. The authors, also, acknowledge Department of Physics and Astronomical Science for providing necessary facility to carry out this work.