The Reactor Anomaly after Daya Bay and RENO

Emilio
Ciuffoli^{*}^{*}*ciuffoli@ihep.ac.cn, Jarah
Evslin^{†}^{†}†jarah@ihep.ac.cn, Hong Li^{‡}^{‡}‡hongli@ihep.ac.cn

1) TPCSF, IHEP, Chinese Academy of Sciences

YuQuan Lu 19(B), Beijing 100049, P.R.China

2) Key Laboratory of Particle Astrophysics, IHEP,

Chinese Academy of Sciences, P.O.Box 918-3, Beijing 100049, P.R.China

3) National Astronomical Observatories, Chinese Academy of Sciences,

Beijing 100012, P.R.China

Abstract

Gallium and short baseline reactor neutrino experiments indicate a short-distance anomalous disappearance of electron antineutrinos which, if interpreted in terms of neutrino oscillations, would lead to a sterile neutrino mass inconsistent with standard cosmological models. This anomaly is difficult to measure at 1 km baseline experiments because its disappearance effects are degenerate with that of . The flux normalization independent measurement of at Daya Bay breaks this degeneracy, allowing an unambiguous differentiation of 1-3 neutrino oscillations and the anomalous disappearance at Double Chooz and RENO. The resulting anomaly is consistent with that found at very short baselines and suggests a downward revision of RENO’s result for . A MCMC global analysis of current cosmological data shows that a quintom cosmology is just compatible at with a sterile neutrino with the right mass to reproduce the reactor anomaly and to a lesser extent the gallium and LSND/MiniBooNE anomalies. However models in which the sterile neutrino acquires a chameleon mass easily satisfy the cosmological bounds and also reduce the tension between LSND and KARMEN.

July 14, 2019

## 1 Introduction

During the past 20 years, the gallium experiments GALLEX [1] and SAGE [2] have observed^{1}^{1}1Here we use the latest calibrations [3]. 145% [4] less neutrinos than expected from radioactive sources inside of the detector. While Lorentz violating [5], CP violating [6, 7] and CPT violating [8] neutrino anomalies have disappeared in the last few months [9, 10, 11], the gallium disappearance anomaly has been corroborated by a new calculation of reactor neutrino fluxes [12]. This new calculation suggests a 3.5% increase in theoretical fluxes from nuclear reactors, leading to a total flux which is about 5.72.3% higher than that observed at short baseline reactor experiments [13]. We will refer to this missing reactor neutrino flux as the reactor anomaly. The reactor and gallium anomalies stand out among neutrino anomalies because they occur in a relatively accessible environment^{2}^{2}2This is in contrast with the low energy solar neutrino deficit [14] and the horizontal flux excess [15]., which can and will be directly probed by future experiments, and yet their simplest interpretation in terms of oscillations of a 1 or more eV sterile neutrino [16] is in conflict with standard cosmology (for a recent analysis see Ref. [17]).

Not all reactor neutrino disappearance is anomalous, reactor neutrinos also disappear as a result of oscillation between the three standard neutrino flavors. At distances below about 100 meters, all such oscillation is well within the experimental errors of current measurements. Above 100 meters 1-3 neutrino oscillation is relevant and above about 2 kilometers 1-2 oscillation is also significant. As the 1-3 mixing angle was poorly measured until quite recently, it has not been possible to disentangle the reactor anomaly from 1-3 oscillation in 1 km baseline reactor experiments. This situation changed 3 months ago when the Daya Bay collaboration measured by comparing neutrino fluxes observed at identical detectors at various distances from an array of reactors [18]. Following the proposal of Ref. [19] this allowed them to determine independently of any assumptions about the overall flux normalization. Recently, with 3 months more data, the Daya Bay collaboration has determined yet more precisely [20].

In this note we will combine the results of the reactor neutrino experiments Double Chooz [21], Daya Bay [18, 20] and RENO [22, 23] to simultaneously estimate the reactor anomaly and . In particular the flux independent mixing angle obtained by Daya Bay allows for an unambiguous separation of standard mixing and anomalous disappearance at the other experiments. As has already been seen in the global fits in Refs. [24, 25], greater reactor anomalies imply smaller mixing angles. Using the reactor anomaly measured at short baseline experiments we will find that the preferred mixing angle is lower than RENO’s result. It is easy to extend our analysis to Palo Verde [26] and Chooz [27] experiments and to gallium experiments; however, as we will see in Fig. 2, this extension does not change our results appreciably, so we will focus on the Daya Bay, RENO and Double Chooz results, confronting them with the reactor anomaly discussed in [13]. We did not use the latest data release from Double Chooz [28] which includes 142 days in which nearly doubled with respect to the first 86 days and is highly dependent upon the analysis used.

We will then consider theoretical models which may explain this anomaly consistently with cosmological constraints. We will see that a quintom dark energy model with a massive sterile neutrino is consistent with these anomalies at about 2. On the other hand a neutrino dark energy model [29] in which a sterile neutrino mass is proportional to the density of its environment [30] easily satisfies cosmological constraints and also reduces the tension between the positive appearance data at LSND and MiniBooNE and the negative results at KARMEN.

After this paper was completed we received the preprints [31, 32] which combine data from the same three experiments. Their analysis differs from ours in the treatments of the reactor anomaly and the total normalization of the theoretical reactor fluxes. Our analysis is 2-dimensional, simultaneously fitting for the anomalous disappearance and 1-3 oscillation induced electron antineutrino disappearance.

## 2 Disappearance Results from RENO

The RENO experiment consists of two detectors, one near and one far, placed near an array of six nuclear reactors. RENO’s results [22] have been posted on the archive in two versions and have been extended at a recent talk [23] at the conference TURN. As they do not provide information about the reactor fluxes, or even about the calibration of the machines as in Ref. [33], no single reference is sufficient to reproduce the claimed results. However, by combining information from, for example, the most recent arXiv version and the TURN talk we were able to reproduce the best mixing angle and flux deficit with an error of about one part in one thousand and to reproduce the value to within the widths of the lines in their graph.

Following the analysis of Ref. [22], we have included all neutrinos with prompt energies up to 12 MeV, however it has been claimed in Ref. [34] that the deficit of neutrinos above 6 MeV is largely unrelated to neutrino oscillation and serves to artificially deflate the best fit mixing angle. A adapting a 6 MeV cut would make the RENO result more compatible with the rate only analysis of the new Double Chooz data [28] which we have not used, increasing our best fit values of both and the reactor anomaly.

Fig. 3 of Ref. [22] provides the oscillation probability, for a given value of , as a function of the weighted baseline , which is the average distance traveled by a neutrino detected at the th detector. Ignoring neutrino oscillations, which is reasonable within the errors at this baseline, the survival probability of neutrinos arriving at the th detector is

(2.1) |

where is the fraction of the neutrino flux observed at the th detector which originated at the th reactor, given in table 1 of version 1 of [22] and also at the talk [23], is the difference between the squared masses of the first and third neutrino flavors and is the fractional energy distribution of neutrinos emitted from the th reactor. This is not equal to the quantity which is plotted in Fig. 3, which is the survival probability at the weighted baseline

(2.2) |

Nonetheless, as we will see the information contained in Fig. 3 is essential to reproduce RENO’s analysis and to extrapolate the observed reactor anomaly.

To determine which values of the anomalous neutrino deficit and best fit RENO’s data, we used the formula (2.1). We used the fractional fluxes reported in table 1 of version 1 of [22] or equivalently page 6 of [23] and the normalization of the expected fluxes from figure 3 of version 2. There are several inequivalent ways of estimating the energy distribution , without making use of the unreleased reactor data. One can add 780 keV to the prompt energy spectrum reported in Fig. 4

(2.3) |

However the statistical errors are large and the quality of the energy calibration is unknown. We found that in analyzing both RENO and Daya Bay data our fits matched those of the experimental groups more closely when we used the sample neutrino energy spectra of Fig. 1.9 of Ref. [35].

To estimate errors we used the pull parameter method of the second version of Ref. [22], with given by their Eq. (2). However as we did not wish to impose any theoretical bias upon the overall reactor flux, we did not restrict its value. We used the systematic errors reported in the second paper. When we optimized the reactor flux within the region that they described, our resulting curve agreed with that shown on the top of RENO’s Fig. 3 to within the widths of their lines.

The flux normalization and are two free parameters. A flux only analysis, such as that performed by RENO, has only two data points, the fluxes observed at the two detectors. Therefore the effective number of degrees of freedom is , reflecting the familiar fact that with two parameters one can generically fit two unknowns. Thus it is of no surprise that the minimum is equal to 0. This is true in their case, in which the reactor flux is optimized and then is fit and it is true in our case, in which we wish to use RENO’s data to fit both the flux deficit and .

Using RENO’s near and far detector fluxes reported in their Fig. 3 to determine the neutrino flux deficit and we obtain essentially the same fit as was obtained by the RENO collaboration, seen as a horizontal stripe in Fig. 1. The stripe is horizontal because, without knowledge of the theoretical flux, the reactor anomaly cannot be estimated.

However, RENO has tentatively reported the theoretical flux. Most recently, in the talk [23] they have provided a preliminary estimate of the anomalous flux deficit. Using their value of they estimate this deficit to be 6%. We have performed a 2-dimensional fit on their new data, providing confidence intervals for both the anomalous flux deficit and . Our result is the ellipse in the left panel of Fig. 2. The center of the ellipse indeed corresponds to the best fit found by the RENO collaboration. Note that a greater anomaly means that less electron neutrinos are expected, because for example they have converted into sterile neutrinos, and so the deficit at both detectors is reduced. This can be compensated at the far detector by decreasing , which explains the inclination of the ellipse and was already evident in the global analyses of Refs. [24, 25].

## 3 Disappearance Results from Double Chooz

Although next year Double Chooz will have both a near and far detector, for now the data contains results from a single detector, which is compared with theoretical fluxes. The Double Chooz collaboration analyzed their data in a number of different ways, arriving at different best fit mixing angles. For the sake of uniformity, we will perform a pure rate analysis, using only the number of neutrinos observed, the theoretical flux and the baseline, even though this leads to large systematic errors [21]. As there is only one detector and both reactors are the same distance from that detector, in this case the spectrum dependence in the survival probability only enters via the energy-weighted quantity

(3.1) |

which is easily extracted from Ref. [21]. Thus in this case we are not required to make any arbitrary assumptions about the isotope fractions in the reactors.

The Double Chooz paper Ref. [21] was written after the upward revision in theoretical fluxes in Ref. [12], in fact the upward revision was motived by its use at Double Chooz. Nonetheless Double Chooz did not normalize its data using the new fluxes, but instead using a hypothetical version of the Bugey4 experiment [36] with modified abundances for the four main fissioning isotopes to fix the flux normalization by hand. However they also reported the theoretical calculation for the flux using Ref. [12], whose comparison with the observed flux we have used to determine their neutrino deficit.

In the analysis we use only a single data point, the total neutrino flux observed at the detector. We fit this data point using 2 parameters, the anomalous flux deficit and . Therefore we have one constraint on two unknowns, and so there is a 1-dimensional curve of solutions with . To reproduce the value of reported by Double Chooz in Ref. [21] one would need to also impose an anomalous flux deficit equal to that measured by Bugey4, which according to [13] was about 6% but according to [21] for the isotope fractions at the Chooz reactors was about 8%. While Bugey4’s flux measurement does represent a typical value of the reactor anomaly and while its errors are relatively small, we will not adopt this strategy. The flux deficit will instead be left as a free parameter in all plots, which can then be compared with the reactor anomaly determined not only at Bugey4 but at all very short baseline reactor experiments and if desired gallium experiments and RENO.

We have performed our analysis without using RENO’s preliminary theoretical flux. Our results are illustrated in Fig. 1. The right panel is a combined fit from the three experiments. The corresponding ellipse is nearly horizontal, implying that the reactor anomaly and are nearly decoupled. This decoupling is a result of the fact that the only information on the absolute normalization of the flux comes from Double Chooz, whose errors are quite large. Therefore the preferred value of the reactor anomaly is about equal to that of Double Chooz, and the preferred value of is nearly a weighted average of that obtained at the three experiments. If in this case one ignores the small correlation between the reactor anomaly and the mixing angle, then the analysis will be reasonably independent of the flux normalizations. With this approximation we would reproduce the analysis of these three experiments performed in Ref. [31].

However if we include the preliminary flux normalization reported by the RENO collaboration in their talk [23] then we find by far the strongest evidence yet for a reactor anomaly, as seen in Fig. 2.

One may try to incorporate the above 1 km baseline fits of the mixing angle and the reactor flux deficit with complimentary data from the flux normalization-independent analysis by Daya Bay [18], the reactor anomaly at short baseline reactor experiments [13], the neutrino deficit observed at gallium experiments [4] and global fits of experimental data more than 10 months old [24, 25].

The Daya Bay data is easily analyzed in this framework. As it depends only upon the mixing angle , the confidence intervals translate into a simple horizontal stripe on our Figs. 1 and 2. Similarly the short distance reactor and gallium anomalies are independent of , and so they correspond to vertical stripes. The global analyses are somewhat more difficult to adapt to our setting, in which the reactor flux normalization is not fixed. However Ref. [24] produces confidence intervals for for both the new reactor fluxes of Ref. [12] and the old normalization, which corresponds to a 3-3.5% lower flux. An adaptation of the global analysis of Ref. [25] to an arbitrary reactor anomaly is more difficult. The mixing angle cited in their abstract includes a fit of short distance reactor data which effectively already changes the flux normalization. Only the mixing angle which they quote in their Eq. (3), in which short baseline reactor data is excluded, uses the flux normalization of Ref. [12]. It is consistent with the analysis of Ref. [24] for this flux choice.

In Fig. 3 we have plotted the 1 and 2 allowed regions obtained by combining the results of Double Chooz, Daya Bay, RENO, Chooz, Palo Verde, very short baseline experiments and gallium experiments with and without RENO’s preliminary theoretical flux reported in the talk [23]. This figure illustrates that kilometer baseline reactor experiments may provide the most convincing evidence yet for a reactor anomaly. For now this result is highly dependent upon preliminary normalization results from RENO, a situation which will be remedied when Daya Bay’s flux normalization analysis is complete.

## 4 Theoretical models of the anomaly

### 4.1 The Sterile Neutrino Solution

The reactor and gallium anomalies are disappearance anomalies. Less electron antineutrinos are observed than are predicted by theoretical models. One consistent explanation is that the models are wrong. In particular, the theoretical calculations generally have uncertainties of order 2.5%, and so even if a deficit of 5% is confirmed the tension with the theoretical calculation will only be 2. In the remainder of this note we will consider the other possibility, that the anomalies indicate new physics. The central question is then, what kind of interaction may be responsible for the disappearance?

The simplest explanation would be an interaction. Perhaps the antineutrinos are simply absorbed as they travel, into the vacuum or the Earth or a dark matter or dark energy field. In this case one would expect the neutrino density to fall exponentially with the baseline. This is not observed. The anomaly at the 100 kilometer distance scales probed by KamLAND is too small to be measured, while it is been observed at the 5% level at numerous experiments with baselines below 100 meters [13]. Therefore the disappearance cannot be caused by simple absorption. The anomaly appears to saturate at a short distance.

How can the anomaly saturate? How can the neutrinos know, after traveling 20 meters, that they have already been absorbed sufficiently and now they should no longer be absorbed? The only apparent causal explanation is that the neutrino beam travels coherently with another beam, which contains the information concerning how many neutrinos have disappeared. The simplest realization of this idea is that the electron neutrinos oscillate into another kind of neutrino which then oscillates back, and so the neutrino density reaches when the ratio of these two kinds of neutrinos reaches a critical value. Either the new neutrino mass is above the electroweak scale or else LEP can exclude that it interacts weakly. Strong and electromagnetic interactions for the new neutrino are already excluded by the fact that it can travel through the Earth as far as, say the KamLand experiment. Therefore the new neutrino is sterile with respect to standard model gauge interactions, although of course dark force interactions are not ruled out.

A new kind of neutrino may sound like a big assumption. However, if any fermion is neutral under the standard model gauge symmetries then Yukawa interactions with neutrinos and for example the Higgs field are marginal and are not forbidden and so are expected. Once such an interaction exists, whatever the nature of the fermion, physicists will call it a sterile neutrino. This is just one possible coupling, that corresponding to the lowest dimensional operator. To distinguish it from other couplings, one needs an energy dependence analysis of the anomaly. Such an analysis is just barely beyond the reach of current Daya Bay data. As soon as they understand their nonlinear detector response away from 2.5 MeV and have accumulated more data, it may be tested.

What properties must the sterile neutrino have? To explain the anomaly is quite easy. To have the right percentage of disappearance one need only fix the mixing angle to be of order degrees. The shortest distance reactor and gallium anomalies indicate that the disappearance has saturated by the time the neutrinos have traveled 10 meters, indicating a sterile neutrino mass above 1 eV. This is all that is needed to be consistent simultaneously with the reactor and gallium anomalies at short baselines and also with limits on neutrino disappearance at longer baselines. If one allows mixing between sterile and muon neutrinos then this can also explain the anomalous oscillations at LSND [6] and MiniBooNE [7, 10], although some tension remains with MiniBooNE’s lowest energy bins.

### 4.2 Cosmological Constraints

The large sterile neutrino mixing angle required by these anomalies implies that sterile neutrinos will be in thermal equilibrium in the early universe (for recent analyses see [37, 38]). After they decouple, the sterile neutrino temperature and therefore energy density may be calculated as a function of time. The energy density impedes the formation of large scale structure (LSS) at scales below the neutrino free streaming length. This can be compensated by increasing the matter fraction . Type Ia supernova data then prefer a lower value of the dark energy equation of state whose extra acceleration compensates for the deceleration caused by the additional matter.

However a larger means that there was less time before recombination. To see this, use the fact that since then the Universe has been matter dominated together with the Friedmann equation to find

(4.1) |

As is fixed by a simple thermodynamic argument together with the CMB temperature, one sees that as increases, decreases. As the primordial plasma lasted for less time , its perturbations do not travel as far before recombination, and so the baryon acoustic oscillation (BAO) scale is smaller than it would have been with massless neutrinos.

This change in absolute size of the BAO scale is not in itself in contradiction with observations as it is not the absolute scale which is measured, but rather the angular scale and redshift depth. The measured scale can appear larger, to preserve agreement with observations of the angular size of the BAO peak [39, 40], if the observed galaxies are closer. As the redshifts of these galaxies are known, they will be closer if the function is increased. Equivalently the BAO angular scale can be recoved if is increased, which at each redshift is equivalent to increasing the Hubble scale

(4.2) |

As it is clear that if increases then increases, decreasing the time that light has traveled from these galaxies. This results in an increased angular size of BAO features. Thus a larger value of makes a smaller BAO scale appear bigger and so makes a large value of compatible with BAO observations. In summary, while reconciling massive neutrinos with supernova requires a smaller value of , baryon acoustic oscillations require a larger value. This tension provides the strongest cosmological bound on neutrino masses.

Some of this tension is relieved by the presence of an additional flavor of relativistic neutrino at high redshift, increasing the number of effective flavors to about 4. Even if the extra flavor is massless, the preferred dark matter density will increase by about 17% [17] in order to preserve the redshift at matter-radiation equality

(4.3) |

which is well determined by the CMB power spectrum. Increasing from 3 to 4 increases the denominator in (4.3) by about 12%, compensating for most of the increase in the numerator. Now when the sterile neutrino masses are turned on, the additional dark matter already somewhat compensates for the free-streaming suppression of medium scale structure.

### 4.3 Satisfying the cosmological constraints

We employ a modified version of CosmoMC [41] to constrain the sterile neutrino mass by performing a Markov Chain Monte Carlo global fitting analysis with the current observations, which includes the 7-year WMAP temperature and polarization power spectrum [42], the Union2.1 sample of type Ia supernova including systematic errors [43], BAO distance ratios from SDSS DR7 galaxies [40], and [44] as recently measured by the Hubble space telescope. Since the neutrino masses and the dark energy equation of state lead to correlated effects, we have compared the sterile neutrino mass limits obtained within different dark energy models: the CDM model and a model with a time-dependent dark energy equation of state

(4.4) |

In the dynamical dark energy models we eliminate the divergences in dark energy perturbations that occur when the dark energy equation of state crosses -1 using the algorithm suggested in Refs. [45].

First, we determined the compatibility of neutrinos of the desired mass with our cosmological model. While the reactor anomaly can be caused by sterile neutrinos as light as eV, such neutrinos would be irrelevant at the short distances probed by the gallium anomaly. Similarly eV neutrinos yield a somewhat better fit for the LSND and MiniBOONE anomalies than eV neutrinos, considering the bounds on the mixing angles from other experiments. Assuming that there are 3 massless neutrinos and 1 (sterile) massive neutrino, we confront the model with current data, and in Fig. 4 we plot the -dimensional frequentist probability distribution of the sterile neutrino mass. By fitting with different dark energy parametrizations, we obtain constraints on the sterile neutrino mass of , and respectively for CDM and for time-evolving dark energy models. In the CDM model, the eV sterile neutrino is compatible within and the eV model, preferred by neutrino experiments, is excluded by more than , while the tension between these observations and 1 eV massive sterile neutrinos can be reduced within the dynamical dark energy models.

To better understand the cosmological consequences of massive sterile neutrinos, we then restricted our attention to 3 flavors of light active neutrinos and a single 1 eV sterile neutrino. We found that cosmology prefers a larger dark matter density and a lower value of than in the massless case with no sterile neutrino. Our mean preferred values and constraints are , which can be compared with the no sterile neutrino case .

In Fig. 5 we plot the cross correlation contours of and with and without sterile neutrinos. The red solid and the black dotted curves bound the 1 and 2-allowed regions with and without the sterile neutrinos respectively. The red contours are shifted towards a higher and a more negative value of as compared with the black contours. This means that a higher value of is preferred at low redshifts, where BAO data dominates, and a lower value at redshifts beyond the reach the SDSS large scale structure surveys, where supernova data dominates. The redshift dependence that we found is much stronger than that found in Ref. [46], reflecting the fact that the Union2.1 supernova dataset that we used contains about twice as many supernova at redshift than the Union2 dataset [47] used in their analysis.

Note that our solution to the cosmological constraints is deep inside of the quintom regime [48], as is far from zero. It may be that it only provides a fair fit to the data because BAO constraints are weak at high redshifts, beyond the distances at which luminous red galaxies have been surveyed. However the steeply sloped dark energy equation of state in such models strongly favor a very large BAO acoustic scale at high redshifts, a prediction which will soon be tested.

As one might expect, more radical solutions can be yet more compatible with the cosmological data, such as a high sterile neutrino initial asymmetry [38] which can prevent sterile neutrinos from thermalizing [49] or a modified theory of gravity [50]. The suppression of small scale structure formation can be compensated by including extended objects in the early universe, such as cosmic strings [51] or giant monopoles [52], obviating the need for an increased dark matter fraction. Nonstandard matter effects may also be responsible for these anomalies, for example in Ref. [53] the authors used sterile neutrinos with masses within cosmological bounds subjected to nonstandard matter effects to explain the LSND and MiniBooNE anomalies.

The tension may also be eliminated altogether if the sterile neutrino couplings are environmentally-dependent [54], for example in a chameleon model in which the sterile neutrino mass is proportional to the density of its environment. It was argued in Ref. [30] that, in the case of the couplings here, such models are consistent with both precision gravity and equivalence principle bounds. LSND, MiniBooNE, reactor anomaly and gallium anomaly experiments all have a significant fraction of the neutrino baseline inside of dense media like the Earth and detector shielding. If the sterile neutrino mass is proportional to the background density and the sterile neutrino mass is of order 2 eV inside of the Earth then in space it will be relativistic at all times since decoupling from the primordial plasma. Thus, so far as big bang nucleosynthesis, structure formation and CMB and baryon oscillations are concerned, the sterile neutrino is massless and all cosmological bounds are easily satisfied. In fact, the addition of a massless neutrino in general improves the cosmological fits compared with with the CDM model.

The fact that about half of the LSND and MiniBooNE baseline is air means that the required neutrino mass will increase to about 2 eV, which is still compatible with constraints from beta decay. However, as the KARMEN baseline has a higher fraction of air, the historical tension between LSND and KARMEN will be reduced. This dependence upon the material in the baseline may be quite easy to test. A minimum amount of shielding is always required around a detector, however this can be thinner in the path followed by the neutrinos. Of course, nuclear power plants themselves always have shielding, and it may be that this shielding is already sufficient for the anomalous oscillation to occur. In this case, gallium experiments may be a more promising setting to test the environmental dependence of sterile neutrino masses.

## 5 Conclusions

Combined with cosmological evidence for an additional relativistic degree of freedom before recombination, the very short distance reactor, gallium and LSND/MiniBooNE neutrino disappearance anomalies provide an ever more substantial argument for a heavy sterile neutrino. While neutrino disappearance has also been observed at baselines ranging from 100 m to 2 km, it has not been so far possible to determine how much of this disappearance is attributable to 1-3 neutrino flavor oscillations and so it has not been possible to use these experiments to support or refute the reactor anomaly hypothesis.

With the flux normalization independent precision measurements of the 1-3 mixing angle at Daya Bay and RENO the degeneracy between the reactor anomaly and has been broken. Therefore, although these experiments provide only moderately precise measurements of the flux normalization themselves, their determinations of may be applied to all neutrino oscillation experiments with baselines longer than 100 meters to determine just how much 1-3 oscillation has occurred, and so what remaining signal is anomalous. This could allow a more precise reanalysis at accelerator experiments such as T2K and MINOS. In this note we applied them to the simpler setting of the short baselines reactor experiments themselves. By combining the disappearance results of Double Chooz, Daya Bay and RENO we were able to simultaneously fit and the anomalous disappearance. We found that these experiments certainly do not exclude a reactor anomaly and indeed lend it some additional support.

Our results are still strongly dependent upon assumptions concerning the theoretic flux at RENO and are bounded by the statistical error at Daya Bay. However both of these shortcomings will remedied by data soon to be released by these experiments, and then we expect that an application of such 2-dimensional fits to the new data will provide a somewhat more precise determination of the reactor anomaly, complimenting the scheduled very short distance measurements.

The theoretical errors in the normalized reactor fluxes remain quite large. And so even if the existence of a reactor anomaly may be convincingly demonstrated in the near future, it will be much more difficult to determine whether it is attributable entirely to error in this calculation or else to new physics. Indeed it is quite plausible that a convincing determination will require inputs from cosmology, which are only weakly dependent upon mixing angles and mass differences but are rapidly providing a more convincing case for a sterile neutrino and potentially, in a few years, a value for at least the total neutrino mass. Therefore we have in this paper also attempted to analyze the consistency of sterile neutrino models with cosmological constraints.

Until last year, CP violating anomalies suggested that there be at least 2 flavors of sterile neutrinos, which were in general strongly disfavored by cosmological constraints. Evidence for CP violation in neutrino anomalies has since weakened appreciably, and so now models with 1 flavor of sterile neutrino provide reasonable fits to experimental data. The best fits are obtained for neutrino masses near eV. As we have reviewed, in a cosmology such massive sterile neutrinos are essentially excluded, and even allowing a constant dark energy equation of state they cannot be made simultaneously compatible with BAO and supernova observations.

However we observe that the various constraints on neutrino masses are relevant at different redshifts. BAO constraints are tight at low redshifts and supernova constraints are tighter at higher redshifts, in particular using the Union2.1 dataset which includes many more supernova at than its predecessors. Therefore, at least for now, a strong redshift dependence in significantly reduces the tension between these cosmological constraints, while making strong predictions for future, higher redshift, large scale structure surveys.

Further we found that the tension may be eliminated entirely by noting that all evidence for massive neutrinos occurs inside of dense media whereas constraints on neutrino masses occur in much lower density environments. Therefore models in which the sterile neutrino mass depends on the background density are consistent with both. These models have the additional attractive feature that they explain the observed anomalous oscillation at LSND despite the negative results at KARMEN. While KARMEN’s total baseline is more than half of that of LSND, LSND has appreciably more dense shielding. The density weighted baseline at LSND is appreciably higher, increasing the amount of oscillation expected at LSND with respect to KARMEN. More importantly, such models are easy to test, in particular at the next generation of small experiments designed to test the reactor anomaly. Changing the ratio of air to shielding between the reactors and detectors when possible, or otherwise changing the density of the shielding will in general affect the neutrino oscillation probability in such models, allowing them to be easily excluded in the near future.

## Acknowledgement

We would like to thank Xinmin Zhang for useful discussions and advice in every phase of this project. We are also grateful to Thierry Lasserre and Sören Jetter for patient discussions and correspondence. We acknowledge the use of the Legacy Archive for Microwave Background Data Analysis (LAMBDA). Support for LAMBDA is provided by the NASA Office of Space Science. The calculation was performed at the Deepcomp7000 of Supercomputing Center, Computer Network Information Center of Chinese Academy of Sciences. JE is supported by the Chinese Academy of Sciences Fellowship for Young International Scientists grant number 2010Y2JA01. EC is supported in part by the NSF of China. HL is supported by the National Science Foundation of China under Grant Nos. 11033005, the 973 program under Grant No. 2010CB83300 and the Chinese Academy of Science under Grant No. KJCX2-EW-W01.

## References

- [1] P. Anselmann et al. [GALLEX. Collaboration], “First results from the Cr-51 neutrino source experiment with the GALLEX detector,” Phys. Lett. B 342 (1995) 440. W. Hampel et al. [GALLEX Collaboration], “Final results of the Cr-51 neutrino source experiments in GALLEX,” Phys. Lett. B 420 (1998) 114.
- [2] D. .N. Abdurashitov, V. N. Gavrin, S. V. Girin, V. V. Gorbachev, T. V. Ibragimova, A. V. Kalikhov, N. G. Khairnasov and T. V. Knodel et al., “The Russian-American gallium experiment (SAGE) Cr neutrino source measurement,” Phys. Rev. Lett. 77 (1996) 4708. J. N. Abdurashitov et al. [SAGE Collaboration], “Measurement of the response of the Russian-American gallium experiment to neutrinos from a Cr-51 source,” Phys. Rev. C 59 (1999) 2246 [hep-ph/9803418]. J. N. Abdurashitov, V. N. Gavrin, S. V. Girin, V. V. Gorbachev, P. P. Gurkina, T. V. Ibragimova, A. V. Kalikhov and N. G. Khairnasov et al., “Measurement of the response of a Ga solar neutrino experiment to neutrinos from an Ar-37 source,” Phys. Rev. C 73 (2006) 045805 [nucl-ex/0512041]. J. N. Abdurashitov et al. [SAGE Collaboration], “Measurement of the solar neutrino capture rate with gallium metal. III: Results for the 2002–2007 data-taking period,” Phys. Rev. C 80 (2009) 015807 [arXiv:0901.2200 [nucl-ex]].
- [3] F. Kaether, W. Hampel, G. Heusser, J. Kiko and T. Kirsten, “Reanalysis of the GALLEX solar neutrino flux and source experiments,” Phys. Lett. B 685 (2010) 47 [arXiv:1001.2731 [hep-ex]].
- [4] C. Giunti and M. Laveder, “Statistical Significance of the Gallium Anomaly,” Phys. Rev. C 83 (2011) 065504 [arXiv:1006.3244 [hep-ph]].
- [5] P. Adamson et al. [ MINOS Collaboration ], “Measurement of neutrino velocity with the MINOS detectors and NuMI neutrino beam,” Phys. Rev. D76 (2007) 072005. [arXiv:0706.0437 [hep-ex]]. T. Adam et al. [ OPERA Collaboration ], “Measurement of the neutrino velocity with the OPERA detector in the CNGS beam,” [arXiv:1109.4897 [hep-ex]].
- [6] A. Aguilar-Arevalo et al. [LSND Collaboration], “Evidence for neutrino oscillations from the observation of anti-neutrino(electron) appearance in a anti-neutrino(muon) beam,” Phys. Rev. D 64 (2001) 112007 [hep-ex/0104049].
- [7] A. A. Aguilar-Arevalo et al. [The MiniBooNE Collaboration], “A Search for electron neutrino appearance at the eV scale,” Phys. Rev. Lett. 98 (2007) 231801 [arXiv:0704.1500 [hep-ex]]. A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration], “Unexplained Excess of Electron-Like Events From a 1-GeV Neutrino Beam,” Phys. Rev. Lett. 102 (2009) 101802 [arXiv:0812.2243 [hep-ex]]. A. A. Aguilar-Arevalo et al. [The MiniBooNE Collaboration], “Event Excess in the MiniBooNE Search for Oscillations,” Phys. Rev. Lett. 105 (2010) 181801 [arXiv:1007.1150 [hep-ex]].
- [8] P. Adamson et al. [MINOS Collaboration], “First direct observation of muon antineutrino disappearance,” Phys. Rev. Lett. 107 (2011) 021801 [arXiv:1104.0344 [hep-ex]].
- [9] A. Zichichi, “Results from LVD-OPERA Combined Analysis: A Time-Shift in the OPERA Setup,” available online at http://agenda.infn.it/getFile.py/access?resId=0&materialId=slides&confId=4896.
- [10] E. D. Zimmerman [MiniBooNE Collaboration], “Updated Search for Electron Antineutrino Appearance at MiniBooNE,” arXiv:1111.1375 [hep-ex].
- [11] P. Adamson et al. [MINOS Collaboration], “An improved measurement of muon antineutrino disappearance in MINOS,” arXiv:1202.2772 [hep-ex].
- [12] T. .A. Mueller, D. Lhuillier, M. Fallot, A. Letourneau, S. Cormon, M. Fechner, L. Giot and T. Lasserre et al., “Improved Predictions of Reactor Antineutrino Spectra,” Phys. Rev. C 83 (2011) 054615 [arXiv:1101.2663 [hep-ex]]. P. Huber, “On the determination of anti-neutrino spectra from nuclear reactors,” Phys. Rev. C 84 (2011) 024617 [Erratum-ibid. C 85 (2012) 029901] [arXiv:1106.0687 [hep-ph]].
- [13] G. Mention, M. Fechner, T. .Lasserre, T. .A. Mueller, D. Lhuillier, M. Cribier and A. Letourneau, “The Reactor Antineutrino Anomaly,” Phys. Rev. D 83 (2011) 073006 [arXiv:1101.2755 [hep-ex]].
- [14] P. C. de Holanda and A. Y. .Smirnov, “Homestake result, sterile neutrinos and low-energy solar neutrino experiments,” Phys. Rev. D 69 (2004) 113002 [hep-ph/0307266]. P. C. de Holanda and A. Y. .Smirnov, “Solar neutrino spectrum, sterile neutrinos and additional radiation in the Universe,” Phys. Rev. D 83 (2011) 113011 [arXiv:1012.5627 [hep-ph]].
- [15] R. Abbasi et al. [IceCube Collaboration], “Measurement of the atmospheric neutrino energy spectrum from 100 GeV to 400 TeV with IceCube,” Phys. Rev. D 83 (2011) 012001 [arXiv:1010.3980 [astro-ph.HE]].
- [16] C. Giunti and M. Laveder, “Implications of 3+1 Short-Baseline Neutrino Oscillations,” Phys. Lett. B 706 (2011) 200 [arXiv:1111.1069 [hep-ph]].
- [17] J. Hamann, S. Hannestad, G. G. Raffelt and Y. Y. Y. Wong, “Sterile neutrinos with eV masses in cosmology: How disfavoured exactly?,” JCAP 1109 (2011) 034 [arXiv:1108.4136 [astro-ph.CO]].
- [18] F. P. An et al. [DAYA-BAY Collaboration], “Observation of electron-antineutrino disappearance at Daya Bay,” Phys. Rev. Lett. 108 (2012) 171803 [arXiv:1203.1669 [hep-ex]].
- [19] L. A. Mikaelyan and V. V. Sinev, “Neutrino oscillations at reactors: What next?,” Phys. Atom. Nucl. 63 (2000) 1002 [Yad. Fiz. 63N6 (2000) 1077] [hep-ex/9908047].
- [20] D. Dwyer, “Daya Bay Results,” presented at Neutrino 2012 in Kyoto. Available at http://neu2012.kek.jp/neu2012/programme.html.
- [21] Y. Abe et al. [DOUBLE-Chooz Collaboration], “Indication for the disappearance of reactor electron antineutrinos in the Double Chooz experiment,” Phys. Rev. Lett. 108 (2012) 131801 [arXiv:1112.6353 [hep-ex]].
- [22] J. K. Ahn et al. [RENO Collaboration], “Observation of Reactor Electron Antineutrino Disappearance in the RENO Experiment,” Phys. Rev. Lett. 108 (2012) 191802 [arXiv:1204.0626 [hep-ex]].
- [23] Y. Kim, “Observation of reactor neutrino disappearance at RENO,” presented at TURN 2012 under Gran Sasso. Available at http://agenda.infn.it/contributionListDisplay.py?confId=4722.
- [24] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo and A. M. Rotunno, “Evidence of from global neutrino data analysis,” Phys. Rev. D 84 (2011) 053007 [arXiv:1106.6028 [hep-ph]].
- [25] T. Schwetz, M. Tortola and J. W. F. Valle, “Where we are on : addendum to ’Global neutrino data and recent reactor fluxes: status of three-flavour oscillation parameters’,” New J. Phys. 13 (2011) 109401 [arXiv:1108.1376 [hep-ph]].
- [26] F. Boehm, J. Busenitz, B. Cook, G. Gratta, H. Henrikson, J. Kornis, D. Lawrence and K. B. Lee et al., “Final results from the Palo Verde neutrino oscillation experiment,” Phys. Rev. D 64 (2001) 112001 [hep-ex/0107009].
- [27] M. Apollonio et al. [Chooz Collaboration], “Search for neutrino oscillations on a long baseline at the Chooz nuclear power station,” Eur. Phys. J. C 27 (2003) 331 [hep-ex/0301017].
- [28] Misaki Ishitsuka, “Double Chooz Results,” presented at Neutrino 2012 in Kyoto. Available at http://neu2012.kek.jp/neu2012/programme.html.
- [29] X. -J. Bi, P. -H. Gu, X. -l. Wang and X. -M. Zhang, “Thermal leptogenesis in a model with mass varying neutrinos,” Phys. Rev. D 69 (2004) 113007 [hep-ph/0311022]. R. Takahashi and M. Tanimoto, “Model of mass varying neutrinos in SUSY,” Phys. Lett. B 633 (2006) 675 [hep-ph/0507142]. R. Takahashi and M. Tanimoto, “Speed of sound in the mass varying neutrinos scenario,” JHEP 0605 (2006) 021 [astro-ph/0601119]. E. Ciuffoli, J. Evslin, J. Liu and X. Zhang, “OPERA and a Neutrino Dark Energy Model,” arXiv:1109.6641 [hep-ph].
- [30] D. B. Kaplan, A. E. Nelson and N. Weiner, “Neutrino oscillations as a probe of dark energy,” Phys. Rev. Lett. 93 (2004) 091801 [hep-ph/0401099].
- [31] M. Tortola, J. W. F. Valle and D. Vanegas, “Global status of neutrino oscillation parameters after recent reactor measurements,” arXiv:1205.4018 [hep-ph].
- [32] G.L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo and A.M. Rotunno, “Global analysis of neutrino masses, mixings and phases: entering the era of leptonic CP violation searches,” arXiv:1205.5254 [hep-ph].
- [33] F. P. An et al. [Daya Bay Collaboration], “A side-by-side comparison of Daya Bay antineutrino detectors,” arXiv:1202.6181 [physics.ins-det].
- [34] T. Lasserre, G. Mention, M. Cribier, A. Collin, V. Durand, V. Fischer, J. Gaffiot and D. Lhuillier et al., “Comment on Phys. Rev. Lett. 108, 191802 (2012): ’Observation of Reactor Electron Antineutrino Disappearance in the RENO Experiment’,” arXiv:1205.5626 [hep-ex].
- [35] X. Guo et al. [Daya-Bay Collaboration], “A Precision measurement of the neutrino mixing angle theta(13) using reactor antineutrinos at Daya-Bay,” hep-ex/0701029.
- [36] Y. Declais, H. de Kerret, B. Lefievre, M. Obolensky, A. Etenko, Y. .Kozlov, I. Machulin and V. Martemyanov et al., “Study of reactor anti-neutrino interaction with proton at Bugey nuclear power plant,” Phys. Lett. B 338 (1994) 383.
- [37] A. Melchiorri, O. Mena, S. Palomares-Ruiz, S. Pascoli, A. Slosar and M. Sorel, “Sterile Neutrinos in Light of Recent Cosmological and Oscillation Data: A Multi-Flavor Scheme Approach,” JCAP 0901 (2009) 036 [arXiv:0810.5133 [hep-ph]].
- [38] S. Hannestad, I. Tamborra and T. Tram, “Thermalisation of light sterile neutrinos in the early universe,” arXiv:1204.5861 [astro-ph.CO].
- [39] D. J. Eisenstein et al. [SDSS Collaboration], “Detection of the baryon acoustic peak in the large-scale correlation function of SDSS luminous red galaxies,” Astrophys. J. 633 (2005) 560 [astro-ph/0501171].
- [40] W. J. Percival, S. Cole, D. J. Eisenstein, R. C. Nichol, J. A. Peacock, A. C. Pope and A. S. Szalay, “Measuring the Baryon Acoustic Oscillation scale using the SDSS and 2dFGRS,” Mon. Not. Roy. Astron. Soc. 381 (2007) 1053 [arXiv:0705.3323 [astro-ph]].
- [41] A. Lewis and S. Bridle, “Cosmological parameters from CMB and other data: A Monte Carlo approach,” Phys. Rev. D 66 (2002) 103511 [astro-ph/0205436].
- [42] E. Komatsu et al. [WMAP Collaboration], “Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation,” Astrophys. J. Suppl. 192 (2011) 18 [arXiv:1001.4538 [astro-ph.CO]].
- [43] N. Suzuki et al., “The Hubble Space Telescope Cluster Supernova Survey: V. Improving the Dark Energy Constraints Above z¿1 and Building an Early-Type-Hosted Supernova Sample,” arXiv:1105.3470 [astro-ph.CO].
- [44] A. G. Riess, L. Macri, S. Casertano, H. Lampeitl, H. C. Ferguson, A. V. Filippenko, S. W. Jha and W. Li et al., “A 3% Solution: Determination of the Hubble Constant with the Hubble Space Telescope and Wide Field Camera 3,” Astrophys. J. 730 (2011) 119 [Erratum-ibid. 732 (2011) 129] [arXiv:1103.2976 [astro-ph.CO]].
- [45] J. -Q. Xia, G. -B. Zhao, B. Feng, H. Li and X. Zhang, “Observing dark energy dynamics with supernova, microwave background and galaxy clustering,” Phys. Rev. D 73 (2006) 063521 [astro-ph/0511625]. G. -B. Zhao, J. -Q. Xia, M. Li, B. Feng and X. Zhang, “Perturbations of the quintom models of dark energy and the effects on observations,” Phys. Rev. D 72 (2005) 123515 [astro-ph/0507482].
- [46] E. Giusarma, M. Archidiacono, R. de Putter, A. Melchiorri and O. Mena, “Sterile neutrino models and nonminimal cosmologies,” Phys. Rev. D 85 (2012) 083522 [arXiv:1112.4661 [astro-ph.CO]].
- [47] N. Suzuki, D. Rubin, C. Lidman, G. Aldering, R. Amanullah, K. Barbary, L. F. Barrientos and J. Botyanszki et al., “The Hubble Space Telescope Cluster Supernova Survey: V. Improving the Dark Energy Constraints Above z=1 and Building an Early-Type-Hosted Supernova Sample,” Astrophys. J. 746 (2012) 85 [arXiv:1105.3470 [astro-ph.CO]].
- [48] B. Feng, M. Li, Y. -S. Piao and X. Zhang, “Oscillating quintom and the recurrent universe,” Phys. Lett. B 634 (2006) 101 [astro-ph/0407432]. B. Feng, X. L. Wang and X. M. Zhang, “Dark energy constraints from the cosmic age and supernova,” Phys. Lett. B 607 (2005) 35 [arXiv:astro-ph/0404224]. X. -F. Zhang, H. Li, Y. -S. Piao and X. -M. Zhang, “Two-field models of dark energy with equation of state across -1,” Mod. Phys. Lett. A 21, 231 (2006) [astro-ph/0501652].
- [49] R. Foot and R. R. Volkas, “Reconciling sterile neutrinos with big bang nucleosynthesis,” Phys. Rev. Lett. 75 (1995) 4350 [hep-ph/9508275].
- [50] H. Motohashi, A. A. Starobinsky and J. ’i. Yokoyama, “Cosmology based on f(R) Gravity admits 1 eV Sterile Neutrinos,” arXiv:1203.6828 [astro-ph.CO].
- [51] R. H. Brandenberger, N. Kaiser, D. N. Schramm and N. Turok, “Galaxy and Structure Formation with Hot Dark Matter and Cosmic Strings,” Phys. Rev. Lett. 59 (1987) 2371. R. H. Brandenberger, A. Mazumdar and M. Yamaguchi, “A Note on the robustness of the neutrino mass bounds from cosmology,” Phys. Rev. D 69 (2004) 081301 [hep-ph/0401239].
- [52] J. Evslin and S. B. Gudnason, “High Q BPS Monopole Bags are Urchins,” arXiv:1111.3891 [hep-th]. J. Evslin and S. B. Gudnason, “Dwarf Galaxy Sized Monopoles as Dark Matter?,” arXiv:1202.0560 [astro-ph.CO].
- [53] G. Karagiorgi, M. H. Shaevitz and J. M. Conrad, “Confronting the short-baseline oscillation anomalies with a single sterile neutrino and non-standard matter effects,” arXiv:1202.1024 [hep-ph].
- [54] G. Karagiorgi, M. H. Shaevitz and J. M. Conrad, “Confronting the short-baseline oscillation anomalies with a single sterile neutrino and non-standard matter effects,” arXiv:1202.1024 [hep-ph].