GUT-scale inflation with sizeable tensor modes
A sizeable tensor-to-scalar ratio, such as recently claimed by BICEP2, would imply a scale of inflation at the typical scale of supersymmetric grand unification. This could be an accident, or strong support for supersymmetric theories. Models of F-term hybrid inflation naturally connect the GUT scale with the inflationary scale, but they also predict the tensor-to-scalar ratio to be unmeasurably small. In this work we analyze a general UV embedding of F-term hybrid inflation into a supergravity theory with a general Kähler potential. The CMB observables are generated during the early phase of inflation, at large inflaton values, where the potential is dominated by Planck-suppressed operators. Tuning the leading higher-order terms can give an inflaton potential with sizeable tensor fluctuations and a field excursion which is still sub-Planckian but close to the Planck scale, as expected from the Lyth bound.
a]Felix Brümmer, a]Valerie Domcke b]and Verónica Sanz Prepared for submission to JCAP
GUT-scale inflation with sizeable tensor modes
SISSA/INFN, Via Bonomea 265, 34136 Trieste, Italy
Department of Physics and Astronomy, University of Sussex, Brighton BN1 9QH, UK
Recently the BICEP2 collaboration reported a discovery of B-mode polarization in the cosmic microwave background . When interpreted as originating from primordial gravitational waves from inflation, this signal corresponds to a tensor-to-scalar ratio of around , and to a scale of inflation around GeV. This claim is currently under intense scrutiny; however, even if the amplitude of the primordial tensor fluctuations were an order of magnitude smaller than reported, this would still point towards a scale of inflation near the scale of supersymmetric grand unification. Hence, a detection of primordial gravity waves immediately has two remarkable consequences: First, it suggests a connection between inflation and SUSY GUTs. Second, assuming single-field slow-roll inflation with a potential whose slope increases monotonically until slow-roll is violated at the end of inflation, it implies super-Planckian field values during inflation , .
These two implications are difficult to reconcile. On the one hand, at such large field values, the theory is generically dominated by uncontrolled higher-dimensional operators. They can be forbidden by symmetries, but these symmetries, against generic expectations, would need to be respected by quantum gravity. Hence inflation might be governed by physics which is not described by supergravity as an effective field theory, and the appearance of the GUT scale might be a pure coincidence. On the other hand, taking the hint for linking inflation and supersymmetric grand unification seriously typically leads to small-field models. Inflation is linked to a GUT-breaking phase transition, with field values rendering these models calculable in an effective supergravity framework. A prototypical example is -term hybrid inflation (FHI) [3, 4]: The inflaton potential is essentially logarithmic, which allows to reach 50 - 60 e-folds within a small field range but at the same time predicts to be very small.
In this paper, we aim to find a compromise between these two seemingly conflicting observations. We take the connection of inflation and supersymmetric grand unification seriously, by studying a variant of FHI which is neither strictly small-field nor strictly large-field inflation. When allowing the inflaton to reach values close to (but not exceeding) the Planck scale, -suppressed operators become important. As was shown in  for general models of inflation, these may significantly affect the shape of the inflaton potential and thus the prediction for . Relaxing all assumptions on the monotonicity of the inflaton potential allows us to realize with sub-Planckian field values, avoiding the bound on the field excursion of  but in accordance with the weaker bound of  and the studies of [7, 8, 9]. See also Ref.  for related work in the context of natural inflation. Simultaneously, the Planck- suppressed operators can account for a sufficiently red-tilted spectral index, , cf. .
In the context of FHI we can obtain suitable potentials by allowing for -symmetry breaking terms in the Kähler potential, suppressed by powers of .111The effects of Planck-suppressed operators from -symmetry preserving non-canonical Kähler terms were recently studied in . However, the result of this study was that remains small even when taking those terms into account: generally for sub-Planckian inflaton values. This is in line with the expectation that the continuous global -symmetry which governs the superpotential of ordinary FHI will be broken by quantum gravity effects. The inflaton potential reduces to that of FHI in the global SUSY limit . During the early inflationary phase which determines the cosmic microwave background (CMB) observables, it is however completely dominated by higher-dimensional operators (whose coefficients must be suitably tuned to obtain the desired values for and without violating the Planck bounds on the running, and the running of the running, of the spectral index). We find that, working with terms up to the order , our model can saturate the Lyth bound , leading to inflaton values for . For such large values of , terms of even higher order are therefore generically not under good control. Our conclusion is that this model can be reconciled with the Planck and BICEP2 data, but only by judiciously choosing the parameters entering the scalar potential up to rather high order in the expansion.
2 F-term hybrid inflation and its supergravity embedding
Recall that the superpotential in -term hybrid inflation is linear in the inflaton superfield ,
There is a global continuous symmetry under which carries charge and , are neutral. During the inflationary phase, the inflaton superpotential reduces to
where is the scale of inflation. With a canonical Kähler potential, and in the limit , the scalar potential is exactly flat at the tree level. At the one-loop level, there is a logarithmic correction coming from and , which owe their mass to the inflaton,
Here is a group-theoretical prefactor and is the renormalization scale. The tensor-to-scalar ratio is small in this model, even when incorporating higher-order terms from supergravity and from a non-canonical Kähler potential respecting .
The -symmetry must, at the latest, be broken at some point after inflation in order to obtain an (almost) Minkowski vacuum with broken supersymmetry. Depending on the model, it may also be broken explicitly by renormalizable inflaton couplings to other fields, as in the model of  which we review in Appendix A. In any case we do not expect quantum gravity effects to respect global symmetries, so it is reasonable to supplement the superpotential Eq. (2.2) by a Kähler potential of the form
where we have assumed to be canonically normalized up to quadratic order. We also allow for a constant term in the superpotential,
Note the absence of terms such as or in . This structure could be the result of the -symmetry being broken in a separate sector with vanishing or at most very small couplings to the inflaton sector. Gravitational physics will still communicate -breaking to the inflaton, but only in the form of -suppressed operators in as in Eq. (2.4). (Of course, some of these operators can partly be absorbed in by a Kähler-Weyl transformation , , but the resulting corrections to are always suppressed by powers of and thus very small, as opposed to, say, a inflaton mass term which would completely upset the model.)
Since all -symmetry breaking corrections are suppressed by powers of , we recover FHI in the rigid limit. On the other hand, for inflaton values of the coefficients and the free parameter222It is tempting to identify with the gravitino mass as expected for TeV-scale supersymmetry, and to therefore impose . Here we will make no such assumption, since SUSY may as well be broken at a higher scale, or there may be other contributions to the gravitino mass after the end of inflation. For instance, in the model of , must be of the order of to cancel the cosmological constant after inflation has ended (see Appendix A for details). We will however assume throughout this paper. can take values which allow the CMB observables to significantly deviate from the FHI predictions, as we will detail in the next Section.
The scalar potential is
where we approximate by the one-loop Coleman-Weinberg potential in the globally supersymmetric limit, cf. Eq. (2.3), discarding terms that are doubly suppressed by both and by a loop factor.
The -symmetry breaking terms in the superpotential and in the Kähler potential break the degeneracy in the phase of the complex inflaton field , turning FHI into a two-field inflation model . However, the real axis remains a self-consistent solution provided that and the are real, and in order to prove that this setup can account for a large tensor-to-scalar ratio, it will be sufficient to focus on this solution. We leave the analysis of the full two-field model to future work. The canonically normalized inflaton field is now
and its scalar potential can be written as
where the coefficients are function of the parameters in the Kähler potential, and we have defined
The expressions for the coefficients of the higher-order terms quickly become very unwieldy in this expansion. For practical purposes it is more convenient to work directly in terms of the coefficients appearing in the scalar potential, as we will do in the following. Any potential in which the -suppressed terms take the general form of Eq. (2.8) can be obtained by choosing and suitably. In Appendix B, we give a translation between the terms in the Kähler and the leading coefficients in Eq. (2.8).
3 An upper bound on the tensor-to-scalar ratio
Our goal is now to find a parameter region in which the potential Eq. (2.8) gives rise to inflation with a relatively large tensor-to-scalar ratio and a relatively small field excursion , such that higher-order terms in the expansion are under control. The free parameters entering Eq. (2.8) are the global SUSY vacuum energy during inflation defined by Eq. (2.9), the superpotential coupling constant , the constant term in the superpotential , and the coefficients of the last term in Eq. (2.8), which we shall truncate at leaving us with the parameters . Moreover, in the following we set (a different value can always be absorbed into a redefinition of and ) and (our results are not sensitive to the precise choice of the renormalization scale).
As it will turn out, truncating the series at the order does not mean that higher terms can safely be neglected for generic values of the coefficients . By definition it is clear that, for any small-field model, eventually higher terms become negligible as long as the sequence of the is well-behaved, but in our case (as should become clear below) this point is only reached at much higher order. In fact we will need to fine-tune the leading coefficients, and setting to zero merely corresponds to a specific choice for this fine-tuning.
To determine a viable parameter region, observe that the slow-roll parameters and , defined by
where and refer to the slow-roll parameters evaluated at , e-folds before the end of inflation. It should be noted that the value of may decrease depending on the foreground dust model which one subtracts, , bringing it into somewhat better agreement with the upper bound from Planck . Future measurements of the B-mode spectrum and possible dust foregrounds will be crucial for a precise determination of .
From now on we assume that the vacuum energy during inflation is dominated by the global SUSY term , i.e. that the second term on the RHS of Eq. (2.8) can be neglected. Then is fixed by the amplitude of the power spectrum ,
Given the structure of the scalar potential Eq. (2.8), the parameter enters only into , but cancels in all the slow-roll parameters and in the slow-roll equation of motion. Thus, after ensuring the correct slow-roll dynamics, we can always determine a posteriori using Eq. (3.3).
We can further restrict the parameter space analytically by confronting our model with the Lyth bound . With no assumptions on the monotonicity properties of the inflaton potential, except that should be approximately constant during the first e-folds which leave their imprint on the CMB observables, the field excursion is bounded from below as (see also )
If we want to succeed in reproducing a large tensor-to-scalar ratio for , we must get at least close to saturating Eq. (3.4). Hence we need a scalar potential which is rather steep for the first e-folds and then quickly becomes very flat to accommodate the remaining e-folds. Achieving this with an analytic single-field inflation potential as in Eq. (2.8) implies that the higher-order derivatives of are generically large (see also ), typically , much larger than preferred by the Planck data . However, since the sign of is not fixed, will have zeros in which the coefficients in Eq. (2.8) conspire so that the higher derivatives are small. Hence without loss of generality we can set the higher-dimensional slow-roll parameters at horizon crossing to an arbitrary value allowed by the Planck data,333Note that the bound on depends on the inclusion of a non-vanishing and . With no bounds given by the Planck collaboration including both additional parameters, we here opt for using the fit including . We expect that including the non-vanishing should not significantly alter the best-fit value, but will possibly enlarge the error bands. thus eliminating two further parameters [17, 15]:
As expected, varying these bounds within the experimental errors does not change the qualitative picture, but does impact the quantitative results somewhat. We find that in order to achieve our goal of large for moderate , the most convenient choice is to set both and to their experimental upper bound. This is in agreement with analyses based on flow equations [18, 19], which indicate that a very flat potential at the end of inflation can be achieved for sufficiently large values of and .
Finally, again exploiting the Lyth bound, we know that the initial value of the inflaton field will need to be sizeable for the -suppressed operators to be relevant, but not too large in order to retain control over the subdominant terms. Fixing and using Eqs. (3.2) to (3.5) to eliminate five of the eight remaining parameters in Eq. (2.8) allows us to restrict the final three free parameters, subject to the consistency condition that e-folds are realized.
Fig. 1 illustrates the result of this analysis, showing the maximal tensor-to-scalar ratio achievable for a given value of in this framework. Starting from the scalar potential in Eq. (2.8), we require e-folds of slow-roll inflation, and to be at their best fit values according to Eqs. (3.2) and (3.3), and to be at their respective upper limit according to Eqs. (3.5). We then perform a parameter scan in the remaining three free parameters, requiring all higher order derivatives to be under control, , with and serving as representative examples for a less and a more conservative bound, respectively. For comparison we also show the Lyth bound, cf. Eq. (3.4), with . We find that our the setup of generalized hybrid inflation discussed in this paper can indeed saturate (and slightly exceed) even the most conservative Lyth bound, allowing for a relatively large tensor-to-scalar ratio for a moderate value of the inflaton field , marginally justifying the expansion in suppressed operators even in the light of the BICEP2 result, at the price of tuning the coefficients of the first few operators.
We should point out here that there are two kinds of observational constraints which might still threaten the validity of these parameter points.444We thank Shaun Hotchkiss for helpful discussions on these issues. Firstly, a too flat potential towards the end of inflation may lead to overproduction of primordial black holes, see e.g.  for an analysis of the resulting bounds. Secondly, the CMB data severely constrains any variation in the power spectrum amplitude over the first few e-folds around , so even the more conservative assumption might turn out not to be conservative enough. We have verified that, while the constraints on primordial black holes do rule out a part of the parameter space, there are nevertheless many valid parameter points left even towards the region of low and low (and in particular, on the curves shown in Fig. 1). As to the second point, addressing it would require a rather involved analysis beyond the slow-roll approximation and beyond the derivative expansion for the potential. A tentative check (which still partly relies on the slow-roll approximation despite briefly becoming in our scenario) seems to indicate that one may need to go to values of even smaller than such as to keep the variation of the power spectrum amplitude under control. However, the situation is not conclusive, and further study is needed to settle this issue.
What characterizes the scalar potentials which lead to large values of for moderate values of ? To achieve a large value of , the potential must feature a rather large first derivative at , when the CMB scales left the horizon, while all higher derivatives should be small to satisfy the Planck constraints. To account for a small value of , this linear behaviour of the potential must transition to a very flat part of the potential which accounts for most of the e-folds at small field values. The total field excursion is minimized if this transition is sharp, typically rendering higher-order derivatives large. This renders our upper bound on the tensor-to-scalar ratio sensitive to the bounds imposed on the higher slow-roll parameters, cf. Fig. 1. It also implies that future measurements constraining these higher derivatives will be crucial to test this class of ‘intermediate-field’ models.555With a sufficient amount of tuning in the parameters, the impact on the higher derivatives can be suppressed/delayed to the (yet unreported) value of the higher-derivatives of the inflationary potential. In this case, a direct comparison with the observed temperature two-point function might prove more restrictive than the usual expansion in , see e.g.  for a related analysis. Moreover, we note that we do not saturate the upper bound on the tensor-to-scalar ratio recently published in , which might naively allow for for field excursion under the slow-roll condition . The reason, as pointed out also by the authors of Ref. , are again the constraints on the higher-order derivatives ( and beyond) which enforce the transition between the two phases (large and small ) to happen gradually.
Fig. 2 shows this behaviour of the scalar potential for some of the parameter combinations saturating the upper bound on in Fig. 1. The left panel demonstrates how a larger tensor-to-scalar ratio impacts the linear part of the potential, leading to larger field excursion. The right panel shows the decomposition of the total scalar potential in terms of its various contributions according to Eq. (2.8). Generically, the individual supergravity contributions are large, and as expected, it is necessary to tune the coefficients in order to obtain the desired shape of the potential. The globally supersymmetric contribution is responsible for ending inflation but is subdominant for most of the inflationary trajectory. To give an example, the parameters for the potential depicted in the right panel of Fig. 2 are , , , , , , and . This is a quite typical example, in the sense that all dimensionless coefficients are and (as determined by and , cf. Eq. (3.3)) is close to the SUSY GUT scale. Fig. 2 also suggests that, if the coefficient which we switched off had instead an value, the behaviour of the potential would be significantly changed at large . As stated above, our motivation for truncating the potential at is not that higher terms are negligibe (which is the case only for in this example) but that nonzero is the minimal set of nonzero coefficients needed.
What could be the physical reason justifying this particular shape of the scalar potential? The two regimes (characterized by an approximately flat potential at small , and an approximately linear potential at large ) cannot be simply identified with the domains where -suppressed terms are respectively negligible and dominant. This can be seen from the right panel of Fig. 2, which shows that various -suppressed terms contribute significantly to flattening the potential at small , together with the Coleman-Weinberg contribution. Furthermore, the transition between these two domains is generically smooth, while to obtain a large tensor-to-scalar ratio with minimal field excursion, one would prefer them to be linked by a sharper bend. One may speculate that the sudden increase in the potential might be triggered by a phase transition, and that the linear regime at large field values might be an exact feature of the UV theory of quantum gravity (for instance, it might be somehow linked to the linear inflaton potentials in axion monodromy inflation in superstring theory ). In the present work, however, we merely show that tuning the parameters of the scalar potential can mimic such effects.
4 Conclusion and Outlook
In this paper we have revisited models of -term hybrid inflation, which are well known to link the inflationary scale with the scale of supersymmetric grand unification. In its minimal form, FHI also predicts an unobservable small tensor-to-scalar ratio. Motivated by the recent results from BICEP2, we have therefore considered a non-minimal extension of FHI by -suppressed operators.
In fact, if the field excursion of the inflaton is of the order of the Planck scale, the scalar potential is dominated by -suppressed terms which can significantly change the predictions for the CMB observables. Since quantum gravity is expected to break global symmetries, these -suppressed operators generically include non-minimal Kähler terms which explicitly break the -symmetry governing minimal FHI. Treating the coefficients of these terms as free parameters, one can tune them to obtain an inflaton potential which is approximately linear during the first e-folds, and then approximately constant for the remaining e-folds. This allows for a sizeable tensor-to-scalar ratio for sub-Planckian values of the inflaton field. The potential reduces to that of FHI only at very small inflaton values, where -suppressed terms are negligible. One thus retains the connection between the GUT scale and the inflationary scale, and more generally the connection to field-theoretical particle physics model building as in FHI models, while the predictions for the inflationary observables depend mostly on Planck-scale physics.
By carefully choosing the coefficients of the higher-dimensional operators, we can obtain models where the inflaton field excursion saturates the Lyth bound. Since this bound already implies a minimal field excursion of the order of for a tensor-to-scalar ratio of , the Kähler potential must be tuned to a fairly high order in the expansion in order to reproduce the inflationary observables correctly. Nevertheless, we find it interesting that one can write down a model with (slightly) sub-Planckian field values which connects the GUT scale with the scale of inflation, and which allows for a sizeable value of .
At the small field values at the end of inflation the supergravity contributions become negligible. We hence expect the subsequent cosmological processes, i.e. the generation of a thermal bath, of a matter-antimatter asymmetry and of dark matter to proceed as in the globally supersymmetric case, albeit with a larger value for the superpotential coupling as well as for the energy density than is usually assumed, see e.g. [23, 24] for recent analyses. In this parameter range both potential gravitino overproduction and the possible formation of cosmic strings at the end of inflation are potentially dangerous and require careful treatment. The investigation of these model-dependent constraints is however beyond the scope of the present paper.666Cosmic strings with a string tension  would be in serious tension with the Planck data . This might be avoided by reducing the cosmic string tension through a coupling to the MSSM Higgs fields  or by considering GUT groups which, against generic expectation , do not produce topological defects, e.g. flipped SU(5). For avoiding non-thermal gravitino overproduction, cf. e.g. . A further interesting question is the impact of the full two-field dynamics in the complex inflaton plane. Additional fields introduce extra friction, allowing for slow-roll on steeper potentials , and for complicated trajectories the Lyth bound on the total length of the trajectory can still allow for small total field excursions. A extreme example of the latter point was recently given in . However, even the introduction of a single -symmetry violating term in -term hybrid inflation can yield non-trivial trajectories, cf. .
Finally, a crucial task will be to verify the BICEP2 signal in an independent experiment. Considering the tension with the Planck data and uncertainties involving the modelling of the dust foreground, upcoming data from ABS, ACTPol, EBEX, Planck, POLARBEAR, Spider and SPT will hopefully provide a clearer picture. As Fig. 1 demonstrates, a value of the tensor-to-scalar ratio somewhat smaller than the current best-fit value would render this scenario less contrived (but also less testable), while a larger value can hardly be accommodated.
We thank Ido Ben-Dayan, Shaun Hotchkiss, and David Seery for useful discussions. This work has been supported in part by the European Union FP7-ITN INVISIBLES (Marie Curie Action PITAN-GA-2011-289442-INVISIBLES) (VD) and by ERC Advanced Grant 267985 “Electroweak Symmetry Breaking, Flavour and Dark Matter” (FB). VS is supported by the Science and Technology Facilities Council (grant number ST / J000477 / 1). This project was initiated at the March 2014 workshop on the “Implications of the 125 GeV Higgs Boson” at LPSC Grenoble.
Appendix A A model of dynamical hybrid inflation with explicit R breaking
An example for a model of -term hybrid inflation with explicit -symmetry breaking was constructed by Dimopoulos, Dvali and Rattazzi (DDR) in , and will be briefly reviewed and slightly generalized here. We will discuss the limit of rigid supersymmetry first. Following [13, 32] we consider supersymmetric gauge theory with pairs of quarks and antiquarks . This theory is asymptotically free because , and it becomes strongly coupled at a scale . We add a singlet and a superpotential
This theory has a symmetry under which carries -charge 2 and the , are neutral. For the quarks decouple at the scale , below which the theory reduces to super-Yang-Mills theory with beta function coefficient and scale . One-loop matching at the scale yields
Gaugino condensation in the super-Yang-Mills theory thus induces an effective superpotential for , valid at scales , which is exactly the inflaton superpotential of FHI Eq. (2.2):
At large the massive quarks and generate a logarithmic one-loop correction to the scalar potential. Together with the superpotential Eq. (A.3) this results in the inflaton potential of FHI,
where is given by Eq. (2.3) with .
In this model the scale of inflation can be identified with the scale of grand unification when is coupled to GUT-symmetry breaking such that its vacuum expectation value provides a mass to the GUT-breaking field . Taking to be in the of ,777More realistically one should perhaps consider other groups than , since the breaking of to the Standard Model at energies below the scale of inflation may produce magnetic monopoles which cannot be inflated away. Also symmetry breaking patterns which produce cosmic strings at the end of inflation are dangerous if the symmetry breaking scale is indeed as high as indicated by the BICEP2 results. We use here for illustration because our main interest is in the physics of inflation, which largely does not depend on the details of GUT breaking. the tree-level superpotential ,
gives rise to the following effective superpotential at low energies:
Here , , and are composite meson and baryon superfields, and is a Lagrange multiplier enforcing the quantum deformed moduli space constraint . This theory has an isolated supersymmetric vacuum at
where we have defined
Note that the -flatness condition for enforces that the and VEVs are proportional. Also note that the last term in Eq. (A.5) explicitly breaks the U(1) -symmetry. There is now no obvious symmetry reason why and terms in should be absent, but setting their coefficients to negligibly small values is of course technically natural thanks to the non-renormalization theorem.
As long as the gauge coupling is small, the dynamics will not significantly affect the dynamically generated superpotential for large Eq. (A.3). If , is stabilised at zero throughout the inflationary phase and the additional logarithmic contribution from to the Kähler potential Eq. (2.4) will be subdominant.
For studying the effects of -suppressed operators on this model, it needs to be embedded into supergravity. In order to end up in a vacuum with (approximately) vanishing cosmological constant after inflation, we add a constant term,
to the superpotential. Thus the DDR model provides an example for a model in which is not given by the gravitino mass after SUSY breaking (in Planck units), but is instead of the order of .
Appendix B Higher-order contributions to the scalar potential
Note that in the main text we are assuming that the vacuum energy during inflation is dominated by the term from global SUSY, i.e. .
-  P. A. R. Ade et al. [BICEP2 Collaboration], arXiv:1403.3985 [astro-ph.CO].
-  L. Boubekeur and D. .H. Lyth, JCAP 0507, 010 (2005) [hep-ph/0502047].
-  E. J. Copeland, A. R. Liddle, D. H. Lyth, E. D. Stewart and D. Wands, Phys. Rev. D 49 (1994) 6410 [astro-ph/9401011].
-  G. R. Dvali, Q. Shafi and R. K. Schaefer, Phys. Rev. Lett. 73 (1994) 1886 [hep-ph/9406319].
-  See, X. Calmet and V. Sanz, arXiv:1403.5100 [hep-ph], and references therein, for a recent study.
-  D. H. Lyth, Phys. Rev. Lett. 78, 1861 (1997) [hep-ph/9606387].
-  I. Ben-Dayan and R. Brustein, JCAP 1009 (2010) 007 [arXiv:0907.2384 [astro-ph.CO]].
-  S. Hotchkiss, A. Mazumdar and S. Nadathur, JCAP 1202, 008 (2012) [arXiv:1110.5389 [astro-ph.CO]].
-  S. Antusch and D. Nolde, arXiv:1404.1821 [hep-ph].
-  A. Hebecker, S. C. Kraus and A. Westphal, Phys. Rev. D 88, 123506 (2013) [arXiv:1305.1947 [hep-th]].
-  M. Bastero-Gil, S. F. King and Q. Shafi, Phys. Lett. B 651, 345 (2007) [hep-ph/0604198].
-  M. Civiletti, C. Pallis and Q. Shafi, arXiv:1402.6254 [hep-ph].
-  S. Dimopoulos, G. R. Dvali and R. Rattazzi, Phys. Lett. B 410 (1997) 119 [hep-ph/9705348].
-  W. Buchmuller, V. Domcke, K. Kamada and K. Schmitz, arXiv:1404.1832 [hep-ph].
-  P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5082 [astro-ph.CO].
-  D. H. Lyth, arXiv:1403.7323 [hep-ph].
-  A. R. Liddle, P. Parsons and J. D. Barrow, Phys. Rev. D 50, 7222 (1994) [astro-ph/9408015].
-  S. H. Hansen and M. Kunz, Mon. Not. Roy. Astron. Soc. 336, 1007 (2002) [hep-ph/0109252].
-  C. Caprini, S. H. Hansen and M. Kunz, Mon. Not. Roy. Astron. Soc. 339, 212 (2003) [hep-ph/0210095].
-  B. J. Carr, K. Kohri, Y. Sendouda and J. ’i. Yokoyama, Phys. Rev. D 81 (2010) 104019 [arXiv:0912.5297 [astro-ph.CO]].
-  P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5076 [astro-ph.CO].
-  L. McAllister, E. Silverstein and A. Westphal, Phys. Rev. D 82 (2010) 046003 [arXiv:0808.0706 [hep-th]]. See also reviews, D. Baumann and L. McAllister, arXiv:1404.2601 [hep-th]. E. Silverstein, arXiv:1311.2312 [hep-th].
-  K. Nakayama, F. Takahashi and T. T. Yanagida, JCAP 1012, 010 (2010) [arXiv:1007.5152 [hep-ph]].
-  W. Buchmuller, V. Domcke, K. Kamada and K. Schmitz, CERNâProceedingsâ2014â001, (CERN, Geneva, 2014), pp. 47-77 [arXiv:1309.7788 [hep-ph]].
-  M. Hindmarsh, Prog. Theor. Phys. Suppl. 190, 197 (2011) [arXiv:1106.0391 [astro-ph.CO]].
-  P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5085 [astro-ph.CO].
-  R. Jeannerot, J. Rocher and M. Sakellariadou, Phys. Rev. D 68, 103514 (2003) [hep-ph/0308134].
-  M. Hindmarsh and D. R. T. Jones, Phys. Rev. D 87, 075022 (2013) [arXiv:1203.6838 [hep-ph]].
-  K. Nakayama, F. Takahashi and T. T. Yanagida, Phys. Lett. B 718, 526 (2012) [arXiv:1209.2583 [hep-ph]].
-  A. R. Liddle, A. Mazumdar and F. E. Schunck, Phys. Rev. D 58, 061301 (1998) [astro-ph/9804177].
-  J. McDonald, arXiv:1404.4620 [hep-ph].
-  S. Dimopoulos, G. R. Dvali, R. Rattazzi and G. F. Giudice, Nucl. Phys. B 510 (1998) 12 [hep-ph/9705307].
-  N. Seiberg, Phys. Rev. D 49 (1994) 6857 [hep-th/9402044].