On the effect of renormalization group improvement on the cosmological power spectrum

Asymptotically safe quantum gravity predicts running gravitational and cosmological constants, while it remains a meaningful quantum field theory because of the existence of a finite number of non-Gaussian ultraviolet fixed points. We have investigated the effect of such running couplings on the cosmological perturbations. We have obtained the improved Mukhanov–Sassaki equation and solved it for two models. The effect of such running of the coupling constants on the cosmological power spectrum is also studied.


Introduction
The quest to construct a covariant renormalizable quantum gravity attracted more attention in recent years. The various motivations to quantize gravity are classified by Kiefer [1] into three categories. First, the unification. After having got a worthy quantum field theory for non-gravitational interactions, unifying quantum theory and general relativity would be a logical wish. The next motivation comes from singularities of general theory of relativity. It seems that these singularities break down the theory and it is expected to be solved with an appropriate quantum theory of gravity. Finally, the last motivation arises from the fact that there are different concepts of time in quantum theory and general relativity. Time in the latter theory is a dynamical object, in contrast to quantum theory, in which it is introduced as an external parameter. Quantum gravity may modify the concept of time in one or both of these theories.
There are various approaches for building a quantum theory of gravity. In [2], these approaches are categorized in three sets. The first category is the covariant line, which is simply the quantum field theory of fluctuations over some background metric. This idea leads to higher derivative thea e-mail: ashojai@ut.ac.ir ories and string theory. The canonical line is the second approach, in which there is no need to introduce the background metric. The Hilbert space can be considered as a representation of operators corresponding to the metric or functions of it. Loop quantum gravity seems to be the latest one in this category. Finally, the name sum over histories is assigned to the third set, which contains versions of Feynman's functional integral quantization for the suggested quantum gravity theory. Discrete approaches and the spin-foam formalism belong to this set.
In any perturbative approach, which is attractive to particle physicists, offering a finite quantity is an important point in a theory. One of the difficulties in perturbative quantization of gravity is the negative (mass) dimension of gravity coupling, i.e. the Newton constant, which leads to divergences at short distances or high energies. Although renormalization of divergences necessarily does not result in a covariant formalism, it would be a vital issue in the suggested quantum gravity theory. In this context, Weinberg's asymptotic safety suggestion in 1979 [3,4], based on the existence of a finite number of non-Gaussian fixed points, for which the renormalization group (RG) flow is attracted by them at infinity, is notable. Much successful research work has been done on the existence of these fixed points (see the references in [5]). The attraction of the RG-flow to these non-Gaussian UV fixed points at short distances protects the theory from divergences in that limit.
In 1998, Reuter suggested the truncated exact renormalization group (ERG) method [6,7] for probing non-Gaussian fixed points, where the trajectories of running essential gauge couplings lie on the finite dimensional subspace of the theory space [8]. It is known that the β-functions give us running couplings of the theory which are screened or anti-screened by loop corrections at short distances. Treating couplings as fields can be found in other theories such as scalar-tensor theories (see the references in [9]), but here couplings run because of quantization in a systematic manner, i.e. RG equa-tion. This evolved constant parameter, i.e. the running coupling constant, can change the behavior of gravitational phenomena like black holes [10,11], galaxy rotation [12], CMB, etc.
In this paper, we investigate the effects of the improvement of gravitational and cosmological constants on the Mukhanov-Sassaki equation (MSE). This equation describes the growth of gauge invariant quantities constructed from quantum perturbations of metric and the inflation scalar field. These perturbations are usually considered as the primary seeds for inhomogeneities of CMB and the structure formation. Therefore investigation of the effects of RG improved couplings on the MSE would be remarkable.
The next section is dedicated to a brief introduction of ERG and various improvement methods. Then the improved MSE (IMSE) would be derived in Sect. 3. In Sect. 4, we obtain the solution of IMSE for two models, one with a scalar field responsible for the inflation, and one with a cosmological constant. Finally, in Sect. 5, the effects of this improvement on the power spectrum are studied.
It should be noted that although we are investigating the improvement of the cosmological power spectrum via the running coupling constants obtained from the asymptotic safe theory, there are other points of view to see the potential impact of renormalization in the power spectrum. For example see [13,14].

Truncated ERG in asymptotically safe gravity
Truncated ERG is one of the several methods for probing non-Gaussian fixed points of gravity theory [5]. In this approach, by truncating the scale-dependent effective action k [g αβ ] up to appropriate interaction terms, other non-effective interaction terms would be ignored. The evolution of these remaining gauge couplings, which cannot be eliminated by a redefinition of the fields, is obtained from the exact renormalization group equation (ERGE). The trajectory of RG flow is where k integrates out all the fluctuations at scale k and connects the admissible fundamental action at k→∞ = S to the conventional effective action at k→0 = . This average-like effective action at tree level describes all gravitational phenomena for each momentum of order k [15]. Indeed, the IR-cutoff R k ( p 2 ), which appears in the definition of k , eliminates the effects of fluctuations of p 2 < k 2 on RG flow and is defined by an arbitrary smooth function where R (0) (ψ) satisfies the conditions R (0) (0) = 1 and R (0) (ψ → ∞) → 0. The exponential form R (0) (ψ) = ψ exp(ψ)−1 is a common chosen form in the literature [16]. Since the multiplicity of couplings in the effective action makes the β-function intricate, truncation would project the RG flow into the finite dimensional subspace spanned by the essential couplings. This method gives a finite number of ordinary differential equations.
The Einstein-Hilbert truncation, is a common truncation for cosmological models. It is shown numerically in [17] that, for small values of the cutoff k → 0, at perturbative regime, the solutions of βfunctions for this model lead to the following power-series for dimensionful couplings where ω = 1 6π [24 2 2 (0)− 1 1 (0)], ν = 1 4π 1 2 (0) and p n (w) is the threshold function, which depends on the IR-cutoff: where g U V * = 0.27 and λ U V * = 0.36 are the non-Gaussian fix point values of the dimensionless couplings g(k) ≡ k 2 G(k) and λ(k) ≡ (k) k 2 , respectively. In the following, we are dealing with the end of the inflation epoch with k m Pl , and thus it is logical to use the perturbative regime. Hence we use Eqs. (3) and (4) as the running couplings for this period.
As is the case for any quantum field theory, running couplings are scaled with momentum k. On the other hand, in effective field theory, k uses this scaling to specify the cutoff point where the fluctuations with momenta smaller than k are ignored. In the space-time picture this cutoff momentum should be related to the inverse of some physical length. Generally speaking, using isotropy and homogeneity of FRW metric, a suitable function k(t, a(t),ȧ(t), ...) seems to be the best choice for scaling parameter. Noting that in the standard cosmology, when the universe is aged t, fluctuations of frequency smaller than 1/t are ignorable, and taking into account the dimensional consideration,s one observes that the cutoff identification is a suitable approximation. Here ξ is a positive constant of order unity. (This may also be justified for the fixed point regime, where (k) ∝ k 2 , and since we know that ∝ H 2 (t). See [18,19]). One may also use the conformal time dt/a(t) instead of t, but the difference is ignorable because of the ignorance of frequencies less than 1/t. As a result, the above cutoff identification is usually used for improvement of the cosmological models in the literature. For more details as regards the cutoff identification, the reader is referred to [15,20].
With this identification, we shall have the time dependent parameters G(k) and (k) whereω ≡ ωξ 2 andν ≡ νξ 4 . The conversion of fundamental units such as c,h and G to variable ones is a debatable issue [21]. In this regard, considering G as a coupling constant, the decision of where and how we should apply the improvement of G 0 to G(x) is an important question.
The first and simplest way to do this can be called the solution improvement, in which the parameter G 0 is replaced by G(x) in any solution of the non-improved theory. The second approach is the equation improvement. This is done at the level of the equations of motion, not the solutions. The difference between these two methods becomes bold for non-vacuum solutions and the latter seems to be more acceptable if the quantum corrections are negligible in the action. Generally speaking the improvement of the equations of motion may lead to solutions different from the former method.
By the third approach which we call the action parameters improvement [22], one means a substitution of G 0 with G(x) in the action, without adding any kinetic term for it. The improved field equations are obtained from this new action and the externally prescribed field G(x) equation comes from the RGE. If one adds some kinetic term for it, there would be no guaranty that the obtained G(x) coincides with the result of RGE. Finding a suitable kinetic term is very intricate.
Here we shall improve the equations of motion, which seems more suitable.

Improved perturbations and MSE
Metric perturbation during the inflation era and its relation to the matter inhomogeneity is an apt quantum mechanical mechanism generating initial seeds of structure formation (for a complete review see Ref. [23]). The perturbed metric has 10 df, but only six of them are physical. The others are unphysical because of gauge (coordinate) transformations. Therefore, introducing gauge invariant quantities helps us to have results that are independent of the chosen coordinate system. Among different gauge invariant quantities, the comoving curvature perturbation R would be an appropriate quantity to analyze the power spectrum of the CMB. The symmetries of perturbation equations under translations make it easier to work with Fourier modes, R q .
The evolution of R q in the non-improved theory follows the MSE [23]: where τ = t t 0 dt a(t ) is the conformal time and z, which is related to the time derivative of the scalar field by z = aφ/H , can be considered as the redshift. Since the Fourier transformation of the 2-point function of the inflation field, i.e. the power spectrum, needs the solution of this differential equation, introducing the initial conditions is unavoidable.
The wave number q is proportional to the inverse of the comoving wavelength of the perturbations, and the conformal time is the comoving horizon, hence the quantity −qτ will be the ratio of causal horizon length to the comoving wavelength of the perturbations. For −qτ 1, the wavelength of the perturbations is larger than the length of the causal horizon and the perturbations exit the horizon before the end of inflation. At the end of inflation, with growing comoving Hubble radius, 1 a H , the wavelength of the perturbations becomes smaller than the comoving Hubble radius and thus the inhomogeneities and their effect on CMB gradually become observable.
Here we are interested in investigating the effects of quantum improvements of G 0 and 0 on these observables. Therefore, at the first step, we have to study briefly the effect of an improvement on the perturbation equations. Then IMSE in the CDM universe will be obtained.

Perturbation equations with G(t) & (t)
The perturbed metric can be written as the sum of the background metricḡ αβ and the perturbation h αβ , g αβ = g αβ + h αβ . We choose the background to be a flat homogeneous and isotropic FLRW metric: with scale factor a(t), and the perturbation decomposed into (spatial) scalar, vector, and tensor modes. Hence, the line element becomes where A, B, E and F are scalars, C i and G i are divergenceless vectors and D i j is a symmetric divergenceless-traceless tensor.
Analogous to the metric decomposition, the perfect fluid energy-momentum tensor,T αβ =pḡ αβ + (p +ρ)ū αūβ is perturbed as where a bar over any quantity represents its non-perturbed value. δp and δρ are pressure and density perturbations. The velocity perturbation is decomposed as to a scalar velocity potential δu and a divergenceless vector δu i . The other parameters in δT i j , i.e. the scalar π S i , the divergenceless vector π V i and the symmetric divergencelesstraceless tensor π T i j , characterize the departure from the energy-momentum tensor of a perfect fluid.
Since the CMB inhomogeneity comes from scalar perturbation, in the following we omit the vector and tensor perturbations and just consider the scalar mode. Separating these terms on both sides of the Einstein perturbed equation, 8π and the perturbed energy-momentum conservation equation gives These are constraints on the perturbation evolutions, by which we shall obtain the IMSE in the next subsection. Note that the perturbation equations differ from the non-improved ones only in changing G 0 to G(t), as is expected in the equation of motion improvement method.

IMSE
Assuming the inhomogeneity of CMB to be produced by the perturbation of the background inflation fieldφ, the action functional would be where M 2 Pl = (8π G 0 ) − 1 2 is the reduced Planck mass. The dynamics of the model is given by the Friedmann equation for the classical background fieldφ: whereρ =ρ ϕ +ρ andp =p ϕ +p are the total density and power, respectively. Note that the dynamics of the background field is assumed to be not affected by the improvement and the effect of RG improvement on the perturbations would be investigated. Eliminating the potential of the scalar field using (21) and (22), the time evolution of the scalar non-perturbed field would be given by In order to get the evolution equation of the gauge invariant quantity, we have to fix the gauge. We take the same gauge as in [23] and consider δϕ q = 0 and B q = 0 to find R q . Hence δu vanishes. Indeed, perturbations of pressure and density in this gauge are δρ ϕ = δp ϕ = − 1 2 Eφ 2 . Therefore the energymomentum conservation (19) leads to It is worth noting that the improvement given by the solutions of RGE does not change Eqs. (14) and (16), but the conservation equation of the background field changes, as given by (24), since the scalar field exchanges energy with the quantum improved part. This energy exchange can be seen more explicitly in Eq. (22) where the dynamics of the scalar field is under the influence of the improvement term, Substituting E =Ȧ H from (16) and using (14), the evolution equation of A can be written in the form As mentioned previously, working with Fourier modes is easier, and thus in what follows we use Fourier transformed equations. It is clear that the scalar perturbations of metric would not be invariant under a transformation of x α → x α + ξ α . The mentioned transformation changes the perturbation A q toÃ q = A q + H (t)ξ 0 . To have a covariant description, usually in the non-improved theory the gauge invariant quantity R q ≡ − A q 2 + H δu q is defined. Since in this special gauge the coefficients of G(t) in Eqs. (14) and (16) vanish, the corrections due to running couplings do not affect the scalar perturbations. The corrections manifest themselves in the conservation equation of the background flow. The Hubble parameter, H (t), has a correction term in R q which leads to a second order perturbation in gauge invariant quantity and can be ignored. Therefore, the non-improved form of the gauge invariant quantity saves its invariance after improvement of the couplings. On the other hand the assumption −qτ 1 at the late time of inflation causes the ignorance of the term H δu q with respect to the scalar metric perturbation, Finally, considering the conformal time as an independent variable, R q evolves as Here the prime denotes the derivative with respect to the conformal time. This is the improved Mukhanov-Sassaki equation (IMSE). Denoting the coefficient of R q by S(τ ): we see that it can be decomposed as

Solutions of IMSE
In this section the solution of the IMSE for two cases will be obtained. First, for a model in which 0 = 0 and remaining zero as time passing. Second, for the -inflation model. We shall see that IMSE can be solved iteratively for both cases.
Considering the well-known exponential potential V (ϕ) = ge −λϕ , where g and λ are arbitrary real constants, the perturbed H = H 0 +ωH 1 and ϕ = ϕ 0 +ωϕ 1 seem to be appropriate solutions for (29) and (30) up to first order. Substituting G(t) from (6) gives H 1 and its time evolution,Ḣ 1 , up to O(1): The independent variable in the IMSE is the conformal time, τ . To determine t (τ ), the scale factor a(t) is necessary. Since H =ȧ (t) a(t) , the scale factor becomes where γ =ω Inverting this, it can be shown that Finally, to obtain the IMSE, we can use the above results to get and Therefore the IMSE is given by Since the right hand side of IMSE is just a small correction, we can solve it iteratively. At zeroth order of iteration, ignoring this term, the solution is Hankel functions H (1) α (−qτ ) and H (2) α (−qτ ). We choose the sub-horizon initial condition as exp(−iqτ ) at large −qτ . This initial condition is the result of applying WKB methods suggesting a plane wave solution for the initial conditions at early times, q a H , as in [23]. With the conditions for the CDM model of the universe, we have the same initial condition as for the non-improved model: , the first Hankel function would be an appropriate choice. After renormalizing the zeroth order solution becomes To obtain the solution of IMSE up to the next order of iteration we have to substitute the zeroth order solution in the right hand side of (40) and solve the equation. Defining in which D 1(α) (q) and D 2(α) (q) are integration constants, the solution up to first order becomes The perturbations outside the horizon are defined by the limit −qτ 1. At this limit, for α > 0, the term On the other hand, F α at −qτ 1 becomes All this leads to the following relation for the perturbations outside the horizon: The parameter α will approach 3/2 on applying the slow roll condition, 1. Since −qτ 1, in the slow roll inflation, the terms (−qτ ) 4α−3 = (−qτ ) 3 and (−qτ ) 4α−1 = (−qτ ) 5 do not have any significant contribution in F α . Beside, if we consider the case in which the non-improved solution at τ 0 satisfies the initial condition R q (0) (τ 0 ) = 0, we obtain D 1(α) (q) = −C α (−qτ 0 ) 2α−1 . The gauge invariant quantity with this initial condition becomes is a constant with the dimension of [L] −1 and C 3 2 = − 4i 6π . Clearly the second term in the parentheses is a consequence of the improvement.
Before continuing, we have to note that here (26), the IMSE, for inflation with exponential potential V (ϕ) = ge −λϕ , was solved. One may wonder how the results depend on the specific inflationary model we used. Therefore, it seems appropriate to take a look at the other models such as the model due to Starobinsky [25] or chaotic inflation [26]. The former, which seems to have a compatible description of the CMB, is the result of a specific f (R) gravity (with f (R) = R + R 2 /6M 2 ). The Starobinsky model in the Einstein frame (obtained by the conformal transformation where For the FRW background this leads to an inflationary scale factor a(t) ∼ e t 2 /12t 2 Pl . Since we are interested in the effects of the improvement at the horizon exit time, and expanding t t 0 + δt (with t 0 much larger than t Pl ), the behavior of the scale factor at times that we are interested in is given by This is a similar behavior to the one obtained by the exponential potential A similar treatment is in order for the chaotic inflation models with potential V (ϕ) = gϕ n . For example the scale factor for typical n = 2 is given as a(t) ∼ e −α(β+t/t Pl ) 2 , leading to the same behavior for the times that we are interested in. This is what is expected, since all the inflationary models are built to purge the effects of the initial conditions, although they have different scale factors at the start of the inflation era.
As a result of the above discussion, one can conclude that although the details of the improved gauge invariant quantity (47) depend on the inflationary model used, the general behavior is the same. Because of this, we confine ourselves to Eq. (47) in the next sections.

Case II: -inflation
For a -inflation model we consider the running gauge coupling, (t) and not the ϕ-field, as the agent that produces exponential expansion in the inflation era. Although this not a good model (because of giving a non-improved spectral index n s = −2), we investigate it to see how the improvement affects such a model.
For simplicity, here we use the IMSE in terms of time t, instead of the cosmological time: in which If we assume that the inflation field only produces the scalar perturbations and does not affect the expansion rate of universe, we can neglect ρ ϕ with respect to ρ . Substituting H = (t) 3 and G(t) and (t) from (6) and (7) in (52), we get In terms of the dimensionless time T = t t Pl this reads where z 1 = α 1 β 2 , z 2 = α 1 γ 2 + β 1 β 2 , z 3 = α 1 δ 2 + β 1 γ 2 + γ 1 β 2 and z 4 = α 1 2 + β 1 δ 2 + γ 1 γ 2 + δ 1 β 2 . Since at the end of inflation we can ignore z 2 , z 3 T −2 , z 4 T −4 and the constant √ 3 0 t Pl . As a result the IMSE in this case is given by where It has to be noted that the term 5 T comes from a d dt ln(a 2 3 ˙ 2 ) and has the highest contribution in this model. In contrast, −a 0 z 1 T is a consequence of the running gravitational coupling.
In a very similar way like the previous subsection the solution to (55) can be obtained: where K is a constant. It is notable that the parameter α in the Hankel function H (1) α (Q q T ) is fixed now. As for the previous model the solution in the next order of iteration can be obtained leading to the following perturbations outside the horizon: where Putting everything in place, the gauge invariant quantity would be In the next section, we shall investigate the effect of these improvements on the gauge invariant quantity on the power spectrum of perturbations for Case I, which is a realistic model.

Improved power spectrum
The freezing of a gauge invariant quantity after exiting from the horizon leads one to propose cosmological observable quantities such as the correlation function, (q), which is the Fourier transform of the 2-point function and gives us some information as regards the inhomogeneities. It is given by the relation Using the results of the previous section we can easily obtain the correlation function for the first case ( 0 = 0) as in which q is normalized to q 0 = 1/τ . The normalized correlation function both in the non-improved and improved cases are plotted in Fig. 1.
One can see in this plot that the effect of improvement coming from the running couplings is to have a little larger correlation for large wave numbers.
The relative correction of the correlation function is shown in Fig. 2 for three different values of the slow roll parameter. It shows that the correction becomes smaller as the slow roll parameter decreases. This correction is about (for q = q 0 ) 5 % for = 0.300, 5 × 10 −2 % for = 0.030 and 5 × 10 −4 % for = 0.003.
In order to obtain the power spectrum, note that the improved form of G(t) (6) quickly approaches G 0 for larger times and thus can be considered constant after inflation and in the radiation and matter dominated era. This means that the non-improved growth of perturbations works also in this case and the transfer function does not have any improvement and is given by the standard one [23]. Therefore the power spectrum is given by where T is the transfer function [23], κ = √ 2q/q E Q , q E Q M h 2 /13.6 Mpc (with M the matter density parameter, and H = h×100 km/s Mpc) is the exit wave number at radiationmatter equality time, and The result is plotted in Fig. 3 and the improved power spectrum is compared with the non-improved one.  As can be seen from this plot, the RG improvement changes slightly the power spectrum for large wave numbers. In order to see how much such an improvement is, the observational data and the obtained results are compared in Fig. 4.
The observational data are from [24], the solid curve is the non-improved power spectrum with the spectral index n s = 0.967 (κ n s T (κ)). The dashed curve is the prediction of the improved model considered in this paper with the slow roll parameter = 0.008. In both curves we choose M h = 0.16.
As is seen from the results, the effect of the improvement is to bring the tail of the graph up and leads to more accurate fit with the observed results.

Conclusions
As the asymptotically safe gravity leads to running gravitational and cosmological couplings, it is a natural question to look for its effect on the cosmological models. The effect of RG improvement of the couplings as a result of asymptotic safety becomes more important at early universe.
Here we have investigated the effect of running couplings on the curvature invariant as the seed for initial perturbations, by deriving the improved Mukhanov-Sassaki equation. We saw that this can be solved iteratively.
The gravitational coupling is almost constant in the radiation and matter dominated eras and thus the transfer function of perturbations to the power spectrum is not improved. We obtained the improved power spectrum and observed that it is slightly improved for large wave numbers. This means that the tail of the predicted power spectrum comes more closer to the observed data.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .