Quantum-Gravitational Effects on Primordial Power Spectra in Slow-Roll Inflationary Models

We review the computation of the power spectra of inflationary gauge-invariant perturbations in the context of canonical quantum gravity for generic slow-roll models. A semiclassical approximation, based on an expansion in inverse powers of the Planck mass, is applied to the complete Wheeler–DeWitt equation describing a perturbed inflationary universe. This expansion leads to a hierarchy of equations at consecutive orders of the approximation and allows us to write down a corrected Schrödinger equation that encodes information about quantum-gravitational effects. The analytical dependence of the correction to the power spectrum on the wavenumber is obtained. Nonetheless, some numerical work is needed in order to obtain its precise value. Finally, it is shown that the correction turns out to be positive, which leads to an enhancement of the power spectrum especially prominent for large scales. We will also discuss whether this correction leads to a measurable effect in the cosmic microwave background anisotropies.


Introduction
The lack of experimental guidance is still one of the main obstacles to advance in the construction and verification of the ultimate theory of quantum gravity [1].There are several indications that the initial phases of the universe might be a good laboratory for such a purpose.In particular, during the inflationary phase of the primordial universe energies close to the Planck scale were present and the conditions for structure formation were set, which makes it plausible that Planck-scale physics might be directly related to the anisotropies of the cosmic microwave background (CMB) and, thus, be observable by our telescopes [2].
In this conference proceedings contribution based on our previous work published in [3,4], we review a computation of the quantum-gravity effects in the power spectra of both gauge-invariant scalar and tensor perturbations.For such a purpose, the perturbations of a homogeneous background universe with an inflaton field will be considered.The gauge-invariant master variables, which encode the complete physical information of the scalar and tensor perturbations, will be constructed in order to deal with the perturbative gauge invariance.The quantization of the whole system (background plus perturbations) will then be carried out following the canonical approach, which will lead to a master Wheeler-DeWitt equation.In order to solve it, a Born-Oppenheimer type of approximation will be applied by performing a power expansion in the inverse of the Planck mass.In this way, the search for a solution of the Wheeler-DeWitt equation will be translated into an infinite set of problems, one at every order of the expansion parameter.Finally, in order to keep the potential of the inflaton field as generic as possible, all the computations will be performed assuming a slow-roll regime.

Hamiltonian Framework for Classical Perturbations of an Inflationary Universe
Perturbation theory is a very useful tool in many areas of physics and, in particular, in gravitational physics.The idea is fairly simple; one typically considers a known solution g µν of the Einstein equations as the background manifold.The perturbations of the metric h µν are then understood as fields propagating on that background, such that the dynamics of h µν is ruled by the linearized version of the Einstein equations.
In particular, the action of general relativity with a scalar matter field φ is given by, π ij and π φ being the conjugate momentum of the spatial metric g ij and of the matter field, respectively.
The Hamiltonian H is a linear combination of constraints, with N being the lapse function and N k the shift vector.While the first-order variation of the action (δS = 0) provides the background equations of motion, the second variation δ 2 S provides an effective action functional for the linear perturbations: In our approach, we intend to quantize all degrees of freedom, thus a complete action for the whole system is needed.This is obtained by considering the sum of both actions, which implies that the total Hamiltonian has the following form, The equations of motion can be obtained by computing the Poisson brackets of different variables with this Hamiltonian or, equivalently, by taking variations of the full action.In particular, it is straightforward to see that this system has eight constraints: The first line corresponds to the background constraints, corrected by terms quadratic in the perturbations, whereas the second line represents expressions linear in the perturbations and it corresponds to the perturbative gauge freedom.At this point, one could follow several routes in order to perform the quantization of this system.On the one hand, one could quantize the full system à la Dirac, without fixing any gauge classically.This approach was followed by Halliwell and Hawking (for a homogeneous and isotropic background) in a well-known paper [5].On the other hand, one could either totally or partially fix the gauge classically and then proceed to the quantization.The problem of this approach is that when extracting physical results, it is difficult to know whether those results are gauge independent and thus physically meaningful.
In order to deal with the perturbative gauge freedom, there is an intermediate route which can be followed, which consists in constructing master gauge-invariant variables.Examples of such variables are Regge-Wheeler [6] and Zerilli [7] variables for the Schwarzschild background; and Mukhanov-Sasaki [8,9] variables for the Friedmann-Robertson-Walker universe.Technically, these variables are constructed by performing an adequate canonical transformation, with the requirement that in the new variables the perturbative constraints are trivialized.The great advantage of this approach is that one does not need to consider the quantization of the perturbative constraints (since they are trivial) and all the results obtained are automatically gauge invariant.
Let us now specialize all of the above theory to the Friedmann-Robertson-Walker universe.By using the conformal time η, the background metric is given as follows, with a being the scale factor.Note that, in this homogeneous case, the shift vector is chosen to be vanishing, N i = 0, and thus the only constraint left from (4) will be the Hamiltonian constraint H + 1 2 δ 2 H = 0.As it is usually done, the perturbation tensor h µν is then decomposed into scalar, vector, and tensor components.The construction of master gauge-invariant quantities leads to three variables that contain the three physical degrees of freedom of the system (two corresponding to the gravitational wave and one to the scalar matter field): the Mukhanov-Sasaki variable v encodes the physical information of the scalar sector, whereas the tensorial sector is by construction gauge invariant and the two polarizations of the gravitational wave will be named as v (λ) = ah (λ) , with λ = +, ×.The Hamiltonian constraint is then decomposed as follows, Here, H 0 is the Hamiltonian of the background minisuperspace model, where π α is the conjugate momentum of α := ln a, V (φ) is the potential of the inflaton field, and we have defined the reduced Planck mass m 2 P := 3 4πG .Furthermore, H S and H T are the Hamiltonians corresponding to the scalar and tensorial sector, respectively.These two are Hamiltonians of parametric harmonic oscillators.In particular, the scalar one can be written as follows, where π k is the conjugate momentum of the Mukhanov-Sasaki variable v k , and the frequency is defined as ω 2 S (η) = k 2 − z z , z being a given function of background objects.The Hamiltonian of the tensor part has exactly the same form, the only difference being that the frequency is given by

Quantization and Semiclassical Approximation
The quantization of the Hamiltonian constraint (7) gives rise to a master Wheeler-DeWitt equation.By imposing an adequate product ansatz for the full wave function, and assuming that different k-modes are independent (which translates to assuming that the action of the Hamiltonian for every k defined below in ( 12) is negligible when acting on another k ), one obtains the following equation for every mode k as well as for both the scalar and tensor sectors: Note that in this equation the minisuperspace variables and the part involving inhomogeneous degrees of freedom for a given k appear on the same footing, and the fiducial volume of the minisuperspace has been absorbed in the minisuperspace variables.The frequency ω should be taken as ω S or ω T , depending on whether we are dealing with the scalar or tensorial sector respectively.To simplify the presentation from here on, if not otherwise stated, we will only consider the scalar sector, which is more involved, and thus ω should be understood as ω S ; in the last Section 6, we will also summarize the results corresponding to the tensorial part.
Following the approach presented in [10,11], and also used in [12,13], the next step of the analysis involves a power expansion in the inverse of the squared Planck mass m 2 P .At leading order, one recovers the Hamilton-Jacobi equation, which encodes the classical dynamics of the system.At the next order, the approximation of quantum field theory on a curved (classical) spacetime is obtained via the following Schrödinger equation: where the matter Hamiltonian operator is, as usual, the one corresponding to a parametric oscillator, At the following order, new physics is obtained and the equation to be solved is a modified Schrödinger equation: where an auxiliary potential has been defined as Note that the quantum-gravitational effects are encoded in the last term of Equation (13).In fact, this term is just a product between ψ (1) k and another term that depends on background quantities and the wave function ψ (0) k .This latter is in general a complex function and thus this correction term is a possible source of unitarity violation.In order to avoid such an undesirable effect and to construct a consistent framework, we will just consider the real part of this term.The appearance of these unitarity-violating terms has been discussed in detail in our previous work [3,4] and will be subject to further investigations.

Gaussian Ansatz and Power Spectra
In summary, we need to solve Equations ( 11) and ( 13) for wave functions ψ k .Since the matter Hamiltonian (12) corresponds to a parametric oscillator, it is thus natural to consider a Gaussian ansatz as follows, for both ψ k .Here, and in the following, the superscript (0) stands for the uncorrected and (1) for the corrected case.Assuming this form for the wave functions, the uncorrected and corrected power spectra are directly given by the real part of the corresponding Gaussian width, Ω (0) k and Ω (1) k , respectively.More precisely, for the scalar perturbations the corresponding power spectrum is given by P (1) where is the slow-roll parameter, which will be defined in (21) below.Furthermore, P S (k) is the usual spectrum that is obtained by quantum field theory methods [14].In this setting, this would be obtained by replacing Ω (1) k by Ω (0) k in this latter equation.Therefore, the relative difference between the two (corrected and uncorrected) spectra is encoded in the object evaluated in the large-scale limit.
Since the quantities of interest are the Gaussian widths, let us present the equations obeyed by them in both the uncorrected and the corrected cases: i i Ω (1) where the prime denotes the derivative with respect to the conformal time.As can be seen, the equations are exactly the same, with the only exception that for the former case the frequencies are the ones defined after Equation ( 9), ω S or ω T , depending on whether one is considering the scalar or tensorial sector, whereas for the latter case the corrected frequencies ω are defined as follows: In summary, in order to get the quantum-gravity correction to the power spectrum, one needs to solve Equations ( 18) and ( 19) and compute then ∆ S as defined in (17).Let us briefly comment about the initial conditions to be considered.The usual initial conditions that are considered in the quantum-field-theoretical approach are given by the so-called Bunch-Davies vacuum, which implies that at the beginning of inflation, when the mode is well inside the horizon, it behaves as a mode on a Minkowski background, without feeling the curvature.In this setting, this implies to choose Ω (0) k = ω at the onset of inflation.For the corrected case, we mimic this initial data and also request that the mode is oscillating with constant frequency and amplitude, which is translated to the initial condition Ω (1) k = ω.The difference with the uncorrected case is that the frequency of the mode, even if constant at the beginning, gets modified.

Slow-Roll Approximation
The numerical value of ∆ S will depend on the inflaton potential V (φ).Instead of fixing this potential, we will consider the slow-roll regime, which is quite generic and includes several types of potentials.Therefore, we define the two slow-roll parameters, and their combination γ := 2 − δ.In these definitions, H is the Hubble rate and the dot stands for the derivative with respect to cosmic time.
By performing an expansion in these parameters, and dropping quadratic and higher-order terms, one can obtain the following form for the scalar and tensorial frequencies [14]: Note that γ = converts the scalar equation into the tensorial one.We will use this fact in the following section to recover the tensorial results from the ones of the scalar sector.
In addition, at this level of approximation, the auxiliary potential V (14) can be written in the following way [4]: where H k is the value of the Hubble factor at the horizon exit, which depends on the mode k.
With these expressions at hand, one can get the form of the corrected frequencies (20) in the slow-roll approximation and thus write our main Equations ( 18) and ( 19) in this regime without any reference to the potential of the inflaton field.

Results
The main Equation ( 19) is quite complicated and one needs to resort to numerical methods in order to solve it in the slow-roll regime.Nevertheless, it is possible to obtain the analytical dependence of the corrections to the power spectrum with numerical coefficients.More precisely, for the scalar sector it can be written as follows, We had to introduce a reference scale k at this point, for which the most reasonable choice in order to give an estimate of the corrections seems to be the pivot scale used in the Planck data analysis.The issue in choosing the most appropriate reference scale is, however, still under debate and was also discussed in [15].
The most important feature of this correction is that it is not scale independent, since it goes as k −3 , and it is thus more relevant for large scales.This is something expected since modes with small wave number are the first to exit the Hubble horizon and therefore experience stronger quantum-gravity effects.Furthermore, since the absolute value of and δ is smaller than 10 −2 , the correction term (24) leads to an enhancement of the power spectrum.One can also analytically compute the corrections for the spectral index n S and its running.In particular, the spectral index is corrected by a negative term, whereas its running gets a positive correction [4].All these corrections are of the order of ∆ S .
Concerning the corrections to the tensorial power spectrum, as commented in the previous section, these can be obtained from the result of the scalar part, just by imposing γ = .In this way, one obtains The form of this correction is qualitatively the same as the scalar one.Therefore, the main features and conclusions drawn above for the scalar sector are also valid for the tensorial case.We can also compute the correction to the tensor-to-scalar ratio r, which leads to r (1) := P (1) Regarding the possible observability of these effects, taking into account data of the Planck mission [16], it is possible to give an estimate of the correction terms ∆ S and ∆ T .Given that the tensor-to-scalar ratio deduced from the Planck data is less than r 0.11, the ratio H 2 inf /m 2 P , where H inf is the Hubble rate during inflation, is bounded to be less than This leads to the following upper bound for our quantum-gravitational correction terms where the ratio between ∆ S and ∆ T is around ∆ S /∆ T ≈ 1.02.For the bound on the correction ∆r := r (1) − r (0) of the tensor-to-scalar ratio, we obtain ∆r In order to see whether these bounds lead to corrections that are too small to be observed with the current experimental data, we calculate the quantities C that are used to describe the CMB anisotropies.They are defined as where we introduced the conformal time at horizon crossing η hor and recombination η rec and used i = 0, 1 to indicate the uncorrected or corrected quantities, respectively.We can, in fact, calculate the correction to the C expressed as ∆C := C This allows us to compare this correction to cosmic variance, which fundamentally restricts the accuracy with which we can measure the CMB anisotropies.The statistical uncertainty arising from cosmic variance can be expressed by such that for = 2, we have ∆C  [4], the order of magnitude that the ratio H 2 k /m 2 P would need to have, such that the quantum-gravitational corrections obtained here overcome cosmic variance and could thus be measurable, can be estimated to be which is way above the bound obtained in (27).Hence, we need to conclude that these corrections are unfortunately too small to be observed with the current CMB data.Nevertheless, our calculations show that concrete predictions can be obtained from approaches to quantum gravity.The question as to whether there are other situations where such effects can be observed, remains presently open.