Phenomenology of MaVaN’s Models in Reactor Neutrino Data

Mass Varying Neutrinos (MaVaN’s) mechanisms were proposed to link the neutrino mass scale with the dark energy density, addressing the coincidence problem. In some scenarios, this mass can present a dependence on the baryonic density felt by neutrinos, creating an effective neutrino mass that depends both on the neutrino and baryonic densities. In this work, we study the phenomenological consequence of MaVaN’s scenarios in which the matter density dependence is induced by Yukawa interactions of a light neutral scalar particle which couples to neutrinos and matter. Under the assumption of one mass scale dominance, we perform an analysis of KamLAND neutrino data which depends on 4 parameters: the two standard oscillation parameters, Δm2 0,21


Introduction
In cosmology, dark energy is a hypothetical form of energy which permeates all space and leads to the accelerated expansion of the universe [1].The current standard model of cosmology has 68% of the total energy of the universe in dark energy [2].Regardless of its nature, dark energy must have a strong negative pressure to explain the accelerated expansion of the universe and was initially motivated by evidence from supernovae [3,4] and data from cosmic microwave background [5,6].
The need for a new mysterious dark energy component can be interpreted as an indication for physics beyond the standard model.In the present time, the density of dark matter and dark energy is similar,  CDM / Λ ∼ 1/3; however this ratio evolves as a function of the scale factor with 1/ 3 .The fact that today these two quantities are of the same order, although they vary dramatically in the history of the universe, became known as the cosmic coincidence problem.The model that intends to explain the dark energy through a cosmological constant needs very specific initial conditions to achieve the correct values for the dark energy density today, giving us clues that this coincidence can be better explained if dark energy has tracked another components of matter since recombination.
Relying on the similarity of scales, (2 × 10 −3 eV) 4 for dark energy and (10 −2 eV) 2 for the neutrino mass split, [7,8] propose to explain the nature of dark energy postulating the neutrino mass as a dynamic quantity, depending on the value of a scalar field .Considering a very slow variation in the parameters involved, one can take the system as always in equilibrium, with the potential following the minima.
The potential  is then considered very flat and its magnitude dependent on the density of cosmological neutrinos.As a result, these mass varying neutrinos (MaVaN's) become heavier with the decrease of its density.The total energy of the fluid (contained in the neutrinos and in the field ), identified as dark energy, may vary smoothly while the density of neutrinos decreases.This would explain not only the origin of dark energy but can also significantly modify the limits of cosmological neutrino mass [7], the relationship between the mass of Advances in High Energy Physics the neutrino and leptogenesis [8], and the expected change in the split of flavors for neutrinos and cosmic background from distant astrophysical sources [9].
It can be shown that interactions of subgravitational force can occur naturally between ordinary matter and the field , which can cause the value of  to vary from its value in vacuum [10].This leads to neutrinos with masses dependent on the density of the medium and new effects in the flavor oscillation, motivating the study of phenomenological effects of this type of model on the Sun [11][12][13], in atmospheric neutrinos [14,15], supernovae [16,17], reactor neutrinos [18,19], and neutrinos propagating on Earth [20].
In this work, we investigate the phenomenology of MaVaN's effects in neutrino oscillation, focusing in reactor neutrino data.By applying a parametrization of the effect already used in another context and with a detailed Earth's crust description, we intend to analyze experimental data with this hypothetical matter-neutrino interaction.
In Section 2, we describe the MaVaN formalism and the particular environment approach used and its effects in neutrino oscillations.In Section 3, we present the results and examine how these modifications can affect the reactor neutrino data.In Section 4, we summarize our conclusions.

MaVaN's Mechanism and Parametrization.
Previous works found limits for the product of the effective neutrinoscalar and matter-scalar Yukawa coupling described in [21].With the assumption of environment homogeneity, the new physics evoked always plays the role of a subleading effect compared to the standard oscillation scenario.
The aim of this work is to find at least one combination of parameters for the new physics that could lead to an acceptable solution to the neutrino oscillation data where such new physics is more than a sub-leading effect.This implies a nonhomogeneity effect of the Earth's crust in the neutrino evolution.
We consider an effective low energy model containing the standard model particles plus a light scalar () of mass   which couples very weakly both to neutrinos (]  ) and the matter fields  = , , .
The Lagrangian takes the form where  0  are the vacuum mass that the neutrinos would have in the presence of the cosmic neutrino background and  ]  and   are, respectively, the effective neutrino-scalar and matter-scalar couplings.We have written a Lagrangian for Dirac neutrinos but equivalently it could be written for Majorana neutrinos.
It has been argued in [22] that, generically, these models contain a catastrophic instability which occurs when neutrinos become nonrelativistic.Since this transition occurs later for hierarchical neutrinos, in what follows we assume the vacuum neutrino masses to be hierarchical: For solar neutrinos of hierarchical masses, the dominant contribution to the neutrino mass is due to the matter background density.Correspondingly, we neglect the contribution to the neutrino mass from the background neutrino density, and we concentrate on the matter density dependence.
With the additional freedom that the new matter density dependence provides, there is no reason to believe that the three neutrino oscillation dynamics factorize into the dynamics of two neutrino subsystems.However, we will assume that this is still the case and study their effect on solar and KamLAND oscillations under the hypothesis of one mass-scale dominance.Under this assumption, we parametrize the evolution equation as where   12 is the 2×2 mixing matrix in vacuum parameterized by the angle  12 and   () are the environment dependence contributions to the neutrino masses; meanwhile  CC () = √ 2    () is the MSW potential proportional to the electron number density   () in the medium.The environment effect is introduced as a dependence of the mass terms with the baryonic matter density with the following parametrization: which has been chosen to reproduce two features that we intend our mass matrix to present a linear growth of mass with baryonic density for small values of this density as suggested in [21] and a saturation of the environmental dependence of neutrino masses for large values of the baryonic density.It is important to stress that such behaviour is an ad hoc choice of our parametrization in a phenomenological approach to the problem, and it is not derived from details of the interaction between neutrinos, baryonic matter, and the scalar field.

A Study of the Environment.
The KamLAND collaboration uses the constant density crust approximation of 2.7 gr/cm 3 in the neutrino evolution analysis.In order to see a more detailed effect than that found in the general literature, we use a specific model of the Earth's crust.The model used was the CRUST 2.0 [23] which is specified in a grid of 2 × 2 degrees on the terrestrial surface.The 2×2 degree model is composed of 360 key 1d-profiles where one of these profiles is assigned to each 2 × 2 degree cell.Each individual profile is a one-dimensional description of 7 layers: ice, water, soft sediments, hard sediments, upper crust, middle crust, and lower crust.For each of these layers, thickness and density are locally given.Any path determined to be above the ice is considered outside the crust and therefore under the influence of typical atmosphere density.
The KamLAND experiment uses as a source of antineutrinos several different reactors.With the geographic location of each reactor and considering a linear path to the detector, this path will describe a density map specific for each source.The neutrino cross specific values of depth as shown in Figure 1, and hence different densities.Examples are shown in Figure 2. To the best knowledge of the authors of this work this model has not been used before in this context in the reactor neutrinos literature.

Effects in Reactor Neutrinos
Oscillations.As we can see in the specific density maps, it is clear that the description leads to a completely nonadiabatic evolution of the neutrino.So the effective angle variation relative to the local density change causes the nondiagonal terms in the evolution equation of the mass eigenstates to become relevant.So one cannot use here the analytic survival probability.
We use the analysis of the probability amplitudes, which allows local calculation for each region of constant density.Using what is called slab approximation in [24], the effective massive neutrinos propagate as plane waves in regions of constant density.
The transition amplitude will be defined in each slab of constant density so The notation   indicates that all the matter-dependent quantities in the  definition must be evaluated with the matter density in the th slab.The final amplitude of a constant density slab is used as initial value at the beginning of the next slab.
In the case of choice, we consider a pure electronic initial neutrino so ] 0 | =0 = ( 1 0 )| 0 .In the two families dynamics we have as the survival probability where  21 is the final transition amplitude between the point of creation and detection.We could have a smooth description of the density profiles based on the stepped profiles generated by the model.However, this would result in a calculation much more expensive computationally since integration would be numerical outside the approximation of (6).

𝜒 2
Analysis.The results of the simulation, that is, the expected number of events for each energy range are analyzed by the method of maximum likelihood.Particularly, the Poisson statistics where the  2 to be minimized is defined by where  is a parameter that can vary, leading to a theoretical flux free normalization.The sum is made on the  energy intervals,  teo  is the event number expected theoretically, and  obs  is the number of observed events given by collaboration, both in the energy interval .

Results and Discussion
We developed a code to simulate the KamLAND neutrino events and, based on the presented parametrization and phenomenology of MaVaN's models, we tested the new physics.We implemented the model with normal hierarchy and a lighter neutrino with zero mass, to avoid the instability already mentioned: that is, leading to an effect only on the diagonal mass eigenstates Therefore, the probability amplitude for each interval of constant density can be written as The effective KamLAND mass split is given by where (  ) is specific density of the interval.
As we can see in (10), we included an effect that does not generate oscillations in the mass eigenstates adding two free parameters that does not interfere with the kinetic term evolution of the neutrino.
To cover then the parameter space, we let the values run between 0.2 < tan 2  12 < 0.8, 6.0 × 10 −5 eV 2 < Δ with Figures 4-7 show the simulation result and its relation with the values previously obtained from the standard model.The description of the data is not improved in any specific bin.There is an improved description of the data in the general distribution of energy, as we can see in Figure 5. Figures 6 and  7 show, respectively, Δ 2 as a function of Δ 2 and Δ 2 as a function of tan  with the other parameters fixed in the best fit point given by (16).
We can see that there is a statistically significant improvement of the fit with the inclusion of terms MaVaN; however, we must be careful when interpreting this result, since we choose a specific parametrization on the behaviour of the MaVaN effects that is expected to at least reproduce the quality of the fit from the standard oscillation mechanism.
The value of  2 is comparable to  0 2 for certain values of density and influences directly the effective mass of the neutrino.It is not a second order effect.Indeed  2 →  0 2 for  ∼ 0.4 g/cm 3 so the difference in effective mass in matter Δ M2 21 (  ) varies widely for neutrinos crossing density areas with this feature.This kind of resonance creates, in this case, a considerably better description of the data, as can be seen in Figure 3.
Such effect is observed only in a small island in parameter space, as seen in Figure 8.When any one of two parameters tends to zero, the effect MaVaN also vanishes and returns to the standard model, which is not excluded by our results within 99% confidence level.A degeneracy occurs because the new parameters tend to compensate each other.We cannot write specific limits for the new physics does to this degeneracy.
In order to quantify the features of this effect inside such island on MaVaN's parameters, we display the global dependence of  2 on  2 ( 02 ) in Figure 9 (Figure 10) after the marginalization over Δ taking into account, for these specific limits, only the island of best fit.The standard model would be included in the intervals of 99% confidence level if we take into account the whole parameter space.
The results and analysis presented here must be viewed as a specific choice of a particular model.We have three sources of arbitrariness in the model, the choice of  1 =  1 = 0 as  3 = 0 and the hypothesis of the hyperbolic tangent   8: Exclusion curves in the parameter space MaVaN to 95%, 99%, and 99.73% confidence levels.We can see that the effect appears only in an island of values and does not exclude the standard model, which becomes dominant when the parameters tend to zero.parametrization.However, we expect that similar effects would be present if we change any of these assumptions.

Conclusions
We studied the phenomenological effects of a specific dependence of the effective neutrino mass with the local density of the propagation medium, considering a model that included the standard model plus a light scalar () that couples weakly to all constituents of the matter.We used a specific parameterization for the new physics and a detailed description of the Earth's crust density, the CRUST 2.0, which was not previously used in the literature in this context.
Assuming that neutrinos masses follow the hierarchy 0 =  0 1 <  0 2 <  0 3 , we analyzed data from KamLAND (20 points −0.01 × 10 −2 eV.This point corresponds to a decrease of Δ 2 /Δd.o.f.= 6.22/2 compared to the minimum where only the standard matter effect was considered.However, we have seen that this strong effect, although statistically significant, appears in an island of specific values, and the results remain consistent with the standard model within 99% confidence level.
We observed that the new physics is more then a subleading effect and clearly dependent on the specific density description used.Such effect appears due to a strong variation of neutrino parameters at lower densities than the value usually used as constant for Earth crust in the literature, impliyng a nonhomogeneity effect of the Earth's crust in the neutrino evolution.
To write specific limits including the null value for the new physics parameters is not possible due to a parameter degeneracy.The standard model is reproduced when any of the parameters that describes the MaVaN effects goes to zero.We believe that a joint analysis with solar neutrinos may clarify this point.These limits will be tested in the future for an analysis with solar neutrinos.

Figure 1 :Figure 2 :
Figure 1: Simplified representation of the creation-detection neutrino path, different cases generate different density maps.

Figure 3 : 2 ) 2 Figure 4 :
Figure 3: Correlation between the effective value of mass after the inclusion of the effect MaVaN and the same term with only the parameters of the standard model.
2  21, tan 2  12 , and  02 ( 2 ).From the results we found the following bounds for 90% and 99.73% CL:1.18 <  2 < 1.42, 1.89 × 10 −2 eV <  02 < 1.92 × 10 −2 eV,Figure 5: Ratio curve of the survival probability of survival with the case of no oscillation.We compared the experimental data (red) with the description of the standard model (black) and the obtained for the MaVaN model (green).Figure 6: Δ 2 as a function of Δ 2 , when the other parameters are fixed in the best fit point, given by (16).Dotted lines indicate, from bottom to top, 90% and 99.73% confidence levels.
Figure7: Δ 2 as a function of tan , when the other parameters are fixed in the best fit point, given by (16).Dotted lines indicate, from bottom to top, 90% and 99.73% confidence levels.