Sensitivity to sterile neutrino mixing using reactor antineutrinos

The reactor antineutrinos are used for the precise measurement of oscillation parameters in the 3-neutrino model, and also used to investigate active-sterile neutrino mixing sensitivity in the 3+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3+1$$\end{document} neutrino framework. In the present work, we study the feasibility of sterile neutrino search with the indian scintillator matrix for reactor anti-neutrino (ISMRAN) experimental set-up using electron antineutrinos (ν¯e\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{\nu }}_e$$\end{document}) produced from reactor as a source. The so-called 3+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3+1$$\end{document} scenario is considered for active-sterile neutrino mixing, which leads to projected exclusion curves in the sterile neutrino mass and mixing angle plane. The analysis is performed considering both the reactor and detector related parameters. It is found that, the ISMRAN set-up can observe the active-sterile neutrino mixing sensitivity for sin22θ14≥0.064\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sin ^{2}2\theta _{14} \ge 0.064$$\end{document} and Δm412=1.0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varDelta m^{2}_{41} = 1.0$$\end{document}eV2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {eV}^2$$\end{document} at 90%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$90\%$$\end{document} confidence level for an exposure of 1 ton-year by using neutrinos produced from the DHRUVA reactor with thermal power of 100MWth\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$100 \hbox {MW}_{{th}}$$\end{document}. It is also observed that, there is a significant improvement of the active-sterile neutrino mixing parameter sin22θ14\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sin ^{2}2\theta _{14}$$\end{document} to ∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document} 0.03 at the same Δm412\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varDelta m^{2}_{41}$$\end{document} by putting the ISMRAN detector set-up at a distance of 20 m from the compact proto-type fast breeder reactor facility with thermal power of 1250 MWth\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {MW}_{{th}}$$\end{document}.


Introduction
The phenomena of neutrino oscillation have been established by several experiments using neutrinos from both natural (atmospheric and solar) and man-made (reactor and accelerator) sources. It shows the mixing between flavor and mass eigenstates, hence established that neutrinos have non-zero masses. Presently, the study of neutrino physics is in the precision era. However, experimental observations from various short baseline (SBL) experiments cannot be explained a e-mail: shiba@barc.gov.in b e-mail: dkmishra@barc.gov.in by 3-neutrino mixing paradigm which requires new additional neutrino called as 'sterile neutrino'. So the concept of this sterile neutrino could explain the results from GALLEX [1,2] and SAGE [3,4] Gallium experiments, find a deficit in neutrino flux while calibrating the detectors with radioactive sources. They have reported that the ratio of numbers of observed to predicted events is 0.88±0.05 [5] and it is known as "Gallium anomaly". The accelerator based SBL experiments such as Liquid Scintillator Neutrino Detector (LSND) [6] at a baseline of ∼ 30 m observed an unexplained excess of electron anti-neutrinos (ν e ) in a muon anti-neutrino beam. The MiniBooNE experiment also observed similar excess in ν μ → ν e mode [7]. The recent MiniBooNE data are consistent with the excess of events reported by the LSND. The significance of the combined analysis of both the experiments is an excess of 6.0σ [8]. There is an anomalous behavior has also been observed in the measurement of the reactor ν e flux and spectrum. The precise energy spectrum of antineutrino flux produced by the reactors are recalculated by Mueller et al. [9] which shows a significantly about 6% higher than experimental measurements at small distance. This discrepancy between the predicted and observed reactor antineutrino flux is known as the "reactor antineutrino anomaly" (RAA) [10]. There are basically two possible explanations for this discrepancy. One is the incomplete reactor models or nuclear data due to underestimated systematics of the measurements of beta spectra emitted after fission [11][12][13] or of the conversion method [9,[14][15][16]. The other explanation is an oscillation of ν e into a fourth light sterile neutrino. Moreover, measurements of the reactor ν e spectra show a discrepancy compared to predictions, particularly at energies of ∼ 5 MeV. The discrepancy in ν e spectra is confirmed by RENO [17], Daya Bay [18], Double Chooz [19], and NEOS [20] collaborations by measuring the reactor ν e energy spectrum. The distortion in energy spectra has been correlated to the reactor power [18], which may be due to the 235 U fuel [21]. In order to verify the existence of active to sterile neutrino oscillation hypothesis as the possible origin of the RAA and, also to clarify the origin of the bump at 5 MeV in the ν e spectra, there are several experiments underway and some will take data soon.
To address the RAA, the SBL experiments are aiming to measure the reactor ν e spectrum at two or more different distances and trying to reconstruct the ν e survival probability both as a function of energy and source to detector distance, L independent of any reactor model prediction. The L dependence is what gives the cleanest signal in the case of the sterile neutrino, and studying the ratio of the spectra measured at two different distances allows to avoid the problem of the theoretical spectrum. The DANSS group has performed the experiment at 3 distances from reactor core varied from 10.7 m to 12 m to find out the active-sterile neutrino mixing by measuring the positron energy spectra. They have observed that the excluded area in the sin 2 2θ 14 −Δm 2 41 (= m 2 4 −m 2 1 ) plane covers a wide range of the sterile neutrino parameters up to sin 2 2θ 14 < 0.01 [22]. Similarly, the STEREO [23] group has measured the antineutrino energy spectrum in six different detector cells covering baselines between 9 and 11 m from the compact core of the ILL research reactor. Their results are compatible with the null oscillation hypothesis and the best fit of the reactor antineutrino anomaly is excluded at 97.5% confidence level. Recently, PROSPECT group has measured the reactor ν e spectra using a movable segmented detector array and their observation disfavors the RAA best fit point at 2.2σ C.L. and constrains significant portions of the previously allowed parameter space at 95% confidence level [24]. This paper presents the results of an investigation on finding a possible mixing of a single sterile neutrino with the three known active neutrinos, vi z. the (3 + 1) model. It is the only allowed active-sterile neutrino mixing scheme [25] under the assumption of four neutrino model. At SBL, the presence of sterile neutrinos with squared mass difference Δm 2 41 ∼ 1eV 2 leads to fast oscillations resulting the reduction of reactor ν e flux, otherwise absent in the standard 3neutrino paradigm. This study quantifies the sensitivity of indian scintillator matrix for reactor anti-neutrino (ISMRAN) experimental set-up in constraining the active-sterile neutrino mixing parameters. In this work, we have considered various reactor (vi z. thermal power, core size, duty cycle, burn up) as well as detector response related parameters (vi z. energy resolution and detection efficiency) and also at several reactor core to detector distance while constraining active-sterile neutrino mixing at an exposure of 1 ton-year.
The outline of the paper is as follows. A detailed description of the ISMRAN detector set-up and the neutrino detection principle is discussed in Sect. 2 and in Sect. 3, respectively. The sterile neutrino oscillation formalism is introduced in Sect. 4. The incorporation of detector resolutions on neutrino induced true events is discussed in Sect. 5. The statistical analysis using both oscillated and without oscillated events based on χ 2 estimation is given in Sect. 6. The sensi-tivity to sterile neutrino mixing at an exposure of 1 ton-year is discussed in Sect. 7. Finally, in Sect. 8, we summarize our findings and discuss the implication of this work.

ISMRAN detector
The ISMRAN experimental set-up is being developed for detecting reactor ν e , searching for possible existence of sterile neutrino and monitoring of reactor power at the DHRUVA reactor facility in Bhabha Atomic Research Centre (BARC), India. The ISMRAN detector set-up will consist of an array of 100 plastic scintillator (PS) bars with weight of about 1 ton [26]. The dimension of each PS bar is 100 cm×10 cm×10 cm wrapped with Gadolinium coated aluminized mylar foils. Each PS bar is coupled with two 3" Photo-multiplier tubes at both ends. The schematic of the detector set-up is shown in Fig. 1. The advantage of the ISMRAN set-up is that, it is compact in size and maneuvered from one place to another easily. Also the segmented detector array can provide the additional position information while recconstructing the neutrino induced events and can improve the active sterile neutrino mixing sensitivity of the ISMRAN detector. To suppress both the natural and reactor related background, detectors are covered by a passive shielding material Lead (for gamma rays) of 10 cm thick and then followed by 10 cm thick of borated polyethylene (for neutrons). The detector is positioned at a distance of ∼ 13 m from the center of a cylindrical reactor core and can be moved closer to the core upto 7 m. The reactor has radius ∼ 1.5m and height ∼ 3.03 m Fig. 1 Schematic of 100 cm × 10 cm × 10 cm plastic scintillators array with shielding materials for reactor antineutrino measurement (defined as an extended source) [27]. The reactor can operate at a maximum thermal power of 100 MW th consuming natural uranium as fuel and producing about 10 19 ν e /s. In future, it is planned to put the detector set-up at proto-type fast breeder reactor (PFBR) facility, IGCAR, Kalpakkam, India [28]. The PFBR has dimension of about 1 m both in radius and height (defined as a compact source), and can operate at a maximum thermal power of 1250 MW th . As the reactor is compact and produces higher thermal power, it is an ideal case to utilize the detector set-up for investigating the active-sterile neutrino mixing. With this experimental set-up, it can be possible to confirm or reject the existence of a light sterile neutrino by measuring the ν e flux and energy spectra. At present a prototype ISMRAN set-up of 1/5-th of the final detector volume which is under operation at DHRUVA reactor facility [26].

ν e Detection principle
The electron-antineutrinos produce from the reactor interact with protons in the PS bars, via the inverse beta decay (IBD) process, The Q value of the above reaction is about −1.80 MeV and hence it limits the detection of antineutrinos. The positron which carries almost all of the available energy, loses it by ionization process in the detector and gets annihilated producing two gammas. The energy loss of the positron constitutes the prompt signal along with the Compton scattered annihilated gammas given by where Eν e is the energy of electron-antineutrino. So from Eq. 2, it is observed that there is a one-to-one correspondence between the positron energy and Eν e . The neutron produced in Eq. 1 carries a few keV's of energy and gets thermalized in collisions with protons in the PS bar. The neutron takes about 180μs in order to gets captured by proton in the PS bar produces gamma ray which is considered as a delayed signal. To further decrease the neutron captured time and improve the detector efficiency, PS bars are wrapped with Gadolinium coated aluminized mylar foil which has very high neutron captured cross-section. Further, the neutron captured time reduces to about 30-40 μs and a cascade of gamma rays produce with total energy ∼ 8 MeV due to Gadolinium. The coincidence of a prompt positron signal and a delayed signal from neutron captured by Gadolinium (Gd) uniquely identifies the IBD event.

Neutrino oscillation probability with 3 + 1 model
The sterile neutrino oscillation probabilities are based on expansion of the three generation Pontecorvo-Maki-Nakagawa-Sakata (PMNS) [29] matrix to 3 + 1 generation, where "3" stands for three active neutrinos (ν e , ν μ , ν τ ) and "1" for a sterile neutrino (ν s ). The neutrino flavors and mass eigenstates are related through ⎛ where U is a unitary mixing matrix. In this analysis the following parametrization for U has been considered where R(θ i j ) are the (complex) rotation matrices, θ i j are the mixing angles with i, j = 1, 2, 3, 4; and the order of rotation angles are considered as given in Ref. [30]. Using the above definition, the flavor change can be described as a function of the mixing matrix elements and masses in terms of the neutrino oscillation probability where α, β correspond to e, μ, τ , s; Δm 2 i j = m 2 i −m 2 j with i > j, L is the source to detector distance in 'meter' and E ν is the energy of neutrinos in 'MeV'. The oscillation probabilities for antineutrinos can be obtained by replacing mixing matrix elements U s with its complex conjugate (U * s). Since Eq. 4 is independent of the CP-violating phases as they are not observable at SBL reactor setups, the third term in Eq. 5 will be zero [30]. For a small value of mixing angle θ 14 and source to detector distance of few meters (< 100 m), the oscillation from 3 × 3 mixing parameters can be neglected. Hence, the electron antineutrino survival probability in Eq. 5 is approximated to The analysis of 3 + 1 generation is reduced to that of two flavor framework with the oscillation parameters Δm 2 41 and sin 2 2θ 14 are given by where U e4 = sin θ 14 .

Simulation procedure
The active-sterile neutrino mixing sensitivity of ISMRAN set-up will be explored at DHRUVA as well as PFBR reactor facilities. The number of neutrinos produced from the reactor depends on the thermal power. It is essential to know the fuel compositions contributing to the thermal power of the reactor. In order to estimate the number of ν e induced events produced in the detector, assumed parametrization for antineutrino flux considered in the analysis is as follows where 'a i ' is the fractional contribution from ith isotope to the reactor thermal power,'b j 's are the constant term used to fit the neutrino spectrum and E ν e is neutrino energy in MeV.
For DHRUVA reactor, we have assumed that the fractional contribution for each isotope to the reactor thermal power as given in Ref. [31] and the list of parameters used to fit the ν e spectrum due to 235 U, 239 Pu and 241 Pu are considered from Ref. [14] and for 238 U is taken from Ref. [9]. The list of parameters used in this analysis are listed in Table 1.
Similarly, for PFBR we have considered the fractional contributions due to 235 U and 239 Pu are 70% and 30%, respectively [32]. We have also considered the neutrino flux variation due to a finite size cylindrical reactor which depends on its radius and height as follows [33], where φ 0 is flux at the center of the reactor core, R and H are the physical radius and height of the cylinder, respectively, J 0 is the zeroth order Bessel function of first kind with r (0 ≤ r ≤ R) and z (0 ≤ z ≤ H ) are the vertex position of the produced neutrinos in the reactor. The leading order interaction cross-section [34] of ν e for the IBD process is given by where E e + = E ν e − (m n − m p ) is the positron energy with neglecting recoil neutron energy and p e + is the positron momentum. It can be mentioned here that the neutrino is a neutral particle and can not be detected directly. In the detector we measure the neutrino induced charged particle, for the present case it is e + . The detector resolution on true positron Here E T and E are true and the measured positron energy, respectively. The detector resolution considered in this study is in the form of σ/E ∼ 20%/ √ E. In the analysis, the neutrino induced events are distributed in terms of positron energy spectrum. We have considered total 80 bins in the e + energy range of 0-8 MeV. The number of events in i-th energy bin after incorporating the detector resolution is given as The index i corresponds to the measured energy bin and N r i corresponds to the number of reconstructed events, k is summed over the true energy of positron and n k is the number of events in k-th true energy bin. Further, K k i being the integral of the detector resolution function over the E bins and is given by The integrations are performed between the lower and upper boundaries of the measured energy (E L i and E H i ) bins. After incorporating detector energy resolution on neutrino induced events, both unoscillated and oscillated event distributions as a function of energy are shown in Fig. 2  the detection efficiency, 80% fiducial volume of the detector, 70% reactor duty cycle and, for an exposure of 1 ton-year while placing ISMRAN set-up at a distance of 13 m from the reactor core. The neutrino oscillation probability from one flavor to another not only depends on precise measurement of the source to detector distance but also on energy of neutrinos. The uncertainty in distance traveled by neutrino should be less than the oscillation wavelength in order to avoid a washout of the oscillation signal. Figure 3 shows the comparison of oscillated to unoscillated event ratios with and without incorporating the detector resolution as well as varying the source to detector distance as a function of L/E ν . Oscillated events are estimated by considering best fit values of active-sterile neutrino mixing parameter at sin 2 2θ 14 = 0.062 and Δm 2 41 = 1.7 eV 2 [25]. In Fig. 3, the black dotted line shows the event ratios in an ideal case i.e. without incorporating the detector response and at a fixed source to detector distance. The red solid line shows the event ratios by considering a cylindrical reactor core where the position of the reactor core is generated using a Monte-Carlo method and point detector. It is observed that with the variation of source to detector distance, the neutrino oscillation probability washed out as compared to fixed path length. Also by varying the path length due to random vertices of reactor core and incorporating the detector resolution on neutrino energy spectrum, the oscillation probability has further washed out as shown by blue dashed line. It is to be mentioned here that rest of the studies are performed using the randomized vertex in the reactor core and 80% fiducial volume of the detector.

Sensitivity estimation of simulated data
In order to quantify the sensitivity of ISMRAN experimental set-up to the active-sterile neutrino mixing parameters θ 14 and Δm 2 41 , we perform the statistical analysis of event distribution for an exposure of 1 ton-year. After incorporating the detector response on the number of events estimated by considering with and without active-sterile neutrino oscillation, the sensitivity to the sterile neutrino mixing parameters has been obtained by calculating the χ 2 . To determine the exclusion limit for a given confidence interval at each value of Δm 2 41 we have scanned over the values of sin 2 2θ 14 to simulate active-sterile neutrino oscillated event spectrum, and determine the boundary of the corresponding χ 2 [e.g., χ 2 = 4.61 for 90% confidence limit(C.L.)]. The χ 2 can be defined as follows [35] where n is the number of energy bins, R ex n , R th n are with oscillated and without oscillated (or theoretically predicted) events, respectively. The R th n carries the information about systematic uncertainties given by with π i n being the strength of the coupling between the pull variable ξ i and R th n . Equation (14) is minimized with respect to pull variables. Four systematic uncertainties such as 3% normalization uncertainty (including reactor total neutrino flux, number of target protons, and detector efficiency), nonlinear energy response of the detector by 1%, uncertainty in energy calibration by 0.5%. We have also considered the possibility of an uncorrelated experimental bin-to-bin systematic error of 2% which could result from insufficient knowledge of some source of background [36].

Results and discussions
The active-sterile neutrino mixing sensitivity depends on various reactor and detector parameters. As mentioned earlier the reactor related parameters are such as the thermal power, its fuel components, duty cycle and the core size. Apart from the reactor parameters, active-sterile neutrino mixing sensitivity also depends on detector mass, its fiducial volume, energy resolution, and detection efficiency. The simulation has been carried out by varying above mentioned parameters while finding the active-sterile neutrino mixing sensitivity as discussed below. The oscillation probabilities of ν e depend on the active-sterile neutrino mixing parameters such as angle and squared mass difference. The sensitivity of both these parameters depends on detector response such as resolution and efficiency. The upper limit for the active-sterile neutrino mixing angle θ 14 for an exposure of 1 ton-year is shown in Fig. 4 at 90% C.L. in the Δm 2 41 −sin 2 2θ 14 plane considering different detector resolutions, σ/E = 5%-25%/ √ E. The analysis is carried out considering reactor thermal power of 100 MW th produced from the extended reactor core and at 70% of its duty cycle. The detector is placed at a distance of about 13 m from the center of the reactor core. The detector has better active-sterile neutrino mixing sensitivity for resolution of σ/E = 5%/ √ E. It is observed that at Δm 2 41 < 0.5eV 2 , active-sterile neutrino mixing sensitivity is independent of detector resolution whereas for higher Δm 2 41 ≥ 0.5 eV 2 , the active-sterile neutrino mixing sensitivity improves by ∼ 23% for the detector resolution of 5% from 25% at Δm 2 41 = 1.0 eV 2 . Here it is to be mentioned that the precision on the Δm 2 41 is controlled by the precise measurement of energy (and also L) for individual events which depends on the resolution of the detector. Further studies are carried out considering the detector resolution of σ/E = 20% (which is the energy resolution of the PS obtained from the measurements) [26]. We have also studied the active-sterile neutrino mixing angle, sin 2 2θ 14 sensitivity by varying the detector efficiencies from 20% to 50% as shown in Fig 5. It is observed that the sensitivity on sin 2 2θ 14 improves with increase of detector efficiency for Δm 2 41 ≤ 4.0eV 2 and has less impact beyond this value. For higher values of Δm 2 41 the oscillation probability washed out. From this study, it is concluded that with better detector response, we will have better sensitivity in both the activesterile neutrino mixing parameters sin 2 2θ 14 and Δm 2 41 . In Detector efficiency 20% 25% 30% 40% 50%

Reactor core to detector distance
The distance between the reactor and the detector is not uniquely defined because of the extended reactor core such as DHRUVA reactor. Figure 6 shows the active-sterile neutrino mixing sensitivity at source to detector distances of 7 m, 10 m, and 13 m. These distances correspond the center to center distance between the reactor core and center of the detector. In our calculation neutrino vertices are generated randomly in the reactor core using MC method and assumed a point detector, for a given energy resolution, thermal power (100 MW th ) and, duty cycle of 70%. The lower limit on source to detector distance (7 m) is based on the closest accessible baseline available to place the detector. It can be observed (1 − U e4 ) 2 , expected from 1 ton-year of the data at different reactor thermal power that at Δm 2 41 = 1.0 eV 2 , the active-sterile neutrino mixing sensitivity sin 2 2θ 14 improves by ∼ 53% for the path length of 7 m from 13 m. In addition, one can maximize the event statistics and experimental sensitivity by placing the detector close to the reactor, however there is a trade-off between distance, other shielding material structures surrounding the reactor core and associated reactor background. The green dashed-dot line shows the sensitivity on mixing parameters by considering extended source as well as detector with centre to centre distance between reactor core and detector is 13 m. Both neutrinos production and their interaction point in the detector are generated on MC basis. Hence the the closest neutrinos can have a path of less than ∼ 11 m and the farthest ones oscillate for more than ∼ 15 m for given both the detector and reactor geometries. It is found that with extended detector, the sensitivity of the ISMRAN further reduces in the range of 0.3 eV 2 ≤ Δm 2 41 ≤ 4.0 eV 2 as compared to case with extended source and point detector placed at distance of 13 m.

Reactor power and duty cycle
The antineutrino flux emitted from the reactor is proportional to its operating thermal power. The DHRUVA research reactor [27] can operate at a maximum thermal power of 100 MW th , where as PFBR power reactor [28] can operate at a maximum thermal power of 1250 MW th which is an order of magnitude higher than research reactor. Figure 7 shows the comparison of exclusion limits on Δm 2 41 -sin 2 2θ 14 plane at various reactor thermal power of 60 MW th , 100 MW th , and 1250 MW th for 1 ton-year of detector exposure time at a distance of 13 m. With the increase in thermal power, there is an increase in event statistics hence increase in sensitivity of the experiment at all Δm 2 41 . Due to the operation of the nuclear reactors below than its maximum thermal output and reactor-off period, the total Reactor duty cycle 50% 70% 90% Fig. 8 The 90% C.L. exclusion limits in the Δm 2 41 − sin 2 2θ 14 plane, where sin 2 2θ 14 = 4U 2 e4 (1 − U e4 ) 2 , expected from 1 ton-year of the data at different duty cycle of the reactor ν e event statistics gets affected. Hence, sensitivity of the sterile neutrino oscillation decreases with lower duty cycle. Figure 8 shows the active-sterile neutrino mixing sensitivity of the detector at reactor duty-cycles of 50%, 70% and 90% for source to detector distance of 13 m and, 100 MW th reactor (extended source) thermal power. It is observed that the active sterile neutrino mixing sensitivity improves with the duty cycle for Δm 2 41 < 2.0 eV 2 and beyond this there is no effect. This is due to the averaging out of oscillation probabilities at higher Δm 2 41 .

Reactor fuel
Above studies are carried out assuming time averaged constant fission fraction contributions of various fuel elements as mentioned in Table 1 to the thermal power. However, the study has also been carried out considering the time evolution of the fissile fraction contribution to the reactor thermal power. The time evolution of nuclear reactor depends on its fuel cycles which can go from about a month to one or two years. In each of the new cycle, a partial or complete fuel is replaced with fresh fuel, which is typically composed by enriched 235 U. At the beginning of each reactor cycle, the ν e flux comes mainly from the fission of 235 U, with a small contribution of 238 U isotope. The neutron flux produced from the fission helps to produce 239 Pu and a small quantity of 241 Pu. Hence, as 235 U is consumed with time, then its contribution to the ν e flux decreases, whereas the contributions from 239 Pu and 241 Pu increase. However, the dominant contribution comes from the 239 Pu, which is comparable with the 235 U towards the end of each cycle [37]. Since our reactor is of CANDU type, we have used the burn-up as given in Ref. [38]. Figure 9 shows the comparison of active-sterile neutrino mixing exclusion limits between fixed fission fraction [31] and with the fission fraction variation due to burn up of the fissile element [38] at reactor thermal power 100 MW th for source to detector path lengths of 7 m and 13 m. It has been observed that the burn up variation of reactor fuel has marginal effect on the active-sterile neutrino mixing sensitivity for all considered values of Δm 2 41 for this assumed fuel cycle. However, burn-up effect may be observed for longer duration of fuel cycle. Also we have shown the effect of fuel cycle which is assumed to be 100 days in our case, on the active-sterile neutrino mixing and it shows a similar sensitivity with respect to time variation reactor burn up.

Backgrounds
The active-sterile neutrino mixing sensitivity has been obtained with inclusion of backgrounds with an assumption of signal to background ratio is one. In the analysis, three different cases of background has been considered such as the default 1/E 2 shape represents the spectral shape provided by accidental backgrounds due to contribution from intrinsic detector radioactivity, a flat distribution in antineutrino energy due to fast neutron backgrounds [39] and the combination of both these backgrounds shown in Fig. 10a. In this study, an associated 10% systematic uncertaintity is considered due to these backgrounds. Figure 10b shows the comparison of ISMRAN detector sensitivity with and without inclusion of different backgrounds. It is observed that with the contribution of both backgrounds, the active-sterile neutrino mixing angle sensitivity is further reduced by ∼ 20 % at Δm 2 41 = 1.0 eV 2 for the case of detector placed at 13 m from the 100 MW th DHRUVA reactor core.

Comparison to the other measurements
The exclusion limits at 90% C.L. on the sin 2 2θ 14 value for each Δm 2 41 obtained from the ISMRAN set-up at two different reactors are shown in Fig. 11. It can be noted here that the analysis has been carried out by smearing both the  [20], the Daya Bay [40], Bugey-3 [41], STEREO [23] and the symbol '+' is the present best fit value from the global analysis [25] are also shown for comparison at 90% C.L. The results from ISM-RAN at a distance of 13 m from DHRUVA reactor core is comparable to the NEOS results at lower Δm 2 41 < 2 eV 2 , at higher Δm 2 41 our results are comparable with the Bugey results and outperform the Daya Bay results for Δm 2 41 > 2 eV 2 . At a distance of 7 m from the reactor core and Δm 2 41 > 1 eV 2 , the results from ISMRAN are comparable with NEOS and Bugey. The ISMRAN has better sensitivity on the activesterile neutrino mixing with respect to STEREO [23]. The exclusion plot from the Daya Bay [40] experiment at lower values of Δm 2 41 = 0.1eV 2 has better sensitivity compared to ISMRAN and also other measurements. However, it is found that the active-sterile neutrino sensitivity of ISMRAN improves substantially if the measurement will be carried out at PFBR facility as shown in Fig. 11b  It is to be noted that NEOS measurements are performed at a distance of ∼ 24 m from the reactor core with thermal power of about 3 GW th . Although the reactor power of PFBR is lower compared to reactor used for NEOS measurements, the results from ISM-RAN at PFBR can give a better sensitivity as compared to other measurements. This is due to the compact core size of the PFBR facility.

Summary
In the near future, results form various SBL experiments using reactor neutrino as a source may resolve the uncertainty for the existence of light sterile neutrino hypothesis as the possible origin of the RAA and in addition it may clarify the origin of the 5 MeV distortion in the ν e energy spectra. The feasibility study on active-sterile neutrino mixing sensitivity is performed with the upcoming ISMRAN experimental set-up for an exposure of 1 ton-year employing ν e produced from the extended core of DHRUVA and compact core of PFBR reactor facility, India. The study is carried out considering both reactor as well as detector related parameters. With varying the source to detector distance of 7 m from 13 m at Δm 2 41 = 1.0 eV 2 , the sensitivity on sin 2 2θ 14 improves by twice. It has been observed that, the burn-up variation of the reactor fuel elements has very less impact on active-sterile neutrino mixing sensitivity. At reactor power of 100 MW th produced from DHRUVA reactor, the experimental set up may see the active-sterile neutrino mixing sensitivity if sin 2 2θ 14 ≥ 0.064 at Δm 2 41 = 1.0 eV 2 . On the otherhand, there is an improvement on the active-sterile neutrino mixing parameter sin 2 2θ 14 to ∼ 0.03 for the same Δm 2 41 by putting the ISMRAN detector set-up at PFBR facility. We have found the limit on active-sterile neutrino mixing parameters is of the same order as that of Bugey. However, in the range of 0.2 eV 2 < Δm 2 41 < 3.0 eV 2 , the present analysis for DHRUVA reactor predicts the same sensitivity limits as that of the results from NEOS. Further, the sensitivity at lower values of Δm 2 41 < 0.2 eV 2 , we may have better sensitivity compared to NEOS. With the ISMRAN set-up, it can be possible to verify the existence of active to sterile neutrino oscillation hypothesis as the possible origin of the RAA and, also to clarify the origin of the bump at 5 MeV in the ν e spectra.