# Rescuing Quadratic Inflation

###### Abstract

Inflationary models based on a single scalar field with a quadratic potential
are disfavoured by the recent Planck constraints on the scalar index, , and the tensor-to-scalar ratio
for cosmological density perturbations, . In this paper we
study how such a quadratic inflationary model can be rescued by postulating additional fields with
quadratic potentials, such as might occur in sneutrino models, which might serve as either
curvatons or supplementary inflatons. Introducing a second scalar field reduces but does not
remove the pressure on quadratic inflation, but we find a sample of three-field models that
are highly compatible with the Planck data on and . We exhibit a specific three-sneutrino example that is
also compatible with the data on neutrino mass difference and mixing angles.

KCL-PH-TH/2013-40, LCTS/2013-27, CERN-PH-TH/2013-293

## 1 Introduction

Inflation is a very promising paradigm for the behaviour of the scale factor in the early Universe, which offers a solution to the cosmological horizon problem, explaining the observed large-scale isotropy and giving rise to small fluctuations in energy density as observed in the cosmic microwave background (CMB) radiation that seeded galactic structure formation, as well as explaining the absence of topological defects. In the simplest models of inflation, the energy density responsible for the accelerated early expansion of the scale factor arises from the expectation value of a scalar field rolling down a potential. As discussed in [1], a huge number of candidate scalar fields appearing in extensions of the standard model of particle physics and in models with modified gravitational sectors have the freedom to match observations: for a sampling of recent models, see [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12].

In the simplest class of single-field models, the energy density in the potential causes an accelerated expansion of the Universe, with the strength of tensor perturbations being directly related to the magnitude of the energy density. Since the inflaton rolls slowly down the potential, the spectrum of scalar perturbations is tilted. Once the value of becomes sub-Planckian, the slow roll ends and the field starts to oscillate and then decays into radiation, reheating the Universe. This simplest model makes surprisingly successful predictions for the main inflationary observables. For example, it provides an almost scale-invariant spectrum of density perturbations and predicts that the non-Gaussianity parameter is relatively small.

The simplest scenario is a single scalar field with a monomial potential such as or . However, models were already under pressure from the results of WMAP [13] and can now be regarded as excluded by the observations of the Planck satellite [14, 15], and quadratic models are now also under severe pressure, since they are unable to fit simultaneously the observed spectral index of scalar perturbations and the low tensor-to-scalar ratio .

For a range of typical values of the number of efolds during inflation, , single-field inflationary models with a potential predict a value for in the range , whereas the Planck data require . For a lower value of , i.e., a smaller amount of expansion between the largest scales leaving the horizon and the end of inflation, it is possible for a quadratic model to give a value for that obeys the Planck constraint, but the quadratic model then makes a prediction for the spectral index of the perturbations that does not respect the Planck constraint , since there is a one-to-one mapping between the number of efolds and the scalar spectral index.

In this work we explore ways to rescue models with quadratic potentials. Our objective is to find a model of quadratic inflation that (i) gives the correct magnitude for the scalar density perturbation , (ii) yields a Planck-compatible value of the scalar spectral index , (iii) predicts a value for the tensor-to-scalar ratio that obeys the Planck constraint, and (iv) does not predict significant non-Gaussianity.

Attaining all these objectives simultaneously requires modifying the mapping between the number of efolds and the scalar spectral index. This may be done by introducing extra fields that act either as secondary inflatons or as curvatons. We restrict our attention to minimal models with multiple potentials. Apart from simplicity, one of our motivations for this restriction is the idea that the inflaton might be a singlet sneutrino in a supersymmetric seesaw model of neutrino masses [16], which is a motivated and minimal extension of the Standard Model. In such a scenario, there must be at least two such sneutrino fields so that two light neutrinos can acquire masses, and probably three so that all three light neutrinos can be massive.

In the next Section we review the simplest case of a single inflationary field, and the above-mentioned reasons why it has difficulties fitting the latest data. In Section 3 we go on to examine models with an additional scalar field acting as a curvaton. Whilst the disagreement with the data can be reduced in such scenarios, they are still in some tension with the data, although not enough to be considered as excluded. Then, in Section 4 we look at models with two inflatons acting in consort accompanied by an curvaton, which we find to be the minimal combination of quadratic fields that obeys completely all the Planck constraints. In Section 5 we show how such a supersymmetric seesaw model can be realized with the three quadratic scalar fields being interpreted as supersymmetric partners of heavy singlet neutrinos, showing how such a model can fit simultaneously the available data on neutrino oscillations. Finally, Section 6 contains a summary and discusses future prospects for quadratic inflationary models.

## 2 Single-Field Quadratic Inflation

We first review the standard lore for single-field inflation with a quadratic potential . The power spectrum of the scalar perturbations is given by [17, 18]

(1) |

where starred quantities correspond to the values of the various parameters at the epoch at which the last scales to re-enter the horizon before the last scattering surface were initially leaving the horizon during inflation. The Hubble parameter is given by the Friedmann equation , where is the reduced Planck mass, GeV. The density perturbation is [17], so using the above equation we have

(2) |

The two slow-roll parameters in this case are:

(3) |

where dashed quantities denote derivatives with respect to .

The scalar spectral index is given in terms of the slow-roll parameters by [19], so in this case it is given by

(4) |

The tensor-to-scalar ratio is defined as the ratio of the power spectrum of the tensor perturbations [20]

(5) |

to the power spectrum of the scalar perturbations, given in (1). In the simple quadratic inflation model the tensor-to-scalar ratio is given by

(6) |

yielding a very restrictive relationship between the tensor-to-scalar ratio and the scalar spectral index:

(7) |

Ultimately, both quantities are fixed by the expectation value of the inflaton corresponding to the scales entering the horizon at the last scattering surface, which is around 40-60 efolds before the end of inflation. In a particular model the exact number of efolds is not a free parameter, but depends upon the rate at which the coherent oscillations of the inflaton decay into radiation after inflation.

To find the time when the largest scale crosses the horizon during inflation, we must solve the equation

(8) |

where is the potential at the end of inflation and is the energy density evaluated at the time of the decay of the longest-lived field in the theory [17]. Solving (8) is quite simple for quadratic inflation, and means that the number of efolds , and therefore the time of horizon crossing , is not an arbitrary parameter but is constrained for a given set of masses and decay rates. A change in the duration of matter or radiation domination of the energy density up until the last decay when the Universe is thermalized for the last time will affect the value of . In the case of single-field (quadratic) inflation, the two model parameters are therefore the inflaton mass and its decay rate , as shown in the first row of Table 1.

Free parameters for single-field quadratic inflation.

Free parameters for the model with a quadratic inflaton and a quadratic curvaton.

Free parameters for the model with two quadratic inflatons and a quadratic curvaton.

## 3 Model with one Inflaton and one Curvaton

We now consider a second model that contains, in addition to the field with a quadratic potential that acts as the inflaton, a second scalar field with a quadratic potential that acts as a the curvaton [21, 22, 23, 24]. This model is therefore described by

(9) |

The curvaton field is frozen at an expectation value during the slow-roll of since . After inflation ends and the Hubble parameter becomes comparable to the curvaton mass, , the curvaton begins to oscillate around the minimum of its potential. When the Hubble parameter falls to the curvaton decay rate, , the curvaton decays to radiation. The parameters of this model are the masses and decay rates of the inflaton and the curvaton, , and , respectively, and the expectation value of the curvaton during inflation , as listed in the second row of Table 1.

Now, although only is responsible for providing the expansion of the Universe during inflation, both fields can contribute to , , and . From [25, 26] we see that the total density perturbation is given by

(10) |

with the individual components evaluated at horizon crossing ^{1}^{1}1All quantities
evaluated at this time have the index *., where we have

(11) |

and is evaluated at the epoch of the last decay, i.e., the epoch at which the curvaton decays, rather than the inflaton, so that

(12) |

where at the epoch of curvaton decay [27, 28, 29]. The power spectrum of the created density perturbation is then

(13) |

and to evaluate the spectral index we simply use the definition

(14) |

where we transformed the derivative using , and note that the expansion is dominated by the effects of at this epoch, which is during slow-roll inflation. The tensor-to-scalar ratio is

(15) |

which may be evaluated using (13) and

(16) |

Once the time has been established in the same way as in Section 2, we can calculate the density perturbation using (10), the spectral index using (14), and the tensor-to-scalar ratio using (15). Employing a Markov Chain Monte Carlo algorithm, we fit the measured values of and and explore the predictions for , as shown in Figs. 1 and 2.

We find in our sampling of its five-parameter space that this model is able to predict a low value for the tensor-to-scalar ratio , well under the constraint coming from Planck. This happens when the curvaton field is significant at the time of the last decay, and dominates the energy density. However, obtaining a low value for is possible only at the expense of an unsuccessful prediction for the scalar spectral index. As is shown in the left panel of Fig. (2), when falls to a value below , grows to a value that is approximately away from the current measurement as .

This conclusion is qualitatively the same as the results displayed in fig. 6 of [30] and fig. 1 of [31]. In both cases it was shown that in a model with one inflaton and one curvaton, a low predicted value for is possible but it leads to a higher value of . However, that analysis was done for a set range of the number of efolds as well as a set value of the curvaton’s contribution to the perturbations. Our results, shown in Figs. 1 and 2, differ from the results in [28], where a Planck-consistent value of was found in a fit to the values of and in a model with one inflaton and one curvaton. The reason for this difference is that in [28] the inflationary scale at horizon crossing was set by hand by assuming , which corresponds to a very low number of efolds, . In contrast, in our work we solve for the field value at , so that our model is self-consistent. We find that has to be in the range as shown in the right plot of Fig. (2).

We find that the values of and depend strongly on the value of the inflaton mass, . As is shown in the left panel of Fig. 3, a lower value of can be obtained with a lower value of , but this corresponds to a higher value of . The same is observed in fig. 6 and 7 of [32]; for an analytic solution, the authors study and in terms of the curvaton’s expectation value and the ratio of the decay rates of the inflaton and the curvaton. They find that the combinations of parameters that result in a Planck-consistent value of are associated with a value of that is significantly higher than the measured one, which is what we find as well. This behaviour is seen also in Fig. 1, where one can see that, although curvaton models fit the Planck data somewhat better than single-field quadratic inflation, they are still disfavoured, as they do not come within the 68% CL region favoured by the Planck + BAO data.

We therefore conclude that, although a model with one inflaton and one curvaton, with both fields contributing to the generation of the perturbations, is able to predict a value for that is within the Planck constraint, this model cannot combine this good prediction with a successful prediction for the spectral index .

This conclusion could be relaxed by assuming additional inflation at some later time in the Universe as in, e.g., thermal inflation [17, 20]. However, it seems quite difficult to obtain enough thermal inflation to reconcile the Planck data on and with the model combining a quadratic inflaton and a quadratic curvaton, since typically one requires 25-30 efolds of later inflation.

## 4 Model with two Inflatons and one Curvaton

In view of this setback, in this Section we further augment the model by including another scalar field, . This is assumed to be a second inflaton field, again with a quadratic term in the potential, so that it becomes

(17) |

In this case, both and roll slowly down their potentials, giving rise to inflation while, as before, the curvaton field remains frozen at its expectation value during inflation and then oscillates and decays into radiation.

Curvaton domination of the energy density in this model can lead to a low value of consistent with the Planck data, while the addition of the second inflaton helps to lower the prediction for , making the model consistent with the Planck value. This is because, in the case of curvaton dominance, the power spectrum of the scalar density perturbations is simply . Then, from the expression (14) for the spectral index we obtain the expression:

(18) |

We see from (18) that, in order to alter the high value of that was obtained in the two-field model we must alter the time derivative of the potential evaluated at horizon crossing. Since the curvaton expectation value does not change during inflation, it is only with the addition of the second inflaton that we can achieve this.

Substituting (17) in (18), we obtain the following expression for the scalar spectral index in terms of our model parameters:

(19) |

In our numerical analysis, we find that the heavier of the two inflaton fields (let us call it ) stops rolling once its expectation value drops down to the Planck scale. It then begins to oscillate, but the Universe is still expanding exponentially due to the other inflaton field. This means that is redshifted, with its expectation value dropping rapidly, and by the end of inflation, i.e., when the expectation value of becomes equal to , the contribution of to the energy density is negligible. It is therefore not necessary to include the decay rate of the heavy inflaton among the relevant parameters of our model listed in the third row of Table 1. The parameter corresponds to the angle in field space at , i.e. the angle between the and direction when the largest scales in the CMB are leaving the horizon. Typically the lighter field will be frozen at this moment, although we do not assume this, since we analyze this stage numerically.

As in the model with one inflaton and one curvaton, we solve equation (8) numerically to find the time during inflation when the largest scales left the horizon and determine the expectation values and at that time. Since we are looking at the cases when the curvaton energy density dominates the Universe at the time of its decay, we can compute the density perturbation from the curvaton component of (10) and the spectral index from (18). The prediction for the tensor-to-scalar ratio is again calculated from (15), but in this case the spectrum of scalar perturbations is given by the second term on the right-hand side of (13):

(20) |

where in the case of curvaton domination will be very close to . This results in the expression

(21) |

for the tensor-to-scalar ratio.

As in the previous Section, we employ a Markov Chain Monte Carlo algorithm to study the predictions of this model for , while fitting the measured values of and . We find that this model with two inflatons and one curvaton is indeed able to accommodate values of that are in agreement with the Planck results. The left panel of Fig. 5 shows the spectral index and the tensor-to-scalar ratio of this 3-field model, and the model is compared to the Planck data in Fig. 7.

Fig. 5 confirms what we expect from (21), namely that a lower value of can be obtained from a lower curvaton expectation value . However, we also find that there is a minimum in the prediction for the value of . As the curvaton expectation value becomes smaller, it becomes more difficult for the curvaton to dominate completely the energy density of the Universe at the time of its decay. Eventually, for GeV the parameter measuring the significance of the curvaton energy density, , becomes smaller than . We see in Fig. 5 that this corresponds to a minimum value of around . For even smaller curvaton expectation values, the curvaton energy density becomes less significant and rapidly becomes smaller. This means that the model becomes more and more like the single-field inflation model. In this case we see that the value of begins to grow, as expected. In Fig. 5 we see that the number of efolds is slightly lower than in the case with one inflaton and one curvaton, and that in order to obtain a good spectral index we require to be greater than two or three times .

It is essential that all the fields in our model decay before nucleosynthesis starts. This requirement places a lower constraint on the temperature of the Universe when the last particle in our model decays into radiation, namely MeV [33].To this end, we calculate the reheating temperature from the total energy density at the time of the last particle’s decay:

(22) |

We find that, in order to get good inflationary parameters, we typically find a low reheating temperature, e.g., - GeV, that resolves the cosmological gravitino problem while still respecting the nucleosynthesis constraint.

Finally, we consider the prediction of this model for the non-Gaussianity parameter . For this model to be consistent with the Planck results, it needs to predict a value in the range . This parameter depends on the sequence of oscillations and decays in the model, and is evaluated at the time of the last decay [29, 34, 35]. In the model with two inflatons and one curvaton, we have seen that the energy density of the heavier inflaton is insignificant from the end of inflation onwards. We also note that for such a model to predict a low value of the tensor-to-scalar ratio that is in agreement with the Planck constraint, the curvaton must dominate the energy density when it decays. Therefore, after the end of inflation, we have the equivalent of a model with one inflaton and one curvaton, with the curvaton energy density being the dominant component of the energy density of the Universe, i.e., at the epoch of last decay. It was shown in [29] that models with one inflaton and one curvaton, with both particles contributing to the perturbations, cannot predict a large when the curvaton dominates before its decay. In Fig. 2 of [29] we see that, for the range of curvaton expectation values that we have in our model and for , the non-Gaussianity parameter is small, . So this model with two inflatons and one curvaton is in agreement with the constraint imposed by the Planck data.

## 5 Three-Sneutrino Inflation

As already mentioned, part of our motivation for studying extended models of quadratic
inflation was the possibility of rescuing sneutrino inflation. In most supersymmetric
seesaw models for the light neutrino masses, there are three heavy singlet supersymmetric partners of
(right-handed) neutrinos ^{2}^{2}2The neutrino oscillation data require only two light neutrinos
to have non-zero masses, which is possible in principle with just two heavy singlet neutrinos in a Type-I seesaw model.,
each with a quadratic effective potential. In this Section we
investigate the possibility of identifying the inflatons and the curvaton of the previous
Section with these sneutrinos.
Typical masses of the singlet (right-handed) neutrino
fields are of order - GeV, without extreme fine-tuning of the Yukawa
couplings ^{3}^{3}3In the Standard Model these range
from order unity for the top quark down to order for the electron, so a wide range of possible values could be
considered.. We look here for models that can explain neutrino masses
and mixing angles as well as provide good parameters for Planck observables.

We start with the following mass matrix for the singlet (right-handed) (s)neutrinos ^{4}^{4}4Supersymmetry
breaking is not important for our analysis.:

(23) |

and a Yukawa matrix:

(24) |

The decay rate of each heavy particle is then given by [36]

(25) |

where .

The usual left-handed neutrino mass matrix depends upon the masses and Yukawa coupllings of the heavy neutrinos, and is given by:

(26) |

By diagonalizing this expression we obtain three eigenvalues, which are the three neutrino masses, and three eigenvectors, which determine the unitary mixing matrix :

(27) |

Here is the leptonic equivalent of the CKM matrix in the quark sector [37, 38], which can be written in the following parametrization:

(28) |

where and . This matrix contains three neutrino mixing angles . In general, the mixing matrix can also contain CP-violating phases, one detectable in neutrino oscillations and two Majorana phases that would affect neutrinoless double- decay. For simplicity, here we discard these phases and assume a real Yukawa matrix for the sneutrinos.

We use the 3 masses to fit the two measured light-neutrino mass-squared differences and . We then choose to fit the three mixing angles , and , as measured in neutrino oscillation experiments [39].

Within this framework, we display one illustrative model with two sneutrino inflatons and one curvaton that leads to predictions for the tensor-to-scalar ratio and the scalar spectral index that are consistent with the Planck results. The parameters of this model are outlined in Table 2, and the corresponding point in parameter space is shown as a black cross in the right panel of Fig. 7.

Parameter | Value |
---|---|

2.6 GeV | |

2.2 GeV | |

2.2 GeV | |

2.2 GeV | |

1.5 GeV | |

2.4 GeV | |

980 GeV | |

1.8 GeV | |

1.2 GeV | |

2.7 GeV | |

1.8 | |

0.9592 | |

0.042 | |

4.85 eV | |

4.92 eV | |

3.1 eV | |

see Fig. 6 | |

see Fig. 6 |

## 6 Conclusions

In this paper we have looked for the simplest possible fit to the Planck data using only fields with quadratic potentials. We first recalled the well-known result that single-field quadratic inflation is under pressure from the Planck data. We then went on to show that, while quadratic inflation with a quadratic curvaton fits the data slightly better, it also is disfavoured in its simplest form. This is because the power spectrum is fixed during the normal single-field inflationary phase, and requiring enough efolds forces one to focus on the region where is too high. We went on to show that with three quadratic potentials we can fit both and . We also found that such a model provides a minimum in the prediction for the value of . Finally, we exhibited a model where the three quadratic fields are identified with three sneutrino fields, two playing the rôles of inflatons and one being a curvaton, displaying one example of a point in parameter space that fits both the neutrino and cosmological data.

Our results show that it is possible to rescue quadratic inflation, and that this does not require a very exotic model. Indeed, the three fields required can be identified with singlet (right-handed) sneutrinos.

Acknowledgments

The authors would like to thank David Mulryne for his very useful comments on the first version of this paper.
The work of J.E. was supported in part by
the London Centre for Terauniverse Studies (LCTS), using funding from
the European Research Council
via the Advanced Investigator Grant 267352.
J.E. and M.F. are grateful for funding from the Science and Technology Facilities Council (STFC).

## References

- [1] J. Martin, C. Ringeval and V. Vennin, arXiv:1303.3787 [astro-ph.CO].
- [2] D. Croon, J. Ellis and N. E. Mavromatos, Physics Letters B 724 (2013) , 165 [arXiv:1303.6253 [astro-ph.CO]].
- [3] J. Ellis, D. V. Nanopoulos and K. A. Olive, Phys. Rev. Lett. 111 (2013) 111301 [arXiv:1305.1247 [hep-th]] and JCAP 1310 (2013) 009 [arXiv:1307.3537 [hep-th]].
- [4] K. Nakayama, F. Takahashi and T. T. Yanagida, JCAP 1308 (2013) 038 [arXiv:1305.5099 [hep-ph]] and arXiv:1311.4253 [hep-ph].
- [5] R. Kallosh and A. Linde, JCAP 1306 (2013) 028 [arXiv:1306.3214 [hep-th]]; JCAP 1307 (2013) 002 [arXiv:1306.5220 [hep-th]].
- [6] W. Buchmuller, V. Domcke and K. Kamada, Phys. Lett. B 726 (2013) 467 [arXiv:1306.3471 [hep-th]].
- [7] F. Farakos, A. Kehagias and A. Riotto, Nucl. Phys. B 876 (2013) 187 [arXiv:1307.1137 [hep-th]].
- [8] D. Roest, M. Scalisi and I. Zavala, JCAP 1311 (2013) 007 [arXiv:1307.4343 [hep-th]].
- [9] E. Kiritsis, JCAP 1311 (2013) 011 [arXiv:1307.5873 [hep-th]].
- [10] J. Ellis and N. E. Mavromatos, Phys. Rev. D 88 (2013) 085029 [arXiv:1308.1906 [hep-th]].
- [11] P. Fre and A. S. Sorin, arXiv:1308.2332 [hep-th] and arXiv:1310.5278 [hep-th].
- [12] T. Li, Z. Li and D. V. Nanopoulos, arXiv:1310.3331 [hep-ph] and arXiv:1311.6770 [hep-ph].
- [13] G. Hinshaw et al. [WMAP Collaboration], Astrophys. J. Suppl. 208 (2013) 19 [arXiv:1212.5226 [astro-ph.CO]].
- [14] P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5062 [astro-ph.CO].
- [15] P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5082 [astro-ph.CO].
- [16] J. R. Ellis, M. Raidal and T. Yanagida, Phys. Lett. B 581 (2004) 9 [hep-ph/0303242].
- [17] A. R. Liddle, D. Lyth, Cosmological Inflation and Large-Scale Structure, (CAMBRIDGE UNIVERSITY PRESS, 2000)
- [18] C. T. Byrnes and D. Wands, Phys. Rev. D 74 (2006) 043529 [astro-ph/0605679].
- [19] D. Wands, N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 66 (2002) 043520 [astro-ph/0205253].
- [20] A. R. Liddle and D. H. Lyth, Phys. Rept. 231 (1993) 1 [astro-ph/9303019].
- [21] K. Enqvist and M. S. Sloth, Nucl. Phys. B 626 (2002) 395 [hep-ph/0109214].
- [22] D. H. Lyth, C. Ungarelli and D. Wands, Phys. Rev. D 67 (2003) 023503 [astro-ph/0208055].
- [23] T. Moroi, T. Takahashi and Y. Toyoda, Phys. Rev. D 72 (2005) 023502 [hep-ph/0501007].
- [24] K. Ichikawa, T. Suyama, T. Takahashi and M. Yamaguchi, Phys. Rev. D 78 (2008) 023513 [arXiv:0802.4138 [astro-ph]].
- [25] D. H. Lyth and D. Wands, Phys. Lett. B 524 (2002) 5 [hep-ph/0110002].
- [26] K. Dimopoulos, G. Lazarides, D. Lyth and R. Ruiz de Austri, Phys. Rev. D 68 (2003) 123515 [hep-ph/0308015].
- [27] K. Enqvist and T. Takahashi, JCAP 0809 (2008) 012 [arXiv:0807.3069 [astro-ph]].
- [28] J. Fonseca and D. Wands, JCAP 1206 (2012) 028 [arXiv:1204.3443 [astro-ph.CO]].
- [29] M. Sueiro and M. Fairbairn, Phys. Rev. D 87 (2013) 4, 043524 [arXiv:1208.0559 [astro-ph.CO]].
- [30] T. Moroi and T. Takahashi, Phys. Rev. D 72 (2005) 023505 [astro-ph/0505339].
- [31] K. Enqvist and T. Takahashi, JCAP 1310 (2013) 034 [arXiv:1306.5958 [astro-ph.CO]].
- [32] J. Meyers and E. R. M. Tarrant, arXiv:1311.3972 [astro-ph.CO].
- [33] E. Kolb, M. Turner The Early Universe, (Westview Press, 1994).
- [34] H. Assadullahi, J. Valiviita and D. Wands, Phys. Rev. D 76 (2007) 103003 [arXiv:0708.0223 [hep-ph]].
- [35] D. Langlois, Prog. Theor. Phys. Suppl. 190 (2011) 90 [arXiv:1102.5052 [astro-ph.CO]].
- [36] J. R. Ellis, Nucl. Phys. Proc. Suppl. 137 (2004) 190 [hep-ph/0403247].
- [37] K. Nakamura et al. JPG 37, 075021 (2010).
- [38] W. Rodejohann, Pramana 72 (2009) 217 [arXiv:0804.3925 [hep-ph]].
- [39] C. Cheung, L. J. Hall and D. Pinner, arXiv:1103.3520 [hep-ph].