Next Article in Journal
Higher Dimensional Lie Algebroid Sigma Model with WZ Term
Next Article in Special Issue
NuFIT: Three-Flavour Global Analyses of Neutrino Oscillation Experiments
Previous Article in Journal
Asymptotic Safety: Swampland or Wonderland?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalised Asymptotic Solutions for the Inflaton in the Oscillatory Phase of Reheating

by
Gabriel Álvarez
1,†,
Luis Martínez Alonso
1,*,† and
Elena Medina
2,†
1
Departamento de Física Teórica, Facultad de Ciencias Físicas, Universidad Complutense, 28040 Madrid, Spain
2
Departamento de Matemáticas, Facultad de Ciencias, Universidad de Cádiz, Puerto Real, 11510 Cádiz, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Universe 2021, 7(10), 390; https://doi.org/10.3390/universe7100390
Submission received: 15 September 2021 / Revised: 14 October 2021 / Accepted: 15 October 2021 / Published: 19 October 2021
(This article belongs to the Special Issue Studying the Universe from Spain)

Abstract

:
We determine generalised asymptotic solutions for the inflaton field, the Hubble parameter, and the equation-of-state parameter valid during the oscillatory phase of reheating for potentials that close to their global minima behave as even monomial potentials. For the quadratic potential, we derive a generalised asymptotic expansion for the inflaton with respect to the scale set by inverse powers of the cosmic time. For the quartic potential, we derive an explicit, two-term generalised asymptotic solution in terms of Jacobi elliptic functions, with a scale set by inverse powers of the square root of the cosmic time. In the general case, we find similar two-term solutions where the leading order term is defined implicitly in terms of the Gauss hypergeometric function. The relation between the leading terms of the instantaneous equation-of-state parameter and different averaged values is discussed in the general case. Finally, we discuss the physical significance of the generalised asymptotic solutions in the oscillatory regime and their matching to the appropriate solutions in the thermalization regime.
PACS:
98.80.Cq; 02.30.Mv

1. Introduction

The reheating period in inflationary cosmology is a transition period connecting the end of inflation to the radiation-dominated era [1,2,3,4,5,6,7,8]. During reheating, the energy of the inflaton field Φ transforms into relativistic particles. A basic model for the inflaton decay is given by the system of Equations [1,9,10]
Φ ¨ + ( 3 H + Γ ) Φ ˙ + V ( Φ ) = 0 ,
ρ ˙ γ = 4 ρ γ H + ( 1 + w Φ ) Γ ρ Φ ,
where dots denote derivatives with respect to the cosmic time t, and V ( Φ ) is the inflaton potential. The inflaton pressure, energy density, and equation of state (eos) parameter are given, respectively, by
p Φ = 1 2 Φ ˙ 2 V ( Φ ) ,
ρ Φ = 1 2 Φ ˙ 2 + V ( Φ ) ,
w Φ = p Φ / ρ Φ ,
Γ is the total inflaton decay rate, ρ γ is energy density of the radiation, and the Hubble parameter H satisfies the Friedmann–Lemaître Equation
H 2 = 1 3 M Pl 2 ( ρ Φ + ρ γ ) ,
with M Pl = 2.435 × 10 18   GeV being the reduced Planck mass.
The model (1)–(2) represents reheating in terms of an effective fluid with total pressure p = p Φ + p γ and density ρ = ρ Φ + ρ γ . The radiation eos p γ = ρ γ / 3 leads to an eos parameter during reheating given by [10]
w re = p ρ = 1 2 Φ ˙ 2 V ( Φ ) + 1 3 ρ γ 1 2 Φ ˙ 2 + V ( Φ ) + ρ γ .
The first phase of reheating can be described as an oscillatory phase around the true minimum of the inflaton potential, wherein the inflaton is not yet effectively coupled to the radiation fields, and the friction term Γ Φ ˙ in Equation (1) can be neglected. During this period, Equations (1) and (6) reduce to
Φ ¨ + 3 H Φ ˙ + V ( Φ ) = 0 ,
H 2 = 1 3 M Pl 2 1 2 Φ ˙ 2 + V ( Φ ) .
Conversely, in the last phase of reheating, the term 3 H Φ ˙ can be neglected with respect to Γ Φ ˙ in (1), leading to the simple equation
Φ ¨ + Γ Φ ˙ + V ( Φ ) = 0 .
During the oscillatory period, where ρ γ ρ Φ , w re reduces to w Φ , whereas, at the end of reheating, ρ ρ γ and w re tends to 1 / 3 , i.e., to the onset of the radiation era.
The aim of this paper is to determine generalised asymptotic solutions for the inflaton field Φ ( t ) and derived quantities valid during the oscillatory phase of reheating for potentials that close to their global minima behave as the even monomial potentials
V ( Φ ) = M 4 p Φ M Pl 2 p , ( p = 1 , 2 , ) .
Recently [8], it has been claimed that the origin of dark matter may reside in the process of reheating with potentials of the form (11) with p > 1 . We also mention that asymptotic solutions similar to those derived in this paper have been recently used to study the pre-inflationary (kinetic) and inflationary stages of inflaton models [11,12].
Turner [9], Mukhanov [5], and Rendall [13] emphasised that, during the oscillatory phase, the model is approximated in some averaged sense by a perfect fluid with a linear equation of state.
Since these generalised asymptotic solutions will be formally valid as the cosmic time t tends to infinity, we raise the question of its physical relevance in comparison with the exact solutions of the full Equations (1), (2) and (6), because the power-law-damped oscillatory approximation will certainly break down when H Γ and ρ γ ρ Φ , and the oscillations become exponentially damped.
However, typical values of the parameters for, e.g., the nondegenerate case p = 1 , are M = 3 × 10 3 M Pl and Γ = 1.375 × 10 9 GeV, and order-of-magnitude estimates and numerical calculations confirm that there is a large time interval (several e-folds) in which the asymptotic solutions in the oscillatory phase will be valid (see, for example, Figure 13 in Ref. [10]).
If we write
Φ ( t ) = 2 3 M Pl φ ( τ ) ,
where
τ = m t ,
and
m = 3 p 2 3 p / 2 M 2 M Pl ,
and substitute Equations (9) and (11) into (8), we find that the system is equivalent to the single non-linear differential equation
φ ( τ ) + ( φ ( τ ) 2 + φ ( τ ) 2 p ) 1 / 2 φ ( τ ) + p φ ( τ ) 2 p 1 = 0 .
Likewise, the Hubble parameter can be written in terms of φ as
H ( t ) = m 3 h ( τ ) ,
where, for later convenience, we introduced the reduced Hubble parameter
h ( τ ) = ( φ ( τ ) 2 + φ ( τ ) 2 p ) 1 / 2 .
Incidentally, by differentiation, the former equation and using (15), we find that
h ( τ ) = φ ( τ ) 2 .
As usual, primes denote derivatives of φ with respect to its argument.
By introducing polar coordinates in the ( φ , φ ) plane, Mukhanov [5] and Rendall [13] find the first two terms of the asymptotic expansion of the inflaton field for the quadratic potential corresponding to p = 1 in Equation (11). (In fact, Equation (24) of Ref. [13] corrects a mistake in Equation (5.45) of Ref. [5].) For general values of p and using different averaging procedures, the authors arrive at a common constant leading term for the averaged value
w Φ p 1 p + 1 , ( t ) ,
i.e., to the above-mentioned averaged linear equation of state. Thus, w Φ 0 for p = 1 and w Φ 1 / 3 for p = 2 , which correspond to a universe dominated by nonrelativistic matter and to a universe dominated by radiation, respectively. It was noted by Turner [9] that the amplitude of the inflaton oscillations decreases as τ 1 / p . However, for p 2 , we are not aware of any explicit results for the asymptotic form of the inflaton field.
In Section 2, we reconsider the nondegenerate harmonic potential discussed by Mukhanov [5] and Rendall [13]. First, we present a method to calculate a generalised asymptotic expansion for the inflaton field
φ ( τ ) n = 1 φ n ( τ ) τ n , ( τ ) ,
where the coefficients φ n are periodic (and therefore bounded) functions. To calculate these expansions we substitute the ansatz (20) into Equation (15) and solve recursively the system of differential equations for the φ n . The requisite that φ n be bounded is implemented by imposing the cancellation of resonant terms in these equations. As a byproduct, we will show how to calculate asymptotic expansions for derived quantities like H and w Φ . A formal proof of the existence of the generalised asymptotic expansion (20) is deferred to Appendix A.
A straightforward generalisation of the ansatz (20) does not work for potentials (11) with p > 1 because the resonant terms cannot be systematically eliminated past a certain order. Therefore, we look for generalised asymptotic formulas with a finite number of terms [14]. In Section 3, we discuss separately the quartic potential. This is a particularly relevant potential because the corresponding oscillatory phase of reheating is radiation-dominated, and the universe would be radiation-dominated since the end of inflation with an ensuing reduction of the uncertainty on the period were the observable perturbations of the universe were produced [15].
We obtain explicit asymptotic formulas for the inflaton field and derived quantities in terms of Jacobi elliptic functions, which have appeared recently in perturbative approaches to inflationary cosmology due to its relation to Hill’s Equation (16). In Section 4, we give a general derivation valid for any value of p, where the leading term of the asymptotic formula for the inflaton field is defined implicitly in terms of Gauss’ hypergeometric function, whereas the second term can be given explicitly in terms of the first. We also obtain the leading term of the instantaneous eos parameter w Φ and show how different averaged values of w Φ coincide with the standard result (19).
In Section 5, we discuss the physical significance of the generalised asymptotic solutions in the oscillatory regime and their matching to the appropriate solutions in the thermalization regime. In particular, we provide values for the cosmic time that marks the end of the oscillatory regime and of the Hubble parameter at that time. (The latter turns out to be independent of the degree 2 p of the potential close to its absolute minimum.) We also illustrate these results with a numerical calculation. The paper ends with some brief remarks on possible developments and applications of these generalised solutions.

2. The Harmonic Potential

In this section, we consider the nondegenerate harmonic potential discussed by Turner [9], Mukhanov [5], and Rendall [13], which corresponds to setting p = 1 in Equation (11), and for which Equation (15) reduces to
φ ( τ ) + ( φ ( τ ) 2 + φ ( τ ) 2 ) 1 / 2 φ ( τ ) + φ ( τ ) = 0 .
As we mentioned in the introduction, we look for a generalised asymptotic expansion of φ ( τ ) with respect to the asymptotic sequence τ n , ( n = 1 , 2 , 3 , ) (see Section 10.3 of [16]),
φ ( τ ) n = 1 φ n ( τ ) τ n , ( τ )
where the φ n ( τ ) are periodic functions. Substituting Equation (22) into Equation (21) and arranging the result by inverse powers of τ (note that all the derivatives of the φ n ( τ ) will also be periodic and, therefore, bounded), we find an infinite system of differential equations for the φ n ( τ ) . The first three, which are enough to show the pattern, are:
φ 1 + φ 1 = 0 ,
φ 2 + φ 2 = ( 2 ( ( φ 1 ) 2 + φ 1 2 ) 1 / 2 ) φ 1 ,
φ 3 + φ 3 = 4 φ 2 2 φ 1 φ 1 2 φ 2 2 φ 1 ( φ 1 ) 2 + φ 2 φ 1 φ 1 + 2 ( φ 1 ) 2 φ 2 φ 1 3 ( ( φ 1 ) 2 + φ 1 2 ) 1 / 2 .
Equation (23) is a harmonic oscillator, whose general solution can be written as
φ 1 ( τ ) = b cos ( τ τ 0 ) .
Substituting this solution into the right hand side of Equation (24), we find
φ 2 + φ 2 = b ( b 2 ) sin ( τ τ 0 ) .
Note that since the natural frequency of the the homogeneous equation φ 2 + φ 2 = 0 is ω = 1 , the inhomogeneous term sin ( τ τ 0 ) is resonant, as would be cos ( τ τ 0 ) . (Terms in sin ( ω ( τ τ 0 ) ) or cos ( ω ( τ τ 0 ) ) with ω 1 appearing below are not resonant.) Therefore, the solution for φ 2 ( τ ) would be unbounded unless we set b = 2 ( b = 0 leads to a trivial solution for φ 1 ( τ ) ), and we have again a harmonic oscillator with general solution
φ 2 ( τ ) = c cos ( τ τ 0 ) + d sin ( τ τ 0 ) .
Substituting this solution into the right-hand side of Equation (25) we find,
φ 3 + φ 3 = 1 + 2 d cos ( τ τ 0 ) cos ( 3 ( τ τ 0 ) ) .
The first term in the right-hand side is again resonant, and the solution for φ 3 ( τ ) would be unbounded unless we set d = 1 / 2 . The coefficient c remains free.
This recursive procedure can be easily programmed in a computer: at each stage, the function φ n ( τ ) is substituted into the right-hand side of the equation for φ n + 1 ( τ ) and resonant terms eliminated. It turns out that all the coefficients depend only on τ 0 and c. The first four coefficients are:
φ 1 ( τ ) = 2 cos ( τ τ 0 ) ,
φ 2 ( τ ) = c cos ( τ τ 0 ) 1 2 sin ( τ τ 0 ) ,
φ 3 ( τ ) = 1 16 8 c 2 + 9 cos ( τ τ 0 ) c 2 sin ( τ τ 0 ) + 1 8 cos ( 3 ( τ τ 0 ) ) ,
φ 4 ( τ ) = c 32 8 c 2 + 27 cos ( τ τ 0 ) 1 192 72 c 2 + 7 sin ( τ τ 0 )
+ 3 c 16 cos ( 3 ( τ τ 0 ) ) 5 32 sin ( 3 ( τ τ 0 ) ) .
In light of these four coefficients one might conjecture the general form
φ n ( τ ) = k = 1 , k odd n c n , k cos ( k ( τ τ 0 ) ) + d n , k sin ( k ( τ τ 0 ) ) ,
with d n , n = 0 for n odd. (In fact, the most efficient way to calculate the expansion is to use this ansatz and to equate coefficients of the inverse powers of τ .) A formal proof of the fact that Equation (35) is indeed true is deferred to Appendix A.
To illustrate the accuracy of this generalised asymptotic solution, in Figure 1 we plot the results of a numerical integration of the differential Equation (21) with initial conditions φ ( 1 ) = 0 , φ ( 1 ) = 1 , as well as the leading term of the asymptotic expansion (22). The numerical solution quickly tends to the asymptotic solution. The dashed lines are the envelopes ± 2 / τ that set the scale for the generalised asymptotic expansion (22). (The factor 2 is precisely the 2 in Equation (30) for φ 1 ( τ ) .) Other initial conditions lead to similar plots and are not represented to avoid clutter in the figure.
As a consequence of the main expansion (22), we can calculate the corresponding expansion for the reduced Hubble parameter (17),
h ( τ ) n = 1 h n ( τ ) τ n ,
whose first terms are
h 1 ( τ ) = 2 ,
h 2 ( τ ) = c + sin ( 2 ( τ τ 0 ) ) ,
h 3 ( τ ) = 1 4 ( 3 + 2 c 2 ) + 3 2 cos ( 2 ( τ τ 0 ) ) + c sin ( 2 ( τ τ 0 ) ) .
In Figure 2, we plot three typical asymptotic trajectories on the phase map of the differential Equation (21). Note that, for sufficiently large times, the approximation of the representative point to the origin is monotonic: in fact, for the quadratic potential, the distance of the representative point to the origin of the phase plane is simply the reduced Hubble constant h ( τ ) 2 / τ , as shown in Equation (37).
Using Equation (16), we find the expansion for the Hubble parameter
H ( t ) 2 3 t + c + sin ( 2 m ( t t 0 ) ) 3 m t 2 + 3 + 2 c 2 + 4 c sin ( 2 m ( t t 0 ) ) + 6 cos ( 2 m ( t t 0 ) ) 12 m 2 t 3 .
Here, m is the value of Equation (14) for p = 1 , i.e., m = 2 M 2 / M Pl , and τ 0 = m t 0 . In this case, the inflaton potential reduces to
V ( Φ ) = 1 2 m 2 Φ 2 ,
i.e., m is the p = 1 inflaton mass. Likewise, using
w Φ ( t ) = 1 2 3 H ˙ H 2 ,
we find the expansion for the equation of state parameter
w Φ ( t ) cos ( 2 m ( t t 0 ) ) + 3 sin ( 2 m ( t t 0 ) ) + sin ( 4 m ( t t 0 ) ) 2 m t + 1 16 m 2 t 2 ( 4 + 23 cos ( 2 m ( t t 0 ) ) + 12 c sin ( 2 m ( t t 0 ) ) + 20 cos ( 4 m ( t t 0 ) ) + 4 c sin ( 4 m ( t t 0 ) ) + 3 cos ( 6 m ( t t 0 ) ) ) + .
Note that these first two terms of the complete asymptotic expansion for w Φ ( t ) are independent of c. Note also that the average of the leading term is trivially zero.

3. The Quartic Potential

In this section, we discuss the quartic potential, which corresponds to setting p = 2 in Equation (15):
φ ( τ ) + ( φ ( τ ) 2 + φ ( τ ) 4 ) 1 / 2 φ ( τ ) + 2 φ ( τ ) 3 = 0 .
As was noted by Turner [9] and numerical calculations confirm, the solutions of this differential equation oscillate with a maximal amplitude that decreases as 1 / τ 1 / 2 . Therefore, one might try a generalised asymptotic solution similar to Equation (22) but with τ replaced by τ 1 / 2 on the right-hand side of that equation. It turns out that resonant terms cannot be systematically eliminated at fourth order, and, as a consequence, there is an undetermined constant at third order. Therefore, we look for a two-term generalised asymptotic solution of the form
φ ( τ ) φ 1 ( σ ) σ + φ 2 ( σ ) σ 2 + o ( 1 / σ 2 ) ,
where
σ = τ 1 / 2 ,
φ 1 ( σ ) and φ 2 ( σ ) are periodic functions of σ . By substituting the ansatz (45) into Equation (44), we find the following equations for φ 1 ( σ ) and φ 2 ( σ ) , which are the analogues of Equations (23) and (24):
φ 1 + 8 ( φ 1 ) 3 = 0 ,
φ 2 + 24 ( φ 1 ) 2 φ 2 = ( 3 ( ( φ 1 ) 2 + 4 φ 1 4 ) 1 / 2 ) φ 1 .
Equation (47) for φ 1 ( σ ) is a nonlinear differential equation whose solutions can be written in terms of the Jacobi elliptic function sn of modulus 1 [17]:
φ 1 ( σ ) = b sn ( 2 b ( σ σ 0 ) | 1 ) .
By substituting this solution into Equation (48), we find
φ 2 + 24 ( φ 1 ) 2 φ 2 = ( 3 2 b 2 ) φ 1 .
Incidentally, that the square root ( ( φ 1 ) 2 + 4 φ 1 4 ) 1 / 2 in Equation (48) reduces to the constant 2 b 2 is just the law of energy conservation for Equation (47) considered as a one-dimensional dynamical system. This fact, which also follows from the explicit form of the solution (49), will be essential in our discussion of the general case in the next section. Equation (50) is a nonhomogeneous linear differential equation with a non-constant coefficient. However, it is immediate to find two linearly independent solutions of the homogenous equation, namely u 1 ( σ ) = φ 1 ( σ ) and u 2 ( σ ) = φ 1 ( σ ) + σ φ 1 ( σ ) , and using the method of variation of constants, to find the general solution:
φ 2 ( σ ) = c φ 1 ( σ ) + d ( φ 1 ( σ ) + σ φ 1 ( σ ) ) + 3 2 b 2 3 σ φ 1 ( σ ) + σ 2 2 φ 1 ( σ ) .
The condition that φ 2 ( σ ) be bounded requires that
d = 0 , b = 3 2 ,
and we obtain explicit formulas for the coefficients φ 1 ( σ ) and φ 2 ( σ ) in the generalised asymptotic solution (45),
φ 1 ( σ ) = 3 2 sn ( 6 ( σ σ 0 ) | 1 ) ,
φ 2 ( σ ) = c φ 1 ( σ )
= 3 c cn ( 6 ( σ σ 0 ) | 1 ) dn ( 6 ( σ σ 0 ) | 1 ) ,
where cn and dn denote Jacobi elliptic functions [17]. Note that, again, the generalised asymptotic solution depends on two constants σ 0 and c, and that since σ = τ 1 / 2 = ( m t ) 1 / 2 , the frequency of the oscillations is not constant.
In Figure 3, we plot the results of a numerical integration of the differential Equation (44) with initial conditions φ ( 1 ) = 0 , φ ( 1 ) = 1 , as well as the leading term of the asymptotic expansion (45). In this case, the dashed lines are the envelopes ± ( 3 / ( 2 τ ) ) 1 / 2 that set the scale for the generalised asymptotic expansion (45), and the factor ( 3 / 2 ) 1 / 2 is the prefactor in Equation (53) for φ 1 ( τ ) .
As in the quadratic case, in Figure 4, we plot three typical asymptotic trajectories on the phase map of the differential Equation (44). Even for arbitrarily large times, the approximation of the representative point to the origin is not monotonic. (In this case, h ( τ ) does not admit the same geometric interpretation as in the quadratic case.)
Again, using Equations (16) and (17), we find the asymptotic formula for the Hubble parameter,
H ( t ) 1 2 t 3 m ( 6 m t ) 3 / 2 sn ( s ( t ) | 1 ) cn ( s ( t ) | 1 ) dn ( s ( t ) | 1 ) + ,
where
s ( t ) = 6 m t 6 m t 0 ,
and using Equation (42), the corresponding formula for the eos parameter,
w Φ ( t ) 1 2 sn ( s ( t ) | 1 ) 4 + .
At this order, w Φ ( t ) is independent of c. In the following section, we show that the leading term of the usual averages of this eos parameter is 1 / 3 .

4. The General Case

In this section, we discuss the generalised asymptotic solution of Equation (15) for any value of p. The results of this section include as particular cases the quadratic potential ( p = 1 ) and the quartic potential ( p = 2 ) studied in the previous sections, but the results obtained there were sharper than the results we obtain in this section. In the quadratic case, we could find explicit formulas and a complete asymptotic expansion. In the quartic case, we could still find explicit formulas although we had to limit ourselves to a two-term asymptotic solution. In this section, φ 1 ( σ ) is defined as the solution of an implicit equation.
The general equation, which we repeat here for convenience, is
φ ( τ ) + ( φ ( τ ) 2 + φ ( τ ) 2 p ) 1 / 2 φ ( τ ) + p φ ( τ ) 2 p 1 = 0 .
In analogy with the quartic case, we look for a two-term generalised asymptotic solution of the form
φ ( τ ) φ 1 ( σ ) σ + φ 2 ( σ ) σ 2 + o ( 1 / σ 2 ) ,
where
σ = τ 1 / p ,
and the conditions on the coefficients are the same as in the previous section. By substituting the ansatz (60) into Equation (59), we find the following equations for φ 1 and φ 2 :
φ 1 + p 3 ( φ 1 ) 2 p 1 = 0 ,
φ 2 + p 3 ( 2 p 1 ) ( φ 1 ) 2 p 2 φ 2 = ( p + 1 ( ( φ 1 ) 2 + p 2 φ 1 2 p ) 1 / 2 ) φ 1 .
(Again, primes denote derivatives of the functions with respect to their arguments, in this case σ .)
In general, Equation (62) cannot be solved in closed form. However, it is readily interpreted as the Newton’s equation of motion with respect to the time σ of a unit-mass particle with position φ 1 under the confining potential U ( φ 1 ) = p 2 ( φ 1 ) 2 p / 2 . The corresponding conserved energy is
E = 1 2 ( φ 1 ) 2 + p 2 2 ( φ 1 ) 2 p ,
and the solution φ 1 ( σ ) is implicitly defined by
σ σ 0 = 0 φ 1 d x 2 E p 2 x 2 p ,
or, equivalently, by
σ σ 0 = φ 1 2 E 2 F 1 1 2 , 1 2 p , 1 + 1 2 p , p 2 2 E ( φ 1 ) 2 p ,
where 2 F 1 denotes the Gauss’ hypergeometric function. Therefore, φ 1 ( σ ) represents a periodic motion on the interval ( 2 E / p 2 ) 1 / 2 p φ 1 ( 2 E / p 2 ) 1 / 2 p .
Equation (63) is, again, a nonhomogeneous linear differential equation with a non-constant coefficient. Two linearly independent solutions of the homogenous equation are u 1 ( σ ) = φ 1 ( σ ) and u 2 ( σ ) = φ 1 ( σ ) + ( p 1 ) σ φ 1 ( σ ) , and, using the method of variation of constants, we find the general solution for φ 2 ( σ ) in terms of φ 1 ( σ ) , namely
φ 2 ( σ ) = c φ 1 ( σ ) + d ( φ 1 ( σ ) + ( p 1 ) σ φ 1 ( σ ) ) + 1 2 E p + 1 σ φ 1 ( σ ) + p 1 2 σ 2 φ 1 ( σ ) .
The condition that φ 2 ( σ ) be bounded is that the terms proportional to σ and σ 2 vanish, which leads to
d = 0 , E = 1 2 ( p + 1 ) 2 ,
and to
φ 2 ( σ ) = c φ 1 ( σ ) .
We may also determine the first terms of the expansion of the reduced Hubble parameter (17),
h ( τ ) h 1 ( σ ) σ p + h 2 ( σ ) σ p + 1 + o ( 1 / σ p + 1 ) ,
where
h 1 ( σ ) = 1 p ( φ 1 ( σ ) 2 + p 2 φ 1 ( σ ) 2 p ) 1 / 2 ,
h 2 ( σ ) = 1 p 2 h 1 ( σ ) ( φ 1 ( σ ) φ 2 ( σ ) + p 3 φ 1 ( σ ) 2 p 1 φ 2 ( σ ) φ 1 ( σ ) φ 1 ( σ ) ) .
As a consequence of (62), (64), (68) and (69) we have
h 1 ( σ ) = p + 1 p ,
h 2 ( σ ) = 1 p ( p + 1 ) φ 1 ( σ ) φ 1 ( σ ) .
Then the leading term of the eos parameter is
w Φ ( t ) p 1 p + 1 2 p ( p + 1 ) 2 h 2 ( σ )
= 1 2 p 2 ( p + 1 ) 2 φ 1 ( σ ) 2 p .
In particular, the first terms of Equation (43) and of Equation (58) are, respectively, the cases p = 1 and p = 2 of Equation (76). That, at this order, w Φ ( t ) is independent of c is, thus, a general result.
As we mentioned in the introduction, different averaging procedures lead to the same leading contribution of the eos parameter (19). Equation (75) leads to the same result if we average over a period of the variable σ as an immediate consequence of this equation and of the periodicity of h 2 ( σ ) in σ (see Equation (74)):
w Φ 1 = 1 T σ 0 T σ w Φ ( σ ) d σ p 1 p + 1 .
It is also easy to derive the same result for other averages used in the literature. For example, Rendall [13] uses the averaging procedure defined by
w Φ 2 = t p Φ ( t ) d t t ρ Φ ( t ) d t .
In the variable τ ,
p Φ ( τ ) = M 4 p 2 3 p ( φ ( τ ) 2 φ ( τ ) 2 p ) ,
ρ Φ ( τ ) = M 4 p 2 3 p ( φ ( τ ) 2 + φ ( τ ) 2 p ) .
Using Equation (18), we find that,
τ φ ( τ ) 2 d τ = h ( τ ) p + 1 p 1 τ ,
where the asymptotic estimate follows from Equations (70) and (73). Furthermore, using the following consequence of Equations (15) and (17),
p φ ( τ ) 2 p = h ( τ ) 1 2 φ ( τ ) 2 φ ( τ ) 2 d d τ φ ( τ ) φ ( τ ) + 1 2 h ( τ ) φ ( τ ) 2 ,
we note that, as τ
h ( τ ) 1 + 1 2 φ ( τ ) 2 = O ( 1 / τ ) ,
φ ( τ ) φ ( τ ) = O ( 1 / τ 1 + 1 / p ) ,
and
τ φ ( τ ) 2 φ ( τ ) 2 d τ = O ( 1 / τ 1 + 2 / p ) .
Therefore,
τ φ ( τ ) 2 p d τ p + 1 p 2 1 τ ,
and from (81) and (86), we obtain that
w Φ 2 p 1 p + 1 .

5. Connection between the Oscillatory Period and the Thermalization Period

As we said in the introduction, the oscillatory phase of reheating is an intermediate period between the end of inflation and the beginning of thermalization. Let us denote by t end the cosmic time corresponding to the end of inflation and by t osc the cosmic time corresponding to the end of the oscillatory period. There are several criteria to determine the value of the inflaton Φ at the end of inflation, which typically give results differing in factors of order unity. For concreteness we quote Equations (2.23) and (2.25) of Ellis et al. [18] for the monomial potentials (11), which, with our conventions, lead to
Φ ( t end ) = 4 p 1 2 2 M Pl ,
Φ ˙ ( t end ) = 4 p 1 2 2 p M 2 p .
However, these (or similar) values cannot be used directly to relate the parameters σ 0 (or t 0 ) and c in the generalised asymptotic solutions (45) and (60) to the parameters of the model M and p for at least two reasons: First, the translational invariance of the inflaton models and the existence of an attractor in the phase plane [5,19] make it physically dubious relating t 0 to a specific moment of physical significance; second, as discussed in Refs [16,20], in general successive terms in generalised asymptotic expansions are not successively smaller at every time: they essentially fix the overall scale of the oscillatory process. (This is the case of Φ ( t ) and Φ ˙ ( t ) .)
However, suitable truncations of the generalised expansions may be free or the ambiguities discussed in the previous paragraph and can be used to link the parameter c to the values of physical magnitudes at t 0 . For concreteness, consider the first two terms of the expansion (40) for H ( t ) in the harmonic case, which we repeat here for convenience:
H ( t ) 2 3 t + c + sin ( 2 m ( t t 0 ) ) 3 m t 2 .
The idea is that the leading term in Equation (90) is not oscillating and, therefore, dominates the following terms at any large time, and that c in the next term is also not multiplied by an oscillating function. Therefore, setting t = t 0 , we find that
c = 3 m t 0 2 H ( t 0 ) 2 m t 0 .
Finally, let us briefly discuss the connection between the generalised asymptotic solutions in the oscillatory period to suitable asymptotic solutions in the thermalization period. To this end, we will denote by t osc the end of the oscillatory period, and by Φ I and Φ II the respective asymptotic forms of the inflaton field.
We consider first the harmonic case. From (12) and (30), we have that the leading term of the asymptotic expansion (22) is given by
Φ I ( t ) = 2 2 3 M Pl cos ( m ( t t 0 ) ) m t .
Furthermore, Equation (10) for the inflaton in the thermalization regime reduces in the harmonic case to
Φ ¨ + Γ Φ ˙ + m 2 Φ = 0 ,
whose exact solution can be written as
Φ ( t ) = A e Γ 2 t cos ( ω ( t t 1 ) ) ,
where
ω = m 2 Γ 2 4 ,
and where A and t 1 are arbitrary constants. In addition, for typical situations m Γ and, therefore, ω m . Thus, we take
Φ II ( t ) = A e Γ 2 t cos ( m ( t t 1 ) ) .
Again, since the expansion (22) is a generalised asymptotic series, by the same arguments used above it is physically meaningful to match the envelopes and the first derivatives of the envelopes of Φ I and Φ II at t = t osc , instead of matching instantaneous values at a certain time. Thus, we have
2 2 3 M Pl m t osc = A e Γ 2 t osc ,
2 2 3 M Pl m t osc 2 = A Γ 2 e Γ 2 t osc .
Hence, we obtain
t osc = 2 Γ ,
and
A = e 3 M Pl M 2 Γ .
Since ω m , we can also set t 1 = t 0 . Moreover, if we approximate the Hubble parameter H by the first term of (40), we obtain the following estimate of H osc H ( t osc ) ,
H osc Γ 3 .
To illustrate graphically the result (99) for t osc , in Figure 5, we present a numerical calculation similar to that of Figure 13 in Ref. [10]. More concretely, we performed a numerical integration of Equations (1)–(6) with the same parameters quoted in the introduction. In Figure 5a, we show the eos parameter (7) as a function of the e-fold number N. The dashed vertical line marks the matching point (99). Perhaps more informative is Figure 5b, in which we show the energy densities of the inflaton (4) and of the radiation as a function of the e-fold number. Again, the dashed vertical line marks the matching point (99), which turns out to be 0.38 e-folds to the right of the point where the energy densities become equal.
For general p, the first term of the asymptotic expansion (22) is given by
Φ I ( t ) = 2 3 M Pl φ 1 ( σ ) σ ,
where σ = ( m t ) 1 / p and φ 1 ( σ ) is a periodic function of σ , which satisfies Equation (62). As a consequence of Equations (64) and (68), φ 1 has a maximal amplitude given by ( 1 + 1 / p ) 1 / p . Then, from (102), it follows that Φ I is a damped oscillatory function of t with envelope given by
Φ I ( env ) ( t ) = 2 3 1 + 1 p 1 / p M Pl ( m t ) 1 / p .
Furthermore, for general p, the Equation (10) for the thermalization regime is a nonlinear Liénard differential equation of the form
Φ ¨ + Γ Φ ˙ + 2 M 4 M Pl 2 p Φ 2 p 1 = 0 ,
which describes a nonlinear oscillator with linear damping. Cveticanin [19] proved that the approximate solution of this equation describes a vibration with decreasing amplitude of the form (cf. Equations (72) and (73) in Ref. [19])
Φ II ( env ) ( t ) = A e Γ p + 1 t .
If we match the envelopes Φ I ( env ) and Φ II ( env ) and their first derivatives at t osc , we obtain at once
t osc = 1 + 1 p 1 Γ ,
and
A = 2 3 M Pl e Γ m 1 / p
= p 3 1 / ( 2 p ) M Pl 1 + 1 / p M 2 / p ( e Γ ) 1 / p .
In addition, if we approximate the Hubble parameter H by the first term of (70), we obtain the following p-independent estimate of H osc
H osc Γ 3 .
Finally, note that the matching of the oscillatory factor for p > 1 , although feasible, is complicated by the fact that the oscillations are not isochronous [19].

6. Conclusions

We analysed the asymptotic properties of the dynamical equation of the inflaton during the oscillatory phase of reheating to extend the existing results for the quadratic case, to provide new explicit results for the quartic case, and implicit results (that include the former as particular cases) for potentials that close to their global minima behave as even monomial potentials.
We paid special attention to the derivation of explicit expressions for the instantaneous eos parameter w Φ , whose averages turn out to be in agreement with the averaged values found in the literature. Several studies assumed that, during part or the whole reheating period, the constant value w rhe w Φ , thus, neglecting the coupling between the inflaton field and the radiation [10,21,22]. A more realistic yet manageable approach is to use a range of values for w rhe [23,24]. For instance, a typical range of values for p = 1 is 0 = w Φ w rhe 0.25 . In this spirit, our results may be useful to formulate constraints on effective eos parameters w rhe based on the instantaneous values w Φ ( t ) rather than on the averaged values w Φ .
Finally, as a line of future development, we mention the possible existence of a finer scale than τ 1 / p that might permit the calculation of full asymptotic expansions for the inflaton, instead of the truncated asymptotic formulas derived in the present paper for the degenerate potentials.

Author Contributions

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

Funding

This research was funded by the Spanish Ministerio de Economía y Competitividad grant number FIS2015-63966-P and PGC2018-094898-B-I00.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

In this appendix, we prove that the ansatz (35) leads to a recurrence procedure that can be carried out to arbitrary order. To this end, first note that (for any p) the main Equation (15) and the definition (17) for h ( τ ) imply that
φ ( τ ) + h ( τ ) φ ( τ ) + p φ ( τ ) 2 p 1 = 0 ,
h ( τ ) = φ ( τ ) 2 .
Conversely, since the system (A1)–(A2) implies that
d d τ h 2 ( ( φ ) 2 + φ 2 p ) = 0 ,
and φ , φ and h tend to zero, it follows that the system (A1)–(A2) and the main Equation (15) are equivalent [13].
Hereafter, we set p = 1 . In addition to the expansions (22) for φ ( τ ) , and (36) for h ( τ ) , it is convenient to introduce an expansion for φ ( τ ) ,
φ ( τ ) = n = 1 ψ n ( τ ) τ n ,
where
ψ n ( τ ) = φ n ( τ ) ( n 1 ) φ n 1 ( τ ) ,
and we set φ 0 ( τ ) 0 . Then, the system (A1)–(A2) leads to the following system of differential equations:
φ n + φ n = ( n 1 ) ( 2 φ n 1 ( n 2 ) φ n 2 ) l = 1 n 1 ψ l h n l ,
h n = ( n 1 ) h n 1 l = 1 n 1 ψ l ψ n l .
For n = 1 , these equations are
φ 1 + φ 1 = 0 ,
h 1 = 0 ,
and their general solutions are
φ 1 ( τ ) = c 1 , 1 cos ( τ τ 0 ) ,
h 1 ( τ ) = f 1 , 0 ,
respectively, where c 1 , 1 and f 1 , 0 are arbitrary real constants.
For n 2 , the solution of (A6)–(A7) can be written as
φ n ( τ ) = c n , 1 cos ( τ τ 0 ) + d n , 1 sin ( τ τ 0 ) + φ n ( p ) ( τ ) ,
h n ( t ) = f n , 0 + h n ( p ) ( τ ) ,
where φ n ( p ) ( τ ) and h n ( p ) ( τ ) are particular solutions of (A6) and (A7), respectively. To avoid secular terms in the functions φ n ( p ) ( τ ) and h n ( p ) ( τ ) , we have to ensure that the inhomogeneous part of (A6) is free of terms proportional to cos ( τ τ 0 ) or sin ( τ τ 0 ) , and that the inhomogeneous part of (A7) does not contain a constant term.
We now prove that this process of eliminating secular terms yields an expansion (22) for φ ( τ ) with coefficients of the form (35), and that the corresponding expansion (36) is of the form
h n ( τ ) = k = 0 , k even n f n , k cos ( k ( τ τ 0 ) ) + g n , k sin ( k ( τ τ 0 ) ) ,
with f n , n = 0 for even n.
Thus, for n = 2 the Equations (A6) and (A7) are
φ 2 + φ 2 = ( 2 h 1 ) φ 1
h 2 = h 1 ( φ 1 ) 2 .
Then, taking into account (A10), the inhomogeneous parts of (A15) and (A16) become
( 2 h 1 ( τ ) ) φ 1 ( τ ) = ( f 1 , 0 2 ) c 1 , 1 sin ( τ τ 0 ) ,
h 1 ( τ ) φ 1 ( τ ) 2 = f 1 , 0 c 1 , 1 2 2 + c 1 , 1 2 2 cos ( 2 ( τ τ 0 ) ) .
Therefore, the particular solutions φ 2 ( p ) ( τ ) and h 2 ( p ) ( τ ) of (A15) and (A16) are free of the resonant terms provided
f 1 , 0 = 2 , c 1 , 1 = 2 .
Then, we obtain
φ 1 ( τ ) = 2 cos ( τ τ 0 ) ,
h 1 ( τ ) = 2 ,
and
φ 2 ( τ ) = c 2 , 1 cos ( τ τ 0 ) + d 2 , 1 sin ( τ τ 0 ) ,
h 2 ( t ) = f 2 , 0 + sin ( 2 ( τ τ 0 ) ) .
The same procedure for n = 3 leads to
c 2 , 1 = f 2 , 0 = c , d 2 , 1 = 1 2 ,
where c is a free parameter.
We proceed by induction and assume that, for any given integer n 3 , we have already obtained the coefficients φ k and h k with k = 1 , , n 1 , of the form (35) and (A14).
The secular terms in Equation (A7) for h n come from
( n 1 ) h n 1 2 ψ 1 ψ n 1 l = 2 n 2 ψ l ψ n l ,
and vanish provided that
2 c n 1 , 1 + ( n 1 ) f n 1 , 0 = R n 2 ,
where R n 2 denotes a sum of terms depending on the coefficients c k , j d k , j with k n 2 .
Next, we note that the secular terms in Equation (A6) for φ n come from the terms
( n 1 ) 2 φ n 1 ( n 2 ) φ n 2 l = 1 n 1 ψ l h n l ,
and are given by
2 ( n 2 ) c n 1 , 1 sin ( τ τ 0 ) d n 1 , 1 cos ( τ τ 0 ) + 2 f n 1 , 0 sin ( τ τ 0 ) S n 2 sin ( τ τ 0 ) T n 2 cos ( τ τ 0 ) ,
where S n 2 and T n 2 denote sums of terms depending on the coefficients c k , j , d k , j , f k , j , g k , j , with k n 2 and on f n 1 , 2 , g n 1 , 2 . Therefore, to eliminate the secular terms in cos ( τ τ 0 ) and sin ( τ τ 0 ) , we must set
2 ( n 2 ) d n 1 , 1 = T n 2 ,
and
2 ( n 2 ) c n 1 , 1 + 2 f n 1 , 0 = S n 2 ,
respectively. Equation (A29) determines d n 1 , 1 in terms of the coefficients c k , j d k , j , f k , j , g k , j with k n 2 and on g n 1 , 2 . Furthermore, Equations (A26) and (A30) are a linear system for c n 1 , 1 and f n 1 , 0 whose coefficient matrix has determinant 2 n ( n 3 ) . Thus, this linear system is undetermined for n = 3 (with a free unknown c = c 2 , 1 = f 2 , 0 ), and it has a unique solution for n > 3 . In this way, it is clear that all the coefficients c n , 1 , d n , 1 , and f n , 0 , with the exception of c 2 , 1 = f 2 , 0 , are determined in terms of c.
Due to our induction hypothesis, from the structure of the right-hand side of the Equation (A7), it is straightforward to deduce that h n is also of the form (A14). For example, if n is odd ( n = 2 r + 1 ), then h n 1 is a linear combination of cos ( 2 k ( τ τ 0 ) ) and sin ( 2 k ( τ τ 0 ) ) with k = 0 , 1 , , r 1 and k = 1 , , r . Similarly, the products ψ l ψ n l with l = 1 2 r in Equation (A7) are linear combinations of cos ( 2 k ( τ τ 0 ) ) and sin ( 2 k ( τ τ 0 ) ) with k = 0 , 1 , , r . Therefore, h n is of the form (A14).
A similar argument proves the ansatz (35) for φ n .

References

  1. Albrecht, A.; Steinhardt, P.J.; Turner, M.S.; Wilczek, F. Reheating an Inflationary Universe. Phys. Rev. Lett. 1982, 48, 1437. [Google Scholar] [CrossRef]
  2. Kofman, L.; Linde, A.; Starobinsky, A.A. Reheating and Inflation. Phys. Rev. Lett. 1994, 73, 3195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Shtanov, Y.; Traschen, J.; Brandenberg, R. Universe reheating after inflation. Phys. Rev. D 1995, 51, 5438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Kofman, L.; Linde, A.; Starobinsky, A.A. Towards the theory of reheating after inflation. Phys. Rev. D 1997, 56, 3258. [Google Scholar] [CrossRef] [Green Version]
  5. Mukhanov, V. Physical Foundations of Cosmology; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  6. Bassett, B.A.; Tsujikawa, S.; Wands, D. Inflation dynamics and reheating. Rev. Mod. Phys. 2006, 78, 537–589. [Google Scholar] [CrossRef] [Green Version]
  7. Baumann, D. TASI lectures on inflation. arXiv 2009, arXiv:0907.5424. [Google Scholar]
  8. Garcia, M.A.G.; Kaneta, K.; Mambrini, Y.; Olive, K.A. Reheating and post-inflationary production of dark matter. Phys. Rev. D 2020, 101, 123507. [Google Scholar] [CrossRef]
  9. Turner, M.S. Coherent scalar-field oscillations in an expanding universe. Phys. Rev. D 1983, 28, 1243–1247. [Google Scholar] [CrossRef]
  10. Martin, J.; Ringeval, C. First CMB constraints on the inflationary reheating temperature. Phys. Rev. D 2010, 82, 023511. [Google Scholar] [CrossRef] [Green Version]
  11. Álvarez, G.; Martínez Alonso, L.; Medina, E.; Vázquez, J.L. Separatrices in the Hamilton-Jacobi formalism of inflaton models. J. Math. Phys. 2020, 61, 043501. [Google Scholar] [CrossRef] [Green Version]
  12. Medina, E.; Martínez Alonso, L. Kinetic dominance and psi series in the Hamilton-Jacobi formulation of inflaton models. Phys. Rev. D 2020, 102, 103517. [Google Scholar] [CrossRef]
  13. Rendall, A.D. Late-time oscillatory behaviour for self-gravitating scalar fields. Class. Quantum Grav. 2007, 24, 667–677. [Google Scholar] [CrossRef] [Green Version]
  14. Miller, P.D. Applied Asymptotic Analysis; AMS: Providence, RI, USA, 2006. [Google Scholar]
  15. 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] [Green Version]
  16. Olver, F.W.J. Asymptotics and Special Functions; A K Peters: Wellesley, MA, USA, 1997. [Google Scholar]
  17. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions; Dover: New York, NY, USA, 1972. [Google Scholar]
  18. Ellis, J.; Garcia, M.A.; Nanopoulos, D.V.; Olive, K.A. Calculations of inflaton decays and reheating: With applications to no-scale inflation models. J. Cosmol. Astropart. Phys. 2015, 2015, 50. [Google Scholar] [CrossRef]
  19. Cveticanin, L. Oscillators with nonlinear elastic and damping forces. Comput. Math. Appl. 2011, 62, 1745. [Google Scholar] [CrossRef] [Green Version]
  20. Wong, R. Asymptotic Approximations of Integrals; SIAM: Philadelphia, PA, USA, 2001. [Google Scholar]
  21. Ueno, Y.; Yamamoto, K. Constraints on α-attractor inflation and reheating. Phys. Rev. D 2016, 93, 083524. [Google Scholar] [CrossRef] [Green Version]
  22. Drewes, M.; Kang, J.U.; Mun, U.R. CMB constraints on the inflaton couplings and reheating temperature in α-attractor inflation. J. High Energy Phys. 2017, 2017, 72. [Google Scholar] [CrossRef] [Green Version]
  23. Cook, J.L.; Dimastrogiovanni, E.; Easson, D.A.; Kraussa, L.M. Reheating predictions in single inflation. J. Cosmol. Astropart. Phys. 2015, 2015, 47. [Google Scholar] [CrossRef]
  24. Muñoz, J.B.; Kamionkowski, M. Equation-of-state parameter for reheating. Phys. Rev. D 2015, 91, 043521. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Numerical integration (blue line) with initial conditions φ ( 1 ) = 0 , φ ( 1 ) = 1 and leading term of the asymptotic expansion (brown line) for the differential Equation (21) corresponding to the quadratic potential. The dashed lines are the envelopes ± 2 / τ that set the scale for the generalised asymptotic expansion (22).
Figure 1. Numerical integration (blue line) with initial conditions φ ( 1 ) = 0 , φ ( 1 ) = 1 and leading term of the asymptotic expansion (brown line) for the differential Equation (21) corresponding to the quadratic potential. The dashed lines are the envelopes ± 2 / τ that set the scale for the generalised asymptotic expansion (22).
Universe 07 00390 g001
Figure 2. Three typical asymptotic trajectories on the phase map of the differential Equation (21) corresponding to the quadratic potential. Note that, for sufficiently large times, the approximation of the representative point to the origin is monotonic.
Figure 2. Three typical asymptotic trajectories on the phase map of the differential Equation (21) corresponding to the quadratic potential. Note that, for sufficiently large times, the approximation of the representative point to the origin is monotonic.
Universe 07 00390 g002
Figure 3. Numerical integration (blue line) with initial conditions φ ( 1 ) = 0 , φ ( 1 ) = 1 and leading term of the asymptotic expansion (brown line) for the differential Equation (44) corresponding to the quartic potential. The dashed lines are the envelopes ± ( 3 / ( 2 τ ) ) 1 / 2 that set the scale for the generalised asymptotic solution (45).
Figure 3. Numerical integration (blue line) with initial conditions φ ( 1 ) = 0 , φ ( 1 ) = 1 and leading term of the asymptotic expansion (brown line) for the differential Equation (44) corresponding to the quartic potential. The dashed lines are the envelopes ± ( 3 / ( 2 τ ) ) 1 / 2 that set the scale for the generalised asymptotic solution (45).
Universe 07 00390 g003
Figure 4. Three typical asymptotic trajectories on the phase map of the differential Equation (44) corresponding to the quartic potential. Even for arbitrarily large times the approximation of the representative point to the origin is not monotonic.
Figure 4. Three typical asymptotic trajectories on the phase map of the differential Equation (44) corresponding to the quartic potential. Even for arbitrarily large times the approximation of the representative point to the origin is not monotonic.
Universe 07 00390 g004
Figure 5. Numerical integration of Equations (1)–(6) with the parameters quoted in the introduction. (a) The eos parameter w re as a function of the e-fold number N. (b) Inflaton and radiation energy densities. In both cases, the dashed vertical line marks the matching point t osc = 2 / Γ .
Figure 5. Numerical integration of Equations (1)–(6) with the parameters quoted in the introduction. (a) The eos parameter w re as a function of the e-fold number N. (b) Inflaton and radiation energy densities. In both cases, the dashed vertical line marks the matching point t osc = 2 / Γ .
Universe 07 00390 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Álvarez, G.; Martínez Alonso, L.; Medina, E. Generalised Asymptotic Solutions for the Inflaton in the Oscillatory Phase of Reheating. Universe 2021, 7, 390. https://doi.org/10.3390/universe7100390

AMA Style

Álvarez G, Martínez Alonso L, Medina E. Generalised Asymptotic Solutions for the Inflaton in the Oscillatory Phase of Reheating. Universe. 2021; 7(10):390. https://doi.org/10.3390/universe7100390

Chicago/Turabian Style

Álvarez, Gabriel, Luis Martínez Alonso, and Elena Medina. 2021. "Generalised Asymptotic Solutions for the Inflaton in the Oscillatory Phase of Reheating" Universe 7, no. 10: 390. https://doi.org/10.3390/universe7100390

APA Style

Álvarez, G., Martínez Alonso, L., & Medina, E. (2021). Generalised Asymptotic Solutions for the Inflaton in the Oscillatory Phase of Reheating. Universe, 7(10), 390. https://doi.org/10.3390/universe7100390

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