Next Article in Journal
Energetic Particles and High-Energy Processes in Cosmological Filaments and Their Astronomical Implications
Next Article in Special Issue
Lorentzian Quantum Cosmology from Effective Spin Foams
Previous Article in Journal
Galaxy Group Ellipticity Confirms a Younger Cosmos
Previous Article in Special Issue
Constraints on the Minimally Extended Varying Speed of Light Model Using Pantheon+ Dataset
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Introduction to the Number of e-Folds in Slow-Roll Inflation

by
Alessandro Di Marco
1,2,*,
Emanuele Orazi
3 and
Gianfranco Pradisi
2,4
1
Istituto Nazionale di Astrofisica, Istituto di Astrofisica e Planetologia Spaziali (INAF-IAPS), Via Fosso del Cavaliere, 100, 00133 Rome, Italy
2
Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Roma Tor Vergata, Via Della Ricerca Scientifica 1, 00133 Roma, Italy
3
Escola de Ciência e and International Institute of Physics, Federal University of Rio Grande do Norte, Campus Universitario, Lagoa Nova 59078-970, RN, Brazil
4
Dipartimento di Fisica, Università di Roma Tor Vergata, Via Della Ricerca Scientifica 1, 00133 Roma, Italy
*
Author to whom correspondence should be addressed.
Universe 2024, 10(7), 284; https://doi.org/10.3390/universe10070284
Submission received: 29 May 2024 / Revised: 24 June 2024 / Accepted: 25 June 2024 / Published: 29 June 2024
(This article belongs to the Special Issue Universe: Feature Papers 2024—'Cosmology')

Abstract

:
In this review, a pedagogical introduction to the concepts of slow-roll inflationary universe and number of e-folds is provided. In particular, the differences between the basic notion of number of e-folds ( N e ), total number of e-folds ( N T ) and number of e-folds before the end of inflation (N) are outlined. The proper application of the number of e-folds before the end of inflation is discussed both as a time-like variable for the scalar field evolution and as a key parameter for computing inflationary predictions.

1. Introduction

The paradigm of cosmological inflation [1,2,3,4,5,6] (for specific reviews, see [7,8,9,10], and for a pedagogical wide-ranging treatment, see [11,12,13,14,15]) offers a comprehensive solution for understanding the origins of the specific initial conditions that underpin the standard Hot Big Bang (HBB) theory. These conditions encompass the flatness of three-dimensional constant-time hypersurfaces, the homogeneity and isotropy of the cosmic microwave background (CMB), the shortage of hypothetical heavy relics, and the significant observable entropy. Moreover, inflation generates adiabatic, Gaussian, and nearly scale-invariant scalar metric perturbations that are responsible both for matter inhomogeneities leading to the formation of the observable large-scale structures (galaxies, clusters of galaxies, superclusters, etc.) and for the primary (relative) temperature anisotropies (∼10−5) of the CMB [16,17,18,19,20,21,22,23,24,25]. Nevertheless, inflation gives rise to tensor perturbations, or primordial gravitational waves (GWs), which could be detectable if the inflationary energy scale is sufficiently high [26,27,28,29,30,31,32,33,34,35,36]. For complete reviews on scalar and tensor perturbations, see [37,38,39]. The inflationary phase can be realized via several mechanisms, and the simplest one is the so-called (single-field) slow-roll inflation [40,41,42] that exists in a large number of varieties [43]. The minimal version of the slow-roll scenario involves a single, homogeneous, neutral, minimally coupled, and canonically normalized (pseudo)scalar field ϕ , known as the inflaton, which typically dominates the stress–energy tensor T μ ν of a reliable local universe’s patch at (or below) the Planck scale (see [44,45,46,47,48,49,50,51,52,53,54,55,56] for details on such a slow-roll onset). On general ground, the inflaton’s scalar potential V ( ϕ ) , referred to as the inflationary potential, exhibits an almost flat region and a global vacuum. Initially, the inflaton field is misaligned from the potential minima and slowly crosses the flat region. Consequently, the potential contribution dominates over the field kinetic term, resembling the presence of a false vacuum or a transient cosmological constant and inducing an almost de Sitter expansion of the universe.
After traversing a total distance Δ ϕ [57,58,59,60,61,62] and reaching a value ϕ end , the slow-roll of the field breaks down and the accelerated phase ends. The inflationary expansion yields a universe that is approximately flat, smooth, empty, and cold, as all pre-existing energy and numerical densities are diluted away. Nevertheless, as ϕ approaches ϕ e n d , the kinetic term becomes important and the inflaton rapidly drops to the true vacuum, around which it begins to oscillate. The inflaton field must also be coupled to the degrees of freedom of the Standard Model (SM) or of some Beyond the Standard Model (BSM) extensions, enabling the transfer of energy density stored in the scalar oscillation to the SM (or BSM) component fields. This phase is called reheating. It produces the relativistic plasma of light SM (or BSM) particles and the corresponding large entropy, and triggers the radiation epoch of the Hot Big Bang (HBB) cosmology (see [63,64,65,66,67,68,69] for perturbative reheating, Refs. [70,71,72,73,74,75,76,77] for non-perturbative reheating, and [78,79,80,81,82] for reviews).
This article reviews the fundamental aspects of the standard inflationary dynamics, focusing on the concept of the number of e-folds, i.e., the quantity that counts the number of exponential expansions of the universe, whose various definitions often give rise to some confusions in the literature. In particular, Section 2 is dedicated to the basics of the inflationary dynamic evolution in time, essentially described by a non-linear second-order differential equation for the homogeneous inflaton field. Section 3 presents a derivation of the alternative Hamilton–Jacobi description of inflation, founded on a system of first-order differential equations. Section 4 introduces the hierarchies of slow-roll parameters, useful both for a reliable description of the slow-roll phase and for a simpler treatment of inflationary observables O ( ϕ ) at some desired perturbative order. Section 5 contains a discussion of the general concept of the number of e-foldings (e-folds), N e , together with a derivation of the related notions of “total inflationary number of e-folds”, N T , and of the fundamental “number of e-folds before the end of inflation”, N. Section 6 provides a detailed and exhaustive derivation of the inflaton equation for the case in which the value of N plays the role of a fundamental time-like variable. Section 7 discusses the computation of the number of e-folds before the end of inflation related to the observable cosmological scales. In Section 8, solutions ϕ ( N ) of the inflaton equations are used to compute numerical estimates of the main inflationary observables O ( ϕ ) in terms of N. In particular, this recipe is applied to some representative classes of inflationary potentials, like monomial potentials, α -attractor models and a class of non-linear Einstein–Cartan gravities, classically equivalent to models with a pseudoscalar slow-rolling inflaton. Finally, Section 9 provides conclusions and discussions.
In this manuscript, the particle natural units ħ = c = 1 and the “mostly minus” Lorentzian metric signature ( + , , , ) are used. In addition, M p = 1 / 8 π G N denotes the reduced Planck mass, where G N is the Newton’s gravitational constant.

2. Standard Slow-Roll Inflationary Dynamics

The cosmological action for the slow-roll inflation can be written by adding to the classical Einsten–Hilbert action, describing the gravitational sector, a matter term whose particle content includes the inflaton (typically a scalar field), together with other SM or BSM constituents. However, on the one hand, the inflaton is expected to dominate over the remaining matter components. On the other hand, as the inflationary evolution proceeds, the contributions due to additional matter degrees of freedom tend to become rapidly and increasingly negligible. Therefore, a good description is given by
S g μ ν , ϕ S EH + S ϕ = = d 4 x g M p 2 2 R + 1 2 μ ϕ μ ϕ V ( ϕ ) ,
where ϕ is the (minimally coupled to gravity) inflaton equipped with the (slow-roll) potential V ( ϕ ) ; R is the Ricci scalar curvature; and g is the determinant of the metric tensor g μ ν , assumed to be the one giving the line element:
d s 2 = d t 2 a ( t ) 2 d l 2 , d l 2 = d r 2 + r 2 d θ 2 + sin 2 θ d φ 2 .
Here, t represents the cosmic time, a ( t ) is the (dimensionless) cosmic scale factor while d l 2 is the (flat) euclidean line element of the metric describing three-dimensional spatial hypersurfaces at constant time in terms of the triplet of comoving coordinates ( r , θ , φ ) . The validity of the assumption in Equation (2) is supported by the fact that, during the slow-roll phase, the hypothethical spatial curvature contribution quickly becomes negligible. The variation of the action in Equation (1) with respect to ϕ provides the Klein–Gordon-like equation for the inflaton field, while the variation with respect to the gravitational metric field g μ ν provides the usual Einstein field equations. The metric Ansatz in Equation (2) determines a couple of differential equations for the cosmic scale factor a ( t ) or for the Hubble rate H = a ˙ / a , sourced, as usual, by the inflaton energy–momentum tensor T μ ν interpretable as that of a perfect fluid with a given energy density ρ ϕ and pressure p ϕ . Indeed, the standard Einstein–Friedmann equations come out to be
a ˙ a 2 = 1 3 M p 2 ρ ϕ
and
a ¨ a = 1 6 M p 2 ρ ϕ + 3 p ϕ ,
where energy density and pressure can be identified with
ρ ϕ = 1 2 ϕ ˙ 2 + V ( ϕ ) , p ϕ = 1 2 ϕ ˙ 2 V ( ϕ ) .
The scalar wave equation for the inflaton reads
ϕ ¨ + 3 H ϕ ˙ + d V d ϕ = 0 ,
and it is important to stress that the “Hubble friction term” containing the first-time derivative of the field is distinctive of the expansion of the universe.
The Einstein–Friedmann–Klein–Gordon (EFKG) system of Equations (3)–(6) represents the standard tool in order to describe an (hypothetical) inflationary phase. In addition, the system allows us to identify the boundary conditions providing a suited accelerated inflating evolution of the universe. In particular, Equation (4) shows that the necessary requirement for a (positive) accelerated expanding phase is the so-called first slow-roll condition [40,41,42]:
ϕ ˙ 2 V ( ϕ ) .
Moreover, the accelerated expansion must be sustained for a sufficiently long period in order to obtain a proper smooth and flat universe, in agreement with observations. This is possible only if the inertial term in the inflaton Equation (6) is negligible compared to the Hubble friction term and the potential term. In light of this, one has to impose a second slow-roll condition:
| ϕ ¨ | | 3 H ϕ ˙ | , | ϕ ¨ | | d V ( ϕ ) d ϕ | .
The two slow-roll conditions can be translated into two corresponding slow-roll parameters:
ϵ ( t ) = H ˙ H 2 , η ( t ) = ϕ ¨ H ϕ ˙ .
The first parameter can be derived just by noticing that a ¨ / a = H ˙ + H 2 , while the second term is constructed by the ratio of the inertial and friction terms of the scalar wave Equation (6). As mentioned, the smallness of the first one ( ϵ 1 ) ensures the realization of the accelerated phase, while the smallness of the second one ( | η | 1 ) ensures that the accelerated phase lasts long enough to sufficiently stretch the universe in a way that is compatible with observations. In general, the solutions of the first Friedmann equation and of the scalar field equation give a complete description of the dynamics, i.e., determine the evolution of the scale factor a ( t ) (or H ( t ) ) and of the inflaton field ϕ ( t ) . The knowledge of the precise expansion rate of the universe is surely very interesting. However, the inflaton trajectory ϕ ( t ) represents the truly fundamental quantity, so that the system of equations is often reduced to the Equation (6) coupled to
H 2 = 1 3 M p 2 1 2 ϕ ˙ 2 + V ( ϕ ) .
They can be combined in a single non-linear second-order differential equation for the inflaton field
ϕ ¨ + 3 M p 1 2 ϕ ˙ 2 + V ( ϕ ) 1 / 2 ϕ ˙ + d V d ϕ = 0 ,
promoted to a proper Cauchy problem by adding a reliable pair of initial conditions for the inflaton field and its first derivative at a given (initial) reference time t i
ϕ ( t i ) = ϕ i , ϕ ˙ ( t i ) 1 .
The time t i can be arbitrary, although the most natural choice for it will be discussed later. In principle, the solution of the Cauchy problem is obtained via numerical integration and provides an inflaton trajectory ϕ ( t ) as a function of the cosmic time. This function shows attractor-like behavior [83,84,85,86] that is compatible with the inflationary phenomena and is extremely useful because it allows us to analyze the behavior of the inflaton field both (i) during the fundamental slow-roll phase and (ii) after the slow-roll phase, namely when the inflaton moves towards the true vacuum and begins to oscillate around it. Eventually, in the oscillating phase, the inflaton equation could be extended by adding phenomenological terms devoted to describe the decay of the inflaton modes into other particles. The related non-perturbative and perturbative reheating phases are fundamental to connect the inflationary epoch to the Friedmann–HBB dynamics in a way that preserves the precise predictions of the Big Bang Nucleosynthesis [78,79,80,81,82].

3. Hamilton–Jacobi Inflationary Dynamics

The standard slow-roll dynamics is founded on a second-order scalar wave equation with the time-like variable given by the usual cosmic time t. However, this scenario is not unique and other possibilities exist. For example, an Hamilton–Jacobi treatment of the inflationary equations based on a system of first-order (non-linear) differential equations in the time-domain has been introduced in [87], while a more convenient version with the scalar field itself as the main time-like variable has been originally proposed in [88]. The starting point is the system of Friedmann equations that can be written just in terms of the Hubble rate:
H 2 = 1 3 M p 2 ρ ϕ ,
H ˙ + H 2 = 1 6 M p 2 ρ ϕ + 3 p ϕ .
Using the expressions of the energy density and pressure in Equation (5), the two Equations (13) and (14) can be rephrased in terms of the inflaton field and combined together in the single non-linear equation
ϕ ˙ 2 = 2 M p 2 H ˙ ,
which relates the time derivative of the Hubble rate to the square of the inflaton time derivative. It should be noticed that substituting H ˙ from Equation (15) inside the time derivative of Equation (13) returns the original Klein–Gordon equation of the inflaton field, thereby confirming, in a sense, the redundancy of the original system of Equations (3)–(6). In addition, Equation (15) also suggests the alternative definition
ϵ ( t ) = H ˙ H 2 = ϕ ˙ 2 2 M p 2 H 2
of the first slow-roll parameter. Moreover, one can substitute Equation (15) in Equation (14), obtaining a single non-linear first-order differential equation for the Hubble rate, directly in terms of the scalar field potential
3 M p 2 H 2 + M p 2 H ˙ = V ( ϕ ) .
Therefore, Equations (15) and (17) determine a time-dependent first-order system of differential equations [87]. One should then trade the cosmic time variable for the scalar field itself to have a measure of the Hubble variation in terms of the scalar field variation. It can be carried out by deriving the inflaton energy density ρ ϕ with respect to time [88] and applying the Klein–Gordon equation to obtain
ρ ˙ ϕ = 3 H ϕ ˙ 2 .
By dividing both members for ϕ ˙ , one obtains
d t d ϕ = 3 H ρ ϕ ,
that implicitly returns the functional relation t ( ϕ ) . By using the derivative of Equation (13) with respect to ϕ , one can finally obtain the linear equation
ϕ ˙ = 2 M p 2 H ( ϕ ) .
It should be noticed that the condition ϕ ˙ 0 is crucial in this treatment. Interestingly, the same result can also be achieved by using the simple relation
d H d t = ϕ ˙ d H d ϕ ,
that always holds in the standard single-field slow-roll regime. Equation (20) provides two different scenarios. The first corresponds to a decreasing scalar field ( ϕ ˙ < 0 ) during the slow-roll. In this case, H ( ϕ ) > 0 and, consequently, the Hubble rate also decreases with time, H ˙ ( t ) < 0 . The second is the opposite one. An increasing scalar field ( ϕ ˙ > 0 ) during the slow-roll corresponds to H ( ϕ ) < 0 , resulting in an increasing-with-time Hubble rate, H ˙ ( t ) > 0 . Notice that a small variation in the time of the field does correspond to a small variation in the Hubble rate with the field, as heuristically expected. Equation (20) can also be inserted into Equation (17), giving rise to the following additional equation in terms of the derivative with respect to ϕ :
V ( ϕ ) = 3 M p 2 H 2 ( ϕ ) 2 M p 4 H 2 ( ϕ ) ,
The first-order Hamilton–Jacobi system of Equations (20) and (22), must be decorated with a set of initial conditions for the Hubble rate and for the scalar field:
H ( ϕ i ) = H i , ϕ ( t i ) = ϕ i .
The Hamilton–Jacobi recipe furnishes an alternative window on the inflationary dynamics and has been used to study several features of the inflationary phase. The first example is due to Salopek and Bondi [89], who both discussed the analogy with the Hamilton–Jacobi formalism of classical mechanics and the non-linear evolution of long-wavelength metric fluctuations in the context of stochastic inflation. These equations are also very useful because they allow us to generate exact inflationary solutions starting from a chosen form of the H ( ϕ ) function, an aspect that has been intensively explored in [90,91,92,93,94].

4. Slow-Roll Parameters, Inflationary Perturbations, and Observables

The previously discussed Hamilton–Jacobi treatment of inflation represents a fundamental framework to introduce a complete hierarchy of field-dependent slow-roll parameters, as discussed by Liddle et al. [42]. In particular, it is possible to define the so-called Hubble slow-roll parameters (HSRPs) as
ϵ H ( ϕ ) = 2 M p 2 H ( ϕ ) H ( ϕ ) 2 ,
β H ( n ) ( ϕ ) = ( 2 M p 2 ) n ( H ) n 1 H ( n + 1 ) H n , n 1 .
In literature, some higher-order HSRPs are also indicated with different symbols, like η H ( ϕ ) = β H ( 1 ) , ξ H 2 ( ϕ ) = β H ( 2 ) , σ H 3 ( ϕ ) = β H ( 3 ) , and so on. It is important to note that the choice ϕ ˙ < 0 naturally implies ϵ H ( ϕ ) > 0 (see Equation (20)). These parameters describe the properties of the inflationary dynamics, while their hierarchical smallness contains the initial conditions required to obtain a slow-roll inflationary phase. In particular, ϵ H ( ϕ ) 1 and | η H ( ϕ ) | 1 correctly reproduce the slow-roll conditions of Section 2. This is not a surprise since, for instance, the ϵ H parameter in Equation (24) can be obtained by simply combining Equation (16) with Equation (20). Needless to say, one can perform the same in terms of the potential, and indeed, one can define a second hierarchy, known as potential slow-roll parameters (PRSPs), as
ϵ V ( ϕ ) = M p 2 2 V ( ϕ ) V ( ϕ ) 2 ,
β V ( n ) ( ϕ ) = ( M p 2 ) n ( V ) n 1 V ( n + 1 ) V n , n 1 .
Again, the higher-order PSRPs are also denoted as η V ( ϕ ) = β ( 1 ) , ξ V 2 ( ϕ ) = β V ( 2 ) , σ V 3 ( ϕ ) = β V ( 3 ) , and so on. However, this set of parameters does not include the necessary conditions for having a slow-roll phase. Rather, it only classifies the flatness of the potential. In principle, the PSRPs are equivalent to the HSRPs if one were to add a constraint, i.e., the condition that the trajectory ϕ ( t ) approaches an attractor solution [42], corresponding to an asymptotic behavior given by the approximation of the scalar wave equation ϕ V / ( 3 H ) . However, the HSRP and PSRP hierarchies are closely correlated rather than independent. Indeed, the Hubble parameters can always be expanded (up to an arbitrary order) in series of the potential parameters. For example, the second-order expansions of the first three HSRPs are as follows:
ϵ H = ϵ V 4 3 ϵ V 2 + 2 3 ϵ V η V + o ( β V 3 ) ,
η H = η V ϵ V + 8 3 ϵ V 2 + 1 3 η V 2 8 3 ϵ V η V + o ( β V 3 ) ,
ξ H 2 = ξ V 2 3 ϵ V η V + 3 ϵ V 2 + o ( β V 3 ) .
Vice versa, the PSRPs admit an exact expression in terms of the HSRPs (see [42] for details), as one can deduce from Equation (22). If needed, one can as well expand the PSRPs in a series of HSRPs. The inflationary stage is notoriously characterized by a decrease in the comoving Hubble horizon (or radius) 1 / a H and by an ubiquitous production of vacuum fluctuations on all comoving subhorizon scales k a H of scalar fields with an effective mass smaller than the Hubble rate. In the simplest case, this involves only the inflaton field driving the expansion. These inflaton fluctuations naturally turn into metric fluctuations around the homogeneous FLRW spacetime [38], resulting in a metric tensor of the form
g μ ν g μ ν FLRW + δ g μ ν ( 1 ) + δ g μ ν ( 2 ) +
In light of this, one can single out three main types of fluctuation modes, i.e., scalar (S), vector (V), and tensor (T) modes. In the first approximation, one can consider small fluctuations around the FLRW metric and simply take into account the linear contributions of the expansion in Equation (31) [38]. The first-order approach provides simple field equations for the three types of modes in the space of comoving k. Such modes result in a complete decoupling from each other at a given scale k. In particular, the solution of the field equations suggests that the vector modes can be neglected, while scalar and tensor modes are naturally stretched beyond the (slowly changing) inflationary Hubble horizon by the accelerated expansion, quickly freezing to superhorizon classical metric perturbations k a H at a conserved value. The related period is usually called horizon crossing and can be referred to any perturbed arbitrary scale k. However, for the sake of convenience, it is better to reserve such a definition for the indication of the specific crossing of an observable perturbed scale, or the period in which the entire package of (human) observable’ perturbed scales are stretched out.
The superhorizon solutions for the amplitude of scalar and tensor modes can be used to construct the so-called power spectra of perturbations, that can be written as
P S ( k ) 1 8 π 2 M p 2 H 2 ϵ H | k = a H ; P T ( k ) 2 π 2 H 2 M p 2 | k = a H ,
where only P S happens to depend upon ϵ H , while P T does not. It is important to stress that the solutions of Equations (32) do not represent the power spectra at horizon crossing. Indeed, they are the asymptotic solutions only written in terms of the horizon crossing phase k a H (for details, see [95]). Moreover, since the evolution of the mode(s) k depends on the scalar field, Equations (32) can be thought of as functions of ϕ at horizon crossing, via H and ϵ H . The solutions of Equations (32) also allow us to conclude that single-field slow-roll inflation generates almost-scale-invariant metric perturbations with a corresponding Gaussian probability density function (PDF) [38].
The power spectra also allow us to derive the so-called inflationary observables, O ( ϕ ) , that inherit the field dependence [95]
n S ( ϕ ) 1 4 ϵ H ( ϕ ) + 2 η H ( ϕ ) ,
n T ( ϕ ) 2 ϵ H ( ϕ ) , r ( ϕ ) 16 ϵ H ( ϕ ) ,
α S ( ϕ ) 2 ξ H 2 ( ϕ ) + 10 ϵ H ( ϕ ) η H ( ϕ ) 8 ϵ H 2 ( ϕ ) ,
α T ( ϕ ) 4 ϵ H ( ϕ ) ( ϵ H ( ϕ ) η H ( ϕ ) ) ,
where n S is the scalar tilt, n T is the tensor tilt, r is the tensor-to-scalar ratio, α S is the running of the scalar tilt, and α T represents the running of tensors. Currently, the parameter n T holds more theoretical than practical interest due to the lack of a (direct or indirect) detection of an inflationary gravitational wave background. These functions describe how the value of the power spectra in Equation (32) changes with (the log of) k [95]. Typically, the inflationary stage predicts tiny inflationary observables confirming the almost-scale-invariant behavior in the k space for the power spectra of Equation (32). The specific dependency of the inflationary observables from the scalar field value ϕ is not simple to derive, because it needs the knowledge of the form of the Hubble rate H ( ϕ ) , which is strongly model-dependent. Fortunately, it is possible to use the expansions of the HSRPs in terms of the PSRPs to convert the H ( ϕ ) dependency to a potential one, obtaining [95]:
n S ( ϕ ) 1 6 ϵ V ( ϕ ) + 2 η V ( ϕ ) ,
n T ( ϕ ) 2 ϵ V ( ϕ ) , r ( ϕ ) 16 ϵ V ( ϕ ) ,
α S ( ϕ ) 2 ξ V 2 ( ϕ ) + 16 ϵ V ( ϕ ) η V ( ϕ ) 24 ϵ V 2 ( ϕ ) ,
α T ( ϕ ) 4 ϵ V ( ϕ ) ( 2 ϵ V ( ϕ ) η V ( ϕ ) ) .
In this way, the inflationary observables are explicitly determined by (the flatness of) the scalar potential, which are functions of ϕ through V ( ϕ ) . The cosmological observations reveal that the amplitude of the scalar perturbations is of order P S 2 × 10 9 [96], the scalar spectral index, measured at the reference comoving scale k 0.002 Mpc 1 , is n S = 0.9649 ± 0.0042 at 68 % of confidence level (CL) [97], while the BICEP/Keck upper limit on the tensor-to-scalar ratio at the reference comoving scale k 0.05 Mpc 1 , results r < 0.036 at 95 % CL [98].
The observables above enable testing a single-field slow-roll inflation model with potential V ( ϕ ) against current observations by varying the field as a parameter. This can be performed by selecting a reasonable sample of scalar field values from the potential plateau. However, this method is approximate, as the chosen field values may not correspond to those at horizon crossing for the relevant scale(s) k. A reliable inflaton trajectory is needed, which requires solving the inflaton equation of motion (11) with proper initial conditions (Equation (12)) for (the horizon crossing of) the scales of interest, such as 0.002 Mpc 1 .
This task is rather delicate because it is not easy to select (specific) time(s) and durations of the desired horizon crossing. Therefore, the introduction of a new time-like variable is necessary. The following three sections will be devoted to presenting this variable, exploring its utility, and discussing its natural values in the context of horizon crossing.

5. Definitions of the Number of e -Folds

In inflationary cosmology, the evolution of the inflating universe’s patch can be described, to a first approximation, by an explicit exponential term. Indeed, by definition, the Hubble rate expression can be interpreted as an ordinary differential equation once equipped with an initial condition for the cosmic scale factor at a generic slow-roll trigger time t i (see Equation (12)):
a ˙ a = H ( t ) , a ( t i ) = a i ,
whose solution obviously reads
a ( t ) = a i e t i t d τ H ( τ ) ( t > t i ) .
Imposing the first slow-roll condition in Equation (7), the Hubble function depends only on the value acquired by the inflaton (potential) energy density during the inflationary epoch
H 2 ( ϕ ) 1 3 M p 2 ρ ϕ , ρ ϕ V ( ϕ ) ,
corresponding to a “transient” ideal fluid with an equation of state (EoS) w ϕ = ρ ϕ / p ϕ 1 . In general, the time integral of the Hubble rate is called the number of e-folds of the expansion from t i to the generic cosmic time t. It measures exactly the rate of exponential growing of the universe during inflation, and can be considered as a function of the upper extreme:
N e ( t ) = t i t d τ H ( τ ) , d N e d t = H ( t ) .
If one assumes V ( ϕ ) as an almost constant plateau (for a reasonable range of the scalar field values) of height M inf 4 , the result of the integration is a linear function of time H inf ( t t i ) , with H inf M inf 2 / M p . Generically, the integral is more complicated, but still, N ( t ) is an increasing function of time, with a positive derivative also causing the cosmic scale factor a ( t ) to monotonically increase. Equation (44) allows us to write down the total number of e-folds of an inflationary period from t i to a final time t end as
N T = t i t end d τ H ( τ ) .
If single-field slow-roll inflation is assumed to start at a time t i immediately following a primordial quantum gravity phase and thought to propagate below an appropriate BSM phase (for instance, a GUT phase), the total number of exponential expansions could amount to millions—an information that may be indicative of the real size of the causal-connected universe, the so-called particle horizon [11,12,13,14]. However, it is commonly believed that the observable perturbed scales (and in some sense, the current observable universe) has been generated at certain times during the last stages of the inflationary expansion. In light of this, one can introduce a new convenient variable, the number of e-folds before the end of inflation
Δ N ( t , t e n d ) = N T N e ( t ) = t i t e n d d τ H ( τ ) t i t d τ H ( τ ) = t t e n d d τ H ( τ ) ,
obtained by subtracting from N T the number of e-folds that precede the beginning of the “interesting phase”, where the observable perturbed scales are generated. This variable is decreasing in time, since
d Δ N ( t , t e n d ) d t = d N e d t < 0 ,
and can be used as a backwords time-like variable, vanishing when t = t end .
The number of e-folds before the end of inflation is often rewritten in terms of the inflaton field ϕ along the flat potential direction and of the first HSRP ϵ H ( ϕ ) . In particular, formally inverting the time profile of ϕ , one can turn the time integral into a field-space integral, obtaining
Δ N ( ϕ , ϕ e n d ) = ϕ ϕ e n d d ϕ ˜ H ( ϕ ˜ ) ϕ ˜ ˙ .
Using Equation (20), assuming ϕ ˙ < 0 and inserting the parameter ϵ H ( ϕ ) , one finally obtains
Δ N ( ϕ , ϕ e n d ) = 1 M p ϕ e n d ϕ d ϕ ˜ 1 2 ϵ H ( ϕ ˜ ) .
In the following, we define
N ( ϕ ) Δ N ( ϕ , ϕ end ) .
The reader should be aware that many authors label the same quantity with other notations, e.g., N e or N .
It should be noticed that Equation (49), under the approximation ϵ H ( ϕ ) ϵ V ( ϕ ) , is equivalent to the inverse solution of the first-order differential equation
d ϕ d N M p 2 ϵ V ( ϕ )
that provides the scalar field profile as a function of N, ϕ ( N ) . As aforementioned, such field trajectory is not an exact solution, as it is rather the slow-roll limit of a more general equation, and is to be discussed in the next section.

6. Inflaton Evolution and Number of e-Folds

In Section 2, the system of Einstein–Friedmann–Klein–Gordon equations was introduced, which provides the function a ( t ) and the global inflaton behavior ϕ ( t ) . Nevertheless, as anticipated at the end of Section 4, the knowledge of the time profile of the field ϕ is not very useful for the numerical estimate of the inflationary observables O ( ϕ ) . These quantities are closely related to the horizon crossing window, so it would be advisable to compute them for a proper range of the number of e-folds N. In this respect, it would be much better to obtain inflaton solutions as explicit functions of the number of e-folds
Observable   O [ ϕ ( t ) ] Horizon   Crossing N O [ ϕ ( N ) ] ,
in a such a way that it is possible to literally sample O ( ϕ ) in terms of N. To this aim, it is useful to rewrite the Klein–Gordon equation with N (in place of the cosmic time) as the basic time-like variable. Its derivation, often overlooked in the current literature (see, however, [43] for a brief description) is now reviewed in detail.
The relation between time and N is contained in the equation N ˙ = H , which results from the combination of Equation (44) and Equation (46).
Given a generic function whose time dependence f ( t ) = f [ N ( t ) ] is through N, it is straightforward to derive the expressions for the first and second derivatives in the form
d f d t = H d f d N ,
d 2 f d t 2 = H 2 ( N ) ϵ ( N ) d f d N + H 2 ( N ) d 2 f d N 2 ,
where the relation
ϵ ( N ) = ϵ ( t ( N ) ) = H ˙ ( t ) H 2 = 1 H d H d N
has been used to insert the first slow-roll parameter ϵ ( N ) . Applying the above relations to ϕ ( N ) inside the Klein–Gordon equation, one obtains
H 2 ( N ) ϵ ( N ) d ϕ d N + H 2 ( N ) d 2 ϕ d N 2 3 H 2 ( N ) d ϕ d N + d V d ϕ = 0 ,
or
d 2 ϕ d N 2 + ϵ ( N ) d ϕ d N 3 d ϕ d N + 1 H 2 ( N ) d V d ϕ = 0 .
It is useful to highlight that Equation (3) gives back H as a function of N in the form
H 2 ( N ) = 2 V ( ϕ ) 6 M p 2 d ϕ d N 2 ,
while Equation (15) in terms of N provides
d ϕ d N 2 = 2 M p 2 ϵ ( N ) .
As a result, the Hubble rate can be written as
H 2 ( N ) = 1 M p 2 V ( ϕ ) 3 ϵ ( N )
in terms of the potential and the “new” slow-roll parameter ϵ ( N ) , introduced in Equation (55). Performing an inversion of sign N N in order to allow for a forward integration over the time-like variable, the Cauchy problem related to the Klein–Gordon equation finally reads
d 2 ϕ d N 2 + ( 3 ϵ ( N ) ) d ϕ d N + ( 3 ϵ ( N ) ) M p 2 ϵ V ( ϕ ) = 0 , ϕ ( N i ) = ϕ i , ϕ ˙ i 1 .
The system is a second-order differential problem. It is non-linear because ϵ ( N ) contains a derivative squared. Moreover, the first derivative term plays the same role as the Hubble friction in the original formulation. The Klein–Gordon equation needs to be numerically integrated to provide the exact inflaton trajectory in terms of the e-folds number. However, the slow-roll limit corresponding to the conditions
ϵ ( N ) 1 , d 2 ϕ d N 2 d ϕ d N , M p 2 ϵ V ( ϕ ) ,
can be used to reduce the equation to the linear one
d ϕ d N + M p 2 ϵ V ( ϕ ) 0 ,
that simply matches the equation met in the previous section (with N N ). Of course, the integration of the linear equation correctly reproduces the (first-order) definition of the number of e-folds before the end of inflation
N ( ϕ ) Δ N ( ϕ , ϕ e n d ) 1 M p ϕ e n d ϕ d ϕ ˜ 1 2 ϵ V ( ϕ ˜ ) ,
where the inflaton value ϕ end is computed from the corresponding slow-roll breaking condition
ϵ V ( ϕ e n d ) 1 .
This treatment suggests that, given a scalar potential V ( ϕ ) , an approximate analytic solution ϕ ( N ) can be simply obtained by inverting Equation (64), if possible. Such a solution allows us to choose an initial N and then to compute a related initial field value for the Cauchy problem in Equation (61).
The numerical integration and the related field function should be given over a reliable range of N associated with the horizon crossing of the reference-observable scales. In principle, nothing prevents assuming arbitrary initial conditions for the slow-rolling field. However, it is certainly better to consider simplified conditions already compatible with the horizon crossing of the observable perturbed scales. In the next section, we present a well-codified strategy to select the value of N related exactly to the epoch in which the relevant observable scales cross the inflationary Hubble horizon.

7. Estimates of the Number of e -Folds before the End of Inflation

The number of e-folds occurring before the end inflation, N, can be estimated in a very simple way in terms of generic properties of the inflationary phase and the features qualifying the postinflationary evolution. In this context, it is interesting to consider a comoving scale k that crosses the Hubble horizon at a certain time t during inflation, k a H , and the comoving scale corresponding to the size of the current Hubble horizon, k 0 a 0 H 0 . The ratio of the two comoving scales
a H a 0 H 0 = k k 0 = λ 0 λ
is naturally constant throughout the evolution of the universe. This ratio can also be split as
k a 0 H 0 = a a end a end a reh a reh a 0 H H 0 ,
where the first partial ratios are referred as the inflationary evolution, the reheating epoch, and the HBB epoch, respectively, while the last fraction represents the ratio between the Hubble scale at the inflationary time t and the current one. In order to obtain the number of e-folds, it is convenient to take the logarithm of the ratio, obtaining
log k a 0 H 0 = log a a end + log a end a reh + log a reh a 0 + log H H 0 = N N reh N HBB + log H H 0 .
Here, as said, N represents the number of e-folds occurring between the generic time t and the end of inflation, t end , while N reh describes the number of e-folds during the reheating phase and N HBB indicates the e-folds related to the standard Big Bang phase. One can thus write the following:
N = log k a 0 H 0 N reh N HBB + log H H 0
Notice that the first contribution is just a numerical term that depends on the chosen comoving scale k. N reh , on the other hand, can be estimated by assuming a global and effective description of the reheating phase. In particular, one can suppose a reheating phase characterized by a single, global fluid with an average equation of state (EoS) w reh and a conserved average energy density described by a proper Boltzmann equation:
ρ ˙ reh + 3 H ( 1 + w reh ) 0 .
The solution of Equation (70) between t end and t reh allows us to express N reh in terms of the average w reh , the energy density at the end of slow-roll inflation ρ end , and the energy density at the end of the reheating ρ reh :
N reh 1 3 ( 1 + w reh ) log ρ end ρ reh .
The third contribution in Equation (69) can be easily extracted from the fundamental definition of an HBB cosmology, i.e., the conservation of the comoving entropy density:
g S ( t reh ) a reh 3 T reh 3 = g S ( t 0 ) a 0 3 T 0 3 ,
where g S is the number of the relativistic degrees of freedom contributing to the entropy density, T reh is the temperature of the relativistic matter close to the reheating scale, and T 0 is the current temperature of the relativistic component of the universe, mainly CMB photons. Remembering that
ρ reh π 2 30 g E ( t reh ) T reh 4 ,
where g S is the number of the relativistic degrees of freedom contributing to the energy density, and by assuming that g S ( t reh ) g E ( t reh ) g reh , one can infer that
N HBB log T 0 1 4 log ρ reh 1 12 ln g r e h + log 43 11 1 3 π 2 30 1 4 .
The substitution of the relations Equation (71), Equation (74), and Equation (43) into Equation (69) leads to
N ζ 1 3 ( 1 + w reh ) ln ρ end ρ reh + 1 4 ln V 2 M p 4 ρ reh
where
ζ = ln k a 0 H 0 + ln T 0 H 0 + ζ 0
and ζ 0 is just a numerical factor
ζ 0 = 1 12 ln g r e h + 1 4 ln 1 9 + ln 43 11 1 3 π 2 30 1 4 .
Equation (75) describes the number of e-folds related to the horizon exit of the mode k in terms of both the properties of the inflationary phase ( V ( ϕ ) , ρ end ) and the properties of the postinflationary phase, relative to reheating and HBB cosmology epochs. In particular, considering for simplicity an almost instantaneous reheating phase and neglecting the subdominant numerical term ζ 0 , one can obtain the following simplified expression for the number of e-folds:
N ln k a 0 H 0 + ln T 0 H 0 + ln M inf M p .
It is important to emphasize that N depends on the chosen “pivot” comoving wave number k. It is worth mentioning that a smaller comoving scale k < a 0 H 0 would obviously give a smaller value of N, reflecting the fact that the smaller the comoving scale crossing the horizon, the lower the number of e-folds before the end of inflation. It is customary to focus on perturbed scales that can be experimentally probed, like that of the Planck mission, k 0.002 Mpc 1 . Remembering that a 0 = 1 , H 0 1.75 × 10 42 GeV and T 0 = 2.3 × 10 13 GeV, one obtains
N 64 + ln M inf M p .
In Equation (79), N has a mild dependency on the inflationary scale. It is indeed possible to demonstrate that it lies in the range 60–50 [99,100] for M inf ( 10 16 10 13 ) GeV. A more complicated postinflationary phase (see, for example, [101,102,103,104]) would induce modifications to the standard expression of Equation (75).
In the next section, the first-order computation of inflationary observables in terms of the variable N for several paradigmatic classes of inflationary models is presented. In addition, numerical predictions for the commonly chosen value of N 60 and their compatibility with current experimental constraints are also discussed.

8. Inflationary Observables and Number of e -Folds

In the realm of the single-field slow-roll inflationary models, the simplest example is surely given by the class of monomial potentials. This class of potentials was firstly introduced contextually to the chaotic inflation scenario [6] and emerge, for instance, in a number of BSM-motivated cosmologies like those within supergravity [105,106,107,108,109,110] or in superstring theories (see, e.g., axion monodromy inflation [111,112,113]). The form of the potential is
V ( ϕ ) = M inf 4 ϕ M p n ,
which is often also written as
V ( ϕ ) = λ n ϕ n , λ n = M inf 4 M p n ,
where the parameter λ n is essentially the coupling constant of the inflaton self-interaction. This class of models falls within the larger class of large-field inflation since the vacuum expectation value of the inflaton field reduces from super-Planckian values ϕ M p to the vacuum value ϕ = 0 with ϕ ˙ < 0 . The most known scenarios are the quadratic case
V ( ϕ ) = 1 2 m ϕ 2 ϕ 2 , m ϕ 2 = 2 M inf 4 M p 2 ,
and the quartic case
V ( ϕ ) = λ ϕ 4 ϕ 4 , λ ϕ = 4 M inf 4 M p 4 ,
originally discussed by Andrei Linde as the simplest examples of chaotic inflation [6]. The k-th derivative of the scalar potential in Equation (81)
V ( k ) ( ϕ ) = M inf 4 M p k Γ ( n + 1 ) Γ ( n k + 1 ) ϕ M p n k , k n ,
allows us to write, for instance, the first three potential parameters:
ϵ V ( ϕ ) = n 2 2 M p ϕ 2 ,
η V ( ϕ ) = n ( n 1 ) M p ϕ 2 ,
ξ V 2 ( ϕ ) = n 2 ( n 1 ) ( n 2 ) M p ϕ 4 .
Thus, one can derive the expression of the scalar power spectrum as
P S ( ϕ ) 1 12 π 2 n 2 M inf M p 4 ϕ M p n + 2
as well as the inflationary observables introduced in Section 4 as explicit functions of the scalar field:
n S ( ϕ ) 1 n ( n + 2 ) M p ϕ 2 ,
r ( ϕ ) 8 n 2 M p ϕ 2 ,
α S ( ϕ ) 2 n 2 ( n + 2 ) M p ϕ 4 ,
α T ( ϕ ) 2 n 3 M p ϕ 4 .
The condition ϕ ˙ < 0 naturally suggests how the scalar spectral index tends to decrease while the tensor-to-scalar ratio has an opposite behavior. As discussed in Section 4, one can roughly deduce the value of the cosmological parameters by considering a reasonable set of inflaton configurations in the slow-roll plateau region of the scalar potential V ( ϕ ) . Nevertheless, a preferred approach is to determine the inflaton values by solving the inflaton equation of motion (61), derived to this aim in Section 6. The associated numerical solution gives an inflaton trajectory ϕ ( N ) that can be used to obtain a final, numerical evolution of the interesting observables with respect to N. Finally, first-order analytic estimates of the cosmological parameters in terms of N can be very useful as well. To achieve them, one can consider the slow-roll limit of the exact Equation (64) for N
N ( ϕ ) 1 n M p 2 1 2 ϕ 2 ϕ e n d 2 ,
and try to invert it explicitly. The result is
ϕ M p 2 2 n N + ϕ e n d M p 2 ,
where the final field value, obtained by the condition ϵ V ( ϕ e n d ) 1 , comes out to be
ϕ e n d M p 2 n 2 2
and leads to
ϕ M p 2 2 n N + n 2 2 .
The fundamental result in Equation (96) provides the expression
P S 1 12 π 2 n 2 M inf 4 M p 4 2 n N + n 2 / 2 1 + n / 2
for the scalar spectrum. It can be inverted in terms of the ratio M inf / M p and, if one applies the current estimate for the scalar spectrum amplitude of the COBE mission, P S cobe 2 × 10 9 , it is possible to constrain the inflationary scale for a given n as
M inf M p 12 π 2 n 2 P S cobe 2 n N + n 2 / 2 1 + n / 2 1 / 4 .
In principle, one can assume that inflation occurs at some specific scale, for example, around 10 15 GeV. This statement, together with the COBE outcome, offers the possibility to constrain the model parameter n. At the same time, the expressions for the set of inflationary parameters in terms of N are
n S ( N ) 1 n ( n + 2 ) 2 n N + n 2 / 2 ,
r ( N ) 8 n 2 2 n N + n 2 / 2 ,
α S 2 n 2 ( n + 2 ) ( 2 n N + n 2 / 2 ) 2 ,
α T 2 n 3 ( 2 n N + n 2 / 2 ) 2 .
To give an example, selecting the massive case ( n = 2 ) that one has at first order leads to
n S ( N ) 1 2 N
r ( N ) 8 N
α S 2 N 2
α T 1 N 2 .
If one assumes that the exits of reference scales correspond to N 60 , then the related inflaton value is of the order of ϕ 15 M p , the predicted reference energy scale is M inf 2 × 10 15 GeV, and, as a consequence, the inflaton mass is estimated to be m ϕ 10 12 GeV. The numerical predictions for the two main inflationary observables result in
N 60 , n S 0.9667 r 0.1333
N 50 , n S 0.9600 r 0.1600
It is worth to stress that the numerical solution of the complete Equation (61) provides slightly different (and clearly more precise) results. For example, for N 60 , one obtains n S 0.9669 and r 0.1322 , while for N 50 , n S 0.9606 and r 0.1573 .
The same recipe can be used to evaluate other scenarios, as, for example, the model with n = 4 . In the single case of N 60 , the scalar field value of the quartic case is ϕ 20 M p and the related energy scale is M inf 10 14 GeV, with a very small self-coupling constant λ ϕ 10 15 10 16 . Unfortunately, both the quadratic and quartic cases are not compatible with the current Planck and BICEP results.
The monomial scenario provides a very simple function N ( ϕ ) that can be easily inverted to give the solution ϕ ( N ) . Moreover, a similar discussion can be performed for the exponential scenario (see [114,115,116]), where an exact solution for ϕ ( N ) and a ( t ) can also be obtained. Unfortunately, this is not always the case. For instance, there are models in which the inversion must proceed via some special functions, preventing a simple analytical interpretation for the solution ϕ ( N ) . An interesting example is given by the class of the so-called α -attractor models. They can be generated in several different ways, although the most advanced version emerges from supergravity [117,118,119,120,121,122,123,124]). The α -attractor potentials are typically divided into two subclasses called E models and T models, both characterized by a flat region (for ϕ > M p ) that comes out to be protected against quantum corrections [124]. In the present review, the focus will be limited only to the E-models, whose scalar potential is
V ( ϕ ) = M inf 4 1 e b ϕ / M p 2 n , b = 2 3 α ,
with the global minimum located at ϕ = 0 . The scalar potential is characterized by a couple of free parameters, n and α . In supergravity, the α parameter is related to the Kähler curvature of the moduli space manifold, where the inflaton superfield is a coordinate [123]. It is important to stress how this family of potentials naturally reproduces well-known inflationary models. For example, the case n = 1 corresponds to the Einstein frame version of the Starobinsky model of inflation [125], while the case α = 1 / 9 is the supergravity Linde–Goncharov model [126,127]. It should also be noticed that for large values of b (or α 1 ), the exponential term tends to be suppressed and the scalar function approximates a step-behavior. On the contrary, for small values of b (or α 1 ), V ( ϕ ) tends to assume a simple quadratic parabolic shape. In order to simplify the treatment, one can simply consider the scenario with n = 1 . In this case, the derivatives of the potential
V ( k ) ( ϕ ) = M inf 4 M p k ( 1 ) k 1 2 b k e b ϕ / M p 1 2 k 1 e b ϕ / M p
allow us to write down the PRSPs. In particular, the first two come out to be
ϵ V ( ϕ ) = 2 b 2 e b ϕ / M p 1 2 ,
η V ( ϕ ) = 2 b 2 e b ϕ / M p 2 e b ϕ / M p 1 2 .
As a consequence, one obtains
P S ( ϕ ) 1 48 π 2 M inf M p 4 e 2 b ϕ / M p b 2 1 e b ϕ / M p 4
and
n S ( ϕ ) 1 4 b 2 e b ϕ / M p + 1 e b ϕ / M p 1 2 ,
r ( ϕ ) 32 b 2 e b ϕ / M p 1 2 .
The solution of Equation (64) gives
N ( ϕ ) 1 2 b 2 e b ϕ / M p e b ϕ e n d / M p 1 2 b ϕ M p ϕ e n d M p ,
where ϕ end , defined by ϵ V ( ϕ e n d ) 1 , results in
ϕ e n d M p 1 b log ( 1 + 2 b ) .
The inversion of the expression in Equation (116) is manifestly not a simple task. One can try to use to the Lambert function W [43] to obtain
ϕ ( N ) M p = 1 b F ( N ) W ( e F ( N ) ) ,
where
F ( N ) = 2 b 2 N + e b ϕ e n d / M p b ϕ end M p .
Unfortunately, this result does not provide a straightforward interpretation of the field evolution. It is thus convenient to adopt an additional approximation, assuming that the plateau inflaton value is larger than the final field value, ϕ ϕ end . It implies that
e b ϕ / M p 2 b 2 N , ϕ ( N ) 1 b ln 2 b 2 N ,
giving rise to the scalar power spectrum
P S N 2 18 α π 2 M inf M p 4 ,
whose inversion for the inflation scale provides
M inf M p 18 α π 2 N 2 P S c o b e .
On the other hand, the inflationary observables decrease with N as
n S ( N ) 1 4 b 2 e b ϕ ( N ) / M p + = 1 2 N ,
r ( N ) 32 b 2 e 2 b ϕ ( N ) / M p + = 12 α N 2 .
The simplest example to discuss is the Starobinsky model corresponding, as mentioned, to n = 1 . Taking the standard horizon crossing option N 60 , one has a scalar field at horizon crossing ϕ 5 M p , an inflationary reference scale M inf 3 × 10 15 GeV, and numerical estimates for scalar tilt and tensor-to-scalar ratio:
N 60 , n S 0.9667 , r 0.0033
N 50 , n S 0.9600 , r 0.0048 .
The numerical investigation would reveal n S 0.9653 and r 0.0034 for N 60 and n S 0.9584 with r 0.0049 for N 50 . Therefore, the Starobinsky case as well as all the α -attractor cases with a small model parameter (i.e., α < 1 ) are extremely compatible with the current observational constraints.
To conclude this very brief overview, it can be interesting to analyze models where the integral giving N ( ϕ ) can be exactly performed, but the resulting expression cannot be explicitly inverted, even resorting to special functions. In this context, one can discuss a recently discovered class of models, based on Einstein-Cartan extension of general relativity [128,129,130]. In particular, it has been shown that one may consider a gravity theory whose action is a combination of the scalar curvature term and a non-linear function of the so-called parity violating Holst invariant term:
R ε μ ν ρ σ R μ ν ρ σ ,
where R μ ν ρ σ is the curvature tensor and ε μ ν ρ σ is the Levi–Civita invariant tensor. R obviously vanishes in a Riemannian geometry, but in a Palatini approach, it provides a dynamical torsion and results classically equivalent (on shell) to a pseudoscalar field, minimally coupled to general relativity, and equipped with a non-trivial potential suitable to drive a single-field slow-roll inflationary phase. Indeed, if ϕ indicates the pseudoscalar field unequivocally emerging from the underlying non-Riemannian geometry, the induced potential comes out to be
V ( ϕ ) = M inf 4 | γ sinh X ( ϕ ) 1 | p p 1 ,
with
X ( ϕ ) = 2 3 ϕ M p + θ γ
and the reference scale of inflation is given by
M inf 4 = p 1 p p / ( p 1 ) 1 ξ 1 p 1 | M p 2 4 γ | p p 1 .
In the previous expressions (for details, see [130] and refeences therein), γ is the so-called Barbero–Immirzi parameter, p 1 is a real parameter, ξ is the coupling of the non-linear term to the metric, and
θ γ = sinh 1 ( γ 1 ) .
The sign of the Barbero–Immirzi parameter determines the direction of the slow-roll phase. Specifically, the slow-roll phase occurs for decreasing values of the inflaton field (i.e., ϕ ˙ < 0 ) for negative values of γ , while it occurs for increasing values of ϕ (i.e., ϕ ˙ > 0 ) for positive values of γ . The first derivative of the scalar potential reads
V ( ϕ ) = M inf 4 M p 2 γ p cosh X ( ϕ ) γ sinh X ( ϕ ) 1 p p 1 3 p 1 γ sinh X ( ϕ ) 1
and the second derivative is
V ( ϕ ) = M inf 4 M p 2 2 γ p γ sinh X ( ϕ ) 1 p p 1 3 ( p 1 ) γ sinh X ( ϕ ) 1 2 γ p 1 ( p sinh 2 X ( ϕ ) + 1 ) sinh X ( ϕ ) .
Thus, the firsts two PSRPs can be calculated and result in
ϵ V ( ϕ ) = γ 2 3 p p 1 2 cosh 2 X ( ϕ ) γ sinh X ( ϕ ) 1 2
and
η V ( ϕ ) = 2 γ p 3 ( p 1 ) γ sinh X ( ϕ ) 1 2 γ p 1 ( p sinh 2 X ( ϕ ) + 1 ) sinh X ( ϕ ) .
The power spectrum assumes the very complicated expression
P S ( ϕ ) M inf 4 18 π 2 M p 4 γ 2 p 1 p 2 γ sinh X ( ϕ ) 1 2 cosh 2 X ( ϕ ) | γ sinh X ( ϕ ) 1 | p p 1 ,
while the two main observables turn out to be
n S ( ϕ ) 1 2 γ 2 3 p p 1 2 cosh 2 X ( ϕ ) γ sinh X ( ϕ ) 1 2 + 4 γ p 3 ( p 1 ) γ sinh X ( ϕ ) 1 2 γ p 1 p sinh 2 X ( ϕ , 1 sinh X ( ϕ )
and
r ( ϕ ) 16 3 γ 2 p p 1 2 cosh 2 X ( ϕ ) γ sinh X ( ϕ ) 1 2 .
As usual, to obtain first-order expressions of the observables and an estimate of the model parameter M inf , the solution of the integral N ( ϕ ) is required. One obtains
N ( ϕ ) = 3 ( p 1 ) 2 p ln | cosh X ( ϕ ) | 1 γ tan 1 sinh X ( ϕ ) | ϕ e n d ϕ = 3 ( p 1 ) 2 p ln | cosh X ( ϕ ) | | cosh X ( ϕ e n d ) | 3 ( p 1 ) 2 γ p tan 1 sinh X ( ϕ ) tan 1 sinh X ( ϕ e n d ) ,
which neatly show how there is no hope to analytically determine ϕ ( N ) by inversion, unless a very specific range of the Barbero–Immirzi parameter ( | γ | 10 1 ) is considered. However, the function in Equation (139) is still very useful to extract possible pairs ( N , ϕ ) , allowing us to derive simple numerical predictions for the cosmological observables O ( ϕ ) . For instance, in the case p = 2 , choosing γ 10 3 and N 60 , one obtains
n S 0.9680 , r 0.003 ,
with a very high inflation reference scale M inf 5.6 × 10 15 GeV. Such predictions result in a very good agreement with the current constraints.
A final observation is in order: there are cases, of course, where the potential function V ( ϕ ) does not allow for an exact calculation of the integral N ( ϕ ) . In a similar situation, it is convenient to proceed by approximating the potential with its slow-roll expression. It happens, for instance, in superstring-inspired models like the ones called fiber-inflation scenarios [131,132,133,134].

9. Summary and Conclusions

The conventional Einstein–Friedmann–Klein–Gordon system of equations offers the simplest framework for exploring the evolution of the single-field slow-roll inflationary scenario in the early universe. Moreover, the solutions for the scale factor a ( t ) and the field trajectory ϕ ( t ) have a straightforward interpretation and provide insights into the overall dynamics. However, the time-dependent profile of the inflaton solution proves to be less practical for tasks like the computation of inflationary observables that are, quite naturally, functions of the scalar field ϕ via potential slow-roll parameters (PRSPs) and typically exhibit an unusable dependence upon time. In light of this, it becomes unavoidable to reformulate the inflaton equation using an alternative and more suitable time-like variable. The number of e-folds before the end of inflation, denoted by N and tracking the times the universe grows exponentially until the end of the slow-roll phase, appears to be the most natural. Rough evaluations suggest that its value from the horizon crossing (of the observable perturbed scales) to the end of inflation typically falls within the range between 50 and 60. Actually, many more possibilities exist, also depending on the properties of the postinflationary phase (see, for example, [101,102,103,104]).
In this paper, we have reviewed the derivation and uses of the inflaton equations based on N as a fundamental time-like variable in an approach that is widely used and appreciated in the current literature, especially in order to obtain deeper insights on the inflationary dynamics. The procedure is characterized by the introduction of a new form of the first slow-roll parameter ϵ , which is slightly different from the traditional one.
Complete numerical solutions of the new equation require initial conditions that are suited to describe the horizon crossing of the relevant cosmological perturbed scales. Alternatively, one can apply the slow-roll approximation, which returns the standard definition of N and facilitates the computation of a first-order, analytical field solution ϕ ( N ) , which is extremely useful to rephrase the dependency of the standard inflationary observables in terms of N.
Finally, it is worth emphasizing another point that underscores the importance of the number of e-folds. In Section 4, the focus was on the theory of linear or first-order metric perturbations. However, as is well known, it is also important to study the behavior of non-linear curvature perturbations since they can provide additional inflation observables like higher-order spectra (bispectrum, trispectrum…), which are useful to estimate non-Gaussianities on CMB scales (see, for instance, [135,136,137,138,139,140]). In those situations, the standard computation approach proves to be quite cumbersome. More suitable non-perturbative methods have thus been introduced, such as the δ N formalism [141,142,143,144,145], based on counting the local number of e-folds of various local patches. The δ N formalism is widely used in the current literature due to its easy-to-use nature (for a recent review, see [146]).

Author Contributions

Conceptualization, A.D.M., E.O. and G.P.; methodology, A.D.M., E.O. and G.P.; software, A.D.M., E.O. and G.P.; validation, A.D.M., E.O. and G.P.; formal analysis, A.D.M., E.O. and G.P.; investigation, A.D.M., E.O. and G.P.; resources, A.D.M., E.O. and G.P.; data curation, A.D.M., E.O. and G.P.; writing—original draft preparation, A.D.M., E.O. and G.P.; writing—review and editing, A.D.M., E.O. and G.P.; visualization, A.D.M., E.O. and G.P.; supervision, A.D.M., E.O. and G.P.; project administration, A.D.M., E.O. and G.P.; funding acquisition, A.D.M., E.O. and G.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study.

Acknowledgments

We would like to thank G. De Gasperis and Y. Mambrini for interesting discussions. E.O. would like to thank the Physics Department of the University of Rome “Tor Vergata” for the kind hospitality while this work was being undertaken. E.O. was partially supported by the “Bando Visiting Professor 2023, University of Rome “Tor Vergata”, Prot. n. 0038224 (26/07/2023), by the project PROMETEO/2020/079 (Generalitat Valenciana), CNPq (Brazil), grant No.314392/2021-1 and by the Spanish Agencia Estatal de Investigacion Grant No. PID2022-138607NB-I00, funded by 14 MCIN/AEI/10.13039/501100011033, FEDER, UE, and ERDF (“A way of making Europe”).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Starobinsky, A.A. A New Type of Isotropic Cosmological Models Without Singularity. Phys. Lett. B 1980, 91, 99–102. [Google Scholar] [CrossRef]
  2. Guth, A.H. The Inflationary Universe: A Possible Solution to the Horizon and Flatness Problems. Phys. Rev. D 1981, 23, 347. [Google Scholar] [CrossRef]
  3. Linde, A.D. A New Inflationary Universe Scenario: A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole Problems. Phys. Lett. 1982, 108B, 389. [Google Scholar] [CrossRef]
  4. Albrecht, A.; Steinhardt, P.J. Cosmology for Grand Unified Theories with Radiatively Induced Symmetry Breaking. Phys. Rev. Lett. 1982, 48, 1220. [Google Scholar] [CrossRef]
  5. Hawking, S.W.; Moss, I.G. Supercooled Phase Transitions in the Very Early Universe. Phys. Lett. 1982, 110B, 35. [Google Scholar] [CrossRef]
  6. Linde, A.D. Chaotic Inflation. Phys. Lett. 1983, 129B, 177. [Google Scholar] [CrossRef]
  7. Linde, A.D. Particle physics and inflationary cosmology. Contemp. Concepts Phys. 1990, 5, 1–362. [Google Scholar]
  8. Linde, A.D. Inflationary Cosmology. Lect. Notes Phys. 2008, 738, 1–54. [Google Scholar] [CrossRef]
  9. Olive, K.A. Inflation. Phys. Rept. 1990, 190, 307. [Google Scholar] [CrossRef]
  10. Baumann, D. INFLATION. Phys. Large Small 2011, 523–686. [Google Scholar] [CrossRef]
  11. Kolb, E.W.; Turner, M.S. The Early Universe. Front. Phys. 1990, 69, 1–547. [Google Scholar] [CrossRef]
  12. Mukhanov, V. Physical Foundations of Cosmology; Cambridge University Press: Cambridge, UK, 2005; ISBN 978-0-521-56398-7. [Google Scholar] [CrossRef]
  13. Weinberg, S. Cosmology; Oxford University Press: Oxford, UK, 2008; 593p. [Google Scholar]
  14. Gorbunov, D.S.; Rubakov, V.A. Introduction to the Theory of the Early Universe: Cosmological Perturbations and Inflationary Theory; World Scientific: Hackensack, NJ, USA, 2011; 489p. [Google Scholar] [CrossRef]
  15. Baumann, D. Primordial Cosmology. PoS 2018, TASI2017, 009. [Google Scholar] [CrossRef]
  16. Mukhanov, V.F.; Chibisov, G.V. Quantum Fluctuations and a Nonsingular Universe. JETP Lett. 1981, 33, 532–535. [Google Scholar]
  17. Hawking, S.W. The Development of Irregularities in a Single Bubble Inflationary Universe. Phys. Lett. B 1982, 115, 295. [Google Scholar] [CrossRef]
  18. Starobinsky, A.A. Dynamics of Phase Transition in the New Inflationary Universe Scenario and Generation of Perturbations. Phys. Lett. B 1982, 117, 175–178. [Google Scholar] [CrossRef]
  19. Guth, A.H.; Pi, S.Y. Fluctuations in the New Inflationary Universe. Phys. Rev. Lett. 1982, 49, 1110–1113. [Google Scholar] [CrossRef]
  20. Linde, A.D. Scalar Field Fluctuations in Expanding Universe and the New Inflationary Universe Scenario. Phys. Lett. B 1982, 116, 335–339. [Google Scholar] [CrossRef]
  21. Bardeen, J.M.; Steinhardt, P.J.; Turner, M.S. Spontaneous Creation of Almost Scale—Free Density Perturbations in an Inflationary Universe. Phys. Rev. D 1983, 28, 679. [Google Scholar] [CrossRef]
  22. Lyth, D.H. Large Scale Energy Density Perturbations and Inflation. Phys. Rev. D 1985, 31, 1792–1798. [Google Scholar] [CrossRef]
  23. Mukhanov, V.F. Gravitational Instability of the Universe Filled with a Scalar Field. JETP Lett. 1985, 41, 493–496. [Google Scholar]
  24. Sasaki, M. Large Scale Quantum Fluctuations in the Inflationary Universe. Prog. Theor. Phys. 1986, 76, 1036. [Google Scholar] [CrossRef]
  25. Stewart, E.D.; Lyth, D.H. A More accurate analytic calculation of the spectrum of cosmological perturbations produced during inflation. Phys. Lett. B 1993, 302, 171–175. [Google Scholar] [CrossRef]
  26. Grishchuk, L.P. Amplification of gravitational waves in an istropic universe. Zh. Eksp. Teor. Fiz. 1974, 67, 825–838. [Google Scholar]
  27. Starobinsky, A.A. Spectrum of relict gravitational radiation and the early state of the universe. JETP Lett. 1979, 30, 682–685. [Google Scholar]
  28. Rubakov, V.A.; Sazhin, M.V.; Veryaskin, A.V. Graviton Creation in the Inflationary Universe and the Grand Unification Scale. Phys. Lett. B 1982, 115, 189–192. [Google Scholar] [CrossRef]
  29. Fabbri, R.; Pollock, M.d. The Effect of Primordially Produced Gravitons upon the Anisotropy of the Cosmological Microwave Background Radiation. Phys. Lett. B 1983, 125, 445–448. [Google Scholar] [CrossRef]
  30. Abbott, L.F.; Wise, M.B. Constraints on Generalized Inflationary Cosmologies. Nucl. Phys. B 1984, 244, 541–548. [Google Scholar] [CrossRef]
  31. Allen, B. The Stochastic Gravity Wave Background in Inflationary Universe Models. Phys. Rev. D 1988, 37, 2078. [Google Scholar] [CrossRef] [PubMed]
  32. Lucchin, F.; Matarrese, S.; Mollerach, S. The Gravitational wave contribution to CMB anisotropies and the amplitude of mass fluctuations from COBE results. Astrophys. J. Lett. 1992, 401, L49. [Google Scholar] [CrossRef]
  33. Turner, M.S.; White, M.J.; Lidsey, J.E. Tensor perturbations in inflationary models as a probe of cosmology. Phys. Rev. D 1993, 48, 4613–4622. [Google Scholar] [CrossRef]
  34. Crittenden, R.; Davis, R.L.; Steinhardt, P.J. Polarization of the microwave background due to primordial gravitational waves. Astrophys. J. Lett. 1993, 417, L13–L16. [Google Scholar] [CrossRef]
  35. Knox, L.; Turner, M.S. Detectability of tensor perturbations through CBR anisotropy. Phys. Rev. Lett. 1994, 73, 3347–3350. [Google Scholar] [CrossRef] [PubMed]
  36. Turner, M.S. Detectability of inflation produced gravitational waves. Phys. Rev. D 1997, 55, R435–R439. [Google Scholar] [CrossRef]
  37. Mukhanov, V.F.; Feldman, H.A.; Brandenberger, R.H. Theory of cosmological perturbations. Part 1. Classical perturbations. Part 2. Quantum theory of perturbations. Part 3. Extensions. Phys. Rept. 1992, 215, 203–333. [Google Scholar] [CrossRef]
  38. Riotto, A. Inflation and the theory of cosmological perturbations. ICTP Lect. Notes Ser. 2003, 14, 317–413. [Google Scholar]
  39. Guzzetti, M.C.; Bartolo, N.; Liguori, M.; Matarrese, S. Gravitational waves from inflation. Riv. Nuovo Cim. 2016, 39, 399–495. [Google Scholar] [CrossRef]
  40. Steinhardt, P.J.; Turner, M.S. A Prescription for Successful New Inflation. Phys. Rev. D 1984, 29, 2162–2171. [Google Scholar] [CrossRef]
  41. Liddle, A.R.; Lyth, D.H. COBE, gravitational waves, inflation and extended inflation. Phys. Lett. B 1992, 291, 391–398. [Google Scholar] [CrossRef]
  42. Liddle, A.R.; Parsons, P.; Barrow, J.D. Formalizing the slow roll approximation in inflation. Phys. Rev. D 1994, 50, 7222–7232. [Google Scholar] [CrossRef]
  43. Martin, J.; Ringeval, C.; Vennin, V. Encyclopædia Inflationaris. Phys. Dark Univ. 2014, 5–6, 75–235. [Google Scholar] [CrossRef]
  44. Coughlan, G.D.; Ross, G.G. Initial Conditions for Inflation. Phys. Lett. B 1985, 157, 151–156. [Google Scholar] [CrossRef]
  45. Linde, A.D. Initial conditions for inflation. Phys. Lett. B 1985, 162, 281–286. [Google Scholar] [CrossRef]
  46. Albrecht, A.; Brandenberger, R.H. On the Realization of New Inflation. Phys. Rev. D 1985, 31, 1225. [Google Scholar] [CrossRef] [PubMed]
  47. Albrecht, A.; Brandenberger, R.H.; Matzner, R. Numerical Analysis of Inflation. Phys. Rev. D 1985, 32, 1280. [Google Scholar] [CrossRef]
  48. Albrecht, A.; Brandenberger, R.H.; Matzner, R. Inflation with Generalized Initial Conditions. Phys. Rev. D 1987, 35, 429. [Google Scholar] [CrossRef] [PubMed]
  49. Kurki-Suonio, H.; Matzner, R.A.; Centrella, J.; Wilson, J.R. Inflation From Inhomogeneous Initial Data in a One-dimensional Back Reacting Cosmology. Phys. Rev. D 1987, 35, 435–448. [Google Scholar] [CrossRef]
  50. Goldwirth, D.S.; Piran, T. Inhomogeneity and the Onset of Inflation. Phys. Rev. Lett. 1990, 64, 2852–2855. [Google Scholar] [CrossRef]
  51. Goldwirth, D.S.; Piran, T. Initial conditions for inflation. Phys. Rept. 1992, 214, 223–291. [Google Scholar] [CrossRef]
  52. Kurki-Suonio, H.; Laguna, P.; Matzner, R.A. Inhomogeneous inflation: Numerical evolution. Phys. Rev. D 1993, 48, 3611–3624. [Google Scholar] [CrossRef]
  53. Iguchi, O.; Ishihara, H. Onset of inflation in inhomogeneous cosmology. Phys. Rev. D 1997, 56, 3216–3224. [Google Scholar] [CrossRef]
  54. Vachaspati, T.; Trodden, M. Causality and cosmic inflation. Phys. Rev. D 1999, 61, 023502. [Google Scholar] [CrossRef]
  55. Easther, R.; Price, L.C.; Rasero, J. Inflating an Inhomogeneous Universe. J. Cosmol. Astropart. Phys. 2014, 8, 41. [Google Scholar] [CrossRef]
  56. East, W.E.; Kleban, M.; Linde, A.; Senatore, L. Beginning inflation in an inhomogeneous universe. J. Cosmol. Astropart. Phys. 2016, 9, 10. [Google Scholar] [CrossRef]
  57. Lyth, D.H. What would we learn by detecting a gravitational wave signal in the cosmic microwave background anisotropy? Phys. Rev. Lett. 1997, 78, 1861. [Google Scholar] [CrossRef]
  58. Efstathiou, G.; Mack, K.J. The Lyth bound revisited. J. Cosmol. Astropart. Phys. 2005, 2005, 008. [Google Scholar] [CrossRef]
  59. Easther, R.; Kinney, W.H.; Powell, B.A. The Lyth bound and the end of inflation. J. Cosmol. Astropart. Phys. 2006, 0608, 004. [Google Scholar] [CrossRef]
  60. Garcia-Bellido, J.; Roest, D.; Scalisi, M.; Zavala, I. Can CMB data constrain the inflationary field range? J. Cosmol. Astropart. Phys. 2014, 1409, 006. [Google Scholar] [CrossRef]
  61. Garcia-Bellido, J.; Roest, D.; Scalisi, M.; Zavala, I. Lyth bound of inflation with a tilt. Phys. Rev. D 2014, 90, 123539. [Google Scholar] [CrossRef]
  62. Marco, A.D. Lyth Bound, eternal inflation and future cosmological missions. Phys. Rev. D 2017, 96, 023511. [Google Scholar] [CrossRef]
  63. Albrecht, A.; Steinhardt, P.J.; Turner, M.S.; Wilczek, F. Reheating an Inflationary Universe. Phys. Rev. Lett. 1982, 48, 1437. [Google Scholar] [CrossRef]
  64. Dolgov, A.D.; Linde, A.D. Baryon Asymmetry in Inflationary Universe. Phys. Lett. 1982, 116B, 329. [Google Scholar] [CrossRef]
  65. Abbott, L.F.; Farhi, E.; Wise, M.B. Particle Production in the New Inflationary Cosmology. Phys. Lett. 1982, 117B, 29. [Google Scholar] [CrossRef]
  66. Turner, M.S. Coherent Scalar Field Oscillations in an Expanding Universe. Phys. Rev. D 1983, 28, 1243. [Google Scholar] [CrossRef]
  67. Shtanov, Y. Scalar-field dynamics and reheating of the universe in chaotic inflation scenario. Ukr. Fiz. Zh. 1993, 38, 1425–1434. [Google Scholar]
  68. Marco, A.D.; Gasperis, G.D.; Pradisi, G.; Cabella, P. Energy Density, Temperature and Entropy Dynamics in Perturbative Reheating. Phys. Rev. D 2019, 100, 123532. [Google Scholar] [CrossRef]
  69. Marco, A.D.; Pradisi, G. Variable inflaton equation-of-state and reheating. Int. J. Mod. Phys. A 2021, 36, 2150095. [Google Scholar] [CrossRef]
  70. Dolgov, A.D.; Kirilova, D.P. On Particle Creation by a Time Dependent Scalar Field. Sov. J. Nucl. Phys. 1990, 51, 172. [Google Scholar]
  71. Traschen, J.H.; Brandenberger, R.H. Particle Production During Out-of-equilibrium Phase Transitions. Phys. Rev. D 1990, 42, 2491. [Google Scholar] [CrossRef]
  72. Kofman, L.; Linde, A.D.; Starobinsky, A.A. Reheating after inflation. Phys. Rev. Lett. 1994, 73, 3195. [Google Scholar] [CrossRef]
  73. Shtanov, Y.; Traschen, J.H.; Brandenberger, R.H. Universe reheating after inflation. Phys. Rev. D 1995, 51, 5438. [Google Scholar] [CrossRef]
  74. Kofman, L.; Linde, A.D.; Starobinsky, A.A. Towards the theory of reheating after inflation. Phys. Rev. D 1997, 56, 3258. [Google Scholar] [CrossRef]
  75. Greene, P.B.; Kofman, L.; Linde, A.D.; Starobinsky, A.A. Structure of resonance in preheating after inflation. Phys. Rev. D 1997, 56, 6175–6192. [Google Scholar] [CrossRef]
  76. Greene, B.R.; Prokopec, T.; Roos, T.G. Inflaton decay and heavy particle production with negative coupling. Phys. Rev. D 1997, 56, 6484. [Google Scholar] [CrossRef]
  77. Greene, P.B.; Kofman, L. Preheating of fermions. Phys. Lett. B 1999, 448, 6–12. [Google Scholar] [CrossRef]
  78. Bassett, B.A.; Tsujikawa, S.; Wands, D. Inflation dynamics and reheating. Rev. Mod. Phys. 2006, 78, 537. [Google Scholar] [CrossRef]
  79. Frolov, A.V. Non-linear Dynamics and Primordial Curvature Perturbations from Preheating. Class. Quant. Grav. 2010, 27, 124006. [Google Scholar] [CrossRef]
  80. Allahverdi, R.; Brandenberger, R.; Cyr-Racine, F.Y.; Mazumdar, A. Reheating in Inflationary Cosmology: Theory and Applications. Ann. Rev. Nucl. Part. Sci. 2010, 60, 27–51. [Google Scholar] [CrossRef]
  81. Amin, M.A.; Hertzberg, M.P.; Kaiser, D.I.; Karouby, J. Nonperturbative Dynamics of Reheating After Inflation: A Review. Int. J. Mod. Phys. D 2014, 24, 1530003. [Google Scholar] [CrossRef]
  82. Lozanov, K.D. Lectures on Reheating after Inflation. arXiv 2019, arXiv:1907.04402. [Google Scholar]
  83. Belinsky, V.A.; Khalatnikov, I.M.; Grishchuk, L.P.; Zeldovich, Y.B. Inflationary stages in cosmological models with a scalar field. Phys. Lett. B 1985, 155, 232–236. [Google Scholar] [CrossRef]
  84. Guo, Z.K.; Piao, Y.S.; Cai, R.G.; Zhang, Y.Z. Inflationary attractor from tachyonic matter. Phys. Rev. D 2003, 68, 043508. [Google Scholar] [CrossRef]
  85. Urena-Lopez, L.A.; Reyes-Ibarra, M.J. On the dynamics of a quadratic scalar field potential. Int. J. Mod. Phys. D 2009, 18, 621–634. [Google Scholar] [CrossRef]
  86. Remmen, G.N.; Carroll, S.M. Attractor Solutions in Scalar-Field Cosmology. Phys. Rev. D 2013, 88, 083518. [Google Scholar] [CrossRef]
  87. Ellis, G.F.R.; Madsen, M.S. Exact scalar field cosmologies. Class. Quant. Grav. 1991, 8, 667–676. [Google Scholar] [CrossRef]
  88. Lidsey, J.E. The Scalar field as dynamical variable in inflation. Phys. Lett. B 1991, 273, 42–46. [Google Scholar] [CrossRef]
  89. Salopek, D.S.; Bond, J.R. Nonlinear evolution of long wavelength metric fluctuations in inflationary models. Phys. Rev. D 1990, 42, 3936–3962. [Google Scholar] [CrossRef]
  90. Muslimov, A.G. On the Scalar Field Dynamics in a Spatially Flat Friedman Universe. Class. Quant. Grav. 1990, 7, 231–237. [Google Scholar] [CrossRef]
  91. Barrow, J.D. Graduated inflationary universes. Phys. Lett. B 1990, 235, 40–43. [Google Scholar] [CrossRef]
  92. Barrow, J.D. New types of inflationary universe. Phys. Rev. D 1993, 48, 1585–1590. [Google Scholar] [CrossRef]
  93. Barrow, J.D. Exact inflationary universes with potential minima. Phys. Rev. D 1994, 49, 3055–3058. [Google Scholar] [CrossRef]
  94. Barrow, J.D.; Parsons, P. Inflationary models with logarithmic potentials. Phys. Rev. D 1995, 52, 5576–5587. [Google Scholar] [CrossRef] [PubMed]
  95. Lidsey, J.E.; Liddle, A.R.; Kolb, E.W.; Copeland, E.J.; Barreiro, T.; Abney, M. Reconstructing the inflation potential: An overview. Rev. Mod. Phys. 1997, 69, 373–410. [Google Scholar] [CrossRef]
  96. Smoot, G.F.; Bennett, C.L.; Kogut, A.; Wright, E.L.; Aymon, J.; Boggess, N.W.; Cheng, E.S.; De Amici, G.; Gulkis, S.; Hauser, M.G.; et al. Structure in the COBE differential microwave radiometer first year maps. Astrophys. J. Lett. 1992, 396, L1–L5. [Google Scholar] [CrossRef]
  97. Akrami, Y.; Arroja, F.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. Planck 2018 results. X. Constraints on inflation. Astron. Astrophys. 2020, 641, A10. [Google Scholar] [CrossRef]
  98. Ade, P.A.; Ahmed, Z.; Amiri, M.; Barkats, D.; Thakur, R.B.; Bischoff, C.A.; Beck, D.; Bock, J.J.; Boenish, H.; Bullock, E.; et al. Improved Constraints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season. Phys. Rev. Lett. 2021, 127, 151301. [Google Scholar] [CrossRef] [PubMed]
  99. Liddle, A.R.; Leach, S.M. How long before the end of inflation were observable perturbations produced? Phys. Rev. D 2003, 68, 103503. [Google Scholar] [CrossRef]
  100. Germán, G.; Quaglia, R.G.; Colorado, A.M.M. Model independent bounds for the number of e-folds during the evolution of the universe. JCAP 2023, 03, 004. [Google Scholar] [CrossRef]
  101. Das, K.; Dutta, K.; Maharana, A. Inflationary Predictions and Moduli Masses. Phys. Lett. B 2015, 751, 195–200. [Google Scholar] [CrossRef]
  102. Cicoli, M.; Dutta, K.; Maharana, A.; Quevedo, F. Moduli Vacuum Misalignment and Precise Predictions in String Inflation. J. Cosmol. Astropart. Phys. 2016, 08, 006. [Google Scholar] [CrossRef]
  103. Maharana, A.; Zavala, I. Postinflationary scalar tensor cosmology and inflationary parameters. Phys. Rev. D 2018, 97, 123518. [Google Scholar] [CrossRef]
  104. Marco, A.D.; Pradisi, G.; Cabella, P. Inflationary scale, reheating scale, and pre-BBN cosmology with scalar fields. Phys. Rev. D 2018, 98, 123511. [Google Scholar] [CrossRef]
  105. Kawasaki, M.; Yamaguchi, M.; Yanagida, T. Natural chaotic inflation in supergravity. Phys. Rev. Lett. 2000, 85, 3572–3575. [Google Scholar] [CrossRef] [PubMed]
  106. Takahashi, F. Linear Inflation from Running Kinetic Term in Supergravity. Phys. Lett. B 2010, 693, 140–143. [Google Scholar] [CrossRef]
  107. Nakayama, K.; Takahashi, F. Running Kinetic Inflation. J. Cosmol. Astropart. Phys. 2010, 11, 009. [Google Scholar] [CrossRef]
  108. Nakayama, K.; Takahashi, F. Higgs Chaotic Inflation in Standard Model and NMSSM. J. Cosmol. Astropart. Phys. 2011, 02, 010. [Google Scholar] [CrossRef]
  109. Nakayama, K.; Takahashi, F. General Analysis of Inflation in the Jordan frame Supergravity. J. Cosmol. Astropart. Phys. 2010, 11, 039. [Google Scholar] [CrossRef]
  110. Harigaya, K.; Ibe, M.; Schmitz, K.; Yanagida, T.T. Chaotic Inflation with a Fractional Power-Law Potential in Strongly Coupled Gauge Theories. Phys. Lett. B 2013, 720, 125–129. [Google Scholar] [CrossRef]
  111. Silverstein, E.; Westphal, A. Monodromy in the CMB: Gravity Waves and String Inflation. Phys. Rev. D 2008, 78, 106003. [Google Scholar] [CrossRef]
  112. McAllister, L.; Silverstein, E.; Westphal, A. Gravity Waves and Linear Inflation from Axion Monodromy. Phys. Rev. D 2010, 82, 046003. [Google Scholar] [CrossRef]
  113. McAllister, L.; Silverstein, E.; Westphal, A.; Wrase, T. The Powers of Monodromy. J. High Energy Phys. 2014, 1409, 123. [Google Scholar] [CrossRef]
  114. Lucchin, F.; Matarrese, S. Power Law Inflation. Phys. Rev. D 1985, 32, 1316. [Google Scholar] [CrossRef] [PubMed]
  115. Halliwell, J.J. Scalar Fields in Cosmology with an Exponential Potential. Phys. Lett. B 1987, 185, 341. [Google Scholar] [CrossRef]
  116. Burd, A.B.; Barrow, J.D. Inflationary Models with Exponential Potentials. Nucl. Phys. B 1988, 308, 929–945. [Google Scholar] [CrossRef]
  117. Kallosh, R.; Linde, A. Universality Class in Conformal Inflation. J. Cosmol. Astropart. Phys. 2013, 1307, 002. [Google Scholar] [CrossRef]
  118. Kallosh, R.; Linde, A. Multi-field Conformal Cosmological Attractors. J. Cosmol. Astropart. Phys. 2013, 1312, 006. [Google Scholar] [CrossRef]
  119. Ferrara, S.; Kallosh, R.; Linde, A.; Porrati, M. Minimal Supergravity Models of Inflation. Phys. Rev. D 2013, 88, 085038. [Google Scholar] [CrossRef]
  120. Kallosh, R.; Linde, A. Superconformal generalizations of the Starobinsky model. J. Cosmol. Astropart. Phys. 2013, 1306, 028. [Google Scholar] [CrossRef]
  121. Kallosh, R.; Linde, A.; Roest, D. Superconformal Inflationary α-Attractors. J. High Energy Phys. 2013, 1311, 198. [Google Scholar] [CrossRef]
  122. Galante, M.; Kallosh, R.; Linde, A.; Roest, D. Unity of Cosmological Inflation Attractors. Phys. Rev. Lett. 2015, 114, 141302. [Google Scholar] [CrossRef]
  123. Kallosh, R.; Linde, A. Escher in the Sky. Comptes Rendus Phys. 2015, 16, 914–927. [Google Scholar] [CrossRef]
  124. Kallosh, R.; Linde, A. Cosmological Attractors and Asymptotic Freedom of the Inflaton Field. J. Cosmol. Astropart. Phys. 2016, 1606, 047. [Google Scholar] [CrossRef]
  125. Maeda, K.i. Inflation as a Transient Attractor in R2 Cosmology. Phys. Rev. D 1988, 37, 858. [Google Scholar] [CrossRef] [PubMed]
  126. Goncharov, A.B.; Linde, A.D. Chaotic Inflation in Supergravity. Phys. Lett. B 1984, 139, 27–30. [Google Scholar] [CrossRef]
  127. Goncharov, A.S.; Linde, A.D. Chaotic inflation of the universe in supergravity. Sov. Phys. JETP 1984, 59, 930–933. [Google Scholar]
  128. Pradisi, G.; Salvio, A. (In)equivalence of metric-affine and metric effective field theories. Eur. Phys. J. C 2022, 82, 840. [Google Scholar] [CrossRef]
  129. Salvio, A. Inflating and reheating the Universe with an independent affine connection. Phys. Rev. D 2022, 106, 103510. [Google Scholar] [CrossRef]
  130. Marco, A.D.; Orazi, E.; Pradisi, G. Einstein– Cartan pseudoscalaron inflation. Eur. Phys. J. C 2024, 84, 146. [Google Scholar] [CrossRef]
  131. Cicoli, M.; Conlon, J.P.; Quevedo, F. General Analysis of LARGE Volume Scenarios with String Loop Moduli Stabilisation. J. High Energy Phys. 2008, 10, 105. [Google Scholar] [CrossRef]
  132. Cicoli, M.; Burgess, C.P.; Quevedo, F. Fibre Inflation: Observable Gravity Waves from IIB String Compactifications. J. Cosmol. Astropart. Phys. 2009, 03, 013. [Google Scholar] [CrossRef]
  133. Cicoli, M.; Muia, F.; Shukla, P. Global Embedding of Fibre Inflation Models. J. High Energy Phys. 2016, 11, 182. [Google Scholar] [CrossRef]
  134. Cicoli, M.; Valentino, E.D. Fitting string inflation to real cosmological data: The fiber inflation case. Phys. Rev. D 2020, 102, 043521. [Google Scholar] [CrossRef]
  135. Acquaviva, V.; Bartolo, N.; Matarrese, S.; Riotto, A. Second order cosmological perturbations from inflation. Nucl. Phys. B 2003, 667, 119–148. [Google Scholar] [CrossRef]
  136. Maldacena, J.M. Non-Gaussian features of primordial fluctuations in single field inflationary models. JHEP 2003, 05, 013. [Google Scholar] [CrossRef]
  137. Bartolo, N.; Komatsu, E.; Matarrese, S.; Riotto, A. Non-Gaussianity from inflation: Theory and observations. Phys. Rept. 2004, 402, 103–266. [Google Scholar] [CrossRef]
  138. Babich, D.; Creminelli, P.; Zaldarriaga, M. The Shape of non-Gaussianities. J. Cosmol. Astropart. Phys. 2004, 08, 009. [Google Scholar] [CrossRef]
  139. Chen, X. Primordial Non-Gaussianities from Inflation Models. Adv. Astron. 2010, 2010, 638979. [Google Scholar] [CrossRef]
  140. Arkani-Hamed, N.; Maldacena, J. Cosmological Collider Physics. arXiv 2015, arXiv:1503.08043. [Google Scholar]
  141. Sasaki, M.; Stewart, E.D. A General analytic formula for the spectral index of the density perturbations produced during inflation. Prog. Theor. Phys. 1996, 95, 71–78. [Google Scholar] [CrossRef]
  142. Lyth, D.H.; Malik, K.A.; Sasaki, M. A General proof of the conservation of the curvature perturbation. J. Cosmol. Astropart. Phys. 2005, 05, 004. [Google Scholar] [CrossRef]
  143. Wands, D.; Malik, K.A.; Lyth, D.H.; Liddle, A.R. A New approach to the evolution of cosmological perturbations on large scales. Phys. Rev. D 2000, 62, 043527. [Google Scholar] [CrossRef]
  144. Dias, M.; Ribeiro, R.H.; Seery, D. The δN formula is the dynamical renormalization group. J. Cosmol. Astropart. Phys. 2013, 10, 062. [Google Scholar] [CrossRef]
  145. Sugiyama, N.S.; Komatsu, E.; Futamase, T. δN formalism. Phys. Rev. D 2013, 87, 023530. [Google Scholar] [CrossRef]
  146. Abolhasani, A.A.; Firouzjahi, H.; Naruko, A.; Sasaki, M. Delta N Formalism in Cosmological Perturbation Theory; WSP: Phoenix, AZ, USA, 2019. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Di Marco, A.; Orazi, E.; Pradisi, G. Introduction to the Number of e-Folds in Slow-Roll Inflation. Universe 2024, 10, 284. https://doi.org/10.3390/universe10070284

AMA Style

Di Marco A, Orazi E, Pradisi G. Introduction to the Number of e-Folds in Slow-Roll Inflation. Universe. 2024; 10(7):284. https://doi.org/10.3390/universe10070284

Chicago/Turabian Style

Di Marco, Alessandro, Emanuele Orazi, and Gianfranco Pradisi. 2024. "Introduction to the Number of e-Folds in Slow-Roll Inflation" Universe 10, no. 7: 284. https://doi.org/10.3390/universe10070284

APA Style

Di Marco, A., Orazi, E., & Pradisi, G. (2024). Introduction to the Number of e-Folds in Slow-Roll Inflation. Universe, 10(7), 284. https://doi.org/10.3390/universe10070284

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop