Inflation from Supersymmetry Breaking

We explore the possibility that inflation is driven by supersymmetry breaking with the superpartner of the goldstino (sgoldstino) playing the role of the inflaton. Moreover, we impose an R-symmetry that allows to satisfy easily the slow-roll conditions, avoiding the so-called $\eta$-problem, and leads to two different classes of small field inflation models; they are characterised by an inflationary plateau around the maximum of the scalar potential, where R-symmetry is either restored or spontaneously broken, with the inflaton rolling down to a minimum describing the present phase of our Universe. To avoid the Goldstone boson and remain with a single (real) scalar field (the inflaton), R-symmetry is gauged with the corresponding gauge boson becoming massive. This framework generalises a model studied recently by the present authors, with the inflaton identified by the string dilaton and R-symmetry together with supersymmetry restored at weak coupling, at infinity of the dilaton potential. The presence of the D-term allows a tuning of the vacuum energy at the minimum. The proposed models agree with cosmological observations and predict a tensor-to-scalar ratio of primordial perturbations $10^{-9}\lesssim r\lesssim 10^{-4}$ and an inflation scale $10^{10}$ GeV $\lesssim H_*\lesssim 10^{12}$ GeV. $H_*$ may be lowered up to electroweak energies only at the expense of fine-tuning the scalar potential.


Introduction
Inflationary models [1] in supergravity 5 suffer in general from several problems, such as fine-tuning to satisfy the slow-roll conditions, large field initial conditions that break the validity of the effective field theory, and stabilisation of the (pseudo) scalar companion of the inflaton arising from the fact that bosonic components of superfields are always even. The simplest argument to see the fine tuning of the potential is that a canonically normalised kinetic term of a complex scalar field X corresponds to a quadratic Kähler potential K = XX that brings one unit contribution to the slow-roll parameter η = V ′′ /V , arising from the e K proportionality factor in the expression of the scalar potential V . This problem can be avoided in models with no-scale structure where cancellations arise naturally due to non-canonical kinetic terms leading to potentials with flat directions (at the classical level). However, such models require often trans-Planckian initial conditions that invalidate the effective supergravity description during inflation. A concrete example where all these problems appear is the Starobinsky model of inflation [3], despite its phenomenological success.
In this work we show that all three problems above are solved when the inflaton is identified with the scalar component of the goldstino superfield 6 , in the presence of a gauged R-symmetry. Indeed, the superpotential is in that case linear and the big contribution to η described above cancels exactly. Since inflation arises in a plateau around the maximum of the scalar potential (hill-top) no large field initial conditions are needed, while the pseudo-scalar companion of the inflaton is absorbed into the Rgauge field that becomes massive, leading the inflaton as a single scalar field present in the spectrum. This model provides therefore a minimal realisation of natural smallfield inflation in supergravity, compatible with present observations, as we show below. Moreover, it allows the presence of a realistic minimum describing our present Universe with an infinitesimal positive vacuum energy arising due to a cancellation between an Fand D-term contributions to the scalar potential, without affecting the properties of the inflationary plateau, along the lines of Ref. [5,6,7].
On general grounds, there are two classes of such models depending on whether the maximum corresponds to a point of unbroken (case 1) or broken (case 2) R-symmetry. The latter corresponds actually to a generalisation of the model we studied recently [6], inspired by string theory [5]. It has the same field content but in a different field basis with a chiral multiplet S ∝ ln X playing the role of the string dilaton. Thus, S has a shift symmetry which is actually an R-symmetry gauged by a vector multiplet and the superpotential is a single exponential. The scalar potential has a minimum with a tuneable vacuum energy and a maximum that can produce inflation when appropriate corrections are included in the Kähler potential. In these coordinates R-symmetry is restored at infinity, corresponding to the weak coupling limit. As we show below, this model can be generalised to a class of models (case 2) where inflation arises at a plateau where (gauged) R-symmetry is spontaneously broken. Small field inflation is again guaranteed consistently with the validity of the effective field theory. In this work, we are mostly focused on a new possibility (case 1) where the maximum is around the origin of X where R-symmetry is restored.
The outline of the paper is the following: In Section 2, we describe the model and the two cases to study. In Section 3, we analyse in detail case 1. We first work out model independent predictions valid under the assumption that inflation ends very rapidly after the inflationary plateau of the scalar potential (section 3.1). We then discuss the possibility for the existence of a minimum nearby the maximum with an infinitesimal tuneable positive vacuum energy and show that additional corrections to the Kähler potential are needed (section 3.2). The notion of nearby is defined in the sense that perturbative expansion around the maximum is valid for the Kähler potential, but not for the slow-roll parameters. We subsequently present an explicit example and work out the cosmological predictions (section 3.3). Finally, in section 3.4, we discuss the anomaly cancellation associated to the U(1) R gauge symmetry. In Section 4 we analyse case 2, firstly in general (section 4.1) and then in a particular example. In Section 5 we derive a stringent constraint on low energy inflation models in a model-independent way. Section 5 can be read independently of the rest of the paper. Our conclusions are presented in Section 6, while Appendix A contains a proof of an identity used in Section 5.

Symmetric versus non-symmetric point
In this work we are interested in supergravity theories containing a single chiral multiplet transforming under a gauged R-symmetry with a corresponding abelian vector multiplet. We assume that the chiral multiplet X (with scalar component X) transforms as: where q is its charge, and ω is the gauge parameter. The Kähler potential is therefore a function of XX, while the superpotential is constrained to be of the form X b : where X is a dimensionless field and κ −1 = m p = 2.4 × 10 15 TeV is the (reduced) Planck mass. For b = 0, the gauge symmetry eq. (1) becomes a gauged R-symmetry. The gauge kinetic function can have a constant contribution as well as a contribution proportional to ln X f (X) = γ + β ln X.
The latter contribution proportional to β is not gauge invariant and can be used as a Green-Schwarz counter term to cancel possible anomalies. We will show however below, in Section 3.4, that the constant β is fixed to be very small by anomaly cancellation conditions and does not change our results. We will therefore omit this term in our analysis below.
We are interested in the general properties of supergravity theories of inflation that are of the above form. Before performing our analysis, a distinction should be made concerning the initial point where slow-roll inflation starts. The inflaton field (which will turn out to be ρ, where X = ρe iθ ) can either have its initial value close to the symmetric point where X = 0, or at a generic point X = 0. The minimum of the potential, however, is always at a nonzero point X = 0. This is because at X = 0 the negative contribution to the scalar potential vanishes and no cancellation between F-term and D-term is possible. The supersymmetry breaking scale is therefore related to the cosmological constant as κ −2 m 2 3/2 ≈ Λ. One could in principle assume that the value of the potential at its minimum is of the supersymmetry breaking scale. However, in this case additional corrections are needed to bring down the minimum of the potential to the present value of the cosmological constant, and we therefore do not discuss this possibility.
In the first case, inflation starts near X = 0, and the inflaton field will roll towards a minimum of the potential at X = 0. On the other hand, in the second case inflation will start at a generic point X = 0. In order to make easier contact with previous literature [5,6], it is convenient to work with another chiral superfield S, which is invariant under a shift symmetry by performing a field redefinition In this case the most general Kähler potential and superpotential are of the form Note that this field redefinition is not valid at the symmetric point X = 0 for the first case.
The first case will be discussed in Section 3, and the second case will be discussed in Section 4.
3 Case 1: Inflation near the symmetric point 3

.1 Slow roll parameters
In this section we derive the conditions that lead to slow-roll inflation scenarios, where the start of inflation is near a local maximum of the potential at X = 0. Since the superpotential has charge 2 under R-symmetry, one has W = 0 as long as R-symmetry is preserved. Therefore, W can be regarded as the order parameter of R-symmetry breaking. On the other hand, the minimum of the potential requires W = 0 and broken R-symmetry. It is therefore attractive to assume that at earlier times R-symmetry was a good symmetry, switching off dangerous corrections to the potential. As similar approach was followed in [8], where a discrete R-symmetry is assumed. Instead, we assume a gauged R-symmetry which is spontaneously broken at the minimum of the potential.
While the superpotential is uniquely fixed in eq. (2), the Kähler potential is only fixed to be of the form K(XX). We expand the Kähler potential as follows where A and f are constants. The gauge kinetic function is taken to be constant since it will be shown in Section 3.4 that the coefficient β in front of the logarithmic term in eq. (3) is fixed to be very small by anomaly cancellation conditions. As far as the scalar potential is concerned, the coefficient γ can be absorbed in other parameters of the theory. We therefore take γ = 1. The scalar potential is given by 7 where and The superpotential is not gauge invariant under the U(1) gauge symmetry. Instead it transforms as 7 We follow the conventions of [9]. Therefore, the U(1) is a gauged R-symmetry which we will further denote as U(1) R . From W X k X R = −r R κ 2 W , where k X R = −iqX is the Killing vector for the field X under the R-symmetry, r R = iκ −2 ξ R with κ −2 ξ R the Fayet-Iliopoulos contribution to the scalar potential, and W X is short-hand for ∂W/∂X, we find A consequence of the gauged R-symmetry is that the superpotential coupling b enters the D-term contribution of the scalar potential as a constant Fayet-Iliopoulos contribution. 8 Note that the scalar potential is only a function of the modulus of X and that the potential contains a Fayet-Iliopoulos contribution for b = 0. Moreover, its phase will be 'eaten' by the U(1) gauge boson upon a field redefinition of the gauge potential similarly to the standard Higgs mechanism. After performing a change of field variables the scalar potential is a function of ρ, Since we assume that inflation starts near ρ = 0, we require that the potential eq. (14) has a local maximum at this point. It turns out that the potential only allows for a local maximum at ρ = 0 when b = 1. For b < 1 the potential diverges when ρ goes to zero. For 1 < b < 1.5 the first derivative of the potential diverges, while for b = 1.5, one has V ′ (0) = 9 4 f 2 + 3 2 q 2 > 0, and for b > 1.5, on has V ′′ (0) > 0. We thus take b = 1 and the scalar potential reduces to Note that in this case the the superpotential is linear W = f X, describing the sgoldstino (up to an additional low-energy constraint) [12]. Indeed, modulo a D-term contribution, the inflaton in this model is the superpartner of the goldstino. In fact, for q = 0 the inflaton reduces to the partner of the goldstino as in Minimal Inflation models [13]. The important difference however is that this is a microscopic realisation of the identification of the inflaton with the sgoldstino, and that the so-called η-problem is avoided (see discussion below). The kinetic terms for the scalars can be written as 9 8 For other studies of inflation involving Fayet-Iliopoulos terms see for example [10], or [11] for more recent work. Moreover, our motivations have some overlap with [8], where inflation is also assumed to start near an R-symmetric point at X = 0. However, this work uses a discrete R-symmetry which does not lead to Fayet-Iliopoulos terms. 9 The covariant derivative is defined as∂ µ X = ∂ µ X − A µ k X R , where k X R = −iqX is the Killing vector for the U(1) transformation eq. (1).
It was already anticipated above that the phase θ plays the role of the longitudinal component of the gauge field A µ , which acquires a mass by a Brout-Englert-Higgs mechanism.
We now interpret the field ρ as the inflaton. It is important to emphasise that, in contrast with usual supersymmetric theories of inflation where one necessarily has two scalar degrees of freedom resulting in multifield inflation [14], our class of models contains only one scalar field ρ as the inflaton. In order to calculate the slow-roll parameters, one needs to work with the canonically normalised field χ satisfying The slow-roll parameters are given in terms of the canonical field χ by Since we assume inflation to start near ρ = 0, we expand where we defined q = f x. Notice that for ρ ≪ 1 the ǫ parameter is very small, while the η parameter can be made small by carefully tuning the parameter A. Any higher order corrections to the Kähler potential do not contribute to the leading contributions in the expansion near ρ = 0 for η and ǫ. Such corrections can therefore be used to alter the potential near its minimum, at some point X = 0 without influencing the slow-roll parameters.
A comment on the η-problem in Supergravity A few words are now in order concerning the η-problem [15]. The η problem in N = 1 supergravity is often stated as follows (see for example [16]): If, for instance, a theory with a single chiral multiplet with scalar component ϕ is taken, then the Kähler potential can be expanded around a reference location ϕ = 0 as where φ is the canonically normalised field φφ = K ϕφ (0)ϕφ, and the ellipses stand for extra terms in the expansion coming from K and W . Following this argument, the mass m φ turns out to be proportional to the Hubble scale and therefore Or otherwise stated, this leading contribution of order 1 to the η-parameter has its origin from the fact that the F-term contribution to the scalar potential contains an exponential factor e K : However, in our model the factor '1' drops out for the particular choice b = 1 in the superpotential 10 , resulting in an inflaton mass m 2 ρ which is determined by the next term A(XX) 2 in the expansion of the Kähler potential, As a result, there are two ways to evade the η-problem: • First, one can obtain a small η by having a small q ≪ f , while A should be of order O(10 −1 ). In this case, the rôle of the gauge symmetry is merely to constrain the form of the Kähler potential and the superpotential, and to provide a Higgs mechanism that eliminates the extra scalar (phase) degree of freedom.
• Alternatively there could be a cancellation between q 2 and 4Af 2 .
Since A is the second term in the expansion of the Kähler potential eq. (7), it is natural to be of order O(10 −1 ) and therefore providing a solution to the η-problem. This will be demonstrated in an example in Section 3.3. Note that the mass of the inflaton given in eqs. (23) is only valid during inflation at small ρ. The mass of the inflaton at its VEV (Vacuum Expectation Value) will be affected by additional corrections that are needed to obtain in particular a vanishing value for the scalar potential at its minimum, which will be discussed in the following sections.

The upper bound on the tensor-to-scalar ratio
Before moving on to the next section, let us focus on the approximation at ρ ≪ 1 where the perturbative expansion of the slow-roll parameters in eqs. (19) is valid, and assume that the horizon exit occurs at the field value ρ * very close to the maximum ρ = 0. In this approximation, eqs. (19) become where the asterisk refers to the value of parameters evaluated at the horizon exit. 10 Note that in hybrid inflation models the η-problem is also evaded by a somewhat similar way, but these models generally include several scalar fields (and superfields) besides the inflaton (see e.g. [17]).
To discuss the upper bound on the tensor-to-scalar ratio, it is convenient to divide the region [ρ = 0, ρ end ] into two regions: one is [0, ρ p ], where the approximation (24) is valid, and the other is the rest [ρ p , ρ end ]. Here ρ end means the inflation end. Note that ρ p < ρ end because the approximation (24) breaks down before the end of inflation where ǫ(ρ end ) = 1 or |η(ρ end )| = 1. In terms of this division, the number of e-folds from the horizon exit to the end of inflation can be approximated by where we introduced Here χ is the canonically normalised field defined by eq. (17). Let us next focus on the region [ρ p , ρ end ]. It is natural to expect the following inequality This is based on the following observation. The right hand side describes a hypothetical situation, as if the slow-roll condition were valid throughout the inflation until its end. But since in the actual inflation the slow-roll condition breaks down in the region [ρ p , ρ end ], the actual number of e-folds in this region will be smaller than that in the hypothetical situation. Adding N pert (ρ * , ρ p ) to the both hand sides of (27) and using (25), we find Using (24) and the definition of the tensor-to-scalar ratio r = 16ǫ * , we obtain the upper bound: To satisfy CMB data, let us choose η = −0.02 and N CMB ≈ 50. Assuming ρ end 1/2, we obtain the upper bound r 10 −4 . Note that this is a little bit lower than the Lyth bound [18] for small field inflation, r 10 −3 . From the upper bound on r, we can also find the upper bound on the Hubble parameter as follows. In general, the power spectrum amplitude A s is related to the Hubble parameter at horizon exit H * by Combining this with the upper bound r 10 −4 and the value A s = 2.2 × 10 −9 by CMB data, we find the upper bound on the Hubble parameter H * 10 9 TeV. In Section 5, we will also find the lower bound r 10 −9 (equivalently H * 10 7 TeV), based on an model-independent argument. This bound can be lowered at the cost of naturalness between parameters in the potential.

de Sitter vacua
In the previous section we showed that for these models the slow-roll parameters can be in eqs. (19). We showed above that it is indeed easy to satisfy this condition, providing a potential solution to the η-problem. In this set-up however, we are ignorant of what happens after the end of inflation. In the following, we will extend the above analysis to include the minimum of the potential, and we require after inflation the scalar field ρ to roll into a 'nearby' minimum of the potential with a tunably small but positive cosmological constant. We show below that, for both cases q = 0 and q = 0, that the slow-roll conditions do not allow for the existence of such a minimum. In other words, while it is easy to safisfy the slow-roll conditions at the maximum of the potential, additional corrections are needed to ensure a vanishing (or tunably small and positive) cosmological constant. These corrections terms modify the potential near the minimum, while the leading contributions to the slow-roll parameters, given by eqs. (19), only depend on the first two terms in the expansion of the Kähler potential K = XX + A(XX) 2 + . . . . Moreover, as is often the case, having a vanishing cosmological constant comes at the cost of fine-tuning parameters in the model. An example of such a correction is proposed in Section 3.3 and compared with the most recent CMB results.
Let us here explain more precisely our definition for a minimum 'nearby' the maximum. In a phenomenologically realistic model, inflation should end before the minimum of the scalar potential. We mentioned in previous sections that the perturbative expansion of the slow-roll parameters in eq. (19) is not valid at the end of inflation. So, it must not be valid at the minimum either. We should then define the existence of a nearby minimum around the maximum of the potential in the 'weaker' sense in which the minimum is in a region where the perturbative expansion of the Kähler potential is valid but not of the slow-roll parameters.
The need for additional corrections for q = 0 In this section we assume q = 0 (the case q = 0 is treated separately below), and we show that a model defined by eqs. (7) does not allow for a tunably small cosmological constant at a nearby minimum when |η| ≪ 1 at the maximum is required.
The scalar potential V = V F given by eq. (15) with q = 0, and is repeated here for convenience The solutions to V(ρ 0 ) = 0 and V ′ (ρ 0 ) = 0 give A plot of this potential is shown in Figure 1. The inflaton starts near the local maximum at ρ = 0 and rolls towards the minimum at ρ = 0.91082. Near the local maximum of the potential at ρ = 0, η is given by eq. (19) while the parameter A is already fixed by requiring a Minkowski vacuum in eq. (33). As a result, we have We conclude that the slow-roll condition |η| ≪ 1 is not consistent with the existence of a tunably small cosmological constant at a nearby minimum for q = 0.
The need for additional corrections for q = 0 The scalar potential is given by eq. (15), and is repeated here for convenience The equation V ′ (ρ) = 0 is solved by The other solutions to V ′ (ρ) = 0 are ρ = 0, which is assumed to be a maximum near the start of inflation, and 1 + ρ 2 + 2Aρ 4 = 0 which corresponds to an AdS (Anti-de Sitter) minimum (if A < 0). We therefore focus on the solution eq. (36) with A > 0 . The condition for a vanishing value of the scalar potential at the minimum V(ρ) = 0 can be combined with eq. (36) to yield For A > 0, this equation can be solved to give ρ = ρ 0 , the value of ρ at the minimum, while V(ρ 0 ) = 0 is satisfied when the relation between parameters f and q is given by where As an example, for A = 1 2 , eq. (37) gives ρ 0 = 0.872008 (the other solutions are not physical). Taking for example q = 1, gives f ≈ 2.928 by eq. (38). The result is plotted in Figure 2. We now show that this is inconsistent with |η| ≪ 1 at the maximum. Eq. (19) can be solved to give This can be combined with eq. (38) to give an expression for η in terms of A and ρ 0 One can now calculate η for every A > 0 (We do not consider A < 0 since this would lead to an AdS vacuum). Numerical analysis however shows that |η| > 1 for all A > 0, and a small η-parameter is impossible in these models when a vanishing cosmological constant is imposed. We conclude that additional corrections to the Kähler potential are necessary in order to obtain a tunably small cosmological constant consistent with |η| ≪ 1. These corrections can be higher order or non-perturbative corrections.

Correction terms that allow for a tunable minimum
We propose corrections to the Kähler potential of the form while the superpotential is fixed by the gauge symmetry with b = 1 (the arguments that excluded b = 1 in Section 3.1 still apply). The resulting scalar potential is given by where F X = ∂ X F and F XX = ∂ X ∂X F . 11 For example, if we choose the Kähler potential can be expanded near X = 0 as where we made a field redefinition X ′ = √ 1 + ξX to bring the Kähler potential into the form of the expansion eqs. (7). 12 It is important to emphasise that the higher order corrections to the Kähler potential do not contribute to the leading order of slow-roll parameters. One can therefore apply eqs. (19) where we used that ρ ′ 2 = (1 + ξ)ρ 2 , and The scalar potential for the parameters is plotted in Figure 3. By using ρ * = 0.04493 and ρ end = 0.32844 one obtains the results in Table 1. The slow-roll parameters ǫ and η during the inflation are shown in Figure 4.  It is important to emphasise that neither the choice of function F in eq. (44), nor the choice for the particular values of the parameters is unique. Different choices of F , or even different choices for the parameters α, ξ and B can lead to similar CMB results as the ones presented in Table 1. Figure 5 shows that our predictions for n s and r are within 2σ C.L. of Planck '15 contours with the number of e-folds N = 57.98.
On the other hand, the gravitino mass is given by m 3/2 = 5.96 × 10 7 TeV and the Hubble scale is 6.80×10 8 TeV. No other parameter set consistent with the CMB data were found with a Hubble scale below 10 8 TeV. Moreover, the tensor-to-scalar ratio satisfies r > 10 −6 . We believe that the inability to find models with gravitino mass of order 10 TeV 11 This can in principle be written in terms of ρ by using F X = 1 2 F ρ e −iθ and F XX = 1 4 F ρρ + 1 4ρ F ρ . 12 Note that we also have W = 1 1+ξ f X ′ , although this is not important for the discussion below.
is independent of the choice of F . We will put stringent conditions on the existence of models with a TeV scale Hubble parameter in Section 5. This will be done in a way not leaning on any ingredients of supergravity, and these results are therefore more general than the scope of this paper. However, in supergravity the gravitino mass is usually of the same order as the Hubble scale, which heavily constraints the supersymmetry breaking scale.
In particular, in Section 5 we show that a necessary, but not sufficient condition to obtain order 10 TeV scale inflation is given in eq. (82), repeated here for convenience of the reader where the derivatives are with respect to the canonically normalised field χ satisfying eq. (17). This can for example be realised by including an extra term in the Kähler potential proportional to Z(XX) 3 . The above constraint then requires Z ≫ 1, which violates our assumption that the Kähler potential can be expanded around the maximum of the potential at the symmetric point X = 0. As a result, the expansions for ǫ and η in eqs. (19) require correction terms of higher order in ρ, and the above analysis is not valid. Moreover, the above condition (50) forces Z to be several orders of magnitude larger than the other parameters, which raises questions about naturalness of such a model. The authors confirm that the inclusion of the term Z(XX) 3 indeed leads to potentials that allow a Hubble scale of order 10 TeV for large Z. The discussion of such models however is postponed for future work, since the inclusion of a very large parameter is not natural, and such a term would reintroduce higher order terms in ρ in the above analysis.
A final comment is in order: In contrast with Case 2 in Section 4, where the D-term and the F-term contributions to the scalar potential are of the same order such that a vanishing cosmological constant can be found by a careful cancellation between the Fterm and the D-term, the model above allows for a vanishing cosmological constant even when q 2 ≪ f 2 . A vanishing cosmological constant can be obtained by carefully tuning the parameters governing the F-term contribution to the scalar potential. However, we emphasise that even in this scenario the gauge symmetry still plays an important rôle: It constrains the form of the Kähler potential and the superpotential, while the U(1) gauge boson 'eats' the phase of X resulting in an inflation scenario with a single scalar degree of freedom. 13

Anomaly cancellation
In this section we show that the coefficient β in eq. (3) is indeed small and can be neglected. First, as was already noted in eq. (11), the superpotential is not gauge invariant under the  Table 1 for a scalar potential in eq. (43) with correction term eq. (44), in the n s − r Planck '15 results for TT,TE,EE,+lowP and assuming ΛCDM + r.
U(1) R gauge symmetry, and therefore is an R-symmetry. A consequence of U(1) R is that scalars and fermions within the same multiplet carry different charges under the symmetry. In particular, the total U(1) R -charges of the chiral fermion χ, the U(1) R -gaugino and the gravitino are given by This results in a contribution to the cubic U(1) 3 R anomaly proportional to C R = Tr[R 3 ] = R 3 χ + R 3 λ + 3R 2 3/2 14 , given by [22] This can however be canceled by a Green-Schwarz mechanism, since the field-dependent term in gauge kinetic function eq. (3) results in a contribution to the Lagrangian of the form Under a gauge transformation, this results in a contribution Anomaly cancelation requires δL 1−loop + δL GS = 0, (55) 14 Note that the contribution of the gravitino is three times that of a gaugino [21].

Contributions from MSSM fermions and mixed anomalies can be treated similarly.
If we instead assume that no extra matter fields are present, we have C R = −5 q 2 3 , and β R = 5 96π 2 q 2 . It follows that β R is very small if q is small. For example, in the previous section we had q = 2.91 ×10 −8 and ρ * = 0.04493, which gives |β R ln ρ * | ≈ 1.8 ×10 −17 ≪ 1, which justifies the approximation f (X) = 1 + β R ln X ≈ 1.
(56) 4 Case 2: Inflation away from the symmetric point In this section we consider the case where inflation starts at a point far away from the symmetric point X = 0. This allows us to make a field redefinition eq. (5), and work with a chiral superfield S (with scalar component s) invariant under a shift symmetry. In this case, the most general superpotential is given by W = κ −3 ae bs , while the Kähler potential should be a function of s +s. The gauge kinetic function is at most linear in s i.e. f (s) = γ + βs.
Since we are focused on small-field inflation, it is natural to assume that the inflaton starts near a (local) maximum of the potential 15 . However, in contrast with Case 1, the maximum of the potential does not correspond to a point of restored R-symmetry and there is no particular reason why the inflation started precisely near this point. This poses a fine-tuning problem that is present in many small-field models of inflation, and we shall not address this further. We will show that, as in the case 1, the minimum of the inflation potential lies near the local maximum in the weaker sense.

Behaviour near the maximum of the potential
Let us start by considering a perturbation around the local maximum of the scalar potential φ 0 , or equivalently making the change of variable s = δs + φ 0 /2. We expand the Kähler potential around the maximum up to the 4th order of the classical fluctuation δs (and δs) K(δs, δs) = κ −2 (Ã +B(δs + δs) +C(δs + δs) 2 +D(δs + δs) 3 +Ẽ(δs + δs) 4 ). (57) We will show below that higher-order terms in the expansion of the Kähler potential do not contribute to the "leading order" terms in slow-roll parameters. Indeed, it was already emphasised in the previous section that only terms up to (XX) 2 are needed in the Kähler potential. We can also express the superpotential and gauge kinetic function as f (δs) = γ + β(δs + φ 0 /2). (59) By a Kähler transformation to absorbÃ followed by the redefinitions the Kähler potential and superpotential can be simplified into K(δs, δs) = κ −2 (C(δs + δs) 2 +D(δs + δs) 3 +Ẽ(δs + δs) 4 ), where we dropped the primes on a, b, γ. The D-term contribution to the scalar potential is given by and the F-term contribution is Let us define φ = s +s such that δφ = δs + δs represents the (classical) fluctuation around the maximum of the potential at φ 0 . Then the potential can be written in terms of δφ as V = 1 2κ 4 a 2 eC δφ 2 +Dδφ 3 +Ẽδφ 4 +bδφ (b + 2Cδφ + 3Dδφ 2 + 4Ẽδφ 3 ) 2 C + 3Dδφ + 6Ẽδφ 2 − 6 We fix γ = 1 and β = 0 for simplicity. The slow-roll parameters ǫ and η are defined in terms of the canonically normalised fluctuation δχ, which is defined by Since we expand the potential around the maximum of the potential, the slow-roll parameter must satisfy ǫ(δφ = 0) = 0. This gives us the following constraint on the parameters, Using the above constraint, we can write the expansion of the slow-roll parameters near the maximum of the potential as and Note that η parameter can be small and negative by carefully fine-tuning four parameters a, b,C andẼ. It is also important to note that the slow-roll expansions are valid only near the maximum of the potential i.e. δφ ≪ 1. They may break down during inflation or at the minimum of the potential. These expansions are useful for showing qualitatively that the η-problem can be avoided. In order to compare any predictions with the CMB data, one needs to use the full expression for ǫ and η.

Example
In order to give an example of this class of models, let us consider the case γ = 1, β = 0, with Kähler potential [6] This model is obviously invariant under the shift symmetry in eq. (4). With φ = s +s = δφ + φ 0 , the potential in terms of fluctuation δφ around the maximum of the potential is given by In this example, the canonically normalised field χ is defined by For this case, we need contributions from both F-term and D-term in order to obtain Minkowski vacua. It was shown in [6,22] that this model is anomaly-free. By choosing, a = 2.34422 × 10 −6 , b = −0.0234, c = 7.10 × 10 −6 , ξ = 0.3023, α = −0.7813, we obtain an appropriate inflationary potential with a flat plateau around the maximum φ 0 ≈ 66.37 illustrated in Figure 6. The horizon exit is at φ * = 64.5315 and inflation ends at φ end = 50.9915. By using eq. (18), the value of the slow-roll parameters at horizon exit are ǫ(φ * ) ≃ 1.30 × 10 −7 and η(φ * ) ≃ −2.01 × 10 −2 . The number of e-folds N, the scalar power spectrum amplitude A s , the spectral index of curvature perturbation n s and the tensor-to-scalar ratio r are calculated and summarised in Table 2, in agreement with Planck '15 data. Figure 7 shows that our predictions for n s and r are within 1σ C.L. of Planck '15 contours with the number of e-folds N ≈ 56.82. During inflation, we can show that a field fluctuation from the maximum of the potential in terms of the canonically normalised field is quite small in Planck units, i.e. κ|δχ * | = κ|χ * − χ 0 | ≈ 0.036 and κ|δχ end | = κ|χ end − χ 0 | ≈ 0.351, but at the minimum of the potential φ min ≈ 21.50 the field fluctuation is large, i.e. κ|δχ min | = κ|χ min −χ 0 | ∼ 1.81. However, one can show that the perturbative expansion of the Kähler potential is still valid and the minimum is nearby the maximum in the weaker sense. Indeed, using the values of parameters a, b, c, ξ and α given above, we can expand the Kähler potential (69) in the fluctuation δφ around the maximum of the scalar potential as It is easy to prove that this expansion makes sense as a perturbative expansion even for |δφ min | = |φ min − φ 0 | ∼ 44.9 (equivalent to κ|δχ min | ∼ 1.81 as given above). Thus, in this example the minimum lies nearby the maximum in the sense defined above. Note however that in contrast with the claim in [6], the potential given in eq. (70) is inconsistent with a supersymmetry breaking scale (and Hubble scale) in the multi-TeV range. It has turned out to be very difficult to find any parameters, consistent with the CMB data, for which the Hubble scale H * is lower than 10 8 TeV. The reason is outlined   Table 2: The theoretical predictions for φ * = 64.5316 and φ end = 50.9915 and the parameters given in Figure 6.
in Section 5: Eq. (76), defining the amplitude A s , fixes ǫ * in terms of H * (for A s fixed by the CMB data), resulting in a very small (∼ 10 −23 ) value for ǫ * . This in turn results in a big number of e-folds occurring near the horizon exit, which rapidly exceeds the required amount of e-folds N CMB between the horizon exit and the end of inflation, which should be of order ∼ 40. Indeed, the potential in eq. (70) does not satisfy the necessary (but not sufficient) condition to obtain order 10 TeV scale inflation 1 κ 3 V ′′′ * V * ≫ 10 6 in eq. (82).

A stringent constraint on TeV Hubble scale
In the previous sections, we have seen that it is difficult to realise inflation models with a TeV scale gravitino mass. Here we characterise this difficulty in a more general way, and independent of the supergravity framework. In particular, we show that a necessary but not sufficient condition to have order 10 TeV inflation is that 1 The CMB observation covers the inflation between the horizon exit and the end of inflation. We first estimate the number of e-folds during this epoch using the following formula [23,24,25,26]: where q 0 , T 0 are the wave number in the physical spatial coordinate and the CMB temperature at present, respectively. Note that this formula is independent of inflation models.
Here we take T 0 /q 0 ∼ 10 28 with q 0 ≃ 0.002 Mpc −1 following [25, pp.325-326]. We assume that H end ≃ H * , and that during the reheating era the energy density evolves as ρ ∝ a −y with 3 ≤ y < 4. The asterisk refers to the horizon exit as in the previous sections. The reheating temperature T reh is assumed to be of the same order as H end . Under the above rough assumptions we find that, for H * ranging from 1 TeV to 10 10 TeV, the number of e-folds N CMB is roughly between 40 and 60. Let us next discuss the condition that an inflaton potential should yield N CMB = 40 -60 for a given Hubble parameter H * . The amplitude A s and tilt n s are given by where we used that the Hubble scale is given by H * = κ V * /3. We choose A s = 2.2 × 10 −9 and n s = 0.96 to satisfy the CMB data. Putting them into eqs. (75) and (76) yields the slow-roll parameters ǫ * and η * as functions of the Hubble constant H * , given by (77) The plot of these expressions is given in Figure 8. The Lyth bound [18] for small field Figure 8: A plot of ǫ * and η * versus the Hubble parameter H * at the horizon exit. The blue line refers to ǫ * and the orange one refers to η * . Note that ǫ * is much smaller than η * for H * 10 10 TeV, which fits our models given in the previous sections.
inflation, r 10 −3 , implies H * 10 10 TeV by eq. (77). From Figure 8 it is clear that if the Hubble scale satisfies H * 10 10 TeV, then the slow-roll parameters at horizon exit satisfy ǫ * ≪ |η * |. 16 In the following we consider the Hubble scale in the range 1 TeV ∼ 10 10 TeV. 17 16 An upper bound on the tensor-to-scalar ratio r 10 −4 from Section 3 implies H * 10 9 TeV. This is consistent with H * 10 10 TeV by the Lyth bound. 17 Note that the maximum Hubble scale H * ∼ 10 10 TeV corresponds to the GUT scale vacuum energy The inflaton potential decreases monotonically through the end of inflation χ end towards the global minimum. For simplicity and without loss of generality, we assume that the field value of the inflaton decreases monotonically during inflation. We can express the number of e-folds between the horizon exit χ * and some field value χ in terms of the inflaton potential as In terms of eq. (78), the number of e-folds between the horizon exit and the end of inflation is given by N(χ end ). Any scalar potential should therefore satisfy N(χ end ) ≈ N CMB in order to be a good candidate for slow-roll inflation satisfying the CMB data. Let us divide the region [χ end , χ * ] into two parts, region A and B. Region A is in the neighbourhood of the horizon exit χ * and the slow-roll parameters in region A are required to satisfy the conditions and the region B is its complement in [χ end , χ * ]. Note that the condition (79) is natural in the neighbourhood of the horizon exit as long as the Hubble parameter satisfies H * 10 10 TeV, as was already observed in Figure 8. The inflaton first passes through region A and enters region B when the conditions (79) are violated. This transition occurs well before the end of inflation where the slow-roll parameters are of order 1. It is shown in the Appendix that if χ is in region A, then the number of e-folds N(χ) can be approximated as where ∆χ = |χ * − χ|. Using this expression, we can show that, for any χ in region A (and therefore satisfy eqs. (79)), the following inequality holds: This inequality can be used to exclude a given potential in the following way: If for any particular inflaton potential there exists some χ in region A for which N A (χ) > N CMB , then it contradicts the inequality (81), and the inflaton potential can be excluded. In other words, if the flat region near the start of inflation is too big, too many e-folds will occur in region A. As a result, the number of e-folds in this region exceeds the required number of e-folds N CMB . This allows us to put a limit on the scalar potential by investigating the size of region A, characterised by κ∆χ.
From eq. (80) we find that, for the Hubble scale H * ∼ 10 p TeV with 0 ≤ p ≤ 10, the choice κ∆χ ≃ 10 p−9 results in a number of e-folds in region A given by N A (χ) ≃ 80, which is sufficiently larger than N CMB ≃ 40 -60. It remains to check whether this χ is still in region A and therefore satisfies eqs. (79). It can be shown 18 that χ with κ∆χ ≃ 10 p−9 lies in region A if the inflaton potential satisfies In summary, we have shown that if eq. (82) is satisfied (for H * = 10 p TeV with 0 ≤ p ≤ 10), then there exists χ in region A which gives N(χ) = 80, which exceeds N CMB , and therefore contradicts the inequality (81). Therefore, the condition in eq. (82) is a sufficient condition to exclude an inflaton potential.
In other words, the condition is a necessary (but not sufficient) condition for any the inflaton potential such that the number of e-folds does not exceed N CMB .
A rough estimate on the necessary condition implies that if H * ≥ 10 7 TeV, a milder condition on the potential 1 κ 3 V ′′′ * V * < 1 should be satisfied. It should therefore be much easier to construct potentials with a Hubble scale above 10 7 TeV. Note that this implies a theoretical (rough) lower limit on the tensor-to-scalar ratio, namely r > 10 −9 (or equivalently, ǫ > 10 −10 ). On the other hand, if the Hubble scale is low, around the TeV scale, then the necessary condition puts a very stringent constraint on the potential which poses questions about the naturalness of the theory since a large (∼ 10 9 ) hierarchy between parameters in the model is necessary.
It is important to emphasise that the above argument holds for all models of (slow-roll) inflation, and is not limited to supergravity models. In supergravity models however, the Hubble scale H * is usually of the same order as the gravitino mass m 3/2 . 19 The above argument therefore puts stringent constraints on the lower limit of the supersymmetry breaking scale. It states that a theory with low energy supersymmetry breaking (m 3/2 of order O(10 TeV)) is very unlikely to simultaneously be a theory for inflation, unless the third derivative of the scalar potential attains very high values near the start of inflation. 18 Combining κ∆χ = 10 p−9 with the assumptions (77) and (82), we find that at the horizon exit. As a result, the condition η(χ) ≈ η(χ * ) satisfied. Similarly, by Taylor expanding √ ǫ in ∆χ, we can show that ǫ(χ) ≪ |η(χ)|. 19 An exception to this occurs in so-called New Inflation models [27], where the Hubble scale is typically much larger than the gravitino mass and a gravitino mass as low as 100 TeV is allowed [28].

Conclusions
In this work we proposed a class of supergravity models of inflation in which the inflaton is the superpartner of the goldstino, carrying a unit charge under a gauged R-symmetry. Inflation occurs around the maximum of the scalar potential where R-symmetry is either restored (case 1) or spontaneously broken (case 2), with the inflaton rolling down towards the electroweak minimum of the supersymmetric Standard Model. Analyticity and Rinvariance imply a linear superpotential which automatically avoids the so-called etaproblem in supergravity, allowing for slow-roll and small field inflation consistently with an effective field theory description. Case 2 is a generalisation of a particular model we studied in a previous publication [6], inspired by type I string theory with moduli stabilisation by internal magnetic fluxes and the inflation identified with the string dilaton. Case 1 offers a new possibility, with the R-symmetry restored around the inflationary plateau, which we analysed in detail in this work: • Assuming first that inflation ends very rapidly after the plateau where small fluctuations are valid around the maximum of the scalar potential, we were able to find simple analytic formulae for the slow-roll parameters and the number of efolds, leading to an upper bound for the tensor-to-scalar ratio r < ∼ 10 −4 and for the inflation scale H * < ∼ 10 12 GeV.
• We then show that additional corrections to the Kähler potential are needed in order to describe a minimum of the scalar potential 'nearby' the maximum with an infinitesimal tuneable positive vacuum energy. The 'nearby' notion is defined in the sense that perturbative expansion around the maximum is valid for the Kähler potential but not for the slow-roll parameters. We subsequently study an example of such corrections and work out the experimental predictions for the cosmological observables.
We finally derived a stringent and model-independent constraint on TeV scale inflation (or in general low scale). The constraint implies that in order to attain a theory with Hubble scale in the multi-TeV range, one has to tune the form of the potential, implying in general a naturalness problem between its parameters.
Note first that the identity can be approximated between χ and χ * by This can be solved analytically by The number of e-folds at χ then becomes This completes the proof of eq. (80).