Hybrid inflation from supersymmetry breaking

We extend a recently proposed framework, dubbed inflation by supersymmetry breaking, to hybrid inflation by introducing a waterfall field that allows to decouple the supersymmetry breaking scale in the observable sector from the inflation scale, while keeping intact the inflation sector and its successful predictions: naturally small slow-roll parameters, small field initial conditions and absence of the pseudo-scalar companion of the inflaton, in terms of one free parameter which is the first order correction to the inflaton K\"ahler potential. During inflation, supersymmetry is spontaneously broken with the inflaton being the superpartner of the goldstino, together with a massive vector that gauges the R-symmetry. Inflation arises around the maximum of the scalar potential at the origin where R-symmetry is unbroken. Moreover, a nearby minimum with tuneable vacuum energy can be accommodated by introducing a second order correction to the K\"ahler potential. The inflaton sector can also play the role of the supersymmetry breaking 'hidden' sector when coupled to the (supersymmetric) Standard Model, predicting a superheavy superparticle spectrum near the inflation scale. Here we show that the introduction of a waterfall field provides a natural way to end inflation and allows for a scale separation between supersymmetry breaking and inflation. Moreover, the study of the global vacuum describing low energy Standard Model physics can be done in a perturbative way within a region of the parameter space of the model.


Introduction
In past works [1][2][3], a framework of intimate connection between supersymmetry and inflation was introduced, dubbed inflation by supersymmetry breaking, relating two theoretical proposals motivated, correspondingly, by particle physics and cosmology.The basic idea is to identify the inflaton with the superpartner of the goldstino, 1 charged under a gauged R-symmetry with the corresponding U (1) field becoming massive by absorbing the (complex) inflaton phase which is the R-goldstone boson.The superpotential is forced to be linear in the inflaton superfield X by symmetry, while the Kähler potential can be expanded around its canonical form in powers of X X with small coefficients representing quantum corrections of the underlying supergravity theory.For a positive sign of the first order correction, neglecting D-term contributions suppressed by a small R-gauge coupling, the scalar potential has a maximum that can realise hilltop inflation with a spectral index of primordial scalar perturbations determined by this correction.Moreover, a second order correction allows for accommodating a nearby minimum with tuneable vacuum energy.The (supersymmetric) Standard Model (SM) can be coupled in a straightforward way with the inflaton part playing also the role of the supersymmetry breaking sector, while the gauge R-symmetry may contain the usual R-parity as a subgroup.The model is very predictive but leads to a superheavy spectrum of superparticles near the inflation scale.
In this work we address the problem of scale separation, aiming to decouple the scale of inflation from the one of low energy supersymmetry breaking but staying always within the successful framework of inflation by supersymmetry breaking.To this end, we generalise the above framework to hybrid inflation [16,17], by introducing a waterfall direction in the scalar potential that opens up at a critical point at the end of inflation along which the waterfall field falls rapidly to a deep global minimum while the inflaton is displaced slightly.The waterfall field is neutral under the R-symmetry and its superpotential can be adjusted so that the global minimum stays infinitesimally close to a supersymmetric one which consists of a valley around a flat direction for the inflaton field.This allows for a perturbative treatment of the whole dynamics around the origin for the inflaton field and around the supersymmetric minimum for the waterfall with the supersymmetry breaking scale being a free parameter independent of the inflation scale.Moreover, without affecting the inflationary predictions of the model, the presence of the waterfall direction provides a way to end inflation efficiently and tune the vacuum energy of the global minimum at a value infinitesimally close to zero.The cancellation occurs between the negative F-term and positive D-term contributions to the scalar potential. 2he outline of the paper is the following.In Section 2, we present the model we study in this work containing the inflation (dubbed hidden) and the observable SM sector.The inflation sector contains the inflaton chiral superfield charged under an abelian vector field that gauges the U (1) R symmetry, as well as the waterfall chiral superfield which is neutral under U (1) R .We give the Kähler potential and superpotential that determine the N = 1 supergravity effective action, assuming a constant (field independent) R-gauge kinetic function.We describe the setup of inflation and recall its predictions.In Section 3, we analyse the vacuum structure of the model, first without and then with the D-term potential contribution.The analysis is perturbative around the origin of the inflaton direction and to the leading order in the supersymmetry breaking scale which can be parametrically small compared to the inflation scale.In Section 4, we impose theoretical and observational constraints and determine the allowed parameter region.In Section 5, we discuss the supersymmetry breaking and the particle spectrum in both hidden and observable sectors of the theory.We end with our conclusions in Section 6.Finally, there are two appendices with technical details on the computation of the parameter space and of fermion masses.

Inflaton sector
We work with four-dimensional N = 1 supergravity theories which contain two chiral multiplets and an abelian vector multiplet that is associated with a gauged U (1) R transformation.We denote the two chiral superfields by X, ϕ, and suppose that X transforms as X → e −iqθ X but ϕ is neutral under the U (1) R .Our model is defined by the following Kähler potential, superpotential and gauge kinetic function, where κ −1 = 2.4 × 10 15 TeV is the reduced Planck mass.The superfields X, ϕ have mass dimension 1, while the parameters A, B, z, γ, f, µ, λ, β R are dimensionless.In the rest of this paper, we set κ = 1.The parameters A, B, z, γ, f, µ are (or can be chosen to be) real numbers, while λ -although complex in general -is assumed to be also real for simplicity.The superpotential has U (1) R charge q.The logarithmic correction in h is introduced for the anomaly cancellation involving the gauged U (1) R (Green-Schwarz mechanism) [22][23][24].The coefficient β R is proportional to q 2 .We will neglect this logarithmic term as the charge squared q will be supposed to be tiny, except when we discuss the U (1) R gaugino mass which is determined by the logarithmic correction.We parametrise the scalar fields X, ϕ, which are the lowest components of the superfields X, ϕ, respectively, as where ρ, υ, ϕ R , φ are real fields.When ρ ̸ = 0, the phase υ is absorbed to make the U (1) R gauge field massive.We will identify ρ with the inflaton.The scalar potential is the sum of the F -and D-term potentials V = V F + V D , where each one is given by,3 (2.5) The potential is invariant under ϕ → φ plus X → X because it is real.It is also invariant under X → X because it is independent of the phase of X.Therefore, ϕ → φ, equivalent to φ → −φ, is a symmetry of the potential.Furthermore, the potential is invariant under ϕ → −ϕ alone because this keeps K and W invariant and the expression (2.5) contains even numbers of ϕ and φ.Combining them yields the invariance under ϕ R → −ϕ R .It is therefore enough to consider the region ϕ R ≥ 0 and φ ≥ 0.

Coupling with supersymmetric Standard Model sector
The model given above can be coupled with the supersymmetric Standard Model (MSSM).The inflaton sector then plays a role of supersymmetry breaking sector.Following [25], we consider two types of the MSSM superpotential: one in which the MSSM superfields are all neutral under U (1) R , and the other in which some of the MSSM superfields are charged under U (1) R in such a manner Figure 1: Inflation followed by waterfall.The figure in the bubble magnifies the detailed shape of the potential during the inflation, which is a shallow valley.φ ± s and φ 0 will be introduced explicitly in Section 3.1.The real part ϕ R is suppressed because the potential is stable in this direction and the trajectory always runs with ϕ R = 0. that the MSSM superpotential has U (1) R charge q.Concretely, we will consider the following two models [25]: where ζ collectively denotes the MSSM matter chiral superfields H u , H d , ū, d, Q, ē, L, and the sum is over ζ.The vacuum expectation value of ζ is zero.In Model 1, the MSSM superfields are neutral under U (1) R .In Model 2, H u and H d are neutral while the other MSSM superfields have R-charge q/2, so that Ŵ2 has R-charge q.Under this charge assignment, the U (1) R group contains the usual R-parity as its subgroup.

Scenario
Our scenario is inflation followed by waterfall (see Figure 1): First, the inflaton is located near the origin (ρ, ϕ) = (0, 0), and then rolls down slowly.During this slow-roll, the potential is stable in ϕ, forming a valley along the ρ-axis.In the meantime the inflaton reaches a "critical" point ρ c , where the potential becomes unstable in some direction in ϕ.The inflaton then deviates from the ρ axis, falls in this unstable direction (waterfall field ϕ is turned on), and reaches a vacuum at (ρ m , ϕ m ).We assume that the inflaton still moves in the ρ direction after ρ c and hence ρ c and ρ m satisfy where the last inequality guarantees perturbative treatment.
Let us analyse the structure near the origin ρ = ϕ = 0.The mass term for ϕ at each ρ is (2.10) The mass functions at the origin ρ = 0 are given by (2.12) Our scenario requires that the potential is stable in ϕ R , φ and unstable in ρ around ρ = 0.For this, we impose that µ > 0 and m I (0) 2 > 0, that is, In accordance with this condition, we introduce a parameter ε by The mass function m 2 I at each non-vanishing ρ behaves as

.15)
As ρ increases from zero, this is positive for a while, but becomes zero at For consistency, we need to require ρ c < 1 and thus ϵ < 1 too.After ρ c , φ becomes tachyonic and the inflaton will likely fall into this direction ("waterfall") due to quantum fluctuations.This requires a dedicated careful investigation that goes beyond the scope of our paper.

Inflation
We assume that inflation ends by ρ c , so that it is of single-field, hilltop type and described by the potential at ϕ = 0, around the origin of ρ.This is in the framework of inflation by supersymmetry breaking [1][2][3]25], as already described in the Introduction.We summarise some results that will be used later.
The inflaton potential V inf = V (ρ, ϕ = 0) is given by (2.17) The slow-roll parameters then read Note that the leading order of η V can be tuned at will, while ϵ V is much smaller, ϵ V ≪ O(η 2 V ).This is an important feature of models in the framework of inflation by supersymmetry breaking [1].
The amplitude of scalar curvature fluctuations A s , tilt n s , and tensor-to-scalar ratio r for the cosmic microwave background (CMB) are given by where ρ * is ρ at the horizon exit.Combining them, we can express the model parameters ρ * , A, f in terms of the CMB parameters n s , r, A s as (2.21) In the limit where U (1) R D-term contributions are neglected when the corresponding gauge coupling is very small (as will be justified in Section 3.2.2),q can be set to zero, leading essentially to one parameter A fixed by the spectral index n s , while the inflation scale f is constrained by the observational data according to (2.21) for q = 0.The horizon exit ρ * and the end of inflation ρ end are related via the number of e-folds ∆N by As mentioned above, ρ end is required to satisfy ρ end ≤ ρ c .In particular, when the equality holds, ρ c is determined by √ r through the expression for ρ * in (2.21).

Vacuum structure
The next task is to find a non-supersymmetric vacuum at (ρ m , ϕ m ).We require that the vacuum energy V m = V (ρ m , ϕ m ) can be tuned to zero.As we have seen above, the inflaton falls in the direction of the imaginary part φ, keeping ϕ R = 0. We may therefore restrict our analysis at ϕ R = 0. Recall that inflation with the subsequent waterfall is supposed to occur for small ρ.We can then deal with the potential perturbatively in ρ.The strategy of our analysis is that we first find the extrema of the potential at ρ = 0 as the unperturbed part, and then study how they evolve when ρ is turned on.In the following analysis, we shall use a new parameter δ instead of λ, defined by4 and also consider it to be small.Indeed, as we will see below, for δ = 0 one obtains a supersymmetric vacuum with a flat direction along ρ; thus, the parameter δ is expected to control the scale of supersymmetry breaking.We will therefore analyse the potential perturbatively both in ρ and in δ.

Potential at ρ = 0
Let us first consider the potential at ρ = 0. Since the D-term potential at ρ = 0 is just a constant in ϕ, it is sufficient to analyse the F-term potential.The F-term potential at ρ = 0 reads Its extrema with zero or small ϕ R are given by They are given explicitly by Notice that they are all pure imaginary.The potential is singular for |ϕ| = z −1/2 .Its denominator is the Kähler metric K X X for ρ = 0, which is positive if |ϕ| < z −1/2 .The extrema with |ϕ| < z −1/2 differ depending on whether the parameters µ, z satisfy µ > 2z or µ < 2z, which are given by: where ε is defined in (2.14).Whether µ > 2z or µ < 2z, the potential is stable in ϕ R .On the other hand, the stability around each extremum of (3.6) in φ is summarised in the following table: The potential for each case is depicted schematically in Figure 2.
Case µ > 2z: The extremum ϕ = 0 is a local minimum in both ϕ R and φ, as already stated in Section 2.1.Next, ϕ − 1 is a saddle point near the origin.Indeed, ϕ − 1 is of order ε, Figure 2: Structure of the F-term potential with ρ = 0 is given for the cases µ > 2z and µ < 2z.The potential is stable (minimum) at ϕ R = 0 in both cases.Only φ-axis (the imaginary axis) is shown.The region φ > z −1/2 is not drawn because the Kähler metric becomes negative.
Next, ϕ ± s are supersymmetric minima as the superpotential, D X W and D ϕ W vanish.Finally, the extremum ϕ 0 is a saddle point that is located near (sandwiched by) ϕ ± s .The value of the F-term potential at ϕ = ϕ 0 is The stability around ϕ 0 is dictated by Case µ < 2z: The extrema ϕ = 0, ϕ − 1 have the same properties as in the case µ > 2z.The extremum ϕ + 1 is a minimum, which becomes where we replaced µ by ε 2 through (2.14).The value of the F-term potential at ϕ + 1 is given by ) and hence the minimum energy is positive.

Potential at nonzero ρ
Before we proceed to detailed analysis of the potential at finite ρ, we give a brief summary on the structure of the potential.In this paper, we only consider the case µ > 2z.The other case µ < 2z will be briefly summarised later.Its main drawback is that the vacuum energy cannot be tuned to zero consistently with our perturbative treatment since its leading order value (3.12) is of order unity.
As ρ is turned on perturbatively, the potential will keep for a while the same structure of the extrema, which we denote by 0, ϕ − 1 (ρ), ϕ 0 (ρ), ϕ ± s (ρ).The potential continues to have the extrema 0, ϕ − 1 (ρ) at values of order of the inflaton potential energy, while the other extrema at ϕ 0 (ρ), ϕ ± s (ρ) have much smaller energy of order δ 4 .Remarkably, as we shall see soon below, when ρ becomes greater than some small value ρ, the three points ϕ 0 (ρ), ϕ ± s (ρ) will merge into a single minimum (being pure imaginary) with energy of order δ 4 .In other words, at ρ, the extremum ϕ 0 (ρ) changes from a ridge to a valley.Recall also that 0 and ϕ − 1 (ρ) merge at ρ c to form the critical point as explained in Section 2.1.Thus, the potential at each ρ > ρ slice has only one maximum at ϕ = 0 and one minimum at ϕ 0 (ρ).In particular, the line {ϕ 0 (ρ) : ρ ≥ ρ} forms a valley and the minimum of the potential along this valley V (ρ, ϕ 0 (ρ)) with respect to ρ is the true minimum, namely the vacuum (ρ m , ϕ 0 (ρ m )). Figure 3 illustrates the potential at three ρ-slices.
As mentioned above and will be shown later, the potential along the line {ϕ 0 (ρ) : ρ ≥ 0} is of order δ 4 , while when δ = 0, the line {ϕ 0 (ρ) : ρ ≥ 0} reduces to the supersymmetric flat direction ϕ(ρ) = i 2/µ.This indicates that δ 2 controls the supersymmetry breaking scale, which is true as will be demonstrated in Section 5. Since the inflation scale does not involve δ 2 , the supersymmetry breaking scale is not correlated with the inflation scale.It will also be shown in Section 3.3 that the tuning of the vacuum energy does not affect the inflation scale.
Summary on the case µ < 2z.As ρ is turned on perturbatively, the potential will keep the same structure of the extrema, which we denote by 0, ϕ ± 1 (ρ).The behaviour of 0, ϕ − 1 (ρ) is the same as in the case µ > 2z.On the other hand, ϕ + 1 (ρ) continues to be a minimum (in other words, the potential forms a valley along ϕ + 1 (ρ) for ρ ≥ 0).Therefore, as in the case µ > 2z, the inflaton falls in the φ direction and reach the minimum at some ϕ + 1 (ρ m ) which gives the true minimum among ϕ + 1 (ρ)'s for ρ ≥ 0. However, the vacuum energy V (ρ, ϕ + 1 (ρ m )) is independent of δ 2 , as already indicated in (3.12).This is a consequence of the fact that in the µ < 2z case, the potential cannot have the supersymmetric flat direction when δ = 0, in contrast with the case µ > 2z.The vacuum energy and the supersymmetry breaking scales are both controlled by µ and z, and also the supersymmetry breaking scales are correlated with the inflation scales in general.We will not therefore consider this case further.

Analysis in the absense of the D-term potential
We first analyse the potential in the case without its D-term part.The D-term contribution will be taken into account in such a manner that it does not destroy the structure of the potential obtained by the analysis without the D-term contribution.
ϕ 0 (ρ) is defined by the extremality condition Let us consider the perturbative expansion of ϕ 0 in δ and ρ.First, ϕ 0 depends on δ, ρ only through δ 2 , ρ 2 because so does the potential.Second, when δ = 0, ϕ 0 (ρ) becomes i 2/µ, which is constant in ρ.As a result, the expansion in δ and ρ must have the form According to ϕ 0 (0) = ϕ 0 given by (3.5), the coefficients φ (0,0) , φ (1,0) are given by Next, let us consider the structure of the potential along ϕ 0 (ρ).As we explained above, the goal is to find a nontrivial minimum of the potential along the valley of φ-minima V F (ρ, ϕ 0 (ρ)) with respect to ρ in a perturbative manner.We therefore need its perturbative expansion in ρ at least up to ρ 4 .It turns out that the leading term of the expansion in δ is of order δ 4 , and the coefficients v F 0 , v F 1 , v F 2 are determined only by φ (0,0) , φ (1,0) , being independent of the other coefficients φ (m,n) .The three coefficients are given by where we introduced γ ′ = γ + (Aµ/2) to absorb A. 5 The stability at the extremum ϕ 0 (ρ) can be seen by expanding the potential there: They imply that ϕ = ϕ 0 (ρ) is a minimum if ρ satisfies (as long as ρ is smaller than 1) (3.20)

Adding the D-term potential
Let us turn on the D-term potential.We regard it as a perturbative correction so that the vacuum structure we found for the F-term potential should not get destroyed.Since V F (ρ, ϕ 0 (ρ)) is of order δ 4 , a natural choice of q may be q ∼ δ 2 , so that V D ∝ q 2 ∼ δ 4 .We therefore set We can expand ϕ 0 (ρ) in δ, ρ in the same way as in the case without the D-term.Moreover, the D-term potential vanishes at ρ = 0. Therefore, the expansion of ϕ 0 (ρ) has the form (3.14) with the coefficients φ (0,0) , φ (1,0) given by (3.15).The value of the potential at ϕ 0 (ρ) for each ρ at the leading order in δ is still of order δ 4 due to our choice of q above: and the coefficients v 0 , v 2 , v 4 depend only on φ (0,0) , φ (1,0) , independent of the other coefficients.The three coefficients are given by where v F 0 , v F 2 , v F 4 are given in (3.17).We can also show that the D-term potential does not affect the stability at ϕ 0 (ρ) given by (3.18) and (3.19) for each ρ.Therefore, ρ is the same as (3.20), and ϕ 0 (ρ) is the minimum for each ρ > ρ.
The vacuum is obtained by minimising V (ρ, ϕ 0 (ρ)) with respect to ρ in the region ρ < ρ < 1, Note that v 0 > 0. For V m to be a minimum, we need Let us come back to the consistency conditions for our model.For our scenario to work, ρ < ρ c < ρ m is needed.Combining this with (2.9) and taking the perturbative nature of the analysis above into account, we require ρ < ρ c < ρ m < 1. (3.26)

Tuning of the vacuum energy
Let us require that the vacuum energy V m is zero, V m = 0, which yields v 2 2 = 4v 0 v 4 .Since v 0 > 0, the positivity of v 4 is automatically guaranteed.Therefore, we only have to impose v 2 < 0 out of the conditions in (3.25).Remarkably, the equation v 2 2 = 4v 0 v 4 is just a first-order polynomial in B, so exactly solvable for B. Let us denote the solution by B = B m , which is a quadratic polynomial in γ ′ , where the coefficients are given by . (3.28) Under this choice of B, the ρ-coordinate ρ m at the Minkowski vacuum reads Its positivity is guaranteed by the negativity of v 2 in (3.25).
To check the validity of the model around the vacuum, we need the trace and determinant of the Kähler metric g i j = K i j at the vacuum.For generic ρ, B, they read tr g = 1 + det g + O(ρ 4 ), (3.30) where g i j = K i j .The trace is automatically positive if the determinant is positive.

Positivity of mass squared
The U (1) R gauge boson is massive around the Minkowski vacuum by Higgs mechanism.When the model is coupled with MSSM, the sparticles in the MSSM sector acquire soft masses.We require that they are not tachyonic.The mass squared of the U (1) R gauge boson around the vacuum is given by where the Kähler metric component K X X satisfies The positivity of K X X is guaranteed as long as det g > 0 and z > 0.
Next, let us consider soft scalar masses.Note that the coupling with the MSSM sector does not change the vacuum energy because the vacuum expectation value of each MSSM field vanishes at the vacuum.Let ζ collectively denote squarks and sleptons.The soft scalar mass squared of ζ is given by where Q ζ = 0 in Model 1 and Q ζ = 1/2 in Model 2 for every ζ, and the bracket ⟨ • • • ⟩ means the vacuum expectation value.This mass squared is not positive definite in general because V F must be negative at the Minkowski vacuum, ⟨V F ⟩ = −⟨V D ⟩ < 0. Its explicit expression is given by

Allowed parameter region
Our model depends on the following parameters A, B, z, γ, f, µ, δ, q.We have already fixed B by tuning the vacuum energy to zero.We will find the allowed region for the other parameters by imposing the observational and theoretical constraints.
The observational constraints come from the consistency of the inflation period of our model with the observational data on CMB [26].This yields constraints on the model parameters through the formulas (2.21).We will use the following values [26]: Actually, the number of e-folds ∆N is constrained from below by a value around 34.2 plus a logarithmic contribution depending on reheating temperature [26], ∆N ≳ 34.2 + ln (T reh /1 TeV), allowing for a wider range of ∆N than our choice above.However, in our setup ∆N appears only through the factor e (1−ns)∆N/2 (see eq.(2.22)).Even if we adopt for example ∆N = 40 ∼ 60, this factor varies only from 2.2 to 3.3, which does not substantially affect the allowed parameter region we find below.
We have also imposed various conditions to guarantee the consistency of our analysis.They are summarised as follows: no ghost around the vacuum, (4.7) Their meanings are summarised as follows: The first condition is needed to justify the perturbative expansion in δ 2 .The second one (4.3) is nothing but (3.26).The third one (4.4)means that the inflation is of the single-field type until its end.The fourth one (4.5)requires our model should be in a perturbative regime as an effective theory.The fifth one (4.6)guarantees that the vacuum is indeed a minimum.The sixth one (4.7)requires the theory around the vacuum to be physical.This condition consists of two parts: one for the kinetic terms of the inflaton and waterfall field, det g i j > 0, and the other for the kinetic term of the fermion orthogonal to the goldstino, which will be given later in (5.6).The seventh one (4.8)means that the soft scalar masses for quarks and leptons are not tachyonic.The analysis of the constraints in what follows consists of two parts: one for upper bounds on δ 2 and the other to find an allowed region for the other parameters.To facilitate the analysis, we introduce x by We also introduced We will use x, ε 2 instead of µ, z.Since x, ε 2 > 0, we have µ < 2/3 and z < 1/3.Moreover, µ > 0 yields ε 2 < 1.

Upper bounds on δ 2 and tensor-to-scalar ratio
We first present upper bounds on δ 2 .Technical details are given in Appendix A. The constraints (4.2) and ρ < ρ c in (4.3) yield the following upper bound on δ 2 / √ r: On the other hand, the following upper bound allows us to neglect the D-term contributions q 2 in the expressions (2.21) for A and f : We will choose δ 2 / √ r so that the two upper bounds are satisfied.As will be demonstrated in the next section, δ 2 controls the supersymmetry breaking scale.Therefore, the upper bound on δ 2 gives a strong suppression of the supersymmetry breaking scales and soft masses.
Since the horizon exit ρ * is proportional to √ r according to (2.21), we can find an upper bound on r from the condition ρ 2 end < 1 from (4.3), 4.2 Allowed parameter region for (x, γ ′ ) Let us analyse the other conditions: concretely, where we dropped A < 1 because it has already been achieved in (4.13), and k ϕ φ is the coefficient of the kinetic term of the fermion orthogonal to the goldstino, which will be given by (5.6) in the next section.We will find the allowed region in the (x, γ ′ )-plane, regarding ε, q, r as parameters put by hand.The difference from the last analysis for δ 2 is that the analysis here is independent of δ 2 .
In this section, we will just summarise the result, delegating technical details to Appendix A.
For each x, the conditions ρ end ≤ ρ c and m 0 (Q) 2 > 0 bound γ ′ from above, while ρ c < ρ m < 1 and the conditions for the correct signs of the kinetic terms bound γ ′ from below.In particular, in Model 1, compatibility of the upper bound on γ ′ from m 0 (Q) 2 > 0 with the value of γ ′ for a given ρ 2 m yields an upper bound on q2 , Since ρ 2 m should be substantially smaller than 1 for the validity of the perturbative analysis, it is natural to suppose q < 1.On the other hand, in Model 2, m 0 (Q) 2 > 0 does not yield an upper bound on γ ′ but just xq 2 < 2. Therefore, Model 2 does not have a counterpart of the condition (4.16).
Here we give examples: Figure 4 shows the allowed region for (q, ε, r) = (1, 0.25, 1.0 × 10 −4 ) in Model 1, and Figure 5 shows the allowed region for (q, ε, r) = (0.5, 0.4, 2.0 × 10 −4 ) in Model 1.There the conditions (4.5) are satisfied on the right of the B m = 1 line.The condition m 0 (Q) 2 > 0 does not contribute to the boundary of the allowed regions in these examples because it is weaker than the other conditions.This is true also in Model 2, so the allowed region for Model 2 for each case is the same.
As seen from the figures, the allowed region is upper bounded by the curve for ρ end = ρ c , which is almost constant in x because according to the second inequality of (A.16), its x-dependence is A/(3 + x) with A being as small as 0.005.Its value is given by where we simplified the expression (2.16) for ρ 2 c by neglecting the terms other than 4γ, ε 4 , which is true as long as ε < 1 is chosen as in the examples above. 6We can see that for a given ε, this upper bound on γ ′ is fixed by the tensor-to-scalar ratio.In particular, when x is chosen from the upper bound line ρ end = ρ c , the tensor-to-scalar ratio fixes γ ′ .
On the other hand, since the curve for ρ m = 1 increases monotonically, it will meet the line for ρ end = ρ c at some x, which gives the maximum value for x, which is given approximately as Figure 4: The green region shows the allowed parameter region for (q, ε) = (1, 0.25).The tensor-toscalar ratio is r = 1.0 × 10 −4 .Model 1 and Model 2 have the same allowed region.
where we used q2 x, q2 x 2 ≪ 2. In this approximation, the positivity of its right hand side, which comes from x > 0, gives a lower bound on r, Here we make comments on possibilities of reducing parameters.We already pointed out that under the identification ρ end = ρ c , the tensor-to-scalar ratio fixes the values of γ ′ .Note that γ ′ = 0 is excluded because x > 0 by definition and hence ρ 2 m < 1 forces γ ′ > 0.Moreover, γ = 0, namely γ ′ = Aµ/2, is also excluded: the condition k ϕ φ > 0 forces ε to be very close to 1, which, however, makes ε 4 dominant in the denominator of ρ 2 c in (2.16) and hence ρ 2 c < 0. On the other hand, one may wonder the possibility of B m = 0, but it is also impossible, as demonstrated in the figures: The condition B m = 0 yields two roots that form two curves in the (x, γ ′ )-plane, and they do not intersect with the allowed region.

Supersymmetry breaking scale
In this section, we discuss supersymmetry breaking around the vacuum and the particle spectrum in both hidden and observable sectors.

Supersymmetry breaking scale
The gravitino mass is given by where the bracket ⟨ • • • ⟩ denotes the vacuum expectation value.The non-vanishing vacuum expectation values of the auxiliary fields in the chiral supermultiplets X, ϕ and in the R-vector multiplet are given by The ratio |⟨F ϕ ⟩/⟨F X ⟩| is of order ρ m , and the ratio |⟨D⟩|/|⟨F X ⟩| is of order xq.Therefore, both ratios are smaller than 1 as long as x and q are taken to be smaller than 1 and supersymmetry breaking is dominated by the F-auxiliary of the inflaton multiplet.Thus, its fermionic component is mainly the goldstino that makes it superparnter of the inflaton, as expected in the framework of "inflation by supersymmetry breaking".It is also manifest from the last equations that δ 2 controls the supersymmetry breaking scale.

Bosonic spectrum
The bosonic particles in the hidden (inflaton) sector are two scalars and one U (1) R vector boson.The mass of the vector boson m A and the masses m w , m i of the two scalars at the vacuum are given by where m A comes from (3.32).The expression inside the square root of m i is positive because v 2 < 0 (3.25).The mass m w corresponds to a scalar linear combination dominated mainly by the waterfall field φ, while m i corresponds to a linear combination dominated mainly by the inflaton ρ.The real part of the waterfall field ϕ R has the mass ≃ 4f e 1/µ µρ m equal to m w up to O(δ 2 ).
In the observable MSSM sector, as we have already discussed, the sfermions acquire a common mass m 0 (Q) given by (3.35), where Q = 0 in Model 1 and Q = 1/2 in Model 2.

Fermionic spectrum
The fermion mass spectrum can be read off from the action in the super-unitary gauge given in (B.8).In our setup, the coefficients in the kinetic terms are given by and the other coefficients are zero.The subscript 'R' refers to U (1) R , and the index a refers to the MSSM gauge groups: U (1), SU (2), SU (3).The indices i, j refer to the MSSM chiral superfields. 7he coefficients β a , together with β R for U (1) R in (2.3), are determined by the cancellation of the anomalies involving U (1) R .In our models they read [25] Model 1: Model 2: where g 1 , g 2 , g 3 are the gauge coupling constants for U (1), SU (2), SU (3), respectively.The fermions in the hidden sector are the U (1) R gaugino and the fermion orthogonal to the goldstino which is eaten by the gravitino.Their masses are given by where m f is equal to m w up to O(δ 2 ).In the supersymmetric limit δ = 0, the fermion mass m f and the waterfall boson mass m w become exactly equal, leading to a massive waterfall supermultiplet.
The fermions in the observable sector are the higgsinos and the gauginos for the MSSM gauge groups.The higgsinos have mass µρ m .The gaugino mass m a for the MSSM gauge group labelled by a is given by (5.10) Notice that we can control the overall scale of m a , m i , m 0 , m λ R by tuning δ 2 .All particles with masses proportional to δ 2 constitute the light spectrum of the theory.Besides the MSSM states, these are the inflaton, the R-gauge boson, the gravitino and the R-gaugino.Since the light states charged under U (1) R have comparable mass and charge of order δ 2 in Planck units, it is interesting to check the weak gravity conjecture postulating that gravity is the weakest force in Nature and black-holes decay without remnants [27].This requires the existence of at least one state whose U (1) charge is no less than its mass.In our case, for Model 1, the lightest charged state is the gravitino which has U (1) R charge q/2 and mass m 3/2 , given by (3.21) and (5.1).Imposing the condition q/2 ≥ m 3/2 one finds q ≥ 2ρ m . (5.11) It turns out that in Model 1, this lower bound contradicts the upper bound (4.16).On the contrary, Model 2 survives because q is just subject to the bound q2 < 2x −1 , which, combined with the lower bound (5.11), yields x < (2ρ 2 m ) −1 .On the other hand in Model 2, MSSM fields with odd R-parity are also charged under U (1) R and in particular electroweak gauginos can be lighter than the gravitino (see the example below), which weakens the constraint (5.11).
We finish our analysis by presenting as benchmark model an explicit example with low mass spectrum in the multi-TeV region within the reach of next generation of particle colliders.It contains a neutralino, mostly electroweak gaugino, as the lightest supersymmetric particle (LSP) which is a standard Dark Matter candidate and predicts the existence of the universal inflaton sector (U (1) R gauge boson and gaugino, as well as the inflaton) among the low energy spectrum, which is the smoking gun of inflation by supersymmetry breaking.Let us consider Model 2 in Figure 4, where the choice (q, ρ 2 m ) = (1, 1/4) satisfies the lower bound (5.11) from the weak gravity conjecture applied to the gravitino.We choose (x, γ ′ ) ≃ (0.1454, 0.09318), which is the intersection of the ρ 2 m = 1/4 and ρ end = ρ c lines.The constraint on δ 2 is δ 2 ≪ 2.4 × 10 −8 .For example, we may choose δ 2 so that the masses of the soft scalars and MSSM gauginos are around 10 ∼ 100 TeV.If we choose δ 2 ∼ 10 −14 and use the values for the coupling constants of the MSSM gauge groups (see for instance [25]), g U (1) ∼ 0. Notice that the waterfall field is not part of the light spectrum; it is superheavy, of the order of the inflation scale.
The formula for the soft scalar masses m 2 0 is valid only for the squarks and sleptons, while the masses for the Higgs fields get an additional contribution from the superesymmetric mass term µ H H u H d in the superpotential (2.8) that introduces an additional mass parameter µ H .It is convenient for the following to introduce μH ≡ µ H /m 3/2 which satisfies: where ϕ m is the vacuum expectation value of the waterfall field at the supersymmetry breaking minimum of the scalar potential (3.24).The Higgs masses and mixing are then given by: where we kept only the dominant contribution, neglecting higher order corrections in the Kähler potential and corrections from U (1) R gauge interactions.Obviously, these contributions are important when µ H is large, i.e. μH of order unity, which may be required in some models for the avoidance of the cosmological moduli problem (see e.g.[28]).

Conclusion
In this work we presented an extension of the simple framework called 'inflation from supersymmetry breaking' by introducing a waterfall field, allowing to stop inflation fast and to decouple parametrically the supersymmetry breaking and inflation scales.The idea behind inflation from supersymmetry breaking consists of identifying the inflaton with the superpartner of the goldstino in the presence of a gauged R-symmetry and the inflation sector with the supersymmetry breaking sector in the context of gravity mediation models.The gauge R-symmetry fixes the inflaton superpotential, absorbs its pseudo-scalar companion into the massive R-gauge boson and makes slow-roll parameters naturally small, determined by a quartic correction to the canonical inflaton Kähler potential.Thus, the minimal field content of the inflaton sector consists of two superfields: one chiral (the inflaton) and one vector (the R-gauge boson).
The simplest model has then 4 parameters: the inflation scale f , the quadratic A and quartic B corrections to the inflaton Kähler metric, as well as the U (1) R gauge coupling normalised at the inflaton charge q which can be chosen sufficiently small so that it plays no role during inflation.For positive and small A, the scalar potential has a maximum at the origin of the inflaton field ρ with a plateau, realising hilltop inflation.The second slow-roll parameter η is determined by A, while the first slow-roll parameter ϵ is much smaller: ϵ ∼ η 2 ρ 2 * < η 2 with ρ * the inflaton value at the horizon exit, suppressing the fraction of primordial gravity waves.The value of A is then fixed by the spectral index of primordial scalar perturbations.At the end of inflation, the inflaton rolls down to the minimum of the potential with a tuneable infinitesimally small value using the parameter B. A microscopic supergravity theory generating the above model and its parameters was provided by a generalisation of the Fayet-Iliopoulos model to a U (1) R [3].As mentioned above the coupling of the inflaton sector to the supersymmetric Standard Model is straightforward with an interesting possibility of having the usual R-parity as a remnant of the underlying gauged R-symmetry.Moreover, supersymmetry breaking is transmitted from the inflaton to the observable sector as in gravity mediation models, predicting however a heavy spectrum of superparticles, near the inflation scale.
The extension of the above model to hybrid inflation consists in introducing a new chiral superfield ϕ (the waterfall) neutral under R-symmetry with a superpotential w(ϕ) which is a small deformation of a superpotential w 0 (ϕ) that admits a supersymmetric vacuum: ∂ ϕ w 0 = w 0 = 0 at ϕ = ϕ 0 .One then writes w = w 0 + δ 2 w 1 with δ the deformation parameter that breaks supersymmetry in the vacuum at a scale δ 2 which is different and can be hierarchically smaller than the inflation scale f .For simplicity, we considered as an example a quartic function invariant under ϕ → −ϕ depending on two couplings that can be traded as δ 2 and a mass µ.The latter controls the critical inflaton value ρ c at which the waterfall direction opens up because a vanishing expectation value becomes unstable (tachyonic).We impose that this happens after the end of inflation (of at least 60 e-folds).We expect that our results are independent of the particular form chosen for the superpotential.Finally, in analogy with the inflaton Kähler metric, we introduced quadratic and quartic (in the inflaton) corrections to the waterfall Kähler metric with corresponding parameters z and γ; these are free parameters of the model, subject to theoretical and experimental constraints.
The requirement of perturbative treatment of both inflation and potential minimum, together with detailed investigation of the parameter space done in Section 4, leads a light spectrum of order δ 2 ≪ f that contains besides the usual MSSM fields (scalar quarks and leptons and SM gauginos) the universal inflaton sector of inflation by supersymmetry breaking: U (1) R gauge boson and gaugino together with the inflaton.The light mass spectrum has the following approximate pattern: the electroweak gauginos tend to be the lightest as long as shown in the explicit example we presented in Section 5.The next lightest group contains the squarks and sleptons, U (1) R gauge boson and gaugino, and the gravitino, having masses of more or less the same order.Then the inflaton comes next a bit heavier.This is a model-independent feature of our framework which connects supersymmetry breaking and inflation in a minimal way.Its embedding in a fundamental theory, such as string theory, remains an open challenging question.

B Fermion masses
Let Φ I collectively denote the chiral superfields X, ϕ and the MSSM ones denoted by ζ, and χ I be the fermionic component of Φ I .The fermions in our models have the following quadratic Lagrangian: In the superunitary gauge, the goldstino η is set to zero: Eliminating χ X from the action by solving this constraint gives where the coefficients in the kinetic terms are given by (B.9) (B.10) 11) and the coefficients in the mass terms are given by where F I = g I J F J .

− 1 2 2 4
ig I J χ J σµ ∂ µ χ I − 1 2 ih * AB λB σµ ∂ µ λ A − 1 2 m IJ χ I χ J − m IA λ A χ I − 1 2 m AB λ A λ B + c.c. , (B.1)where the labels i, j, • • • refer to the chiral superfields except X.The labels A, B, • • • refer to the gauged U (1) R and the MSSM gauge groups U (1), SU (2), SU (3), which we denote by G A .The labels a, b, • • • refer to the MSSM gauge groups only, which we denote by G a .The gauge kinetic functions are chosen to beh AB = h A δ AB , h A = 1 + β A ln X, (B.2)where the coefficients β A are determined by the cancellation of the anomalies involving U (1) R .The mass coefficients are given bym IJ = W IJ + K IJ W + K I ∇ J W + K J ∇ I W − K I K J W − Γ K IJ ∇ K W, i∂ I h AB D B , (B.4) m AB = 1 2 F I ∂ I h AB , (B.5)whereK IJ = ∂ I ∂ J K and Γ K IJ = g LK ∂ J g I L. X I A isthe Killing vector for Φ I , defined by splitting off transformation parameters from the infinitesimal transformation of Φ I under gauge group G A .In our setup m aI = 0 holds since ϕ is neutral under U (1) R and the MSSM fields have vanishing expectation values at the vacuum.The auxiliary fields are defined by F I = −e K/2 g I J ∇ J W , D A = iX I A (K I + W −1 W I ), D A = (h −1 ) AB D B .(B.6) 5, g SU (2) ∼ 0.59, g SU (3) ∼ 0.72,