On the Cutoff Scale Identification of FLRW Cosmology in Asymptotically Safe Gravity


Abstract

We examine Friedmann–Lemaı̂tre–Robertson–Walker cosmology, incorporating quantum gravitational corrections through the functional renormalization group flow of the effective action for gravity. We solve the Einstein equation with quantum improved coupling perturbatively including the case with non-vanishing classical cosmological constant (CC) which was overlooked in the literatures. We discuss what is the suitable identification of the momentum cutoff k with time scale, and find that the choice of the Hubble parameter is suitable for vanishing CC but not so for non-vanishing CC. We suggest suitable identification in this case. The energy-scale dependent running coupling breaks the time translation symmetry and then introduces a new physical scale.

1 Introduction↩︎

The effects of quantum gravity are crucial for establishing a consistent cosmological model spanning from the early to late universe. Developing a consistent and predictive theory of quantum gravity to achieve this is a significant challenge in theoretical physics. The perturbative approach in quantum gravity has been limited due to the negative mass dimension of the gravitational constant. However, in recent years, the asymptotically safe gravity has emerged as a promising scenario. The presence of a nontrivial fixed point renders the theory safe from UV divergences [1], [2]. The use of the functional renormalization group (FRG) within the context of asymptotically safe gravity has provided new avenues to explore phenomena in black hole physics and cosmology, extending from the UV scale to the IR scale [1], [3][5]. This approach enables the examination of consistent cosmological models across the entire UV to IR range.

Building on the framework of the FRG, the central element in the asymptotically safe gravity is the effective average action, Γk[gμν]. This action is constructed to describe gravitational phenomena at a momentum scale k, accounting for the effects of quantum loops. To define the effective average action, a regulator term Rk is introduced, which suppresses the contributions from momentum modes p below k. As a result, quantum fluctuations below the FRG scale k are excluded from the effective action, while modes with p>k are fully integrated out. More explicitly the FRG equation is (1)dΓkdt=12Tr[Γk(2)+Rk]1dRkdt, where Γk(2) is the second variation of the effective average action and t=log(k/k0) with k0 being an arbitrary initial value.

However, solving fully the FRG equation in a realistic theory is challenging. To make progress, the Einstein-Hilbert truncation is often employed. The action is defined as S=d4xg116πG(R2Λ), where G is the Newton constant and Λ is the cosmological constant. The FRG ensures that the effective average action interpolates between the fundamental action, devoid of quantum corrections in the UV regime, and the quantum effective action as k0. This interpolation makes the coupling constants scale dependent, and their evolution is governed by the FRG equation [1], [6], [7]. We show in the appendix 6 that the FRG equations yield scale-dependent expressions for these coupling constants: (2)G(k)=G0[1ωG0k2+ω1G02k4+O(G03k6)],Λ(k)=Λ0[1μG0k2+μ1G02k4+O(G03k6)]+G0k4[ν+ν1G0k2+O(G02k4)], where G0 and Λ0 are the Newton coupling and cosmological constant at k=0. The dimensionless parameters ω,ω1 etc. depend explicitly on the choice of regularization scheme, and we give the results for two schemes, the optimized and exponential cutoffs there. However the results are qualitatively similar, and both G(k) and Λ(k) may be written compactly to allow the Λ0=0 solution by simply setting Λ0=0 in 2 . The solutions correspond to two distinct branches. The first branch, known as the Type IIa trajectory or the separatrix, assumes Λ0=0 [7][9]. It extends to k0, smoothly approaching the Gaussian fixed point in the IR regime. In contrast, the second branch, referred to as Type IIIa trajectories, corresponds to Λ00. This branch may not extend to k0 due to the emergence of a singularity in the IR regime, limiting its validity to a finite momentum scale. However, as discussed in [10], the nonlocal extension of the Einstein-Hilbert truncation can potentially resolve the infrared singularity in the RG trajectories with a positive cosmological constant. It is crucial to note that the solutions for Λ0=0 and Λ00 are derived separately, and we will see that these give distinct quantum corrections for the two branches. Nevertheless, the coefficients may be written uniformly like 2 . We note that the quantum correction to the cosmological constant begins at the k4 order for Λ0=0, while it starts at the k2-order for Λ00, which differs from previous reports in the literature [8], [11]. This distinction shows the need to consider quantum-improved cosmological solutions for the two cases separately.

Quantum effects can be incorporated into cosmological studies at various levels. In this work, we incorporate quantum corrections into the Einstein equations by replacing the ordinary gravitational constant and cosmological constant with their scale-dependent counterparts. This approach allows us to explore the impact of quantum gravity on cosmology. A crucial step in this setup is the identification of the IR momentum scale k with a suitable physical cutoff. Unlike black hole physics, where the cutoff scale can often be determined from physical mechanisms, cosmology lacks a direct method for such identification. Existing literature addresses this challenge by expressing k in terms of physical quantities such as particle momenta or the spacetime curvature. Most natural and common choice in cosmology is to identify k with the Hubble parameter, which is an almost unique physical scale. Various studies have explored Friedmann–Lemaı̂tre–Robertson–Walker (FLRW) cosmology choosing different cutoff scales [12][21]. In this work, we aim to determine a suitable cutoff scale motivated by the classical solutions at late times. In general, different choices for the cutoff scale can lead to varying cosmic evolution scenarios. This raises a critical question: how can we identify which cutoff choices are more viable? To address this, we investigate the quantum-corrected late-time behaviour of the Hubble parameter and the scale factor. Since quantum corrections are expected to be small at late times, the quantum-corrected solutions should closely resemble the classical solutions with minor quantum improvements. This implies that the perturbations introduced by quantum corrections must not be so large as to significantly alter the form of the solutions at late times. Through this analysis, we can identify a way to assess the viability of different cutoff scale choices.

We study the FLRW universe using the FRG-improved Einstein equations. This approach leads to a set of differential equations involving the Hubble parameter, energy density, the scale-dependent Newton’s gravitational constant, the cosmological constant, and a cutoff function. An important input for defining these equations is the requirement that the left-hand side of the Einstein equations is covariantly conserved, resulting in a modified continuity equation. This equation establishes a relationship between the matter density, the time-dependent gravitational constant, and the cosmological constant. Many cosmological studies have employed this modified continuity equation approach [16], [20], [22], [22][24]. It is important to note that this method differs from the consistency approach, which enforces the covariant conservation of both the energy-momentum tensor and the left-hand side of the Einstein equations [6], [12], [25]. In contrast, the modified continuity equation approach focuses on the conservation of the left-hand side of the Einstein equations only. Using this framework, we investigate the cosmological evolution for different cutoff choices across two trajectory branches: one with Λ0=0 and the other with Λ00. For the Λ0=0 branch, we consider cutoff identifications where k is either the classical Hubble parameter or a function of the quantum Hubble parameter, and find that the choice of the classical Hubble parameter gives reasonable cosmology, and another choice of quantum Hubble parameter gets some constraint. For the Λ00 branch, the Hubble parameter becomes constant in late time, so it is not suitable to describe the time development. We find that the choice of the inverse cosmic time successful for Λ0=0 is not suitable. Instead, we consider a time-dependent cutoff scale proportional to the inverse of the classical scale factor, k1/acl(t), which is suggested by the scalar curvature for compact space. We find that this choice ensure consistent solutions with the classical Hubble parameter and scale factor.

This paper is organized as follows. In sect. 2, we present the improved field equations within the FRG framework. Section 3 focuses on FLRW cosmology with Λ0=0. We first summarize the classical solution for Λ0=0 in subsect. 3.1. In subsect. 3.2, we use the identification of the momentum scale with the classical Hubble parameter Hcl and find reasonable behavior of the late-time cosmology with quantum corrections. In subsect. 3.3, we modify the identification to use the quantum Hubble parameter, discuss how the solution is different. In sect. 4, we consider FLRW cosmology with Λ00. We then summarize the classical solution in subsect. 4.1. In this case, we show in subsect. 4.2 that the same choice of the identification as the Λ0=0 case, namely the inverse of the cosmic time does not give good late-time behavior, which does not reproduce the classical behavior in the classical limit. So it is not a viable identification. Then in subsect. 4.3, we show that good identification is obtained by using the scale factor, ensuring consistent solutions at late times. Finally, in sect. 5, we summarize our results.

2 The Improved Field Equation↩︎

In this section, the FRG improvement is applied at the “equation level” by replacing the ordinary Newtonian gravitational constant and cosmological constant with the scale-dependent running coupling constants G(k) and Λ(k) in the Einstein equation. We take the viewpoint that these running coupling constants capture quantum effects at the leading order. The quantum improved Einstein equation can thus be written as (3)Rμν12Rgμν=Λ(t)gμν+8πG(t)Tμν . In writing down this equation, G(k) and Λ(k) are replaced with G(t) and Λ(t) to facilitate the study of cosmology, where the infrared cutoff is identified as k=k(t) with the cosmic time t.

We consider a spatially flat, homogeneous, and isotropic cosmological model described by the FLRW spacetime: (4)ds2=dt2+a2(t)(dx2+dy2+dz2), where a(t) is the scale factor. The associated geometric quantities for this spacetime are given by (5)Rtt=3a¨a,Rii=aa¨+2a˙2,R=6(a¨a+a˙2a2) . Moreover, The cosmic matter is assumed to be a perfect fluid, described by the energy-momentum tensor (6)Tμν=diag(ρ,p,p,p), where ρ is the energy density, p is the pressure and they are related by the equation of state (7)p=wρ. From these, the first modified Friedmann equation can be written as (8)H2=a˙2a2=8π3G(t)ρ+Λ(t)3 . In addition, the left-hand side of the Einstein equation is covariantly conserved, expressed as Dμ(RμνRgμν/2)=0, which is a mathematical identity in Riemannian geometry known as the Bianchi identity. Consequently, the entire right-hand side of the Einstein equation must also be covariantly conserved Dμ(Λgμν+8πGTμν)=0. This leads to the continuity equation: (9)ρ˙+3H(p+ρ)=8πρG˙+Λ˙8πG(t). Then there are two possible standpoints at this stage.

(i) Consistency condition

In the absence of quantum effects, in particular if the Newton coupling G and cosmological constant Λ are independent of time, the energy-momentum tensor is conserved DμTμν=0. So it may appear natural to require this. This leads to the ordinary continuity equation, and both the LHS and RHS of 9 should vanish [6], (10)ρ˙+3H(p+ρ)=0,8πρG˙+Λ˙=0. However in this case, the second equation gives ρ=Λ˙8πG˙. As we will discuss later, if we identify the momentum scale k with inverse of the cosmic time t, this is translated into ρ=kΛ8πkG. We see from the low energy expansions in Eq. 2 , the leading term of kG is negative (for ω>0). So the sign of the energy density ρ depends entirely on the sign of the coefficient of the leading quantum correction to the cosmological constant. For Λ00, it is negative, giving negative ρ. For Λ0=0, it is positive (for ν>0), resulting in positive ρ. The case of Λ0=0 is often considered, and discussed in detail especially in Ref. [6]. However, in our case, there is another possibility, which we employ to identify the cutoff scale in both cases.

(ii) Modified continuity equation

When we consider the quantum effects, G and Λ depend on the cosmic time. This arises because of the quantum effects of gravity. It is then expected that they also contribute to the energy-momentum. It must be G in front of the energy-momentum tensor and Λ that give such contribution in the present formulation. This suggests that we should include the time derivative of G(t) and Λ(t) in Eq. 9 as the energy-momentum from the gravity, instead of making its both sides separately vanish. This modified equation has been employed in various studies of cosmology within the asymptotic safety framework [16], [20], [22], [22][24]. This is weaker condition than (i) and gives interesting cosmology for Λ00. In this paper, we take this condition.

Let us write the modified continuity equation as (11)8πt[G(t)ρ+Λ(t)8π]=24π(1+w)HG(t)ρ. where the equation of state 7 has been used. Substituting G(t)ρ from Eq. 8 into the above equation, we obtain (12)H˙=12(3+3w)[H213Λ(t)]=α(H2Λ3), where α3(1+w)/2. Furthermore, the energy density ρ(t) can be expressed in terms of the solution for H(t) and the cosmological constant Λ(t): (13)ρ=38πG(H2Λ3). We shall solve the system of differential equations given by 12 and 13 for H(t) and ρ(t), under the condition that H(t)0 for the prescribed form of G(t) and Λ(t). It is important to note that to solve Eq. 12 for the Hubble parameter H(t), only the form of the external function Λ(t) is required. The running gravitational constant G(t) only appears in Eq. 13 when solving for ρ(t) in terms of H(t) and Λ(t). To determine Λ(t), we have to make identification of the infrared momentum scale k with the physical cutoff scale. Existing literature often expresses k in terms of scales, such as particle momenta, field strengths, or the curvature of spacetime [6]. However, there is no satisfactory principles to directly correlate the FRG scale k with the physical scale in cosmology, similar to the consistent thermodynamics for quantum improved black holes [26], [27]. For the FLRW universe, homogeneity and isotropy imply that the only physical sacle is the Hubble parameter which is determined by the cosmic time. Thus we have (14)G(t)G(k=k(t)),Λ(t)Λ(k=k(t)) . A natural choice of the cutoff identification is in terms of the Hubble parameter H(t). We still have the possibility of using classical or quantum Hubble parameters. In the next section, we first consider FLRW cosmology for Λ0=0 using the identification with the classical Hubble parameter, and then the quantum Hubble. In the latter case, we consider a general form of k(H)=ξG0(β1)/2Hβ, where β is constant.

3 FLRW Cosmology for Λ0=0↩︎

In this section, we discuss the classical and quantum-corrected late-time solutions for H(t), which, in turn, provide the solution for the scale factor a(t) for Λ0=0. This corresponds to a unique trajectory that reaches the Gaussian fixed point as k0, with Λ0=0; this is commonly referred to in the literature as the Type IIa or separatrix trajectory [7], [8]. The second branch is discussed in sect. 4.

3.1 Classical solution↩︎

The classical solution for the Hubble parameter, denoted as Hcl(t), in a single component universe specifically for radiation (w=1/3) and matter (w=0) can be obtained from Eq. 12 by setting G(t)=G0 and Λ(t)=Λ0. The classical solution for the Hubble parameter with Λ0=0 is given by (15)Hcl(t)=1α(tc), where c is an integration constant. Because the system has the time translation symmetry in Einstein equation, the free time parameter c does not carry physical significance and we can set c=0 or absorb it into the definition of the cosmic time t: (16)Hcl(t)=1αt By this, we fix the time origin in the system. The classical scale factor acl(t) is then derived by integrating the Hubble parameter as (17)acl(t)=a0t1/α . The energy density can then be expressed using the classical Hubble parameter solution from Eq. 15 as follows: (18)ρcl(t)=(3/8πG0)/(αt)2.

3.2 The identification with the classical Hubble parameter↩︎

We begin by discussing the quantum-corrected solutions and their implications. The unique physical scale in this case is the Hubble parameter. Here we consider the classical Hubble parameter in the identification. So we take (19)k=ξHcl(t)=ξtaclα(t), where ξ and ξ(=ξ/α) are dimensionless parameters of order 1. By substituting this cutoff scale into Eq. 2 , the running gravitation coupling and the cosmological coupling can be written as a function of time: (20)G(t)=G0[1ω~G0t2+ω~1G02t4+O(G03t6)],Λ(t)=t2[ν~G0t2+ν~1G02t4+O(G03t6)], where ω~ωξ2,ω~1ω1ξ4 and ν~νξ4,ν~1=ν1ξ6. Note that the time translation symmetry is no longer there by these parameters.

Substituting the cosmological constant 20 into 12 , and taking the ansatz for H(t) as (21)H(t)=1αt(1+ct+c1t2+c2t3+c3t4+) , we can determine the the above constants by comparing terms order by order in 1/t. We then obtain (22)H(t)=1αt[(1+ct+c2t2+c3t3+c4t4+)α2ν~G03t2(1+ct+c2t2+c23t2α2ν~G09t2+)α2ν~1G029t4(1+)], where c is a new constant parameter. The appearance of one undetermined parameter c is expected, as Eq. 12 is a first-order differential equation, which naturally provides one integration constant. We see the terms in the first bracket sums up to 1/(tc), quite similar to the classical solution. However we cannot remove this parameter by absorbing it to the cosmic time t because it changes the form of the cosmological constant term, in contrast to classical case. But in the absence of the quantum corrections (ν~=ν~1=0), c becomes identical with the classical parameter.

We now show that this new parameter c has physical meaning. This point was not discussed in [24]. The solution, up to fourth order in 1/t, can be re-expressed as (23)H(t)=1αt[1α2ν~G03tct(tct1)(α2ν~G03tct)2(tct1)(α2ν~G03tct)3(tct1)+O(t4)], where the critical time scale tc is defined by tc=α2ν~G0/3c . The quantum parameter ν~, which defines the critical time, is always positive for both the optimized and exponential cutoffs: specifically, ν~=ξ4/8π for the former case, as defined in Eq. 55 and ν~=ξ4ζ(3)/2π for the latter in Eq. 54 . Therefore, this positive critical time scale indicates that the leading quantum correction has an attractive effect (AdS-phase) for which the expansion rate of the FLRW universe slows down compared with the speed of classical expansion, H(t)<Hcl(t)=1/αt, for t<tc. Conversely, it becomes repulsive (dS-phase) and makes the expansion rate bigger than the classical expansion rate when t>tc. The factor (tc/t1) is not universal, but appears up to t3-order and the remaining terms are much smaller. So the critical time scale, which distinguishes the AdS and dS phases, is qualitatively correct.

The value of the parameter c can be estimated by physical observation data. We can get approximated idea about the value of c is by setting the Hubble parameter H(t0)=H0 at present time t0, with α=3/2 (w=0 matter dominated era) in the solution H(t0)=1/(α(t0c)) with the leading correction from Eq. 22 . From the condition that this should be the current time Hubble parameter H0 at t0, we get (24)c=t01αH0,t0=(1.373±0.012)×1010years,1/H0=1.45×1010years. This means that in general c is a huge positive number of the same order or one order smaller than t0, and thus the value of tc is very small. This implies the quantum correction always provides repulsive effect in late time of FLRW cosmology ttc. Note that the late time solution 22 is valid for t>c, implying 24 is a consistent condition. The quantum corrected scale factor a(t) can be calculated integrating the Hubble parameter H(t)=a˙(t)/a(t) as follows: (25)a(t)=a0t1/α[1cαt(α1)c22α2t2+αν~G06t2+O(t3)]=a0t1/α[1αν~G03tctα2(α1)ν~2G0218tc2t2+αν~G06t2+O(t3)] . Substituting Λ(t) from Eq. 20 and the Hubble parameter H(t) from Eq. 22 in Eq. 13 , we obtain the quantum corrected energy density to be (26)ρ(t)=38πα2G0t2[1+2ct+3c2t2+ω~G0α2ν~G0t2+O(t3)] . Note that, by setting ν~=ω~=0 and c=0, we recover the classical result of energy density in FLRW cosmology. Here we can have a further physical understanding for the existence of the critical time tc. The leading-order quantum correction to the cosmological constant, Λν~G0/t4>0, introduces a repulsive force. However, the integration constant c increases the energy density ρ(t) which provides an attractive force. The entire effect is the competition of these two opposite contributions, with the critical time tc characterizing the balance point between them, as can be seen in Eq. 23 .

3.3 The identification by the Hubble parameter with quantum correction↩︎

Here we consider some general type of cutoff identification as mentioned before, defined as a function of the quantum Hubble parameter H(t) in the form k(H)=ζG0(β1)/2Hβ, where ζ is a dimensionless constant. To ensure that ζ in the cutoff scale is dimensionless, unlike in [24], we introduce G0 within the cutoff scale. Three specific values of β=1/4,3/4,1 were considered in [24], and we study which of these values are viable. By substituting this cutoff choice into 2 with Λ0=0, we obtain the leading correction to the cosmological constant as (27)Λ=3γH4β+, where γ=ζ4νG02β1/3. Putting 27 in 12 , the Einstein equation 12 for H(t) becomes (28)dHdt=α(H2γH4β) . Using the relation H1(dH/dt)=a(dH/da), we can rewrite the above equation in the form (29)H24βH(γH24β)dH=αdaa . Integrating the above equation, we obtain the solution of H(t) in terms of scale factor a(t) as (30)H=aα(γaα(4β2)C0)1/(4β2). where C0 is the integration constant. To determine this, we require the scale factor a(t0)=1 and H(t0)=H0 at present time, yielding (31)C0=γH024β . By substituting the integration constant C0 into Eq. 30 , we get the final form of the quantum corrected Hubble parameter H(t): (32)H=H0aα[1+γH04β2(aα(4β2)1)]1/(4β2) . We will examine the asymptotic behaviour of the Hubble parameter at t. Our aim is to determine if there exists a critical value of the cutoff scale exponent β at which the entropy transitions from constant to divergent occurs at t, as previously discussed [24]. In the limit of infinitely large time, the scale factor a becomes for spatially flat FLRW universe. We then observe from Eq. 32 two distinct behaviours of H(t) depending on the value of β: H{H~0aα,β>1/2,γ~,β<1/2, where H~0=H0(1γH04β2)1/(4β2) and γ~=γ1/(4β2). Furthermore, the scale factor a(t) takes the form (33)a(t){(αH~0)1/αt1/α,β>1/2,eγ~t,β<1/2,. From the asymptotic behaviour of H(t) and the scale factor a(t), it is clear that β=1/2 is a critical value which changes the scale factor behaviour from a power law (β>1/2) to exponential (β<1/2).

In order to find the higher order terms in the solution of the Hubble parameter H(t) at late time, we expand H(t) as (34)H(t)={H~0aα[1+Δaα(4β2)]1/(4β2)=H~0aα[1Δ4β2aα(4β2)+],β>1/2γ~[1+Δ1aα(24β)]1/(4β2)=γ~[11Δ(4β2)aα(24β)+],β<1/2, where Δ=(H~0/H0)4β21. Integrating the above equation, we obtain the late time behaviour of the scale factor a(t) with leading order correction as (35)a(t)={(αH~0)1/αt1/α[1+Δ(4β1)α(4β2)2(αH~0)(4β2)t(4β2)+],β>1/2exp[γ~t1Δα(24β)2eαγ~(24β)t+],β<1/2. The exponential solution for the cutoff choice with β<1/2 suggests that this choice is not viable. For β<1/2, large momentum scales k cause significant perturbations, leading to exponential expansion which is completely different from the classical solution. On the other hand, for β>1/2, the behaviour follows a power law similar to the classical case, as the perturbations remain small enough to maintain this expected behaviour with the small quantum corrections. Thus, the cutoff identification is only valid for β>1/2 for late time.

For the sake of completeness, let us examine the cosmic evolution when the exponent of the momentum cutoff scale is exactly β=1/2. This gives k=ζG01/4H1/2=ζmPlH1/2 and the leading order cosmological constant becomes Λ=ζ4νH2. The Friedmann equation then takes the form (36)dHdt=α(1ζ4ν3)H2 . Solving this equation, we find the Hubble parameter as a function of cosmic time: H(t)=1α(1ζ4ν/3)(tc), where c is an arbitrary integration constant. Integrating H(t), the scale factor can be written as a(t)=a0(tc)1α(1ζ4ν/3) . We observe that the solution maintains a power-law behaviour, similar to the classical case, unlike the solutions for β<1/2. Although the cutoff scale k is large, scaling with the Planck mass, the running cosmological constant does not depend on the Planck mass and remains small. Thus, β1/2 is a viable cutoff choice and β=1/2 represents the critical value, below which β are not good choices for late time.

4 FLRW cosmology for Λ0>0↩︎

In this section, we begin by summarizing the classical solutions for H(t), a(t) and ρ(t) for the branch of trajectories where Λ=Λ0>0.

4.1 Classical solution↩︎

To calculate the classical solution for the Hubble parameter Hcl(t) in a universe dominated by a cosmological constant, the classical Einstein equation 12 with Λ=Λ0>0 yields: (37)Hcl(t)=Λ0/31+c¯exp(2αΛ0/3t)1c¯exp(2αΛ0/3t), where c¯ is the integration constant. This exact solution has two branches, depending on the sign of the free parameter c¯. Again using the time translation symmetry to shift the cosmic time t, we can set |c¯|=1. Depending on the sign of c¯, the classical solution then becomes (38)Hcl(t)={Λ0/3tanh(αΛ0/3t),sign(c¯)=1,Λ0/3coth(αΛ0/3t),sign(c¯)=1. At late times, both branches asymptotically approach a constant value of Λ0/3, consistent with the expected behaviour of a cosmological constant dominated universe. However, we select the branch with sign(c¯)=1, as it provides a continuous solution at t=0, unlike the sign(c¯)=1 branch, which has a discontinuity at t=0. For sign(c¯)=1, the exact late-time behaviour is written as (39)Hcl(t)=Λ0/3(1+2j=1(1)je2jαΛ0/3t) . The classical scale factor for the sign(c¯)=1 branch in Eq. 38 is (40)acl(t)=a0cosh1/α(αΛ0/3t)a0exp(Λ0/3t) . Finally, the energy density is obtained from Eq. 13 as (41)ρcl(t)=(Λ0/8πG0)cosh2(αΛ0/3t)(Λ0/2πG0)exp(2αΛ0/3t)0,. Thus, for Λ>0, acl at very late times (t) grows exponentially and ρcl effectively vanishes.

4.2 The identification using the cosmic time↩︎

We now analyze how the cosmic evolution of the FLRW universe differs in the different branches of FRG trajectories with Λ0>0 compared to the branch with Λ0=0.

The first problem is how to determine the identification of the momentum scale. If we follow the standard common sense, it is natural to use the Hubble scale. However it becomes constant in the late time and is not suitable for identification with the momentum scale k. So let us consider the same cutoff scale identification k=ξt, which was successful for Λ0=0. By inserting this cutoff scale in Eq. 2 , we can express the running gravitational coupling and cosmological coupling as functions of time as (42)G(t)=G0[1ω~G0t2+ω~1G02t4+O(G03t6)],Λ(t)=Λ0[1μ~G0t2+μ~1G02t4+O(G03t6)]+t2[ν~G0t2+ν~1G02t4+O(G03t6)], where ω~ωξ2,ω~1ω1ξ4,μ~μξ2,μ~1μ1ξ4 and ν~νξ4,ν~1=ν1ξ6.

We can use the Einstein equation 12 to obtain the quantum-corrected solution for H(t) as the inverse power series in the cosmic time: H(t)=Λ03[1μ~G021t23μ~G02αΛ01t3+(ν~G02Λ0μ~G0(18+α2μ~G0Λ0)8α2Λ0+μ~1G022)1t4+3Λ03(ν~G0αμ~G0(9+α2μ~G0Λ0)2α3+μ~1G02Λ0α)1t5+O(t6)], which leads to the scale factor a(t)=a0exp[Λ03t+Λ03μ~G02t1+μ~G04αt23Λ0(ν~G018μ~G0(18+α2μ~G0Λ0)72α2+μ~1G02Λ018)t3+O(t4)]=a0exp(Λ0/3t)[1+Λ03μ~G02t1+μ~G04(1α+μ~G0Λ06)t2+O(t3)]. By setting the quantum parameters to zero, we should recover the leading-order classical term of the Hubble parameter. However, the classical solution shows the exponential expansion, which can never reproduced from the power expansion in t. Another important observation is that, unlike the Λ0=0 case, this solution lacks an undetermined constant. This absence is unusual because the Einstein equation is a first-order differential equation for H(t), which typically permits a constant of integration. All of these suggest that the cutoff identification k=ξ/t is not a physically viable choice.

4.3 The identification using the scale factor↩︎

In pursuit of a consistent identification of the cutoff scale, we notice the scale factor is another physical quantity that determines the scale of our universe. So we propose to use it for the identification. Define τ=12αΛ0/3exp(2αΛ0/3t)acl2α(t), where we have used the cosmological constant to introduce dimension, and we use this for identification (43)k=ξτacl2α(t), where ξ is a constant of order 1. Here we have put the power α on the scale factor in analogy with 19 . The Einstein equation 12 takes the form (44)Λ0/3τdHdτ=12(H2Λ/3) . In the literature, this cutoff was employed in [6] within the consistency condition approach, where the energy-momentum tensor is conserved along with the right-hand side of the Einstein equation. In that approach, however, it was shown that no consistent solution exists for a flat FLRW universe, except in the presence of exotic matter.

Here we demonstrate that this cutoff choice is a viable option for the case Λ0>0 if we use the modified continuity equation, in contrast to the consistency condition approach. It is also important to note that we select Λ0 to match the dimension of the momentum of the cutoff scale k, unlike the Λ0=0 case, where G0 is chosen. This is because the quantum effects in the solution is expected to become strong around the order of Planck scale: If we had used G0 with positive power (for dimensional reason) in τ, the Planck scale in the quantum effects would cancel out from Eq. 2 , in contradiction to this expectation.

With the choice of cutoff 43 , the scale-dependent cosmological constant 2 becomes (45)Λ(τ)=Λ0[1μ~G0τ2+μ~1G02τ4+O(G03τ6)]+τ2[ν~G0τ2+ν~1G02τ4+O(G03τ6)] . By substituting this into Eq. 44 , we obtain the solution for the Hubble parameter in terms of τ as (46)H(τ)=Λ0/3[1+2cτ+2c2τ2+2c3τ3+2c4τ4+μ~G02τ2+μ~G0c2τ3+2μ~G0c23τ4+(μ~24μ~1)G0224τ4ν~G06Λ0τ4+O(τ5)] . As discussed before for Λ0=0, the time translation symmetry is broken here due to the time dependence of the cosmological term in the Einstein equation. Consequently, the new parameter c plays a significant role in obtaining the quantum-corrected solution for the Hubble parameter. However, as mentioned for the Λ0=0 case, to recover the classical solution, we should set the quantum parameters μ~, ν~ and μ~1 to zero. This allows us to absorb the free parameter within time or set it to zero. This point will become clearer in the following steps.

If we set μ~, ν~ and μ~1 to zero, the expansion of the Hubble parameter 46 becomes (47)Hcl=Λ03[1+2cτ+2c2τ2+2c3τ3+2c4τ4+]=Λ03(1+cτ)(1+cτ+c2τ2+c3τ3+c4τ4+)=Λ031+c/τ1c/τ . Redefining a new free parameter c (48)c=ln(2αΛ0/3|c|)2αΛ0/32αΛ0/3c=sign(c)e2αΛ0/3c, we obtain the exact solution for the Hubble parameter: Hcl=Λ031+sign(c)e2αΛ0/3(tc)1sign(c)e2αΛ0/3(tc)={Λ03tanh(αΛ0/3(tc)),sign(c)=1,Λ03coth(αΛ0/3(tc)),sign(c)=1. Here the parameter c can be either absorbed into the cosmic time or set to c=0. However, in the case of the quantum-corrected solution, the free parameter c, or more precisely |c|, acquires a physical meaning due to quantum corrections through Λ(t). This becomes evident when re-expressing the quantum-corrected solution as (49)H(τ)=Λ0/3[1+2|c|τ(τcτ+sign(c))+2|c|2τ2(τcτsign(c)+1)+] , where the parameter |c| gives a critical scale determined by τc=μ~G0/4|c| . If c>0, which implies sign(c)=1, the quantum correction enhances the Hubble parameter eternally. On the other hand, if c<0, meaning sign(c)=1, the infintely large late-time expansion rate remains constant at Λ0/3 as the quantum corrections become negligible. However, the value of the constant is still affected by the quantum correction terms at late time. Since the leading-order quantum correction is stronger than the subleading orders, we observe from 49 that late-time cosmic expansion is suppressed due to quantum corrections compared to the very late-time classical behavior when τ>τc, and the opposite occurs for τ<τc.

In reality, based on Eq. 48 , |c| is found to be very large since the leading-order parameter c is of the order of the present time, making τc unreasonably small. As a result, the τ>τc case dominates. While the quantum corrections are small, they produce an effect opposite to what is expected from classical behavior. Specifically, while the present universe is accelerating, the quantum-corrected Hubble parameter appears to be smaller than the classical Hubble parameter. Even in the case of Λ0=0, the quantum-corrected Hubble parameter increases, contrary to the classical expectation of a decelerating universe.

Since the Hubble parameter is defined in terms of τ as H(τ)=a1(da/dτ)(dτ/dt), the scale factor can be expressed as a(τ)=τ1/2α[1cατ(4c2+μ~G08αc22α2)1τ2+O(τ3)] . The energy density can be computed straightforwardly from Eq. 13 as ρ(τ)=Λ0c2πG0τ+Λ0(4c2+μ~G0)4πG0τ2+3Λ0c(4c2+μ~G0+4ω~G0/3)8πG0τ3+O(τ4). As expected, the entire contribution to the energy density at late times arises from quantum corrections. For c<0, the critical scale τc emerges from the competition between the first term and the term Λ0μ~G0/(4πG0τ2) which appears in the parenthesis of the second term.

It is worth noting that the quantum-corrected solution for the Hubble parameter under this cutoff choice addresses the issue of requiring an additional integration constant. Furthermore, this solution reproduces exponential subleading terms in H(t) for late times, in contrast to the power expansion terms obtained with the k=ξ/t choice. Therefore, it can be concluded from the quantum-corrected solution for the Hubble parameter that the late-time behaviour aligns with the classical solution, with minor quantum improvements. This demonstrates the viability of the proposed cutoff choice, unlike the earlier k=ξ/t choice.

5 Summary and Conclusions↩︎

In this work, we have studied a spatially flat FLRW universe in late time, taking quantum gravitational effects into account, within the framework of asymptotically safe gravity. The FRG flow of the effective average action in asymptotically safe gravity has resulted in the running of the gravitational constant and the cosmological constant under the Einstein-Hilbert truncation. We have analyzed the cosmic evolution of the FLRW cosmology at late times by improving the Einstein equations to include the scale dependence of Newton’s constant and the cosmological constant. Specifically, we have considered the conservation of the entire right-hand side of the Einstein equations, which has led to a modified continuity equation. The main focus of this work has been to identify a viable cutoff scale, k, that can yield a consistent FLRW cosmology.

In this context, we have examined cosmological solutions for two branches of FRG trajectories. The first branch has involved a power series expansion of the dimensionful cosmological constant, Λ(t), beginning with a k4 term, with Λ0=0. This trajectory has approached the Gaussian fixed point as k0 and has been referred to in the literature as the “separatrix" [7], [8]. The second branch, which may be realized in nature, starts with a k2 term in the expansion of the cosmological constant and is valid up to a certain non-zero momentum scale, where the β-function becomes singular as the dimensionless cosmological constant λ(k) comes to a”singular line” near 1/2.

For studying cosmology in these two branches, we have considered the cutoff scale k in terms of either classical or quantum Hubble parameter for the Λ0=0 case. We have chosen kHcl(t), which measures the curvature of the FLRW spacetime. Before incorporating quantum corrections through the Einstein equations, there was the time translation symmetry, implying that the integration constant in the classical solution does not carry any physical significance. However, in the quantum-corrected solution, a new constant comes into the solution with physical significance through the critical time scale tc which defines two distinct phases of the quantum-corrected Hubble parameter: the repulsive phase (dS phase) and the attractive phase (AdS phase). In the attractive phase, the leading quantum correction has introduced an attractive effect, causing the expansion rate of the FLRW universe to slow down compared to the classical expansion rate, H(t)<Hcl(t)=1/(αt) for t<tc. In the repulsive phase (dS phase), it has enhanced the expansion rate beyond the classical rate when t>tc.

Considering another possibility that the cutoff scale k may have a simple functional relationship with the quantum Hubble parameter rather than with the classical Hubble parameter, we examined the cutoff choice k(H)=ζG0(β1)/2Hβ, where β is a constant. In earlier work [24], three specific values of β namely, β=1/4,3/4 and 1 were analyzed, and the corresponding entropy generation due to quantum effects was investigated. It was observed that entropy generation diverges at late times for β=1/4, indicating that this cutoff choice is not suitable for studying cosmology. On the other hand, for β=3/4 and β=1, entropy generation converges to a constant and thus these viable cutoff choices are viable. This suggests the existence of a critical value of β that determines the appropriateness of a cutoff choice. By studying the asymptotic behaviour of the Hubble parameter at large times, we have demonstrated that the critical value is β=1/2. This value marks a transition in the scale factor’s behaviour from power-law growth (β1/2) to exponential growth (β<1/2). From the exponential behaviour of the scale factor for β<1/2, we infer that such cutoff choices are not viable, as the large momentum scale k introduces significant perturbations that entirely disrupt the classical power-law behaviour for the Λ0=0 trajectory. On the other hand, for β1/2, the behaviour follows a power law similar to the classical case. Here the perturbations remain small enough to preserve this expected behaviour with minor quantum corrections. Therefore, the cutoff identification is valid for β1/2 making it suitable for cosmological studies.

For trajectories with Λ0>0, the Hubble parameter is not suitable for the identification because it becomes constant quickly. Following the successful case for Λ0=0, we tried to use the cutoff scale k=ξ/t in terms of the cosmic time, but found that it does not provide a consistent solution because the quantum corrections to the Hubble parameter introduce power expansion terms that become more significant than the classical subleading exponential terms and cannot reproduce the classical solution in the classical limit. Moreover, unlike the Λ0=0 case, this solution with this cutoff lacks an undetermined constant term, despite the fact that the Einstein equation is a first-order differential equation for H(t). All of these suggest that the cutoff identification k=ξ/t may not represent a physically viable choice. In contrast, we can resolve this issue with k=ξ/τ which indicates that the quantum solution resembles the classical solution, along with quantum improvements in the same form. We expect that the perturbations introduced by quantum corrections should not be too large to significantly alter the form of the solutions at late times. Through this analysis, we have demonstrated a method to evaluate the viability of different cutoff scale choices for two branches of the trajectories Λ0=0 and Λ0>0.

We would like to thank Akihiro Ishibashi for valuable discussions at the early stage of this work. The work of C.M.C. was supported by the National Science and Technology Council of the R.O.C. (Taiwan) under the grant NSTC 113-2112-M-008-027. The work of R.M. was supported by the National Science and Technology Council of the R.O.C. (Taiwan) under the grant NSTC 113-2811-M-008-046. The work of N.O. was supported in part by the Grant-in-Aid for Scientific Research Fund of the JSPS (C) No. 20K03980.

6 Perturbative Solution of Running Couplings↩︎

The FRG equations for running couplings G(k) and Λ(k) for Einstein-Hilbert truncation are [1] kkG(k)=ηNG(k),kkΛ(k)=ηNΛ+k4G(k)2π[10Φ21(2λ)8Φ21(0)5ηNΦ~21(2λ)], where ηN is the anomalous dimension of gR: ηN=gB1(λ)1gB2(λ), in which the dimensionless running couplings are defined as g(k)=G(k)k2,λ(k)=Λ(k)/k2 and B1(λ),B2(λ) are abbreviations of B1(λ)=13π[5Φ11(2λ)18Φ22(2λ)4Φ11(0)6Φ22(0)],B2(λ)=16π[5Φ~11(2λ)18Φ~22(2λ)]. The threshold functions that appeared here are defined as Φnp(w)=1Γ(n)0dzzn1R(0)(z)zzR(0)(z)[z+R(0)(z)+w]p,Φ~np(w)=1Γ(n)0dzzn1R(0)(z)[z+R(0)(z)+w]p, which depend on the choice the cutoff function R(0)(z).

The FRG equations can be further reexpressed as (50)kkG(k)=k2G2(k)B1(λ)1k2G(k)B2(λ),(51)kkΛ(k)=k2G(k){k4G(k)[A1(λ)B2(λ)+A2(λ)B1(λ)]A1(λ)k2B1(λ)Λ(k)}1k2G(k)B2(λ), with two additional abbreviations A1(λ)=1π[5Φ21(2λ)4Φ21(0)],A2(λ)=52πΦ~21(2λ).

For small value of k, appropriate to late time behavior, the running couplings are assumed as (52)G(k)=j=0Gjkj,Λ(k)=j=0Λjkj. In what follows, we solve the FRG equations perturbatively.

6.1 Exponential Cutoff↩︎

In [6], the low energy expansions of running couplings were derived with the exponential cutoff R(0)(z)=zexp(z)1.

In order to solve the FRG equations perturbatively, one should expand the factor (Z2Λ(k)/k2)p with Z=z+R(0) in power of k. For Λ0=0, we write (Z2Λ(k)k2)p=Zp[1+2pΛ4Zk2+(2pΛ6Z+2p(p+1)Λ42Z2)k4+], and then the threshold functions Φnp(2λ) can be expanded as Φnp(2Λ(k)k2)=Φnp(0)+2pΛ4Φnp+1(0)k2+[2pΛ6Φnp+1(0)+2p(p+1)Λ42Φnp+2(0)]k4+, and similarly for Φ~np(2λ) as Φ~np(2Λ(k)k2)=Φ~np(0)+2pΛ4Φ~np+1(0)k2+[2pΛ6Φ~np+1(0)+2p(p+1)Λ42Φ~np+2(0)]k4+. By solving the FRG equation order by order, we determine the leading terms of the solutions 52 as G2=Φ11(0)24Φ22(0)6πG02,G4=2[Φ11(0)24Φ22(0)]2[Φ11(0)24Φ22(0)][5Φ~11(0)18Φ~22(0)]+3Φ21(0)[5Φ12(0)36Φ23(0)]72π2G03,Λ4=Φ21(0)4πG0,Λ6=[Φ11(0)24Φ22(0)][3Φ11(0)10Φ~21(0)]+30Φ21(0)Φ22(0)72π2G02, with the specific values (53)Φ10(0)=π23,Φ11(0)=π26,Φ21(0)=2ζ(3),Φ12(0)=1,Φ22(0)=1,Φ23(0)=12,Φ~11(0)=1,Φ~21(0)=1,Φ~12(0)=ln2,Φ~22(0)=12. The solutions agree with the results in [6] as, in the notation of 2 : G2=ωG02,G4=ω1G03,Λ4=νG0,Λ6=ν1G02, (54)ω=4π(1π2144),ν=ζ(3)2π,ω1=ω2ω3π13ν6π,ν1=ων+5ω6π+5ν3π.

The expansion of the factor (Z2Λ(k)/k2)p for Λ00 is written as (Z2Λ(k)k2)p=(1)p2pΛ0pk2p(1)ppΛ12pΛ0p+1k2p+1+(1)p2p+1(p(p+1)Λ122pΛ0Λ2Λ0p+2+pZΛ0p+1)k2p+2+(1)p2p+1(p(p+1)(p+2)Λ13+6pΛ0Λ36p(p+1)Λ0Λ1Λ23Λ0p+3+p(p+1)Λ1ZΛ0p+2)k2p+3+, and then it leads to the expansion of Φnp(2λ) Φnp(2Λ(k)k2)=(1)p2pΛ0pΦn0(0)k2p(1)ppΛ12pΛ0p+1Φn0(0)k2p+1+(1)p2p+1(p(p+1)Λ122pΛ0Λ2Λ0p+2Φn0(0)+pΛ0p+1Φn1(0))k2p+2+(1)p2p+1(p(p+1)(p+2)Λ13+6pΛ0Λ36p(p+1)Λ0Λ1Λ23Λ0p+3Φn0(0)+p(p+1)Λ1Λ0p+2Φn1(0))k2p+3+, and similarly for Φ~np(2λ). The FRG equations lead to 0=G1=G3=,0=Λ1=Λ3= and G2=2Φ11(0)+3Φ22(0)3πG02,G4=[2Φ11(0)+3Φ22(0)]29π2G035Φ10(0)24πG02Λ0,Λ2=2Φ11(0)+3Φ22(0)3πG0Λ0,Λ4=[2Φ11(0)+3Φ22(0)]29π2G02Λ05Φ10(0)+24Φ21(0)24πG0. With the specific values 53 we have G(k)=G0[19+π29πG0k2+[(9+π29π)25π721G0Λ0]G02k4+O(k6)],Λ(k)=Λ0[19+π29πG0k2+(9+π29π)2G02k4+O(k6)]+G0k4[5π2+144ζ(3)72π+O(k2)]. This result gives the coefficients in Eq. 2 for the exponential cutoff.

6.2 Optimized Cutoff↩︎

For the optimized cutoff [7], the cutoff function is chosen as R(0)(z)=(1z)Θ(1z), then the integration of threshold functions can be carried out Φnp(w)=(1+w)pΓ(n+1),Φ~np(w)=(1+w)pΓ(n+2). The four abbreviations in FRG equations 50 and 51 can be simplified A1=k2+8Λ2π(k22Λ),A2=5k212π(k22Λ),B1=11k418k2Λ+28Λ23π(k22Λ)2,B2=k4+10k2Λ12π(k22Λ)2.

For the case Λ0=0, the leading orders of equations lead to Λ1=Λ2=Λ3=0 and the solutions are (55)G(k)=G0[1116πG0k2+21772π2G02k4+O(G03k6)],Λ(k)=G0k4[18π+754π2G0k2+O(G02k4)].

The solutions for the case Λ00 are G(k)=G0[176πG0k2+(4936π2524πG0Λ0)G02k4+O(G03k6)],Λ(k)=Λ0[176πG0k2+4936π2G02k4(2972πG02Λ02+343216π3)G03k6+O(G04k8)]+G0k4[1724π+119144π2G0k2+O(G02k4)].

Note that these two cases are “disconnected”, i.e. the coefficients in the case Λ0=0 are not the Λ00 limit of the related coefficients in the case Λ00.

References↩︎

[1]
M. Reuter, “Nonperturbative evolution equation for quantum gravity,” Phys. Rev. D 57(1998) 971 [arXiv:hep-th/9605030 [hep-th]].
[2]
W. Souma, “Nontrivial ultraviolet fixed point in quantum gravity,” Prog. Theor. Phys. 102(1999) 181 [hep-th/9907027].
[3]
R. Percacci, “An Introduction to Covariant Quantum Gravity and Asymptotic Safety", 100 Years of General Relativity Vol. 3 (World Scientific Publishing, Singapore, 2017).
[4]
A. Eichhorn, “An asymptotically safe guide to quantum gravity and matter,” Front. Astron. Space Sci. 5(2019) 47 [arXiv:1810.07615 [hep-th]].
[5]
M. Reuter and F. Saueressig, “Quantum Gravity and the Functional Renormalization Group: The Road towards Asymptotic Safety,” Cambridge monographs on Mathematical physics.
[6]
A. Bonanno and M. Reuter, “Cosmology of the Planck era from a renormalization group for quantum gravity,” Phys. Rev. D 65(2002) 043508 [arXiv:hep-th/0106133 [hep-th]].
[7]
A. Codello, R. Percacci and C. Rahmede, “Investigating the Ultraviolet Properties of Gravity with a Wilsonian Renormalization Group Equation,” Annals Phys. 324(2009) 414 [arXiv:0805.2909 [hep-th]].
[8]
M. Reuter and F. Saueressig, “Renormalization group flow of quantum gravity in the Einstein-Hilbert truncation,” Phys. Rev. D 65(2002) 065016 [arXiv:hep-th/0110054 [hep-th]].
[9]
H. Kawai and N. Ohta, “Wave Function Renormalization in Asymptotically Safe Quantum Gravity,” [arXiv:2412.08808 [hep-th]].
[10]
P. F. Machado and F. Saueressig, “On the renormalization group flow of f(R)-gravity,” Phys. Rev. D 77(2008), 124045 [arXiv:0712.0445 [hep-th]].
[11]
M. Reuter and H. Weyer, “Quantum gravity at astrophysical distances?,” JCAP 12(2004) 001 [arXiv:hep-th/0410119 [hep-th]].
[12]
A. Bonanno and M. Reuter, “Cosmology with selfadjusting vacuum energy density from a renormalization group fixed point,” Phys. Lett. B 527(2002) 9 [arXiv:astro-ph/0106468 [astro-ph]].
[13]
A. Bonanno and M. Reuter, “Cosmological perturbations in renormalization group derived cosmologies,” Int. J. Mod. Phys. D 13(2004) 107 [arXiv:astro-ph/0210472 [astro-ph]].
[14]
E. Bentivegna, A. Bonanno and M. Reuter, “Confronting the IR fixed point cosmology with high redshift supernova data,” JCAP 01(2004), 001 [arXiv:astro-ph/0303150 [astro-ph]].
[15]
M. Reuter and F. Saueressig, “From big bang to asymptotic de Sitter: complete cosmologies in a quantum gravity framework,” JCAP 09(2005) 012 [arXiv:hep-th/0507167].
[16]
A. Bonanno, and M. Reuter, “Entropy signature of the running cosmological constant,” JCAP 0708 (2007) 024 [arXiv:0706.0174 [hep-th]].
[17]
A. Bonanno, M. Reuter, “Entropy Production during Asymptotically Safe Inflation", Entropy 13(2011) 274 [arXiv:1011.2794].
[18]
A. Bonanno and F. Saueressig, “Asymptotically safe cosmology - a status report", Comptes Rendus Physique 18, (2017) 254 [arXiv:1702.04137 [hep-th]].
[19]
R. Moti and A. Shojai, “On the cutoff identification and the quantum improvement in asymptotically safe gravity", Phys. Lett. B 793(2019) 313.
[20]
M. Hindmarsh, D. Litim, C. Rahmede, “Asymptotically safe cosmology,” JCAP 07(2019) 019.
[21]
A. Platania, “From renormalization group flows to cosmology", Front.in Phys. 8(2020) 188 [arXiv:2003.13656 [gr-qc]].
[22]
M. Reuter and H. Weyer, “Renormalization group improved gravitational actions: A Brans-Dicke approach,” Phys. Rev. D 69(2004) 104022 [arXiv:hep-th/0311196 [hep-th]].
[23]
M. Reuter and H. Weyer, “Running Newton constant, improved gravitational actions, and galaxy rotation curves,” Phys. Rev. D 70(2004) 124028 [arXiv:hep-th/0410117 [hep-th]].
[24]
R. Mandal, S. Gangopadhyay and A. Lahiri, “Cosmology with modified continuity equation in asymptotically safe gravity,” Eur. Phys. J. Plus 137(2022) 10 [arXiv:2010.09716 [gr-qc]].
[25]
R. Mandal, S. Gangopadhyay and A. Lahiri, “Cosmology of Bianchi type-I metric using renormalization group approach for quantum gravity,” Class. Quant. Grav. 37(2020) 065012 [arXiv:1906.08674 [gr-qc]].
[26]
C.-M. Chen, Y. Chen, A. Ishibashi, N. Ohta and D. Yamaguchi, “Running Newton Coupling, Scale Identification and Black Hole Thermodynamics,” Phys. Rev. D 105(2022) 106026 [arXiv:2204.09892 [hep-th]].
[27]
C.-M. Chen, Y. Chen, A. Ishibashi and N. Ohta, “Quantum improved regular Kerr black holes,” Chin. J. Phys. 92(2024) 766 [arXiv:2308.16356 [hep-th]].