-
In general relativity, the event horizon of a black hole is defined as the frontier of a region of space-time, where gravity is so strong that nothing, not even light, that enters that region can ever escape. The characteristic features of black holes, including event horizons, accretion and outflow processes, singularities, and Hawking radiation, as well as the possibility of observation of their immediate environment with angular resolution comparable to the event horizon, using facilities such as the Event Horizon Telescope [1] and the GRAVITY instrument [2], make black holes useful laboratories to test general relativity in the strong field regime and, hopefully, to test ideas for quantum gravity.
The accretion process onto compact objects is relevant as it is believed to provide the energy of Quasars, accreting supermassive black holes, X-ray binaries, and Gamma-ray bursts. The simplest situation, which has reached a paradigmatic status in studies on accretion astrophysics, is spherically symmetric accretion onto a static compact object of an inviscid gas, which is at rest at infinity. In the framework of Newtonian gravity the problem was analyzed in a seminal paper by Bondi [3]. The relativistic generalization of the Bondi accretion problem was performed in [4], which was followed by numerous studies analyzing spherical accretion onto a wide variety of astrophysical objects, including black holes [5-18].
The problem of the stability of spherically symmetric accretion under linearized perturbations has been studied extensively in the literature both in the context of Newtonian gravity [19-22] and in the classical general relativistic realm [23-31]. In both of these regimes, no evidence has been found of the development of instabilities. In [19], it has been shown that in the Newtonian regime and neglecting viscosity, the perturbation has a constant amplitude for the conserved flow so that stability is ensured. In general relativity, the background solution is likewise stable with the amplitude of the perturbation being damped in time, that is, general relativistic effects enhance the stability [23, 30, 31]. A mechanism that favors stability in this case is the coupling of the infalling flow with the geometry of space-time, which acts in the manner of a dissipative effect.
In this work, we consider quantum gravity corrections to the accretion, and to its stability, onto an improved Schwarzschild black hole in the framework of the asymptotic safety scenario for quantum gravity. In this scenario, it is conjectured that gravity constitutes a fundamental theory at the non-perturbative level [32]. A basic ingredient in this construction is a non-Gaussian fixed point (NGFP) of the gravitational renormalization group (RG) flow, which controls the behavior of the theory at trans-Planckian energies, where the physical degrees of freedom interact predominantly antiscreening, which in turn renders physical quantities safe from unphysical divergences at all scales [33]. This means that the gravitational field itself is the principal carrier of relevant classical and quantum degrees of freedom, with the renormalization flow relating its infrared and ultraviolet physics. On the condition that the set of RG trajectories approaching the fixed point in the UV are parameterized by a finite number of relevant (i.e., physically motivated) running couplings, asymptotic safety is as predictive as the standard perturbatively renormalizable quantum field theory. Notwithstanding asymptotic safety defines a consistent and predictive quantum theory for gravity within the framework of quantum field theory, it remains a prediction in the sense that a rigorous existence proof for the NGFP is still lacking. There is, however, substantial evidence for the existence of such a fixed point suitable for the Weinberg's asymptotic safety scenario [34-50].
Quantum corrections to the structure of black hole spacetimes within the asymptotic safety program was initiated in [51] by constructing the so-called "RG-improved" solution to the Schwarzschild metric. Following up on these initial studies, the RG-improved Kerr metric has been discussed in [52], and the improved Schwarzschild-(A)dS and Reissner-Nordström-(A)dS metrics have been studied in [53] and [54], respectively. The inclusion of a cosmological constant is relevant, as it has been shown in [53] that this constant determines the short-distance structure of the RG-improved black holes. Further aspects of these quantum compact objects including thermodynamical properties, evaporation processes, and mini-black hole production, have been analyzed in [55-59], while quantum corrections to the accretion process onto a Schwarzschild black hole have been discussed in [60, 61].
In this paper, we review the results reported by the authors of [61], where the study of the accretion onto an improved asymptotically safe Schwarzschild metric with high-derivative terms was conducted; we furrther discuss the stability of the accretion process and contrast it with the result from GR. The interest in the improved Schwarzschild space-time with higher-derivatives in the infrared limit arises from the fact that, as shown in [62], it describes a black hole that never completely evaporates. Consequently, it makes an excellent dark matter candidate. Because our interest is to analyze the most general aspects of the problem, we will neglect viscous effects, heat transport, fluid self-gravity, and effects associated to the back-reaction of the fluid on the geometry. As done in [61], we describe the accretion process of isothermal fluids using the Hamiltonian dynamical formalism and show that, contrary to what claimed in that work, there are physical solutions for the accretion of an ultra-relavitivistic fluid in AS. Regarding the analysis of the stability, instead of adopting the approach of perturbing a scalar potential whose gradient is prescribed to be the velocity of the ideal fluid [23], we use a perturbation scheme based on the continuity equation [19, 20, 30, 31]. The stationary solution of the continuity equation gives a first integral, which within a constant geometric factor, is the matter accretion rate. Our perturbative procedure entails the perturbation of this constant of motion. This quantity, being a flux of mass can be observed and, in principle, can be measured using the present day and the projected observational instruments [1, 2].
This paper is organized as follows. In the next section, we present the RG-improved Schwarzschild metric. In Sec. III, we review the mathematics of the accretion process and present the Hamiltonian dynamical system formalism. In Sec. IV, we apply the formalism to the spherically symmetric accretion of isothermal fluids and analyze the quantum gravity effects by comparing with the GR description. In Sec. V, we discuss the implementation of the perturbation scheme based on the continuity equation, and analyze the quantum gravity corrections to the stability of the accretion by studying the behavior of standing wave and traveling wave perturbations to the mass accretion rate. Finally, in the last section, we present our conclusions.
-
We start by recalling that the classical Schwarzschild line element for a black hole with mass
M , in unitsc=G0=1 , is written asds2=−f0(r)dt2+f0(r)−1dr2+r2(dθ2+sin2θdφ2),
(1) with the metric coefficient
gtt given by:f0(r)=1−2M/r . The event horizon or Schwarzschild radius comes from solvingf0(r)=0 , which yieldsrhS=2M .In the high-derivative gravity theory, which includes the Ricci scalar square, the Ricci tensor square, the Kretschmann scalar, and running gravitational couplings in the effective action [62], the quantum gravity corrections to the metric are accounted for by promoting the classical gravitational Newton constant through a running coupling that evolves under the equations of the gravitational RG-flow. Thus, in the low energy limit the RG-improvement to the metric in Eq. (1) is obtained by doing
f0(r)→f(r)=1−2Mr(1−ξr2),
(2) where
ξ is a parameter with dimensions of length squared associated to the scale identification between the momentum scalep and the radial coordinater , which in the IR regime, takes the formp∼1/r . Notably, the quantum effects on the space-time geometry are all encoded inξ , such that forξ=0 the classical metric coefficientf0(r) is recovered. The condition for the radius of the new event horizon:f(r)=0 , takes the form of a generic cubic equationr3−2Mr2+2M3˜ξ=0,
(3) where the dimensionless parameter
˜ξ=ξ/M2 has been introduced for convenience. The only real solutionrhIR to Eq. (3) isrhIRrhS=16(2+43√8−27˜ξ+3√3√˜ξ(27˜ξ−16)+3√8−27˜ξ+3√3√˜ξ(27˜ξ−16)).
(4) The left panel in Fig. 1 shows the quotient
q≡rhIR/rhS as a function of˜ξ . For˜ξ=0 , we haverhIR=rhS as expected, while for˜ξ greater than the critical value˜ξc=16/27 , there is no horizon at all, and a naked singularity arises. This implies that each value of˜ξ in the range0⩽˜ξ⩽˜ξc selects a critical massMc=√27ξ/16 in such a way that forM>Mc there are two horizons: one inner Cauchy horizon and one outer event horizon, forM<Mc a naked singularity develops, while whenM=Mc the two horizons merge. The right panel in Fig. 1 illustrates this situation forξ=0.5 that is, forMc=0.918 . Because our aim is to discuss quantum gravity effects on the accretion onto a Schwarzschild black hole, naked space-time singularities will be ignored.Figure 1. (color online) Left panel: Dependence of
q≡rhIR/rhS on˜ξ=ξ/M2 . Rigt panel: plots of the improved metric coefficientf(r) forM<Mc (blue),M=Mc (red), andM>Mc (green), withMc=0.918 corresponding toξ=0.5 . The dashed line correspond to the classicalf0(r) forM=Mc .For future reference, we note that for
M>Mc , the outer horizon of the improved solution, upon expanding to the leading order inξ , acquires the form:rhIR= rhS−ξ/(2M) . This result exhibits the typical shifting of the horizon of RG-improved black hole solutions towards smaller values with respect to its classical counterpart. -
Here, we describe the accretion flow by adopting the Hamiltonian procedure. The formulation of the problem of accretion onto compact objects as a Hamiltonian dynamical system was proposed for the first time in [63-65].
We first recall that the phenomenon of accretion is based on two conservation laws: the continuity equation, which expresses the conservation of the particle number, and the energy-momentum conservation. These requirements are given, respectively, by the covariant derivatives
∇μ(ηuμ)=0and∇μTμν=0,
(5) where
η is the particle number density of the fluid, anduμ is the four-velocity of the fluid normalized asuμuμ=−1 .By omitting effects related to viscosity or heat transport and assuming that the fluid's energy density is sufficiently small, such that its self-gravity can be ignored, the accreting matter can be approximated as a perfect fluid described by the energy-momentum tensor
Tμν=(ϵ+p)uμuν+pgμν,
(6) where
ϵ andp are the proper energy density and the proper pressure of the fluid, respectively.For purely radial flow in the equatorial plane, the normalization condition
uμuμ=−1 producesut=−fut=√f+(ur)2
(7) where
ur denotes the radial flow velocity, andf=f(r) is the metric coefficient in the spherically symmetric line element.We assume as infalling matter a test fluid, which is both isothermal and isentropic, that is a test fluid with a barotropic equation of state
h=h(η) , whereh is the specific enthalpyh=ϵ+pη.
(8) This fluid is described by the thermodynamical equations
dp=ηdh,dϵ=hdη,
(9) such that the square of the local speed of sound is written as
a2=∂p∂ϵ|ds=0=dpdϵ=ηhdhdη=dlnhdlnη.
(10) Then, as it is well known, for the spherically symmetric stationary accretion ,the conservation laws in Eq. (3) yield
ηurr2=C1,
(11) and
hut=C2,
(12) where
C1 andC2 are arbitrary integration constants.To describe the accretion process as a two-dimensional Hamiltonian dynamical system, it is useful to define the three-velocity
v of a fluid element as measured by a locally static observer asv=dl/dτ , wheredl=dr/√f anddτ=√fdt are the proper radial distance and the proper time, respectively. It is important to remark thatv is defined outside the horizon, that is for an observer for whom the time coordinate is timelike [63-65]. Usingur=dr/dτ ,ut=dt/dτ and Eq. (7) we obtainv2=(ur)2f+(ur)2=(ur)2(ut)2,
(13) from which
(ur)2=fv21−v2,and(ut)2=f1−v2.
(14) As the Hamiltonian for the flow fluid, we can choose either of the two integrals of motion
C1,C2 , or even any combination of them. We chooseC22 as our HamiltonianH and fix the dynamical variables of the system to ber andv . Using Eq. (14), we then haveH(r,v)=h2(r,v)f(r)1−v2,
(15) and the dynamical system is defined by
˙r=∂H∂v,˙v=−∂H∂r,
(16) where the dot denotes time derivative. After employing Eq. (10), these two equations become
˙r=2fh2v(1−v2)2(v2−a2),
(17) ˙v=−h2r(1−v2)[(1−a2)rdfdr−4fa2].
(18) The critical points of the dynamical system, which coincides with the sonic points of the fluid flow, are the points
(rc,vc) , where both˙r and˙v are zero. Assuming thath≠0 for all values ofr>rhIR where alsof≠0 , it follows that at the critical pointsv2c=a2c,rc(1−a2c)fc,rc=4fca2c,
(19) where
fc=f(r)|rc andfc,rc=(df/dr)|rc . -
In this section, according to [66, 67], we will obtain an adequate expression for
h(r,v) in Eq. (15) when the accreted fluid is isothermal. This in turn provides a neat expression for the HamiltonianH(r,v) .The barotropic equation of state for an isentropic fluid can be expressed as
ϵ=ϵ(η) . Moreover, for an isothermal fluid, the equation of state readsp=kϵ , where the constantk is the parameter state:0<k<1 . Note that the definition of the speed of sound in Eq. (10) leads toa2=k , and so the speed of sound remains constant through the accretion process. Now, from Eqs. (9) we haveh=dϵdη,
(20) and
dpdη=ηdhdη=d2ϵdη2,
(21) which, upon integration, produces
ηdϵdη−ϵ(η)=kϵ(η),
(22) where we used
p(η)=kϵ(η) . Integration of this last equation yieldsϵ(η)=Cηk+1=ϵcηk+1cηk+1,
(23) where the integration constant
C has been chosen such that Eqs. (8) and (20) give rise to the same expression forh h(r,η)=(k+1)ϵcηc(ηηc)k,
(24) An expression for
η/ηc can be obtained from the constant of integrationC1 in Eq. (11). In fact, using Eqs. (14) and (19), we can writeC21=r4η2fv21−v2=r4cη2cfcv2c1−v2c=r5cη2cfc,rc4,
(25) such that
(ηηc)k=(r5cfc,rc41−v2r4fv2)k/2.
(26) Substituting this into Eq. (24), we obtain
h2=K(1−v2r4fv2)k,
(27) with the constant
K given byK=(r5cfc,rc/4)k[(k+1)ϵc/ηc]2 .Redefining the Hamiltonian as
H=H/K and using Eq. (26),H finally acquires the formH(r,v)=[f(r)]1−kr4kv2k(1−v2)1−k.
(28) It must be stressed that, due to the definition of the three-velocity
v , this expression forH is valid for an observer outside the horizon. This applies both for the classical Schwarzschild metric and for the improved AS metric.We now focus on four types of isothermal accretion: ultra-stiff fluid (
k=1 ), ultra-relativistic fluid (k=1/2 ), radiation fluid (k=1/3 ), and sub-relativistic fluid (k=1/4 ). For the purpose of comparison with the analysis in [61], we will takeξ=0.5 andM=1>Mc=0.918 in the expressions for the quantumH(AS) and the classicalH(GR) Hamiltonians. Even though this amounts to consider a black hole solution with two horizons, the quantum effects on the accretion we will describe below are present for all values of the parameterξ in the range0⩽ξ⩽(16/27)M2 , with the orbits in the contour plot of the quantum Hamiltonian going continuously towards the orbits of its classical counterpart in the limitξ→0 , as we will show ahead. -
The equation of state for ultra-stiff fluid reads
p=ϵ , which impliesa2=1=v2c . Eq. (19) immediately yieldsfc=0 , and the Hamiltonian in AS coincides with the one in GRH(AS)=H(GR)=1r4v2,
(29) with the only difference being in the location of the horizon that, in both cases, coincides with the location of the critical point:
r(AS)c=rhIR andr(GR)c=rhS .Contour plots of
H(AS) (continuous lines) and ofH(GR) (dashed line) are shown in Fig. 2, where we have only included selected orbits associated to physical flow (|v|<1 ). Here, we retrieve the findings in [61] concerning the contour lines forH(AS) . Notably, both for infalling matter (−1<v<0 ) and for fluid outflow (0<v<1 ), and for fixed values of the radial coordinater withr>rhS , the flow velocity is consistently larger in GR than in AS.Figure 2. (color online) Contour lines of
H(AS)=H(AS)c (black) andH(AS)=H(AS)c±0.06 (red) for accretion of ultra-stiff fluid (k=1 ). The green dashed curve is the contour line ofH(GR)=H(GR)c±0.06 .H(AS,GR)c stands forH(AS,GR)(r)|rc . For all the orbits:M=1 ,ξ=0.5 . The vertical dashed lines locate the horizonsrhIR (left) andrhS (right). The horizontal dashed lines locate the speed of sounda=√k .We note that we have restricted the comparison between classical and quantum accretion and outflow to the region
r>rhS in the plots, because this is the region for which the comparison makes sense. -
In the case of an ultra-relativistic fluid we have
a2=1/2=v2c , and the Hamiltonians readsH(AS)=√1−2Mr+2Mξr3r2v√1−v2,H(GR)=√1−2Mrr2v√1−v2.
(30) For the critical radii, we obtain
r(AS)c≃52M−1425ξM,r(GR)c=52M,
(31) where
r(AS)c is given at leading order inξ .In Fig. 3, the contour plots of
H(AS) (continuous lines) and ofH(GR) (dashed line) are shown. For readability of the figure, we have not depicted the contour lineH(GR)=H(GR)c . The important fact here is the existence of orbits associated to physical flow (|v|<1 ) in AS, which includes the subsonic regime (−vc<v<vc ), supersonic regime (−v<−vc andv>vc ), and transonic regime associated to the solution passing through the sonic point. This is the opposite to the findings by the authors in [61], where no physical significance was established for accretion of an ultra-relativistic fluid in AS. Two quantum gravity effects can be identified: a shifting of the orbits towards the black hole and an increase of the maximum flow velocity, which occurs atr=rc , as compared with the case in GR.Figure 3. (color online) Contour lines of
H(AS)=H(AS)c (black) andH(AS)=H(AS)c±0.06 (red) for accretion of ultra-relativistic fluid:k=1/2 . The green dashed curve is the contour line ofH(GR)=H(GR)c±0.06 . For all the orbits:M=1 ,ξ=0.5 . The vertical dashed lines locate the horizonsrhIR (left) andrhS (right). The horizontal dashed lines locate the speed of sounda=√k . -
For a radiation fluid, the critical radii are given by
r(AS)c≃3M−59ξM,r(GR)c=3M,
(32) while the AS and RG Hamiltonians takes respectively the form
H(AS)=(1−2Mr+2Mξr3)2/3r4/3v2/3(1−v2)2/3,H(GR)=(1−2Mr)2/3r4/3v2/3(1−v2)2/3.
(33) The orbits of
H(AS) (continuous lines) and ofH(GR) (dashed line) are shown in Fig. 4. Here, we also recover the results in [61] concerning the contour lines forH(AS) . Notably, the overall behavior of the accretion and the outflow are the same as in the casek=1/2 , which implies the same effects coming from quantum gravity discussed above.Figure 4. (color online) Contour lines of
H(AS)=H(AS)c (black) andH(AS)=H(AS)c±0.06 (red) for accretion of radiation:k=1/3 . The green dashed curve is the contour line ofH(GR)=H(GR)c±0.06 . For all the orbits:M=1 ,ξ=0.5 . The vertical dashed lines locate the horizonsrhIR (left) andrhS (right). The horizontal dashed lines locate the speed of sounda=√k . -
For sub-relativistic fluid, we have
r(AS)c≃72M−2649ξM,r(GR)c=72M
(34) and
H(AS)=(1−2Mr+2Mξr3)3/4rv1/2(1−v2)3/4,H(GR)=(1−2Mr)3/4rv1/2(1−v2)3/4.
(35) The contour plots for
H(AS) (continuous lines) and ofH(GR) (dashed line) in Fig. 5 show that, once again, the overall behavior of the contour lines is the same as in the two previous cases, with the same global quantum gravity effects on the accretion and outflow processes. Notwithstanding, we observe that for large values ofr (r>rc) , the quantum effects are milder for sub-relativistic matter than for ultra-relativistic and radiation fluids.Figure 5. (color online) Contour lines of
H(AS)=H(AS)c (black) andH(AS)=H(AS)c±0.06 (red) for accretion of sub-relativistic fluid:k=1/4 . The green dashed curve is the contour line ofH(GR)=H(GR)c±0.06 . For all the orbits:M=1 ,ξ=0.5 . The vertical dashed lines locate the horizonsrhIR (left) andrhS (right). The horizontal dashed lines locate the speed of sounda=√k .At this point, it is important to show more explicitly how quantum effects modify accretion, that is, the dependence of the accretion flow on the parameter
ξ . This is displayed in Fig. 6 for the case of accretion of ultra-relativistic fluid (k=1/2 ) and for subsonic flow. With the increase in the value ofξ (from right to left in the figure), the effects of shifting of the orbit and of increasing of the maximum flow velocity, already mentioned, become increasingly pronounced, whereas in the opposite limit,ξ→0 , the classical orbit is retrieved. -
Let us start, following [30, 31], by writing the equations describing the accretion process in a more appropriate manner for the stability analysis.
First, the explicit form of the continuity equation reads
∂t(ηut)+r−2∂r(ηurr2)=0.
(36) To solve for
η andur , we use the equation for energy-momentum conservation and the thermodynamic equation for mass-energy conservationd(ϵ/η)+Pd(1/η)=Tds . Along with the condition of constant entropyds=0 , and with the square of the speed of sound given by Eq. (10), the condition of energy-momentum conservation readsut∂tut+ur∂rut+f(∂rf)urut+a2η[((ut)2−f)∂tη+urut∂rη]=0.
(37) Once the expression for the speed of sound
a of the particular test fluid is known and after fixing the time derivatives to zero, the solution to the system of coupled equations (36) and (37) provides us with the stationary fieldsur(r) ,η(r) , anda(r) and allow, in principle, the quantitative study of the spherical accretion onto a static spherically symmetric black hole. Assuming that the flow is smooth at all points of space-time, the stationary solutions of Eqs. (36) and (37) are written as4πˉμηurr2=˙m,
(38) and
ddr[ln(fut)]+a2ηdηdr=0,
(39) respectively, where the constant of motion
˙m is to be identified as the matter flow rate.The use of the perturbation scheme based on the continuity equation starts by considering linear perturbations to the stationary velocity and particle number density fields
ur(r,t)=ur(r)+ur′(r,t) andη(r,t)=η(r)+η′(r,t) , where primed quantities stand for small time-dependent perturbations. We now define the variableψ=ηurr2 , whose stationary valueψ(r) coincides, except for the constant4πˉμ , with the steady mass accretion rate given by Eq. (38).The first order perturbation to
ψ around stationary values isψ′(r,t)=[ur(r)η′(r,t)+η(r)ur′(r,t)]r2,
(40) and Eq. (36) acquires the form
ut∂tη′+ηurf2ut∂tur′=−1r2∂rψ′.
(41) In general, the first order perturbation to the square of the speed of sound is given by
a′2=a2+da2dηη′,
(42) The time evolution of
η′ andur′ follows from Eqs. (40) and (41) and is given, respectively, by∂tη′=−1r2(urf∂tψ′+fut∂rψ′),
(43) ∂tur′=futηr2(ut∂tψ′+ur∂rψ′).
(44) Eq. (37) is written in terms of perturbed quantities in the form
ut[ura2η∂tur′+∂tη′]+∂r(urur′)+2urur′a2η∂rη+(fut)2∂r(a2ηη′)=0.
(45) Taking the time derivative of this last equation and substituting for the time derivatives of
η′ andur′ from Eqs. (43) and (44), we obtain the differential equation obeyed by the perturbation to the mass accretion rateψ′ ∂t(ηhtt∂tψ′+ηhtr∂rψ′)+∂r(ηhrt∂tψ′+ηhrr∂rψ′)=(1−2a2)(hrt∂tψ′+hrr∂rψ′)dηdr,
(46) where the coefficients
hαβ are given byhtt=ur(fut)f2[(fut)2+ur−(ur)2a2],htr=hrt=(ur)2(fut)2f(1−a2),hrr=ur(fut)[(ur)2−(fut)2a2].
(47) -
Concerning the study of perturbations in the form of a standing wave, it is has been already pointed out in the literature that, because an event horizon instead of a physical surface occurs in a black hole, difficulties appear in fixing an appropriate inner boundary condition. Regularity of the flow at the black hole horizon singles out a unique solution, the Bondi one, which turns out to be transonic after crossing of the sonic point. However, the standing wave perturbation must vanish even in the supersonic regime, but there is no physical mechanism that allows to impose such a constraint. Because the standing wave analysis requires cancellation at the boundaries and the continuity of the solution, we must restrict ourselves to completely subsonic flows even though these may not entirely be representative of the precise manner of the infalling process. Consequently, in this section, we study the stability of a subsonic flow by assuming the trial standing wave perturbation
ψ′(r,t)=ζ(r)exp(−iωt),
(48) which, when substituted into Eq. (46), provides
ω2httζ2+iω{ddr(htrζ2)−2hrtζ2ddr[ln(fut)]}+hrr(fut)2dζdrddr[ζ(fut)2]−ddr(hrrζdζdr)=0.
(49) Integrating Eq. (49) over the radial coordinate with the integrated terms vanishing at the boundaries, a dispersion relation for
ω is obtainedAω2−2iBω+C=0
(50) where
A=∫httζ2dr,B=∫hrtζ2ddr[ln(fut)]dr,C=∫hrr(fut)2dζdrddr[ζ(fut)2]dr.
(51) The roots of Eq. (50) are given by
ω=iBA±i√B2A2+CA.
(52) Clearly, the sign of the discriminant of the relation dispersion determines the stability of the standing wave. Now, from Eq. (47), it follows that
hrt>0 . Further, taking into account that for infalling matter(dη/dr)<0 , and referring to Eq. (39), we can conclude thatB>0 . Again, becauseut=√f+(ur)2/f and for infalling fluidur<0 , we havehtt<0 , which meansA<0 . Finally, it is easily verified that for subsonic flow,hrr>0 . Consequently,(B/A)<0 and(C/A)<0 , which implies whether an oscillatory and damped in time perturbation when|C/A|>(B/A) , or an overdamped perturbation if|C/A|<(B/A) . Thus, in any case, the stationary solution will be stable.Our goal, however, is to determine if quantum gravity effects enhance or diminish the dissipative effect associated to the coupling of the flow with the geometry of space-time. Nevertheless, since the integrands in the coefficients in Eq. (51) depend in a complicated manner on the classical and quantum functions
f0(r) andf(r) , the solution forω is rather impractical for our purpose in the sense that it requires a considerable numerical effort after choosing a suitable mathematical distribution for modeling the amplitude of the standing wave. In the next section, we will see that this aim can be achieved more easily by analyzing a disturbance in the form of traveling wave. -
In this case, the perturbation will be modeled as a high-frequency traveling wave with a wavelength significantly smaller than the horizon radius of the black hole [19]. Then, the spatial part
ζ(r) of the perturbationψ′(r,t) is written as a power series inω in the form [30, 31]ζω(r)=exp[∞∑l=−1ω−lkl(r)].
(53) After replacing into Eq. (49), the coefficients of
ω2 andω can be collected and equated, each to zero. In this manner, two first order differential equations are obtained fork−1 andk0 , respectively. Settingthe coefficients ofω0 likewise to zero, a second order differential equation is obtained fork1 , in terms ofk−1 andk0 . The solutions fork−1 andk0 readk−1=i∫(hrr)−1[htr±√(htr)2−hrrhtt]dr,
(54) and
k0=ln{(fut)2[√(htr)2−hrrhtt]−1}1/2,
(55) while the differential equation for
k1 is2(hrrdk−1dr−ihtr)dk1dr+ddr(hrrdk0dr)+hrrdk0drddr[k0−2ln(fut)]=0.
(56) The power series can be truncated after these first three terms, as indicated by their asymptotic behavior given by
k−1∼r ,k0∼lnr , andk1∼r−1 , which likewise shows that the self-consistency requirementω−l|kl(r)|≫ω−(l+1)|kl+1(r)| for the convergence of the series is fulfilled. Thus, we can writeζω(r)≈exp[ω1k−1(r)+k0(r)+ω−1k1(r)],
(57) It is then clear that
k−1 andk1 only contribute to the phase of the traveling wave, and that the most important contribution to the amplitude of the perturbationψ′(r,t) comes fromk0 . A direct calculation yields|ζω(r)|=χ|exp[k0(r)]|=χ|(fut)2√(htr)2−hrrhtt|1/2=χ[(ut)2(ur)2a2]1/4=χ(1v2a2)1/4,
(58) where
χ is an arbitrary small real constant, and we have used Eq. (13). It must be noted that, due to the fact that for a spherically symmetric accretion,v never equals to zero, and sincea is also different from zero,ζω never diverges, such that the background solution is stable.The dependence of
|ζω(r)| on the metric coefficientf makes the effect of the space-time geometry on the perturbation evident and determines the coupling between the accretion flow and the curvature of space-time. The differences between these effects in AS and in GR for the accretion of isothermal fluids are already evident from Eq. (58), because a higher value of the three-velocity of the accretionv is always accompanied by a lower value of|ζω(r)| and vice versa. The variation of|ζω(r)| with the radial distance to the black hole can be calculated by writing the Hamiltonian for each type of isothermal fluid as a polynomial equation in the three-velocityv . After fixing an appropriate value for the Hamiltonian, this equation can be solved, and the result can be inserted into Eq. (58). These polynomial equations for each value of the state parameterk , except fork=1 for which the expression forv2 is immediate, readw2−w+f(r)H2(r)r4=0(k=1/2),
(59) where
w=v2 ;v3−v+f(r)H3/2(r)r2=0(k=1/3),
(60) and
x4−x−f(r)H4/3(r)r4/3=0(k=1/4),
(61) where
x=v2/3 .From the solutions to previous equations for each value of
k , we selected only the one corresponding to infalling or outflowing matter with velocity approaching zero at the infinite spatial. By replacing the chosen solutions into Eq. (58), we can construct the plots shown in Fig. 7 for subsonic flow (similar plots can also be done for transonic and supersonic flows). In this figure, the upper left panel shows that for an ultra-stiff fluid, the amplitude of the perturbation is larger in AS as compared with the amplitude in GR for all values ofr (r>rhS ) that is, for the same radial distance a locally static observer measures a larger amplitude in AS than in GR, thus leading to a lower stability of the accretion process in AS. This indicates that, for ultra-stiff fluid, the coupling of the infalling matter with the space-time geometry is weaker in AS than in GR, a result which is undoubtedly associated to the anti-screening character of the gravitational interaction in AS. For the remaining cases of isothermal accretion, namelyk=1/2 ,k=1/3 , andk=1/4 , the upper right panel, lower left panel, and lower right panel in Fig. 7, respectively, show that the amplitude of the perturbation in AS is detracted from its value in GR ifrhS<r<rcross and it is instead enhanced ifr>rcross , where the value ofrcross for each type of isothermal fluid is simply obtained by solving the equation forr , which equals the quantum amplitude with the classical one. This means that the coupling of the fluid with the space-time curvature, which acts on the perturbation in the manner of a dissipative effect, is stronger in AS than in GR forr<rcross and is weaker forr>rcross . Notably, for sub-relativistic fluid, the main effect is the reduction of the amplitude of the traveling wave perturbation with the consequent increase in the coupling of the fluid with the curvature of space-time. This result has been already reported in [68].Figure 7. (color online) Behavior of amplitude of traveling wave perturbation
|ζω(r)| in AS (red) versus GR (green) for accretion of ultra-stiff fluid (upper left panel), ultra-relativistic fluid (upper right panel), radiation fluid (lower left panel), and sub-relativistic fluid (lower right panel). The vertical dashed line stands for the horizon radiusrhS of the classical Schwarzschild black hole. The values of the parameters areHAS=HASc+0.01 ,HGR=HGRc+0.01 ,M=1 ,ξ=0.5 ,χ=0.1 -
In this work, we reviewed the problem previously analyzed in Ref. [61] regarding the steady accretion onto a R-G improved Schwarzschild black hole in AS theory with higher derivatives. In this sense, the accretion process is considered as a quantum correction to the known spherically symmetric accretion in general relativity. As in [61], we used an isothermal fluid as our test fluid, and we have specialized to the cases: ultra-stiff fluid (
k=1 ), ultra-relativistic fluid (k=1/2 ), radiation fluid (k=1/3 ), and sub-relativistic fluid (k=1/4 ). Because our interest has been to describe the most general aspects of the problem, we neglected viscous effects, heat transport, fluid self-gravity, and effects associated to the back-reaction of the fluid on the geometry. By using the Hamiltonian dynamical system procedure introduced in [63-65], we have recovered the results by the authors in [61] for the accretion of ultra-stiff, radiation, and sub-relativistic fluids. However, in contrast to what was established in that study, we found that accretion of ultra-relativistic fluids are a physical possibility in AS. In fact, for this type of fluid, as for the others considered here, we found not only a shifting of the orbits of the Hamiltonian in the(v,r) plane towards the central object and an enhancement of the flow velocity near the black hole as compared to what happens in GR, but also subsonic, supersonic, and transonic regimes. This is an important result, as there seems to be no physical reason for the no physical significance of the accretion of an ultra-relativistic fluid in AS. As expected, the quantum effects on the accretion are fully governed by the parameterξ in the RG-improved metric coefficientf(r) in Eq. (2) with the classical behavior recovered in the limitξ→0 .We also analyzed the issue of the stability of the accretion and compared the results with the known stability of the accretion flow in classical general relativity. Using a perturbative procedure based on the continuity equation, the mass accretion rate of the infalling matter, considered as an isothermal fluid, has been subjected to small linear perturbations in the form of a standing and a traveling wave. For the standing wave perturbation, we conclude that, as in the general relativistic realm, the disturbance is damped in time, such that the accretion is stable. As for the traveling wave perturbation a simple criterion, based on the behavior of the coefficient
k0 in the series expansion of the spatial part of the perturbation, has allowed us to compare between the quantum and the classical frameworks, and to reach to the conclusion that, where quantum gravity effects will have to be accounted for, the amplitude of the traveling wave perturbation for ultra-stiff fluid becomes enhanced as compared with the classical one. This indicates that the coupling of this type of fluid with the space-time curvature is weaker in AS than in GR. For ultra-relativistic, radiation, and sub-relativistic fluids, the amplitude of the perturbation is reduced or enhanced depending on whether the local observer is located in the immediate neighborhood of the black hole horizon or not, respectively. The value of the radial coordinatercross at which this transition takes place can be easily calculated by solving the equation forr , which equates the classical amplitude with the quantum one. This means that the coupling of the accreting fluid with the geometry of the space-time, which is responsible for the damping of the perturbation, is stronger in the quantum case than in the classical one forr<rcross and weaker forr>rcross . Because the disturbed quantity is the mass accretion rate, a physical quantity that is measurable in principle, it would be highly interesting to study the possible observable effects that could be used to test asymptotic safety in the future. -
We acknowledge financial support from COLCIENCIAS through the project with code 50754 associated to the CONVOCATORIA 757 PARA DOCTORADOS NACIONALES 2016. We also acknowledge financial support from Universidad Nacional de Colombia through the project with code 36031 within the program for the institutional strengthening of scientific research.
Black holes in asymptotic safety with higher derivatives: accretion and stability analysis
- Received Date: 2021-02-23
- Available Online: 2021-07-15
Abstract: We review steady spherically symmetric accretion onto a renormalization group improved Schwarzschild space-time, which is a solution to an asymptotically safe theory (AS) containing high-derivative terms. We use a Hamiltonian dynamical system approach for the analysis of the accretion of four types of isothermal test fluids: ultra-stiff fluid, ultra-relativistic fluid, radiation fluid, and sub-relativistic fluid. An important outcome of our study is that, contrary to the claim in a recent work, there are physical solutions for the accretion of an ultra-relativistic fluid in AS, which include subsonic, supersonic, and transonic regimes. Furthermore, we study quantum corrections to the known stability of the accretion in general relativity (GR). To this end, we use a perturbative procedure based on the continuity equation with the mass accretion rate being the perturbed quantity. Two classes of perturbations are studied: standing and traveling waves. We find that quantum gravity effects either enhance or diminish the stability of the accretion depending on the type of test fluid and the radial distance to the central object.