-
Einstein's theory of general relativity (GR) is an effective gravitational theory applicable to physical scales much greater than the Planck scale but smaller than the Hubble scale. However, near the Planck scale, the GR theory is ammended in terms of the action by certain averaging effects such as those resulting from the renormalization group [1]. The renormalization group improved (RGI) approach is based on the UV complete and non-perturbative renormalization to quantum gravity, which is manifested in the asymptotically safe gravity as well [2, 3]. Consequently, quantum gravity effects are manifested in the quantum Schwarzschild black hole (BH) spacetime as various correction terms of the form
r−n , where n denotes a positive real number [4].The gravitational effects attributed to these quantum corrections are experienced by test particles located near the BHs. The RGI Schwarzschild BH has 2, 1, or 0 horizons, provided the ADM mass of the BH is greater, equal to, or less than a critical mass, respectively. Within the semi-classical approach to quantum gravity, the Hawking evaporation of a BH leads to the formation of a naked singularity, whereas RGI effects prevent the complete evaporation of a BH when the ADM mass of the BH reaches the critical mass, thus leaving behind a cold remnant. It should be noted that the free parameters of the RGI BH have been constrained recently using data from various astronomical sources involving gravitational lensing by the BHs at the centers of the M87 and Milky Way galaxies [5, 6]. The bound orbits and epicyclic frequencies of massive particles around RGI BHs have been recently reported in Ref. [7].
Astrophysical BHs are surrounded by plasma and magnetic fields. This has been recently confirmed based on astronomical observations of the M87 galactic center by the Event Horizon Telescope (EHT) collaboration using polarized synchrotron radiation probes [8, 9]. In particular, the observations reveal the presence of hot plasma with a temperature on the order of ten billion kelvin, along with a magnetic field with a strength on the order of 30 G. The M87 central BH is also an active BH candidate surrounded by an accretion disk, and the mass accretion rate is approximately
(3−20)×10−4 solar mass per year. Theoretical models depicting the dynamics of non-magnetized plasma near BHs with a spherical or axial symmetry have been developed by Perlick et al. [10–14]. These authors not only studied gravitational lensing and BH shadows in the presence of plasma but also investigated plasma accretion in the radial infall of the BH. In the literature, plasma effects on the shadows of diverse modified gravity BHs have been vigorously investigated for phenomenological reasons [15–23], along with gravitational weak lensing in the presence of plasma [24–34]. Previously, gravitational lensing in the weak and strong field regimes in the vacuum surrounding BHs has been reported by some authors in Refs. [35–48]. More recently, new observations made by the EHT team on the shadow of the Sgr A∗ BH [49] have also spurred interest in further analytical studies on BH and wormhole shadow models [50–56].In this study, we aim to extend the previous phenomenological studies on RGI BHs by analyzing the impact of the presence of plasma, along with the effects of quantum corrections, on the motion of photons near a BH. The outline of this article is as follows: In Sec. II, we outline the equations of motion in a plasma medium. In Sec. III, we investigate the shadow images of an RGI BH. In Sec. IV, we compute the deflection angle of light within the weak limit. In Sec. V, we evaluate the shadow images casted by the infalling gas on the BH. Finally, in Sec. VI, we comment on our results. Note that throughout the paper, we use a system of geometric units wherein
G=1=c . -
The spacetime metric describing a static and spherically symmetric RGI BH is given as follows [4]:
ds2=−f(r)dt2+1f(r)dr2+r2(dθ2+sin2θdϕ2),
(1) where
f(r)=1−2Mr(1+ΩM2r2+ΩγM3r3)−1.
(2) Both the Ω and γ parameters denote the quantum effect and are completely free based on the RGI theory. As mentioned earlier in the Introduction, these parameters have been recently constrained, i.e.,
0.02≤γ≤0.22 and0.165≤Ω≤9.804 , and the constraints can be further refined using future astronomical observations acquired from closer inspections of BHs. However, from a theoretical point of view, Ω should be an extremely small quantity. In the astrophysical observations of supermassive BHs withM≫Ω , quantum effects are therefore extremely small, i.e.,Ω≪1 . However, to identify some interesting effects, in the present study, we assume Ω as a free parameter to extrapolate a stronger quantum effect. In our view, current astronomical observations are not precise enough; for example, a huge uncertainty exists in the BH mass, which is why we cannot constrain such quantum effects to a satisfactory level. The Hamiltonian representing the propagation of photons in the plasma medium has the form [57]H(xα,pα)=12[gαβpαpβ−(n2−1)(pβuβ)2],
(3) where
xα denote the spacetime coordinates,pα anduβ denote the four-momentum and four-velocity vectors of the photon, respectively, and n is the refractive index (n=ω/k , where k is the wave number). The refractive index is expressed as follows [58]:n2=1−ω2pω2,
(4) in terms of the plasma frequency
ω2p(xα)=4πe2N(xα)/me (e andme are the electron charge and mass, respectively, and N is the number density of the electrons); the photon frequencyω(xα) is defined byω2=(pβuβ)2 withω(r)=ω0√f(r),ω0=const.
(5) The lapse function is defined such that
f(r)→1 asr→∞ andω(∞)=ω0=−pt, which represents the energy of the photon at spatial infinity [13]. Moreover, the plasma frequency must be sufficiently smaller than the photon frequency(ω2p≪ω2) , thereby allowing the BH shadow to be differentiated from that in the vacuum case. In other words, the natural frequency of oscillation of electrons within the plasma is much smaller than the frequency of light passing through the medium. Using Eqs. (3) and (4), the Hamiltonian for light rays in the plasma medium has the following form:H=12[gαβpαpβ+ω2p].
(6) The components of the four-velocity vector for photons in the equatorial plane
(θ=π/2,pθ=0) are given by˙t≡dtdλ=−ptf(r),
(7) ˙r≡drdλ=prf(r),
(8) ˙ϕ≡dϕdλ=pϕr2,
(9) where we use the following relationship:
˙xα=∂H/∂pα . From Eqs. (8) and (9), we obtain a governing equation for the phase trajectory of lightdrdϕ=f(r)r2prpϕ.
(10) Using the constraint
H=0 , we can rewrite the above equation as [13]drdϕ=√r2f(r)√h2(r)ω20p2ϕ−1,
(11) where we define [13]
h2(r)=r2[1f(r)−ω2p(r)ω20].
(12) The radius of a circular orbit of light, particularly the one that forms a photon sphere of radius
rp , is determined by solving the following equation [13]:d(h2(r))dr|r=rp=0.
(13) After substituting Eq. (12) into (13), one can write the algebraic equation for
rp in the presence of a plasma medium asr2p−3rpM(rp−2M)2−Ω(2γM4+2M4+M3rp)rp(2M−rp)3=[ω2p(rp)ω20+rω′p(rp)ωp(rp)ω20],
(14) where prime denotes the derivative with respect to radial coordinate r. Evidently, the roots of Eq. (14) cannot be obtained analytically for most choices of
ωp(r) ; however, we will consider a few simplified cases below. -
For the special case of a homogeneous plasma medium with a constant plasma frequency throughout the medium, i.e.
ω2p=const. , Eq. (14) can be solved numerically and is depicted in Fig. 1. The displayed figures indicate that the photon radius decreases with increasing parameters γ and Ω. -
Further, we explore photon spheres in the presence of an inhomogeneous plasma medium, where the plasma frequency must satisfy a simple power-law of the form [59, 60]
ω2p(r)=z0rq,
(15) where
z0 and q are free parameters. To analyze the primary features of the power-law model, we restrict ourselves to the following case:q=1 andz0 as a constant [59]. Using Eqs. (14) and (15), we obtain the radius of the photon sphere based on a numerical scheme for the inhomogeneous plasma medium, as illustrated in Fig. 2. One can observe that this profile of the photon radius is approximately similar to that displayed in Fig. 1. This suggests that testing and distinguishing the homogeneous plasma from inhomogeneous plasma around BHs based on their shadows may be quite challenging. -
In this section, we investigate the radius of the shadow of an RGI Schwarzschild spacetime metric in a plasma medium. The angular radius,
αsh , of the BH shadow is defined by a geometric approach, which results in the following [13, 61]:sin2αsh=h2(rp)h2(ro)=r2p[1f(rp)−ω2p(rp)ω20]r2o[1f(ro)−ω2p(ro)ω20],
(16) where
ro andrp represent the locations of the observer and the photon sphere, respectively. Notably, if the observer is located at a sufficiently large distance from the BH, one can approximate the radius of the BH shadow using Eq. (16), as follows [13]:Rsh≃rosinαsh=√r2p[1f(rp)−ω2p(rp)ω20];
this is based on the fact that
h(r)→r , which follows from Eq. (12), at spatial infinity for both the models of plasma. For vacuumωp(r)≡0 , we recover the radius of the Schwarzschild BH shadow,Rsh=3√3M , whenrp=3M . The radius of the BH shadow is depicted for different parameters in Fig. 3 for a homogeneous plasma medium with a fixed plasma frequency, and Fig. 4 presents the case for a power-law model with a plasma frequency ofω2p(r)=z0/r . Both sets of figures illustrate that the shadow radius decreases much more steeply with Ω but remains almost constant with the variation in γ. -
For simplicty, to perform weak lensing, we focus on a weak-field approximation given by
gαβ=ηαβ+hαβ ,
(17) where
ηαβ andhαβ represent a Minkowski spacetime geometry and a minute perturbation in the background of the Minkowski spacetime, respectively. For the foregoing, the general properties are given as follows [62]:ηαβ=diag(−1,1,1,1),hαβ≪1,hαβ→0underxi→∞ ,gαβ=ηαβ−hαβ, hαβ=hαβ.
(18) We now explore the effect of the plasma in the environment surrounding the BH on the gravitational deflection angle. For the plasma medium, the basic equation for the deflection angle is written as [24, 62]
ˆαk=12∫∞−∞(h33+h00ω2−KeN(xi)ω2−ω2e),kdz ,
(19) with
N(xi) representing the concentration of charged particles in the plasma medium and ω andωe respectively referring to the photon and plasma frequencies. Note thatKe=4πe2/me defines the constant value of a plasma particle. The above equation [62] can be rewritten as follows:ˆαb=12∫∞−∞br(dh33dr+11−ω2e/ω2dh00dr−Keω2−ω2e dNdr)dz.
(20) Note that in further calculations,
ˆab can assume both negative and positive values, depending on the light position and on its motion toward the compact object or away from the compact object.On expanding, for greater r, the BH metric can be expressed as
ds2=ds20+(2Mr−2M4γΩr4−2M3Ωr3)dt2+(2Mr−2M4γΩr4−2M3Ωr3)dr2,
(21) where
ds20 describes the line element for the Minkowski spacetime, which is defined byds20=−dt2+dr2+r2(dθ2+sin2θdϕ2).
(22) For further calculations associated with Eq. (20), we rewrite the
haβ components in Cartesian coordinates for simplicity:h00=(2Mr−2M4γΩr4−2M3Ωr3),hik=(2Mr−2M4γΩr4−2M3Ωr3)nink,h33=(2Mr−2M4γΩr4−2M3Ωr3)cos2χ.
(23) Here,
cos2χ=z2/(b2+z2) andr2=b2+z2 are new notations. In creating these notations, we write the first derivatives ofh00 andh33 as follows:dh33dr=2z2(6γM4Ω+5M5rΩ−3Mr3)r7,dh00dr=8γM4Ωr5+6M3Ωr4−2Mr2.
(24) Further, the deflection angle considered can be obtained based on the following parts [34]:
ˆαb=ˆα1+ˆα2+ˆα3,
(25) with
ˆα1=12∫∞−∞brdh33drdz,ˆα2=12∫∞−∞br(11−ω2e/ω2dh00dr)dz,ˆα3=12∫∞−∞br(−Keω2−ω2e dNdr)dz .
(26) From the above equations,
ˆα1 ,ˆα2 , andˆα3 respectively refer to the contributions of gravity and the homogeneous and inhomogeneous plasma medium to the deflection angle. To consider the impact of the plasma medium on the deflection angle, we exploit Eq. (25). -
Let us consider the impact of a uniform plasma medium on the deflection angle given by Eq. (25) [26]:
ˆαuni=ˆαuni1+ˆαuni2+ˆαuni3,
(27) where
ˆauni1 andˆauni2 respectively refer to the contributions of uniform plasma, andˆauni3=0 represents the contribution of the uniform plasma distribution. Following Eqs. (24), (25), and (27), the deflection angles of photons around static BHs with spherical symmetry in the presence of RGI gravity and uniform plasma can be defined byˆαuni=2Mb−3πγM4Ω8b4−4M3Ω3b3+(2Mb−3πγM4Ω2b4−4M3Ωb3)11−(w2p/w2).
(28) In Fig 5, the deflection angle of light travelling in a uniform density plasma medium about an RGI BH is plotted against the impact parameter and RGI theory parameters. The figures demonstrate that to achieve a large deflection angle of light in the presence of strong gravity, the RGI parameters must be extremely small.
-
Further, we investigate the impact of non-uniform plasma on the deflection angle of photons around a BH in RGI gravity. Note that in a previous study, a Singular Isothermal Sphere (SIS) has been proposed to describe the non-uniform plasma medium distribution (see for example [62]). Thus, the plasma concentration in the SIS medium is defined by [24, 62]
N(r)=ρ(r)kmp,
(29) with the following plasma density
ρ(r)=σ2ν2πr2 , whereσν refers to the dispersion velocity. Thus, Eq. (25) can be written asˆαSIS=ˆα1SIS+ˆα2SIS+ˆα3SIS,
(30) where the first two terms
ˆa1SIS andˆa2SIS reflect the contributions of the gravity and plasma effects to the deflection angle, respectively, whereas the last termˆa3SIS reflects the contribution of the density gradient of the plasma medium. Relying on Eqs. (24), (25), and (30), the deflection angles of the photons in non-uniform plasma assume the following form:ˆαSIS=4Mb−5γM6Ωw2cb6w2−64M5Ωw2c5πb5w2−15πγM4Ω8b4+16M3w2c3πb3w2−16M3Ω3b3+2M2w2cb2w2,
(31) where we define
w2c=σ2νKe2kmpR2s .
(32) In Fig. 6, the deflection angle of light in an SIS with inhomogeneous density plasma about an RGI BH is plotted against the impact parameter and RGI theory parameters. The figures demonstrate that to achieve a large deflection angle of light in the presence of strong gravity, the RGI parameters must be extremely small.
Figure 6. (color online) SIS: dependence of the deflection angle on the impact parameter b and RGI theory parameters
{γ,Ω} .We now consider a comparison of the deflection angles between uniform and non-uniform plasma media. The behavior of the deflection angle for both the uniform and non-uniform plasma medium cases is depicted in Fig. 7. As can be observed from Fig. 7, the behavior evidently indicates that the deflection angle of the photon beam around the BH has larger values compared to that in the non-uniform case. This is because the photon beam is scattered under the influence of the non-uniform plasma. Thus, the deflection angle in the uniform plasma medium is greater than that in the non-uniform plasma medium.
-
Further, we aim to evalute the effect of the plasma medium on the optical images of the RGI BH. To achieve this goal, we use the backward raytracing technique [50, 63–68]. In this technique, for the given image, one has to compute the specific intensity,
Iν0 , observed at some distance from the BH [64]Iobs(νobs,X,Y)=∫γg3j(νe)dlprop,
(33) where
g=νobs/νe is the redshift factor, andνe denotes the photon frequency measured in the rest-frame of the emitter. For the total flux, we have [64, 69]Fobs(X,Y)=∫γIobs(νobs,X,Y)dνobs.
(34) Notably, for radiating gas (in free fall), one has to use the four-velocity components [64]
ute=1f(r),ure=−√1−f(r),uθe=uϕe=0.
(35) In addition, we must consider the relation between the radial and time components of the photon four-velocity vector, given by
kr=±kt√1f(r)(1f(r)−b2r2−ω2pω20),
(36) which is derived from
H=0 . Thus, we determine the presence of an extra term owing to the plasma medium compared to the corresponding form in [64]. The signs+(−) indicate cases when the photon approaches or recedes from the BH. Furthermore, the impact parameter b is also modified in the plasma medium:b=r√1f(r)−ω2p(r)ω20.
(37) For the redshift function
g , we use [64]g=kαuαokβuβe.
(38) A radial profile is assumed with the
1/r2 law and the specific emissivityj(νe)∝δ(νe−ν⋆)r2,
(39) where δ is the Dirac delta function. For the proper length, we have
dlprop=kαuαedλ=−ktg|kr|dr.
(40) For the total flux, we can therefore write [64]
Fobs(X,Y)∝−∫γg3ktr2krdr.
(41) Let us finally consider one additional model consisting of radiating gas at rest (static gas model). In this case, for the redshift factor, we have
g=f(r)1/2 . Assuming a radial profile with1/r2 , the proper length can be written as follows:dlprop=√f(r)−1dr2+r2dϕ2.
(42) For the specific intensity observed, we have
Iobs(νobs)=∫γf(r)3/2r2√f(r)−1+r2(dϕdr)2dr,
(43) where
drdϕ=±r√f(r)(h(r)2b2−1)
(44) is determined based on the equation of motion obtained for the light ray moving along the equatorial plane. As the function
h(r) and the impact parameter encode the plasma effect, the intensity should change. In addition to this, the RGI parameters γ and Ω also exert some effects. In Figs. 8, 9, and 10, we illustrate the shadow images of the BH with/without a plasma medium.Figure 8. (color online) Shadow images with an infalling and static gas model in the RGI BH without a plasma medium with
Ω=0.1 andγ=0.5 . -
In this study, we investigate the effect of plasma on the optical properties of an RGI Schwarzschild BH and discover the presence of relevant significant effects on the photon sphere radius of the BH shadow and gravitational weak lensing. Here, the orbits of photons are obtained using the general method and are presented in Figs. 1 and 2. The plots indicate that the photon orbits come closer to the central object with an increase in γ or Ω. We also discuss the shadow cast by the BH in the presence of RGI gravity and notice the effects of the model parameters γ and Ω on the BH shadow. The results reveal that with increasing values of the parameters γ and Ω, the radius of the shadow of the BH decreases, and this can be visually confirmed from Figs. 3 and 4. Additionally, we consider the deflection angle of light rays around a BH in RGI gravity in the presence of plasma (different distributions: uniform and non-uniform cases) and analytically obtain an equation for the deflection angle in the weak-field regime. In this part of our results, we discuss and analyze the effect of the uniform and non-uniform plasma on the deflection angle for different values of the parameters γ and Ω.
We also demonstrate the effect mechanism of the plasma medium on the motion of light owing to infalling gas or gas at rest and obtain corresponding optical images. Although the final image is found to depend on the specific model, the photon ring is a universal characteristic and therefore interesting from an observational point of view. From Figs. 8–10,it can be observed that the intensity of the photon ring predicted by these models is significantly higher for the static gas model compared to the infalling model. As another observation, one can infer that in the plasma medium, the intensity is lower compared to that in the absence of plasma.
-
F.A. acknowledges support from the Inha University, Tashkent.
Quantum effects on the black hole shadow and deflection angle in the presence of plasma
- Received Date: 2022-11-26
- Available Online: 2023-03-15
Abstract: In this study, the optical properties of a renormalization group improved (RGI) Schwarzschild black hole (BH) are investigated in a plasma medium. Beginning with the equations of motion in a plasma medium, we aim to present the modifications in the shadow radius of the RGI BH. To this end, we compute the deflection angle of light in the weak gravity regime for uniform and non-uniform plasma media. Importantly, owing to the plasma media, we discover that the equations of motion for light obtained from the radiating and infalling/rest gas have to be modified. This, in turn, changes and modifies the expression for the intensity observed far away from the BH. Finally, we obtain the shadow images for the RGI BH for different plasma models. Although quantum effects change the background geometry, such effects are minimal, and practically detecting these effects using the current technology based on supermassive BH shadows is impossible. The parameter Ω encodes the quantum effects, and in principle, one expects such quantum effects to play significant roles only for very small BHs. However, the effects of plasma media can play an important role in the optical appearance of BHs, as they affect and modify the equations of motion.