Processing math: 85%

Photon motion and weak gravitational lensing in black-bounce spacetime

  • The effect of spacetime curvature on photon motion may offer an opportunity to propose new tests on gravity theories. In this study, we investigate and focus on the massless (photon) particle motion around black-bounce gravity. We analyze the horizon structure around a gravitational compact object described by black-bounce spacetime. The photon motion and the effect of gravitational weak lensing in vacuum and plasma are discussed, and the shadow radius of the compact object is also studied in black-bounce spacetime. Additionally, the magnification of the image is studied using the deflection angle of light rays.
  • The effect of a curved spacetime on wave propagation is an important subject in modern relativistic astrophysics. The first ever observation of gravitational waves by the LIGO-Virgo collaboration [1, 2] was the starting point for a new era of gravitational wave observation. Later, the observation of images of M87 [3, 4] and SgrA* [5] motivated specialists to use these observations to test gravity theories in the strong field regime.

    Although well tested in the weak [6] and strong field regimes [2, 3], Einstein's general relativity still encounters several fundamental problems. One is the existence of the physical (curvature) singularity at the origin of most solutions within the theory. There is a belief that an undiscovered new quantum gravity theory will resolve this issue. One of the main aims of introducing alternative theories of gravity and modified versions of general relativity is to move in the direction of constructing a quantum gravity theory.

    Owing to the strong gravitational field, the spacetime around massive objects is curved. The propagation of an electromagnetic wave is affected by this curvature, and an observer at infinity can detect light deflection due to gravity. However, because of the strong gravitational field around black holes, the part of the light rays propagating near a black hole will be captured by the central object. Consequently, the observer at infinity will detect a black spot on the celestial plane near the position of the black hole. The original idea of the possible observation of this black spot, usually referred to as a black hole shadow, was proposed by Synge [7] and later developed by Luminet [8] and Bardeen [9]. The shadow of a black hole for different solutions within general relativity and/or modified/alternative theories of gravity has subsequently been investigated by the various authors [1038].

    The deflection of light, also known as gravitational lensing, around a gravitating object can be used to either study the light source or test the gravitational field around the lensing object. The effect of gravitational lensing was first used to test general relativity by Eddington in 1919 during the solar eclipse. Later, the gravitational lensing around compact objects within different gravity models was investigated by various authors (see, for example, [3949]). The existence of a plasma environment is one of the main markers of most black holes and compact objects. The light propagating through the plasma is also affected by the latter. The effects of both uniform and nonuniform plasma on the optical properties of black holes were investigated in Refs. [5058]). Moreover, using gravitational lensing, one may study the effects of different plasma configurations on the light deflection angle (see, for example, [5980]).

    A regularization method was proposed by Simpson and Visser [81], in which the expression 1/r in the Schwarzschild metric is replaced with 1/r2+a2. As a result, one may obtain a globally regular family of metrics, including black-bounce geometry and a traversable wormhole. A particularly interesting feature of these solutions is that they may be considered a one-parameter extension of the Schwarzschild solution. The generalization of Reissner-Nordström spacetime using this technique has been reported in Ref. [82]. Later, the authors of Ref. [83] extended them and proposed a number of new solutions for black-bounce spacetime. Recently, the effect of retrolensing by two-photon spheres within black-bounce spacetime was studied in detail in Ref. [84]. In this paper, we extend the analysis of this spacetime by studying massless test particle motion around a black-bounce compact object. Moreover, we also study the effect of the gravitational lensing and shadow of the compact object described by black-bounce spacetime.

    The paper is arranged as follows. First, we review the spacetime metric in Sec. II. We study the photon motion around the black bounce compact object in Sec. III. Then, the shadow of the black-bounce object is studied in Sec. IV in the presence of plasma. Gravitational lensing in the presence of plasma in black bounce spacetime is then investigated in Sec. V. We conclude our study in Sec. VI. Throughout, we use the unit system G=c=1.

    Non-rotating spacetime in black bounce gravity has the following linear element in spherical coordinates [83, 84]:

    ds2=f(r)dt2+dr2f(r)+γ(r)(dθ2+sin2θdϕ2) ,

    (1)

    where the metric functions f(r) and γ(r) are defined as

    f(r)=12Mr2(a2+r2)3/2 ,

    (2)

    γ(r)=r2+a2 ,

    (3)

    where M is the mass of the compact object, and a is a nonnegative constant [83]. To understand the nature of the gravitational potential of the source, we may introduce the effective mass as Meff=M(r/r2+a2)3.

    Spacetime metric (1) describes a regular black hole when the condition 0<a/m<43/9 is satisfied. To obtain more information about the spacetime structure, we can analyze the function f(r). The dependence of f(r) on radial coordinates is shown in Fig. 1. From the plot, we can easily see that there is an extreme value of a/M when the graph has one r/M-intercept. The extreme value of the parameter aext=43/9. The values of parameter a that are less than its extreme value correspond to a black hole with two horizons (Fig. 2 shows the existence of two roots for a set of parameters a values, which corresponds to the Cauchy horizon (smaller root) and event horizon (larger root)), whereas the values of parameter a that are greater than its extreme value correspond to a naked singularity. The dependence of the radius of the horizon corresponding to the condition f(r)=0 is represented in Fig. 2. From the plot, we can see that with an increase in parameter a in black-bounce spacetime, the outer horizon decreases until the value M is reached.

    Figure 1

    Figure 1.  (color online) Radial dependence of the lapse function for different values of the parameter a/M.

    Figure 2

    Figure 2.  (color online) Horizon radius of a black hole as a function of a/M.

    Consider the motion of photons around a gravitational object in the framework of black-bounce gravity. To describe the dynamics, we can use the Hamilton-Jacobi equation in the following form [85]:

    gμνdSdxμdSdxν=0,

    (4)

    where gμν is the metric tensor of spacetime metric (1), and S is the action. Using the separation of variables method and the symmetry of spacetime metric (1), the action S can be expressed as the following form:

    S=Et+Lφ+Sr(r)+Sθ(θ),

    (5)

    where E and L can be described as the energy and angular momentum of a massless particle, respectively. The dynamics of particles with zero rest mass are defined by null geodesics of spacetime. Using Eqs. (4) and (5), we can obtain the equations of motion of massless particles at the equatorial plane (θ=0.5π, ˙θ=0) as

    dtdλ=Ef(r),

    (6)

    drdλ=±R(r),

    (7)

    dϕdλ=Lγ(r),

    (8)

    with

    R(r)=E2f(r)γ(r)L2,

    (9)

    where λ is the affine parameter, and we consider θ=0.5π. Accordingly, the effective potential for the radial motion of photons takes the following form:

    Veff=f(r)L2γ(r).

    (10)

    The study of circular orbits r reflects significant interest while considering the photon (massless particle) motion. Photon circular orbits rph are determined using the following conditions:

    Veff=0.

    (11)

    Using the above conditions, we can obtain the equation to determine the radius of circular photon orbits rph as

    (a2+r2)3/2+2Ma23Mr2=0.

    (12)

    We can obtain the radius of the photon sphere (rph) by solving Eq. (12) with respect to the radial coordinate,

    rph=a2(10M2D+1)+D+9M4D+3M2 ,

    (13)

    where D is defined as

    D=(25a4M2245a2M4+552a3M25a24M2+27M6)1/3.

    (14)

    The dependence of the photon orbit rph on the parameter a is shown in Fig. 3. From the plot, we can conclude that with an increase in parameter a, the radius of the photon orbit rph decreases.

    Figure 3

    Figure 3.  (color online) Dependence of the radius of photon circular orbits rph/M on parameter a.

    The detailed study of photon motion around a black hole leads to a phenomenon known as a black hole shadow [16]. The shadow cast by a black hole is a result of photon capture by the central object, and the observer will see a black spot on the bright background [16]. In this section, we explore the radius (Rsh) of the shadow of a black hole in the presence of a plasma medium [55]. The boundary of the shadow is fully defined by the photon (massless particle) trajectory governed by the equation of motion in the plasma medium. To analyze the shadow of the black hole in black-bounce spacetime in a plasma medium, we should first study the photon equation in a plasma environment.

    The Hamiltonian for photons around a black hole in the presence of an electron plasma medium can be expressed as [55, 85]

    H=12[gμνpμpν(n21)p2tf(r)]=0,

    (15)

    with

    n2=1ω2eω2(xi),

    (16)

    where ω(xi) and ωe are the photon and plasma frequencies, respectively, and n is the refractive index of the medium. We denote the value of the photon frequency observed at infinity as ω0:

    ω()=ω0=pt.

    (17)

    The components of the four-velocity (˙t, ˙ϕ, ˙r) for the photons (massless particle) in the equatorial plane (θ=0.5π, pθ=0) are given by

    dtdλ=ptf(r),

    (18)

    drdλ=prf(r),

    (19)

    dϕdλ=pϕr2+a2,

    (20)

    where we consider the relationship dxμ/dλ=H/pμ. From Eqs. (19) and (20), we can obtain an expression for the phase trajectory of light (photon) rays,

    drdϕ=grrprgϕϕpϕ

    (21)

    Using the constraint H=0, we can define the above equation as [50]

    drdϕ=grrgϕϕh2(r)ω20p2ϕ1,

    (22)

    where the following definition is introduced:

    h2(r)gttgϕϕ1gϕϕω2e(r)ω20.

    (23)

    The radius of a circular orbit of photons rph, particularly the one that forms the photon sphere of radius rph, is determined by the following equation [50]:

    d(h2(r))dr|r=rph=0.

    (24)

    Substituting Eq. (24) into Eq. (23), we can define an algebraic expression for rph in the plasma environment as

    (a2+r2)f(r)f(r)+2f(r)ωe(r)ω20[(a2+r2)ωe(r)+rωe(r)]=2r,

    (25)

    where the prime stands for the derivative with respect to r. In the case of homogeneous plasma (ω2e(r)=const.), Eq. (25) can be simplified to

    (a2+r2)f(r)f(r)+2rf(r)ω2eω20=2r.

    (26)

    As mentioned above, the solution to this equation with respect to r leads to the radius of the photon sphere rph of the black hole in a uniform plasma environment.

    Figure 4 gives information on the dependence of the radius of the photon sphere rph on the black-bounce parameter a and plasma parameters. It is clear from the plots that the radius of the photon sphere rph decreases with an increase in parameter a. Note that the presence of plasma causes an increase in the value of the photon sphere radius, as shown in Fig. 4.

    Figure 4

    Figure 4.  (color online) Dependence of the radius of the photon sphererph/M on plasma frequency (top) and black-bounce parameter a (bottom).

    In this section, we study the radius (Rsh) of the shadow of the black hole in the presence of uniform plasma. The angular radius αsh of the black hole shadow is defined by [7, 50]

    sin2αsh=h2(rsh)h2(ro)=(a2+r2ph)(1f(rph)ω2e(rph)ω20)(a2+r2o)(1f(ro)ω2e(ro)ω20),

    (27)

    where ro and rph are the radial positions of the observer and photon sphere (introduced in the previous section), respectively. If the observer is located at a sufficiently large distance from the black hole, we can approximate the radius of the black hole shadow (for the case of uniform plasma) using Eq. (27) [50]:

    Rshrosinαsh=(a2+r2ph)(1f(rph)ω2eω20).

    (28)

    For the case of no plasma in the vicinity of the black-bounce black hole, the shadow radius can be represented by

    Rshrosinαsh=a2+r2phf(rph).

    (29)

    Figure 5 describes the dependence of the radius of black hole shadow on parameter a and plasma parameters. In the plots, we define Rsh in terms of Meff because it is known from EHT data that the masses of M87* and Sgr A* are given by the gravitational potential. Based on these data, we can conclude that the shadow radius is diminished due to the influence of both plasma and the black-bounce parameter a.

    Figure 5

    Figure 5.  (color online) Dependence of the radius of the black hole shadow on parameter a (top) and plasma frequency (bottom).

    Now, we consider the rough assumption that the supermassive black holes M87* and Sgr A* are the black-bounce black holes that we are studying. Then, we can theoretically obtain the constraints from the EHT observation results.

    The angular diameter of the image of the supermassive black hole M87* is θ=42±3 μas at the 1σ confidence level, and the mass of M87* and the distance from the solar system are M=6.5×109M and D= 16.8 Mpc, respectively [3, 4]. The same data for Sgr A* are as follows: θ=48.7±7 μas, M4×106M, and D8 kpc [5]. The shadow radius caused by the black hole per unit mass can be expressed as the following equation:

    dsh=DθMeff,

    (30)

    where Meff=M(r/r2+a2)3 from black-bounce spacetime. Here, we must use the effective mass because the masses of M87 and Sgr A* are given by the gravitational potential.

    Using Eq. (30), we can easily calculate the diameter of the shadow images. The diameter of the black hole shadow is dsh=(11±1.5)Meff for M87* and dsh=(9.5±1.4)Meff for Sgr A*. From these data, we can now obtain the constraints on the black-bounce and plasma parameters for the supermassive black holes Sgr A* and M87*. The results are presented in Fig. 6 for the case of no plasma around the black hole and in Fig. 7 for the black hole in plasma.

    Figure 6

    Figure 6.  (color online) Constraints on the black-bounce parameter a with no plasma around the black hole. The brown and black dashed lines describe the limits of the measured shadow radii of M87* and Sgr A*, respectively.

    Figure 7

    Figure 7.  (color online) Shadow diameter of the black-bounce black hole as a function of the parameter a and plasma parameter. The black curve corresponds to dsh=9.5Meff, which is the lower limit for M87*. The blue curve corresponds to dsh=10.9Meff, which is the upper limit for Sgr A*.

    Now, we can move on to another part of our study: gravitaional lensing in spacetime metric (1) in a weak-field case explicated via [59]

    gαβ=ηαβ+ξαβ ,

    (31)

    where ηαβ is the general metric of Minkowski spacetime, and ξαβ is a perturbation of the flat spacetime metric, which can be described by the spacetime metric theory of gravity and is written as [59]

    ηαβ=diag(1,1,1,1),ξαβ1,ξαβ0underxi ,gαβ=ηαβξαβ,   ξαβ=ξαβ .

    (32)

    Here, we can study the effects of the plasma medium on the deflection angle ˆαk in the gravitational weak field of the black hole. The general expression for the deflection angle in plasma is [59, 64]

    ˆαi=12(ξ33+ξ00ω2KeN(xi)ω2ω2e),idz,  i=1,2

    (33)

    where N(xi) is the number density of particles in a plasma medium around a compact object, Ke=4πe2/me is constant (e and me are the electron charge and mass), and ω and ωe are the photon (massless particle) and plasma frequencies, respectively [59]. Using Eqs. (32) and (33), we can rewrite Eq. (33) for the deflection angle in the following form [59]:

    ˆαb=12br(dξ33dr+11ω2e/ω2dξ00drKeω2ω2e dNdr)dz,

    (34)

    where b is the impact parameter of light rays. It is useful to note that the values of ˆαb can be both negative and positive [59].

    In the weak-field case for far distances from the black hole, the black-bounce spacetime metric (1) can be written as

    ds2=ds20+2Mr2(a2+r2)3/2dt2+2Mr2(a2+r2)3/2dr2+a2(dθ2+sin2θdϕ2),

    (35)

    where ds20=dt2+dr2+r2(dθ2+sin2θdϕ2) is the metric element in the Minkowski spacetime metric, and we can use the notation Rs=2M for further calculations.

    To analyze the deflection angle (ˆαb) of light rays around the black hole in the plasma medium using Eq. (34), we can rewrite the required components (ξ00 and ξ33) in Cartesian coordinates as

    ξ00=Rsr2(a2+r2)3/2,ξ33=Rsr2(a2+r2)3/2cos2χ+a2r2sin4χ ,

    (36)

    where cos2χ=z2/(b2+z2) and r2=b2+z2 are introduced, for example, in [64]. Here, z is the coordinate aligned with the optical axis. We can obtain the derivative of ξ00 and ξ33 with respect to the radial coordinate as

    dξ33dr=3rRsz2(a2+r2)5/22a2(r44r2z2+3z4)r7,dξ00dr=Rs(2a2rr3)(a2+r2)5/2.

    (37)

    For simplicity, the expression for the deflection angle ˆαb can be expanded as [63]

    ˆαb=ˆα1+ˆα2+ˆα3,

    (38)

    with

    ˆα1=12brdξ33drdz,ˆα2=12br(11ω2e/ω2dξ00dr)dz,ˆα3=12br(Keω2ω2e dNdr)dz,

    (39)

    where the notations ˆα1, ˆα2, and ˆα3 represent the contributions to the deflection angle due to gravity, uniform plasma, and non-uniform plasma, respectively. In this study, we use Eqs. (38) and (39) to investigate the effect of plasma on the deflection angle in gravitational weak lensing. In the following sections, we consider each case in detail.

    Here, we can test the deflection angle of photons (light rays) around the black hole in the presence of a uniform plasma environment using expression (38) written in the following form:

    ˆαuni=ˆαuni1+ˆαuni2+ˆαuni3,

    (40)

    We can claim that the third term of this equation will vanish due to the uniform distribution of the plasma environment. Using Eqs. (38), (37), and (40), we can obtain an equation for the deflection angle of light rays (photon) around the black hole in black-bounce spacetime surrounded by uniform plasma:

    ˆαuni=bRsa2+b2+3πa216b2bRs(a2b2)(a2+b2)2ω2ω2ω2e.

    (41)

    In this equation, the second term appears owing to the property of black-bounce spacetime itself.

    Using expression (41), we can plot the deflection angle against the impact b, plasma frequency, and a. Fig. 8 shows the dependence of the deflection angle ˆαuni of light rays (photon) around the black hole on the impact parameter b for different values of the parameter a and plasma parameter ω2e/ω2. The dependence of the deflection angle ˆαuni plasma frequency and parameter a for fixed values of the impact parameter is shown in Fig. 9. The graphs of the defection angle ˆαuni on the impact parameter b and plasma parameter are reliable, as expected. The presence of uniform plasma leads to an increase in the value of the deflection angle of light rays, whereas the angle of deflection decreases dramatically with increasing impact parameter b. However, interestingly, the deflection angle remains almost unchanged with the influence of the black-bounce parameter a. There is only a slight increase in the value of the deflection angle with an increase in parameter a, as shown in Figs. 8 and 9.

    Figure 8

    Figure 8.  (color online) Dependence of the deflection angle ˆαuni on the impact parameter b for different values of the parameter a/M (left panel) and plasma medium ω2e/ω2 (right panel).

    Figure 9

    Figure 9.  (color online) Dependence of the deflection angle ˆαuni on the parameter a/M (left panel) and plasma parameter ω2e/ω2 (right panel) for a fixed value of the impact parameter b/M=30.

    Here, we investigate the deflection angle ˆαuni of light rays (photons) around the black hole in the presence of non-uniform plasma and in the framework of black-bounce spacetime with a nonzero parameter a. We can consider the singular isothermal sphere (SIS) medium as the non-uniform plasma distribution [59, 64]. The plasma number density of the SIS environment is written in the following form [59, 63]:

    N(r)=ρ(r)kmp ,

    (42)

    where mp is the proton mass, and k is a non-dimensional coefficient responsible for the dark matter contribution in the SIS model and with

    ρ(r)=σ2ν2πr2 ,

    (43)

    where σν and ρ(r) are the velocity of dispersion and the density of plasma around the compact object, respectively [59]. Here, Eq. (38) can be written in the non-uniform plasma case as [63]

    ˆαSIS=ˆαSIS1+ˆαSIS2+ˆαSIS3 ,

    (44)

    where ˆαSIS1, ˆαSIS2, and ˆαSIS3 represent the deflection angle due to the gravity of the black hole, due to plasma, and due to the density of the non-uniform plasma medium around the compact object in black-bounce spacetime, respectively. From Eqs. (38), (37), and (44), we can obtain an equation for the deflection angle of light rays (photons) around the black hole in black-bounce surrounded by non-uniform plasma:

    ˆαSIS=3πa216b2+2b3Rs(a2+b2)2+R2sω2c2b2ω2+2bR3sω2c(5b27a2)15πω2(a2+b2)3,

    (45)

    where the new notation for the non-uniform case is

    ω2c=σ2νKe2kmpR2s.

    (46)

    In Fig. 10, we show the plasma effect of non-uniform plasma and the parameter a on the deflection angle of photons around the black hole in black-bounce spacetime. The dependence of the deflection angle ˆαSIS on non-uniform plasma and the parameter a for fixed values of the impact parameter b and the non-uniform plasma parameter is clearly shown in Fig. 11. The plots reveal that the deflection angle increases gradually with increasing plasma parameter, whereas there is a small change in the value of the deflection angle due to the influence of parameter a.

    Figure 10

    Figure 10.  (color online) Dependence of the deflection angle ˆαSIS on the impact parameter b for different values of the parameter a/M (left panel) and non-uniform plasma parameter (right panel).

    Figure 11

    Figure 11.  (color online) Dependence of the deflection angle ˆαSIS on the parameter a/M (left panel) and non-uniform plasma parameter (right panel) for a fixed value of the impact parameter b/M=30.

    We can easily compare the two cases for analysis (uniform and non-uniform) with the help of Fig. 12. It is evident that the deflection angle of light rays (photons) in the presence of a uniform plasma medium is larger than that of the non-uniform plasma case.

    Figure 12

    Figure 12.  (color online) Dependence of the deflection angle ˆαb on the impact parameter b (left panel) and parameter a/M (right panel).

    Now, we can study the brightness of an image source as a consequence of gravitational lensing in the presence of uniform plasma using a lens equation in the form [64, 76]

    θDs=βDs+ˆαbDds ,

    (47)

    where Ds and Dds are the distances from the distant source to the observer and the lens object, respectively, and θ and β represent the angular positions of the image and source, respectively. We use the following relation between the impact parameter b and angle θ:

    b=Ddθ,

    (48)

    where Dd is the distance from the observer to the lens object. Using these expressions, we can rewrite Eq. (47) as [64, 76]

    β=θDdsDsF(θ)Dd1θ ,

    (49)

    with

    F(θ)=|αb|b=|αb(θ)|Ddθ.

    We can write the expression for Einstein's angle θplE in the presence of uniform plasma around a black-bounce black hole using Eq. (41),

    θplE=DdsDsˆαuni(b).

    (50)

    We numerically determine the angle by solving this equation to analyze the dependence of the magnification on the black hole parameter.

    It is well-known that the general equation to compute the magnification of an image source can be expressed in the following form [59]:

    μΣ=ItotI=k|(θkβ)(dθkdβ)| ,k=1,2,...,s ,

    (51)

    where s denotes the total number of images of the source, and Itot can be explained as the value of the total increased brightness due to multiple images of the source with brightness I.

    From Eq. (51), we can obtain the expression for the magnification of the image source as

    μpltot=μpl++μpl=x2uni+2xunix2uni+4 ,

    (52)

    where

    xuni=β(θplE)uni.

    (53)

    The magnification of the image source is written as

    (μpl+)uni=14(xunix2uni+4+x2uni+4xuni+2) ,

    (54)

    (μpl)uni=14(xunix2uni+4+x2uni+4xuni2) .

    (55)

    Using Eq. (52), we present the dependence of the total magnification on the plasma parameter for different values of parameter a in Fig. 13. From these plots, we can conclude that the total magnification of the images increases owing to the influence of uniform plasma, and there is a slight change in the total magnification with an increase in the parameter a.

    Figure 13

    Figure 13.  (color online) Dependence of the total magnification μtot on the parameter a for different values of the plasma parameter and on the plasma parameter for different values of parameter a.

    In this study, particle orbits around black holes in black bounce gravity are studied in detail. The results obtained can be summarized as follows:

    ● We calculate the horizon radius for various values of the black bounce parameter a. The results show that the radius of the horizon rh decreases in the presence of the parameter a in black-bounce spacetime.

    ● We analyze the photon orbits around a black-bounce compact object. It is shown that orbits shift toward the central object because of black bounce gravity.

    ● We obtain and discuss the observable shadow of the black hole described by the black-bounce spacetime metric in plasma. We reveal that the shadow radius (Rsh) of the black hole decreases with an increase in both the black-bounce parameter a and plasma frequency (Fig. 5). The tendency of the black hole shadow size to decrease with an increase in parameter a can be explained in a similar way to that of photon deflection. The effect of parameter a decreases the gravitational effect of the central object on photons, and we may observe a decrease in the size of the shadow in the black-bounce spacetime metric. The negative effect of parameter a on photon orbits may be useful in obtaining constraints using current and future observation of the black hole shadow with event horizon telescopes.

    ● We also study the effect of weak gravitational lensing around the black hole in black bounce spacetime in the presence of a plasma medium around a compact object. Here, uniform and non-uniform cases are considered and discussed in detail (see. Figs. 812).

    ● For the non-uniform case, we find that the plasma frequency depends on 1/r. The deflection angle around compact objects increases with increasing value of the parameter a. We also compare the uniform and non-uniform cases and find that the deflection angle is larger in uniform plasma than in non-uniform plasma.

    ● Additionally, we consider the magnification of the image using the deflection angle of light rays, which is clearly represented in Fig. 13.

    [1] B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016) doi: 10.1103/PhysRevLett.116.061102
    [2] B. Abbott, R. Abbott, T. Abbott et al., Phys. Rev. Let. 116 (2016) doi: 10.1103/physrevlett.116.241102
    [3] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. 875, L1 (2019), arXiv:1906.11238[astro-ph.GA doi: 10.3847/2041-8213/ab0ec7
    [4] K. Akiyama et al., Astrophys. J. 875, L6 (2019), arXiv:1906.11243[astro-ph.GA doi: 10.3847/2041-8213/ab1141
    [5] K. Akiyama et al., Astrophys. J. Lett 930, L12 (2022) doi: 10.3847/2041-8213/ac6674
    [6] C. M. Will, Living Reviews in Relativity 9, 3 (2006), arXiv:gr-qc/0510072[gr-qc doi: 10.12942/lrr-2006-3
    [7] J. L. Synge, Mon. Not. R. Astron. Soc. 131, 463 (1966) doi: 10.1093/mnras/131.3.463
    [8] J. P. Luminet, Astron.Astrophys. 75, 228 (1979)
    [9] J. M. Bardeen, in Black Holes ((Les Astres Occlus, 1973) pp. 215–239
    [10] H. Falcke, F. Melia, and E. Agol, Astrophys. J. 523, L18 (2000), arXiv:9912263 doi: 10.1086/312423
    [11] C. Bambi and K. Freese, Phys. Rev. D 79, 043002 (2009), arXiv:0812.1328[astro-ph doi: 10.1103/PhysRevD.79.043002
    [12] K. Hioki and K.-I. Maeda, Phys. Rev. D 80, 024042 (2009) doi: 10.1103/PhysRevD.80.024042
    [13] A. Abdujabbarov, F. Atamurotov, Y. Kucukakca et al., Astrophys. Space. Sci. 344, 429 (2013), arXiv:1212.4949[physics.gen-ph doi: 10.1007/s10509-012-1337-6
    [14] L. Amarilla and E. F. Eiroa, Phys. Rev. D 85, 064019 (2012) doi: 10.1103/PhysRevD.85.064019
    [15] A. Abdujabbarov, M. Amir, B. Ahmedov et al., Phys. Rev. D 93, 104004 (2016), arXiv:1604.03809[gr-qc doi: 10.1103/PhysRevD.93.104004
    [16] F. Atamurotov, A. Abdujabbarov, and B. Ahmedov, Phys. Rev. D 88, 064004 (2013) doi: 10.1103/PhysRevD.88.064004
    [17] N. Tsukamoto, Phys. Rev. D. 97, 064021 (2018), arXiv:1708.07427[gr-qc doi: 10.1103/PhysRevD.97.064021
    [18] A. Abdujabbarov, F. Atamurotov, N. Dadhich et al., Eur. Phys. J. C 75, 399 (2015), arXiv:1508.00331[gr-qc doi: 10.1140/epjc/s10052-015-3604-5
    [19] R. C. Pantig, P. K. Yu, E. T. Rodulfo et al., Annals of Physics 436, 168722 (2022), arXiv:2104.04304[gr-qc doi: 10.1016/j.aop.2021.168722
    [20] S.-W. Wei and Y.-X. Liu, Eur. Phys. J. Plus 136, 436 (2021), arXiv:2003.07769[gr-qc doi: 10.1140/epjp/s13360-021-01398-9
    [21] M. Ghasemi-Nodehi, M. Azreg-Aïnou, K. Jusufi et al., Phys. Rev. D 102, 104032 (2020), arXiv:2011.02276[gr-qc doi: 10.1103/PhysRevD.102.104032
    [22] P. Kocherlakota et al. (EHT Collaboration), Phys. Rev. D. 103, 104047 (2021), arXiv:2105.09343 doi: 10.1103/PhysRevD.103.104047
    [23] P.-Z. He, Q.-Q. Fan, H.-R. Zhang et al., Eur. Phys. J. C 80, 1195 (2020), arXiv:2009.06705[gr-qc doi: 10.1140/epjc/s10052-020-08707-z
    [24] A. de Vries, Classical and Quantum Gravity 17, 123 (2000) doi: 10.1088/0264-9381/17/1/309
    [25] A. A. Abdujabbarov, L. Rezzolla, and B. J. Ahmedov, Mon. Not. R. Astron. Soc. 454, 2423 (2015), arXiv:1503.09054[gr-qc doi: 10.1093/mnras/stv2079
    [26] A. Grenzebach, V. Perlick, and C. Lämmerzahl, Phys. Rev. D 89, 124004 (2014), arXiv:1403.5234[gr-qc doi: 10.1103/PhysRevD.89.124004
    [27] X. Hou, Z. Xu, and J. Wang, Mon. Not. R. Astron. Soc. (2021), arXiv:2103.11417[gr-qc
    [28] V. Perlick, O. Y. Tsupko, and G. S. Bisnovatyi-Kogan, Phys. Rev. D. 97, 104062 (2018), arXiv:1804.04898[grqc doi: 10.1103/PhysRevD.97.104062
    [29] P. V. P. Cunha, N. A. Eiró, C. A. R. Herdeiro et al., JCAP 2020, 035 (2020), arXiv:1912.08833[gr-qc doi: 10.1088/1475-7516/2020/03/035
    [30] M. Afrin, R. Kumar, and S. G. Ghosh, Mon. Not. R. Astron. Soc. 504, 5927 (2021), arXiv:2103.11417 doi: 10.1093/mnras/stab1260
    [31] P. Bambhaniya, D. Dey, A. B. Joshi et al., Phys. Rev. D 103, 084005 (2021), arXiv:2101.03865 doi: 10.1103/PhysRevD.103.084005
    [32] U. Papnoi, F. Atamurotov, S. G. Ghosh et al., Phys. Rev. D 90, 024073 (2014), arXiv:1407.0834[gr-qc doi: 10.1103/PhysRevD.90.024073
    [33] P. V. Cunha, C. A. Herdeiro, B. Kleihaus et al., Phys. Lett. B 768, 373 (2006) doi: 10.1016/j.physletb.2017.03.020
    [34] F. Atamurotov, U. Papnoi, and K. Jusufi, Classical and Quantum Gravity 39, 025014 (2022), arXiv:2104.14898[gr-qc doi: 10.1088/1361-6382/ac3e76
    [35] F. Atamurotov, I. Hussain, G. Mustafa et al., Eur. Phys. J. C 82, 831 (2022), arXiv:2209.01652[gr-qc doi: 10.1140/epjc/s10052-022-10782-3
    [36] F. Sarikulov, F. Atamurotov, A. Abdujabbarov et al., Eur. Phys. J. C 82, 771 (2022) doi: 10.1140/epjc/s10052-022-10711-4
    [37] U. Papnoi and F. Atamurotov, Physics of the Dark Universe 35, 100916 (2022), arXiv:2111.15523[gr-qc doi: 10.1016/j.dark.2021.100916
    [38] G. Mustafa, F. Atamurotov, I. Hussain et al., Chin. Phys. C 46, 125107 (2022), arXiv:2207.07608[gr-qc doi: 10.1088/1674-1137/ac917f
    [39] V. Bozza, S. Capozziello, G. Iovane et al., General Relativity and Gravitation 33, 1535 (2001), arXiv:gr-qc/0102068[gr-qc doi: 10.1023/A:1012292927358
    [40] Bozza, Phys. Rev. D. 66, 103001 (2002) doi: 10.1103/PhysRevD.66.103001
    [41] V. Bozza, F. D. Luca, and G. Scarpetta, Phys. Rev. D 74, 063001 (2006) doi: 10.1103/PhysRevD.74.063001
    [42] E. F. Eiroa and D. F. Torres, Phys. Rev. D. 69, 063004 (2004) doi: 10.1103/PhysRevD.69.063004
    [43] K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D. 62, 084003 (2000) doi: 10.1103/PhysRevD.62.084003
    [44] K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D. 65, 103004 (2002) doi: 10.1103/PhysRevD.65.103004
    [45] K. S. Virbhadra, Phys. Rev. D. 79, 083004 (2009) doi: 10.1103/PhysRevD.79.083004
    [46] S. U. Islam, R. Kumar, and S. G. Ghosh, JCAP 09, 030 (2020) doi: 10.1088/1475-7516/2020/09/030
    [47] C.-Y. Wang, Y.-F. Shen, and Y. Xie, JCAP 2019, 022 (2019), arXiv:1902.03789[gr-qc doi: 10.1088/1475-7516/2019/04/022
    [48] X. Lu and Y. Xie, Eur. Phys. J. C 79, 1016 (2019) doi: 10.1140/epjc/s10052-019-7537-2
    [49] Y.-X. Gao and Y. Xie, Phys. Rev. D 103, 043008 (2021) doi: 10.1103/PhysRevD.103.043008
    [50] V. Perlick, O. Y. Tsupko, and G. S. Bisnovatyi-Kogan, Phys. Rev. D 92, 104031 (2015), arXiv:1507.04217[gr-qc doi: 10.1103/PhysRevD.92.104031
    [51] V. Perlick and O. Y. Tsupko, Phys. Rev. D. 95, 104003 (2017), arXiv:1702.08768[gr-qc doi: 10.1103/PhysRevD.95.104003
    [52] A. Chowdhuri and A. Bhattacharyya, Phys. Rev. D 104, 064039 (2021), arXiv:2012.12914[gr-qc doi: 10.1103/PhysRevD.104.064039
    [53] F. Atamurotov, B. Ahmedov, and A. Abdujabbarov, Phys. Rev. D 92, 084005 (2015), arXiv:1507.08131[gr-qc doi: 10.1103/PhysRevD.92.084005
    [54] J. Badía and E. F. Eiroa, Phys. Rev. D 104, 084055 (2021), arXiv:2106.07601[gr-qc doi: 10.1103/PhysRevD.104.084055
    [55] F. Atamurotov, K. Jusufi, M. Jamil et al., Phys. Rev. D 104, 064053 (2021), arXiv:2109.08150[gr-qc doi: 10.1103/PhysRevD.104.064053
    [56] G. Z. Babar, A. Z. Babar, and F. Atamurotov, Eur. Phys. J. C 80, 761 (2020), arXiv:2008.05845[gr-qc doi: 10.1140/epjc/s10052-020-8346-3
    [57] M. Fathi, M. Olivares, and J. R. Villanueva, arXiv: 2104.07721
    [58] E. Ghorani, B. Puliçe, F. Atamurotov et al., Eur. Phys. J. C 83, 318 (2023), arXiv:2304.03660[gr-qc doi: 10.1140/epjc/s10052-023-11490-2
    [59] G. S. Bisnovatyi-Kogan and O. Y. Tsupko, Mon. Not. R. Astron. Soc 404, 1790 (2010), arXiv:1006.2321[astro-ph.CO doi: 10.1111/j.1365-2966.2010.16290.x
    [60] A. Rogers, Mon. Not. R. Astron. Soc. 451, 17 (2015) doi: 10.1093/mnras/stv903
    [61] F. Atamurotov, A. Abdujabbarov, and W.-B. Han, Phys. Rev. D 104, 084015 (2021) doi: 10.1103/PhysRevD.104.084015
    [62] A. Abdujabbarov, B. Ahmedov, N. Dadhich et al., Phys. Rev. D 96, 084017 (2017) doi: 10.1103/PhysRevD.96.084017
    [63] F. Atamurotov, A. Abdujabbarov, and J. Rayimbaev, Eur. Phys. J. C. 81, 118 (2021) doi: 10.1140/epjc/s10052-021-08919-x
    [64] G. Z. Babar, F. Atamurotov, and A. Z. Babar, Physics of the Dark Universe 32, 100798 (2021) doi: 10.1016/j.dark.2021.100798
    [65] G. Z. Babar, F. Atamurotov, S. Ul Islam et al., Phys. Rev. D 103, 084057 (2021), arXiv:2104.00714[gr-qc doi: 10.1103/PhysRevD.103.084057
    [66] A. Abdujabbarov, B. Toshmatov, J. Schee et al., Int. J. Mod. Phys. D 26, 1741011 (2017) doi: 10.1142/S0218271817410115
    [67] F. Atamurotov, M. Jamil, and K. Jusufi, Chin. Phys. C 47, 035106 (2023), arXiv:2212.12949[gr-qc doi: 10.1088/1674-1137/acaef7
    [68] W. Javed, I. Hussain, and A. Övgün, Eur. Phys. J. Plus 137, 148 (2022), arXiv:2201.09879[gr-qc doi: 10.1140/epjp/s13360-022-02374-7
    [69] A. Hakimov and F. Atamurotov, Astrophys. Space Sci. 361, 112 (2016) doi: 10.1007/s10509-016-2702-7
    [70] F. Atamurotov, S. Shaymatov, P. Sheoran et al., JCAP 2021, 045 (2021), arXiv:2105.02214[gr-qc doi: 10.1088/1475-7516/2021/08/045
    [71] F. Atamurotov, S. Shaymatov, and B. Ahmedov, Galaxies 9, 54 (2021) doi: 10.3390/galaxies9030054
    [72] F. Atamurotov, D. Ortiqboev, A. Abdujabbarov et al., Eur. Phys. J. C 82, 659 (2022) doi: 10.1140/epjc/s10052-022-10619-z
    [73] F. Atamurotov and S. G. Ghosh, Eur. Phys. J. Plus 137, 662 (2022) doi: 10.1140/epjp/s13360-022-02885-3
    [74] F. Atamurotov, M. Alloqulov, A. Abdujabbarov et al., Eur. Phys. J. Plus 137, 634 (2022) doi: 10.1140/epjp/s13360-022-02846-w
    [75] F. Atamurotov, F. Sarikulov, V. Khamidov et al., Eur. Phys. J. Plus 137, 567 (2022) doi: 10.1140/epjp/s13360-022-02780-x
    [76] V. S. Morozova, B. J. Ahmedov, and A. A. Tursunov, Astrophys. Space Sci. 346, 513 (2013) doi: 10.1007/s10509-013-1458-6
    [77] F. Atamurotov, F. Sarikulov, A. Abdujabbarov et al., Eur. Phys. J. Plus 137, 336 (2022) doi: 10.1140/epjp/s13360-022-02548-3
    [78] F. Atamurotov, I. Hussain, G. Mustafa et al., Chin. Phys. C 47, 025102 (2023) doi: 10.1088/1674-1137/ac9fbb
    [79] F. Atamurotov, H. Alibekov, A. Abdujabbarov et al., Symmetry 15, 848 (2023) doi: 10.3390/sym15040848
    [80] A. Ditta, X. Tiecheng, S. Mumtaz et al., arXiv:2303.05438
    [81] A. Simpson and M. Visser, JCAP 2019, 042 (2019), arXiv:1812.07114 doi: 10.1088/1475-7516/2019/02/042
    [82] E. Franzin, S. Liberati, J. Mazza et al., JCAP 2021, 036 (2021), arXiv:2104.11376[grqc doi: 10.1088/1475-7516/2021/07/036
    [83] F. S. N. Lobo, M. E. Rodrigues, M. V. d. S. Silva et al., Phys. Rev. D 103, 084052 (2021), arXiv:2009.12057[gr-qc doi: 10.1103/PhysRevD.103.084052
    [84] N. Tsukamoto, Phys. Rev. D 105, 084036 (2022), arXiv:2202.09641[gr-qc doi: 10.1103/PhysRevD.105.084036
    [85] J. L. Synge, Relativity: The General Theory, (North-Holland, Amsterdam, 1960)
  • [1] B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016) doi: 10.1103/PhysRevLett.116.061102
    [2] B. Abbott, R. Abbott, T. Abbott et al., Phys. Rev. Let. 116 (2016) doi: 10.1103/physrevlett.116.241102
    [3] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. 875, L1 (2019), arXiv:1906.11238[astro-ph.GA doi: 10.3847/2041-8213/ab0ec7
    [4] K. Akiyama et al., Astrophys. J. 875, L6 (2019), arXiv:1906.11243[astro-ph.GA doi: 10.3847/2041-8213/ab1141
    [5] K. Akiyama et al., Astrophys. J. Lett 930, L12 (2022) doi: 10.3847/2041-8213/ac6674
    [6] C. M. Will, Living Reviews in Relativity 9, 3 (2006), arXiv:gr-qc/0510072[gr-qc doi: 10.12942/lrr-2006-3
    [7] J. L. Synge, Mon. Not. R. Astron. Soc. 131, 463 (1966) doi: 10.1093/mnras/131.3.463
    [8] J. P. Luminet, Astron.Astrophys. 75, 228 (1979)
    [9] J. M. Bardeen, in Black Holes ((Les Astres Occlus, 1973) pp. 215–239
    [10] H. Falcke, F. Melia, and E. Agol, Astrophys. J. 523, L18 (2000), arXiv:9912263 doi: 10.1086/312423
    [11] C. Bambi and K. Freese, Phys. Rev. D 79, 043002 (2009), arXiv:0812.1328[astro-ph doi: 10.1103/PhysRevD.79.043002
    [12] K. Hioki and K.-I. Maeda, Phys. Rev. D 80, 024042 (2009) doi: 10.1103/PhysRevD.80.024042
    [13] A. Abdujabbarov, F. Atamurotov, Y. Kucukakca et al., Astrophys. Space. Sci. 344, 429 (2013), arXiv:1212.4949[physics.gen-ph doi: 10.1007/s10509-012-1337-6
    [14] L. Amarilla and E. F. Eiroa, Phys. Rev. D 85, 064019 (2012) doi: 10.1103/PhysRevD.85.064019
    [15] A. Abdujabbarov, M. Amir, B. Ahmedov et al., Phys. Rev. D 93, 104004 (2016), arXiv:1604.03809[gr-qc doi: 10.1103/PhysRevD.93.104004
    [16] F. Atamurotov, A. Abdujabbarov, and B. Ahmedov, Phys. Rev. D 88, 064004 (2013) doi: 10.1103/PhysRevD.88.064004
    [17] N. Tsukamoto, Phys. Rev. D. 97, 064021 (2018), arXiv:1708.07427[gr-qc doi: 10.1103/PhysRevD.97.064021
    [18] A. Abdujabbarov, F. Atamurotov, N. Dadhich et al., Eur. Phys. J. C 75, 399 (2015), arXiv:1508.00331[gr-qc doi: 10.1140/epjc/s10052-015-3604-5
    [19] R. C. Pantig, P. K. Yu, E. T. Rodulfo et al., Annals of Physics 436, 168722 (2022), arXiv:2104.04304[gr-qc doi: 10.1016/j.aop.2021.168722
    [20] S.-W. Wei and Y.-X. Liu, Eur. Phys. J. Plus 136, 436 (2021), arXiv:2003.07769[gr-qc doi: 10.1140/epjp/s13360-021-01398-9
    [21] M. Ghasemi-Nodehi, M. Azreg-Aïnou, K. Jusufi et al., Phys. Rev. D 102, 104032 (2020), arXiv:2011.02276[gr-qc doi: 10.1103/PhysRevD.102.104032
    [22] P. Kocherlakota et al. (EHT Collaboration), Phys. Rev. D. 103, 104047 (2021), arXiv:2105.09343 doi: 10.1103/PhysRevD.103.104047
    [23] P.-Z. He, Q.-Q. Fan, H.-R. Zhang et al., Eur. Phys. J. C 80, 1195 (2020), arXiv:2009.06705[gr-qc doi: 10.1140/epjc/s10052-020-08707-z
    [24] A. de Vries, Classical and Quantum Gravity 17, 123 (2000) doi: 10.1088/0264-9381/17/1/309
    [25] A. A. Abdujabbarov, L. Rezzolla, and B. J. Ahmedov, Mon. Not. R. Astron. Soc. 454, 2423 (2015), arXiv:1503.09054[gr-qc doi: 10.1093/mnras/stv2079
    [26] A. Grenzebach, V. Perlick, and C. Lämmerzahl, Phys. Rev. D 89, 124004 (2014), arXiv:1403.5234[gr-qc doi: 10.1103/PhysRevD.89.124004
    [27] X. Hou, Z. Xu, and J. Wang, Mon. Not. R. Astron. Soc. (2021), arXiv:2103.11417[gr-qc
    [28] V. Perlick, O. Y. Tsupko, and G. S. Bisnovatyi-Kogan, Phys. Rev. D. 97, 104062 (2018), arXiv:1804.04898[grqc doi: 10.1103/PhysRevD.97.104062
    [29] P. V. P. Cunha, N. A. Eiró, C. A. R. Herdeiro et al., JCAP 2020, 035 (2020), arXiv:1912.08833[gr-qc doi: 10.1088/1475-7516/2020/03/035
    [30] M. Afrin, R. Kumar, and S. G. Ghosh, Mon. Not. R. Astron. Soc. 504, 5927 (2021), arXiv:2103.11417 doi: 10.1093/mnras/stab1260
    [31] P. Bambhaniya, D. Dey, A. B. Joshi et al., Phys. Rev. D 103, 084005 (2021), arXiv:2101.03865 doi: 10.1103/PhysRevD.103.084005
    [32] U. Papnoi, F. Atamurotov, S. G. Ghosh et al., Phys. Rev. D 90, 024073 (2014), arXiv:1407.0834[gr-qc doi: 10.1103/PhysRevD.90.024073
    [33] P. V. Cunha, C. A. Herdeiro, B. Kleihaus et al., Phys. Lett. B 768, 373 (2006) doi: 10.1016/j.physletb.2017.03.020
    [34] F. Atamurotov, U. Papnoi, and K. Jusufi, Classical and Quantum Gravity 39, 025014 (2022), arXiv:2104.14898[gr-qc doi: 10.1088/1361-6382/ac3e76
    [35] F. Atamurotov, I. Hussain, G. Mustafa et al., Eur. Phys. J. C 82, 831 (2022), arXiv:2209.01652[gr-qc doi: 10.1140/epjc/s10052-022-10782-3
    [36] F. Sarikulov, F. Atamurotov, A. Abdujabbarov et al., Eur. Phys. J. C 82, 771 (2022) doi: 10.1140/epjc/s10052-022-10711-4
    [37] U. Papnoi and F. Atamurotov, Physics of the Dark Universe 35, 100916 (2022), arXiv:2111.15523[gr-qc doi: 10.1016/j.dark.2021.100916
    [38] G. Mustafa, F. Atamurotov, I. Hussain et al., Chin. Phys. C 46, 125107 (2022), arXiv:2207.07608[gr-qc doi: 10.1088/1674-1137/ac917f
    [39] V. Bozza, S. Capozziello, G. Iovane et al., General Relativity and Gravitation 33, 1535 (2001), arXiv:gr-qc/0102068[gr-qc doi: 10.1023/A:1012292927358
    [40] Bozza, Phys. Rev. D. 66, 103001 (2002) doi: 10.1103/PhysRevD.66.103001
    [41] V. Bozza, F. D. Luca, and G. Scarpetta, Phys. Rev. D 74, 063001 (2006) doi: 10.1103/PhysRevD.74.063001
    [42] E. F. Eiroa and D. F. Torres, Phys. Rev. D. 69, 063004 (2004) doi: 10.1103/PhysRevD.69.063004
    [43] K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D. 62, 084003 (2000) doi: 10.1103/PhysRevD.62.084003
    [44] K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D. 65, 103004 (2002) doi: 10.1103/PhysRevD.65.103004
    [45] K. S. Virbhadra, Phys. Rev. D. 79, 083004 (2009) doi: 10.1103/PhysRevD.79.083004
    [46] S. U. Islam, R. Kumar, and S. G. Ghosh, JCAP 09, 030 (2020) doi: 10.1088/1475-7516/2020/09/030
    [47] C.-Y. Wang, Y.-F. Shen, and Y. Xie, JCAP 2019, 022 (2019), arXiv:1902.03789[gr-qc doi: 10.1088/1475-7516/2019/04/022
    [48] X. Lu and Y. Xie, Eur. Phys. J. C 79, 1016 (2019) doi: 10.1140/epjc/s10052-019-7537-2
    [49] Y.-X. Gao and Y. Xie, Phys. Rev. D 103, 043008 (2021) doi: 10.1103/PhysRevD.103.043008
    [50] V. Perlick, O. Y. Tsupko, and G. S. Bisnovatyi-Kogan, Phys. Rev. D 92, 104031 (2015), arXiv:1507.04217[gr-qc doi: 10.1103/PhysRevD.92.104031
    [51] V. Perlick and O. Y. Tsupko, Phys. Rev. D. 95, 104003 (2017), arXiv:1702.08768[gr-qc doi: 10.1103/PhysRevD.95.104003
    [52] A. Chowdhuri and A. Bhattacharyya, Phys. Rev. D 104, 064039 (2021), arXiv:2012.12914[gr-qc doi: 10.1103/PhysRevD.104.064039
    [53] F. Atamurotov, B. Ahmedov, and A. Abdujabbarov, Phys. Rev. D 92, 084005 (2015), arXiv:1507.08131[gr-qc doi: 10.1103/PhysRevD.92.084005
    [54] J. Badía and E. F. Eiroa, Phys. Rev. D 104, 084055 (2021), arXiv:2106.07601[gr-qc doi: 10.1103/PhysRevD.104.084055
    [55] F. Atamurotov, K. Jusufi, M. Jamil et al., Phys. Rev. D 104, 064053 (2021), arXiv:2109.08150[gr-qc doi: 10.1103/PhysRevD.104.064053
    [56] G. Z. Babar, A. Z. Babar, and F. Atamurotov, Eur. Phys. J. C 80, 761 (2020), arXiv:2008.05845[gr-qc doi: 10.1140/epjc/s10052-020-8346-3
    [57] M. Fathi, M. Olivares, and J. R. Villanueva, arXiv: 2104.07721
    [58] E. Ghorani, B. Puliçe, F. Atamurotov et al., Eur. Phys. J. C 83, 318 (2023), arXiv:2304.03660[gr-qc doi: 10.1140/epjc/s10052-023-11490-2
    [59] G. S. Bisnovatyi-Kogan and O. Y. Tsupko, Mon. Not. R. Astron. Soc 404, 1790 (2010), arXiv:1006.2321[astro-ph.CO doi: 10.1111/j.1365-2966.2010.16290.x
    [60] A. Rogers, Mon. Not. R. Astron. Soc. 451, 17 (2015) doi: 10.1093/mnras/stv903
    [61] F. Atamurotov, A. Abdujabbarov, and W.-B. Han, Phys. Rev. D 104, 084015 (2021) doi: 10.1103/PhysRevD.104.084015
    [62] A. Abdujabbarov, B. Ahmedov, N. Dadhich et al., Phys. Rev. D 96, 084017 (2017) doi: 10.1103/PhysRevD.96.084017
    [63] F. Atamurotov, A. Abdujabbarov, and J. Rayimbaev, Eur. Phys. J. C. 81, 118 (2021) doi: 10.1140/epjc/s10052-021-08919-x
    [64] G. Z. Babar, F. Atamurotov, and A. Z. Babar, Physics of the Dark Universe 32, 100798 (2021) doi: 10.1016/j.dark.2021.100798
    [65] G. Z. Babar, F. Atamurotov, S. Ul Islam et al., Phys. Rev. D 103, 084057 (2021), arXiv:2104.00714[gr-qc doi: 10.1103/PhysRevD.103.084057
    [66] A. Abdujabbarov, B. Toshmatov, J. Schee et al., Int. J. Mod. Phys. D 26, 1741011 (2017) doi: 10.1142/S0218271817410115
    [67] F. Atamurotov, M. Jamil, and K. Jusufi, Chin. Phys. C 47, 035106 (2023), arXiv:2212.12949[gr-qc doi: 10.1088/1674-1137/acaef7
    [68] W. Javed, I. Hussain, and A. Övgün, Eur. Phys. J. Plus 137, 148 (2022), arXiv:2201.09879[gr-qc doi: 10.1140/epjp/s13360-022-02374-7
    [69] A. Hakimov and F. Atamurotov, Astrophys. Space Sci. 361, 112 (2016) doi: 10.1007/s10509-016-2702-7
    [70] F. Atamurotov, S. Shaymatov, P. Sheoran et al., JCAP 2021, 045 (2021), arXiv:2105.02214[gr-qc doi: 10.1088/1475-7516/2021/08/045
    [71] F. Atamurotov, S. Shaymatov, and B. Ahmedov, Galaxies 9, 54 (2021) doi: 10.3390/galaxies9030054
    [72] F. Atamurotov, D. Ortiqboev, A. Abdujabbarov et al., Eur. Phys. J. C 82, 659 (2022) doi: 10.1140/epjc/s10052-022-10619-z
    [73] F. Atamurotov and S. G. Ghosh, Eur. Phys. J. Plus 137, 662 (2022) doi: 10.1140/epjp/s13360-022-02885-3
    [74] F. Atamurotov, M. Alloqulov, A. Abdujabbarov et al., Eur. Phys. J. Plus 137, 634 (2022) doi: 10.1140/epjp/s13360-022-02846-w
    [75] F. Atamurotov, F. Sarikulov, V. Khamidov et al., Eur. Phys. J. Plus 137, 567 (2022) doi: 10.1140/epjp/s13360-022-02780-x
    [76] V. S. Morozova, B. J. Ahmedov, and A. A. Tursunov, Astrophys. Space Sci. 346, 513 (2013) doi: 10.1007/s10509-013-1458-6
    [77] F. Atamurotov, F. Sarikulov, A. Abdujabbarov et al., Eur. Phys. J. Plus 137, 336 (2022) doi: 10.1140/epjp/s13360-022-02548-3
    [78] F. Atamurotov, I. Hussain, G. Mustafa et al., Chin. Phys. C 47, 025102 (2023) doi: 10.1088/1674-1137/ac9fbb
    [79] F. Atamurotov, H. Alibekov, A. Abdujabbarov et al., Symmetry 15, 848 (2023) doi: 10.3390/sym15040848
    [80] A. Ditta, X. Tiecheng, S. Mumtaz et al., arXiv:2303.05438
    [81] A. Simpson and M. Visser, JCAP 2019, 042 (2019), arXiv:1812.07114 doi: 10.1088/1475-7516/2019/02/042
    [82] E. Franzin, S. Liberati, J. Mazza et al., JCAP 2021, 036 (2021), arXiv:2104.11376[grqc doi: 10.1088/1475-7516/2021/07/036
    [83] F. S. N. Lobo, M. E. Rodrigues, M. V. d. S. Silva et al., Phys. Rev. D 103, 084052 (2021), arXiv:2009.12057[gr-qc doi: 10.1103/PhysRevD.103.084052
    [84] N. Tsukamoto, Phys. Rev. D 105, 084036 (2022), arXiv:2202.09641[gr-qc doi: 10.1103/PhysRevD.105.084036
    [85] J. L. Synge, Relativity: The General Theory, (North-Holland, Amsterdam, 1960)
  • 加载中

Figures(13)

Get Citation
Furkat Sarikulov, Farruh Atamurotov, Ahmadjon Abdujabbarov and Vokhid Khamidov. The photon motion and weak gravitational lensing in a black-bounce spacetime[J]. Chinese Physics C. doi: 10.1088/1674-1137/acedf2
Furkat Sarikulov, Farruh Atamurotov, Ahmadjon Abdujabbarov and Vokhid Khamidov. The photon motion and weak gravitational lensing in a black-bounce spacetime[J]. Chinese Physics C.  doi: 10.1088/1674-1137/acedf2 shu
Milestone
Received: 2023-05-10
Article Metric

Article Views(1273)
PDF Downloads(36)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Photon motion and weak gravitational lensing in black-bounce spacetime

  • 1. Ulugh Beg Astronomical Institute, Astronomy St 33, Tashkent 100052, Uzbekistan
  • 2. New Uzbekistan University, Mustaqillik ave. 54, 100007 Tashkent, Uzbekistan
  • 3. Akfa University, Milliy Bog’ Street 264, Tashkent 111221, Uzbekistan
  • 4. Inha University in Tashkent, Ziyolilar 9, Tashkent 100170, Uzbekistan
  • 5. Institute of Fundamental and Applied Research, National Research University TIIAME, Kori Niyoziy 39,Tashkent 100000, Uzbekistan
  • 6. National University of Uzbekistan, Tashkent 100174, Uzbekistan
  • 7. Tashkent State Technical University, Tashkent 100095, Uzbekistan
  • 8. Tashkent University of Information Technologies named after Muhammad al Khwarizmi, Amir Temur 108, Tashkent 100014, Uzbekistan

Abstract: The effect of spacetime curvature on photon motion may offer an opportunity to propose new tests on gravity theories. In this study, we investigate and focus on the massless (photon) particle motion around black-bounce gravity. We analyze the horizon structure around a gravitational compact object described by black-bounce spacetime. The photon motion and the effect of gravitational weak lensing in vacuum and plasma are discussed, and the shadow radius of the compact object is also studied in black-bounce spacetime. Additionally, the magnification of the image is studied using the deflection angle of light rays.

    HTML

    I.   INTRODUCTION
    • The effect of a curved spacetime on wave propagation is an important subject in modern relativistic astrophysics. The first ever observation of gravitational waves by the LIGO-Virgo collaboration [1, 2] was the starting point for a new era of gravitational wave observation. Later, the observation of images of M87 [3, 4] and SgrA* [5] motivated specialists to use these observations to test gravity theories in the strong field regime.

      Although well tested in the weak [6] and strong field regimes [2, 3], Einstein's general relativity still encounters several fundamental problems. One is the existence of the physical (curvature) singularity at the origin of most solutions within the theory. There is a belief that an undiscovered new quantum gravity theory will resolve this issue. One of the main aims of introducing alternative theories of gravity and modified versions of general relativity is to move in the direction of constructing a quantum gravity theory.

      Owing to the strong gravitational field, the spacetime around massive objects is curved. The propagation of an electromagnetic wave is affected by this curvature, and an observer at infinity can detect light deflection due to gravity. However, because of the strong gravitational field around black holes, the part of the light rays propagating near a black hole will be captured by the central object. Consequently, the observer at infinity will detect a black spot on the celestial plane near the position of the black hole. The original idea of the possible observation of this black spot, usually referred to as a black hole shadow, was proposed by Synge [7] and later developed by Luminet [8] and Bardeen [9]. The shadow of a black hole for different solutions within general relativity and/or modified/alternative theories of gravity has subsequently been investigated by the various authors [1038].

      The deflection of light, also known as gravitational lensing, around a gravitating object can be used to either study the light source or test the gravitational field around the lensing object. The effect of gravitational lensing was first used to test general relativity by Eddington in 1919 during the solar eclipse. Later, the gravitational lensing around compact objects within different gravity models was investigated by various authors (see, for example, [3949]). The existence of a plasma environment is one of the main markers of most black holes and compact objects. The light propagating through the plasma is also affected by the latter. The effects of both uniform and nonuniform plasma on the optical properties of black holes were investigated in Refs. [5058]). Moreover, using gravitational lensing, one may study the effects of different plasma configurations on the light deflection angle (see, for example, [5980]).

      A regularization method was proposed by Simpson and Visser [81], in which the expression 1/r in the Schwarzschild metric is replaced with 1/r2+a2. As a result, one may obtain a globally regular family of metrics, including black-bounce geometry and a traversable wormhole. A particularly interesting feature of these solutions is that they may be considered a one-parameter extension of the Schwarzschild solution. The generalization of Reissner-Nordström spacetime using this technique has been reported in Ref. [82]. Later, the authors of Ref. [83] extended them and proposed a number of new solutions for black-bounce spacetime. Recently, the effect of retrolensing by two-photon spheres within black-bounce spacetime was studied in detail in Ref. [84]. In this paper, we extend the analysis of this spacetime by studying massless test particle motion around a black-bounce compact object. Moreover, we also study the effect of the gravitational lensing and shadow of the compact object described by black-bounce spacetime.

      The paper is arranged as follows. First, we review the spacetime metric in Sec. II. We study the photon motion around the black bounce compact object in Sec. III. Then, the shadow of the black-bounce object is studied in Sec. IV in the presence of plasma. Gravitational lensing in the presence of plasma in black bounce spacetime is then investigated in Sec. V. We conclude our study in Sec. VI. Throughout, we use the unit system G=c=1.

    II.   SPACETIME METRIC
    • Non-rotating spacetime in black bounce gravity has the following linear element in spherical coordinates [83, 84]:

      ds2=f(r)dt2+dr2f(r)+γ(r)(dθ2+sin2θdϕ2) ,

      (1)

      where the metric functions f(r) and γ(r) are defined as

      f(r)=12Mr2(a2+r2)3/2 ,

      (2)

      γ(r)=r2+a2 ,

      (3)

      where M is the mass of the compact object, and a is a nonnegative constant [83]. To understand the nature of the gravitational potential of the source, we may introduce the effective mass as Meff=M(r/r2+a2)3.

      Spacetime metric (1) describes a regular black hole when the condition 0<a/m<43/9 is satisfied. To obtain more information about the spacetime structure, we can analyze the function f(r). The dependence of f(r) on radial coordinates is shown in Fig. 1. From the plot, we can easily see that there is an extreme value of a/M when the graph has one r/M-intercept. The extreme value of the parameter aext=43/9. The values of parameter a that are less than its extreme value correspond to a black hole with two horizons (Fig. 2 shows the existence of two roots for a set of parameters a values, which corresponds to the Cauchy horizon (smaller root) and event horizon (larger root)), whereas the values of parameter a that are greater than its extreme value correspond to a naked singularity. The dependence of the radius of the horizon corresponding to the condition f(r)=0 is represented in Fig. 2. From the plot, we can see that with an increase in parameter a in black-bounce spacetime, the outer horizon decreases until the value M is reached.

      Figure 1.  (color online) Radial dependence of the lapse function for different values of the parameter a/M.

      Figure 2.  (color online) Horizon radius of a black hole as a function of a/M.

    III.   PHOTON MOTION
    • Consider the motion of photons around a gravitational object in the framework of black-bounce gravity. To describe the dynamics, we can use the Hamilton-Jacobi equation in the following form [85]:

      gμνdSdxμdSdxν=0,

      (4)

      where gμν is the metric tensor of spacetime metric (1), and S is the action. Using the separation of variables method and the symmetry of spacetime metric (1), the action S can be expressed as the following form:

      S=Et+Lφ+Sr(r)+Sθ(θ),

      (5)

      where E and L can be described as the energy and angular momentum of a massless particle, respectively. The dynamics of particles with zero rest mass are defined by null geodesics of spacetime. Using Eqs. (4) and (5), we can obtain the equations of motion of massless particles at the equatorial plane (θ=0.5π, ˙θ=0) as

      dtdλ=Ef(r),

      (6)

      drdλ=±R(r),

      (7)

      dϕdλ=Lγ(r),

      (8)

      with

      R(r)=E2f(r)γ(r)L2,

      (9)

      where λ is the affine parameter, and we consider θ=0.5π. Accordingly, the effective potential for the radial motion of photons takes the following form:

      Veff=f(r)L2γ(r).

      (10)

      The study of circular orbits r reflects significant interest while considering the photon (massless particle) motion. Photon circular orbits rph are determined using the following conditions:

      Veff=0.

      (11)

      Using the above conditions, we can obtain the equation to determine the radius of circular photon orbits rph as

      (a2+r2)3/2+2Ma23Mr2=0.

      (12)

      We can obtain the radius of the photon sphere (rph) by solving Eq. (12) with respect to the radial coordinate,

      rph=a2(10M2D+1)+D+9M4D+3M2 ,

      (13)

      where D is defined as

      D=(25a4M2245a2M4+552a3M25a24M2+27M6)1/3.

      (14)

      The dependence of the photon orbit rph on the parameter a is shown in Fig. 3. From the plot, we can conclude that with an increase in parameter a, the radius of the photon orbit rph decreases.

      Figure 3.  (color online) Dependence of the radius of photon circular orbits rph/M on parameter a.

    IV.   BLACK HOLE SHADOW IN THE PRESENCE OF PLASMA
    • The detailed study of photon motion around a black hole leads to a phenomenon known as a black hole shadow [16]. The shadow cast by a black hole is a result of photon capture by the central object, and the observer will see a black spot on the bright background [16]. In this section, we explore the radius (Rsh) of the shadow of a black hole in the presence of a plasma medium [55]. The boundary of the shadow is fully defined by the photon (massless particle) trajectory governed by the equation of motion in the plasma medium. To analyze the shadow of the black hole in black-bounce spacetime in a plasma medium, we should first study the photon equation in a plasma environment.

    • A.   Photon motion in plasma

    • The Hamiltonian for photons around a black hole in the presence of an electron plasma medium can be expressed as [55, 85]

      H=12[gμνpμpν(n21)p2tf(r)]=0,

      (15)

      with

      n2=1ω2eω2(xi),

      (16)

      where ω(xi) and ωe are the photon and plasma frequencies, respectively, and n is the refractive index of the medium. We denote the value of the photon frequency observed at infinity as ω0:

      ω()=ω0=pt.

      (17)

      The components of the four-velocity (˙t, ˙ϕ, ˙r) for the photons (massless particle) in the equatorial plane (θ=0.5π, pθ=0) are given by

      dtdλ=ptf(r),

      (18)

      drdλ=prf(r),

      (19)

      dϕdλ=pϕr2+a2,

      (20)

      where we consider the relationship dxμ/dλ=H/pμ. From Eqs. (19) and (20), we can obtain an expression for the phase trajectory of light (photon) rays,

      drdϕ=grrprgϕϕpϕ

      (21)

      Using the constraint H=0, we can define the above equation as [50]

      drdϕ=grrgϕϕh2(r)ω20p2ϕ1,

      (22)

      where the following definition is introduced:

      h2(r)gttgϕϕ1gϕϕω2e(r)ω20.

      (23)

      The radius of a circular orbit of photons rph, particularly the one that forms the photon sphere of radius rph, is determined by the following equation [50]:

      d(h2(r))dr|r=rph=0.

      (24)

      Substituting Eq. (24) into Eq. (23), we can define an algebraic expression for rph in the plasma environment as

      (a2+r2)f(r)f(r)+2f(r)ωe(r)ω20[(a2+r2)ωe(r)+rωe(r)]=2r,

      (25)

      where the prime stands for the derivative with respect to r. In the case of homogeneous plasma (ω2e(r)=const.), Eq. (25) can be simplified to

      (a2+r2)f(r)f(r)+2rf(r)ω2eω20=2r.

      (26)

      As mentioned above, the solution to this equation with respect to r leads to the radius of the photon sphere rph of the black hole in a uniform plasma environment.

      Figure 4 gives information on the dependence of the radius of the photon sphere rph on the black-bounce parameter a and plasma parameters. It is clear from the plots that the radius of the photon sphere rph decreases with an increase in parameter a. Note that the presence of plasma causes an increase in the value of the photon sphere radius, as shown in Fig. 4.

      Figure 4.  (color online) Dependence of the radius of the photon sphererph/M on plasma frequency (top) and black-bounce parameter a (bottom).

    • B.   Black hole shadow in uniform plasma

    • In this section, we study the radius (Rsh) of the shadow of the black hole in the presence of uniform plasma. The angular radius αsh of the black hole shadow is defined by [7, 50]

      sin2αsh=h2(rsh)h2(ro)=(a2+r2ph)(1f(rph)ω2e(rph)ω20)(a2+r2o)(1f(ro)ω2e(ro)ω20),

      (27)

      where ro and rph are the radial positions of the observer and photon sphere (introduced in the previous section), respectively. If the observer is located at a sufficiently large distance from the black hole, we can approximate the radius of the black hole shadow (for the case of uniform plasma) using Eq. (27) [50]:

      Rshrosinαsh=(a2+r2ph)(1f(rph)ω2eω20).

      (28)

      For the case of no plasma in the vicinity of the black-bounce black hole, the shadow radius can be represented by

      Rshrosinαsh=a2+r2phf(rph).

      (29)

      Figure 5 describes the dependence of the radius of black hole shadow on parameter a and plasma parameters. In the plots, we define Rsh in terms of Meff because it is known from EHT data that the masses of M87* and Sgr A* are given by the gravitational potential. Based on these data, we can conclude that the shadow radius is diminished due to the influence of both plasma and the black-bounce parameter a.

      Figure 5.  (color online) Dependence of the radius of the black hole shadow on parameter a (top) and plasma frequency (bottom).

      Now, we consider the rough assumption that the supermassive black holes M87* and Sgr A* are the black-bounce black holes that we are studying. Then, we can theoretically obtain the constraints from the EHT observation results.

      The angular diameter of the image of the supermassive black hole M87* is θ=42±3 μas at the 1σ confidence level, and the mass of M87* and the distance from the solar system are M=6.5×109M and D= 16.8 Mpc, respectively [3, 4]. The same data for Sgr A* are as follows: θ=48.7±7 μas, M4×106M, and D8 kpc [5]. The shadow radius caused by the black hole per unit mass can be expressed as the following equation:

      dsh=DθMeff,

      (30)

      where Meff=M(r/r2+a2)3 from black-bounce spacetime. Here, we must use the effective mass because the masses of M87 and Sgr A* are given by the gravitational potential.

      Using Eq. (30), we can easily calculate the diameter of the shadow images. The diameter of the black hole shadow is dsh=(11±1.5)Meff for M87* and dsh=(9.5±1.4)Meff for Sgr A*. From these data, we can now obtain the constraints on the black-bounce and plasma parameters for the supermassive black holes Sgr A* and M87*. The results are presented in Fig. 6 for the case of no plasma around the black hole and in Fig. 7 for the black hole in plasma.

      Figure 6.  (color online) Constraints on the black-bounce parameter a with no plasma around the black hole. The brown and black dashed lines describe the limits of the measured shadow radii of M87* and Sgr A*, respectively.

      Figure 7.  (color online) Shadow diameter of the black-bounce black hole as a function of the parameter a and plasma parameter. The black curve corresponds to dsh=9.5Meff, which is the lower limit for M87*. The blue curve corresponds to dsh=10.9Meff, which is the upper limit for Sgr A*.

    V.   WEAK GRAVITATIONAL LENSING IN THE PRESENCE OF PLASMA
    • Now, we can move on to another part of our study: gravitaional lensing in spacetime metric (1) in a weak-field case explicated via [59]

      gαβ=ηαβ+ξαβ ,

      (31)

      where ηαβ is the general metric of Minkowski spacetime, and ξαβ is a perturbation of the flat spacetime metric, which can be described by the spacetime metric theory of gravity and is written as [59]

      ηαβ=diag(1,1,1,1),ξαβ1,ξαβ0underxi ,gαβ=ηαβξαβ,   ξαβ=ξαβ .

      (32)

      Here, we can study the effects of the plasma medium on the deflection angle ˆαk in the gravitational weak field of the black hole. The general expression for the deflection angle in plasma is [59, 64]

      ˆαi=12(ξ33+ξ00ω2KeN(xi)ω2ω2e),idz,  i=1,2

      (33)

      where N(xi) is the number density of particles in a plasma medium around a compact object, Ke=4πe2/me is constant (e and me are the electron charge and mass), and ω and ωe are the photon (massless particle) and plasma frequencies, respectively [59]. Using Eqs. (32) and (33), we can rewrite Eq. (33) for the deflection angle in the following form [59]:

      ˆαb=12br(dξ33dr+11ω2e/ω2dξ00drKeω2ω2e dNdr)dz,

      (34)

      where b is the impact parameter of light rays. It is useful to note that the values of ˆαb can be both negative and positive [59].

      In the weak-field case for far distances from the black hole, the black-bounce spacetime metric (1) can be written as

      ds2=ds20+2Mr2(a2+r2)3/2dt2+2Mr2(a2+r2)3/2dr2+a2(dθ2+sin2θdϕ2),

      (35)

      where ds20=dt2+dr2+r2(dθ2+sin2θdϕ2) is the metric element in the Minkowski spacetime metric, and we can use the notation Rs=2M for further calculations.

      To analyze the deflection angle (ˆαb) of light rays around the black hole in the plasma medium using Eq. (34), we can rewrite the required components (ξ00 and ξ33) in Cartesian coordinates as

      ξ00=Rsr2(a2+r2)3/2,ξ33=Rsr2(a2+r2)3/2cos2χ+a2r2sin4χ ,

      (36)

      where cos2χ=z2/(b2+z2) and r2=b2+z2 are introduced, for example, in [64]. Here, z is the coordinate aligned with the optical axis. We can obtain the derivative of ξ00 and ξ33 with respect to the radial coordinate as

      dξ33dr=3rRsz2(a2+r2)5/22a2(r44r2z2+3z4)r7,dξ00dr=Rs(2a2rr3)(a2+r2)5/2.

      (37)

      For simplicity, the expression for the deflection angle ˆαb can be expanded as [63]

      ˆαb=ˆα1+ˆα2+ˆα3,

      (38)

      with

      ˆα1=12brdξ33drdz,ˆα2=12br(11ω2e/ω2dξ00dr)dz,ˆα3=12br(Keω2ω2e dNdr)dz,

      (39)

      where the notations ˆα1, ˆα2, and ˆα3 represent the contributions to the deflection angle due to gravity, uniform plasma, and non-uniform plasma, respectively. In this study, we use Eqs. (38) and (39) to investigate the effect of plasma on the deflection angle in gravitational weak lensing. In the following sections, we consider each case in detail.

    • A.   Uniform plasma with ω2e=const

    • Here, we can test the deflection angle of photons (light rays) around the black hole in the presence of a uniform plasma environment using expression (38) written in the following form:

      ˆαuni=ˆαuni1+ˆαuni2+ˆαuni3,

      (40)

      We can claim that the third term of this equation will vanish due to the uniform distribution of the plasma environment. Using Eqs. (38), (37), and (40), we can obtain an equation for the deflection angle of light rays (photon) around the black hole in black-bounce spacetime surrounded by uniform plasma:

      ˆαuni=bRsa2+b2+3πa216b2bRs(a2b2)(a2+b2)2ω2ω2ω2e.

      (41)

      In this equation, the second term appears owing to the property of black-bounce spacetime itself.

      Using expression (41), we can plot the deflection angle against the impact b, plasma frequency, and a. Fig. 8 shows the dependence of the deflection angle ˆαuni of light rays (photon) around the black hole on the impact parameter b for different values of the parameter a and plasma parameter ω2e/ω2. The dependence of the deflection angle ˆαuni plasma frequency and parameter a for fixed values of the impact parameter is shown in Fig. 9. The graphs of the defection angle ˆαuni on the impact parameter b and plasma parameter are reliable, as expected. The presence of uniform plasma leads to an increase in the value of the deflection angle of light rays, whereas the angle of deflection decreases dramatically with increasing impact parameter b. However, interestingly, the deflection angle remains almost unchanged with the influence of the black-bounce parameter a. There is only a slight increase in the value of the deflection angle with an increase in parameter a, as shown in Figs. 8 and 9.

      Figure 8.  (color online) Dependence of the deflection angle ˆαuni on the impact parameter b for different values of the parameter a/M (left panel) and plasma medium ω2e/ω2 (right panel).

      Figure 9.  (color online) Dependence of the deflection angle ˆαuni on the parameter a/M (left panel) and plasma parameter ω2e/ω2 (right panel) for a fixed value of the impact parameter b/M=30.

    • B.   Non-uniform plasma with singular isothermal sphere medium

    • Here, we investigate the deflection angle ˆαuni of light rays (photons) around the black hole in the presence of non-uniform plasma and in the framework of black-bounce spacetime with a nonzero parameter a. We can consider the singular isothermal sphere (SIS) medium as the non-uniform plasma distribution [59, 64]. The plasma number density of the SIS environment is written in the following form [59, 63]:

      N(r)=ρ(r)kmp ,

      (42)

      where mp is the proton mass, and k is a non-dimensional coefficient responsible for the dark matter contribution in the SIS model and with

      ρ(r)=σ2ν2πr2 ,

      (43)

      where σν and ρ(r) are the velocity of dispersion and the density of plasma around the compact object, respectively [59]. Here, Eq. (38) can be written in the non-uniform plasma case as [63]

      ˆαSIS=ˆαSIS1+ˆαSIS2+ˆαSIS3 ,

      (44)

      where \hat{\alpha}_{\rm SIS1}, \hat{\alpha}_{\rm SIS2}, and \hat{\alpha}_{\rm SIS3} represent the deflection angle due to the gravity of the black hole, due to plasma, and due to the density of the non-uniform plasma medium around the compact object in black-bounce spacetime, respectively. From Eqs. (38), (37), and (44), we can obtain an equation for the deflection angle of light rays (photons) around the black hole in black-bounce surrounded by non-uniform plasma:

      \hat\alpha_{\rm SIS} = \frac{3 \pi a^2}{16 b^2}+\frac{2 b^3 R_s}{\left(a^2+b^2\right)^2} +\frac{R_s^2 \omega_c^2}{2 b^2 \omega^2} + \frac{2 b R_s^3 \omega_c^2 \left(5 b^2-7 a^2\right)}{15 \pi \omega^2 \left(a^2+b^2\right)^3},

      (45)

      where the new notation for the non-uniform case is

      \omega_c^2= \frac{\sigma^2_{\nu} K_e}{2 k m_p R^2_{s}}.

      (46)

      In Fig. 10, we show the plasma effect of non-uniform plasma and the parameter a on the deflection angle of photons around the black hole in black-bounce spacetime. The dependence of the deflection angle \hat{\alpha}_{\rm SIS} on non-uniform plasma and the parameter a for fixed values of the impact parameter b and the non-uniform plasma parameter is clearly shown in Fig. 11. The plots reveal that the deflection angle increases gradually with increasing plasma parameter, whereas there is a small change in the value of the deflection angle due to the influence of parameter a.

      Figure 10.  (color online) Dependence of the deflection angle \hat \alpha_{\rm SIS} on the impact parameter b for different values of the parameter a/M (left panel) and non-uniform plasma parameter (right panel).

      Figure 11.  (color online) Dependence of the deflection angle \hat\alpha_{\rm SIS} on the parameter a/M (left panel) and non-uniform plasma parameter (right panel) for a fixed value of the impact parameter b/M=30 .

      We can easily compare the two cases for analysis (uniform and non-uniform) with the help of Fig. 12. It is evident that the deflection angle of light rays (photons) in the presence of a uniform plasma medium is larger than that of the non-uniform plasma case.

      Figure 12.  (color online) Dependence of the deflection angle \hat \alpha_{b} on the impact parameter b (left panel) and parameter a/M (right panel).

    • C.   Magnification

    • Now, we can study the brightness of an image source as a consequence of gravitational lensing in the presence of uniform plasma using a lens equation in the form [64, 76]

      \begin{array}{*{20}{l}} \theta D_s=\beta D_s +\hat{\alpha}_{b} D_{ds}\ , \end{array}

      (47)

      where D_s and D_{ds} are the distances from the distant source to the observer and the lens object, respectively, and θ and β represent the angular positions of the image and source, respectively. We use the following relation between the impact parameter b and angle θ:

      \begin{array}{*{20}{l}} b=D_d\theta, \end{array}

      (48)

      where D_d is the distance from the observer to the lens object. Using these expressions, we can rewrite Eq. (47) as [64, 76]

      \beta = \theta - \frac{D_{ds}}{D_s}\frac{F(\theta)}{D_d} \frac{1}{\theta} \ ,

      (49)

      with

      F(\theta)=|\alpha_b| b=|\alpha_b (\theta)| D_d \theta.

      We can write the expression for Einstein's angle \theta_E^{pl} in the presence of uniform plasma around a black-bounce black hole using Eq. (41),

      \theta_E^{\rm pl}=\frac{D_{ds}}{D_{s}}\hat\alpha_{\rm uni}(b).

      (50)

      We numerically determine the angle by solving this equation to analyze the dependence of the magnification on the black hole parameter.

      It is well-known that the general equation to compute the magnification of an image source can be expressed in the following form [59]:

      \mu_\Sigma=\frac{I_{\rm tot}}{I_{*}}=\sum\limits_k \left| \left(\frac{\theta_k}{\beta}\right) \left(\frac{{\rm d}\theta_k}{{\rm d}\beta}\right) \right| \ , \quad k=1,2, ..., s\ ,

      (51)

      where s denotes the total number of images of the source, and I_{\rm tot} can be explained as the value of the total increased brightness due to multiple images of the source with brightness I_{*} .

      From Eq. (51), we can obtain the expression for the magnification of the image source as

      \mu_{\rm tot}^{\rm pl}=\mu_+^{\rm pl}+\mu_-^{\rm pl}=\frac{x_{\rm uni}^2+2}{x_{\rm uni}\sqrt{x_{\rm uni}^2+4}}\ ,

      (52)

      where

      x_{\rm uni}=\frac{\beta}{(\theta_E^{\rm pl})_{\rm uni}} .

      (53)

      The magnification of the image source is written as

      \left(\mu_+^{\rm pl}\right)_{\rm uni}=\frac{1}{4}\Bigg(\frac{x_{\rm uni}}{\sqrt{x_{\rm uni}^2+4}}+\frac{\sqrt{x_{\rm uni}^2+4}}{x_{\rm uni}}+2\Bigg)\ ,

      (54)

      \left(\mu_-^{\rm pl}\right)_{\rm uni}=\frac{1}{4}\Bigg(\frac{x_{\rm uni}}{\sqrt{x_{\rm uni}^2+4}}+\frac{\sqrt{x_{\rm uni}^2+4}}{x_{\rm uni}}-2\Bigg)\ .

      (55)

      Using Eq. (52), we present the dependence of the total magnification on the plasma parameter for different values of parameter a in Fig. 13. From these plots, we can conclude that the total magnification of the images increases owing to the influence of uniform plasma, and there is a slight change in the total magnification with an increase in the parameter a.

      Figure 13.  (color online) Dependence of the total magnification \mu_{\rm tot} on the parameter a for different values of the plasma parameter and on the plasma parameter for different values of parameter a.

    VI.   CONCLUSIONS
    • In this study, particle orbits around black holes in black bounce gravity are studied in detail. The results obtained can be summarized as follows:

      ● We calculate the horizon radius for various values of the black bounce parameter a. The results show that the radius of the horizon r_{h} decreases in the presence of the parameter a in black-bounce spacetime.

      ● We analyze the photon orbits around a black-bounce compact object. It is shown that orbits shift toward the central object because of black bounce gravity.

      ● We obtain and discuss the observable shadow of the black hole described by the black-bounce spacetime metric in plasma. We reveal that the shadow radius (R_{\rm sh}) of the black hole decreases with an increase in both the black-bounce parameter a and plasma frequency (Fig. 5). The tendency of the black hole shadow size to decrease with an increase in parameter a can be explained in a similar way to that of photon deflection. The effect of parameter a decreases the gravitational effect of the central object on photons, and we may observe a decrease in the size of the shadow in the black-bounce spacetime metric. The negative effect of parameter a on photon orbits may be useful in obtaining constraints using current and future observation of the black hole shadow with event horizon telescopes.

      ● We also study the effect of weak gravitational lensing around the black hole in black bounce spacetime in the presence of a plasma medium around a compact object. Here, uniform and non-uniform cases are considered and discussed in detail (see. Figs. 812).

      ● For the non-uniform case, we find that the plasma frequency depends on 1/r . The deflection angle around compact objects increases with increasing value of the parameter a. We also compare the uniform and non-uniform cases and find that the deflection angle is larger in uniform plasma than in non-uniform plasma.

      ● Additionally, we consider the magnification of the image using the deflection angle of light rays, which is clearly represented in Fig. 13.

Reference (85)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return