arxiv: v1 [gr-qc] 24 Jan 2019

Transcription

1 A moment approach to compute quantum-gravity effects in the primorial universe Davi Brizuela 1 an Unai Muniain Fisika Teorikoa eta Zientziaren Historia Saila, UPV/EHU, 644 P.K., Bilbao, Spain arxiv: v1 [gr-qc] 4 Jan 019 Abstract. An approach to compute quantum-gravity corrections to the scalar an tensorial power spectra of the inflationary perturbations is presente. The analysis of the Wheeler-DeWitt equation is performe by a ecomposition of the wave function into its infinite set of moments, which must obey certain system of (first-class) constraints. Consiering a semiclassical approximation, the system is truncate at secon orer in moments an an appropriate gauge-fixing conition is introuce, which allows us to interpret the scale factor of the universe as an internal time. The evolution of the ifferent fluctuations an correlations is then explicitly consiere for a e Sitter universe. An approximate analytical solution is obtaine for the corrections of the power spectra, which prouces an enhancement of power for large scales. Remarkably, the result is in agreement with previous stuies in the literature that mae use of very ifferent semiclassical approximations. Finally, the numerical implementation of the system is also consiere to verify the valiity of the analytical solution. 1 Introuction Even if uring the last ecaes a large amount of approaches have been attempte to construct a consistent theory that escribes the quantum behavior of the gravitational interaction [1], the ultimate theory is yet unknown. In particular, the lack of experimental ata is one of the issues that makes progress in this research area complicate. Since quantum gravity effects are suppresse by the value of the Planck mass, they are inee very ifficult to be measure. Nonetheless, these effects coul have consierable consequences on some astrophysical scenarios, such as the primorial inflationary universe. Although the raiation from such early times is not irectly observable, quantum fluctuations present at those eras result into anisotropies on the Cosmic Microwave Backgroun (CMB), for which quite accurate experimental ata is alreay available []. In this context, it is very usual to apply the approximation of quantum fiel theory (QFT) on fixe classical backgrouns to stuy the evolution of such fluctuations. In this approximation the matter an geometric perturbations are treate as quantum entities, whereas the backgroun universe is consiere to follow a completely classical behavior. In fact, very goo results are obtaine by following this metho (see, e.g., the reviews [3 5] or the textbook [6]), which agree with the commente ata. Nevertheless it is expecte that, since these quantum fluctuations were originate at the very beginning of the universe, the quantum behavior of the backgroun might inee prouce some nontrivial (quantum-gravity) corrections to the stanar results of QFT on classical backgrouns. In the present work we will analyze corrections to the power spectra of the gauge-invariant scalar an tensorial inflationary perturbations base on the canonical approach to quantum gravity that leas to the Wheeler-DeWitt equation. Even if this theory might not be the funamental theory of quantum gravity, it is a very conservative approach an it is expecte to be correct at least at a semiclassical level. There are several works in the literature that have alreay analyze these corrections, by using ifferent Born-Oppenheimer-like approximations. On the one han, in [7 10] an expansion with respect to the Newton s gravitational constant G N was consiere. In particular, the term of orer O(G 0 N ) provies exactly the QFT approximation an quantum-gravity corrections appear at next orer O(G N ). On the other han, in [11, 1] a ecomposition closer to the stanar Born- Oppenheimer approximation use in molecular physics was followe, by separating the wave function of the fast egrees of freeom, which correspon to the perturbations, an of the slow ones, given by the backgroun universe. In this case, the fast Hamiltonian reprouces the result obtaine within the QFT approximation, while quantum corrections are encoe in the slow sector. Both approaches have agree qualitatively in an enhancement of power for large scales. 1 aress: avi.brizuela@ehu.eus aress: unaimuni@gmail.com 1

2 The approach that will be use here is very ifferent from these previous ones. We will use a formalism base on quantum moments for constraine systems, as introuce in [13 15], which oes not require a classical eparametrization. More precisely, the wave function will be first ecompose into its infinite moments, such as its expectation values, fluctuations an correlations. The Wheeler-DeWitt equation will then be rewritten as an infinite set of (first-class) constraint equations for these moments. Since a moment of orer n scales as n/, for semiclassical peake states one can safely truncate the system at a given finite orer an get reliable results. In our case, in orer to obtain the main quantum-gravity corrections to the inflationary power spectra, it will be enough to truncate the equations at secon orer in moments, or orer. After the truncation, a gauge conition will be impose to solve the constraint system, that will allow us to interpret the scale factor as an internal time. The obtaine equations of motion will be solve analytically, making use of certain approximation. In aition, the numerical implementation of the system will also be consiere, which will valiate the approximate analytical solution. The well-known results of the QFT on classical backgroun approximation will be easily recovere just by ropping the quantum moments corresponing to the backgroun egrees of freeom. Finally, let us briefly mention that in the context of loop quantum cosmology, which is also formulate by following a canonical quantization of the system, there has been an intense activity uring the last years to obtain similar quantum-gravity corrections for the inflationary power spectra [16 4]. Even if here we will not consier this theory, we point out that the metho use an evelope on this paper coul also be applie to that framework. The rest of the paper is organize as follows. In Section we will briefly review the classical canonical escription of an inflationary universe, which will be use to carry out the quantization in Section 3. In this section two ifferent quantizations will be presente: the one corresponing to the approximation of QFT of classical backgroun, an the full quantization of the system, that will lea to the Wheeler-DeWitt equation. In Section 4 we will apply the formalism of quantum moments for both quantization schemes an the expression of the power spectra will be obtaine. In Section 5 we will consier the particular case of a e Sitter universe. The usual result of QFT for the power spectrum will be erive an corrections for that formula will be obtaine by solving analytically an numerically the equations for the moments corresponing to the Wheeler-DeWitt quantization. Our main results an conclusions will be summarize in Section 6. The canonical formulation of the classical moel In this section we briefly summarize the stanar canonical formulation of linear perturbations of a homogeneous an isotropic universe..1 Backgroun homogeneous an isotropic universe Inflationary ynamics is very well escribe by consiering the perturbations of a scalar matter fiel φ that propagates on a homogeneous an isotropic backgroun. This backgroun is escribe by the Friemann-Lemaitre- Robertson-Walker metric ( ) r s = N(t) t + a(t) 1 kr + r (θ + sin θϕ ). (1) In this expression there are only two time-epenent variables: the scale factor a(t) an the lapse function N(t). Assuming that the matter fiel is minimally couple with gravity an its self-interaction is given by a potential V (φ), the action for the backgroun is given as [5] [ S[a, N, φ] = L 3 t 1 ( ) ] aȧ G N kn + a3 N φ Na 3 V (φ), () where we have efine the reuce gravitational coupling constant as G = 4πG N, (3) 3 G N being Newton s constant, an L 3 stans for the integrate spatial volume. For the spatially compact cases L 3 is the total finite spatial volume of the universe; whereas for the flat case (k = 0), which will be consiere in the

3 rest of the paper, it is in principle infinite. Nonetheless, in orer to regularize the above expression, it is possible to absorb this volume by rescaling the ifferent objects as follows, a a L, t Lt, N N L. (4) After applying this rescaling, an by performing a Legenre transformation, the Hamiltonian of this system is obtaine in a straightforwar way, ( ) H 0 = N Gπ a a + π φ a 3 + a3 V (φ). (5) As it is well known, the lapse function is not a ynamical variable, but a Lagrange multiplier. The variation of the action with respect to it provies the constraint H 0 = 0, (6) which implies that the Hamiltonian of the system must be vanishing. From this point on, we will choose the lapse as N = a, which efines the conformal time t = η.. Cosmological perturbations The perturbe line element has the general form, s = a (η){ (1 A)η + ( i B)x i η + [(1 ψ)δ ij + i j E + h ij ]x i x j }. (7) The perturbative variables can be classifie in three ifferent sectors: the scalar, vectorial an tensorial ones. At linear level ifferent sectors evolve inepenently an they can thus be treate separately. The functions A, B, E an ψ, in combination with the perturbation of the fiel δφ, are scalar perturbations. Vector perturbations have not been consiere when writing the above form of the metric since, in the absence of vectorial matter, they are purely gauge. Furthermore, this type of perturbations are known to rapily ecay uring the accelerate expansion of the universe. Finally, h ij is the tensorial part that escribes the gravitational waves. Tensorial perturbations are invariant uner gauge transformations, but the scalar perturbations are not. Nevertheless, they can be combine in a unique gauge-invariant quantity that encoes the complete physical information of the scalar sector. This quantity is calle the Mukhanov-Sasaki master variable v(η, x), an is efine as [6] v(η, x) := a a [a δφ + a φ (B E ) + aφ (A + B E )], (8) where the prime enotes the erivative with respect to the conformal time η. Since the equations of motion are linear, it is useful to introuce the Fourier transform as 1 v(η, x) = 3 k v (π) 3/ k (η) e ik x. (9) R 3 In orer to obtain the ynamics of the variable v, the action of the system S is expane up to secon orer in the perturbations, δ S. In terms of the Fourier components the scalar part of the action then reas 3 δ S = 1 ( [ ]) (a ɛ) η 3 k v kv k + v k vk a k, (10) ɛ where ɛ is known as the first slow-roll parameter an is efine as ɛ := H ah, (11) 3 Since the time variable η has been rescale above (4) to absorb the volume of the universe L 3, one also nees to rescale the perturbative variable v k (for both scalar an tensorial sectors) an the wave number k as v k L v k an k k/l so that L oes not appear in these expressions. From this point on, all the variables are rescale as commente. In this way, η an k are imensionless, whereas a has imensions of length an vk has imensions of action. For a more explicit iscussion about this rescaling transformation see, e.g., Section III of [7]. 3

4 H := a /a being the Hubble parameter. From this action, it is straightforwar to obtain the canonical momentum for each k, π k = L v k = v k. (1) Performing then a Legenre transformation the Hamiltonian can be finally euce: H S = 1 ( 3 k π k πk + v k vk [k (a ɛ) ]) a. (13) ɛ This Hamiltonian correspons to a system of ecouple harmonic oscillators, with a time-epenent frequency given by ω S(η, k) := k (a ɛ) a ɛ. (14) Regaring tensorial perturbations, as commente above, they are escribe by the tensor h ij. Besies being symmetric, this tensor is also traceless h i i = 0 an ivergence-free i h ij = 0. Therefore it just contains two physical egrees of freeom, which correspon to the two polarizations of the gravitational wave: h + an h. In particular, if the Cartesian z coorinate is chosen to be the irection of the propagation of the transverse wave, these polarization terms are given by h + = h xx = h yy an h = h xy = h yx. The erivation of the Hamiltonian for the tensorial part is completely analogous to the process explaine above for scalars. In orer to write it own, let us efine the master variable v s k := a 3G hs k, (15) where s stans for the two polarizations, + an. After expaning the action, transforming to Fourier components an performing the Legenre transformation, the tensorial Hamiltonian is obtaine: H T = 3 k ( ) 1 πs kπk s + ω T (η, k) v s kvk s. (16) s Note that for these tensor moes the sum for the two polarizations s must be consiere. The Hamiltonian correspons again to a system of ecouple harmonic oscillators, but in this case the frequency is given by ω T := k a a. (17) Even if both v k an π k are complex, in the Hamiltonian they always appear in quaratic combination with their corresponing complex conjugate. In orer to procee to the quantization, for each of the variables one in principle shoul efine a new set of two real variables, for instance by taking their real an imaginary parts as one, e.g., in [7]. However, since these variables never appear outsie the commente quaratic combination, an the result will thus be unaffecte by this consieration, in orer to enlighten the notation we will just use the variables v k an π k as if they were real. 3 The quantization of the moel In this section two ifferent quantizations of the classical system presente above will be consiere. The first one will be the stanar one given by quantum fiel theory on fixe classical backgrouns. In this case only the perturbative variables v k will be quantize, whereas the backgroun will be left classical. In the secon quantization the full moel (backgroun plus perturbations) will be quantize. 3.1 Quantum fiel theory on classical backgrouns In this approach each moe of the tensorial an scalar perturbations can be consiere inepenently, an one just promotes to an operator their corresponing Hamiltonian, either (13) or (16), in orer to obtain the functional Schröinger equation, ( ˆπ Ĥ k Ψ(η, v k ) := k + ω (η, k)ˆv k ) 4 Ψ(η, v k ) = i Ψ(η, v k). (18) η

5 In this representation the momentum operator acts as a erivative with respect to v k, that is ˆπ k := i / v k, whereas the position operator acts by multiplication. The epenence on the backgroun universe is entirely insie the classical function ω (η, k), which has the form given by (14) or (17) epening whether v k correspons respectively to a scalar or tensorial moe. 3. Full quantization of the moel: the Wheeler-DeWitt master equation In orer to perform a quantization of the full moel, one nees to consier the total action that escribes completely the backgroun an perturbative ynamics. This is given by S total = S + δ S. (19) Starting from this action for a general lapse N(t), the total Hamiltonian is obtaine, which is just the sum of the backgroun an perturbative Hamiltonians: H total = N(H 0 + δ H) = N(H 0 + H S + H T ). (0) In this expression, the lapse N multiplies the whole Hamiltonian an the backgroun constraint H 0 = 0 is not satisfie anymore but it gets quaratic corrections in the perturbative variables. The Hamiltonian constraint now reas H 0 + H S + H T = 0. (1) More explicitly, by choosing the lapse N(t) = a, which correspons to the conformal time η, the classical Hamiltonian constraint takes the following form, H total = G π a + π φ a + a4 V (φ) k ( πk + ω Svk ) k ( (π s k ) + ωt (vk) s ) = 0. () s Once this Hamiltonian is at han, the quantization proceure is performe by promoting ifferent variables into operators that satisfy the canonical commutation relations, [â, ˆπ a ] = i, [ ˆφ, ˆπ φ ] = i, [ˆv k, ˆπ p ] = i δ(k p), (3) where v k stans now for both scalar an tensorial components. The action of these operators on a wave function on the position representation Ψ(a, φ, {v k }) is as follows, âψ = aψ, ˆπ a Ψ = i Ψ a, ˆφΨ = φψ, ˆπφ Ψ = i Ψ φ ˆv k Ψ = v k Ψ, ˆπ k Ψ = i Ψ v k. The master Wheeler-DeWitt equation is then obtaine just by requesting that the action of the Hamiltonian operator on the wave function is vanishing, Ĥ total Ψ = 0. (4) This equation is very ifficult to eal with an, following [10], one can further simplify it by assuming a prouct ansatz for the wave function Ψ an ropping cross terms that involve ifferent moes. In this way, one obtains an iniviual Wheeler-DeWitt equation for each moe, ( ) ĤΨ k := Gˆπ a + ˆπ φ â + â4 V ( ˆφ) + ˆπ k + ω (η, k)ˆv k Ψ k = 0. (5) Physically, consiering this last equation instea of the full one (4) correspons to assuming just the presence of one perturbative moe in the classical moel. Therefore, for a given moe k, we are neglecting possible effects from other moes on its ynamics, but only analyzing how the quantum behavior of the backgroun will affect its evolution. This is, inee, the main goal of the present work. The equation above is vali for both a scalar or a tensorial moe an that is why we have not written any of the subinices S or T in the frequency ω. In fact, from this point on, all expressions will be vali for both scalar an tensor moes, an thus we will use the symbol ω inistinctly. Finally, note that in this quantization proceure the frequency ω has been left as a classical function of time, given by its form (14) or (17). This is a usual approximation in the literature, see e. g. [7,11], since the quantization of ω turns out to be problematic. In particular, besies the quite intricate form of ω in terms of funamental variables, the main ifficulty is that inverse of the momenta appear in its enominator which, upon quantization, woul imply nonlocal operators. 5

6 4 Quantum moments an the power spectra for the perturbations In this section we will analyze the two quantizations performe in the previous section by ecomposing the corresponing wave function into its infinite set of moments. 4.1 Quantum moments for quantum fiel theory on classical backgrouns In orer to analyze the system given in the approximation of QFT on curve backgrouns that is escribe by the Schröinger equation (18), we will efine the following quantum moments, (v n kπ m k ) := (ˆv k v k ) n (ˆπ k π k ) m W eyl, (6) where the Weyl subscript means that the expectation value is evaluate for the totally symmetric orere set of operators, an the expectation values v k := ˆv k an π k := ˆπ k have been efine. This is inee a slight abuse of notation, as we are using the same symbols to enote the classical variables an these expectation values. Nonetheless, we consier it is not worth introucing a more complicate notation since the istinction will be clear from the context an, in particular, from this point on v k an π k will always stan for the expectation values of their corresponing operators, an not for the classical variables. The value n + m will be referre as the orer of the quantum moment. In aition, as we will be ealing with each k moe separately, we will consier that v k an π k are nonfiel quantum variables, satisfying [ˆv k, ˆπ k ] = i, an formally regularize the Dirac eltas as δ(0) = 1. Following this efinition, the wave function is replace by its infinite set of quantum moments. In fact, the knowlege of all moments (incluing the expectation values v k an π k ) is completely equivalent to the full expression of the wave function. The avantage of this formalism is that one works irectly with observable quantities. Furthermore, a moment of orer n has the imensions of n. Therefore, for peake semiclassical states ( 1), one can make the assumption that higher-orer moments are negligible, truncate the infinite set of moments an obtain physically reliable results. The ynamics of these quantum moments will be given by the expectation value of the Hamiltonian Ĥk, which can be written in terms of the above moments by performing an expansion aroun the expectation values v k an π k, Ĥk(ˆv k, ˆπ k ) = H k (v k, π k ) + n+m 1 n+m H k n!m! vk n (vkπ n πm k m ) = 1 ( π k k + ( π k ) + ω vk + ω ( v k ) ). (7) (Here, an in the rest of the paper, in orer to follow the stanar notation, the fluctuations have been enote as ( X) := (X ).) As can be seen, in the general case the expectation value of the Hamiltonian woul be given by an infinite sum in moments. Nonetheless, when the Hamiltonian is harmonic (efine as containing only up to quaratic combinations of basic variables), as it is the case, only up to secon-orer moments appear in the expression of its expectation value. These harmonic Hamiltonians are of a very special kin an their properties in terms of this moment formalism are explaine in etail, for instance, in [8]. One of their properties is that the equations of motion for ifferent orers in moments ecouple. In this paper we will be intereste in computing the power spectra for the perturbations which, as will be explaine below, are given by the fluctuation of the perturbative master variable ( v k ). Therefore, in orer to obtain the exact evolution of that variable given by this harmonic Hamiltonian, it will be enough to consier the system up to secon orer. The equations of motion of all variables are obtaine by computing the Poisson brackets with the above Hamiltonian. These brackets can be efine in terms of the commutators by the stanar relation { ˆX, Ŷ } = i/ [ ˆX, Ŷ ]. In particular, the evolution for the expectation values an secon-orer moments are given by, v k = π k, (8) π k = ω v k, (9) ( ( vk ) ) = (v k π k ), (30) ( (v k π k )) = ( π k ) ω ( v k ), (31) ( ( πk ) ) = ω (v k π k ). (3) For future convenience, we note that the equations of motion for secon-orer moments can be rewritten as a unique equation for the fluctuation ( v k ), ( ( vk ) ) + 4 ω ( ( v k ) ) + 4 ω ω ( v k ) = 0. (33) 6

7 As can be seen, the expectation values evolve exactly as in the classical theory. In orer to solve these equations is usual to fin a region where the frequency ω tens to a constant. In inflationary scenarios this happens at the beginning of inflation (kη ) since the moe is well-insie the horizon an it oes not feel the effects of the curvature. In that limit the frequency tens to the wave number ω k an the moe behaves as a free moe evolving on a Minkowski backgroun. The state of the moe is then assume to be stationary. This stationarity conition can be implemente just by setting equal to zero all time erivatives in the equations above (8) (3); which implies that the expectation values must be vanishing, as well as the correlation, v k = π k = 0, (34) (v k π k ) = 0. (35) In aition, from equation (31), the following relation between fluctuations must hol, ( π k ) = k ( v k ), (36) where the frequency has been replace by its constant value k. Note that, no matter what is the evolution of the frequency, the stationarity conitions for the expectation values (34) will be obeye uring the whole evolution. An, as long as the frequency ω is kept constant, the conitions (35)-(36) on secon-orer moments will also be obeye. But, the evolution of the universe makes the frequency time-epenent an thus, when the moe begins to feel the effects of the curvature, these conitions will cease to be vali an the values of the fluctuations an the correlation will not be constant anymore. Finally the above initial conitions (34) (36) are just 4 relations for our 5-variable system (8) (3). Therefore, one freeom is left to be fixe, which is relate to the energy of the initial state at the beginning of inflation. It is usual, though not necessary an in fact several works can be foun in the literature analyzing an excite initial state, see for instance [9 31], to impose that the state is on its funamental (less-energetic) level. This choice is usually known as the Bunch-Davies vacuum [3]. Replacing conitions (34) (36) in the Hamiltonian (7) an taking into account that the Heisenberg uncertainty relation, ( v k ) ( π k ) ( (v k π k )) 4, (37) must be obeye, it is easy to obtain the initial value for the fluctuations of the basic variables, ( π k ) = k, ( v k) = k. (38) This obviously correspons to the funamental state of the harmonic oscillator an its energy is k/. As a sie remark, let us mention that in this framework this is the only place where the Planck constant enters the computations (it oes not appear in the Hamiltonian (7) or in the equations of motion (8) (3)). Therefore, the appearance of the Planck constant in the result will be just ue to this choice of initial state. Finally, let us note that the ynamics given by the present Hamiltonian (7) conserves the combination that appears on the left-han sie of the uncertainty relation (37), which is another generic property of the harmonic Hamiltonians [8]. This happens because the evolution of the uncertainty relation is given by thir-orer erivatives of the Hamiltonian an, thus, if those are vanishing this combination is conserve. Therefore, any state saturating initially the uncertainty relation will saturate it uring the whole evolution. Even so, in the general case, each moment inepenently will not be constant an thus the state will be eforme through evolution. 4. Quantum moments for the Wheeler-DeWitt equation In the full quantum moel, the backgroun egrees of freeom, escribe by the couples (a, π a ) an (φ, π φ ), are also promote to operators. Therefore, in orer to escribe the complete system, their corresponing moments have also to be taken into account. The notation will be similar to the one use in the previous section. We will only eal with up to secon-orer quantities. Therefore the fluctuation of a given operator ˆX will be enote as ( X), whereas the correlation between two operators ˆX an Ŷ, will be enote as (XY ), where the operators will be symmetrically orere an thus (XY ) = (Y X). 7

8 As oppose to the case of QFT on classical backgrouns of the previous section, in the full quantum picture there is not an evolution equation of the form (18), rather the constraint equation (5) governs the ynamics of the system. Furthermore, in this case the Hamiltonian is not harmonic an ifferent orers in moments couple. In particular, higher-orer moments enter into the equations of motion for the expectation values, an thus these o not follow the classical trajectories anymore. As a first approximation, we will assume a peake semiclassical state an introuce a truncation of the system at secon orer in moments. In orer to rewrite the Wheeler-DeWitt equation (5) in terms of the moments we will follow [33] an request that, since the action of the Hamiltonian on any physical wave function must be vanishing, the expectation value of the Hamiltonian constraint multiplie from the left-han sie by any prouct of basic operators shoul be vanishing as well. More precisely, the following set of infinite relations shoul be obeye (ˆπ a π a ) n1 (â a) n (ˆπ φ π φ ) n3 ( ˆφ φ) n4 (ˆπ k π k ) n5 (ˆv k v k ) n6ĥ = 0. (39) This is in fact a tower of constraints that shoul be solve orer by orer. As commente, we will truncate the system at secon orer. In this way, only 7 constraints will have to be satisfie, which are the following: C := Ĥ = G ( π a + ( π a ) ) + 1 a 4 ( a π φ + a ( π φ ) 4aπ φ (aπ φ ) + 3π φ( a) ) + a V (φ) ( a + 6( a) ) + 8a 3 V (φ) (aφ) + a 4 V (φ)( φ) + π k + ( π k ) + ω (vk + ( v k ) ) = 0, (40) C πa := 1 (ˆπ a π a )Ĥ = Gπ a( π a ) + π φ a (π aπ φ ) + 1 ( 4a 6 a 3 V (φ) πφ ) ( (πa a) i ) + a 4 V (φ) (π a φ) + π k (π a π k ) + ω v k (π a v k ) = 0, (41) C a := 1 (â a)ĥ = G π a( (π a a) + i ) + π φ a (aπ φ) + 1 ( 4a 6 a 3 V (φ) πφ ) ( a) + a 4 V (φ) (aφ) + π k (aπ k ) + ω v k (av k ) = 0, (4) C πφ := 1 (ˆπ φ π φ )Ĥ = Gπ a (π a π φ ) + π φ a ( π φ) + 1 ( 4a 6 a 3 V (φ) πφ ) (aπφ ) + a4 V (φ) ( (π φ φ) i ) + π k (π φ π k ) + ω v k (π φ v k ) = 0, (43) C φ := 1 ( ˆφ φ)ĥ = Gπ a (π a φ) + π φ a ( (π φφ) + i ) + 1 ( 4a 6 a 3 V (φ) πφ ) (aφ) + a 4 V (φ) ( φ) + π k (φπ k ) + ω v k (φv k ) = 0 (44) C πk := 1 (ˆπ k π k )Ĥ = Gπ a (π a π k ) + π φ a (π φπ k ) + 1 ( 4a 6 a 3 V (φ) πφ ) (aπk ) + a 4 V (φ) (φπ k ) + π k ( π k ) + ω v k ( (π k v k ) i ) = 0, (45) C vk := 1 (ˆv k v k )Ĥ = Gπ a (π a v k ) + π φ a (π φv k ) + 1 ( 4a 6 a 3 V (φ) πφ ) (avk ) + a 4 V (φ) (φv k ) + π k ( (π kv k ) + i ) + ω v k ( v k ) = 0. (46) Note that, for convenience, the Hamiltonian operator has been assume to be Weyl orere. In orer to impose such orering, for both the inverse of the scale-factor operator â an the generic potential of the fiel V ( ˆφ) that appear insie (5), a Taylor expansion aroun their corresponing expectation values has been consiere. In the Appenix A all Poisson brackets between these ifferent constraints are explicitly compute up to linear orer in moments or orer. As can be seen, all of them are a linear combination of constraints an thus vanish on shell. Therefore, all the above constraints are first class. In summary, we have a 7-imensional space, compose by the 6 expectation values (a, π a, φ, π φ, v k, π k ), their 6 corresponing fluctuations, an their 15 inepenent correlations. The evolution of the system will be escribe by curves on that space an we have the gauge freeom to choose the time that will parameterize those curves. In principle there are no restrictions for such a choice an any combination of the expectation values, or higher-orer moments, coul o the job. Nonetheless, consiering the physical system we are escribing, it is natural to choose one of the backgroun variables as time, so that the eparametrize system still makes sense for the purely backgroun case without perturbations. There are 4 such variables (a, π a, φ, π φ ) but, in the next section, this formalism will be applie to the e Sitter universe, where 8

9 the matter fiel φ is a constant an not a ynamical variable anymore. (In such a case, an evolution in φ woul escribe the nonphysical evolution between ifferent e Sitter universes.) Hence, we will impose a gauge so that a will be our internal time variable an its canonical conjugate momentum π a will play the role of a physical Hamiltonian. Let us explain in more etail how this is actually erive from the above constraint equations. Since we have ecie to use a as our time variable, the above system of 7 equations must be unerstoo to be solve for all 7 objects relate to its momentum π a ; namely, π a itself, its fluctuation ( π a ), an its 5 correlations with the rest of the variables (Xπ a ) X πa. In particular, we note that this system is linear in the commente fluctuations an correlations. Therefore, it is straightforwar to obtain the solution for the commente seconorer moments (Xπ a ) from (41) (46) an replace it in the expectation value of the Hamiltonian constraint (40). In this way, one ens up with a fourth-orer polynomial equation for the momentum π a : G π 4 a + GA π a + i GA 1 π a + A 0 = 0, (47) where the coefficients A i are real functions of expectation values an moments unrelate to π a, that is, moments of the form (XY ) X πa,y π a. In particular, the only imaginary term appears as a multiplicative factor of the linear term in π a, an its corresponing coefficient A 1 has the form, A 1 = π φ a 3 + 4a3 V (φ). (48) Note that this coefficient is given in terms of expectation values an oes not contain any higher-orer moment. In fact, it is proportional to the time (a-) erivative of the classical constraint H 0. That is, if the Hamiltonian was time (a-) inepenent, this term woul not be present an the equation to solve for π a woul be real. The appearance of imaginary parts for the time variable is a generic feature of this approach, as has been analyze in [34 37]. Nonetheless, as will be explaine below, in our case the evolution of the universe will make the imaginary terms to be ecreasing an thus will be safe to rop them. Furthermore, the other two coefficients of the polynomial equation are given by, A 0 = a (aφ)v (φ) [ πφ 4a 6 V (φ) ] ( a) [ π a 6 φ 4a 6 V (φ) ] a 8 ( φ) V (φ) + (aπ φ) [ π 3 a 5 φ 8a 6 π φ V (φ) ] π φ a 4 ( π φ) a 4 π k (φπ k )V (φ) a 4 ω v k (φv k )V (φ) π φ [ πk a (π φ π k ) ω v k (π φ v k ) ] a π φ (φπ φ )V (φ) + a 3 [π k (aπ k ) + ω v k (av k )] [ π φ 4a 6 V (φ) ] π k( π k ) π k ω v k (v k π k ) ω 4 v ( v k ), A = π φ a + a4 V (φ) + πk + ω vk + 3π φ a 4 ( a) + 1 a ( π φ) 4π φ a 3 (aπ φ) + 1a V (φ)( a) + a 4 ( φ) V (φ) + 8a 3 (aφ)v (φ) + ω ( v k ) + ( π k ). Note that if A 0, A 1 A, the momentum woul take the form π a A /G. That is in fact the solution one woul obtain from the expectation value of the Hamiltonian constraint (40) by consiering that the system was on a momentum eigenstate an thus its fluctuation was vanishing ( π a ) = 0. Therefore, in this sense, the coefficients A 0 an A 1 contain the effects of consiering the quantum behavior of this variable. In fact, in our case, the term A 1 can be neglecte. We will consier the e Sitter case with an initial Bunch- Davies vacuum an, in principle, it is expecte that the system will not epart too much from its classical behavior. Therefore, the Hubble factor H woul be approximately constant an the scale factor will be given by an exponential of the cosmological time t as a e Ht. The rest of the variables classically go as π a e Ht, π φ 0, V (φ) H, v k π k 0. With this information at han, we obtain that the two first terms of the equation (47) behave as π 4 a A π a e 8Ht, (49) whereas the linear term also increases in time, but much slower, A 1 π a e 5Ht. Therefore, in this case, one can neglect this last term. As a sie remark, we will keep A 0, even if its approximate behavior is given by A 0 e 6Ht ( a) an, once a is chosen as time an its fluctuations requeste to be vanishing, it will be quite small. 9

10 In summary, by ropping the linear term A 1 = 0, relation (47) is a quaratic equation for πa, which can reaily be solve, πa = 1 ( ) A ± A G + 4A 0. (50) In orer to choose the correct sign insie the square root one can replace this solution, along with the form for the correlations (Xπ a ), obtaine from (4) (46), in (41) an solve that relation for the fluctuation ( π a ). It can be checke that, by choosing the minus sign in the expression above, an linearizing the result in moments, the fluctuation ( π a ) turns out to be ivergent; whereas for the plus sign the form of this fluctuation is given by an expression linear in other secon-orer moments. Therefore, in orer to be consistent with the semiclassical expansion performe in this formalism, it is necessary to choose the plus sign. In aition, the global sign that comes from taking the square root of both sies in (50) is chosen as positive so that the time flows in the increasing irection of a. In this way, one ens up with the solution, π a = 1 1/ (A + A 0) + 4A =: H a. (51) G At this point we have alreay solve the constraints for the moments relate to the momentum π a, but we have not impose any gauge an thus there is still some freeom left. In orer to iscuss an fix this gauge freeom, let us rewrite the constraints (40) (46) in the following form: C 1 := π a H a = 0, (5) C := (π a a) f = 0, (53) C 3 := ( π a ) f 3 = 0, (54) C 4 := (π a π φ ) f 4 = 0, (55) C 5 := (π a φ) f 5 = 0, (56) C 6 := (π a π k ) f 6 = 0, (57) C 7 := (π a v k ) f 7 = 0, (58) where the functions f i only epen on the expectation values an moments not relate with π a. Note that for writing C 1 a fourth-orer polynomial equation has been solve an the corresponing solution chosen as explaine above. Nonetheless, the rest of the constraints C C 7 correspons to linear combinations of constraints (41) (46) an are thus completely equivalent to those. Since that system of constraints was proven to be first-class up to linear orer with respect to secon-orer quantum moments, the same applies to the last system C 1 C 7. Therefore, we have a system of 7 first-class constraints, an following the stanar approach to constraint systems, the generalize Hamiltonian is given as C total = 7 µ i C i = 0, (59) i=1 with arbitrary Lagrange multipliers µ i. Now, we nee to impose a gauge an choose the Lagrange multipliers µ i that keeps the gauge conitions invariant. The above constraints are meant to be solve for the 7 moments relate to π a, but there are still other 6 variables relate to the egree of freeom (a, π a ); namely the expectation value a itself an its fluctuation ( a), as well as its correlations (aπ φ ), (aφ), (aπ k ), (av k ). These are the 6 gauge parameters of our system, an in orer to reuce them just to one, it is necessary to impose 5 gauge conitions. Since we want to interpret a as time, it is natural to request that it is a parameter with vanishing fluctuations an correlations. Therefore, the most natural gauge choice is given by the conitions, φ 1 = ( a) = 0, (60) φ = (aπ φ ) = 0, (61) φ 3 = (aφ) = 0, (6) φ 4 = (aπ k ) = 0, (63) φ 5 = (av k ) = 0. (64) 10

11 It can be shown that, once constraints an gauge conitions are impose, the gauge conitions φ i commute among themselves an with the physical variables (expectation values an all moments that o not involve any π a or a). Therefore, they also commute with C 1, so that it remains first class. In aition, the constraint C also remains first class but it acts trivially on the physical variables, so it oes not generate any flow on the surface efine by the constraints an gauge conitions. This means that the Hamiltonian (59) with Legenre multipliers µ 1 = 1 an µ i = 0 for i =,..., 7, conserves these gauge conitions throughout evolution. Furthermore, in the Appenix B we show that taking the constraints C 3 C 7 an the gauge-fixing conitions (60) (64) as a new whole set of constraints, their commutation matrix that is forme by their Poisson brackets an it is necessary to construct the Dirac brackets, is invertible in the general case. This proves the aequacy of the chosen gauge. In summary, we en up with the physical Hamiltonian (51), where gauge conitions (60) (64) shoul be impose. Furthermore, a plays the role of the evolution parameter an the physical variables are the 4 expectation values φ, π φ, v k an π k, together with their 4 fluctuations an the 6 correlations between them. 4.3 Power spectra The power spectrum for the variable v can be efine as, P v (k) = k3 π ( v k). (65) Nonetheless, since this quantity is not observe irectly at the CMB, it is convenient to write the power spectra of quantities relate to measurable magnitues. On the one han, scalar perturbations are present in temperature anisotropies on the CMB. They are relate to the comoving curvature perturbation R := H φ v k. (66) Therefore, from (65), the power spectrum of scalar perturbations can be written as P S (k) := P R (k) = H φ P v(k) = H k 3 φ π ( v k) = 3G ɛ k 3 π a ( v k). (67) In the last equality we have use the relation between the slow-roll parameter ɛ an the erivative of the scalar fiel φ that can be obtaine from Friemann equations. Tensor perturbations, on the other han, are expresse with the quantity h k := 3G v k. (68) a Due to the fact that tensor perturbations have two polarizations, the tensorial power spectrum will have a factor of with respect to the power spectrum of h k, P T (k) := P h (k) = 1G a P v(k) = 1 Gk3 π a ( v k). (69) These power spectra will be evaluate at late times ( kη 0 or equivalently a ) [38]. As can be seen both scalar (67) an tensorial (69) power spectra are proportional to the imensionless ratio, P := Gk3 ( v k ) a. (70) Therefore, in the next section we will analyze the behavior of this object for the case of a e Sitter universe. 5 e Sitter universe 5.1 Classical e Sitter backgroun Since we are intereste in computing quantum corrections to the e Sitter power spectra, we will first obtain the result in the approximation of QFT on classical backgrouns. This is a well-known result but, for clarity, we think it is interesting to make it explicit in terms of this moment formalism. 11

12 The e Sitter universe is characterize by a constant value of the Hubble factor H. This universe can be recovere from the moel presente above just by imposing a constant value for the potential V (φ) = H /(G), which fixes the energy of inflation. The evolution of the scale factor in terms of the conformal time η is given by a(η) = 1 Hη. (71) Since in the e Sitter case the slow-roll parameter ɛ vanishes, the frequency of the scalar sector (14) coincies with the one of the tensorial sector (17), an takes the following form, ωs = ωt = a a = k ( η = k 1 ), (7) ξ where we have introuce the imensionless time ξ := kη in orer to simplify expressions below. By replacing this form of the frequency in the equation of the fluctuation ( v k ) (33), one gets ξ 3 3 ( v k ) ξ 3 + 4ξ(ξ ) ( v k) + 8( v k ) = 0. (73) ξ This equation can be solve analytically, which yiels the following form of the fluctuation of the perturbative master variable, ( v k ) = ξ ( A(J 3/ (ξ)) + BJ 3/ (ξ)y 3/ (ξ) + C(Y 3/ (ξ)) ), (74) where J 3/ an Y 3/ are the first- an secon-type Bessel functions respectively, which can be written in terms of elementary trigonometric functions, ( ) sin(ξ) J 3/ (ξ) = cos(ξ), (75) πξ ξ ( ) cos(ξ) Y 3/ (ξ) = + sin(ξ). (76) πξ ξ In orer to fix the integration constants A, B an C, we consier the Bunch-Davies vacuum (34) (38) as initial conition which, as explaine above, correspons to the state with the minimum energy at the beginning of inflation (ξ ). With this initial state, the evolution of the three secon-orer moments is then given by ( v k ) (ξ) = 1 + ξ k ξ, (77) ( π k ) (ξ) = k ξ 4 ξ + 1 ξ 4, (78) (v k π k )(ξ) = ξ 3. (79) In particular, from these expressions, it can be explicitly seen that the uncertainty principle (37) is saturate for all ξ. For later convenience, since in the quantum-backgroun approach the time parameter will be the scale factor a, we reparametrize the fluctuation of the perturbative master variable by using the expression (71), ( v k ) (a) = a H k 3 + k. (80) In the late-time limit a, the secon term above is then suppresse an the stanar result for the power spectrum (70) is obtaine, P = Gk3 a ( v k) = H, (81) where we have efine the imensionless Hubble parameter H := l p H, with the reuce Planck length l p := G = (4π GN /3) 1/. One of the main properties of this power spectrum is its scale invariance, since it oes not epen on the wave number k. 1

13 5. Quantum e Sitter backgroun The quantum-gravity corrections for the power spectrum will be now analyze, where the ynamics of the variables is governe by the Hamiltonian (51). The amount of physical variables is now 14: the 4 expectation values (φ, π φ, v k, π k ), their 4 fluctuations an the 6 correlations between each couple. The evolution is then given by the solutions of the system of 14 ifferential equations of the form X/a = {X, H a }. In the general case the equations turn out to be quite complicate but, as will be shown below, in the e Sitter case it will be possible to fin an approximate analytic solution. As in the previous case we impose the e Sitter conitions V = H /(G) an ω = k a H. In aition, for the expectation values v k, π k an π φ, we will assume that initially they take the same value as in the previous section, that is, v k = 0, π k = 0 an π φ = 0. The first two values come from the stationary-state analysis iscusse in Section 4.1, whereas the last one is motivate from the classical behavior of this variable. It turns out that our quantum equations of motion preserve those values. That is, once replace all the commente conitions, one gets that at this orer the expectation values φ, π φ, v k, an π k are constants of motion. Therefore, concerning these variables, the stationary point an the classical e Sitter behavior are respecte. This is a nontrivial result, since the mentione equations of motion might well have been moifie by ifferent moments in such a way that, for instance, (v k = 0, π k = 0) was not a stationary point of the system anymore. In summary, once impose the above conitions, out of the 14 equations of motion, 5 are trivial. The expectation values of the fiel φ an of the perturbative variables v k, as well as of their momenta π φ an π k, are constants of motion: φ a = 0, (8) a π φ = 0, (83) a v k = 0, (84) a π k = 0. (85) Furthermore, the fluctuation of the momentum of the fiel is also conserve through evolution, Finally, the other 9 nontrivial equations take the following form: (86) a ( π φ) = 0. (87) GH a a (( φ) ) = a (φπ φ) (88) GH a a ( (φπ φ)) = 1 a ( π φ), (89) GH a a ( (φv k)) = (φπ k ) + 1 a (π φv k ), (90) GH a a ( (φπ k)) = (φv k ) ( a H k ) + 1 a (π φπ k ), (91) GH a a ( (π φv k )) = (π φ π k ), (9) GH a a ( (π φπ k )) = (π φ v k ) ( k a H ), (93) GH a a (( v k) ) = (v k π k ), (94) GH a a (( π k) ) = (v k π k ) ( k a H ), (95) GH a a ( (v kπ k )) = ( v k ) ( a H k ) + ( π k ), (96) 13

14 with the physical Hamiltonian given by, H a = 1 [ a 4 G G H + ( v k ) ( k a H ) + ( π k ) + 1 ] 1/ a ( π φ). (97) This is a highly couple an nonlinear system of ifferential equations, in particular ue to the form of this physical Hamiltonian, an analytical solutions are extremely ifficult to get. Looking at the equations, we note that it woul be natural to efine a time τ, in such a way that GH a τ = a, so that the left-han sie of all equations woul be just the erivative of the corresponing variable with respect to this time. The limit of QFT on classical backgrouns is easily recovere from the system above just by ropping all moments that involve any backgroun object. In particular, equations (88) (93) are trivially obeye in that limit. Furthermore, to obtain that limit, the physical Hamiltonian shoul also be taken as a classical variable (since it is the momentum π a ), so one shoul rop all moments from its efinition, which leas to H a = a H/G. In this case the natural time τ is the stanar conformal time τ = η = 1/(Ha) an equations (94) (96) for the fluctuations of the perturbative variables reuce to their form (30) (3). It is interesting to note that, since the fluctuation of the momentum of the fiel ( π φ ), which appears insie the Hamiltonian, is a constant of motion, equations (94) (96) for the moments relate to the perturbative egree of freeom (v k, π k ) ecouple from the rest. In particular, it is possible to write the generalization of equation (73) for the fluctuation of the perturbative master variable as follows, Ha 3 ( v k ) H a ( v k ) a 3 + 3H a a a + [ 4 ω G + H a H a a + ( Ha a ) ] ( v k ) 8 a G ah ( v k ) = 0. (98) Nevertheless, contrary to the QFT version of this equation (73), this full-quantum equation cannot be solve inepenently from the rest, since it is still couple to the fluctuation of the momentum ( π k ) an their correlation (v k π k ) through the Hamiltonian H a. Finally, if the epenence of the Hamiltonian on a was known, one coul write that equation in terms of the time τ in a very compact way, 3 τ 3 ( v k) + 4ω τ (ω( v k) ) = 0. (99) Although this might seem a simple equation to eal with, all the complexity is insie the frequency ω = ω(τ) an analytic solutions are very har to fin Approximate analytical solution for the power spectrum Since one expects to fin small changes to the result of the previous section, a first approximation to obtain the evolution of the fluctuation ( v k ) woul be to assume the behavior erive in the approximation of QFT on classical backgrouns to compute the evolution of the physical Hamiltonian H a (a). Plugging this form into (98) will ecouple this equation from the rest an one can then solve it to obtain the approximate ynamics of the fluctuation ( v k ). Therefore, replacing the form of the moments (77)-(79) in the efinition of the Hamiltonian (97) one gets, H a = a H G 1 + H k 3 (ξ4 ξ 1) + H 4 k 6 σ ξ 6. (100) In orer to absorb some of the funamental constants in this expression we have mae use of the imensionless Hubble parameter H := l p H an the time parameter ξ := k/(ah), as alreay efine above; whereas the constant fluctuation of the momentum ( π φ ) has been escribe by the quantity σ as follows: ( π φ ) =: G σ. In this way, it is shown very explicitly that H is the parameter that measures the strength of quantum-gravity effects. Since this parameter is known to be quite small [], one can perform a power expansion of the square root above an write the Hamiltonian in the following way, H a = a H G = a H G n=0 [ ( H k 3 ) n f n (ξ) (101) 1 + H 4k 3 ( ξ 4 ξ 1 ) H 4 3k 6 [ 4ξ 8 8ξ 6 ( σ + 1 ) + 4ξ + 1 ] + O ( )] H6 k 9, 14

15 where the series only contains even powers of H/k 3/. It is remarkable that the quantity σ only enters at orer H 4, which means that the actual state of the fluctuations of the matter fiel is not the most ominant part among the quantum-gravity corrections. One can then assume a similar power expansion for the fluctuation ( v k ), aroun its value in the QFT approximation (77), as follows ( v k ) = k [ ξ + 1 ξ + n=1 ( ) n H k 3 δ n (ξ)], (10) with unknown functions δ n (ξ). Replacing the last two expressions (101) (10) in equation (98), an taking into account the efinition of the imensionless time ξ = k/(ah), one can solve for the ifferent functions δ n (ξ) iteratively at each orer in H. The power spectrum can be compute just by replacing the above form of ( v k ) in its efinition (70) an it takes the form P = Gk3 a ( v k) = H [ 1 + ξ + n=1 ( H k 3 ) n ξ δ n (ξ)]. (103) In orer to obtain its exact numerical value, the ifferent δ n functions must be obtaine as inicate above. Nonetheless, from this expression one can alreay see the explicit epenence of the power spectrum on the ifferent parameters, like the Hubble factor H an the wave number k. In particular we see that the quantumgravity corrections are given by a power series in the parameter H /k 3, which makes it scale epenent an more relevant for large scales, which are relate with small values of the moe number k. We woul like to emphasize that the very same form for the first term of this series has been obtaine in several other stuies with very ifferent approximation methos, see for instance [7, 10, 11]. In orer to know whether these correction terms imply an enhancement or a suppression of power, we nee to stuy the behavior of the δ n functions in the asymptotic limit of super-hubble scales. Since the expansion parameter is very small, we will just focus on the first correction term given by δ 1, which obeys the following equation, ξ 5 δ 1 (ξ) + 4ξ 3 (ξ )δ 1(ξ) + 8ξ δ 1 (ξ) + 4ξ + 6 = 0. (104) The analytical solution of this equation is obtaine an the constants are fixe by requesting that the solution is not oscillating at the beginning of inflation (ξ ). In this way, we fin the following form for the first corrective term to the evolution of the fluctuation of the perturbative variable ( v k ) : δ 1 = 1 1ξ { Ci(ξ) [( ξ 1 ) cos(ξ) ξ sin(ξ) ] [( ξ 1 ) sin(ξ) + ξ cos(ξ) ] (π Si(ξ)) }, (105) Si an Ci being respectively the sine an cosine integral functions. The asymptotic limit for super-hubble scales (ξ 0) for the coefficient that appears multiplying the correction term of orer H /k 3 in (103) is given by, ξ δ 1 = 1 6 (11 4γ E 4 ln ξ) + O(ξ), (106) γ E being the Euler-Mascheroni constant. This expression still contains certain time epenence in the last logarithmic term, which slowly iverges. In orer to give an estimate, it is usual to evaluate this expression at horizon crossing (ξ = 1), which gives the following numerical value, ξ δ (107) Hence, up to the consiere orer, the power spectrum takes the form P H [ 1 + H ] k 3 + O( H 4 ), (108) where the above numerical value has been roune up to one. This is one of the main results of this paper an implies an enhancement of the power spectrum for large scales. In fact, the numerical value of the corrections obtaine in this analysis is in agreement with the one compute in [7], where a very ifferent semiclassical approximation scheme was use. This is certainly an important inication of the robustness of the presente result. 15