# Constraints on neutrino decay lifetime using long-baseline charged and neutral current data

###### Abstract

We investigate the status of a scenario involving oscillations and decay for charged and neutral current data from the MINOS and T2K experiments. We first present an analysis of charged current neutrino and anti-neutrino data from MINOS in the framework of oscillation with decay and obtain a best fit for non-zero decay parameter . The MINOS charged and neutral current data analysis results in the best fit for eV, and zero decay parameter, which corresponds to the limit for standard oscillations. Our combined MINOS and T2K analysis reports a constraint at the 90% confidence level for the neutrino decay lifetime s/eV. This is the best limit based only on accelerator produced neutrinos.

###### keywords:

neutrino, neutrino decay, neutrino oscillation###### Pacs:

14.60.Lm, 14.60.St, 14.60.Pq^{†}

^{†}journal: \biboptions

numbers,square,comma,sort&compress

## 1 Introduction

Recent results from reactor experiments, Double Chooz Abe et al. (2012), Daya Bay An et al. (2014) and Reno Ahn et al. (2012), complete the picture that three active neutrinos oscillate with two known non-zero mass differences ( and ), three mixing angles (, and ) and an unknown CP phase () Pontecorvo (1957); Maki et al. (1962). For a detailed description of neutrino oscillations see Ref. Gonzalez-Garcia and Nir (2003).

In this picture the large statistics of atmospheric neutrinos of the Super-Kamiokande Abe et al. (2013) and IceCube Aartsen et al. (2013) experiments show us that the deficit in the muon events can be understood as the result of oscillation . Other experiments also show a strong signal for disappearance. MINOS, for instance, fixed very precisely the scale of oscillations at the value eV and Adamson et al. (2013). The reactor experiments Abe et al. (2012); An et al. (2014); Ahn et al. (2012) show evidence for electron neutrino disappearance. For instance, the Daya Bay experiment show a signal for oscillations with a scale of eV (the parameter is the properly averaged quantity between and Minakata et al. (2006)) and with amplitude of An et al. (2014). From solar neutrino experiments Gonzalez-Garcia and Nir (2003) and the reactor experiment KamLAND Gando et al. (2013) there is evidence for (anti-)electron neutrino disappearance. These two signals of oscillations can be explained by eV and the large mixing angle Gando et al. (2013), which are associated to what is called the Large Mixing Angle (LMA) solution for the solar neutrino anomaly.

Now that it is established that neutrinos are massive, this also implies that they could decay. The idea of neutrino decay is as old as the idea of neutrino masses and mixing. The decays of the neutrino in the Standard Model and are already too constrained and will not be discussed here (see Ref. Mirizzi et al. (2007)). An interesting possibility is the scenario where the neutrino decays into another neutrino and a scalar (or Majoron): decays, where can be a scalar or pseudoscalar massless boson Zatsepin and Smirnov (1978); Lindner et al. (2001). These non-radiative decays are from two types: (I) invisible decays, where neutrinos decay into non-observable final states Bandyopadhyay et al. (2001); Kelso and Hooper (2010); Lunardini (2010); Ando (2003); Fogli et al. (2004); Barger et al. (1999); Gonzalez-Garcia and Maltoni (2008); Palomares-Ruiz et al. (2005); Adamson et al. (2010); and (II) visible decays, where the final products contain active neutrinos Raghavan et al. (1988); Berezhiani et al. (1992); Das and Pulido (2011); Eguchi et al. (2004); Aharmim et al. (2004); Berezhiani and Smirnov (1989); Lindner et al. (2002); Ando (2004).

We can parametrize the decay by the ratio of the lifetime parameter and the mass for each of the mass eigenstates . The role of invisible neutrino decay was investigated for the solar neutrino anomaly Bandyopadhyay et al. (2001); Joshipura et al. (2002); Beacom and Bell (2002) and showed no evidence for the dominance of the decay scenario. From these we can constrain values of the decay parameter, at 90% C.L., where and are respectively the lifetime and the highest mass eigenstate in a two generation scenario Bandyopadhyay et al. (2001).

In the visible decay scenario, we can search for conversion using a pure source such as the Sun Raghavan et al. (1988). The null results from the solar appearance impose a constrain on the decay parameter: s/eV from the KamLAND experiment Eguchi et al. (2004).

Neutrinos produced in supernovas are interesting to investigate for the presence of decays due to the large distance traveled by the neutrino. From the observation of electron neutrinos from SN1987A Hirata et al. (1988); Bionta et al. (1987) we should have a lower limit in neutrino lifetime. Otherwise we could not see any signal. For larger values of the mixing angle, such as the current LMA solution for the solar neutrino anomaly, no constraint is possible Frieman et al. (1988). Other possibilities for neutrinos coming from a supernova include the neutrino decay catalyzed by a very dense media Berezhiani and Smirnov (1989); Berezhiani et al. (1992) where the matter effects can increase the decay rate of neutrinos. The diffuse supernova neutrinos (neutrinos coming from all past supernova explosions) Lunardini (2010), can provide very robust sensitivity in the range of s/eV Ando (2003); Fogli et al. (2004).

Astrophysical neutrino sources megaparsec away can generate all neutrino flavors. Due to the long distance from the sources we are in the limit , where all dependence on the lifetime parameter fades away. However, if we have a precise determination of the ratio of flavors of these neutrinos we can discriminate the case with and without decay Meloni and Ohlsson (2007); Maltoni and Winter (2008); Baerwald et al. (2012); Dorame et al. (2013).

Concerning the accelerator and atmospheric neutrinos we can test the decay scenario of the third generation of neutrino mass eigenstates investigating how the decay parameter changes the survival probability Barger et al. (1999); Gonzalez-Garcia and Maltoni (2008). The MINOS experiment made a search for the decaying neutrino and constrained the lifetime to s/eV at 90% C.L., using both neutral and charged current events Adamson et al. (2010). The combined analysis of Super-Kamiokande atmospheric neutrinos with K2K and MINOS accelerator neutrinos show a 90% C.L. lower bound value of s/eV Gonzalez-Garcia and Maltoni (2008).

This article is organized as follows, in Section 2 we discuss our neutrino decay scenario. Next, we introduce the analysis developed for the neutrino charged and neutral current data of the MINOS experiment in Section 3. We then present our bounds for the neutrino lifetime based on a charged current analysis (Section 4.1) and on a combined charged and neutral current analysis (Section 4.2) using MINOS data. In Section 4.3 we discuss the sensitivity of the T2K experiment for the neutrino decay scenario and present the constraints on neutrino lifetime for analyses using T2K data only and the combined MINOS and T2K data. Section 5 presents a discussion on the relation of this scenario under Majoron models.

## 2 Decay model for neutrinos

We are going to introduce the neutrino evolution equation in which the neutrino can decay. This is made by putting an imaginary part related to the neutrino lifetime, which is the ratio , in the evolution equation. We are going to assume the decay of the heaviest state, , where both final products are invisible. The two-generation system is considered, which is adequate to describe the muon neutrino and anti-neutrino data from MINOS. The evolution equation is

(1) |

where the state , and E is the neutrino energy. is the usual rotation matrix,

(2) |

where and . The same evolution equation applies for anti-neutrinos as well.

From Eq. (1) we obtain the muon neutrino survival probability as

(3) | ||||

where is the distance traveled by the neutrinos. We can notice that a non-zero decaying parameter changes only the amplitudes: the constant amplitude (first term of the equation above), and the oscillation amplitude (second term). Both amplitudes are damped but the oscillation phase does not change. In the two- standard oscillation probability formula we have the symmetry , but in Eq. (3) the symmetry is broken, and then we should scan the parameter space of the variable in the range . This broken symmetry will appear in our plots later.

The limiting case where the oscillations are induced only by decay, , can also be tested. In this case the probability assumes the simple form

(4) |

where now we have two free parameters, the mixing amplitude and the decay parameter . Even at this limit we also have an asymmetry between and .

In the standard neutrino oscillation scenario, the sum of the probabilities over the active states is equal to unity. Then the spectrum of neutral current (NC) events is not effected by active oscillation, which means that the expected number of NC events is the same with or without oscillations. But in the extended scenario involving sterile neutrinos the sum of probabilities, , where , obviously does not sum up to 1. This is discussed in the general context of non-unitary neutrino evolution in Ref. Berryman et al. (2014).

We can compute the conversion probability for the two- oscillations with decay scenario,

(5) | ||||

which implies that

(6) |

Thus, from the Eq. (6) we observe that there will be an effect on the neutral current interaction events under the oscillations with decay model.

## 3 Analysis of MINOS Charged and Neutral Current Data

We have performed a combined analysis using the published data of charged and neutral current MINOS analyses. MINOS is a long-baseline neutrino experiment Michael et al. (2008) using two detectors and exposed to a neutrino beam produced at Fermilab. The NuMI beam line is a two-horn-focused neutrino beam that can be configured to produce muon neutrinos or anti-neutrinos. The Near Detector is located at Fermilab, around 1 km from the NuMI target and the Far Detector is 735 km far from the target.

The first data set used in our analysis comprises the charged current (CC) contained-vertex neutrino disappearance data Adamson et al. (2013) using the enhanced beam with exposure of protons on target (POT), with 23 points of non-equally divided bins of energy up to 14 GeV; the second is the CC contained-vertex anti-neutrino disappearance data Adamson et al. (2013) from the enhanced beam with POT, using 12 points for energies up to 14 GeV as well; and the third is the neutral current data Adamson et al. (2011a), based on POT. The spectrum of NC events is described as a function of a reconstruted energy, . We then use the information from a previous MINOS analysis Adamson et al. (2010) to separate the NC data into two bins: (a) events with GeV and (b) events with GeV, which have median neutrino energies of 3.1 and 7.9 GeV, respectively.

For the calculation we use the following function

(7) |

in bins of energy, where is the prediction for the theoretical model () that we are using (e.g. oscillation plus decay scenario), based on the no-oscillation events (), and including the background contribution () adjustable by a parameter . is the data from one of the MINOS analyses and is the total error. The background events of the NC data set come from misidentified CC events. The numbers , and were read off from References Adamson et al. (2013) and Adamson et al. (2011a).

The total error used for both neutrino and anti-neutrino CC data sets is given by

(8) |

where is the total error of the data, and are the statistical and systematic error of the prediction, respectively. We have , where is the lower (upper) error bar from the data; and Adamson et al. (2008). For the NC data set we do not have the total error of the data events. We then calculate it summing in quadrature the statistical and systematic errors. The former is the square root of the number of data events, and the latter is an estimate based on the systematic error of the extracted expectation of events.

We scan the parameter region in the variables , and and for the parameter of our function . The is then marginalized over the nuisance parameter to obtain the effective as a function of the other parameters

To have a guess about the range of the decay parameter we can probe, we will assume the argument of the exponential term in Eq. (3) is of the order of unity, where we can get the most sensitivity. Using the MINOS distance and the energy range of GeV, we get s/eV.

Model | /d.o.f. | ||||||

No oscillation | our MINOS CC | 324.70/35 | |||||

our T2K CC | 185.92/27 | ||||||

Standard oscillation | MINOS Ref. Adamson et al. (2013) | NA | 2.41 | 0.95 | |||

our MINOS CC | 19.80/31 | 2.38 | 0.92 | ||||

T2K Ref. Abe et al. (2014) | NA | 2.51 | 1.00 | 0.51 | |||

our T2K CC | 12.66/25 | 2.44 | 0.97 | 0.58 | |||

Pure decay | MINOS Ref. Adamson et al. (2010) | 76.4/40 | 0.96 | 0.60 | |||

our MINOS CC | 51.18/31 | 0.00 | 1.00 | ||||

our MINOS CC+NC | 69.73/32 | 0.38 | 0.89 | ||||

our T2K CC | 23.11/25 | 0.86 | 0.69 | ||||

our MINOS+T2K | 117.28/59 | 0.08 | 0.98 | ||||

Osc. with decay | MINOS Ref. Adamson et al. (2010) | 47.5/39 | NA | 1.00 | 0.50 | ||

our MINOS CC | 19.28/30 | 2.30 | 0.93 | 0.63 | |||

our MINOS CC+NC | 19.88/31 | 2.34 | 0.96 | 0.60 | |||

our T2K CC | 7.33/24 | 2.46 | 1.00 | 0.49 | |||

and | |||||||

our MINOS+T2K | 32.48/58 | 2.34 | 0.97 | 0.42 |

## 4 Results

### 4.1 MINOS Charged Current Analysis

To test the correctness of our analysis, we first consider the standard oscillation scenario, corresponding to the limit in Eq. (3). The results shown here are from neutrino and anti-neutrino charged current data obtained from Ref. Adamson et al. (2013). However, we also used our calculation on the neutrino charged current data from Ref. Adamson et al. (2011b). For both data sets we obtain good concordance with the MINOS allowed region and best fit.

In Table 1 we show the consistency between the best fit parameters obtained by our standard oscillation analysis (accelerator data only) and MINOS Adamson et al. (2013) (which includes accelerator and atmospheric data). We obtain that the best fit of the oscillation parameters is for non-maximal mixing angle, and for eV. The of the best fit point is 19.80 for 31 degrees of freedom, which corresponds to a G.O.F. (Goodness of Fit) of 94.0%. The for the no-oscillation hypothesis compared to the standard oscillation scenario,

is equal to 304.9, which means that we can exclude the no-oscillation hypothesis with remarkable precision.

The allowed regions in the ( – ) plane for neutrinos and anti-neutrinos under the standard oscillations model is shown in Fig. 1. The anti-neutrino data is compatible with the neutrino data but does not contribute significantly to improve the region of parameters in the combined CC analysis due to its small statistics compared to the neutrino one. Our results show that we are able to reproduce the allowed regions of square mass difference and mixing angle of the MINOS analyses.

Having confirmed the correctness of the analysis, we use it to test the oscillation with decay scenario. The procedure is the same as in the standard oscillation case, but now we have three free parameters, , and . For the combined (neutrino and anti-neutrino) CC analysis under the oscillation with decay framework we find a finite value for the best fit of (see Table 1),

(9) |

The effect of a finite value of results in a slightly lower value of than the one obtained for standard oscillation. The obtained is 19.28 for 30 degrees of freedom, which corresponds to a G.O.F. of 93.4%. The pure decay hypothesis is also considered in our CC analysis resulting in a of 51.18 for 31 degrees of freedom, excluding this hypothesis at the 5.6 standard deviation level.

The projections of for each parameter are obtained minimizing the function accordingly to the parameters. We then use the function , where is the global minimum value for the oscillation with decay scenario. Next, we present the projections of as a function of , and .

The one-dimensional projection for the parameter, shown in Fig. 2 (dotted curve), allows us to get a lower bound on or equivalently an upper bound on . The allowed values for the neutrino decay lifetime is

(10) |

at the 90% C.L. (Table 1). We can compare the value we obtain for the oscillation with decay model to the one from MINOS Adamson et al. (2010), which used both charged and neutral current data in their analysis, but with lower statistics than the extracted data used here. We see that our combined CC result does not improve the MINOS lower limit which is s/eV.

An interesting behavior appears in Fig. 3 (left), where we compare the allowed values for with and without decay. When we include a finite value for the parameter, smaller values of are allowed for a certain confidence level. Then, under this scenario, we can explain the muon disappearance signal seen by MINOS as due partially to oscillation and partially to decay, since the limit in Eq. (3) corresponds to the pure decay model.

The broken symmetry, mentioned after Eq. (3), is manifested in Fig. 3 (right), where the curve is broader for larger values of . This effect can be understood investigating the coefficient on the first term of Eq. (3). For the decay term is more relevant, and for the opposite case the decay contribution is suppressed.

We also present two-dimensional projections of the allowed three-dimensional region after normalization with respect to the undisplayed parameter. Figure 4 shows the 90% C.L. allowed region in the ( – ) and ( – ) planes for the oscillation with decay scenario (dotted curves) and our best fit point (circles) for the MINOS CC data only. We can observe that for smaller values of , when the effects of the decay are larger, the contours allow smaller values of and higher values of (the same behaviour shown in Fig. 3).

The 90% C.L. allowed regions for the oscillation parameters, and , is presented in Fig. 5 both with and without decay. In agreement with the information shown in the discussion of Figs. 3 and 4 we can see that the decay allows a region of smaller values of and larger values of than the standard oscillation model.

The best fit points are also shown in Fig. 5 for both models (with and without decay). The standard oscillation analysis obviously shows two possible values of due to the symmetry . Since the oscillation with decay model does not manifest such a symmetry we find that the best fit point for this scenario is in the octant, with a value of .

### 4.2 MINOS Combined Charged and Neutral Current Analysis

In this section we present the results obtained for the MINOS combined charged and neutral current data. We expect an effect on the neutrino lifetime due to the inclusion of the neutral current data in the analysis as shown in Eq. 6. First, we calculate the ratio

(11) |

for each bin of reconstructed energy extracted from the NC data, where , , and are the number of data, background and expected events, respectively. The ratios we obtain are in good agreement with the ones from MINOS Adamson et al. (2011a).

Under the oscillation with decay model the best fit point for the combined CC and NC analysis can be found in Table 1. The and values are consistent with the ones for standard oscillations. We obtain a neutrino lifetime , which corresponds to the case of no decay. If we compare this result to the best fit of the CC analysis, the decay effect becomes less relevant due to the inclusion of NC data, implying that the best fit value of increases (from eV to eV).

The of the best fit for the combined CC and NC analysis is 19.88 for 31 degrees of freedom resulting in a G.O.F. of 93.8%. The pure decay model is also used in our analysis, being excluded at the 7.1 standard deviation level. Our best fit value of for the pure decay scenario is consistent to the one from MINOS Adamson et al. (2010).

The one-dimensional projection for the parameter in the combined CC and NC analysis is shown in Fig. 2 (solid curve). We find a lower limit of

(12) |

at the 90% C.L., which improves the MINOS Adamson et al. (2010) limit. This constraint is the best so far for long-baseline experiments.

The two-dimensional projections in Figs. 4 and 5 show the effect of including the NC data into our analysis. From Eq. (6) we know that the effect of decay for the NC sample is stronger for large values of than for small ones. And since the combined CC and NC analysis does not suggest evidence for decay (), this explains the shrink of the large values of in the allowed region of Fig. 4 (bottom). The lower values of in the allowed region at the top plot of Fig. 4 was cut due to its correlation with the larger values of .

In Fig. 5 we can see that the inclusion of NC data into our analysis decreases the asymmetry between and which is present in the CC only analysis. For the combined CC and NC analysis, larger values of are also excluded by the same reason pointed out in the discussion of the Fig. 4.

The extracted MINOS spectrum data for CC neutrinos and anti-neutrinos are shown in Fig. 6. In both plots we show the curves for the no-oscillation hypothesis, the best fit of standard oscillations, and the background from NC events. In addition, we also show the curves for the oscillation with decay scenario at the best fit for and and using the 90% C.L. value of the parameter for the CC data only ( s/eV) and for the combined CC with NC data ( s/eV). The reason for using the 90% C.L. values of is that the curves for the best fit parameters show no significant difference if compared to the standard oscillations curve. We can observe that the curve using the higher value of , which corresponds to smaller effect of decay, is closer to the one for standard oscillations. We also observe that the inclusion of the decay scenario into the oscillation model worsens the fit for both neutrino and anti-neutrino analyses.

### 4.3 Combined MINOS and T2K Analysis

We investigate the impact on our analysis of the charged current data from the T2K experiment Abe et al. (2014). T2K is a long-baseline neutrino experiment which uses two detectors exposed to a neutrino beam produced at J-PARC ring. The Near Detector is located around 280 m downstream of the neutrino production target and the off-axis Far Detector is located 295 km far from the target. The data we used here has an integrated POT of from the disappearance analysis.

Figure 7 shows the survival probability as a function of the neutrino energy using T2K data for some values of neutrino lifetime under the oscillation with decay scenario. Considering the 90% C.L. limit for from our MINOS CC and NC analysis (dashed curve) we observe an effect on the probability in the first two T2K data points when comparing to the standard oscillation probability (solid curve). The region at the oscillation minimum, where the errors are smaller than in the other regions, shows very small difference for the scenarios with and without decay. Since T2K data is still dominated by statistical errors their sensitivity for the oscillation with decay scenario would certainly benefit with increasing statistics.

Nevertheless, we use the charged current data from the T2K experiment and perform a similar analysis as made for the MINOS experiment, shown in Eq. (7). The number of data () and expected no-oscillation () events were read off respectively from References Abe et al. (2014) and Giganti et al. (2014). Using the T2K data only we have found that the no-oscillation scenario is highly disfavored as shown in Table 1. Under the standard oscillation scenario, we could reproduce the 90% C.L. allowed region and obtain best fit values consistent to the T2K ones (Table 1). The scenario of pure decay for the T2K data is disfavoured compared to standard oscillations, but still being allowed with a G.O.F. of 57.1%. In the more general scenario of decay and oscilations we have found that T2K data prefers non-zero values of decay parameter, s/eV with eV and . This data favours the oscillation with decay scenario with a /d.o.f. = 7.33/24 over the standard oscillation scenario, /d.o.f. = 12.66/25, mainly due to the energy bin centered at 0.35 GeV (second data point in Fig. 7).

We also test the combined MINOS and T2K data analysis under the pure decay scenario, resulting in its exclusion at the 10.5 standard deviation level. The analysis of the combined data considering the oscillation with decay model results in a finite best fit value for the neutrino lifetime, s/eV, with the oscillation parameters eV and (Table 1).

Figure 2 shows the projections of as a function of for (I) the T2K only and (II) the combined MINOS and T2K analyses under the oscillation with decay scenario. We already observed in that figure that the solid curve, from the MINOS combined charged and neutral current analysis, shows a non-finite best fit value for , as described previously. The dash-dotted curve which is from the T2K data only analysis shows a 90% C.L. range of the neutrino lifetime, s/eV. This result is consistent to the closed allowed regions at 90% C.L. of versus and versus for the oscillation with decay scenario (dash-dotted curve) shown in Fig. 4.

In Fig. 2 we also see the curve (dashed) from the combined MINOS and T2K analysis, that results in a limit of

(13) |

at the 90% C.L., which does not improve the contraint obtained by our MINOS analysis. The results due to the inclusion of T2K in our analysis are in agreement with what we expected from the discussion of Fig. 7.

## 5 Majoron models for neutrino decay

The phenomenological model of neutrino decay that we defined by Eq. (1) can be accommodated in well motivated Majoron models. In general grounds the effective Lagrangian for neutrino decay can be written as

(14) |

where g and h are the Majoron-neutrino couplings in mass basis with scalar (pseudo-scalar) massless boson. We can compute the neutrino lifetime for the decay , where is a neutrino eigenstate. Assuming that the third mass eigenstate decay into a lightest state, we have the lifetime in the laboratory system, for , given by

(15) |

We can use the Eq. (15) to compute the decay parameter as

(16) |

If we use the constraint obtained in our CC (CC+NC) analysis, at 90% C.L., into the Eq. (16) we can get an upper limit for Majoron-neutrino coupling constant of the order of

(17) |

The effective Lagrangian for neutrino decay shown in Eq. (14) can be embedded in different extensions of Standard Model. The general trend is to have the inclusion of new scalar particles in different representations with non-universal couplings between the different families and also the addition of new sterile neutrino states. For instance, Ref. Dorame et al. (2013) presents a model with an symmetry that included as new fields one extra singlet scalar boson and three right-handed neutrinos. This model has the assignments that are family-dependent and therefore the Majoron-neutrino couplings and the mass basis are not proportional to each other, making possible the neutrino decay Dorame et al. (2013).

Another example is the model with an symmetry that can have a see-saw scale at TeV energies, when we include two scalars with a singlet and doublet scalar added to this enlarged gauge symmetry Cogollo et al. (2008). Due to anomaly cancellation this model has already family-dependent couplings that can induce neutrino decays. In both examples we can accommodate the Majoron-neutrino couplings to be below the upper limit obtained in this analysis.

## 6 Conclusions

The present scenario of standard neutrino oscillation observed by different experiments shows a strong case for non-zero neutrino masses and mixing. A direct consequence of non-zero neutrino masses is that the neutrino can decay, but usually with much longer lifetimes for decays like or . The only expectation to observe sizeable effects from neutrino decays is from , where the decay products are sterile states and the scalar or pseudo-scalar massless boson .

We use the more updated charged and neutral current data of the MINOS experiment for accelerator produced neutrinos and anti-neutrinos to study the impact of a non-zero decay parameter, , in the channel. Using the charged current data only we have found that the best fit scenario is for a non-zero parameter, corresponding to a neutrino lifetime s/eV.

We have found that the 90% C.L. range for is eV and for the mixing angle is in the oscillation with decay scenario. We could understand the effect of the non-zero decaying parameter for the CC data analysis, implying smaller values of and larger values of than the ones in standard oscillation scenario.

Based on the MINOS combined charged and neutral current analysis the best fit for the oscillation with decay model indicates a zero value for the parameter, that is equivalent to the standard oscillation model. Using the combined analysis we could improve the previous constraint on the allowed neutrino decay lifetime s/eV at the 90% C.L., which is the best limit based only on accelerator produced neutrino data. The inclusion of NC data, which is compatible to standard oscillations, makes the decay scenario more constrained than in the CC analysis only.

We showed that the effect due to the T2K data into our analysis of the MINOS data would be small. In fact, the combined MINOS and T2K analysis did not improve the 90% C.L. limit on the . However, the combined analysis results in a range of finite values for the neutrino lifetime with low significance.

The phenomenological model used in our analysis can be connected with the Majoron models that include also sterile neutrino states. Our constraint on the neutrino lifetime can be translated into the Majoron-neutrino coupling upper bound as using the combined charged and neutral current analysis.

## Acknowledgments

We would like to thank C. Castromonte for helping on computational tasks and M. Goodman for valuable discussion and careful reading of the manuscript. A. L. G. G. was supported by CAPES during his Msc. degree project. O. L. G. P. thanks the support of funding grant 2012/16389-1, São Paulo Research Foundation (FAPESP) and ICTP for hospitality. And R. A. G. thanks the support of Fulbright/CAPES, and the hospitality at Argonne National Laboratory.

## References

- Abe et al. (2012) Y. Abe et al. (Double Chooz Collaboration), Phys. Rev. D86, 052008 (2012), arXiv:1207.6632.
- An et al. (2014) F. P. An et al. (Daya Bay Collaboration), Phys. Rev. Lett. 112, 061801 (2014), arXiv:1310.6732.
- Ahn et al. (2012) J. K. Ahn et al. (RENO Collaboration), Phys. Rev. Lett. 108, 191802 (2012), arXiv:1204.0626.
- Pontecorvo (1957) B. Pontecorvo, Sov. Phys. JETP 6, 429 (1957).
- Maki et al. (1962) Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962).
- Gonzalez-Garcia and Nir (2003) M. C. Gonzalez-Garcia and Y. Nir, Rev. Mod. Phys. 75, 345 (2003), hep-ph/0202058.
- Abe et al. (2013) K. Abe et al. (Super-Kamiokande Collaboration), Phys. Rev. Lett. 110, 181802 (2013), arXiv:1206.0328.
- Aartsen et al. (2013) M. Aartsen et al. (IceCube Collaboration), Phys. Rev. Lett. 111, 081801 (2013), arXiv:1305.3909.
- Adamson et al. (2013) P. Adamson et al. (MINOS Collaboration), Phys. Rev. Lett. 110, 251801 (2013), arXiv:1304.6335.
- Minakata et al. (2006) H. Minakata, H. Nunokawa, S. J. Parke, and R. Zukanovich Funchal, Phys. Rev. D74, 053008 (2006), hep-ph/0607284.
- Gando et al. (2013) A. Gando et al. (KamLAND Collaboration), Phys. Rev. D88, 033001 (2013), arXiv:1303.4667.
- Mirizzi et al. (2007) A. Mirizzi, D. Montanino, and P. D. Serpico, Phys. Rev. D76, 053007 (2007), arXiv:0705.4667.
- Zatsepin and Smirnov (1978) G. T. Zatsepin and A. Y. Smirnov, Yad. Fiz. 28, 1569 (1978).
- Lindner et al. (2001) M. Lindner, T. Ohlsson, and W. Winter, Nucl. Phys. B607, 326 (2001), hep-ph/0103170.
- Bandyopadhyay et al. (2001) A. Bandyopadhyay, S. Choubey, and S. Goswami, Phys. Rev. D63, 113019 (2001), hep-ph/0101273.
- Kelso and Hooper (2010) C. Kelso and D. Hooper, arXiv:1009.5996 (2010).
- Lunardini (2010) C. Lunardini, arXiv:1007.3252 (2010).
- Ando (2003) S. Ando, Phys. Lett. B570, 11 (2003), hep-ph/0307169.
- Fogli et al. (2004) G. L. Fogli, E. Lisi, A. Mirizzi, and D. Montanino, Phys. Rev. D70, 013001 (2004), hep-ph/0401227.
- Barger et al. (1999) V. D. Barger, J. G. Learned, S. Pakvasa, and T. J. Weiler, Phys. Rev. Lett. 82, 2640 (1999), astro-ph/9810121.
- Gonzalez-Garcia and Maltoni (2008) M. C. Gonzalez-Garcia and M. Maltoni, Phys. Lett. B663, 405 (2008), arXiv::0802.3699.
- Palomares-Ruiz et al. (2005) S. Palomares-Ruiz, S. Pascoli, and T. Schwetz, JHEP 0509, 048 (2005), hep-ph/0505216.
- Adamson et al. (2010) P. Adamson et al. (MINOS Collaboration), Phys. Rev. D81, 052004 (2010), arXiv:1001.0336.
- Raghavan et al. (1988) R. S. Raghavan, X.-G. He, and S. Pakvasa, Phys. Rev. D38, 1317 (1988).
- Berezhiani et al. (1992) Z. G. Berezhiani, G. Fiorentini, M. Moretti, and A. Rossi, Z. Phys. C54, 581 (1992).
- Das and Pulido (2011) C. R. Das and J. Pulido, Phys. Rev. D83, 053009 (2011), arXiv:1007.2167.
- Eguchi et al. (2004) K. Eguchi et al. (KamLAND Collaboration), Phys. Rev. Lett. 92, 071301 (2004), hep-ex/0310047.
- Aharmim et al. (2004) B. Aharmim et al. (SNO Collaboration), Phys. Rev. D70, 093014 (2004), hep-ex/0407029.
- Berezhiani and Smirnov (1989) Z. G. Berezhiani and A. Y. Smirnov, Phys. Lett. B220, 279 (1989).
- Lindner et al. (2002) M. Lindner, T. Ohlsson, and W. Winter, Nucl. Phys. B622, 429 (2002), astro-ph/0105309.
- Ando (2004) S. Ando, Phys. Rev. D70, 033004 (2004), hep-ph/0405200.
- Joshipura et al. (2002) A. S. Joshipura, E. Masso, and S. Mohanty, Phys. Rev. D66, 113008 (2002), hep-ph/0203181.
- Beacom and Bell (2002) J. F. Beacom and N. F. Bell, Phys. Rev. D65, 113009 (2002), hep-ph/0204111.
- Hirata et al. (1988) K. S. Hirata, T. Kajita, M. Koshiba, M. Nakahata, Y. Oyama, et al., Phys. Rev. D38, 448 (1988).
- Bionta et al. (1987) R. M. Bionta, G. Blewitt, C. B. Bratton, D. Casper, A. Ciocio, et al., Phys. Rev. Lett. 58, 1494 (1987).
- Frieman et al. (1988) J. A. Frieman, H. E. Haber, and K. Freese, Phys. Lett. B200, 115 (1988).
- Meloni and Ohlsson (2007) D. Meloni and T. Ohlsson, Phys. Rev. D75, 125017 (2007), hep-ph/0612279.
- Maltoni and Winter (2008) M. Maltoni and W. Winter, JHEP 0807, 064 (2008), arXiv:0803.2050.
- Baerwald et al. (2012) P. Baerwald, M. Bustamante, and W. Winter, JCAP 1210, 020 (2012), arXiv:1208.4600.
- Dorame et al. (2013) L. Dorame, O. G. Miranda, and J. W. F. Valle, Frontiers in Physics 1 (2013), ISSN 2296-424X, arXiv:1303.4891, URL http://www.frontiersin.org/high-energy_and_astroparticle_physics/10.3389/fphy.2013.00025/abstract.
- Berryman et al. (2014) J. M. Berryman, A. de Gouvea, D. Hernandez, and R. L. N. Oliveira (2014), arXiv:1407.6631.
- Michael et al. (2008) D. G. Michael et al. (MINOS Collaboration), Nucl. Instrum. Meth. A596, 190 (2008), arXiv:0805.3170.
- Adamson et al. (2011a) P. Adamson et al. (MINOS Collaboration), Phys. Rev. Lett. 107, 011802 (2011a), arXiv:1104.3922.
- Adamson et al. (2008) P. Adamson et al. (MINOS Collaboration), Phys. Rev. D77, 072002 (2008), arXiv:0711.0769.
- Abe et al. (2014) K. Abe et al. (T2K Collaboration), Phys. Rev. Lett. 112, 181801 (2014), arXiv:1403.1532.
- Adamson et al. (2011b) P. Adamson et al. (MINOS Collaboration), Phys. Rev. Lett. 106, 181801 (2011b), arXiv:1103.0340.
- Giganti et al. (2014) C. Giganti et al. (for the T2K Collaboration), NOW Conference (2014), URL http://www.ba.infn.it/~now/now2014/web-content/TALKS/cWed/Plen/NOW_2014_T2K.pdf.
- Cogollo et al. (2008) D. Cogollo, H. Diniz, C. de S. Pires, and P. Rodrigues da Silva, Eur. Phys. J. C58, 455 (2008), arXiv:0806.3087.