Inflation from Supergravity with Gauged R-symmetry in de Sitter Vacuum

We study the cosmology of a recent model of supersymmetry breaking, in the presence of a tuneable positive cosmological constant, based on a gauged shift symmetry of a string modulus that can be identified with the string dilaton. The minimal spectrum of the `hidden' supersymmetry breaking sector consists then of a vector multiplet that gauges the shift symmetry of the dilaton multiplet and when coupled to the MSSM leads to a distinct low energy phenomenology depending on one parameter. Here we study the question if this model can also lead to inflation by identifying the dilaton with the inflaton. We find that this is possible if the K\"ahler potential is modified by a term that has the form of NS5-brane instantons, leading to an appropriate inflationary plateau around the maximum of the scalar potential, depending on two extra parameters. This model is consistent with present cosmological observations without modifying the low energy particle phenomenology associated to the minimum of the scalar potential.


Introduction
A fundamental theory of Nature, such as string theory, should be able to describe at the same time particle physics and cosmology, which are phenomena that involve very different scales from the microscopic four-dimensional (4d) quantum gravity length of 10 −33 cm to large macroscopic distances of the size of the observable Universe ∼10 28 cm, spanned a region of about 60 orders of magnitude. In particular, besides a e-mail: antoniadis@itp.unibe.ch b e-mail: dma3ac2@gmail.com c e-mail: hiroshi.isono81@gmail.com d e-mail: rob.knoops@cern.ch the 4d Planck mass, there are three very different scales with very different physics corresponding to the electroweak, dark energy and inflation. These scales might be related via the scale of the underlying fundamental theory, such as string theory, or they might be independent in the sense that their origin could be based on different and independent dynamics.
In this work, we make an attempt in this direction by connecting the scale of inflation with the electroweak and supersymmetry breaking scales within the same effective field theory, which at the same time allows the existence of an infinitesimally small (tuneable) positive cosmological constant describing the present dark energy of the Universe. To this end, we use a simple model of supersymmetry breaking in a tuneable metastable de Sitter vacuum, independent of the scale of supersymmetry breaking that can be in the TeV region [1][2][3][4]. The model is based on a shift symmetry of a string modulus S (which we identify with the string dilaton) along its imaginary (axionic) component, which is gauged by a U (1) vector multiplet. The latter can be for instance a linear combination of the ordinary baryon and lepton number, containing the matter parity which guarantees a dark matter candidate [5].
The superpotential W is completely fixed to a simple exponential (W = ae bS , in a Kähler basis where the shift symmetry is an R-symmetry) depending on two parameters, while a third parameter arises from the U (1) gauge coupling. On the other hand, the Kähler potential K is an arbitrary function of S +S. However, since one is interested in a vacuum where the string coupling is weak (and therefore S large), in order to study the supersymmetry phenomenology, one can restrict K to its tree-level logarithmic form, K = −p log(S +S) with p = 1 or 2, depending on whether S is associated to only the D9 or to the D9 and D5 gauge couplings in type I string theory [6]. This model has the necessary ingredients to be obtained as a remnant of moduli stabilisation within the framework of internal magnetic fluxes in type I string theory, turned on along the compact directions for several abelian factors of the gauge group. All geometric moduli can in principle be fixed in a supersymmetric way, while the shift symmetry is associated to the 4d axion and its gauging is a consequence of anomaly cancellation [7][8][9].
The resulting scalar potential has a metastable minimum with a tuneable positive vacuum energy that can be made infinitesimally small and one is left with one free parameter which fixes the scale of supersymmetry breaking. This is due to a tuning between the D-and the F-term contributions to supersymmetry breaking that can have opposite signs in supergravity. Coupling this model to the observable sector (MSSM) is straightforward but requires the introduction of an additional parameter in order to address the problem of anomaly cancellation or of tachyonic scalar masses for U (1) neutral matter fields [3,4].
The main question we address in this work is whether the same scalar potential can provide inflation with the dilaton playing also the role of the inflaton at an earlier stage of the evolution of the Universe. We show that this is possible if one modifies the Kähler potential by a correction that plays no role around the minimum, but creates an appropriate plateau around the maximum. In general, the Kähler potential receives perturbative and non-perturbative corrections that vanish in the weak coupling limit. After analysing all such corrections, we find that only those that have the form of (Neveu-Schwarz) NS5-brane instantons can lead to an inflationary period compatible with cosmological observations. The scale of inflation turns out then to be of the order of low energy supersymmetry breaking, in the TeV region. On the other hand, the predicted tensor-to-scalar ratio is too small to be observed.
An interesting property of this model is that the inflaton is a component of the goldstino superpartner and shares some of the features of the models proposed in Refs. [10,11], 1 although it does not belong to the same class, since there is a D-term contribution to supersymmetry breaking. As a result, the gravitino mass does not satisfy the property to be much lower than the supersymmetry breaking scale.
The outline of the paper is the following. In Sect. 2, we give a brief review of the model and compute the slow-roll parameters of the scalar potential to show that without modifying the Kähler potential, it does not give rise to inflation (Sect. 2.1). We then present its general extension, parametrising appropriately the quantum corrections to the Kähler potential (Sect. 2.2). In Sect. 3, we perform the analysis for a correction that has the form of NS5-brane instantons (depending on two additional parameters) for both p = 2 (Sect. 3.1) and p = 1 (Sect. 3.2) cases. In particular, we compute the corresponding slow-roll parameters and fix the two parameters of the correction, so that our model reproduces the spectral index and amplitude of density fluctuations in agreement with the cosmological data of Planck '15. We also extract the predictions for the inflation scale, the number of e-foldings and the tensor-to-scalar ratio. In Sect. 4, we study the impact of the correction to the Kähler potential for the low energy superparticle spectrum that we find to be of order of 10%. Finally, Sect. 5 contains our concluding remarks. For certain values of the parameters, the model allows for a tuneably small and positive value for the cosmological constant. Its anomaly cancellation conditions are discussed in [3], while the resulting low energy spectrum is discussed in [4,5]. In this section we recall the main properties of the model. We use the conventions of [15].
The model consists of one chiral multiplet S, whose scalar component s is invariant under a gauged shift symmetry 2 where θ is the gauge parameter, and c is a constant. The Kähler potential, superpotential and gauge kinetic function are given by (in appropriate Kähler coordinates where the superpotential is constant) where κ −1 = m p = 2.4 × 10 15 TeV is the inverse of the (reduced) Planck mass. The scalar potential is given by where the Greek indices α, β label the chiral multiplets in the theory, and the capital Roman letters A, B label the different gauge groups. In Eq. (3), the Kähler covariant derivative of the superpotential is and the moment maps P A are given by where k α A are the Killing vectors, and r A is the Fayet-Iliopoulos contribution satisfying W α k α A = −κ 2 r A W . In the above model k s = −ic is the Killing vector associated with the shift symmetry Eq. (1). As a result, the scalar potential is given by (with φ = s +s) For b > 0, the potential always admits a supersymmetric AdS (anti-de Sitter) vacuum at φ = b/ p, while for b = 0 supersymmetry is broken in AdS space. We therefore focus on b < 0. In this case, the potential admits a supersymmetry breaking dS (de Sitter) vacuum for p < 3. For example, for p = 2, β = 1 and γ = 0, 3 the scalar potential Eq. (6) reduces to A vanishing cosmological constant can be found by tuning the parameters of the model. Solving the equations V(φ min ) = 0 and dV(φ min )/dφ = 0 gives where l 0 is the root of the polynomial −x 5 + 7x 4 − 10x 3 − 22x 2 + 40x + 8 close to −0.18, and In Eq. (9), A(l 0 ) is given by A non-zero cosmological constant can be found if The gravitino mass parameter is given by Equation (8) fixes the VEV (Vacuum Expectation Value) of φ as a function of the parameter b. The field φ can be interpreted as the dilaton 4 related to the string coupling constant g s by φ min = 2/g s . A string coupling constant in the perturbative region, for example φ min = 10, can therefore be obtained for example by the choice b = −0.01820. Since the parameters a and c are related by Eq. (9), this leaves only one free parameter, which can be tuned to obtain an O(10 TeV) gravitino mass. For example, the parameter choice c = 0.61 × 10 −13 results in m 3/2 = 12.83 TeV. However, we will show below that this model does not allow slow-roll inflation to occur. The kinetic terms in the Lagrangian for the scalar φ are given by The canonically normalised field χ therefore satisfies χ = κ −1 p 2 log φ. The slow-roll parameters are given by It can be shown that, when the conditions (8) and (9) are satisfied, the slow-roll parameters depend only on ρ = −bφ, , and a plot of the slow-roll parameter η as a function of ρ (right). The slow-roll condition |η| 1 is not satisfied for any value of the parameters a, b, c Moreover, the scalar potential is given in terms of ρ where A(l 0 ) ≈ −50.66 as in Eq. (9). In Fig. 1, a plot is shown of κ 4 V(ρ) |b| 3 c 2 as a function of ρ. The minimum of the potential is at ρ min ≈ 0.1832 (see Eq. (8)), while the potential has a local maximum at ρ max ≈ 0.4551. A plot of the slow-roll parameter η (also in Fig. 1) shows that |η| 1 is not satisfied. This result holds for any parameters a, b, c satisfying Eqs. (8) and (9).
A similar analysis to the one above can be performed for p = 1. In contrast with the case p = 2, vacua with a vanishing cosmological constant can also be found when β = 0 and γ = 0. Note, however, that for β = 0 the Lagrangian contains a Green-Schwarz term which is not gauge invariant. In principle, such a term can be present if the model contains additional fields that are charged under the U (1) gauge symmetry (see Eq. (1)), in order to cancel their contribution to a possible cubic U (1) 3 anomaly via a Green-Schwarz mechanism. However, it turns out (see for example [5]) that even in this case β is very small. We therefore take β = 0 for simplicity. 5 Moreover, as long 5 The Green-Schwarz term Eq. (17) is always present in the case p = 2 above, since in this case one can only find vacua with a vanishing as the scalar potential is concerned, γ can be absorbed in other parameters of the model, and we therefore take γ = 1.
For γ = 1 and β = 0, vacua with a vanishing cosmological constant can be found if bφ min = l 0 ≈ −0.233153, and a 2 bc 2 = A(l 0 ) ≈ −2.783259. However, a similar analysis as the one above shows that also in this case the slow-roll condition η 1 cannot be satisfied.

Extensions of the model that satisfy the slow-roll conditions
In the previous section we showed that the slow-roll conditions can not be satisfied in the minimal versions of the model. In this section we modify the above model by modifying the Kähler potential. While the superpotential is uniquely fixed (up to a Kähler transformation 6 ), the Kähler potential admits corrections that can always be put in the form while the superpotential, the gauge kinetic function and moment map are given by Footnote 5 continued cosmological constant when β > 0. This is the reason why we focus on p = 1 below. 6 For example, by performing the Kähler transformation with J (z) = −κ −2 bs, the linear contribution in the Kähler potential can be absorbed into the superpotential, which becomes of the form W (s) = κ −3 a exp (bs).
Identifying Re(s) with the inverse string coupling, the function F may contain perturbative contributions that can be expressed as power series of 1/(s +s), as well as nonperturbative corrections which are exponentially suppressed in the weak coupling limit. The later can be either of the form e −δ(s+s) for δ > 0 in the case of D-brane instantons, or of the form e −δ(s+s) 2 in the case of (Neveu-Schwarz) NS5-brane instantons (since the closed string coupling is the square of the open string coupling). We have considered a generic contribution of these three different types of corrections and we found that only the last type of contributions can lead to an inflationary plateau providing sufficient inflation. The other corrections can be present but do not modify the main properties of the model (as long as weak coupling description holds). In the following section, we analyse in detail a function F describing a generic NS5-brane instanton correction to the Kähler potential.

p = 2 case
We now consider the case with where b < 0 and α < 0. F(φ) vanishes asymptotically at large φ. In this case, we obtain and There are four parameters in this model namely α, ξ , b and c. The first two parameters α and ξ control the shape of the potential. There are some regions in the parameter space of α and ξ that the potential satisfies the slow-roll conditions i.e. 1 and |η| 1. In order to obtain the potential with flat plateau shape which is suitable for inflation and in agreement with Planck '15 data, we choose 7 α = −4.841115384560439 and ξ = 0.025350051999999998. (24) Note that in the case of ξ = 0 and b < 0, we can find the Minkowski minimum by solving the equations V(φ min ) = 0 and dV(φ min )/dφ = 0, where φ min = s min +s min is the value of φ at the minimum of the potential. In the case of ξ = 0, we cannot solve the equations analytically and the relations (8), (9) are not valid. We can always assume that they are modified into bφ min = l(ξ, α) and a 2 bc 2 = −50.66016761885055 × λ(ξ, α, ) 2 , where λ takes positive values and satisfies |λ − 1| 1. For any given value of parameters ξ , α and the cosmological constant , one can numerically fix the value of l and λ. By fine-tuning the cosmological constant to be very close to zero, we can numerically solve the equations V = 0 and dV/dφ = 0 for the value of l and λ in (25) as From Eq. (25), we can see that the third parameter, b, controls the vacuum expectation value φ min . This can be shown in Fig. 2 where we compare the scalar potential for different values of b. Motivated by string theory, we have the identification φ ∼ 1/g s . We can choose the value of the parameter b such that φ min is of the order of 10 to make sure that we are in the perturbative regime in g s . The last parameter, c, controls the overall scale of the potential but does not change its minimum and its shape. This can be shown in Fig. 3. In the following, we will fix b and c by using the cosmological data.
In order to compare the predictions of our models with Planck '15 data, we choose the following boundary conditions: The initial conditions are chosen very near the maximum on the (left) side, so that the field rolls down towards the electroweak minimum. Any initial condition on the right of the maximum may produce also inflation, but the field will roll towards the SUSY vacuum at infinity. The results are therefore very sensitive to the initial conditions (Eqs. (28) and (29)) of the inflaton field. The slow-roll parameters are given in Eq. (14). The total number of e-folds N can be determined by Note that we choose |η(χ end )| = 1. We can compare the theoretical predictions of our model to the experimental results via the power spectrum of scalar perturbations of the CMB, namely the amplitude A s and tilt n s , and the relative strength of tensor perturbations, i.e. the tensor-to-scalar ratio r . In terms of slow-roll parameters, these are given by where all parameters are evaluated at the field value χ int . In order to satisfy the Planck '15 data, we choose the parameters b = −0.0182, c = 0.61 × 10 −13 . The slow-roll parameters and η during the inflation are shown in Fig. 4. The value of the slow-roll parameters at the beginning of inflation are (φ int ) 1.86 × 10 −24 and η(φ int ) −1.74 × 10 −2 . (34) The total 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 1, in agreement with Planck '15 data [18].   shows that our predictions for n s and r are within 1σ C.L. of Planck '15 contours with the total number of e-folds N ≈ 1075. Note that N is the total number of e-folds from φ int to φ end . However, the number of e-folds associated with the CMB observation corresponds to a period between the time of horizon crossing and the end of inflation, which is much smaller than 1075. According to the general formula in [18], the number of e-folds between the horizon crossing and the end of inflation is roughly estimated to be around 50-60.
We would like to remark that the parameter c also controls the gravitino mass at the minimum of the potential around O(10) TeV. Indeed, the gravitino mass is written as For b = −0.0182, we get φ min ≈ 9.91134 and for the gravitino mass at the minimum of the potential The Hubble parameter during inflation (evaluated at φ * = φ int ) is This shows that our predicted scale for inflation is of the order of TeV. The mass of gravitino during the inflation m * 3/2 = 4.15 TeV is higher than the inflation scale, and the gauge boson mass is M * A μ = 3.12 TeV. 8 In fact, the gauge boson acquires a mass due to a Stueckelberg mechanism by eating the imaginary component of s, where its mass at the minimum of the potential is given by As a result, the model essentially contains only one scalar field Re(s), which is the inflaton. This is in contrast with other supersymmetric models of inflation, which usually contain at least two real scalars [23]. 9 8 The gauge boson mass is given by m Aμ = 2g ss c 2 /Re(s). 9 This is because a chiral multiplet contains a complex scalar.

p = 1 case
In this case, we obtain and The potential has similar properties to the p = 2 case, although it may give different phenomenological results at low energy. For ξ = 0 and b < 0, the Minkowski minimum satisfies the following relations [2]: However, similar to the previous case, the above relations are not valid when ξ = 0 and we assume that they are modified into bφ min = l(ξ, α) and bc 2 a 2 = −0.35929085159984514 × λ(ξ, α, ) −2 .
By choosing α = −0.781 and ξ = 0.3023 and tuning the cosmological constant to be very close to zero, we can numerically fix l = −0.562536 and λ = 1.2937645 for this case. The gravitino mass for p = 1 case can be written as By choosing the parameters b = −0.0234, c = 1 × 10 −13 , the gravitino mass at the minimum of the potential is  As summarised in Table 2, the predictions for the p = 1 case are similar to those of p = 2, in agreement with Planck '15 data with the total number of e-folds N ≈ 888. In this case, the Hubble parameter during inflation is Note that, for the p = 1 case, the mass of the gauge boson is M * A μ = 6.78 TeV, and the mass of the gravitino during inflation is m * 3/2 = 4.72 TeV, a bit smaller than but comparable with the Hubble scale. It would be interesting to investigate how much this would affect the power spectrum and bispectrum, inspired by a class of inflation models called quasisingle field inflation [19][20][21][22][23][24][25][26][27], which may contain also a light inflaton as well as massive fields of Hubble scale mass, which can amplify non-Gaussianities. Moreover, the U (1) mass is close to the Hubble scale (although slightly above it) so that it may produce observable effects in the E-mode polarisation spectrum.

SUGRA spectrum
The above model can be coupled to MSSM-like fields ϕ. In this case the multiplet containing the inflaton Re(s) is considered to be a "hidden sector" field which is responsible for breaking supersymmetry by both F-and D-terms, as described above. The supersymmetry breaking is then communicated to the visible sector (MSSM) via gravity mediation. We consider the following Kähler potential and superpotential: where K(s +s) is given by Eq. (18), the hidden sector superpotential W h (s) is given by Eq. (19), and W MSSM is the MSSM superpotential, which only depends on the MSSM fields ϕ. The soft supersymmetry breaking terms can be calculated as follows: Here, m 2 0 is the scalar soft mass squared. All trilinear couplings are the same and given by A 0ŷi , whereŷ i are the Yukawa couplings of the rescaled MSSM superpotential W MSSM = e κ 2 K/2 W MSSM . The Bμ-term parameter is given by B 0μ , whereμ = e K/2 μ.
For p = 2 the Lagrangian contains a Green-Schwarz term Eq. (17), and the theory is not gauge invariant (without the inclusion of extra fields that are charged under the U (1)). We therefore focus on p = 1. The soft terms can be written in terms of the gravitino mass (see Eq. (35)) where Using the parameters presented in Sect. 3.2, we find m 3/2 = 18.36 TeV and C = 1.53. For ξ = 0 the model reduces to the one analysed in [4], where one has C = 1.52 and m 3/2 = 17.27 TeV (with φ min = 9.96). As in [4], the scalar soft mass is tachyonic. This can be solved either by introducing an extra Polonyi-like field, or by allowing a non-canonical Kähler potential for the MSSM-like fields ϕ. The resulting low energy spectrum is expected to be similar to the one described in [4]. We do not perform this analysis, but we only summarise their results.
Since the tree-level contribution to the gaugino masses vanishes, their mass is generated at one-loop by the so-called 'Anomaly Mediation' contribution [28][29][30]. As a result, the spectrum consists of very light neutralinos (O(10 2 ) GeV), of which the lightest (a mostly Bino-like neutralino) is the LSP dark matter candidate, slightly heavier charginos and a gluino in the 1-3 TeV range. The squarks are of the order of the gravitino mass (∼10 TeV), with the exception of the stop squark, which can be as light as 2 TeV.

Conclusions
In this paper, we have studied the cosmology of a simple model of supersymmetry breaking, based on a single chiral multiplet (the string dilaton) with a gauged shift symmetry (R-symmetry in an appropriate field basis) leading to a de Sitter vacuum with a tuneable cosmological constant [1,2]. When coupled to the MSSM, the model have been shown to yield an interesting low energy pattern of supersymmetry breaking, depending on one parameter, distinct from other proposals [4,5]. By modifying the dilaton Kähler potential with a non-perturbative correction that could arise from NS5brane instantons (depending on two extra parameters), we show that the scalar potential can acquire an inflationary plateau producing sufficient inflation, consistent with cosmological observations, without altering the low energy particle spectrum around the minimum.
The predicted tensor-to-scalar ratio is unfortunately too small to be detected. On the other hand, our model may lead to measurable non-gaussianities and/or E-mode polarisation effects that deserve further study.
Our model provides therefore an interesting example of connecting the inflation sector with the 'hidden' sector of supersymmetry breaking by identifying the inflaton with the dilaton which is also the scalar partner of the goldstino (sgoldstino). It does not belong, however, to the class of models studied in the past [10,11], since there is also a D-term component in the order parameter of supersymmetry breaking. It would be interesting to analyse in detail the generalised class of such models. Another interesting question is the possible implementation of this model in string theory. As mentioned in the introduction, the framework of moduli stabilisation by internal magnetic fields for several abelian factors of the gauge group along the compact directions of the compactified manifold, combined with non-criticality, seems to be a promising direction to explore.