Study of Non-Standard Charged-Current Interactions at the MOMENT experiment

MuOn-decay MEdium baseline NeuTrino beam experiment (MOMENT) is a next-generation accelerator neutrino experiment looking for more physics study. We try to simulate neutrino oscillations confronting with Charged-Current\&Non-Standard neutrino Interactions(CC-NSIs) at MOMENT. These NSIs could alter neutrino production and detection processes and get involved in neutrino oscillation channels. We separate a perturbative discussion of oscillation channels at near and far detectors, and analyze parameter correlations with the impact of CC-NSIs. Taking $\delta_{cp}$ and $\theta_{23}$ as an example, we find that CC-NSIs can induce bias in precision measurements of standard oscillation parameters. In addition, a combination of near and far detectors using Gd-doped water cherenkov technology at MOMENT is able to provide good constraints of CC-NSIs happening at the neutrino production and detection processes.

Neutrinos are massless in Standard Model (SM) and phenomenon of neutrino oscillations is new physics beyond SM. It is required to generate massive neutrinos by extending SM, such popular candidates as seesaw models, supersymmetry models, extra-dimension models and the like. With more particle contents in new physics models, it might contain the sub-leading effects induced by non-standard neutrino interactions (NSIs) in nature. Effective operators towards this direction have been adopted to link neutrino mass models and NSIs [25][26][27][28]. Though there are viable models for sizable NSIs associated with neutral-current interactions [29][30][31][32][33], it is still a tough task to model CC-NSIs surviving the overwhelming constraints from precision measurements of charged lepton properties. Several studies have been conducted on NSIs from the experimental and model-building point of view [34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51]. A review of NSIs is given in detail in Ref. [52,53]. With the help of an effective field theory, we can generally integrate out the mediator/propagator in the Feynman diagram and keep four fermions contact with each other. New physics scale is then embedded into the effective coupling constant αβ α β where α/β or α β are the related fermion flavours. In theory the higher the new physics scale, the harder it is to reach a small effective coupling constant. We have reached an era of precision measurements of neutrino mixing parameters after an establishment of neutrino oscillation. It is promising for us to develop better neutrino detectors to search for sub-leading NSIs in the current and next-generation neutrino oscillation experiments as a complementary to the new physics search with the high intensity machine at the collider.
The MuOn-decay MEdium baseline NeuTrino beam experiment (MOMENT) is a next-generation accelerator neutrino experiment proposed for discovery of leptonic CP violation [54]. The atmospheric neutrino flux is a serious hindrance to the study of CP violation at MOMENT. Backgrounds caused by the atmospheric neutrinos exceeds the oscillation signal events significantly at O(100) MeV. Neutrino beams from such a continuous proton accelerator provide high luminosity fluxes but result in a loss of timing information which is traditionally used to suppress atmospheric neutrino backgrounds in the accelerator neutrino oscillation experiment with the pulsed proton beam facility. A new detector technology, however, might overcome the barrier and offer precision tests of θ 23 . The new detection technology might also lead to a discovery of the CP violating phase in the framework of 3-flavour neutrino oscillations [55], which complement the study at T2K and NOνA to solve the degeneracy problem and exclude CP conserved phase at a relatively high confidence level. In addition, a comprehensive study of the bounds on NSI parameters has been carried out. The bounds on NSI parameters governing the neutrino productions and detections are about one order of magnitude stronger than those related to neutrino propagation in matter, taking the current bounds on ud and µe as an example [56]: A special case study was performed at the Daya Bay reactor neutrino experiment to constrain NSI parameters where neutrinos are produced by beta decays and detected by inverse beta decays [57]. The relevant ee sector of µe got an upper limit of O(10 −3 ). With the help of a perturbation theory, neutrino oscillation probabilities in the presence of source/detector and matter NSIs can be found in the reference [58], which is motivating further study and optimization of new experimental proposals to pin down the current bounds. The first glimpse of NSI effects during neutrino propagation in matter at MOMENT has been shown in the reference [59]. It has discussed the sensitivity of neutralcurrent NSIs by means of accelerator neutrino oscillations in matter. However, the impact of source and detector NSIs associated with charged-current interactions has not been discussed. Within a theoretical model predicting new neutrino interactions, it is natural and fair for neutrinos to carry new charged-current and neutral-current interactions during the production process, the propagation process and the detection process. Furthermore, the current neutrino experiment T2K and NOνA are based on the superbeam neutrino production where neutrinos come from pion decays. We have to stress that NSIs associated with muon decays are very different from those happening at pion decays if we take the new physics into account. Therefore, it is necessary to bring source/detector NSIs for a complete analysis at MOMENT where neutrinos are produced by muon decays.
In this work, we explore the charged current NSIs effects at MOMENT. We focus on the precision measurement of standard neutrino mixing parameters and constraints of NSI parameters in the presence of non-standard chargedcurrent interactions at the source and detector. The paper is organized as follows: we discuss neutrino oscillation channels at a short and long distance in Section II. In Section III, we describe our implementations of MOMENT and details in simulation. In Section IV, we show the impacts of NSIs on precision measurements of standard neutrino parameters and present the correlations and constraints of NSI parameters within production and detection at MOMENT, and compare the expected results with current bounds. The summary follows in Section V.

II. DISCUSSION OF NEUTRINO OSCILLATION CHANNELS
The formalism of NSI is a general way of studying the impacts of new physics in neutrino oscillations. Without dealing with the matter NSIs due to the short baseline, we start with the neutrino production and detection processes involving non-standard interactions. These processes are often related to the charged lepton and called the chargedcurrent-like NSIs. The neutrinos at the MOMENT experiment are produced by the muon decay processes µ − → e − +ν e + ν µ and µ + → e + + ν e +ν µ and are detected mainly through quasielastic charged-current interactions: ν + n → p + − andν + p → n + + (Here denotes e or µ) in the neutrino detector. The CC-NSIs imposed on the production and detection are two different types: the NSIs involved in the muon decay production process are related to charged leptons, while the NSIs involved in the detection process are associated with quarks. For simplicity, we have restricted the operators to (V − A)(V − A) Lorentz structure and neglected NSIs including right handed neutrinos, where the process is helicity suppressed. One may ask why we ignore other potential Lorentz structures, since the most general way of constructing the four-fermion interactions could come with the current of (V ± A)(V ± A), (S ± P )(S ± P ) and T T , where V stands for vector couplings, A for axial-vector couplings, S for scalar couplings, P for pseudo-scalar couplings and T for the tensor couplings. There have been several attempts in the literature towards the chirality discussion of NSIs (see e.g. [60][61][62]). Except (V − A)(V − A), other structures are either helicity suppressed or very small due to their contributions by higher order corrections. Therefore, the most interesting CC-NSIs could be parameterized as the effective four-fermion Lagrangians at the detector: where the superscript "d" represents the interactions at the detector, G F is the Fermi constant, f (f ) represents d (u) quark, α (α = e, µ, τ ) is neutrino index and β (β = e, µ) is lepton index. The NSIs at the detector are parametrized by d αβ which give the strength of NSIs relative to G F . The detector projects the neutrino wave function not only onto the standard weak eigenstates, but onto a combination of them: In a muon-decay accelerator neutrino experiment, the neutrinos are produced by the muon decay process µ + → e + + ν e +ν µ and the charge conjugated process in the Standard Model. The effective lagrangian involving NSIs in production processes can be expressed as: where the superscript "s" represents the interactions at the source and two fermion fields "f " and "f " stand for a neutrino or a charged lepton in order to avoid the confusion of the neutrino flavour index and the index related to the γ matrix. For simplicity, we assume that the dominant NSI processes can interfere coherently with standard oscillations. Here we only consider two separate cases (take the µ + decay as an example): (1) µ + → e + + ν α +ν µ for NSIs s eα with any flavour (α = e, µ, τ ) by fixingν µ , where we assume f = µ, f = ν µ , γ = e.
Similarly, the neutrino flavour states produced at the source can be written as superpositions of pure flavour eigenstates: Thus, the oscillation probability is given by: Here the Hamiltonian takes the form of H = U diag(E + In principle, accelerator neutrinos pass through the earth matter until reaching a detector and matter effects change the oscillation probability. However, MOMENT is a medium baseline experiment and the matter effects are relatively small. Meanwhile, the main topic in the current study is delivered to NSIs happening at the production and detection processes rather than NSIs in matter. Therefore, we will only, for sake of simplicity, display the probabilities perturbatively in vacuum for related appearance and disappearance channels and try to extract useful information for physics performance study. Of course, matter effects are taken into account simultaneously in the complete simulation for physics performance. Since near and far detectors will be used in the simulation later, we will discuss the oscillation channels at a short and far distance separately.

A. Oscillation channels at a near detector
Here a near detector means detecting neutrinos at a distance of O(100) meters. In the standard oscillation frame without non-standard interactions, ν µ (ν µ ) and ν e (ν e ) can not develop neutrino oscillation patterns in such a short distance and their probabilities are equal to 1. However, NSIs are able to generate zero-distance effects so that the disappearance probabilities are allowed to be larger than 1, equal to 1, or smaller than 1. After dropping the terms O( 2 ), we approximate the probabilities as Eq. (8) and Eq. (9): It is easy to see that P (ν s e →ν d e )/P (ν s µ →ν d µ ) deviates from unity with some constant terms in the presence of relevant NSI parameters s ee and d ee ( s µµ and d µµ ). If neutrinos are produced with charged lepton decays and detected by identifying the same charged leptons, the contribution of s ee to the probability is equivalent to d ee and then the sensitivity to these two parameters should be the same at the near detector.
Similarly, the appearance channels ν µ → ν e (ν e → ν µ ) must remain zero in the standard 3-flavor neutrino scenario. After dropping the O( 3 ) and O( 4 ) terms, the expressions for ν e and ν µ appearance probabilities including CC-NSIs can be written as Eq. (10) and Eq. (11), respectively: Here the ν e (ν µ ) appearance probability would depend on s µe and d µe ( s eµ and d eµ ) after we introduce the NSIs at the neutrino source and detector. The probability of each oscillation channel and its conjugate partner shares the same form at the near detector. In fact, there are only four effective channels even though eight neutrino oscillation channels can get involved in the MOMENT experiment. It is a discovery of new physics to observe zero-distance effects at near detectors for disappearance or appearance channels.

B. Oscillation channels at a far detector
Oscillation patterns get more complicated as soon as we consider channels suitable for the far detector at MOMENT. In the standard framework describing three neutrino mixings, the probability of ν µ → ν e channel is calculated by a simple change of sign of the sin δ term in a T-reversed channel of ν e → ν µ . Due to CC-NSIs, ν µ → ν e and ν e → ν µ probabilities are not so obvious any more. We perturbatively derive the explicit expressions of their probabilities in vaccum as given in Eq. (12) and Eq. (16), considering α = νe→νµ ). For an oscillation channel of ν µ → ν e , the probability can be written as: with In a similar way, the probability of ν e → ν µ can be expressed as: with Here s ij = sin θ ij , c ij = cos θ ij , s 2×ij = sin 2θ ij , c 2×ij = cos 2θ ij , ∆ 31 = Our derived oscillation probabilities are consistent with the reference [58]. We can immediately read the impacts of NSI parameters on appearance channels in order: • It is clear that is approximately equivalent to enlarging or compressing the amplitude of sin 2 ∆ 31 . Meanwhile, the term of d µe ( s eµ ) will change the maximal position of its oscillation probability.
• Since these NSI parameters are entangled with standard mixing parameters, NSIs would interfere with precision measurements of the standard CP-violating phase δ and θ 23 which is manifested by related terms in Eq. (12) and Eq. (16). These two channels are sensitive to the octant of θ 23 beacuse the leading order of Eq. (12) and Eq. (16) depend on sin 2 θ 23 . The presence of NSIs may induce a wrong determination of θ 23 . Fig. 1 shows a specific situation of the degeneracy caused by NSIs at the probability level. The continuous line and the dashed line show the case of θ 23 = 41.6 • and θ 23 = 48.4 • without NSI, respectively. We can see a clear separation between two scenarios. However, introducing NSIs can shift the amplitude of the probability and fake the contribution of θ 23 . For instance, the dotted line shows the special case with | s µe | = 0.02 and φ s µe = π/2. It is clear that the dotted oscillation curve almost coincides with the continuous curve. This indicates that we may get the wrong measurement of θ 23 , if there is indeed CC-NSIs but we ignore them by only fitting the data within the standard neutrino oscillation framework.
• An assumption of real NSI parameters with all NSI phases φ = 0 could leads to a substantial simplification in each channel, since those terms proportional to sin φ s/d αβ would vanish. One can observe that the sensitivity to d τ e ( s eτ ) from the oscillation channel ν µ → ν e (ν e → ν µ ) would be tiny for δ = 3π/2, π/2. Moreover, we can obtain a linear correlation of s µe and d µe (or s eµ and d eµ ). In the case of δ = 0/π, sin δ term will vanish and we can get the linear correlation of s µe , d µe and d τ e ( s eµ , d eµ and s eτ ). Futhermore, if θ 23 is close to 45 degrees, the contributions from d µe ( s eµ ) tend to disappear. In brief, the sensitivities to the related CC-NSI parameters at a far detector vary with the standard parameters θ 23 and δ.
Similar to the discussion above, ν µ and ν e disappearance probabilities are given as Eq. (20) and Eq. (21): IG. 2. The oscillation probabilities of P (νµ → νe) and P (νe → νµ) with/without CC-NSIs at a far detector given the baseline of 150 km. In both panels, the shaded region in cyan is due to the variation of relevant NSI parameters as is shown in the legend. The dashed red line represents the case without NSI parameters. Here we assume the normal mass hierarchy and take the following inputs: ∆m 2 21 = 7.5 × 10 −5 eV 2 , θ12 = 33.56 • , ∆m 2 31 = 2.524 × 10 −3 eV 2 and θ13 = 8.46 • .
Apart from standard neutrino oscillations, major contributions come from the terms proportional to . The ν µ → ν µ (ν µ →ν µ ) is an important channel to measure θ 23 , which is expected to judge whether θ 23 is maximal or not. The channel of ν e → ν e (ν e →ν e ) is good at precision measurements of θ 13 . Without NSIs, ν e (ν e ) disappearance channel has no dependence of the standard CP phase. After introducing NSI parameters, however, even the standard CPviolating phase would appear in the ν e (ν e ) disappearance probability. Under the assumption of φ = 0, we will have a rather simplified correlation of standard neutrino oscillation and NSI parameters. A few comments for this special case are given as follows: • The oscillation channel ν µ → ν µ (ν µ →ν µ ) is affected by | On the contrary, these four parameters can be well constrained in appearance channels.
In Fig. 2, we make a comparison of the probabilities P (ν µ → ν e ) and P (ν e → ν µ ) with/without NSIs. We turn on nonvanishing parameters of | s µe | = | d µe | = | d τ e | = 0.03 for an illustration. The amplitude of oscillation pattern is shifted significantly by these NSI impacts as we point out from the perturbative approximation. We display the probability curve of standard three-flavour neutrino oscillation and the probability bands originating from the variations of NSI parameters. The shaded region in Fig. 2 shows changes of oscillation probabilities when the standard CP-violating phase is 3π/2 and the NSI phases vary in [−π, π]. It is easy to see the difference between P (ν µ → ν e ) and P (ν e → ν µ ) without NSI effects, where P (ν µ → ν e ) − P (ν e → ν µ ) becomes negative for δ ∈ (0, π) and P (ν µ → ν e ) is larger than P (ν e → ν µ ) for δ ∈ (π, 2π); with δ = 0/π, there is no difference between P (ν µ → ν e ) and P (ν e → ν µ ). Once we turn on NSI parameters, however, the relationship between the probabilities of these two channels will be changed significantly as can be seen from the shaded regions. After we introduce the characteristics of MOMENT and simulation details in the next section, we will be further convinced by an analysis of event rates as shown in Fig. 3.

III. CHARACTERISTICS OF MOMENT AND THE DETAILS OF SIMULATION
Continuous high-energy proton beam at 1.5 GeV with a power of 15 MW will be produced in the accelerator facility. A choice of target stations with fluidised tungsten is still under optimization to generate charged mesons most of which are pions and kaons. We can expect 1.1 × 10 24 Proton On Target (POT) per year. A magnetic solenoid will be deployed to make the pion beam focused and selected. The curvature of the solenoid helps with selecting muons from pion decays, followed by a straight tunnel to prepare neutrinos from muon decays. The neutrino fluxes are provided by the accelerator working group in MOMENT for our physics performance study [68]. We intend to extract more information from eight oscillation channels using the muon-decay neutrino beams in the simulation study: ν e → ν e , ν e → ν µ , ν µ → ν e , ν µ → ν µ and their conjugate partners. Since we have to conduct flavour and charge identifications to distinguish secondary particles, we consider the new technology using Gd-doping water to separate both Cherenkov and coincident signals from capture of thermal neutrons. Muon taggings can be efficiently obtained by daughter electrons together with pulse shape discrimination of waveforms. We follow the detector description from a sophisticated study in the CERN-MEMPHYS project [69] and update the related new technology with regard to Gd-doping water. The major backgrounds for MOMENT come from atmospheric neutrinos. We believe that they can be suppressed by the beam direction and proper modelling background spectra within the beam-off period, which is to be extensively studied in detector simulations in future study.
Charged-current interactions are used to identify neutrino signals: ν e + n → p + e −ν µ + p → n + µ + ν e + p → n + e + ν µ + n → p + µ − A few remarks for signals and backgrounds are given as follows: • Gd doping into pure water could be used to discriminate electron neutrinos and antineutrinos by whether there is a capture of the scattered thermal neutron or not. Neutron capture on Gd emits the 8 MeV gamma rays. Theν e (ν µ ) signal is reconstructed by tagging the neutron in coincidence with the positron to suppress most of backgrounds associated with single events. While water Cherenkov detection is not significantly changed, ν e signals come from the Cherenkov ring created by ν e elastic scattering with electrons.
• In a Water Cherenkov detector, electron and muon-flavour neutrinos are well separated by event reconstructions, where the former type creates electron showers and the latter type leads to muon tracks. Sometimes, low-energy muons decay and cause flavour misidentifications. Here we ignore the flavour misidentification in the simulation.
• The imperfection of detectors leads to misidentifications of charge-current interactions for neutrino signals. Here we suppose their effects are negligible.
• It is also possible for neutral-current interactions with accidental single events to be identified as the coincident signal. A neutron knocks a nucleus off of an oxygen, resulting in excited states and photons from de-excitations will mimickν e signals oscillated from accelerator neutrino beams. We expect them to be rather small and assign an extremely small background over signal ratio in our simulation. Table II lists the simulation details about the neutrino detector. A baseline of 150 km is assigned in the current proposal based on the neutrino beam energy range [54]. We assume a near detector with a fiducial mass of 100 t and a far detector with a fiducial mass of 500 kton. The running time is 5 years for each polarity. In the massive water cherenkov detector, we follow electron and muon selection efficiencies given in Ref. [70,71]. With regard to the normalization error on signals and the normalization error on backgrounds, we assume they are at the level of 5%. As for the atmospheric backgrounds, they could be suppressed via sending the neutrino beam in short bunches with a suppression factor of 2.2 × 10 −3 [72]. The cross section for quasi-elastic interactions is taken from the reference [73]. The values of the standard neutrino oscillation parameters are taken from the latest nu-fit results [12]. Table I shows the central values and their uncertainties in the present work. Unless otherwise mentioned, we expect a determination of mass hierarchy without NSIs before running MOMENT and assume the normal mass hierarchy in our simulation without a loss of generality, i.e. ∆m 2 31 > 0. A list of assumptions for near and far detectors are given in Table II. We present numerical results by simulating the neutrino oscillation signals and backgrounds using GLoBES [74,75]. Similar to the probability-level analysis given in Fig. 2, we present event rates of ν µ → ν e and ν e → ν µ channels versus the neutrino energy in Fig. 3. The event spectra are shifted significantly after we consider CC-NSI effects. The shaded region highlights the large variation from the new CP phases caused by CC-NSIs even if we fix their strength of NSI couplings. It is then straightforward to discuss physics performance for the MOMENT experiment with simulated event spectra. We calculate the event rates by defining the true values (central values) for standard oscillation parameters and fit with/without NSI impacts to extract useful information. We compute the χ 2 using the following approach: where the index j denotes the channel number and i denotes the bin number, ρ true and ρ test represent the standard vectors of true and test values, respectively, s is the vector of systematics related to the neutrino beam and detector, true and test are the non-standard vectors of true values and test values, and N ij is the expected event for the j-th channel and i-th bin. The second term corresponds to the contribution to χ 2 from the external inputs which are based on results from previous experiments. σ ρα is the external error (input error) in GLoBES imposed on the central values. Similarly, the third term represents the treatment of systematic errors implemented on the χ 2 .

A. Impacts on precision measurements of standard mixing parameters by CC-NSIs
The CKM mixing matrix is well measured in the quark sector at the sub-percent level [76], while mixing parameters in the lepton sector are far away from such a precision. It is very likely for the next-generation experiment like MOMENT to achieve the goal of doing precision measurements. In Fig. 1, we have showed that CC-NSIs may induce a bias in precision measurements of θ 23 at the probability level. In this section, we take the Dirac CP-violating phase δ cp and θ 23 as an illustration to show the impacts from CC-NSIs after the simulation. The true value of θ 23 is taken as 41.6 • . We choose two δ values with δ = π/2 and δ = 3π/2 to simulate all oscillation channels at MOMENT and fit  the neutrino spectra with/without NSIs. Fig. 4 demonstrates the numerical results. The true values of the standard oscillation parameters are shown by a red point in each panel. In all sub-figures (a), (b) and (c), we have considered uncertainties of standard mixing parameters. Panels (a1) and (a2) show the determination of δ cp and θ 23 in the case of the standard neutrino oscillation without NSIs. By running MOMENT, we can determine the mixing angle θ 23 with an error bar of one degree at the 3σ confidence level, while the precision for δ cp is good enough. In sub-figure (b) and (c), NSIs happening at the source and detector are turned on. All the corresponding CC-NSI phases can vary within (0, 2π). In the panel (b1) and (b2), we only consider the CC-NSIs ( s eµ , s eτ and d eµ , and their marginalization ranges are allowed within the current bounds given in Eq. (1) and Eq. (2).) which are related to ν µ (ν µ ) appearance channels. The panels of (b1) and (b2) show the enlarged uncertainties in parameter fittings. Especially, a degeneracy pops up in the measurement of θ 23 -δ cp for the case of δ = 3π/2, while it is still safe for the case of δ = π/2. Furthermore, we go to panels (c1)/(c2) by turning on those CC-NSIs related to the ν e (ν e ) and ν µ (ν µ ) appearance channels ( s eµ , s eτ , d eµ , s µe , d µe , d τ e and their marginalization ranges are within current bounds given in Eq. (1) and Eq. (2)). As we have discussed in Section II, the ν e appearance channel is affected by the parameters: s µe , d µe , d τ e , while the ν µ appearance channel is mainly determined by the parameters s eµ , s eτ , and d eµ . This feature can be understood after a closer look at the Eq. (16). If δ is equal to π/2, the peak of the probability in ν e → ν µ channel is much larger than the case of δ = 3π/2. In turn, the event rate in the detector for δ = π/2 is much higher. The corresponding fitted results are much better. Therefore, CC-NSI parameters destroy the precise determinations of standard mixing parameters. As can be seen from panels (c1) and (c2), the degeneracy even shows up at high confidence levels when we consider all the relevant NSI parameters. We might get into the wrong best-fit region if we neglect the CC-NSIs from the new physics. A combination of different neutrino oscillation experiments might resolve such an ambiguity and finish the task of precision measurements of neutrino mixing parameters to the same level as quark mixing parameters. Or we need a more powerful machine, such as a neutrino factory. We have introduced the NSIs by integrating the potential heavy propagator from the new physics scale based on the effective theory. Each NSI parameter has a magnitude which tells us the strength of new couplings and its associated phase to bridge the CP violating story. Therefore, it is convenient to adopt two methods: one is to take the NSI parameters as real, which corresponds to the strength of the coupling constant or switch off the NSI-induced CP violation phases; the other is to keep general assumptions given complex NSI parameters. In this section we discuss the constraints on source and detector NSIs from the far and near detector, respectively. For the former case, Table III demonstrates the sensitivity of MOMENT in constraining the NSI parameters using the single-parameter-fit at 90% C.L. The results are obtained by using all the neutrino and antineutrino oscillation channels with both near and far detectors. In a comparision of the current bounds and the expected limits from our simulation, we find that the MOMENT experiment with running time of 5+5 years has a potential to improve the constraints for the CC NSIs. In Table III, one can observe that most of the bounds for NSIs are improved except for s eτ , s µτ , d τ e and d τ µ which are marked by red fonts. These results also confirm our observations in Section II: • NSI parameters s/d ee and s/d µµ , at the leading order of NSI terms in the disappearance channels can be well constrained with a combination of near and far detectors.
• s eτ , s µτ , d τ e and d τ µ cannot be well constrained by the near detector due to their negligible effects. Their constraints by the far detector are even weaker than the current bounds since their contributions to the dominant order of O( s 13 ) in Eq. (12) and Eq. (16) would be zero when δ = 270 degrees.
• The constraints on s eµ , s µe , d eµ , d µe are mainly from the appearance channels. A combination of near and far detectors have stronger constraints on s eµ and d eµ than s µe and d µe . Fig. 3 can be used to explain this property: the channel ν e → ν µ has more events compared to the channel ν µ → ν e . In addition, it should be noted that ν e andν e disappearance channels are also sensitive to s eµ and d µe .
In the following part, we will focus on the exclusion curve of the amplitude versus its corresponding phase for each NSI parameter. Based on the previous discussions in Section. II, we neglect certain parameters which have trivial influence on the probabilities. We generate the event spectra with the true central values of standard oscillation parameters given in Table I. Then we turn on one NSI parameter and scan its amplitude and phase to fit the data. At the near detector, we pay more attention to these parameters:  Fig. 5 and Fig. 6, respectively. In Fig. 5, we show the excluded parameter space without any color at a near detector. One can observe that ν µ disappearance channel almost has the same performance with the ν e disappearance channel.  ). This is because the term containing the phases φ s µe and φ d µe plays an important role in two amplitudes as can be seen from Eq. (10) and Eq. (11). Appearance channels can well constrain the magnitude of related CC-NSI parameters while they barely have impacts on their phases. On the other hand, disappearance channels exclude a large parameter space allowed in the current experimental bounds. There is a strong correlation between the coupling strength and its phases in neutrino oscillation experiments.
In Fig. 6, we switch to the sensitivities of CC-NSI parameters at a far detector. The colorful regions are allowed after we run a far detector at the MOMENT experiment. Compared to the current bounds marked by the dashed  red lines, we obtain good constraints on most of NSI parameters at MOMENT, especially for s µe and d eµ . Here we list the features about the exclusion curves at the far detector for MOMENT: • The oscillation channels ν µ → ν e (ν µ →ν e ) and ν e → ν µ (ν e →ν µ ) are T-conjugate inverse of each other, leading to the symmetry between their NSIs: s µe and d eµ have equal contributions in the probability of P F D νµ→νe and P F D νe→νµ , respectively. Similarly, the pair of d τ e and s eτ and the pair of d µe and s eµ follow the same way. Therefore, these pairs of parameters have similar behaviour in Fig. 6, which can be manifested by Eq. (14) and Eq. (18) in Section. II. The dependence of constraints for | s µe | (| d eµ |) on the corresponding phase is not strong, because | s µe | (| d eµ |) depends on both terms of sin(δ + φ s µe ) (or sin(δ − φ d eµ )) and cos(δ + φ s µe ) (or cos(δ − φ d eµ )), which complement each other when varying the phase φ s µe (φ d eµ ). Although | d eµ | (| s eµ |) also depends on both terms of cos(δ + φ d µe ) (or cos(δ − φ s eµ )) and sin(δ + φ d µe ) (or sin(δ − φ s eµ )), the former term is suppressed by the coefficient c 2×23 . Thus, the exclusion curve of | d µe | − φ d µe (| s eµ | − φ s eµ ) mainly varies with the sine term. However, the constraint on | d τ e | (| s eτ |) only depends on the term of cos(δ + φ d τ e ) (or cos(δ − φ s eτ )) so that the limit to | d τ e | (| s eτ |) will be extremely weak for φ d τ e (φ s eτ )=0, ±π and be the best for φ d τ e (φ s eτ )=±π/2.
• Since we can take advantage of the ν e (ν e ) channel in MOMENT, obtaining the sensitivities of s eτ , s eµ and d eµ would be accessible. It is noted that these parameters can not be constrained well in superbeam experiments.
• The sensitivities to s eµ and s eτ are mainly extracted from the ν µ andν µ appearance channels. It should be noted, however, that the ν e andν e disappearance channels can help to enhance the constraints of them. Similarly, with contributions from ν e andν e disappearance channels, the sensitivities to d µe and d τ e will be improved.
• The constraints to NSI parameters s ee , s µµ , s µτ , d τ µ are extracted from the disappearance channels. The exclusion curves for s/d ee and s/d µµ are similar to each other at the near detector, since they are entangled with the same term in oscillation probabilities. In addition, there are some symmetric relationships between the source and detector NSI parameters, such as the pair of s µτ and d τ µ , the pair of s ee and d ee , the pair of s µµ and d µµ . This is not surprising at all, since the pair of effective coupling constants will be the same for neutrinos produced and detected related to the same charged leptons.
The far detector has a good sensitivity to NSI parameters, especially for d eµ and s µe . The numerical results of the correlations between the amplitudes and phases can be interpreted with the previous probability-level discussions. Almost all NSI-induced phases change the exclusion limits severely except the e-mu sector. Meanwhile, limits on other sectors are not as good as those on the e-mu stamped CC-NSIs. Therefore, MOMENT using muon-decay beams has its unique capability of improving the constraints on d eµ and s µe .

V. SUMMARY AND CONCLUSIONS
New Physics beyond SM might cause non-standard neutrino interactions and leave imprints on the neutrino oscillation. The next-generation accelerator neutrino experiment MOMENT intends to produce the powerful neutrino beam with an energy of O(100) MeV by muon decays and leaves plenty of room for detector selections and physics study. At this energy range, quasi-elastic neutrino interactions dominate the detection process and backgrounds from π 0 are highly suppressed. Compared with traditional superbeams from charged meson decays where intrinsic backgrounds have to be alleviated by the off-axis technology like T2K and NOνA, beams from muon decays are cleaner neutrino sources and good at a detection of new physics. CC-NSIs happening at neutrino productions and detections point to the new phenomenon, where a neutrino produced or detected together with the charged lepton will not necessarily share the same flavour, and flavour conversion is present already at the interaction level and "oscillations" can occur at zero distance. With the capability of flavour and charge identifications, we have an opportunity to use eight appearance and disappearance oscillation channels in the physics study. We have chosen the advanced neutrino detector using the Gd-doped Water Cherenkov technology and studied neutrino oscillations confronting with CC-NSIs at the MOMENT experiment.
In order to understand the relevant behavior from NSIs, we have perturbatively derived oscillation probabilities including CC-NSIs at a short and far distance, and tried to analyze parameter correlations of standard neutrino mixing parameters with NSI parameters. We have investigated impacts of the charged current NSIs at the neutrino oscillation probabilities, selected the following dominanting CC-NSI parameters A near detector at MOMENT is good at detecting the zero-distance effects induced by NSIs while the oscillation pattern would have not been developed in the standard neutrino oscillation paradigm. With near and far detectors, we have found that CC-NSIs can induce bias in precision measurements of standard mixing parameters. Taking δ cp and θ 23 as an example, we have found degeneracies after introducing CC-NSI parameters. With a non-maximal θ 23 , its degeneracy with the standard CP phase δ CP gets much worse if CC-NSIs appear in the neutrino production and detection processes. The current bounds on NSI parameters governing the neutrino productions and detections are about one order of magnitude stronger than those related to neutrino propagation in matter. Our study has shown that a combination of near and far detectors at MOMENT is able to provide lower bounds on CC-NSIs where a factor of about two can be envisaged for most of parameters compared with the current experimental bounds, as shown in Table III. We have found strong correlations of NSIs and constrained NSI parameters using a combination of near and far detectors at MOMENT.
The feasibility of physics performance at MOMENT strongly depends on inputs of the accelerator facility and the advanced neutrino detection technology. In the future, results will be further improved by tuning the beam energy and optimizing the baseline. We hope that our study will boost the the research and development activities for MOMENT.