Some Clarifications on the Duration of Inflation in Loop Quantum Cosmology

TL;DR

This paper clarifies how various assumptions influence the predicted duration of inflation in loop quantum cosmology, emphasizing that the quantum regime's background dynamics have limited impact on inflation length.

Contribution

It provides elementary clarifications on the assumptions affecting inflation duration predictions in loop quantum cosmology, highlighting the limited role of quantum background dynamics.

Findings

Inflation duration is not significantly affected by quantum regime background dynamics.

The number of e-folds in inflation is constrained by specific hypotheses.

Clarifications help understand the robustness of inflation predictions in loop quantum cosmology.

Abstract

The prediction of a phase of inflation whose number of e-folds is constrained is an important feature of loop quantum cosmology. This work aims at giving some elementary clarifications on the role of the different hypotheses leading to this conclusion. We show that the duration of inflation does not depend significantly on the modified background dynamics in the quantum regime.

Figures

Equations

$$H^2=\frac{\rho }{3}\left(1-\frac{\rho }{{\rho }_{\text{B}}}\right),$$

$$\ddot{\varphi }+3H\dot{\varphi }+m^2\varphi =0.$$

$$\mu (N)=\mathcal{C}m{N}^{-\frac{1}{2}}\exp (-3N)\{1+1/(6N)\},$$

$${\varphi }_{\text{I}}={\varphi }_{\text{0}}+\mathrm{sgn}({\dot{\varphi }}_{\text{0}})\sqrt{(2/3)}\mathrm{Arcsinh}\left(\Gamma \sqrt{2/\ln (z)}\right),$$

$$\mu (N)=1-m\sqrt{(2N/{\rho }_{\text{0}})},$$

$${\varphi }_{\text{B}}^{\mathrm{GT}}\equiv \sqrt{(2/3)}\ln \left(2\Gamma /\sqrt{\ln \Gamma }\right).$$

$$x^{\text{GR}}$$

$$x^{\text{LQC}}$$

Full text

Loop quantum gravity (LQG) is a nonperturbative and background-independent quantization of general relativity (GR). It relies on the Sen-Ashtekar-Barbero variables, that is SU(2) valued connections and conjugate densitized triads. The quantization is obtained using holonomies of the connections and fluxes of the densitized triads. Loop quantum cosmology (LQC) is an effective theory based on a symmetry reduced version of LQG. In LQC, the big bang is believed to be replaced by a bounce due to repulsive quantum geometrical effects (see Ashtekar and Singh (2011) for a review). For the flat homogeneous and isotropic background cosmology that we consider in this work, the effective LQC-modified Friedmann equation is

[TABLE]

where $H\equiv(\dot{a}/a)$ is the Hubble parameter, $\rho$ is the total energy density and ${\rho}_{\scriptscriptstyle{\textrm{B}}}$ is the critical density at the bounce (expected to be of the order of the Planck density). The dot refers to a coordinate time derivative.

We assume that the dominating energy component in the early universe is a scalar field $\phi$, with potential $V=\tfrac{1}{2}m^{2}\phi^{2}$. The total energy density can be written as $\rho=\tfrac{1}{2}\dot{\phi}^{2}+V$. Based on cosmic microwave background (CMB) measurements and under most reasonable assumptions for the length of observable inflation (between horizon exit of the pivot scale and the end of the inflationary phase), one obtains $m\simeq 10^{-6}{m}_{\scriptscriptstyle{\textrm{Pl}}}$. The equation of motion for the scalar field is

[TABLE]

There are different ways to statistically estimate the duration of inflation in this framework.

At a fixed energy density, ${\rho}_{\scriptscriptstyle{\textrm{0}}}$, one can first ask the following question: for a given number of e-folds $N$, what is the fraction of trajectories, i.e. solutions to Eq.(2), that lead to a phase of slow-roll inflation lasting more than $N$ e-folds? It should be noticed that the set of trajectories can be parametrized by $\{{a}_{\scriptscriptstyle{\textrm{0}}},{\phi}_{\scriptscriptstyle{\textrm{0}}}\}$. As the energy density has been fixed, the initial time derivative of the scalar field, ${\dot{\phi}}_{\scriptscriptstyle{\textrm{0}}}$, is determined in terms of ${\rho}_{\scriptscriptstyle{\textrm{0}}}$ and ${\phi}_{\scriptscriptstyle{\textrm{0}}}$. This also implies that ${\phi}_{\scriptscriptstyle{\textrm{0}}}$ can only take values within a finite interval, ranging from $-(\sqrt{2{\rho}_{\scriptscriptstyle{\textrm{0}}}}/m)$ to $(\sqrt{2{\rho}_{\scriptscriptstyle{\textrm{0}}}}/m)$. In a flat universe, the value of the scale factor has no physical meaning. The number of e-folds of inflation depends on ${\phi}_{\scriptscriptstyle{\textrm{0}}}$ but not on ${a}_{\scriptscriptstyle{\textrm{0}}}$: $N=N({\phi}_{\scriptscriptstyle{\textrm{0}}};m,{\rho}_{\scriptscriptstyle{\textrm{0}}})$. So the fraction of trajectories that achieve a phase of inflation lasting more than $N$ e-folds can be written as $\mu=(m\Delta{\phi}_{\scriptscriptstyle{\textrm{0}}})/(2\sqrt{2{\rho}_{\scriptscriptstyle{\textrm{0}}}})$, where $\Delta{\phi}_{\scriptscriptstyle{\textrm{0}}}$ is the range of initial values of the scalar field that yields the required inflationary phase. It is then necessary to evaluate $\mu$ as a function of $N$. There are two cases in which this can be done analytically: (i) at low energy, ${\rho}_{\scriptscriptstyle{\textrm{0}}}\ll m^{2}$, and (ii) at high energy ${\rho}_{\scriptscriptstyle{\textrm{0}}}\gg m^{2}$. At low energy, the calculation of Gibbons and Turok of the probability for inflation can be used to show that Gibbons and Turok (2008)

[TABLE]

where $\mathcal{C}$ is a numerical factor that does not depend on $m$ or ${\rho}_{\scriptscriptstyle{\textrm{0}}}$. For $N\simeq 60$ e-folds, as required to explain the CMB temperature anisotropy, this leads to $\mu(N)\ll 1$. At high energy, one reaches the opposite conclusion. In this case, one can compute $\Delta{\phi}_{\scriptscriptstyle{\textrm{0}}}$ as follows. For a massive quadratic potential the total number of e-folds of inflation can be expressed in terms of the amplitude of the scalar field at the start of the inflationary phase, ${\phi}_{\scriptscriptstyle{\textrm{I}}}$, as $N\approx({\phi}_{\scriptscriptstyle{\textrm{I}}}^{2}/4)$. In turn, ${\phi}_{\scriptscriptstyle{\textrm{I}}}$ can be expressed in terms of the initial value of the scalar field as

[TABLE]

with $z\equiv 8\Gamma^{2}\exp(\sqrt{6}{\phi}_{\scriptscriptstyle{\textrm{0}}})$ and $\Gamma\equiv\sqrt{3{\rho}_{\scriptscriptstyle{\textrm{0}}}}/m$. This formula for the amplitude of the scalar field at the start of inflation is valid in LQC, with the modified Friedmann equation given by Eq. (1). For the standard flat FLRW dynamics, without LQC modifications, the analytical calculations suggest that at the start of inflation the scalar field reaches a maximum value given by (4) minus $(\ln 2/\Gamma)$. In both cases, we find that the range of values of ${\phi}_{\scriptscriptstyle{\textrm{0}}}$ that do not yield an inflationary phase longer than $N$ e-folds is an interval of size $4\sqrt{N}$ centered on ${\phi}_{\scriptscriptstyle{\textrm{0}}}=0$. Hence,

[TABLE]

and $\mu(60)\simeq 0.99999$, which means that all but a tiny fraction of the possible trajectories do not go through a long inflationary phase. It might be tempting to interpret $\mu$ as a probability measure. This is however not that simple. The phase space of the flat FLRW universe presents a serious ambiguity: the Liouville measure is proportional to the scale factor and the scale factor can be rescaled arbitrarily. In addition, as explained just before, $\mu$ depends on the choice of the surface of initial data. More importantly, the fundamental question to ask is: is there a variable on which a flat (or at least known) probability distribution function (PDF) can be assigned? Assuming implicitly that the initial values of the field should have a flat PDF is arbitrary.

In Ashtekar and Sloan (2011) it was argued that the two first issues mentioned above can be solved in LQC. It was indeed claimed that the scale factor can be rigorously factored out of the Liouville measure, and that the bounce provides a preferred choice for the surface of initial data. In this study, following Linsefors and Barrau (2013, 2015) we choose a different perspective. We decide, the other way round, to set initial conditions in the remote past of the contracting branch, when the Universe is classical and well understood (${\rho}_{\scriptscriptstyle{\textrm{0}}}\ll{\rho}_{\scriptscriptstyle{\textrm{B}}}$). This is not only technically justified but also conceptually necessary if the bounce has to be taken seriously in a causal way. Still, we naturally choose a time which is close enough to the bounce so that it is reasonable to assume a scalar field as the main component of the Universe. The phase of the oscillations of the scalar field in the contracting branch is an obvious variable to which a flat PDF can be assigned Linsefors and Barrau (2013). In addition, the key point is that this PDF is preserved over time (as long as one remains in the classical phase when the field oscillates). The numerical analysis of Linsefors and Barrau (2013) shows that at fixed ${\rho}_{\scriptscriptstyle{\textrm{0}}}$, nearly all possible initial values for the scalar field, ${\phi}_{\scriptscriptstyle{\textrm{0}}}$, yield an inflationary phase whose number of e-folds is peaked around $N=142$ e-folds (with ${\rho}_{\scriptscriptstyle{\textrm{B}}}=0.41{m}_{\scriptscriptstyle{\textrm{Pl}}}^{4}$).

The procedure to derive this result is simple:

• •

Consider an initial energy density ${\rho}_{\scriptscriptstyle{\textrm{0}}}={\rho}_{\scriptscriptstyle{\textrm{Pl}}}/\alpha^{2}$, with $\alpha$ large large enough so that the evolution starts in the remote past of the contracting phase.

• •

Choose an initial value for the scalar field and its time derivative by a random sampling of the phase ${\theta}_{\scriptscriptstyle{\textrm{0}}}$ between [math] and $2\pi$, where ${\theta}_{\scriptscriptstyle{\textrm{0}}}$ is defined such as ${\phi}_{\scriptscriptstyle{\textrm{0}}}=\sqrt{\frac{2}{3}}\frac{\Gamma}{\alpha}\sin{\theta}_{\scriptscriptstyle{\textrm{0}}}$.

• •

Solve the dynamics, across the bounce, until the end of slow-roll inflation in the expanding branch.

• •

For each ${\theta}_{\scriptscriptstyle{\textrm{0}}}$, collect the corresponding number of e-folds.

Finally, one can produce the associated histogram which, in a probabilistic interpretation, is the PDF for the number of e-folds. This is illustrated on the right panel of Fig. 5 where we also present the PDFs for several initial energy densities corresponding to different values of $\alpha\equiv\sqrt{{\rho}_{\scriptscriptstyle{\textrm{Pl}}}/{\rho}_{\scriptscriptstyle{\textrm{0}}}}$ in order to show that for large values of $\alpha$ the PDF becomes independent of the initial energy density, as explained analytically in Linsefors and Barrau (2013). Interestingly, the peakedness of the PDF can be understood as follows. The calculation Gibbons and Turok is often considered controversial in standard cosmology because they somehow set “initial conditions” for the final state. However, in the case of a bouncing Universe it implies that almost none of all the possible trajectories, starting at low energy in the contracting branch, have a significant phase of pre-bounce exponential contraction, that is of so-called deflation. A trajectory with deflation in the contracting phase leading to $({\phi}_{\scriptscriptstyle{\textrm{B}}},{\dot{\phi}}_{\scriptscriptstyle{\textrm{B}}})$ can be identified with a trajectory with inflation in the expanding phase with $({\phi}_{\scriptscriptstyle{\textrm{B}}},-{\dot{\phi}}_{\scriptscriptstyle{\textrm{B}}})$. Equation (4) can be used to calculate the value of the scalar field at the bounce corresponding to the trajectory with no deflation. One simply has to solve Eq. (4) with respect to ${\phi}_{\scriptscriptstyle{\textrm{B}}}$ for ${\phi}_{\scriptscriptstyle{\textrm{I}}}=0$ and ${\dot{\phi}}_{\scriptscriptstyle{\textrm{B}}}<0$. In the limit of large $\Gamma$, the solution is well approximated by

[TABLE]

This can then be inserted back into Eq. (4), with ${\dot{\phi}}_{\scriptscriptstyle{\textrm{B}}}>0$, in order to obtain the value of the field at the start of inflation in the expanding phase and the corresponding number of e-folds of inflation. With the standard values for $m$ and ${\rho}_{\scriptscriptstyle{\textrm{B}}}$, this calculation yields $N=142$, in excellent agreement with the numerics (Fig. 3). Moreover, a closer look at Gibbons and Turok’s PDF for the number of e-folds suggests that most trajectories starting in the remote past have less than one e-folds of deflation, see Fig. 5. This means that nearly all trajectories end up with a value of ${\phi}_{\scriptscriptstyle{\textrm{B}}}$ that belongs to an interval of size $\Delta{\phi}_{\scriptscriptstyle{\textrm{B}}}\approx 4$ centered around ${\phi}_{\scriptscriptstyle{\textrm{B}}}^{\scriptscriptstyle{\mathrm{GT}}}$. In terms of number of e-folds this translates into $\Delta N\approx 4\sqrt{N}$, also in agreement with the numerics as can be seen on Fig. 3.

We shall now investigate to which extent the specific modified dynamics is responsible for the peakedness of the probability density function of the number of e-folds in loop quantum cosmology. The argument we have developed in the previous section did not refer to the modified LQC dynamics. It was essentially based on Gibbons and Turok’s analysis combined with the presence of the bounce at Planckian energy density. It can therefore already be guessed that the peakedness does not depend strongly on the LQC modification to the Friedmann equation. To address this question in more details, we consider an artificial bouncing cosmological scenario where the Friedmann equation is left unchanged even at Planckian energy. In this “GR-like” cosmological scenario, initial conditions for a given trajectory are set in the remote past of the contracting branch at the same energy density and with the same values of ${\phi}_{\scriptscriptstyle{\textrm{0}}}$ and ${\dot{\phi}}_{\scriptscriptstyle{\textrm{0}}}$ than for a trajectory which follows the LQC dynamics (as previously considered). The dynamics is divided into two parts: the contracting branch with a negative Hubble parameter and the expanding branch with a positive Hubble parameter. The evolution, starting in the contracting branch, is artificially stopped when the energy density reaches the LQC critical energy density ${\rho}_{\scriptscriptstyle{\textrm{B}}}$. The values of ${\phi}_{\scriptscriptstyle{\textrm{B}}}$ and ${\dot{\phi}}_{\scriptscriptstyle{\textrm{B}}}$ are collected and used as initial conditions for the dynamics in the expanding branch where the initial Hubble parameter is now positive. At the junction between both phases, the Hubble parameter and therefore $\ddot{\phi}$ are discontinuous but $a$, $\phi$ and $\dot{\phi}$ are continuous, as illustrated in Fig. 1, 2 and 6.

The numerical result for the PDF of the number of e-folds in the GR-like scenario is plotted against the LQC prediction on the left panel of Fig. 3. The PDF has the same width and shape than in LQC. This confirms that the peakedness does not depend strongly on the specific LQC modified dynamics, and suggests that this feature would remain in case one incorporates additional quantum gravity corrections to the LQC effective equations. Nevertheless, the number of e-folds corresponding to the peak of the PDF is slightly different in the GR-like scenario than in LQC. This can be explained as follows.

First, it should be noticed that the difference between the GR-like scenario and LQC becomes significant when $\rho\approx{\rho}_{\scriptscriptstyle{\textrm{Pl}}}$. Second, as shown in the previous sections, the fraction of trajectories that have a significant phase of deflation in the contracting branch is tiny. This means that at high energy density, the dynamics of most trajectories is largely kinetic energy dominated. In simplistic terms, deflation can not bring the amplitude of the scalar field to large values because it stops nearly immediately. To inverstigate the difference between the GR-like scenario and LQC, the equation of motion of the scalar field at high energy, and for kinetic energy domination, need to be studied. It is natural to introduce $x\equiv\phi/\sqrt{2{\rho}_{\scriptscriptstyle{\textrm{B}}}}$ and $y\equiv\dot{\phi}/\sqrt{2{\rho}_{\scriptscriptstyle{\textrm{B}}}}$, the so-called potential and kinetic energy parameters. As the duration of inflation depends on the scalar field amplitude, it is sufficient to focus on the potential energy parameter. It is easy to show that in the interesting regime

[TABLE]

In LQC the scalar field is boosted by a short phase of super-inflation, $\dot{H}>0$, during which its amplitude accumulates a surplus of $(\ln 2/\Gamma)$ compared to the standard FLRW dynamics. This yields a difference of $\sqrt{8/3}\ln 2\sqrt{{N}_{\scriptscriptstyle{\textrm{GR}}}}$ between the number of e-folds in both scenarios. With the standard numerical values of $m$ and ${\rho}_{\scriptscriptstyle{\textrm{B}}}$, one gets ${N}_{\scriptscriptstyle{\textrm{LQC}}}-{N}_{\scriptscriptstyle{\textrm{GR}}}\simeq 13$, in agreement with the numerical results.

One can also study the differences between the primordial power spectra of cosmological perturbations in the GR-like scenario and in LQC. As a toy model to focus on the difference between both background dynamics, we set initial conditions for perturbations at an energy density corresponding to the bounce energy density, choosing the Bunch-Davies vacuum as the initial state. The resulting power spectra are shown on Fig. 4 and compared with the usual slow-roll inflation expectation (dotted lines). Such spectra (and their variants including ore subtle LQC effects) are the main observables associated with loop quantum cosmology.

The length of inflation is crucial because it determines the location of the window of wavenumbers relevant for the cosmic microwave background anisotropy measurements, with respect to the characteristic LQC scale ${k}_{\scriptscriptstyle{\textrm{LQC}}}\equiv{a}_{\scriptscriptstyle{\textrm{B}}}\sqrt{{\rho}_{\scriptscriptstyle{\textrm{B}}}}$. On infrared scales, the mode functions remain in the Bunch-Davies state with $\mathcal{P}(k)\propto k^{2}$. We see that in both LQC and the GR-like scenario the power spectra agree with the slow-roll expectations in the ultraviolet regime. Oscillations are present in both scenarios in the range $10^{-3}<k/{k}_{\scriptscriptstyle{\textrm{LQC}}}<1$. The amplitude is larger in the GR-like scenario than in LQC, however the period of the oscillations does not seem to be affected by the specific modified LQC dynamics. This shows that oscillations in themselves are a bounce feature but not a specific LQC feature. This motivates the search for complementary probes such as primordial non-gaussianity Agullo (2015). For a more detailed comparison of the different kinds of power spectra expected in LQC under different assumptions for the mode propagation, see Bolliet et al. (2015).

The most reliable result of loop quantum cosmology is the modified Friedmann equation describing the background dynamics. It receives a quadratic correction in density which prevents the Universe from collapsing into a singularity. In this article, we have investigated the influence of this modified dynamics on the duration of inflation. The conclusion is that the modification of the Friedmann equation has a very small (but non vanishing) impact on the duration of inflation. The key role played by LQC in “predicting” inflation – or more precisely the duration of inflation – is not due to the modified dynamics in itself. It is grounded in two different aspects. First, LQC sets the energy scale. This is the fundamental point. As far as the Universe is assumed to be filled by a massive scalar field, inflation happens naturally if the energy scale “before” inflation is high enough, $\rho\gg m^{2}$. But whereas starting at the Planck energy density in GR is somehow arbitrary, in LQC the bounce energy density can be calculated (modulo some hypotheses) and derived from the full theory, providing a natural energy scale. This is the first important aspect. Second, LQC selects favored conditions at the bounce, see formula (6), corresponding to a favored duration of inflation $N\simeq 145$, for ${\rho}_{\scriptscriptstyle{\textrm{B}}}=0.41{m}_{\scriptscriptstyle{\textrm{Pl}}}^{4}$. This is an interesting prediction rooted in the existence of a pre-bounce phase where a natural variable to which a known PDF can be assigned was identified. This cannot be produced in standard cosmology and is specific to bouncing models.

I Acknowledgement

The work of BB was supported by a grant from ENS Lyon and a Fulbright Grant from the Franco-American commission.