Inflaton decay in No-Scale Supergravity and Starobinsky-like models

We consider the decay of the inflaton in Starobinsky-like models arising from either an $R+R^2$ theory of gravity or $N=1$ no-scale supergravity models. If Standard Model matter is simply introduced to the $R + R^2$ theory, the inflaton (which appears when the theory is conformally transformed to the Einstein frame) couples to matter predominantly in Standard Model Higgs kinetic terms. This will typically lead to a reheating temperature of $\sim 3 \times 10^9$~GeV. However, if the Standard Model Higgs is conformally coupled to curvature, the decay rate may be suppressed and vanishes for a conformal coupling $\xi = 1/6$. Nevertheless, inflaton decays through the conformal anomaly leading to a reheating temperature of order $10^8$~GeV. The Starobinsky potential may also arise in no-scale supergravity. In this case, the inflaton decays if there is a direct coupling of the inflaton to matter in the superpotential or to gauge fields through the gauge kinetic function. We also discuss the relation between the theories and demonstrate the correspondence between the no-scale models and the conformally coupled $R+R^2$ theory (with $\xi = 1/6$).


Introduction
A key aspect of the difference between the 'old' model of inflation, based on a firstorder phase transition [1], and the 'new' inflationary Universe [2] is the ability to recover from the era of exponential expansion, or a graceful exit from inflation [3].In old inflation, bubble percolation was incompatible with resolving the cosmological problems tackled by inflation.In new inflation (and in all subsequent 'newer' inflationary models), a graceful exit is a built-in feature.Indeed, a graceful exit is possible regardless of whether inflation takes place at small field values rolling toward large field values (hilltop inflation), or at large field values rolling to the origin in field space (plateau models).The evolution of the field driving inflation, the inflaton, 1 is governed by its equation of motion and begins with a period of slow-roll characterized by exponential expansion, a ∼ e H I t , where a is the cosmological scale factor and H I is a nearly constant Hubble parameter during inflation.Slow-roll naturally ends with an oscillatory phase where the inflaton oscillates about a minimum, and a matter-dominated phase, a ∼ t 2/3 begins.
However, at some point, inflaton oscillations must give way to a radiation-dominated Universe.This transition is most easily accomplished through the decay of the inflaton [4,5], leading to a period of reheating.Provided that the inflaton decay products scatter and thermalize [6], the standard radiation-dominated era of the early Universe is produced.In this review, dedicated to the accomplishments of Professor Dimitri Nanopoulos, we revisit the possibilities for inflaton decay in Starobinsky-like models, and explore these same models when derived from no-scale supergravity [7].
The Starobinsky model [8] predates the original models of inflation based on a GUT phase transition, and was an attempt to find a solution for a singularity-free Universe.The model is based on an R + R 2 theory of gravity, which is conformally equivalent to a theory that combines Einstein gravity with a scalar field, potentially serving as the inflaton [9].Furthermore, the first calculation of a nearly flat scale-independent density fluctuation spectrum was done in this context [10].
Shortly after the advent of the new inflationary model, it was recognized that supersymmetry may play an important role [11] in preventing an inflationary hierarchy problem, which requires the inflaton mass to be ≪ M P , where M P = 2.4 × 10 18 GeV is the reduced Planck mass.Supergravity is the extension of supersymmetry that incorporates gravitational interactions and it is natural to construct inflationary models consistent with supergravity.However, minimal supergravity models are plagued with the so-called η problem [12], in which scalars tend to obtain large (Hubble scale) masses, preventing an extended period of inflation.This problem is alleviated [13] in no-scale supergravity models [7], formulated on a maximally symmetric field-space manifold and have been derived as the effective low energy theory from string theory [14].Of particular interest here is the class of inflationary models based on no-scale supergravity, which yield a model related to the Starobinsky model [15][16][17].
In this contribution, we will focus on inflaton decay in Starobinsky-like models and related models based on no-scale supergravity.Some of these results may be applicable to T-model attractors [18] as well.We begin by reviewing the Starobinsky model, augmented with a matter sector that, for simplicity, consists only of the Standard Model Higgs boson (or the supersymmetric Higgs, as the case may be).We also consider the possibility that the Higgs boson is coupled to curvature (with a coupling ξ), where the decay rate of the inflaton vanishes when ξ = 1/6.In this case, the dominant decay mechanism comes from the coupling to T µ µ through the conformal anomaly [19].We also consider Starobinsky-like models derived in no-scale supergravity and examine the prospects for decay in these models.In general, inflaton decays in no-scale models are highly suppressed [20,21].In the absence of a direct (superpotential) coupling of the inflaton to the Standard Model, decay to gauge bosons through the gauge kinetic function may be the dominant channel for inflaton decay.For a review on building inflationary models in no-scale supergravity, see Ref. [22].
In what follows, we begin with the Starobinsky model as a model of inflation and compare the inflationary observables-the tilt in the perturbation spectrum and the tensorto-scalar ratio-to CMB results.We consider several possibilities for coupling the Standard Model to the R + R 2 theory and calculate the inflaton decay rate.Even in the Einstein frame, fields are generally not canonically normalized and we discuss the procedure for determining canonical field coordinates.We also discuss inflaton decay through the conformal anomaly.In Section 3, we discuss the derivation of a Starobinsky-like model in the context of no-scale supergravity.We also review the possibilities for inflaton decay in this context.In this case as well, fields are generally non-canonical, and we briefly review the transformation to canonical field coordinates.Finally, in Section 4, we relate the formulations of the inflationary model in the contexts of R + R 2 gravity and no-scale supergravity.There is a close correspondence between the two when the Higgs is conformally coupled with ξ = 1/6.Our conclusions are presented in Section 5.

The Starobinsky model and inflaton decays
The Starobinsky model of inflation [8] can be formulated by adding an R 2 term to the conventional linear Einstein-Hilbert action: where α is a constant.Here, we work in units where M P = 1.Upon introducing a Lagrange multiplier −αR 2 → 2αΦR + αΦ 2 , this action may be rewritten in the form Then, following the conformal transformation [9,23] gµν = e 2Ω g µν = (1 + 2αΦ)g µν , we can rewrite the action in the Einstein frame as Note the appearance of an additional scalar degree of freedom, with a non-canonical kinetic term and potential V(Φ).By making a field redefinition, ϕ ≡ √ 3/2 ln(1 + 2αΦ), Eq. ( 4) may be written as: When α = 1/6M 2 , the potential is now seen with the well-known form of the Starobinsky potential, This potential is shown in Fig. 1.
As an inflationary model, the potential (6) leads to observables in excellent agreement with observations by Planck and other CMB experiments [24][25][26].The comparison to experiment is easily made by first computing the slow-roll parameters, ϵ and η for single field inflationary models, where the prime denotes a derivative with respect to the inflaton field ϕ.Of particular interest is the tilt in the spectrum of scalar perturbations, n s [24], = 0.965 ± 0.004 (68% CL) , (9) and the tensor-to-scalar ratio, r [25], r ≃ 16ϵ * < 0.036 (95% CL) .(10) The overall inflationary scale, M, is set by the amplitude of scalar perturbations, A s [24], The number of e-folds, N i , of inflation between the initial and final values of the inflaton field is given by the formula The number of e-folds between when the pivot scale k * exits the horizon and the end of inflation is denoted by N * .Inflation ends when a = a end and ä = 0. Typical values of N * are in the range ∼ 50-60, depending on the mechanism that ends inflation [24,27,28].For the Starobinsky model, it is straightforward to compute the inflationary observable listed above [15]: For the Starobinsky potential, ϕ end = ϕ(a end ) = 0.62M P [29,30] and N * = 55 for ϕ * = 5.35M P .In this case, we find M = 1.25 × 10 −5 M P ≃ 3 × 10 13 GeV, n s = 0.965, and r = 0.0035.We show in Fig. 2 taken from Ref. [30], the 68% and 95% CL regions of the (n s , r) plane allowed by the Planck and Keck/BICEP2 data [24,25].The (nearly) horizontal line (labeled α = 1) 2 corresponds to the predictions of the Starobinsky model (6) for N * ∼ 46 to ∼ 52 e-folds, corresponding to limits from reheating (see below) 100 GeV < T RH < 10 10 GeV, the latter corresponding to the upper limit on T RH from the limit on gravitino production and the relic density of a 100 GeV lightest supersymmetric particle.The other nearly vertical curves correspond to (from left to right), T RH = 2 MeV (solid), N * = 50 (dotted), T RH ≃ 5 × 10 14 GeV, obtained when a perturbative decay coupling is O(1) 8π m ϕ ) (solid), and N * = 60 (dotted).For the α-Starobinsky model [16,31,32], we shade the region in green that respects the constraint T reh > T EW and the relic density constraint of m LSP ≃ 100 GeV.For more details, see [30].
As described above, the α-Starobinsky model will inflate and produce density fluctuation in agreement with the experiment.When ϕ < ϕ end , the inflaton begins a period of oscillations, and the energy density in these oscillations scales as a −3 , typical of a matterdominated universe.These oscillations continue until inflaton decays occur sufficiently fast to produce radiation.How this occurs depends on how the Standard Model is introduced in the context of inflationary theory and how the inflaton couples to the Standard Model (SM).
We first consider a simple model by adding a complex scalar field H to the R + R 2 action: After performing the same conformal transformation (3) and the field redefinition ϕ ≡ √ 3/2 ln(1 + 2αΦ), the above action can be rewritten as Expanding the second line around the minimum, ⟨ϕ⟩ = 0, we obtain a canonical kinetic term for the complex scalar, ∂ µ H * ∂ µ H, as well as the couplings of the inflaton to H: The dominant contribution to the decay ϕ → HH comes from the first term (the 2nd term gives a contribution which is suppressed by the Higgs mass) and the rate is given by where N H = 4 is the number of the real scalar degrees of freedom of the SM Higgs doublet.
It is often useful to perform a field redefinition to obtain canonical fields in a particular background using the Riemann normal coordinates (RNC).The term ϕ∂ µ H * ∂ µ H can be eliminated and the coupling of the inflaton to Higgs is purely a potential interaction.The RNC fields can be obtained quite generally.Consider a Lagrangian of the form where g I J (Φ) is a symmetric metric and V(Φ) is the potential.Here we use the Latin indices I, J, ..., for the flavor eigenbasis.To compute the interactions around the minimum, we introduce the field expansion Φ I = v I + ϕ I , where v I is the vacuum expectation value and ϕ I is the field fluctuation.Following the discussions in [33][34][35], the symmetrized covariant derivatives of the potential can be expressed as where the covariant derivatives are symmetrized and include all permutations multiplied by a factor of symmetry 1/n!.Here we take the derivatives with respect to field Φ I , evaluated at their VEV v I .The Christoffel symbols are given by Using these definitions, we find that the canonically normalized mass matrix at the minimum is given by To canonically normalize the fields and transform from the flavor to mass eigenbasis, we introduce a tetrad that flattens the metric Here, the lowercase letters correspond to the mass eigenbasis, and the indices may be raised or lower using the Kronecker delta δ ij or δ ij .The inverse tetrad can be computed using the expression e i I = δ ij e J j g J I .In Riemann normal coordinates, there are no cubic derivative interactions, which significantly simplifies the computation.Thus, in this case, the cubic interaction can arise only from the potential terms.From the action (17), we find that the field metric of the kinetic terms 1  2 g I J (Φ)∂ µ Φ I ∂ µ Φ J can be expressed as where Φ = (ϕ, Re H, Im H).At the minimum, the VEVs of these fields are v = (0, 0, 0).We find that the tetrad at the minimum is given by Finally, we compute the three-particle scattering amplitude using Eq. ( 21) that includes a symmetry factor for two identical final states in the amplitude.To compute the scattering amplitude in the mass eigenbasis, we use and find These couplings lead to the decay rate (19).For more details on how to use the RNC to compute the scattering amplitudes, see Refs.[33][34][35].
Another alternative description of the model entails a conformal transform into the conformal frame where all the scalar fields conformally couple to the gravity.This is achieved by the conformal transformation with the inflaton ϕ and the conformal factor satisfying Then the above action can be rewritten as [37] A In this frame, the scalar fields are canonical and conformally coupled to gravity without any kinetic mixing, and the decay is induced by the cubic interaction ϕ c |H c | 2 from the potential.Focusing on the Standard Model Higgs with V(H) = m 2 |H| 2 + λ|H| 4 , the decay is predominantly induced by the last term in Eq. ( 30), and we recover the decay rate given in Eq. ( 19), which is expected due to frame independence.One may also introduce a non-minimal coupling to gravity ξRH * H: The conformal transformation in this case is In the Einstein frame, this is equivalent to where now ϕ = √ 3/2 ln(1 + 2αΦ − 2ξ|H| 2 ).To compute the decay rate in this frame, one can again expand the couplings around the minimum, which gives rise to the following dominant trilinear couplings: The rate is then 3 3 See also [36] for discussions of inflaton decay in these models.
Alternatively, we may move to the conformal frame by with the inflaton ϕ and the conformal factor satisfying and we obtain [37] A This form of action makes it manifest that the interaction between the inflaton and the Higgs is controlled by the combination 1 − 6ξ.The decay is again induced by the cubic interaction ϕ c |H c | 2 from the last term in the potential, 4 and not surprisingly the rate is given by Eq. ( 35).Note that a very large negative value for ξ forces the Higgs to participate in the inflationary dynamics, making it the so-called Higgs-R 2 inflation model [39].In this case, the (p)reheating dynamics can be non-perturbative [40], and we do not consider this regime in the following.
The above formula shows that, if the Higgs conformally couples to gravity, with ξ = 1/6, the inflaton does not decay into the Higgs.In this case, the inflaton predominantly decays into Standard Model gauge bosons through the trace anomaly [19].In general, the coupling between the inflaton and Standard Model particles in the Starobinsky model is expressed to the lowest order as where T µ µ is the trace of the stress-energy tensor.One can see that Eqs. ( 18) and (34) indeed take this form.If the matter sector is (classically) conformal, and assuming the absence of direct coupling between the gauge sector and the gravity beyond the minimal one in the original frame, T µ µ is predominantly sourced by the trace anomaly.The origin of the trace anomaly depends on the regularization scheme (although the final result itself does not) [41].For instance, in dimensional regularization, the inflaton couples to the gauge bosons as ϵϕF µν F µν , with the spacetime dimension d = 4 − 2ϵ, after the conformal transformation.This factor of ϵ is compensated by the 1/ϵ pole arising from the wavefunction renormalization (and self-interactions in the non-abelian case), resulting in a finite term in the limit ϵ → 0 proportional to βϕF µν F µν , where β is the beta function counting both heavy and light degrees of freedom.In addition, the inflaton directly couples to heavy degrees of freedom through their mass terms.After integrating out heavy degrees of freedom, the contribution from the mass term cancels with the one from the wavefunction renormalization, leaving only light degrees of freedom, as demonstrated in [41].As a result, the trace anomaly is given by where i runs over all the different gauge groups, a runs over the generators of each gauge group, and α i is the fine structure constant of the gauge group i.The beta function is given by where N A i is the number of the gauge bosons in the gauge group i, and we again explicitly include the Planck scale.To evaluate this expression, we use the couplings at the scale of M/2 after the renormalization group running.In the case of the Standard Model, by running the couplings up to two-loop with SARAH [42] and taking the input values at the electroweak scale following [43], we obtain α 1 (M/2) ≃ α 2 (M/2) ≃ 0.024 and α 3 (M/2) ≃ 0.027, where M ≃ 3 × 10 13 GeV.
There are of course numerous other ways to introduce the SM Lagrangian.For example, one can introduce it directly in the Einstein frame leading to an action similar to the one in Eq. ( 17) without the coupling of the inflaton to the Higgs kinetic and potential terms.In this case, there is no coupling of the inflaton to matter, though to achieve this starting from the Jordan frame is quite contrived.Hence we restrict ourselves to the above simple possibilities.
If there is a non-negligible decay width for the inflaton, its decays will start populating the radiation bath and initially there will be a rapid increase in radiation temperature followed by a slow redshift with ρ r ∝ a −3/2 [44][45][46].We define the reheating temperature corresponding to ρ r (a RH ) = ρ ϕ (a RH ), or g RH π 2  30 where g RH is the number of degrees of reheating at T RH .Subsequently, the radiation will redshift normally as ρ r ∝ a −4 .In the SM, g RH = 427/4 and for M = 3 × 10 13 GeV, we find where the former comes from ϕ → HH while the latter from ϕ → AA.Thus, in the Starobinsky model, there is always a significant source of reheating.Before concluding this section, we note that many of the above arguments can be applied to T-model attractors [18].Although originally formulated using an action with an O(1,1) symmetry, these models can also be simply derived from a Jordan frame with an action [47] A with f (ϕ) = 1 − ϕ 2 /6.After a transformation to the Einstein frame we have Taking V( f ) = λ f 2 ϕ k and making the field redefinition ϕ = √ 6 tanh(χ/ √ 6), we obtain Possible decay modes of the inflaton χ, can be obtained through the induced couplings to the Higgs bosons (once introduced) or through the gauge anomaly.

No-scale inflationary models and inflaton decay
The bosonic sector of an N = 1 supergravity theory is specified by a Kähler potential, K(ϕ i , ϕ * j ), which determines the field-space geometry of the chiral scalar fields in the theory, a holomorphic function of these fields, W(ϕ i ), which determines the interactions between these fields and their fermionic partners, and a gauge kinetic function, f αβ (ϕ i ).Taken together, the bosonic Lagrangian can be written as where the first term is the minimal Einstein-Hilbert term of general relativity and in the second term K j i ≡ ∂ 2 K/∂ϕ i ∂ϕ * j is the field-space metric.The effective scalar potential, where and G i ≡ ∂G/∂ϕ i , G j ≡ ∂G/∂ϕ * j , and G −1 i j is the inverse of the matrix of second derivatives of G.In addition, there are also D-term contributions for gauge non-singlet chiral fields.For a review of local supersymmetry, see Ref. [48].Minimal supergravity (mSUGRA) is characterized by a Kähler potential of the form in which case the effective potential (49) can be written in the form where W i ≡ ∂W/∂ϕ i .Unlike the scalar potential in globally supersymmetric models, where V = |W i | 2 , the minimal supergravity potential is not positive semi-definitive.Indeed, the negative term ∝ |W| 2 in (52) generates in general minima with V ∼ −O(m 2 3/2 M 2 P ) [49], where m 2 3/2 = e G is the gravitino mass.In addition, as noted earlier, it is difficult to generate flat directions suitable for inflation, as scalar fields tend to pick up order H 2 masses [12].
In contrast, such difficulties are absent in no-scale supergravity.The simplest such (single field) theory is defined by [7,50] This describes a maximally symmetric field space with constant curvature R = 2/3 and for W = 0, V = 0.The theory is generalized by adding additional scalar fields with [51,52] This Kähler potential still describes a maximally symmetric field-space with curvature R = (N + 1)/3, where N − 1 is the number of 'matter' fields, ϕ i .In this case, the Lagrangian becomes where the effective scalar potential can be written as When W T = 0, the potential takes a form related to that in global supersymmetry, with a proportionality factor of e 2K/3 , where K is the canonically-redefined modulus.Large mass terms are not generated [53] in this case, and the η-problem is avoided [13].
Interestingly, no-scale supergravity can be used as the framework to construct models of inflation [54][55][56][57].In particular, by adopting a very simple (Wess-Zumino) form for the superpotential of a single matter field, ϕ, 6 of the form [15] the Starobinsky model potential ( 6) is obtained once a field redefinition is made to a canonically normalized field χ, and T is stabilized with a vacuum expectation value assumed here to be ⟨T⟩ = 1 2 .Decomposing χ into its real and imaginary parts: χ = (x + iy)/ √ 2, the potential is minimized for y = 0, and in the real direction we obtain the Starobinsky potential with the identification of ϕ in Eq. ( 6) with x here.
Note that the choice of superpotential in Eq. ( 58) is not unique.There is another well-studied example, defined by [58][59][60][61][62] In this case, the Starobinsky model potential ( 6) is obtained once a field redefinition is made to a canonically normalized field, 6 To construct a Starobinsky-like model, at least two fields, T and ϕ are needed [16].
where t is real (t here plays the role of the inflaton ϕ in Eq. ( 6)) and ⟨ϕ⟩ = 0. Indeed there are multiple classes of such models all related by the underlying SU(2,1)/SU(2)×U(1) no-scale symmetry [17].
Inflation in all of these models is indistinguishable from the original Starobinsky model, and the inflationary observables are the same as those given in Eqs. ( 11)- (15).However, as in the discussion of the previous section, the matter sector must also be considered in order to achieve reheating and a radiation-dominated universe.In this case, there are a few possibilities.If we are considering only the possibility for decays to the Higgs boson, the superpotential must include a µ-term (in the minimal supersymmetric Standard Model (MSSM)), The Higgs kinetic terms may arise from the Kähler potential either as untwisted fields or as twisted fields For the case of untwisted fields and the inflaton is ϕ (as opposed to T), it was shown that if the superpotential dependence on the inflaton does not extend beyond W I , from the canonical mass matrix, it is straightforward to see that there are no decay terms for the inflaton to scalars or fermions [20,21].This result can be readily seen if one performs the necessary field redefinitions to canonical fields.By defining Kähler normal coordinates (KNC), we can rewrite the Lagrangian and read off all possible couplings of the inflaton [34,35].
We follow the same procedure as for the RNC.We consider a general theory with N massive complex scalar fields denoted by Z I with general two-derivative interactions and a potential V(Z, Z), given by the Lagrangian7 where K I J is a metric tensor that is Hermitian.Here the use the holomorphic (unbarred) indices on the left and anti-holomorphic (barred) indices on the right, with K I J K M J = δ I M .We expand the complex fields as where w I is the complex scalar field VEV and z I is the field fluctuation.We introduce the complex tetrads that flatten the Kähler metric When we use the complex notation, the Greek indices are used for the mass eigenbasis, and the indices can be raised by using the Kronecker delta δ α β and δ α β.The inverse complex tetrad e α I (w, w) can be computed using the relation e α I (w, w) = δ α βe J β K I J .We introduce the covariant derivatives acting on the scalar potential V(Z, Z): where these covariant derivatives are not symmetrized and we take the derivatives with respect to the complex fields Z i and Z Ī .Following Ref. [63], the non-zero Christoffel symbols are given by and the scalar mass matrix can be expressed as From the Kähler potential (63), we find the following Kähler metric At the minimum, the VEVs of the complex fields Z = (T, ϕ, H u , H d ) are given by w = ( 1 2 , 0, 0, 0), and the complex tetrad can be expressed as If we compute the effective scalar potential (49) with the Cecotti superpotential (60) combined with W SM = µH u H d , and use the covariant derivatives (68) in the mass eigenbasis, we find that the inflaton fluctuation δT contains the following trilinear couplings to MSSM Higgs fields:8 where we have ignored a coupling of T to ϕ 2 as it does not lead to reheating.If instead we use the Wess-Zumino superpotential (58), the inflaton fluctuation δϕ does not have any trilinear couplings.Note these couplings are proportional to µ.The decay rate is given by [21,[64][65][66][67]] for decays into the eight real scalar Higgs fields.For µ ∼ O(1) TeV, the reheating temperature is T RH < .01eV.In the case of high scale supersymmetry (with µ ≳ m ϕ ), the reheating temperature can be significantly higher.
If we repeat the same procedure for the Kähler potential with twisted Higgs fields (64), we find that the Kähler metric is given by Using the same VEVs at the minimum as before, we find the complex tetrad In this case for the Cecotti superpotential combined with W SM = µH u H d , we find the following trilinear couplings to the MSSM Higgs for the inflaton fluctuation δT and for the Wess-Zumino model, there are no trilinear couplings for the inflaton fluctuation δϕ.These couplings, though 3 times larger, still provide a dismal amount of reheating.
Note that there are of course many other potential couplings of the inflaton to MSSM scalars.For example, there is a three-body decay to a Higgs, stop, and antistop, with a decay rate µ 2 y 2 t M/2048π 3 M 2 P giving a reheating temperature of order 5 MeV [21,22].Four-body decay rates into stops (and anti-stops) are more promising and provide a rate 9y 4  t M 3 /262144π 5 M 2 P , giving a reheating temperature of order 10 7 GeV, provided the MSSM fields are in the twisted sector.For untwisted matter fields, this four-body rate vanishes.Including MSSM fermions, three-body decays into Higgs, top and antitop lead to the rate 3y 2 t M 3 /2048π 3 M 2 P for twisted fields, corresponding to T RH ∼ 10 9 GeV.For more information, see, [21,22].
It is also possible that the inflaton does couple directly to SM fields if, for example, there is a superpotential coupling such as [68,69] to the up-like Higgs and lepton doublets.When considered with the superpotential (58) there is an interaction term −My ν H * u L * ϕ 1 , and the inflaton decay width is given by where we have neglected the masses of the final-state particles.There is also a decay to fermions ϕ 1 → H 0 u ν, H + u fL with a rate equal to that in (80).Using Eq. ( 43) for T RH and the MSSM value g RH = 915/4 with M = 3 × 10 13 GeV, T RH = 4.8 × 10 14 y ν GeV . (81) Clearly, this is a very efficient way to reheat the Universe, if such a coupling exists.
It is also possible for the inflaton to decay to gauge bosons (and gauginos) if the gauge kinetic function f αβ has a non-trivial dependence on the inflaton [20,21,70,71].In this case, the decay rate is given by 9 where N G = 12 in the Standard Model (assuming a universal coupling of the inflaton to all gauge bosons), and d g is given by This leads to a reheating temperature of Finally, we note that the T-model potential for inflation found in Eq. ( 47), can also be derived in no-scale supergravity.For Wess-Zumino-like models where the inflaton is ϕ, a choice of superpotential [45] Alternatively, choosing yields the same potential when T is associated with the inflaton.

Relating
No-scale supergravity and R + R 2 It is not coincidental that the Starobinsky potential can be derived from an R + R 2 theory of gravity as well as no-scale supergravity [72,73].The standard formulation of N = 1 supergravity is in the Jordan frame, with L ⊃ − 1  6 ΦR and only after a conformal transformation with e 2Ω = 1 3 Φ do we arrive at the theory defined in the Einstein frame with a Kähler potential given by K = −6Ω = −3 ln(Φ/3) . (87) More specifically, starting with [74] where 9 The numerical prefactor of this rate differs from the prefactor of the rate computed in [71] by a factor of 2. The reason is the definition of d g in terms of a derivative with respect to the complex field ϕ and not the physical inflaton, which is the canonically normalized real part, √ 2 Re ϕ.
and g is a holomorphic function of the scalars ϕ (indices on the scalars ϕ i have been dropped for clarity) and J = 6Ω = 3 ln(Φ/3) is a function of ϕ, ϕ * .This expression only includes the purely scalar part of the Lagrangian and the gravitational curvature.Collecting the kinetic terms in Eqs. ( 89) and (90), we have Setting K = −J, we see that this corresponds to the first two terms in Eq. ( 48).The remaining terms in L aux can be reorganized to give the scalar potential given in Eq. ( 49), with g = 2W.
The simplest no-scale theory defined in Eq. ( 53) is obtained when Φ = 3(T + T * ) or 3(T + T * − ϕϕ * ) in the more general theory with additional chiral fields as in Eq. ( 54).
Notice now that the Starobinsky model can be matched to the real part of the supergravity theory with the identification of (1 + 2αΦ) from Eq. (3) to Φ/3 in Eq. ( 87).The latter can then be identified with T + T * to obtain the Kähler potential in Eq. ( 53).
If we include matter fields in the Starobinsky model as in Eq. ( 31), where H is a representative example of the matter fields ϕ i , we obtain the Kähler potential in Eq. ( 54) if ξ = 1/6.In fact, this was to be expected as we have seen that the supergravity couplings of the inflaton lead to vanishing decay rates from kinetic terms which is precisely the case in the R + R 2 model when matter fields are conformally coupled with ξ = 1/6 as in Eq. (35).
We have seen that the conformal transformation to the Einstein frame in the R + R 2 model gives rise to a potential which becomes the Starobinsky potential when we set ϕ = √ 3/2 ln(1 + 2αΦ).When matter fields are included, the potential is and now we must set ϕ = √ 3/2 ln(1 + 2αΦ − 2ξ|H| 2 ).The second term here again closely resembles the scalar potential in no-scale supergravity given by Eq. ( 56) when we make the identification of V(H) to V, 1 + 2αΦ with (T + T * ), ξ = 1/6, and H to a more general ϕ i .In the supergravity case of course V must be determined by the superpotential as in Eq. (57).To completely match the scalar potentials in the R 2 and no-scale theories, additional potential interactions need to be added to the action in Eq. ( 31) once a superpotential has been specified.
We have seen that matter fields introduced in the no-scale framework appear as conformally coupled fields with ξ = 1/6 accounting for the lack to decay channels beyond those proportional to powers of the scalar masses which breaks the conformal symmetry.Previously, we have also seen that inflaton decays are expected to occur though the trace anomaly, and we expect the same to be true in the context of supergravity.Indeed, the explicit factor of e 2 3 K in Eq. ( 56) is nothing other than the conformal factor e −4Ω = e −2 2  3 ϕ in Eq. ( 17), with K = −6Ω.The classical couplings of K to scalar potential and kinetic terms may also be found from the coupling as in Eq. ( 39).We expect therefore that at the quantum level, there is a coupling of K to T µ µ | anom in Eq. (40).Thus up to a numerical factor, we further expect that this coupling will induce a decay of the inflaton [75] as in Eq. ( 42). 10 This coupling was also considered in the context of anomaly mediated supersymmetry breaking [76].

Summary
Our view of inflation has evolved significantly since the original first order GUT phase transition proposed by Guth [1].The exit from accelerated expansion occurs naturally in slow-roll inflation.The problem of reheating is intimately connected with the detailed model of inflation and how the inflaton couples to the Standard Model.In addition, far from being simply an abstract construct employed to solve cosmological problems, inflation has become a theory with testable experimental predictions.Even simple models of inflation typically make three predictions: The overall curvature, Ω = 1; the tilt of the scalar anisotropy spectrum, n s ≲ 1, and the scalar-to-tensor ratio, r.From experiment, we have Ω = 0.9993 ± 0.0019 [24]; n s = 0.9649 ± 0.0042 [24]; and r < 0.036 95% CL [25] or r < 0.032 95% CL [26].These values (and limits) can be compared for example to the predictions of the Starobinsky model [8] which gives Ω = 1, n s = 0.965, and r = 0.0035.It should be noted that n s (and to a lesser extent r) depends on the number of e-folds which in turn depends on the reheating temperature [30].
For this paradigm to work, reheating and the production of Standard Model particles must occur.Adding couplings of the inflaton to the SM can in some cases become problematic if they distort the inflaton potential and spoil the positive aspects of the inflationary expansion.In this work, we examined the question of inflaton decay and reheating in two related frameworks for inflation.When formulated as a modified theory of gravity with a Lagrangian given by 1  2 (−R + αR 2 ) as in the Starobinsky model [8], we have seen that subsequent to the conformal transformation to the Einstein frame, couplings to Standard Model fields are automatically generated leading to a decay width proportional to m 3 ϕ /M 2 P and a reheating temperature of order 10 9 GeV.This is the case so long as the Standard Model fields are not conformally coupled to curvature with conformal coupling ξ = 1/6 in which case the coupling of the inflaton to the Higgs field vanishes.However, even in this case, at the quantum level, there is a coupling of the inflaton to gauge fields through the trace anomaly leading to decay and a reheating temperature of order 10 8 GeV.
We have also considered the inflationary models constructed in the framework of no-scale supergravity [7].Once a relatively simple superpotential is specified [15,[58][59][60][61][62] (as in Eq. (58) or Eq. ( 60)), a Starobinsky-like potential is generated yielding the same predictions for the inflationary observables.In this case, unless the inflaton is directly coupled to the matter (e.g., by associating the inflaton with the right-handed sneutrino [68,69]), inflaton decay is highly suppressed.This can be understood when relating the no-scale supergravity models to the R + R 2 models as was done in the previous section, and one can see that the two theories are related when Standard Model fields are in fact conformally coupled with ξ = 1/6 in the R + R 2 theory.We expect that in this case too, inflaton decay through the trace anomaly is possible.Alternatively, coupling the inflaton to gauge fields through the gauge kinetic function may also lead to inflaton decay and reheating.
We expect the next significant test of these models will be available in the next round of CMB experiments which can probe the tensor to scalar ratio down to 0.001 and should either confirm or exclude the type of models discussed here which predict r ∼ 0.003.