This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

thanks: [email protected]thanks: Corresponding author.
[email protected]

Probing a polymerized black hole with the frequency shifts of photons

Qi-Ming Fu Key Laboratory of Cosmology and Astrophysics (Liaoning Province), College of Sciences, Northeastern University, Shenyang 110819, China Institute of Physics, Shaanxi University of Technology, Hanzhong 723000, China    Xin Zhang Key Laboratory of Cosmology and Astrophysics (Liaoning Province), College of Sciences, Northeastern University, Shenyang 110819, China Key Laboratory of Data Analytics and Optimization for Smart Industry (Ministry of Education), Northeastern University, Shenyang 110819, China National Frontiers Science Center for Industrial Intelligence and Systems Optimization, Northeastern University, Shenyang 110819, China
Abstract

As is well-known, black hole plays an important role in probing the quantum effects of gravity in the strong-field regime. In this paper, we focus on a polymerized black hole and investigate the quantum effects on the redshift, blueshift and gravitational redshift of photons emitted by massive particles revolving around this black hole. With a general relativistic formalism, we obtain two elegant and concise expressions for the mass of the black hole and the quantum parameter in terms of few direct observables, and find that all the frequency shifts of photons decrease with the quantum parameter. By expanding the total frequency shift for large orbit of the emitter, low peculiar redshift, and small azimuthal angle, we find that the leading contributions to the frequency shift are originated from the kinematic shift and peculiar shift, and the quantum effects make a negligible contribution. In addition, we study the effects of the plasma on the redshift/blueshift for the first time, and conclude that the presence of plasma results in the decrease of the redshift/blueshift.

I Introduction

Although general relativity has been proven to be very successful in modern physics, it may not be the final theory of gravity because of its prediction of spacetime singularity inside a black hole or at the start of the Universe, which indicates that the theory invalidates there. It is widely believed that the singularity can be avoided by incorporating the quantum effects. Many efforts have been devoted to address these issues by resorting to the quantum theory of gravity. One of the most successful attempts is to apply an effective approach developed in loop quantum gravity to the big-bang singularity in cosmological scenarios, where the singularity is replaced by a quantum bounce bridging a contracting branch with an expanding one [1, 2]. The key point is the phase space quantization usually called polymerization [3]. This quantization procedure introduces a natural cutoff called the polymer scale, where the quantum effects become relevant close to this scale.

In recent years, similar technique has been extended to the context of black holes, and a lot of polymerized black hole models have been constructed [4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. Most of these models focused on the spherically symmetric cases. The common feature of these types of models is that the singularity inside a black hole is replaced by a spacelike transition hypersurface, which smoothly connects a trapped region (the black hole region) to an antitrapped region (the white hole region), and the spacetime is regular everywhere. Besides, the resulting spacetime with the quantum corrections can still be described by an explicit metric. Thus, one can investigate such nonsingular models by exploiting the already well-developed mechanism for the semiclassical black holes.

In the last couple of years, various aspects of the polymerized black holes have been studied in the literature. For example, the role of Dirac observables in polymerized black hole models was investigated in Ref. [19], and the interior structure, quasinormal spectra of various perturbation modes were studied in Ref. [20]. The thermodynamic properties of the polymerized black holes were analyzed in Ref. [21], and the effects of the quantum corrections on the relevant thermodynamic quantities also were considered. In addition, the Hawking radiation emitted by polymerized black holes for the fields with different spins were investigated in Refs. [22, 23], and the authors also investigated the primordial black hole evaporation signals with the Hawking radiation rates. Recently, the strong gravitational lensing of the polymerized black holes was investigated in Ref. [24], and the solar system experiments of these black holes were considered in Ref. [25], from which the constraints on the polymer scale were derived as well. Other phenomenological investigations and observational constraints on the polymerized black holes were explored in Refs. [26, 27, 28, 29, 30, 31, 32, 33, 34] and references therein. In this paper, we mainly concentrate on a particular polymerized black hole introduced in Refs. [18, 19] and investigate the frequency shifts of the photons emitted by massive objects orbiting this black hole.

On the other hand, the authors of Ref. [35] introduced a general relativistic method that can express the mass and spin of a Kerr black hole in terms of direct observables, i.e., the radius of the massive particles orbiting the black hole and the frequency shifts experienced by the photons emitted by these revolving particles. Subsequently, this general relativistic mechanism was applied to other black holes, such as, Myers-Perry black hole [36], Kerr-Newman-(anti) de Sitter black hole [37], and Plebanski-Demianski black hole [38]. Besides, the redshift and blueshift of photons emitted by timelike geodesic particles near compact objects, such as, boson stars, as well as the Schwarzschild and Reissner-Nordstro¨\ddot{\text{o}}m black holes, were investigated in Ref. [39]. In Refs. [40, 41], a similar formalism was used to obtain the relations between the parameters describing a class of regular black holes and the redshift observables. Other investigations on the frequency shifts of photons can be found in Refs. [42, 43, 44, 45, 46, 47].

Although the polymerized black hole has been extensively studied in the literature, the frequency shifts of photons emitted by massive particles orbiting this black hole have not been analyzed yet. With the general relativistic formalism, one can obtain closed formulas for the mass and the quantum parameter of the polymerized black hole in terms of the frequency shifts of the emitted photons. Thus, this mechanism not only provides us with valuable information about the black hole with few direct observables, but also can be used as a powerful tool to probe the quantum effects. On the other hand, most astrophysical black holes in nature are surrounded by a plasma medium and the frequency of photons will be influenced by the plasma [48, 49, 50]. Thus, it is necessary to study the polymerized black hole with the general relativistic formalism and it is also of interest to know what the effects are of the plasma on the frequency shifts.

In this paper, we focus on exploring the quantum effects through the frequency shifts of photons for a particular polymerized black hole model, and give a brief analysis about the effects of the plasma on the frequency shifts. This paper is organized as follows: In Sec. II, we first give a brief review of the polymerized black hole and then derive the equations of motion for the orbiting massive particles and the photons emitted by these particles. In Sec. III, we investigate the frequency shifts of the photons and derive the closed formulas for the mass and quantum parameter of the black hole. The effects of the plasma on the frequency shifts are analyzed in Sec. IV. Section V contains the conclusion.

II Particles and photons in a polymerized black hole spacetime

The static and spherically-symmetric metric modified by the effective loop quantum gravity can be expressed as [18, 19, 26],

ds2=𝒜(y)dτ2+dy2𝒜(y)+𝒞2(y)dΩ22,\displaystyle ds^{2}=-\mathcal{A}(y)d\tau^{2}+\frac{dy^{2}}{\mathcal{A}(y)}+\mathcal{C}^{2}(y)d\Omega_{2}^{2},~{} (1)

where the radial coordinate y(,)y\in(-\infty,\infty), and the metric functions are given by

𝒜(y)\displaystyle\mathcal{A}(y) =\displaystyle= (112Aλ(1+y2))1+y2𝒞2(y),\displaystyle\left(1-\frac{1}{\sqrt{2A_{\lambda}(1+y^{2})}}\right)\frac{1+y^{2}}{\mathcal{C}^{2}(y)}, (2)
𝒞2(y)\displaystyle\mathcal{C}^{2}(y) =\displaystyle= Aλ1+y2MB2(y+1+y2)6+MW2(y+1+y2)3.\displaystyle\frac{A_{\lambda}}{\sqrt{1+y^{2}}}\frac{M_{B}^{2}(y+\sqrt{1+y^{2}})^{6}+M_{W}^{2}}{(y+\sqrt{1+y^{2}})^{3}}. (3)

Here MBM_{B} and MWM_{W} represent the masses of an asymptotically Schwarzschild black hole and a white hole, respectively [26], and AλA_{\lambda} is a dimensionless and non-negative parameter defined by Aλ(λk/(MBMW))2/3/2A_{\lambda}\equiv(\lambda_{k}/(M_{B}M_{W}))^{2/3}/2, where λk\lambda_{k} is a quantum parameter originated from holonomy modifications and should be thought of being of Planckian order [18, 19]. Besides, according to its definition, the parameter AλA_{\lambda} is extremely small for astrophysical black holes due to the large suppression of the black hole mass. In this paper, we focus on the case of MB=MW=MM_{B}=M_{W}=M corresponding to a symmetric bounce, i.e., the spacetime is symmetric under yyy\rightarrow-y. Without loss of generality, we focus on the positive branch with y0y\geqslant 0. By redefining two new coordinates r8AλMyr\equiv\sqrt{8A_{\lambda}}My and tτ/(8AλM)t\equiv\tau/(\sqrt{8A_{\lambda}}M), the metric (1) can be rexpressed as [26, 24]

ds2=A(r)dt2+B(r)dr2+C(r)(dθ2+sin2θdϕ2),\displaystyle ds^{2}=-A(r)dt^{2}+B(r)dr^{2}+C(r)(d\theta^{2}+\sin^{2}\theta d\phi^{2}),~{} (4)

where

A(r)\displaystyle A(r)\!\! =\displaystyle= 1B(r)=8AλM2+r2(8AλM2+r22M)2AλM2+r2,\displaystyle\!\!\frac{1}{B(r)}\!=\!\frac{\sqrt{8A_{\lambda}M^{2}+r^{2}}\big{(}\sqrt{8A_{\lambda}M^{2}+r^{2}}-2M\big{)}}{2A_{\lambda}M^{2}+r^{2}},~{}~{}~{}~{} (5)
C(r)\displaystyle C(r)\!\! =\displaystyle= 2AλM2+r2.\displaystyle\!\!2A_{\lambda}M^{2}+r^{2}~{}. (6)

The event horizon of this black hole is determined by A(r)=0A(r)=0, which admits the solution r=2M12Aλr=2M\sqrt{1-2A_{\lambda}}. Hence the permissible range for the quantum parameter AλA_{\lambda} in this black hole spacetime is 0Aλ120\leqslant A_{\lambda}\leqslant\frac{1}{2}.

The Lagrangian for a massive particle moving around this black hole is

2=A(r)t˙2+B(r)r˙2+C(r)(θ˙2+sin2θϕ˙2),\displaystyle 2\mathcal{L}=-A(r)\dot{t}^{2}+B(r)\dot{r}^{2}+C(r)(\dot{\theta}^{2}+\sin^{2}\theta\dot{\phi}^{2}),~{} (7)

where x˙μ=Uμdxμdτ\dot{x}^{\mu}=U^{\mu}\equiv\frac{dx^{\mu}}{d\tau}, UμU^{\mu} denotes the four velocity of the massive particle, and τ\tau is an affine parameter along the geodesic. Since the Lagrangian only depends on rr and θ\theta, there appear two conserved quantities along the geodesic,

pt\displaystyle p_{t} =\displaystyle= t˙=A(r)t˙=A(r)Ut=E,\displaystyle\frac{\partial\mathcal{L}}{\partial\dot{t}}=-A(r)\dot{t}=-A(r)U^{t}=-E,~{} (8)
pϕ\displaystyle p_{\phi} =\displaystyle= ϕ˙=C(r)sin2θϕ˙=C(r)sin2θUϕ=L,\displaystyle\frac{\partial\mathcal{L}}{\partial\dot{\phi}}=C(r)\sin^{2}\theta\dot{\phi}=C(r)\sin^{2}\theta U^{\phi}=L,~{} (9)

which correspond to the energy and angular momentum of the massive particle with unit rest mass, respectively. From the normalized condition UμUμ=1U_{\mu}U^{\mu}=-1, the geodesic equation for the massive particle at the equatorial plane (Uθ=0U^{\theta}=0) can be derived as

(Ur)2+V(r)=0,\displaystyle(U^{r})^{2}+V(r)=0, (10)

with the effective potential

V(r)=1B(1E2A+L2C).\displaystyle V(r)=\frac{1}{B}\left(1-\frac{E^{2}}{A}+\frac{L^{2}}{C}\right). (11)

For circular orbits, i.e., Ur=0U^{r}=0, the effective potential must have an extremum indicating two conditions V(r)=0V(r)=0 and V(r)=0V^{\prime}(r)=0, where the prime stands for the derivative with respect to rr. From these two conditions, one can obtain

E2\displaystyle E^{2} =\displaystyle= A2CACCA,\displaystyle\frac{A^{2}C^{\prime}}{AC^{\prime}-CA^{\prime}},~{} (12)
L2\displaystyle L^{2} =\displaystyle= C2AACCA.\displaystyle\frac{C^{2}A^{\prime}}{AC^{\prime}-CA^{\prime}}.~{} (13)

Besides, for the circular orbits to be stable, the following stable condition should also be satisfied:

V′′=C(2AAC+C(AA′′2A2))ACAC′′ABC(ACCA)>0.\displaystyle~{}V^{\prime\prime}\!\!=\!\!\frac{C^{\prime}\!\left(2AA^{\prime}C^{\prime}\!+\!C\left(AA^{\prime\prime}\!-\!2A^{\prime 2}\right)\right)\!-\!ACA^{\prime}C^{\prime\prime}}{ABC\left(AC^{\prime}-CA^{\prime}\right)}>0. (14)

Inserting Eqs. (12) and (13) into Eqs. (8) and (9), the four-velocity for the circular and equatorial orbiting massive particles can be calculated as

Ut\displaystyle U^{t} =\displaystyle= CACCA,\displaystyle\sqrt{\frac{C^{\prime}}{AC^{\prime}-CA^{\prime}}},~{} (15)
Uϕ\displaystyle U^{\phi} =\displaystyle= AACCA.\displaystyle\sqrt{\frac{A^{\prime}}{AC^{\prime}-CA^{\prime}}}.~{} (16)

On the other hand, the energy and angular momentum for the photons moving along null geodesics can be obtained in a similar way,

Eγ=Akt,Lγ=Ckϕ,\displaystyle E_{\gamma}=Ak^{t},\quad L_{\gamma}=Ck^{\phi},~{} (17)

with kμ=(kt,kr,kθ,kϕ)k^{\mu}=(k^{t},k^{r},k^{\theta},k^{\phi}) the four momentum of the photons, which satisfies kμkμ=0k^{\mu}k_{\mu}=0. Then, the equation of motion for the photons at the equatorial plane is

A(kt)2+B(kr)2+C(kϕ)2=0.\displaystyle-A(k^{t})^{2}+B(k^{r})^{2}+C(k^{\phi})^{2}=0.~{} (18)

After inserting Eq. (17) into the above equation, the radial component of the four momentum can be solved as

(kr)2=ALγ2+CEγ2ABC.\displaystyle~{}(k^{r})^{2}=\frac{-AL_{\gamma}^{2}+CE_{\gamma}^{2}}{ABC}. (19)

Introducing an auxiliary bidimensional and geometrical vector defined by κ2(kr)2+C(kϕ)2\kappa^{2}\equiv(k^{r})^{2}+C(k^{\phi})^{2}, where the components krk^{r} and kϕk^{\phi} satisfy the following decomposition [47]:

kr\displaystyle k^{r} =\displaystyle= κcosϕ,\displaystyle\kappa\cos\phi,~{} (20)
Ckϕ\displaystyle\sqrt{C}k^{\phi} =\displaystyle= κsinϕ,\displaystyle\kappa\sin\phi, (21)

with 0ϕ2π0\leqslant\phi\leqslant 2\pi, the auxiliary vector can be explicitly expressed as

κ2=ALγ2+CEγ2ABC+C(LγC)2.\displaystyle~{}\kappa^{2}=\frac{-AL_{\gamma}^{2}+CE_{\gamma}^{2}}{ABC}+C\left(\frac{L_{\gamma}}{C}\right)^{2}. (22)

Besides, after inserting Eq. (20) into Eq. (19), another expression for the auxiliary vector can be obtained as

κ2=ALγ2+CEγ2ABCcos2ϕ.\displaystyle~{}\kappa^{2}=\frac{-AL_{\gamma}^{2}+CE_{\gamma}^{2}}{ABC\cos^{2}\phi}. (23)

Combining these two expressions (22) and (23) for the auxiliary vector together, one can obtain an equation for the impact parameter bEγLγb\equiv\frac{E_{\gamma}}{L_{\gamma}},

b2A(sin2ϕ+Bcos2ϕ)Csin2ϕ=0,\displaystyle b^{2}A(\sin^{2}\phi+B\cos^{2}\phi)-C\sin^{2}\phi=0, (24)

which admits the following solution

b=CsinϕA(sin2ϕ+Bcos2ϕ).\displaystyle b=-\frac{\sqrt{C}\sin\phi}{\sqrt{A(\sin^{2}\phi+B\cos^{2}\phi)}}.~{} (25)

The maximum or minimum of the impact parameter is reached when kr=0k^{r}=0 corresponding to ϕ=π/2\phi=-\pi/2 or π/2\pi/2, respectively.

III Redshift/Blueshift of emitted photons

The frequency shift of a photon emitted by a massive particle moving around a black hole and then received by a detector can be defined as [35]

1+z=ωeωd=kμUμ|ekμUμ|d,\displaystyle~{}1+z=\frac{\omega_{\text{e}}}{\omega_{\text{d}}}=\frac{-k_{\mu}U^{\mu}|_{\text{e}}}{-k_{\mu}U^{\mu}|_{\text{d}}}, (26)

where the lower indexes e and d represent the emitter and detector, respectively. For the case of a circular and equatorial orbiting massive particle, i.e., Ur=Uθ=0U^{r}=U^{\theta}=0, Eq. (26) can be reduced to

1+z=(EγUtLγUϕ)|e(EγUtLγUϕ)|d=UetbeUeϕUdtbdUdϕ.\displaystyle~{}1+z=\frac{(E_{\gamma}U^{t}-L_{\gamma}U^{\phi})|_{\text{e}}}{(E_{\gamma}U^{t}-L_{\gamma}U^{\phi})|_{\text{d}}}=\frac{U^{t}_{\text{e}}-b_{\text{e}}U^{\phi}_{\text{e}}}{U^{t}_{\text{d}}-b_{\text{d}}U^{\phi}_{\text{d}}}. (27)

There are mainly two contributions to the frequency shift defined above, i.e., the gravitational effect and the kinematic effect or Doppler effect. The impact parameter for a radially emitted photon is usually called the central impact parameter bcb_{\text{c}}, which vanishes for a spherically-symmetric and static black hole, and the only contribution to the frequency shift of this photon is the gravitational effect. Thus, the gravitational redshift can be conveniently defined as [39]

1+zg=UetbcUeϕUdtbcUdϕ=UetUdt.\displaystyle~{}1+z_{\text{g}}=\frac{U^{t}_{\text{e}}-b_{\text{c}}U^{\phi}_{\text{e}}}{U^{t}_{\text{d}}-b_{\text{c}}U^{\phi}_{\text{d}}}=\frac{U^{t}_{\text{e}}}{U^{t}_{\text{d}}}. (28)

Then, the kinematic redshift zkinz_{\text{kin}} originated from the kinematic effect of the emitter can be defined by subtracting the gravitational redshift zgz_{\text{g}} from the frequency shift zz [39],

zkinzzg=UetbeUeϕUdtbdUdϕUetUdt.\displaystyle~{}z_{\text{kin}}\equiv z-z_{\text{g}}=\frac{U^{t}_{\text{e}}-b_{\text{e}}U^{\phi}_{\text{e}}}{U^{t}_{\text{d}}-b_{\text{d}}U^{\phi}_{\text{d}}}-\frac{U^{t}_{\text{e}}}{U^{t}_{\text{d}}}. (29)

Besides, if there is a relative motion of the black hole with respect to a distant observer or detector, there should be an additional correction to the frequency shifts of photons coming from the special relativistic boost, which is defined by [47]

1+zboost=γ(1+β),γ=(1β2)1/2,β=v0c,\displaystyle 1+z_{\text{boost}}=\gamma(1+\beta),\quad\gamma=(1-\beta^{2})^{-1/2},\quad\beta=\frac{v_{0}}{c},~{}~{} (30)

where v0=z0cv_{0}=z_{0}c represents the radial peculiar velocity of the black hole with respect to the observer, and z0z_{0} is called the peculiar redshift which codifies the motion of the black hole receding (positive z0z_{0}) or approaching (negative z0z_{0}) with respect to the observer. Then, by considering Eqs. (15), (16) and (25), the total frequency shift can be expressed as [47]

ztot\displaystyle~{}z_{\text{tot}} =\displaystyle= (1+z)(1+zboost)1\displaystyle(1+z)(1+z_{\text{boost}})-1 (31)
=\displaystyle= (1+zkin+zg)γ(1+β)1\displaystyle(1+z_{\text{kin}}+z_{\text{g}})\gamma(1+\beta)-1
=\displaystyle= 1+z01z0(C(re)A(re)C(re)C(re)A(re)\displaystyle\sqrt{\frac{1+z_{0}}{1-z_{0}}}\bigg{(}\sqrt{\frac{C^{\prime}(r_{\text{e}})}{A(r_{\text{e}})C^{\prime}(r_{\text{e}})-C(r_{\text{e}})A^{\prime}(r_{\text{e}})}}
+\displaystyle+ sinϕC(re)A(re)A(re)C(re)C(re)A(re)A(re)(sin2ϕ+B(re)cos2ϕ))1,\displaystyle\frac{\sin\phi\sqrt{\frac{C(r_{\text{e}})A^{\prime}(r_{\text{e}})}{A(r_{\text{e}})C^{\prime}(r_{\text{e}})-C(r_{\text{e}})A^{\prime}(r_{\text{e}})}}}{\sqrt{A(r_{\text{e}})(\sin^{2}\phi+B(r_{\text{e}})\cos^{2}\phi)}}\bigg{)}-1,

where rer_{\text{e}} is the orbital radius of the revolving massive particles, and the static observer is assumed to be located at infinity with the four velocity Uμ|d=(1,0,0,0)U^{\mu}|_{\text{d}}=(1,0,0,0).

By substituting the metric functions (5) and (6) into Eq. (31), one can obtain the total frequency shift for this particular polymerized black hole, which is not given here due to its length. However, for large orbit of the emitter MreM\ll r_{\text{e}}, low peculiar redshift (z01z_{0}\ll 1), and small azimuthal angle ϕ0\phi\approx 0, the total frequency shift can be expanded as

ztot\displaystyle z_{\text{tot}} \displaystyle\approx 32M~+z0+ϕM~1/2+27M~28+3M~z02+12z02\displaystyle\frac{3}{2}\tilde{M}+z_{0}+\phi\tilde{M}^{1/2}+\frac{27\tilde{M}^{2}}{8}+\frac{3\tilde{M}z_{0}}{2}+\frac{1}{2}z_{0}^{2} (32)
+\displaystyle+ 32M~3/2ϕ+M~1/2z0ϕ6AλM~2,\displaystyle\frac{3}{2}\tilde{M}^{3/2}\phi+\tilde{M}^{1/2}z_{0}\phi-6A_{\lambda}\tilde{M}^{2},

with M~M/re\tilde{M}\equiv M/r_{\text{e}}. Obviously, the first two terms are the leading contribution to the frequency shift. The first term is originated from the so-called Doppler or kinematic shift, i.e., the motion of the emitter revoloving around the black hole. The second term is originated from the peculiar motion, i.e., the emitter and the black hole as a whole entity approaching or receding from us. The third term is the subleading contribution to the frequency shift which is also originated from the kinematic shift. The fifth and eighth terms in this expansion are originated from the combination of the kinematic shift and peculiar redshift. The parameter AλA_{\lambda} first appears in the ninth term of this series expansion, which encodes the effects of the quantum corrections and indicates a very subtle contribution to the total frequency shift. In addition, AλA_{\lambda} is extremely small for astrophysical black holes also signifying that the quantum effects on the frequency shift can be negligible.

Before proceeding, we should note that there exist three special azimuthal angles, i.e., ϕ=0,±π2\phi=0,\pm\frac{\pi}{2}. Here ϕ=0\phi=0 is assumed to be the opposite direction of the line of sight of the observer, and the observer is on the right-hand side of the black hole. Obviously, the impact parameter vanishes in this case and describes a photon emitted by an anticlockwise orbiting particle when moving at the location with the azimuthal angle ϕ=0\phi=0. At this moment, the three velocity of the orbiting particle is perpendicular with respect to the line of sight of the distant observer, and the only contribution to the frequency shift of photons is the gravitational effect, which is defined by ztot3ztot|ϕ0z_{\text{tot}_{3}}\equiv z_{\text{tot}}|_{\phi\rightarrow 0}. Then, Eq. (25) shows that ϕ=±π2\phi=\pm\frac{\pi}{2} corresponds to the minimum and maximum of the impact parameter, respectively. The minimal impact parameter describes the photon emitted by the anticlockwise orbiting particle when arriving at the location with the azimuthal angle ϕ=π2\phi=\frac{\pi}{2}. At this moment, the initial three velocity of the photon is exactly opposite to the instantaneous orbital velocity of the particle. Thus, the redshift of the photon reaches the maximum defined by ztot1ztot|ϕπ/2z_{\text{tot}_{1}}\equiv z_{\text{tot}}|_{\phi\rightarrow\pi/2}. On the other hand, the maximal impact parameter describes the photon emitted by the revolving particle when reaching at the location with the azimuthal angle ϕ=π2\phi=-\frac{\pi}{2} and the initial three-velocity of the photon is the same as the instantaneous orbital velocity of the particle, which corresponds to the case of the maximal blueshift of the photon defined by ztot2ztot|ϕπ/2z_{\text{tot}_{2}}\equiv z_{\text{tot}}|_{\phi\rightarrow-\pi/2}.

Then, by defining the following variables,

R\displaystyle R \displaystyle\equiv 1+ztot1,\displaystyle 1+z_{\text{tot}_{1}}, (33a)
S\displaystyle S \displaystyle\equiv 1+ztot2,\displaystyle 1+z_{\text{tot}_{2}}, (33b)
T\displaystyle T \displaystyle\equiv (1+ztot3)2,\displaystyle\left(1+z_{\text{tot}_{3}}\right)^{2}, (33c)

one can obtain two elegant closed formulas for the mass of the black hole and the quantum parameter in terms of these variables by inverting the relations (33):

M\displaystyle M\!\! =\displaystyle= re(1+ζT2ζRS)3ζT15ζT4RS(2ζ+T),\displaystyle\!\!r_{\text{e}}\left(1+\frac{\zeta}{T}-\frac{2\zeta}{RS}\right)\sqrt{\frac{3\zeta T}{15\zeta T-4RS(2\zeta+T)}},~{} (34)
Aλ\displaystyle A_{\lambda}\!\! =\displaystyle= TR2S2(RS(2ζ+T)3ζT)6ζ(RS(ζ+T)2ζT)2,\displaystyle\!\!\frac{TR^{2}S^{2}(RS(2\zeta+T)-3\zeta T)}{6\zeta(RS(\zeta+T)-2\zeta T)^{2}}, (35)

with ζ1+z01z0\zeta\equiv\frac{1+z_{0}}{1-z_{0}}. The function MM given by Eq. (34) and the radius of the emitter rer_{\text{e}} are expressed in geometrized units (G=c=1G=c=1) and scaled by pMpM_{\odot}, where MM_{\odot} is the solar mass and pp is an arbitrary factor of proportionality, such as, p=4.3×106p=4.3\times 10^{6} for SgrA*. These two formulas tell us that the mass of the black hole and the quantum parameter can be completely and independently given by few direct observables, i.e., the redshift ztot1z_{\text{tot}_{1}}, blueshift ztot2z_{\text{tot}_{2}}, and gravitational redshift ztot3z_{\text{tot}_{3}} of photons, as well as the orbital radius rer_{\text{e}} of the massive particle. It should be noted that the peculiar redshift z0z_{0} is not a measurable quantity, but can be statistically estimated with the help of the direct observables and Eq. (31). Obviously, this general relativistic formalism described above can be used to predict the property of the compact object and probe the quantum effects precisely. To be specific, the compact object is a Schwarzschild black hole for Aλ=0A_{\lambda}=0, a regular black hole for 0<Aλ<120<A_{\lambda}<\frac{1}{2}, and a wormhole for Aλ=12A_{\lambda}=\frac{1}{2}. One can see Ref. [24] for more details.

In the following, we will give a brief analysis about the influences of the quantum effects and peculiar velocity on the frequency shifts of photons. Figure 1 shows that both of the redshift ztot1z_{\text{tot}_{1}} and the gravitational redshift ztot3z_{\text{tot}_{3}} decrease with AλA_{\lambda}. Besides, the blueshift also decreases with AλA_{\lambda} although the value of ztot2z_{\text{tot}_{2}} increases with it, since ωd\omega_{\text{d}} is in the denominator and the smaller ωd\omega_{\text{d}} means the larger ztot2z_{\text{tot}_{2}} but the smaller blueshift. In a word, the quantum effects reduce the frequency shifts of photons. The possible reason is that the quantum effects usually result in a weakening of the gravitational force [51]. What’s more, Fig. 1 shows that the frequency shifts of the photons also decrease with the orbital radius rer_{\text{e}} for fixed AλA_{\lambda} because the larger rer_{\text{e}} is, the weaker the gravitational field becomes. Finally, Fig. 1 shows the frequency shifts versus the peculiar redshift or peculiar velocity. Obviously, the redshift ztot1z_{\text{tot}_{1}} and gravitational redshift ztot3z_{\text{tot}_{3}} increase with the radial receding peculiar velocity (corresponding to positive z0z_{0}) of the black hole with respect to the detector while the blueshift decreases with it although the value of ztot2z_{\text{tot}_{2}} increases. Besides, there exists a zero point for the blueshift since the kinematic effects is counteracted by the peculiar redshift at a particular receding velocity and similar zero point also respectively appears for the redshift and gravitational redshift at a particular approaching velocity (corresponding to negative z0z_{0}).

Refer to caption
Refer to caption
Refer to caption
Figure 1: The top panel plots the frequency shifts versus the quantum parameter AλA_{\lambda} with other dimensionless parameters set to M~=1/7\tilde{M}=1/7 and z0=0z_{0}=0. The middle panel shows the frequency shifts versus the dimensionless orbiting radius re/Mr_{\text{e}}/M with z0=0z_{0}=0 and Aλ=1/3A_{\lambda}=1/3. The bottom panel shows the frequency shifts versus the peculiar redshift z0z_{0} with M~=1/11\tilde{M}=1/11 and Aλ=1/3A_{\lambda}=1/3. It should be noted that the parameter AλA_{\lambda} can be the order of unity only for black holes with Planck mass scale, and is extremely small for astrophysical black holes. Besides, all parameters are set to ensure E2>0E^{2}>0, L2>0L^{2}>0 and the stable condition (14) and rer_{\text{e}} is larger than the radius of the innermost stable circular orbit (ISCO).

IV The presence of plasma

In this section, we consider a more practical scenario, i.e., the black hole surrounded by a nonmagnetized, cold plasma, and investigate the effects of the plasma on the redshift/blueshift. The plasma is assumed as a static, inhomogeneous medium and the Hamiltonian for the photon in this plasma is given by [48, 49]

H=12[gμνkμkν+ωp(r)2],\displaystyle H=\frac{1}{2}\big{[}g^{\mu\nu}k_{\mu}k_{\nu}+\omega_{\text{p}}(r)^{2}\big{]}, (36)

where ωp\omega_{\text{p}} stands for the electron plasma frequency. By introducing the refractive index of the plasma medium

n2=1ωp2ω2,\displaystyle n^{2}=1-\frac{\omega_{\text{p}}^{2}}{\omega^{2}}, (37)

where ω\omega stands for the frequency of photons measured by a static observer with respect to the plasma, the Hamiltonian can be rewritten as

H=12[gμνkμkν(n21)(kμVμ)2]=12geffμνkμkν,\displaystyle~{}H=\frac{1}{2}\big{[}g^{\mu\nu}k_{\mu}k_{\nu}-(n^{2}-1)(k_{\mu}V^{\mu})^{2}\big{]}=\frac{1}{2}g^{\mu\nu}_{\text{eff}}k_{\mu}k_{\nu},~{}~{}~{} (38)

where VμV^{\mu} denotes the four velocity of the plasma, and geffμνgμν(n21)VμVνg^{\mu\nu}_{\text{eff}}\equiv g^{\mu\nu}-(n^{2}-1)V^{\mu}V^{\nu}. Equation (38) shows that photons travel along the geodesic in the effective metric geffμνg^{\mu\nu}_{\text{eff}}. Thus, for the case of a circular and equatorial orbiting massive particle, Eq. (26) can still be formally expressed as

1+z=(EγUtLγUϕ)|e(EγUtLγUϕ)|d=Uetb^eUeϕUdtb^dUdϕ,\displaystyle~{}1+z=\frac{(E_{\gamma}U^{t}-L_{\gamma}U^{\phi})|_{\text{e}}}{(E_{\gamma}U^{t}-L_{\gamma}U^{\phi})|_{\text{d}}}=\frac{U^{t}_{\text{e}}-\hat{b}_{\text{e}}U^{\phi}_{\text{e}}}{U^{t}_{\text{d}}-\hat{b}_{\text{d}}U^{\phi}_{\text{d}}}, (39)

with b^LγEγ\hat{b}\equiv\frac{L_{\gamma}}{E_{\gamma}}. Although the definition is the same as the case of the absence of plasma, the impact parameter now incorporates the effects of plasma as explained in the following.

Considering the following plasma frequency with a radial power-law number density [48, 49]

ωp(r)2=4πe2mN(r),N(r)=N0rh,\displaystyle\omega_{\text{p}}(r)^{2}=\frac{4\pi e^{2}}{m}N(r),\quad N(r)=\frac{N_{0}}{r^{h}}, (40)

where N0N_{0} is a constant, and ee and mm respectively represent the electron charge and mass, the refractive index of the plasma medium can be expressed as

n2=1A(r)krh,\displaystyle n^{2}=1-A(r)\frac{k}{r^{h}}, (41)

with k4πe2mω02N0k\equiv\frac{4\pi e^{2}}{m\omega_{0}^{2}}N_{0} and ω0\omega_{0} the frequency of the photon measured by a detector at infinity. For simplicity, we focus on the case of h=1h=1. Within this effective geometry, the energy and angular momentum for photons moving along the geodesics can be obtained as Eγ=An2ktE_{\gamma}=\frac{A}{n^{2}}k^{t} and Lγ=CkϕL_{\gamma}=Ck^{\phi}. Then, the equation of motion for photons at the equatorial plane can be given by

n2AEγ2+B(kr)2+Lγ2C=0.\displaystyle-\frac{n^{2}}{A}E_{\gamma}^{2}+B(k^{r})^{2}+\frac{L_{\gamma}^{2}}{C}=0. (42)

Following the same procedure in the last section, the impact parameter can be calculated as

b^=CnsinϕA(sin2ϕ+Bcos2ϕ).\displaystyle~{}\hat{b}=-\frac{\sqrt{C}n\sin\phi}{\sqrt{A(\sin^{2}\phi+B\cos^{2}\phi)}}. (43)

Obviously, the impact parameter incorporates the refractive index nn of the plasma medium. Inserting the impact parameter into Eq. (39), the total frequency shift with the presence of plasma can be expressed as

z^tot\displaystyle\hat{z}_{\text{tot}} =\displaystyle= 1+z01z0(C(re)A(re)C(re)C(re)A(re)\displaystyle\sqrt{\frac{1+z_{0}}{1-z_{0}}}\bigg{(}\sqrt{\frac{C^{\prime}(r_{\text{e}})}{A(r_{\text{e}})C^{\prime}(r_{\text{e}})-C(r_{\text{e}})A^{\prime}(r_{\text{e}})}} (44)
+\displaystyle+ n(re)sinϕC(re)A(re)A(re)C(re)C(re)A(re)A(re)(sin2ϕ+B(re)cos2ϕ))1,\displaystyle\frac{n(r_{\text{e}})\sin\phi\sqrt{\frac{C(r_{\text{e}})A^{\prime}(r_{\text{e}})}{A(r_{\text{e}})C^{\prime}(r_{\text{e}})-C(r_{\text{e}})A^{\prime}(r_{\text{e}})}}}{\sqrt{A(r_{\text{e}})(\sin^{2}\phi+B(r_{\text{e}})\cos^{2}\phi)}}\bigg{)}-1,

where the static observer is also assumed to be located at infinity with the four-velocity Uμ|d=(1,0,0,0)U^{\mu}|_{\text{d}}=(1,0,0,0).

Similarly, one can also define three variables

R^\displaystyle\hat{R} \displaystyle\equiv 1+z^tot1,\displaystyle 1+\hat{z}_{\text{tot}_{1}}, (45a)
S^\displaystyle\hat{S} \displaystyle\equiv 1+z^tot2,\displaystyle 1+\hat{z}_{\text{tot}_{2}}, (45b)
T^\displaystyle\hat{T} \displaystyle\equiv (1+z^tot3)2,\displaystyle\left(1+\hat{z}_{\text{tot}_{3}}\right)^{2}, (45c)

where z^tot1z^tot|ϕπ/2\hat{z}_{\text{tot}_{1}}\equiv\hat{z}_{\text{tot}}|_{\phi\rightarrow\pi/2}, z^tot2z^tot|ϕπ/2\hat{z}_{\text{tot}_{2}}\equiv\hat{z}_{\text{tot}}|_{\phi\rightarrow-\pi/2}, and z^tot3z^tot|ϕ0\hat{z}_{\text{tot}_{3}}\equiv\hat{z}_{\text{tot}}|_{\phi\rightarrow 0}. Then, two elegant closed formulas for the mass of the black hole and the quantum parameter can be derived from these variables by inverting the relations Eq. (45):

M\displaystyle M =\displaystyle= nre(T^n2(T^ζ)+(ζ+T^)(R^S^T^))T^(R^S^T^(1n2))3ζT^T^n2(7ζ4T^)4(2ζ+T^)(R^S^T^),\displaystyle\frac{nr_{\text{e}}\left(\hat{T}n^{2}(\hat{T}-\zeta)+(\zeta+\hat{T})(\hat{R}\hat{S}-\hat{T})\right)}{\hat{T}\left(\hat{R}\hat{S}-\hat{T}\left(1-n^{2}\right)\right)}\sqrt{\frac{3\zeta\hat{T}}{\hat{T}n^{2}(7\zeta-4\hat{T})-4(2\zeta+\hat{T})(\hat{R}\hat{S}-\hat{T})}}, (46)
Aλ\displaystyle A_{\lambda} =\displaystyle= T^(R^S^T^(1n2))2(T^n2(T^ζ)+(2ζ+T^)(R^S^T^))6ζ(T^n3(T^ζ)+n(ζ+T^)(R^S^T^))2.\displaystyle\frac{\hat{T}\left(\hat{R}\hat{S}-\hat{T}\left(1-n^{2}\right)\right)^{2}\left(\hat{T}n^{2}(\hat{T}-\zeta)+(2\zeta+\hat{T})(\hat{R}\hat{S}-\hat{T})\right)}{6\zeta\left(\hat{T}n^{3}(\hat{T}-\zeta)+n(\zeta+\hat{T})(\hat{R}\hat{S}-\hat{T})\right)^{2}}. (47)

In the following, we make a brief analysis about the effects of the plasma on the redshift/blueshift. To make the analysis clearly, we assume fixed ω0\omega_{0} or ωd\omega_{\text{d}}, which is the natural frequency of the photon measured by a static observer located at infinity, and investigate the effects of the plasma on the frequency shifts through the emitting frequency ωe\omega_{\text{e}}. Figure 2 shows that both of the redshift and blueshift decrease with kk, and there is a bound for the range of kk since the photons cannot propagate in the region with ω0<ωp(r)A(r)\omega_{0}<\omega_{\text{p}}(r)\sqrt{A(r)}. Besides, it should be noted that the plasma has no effects on z^tot3\hat{z}_{\text{tot}_{3}}. The reason is that the impact parameter vanishes when the photon propagate radially along the line of sight of the observer, and Eq. (39) shows that the influence of the plasma vanishes when b^0\hat{b}\rightarrow 0.

Refer to caption
Refer to caption
Figure 2: The redshift (top panel) and blueshift (bottom panel) versus the parameter kk with vanishing peculiar redshift z0=0z_{0}=0. The other dimensionless parameters are set to M~=1/10\tilde{M}=1/10 and Aλ=1/3A_{\lambda}=1/3. All parameters satisfy E2>0E^{2}>0, L2>0L^{2}>0 and the stable condition (14).

V Conclusion

In this paper, we investigated the frequency shifts of photons emitted by massive particles revolving around a particular polymerized black hole introduced in Refs. [18, 19], which is characterized by two parameters, the black hole mass MM and the non-negative parameter AλA_{\lambda} originated from quantum corrections. It is found that these two parameters can be completely and independently given by few direct astrophysical observables related to the frequency shifts of the photons. Then, we gave a brief analysis about the influences of the quantum effects on the frequency shifts. Besides, since plasma will change the behavior of the photon and most real black holes are surrounded by plasma, we studied the effects of the plasma on the frequency shifts of photons.

To be specific, for the case of a circular and equatorial orbiting massive particle, we computed the total frequency shift of the photons emitted by the massive particle. By expanding the total frequency shift for large orbits of the emitter, low peculiar redshifts, and small azimuthal angle, we find that the leading contributions to the frequency shift originate from the kinematic shift and peculiar shift, and the quantum effects make a negligible contribution to it since AλA_{\lambda} first appears in the ninth term of the series expansion and the quantum parameter is largely suppressed by the masses of astrophysical black holes.

The total frequency shift ztotz_{\text{tot}} includes three special cases, which respectively correspond to the max redshift ztot1z_{\text{tot}_{1}} for ϕ=π2\phi=\frac{\pi}{2}, the max blueshift ztot2z_{\text{tot}_{2}} for ϕ=π2\phi=-\frac{\pi}{2}, and the gravitational redshift ztot3z_{\text{tot}_{3}} for ϕ=0\phi=0. By solving an inverse problem, we obtained two closed analytic expressions for the black hole mass and quantum parameter in terms of the redshift ztot1z_{\text{tot}_{1}}, blueshift ztot2z_{\text{tot}_{2}} and gravitational redshift ztot3z_{\text{tot}_{3}}, as well as the orbital radius rer_{\text{e}} of the massive particle and the peculiar redshift z0z_{0} which codifies the peculiar motion of the black hole. Namely, the mass and the quantum parameter describing this particular polymerized black hole can be completely determined by few direct observables. Besides, we found that all frequency shifts of photons decrease with the quantum parameter AλA_{\lambda} since the quantum effects weaken the strength of the gravitational field. The frequency shifts also decrease with rer_{\text{e}} since the larger orbital radius means a weaker gravitational field. In addition, we investigated the frequency shifts versus the peculiar redshift and found that the redshift and gravitational redshift increase with the receding peculiar velocity of the black hole while the blueshift decreases with it. There exists a zero point for the blueshift where the kinematic effect is counteracted by the peculiar redshift at a particular receding velocity and similar zero point also respectively appears for the redshift and gravitational redshift at a particular approaching velocity.

Then, we investigated the refractive plasma effects on the frequency shifts of photons. After introducing an index of refraction of the plasma medium, the photons can still travel along the geodesic of an effective metric. Then, we calculated the frequency shifts of the photons following the same way as in the case of the absence of the plasma, and found that the black hole mass and quantum parameter also can be analytically solved in terms of the redshift z^tot1\hat{z}_{\text{tot}_{1}}, blueshift z^tot2\hat{z}_{\text{tot}_{2}} and gravitational redshift z^tot3\hat{z}_{\text{tot}_{3}}. Besides, we analyzed the dispersion effects of the plasma and found that the frequency shifts always decrease with the plasma density.

Finally, we stress that this general relativistic formalism has been applied to a lot of Schwarzschild black holes hosted at the core of different active galactic nuclei, such as, NGC 4258 [52], TXS-2226-184 [53], J1346+5228 [54], and so on, where the authors estimate the mass-to-distance ratio of these black holes with a Bayesian fitting method by using the large data set of the frequency shifts of photons emitted by water masers orbiting them. In our forthcoming work, we will estimate the mass of the polymerized black hole hosted at different active galactic nuclei and evaluate the quantum effects with the general relativistic formalism.

Acknowledgements.
This work was supported by the National Natural Science Foundation of China (Grants Nos. 11975072, 11875102 and 11835009), the National SKA Program of China (Grants Nos. 2022SKA0110200 and 2022SKA0110203), the China Postdoctoral Science Foundation (Grant No. 2021M700729), and Shaanxi Provincial Education Department (Grant No. 21JK0556).

References

  • [1] A. Ashtekar, T. Pawlowski, and P. Singh, Phys. Rev. D 74, 084003 (2006).
  • [2] A. Ashtekar and P. Singh, Classical Quantum Gravity 28, 213001 (2011).
  • [3] A. Corichi, T. Vukasinac, and J. A. Zapata, Phys. Rev. D 76, 044016 (2007).
  • [4] A. Ashtekar and M. Bojowald, Classical Quantum Gravity 23, 391 (2006).
  • [5] L. Modesto, Classical Quantum Gravity 23, 5587 (2006).
  • [6] M. Campiglia, R. Gambini, and J. Pullin, Classical Quantum Gravity 24, 3649 (2007).
  • [7] C. G. Böhmer and K. Vandersloot, Phys. Rev. D 76, 104030 (2007).
  • [8] D. W. Chiou, Phys. Rev. D 78, 064040 (2008).
  • [9] A. Corichi and P. Singh, Classical Quantum Gravity 33, 055006 (2016).
  • [10] J. Olmedo, S. Saini, and P. Singh, Classical Quantum Gravity 34, 225011 (2017).
  • [11] N. Bodendorfer, F. M. Mele, and J. Münch, Classical Quantum Gravity 36, 187001 (2019).
  • [12] N. Bodendorfer, F. M. Mele, and J. Münch, Classical Quantum Gravity 36, 195015 (2019).
  • [13] J. G. Kelly, R. Santacruz, and E. Wilson-Ewing, Phys. Rev. D 102, 106024 (2020).
  • [14] V. Faraoni and A. Giusti, Symmetry 12, 1264 (2020).
  • [15] W. C. Gan, N. O. Santos, F. W. Shu, and A. Wang, Phys. Rev. D 102, 124030 (2020).
  • [16] M. Bojowald, Phys. Rev. D 102, 046006 (2020).
  • [17] F. Sartini and M. Geiller, Phys. Rev. D 103, 066014 (2021).
  • [18] N. Bodendorfer, F. M. Mele, and J. Münch, Phys. Lett. B 819, 136390 (2021).
  • [19] N. Bodendorfer, F. M. Mele, and J. Münch, Classical Quantum Gravity 38, 095002 (2021).
  • [20] M. Bouhmadi-López, S. Brahma, C. Y. Chen, P. Chen, and D. Yeom, J. Cosmol. Astropart. Phys. 07 (2020) 066.
  • [21] F. M. Mele, J. Münch, and S. Pateloudis, J. Cosmol. Astropart. Phys. 02 (2022) 011.
  • [22] A. Arbey, J. Auffinger, M. Geiller, E. R. Livine, and F. Sartini, Phys. Rev. D 103, 104010 (2021).
  • [23] A. Arbey, J. Auffinger, M. Geiller, E. R. Livine, and F. Sartini, Phys. Rev. D 104, 084016 (2021).
  • [24] Q. M. Fu and X. Zhang, Phys. Rev. D 105, 064020 (2022).
  • [25] Y. L. Liu, Z. Q. Feng, and X. D. Zhang, Phys. Rev. D 105, 084068 (2022).
  • [26] S. Brahma, C. Y. Chen, and D. H. Yeom, Phys. Rev. Lett. 126, 181301 (2021).
  • [27] J. B. Achour, S. Brahma, and J. Uzan, J. Cosmol. Astropart. Phys. 03 (2020) 041.
  • [28] J. B. Achour and J. P. Uzan, Phys. Rev. D 102, 124041 (2020).
  • [29] J. B. Achour, S. Brahma, S. Mukohyama, and J. P. Uzan, J. Cosmol. Astropart. Phys. 09 (2020) 020.
  • [30] K. Blanchette, S. Das, S. Hergott, and S. Rastgoo, Phys. Rev. D 103, 084038 (2021).
  • [31] J. Münch, Phys. Rev. D 104, 046019 (2021).
  • [32] Y. C. Liu, J. X. Feng, F. W. Shu, and A. Wang, Phys. Rev. D 104, 106001 (2021).
  • [33] R. K. Walia, arXiv: 2207.02106.
  • [34] S. J. Yang, Y. P. Zhang, S. W. Wei, and Y. X. Liu, J. High Energy Phys. 04 (2022) 066.
  • [35] A. Herrera-Aguilar and U. Nucamendi, Phys. Rev. D 92, 045024 (2015).
  • [36] M. Sharif and S. Iftikhar, Eur. Phys. J. C 76, 404 (2016).
  • [37] G. V. Kraniotis, Eur. Phys. J. C 81, 147 (2021).
  • [38] D. Ujjal, Chin. J. Phys. 70, 213 (2021).
  • [39] R. Becerril, S. Valdez-Alvarado, and U. Nucamendi, Phys. Rev. D 94, 124024 (2016).
  • [40] L. A. López and J. C. Olvera, Eur. Phys. J. Plus 136, 64 (2021).
  • [41] R. Becerril, S. Valdez-Alvarado, U. Nucamendi, P. Sheoran, and J. M. Dávila, Phys. Rev. D 103, 084054 (2021).
  • [42] U. Rashmi, N. Hemwati, and K. D. Purohit, Classical Quantum Gravity 35, 025003 (2017).
  • [43] P. Sheoran, A. Herrera-Aguilar, and U. Nucamendi, Phys. Rev. D 97, 124049 (2018).
  • [44] S. Komarov, A. Gorbatsievich, and A. Tarasenko, Gen. Relativ. Gravit. 50, 132 (2018).
  • [45] L. A. López and N. Bretón, Astrophys. Space Sci. 366, 55 (2021).
  • [46] A. Herrera-Aguilar, R. Lizardo-Castro, and U. Nucamendi, Astron. Nachr. 342, 198 (2021).
  • [47] P. Banerjee, A. Herrera-Aguilar, M. Momennia, and U. Nucamendi, Phys. Rev. D 105, 124037 (2022).
  • [48] A. Rogers, Mon. Not. R. Astron. Soc. 451, 17 (2015).
  • [49] V. Perlick, O. Yu. Tsupko, and G. S. Bisnovatyi-Kogan, Phys. Rev. D 92, 104031 (2015).
  • [50] F. Atamurotov, B. Ahmedov, and A. Abdujabbarov, Phys. Rev. D 92, 084005 (2015).
  • [51] R. R. Caldwell and D. Grin, Phys. Rev. Lett. 100, 031301 (2008).
  • [52] U. Nucamendi, A. Herrera-Aguilar, R. Lizardo-Castro, and O. L. Cruz, Astrophys. J. Lett. 917, L14 (2021).
  • [53] A. Villalobos-Ramirez, O. Gallardo-Rivera, A. Herrera-Aguilar, and U. Nucamendi, Astron. Astrophys. 662, L9 (2022).
  • [54] A. Villalobos-Ramírez, A. González-Juárez, M. Momennia, and A. Herrera-Aguilar, arXiv: 2211.06486.