Clocking the End of Cosmic Inflation
Abstract
Making observable predictions for cosmic inflation requires determining when the wavenumbers of astrophysical interest today exited the Hubble radius during the inflationary epoch. These instants are commonly evaluated using the slow-roll approximation and measured in e-folds , in reference to the e-fold at which inflation ended. Slow roll being necessarily violated towards the end of inflation, both the approximated trajectory and are determined at, typically, one or two e-folds precision. Up to now, such an uncertainty has been innocuous, but this will no longer be the case with the forthcoming cosmological measurements. In this work, we introduce a new and simple analytical method, on top of the usual slow-roll approximation, that reduces uncertainties on to less than a tenth of an e-fold.
1 Introduction
Cosmic Inflation is a phase of accelerated expansion of the primordial universe which addresses various puzzles of the Big-Bang model as, for instance, the so-called horizon problem and the smallness of the spatial curvature today [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. Inflation also provides a convincing and simple physical explanation for the origin of cosmic structures: they are seeded by vacuum quantum fluctuations of both the metric and a yet unknown scalar degree of freedom [13].
In its simplest incarnation, both the accelerated expansion and the quantum fluctuations are the outcome of a self-gravitating scalar field , named the inflaton, slowly rolling down its potential energy . This class of scenarios is a populated landscape counting hundreds of models, all of them making definite predictions which can be confronted by cosmological observations [14, 15, 16]. As of today, of the proposed scenarios in this class have been ruled-out by Cosmic Microwave Background (CMB) and Large Scale Structure measurements [17, 18, 19]. Still, of the remaining models are compatible with the data. With the deployment of ground-based CMB-S4 polarization telescopes [20, 21, 22], the soon-to-be released Euclid satellite data [23, 24], unprecedented galaxy surveys [25], and the search for -mode polarization from space by the LiteBIRD satellite [26], one should reasonably expect many of the remaining models to be disambiguated and tested. This necessitates, however, that theoretical predictions are made at the required accuracy [27, 28].
In a brute-force manner, one can simply solve the field and gravitational evolution numerically. However, the underlying gravity theory is General Relativity, and, as of today, only parts of the inflationary evolution can be solved without approximation [29, 30, 31, 32, 33, 34, 35, 36]. Considering the regime in which the quantum fluctuations do not dominate the dynamics [37, 38, 39, 40, 41, 42], one can alternatively solve for linear and non-linear perturbations numerically, around a homogeneous background, without any other approximations [43, 44, 45, 46, 47, 48, 49, 50, 51]. These methods are accurate as long as the gravitational effects remain small, but they are computationally too much demanding when dealing with hundreds of different models [52, 53].
There exists, however, a model-free and perturbative treatment for the single-field scenarios in which the assumption of “slow roll” can be made. This approach, initiated in Ref. [1] for the tensor modes, has been extended to scalar perturbations in Refs. [54, 55] and generalized to higher orders in [56, 57, 58, 59, 60, 61, 62, 63]. It has also found applications out of the original context and can be extended to other classes of inflationary models [64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85]. In modern terminology, the slow-roll approximation introduces the Hubble-flow functions defined by [59]
(1.1) |
where denotes the Hubble parameter during inflation and is the number of e-folds, being the Friedmann-Lemaître-Robertson-Walker (FLRW) scale factor. For a quasi-de Sitter accelerated expansion, is nearly constant and all the Hubble-flow functions are expected to be small. It is therefore possible to solve the linearized Einstein’s equations for both the tensor and scalar perturbations by performing a consistent expansion in terms of these functions. Analytical solutions have currently been derived up to third order [63] and they allow us to calculate the primordial power spectra of the comoving curvature perturbation and of the primordial gravitational waves . For instance, keeping only the first order terms, one gets, for the power spectrum of the curvature perturbations,
(1.2) |
where the constant and is a wavenumber around which the expansion is made (an observer choice). For the wavenumbers probed by the Cosmic Microwave Background anisotropies, one usually takes to be in the middle of the observable range of modes. All the other “starred” quantities in Eq. 1.2 refer to the Hubble-flow functions evaluated at a given e-fold number , i.e., and . This e-fold number is the time at which the physical pivot wavenumber exited the Hubble radius during inflation, namely the solution of111For practical reasons, is usually defined in terms of the conformal time by , which coincides with Hubble radius crossing at leading order in the Hubble-flow functions. For the present discussion, these differences will not play a role, but they are important when considering higher order terms [63].
(1.3) |
As a result, even if the accuracy at which Eq. 1.2 is derived is under control, another source of uncertainties in making observable predictions comes from our ability to determine a precise value for , and, as we will see, for . This is the main focus of this paper.
Before entering into details, let us further express Eq. 1.3 in terms of observable quantities. The physical pivot wavenumber is measured today, for a scale factor given by , in terms of which Eq. 1.3 reads
(1.4) |
We have made explicit , the redshift at which inflation ended. It depends on the universe history after inflation and, in particular, it is sensitive to the so-called reheating era. Following Refs. [52, 86], one can conveniently absorb all the kinematic effects associated with this era into the reheating parameter defined by
(1.5) |
Here and stand for the energy density of the universe at the end of inflation and at beginning of the radiation era (the end reheating), respectively. In terms of , one has
(1.6) |
where is a measure of the change of number of entropic () and energetic () relativistic degrees of freedom between the beginning of the radiation era and today [87]. For instance, one has for the Standard Model [88]. One can further expand Eqs. 1.4 to 1.6 for single-field inflationary models by making use of the Friedmann-Lemaître equations for a self gravitating scalar field . As shown in Section 2, they allow us to express the Hubble parameter during inflation as
(1.7) |
At the end of inflation, one therefore has
(1.8) |
Plugging this expression into Eq. 1.6, taking the logarithm of Eq. 1.4, one finally gets
(1.9) |
where is defined by
(1.10) |
the absolute value of which giving the typical number of e-folds of decelerated expansion after inflation. In all practical situations, one has and, for consistency, we have kept only the leading order terms in while deriving Eq. 1.9. Moreover, in the last term, we have made explicit the quantity , which is a very well measured observable as [17].
Let us now explain how to determine the value of under the hypothesis that an inflationary model, given by its potential , is specified. Any possible reheating history is associated with definite values for (and ). For instance, a radiation-like, or an instantaneous reheating, are both associated with . In this situation, Eq. 1.9 ends up being a simple algebraic equation for provided one has an explicit expression for , and . In other words, one must determine the field trajectory to solve Eq. 1.9 at given reheating history. Usually, this cannot be made exactly and one has to resort to an exact numerical integration, or, to some slow roll approximation to evaluate the field trajectory . Notice that, it is also possible to interpret Eq. 1.9 as an algebraic equation on , but this still requires determining the field trajectory in order to evaluate .
The fastest and most practical method used to determine the field trajectory is the slow-roll approximation. As we show in Section 2.2, it induces errors, which have been, up to now, not a concern. Indeed, most of the theoretical unknowns in Eq. 1.9 are actually associated with the reheating, namely the values of ( does not have significant effects provided the number of relativistic degrees of freedom does not take exponentially large values [89]). The actual value of is unknown by orders of magnitude and, in principle, it is allowed to vary from a lower bound as small as Big-Bang Nucleosynthesis to which can be as large as . Under very reasonable assumptions, one can show that (see Ref. [90, 52]).
This justifies why questioning the accuracy at which is evaluated was not a concern. However, as shown in Ref. [19], the current cosmological data are now constraining the reheating era and models having exactly the same accelerated inflationary phase but differing only by their reheating histories, i.e., predicting different values of , can now be disambiguated. From another point of view, even for inflationary scenarios not specifying the reheating, the current data allow us to determine the favoured values of . Any uncertainty in the determination of will then bias the constraints on . As such, it is becoming relevant to improve the accuracy at which the function can be determined.
The paper is organized as follows. In Section 2, we recap how to obtain the field trajectory within a FLRW metric and detail the slow-roll method commonly used to approximate the solution. In particular, we use a numerical integration to discuss the amplitude and the origin of the uncertainties made by using the slow-roll approximated trajectory instead of the exact one. In Section 3, we present new analytical results and an exact expansion of the trajectory which allow us to propose a “velocity correction” to the traditional slow-roll. We show that such a correction reduces the uncertainties by an order of magnitude. Section 4 is dedicated to the problem of determining the field value at which inflation ends, which is another (small) source of errors on the determination of . We present various analytical approaches to address this issue and test them within various inflationary scenarios. Here as well, we show that our method reduces the uncertainties on by an order of magnitude. Our conclusion are presented in Section 5.
2 Basics on the field trajectory
2.1 Equations of motion
In the following, we assume a minimally coupled single scalar field within a FLRW metric. The Friedmann-Lemaître and Klein-Gordon equations read
(2.1) | ||||
(2.2) | ||||
(2.3) |
where a dot denotes differentiation with respect to the cosmic time and . In terms of the number of e-fold , these equations can be decoupled. From Eqs. 1.1, 2.1 and 2.2, one has
(2.4) |
and the first Hubble-flow function measures the kinetic energy of the field when time is counted in e-fold. In order to simplify the notations, let us introduce the “field velocity” in e-fold as
(2.5) |
such that . Expressing Eqs. 2.1 and 2.2 in e-fold time, one obtains Eq. 1.7 for the Hubble parameter, which can be finally plugged into Eq. 2.3 to obtain a decoupled equation of motion for
(2.6) |
From now on, we will be working in Planck units with such that, making use of Eq. 2.4, the previous equation simplifies to
(2.7) |
As discussed in Ref. [18], this equation is similar to the one of a relativistic particle in presence of friction and accelerated by an external force created by the potential , the value giving the maximal possible speed for the field . Indeed, positivity of Eq. 1.7 enforces that all field trajectories must satisfy , i.e., .
2.2 Slow-roll trajectory
There is no known analytical solution of Eq. 2.7 for an unspecified potential , although various approximated solutions, in different regimes, have been derived [18] (see, however, Appendix A). In the slow-roll regime we are interested in, one can remark that the field acceleration can also be expressed in terms of . From Eqs. 2.5 and 1.1, one has
(2.8) |
in terms of which Eq. 2.7 reads
(2.9) |
Assuming a slowly rolling field evolution implies that all the are small and, at leading order, Eq. 2.9 can be approximated by
(2.10) |
which has the solution
(2.11) |
Here the prime stands for the derivative with respect to field value. Another way to interpret this slow-roll trajectory is to remark that the acceleration term of Eq. 2.7 is ignored, which means that we are only considering the friction dominated regime. In fact, were the force term on the right-hand-side be constant, Eq. 2.10 would give the exact terminal velocity of Eq. 2.7, and, Eq. 2.11 would be the exact attractor solution. In the general case, however, there is a small drift sourced by the non-constancy of the force term and this induces differences between and the exact attractor solution of Eq. 2.7.
2.3 Characterizing the end of inflation
By definition, inflation stands for an accelerated expansion of the spacetime, i.e., . From Eq. 2.4, using , one has and the condition for acceleration translates into . In the vanilla single-field inflationary models, the accelerated expansion ends by itself with a so-called “graceful exit”: the potential becomes steeper, and the field accelerates up to the point at which
(2.12) |
Translated into velocities, one therefore has , the sign being related to the direction in which inflation proceeds. Indeed, depending on the shape of the potential, either the field increases during inflation and , or it decreases and . It is also possible that inflation ends by another mechanism than a graceful exit, as for instance by a tachyonic instability triggered by an extra field, as in the prototypical hybrid inflation model [91]. In that situation, is no longer set by the condition . Instead, it may be viewed as an additional model parameter. The determination of in these situations has, therefore, nothing to do with the inflationary dynamics, and we will not consider these cases. Let us however stress that a tachyonic instability is relevant only if it affects the inflaton while , otherwise it would rather be interpreted as an event belonging to the reheating era.
Solving Eq. 2.12 for is problematic in various aspects. It is a condition on the first Hubble-flow function, or equivalently, on the field velocity , whose e-fold dependency would require to solve Eq. 2.7 exactly. Without knowing the exact solution, the best one can do is to use the approximation of Eq. 2.10 and solve instead of . However, by doing so, we break our working hypothesis that the functions have to be small, as the end of inflation is indeed manifestly violating slow roll.
In spite of this, in essentially all works on slow-roll inflation, is the value actually used for . As we will demonstrate in Section 3, this is quite a good approximation because turns out to be the leading order solution of yet another expansion of the field trajectory valid even when slow roll is violated. Another more intuitive explanation justifying the extrapolation of to the end of inflation is to remark that when slow-roll is violated inflation cannot be sustained for a long time, typically not more than e-fold. As such, one cannot make a larger error than that on by using instead.
In practice, solving consists in finding the root of the algebraic equation
(2.13) |
and injecting it into the slow-roll trajectory of Eq. 2.11, i.e., .
2.4 Assessing slow-roll accuracy
In this section, we compare, for various potentials, the slow-roll approximated trajectory defined by Eq. 2.14, to an exact numerical integration of Eq. 2.7 complemented by a root finding algorithm to numerically determine , the solution of .
As discussed in the previous sections, the errors made by using instead of come from both the approximation and . In order to separate both, let us define a semi-numerical solution, built upon the slow-roll trajectory, based on the exact field value for the end of inflation
(2.15) |
All these functions take as input the field value and return some approximated number of e-folds. Once we have (numerically) integrated the field trajectory exactly, we have at our disposal the functions , as well as the value of and . From these, we can numerically determine the exact functions and .
Starting from some initial conditions, at and , a first approach would be to compare the exact solution to its slow-roll approximated version . Equivalently, one could also compare to . Let us first remark that the field value plays no role as, for a given potential, and once on the attractor, the trajectory is universal and always ends in the same manner. Intuitively, one expects , as well as , to be small deep in slow-roll while growing towards the end of inflation and this is exactly what happens. However, we have chosen not to show these trajectories in the following. Indeed, as discussed at length in the introduction, the observable quantity entering the reheating equation is , in reference to the end of inflation. As such, any errors damaging the actual value of , as the ones building up close to the end of inflation, will be necessarily folded into all the values of , even if and match well in those regions. Hence, it is actually much more informative to compare , to the exact , all of these quantities being sensitive to the end of inflation.
Last but not least, the values taken by are also not very much informative as observable predictions are mostly sensitive to e-fold numbers. Knowing the exact trajectory , we can easily trade for and discuss all error made in terms of the latter quantity. This is relevant because the prototypical value of and minimizing errors is particularly important around these figures rather than towards the end of inflation, or, much earlier.
In Fig. 1, we have therefore represented, for various models, as a function of , the absolute errors in the number of e-folds made by using the slow-roll approximated trajectories instead of the exact one. The red curve in these plots shows whereas the blue curve is for the semi-analytical trajectory . The differences between the blue and red curves are thus coming from the uncertainties in determining the field value at which inflation ends.
The six models considered in Fig. 1 are a few representative of the ones discussed in the Encyclopædia Inflationaris paper of Ref. [14]. We have picked up a quadratic large field inflation model (), having a potential , a quartic small field inflation model () with where the vacuum expectation value is super-Planckian, Starobinsky Inflation () having , a quadratic T-model inflation () with , an exponential supersymmetric inflation model () with and a pseudo-natural inflation model () having a potential with . Not all of these models are compatible with current cosmological data, for instance is strongly disfavoured whereas has its parameters purposely chosen to violate slow-roll ( during inflation). The small field scenario is a model which is compatible with the data whereas , and are different incarnations of the so-called plateau-type models and belong to most favoured scenarios [19].
Let us remark in Fig. 1 that, for all models, the errors generated by the slow-roll trajectory of Eq. 2.11 grow with the number of e-folds before the end of inflation. This growth is precisely due to the small terms omitted by making the assumption and confirms that is slightly off-track the true slow rolling attractor velocity. One can also notice that the error jumps quite fast close to whereas, up to one model (), it increases like a logarithm at larger values of . This is due to the fact that slow-roll is most violated towards the end of inflation and the assumption is quite wrong in these regions. For (lower right panel), the errors seem to increase linearly with , as opposed to a logarithm-like growth. The reason being that, as aforementioned, it is far from slow roll also during inflation ().
Finally, these plots confirm that for the fiducial value , the typical errors on the trajectory are e-folds. Only for the extreme slow-roll violating model , one gets a larger, but still reasonable error.
3 Correcting slow-roll
In this section, we address the main source of error eroding the traditional slow-roll trajectory: being slightly off-track the attractor solution.
3.1 Integral constraints
Let us first show that, even though the slow-roll trajectory is not right on the attractor, the shift with respect to the exact solution is actually bounded. One can define the absolute error
(3.1) |
which can be viewed as a function by formally making use of the exact field trajectory for . From the definition of in Eq. 2.5 (still in Planck units) one can rewrite Eq. 2.7 in terms of as
(3.2) |
This equation can actually be integrated by separating variables and isolating as
(3.3) |
The left-hand side can be integrated exactly and, using , one gets the integral constraint
(3.4) |
In the slow roll regime and the integrated error made between and is . Notice the negative sign, which shows that, for , one has and the approximated trajectory is slightly advanced compared to the exact one (see also Fig. 1).
We can also derive a second integral constraint by integrating with respect to the number of e-fold. Starting again from Eq. 2.7 and separating the variables and , one has
(3.5) |
Again, the left-hand side can be integrated exactly while the right hand side can be expressed in terms of by using Eq. 2.5. One obtains another integral constraint
(3.6) |
where is the relative error between the slow-roll and exact trajectory. The sign is for , depending on which direction inflation proceeds. Provided , the relative integrated error is bounded and reads . We recover that, for , and .
Both Eqs. 3.4 and 3.6 are finite and shows that both the absolute error and the relative error are under control, deep in the slow roll regime as well as at the end of inflation when approaches . This suggests using either or as a small parameter against which the exact solution of Eq. 2.7 can be expanded.
3.2 New expansion for the field trajectory
From the exact field velocity , the true number of e-fold is given, up to a constant, by
(3.7) |
Instead, the traditional slow-roll approximation replaces it with Eq. 2.11. Let us use Eq. 3.7 to express the exact field trajectory as
(3.8) |
where we have artificially introduced the known function defined in Eq. 2.10. From the definition of in Eq. 3.1, one has
(3.9) |
which can be plugged into Eq. 3.8 to get
(3.10) |
the first term giving back the traditional slow-roll approximation. The second term can be further expanded by remarking that
(3.11) |
giving a new and exact expansion for the field trajectory
(3.12) |
Let us notice the first term, which is as the field value here has to be the exact one. Quite importantly, this expansion is not based on the usual slow-roll expansion, one does not need to assume . Instead, the “small parameter” is the relative error functional, , which is ensured to be under control thanks to the integral constraints derived earlier. The benefit of having expanded the exact trajectory as in Eq. 3.12 is that, as we show in the next section, the second term, which acts as a first correction, is exactly calculable.
Strictly speaking, the expansion of Eq. 3.11 is converging only if the relative error . Although this is ensured for most of the inflationary trajectory, in slow roll, it may exceed unity very close to the end of inflation for . However, in order to satisfy the integral constraint of Eq. 3.6, the field domain over which this happens must be small (in Planck units). Similarly, Eq. 3.4 gives a constraint on the absolute error over time, this one cannot not be of order unity for more than a fraction of an e-fold. Although these cases are not of immediate interest when considering , it is interesting to remark that for , one has and another expansion can be performed
(3.13) |
Plugging this expression into Eq. 3.8, one gets another exact expansion
(3.14) |
which shows that, up to the field value at which inflation ends, the usual slow-roll approximated trajectory remains the leading order term.
3.3 Velocity correction
Let us now assume that we are in the field domain for which . From Eq. 3.2, one has, exactly
(3.15) |
such that
(3.16) |
where use has been made of . As such, the first correction appearing in the expansion of Eq. 3.12 is a simple velocity correction. Even in slow-roll, for , we see that this term matters. As a matter of fact, it grows logarithmically when becomes small, and it cancels most of the errors associated with .
We can also rewrite Eq. 3.12 as
(3.17) |
to render explicit the deviations with respect to . The third term appears because assumes inflation to end at instead of the exact value . We will discuss this issue in Section 4.
An additional issue with Eq. 3.17 is that, in principle, one does not know . However, we are interested in far from the end of inflation, where slow roll is verified. As such, it is perfectly justified to evaluate within the second term of Eq. 3.17. If a higher precision is needed, it is always possible to account for higher-derivative perturbative corrections by determining in terms of [92, 19]. In fact, it is interesting to compare our velocity refinement, the second term of Eq. 3.17, to these higher-derivative corrections. As shown in Ref. [93], including the next-to-leading order term to determine yields a correction to in . This can be compared to first term in the right-hand side of Eq. 3.16, i.e., . Deep in slow roll, using , it approximates to and this improves the next-to-leading order higher-derivative correction by a constant shift of e-fold. Let us also stress that, as opposed to perturbative higher-order terms, Eq. 3.16 is exact and this is why we can safely incorporate into the velocity correction some effects coming from the end of inflation (the terms involving ).
In order to check the accuracy of Eq. 3.17, we have plotted in Fig. 2 the absolute error (blue curve) where
(3.18) |
Let us stress that we have traded for in this expression. Compared to the traditional slow-roll trajectory (red), it is evident that the velocity corrections erase the logarithmic error growth with respect to by a factor . In view of such a success, one may be tempted in trying to evaluate the higher order terms of Eq. 3.17 similarly. However, because they involve powers of , using Eq. 3.15 does not allow for an exact integration, even though some parts can still be estimated. We have also tried to calculate exactly these terms starting from Eq. 2.7 and this has allowed us to derive a new exact formula, presented in Appendix A, but which would require the knowledge of , which is not usually the case (see, however, Section 4.2).
4 Pinpointing the end of inflation
As explained in Section 2.3, the uncertainties associated with come from the trading between and at the end of inflation, when slow-roll is manifestly violated. A first approach may be to use Eq. 3.9 close to the end of inflation and making use of Eq. 3.15 while replacing by in the derivative. However, we are now in a regime in which is not necessarily close to and could also exceed unity. We have checked that the relative error indeed exceeds unity at the end of inflation for two of the tested models: and . A similar approach has been discussed in Refs. [94, 95] and it has been shown to reduce the error in determining the end of inflation down to for some specific potentials ( and -like, which both satisfy ). In the following, we will present other methods, performing only slightly better ( uncertainties on ) but designed to be robust for all slow-roll models.
4.1 Constrained extrapolations
For a given potential, one has a perfect knowledge of and , and we would like to have an accurate determination of close to . As such, we could circumvent the determination of by trying to approximate instead the integral constraints of the error function towards the end of inflation.
For instance, using a trapezoidal approximation for the integrals, one can approximate the first constraint Eq. 3.4 at the end of inflation as
(4.1) |
Similarly, one may approximate the second constraint Eq. 3.6, as
(4.2) |
Under these approximations, Eqs. 4.1 and 4.2 only involve two unknown and independent variables and . Indeed, the error function in both points can be written as
(4.3) |
where, as before, . This algebraic system can be solved numerically to obtain an estimate of the field value at the end of inflation . Let us stress that, according to the previous discussion, the domain in which these equations are solved requires and, either or , depending on whether inflation proceeds at decreasing or increasing field values, respectively.
Let us notice that the method could also be accommodated with other functional shapes to model in the constraint integrals. We have tested some power-law and exponential extrapolations, but they do not perform better than the simple trapezoidal rule presented here. In the next section, we discuss a similar method anchored on a family of exact field trajectories.
4.2 Matching to Mukhanov inflation
Mukhanov inflation is one of the very few inflationary models for which the exact field trajectory is analytically known, and, the only one which exhibits a graceful exit [96].
Without giving details, the potential is parametrized by two constants and and reads, in Planck units [14]
(4.4) |
The expression for is analytically known and reads
(4.5) |
The field value at which Mukhanov inflation ends is obtained by solving and reads222Another solution exists for as the potential develops a maximum located at and there is a symmetry in these cases.
(4.6) |
These functions being exact solution of Eq. 2.7, they automatically satisfy all the integral constraints discussed in Section 3.1.
The present problem is to estimate knowing only and one can use Mukhanov inflation as a proxy. For instance, one can determine the value of and such that . If the slow-roll trajectories of Mukhanov inflation and of the model under scrutiny are close, then should also be their respective exact trajectories. As such, one can use as an approximation of the unknown .
In order to determine the two parameters and , one needs two equations, which cannot be provided by the integral constraints of Section 3.1 as they are already verified. Let us match the first and second derivatives of the potential’s logarithm towards the end of inflation by imposing
(4.7) |
Here, refers to the slow-roll approximated value, precisely obtained by solving , while the slow-roll velocity of Mukhanov inflation has still to be determined. From Eq. 4.4, one has
(4.8) |
where
(4.9) |
Taking the derivative of Eq. 4.8 with respect to gives
(4.10) |
Let us mention that the slow-roll approximated functions of Eqs. 4.8 and 4.10 are also related to the so-called potential slow-roll parameters and .
For a given inflationary potential , plugging Eqs. 4.8 and 4.10 into Eq. 4.7 gives a set of two algebraic equations for and that has to be solved numerically. Once the value of and are determined, the exact field value at which Mukhanov inflation ends is given by Eq. 4.6 and this will be taken as the estimator of the unknown .
4.3 End point correction
In Table 1, we give the exact numerical value of , the slow-roll approximated value , the trapezoid-approximated value and the Mukhanov-approximated , for the six inflationary models presented in Section 2.4. For the quite extreme case , we have not reported the value of as we have found more than one numerical solution for Eqs. 4.1 and 4.2 thereby preventing an easy determination of the best estimator. Let us notice that the value of , even if close to the exact one, lies in a domain for which the potential is not defined () and, as such, it is certainly not really useful. Let us stress, again, that numerically solving algebraic equations, such as Eqs. 4.1 and 4.2, or, Eq. 4.7, is orders of magnitude faster than numerically integrating Eq. 2.7 all along inflation. As can be seen in this table, the relative error in reference to is reduced by an order of magnitude using either the trapezoidal approximation or the Mukhanov inflation matching method instead of the traditional slow-roll value.
One can now define an additional correction to slow roll, including both the velocity correction of Section 3.3 and a better estimation of the end point field value, as
(4.11) | ||||
and a similar expression for the Mukhanov-approximated value
(4.12) |
The performance of and have been represented in Fig. 2, as green and magenta curves, compared to the traditional slow roll (red) and to the velocity corrected trajectory (blue). There is some improvement for while for the other model, using a more accurate field value for the end of inflation produces a slight overshoot for the total correction. As can be checked in Eq. 3.17, the next term () in the expansion has an opposite sign to the third one, , which is the error induced by using instead . As such, they somehow compensate and, for some models, the best option is to keep the simple velocity correction with in the trajectory. Notice, however, that a more accurate value for is always beneficial for estimating while the small overshooting is never large enough to spoil the velocity correction.
5 Conclusion
In this work, we have proposed new methods to improve the analytical observable predictions of the slow-roll single field models of inflation. Complementing most of the works in the literature that have been focused on the generation of cosmological perturbations during inflation, we have been focused here on a quite neglected aspect which concerns the accuracy at which the background field trajectory can be determined. As explained in the introduction, determining with precision the relation is crucial to correctly map wavenumbers today to wavenumbers during inflation. Moreover, because the reheating era lies in between the standard hot Big-Bang model eras and Cosmic Inflation, any uncertainties on will bias any inference made onto the kinematics of the reheating era.
One of the main results of this work is the exact expansion of Eq. 3.12, which has allowed us to present a simple and practical velocity correction to the usual slow-roll trajectory, as defined in Eq. 3.18. Adding this correction is trivial and, as shown in Fig. 2, immediately kills the absolute error on by an order of magnitude, for all the tested models. We have also discussed additional improvements to better determine the field value at which inflation ends. This end-point correction does not necessarily perform better than the velocity correction alone as it breaks some fortuitous compensation of some neglected higher order terms. Nonetheless, it never degrades significantly the velocity-corrected trajectory and always allows for a more accurate determination of .
Various other new results have been obtained along the course of searching for slow-roll improvements, such as the derivation of new integral constraints in Section 3.1, and, an exact, but still formal, new solution for the field trajectory when the functional is known, see Appendix A. Our work could be improved in various directions, such as estimating the next terms in the expansion of Eq. 3.18, or, devising more involved methods to determine . However, one should keep in mind that more involved analytical methods are relevant only if they remain simpler, and numerically much faster, than bruteforcely integrating Eq. 2.7. The present work may be precisely filling this niche.
Acknowledgements
We would like thank J. Martin and V. Vennin for enlightening discussions and for providing useful comments on the manuscript. This work is supported by the ESA Belgian Federal PRODEX Grant , the Wallonia-Brussels Federation Grant ARC . B. B. is publishing in the quality of ASPIRANT Research Fellow of the FNRS.
Appendix A Formal solution
Starting from Eq. 2.7, instead of trying to use the field value as a variable, one may switch to the number of e-folds and express the right-hand side as
(A.1) |
where the potential is now viewed as a . Plugging Eq. A.1 into Eq. 2.7, multiplying both sides by one gets a differential equation for
(A.2) |
This differential equation is a non-homogeneous Riccati equation and can be solved analytically [97]. Let us define the new “boost” function
(A.3) |
In terms of , Eq. A.2 considerably simplifies into
(A.4) |
This is a non-homogeneous linear differential equation which admits the exact solution
(A.5) |
where, as before, . This solution extends an approximated one derived in Ref. [18] under the “non-relativistic” assumption ().
References
- [1] A. A. Starobinsky, Spectrum of relict gravitational radiation and the early state of the universe, JETP Lett. 30 (1979) 682–685.
- [2] A. A. Starobinsky, A New Type of Isotropic Cosmological Models Without Singularity, Phys. Lett. B 91 (1980) 99–102.
- [3] A. H. Guth, The Inflationary Universe: A Possible Solution to the Horizon and Flatness Problems, Phys. Rev. D 23 (1981) 347–356.
- [4] A. D. Linde, A New Inflationary Universe Scenario: A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole Problems, Phys. Lett. B 108 (1982) 389–393.
- [5] A. Albrecht and P. J. Steinhardt, Cosmology for Grand Unified Theories with Radiatively Induced Symmetry Breaking, Phys. Rev. Lett. 48 (1982) 1220–1223.
- [6] A. D. Linde, Chaotic Inflation, Phys. Lett. B 129 (1983) 177–181.
- [7] V. F. Mukhanov and G. V. Chibisov, Quantum Fluctuations and a Nonsingular Universe, JETP Lett. 33 (1981) 532–535.
- [8] V. F. Mukhanov and G. V. Chibisov, The Vacuum energy and large scale structure of the universe, Sov. Phys. JETP 56 (1982) 258–265.
- [9] A. A. Starobinsky, Dynamics of Phase Transition in the New Inflationary Universe Scenario and Generation of Perturbations, Phys. Lett. B 117 (1982) 175–178.
- [10] A. H. Guth and S. Y. Pi, Fluctuations in the New Inflationary Universe, Phys. Rev. Lett. 49 (1982) 1110–1113.
- [11] S. W. Hawking, The Development of Irregularities in a Single Bubble Inflationary Universe, Phys. Lett. B 115 (1982) 295.
- [12] J. M. Bardeen, P. J. Steinhardt and M. S. Turner, Spontaneous creation of almost scale-free density perturbations in an inflationary universe, Phys. Rev. D 28 (1983) 679–693.
- [13] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Theory of cosmological perturbations. Part 1. Classical perturbations. Part 2. Quantum theory of perturbations. Part 3. Extensions, Phys. Rept. 215 (1992) 203–333.
- [14] J. Martin, C. Ringeval and V. Vennin, Encyclopædia Inflationaris, Phys. Dark Univ. 5-6 (2014) 75–235, [1303.3787].
- [15] D. Baumann and L. McAllister, Inflation and String Theory. Cambridge Monographs on Mathematical Physics. Cambridge University Press, 5, 2015, 10.1017/CBO9781316105733.
- [16] V. Vennin, K. Koyama and D. Wands, Encyclopædia curvatonis, JCAP 11 (2015) 008, [1507.07575].
- [17] Planck collaboration, Y. Akrami et al., Planck 2018 results. X. Constraints on inflation, Astron. Astrophys. 641 (2020) A10, [1807.06211].
- [18] D. Chowdhury, J. Martin, C. Ringeval and V. Vennin, Inflation after Planck: Judgement Day, Phys. Rev. D 100 (2019) 083537, [1902.03951].
- [19] J. Martin, C. Ringeval and V. Vennin, Cosmic Inflation at the Crossroads, 2404.10647.
- [20] CMB-S4 collaboration, K. N. Abazajian et al., CMB-S4 Science Book, First Edition, 1610.02743.
- [21] Simons Observatory collaboration, P. Ade et al., The Simons Observatory: Science goals and forecasts, JCAP 02 (2019) 056, [1808.07445].
- [22] M. Mallaby-Kay et al., The Atacama Cosmology Telescope: Summary of DR4 and DR5 Data Products and Data Access, Astrophys. J. Supp. 255 (2021) 11, [2103.03154].
- [23] Euclid collaboration, R. Scaramella et al., Euclid space mission: a cosmological challenge for the next 15 years, IAU Symp. 306 (2014) 375–378, [1501.04908].
- [24] Euclid collaboration, S. Ilić et al., Euclid preparation. XV. Forecasting cosmological constraints for the Euclid and CMB joint analysis, Astron. Astrophys. 657 (2022) A91, [2106.08346].
- [25] LSST Science, LSST Project collaboration, P. A. Abell et al., LSST Science Book, Version 2.0, 0912.0201.
- [26] LiteBIRD collaboration, E. Allys et al., Probing Cosmic Inflation with the LiteBIRD Cosmic Microwave Background Polarization Survey, PTEP 2023 (2023) 042F01, [2202.02773].
- [27] C. Ringeval, Fast Bayesian inference for slow-roll inflation, Mon. Not. Roy. Astron. Soc. 439 (Apr., 2014) 3253–3261, [1312.2347].
- [28] J. Martin, C. Ringeval and V. Vennin, Shortcomings of New Parametrizations of Inflation, Phys. Rev. D 94 (2016) 123521, [1609.04739].
- [29] H. Kurki-Suonio, P. Laguna and R. A. Matzner, Inhomogeneous inflation: Numerical evolution, Phys. Rev. D 48 (1993) 3611–3624, [astro-ph/9306009].
- [30] W. E. East, M. Kleban, A. Linde and L. Senatore, Beginning inflation in an inhomogeneous universe, JCAP 09 (2016) 010, [1511.05143].
- [31] K. Clough, R. Flauger and E. A. Lim, Robustness of Inflation to Large Tensor Perturbations, JCAP 05 (2018) 065, [1712.07352].
- [32] J. C. Aurrekoetxea, K. Clough, R. Flauger and E. A. Lim, The Effects of Potential Shape on Inhomogeneous Inflation, JCAP 05 (2020) 030, [1910.12547].
- [33] C. Joana and S. Clesse, Inhomogeneous preinflation across Hubble scales in full general relativity, Phys. Rev. D 103 (2021) 083501, [2011.12190].
- [34] C. Joana, Gravitational dynamics in Higgs inflation: Preinflation and preheating with an auxiliary field, Phys. Rev. D 106 (2022) 023504, [2202.07604].
- [35] M. Elley, J. C. Aurrekoetxea, K. Clough, R. Flauger, P. Giannadakis and E. A. Lim, Robustness of inflation to kinetic inhomogeneities, 2405.03490.
- [36] C. Joana, Beginning inflation in conformally curved spacetimes, 2406.00811.
- [37] A. A. Starobinsky and J. Yokoyama, Equilibrium state of a selfinteracting scalar field in the De Sitter background, Phys. Rev. D50 (1994) 6357–6368, [astro-ph/9407016].
- [38] A. A. Starobinsky, Stochastic de Sitter (inflationary) Stage in the Early Universe, in Field Theory, Quantum Gravity and Strings (H. J. de Vega & N. Sánchez, ed.), vol. 246 of Lecture Notes in Physics, Berlin Springer Verlag, p. 107, 1986. DOI.
- [39] V. Vennin and A. A. Starobinsky, Correlation Functions in Stochastic Inflation, Eur. Phys. J. C 75 (2015) 413, [1506.04732].
- [40] K. Ando and V. Vennin, Power spectrum in stochastic inflation, JCAP 04 (2021) 057, [2012.02031].
- [41] B. Blachier, P. Auclair, C. Ringeval and V. Vennin, Spatial curvature from super-Hubble cosmological fluctuations, Phys. Rev. D 108 (2023) 123510, [2302.14530].
- [42] K. Tokeshi and V. Vennin, Why does inflation look single field to us?, 2310.16649.
- [43] D. S. Salopek and J. R. Bond, Nonlinear evolution of long-wavelength metric fluctuations in inflationary models, Phys. Rev. D 42 (Dec, 1990) 3936–3962.
- [44] J. A. Adams, B. Cresswell and R. Easther, Inflationary perturbations from a potential with a step, Phys. Rev. D64 (2001) 123514, [astro-ph/0102236].
- [45] C. Ringeval, P. Brax, C. van de Bruck and A.-C. Davis, Boundary inflation and the wmap data, Phys.Rev. D73 (2006) 064035, [astro-ph/0509727].
- [46] A. Makarov, On the accuracy of slow-roll inflation given current observational constraints, Phys. Rev. D72 (2005) 083517, [astro-ph/0506326].
- [47] M. J. Mortonson, H. V. Peiris and R. Easther, Bayesian Analysis of Inflation: Parameter Estimation for Single Field Models, Phys.Rev. D83 (2011) 043505, [1007.4205].
- [48] L. C. Price, H. V. Peiris, J. Frazer and R. Easther, Designing and testing inflationary models with Bayesian networks, JCAP 02 (2016) 049, [1511.00029].
- [49] D. Seery, CppTransport: a platform to automate calculation of inflationary correlation functions, 1609.00380.
- [50] D. Werth, L. Pinol and S. Renaux-Petel, CosmoFlow: Python Package for Cosmological Correlators, 2402.03693.
- [51] A. Caravano, K. Inomata and S. Renaux-Petel, The Inflationary Butterfly Effect: Non-Perturbative Dynamics From Small-Scale Features, 2403.12811.
- [52] J. Martin and C. Ringeval, Inflation after WMAP3: Confronting the Slow-Roll and Exact Power Spectra to CMB Data, JCAP 08 (2006) 009, [astro-ph/0605367].
- [53] R. Easther and H. V. Peiris, Bayesian Analysis of Inflation II: Model Selection and Constraints on Reheating, Phys.Rev. D85 (2012) 103533, [1112.0326].
- [54] V. F. Mukhanov, L. A. Kofman and D. Y. Pogosian, Cosmological Perturbations in the Inflationary Universe, Phys. Lett. B 193 (1987) 427–432.
- [55] V. F. Mukhanov, Quantum Theory of Gauge Invariant Cosmological Perturbations, Sov. Phys. JETP 67 (1988) 1297–1302.
- [56] E. D. Stewart and D. H. Lyth, A more accurate analytic calculation of the spectrum of cosmological perturbations produced during inflation, Phys. Lett. B302 (1993) 171–175, [gr-qc/9302019].
- [57] E. D. Stewart, The Spectrum of density perturbations produced during inflation to leading order in a general slow roll approximation, Phys. Rev. D 65 (2002) 103508, [astro-ph/0110322].
- [58] J.-O. Gong and E. D. Stewart, The Density perturbation power spectrum to second order corrections in the slow roll expansion, Phys. Lett. B 510 (2001) 1–9, [astro-ph/0101225].
- [59] D. J. Schwarz, C. A. Terrero-Escalante and A. A. Garcia, Higher order corrections to primordial spectra from cosmological inflation, Phys. Lett. B 517 (2001) 243–249, [astro-ph/0106020].
- [60] S. M. Leach, A. R. Liddle, J. Martin and D. J. Schwarz, Cosmological parameter estimation and the inflationary cosmology, Phys. Rev. D 66 (2002) 023515, [astro-ph/0202094].
- [61] J. Choe, J.-O. Gong and E. D. Stewart, Second order general slow-roll power spectrum, JCAP 07 (2004) 012, [hep-ph/0405155].
- [62] D. J. Schwarz and C. A. Terrero-Escalante, Primordial fluctuations and cosmological inflation after WMAP 1.0, JCAP 0408 (2004) 003, [hep-ph/0403129].
- [63] P. Auclair and C. Ringeval, Slow-roll inflation at N3LO, Phys. Rev. D 106 (2022) 063512, [2205.12608].
- [64] J. Martin and D. J. Schwarz, Wkb approximation for inflationary cosmological perturbations, Phys. Rev. D67 (2003) 083512, [astro-ph/0210090].
- [65] S. Habib, K. Heitmann, G. Jungman and C. Molina-Paris, The inflationary perturbation spectrum, Phys. Rev. Lett. 89 (2002) 281301, [astro-ph/0208443].
- [66] S. Habib, A. Heinen, K. Heitmann, G. Jungman and C. Molina-Paris, Characterizing inflationary perturbations: The uniform approximation, Phys. Rev. D70 (2004) 083507, [astro-ph/0406134].
- [67] R. Casadio, F. Finelli, M. Luzzi and G. Venturi, Improved wkb analysis of cosmological perturbations, Phys. Rev. D71 (2005) 043517, [gr-qc/0410092].
- [68] R. Easther and J. T. Giblin, The Hubble slow roll expansion for multi field inflation, Phys. Rev. D72 (2005) 103505, [astro-ph/0505033].
- [69] F. Di Marco and F. Finelli, Slow-roll inflation for generalized two-field Lagrangians, Phys. Rev. D71 (2005) 123502, [astro-ph/0505198].
- [70] R. Casadio, F. Finelli, M. Luzzi and G. Venturi, Higher order slow-roll predictions for inflation, Phys. Lett. B625 (2005) 1–6, [gr-qc/0506043].
- [71] X. Chen, M.-x. Huang, S. Kachru and G. Shiu, Observational signatures and non-Gaussianities of general single field inflation, JCAP 0701 (2007) 002, [hep-th/0605045].
- [72] T. Battefeld and R. Easther, Non-Gaussianities in Multi-field Inflation, JCAP 0703 (2007) 020, [astro-ph/0610296].
- [73] W. H. Kinney and K. Tzirakis, Quantum modes in DBI inflation: exact solutions and constraints from vacuum selection, Phys. Rev. D77 (2008) 103517, [0712.2043].
- [74] S. Yokoyama, T. Suyama and T. Tanaka, Primordial Non-Gaussianity in Multi-Scalar Slow-Roll Inflation, JCAP 0707 (2007) 013, [0705.3178].
- [75] L. Lorenz, J. Martin and C. Ringeval, K-inflationary Power Spectra in the Uniform Approximation, Phys.Rev. D78 (2008) 083513, [0807.3037].
- [76] K. Tzirakis and W. H. Kinney, Non-canonical generalizations of slow-roll inflation models, JCAP 01 (2009) 028, [0810.0270].
- [77] N. Agarwal and R. Bean, Cosmological constraints on general, single field inflation, Phys. Rev. D79 (2009) 023503, [0809.2798].
- [78] T. Chiba and M. Yamaguchi, Extended Slow-Roll Conditions and Primordial Fluctuations: Multiple Scalar Fields and Generalized Gravity, JCAP 0901 (2009) 019, [0810.5387].
- [79] K. Ichikawa, T. Suyama, T. Takahashi and M. Yamaguchi, Non-Gaussianity, Spectral Index and Tensor Modes in Mixed Inflaton and Curvaton Models, Phys. Rev. D78 (2008) 023513, [0802.4138].
- [80] D. Langlois, S. Renaux-Petel, D. A. Steer and T. Tanaka, Primordial perturbations and non-Gaussianities in DBI and general multi-field inflation, Phys.Rev. D78 (2008) 063523, [0806.0336].
- [81] A. De Felice and S. Tsujikawa, Conditions for the cosmological viability of the most general scalar-tensor theories and their applications to extended Galileon dark energy models, JCAP 1202 (2012) 007, [1110.3878].
- [82] J. Martin, C. Ringeval and V. Vennin, K-inflationary Power Spectra at Second Order, JCAP 06 (2013) 021, [1303.2120].
- [83] J. Beltran Jimenez, M. Musso and C. Ringeval, Exact Mapping between Tensor and Most General Scalar Power Spectra, Phys. Rev. D88 (2013) 043524, [1303.2788].
- [84] A. Karam, T. Pappas and K. Tamvakis, Frame-dependence of higher-order inflationary observables in scalar-tensor theories, Phys. Rev. D 96 (2017) 064036, [1707.00984].
- [85] E. Bianchi and M. Gamonal, Primordial power spectrum at N3LO in effective theories of inflation, 2405.03157.
- [86] J. Martin and C. Ringeval, First CMB Constraints on the Inflationary Reheating Temperature, Phys. Rev. D 82 (2010) 023511, [1004.5525].
- [87] C. Ringeval, T. Suyama and J. Yokoyama, Magneto-reheating constraints from curvature perturbations, JCAP 09 (2013) 020, [1302.6013].
- [88] M. Hindmarsh and O. Philipsen, WIMP dark matter and the QCD equation of state, Phys. Rev. D71 (2005) 087302, [hep-ph/0501232].
- [89] G. Dvali and M. Redi, Phenomenology of Dark Sectors, Phys. Rev. D 80 (2009) 055001, [0905.1709].
- [90] C. Ringeval, The exact numerical treatment of inflationary models, Lect. Notes Phys. 738 (2008) 243–273, [astro-ph/0703486].
- [91] A. D. Linde, Hybrid inflation, Phys.Rev. D49 (1994) 748–754, [astro-ph/9307002].
- [92] A. R. Liddle, P. Parsons and J. D. Barrow, Formalizing the slow roll approximation in inflation, Phys. Rev. D 50 (1994) 7222–7232, [astro-ph/9408015].
- [93] V. Vennin, Horizon-Flow off-track for Inflation, Phys. Rev. D 89 (2014) 083526, [1401.2926].
- [94] J. Ellis, M. A. G. Garcia, D. V. Nanopoulos and K. A. Olive, Calculations of Inflaton Decays and Reheating: with Applications to No-Scale Inflation Models, JCAP 07 (2015) 050, [1505.06986].
- [95] L. Iacconi, M. Fasiello, J. Väliviita and D. Wands, Novel CMB constraints on the parameter in alpha-attractor models, JCAP 10 (2023) 015, [2306.00918].
- [96] V. Mukhanov, Quantum Cosmological Perturbations: Predictions and Observations, Eur. Phys. J. C 73 (2013) 2486, [1303.3925].
- [97] I. S. Gradshteyn, I. M. Ryzhik, A. Jeffrey and D. Zwillinger, Table of Integrals, Series, and Products. 2007.