Eur. Phys. J. C manuscript No. (will be inserted by the editor) Sensitivity to sterile neutrino mixing using reactor antineutrinos S. P. Behera a, D. K. Mishra b, and L. M. Pant Nuclear Physics Division, Bhabha Atomic Research Centre, Mumbai - 400085, India Received: date / Accepted: date arxiv:90.04746v [hep-ph 5 Jan 09 Abstract. 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+ 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) produced from reactor as a source. The so-called 3+ 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 sin θ 4 0.064 and m 4 =.0 ev at 90% confidence level for an exposure of ton-year by using neutrinos produced from the DHRUVA reactor with thermal power of 0 MW th. It is also observed that, there is a significant improvement of the active-sterile neutrino mixing parameter sin θ 4 to 0.03 at the same m 4 by putting the ISMRAN detector set-up at a distance of 0 m from the compact proto-type fast breeder reactor (PFBR) facility with thermal power of 50 MW th. PACS. XX.XX.XX No PACS code given 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 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 [ and SAGE [ 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 [3 and it is known as Gallium anomaly. The accelerator based SBL experiments such as Liquid Scintillator Neutrino Detector (LSND) [4 at a baseline of 30 m observed an unexplained excess of electronanti-neutrinos(ν e )inamuonanti-neutrinobeam. The MiniBooNE experiment also observed similar excess inν µ ν e mode[5. TherecentMiniBooNEdataareconsistent 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σ [6. There is an anomalous a e-mail:shiba@barc.gov.in b dkmishra@barc.gov.in 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. [7 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) [8. 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 [9 or of the conversion method [7, 4. 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 [5, Daya Bay [6, and Double Chooz [7 collaborations by measuring the reactor ν e energy spectrum. The distortion in energy spectra has been correlated to the reactor power [6, which may be due to the 35 U fuel [8. 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 sur-
S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos vival 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.7 m to 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 θ 4 m 4 plane covers a wide range of the sterile neutrino parameters up to sin θ 4 < 0.0 [9. Similarly, the STEREO [0 group has measured the antineutrino energy spectrum in six different detector cells covering baselines between 9 and meters 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.σ C.L. and constrains significant portions of the previously allowed parameter space at 95% confidence level [. This paper presents the results of an investigation on finding a possible mixing of a single sterile neutrino with the 3 known active neutrinos, viz. the (3 + ) model. It is the only allowed active-sterile neutrino mixing scheme[ under the assumption of 4 neutrino model. At SBL, the presence of sterile neutrinos with squared mass difference m 4(= m 4 m ) ev leads to fast oscillations resulting the reduction of reactor ν e flux, otherwise absent in the standard 3-neutrino 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 (viz. thermal power, core size, duty cycle, burn up) as well as detector response related parameters (viz. energy resolution and detection efficiency) and also at several reactor core to detector distance while constraining active-sterile neutrino mixing at an exposure of 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 Sec. and in Sec. 3, respectively. The sterile neutrino oscillation formalism is introduced in Sec. 4. The incorporation of detector resolutions on neutrino induced true events is discussed in Sec. 5. The statistical analysis using both oscillated and without oscillated events based on χ estimation is given in Sec. 6. The sensitivity to sterile neutrino mixing at an exposure of ton-yr is discussed in Sec. 7. Finally, in Sec. 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 0 plastic scintillator (PS) bars with weightofaboutton[3.thedimensionofeachpsbaris 0 cm cm cm wrapped with Gadolinium coated aluminizedmylarfoils.eachpsbariscoupledwithtwo3 Photo-multiplier tubes at both ends. The schematic of the detector set-up is shown in Fig.. 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 cm thick and then followed by cm thick of borated polyethylene (for neutrons). The detector is positioned at adistanceof 3mfromthecenterofacylindricalreactor core and can be moved closer to the core upto 7 m. The reactor has radius.5 m and height 3.03 m (defined as an extended source) [4. The reactor can operate at a maximum thermal power of 0 MW th consuming natural uranium as fuel and producing about 9 ν e /s. In future, it is planned to put the detector set-up at proto-type fast breeder reactor (PFBR) facility, IGCAR, Kalpakkam, India [5. The PFBR has dimension of about m both in radius and height (defined as a compact source), and can operate at a maximum thermal power of 50 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 proto-type ISM- RAN set-up of /5-th of the final detector volume which is under operation at DHRUVA reactor facility [4. 3 ν e DETECTION PRINCIPLE The electron-antineutrinos produce from the reactor interact with protons in the PS bars, via the Inverse Beta Decay (IBD) process, ν e +p n+e + () TheQ-valueoftheabovereactionisabout.80MeVand 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 E prompt = E νe +Q+m e c, () where E νe is the energy of electron-antineutrino. So from Eq., it is observed that there is a one-to-one correspondence between the positron energy and E νe. The neutron
S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos 3 where R(θ ij ) are the (complex) rotation matrices, θ ij are the mixing angles with i,j =,, 3, 4; and the order of rotation angles are considered as given in Ref. [7. 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 P αβ = δ αβ 4 ( m ) Re(U αi UβiU αju βj )sin ij L 4E i>j ν + ( m ) Im(U αi Uβi U αj U βj)sin ij L, E i>j ν (5) Fig.. Schematic of 0 cm cm cm plastic scintillators array with shielding materials for reactor antineutrino measurement. produced in Eq. carries a few kev s of energy and gets thermalized in collisions with protons in the PS bar. The neutron takes about 80µs in order to gets captured by protonin the PS barproducesgammaraywhich 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. 4 NEUTRINO OSCILLATION PROBABILITY WITH 3 + MODEL The sterile neutrino oscillation probabilities are based on expansion of the 3 generation Pontecorvo-Maki-Nakagawa- Sakata (PMNS) [6 matrix to 3+ generation, where 3 stands for three active neutrinos (ν e, ν µ, ν τ ) and for a sterile neutrino (ν s ). The neutrino flavors and mass eigenstates are related through ν e ν µ ντ ν s U e U e U e3 U e4 U = µ U µ U µ3 U µ4 U τ U τ U τ3 U τ4 U s U s U s3 U s4 ν ν ν 3 ν 4, (3) where U is a unitary mixing matrix. In this analysis the following parametrization for U has been considered U = R(θ 34 )R(θ 4 )R(θ 3 )R(θ 4 )R(θ 3 )R(θ ), (4) where α, β correspond to e, µ, τ, s; m ij = m i m 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 CPviolating phases as they are not observable at SBL reactor setups, the third term in Eq. 5 will be zero [7. For a small value of mixing angle θ 4 and source to detector distanceoffew meters (< 0m), the oscillationfrom 3 3 mixing parameters can be neglected. Hence, the electron antineutrino survival probability in Eq. 5 is approximated to ( ).7 m P νeν e (E,L) sin θ 4 sin 4 L, (6) E ν The analysis of 3+ generation is reduced to that of two flavor framework with the oscillation parameters m 4 and sin θ 4 are given by m 4 = m 4 m ; sin θ 4 = 4 U e4 ( U e4 ), (7) where U e4 = sinθ 4. 5 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 f(e νe ) = 4 i= 0 ( 6 a i exp j= 0 b j E j ν e ), (8) where a i is the fractional contribution from ith isotope to the reactor thermal power, b j s are the constant term
4 S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos Table. Fractional contribution of each element to reactor power and parameters used to fit the neutrino spectrum Element a b 0 b b b 3 b 4 b 5 35 U 0.58 4.367-4.577. -0.594 0.0686-0.00777 39 Pu 0.30 4.757-5.39.63-0.6596 0.078-0.003536 4 Pu 0.05.99 -.88.78-0.3343 0.03905-0.00754 38 U 0.07 4.833.97 -.83-6.76.33 -.536 used to fit the neutrino spectrum and E νe is neutrino energyinmev.fordhruvareactor,wehaveassumedthat the fractional contribution for each isotope to the reactor thermal power as given in Ref. [8 and the list of parameters used to fit the ν e spectrum due to 35 U, 39 Pu and 4 Pu are considered from Ref. [ and for 38 U is taken from Ref. [7. The list of parameters used in this analysis are listed in Table. Similarly, for PFBR we have considered the fractional contributions due to 35 U and 39 Pu are 70% and 30%, respectively [9. We have also considered the neutrino flux variation due to a finite size cylindrical reactor which depends on its radius and height as follows [30, φ = φ 0 J 0 (.405r/R) cos(πz/h) (9) 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 [3 of ν e for the inverse IBD process is given by σ IBD = 0.095 4 cm (E e p e /MeV ), () 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 energy(kinetic) spectrum is incorporated assuming a standard Gaussian form of the energy resolution: R(E,E T ) = πσ exp( (E E T) σ ). () 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 0%/ 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 N r i = k K k i (E k T)n k () 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 K k i = EHi E Li de e ( E πσ E k T E ) σ E (3) The integrations are performed between the lower and upper boundaries of the measured energy (E Li and E Hi ) 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. assuming 5% of the detection efficiency, 80% fiducialvolumeofthedetector,70%reactordutycycleand,for an exposure of ton-year while placing ISMRAN set-up at a distance of 3 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 θ 4 = 0.06 and m 4 =.7 ev [. In Fig. 3, the black dotted line shows the event ratiosin 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.
S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos 5 Counts 500 400 300 00 0 w/o oscillation with oscillation 3 4 5 6 7 8 9 Energy (MeV) Fig.. Simulated event distribution without and with activesterile neutrino oscillation after incorporating detector response. Oscillated events are estimated with m 4 =.7 ev and sin θ 4 = 0.06. Oscillated/Unoscillated.00 0.95 0.90 w/o det. reso, fixed L w/o det. reso., var L with det. reso., var L 4 5 6 7 L/E (m/mev) Fig. 3. Ratios of oscillated to unoscillated event distribution as a function L/E. Oscillated events are estimated with m 4 =.7 ev and sin θ 4 = 0.06. The black dotted line shows the ratios in ideal case i.e. without incorporating the detector response and at a fixed path length of 3 m. The red solid line shows the events ratios by considering a cylindrical reactor core (parameters mentioned earlier) without incorporating the detector resolution and the same with detector resolution is shown in blue dashed line. 6 SENSITIVITY ESTIMATION OF SIMULATED DATA In order to quantify the sensitivity of ISMRAN experimental set-up to the active-sterile neutrino mixing parameters θ 4 and m 4, we perform the statistical analysis of event distribution for an exposure of 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 χ. To determine the exclusion limit foragivenconfidence interval at each value of m 4 we have scanned over the values of sin θ 4 to simulate active-sterile neutrino oscillated event spectrum, and determine the boundary of the corresponding χ (e.g. χ = 4.6 for 90% confidence limit(c.l.)). The χ can be defined as follows [3 χ = N ( R th n=0 n Rn ex σ(rn ex) ) (4) wherenisthenumberofenergybins,rn ex,rth n arewithoscillated and without oscillated (or theoretically predicted) events,respectively.thern th carries the information about systematic uncertainties given by ( Rn th = th R n + k ) πn i ξ i +O(ξ ) (5) i=0 with πn i being the strength of the coupling between the pull variable ξ i and R th n. Equation (4) 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 %, uncertainty in energy calibration by 0.5%. We have also considered the possibility of an uncorrelated experimental bin-to-bin systematic error of % which could result from insufficient knowledge of some source of background [33. 4 [ev Detector resolution 5% % 5% 0% 5% 0.0 0. sin θ 4 Fig. 4. The 90% C.L. exclusion limits in the m 4 sin θ 4 plane, where sin θ 4 = 4U e4( U e4), expectedfrom ton-yr of the data at different PS detector resolution. 7 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
6 S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos 4 [ev Detector efficiency 0% 5% 30% 40% 50% 4 [ev Reactor power Ext. source 60MW Ext. source 0MW PFBR 0MW PFBR 50MW 0.0 0. sin θ 4 Fig. 5. The 90% C.L. exclusion limits in the m 4 sin θ 4 plane, where sin θ 4 = 4U e4( U e4), expectedfrom ton-yr of the data at various detector efficiencies. 0.0 0. sin θ 4 Fig. 7. The 90% C.L. exclusion limits in the m 4 sin θ 4 plane, where sin θ 4 = 4U e4( U e4), expectedfrom ton-yr of the data at different reactor thermal power. 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 activesterile neutrino mixing sensitivity as discussed below. 4 [ev Distance from reactor core 7m m 3m 0.0 0. sin θ 4 Fig. 6. The 90% C.L. exclusion limits in the m 4 sin θ 4 plane, where sin θ 4 = 4U e4( U e4), expected from tonyr of the data at different source to detector path lengths with 0 MW th reactor power. 7. DETECTOR RESPONSE The oscillation probabilities of ν e depend on the activesterile 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 θ 4 for an exposure of ton-yr is shown in Fig. 4 at 90% C.L. in the m 4 - sin θ 4 plane considering different detector resolutions, σ/e = 5% 5%/ E. The analysis is carried out considering reactor thermal power of 0 MW th produced from the extended reactor core and at 70% of its duty cycle. The detector is placed at a distance of about 3 m from the center of the reactorcore. The detector has better active-sterile neutrino mixing sensitivityforresolutionofσ/e =5%/ E.Itisobservedthat at m 4 < 0.5 ev, active-sterile neutrino mixing sensitivity is independent of detector resolution whereas for higher m 4 0.5 ev, the active-sterile neutrino mixing sensitivity improves by 3% for the detector resolution of 5% from 5% at m 4 =.0 ev. Here it is to be mentioned that the precision on the m 4 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 detectorresolutionofσ/e =0%(which isthe energyresolution of the PS obtained from the measurements) [3. We have also studied the active-sterile neutrino mixing angle,sin θ 4 sensitivitybyvaryingthedetectorefficiencies from 0%to 50% asshown in Fig 5. It is observedthat the sensitivity on sin θ 4 improves with increase of detector efficiency for m 4 4.0 ev and has less impact beyond this value. For higher values of m 4 the oscillation probability washed out. From this study, it is concluded that with better detector response, we will have better sensitivity in both the active-sterile neutrino mixing parameters
S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos 7 sin θ 4 and m 4. In the subsequent analyses, we have used energy resolution as σ/e = 0%/ E and efficiency of 5% unless otherwise stated. 7. 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 7m, m, and3m. These distancescorrespondsthe 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 (0 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 that at m 4 =.0 ev, the active-sterile neutrino mixing sensitivity sin θ 4 improves by 53% for the path length of 7 m from 3 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 3 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 m and the farthest ones oscillate for more than 5 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 m 4 4.0 ev as compared to case with extended detector and point detector placed at distance of 3 m. 7.3 REACTOR POWER AND DUTY CYCLE The antineutrino flux emitted from the reactor is proportional to its operating thermal power. The DHRUVA research reactor [4 can operate at a maximum thermal power of 0 MW th, where as PFBR power reactor [5 can operate at a maximum thermal power of 50 MW th which is an order of magnitude higher than research reactor. Figure 7 shows the comparison of exclusion limits on m 4 - sin θ 4 plane at various reactor thermal power of 60 MW th, 0 MW th, and 50 MW th for ton-yr of detector exposure time at a distance of 3 m. With the increase in thermal power, there is an increase in event statistics hence increase in sensitivity of the experiment at all m 4. Due to the operation of the nuclear reactors below than its maximum thermal output and reactor-off period, the total ν e eventstatistics gets affected. Hence, sensitivity of the sterile neutrino oscillation decreases with lower duty 4 [ev Reactor duty cycle 50% 70% 90% 0.0 sin 0. θ 4 Fig. 8. The 90% C.L. exclusion limits in the m 4 sin θ 4 plane, where sin θ 4 = 4U e4( U e4), expectedfrom ton-yr of the data at different duty cycle of the reactor. 4 [ev 7m ave. fraction 7m burnup with time 3m ave. fraction 3m burnup with time 3m burnup with refill 0.0 0. sin θ 4 Fig. 9. The comparison of 90% C.L. exclusion limits in the m 4 sin θ 4 plane, where sin θ 4 = 4U e4( U e4), expected from ton-yr of the data at different reactor fuel evolution. 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 3 m and, 0 MW th reactor (extended source) thermal power. It is observed that the active sterile neutrino mixing sensitivity improves with the duty cycle for m 4 <.0 ev and beyond this there is no effect. This is due to the averaging out of oscillation probabilities at higher m 4.
8 S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos 7.4 REACTOR FUEL Above studies are carried out assuming time averaged constant fission fraction contributions of various fuel elements as mentioned in Table 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 composedby enriched 35 U. At the beginning ofeach reactor cycle, the ν e flux comes mainly from the fission of 35 U, with a small contribution of 38 U isotope. The neutron flux produced from the fission helps to produce 39 Puandasmallquantityof 4 Pu.Hence,as 35 Uisconsumed with time, then its contribution to the ν e flux decreases, whereas the contributions from 39 Pu and 4 Pu increase. However, the dominant contribution comes from the 39 Pu, which is comparable with the 35 U towards the endofeachcycle[34.sinceourreactorisofcandutype, we have used the burn-up as given in Ref. [35. Figure 9 shows the comparison of active-sterile neutrino mixing exclusion limits between fixed fission fraction [8 and with the fission fraction variation due to burn up of the fissile element [35 at reactor thermal power 0 MW th for source to detector path lengths of 7 m and 3 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 4 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 0 days in our case, on the active-sterile neutrino mixing and it shows a similar sensitivity with respect to time variation reactor burn up. Counts 4 [ev 800 600 400 00 Inverse Beta signal Flat background /E background (a) 3 4 5 6 7 8 Energy (MeV) (b) 0.0 0. sin θ 4 Different background No background Flat /E Flat + /E Fig.. Signal and different background energy spectra (in panel (a)). The comparison of 90% C.L. exclusion limits in the m 4 sin θ 4 plane, where sin θ 4 = 4U e4( U e4), expected from ton-yr of the data for various background shapes. 7.5 BACKGROUNDS The active-sterile neutrino mixing sensitivity has been obtained with inclusion of backgrounds with an assumption of signal to background ratio is. In the analysis, three different cases of background has been considered such as the default /E 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 [36 and the combination of both these backgrounds shown in Fig. (a). In this study, an associated % systematic uncertaintity is considered due to these backgrounds. Figure (b) 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 0%at m 4 =.0eV for the case of detector placed at 3 m from the 0 MW th DHRUVA reactor core. 7.6 COMPARISON TO THE OTHER MEASUREMENTS The exclusion limits at 90% C.L. on the sin θ 4 value for each m 4 obtained from the ISMRAN set-up at two different reactors are shown in Fig.. It can be noted here that the analysis has been carried out by smearing both the extended source and detector volumes. Neutrinos production in the reactor core and their interaction in the detector are generated on MC basis. The upper panel shows the results obtained assuming the detector set-up placed at distance of 7 m and 3 m from the core of DHRUVA reactor and the lower panel shows sensitivity of the detector by placing at different distances in the PFBR facility. Measurements from other experiments such as, the NEOS[37, the Daya Bay[38, Bugey-3[39, STEREO[0 and the symbol + is the present best fit value from the global analysis [ are also shown for comparison at 90% C.L. The results from ISMRAN at a distance of3mfrom DHRUVA reactor core is comparable to the NEOS results at lower m 4 < ev, at higher m 4 our results are comparable with the Bugey results and outperform the
S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos 9 4 [ev 90% C.L. with ISMRAN 7m DHRUVA 3m DHRUVA 4 [ev 90% C.L. with ISMRAN 3m PFBR 0m PFBR 5m PFBR 30m PFBR sin θ 4 (a) 90% C.L. NEOS Bugey Daya Bay Global best fit (b) 0.0 0. Fig.. The comparison of 90% C.L. exclusion limits in the m 4 sin θ 4 plane, where sin θ 4 = 4U e4( U e4) between ISMRAN and other experiments. Daya Bay results for m 4 > ev. At a distance of 7 m from the reactor core and m 4 > ev, the results from ISMRAN are comparable with NEOS and Bugey. The ISMRAN has better sensitivity on the active-sterile neutrino mixing with respect to STEREO [0. The exclusion plot from the Daya Bay [38 experiment at lower values of m 4 = 0. ev has better sensitivity compared to ISMRAN and also other measurements. However, it is found that the active-sterile neutrino sensitivity of ISM- RAN improves substantially if the measurement will be carried out at PFBR facility as shown in Fig. (b). The ISMRAN results at a distance of 0 m from the reactor core are comparable to NEOS and Bugey at all values of m 4 and exclusion limits are better for m 4 < ev. It is to be noted that NEOS measurements are performed at a distance of 4 m from the reactor core with thermal power of 3 GW th. Although the reactor power of PFBR is lower compared to reactor used for NEOS measurements, the results from ISMRAN at PFBR can give a better sensitivity as compared to other measurements. This is due to the compact core size of the PFBR facility. 8 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 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 7m from 3m at m 4 =.0 ev, the sensitivity on sin θ 4 improves by twice. It has been observed that, the burn-up variation of the reactor fuel elements has very less impact on active-sterile neutrinomixingsensitivity.atreactorpowerof0mw th produced from DHRUVA reactor, the experimental set up may see the active-sterile neutrino mixing sensitivity if sin θ 4 0.064 at m 4 =.0 ev. On the otherhand, there is an improvement on the active-sterile neutrino mixing parameter sin θ 4 to 0.03 for the same m 4 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. ev < m 4 < 3.0 ev, the present analysis for DHRUVA reactor predicts the same sensitivity limits as that of the results from NEOS. Further, the sensitivity at lowervalues of m 4 < 0.eV, 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. ACKNOWLEDGMENTS We thank A. K. Mohanty, V. M. Datar and Anushree Ghosh for their helpful suggestions and useful discussions. We also thank ISMRAN group members for useful discussion. References. P. Anselmann et al. [GALLEX Collaboration, Phys. Lett. B 34, 440 (995), Phys. Lett. B 40, 4 (998).. D. N. Abdurashitov et al., Phys. Rev. Lett. 77, 4708 (996); J. N. Abdurashitov et al. [SAGE Collaboration, Phys. Rev. C 59, 46 (999). 3. J. N. Abdurashitovet al., Phys. Rev.C 73, 045805 (006).
S. P. Behera et al.: Sensitivity to sterile neutrino mixing using reactor antineutrinos 4. A. Aguilar-Arevalo et al. [LSND Collaboration, Phys. Rev. D 64, 07 (00). 5. A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration, Phys. Rev. Lett., 680 (03). 6. A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration, Phys. Rev. Lett., no., 80 (08). 7. T. A. Mueller et al., Phys. Rev. C 83, 05465 (0). 8. G. Mention, M. Fechner, T. Lasserre, T. A. Mueller, D. Lhuillier, M. Cribier and A. Letourneau, Phys. Rev. D 83, 073006 (0). 9. F. von Feilitzsch, A. A. Hahn, K. Schreckenbach, Phys. Lett. B 8, 6 (98).. K. Schreckenbach, G. Colvin, W. Gelltely, F. von Feilitzsch, Phys. Lett. B 60, 35 (985).. A. A. Hahn, K. Schreckenbach, W. Gelletly, F. von Feilitzsch, G. Colvin, B. Krusche, Phys. Lett. B 8, 365 (989).. P. Huber, Phys. Rev. C 84, 0467 (0). 3. A. C. Hayes, P. Vogel, Ann. Rev. Nucl. Part. Scien. 66, 9 (06). 4. P. Huber, Nucl. Phys. B 908, 68 (06). 5. J. H. Choi et al. [RENO Collaboration, Phys. Rev. Lett. 6, 80 (06). 6. F. P. An et al. [Daya Bay Collaboration, Phys. Rev. Lett. 6, 0680 (06). 7. Y. Abe et al. [Double Chooz Collaboration, JHEP 60, 63 (06). 8. F. P. An et al. [Daya Bay Collaboration, Phys. Rev. Lett. 8, 580 (07). 9. I. Alekseev et al. [DANSS Collaboration, arxiv:804.04046 [hep-ex. 0. H. Almazn et al. [STEREO Collaboration, arxiv:806.0096 [hep-ex.. J. Ashenfelter et al. [PROSPECT Collaboration, arxiv:806.0784 [hep-ex.. S. Gariazzo, C. Giunti, M. Laveder and Y. F. Li, JHEP 706, 35 (07). 3. D. Mulmule et al., arxiv:806.044 [physics.ins-det. 4. S. K. Agarwal, C. G. Karhadkar, A. K. Zope and Kanchhi Singh, Nuclear Engineering and Design 36, 747 (006). 5. S. C. Chetal et al., Nuclear Engineering and Design 36, 85 (006). 6. Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 8, 870 (96). 7. A. Palazzo, JHEP 3, 7 (03). 8. L. Zhan, Y. Wang, J. Cao and L. Wen, Phys. Rev. D 78, 3 (008). 9. M. Schneider, Science & Global Security, 7, 36 (009). 30. Nuclear Reactor Engineering, Vol., Samuel Glasstone, Alexander Sesonske, Von-Nostrand, 994. 3. P. Vogel and J. F. Beacom, Phys. Rev. D 60, 053003 (999). 3. M. C. Gonzalez-Garcia and M. Maltoni, Phys. Rev. D 70, 0330 (004). 33. P. Huber, M. Lindner, T. Schwetz and W. Winter, Nucl. Phys. B 665, 487 (003). 34. C. Giunti, X. P. Ji, M. Laveder, Y. F. Li and B. R. Littlejohn, JHEP 7, 43 (07). 35. M.W. Francis, C.F. Weber, M.T. Pigni and I.C. Gauld, Report ORNL/TM-04/464, Oak Ridge National Laboratory, Oak Ridge, Tennessee 3783-683, U.S.A., September 04. 36. K.M. Heeger,B. R.Littlejohn, H.P.MummandM. N.Tobin, Phys. Rev. D 87, 073008 (03). 37. Y. J. Ko et al., Phys. Rev. Lett. 8, 80 (07). 38. P. Adamson et al. [Daya Bay and MINOS Collaborations, Phys. Rev. Lett. 7, 580 (06). 39. Y. Declais et al., Nucl. Phys. B 434, 503 (995).