Processing math: 99%

Quasi-normal modes, emission rate, and shadow of charged AdS black holes with perfect fluid dark matter

  • In this study, we comprehensively investigated charged AdS black holes surrounded by a distinct form of dark matter. In particular, we focused on key elements including the Hawking temperature, quasi-normal modes (QNMs), emission rate, and shadow. We first calculated the Hawking temperature, thereby identifying critical values such as the critical radius and maximum temperature of the black hole, essential for determining its phase transition. Further analysis focused on the QNMs of charged AdS black holes immersed in perfect fluid dark matter (PFDM) within the massless scalar field paradigm. Employing the Wentzel-Kramers-Brillouin (WKB) method, we accurately derived the frequencies of these QNMs. Additionally, we conducted a meticulous assessment of how the intensity of the PFDM parameter α influences the partial absorption cross sections of the black hole, along with a detailed study of the frequency variation of the energy emission rate. The pivotal role of geodesics in understanding astrophysical black hole characteristics is highlighted. Specifically, we examined the influence of the dark matter parameter on photon evolution by computing the shadow radius of the black hole. Our findings distinctly demonstrate the significant impact of the PFDM parameter α on the boundaries of this shadow, providing crucial insights into its features and interactions. We also provide profound insights into the intricate dynamics between a charged AdS black hole, novel dark matter, and various physical phenomena, elucidating their interplay and contributing valuable knowledge to the understanding of these cosmic entities.
      PCAS:
  • The General Theory of Relativity, introduced in 1915, is the foremost comprehensive explanation of gravitational phenomena. Notably, Einstein's field equations, formulated within this theory, predicted the existence of one of the most enigmatic entities in the cosmos: the black hole. Originating from the collapse and subsequent extinguishment of massive stars, black holes manifest as exceedingly compact regions with gravitational forces so intense that even light cannot escape their grasp. The oscillations of black holes have been observed to emit gravitational radiation, carrying pivotal information about their internal characteristics [1, 2]. Observations of supermassive black holes, such as those in M87* [35] and Sagittarius A* [6, 7], provide compelling evidence supporting the existence of these enigmatic entities within our universe. The formation of the shadow of a black hole primarily arises from the gravitational deflection of light [8, 9]. Accurate measurement and analysis of this shadow offer valuable insights into the gravitational field surrounding the black hole [1014].

    The study of black hole oscillations when subjected to external perturbations is central in the analysis of the stability of the spacetime of black holes. This study primarily examined the evolution of the background field or black hole mergers. When an external field is present, it triggers the generation of QNMs in scalar [15], electromagnetic [16], and gravitational fields [17] surrounding the black hole. These QNMs, characterized by complex frequencies, dominate the gravitational waves emitted during perturbation processes. The real component of the frequency represents the oscillation frequency exhibited by the black hole, while the imaginary part denotes its decay rate. The study of these QNMs offers valuable insights into properties within the event horizon radius of a black hole [18]. Gravitational waves are emitted during the merger phase when black holes collide. The first detection of these waves from black hole mergers was announced by LIGO/VIRGO in 2016 [19]. The phenomenon known as ringdown in gravitational wave signals is intriguing, providing essential insights into the complexities of black hole physics. Ring decay refers to the rapid attenuation of oscillations observed towards the end of the waveform, representing QNMs exhibited by remnants of black holes during this phase. These QNMs can be conceptualized as characteristic vibrations or resonances arising when perturbations disrupt the event horizon of a black hole. QNMs in black holes have attracted significant attention from the academia. In Ref. [20], the authors explored the pseudo-spectra of non-horizon exotic compact objects (ECOs) with reflective surfaces situated near Schwarzschild radii. Their findings revealed that the QNMs of ECOs are influenced by their spectral instability, offering valuable insights into the dynamic behavior of peculiar entities surrounding black holes and the mechanism behind gravitational wave signals. Moreover, in Ref. [21], researchers discovered that, despite the distinct QNM spectra exhibited by black holes, very compact objects with rings demonstrate a similar decay phase in their rings. Therefore, precise observation of the late-time decay signal becomes crucial to eventually distinguish differences in QNM spectra. The WKB method is widely acknowledged as one of the primary mathematical tools for studying QNMs. Konoplya and Zhidenko proposed a more precise sixth-order WKB method to accurately compute the QNMs and grey factors of various black holes [2229].

    The primary objective of this study was to explore the Hawking temperature, QNMs, shadow formation, and emission rate of charged AdS black holes encompassed by perfect fluid dark matter. In the context of the standard cosmological model, astronomical observations have led to conclude that dark energy constitutes approximately 73% of the cosmos, while dark matter accounts for approximately 23%, with baryonic matter making up a mere 4% [30]. Theoretical models, such as the quintessence and quintom models, have been proposed to explain these empirical findings [3133]. Interestingly, the quintessence model focuses on elucidating the relationship between the pressure of dark energy and its energy density. Understanding the interplay between quintessence and black holes holds significant importance in the realm of cosmology. Building upon Kiselev's groundbreaking work [34], this study further investigated the Schwarzschild black hole solution enveloped by quintessence. Extensive research has been carried out on various black holes immersed in quintessence matter [3539]. Although real dark matter or dark energy could not satisfy the perfect fluid condition, perfect fluid is a good approximation to describe their behaviour. PFDM has attracted considerable research attention as a novel dark matter model owing to its efficacy in explaining the asymptotically flat rotation curve observed in spiral galaxies [4042]. Despite manifesting gravitational effects across diverse systems, the true nature of PFDM remains enigmatic [4345]. Another important difference between black holes in AdS space and flat or de Sitter space-time is that the horizon of black holes in flat or de Sitter space must have a spherical structure, while the horizon topology of black holes in AdS spacetime could be a zero or exhibits a negative constant curvature surface, except for the case with a positive constant curvature surface. These types of black holes with zero or negative constant curvature horizon have been studied in the literature [4659].

    Recently, Abbas and Ali conducted a detailed analysis on the thermodynamic properties concerning the stability of phase transitions in charged AdS black holes immersed in PFDM [60]. In this study, we aimed to provide pivotal insights into the horizon structure of charged AdS black holes hosting PFDM through a comprehensive analysis of their stability. Our investigation explored the Hawking temperature and QNMs of these black holes and evaluated the energy emission rate akin to Hawking radiation using grey-body factors. Additionally, we aimed at computing the shadow radius of the black hole to offer valuable insights for astronomical observations focused on charged AdS black holes featuring PFDM.

    The structure of this paper is as follows. Section II provides a concise review of the solution representing a charged AdS black hole immersed in PFDM. This section also includes an analysis of the critical radius and maximum temperature associated with the phase transition of the black hole. In Sec. III, our focus shifts to analyzing the real and imaginary components within the spectrum of QNM frequencies pertaining to the black hole. Additionally, we examine the grey-body factor of the black hole and the energy emission rate akin to the Hawking radiation of the black hole. In Sec. IV, we elaborate on the shadow radius of the black hole. This section aims to provide insights into the observational aspects related to the shadow of the black hole. Lastly, Sec. V concludes the paper by summarizing our findings and the implications drawn from the outcomes obtained.

    The electromagnetic and PFDM fields are widely acknowledged to exhibit minimal coupling with gravity and the cosmological constant [61]. The solution for a charged AdS black hole with PFDM results in an expression for the metric as follows [60]:

    N(r)=12mr+Q2r213Λr2+αrlnr|α|.

    (1)

    Based on this metric, we can explore the distinct properties and traits of black holes. As the strength parameter α tends toward zero, the black hole undergoes a transition into an RN AdS black hole. This specific type of black hole retains both the charge Q and mass m, existing within an antisymmetric space. Furthermore, in a scenario where the charge Q is absent, the black hole transforms into a Schwarzschild AdS black hole.

    The subsequent step entails evaluating the Hawking temperature of the black hole, a crucial parameter that determines its inherent thermodynamic properties:

    TH=14πdN(r)dr|r=rh.

    (2)

    We calculate the outer horizon by setting the lapse function to be zero:

    N(rh)=12Mrh+Q2r2h13Λr2h+αrhlnrh|α|=0.

    (3)

    According to Eq. (3), the mass can be expressed as a function of the horizon radius as

    M=12[r3hˉl2+Q2rh+αln(rh|α|)+rh],

    (4)

    where we use Λ=3/ˉl2. Therefore, the Hawking temperature is expressed by the following equation:

    TH=ˉl2(rh(α+rhQ2))+3r4h4πˉl2r3h.

    (5)

    Figure 1 shows the Hawking temperature versus the horizon radius. The analysis demonstrates that the evaporation of the black hole gives rise to a residual mass that emerges in the final stage, contingent upon the α parameter. By setting Eq. (5) to zero, we can calculate this residual mass, which is listed for various α values alongside their associated radii (rrem) in Table 1. Additionally, the critical radius (rcr) and maximum temperature (Tmax), indicative of the phase transition, have also been determined.

    Figure 1

    Figure 1.  (color online) Hawking temperature for Q=1, ˉl=20 and various values of α.

    Table 1

    Table 1.  Q=1, ˉl=20 and various values of parameter α.
    αrcrrremMremTmax
    0.101.650920.9482171.114950.0344216
    0.151.602710.9250661.140470.0359256
    0.201.556130.9025091.156870.0375213
    0.251.511140.8805431.166340.0392136
    0.301.467710.8591651.170160.0410079
    0.351.425810.8383711.169180.0429096
    0.401.385390.8181551.164010.0449242
    0.451.346430.7985121.15510.0470577
    0.501.308870.7794331.142790.0493158
    DownLoad: CSV
    Show Table

    The perturbation caused by the scalar field interacting with the black hole is mathematically expressed through the covariant Klein-Gordon equation in curved space-time,

    1gμ(ggμννΨωlm)=0,

    (6)

    where we use the metric

    ds2=N(r)dt2+N1(r)dt2+r2dΩ2,

    (7)

    with N(r) being determined by Eq. (1).

    Next, we proceed by assuming that the wave function within this spherically symmetric system can be decomposed into the following form:

    Ψωlm(r,t)=Rωl(r)rYlm(θ,φ)eiωt,

    (8)

    where Ylm(θ,φ) represents the spherical harmonics, and ω denotes the complex QNM frequency. Upon substituting the decomposed wave function Ψωlm(r,t) into Eq. (6), we derive the radial wave function Rωl(r) as follows:

    N(r)ddr(N(r)dRωldr)+(ω2Veff.)Rωl=0.

    (9)

    Here, Veff.(r), which is called the effective potential, has the following form:

    Veff.(r)=N(r)l(l+1)r2+N(r)rdNdr.

    (10)

    Then, by employing a tortoise coordinate defined as dr=dr/N(r) in Eq. (9), we transform the equation into a Schrödinger-like form:

    [d2dr2+ω2Veff.]Rωl(r)=0.

    (11)

    Various approaches have been employed to explore the wave function and scattering properties [6267]. One frequently used method for QNM analysis, initially suggested by Schutz and Will [68], is the WKB method [22, 23]. Notably, higher-order corrections for this method have been proposed by Konoplya [24, 25]. In our investigation, we aimed to assess the influence of the parameter α associated with a charged AdS black hole featuring PFDM on its QNMs. These QNMs can be computed by solving the wave equation presented in Eq. (11) using the sixth-order WKB approximation method. Within this approach, the QNMs are determined based on boundary conditions that mandate outgoing waves extending to infinity and incoming waves reaching the event horizon exclusively. The expression for the QNMs is given by

    i(ω2V0)2V06j=2Ωj=n+12.

    (12)

    Here, V0 denotes the second derivative of the potential evaluated at its maximum point r0, while Ωj represents correction terms dependent on the effective potential and its derivatives at this maximum. The QNMs are denoted by ω, having a form ω=ωRiωI. Figures 25 show the real and imaginary parts of the QNMs, computed using the sixth-order WKB method, for various values of the parameter α. These figures are categorized based on the multipole number l, as shown in Figs. 2, 3, 4, and 5, corresponding to l=0,1,2, and 3, respectively. The results presented in the tables indicate a trend in which both the real frequency and magnitude of the imaginary components of the QNMs increase with higher values of the α parameter. Our findings suggest that higher α values lead to higher propagating frequencies and faster damping of scalar waves.

    Figure 2

    Figure 2.  (color online) Variation of real ωR and imaginary ωR QNMs with α, using M=1, Q=0.2, ˉl=20, l=0, and its related monopole nl.

    Figure 3

    Figure 3.  (color online) Variation of real ωR and imaginary ωR QNMs with α, using M=1, Q=0.2, ˉl=20, l=1, and its related monopole nl.

    Figure 4

    Figure 4.  (color online) Variation of real ωR and imaginary ωR QNMs with α, using M=1, Q=0.2, ˉl=20, l=2, and its related monopole nl.

    Figure 5

    Figure 5.  (color online) Variation of real ωR and imaginary ωR QNMs with α, using M=1, Q=0.2, ˉl=20, l=3, and its related monopole nl.

    Furthermore, the WKB method has demonstrated its efficacy in analyzing wave scattering and absorption cross-sections within black hole spacetimes. The grey-body factor, a critical element in determining the tunneling probability through the potential of a black hole [6971], relies on the probability of either incoming waves being absorbed or outgoing waves reaching infinity. This factor can be determined via the WKB method by imposing appropriate boundary conditions expected near the horizon and at infinity [27]:

    Rωl={eiωr+Reiωrif r(rrh),Teiωrif r+(r),

    (13)

    The reflection coefficient is symbolized as R, while the transmission coefficient is defined as T. When employing the WKB method to investigate the influence of the PFDM parameter α on the scattering problem [24, 27, 72], it can be expressed more precisely as follows:

    |R|2=11+e2iπΥ,

    (14)

    |T|2=11+e+2iπΥ=1|R|2,

    (15)

    where Υ is the phase factor; Υ can be expressed as [73]

    Υ=i(ω2V0)2V06j=2Ωj(Υ).

    (16)

    The coefficients Ωj(Υ) represent functions of the effective potential, its derivatives, and the quantity Υ. With this information, we can calculate the grey-body factor using the following expression:

    |Al|2=1|Rl|2=11+e+2iπΥ,

    (17)

    Figure 6 presents the grey-body factors of the scalar field computed using the WKB method for different values of the dark matter parameter α (l=1). The graphical representation demonstrates that an increase in the α value results in a reduction of the grey-body factors. This decrease means that a smaller portion of the scalar field can traverse through the potential barrier. The computation of the partial absorption cross-section involves utilizing the transmission coefficient [7375], which is defined as follows:

    Figure 6

    Figure 6.  (color online) Grey-body factors of the scalar field obtained by the WKB method for M=1, Q=0.2, and ˉl=20. In the upper panel, l=1 and α=0.02,0.04,0.06,0.08; the lower panel shows plots for α=0.02 and l=1,2,3.

    σl=π(2l+1)ω2|Tl(ω)|2,

    (18)

    and the total absorption cross-section is derived by aggregating the partial absorption cross-sections, where l denotes the number of modes and ω denotes the frequency. Hence, it can be expressed as follows:

    σabs=lσl.

    (19)

    Figure 7 depicts the partial absorption cross-sections as a function of the frequency for l=1, showing variations for different values of the dark matter parameter.

    Figure 7

    Figure 7.  (color online) Partial absorption cross-section of the scalar field obtained from the WKB method for M=1 and Q=0.2, and ˉl=20. The upper panel presents plots for l=1 and various values of α, whereas the lower panel shows plots for α=0.02 and different l values.

    The plots show that higher values of the dark matter parameter correspond to lower partial absorption cross-sections. Furthermore, we analyzed the behavior of the energy emission rate with respect to the frequency.

    Near black hole horizons, certain particle pairs can be generated by quantum fluctuations. Pairs with positive energy can escape from the black hole via tunneling, leading to the phenomenon known as Hawking radiation. This process is responsible for the eventual evaporation of black holes within a specific timeframe. It is important to mention that, in many scenarios, the final mass of the black hole is not considered zero owing to the existence of a residual mass in the final phase. This residual mass holds significant importance because it ensures that not all information within the black hole is lost.

    According to Ref. [76], the energy emission rates are calculated as

    d2Edωdt=2π2σω3eωT1,

    (20)

    where ω symbolizes the frequency of a photon, and T denotes the Hawking temperature corresponding to the outer event horizon.

    Figure 8 shows the emission rate versus the frequency ω for different values of α. Note that the maximum emission rate can be determined by setting the derivative of the emission rate with respect to the frequency to zero. Additionally, according to the emission rate equation, when frequencies tend towards zero or infinity, the emission rate diminishes. Figure 8 clearly shows that, as the parameter α increases, the peak of the emission rate, representing the maximum energy emission, decreases. This escalation in α also leads to a shift of the peak of the emission rate towards lower frequencies.

    Figure 8

    Figure 8.  (color online) Emission rate for M=1 , Q=1, and ˉl=20 for various values of α.

    Considering the initial image of a black hole revealed by the Event Horizon Telescope [3], exploring shadows and geodesics becomes profoundly significant for understanding the physical attributes of black holes [7791]. This study aimed to investigate the equation of null geodesics and analyze the influence of the dark matter parameter on photon evolution. Assuming a spherically symmetric spacetime metric as that presented in Eq. (7), the Hamilton-Jacobi action is employed following the methodology proposed in [92]:

    Sτ=12gμνSτμSτν.

    (21)

    In this context, the Jacobi action S and the arbitrary affine parameter τ are utilized. We assume the variables to be separable, allowing us to obtain

    S=12m2τEt+Lϕ+Sr(r)+Sθ(θ).

    (22)

    We express Sr(r) and Sθ(θ) as functions solely dependent on r and θ, respectively. Additionally, we consider a photon with zero mass (m=0) and introduce the constants of motion—energy denoted as E and angular momentum as L—along the trajectory of the photon [92]. We can derive the null geodesic equations using Eqs. (21) and (22) as follows:

    dtdτ=EN(r),drdτ=R(r)r2,dθdτ=±Q(θ)r2,dφdτ=Lcsc2θr2,

    (23)

    where R(r) and Q(θ) are defined as

    R(r)=E2r4(K+L2r2N(r))

    (24)

    Q(θ)=KL2cot2θ,

    (25)

    and K is the Carter constant [93]. In this study, we only consider the equatorial plane (θ=π/2). Now, we concentrate on the radial equation:

    (drdτ)2+Veff.(r)=0.

    (26)

    Here, the effective potential is expressed as

    Veff.(r)=(L2+KN(r)r2E2)

    (27)

    and two impact parameters are introduced

    ξ=LE,η=KE2.

    (28)

    The critical radius of a photon can be determined by applying the unstable condition to the effective potential:

    Veff.=dVeff.dr=0,

    (29)

    which leads to

    2rN(r)N(r)\uparrowvertr=rc2r(2Q2r3+2mr2+αr22rΛ3αln(r|α|)r2)12mr+Q2r213Λr2+αrlnr|α|\uparrowvertr=rc=0.

    (30)

    By applying both Eqs. (27) and (28) to the above expression, we obtain [92]

    ξ2+η=r2cN(rc).

    (31)

    The shadow radius for an AdS black hole can be expressed as

    Rsh=rcN(ro)N(rc)=rc12mro+Q2r2o13Λr2o+αrolnro|α|12mrc+Q2r2c13Λr2c+αrclnrc|α|,

    (32)

    where rcrepresents the radius of the photon sphere obtained from Eq. (30), and ro denotes the distance between the observer and black hole. By employing this formula, we can now establish constraints on the values of Λ and α based on the EHT data. Next, we examine the observational restriction on the parameters Λ and α by analyzing the data acquired from SgrA in connection with the black hole shadow. According to the findings reported in Ref. [6] regarding SgrA, the observational data include the distance of SgrA from Earth, which is determined to be D=8277±33pc, and the mass of the black hole, which is calculated to be MSgrA=4.3±0.013×106M [94, 95].

    Figure 9 shows the shadow shape for different values of α. Note that, when α decreases, the radius of the shadow increases. This indicates that the shadow radius is reduced as the effect of the PFDM parameter increases.

    Figure 9

    Figure 9.  (color online) Shadow shape for M=1, Q=1, and ˉl=20 for various values of α.

    Figure 10 shows the shadow shape for different values of ˉl. Note that an increase in ˉl results in a reduction in the radius of the shadow. However, as the value of ˉl increases, the discrepancy in the shadow radius becomes negligible, albeit not entirely absent.

    Figure 10

    Figure 10.  (color online) Shadow shape for M=1, Q=0.1 and α=0.01for various values of ˉl.

    The uncertainty at 1σ and 2σ levels is illustrated in Fig. 11, which depicts the variation of the shadow radius with the PFDM parameter α and cosmological constant Λ for SgrA.

    Figure 11

    Figure 11.  (color online) The upper panel shows the shadow radius of the black hole as a function of the parameter α for three Λ values: 0.007, 0.005, and 0.003, represented by pink, dark red, and gray lines, respectively. In the lower panel, the shadow radius is plotted as a function of Λ, with α values of 0.01, 0.05, and 0.1, represented by pink, dark red, and gray lines, respectively. The black hole mass is set to M=1. The dark cyan regions indicate α and Λ values inconsistent with observations of stellar dynamics for SgrA. The white and light green shaded areas correspond to the EHT horizon-scale images of Sgr A* at 1σ and 2σ confidence levels, respectively. In both panels, these shaded regions represent α and Λ values consistent with the averaged Keck and VLTI mass-to-distance ratio priors for SgrA.

    We list the upper and lower limits of the acceptable parameter values in Tables 2 and 3.

    Table 2

    Table 2.  Acceptable values for α can be determined based on three different values of Λ=0.007,0.005, and 0.003, represented in Fig. 11 by pink, dark red, and gray lines, respectively. These values correspond to the shadow radius of the black hole that matches the EHT horizon-scale image of SgrAwithin 1σ and 2σ confidence levels.
    α1σ2σ
    UpperLowerUpperLower
    0.050.12
    Λ=0.0072.972.263.341.88
    14.6214.9014.4615.02
    0.030.10
    Λ=0.0053.052.333.411.96
    14.3814.6414.2314.75
    0.020.08
    Λ=0.0033.132.413.492.04
    14.1414.3814.0014.48
    DownLoad: CSV
    Show Table

    Table 3

    Table 3.  Acceptable values for Λ can be determined based on three different values of α=0.01,0.05, and 0.1, represented in Fig. 11 by pink, dark red, and gray lines, respectively. These values correspond to the shadow radius of the black hole that matches the EHT horizon-scale image of SgrAwithin 1σ and 2σ confidence levels.
    Λ1σ2σ
    UpperLowerUpperLower
    α=0.010.016040.001470.02585
    α=0.050.024760.007130.036650.00031
    α=0.10.034300.013340.048430.00523
    DownLoad: CSV
    Show Table

    We conclude that, considering various potential ranges of α and Λ at both 1σ and 2σ uncertainty levels, SgrA could be a charged AdS-PFDM black hole within the current precision of astrophysical observations.

    In summary, this study focused on charged AdS black holes with PFDM as the primary research subject. Our exploration involved probing the Hawking temperature of this type of black holes and analyzing the massless scalar field equation. Additionally, we employed the WKB method to determine the frequencies of the QNMs associated with the black hole while examining its partial absorption cross-section. Using the grey-body factor, we aimed at understanding how the energy emission rate varies according to frequency. Lastly, we aimed at determining the shadow radius of this particular type of black hole. Our findings encompass crucial insights into these aspects of charged AdS black holes with PFDM.

    (i) During the final phase of black hole evaporation, a residual mass is observed, with a maximum remnant mass of 1.17016(Q=20,l=1,α=0.3). The Hawking temperature increases with the dark matter parameter α. Additionally, we provide insights into the critical radius and maximum temperature associated with the phase transition of black holes.

    (ii) Utilizing the WKB method, we calculated the QNM frequency of the black hole. The results reveal that an increase in PFDM parameter α enhances both the real and absolute values of the imaginary part of the QNMs. Higher PFDM values lead to higher frequencies for propagating scalar waves and a faster damping rate.

    (iii) We observed a decrease in the partial absorption cross-section as the PFDM parameter α increases. The energy emission rate diminishes as the frequency approaches zero and infinity, causing the maximum energy emission to decrease with higher parameter values. Consequently, the peak of the emission rate shifts towards lower frequencies in this scenario.

    (iv) Figure 11 shows constraints on α and Λ based on EHT data. The upper panel displays two upper bounds and three lower bounds for both 1σ and 2σ levels. Additionally, the upper panel in Fig. 11 shows changes in the behavior of the shadow radius before and after two extremum points. Specifically, before α=0.68 and after α=9.73, the shadow size decreases with increasing values of α; in the range (0.68,9.73), the shadow size increases with increasing values of α. In the lower panel, both upper and lower bounds are presented for the lines corresponding to α=0.05, 0.1, while for the line corresponding to α=0.01 at the 2σ level, only one upper bound exists without a lower bound. However, at the 1σ level, both upper and lower bounds are provided. This finding may be valuable for astronomical observations of AdS black holes filled with PFDM.

    We are grateful to the referees for their insightful comments and suggestions, which have allowed us to improve this paper significantly.

    [1] W. Israel, Phys. Rev 164, 1776 (1967) doi: 10.1103/PhysRev.164.1776
    [2] B. Carter, Phys. Rev. Lett. 26, 331 (1971) doi: 10.1103/PhysRevLett.26.331
    [3] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 875, L1 (2019) doi: 10.3847/2041-8213/ab0ec7
    [4] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 875, L4 (2019) doi: 10.3847/2041-8213/ab0e85
    [5] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 875, L6 (2019) doi: 10.3847/2041-8213/ab1141
    [6] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 930, L12 (2022) doi: 10.3847/2041-8213/ac6674
    [7] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 930, L14 (2022) doi: 10.3847/2041-8213/ac6429
    [8] J. M. Bardeen, W. H. Press, and S. A. Teukolsky, Astrophys. J. 178, 347 (1972) doi: 10.1086/151796
    [9] X. J. Gao, J. M. Chen, H. Zhang et al., Phys. Lett. B 822, 136683 (2021) doi: 10.1016/j.physletb.2021.136683
    [10] Y. Mizuno, Z. Younsi, C. M. Fromm et al., Nature Astron. 2, 585 (2018) doi: 10.1038/s41550-018-0449-5
    [11] V. Perlick and O.Y. Tsupko, Phys. Rep. 947, 39 (2022) doi: 10.1016/j.physrep.2021.10.004
    [12] X. J. Gao, X. K. Yan, Y. Yin et al., Eur. Phys. J. C 83, 281 (2023) doi: 10.1140/epjc/s10052-023-11414-0
    [13] D. Anish, A. Saha, and S. Gangopadhyay, Class. Quant. Grav. 38, 065015 (2021) doi: 10.1088/1361-6382/abd95b
    [14] C. K. Qiao and M. Zhou, JCAP 12, 005 (2023) doi: 10.1088/1475-7516/2023/12/005
    [15] R. Moderski and M. Rogatko, Phys. Rev. D 64, 044024 (2001) doi: 10.1103/PhysRevD.64.044024
    [16] M. Bouhmadi-Lopez, S. Brahma, C. Y. Chen et al., arXiv: 1902.07874 [gr-qc]
    [17] C. Y. Chen and P. Chen, Phys. Rev. D 99, 104003 (2019) doi: 10.1103/PhysRevD.99.104003
    [18] C. V. Vishveshwara, Nature 227, 936 (1970) doi: 10.1038/227936a0
    [19] B. P. Abbott et al. (LIGO Scientific, Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016) doi: 10.1103/PhysRevLett.116.061102
    [20] V. Boyanov, K. Destounis, R. Panosso Macedo et al., Phys. Rev. D 107, 064012 (2023) doi: 10.1103/PhysRevD.107.064012
    [21] V. Cardoso, E. Franzin, and P. Pani, Phys. Rev. Lett. 116, 171101 (2016) doi: 10.1103/PhysRevLett.116.171101
    [22] R. A. Konoplya and A. Zhidenko, Rev. Mod. Phys. 83, 793 (2011) doi: 10.1103/RevModPhys.83.793
    [23] R. A. Konoplya, Phys. Rev. D 66, 044009 (2002) doi: 10.1103/PhysRevD.66.044009
    [24] R. A. Konoplya, A. Zhidenko, and A F Zinhailo, Class. Quant. Grav. 36, 155002 (2019) doi: 10.1088/1361-6382/ab2e25
    [25] R. A. Konoplya, Phys. Rev. D 68, 024018 (2003) doi: 10.1103/PhysRevD.68.024018
    [26] R. A. Konoplya and A. Zinhailo, Phys. Rev. D 99, 104060 (2019) doi: 10.1103/PhysRevD.99.104060
    [27] R. A. Konoplya, Phys. Lett. B 804, 135363 (2020) doi: 10.1016/j.physletb.2020.135363
    [28] R. A. Konoplya and A. Zhidenko, JCAP 08, 008 (2023) doi: 10.1088/1475-7516/2023/08/008
    [29] R. A. Konoplya and A. Zhidenko, Phys. Rev. D 107, 044009 (2023) doi: 10.1103/PhysRevD.107.044009
    [30] D. Kastor et al., Class. Quant. Grav. 27, 235014 (2010) doi: 10.1088/0264-9381/27/23/235014
    [31] Ė. Komatsu et al. (WMAP Collaboration), Astrophys. J. Suppl. 192, 18 (2011) doi: 10.1088/0067-0049/192/2/18
    [32] R. Pantig and A. Övgün , Eur. Phys. J. C 82, 391 (2022) doi: 10.1140/epjc/s10052-022-10319-8
    [33] A. G. Riess, L. G. Strolger, J. Tonry et al., Astrophys. J. 607, 665 (2004) doi: 10.1086/383612
    [34] V. V. Kiselev, Class. Quantum Grav 20, 1187 (2003) doi: 10.1088/0264-9381/20/6/310
    [35] S. B. Chen, B. Wang, and R. Su, Phys. Rev. D 77, 124011 (2008) doi: 10.1103/PhysRevD.77.124011
    [36] G. Q. Li, Phys. Lett. B 735, 256 (2014) doi: 10.1016/j.physletb.2014.06.047
    [37] A. Younas, M. Jamil, and S. Hussain, Phys. Rev. D 92, 084042 (2015) doi: 10.1103/PhysRevD.92.084042
    [38] H. Chen, B. C. Lütfüolu, H. Hassanabadi et al., Phys. Lett. B 827, 136994 (2022) doi: 10.1016/j.physletb.2022.136994
    [39] B. Hamil and B.C. Lütfüoğlu , Nucl. Phys. B 990, 116191 (2023) doi: 10.1016/j.nuclphysb.2023.116191
    [40] G. Siddhartha, T. Matos, D. Nunez et al., arXiv: astro-ph/0003105
    [41] H. X. Zhang, Y. Chen, T. C. Ma et al., Chin. Phys. C 45, 055103 (2021) doi: 10.1088/1674-1137/abe84c
    [42] S. Qi, Q. Li, Y. Zhang et al., Mod. Phys. Lett. A 38, 2350102 (2023) doi: 10.1142/S021773232350102X
    [43] K. Saurabh and K. Jusufi, Eur. Phys. J. C 81, 490 (2021) doi: 10.1140/epjc/s10052-021-09280-9
    [44] S. Shaymatov, B. Ahmedov, and M. Jamil, Eur. Phys. J. C 81, 588 (2021) doi: 10.1140/epjc/s10052-021-09398-w
    [45] A. Das, A. R. Chowdhury, and S. Gangopadhyay, Class. Quant. Grav. 41, 015018 (2024) doi: 10.1088/1361-6382/ad0fb7
    [46] J. P.S. Lemos, Class. Quant. Grav. 12, 1081 (1995) doi: 10.1088/0264-9381/12/4/014
    [47] C.-G. Huang and C.-B. Liang, Phys. Lett. A 201, 27 (1995) doi: 10.1016/0375-9601(95)00229-V
    [48] Jose P.S. Lemos and Vilson T. Zanchin, Phys. Rev. D 54, 3840 (1996) doi: 10.1103/PhysRevD.54.3840
    [49] R.-G. Cai and Y.-Z. Zhang, Phys. Rev. D 54, 4891 (1996) doi: 10.1103/physrevd.54.4891
    [50] R.-G. Cai, J.-Y. Ji, and K.-S. Soh, Phys. Rev. D 57, 6547 (1998) doi: 10.1103/PhysRevD.57.6547
    [51] S. Aminneborg, I. Bengtsson, S. Holst et al., Class. Quant. Grav. 13, 2707 (1996) doi: 10.1088/0264-9381/13/10/010
    [52] R. B. Mann, Class. Quant. Grav., 14, L109 (1997) doi: 10.1088/0264-9381/14/5/007
    [53] L. Vanzo, Phys. Rev. D 56, 6475 (1997) doi: 10.1103/PhysRevD.56.6475
    [54] D. R. Brill, J. Louko, and P. Peldan, Phys. Rev. D 56, 3600 (1997) doi: 10.1103/PhysRevD.56.3600
    [55] D. Klemm, V. Moretti, and L. Vanzo, Phys. Rev. D 57, 6127 (1998) doi: 10.1103/PhysRevD.57.6127
    [56] R.-G. Cai and K.-S. Soh, Phys. Rev. D 59, 044013 (1999) doi: 10.1103/PhysRevD.59.044013
    [57] D. Birmingham, Class. Quant. Grav. 16, 1197 (1999) doi: 10.1088/0264-9381/16/4/009
    [58] R.-G. Cai, Phys. Rev. D 65, 084014 (2002) doi: 10.1103/PhysRevD.65.084014
    [59] M.-Y. Zhang, H. Chen, H. Hassanabadi et al., Eur. Phys. J. C 83, 773 (2023) doi: 10.1140/epjc/s10052-023-11933-w
    [60] G. Abbas and R. H. Ali, Eur. Phys. J. C 83, 407 (2023) doi: 10.1140/epjc/s10052-023-11580-1
    [61] M. H. Li and K. C. Yang, Phys. Rev. D 86, 123015 (2012) doi: 10.1103/PhysRevD.86.123015
    [62] V. Ferrari and B. Mashhoon, Phys. Rev. D 30, 295 (1984) doi: 10.1103/PhysRevD.30.295
    [63] E. W. Leaver, Proc. Roy. Soc. Lond. A 402, 285 (1985) doi: 10.1098/rspa.1985.0119
    [64] N. Heidari and H. Hassanabadi, Phys. Lett. B 839, 137814 (2023) doi: 10.1016/j.physletb.2023.137814
    [65] L. Dong, Y. Yang, A. Övgün et al., Eur. Phys. J. C 83, 565 (2023) doi: 10.1140/epjc/s10052-023-11739-w
    [66] A. Baruah, A. Övgün, and A. Deshamukhya , Ann. Phys. 455, 169393 (2023) doi: 10.1016/j.aop.2023.169393
    [67] E. Berti, V. Cardoso, and A. O. Starinets, Class. Quant. Grav. 26, 163001 (2009) doi: 10.1088/0264-9381/26/16/163001
    [68] B. F. Schutz and C. M. Will, Astrophys. J. 291, L33 (1985) doi: 10.1086/184453
    [69] i. Sakalli and S. Kanzi, Turk. J. Phys. 46, 51 (2022) doi: 10.55730/1300-0101.2691
    [70] L. C. Crispino, A. Higuchi, E. S. Oliveira et al., Phys. Rev. D 87, 104034 (2013) doi: 10.1103/PhysRevD.87.104034
    [71] V. Cardoso, M. Cavaglia, and L. Gualtieri, Phys. Rev. Lett. 96, 071301 (2006) doi: 10.1103/PhysRevLett.96.071301
    [72] J. A. V Campos, M. A. Anacleto, F.A. Brito et al., Sci. Rep. 12, 8516 (2022) doi: 10.1038/s41598-022-12343-w
    [73] S. Iyer and C. M. Will, Phys. Rev. D 35, 3621 (1987) doi: 10.1103/PhysRevD.35.3621
    [74] S. Dey and S. Chakrabarti, Eur. Phys. J. C 79, 504 (2019) doi: 10.1140/epjc/s10052-019-7004-0
    [75] H. Chen, H. Hassanabadi, B. C. Lütfüolu et al., Gen. Rel. Grav. 54, 143 (2022) doi: 10.1007/s10714-022-03037-9
    [76] S. Hawking, Commun. Math. Phys. 43, 199 (1975) doi: 10.1007/BF02345020
    [77] K. Jusufi, Phys. Rev. D 101, 124063 (2020) doi: 10.1103/PhysRevD.101.124063
    [78] H. Falcke, F. Melia, and E. Agol, Astrophys. J. 528, L13 (1999) doi: 10.1086/312423
    [79] F. Atamurotov, A. Abdujabbarov, and B. Ahmedov, Phys. Rev. D 88, 064004 (2013) doi: 10.1103/PhysRevD.88.064004
    [80] L. Amarilla and E. F. Eiroa, Phys. Rev. D 85, 064019 (2012) doi: 10.1103/PhysRevD.85.064019
    [81] P. C. Li, M. Y. Guo, and B. Chen, Phys. Rev. D 101, 084041 (2020) doi: 10.1103/PhysRevD.101.084041
    [82] I. Cimdiker, D. Demir, and A. Övgün , Phys. Dark Univ. 34, 100900 (2021) doi: 10.1016/j.dark.2021.100900
    [83] A. F. Zakharov, F. Paolis, G. Ingrosso et al., New Astron. Rev. 56, 64 (2012) doi: 10.1016/j.newar.2011.09.002
    [84] S. Vagnozzi, R. Roy, Y. D. Tsai et al., Class. Quant. Grav. 40, 16 (2023) doi: 10.1088/1361-6382/acd97b
    [85] R. Kumar and S. G. Ghosh, Astrophys. J. 892, 78 (2020) doi: 10.3847/1538-4357/ab77b0
    [86] B. Cuadros-Melgar, R. Fontana, and J. Oliveira, Phys. Lett. B 811, 135966 (2020) doi: 10.1016/j.physletb.2020.135966
    [87] R. Pantig, P. K. Yu, E. T. Rodulfo et al., Ann. Phys. 436, 168722 (2022) doi: 10.1016/j.aop.2021.168722
    [88] G. Mustafa, F. Atamurotov, I. Hussain et al., Chin. Phys. C 46, 125107 (2022) doi: 10.1088/1674-1137/ac917f
    [89] B. Hamil and B. C. Lütfüoğlu, Phys. Dark Univ. 42, 101293 (2023) doi: 10.1016/j.dark.2023.101293
    [90] Y. Yang, D. Liu, A. Övgün et al., Eur. Phys. J. C 84, 63 (2024) doi: 10.1140/epjc/s10052-024-12412-6
    [91] Y. Yang, D. Liu, A. Övgün et al., Phys. Rev. D 109, 024002 (2024) doi: 10.1103/PhysRevD.109.024002
    [92] K. Jusufi, Phys. Rev. D 101, 084055 (2020) doi: 10.1103/PhysRevD.101.084055
    [93] B. Carter, Phys. Rev. 174, 1559 (1968) doi: 10.1103/PhysRev.174.1559
    [94] E. González, K. Jusufi, G. Leon et al., Phys. Dark Univ. 42, 101304 (2023) doi: 10.1016/j.dark.2023.101304
    [95] S. Capozziello, S. Zare, D. F. Mota et al., J. Cosmol. Astropart. Phys. 05, 027 (2023) doi: 10.1088/1475-7516/2023/05/027
  • [1] W. Israel, Phys. Rev 164, 1776 (1967) doi: 10.1103/PhysRev.164.1776
    [2] B. Carter, Phys. Rev. Lett. 26, 331 (1971) doi: 10.1103/PhysRevLett.26.331
    [3] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 875, L1 (2019) doi: 10.3847/2041-8213/ab0ec7
    [4] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 875, L4 (2019) doi: 10.3847/2041-8213/ab0e85
    [5] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 875, L6 (2019) doi: 10.3847/2041-8213/ab1141
    [6] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 930, L12 (2022) doi: 10.3847/2041-8213/ac6674
    [7] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 930, L14 (2022) doi: 10.3847/2041-8213/ac6429
    [8] J. M. Bardeen, W. H. Press, and S. A. Teukolsky, Astrophys. J. 178, 347 (1972) doi: 10.1086/151796
    [9] X. J. Gao, J. M. Chen, H. Zhang et al., Phys. Lett. B 822, 136683 (2021) doi: 10.1016/j.physletb.2021.136683
    [10] Y. Mizuno, Z. Younsi, C. M. Fromm et al., Nature Astron. 2, 585 (2018) doi: 10.1038/s41550-018-0449-5
    [11] V. Perlick and O.Y. Tsupko, Phys. Rep. 947, 39 (2022) doi: 10.1016/j.physrep.2021.10.004
    [12] X. J. Gao, X. K. Yan, Y. Yin et al., Eur. Phys. J. C 83, 281 (2023) doi: 10.1140/epjc/s10052-023-11414-0
    [13] D. Anish, A. Saha, and S. Gangopadhyay, Class. Quant. Grav. 38, 065015 (2021) doi: 10.1088/1361-6382/abd95b
    [14] C. K. Qiao and M. Zhou, JCAP 12, 005 (2023) doi: 10.1088/1475-7516/2023/12/005
    [15] R. Moderski and M. Rogatko, Phys. Rev. D 64, 044024 (2001) doi: 10.1103/PhysRevD.64.044024
    [16] M. Bouhmadi-Lopez, S. Brahma, C. Y. Chen et al., arXiv: 1902.07874 [gr-qc]
    [17] C. Y. Chen and P. Chen, Phys. Rev. D 99, 104003 (2019) doi: 10.1103/PhysRevD.99.104003
    [18] C. V. Vishveshwara, Nature 227, 936 (1970) doi: 10.1038/227936a0
    [19] B. P. Abbott et al. (LIGO Scientific, Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016) doi: 10.1103/PhysRevLett.116.061102
    [20] V. Boyanov, K. Destounis, R. Panosso Macedo et al., Phys. Rev. D 107, 064012 (2023) doi: 10.1103/PhysRevD.107.064012
    [21] V. Cardoso, E. Franzin, and P. Pani, Phys. Rev. Lett. 116, 171101 (2016) doi: 10.1103/PhysRevLett.116.171101
    [22] R. A. Konoplya and A. Zhidenko, Rev. Mod. Phys. 83, 793 (2011) doi: 10.1103/RevModPhys.83.793
    [23] R. A. Konoplya, Phys. Rev. D 66, 044009 (2002) doi: 10.1103/PhysRevD.66.044009
    [24] R. A. Konoplya, A. Zhidenko, and A F Zinhailo, Class. Quant. Grav. 36, 155002 (2019) doi: 10.1088/1361-6382/ab2e25
    [25] R. A. Konoplya, Phys. Rev. D 68, 024018 (2003) doi: 10.1103/PhysRevD.68.024018
    [26] R. A. Konoplya and A. Zinhailo, Phys. Rev. D 99, 104060 (2019) doi: 10.1103/PhysRevD.99.104060
    [27] R. A. Konoplya, Phys. Lett. B 804, 135363 (2020) doi: 10.1016/j.physletb.2020.135363
    [28] R. A. Konoplya and A. Zhidenko, JCAP 08, 008 (2023) doi: 10.1088/1475-7516/2023/08/008
    [29] R. A. Konoplya and A. Zhidenko, Phys. Rev. D 107, 044009 (2023) doi: 10.1103/PhysRevD.107.044009
    [30] D. Kastor et al., Class. Quant. Grav. 27, 235014 (2010) doi: 10.1088/0264-9381/27/23/235014
    [31] Ė. Komatsu et al. (WMAP Collaboration), Astrophys. J. Suppl. 192, 18 (2011) doi: 10.1088/0067-0049/192/2/18
    [32] R. Pantig and A. Övgün , Eur. Phys. J. C 82, 391 (2022) doi: 10.1140/epjc/s10052-022-10319-8
    [33] A. G. Riess, L. G. Strolger, J. Tonry et al., Astrophys. J. 607, 665 (2004) doi: 10.1086/383612
    [34] V. V. Kiselev, Class. Quantum Grav 20, 1187 (2003) doi: 10.1088/0264-9381/20/6/310
    [35] S. B. Chen, B. Wang, and R. Su, Phys. Rev. D 77, 124011 (2008) doi: 10.1103/PhysRevD.77.124011
    [36] G. Q. Li, Phys. Lett. B 735, 256 (2014) doi: 10.1016/j.physletb.2014.06.047
    [37] A. Younas, M. Jamil, and S. Hussain, Phys. Rev. D 92, 084042 (2015) doi: 10.1103/PhysRevD.92.084042
    [38] H. Chen, B. C. Lütfüolu, H. Hassanabadi et al., Phys. Lett. B 827, 136994 (2022) doi: 10.1016/j.physletb.2022.136994
    [39] B. Hamil and B.C. Lütfüoğlu , Nucl. Phys. B 990, 116191 (2023) doi: 10.1016/j.nuclphysb.2023.116191
    [40] G. Siddhartha, T. Matos, D. Nunez et al., arXiv: astro-ph/0003105
    [41] H. X. Zhang, Y. Chen, T. C. Ma et al., Chin. Phys. C 45, 055103 (2021) doi: 10.1088/1674-1137/abe84c
    [42] S. Qi, Q. Li, Y. Zhang et al., Mod. Phys. Lett. A 38, 2350102 (2023) doi: 10.1142/S021773232350102X
    [43] K. Saurabh and K. Jusufi, Eur. Phys. J. C 81, 490 (2021) doi: 10.1140/epjc/s10052-021-09280-9
    [44] S. Shaymatov, B. Ahmedov, and M. Jamil, Eur. Phys. J. C 81, 588 (2021) doi: 10.1140/epjc/s10052-021-09398-w
    [45] A. Das, A. R. Chowdhury, and S. Gangopadhyay, Class. Quant. Grav. 41, 015018 (2024) doi: 10.1088/1361-6382/ad0fb7
    [46] J. P.S. Lemos, Class. Quant. Grav. 12, 1081 (1995) doi: 10.1088/0264-9381/12/4/014
    [47] C.-G. Huang and C.-B. Liang, Phys. Lett. A 201, 27 (1995) doi: 10.1016/0375-9601(95)00229-V
    [48] Jose P.S. Lemos and Vilson T. Zanchin, Phys. Rev. D 54, 3840 (1996) doi: 10.1103/PhysRevD.54.3840
    [49] R.-G. Cai and Y.-Z. Zhang, Phys. Rev. D 54, 4891 (1996) doi: 10.1103/physrevd.54.4891
    [50] R.-G. Cai, J.-Y. Ji, and K.-S. Soh, Phys. Rev. D 57, 6547 (1998) doi: 10.1103/PhysRevD.57.6547
    [51] S. Aminneborg, I. Bengtsson, S. Holst et al., Class. Quant. Grav. 13, 2707 (1996) doi: 10.1088/0264-9381/13/10/010
    [52] R. B. Mann, Class. Quant. Grav., 14, L109 (1997) doi: 10.1088/0264-9381/14/5/007
    [53] L. Vanzo, Phys. Rev. D 56, 6475 (1997) doi: 10.1103/PhysRevD.56.6475
    [54] D. R. Brill, J. Louko, and P. Peldan, Phys. Rev. D 56, 3600 (1997) doi: 10.1103/PhysRevD.56.3600
    [55] D. Klemm, V. Moretti, and L. Vanzo, Phys. Rev. D 57, 6127 (1998) doi: 10.1103/PhysRevD.57.6127
    [56] R.-G. Cai and K.-S. Soh, Phys. Rev. D 59, 044013 (1999) doi: 10.1103/PhysRevD.59.044013
    [57] D. Birmingham, Class. Quant. Grav. 16, 1197 (1999) doi: 10.1088/0264-9381/16/4/009
    [58] R.-G. Cai, Phys. Rev. D 65, 084014 (2002) doi: 10.1103/PhysRevD.65.084014
    [59] M.-Y. Zhang, H. Chen, H. Hassanabadi et al., Eur. Phys. J. C 83, 773 (2023) doi: 10.1140/epjc/s10052-023-11933-w
    [60] G. Abbas and R. H. Ali, Eur. Phys. J. C 83, 407 (2023) doi: 10.1140/epjc/s10052-023-11580-1
    [61] M. H. Li and K. C. Yang, Phys. Rev. D 86, 123015 (2012) doi: 10.1103/PhysRevD.86.123015
    [62] V. Ferrari and B. Mashhoon, Phys. Rev. D 30, 295 (1984) doi: 10.1103/PhysRevD.30.295
    [63] E. W. Leaver, Proc. Roy. Soc. Lond. A 402, 285 (1985) doi: 10.1098/rspa.1985.0119
    [64] N. Heidari and H. Hassanabadi, Phys. Lett. B 839, 137814 (2023) doi: 10.1016/j.physletb.2023.137814
    [65] L. Dong, Y. Yang, A. Övgün et al., Eur. Phys. J. C 83, 565 (2023) doi: 10.1140/epjc/s10052-023-11739-w
    [66] A. Baruah, A. Övgün, and A. Deshamukhya , Ann. Phys. 455, 169393 (2023) doi: 10.1016/j.aop.2023.169393
    [67] E. Berti, V. Cardoso, and A. O. Starinets, Class. Quant. Grav. 26, 163001 (2009) doi: 10.1088/0264-9381/26/16/163001
    [68] B. F. Schutz and C. M. Will, Astrophys. J. 291, L33 (1985) doi: 10.1086/184453
    [69] i. Sakalli and S. Kanzi, Turk. J. Phys. 46, 51 (2022) doi: 10.55730/1300-0101.2691
    [70] L. C. Crispino, A. Higuchi, E. S. Oliveira et al., Phys. Rev. D 87, 104034 (2013) doi: 10.1103/PhysRevD.87.104034
    [71] V. Cardoso, M. Cavaglia, and L. Gualtieri, Phys. Rev. Lett. 96, 071301 (2006) doi: 10.1103/PhysRevLett.96.071301
    [72] J. A. V Campos, M. A. Anacleto, F.A. Brito et al., Sci. Rep. 12, 8516 (2022) doi: 10.1038/s41598-022-12343-w
    [73] S. Iyer and C. M. Will, Phys. Rev. D 35, 3621 (1987) doi: 10.1103/PhysRevD.35.3621
    [74] S. Dey and S. Chakrabarti, Eur. Phys. J. C 79, 504 (2019) doi: 10.1140/epjc/s10052-019-7004-0
    [75] H. Chen, H. Hassanabadi, B. C. Lütfüolu et al., Gen. Rel. Grav. 54, 143 (2022) doi: 10.1007/s10714-022-03037-9
    [76] S. Hawking, Commun. Math. Phys. 43, 199 (1975) doi: 10.1007/BF02345020
    [77] K. Jusufi, Phys. Rev. D 101, 124063 (2020) doi: 10.1103/PhysRevD.101.124063
    [78] H. Falcke, F. Melia, and E. Agol, Astrophys. J. 528, L13 (1999) doi: 10.1086/312423
    [79] F. Atamurotov, A. Abdujabbarov, and B. Ahmedov, Phys. Rev. D 88, 064004 (2013) doi: 10.1103/PhysRevD.88.064004
    [80] L. Amarilla and E. F. Eiroa, Phys. Rev. D 85, 064019 (2012) doi: 10.1103/PhysRevD.85.064019
    [81] P. C. Li, M. Y. Guo, and B. Chen, Phys. Rev. D 101, 084041 (2020) doi: 10.1103/PhysRevD.101.084041
    [82] I. Cimdiker, D. Demir, and A. Övgün , Phys. Dark Univ. 34, 100900 (2021) doi: 10.1016/j.dark.2021.100900
    [83] A. F. Zakharov, F. Paolis, G. Ingrosso et al., New Astron. Rev. 56, 64 (2012) doi: 10.1016/j.newar.2011.09.002
    [84] S. Vagnozzi, R. Roy, Y. D. Tsai et al., Class. Quant. Grav. 40, 16 (2023) doi: 10.1088/1361-6382/acd97b
    [85] R. Kumar and S. G. Ghosh, Astrophys. J. 892, 78 (2020) doi: 10.3847/1538-4357/ab77b0
    [86] B. Cuadros-Melgar, R. Fontana, and J. Oliveira, Phys. Lett. B 811, 135966 (2020) doi: 10.1016/j.physletb.2020.135966
    [87] R. Pantig, P. K. Yu, E. T. Rodulfo et al., Ann. Phys. 436, 168722 (2022) doi: 10.1016/j.aop.2021.168722
    [88] G. Mustafa, F. Atamurotov, I. Hussain et al., Chin. Phys. C 46, 125107 (2022) doi: 10.1088/1674-1137/ac917f
    [89] B. Hamil and B. C. Lütfüoğlu, Phys. Dark Univ. 42, 101293 (2023) doi: 10.1016/j.dark.2023.101293
    [90] Y. Yang, D. Liu, A. Övgün et al., Eur. Phys. J. C 84, 63 (2024) doi: 10.1140/epjc/s10052-024-12412-6
    [91] Y. Yang, D. Liu, A. Övgün et al., Phys. Rev. D 109, 024002 (2024) doi: 10.1103/PhysRevD.109.024002
    [92] K. Jusufi, Phys. Rev. D 101, 084055 (2020) doi: 10.1103/PhysRevD.101.084055
    [93] B. Carter, Phys. Rev. 174, 1559 (1968) doi: 10.1103/PhysRev.174.1559
    [94] E. González, K. Jusufi, G. Leon et al., Phys. Dark Univ. 42, 101304 (2023) doi: 10.1016/j.dark.2023.101304
    [95] S. Capozziello, S. Zare, D. F. Mota et al., J. Cosmol. Astropart. Phys. 05, 027 (2023) doi: 10.1088/1475-7516/2023/05/027
  • 加载中

Figures(11) / Tables(3)

Get Citation
Hao Chen, Shi-Hai Dong, Sepehr Hassanabadi, Narges Heidari and Hassan Hassanabadi. Quasi-normal modes, Emission rate and Shadow of the charged AdS black hole with perfect fluid dark matter[J]. Chinese Physics C. doi: 10.1088/1674-1137/ad4af4
Hao Chen, Shi-Hai Dong, Sepehr Hassanabadi, Narges Heidari and Hassan Hassanabadi. Quasi-normal modes, Emission rate and Shadow of the charged AdS black hole with perfect fluid dark matter[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ad4af4 shu
Milestone
Received: 2024-03-13
Article Metric

Article Views(1610)
PDF Downloads(57)
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:

Quasi-normal modes, emission rate, and shadow of charged AdS black holes with perfect fluid dark matter

    Corresponding author: Hao Chen, haochen1249@yeah.net (Corresponding author)
  • 1. School of Physics and Electronic Science, Zunyi Normal University, Zunyi 563006, China
  • 2. Research Center for Quantum Physics, Huzhou University, Huzhou 313000, China
  • 3. Centro de Investigación en Computación, Instituto Politécnico Nacional, UPALM, Mexico City 07700, Mexico
  • 4. Department of Physics, University of Tehran, North Karegar Avenue, Tehran 14395-547, Iran
  • 5. Faculty of Physics, Shahrood University of Technology, Shahrood, Iran
  • 6. Department of Physics, University of Hradec Králové, Rokitanského 62, 500 03 Hradec Králové, Czechia

Abstract: In this study, we comprehensively investigated charged AdS black holes surrounded by a distinct form of dark matter. In particular, we focused on key elements including the Hawking temperature, quasi-normal modes (QNMs), emission rate, and shadow. We first calculated the Hawking temperature, thereby identifying critical values such as the critical radius and maximum temperature of the black hole, essential for determining its phase transition. Further analysis focused on the QNMs of charged AdS black holes immersed in perfect fluid dark matter (PFDM) within the massless scalar field paradigm. Employing the Wentzel-Kramers-Brillouin (WKB) method, we accurately derived the frequencies of these QNMs. Additionally, we conducted a meticulous assessment of how the intensity of the PFDM parameter α influences the partial absorption cross sections of the black hole, along with a detailed study of the frequency variation of the energy emission rate. The pivotal role of geodesics in understanding astrophysical black hole characteristics is highlighted. Specifically, we examined the influence of the dark matter parameter on photon evolution by computing the shadow radius of the black hole. Our findings distinctly demonstrate the significant impact of the PFDM parameter α on the boundaries of this shadow, providing crucial insights into its features and interactions. We also provide profound insights into the intricate dynamics between a charged AdS black hole, novel dark matter, and various physical phenomena, elucidating their interplay and contributing valuable knowledge to the understanding of these cosmic entities.

    HTML

    I.   INTRODUCTION
    • The General Theory of Relativity, introduced in 1915, is the foremost comprehensive explanation of gravitational phenomena. Notably, Einstein's field equations, formulated within this theory, predicted the existence of one of the most enigmatic entities in the cosmos: the black hole. Originating from the collapse and subsequent extinguishment of massive stars, black holes manifest as exceedingly compact regions with gravitational forces so intense that even light cannot escape their grasp. The oscillations of black holes have been observed to emit gravitational radiation, carrying pivotal information about their internal characteristics [1, 2]. Observations of supermassive black holes, such as those in M87* [35] and Sagittarius A* [6, 7], provide compelling evidence supporting the existence of these enigmatic entities within our universe. The formation of the shadow of a black hole primarily arises from the gravitational deflection of light [8, 9]. Accurate measurement and analysis of this shadow offer valuable insights into the gravitational field surrounding the black hole [1014].

      The study of black hole oscillations when subjected to external perturbations is central in the analysis of the stability of the spacetime of black holes. This study primarily examined the evolution of the background field or black hole mergers. When an external field is present, it triggers the generation of QNMs in scalar [15], electromagnetic [16], and gravitational fields [17] surrounding the black hole. These QNMs, characterized by complex frequencies, dominate the gravitational waves emitted during perturbation processes. The real component of the frequency represents the oscillation frequency exhibited by the black hole, while the imaginary part denotes its decay rate. The study of these QNMs offers valuable insights into properties within the event horizon radius of a black hole [18]. Gravitational waves are emitted during the merger phase when black holes collide. The first detection of these waves from black hole mergers was announced by LIGO/VIRGO in 2016 [19]. The phenomenon known as ringdown in gravitational wave signals is intriguing, providing essential insights into the complexities of black hole physics. Ring decay refers to the rapid attenuation of oscillations observed towards the end of the waveform, representing QNMs exhibited by remnants of black holes during this phase. These QNMs can be conceptualized as characteristic vibrations or resonances arising when perturbations disrupt the event horizon of a black hole. QNMs in black holes have attracted significant attention from the academia. In Ref. [20], the authors explored the pseudo-spectra of non-horizon exotic compact objects (ECOs) with reflective surfaces situated near Schwarzschild radii. Their findings revealed that the QNMs of ECOs are influenced by their spectral instability, offering valuable insights into the dynamic behavior of peculiar entities surrounding black holes and the mechanism behind gravitational wave signals. Moreover, in Ref. [21], researchers discovered that, despite the distinct QNM spectra exhibited by black holes, very compact objects with rings demonstrate a similar decay phase in their rings. Therefore, precise observation of the late-time decay signal becomes crucial to eventually distinguish differences in QNM spectra. The WKB method is widely acknowledged as one of the primary mathematical tools for studying QNMs. Konoplya and Zhidenko proposed a more precise sixth-order WKB method to accurately compute the QNMs and grey factors of various black holes [2229].

      The primary objective of this study was to explore the Hawking temperature, QNMs, shadow formation, and emission rate of charged AdS black holes encompassed by perfect fluid dark matter. In the context of the standard cosmological model, astronomical observations have led to conclude that dark energy constitutes approximately 73% of the cosmos, while dark matter accounts for approximately 23%, with baryonic matter making up a mere 4% [30]. Theoretical models, such as the quintessence and quintom models, have been proposed to explain these empirical findings [3133]. Interestingly, the quintessence model focuses on elucidating the relationship between the pressure of dark energy and its energy density. Understanding the interplay between quintessence and black holes holds significant importance in the realm of cosmology. Building upon Kiselev's groundbreaking work [34], this study further investigated the Schwarzschild black hole solution enveloped by quintessence. Extensive research has been carried out on various black holes immersed in quintessence matter [3539]. Although real dark matter or dark energy could not satisfy the perfect fluid condition, perfect fluid is a good approximation to describe their behaviour. PFDM has attracted considerable research attention as a novel dark matter model owing to its efficacy in explaining the asymptotically flat rotation curve observed in spiral galaxies [4042]. Despite manifesting gravitational effects across diverse systems, the true nature of PFDM remains enigmatic [4345]. Another important difference between black holes in AdS space and flat or de Sitter space-time is that the horizon of black holes in flat or de Sitter space must have a spherical structure, while the horizon topology of black holes in AdS spacetime could be a zero or exhibits a negative constant curvature surface, except for the case with a positive constant curvature surface. These types of black holes with zero or negative constant curvature horizon have been studied in the literature [4659].

      Recently, Abbas and Ali conducted a detailed analysis on the thermodynamic properties concerning the stability of phase transitions in charged AdS black holes immersed in PFDM [60]. In this study, we aimed to provide pivotal insights into the horizon structure of charged AdS black holes hosting PFDM through a comprehensive analysis of their stability. Our investigation explored the Hawking temperature and QNMs of these black holes and evaluated the energy emission rate akin to Hawking radiation using grey-body factors. Additionally, we aimed at computing the shadow radius of the black hole to offer valuable insights for astronomical observations focused on charged AdS black holes featuring PFDM.

      The structure of this paper is as follows. Section II provides a concise review of the solution representing a charged AdS black hole immersed in PFDM. This section also includes an analysis of the critical radius and maximum temperature associated with the phase transition of the black hole. In Sec. III, our focus shifts to analyzing the real and imaginary components within the spectrum of QNM frequencies pertaining to the black hole. Additionally, we examine the grey-body factor of the black hole and the energy emission rate akin to the Hawking radiation of the black hole. In Sec. IV, we elaborate on the shadow radius of the black hole. This section aims to provide insights into the observational aspects related to the shadow of the black hole. Lastly, Sec. V concludes the paper by summarizing our findings and the implications drawn from the outcomes obtained.

    II.   THERMODYNAMICS OF A CHARGED ADS BLACK HOLE WITH PFDM
    • The electromagnetic and PFDM fields are widely acknowledged to exhibit minimal coupling with gravity and the cosmological constant [61]. The solution for a charged AdS black hole with PFDM results in an expression for the metric as follows [60]:

      N(r)=12mr+Q2r213Λr2+αrlnr|α|.

      (1)

      Based on this metric, we can explore the distinct properties and traits of black holes. As the strength parameter α tends toward zero, the black hole undergoes a transition into an RN AdS black hole. This specific type of black hole retains both the charge Q and mass m, existing within an antisymmetric space. Furthermore, in a scenario where the charge Q is absent, the black hole transforms into a Schwarzschild AdS black hole.

      The subsequent step entails evaluating the Hawking temperature of the black hole, a crucial parameter that determines its inherent thermodynamic properties:

      TH=14πdN(r)dr|r=rh.

      (2)

      We calculate the outer horizon by setting the lapse function to be zero:

      N(rh)=12Mrh+Q2r2h13Λr2h+αrhlnrh|α|=0.

      (3)

      According to Eq. (3), the mass can be expressed as a function of the horizon radius as

      M=12[r3hˉl2+Q2rh+αln(rh|α|)+rh],

      (4)

      where we use Λ=3/ˉl2. Therefore, the Hawking temperature is expressed by the following equation:

      TH=ˉl2(rh(α+rhQ2))+3r4h4πˉl2r3h.

      (5)

      Figure 1 shows the Hawking temperature versus the horizon radius. The analysis demonstrates that the evaporation of the black hole gives rise to a residual mass that emerges in the final stage, contingent upon the α parameter. By setting Eq. (5) to zero, we can calculate this residual mass, which is listed for various α values alongside their associated radii (rrem) in Table 1. Additionally, the critical radius (rcr) and maximum temperature (Tmax), indicative of the phase transition, have also been determined.

      Figure 1.  (color online) Hawking temperature for Q=1, ˉl=20 and various values of α.

      αrcrrremMremTmax
      0.101.650920.9482171.114950.0344216
      0.151.602710.9250661.140470.0359256
      0.201.556130.9025091.156870.0375213
      0.251.511140.8805431.166340.0392136
      0.301.467710.8591651.170160.0410079
      0.351.425810.8383711.169180.0429096
      0.401.385390.8181551.164010.0449242
      0.451.346430.7985121.15510.0470577
      0.501.308870.7794331.142790.0493158

      Table 1.  Q=1, ˉl=20 and various values of parameter α.

    III.   QNMs AND EMISSION RATE OF THE CHARGED AdS BLACK HOLE WITH PFDM
    • The perturbation caused by the scalar field interacting with the black hole is mathematically expressed through the covariant Klein-Gordon equation in curved space-time,

      1gμ(ggμννΨωlm)=0,

      (6)

      where we use the metric

      ds2=N(r)dt2+N1(r)dt2+r2dΩ2,

      (7)

      with N(r) being determined by Eq. (1).

      Next, we proceed by assuming that the wave function within this spherically symmetric system can be decomposed into the following form:

      Ψωlm(r,t)=Rωl(r)rYlm(θ,φ)eiωt,

      (8)

      where Ylm(θ,φ) represents the spherical harmonics, and ω denotes the complex QNM frequency. Upon substituting the decomposed wave function Ψωlm(r,t) into Eq. (6), we derive the radial wave function Rωl(r) as follows:

      N(r)ddr(N(r)dRωldr)+(ω2Veff.)Rωl=0.

      (9)

      Here, Veff.(r), which is called the effective potential, has the following form:

      Veff.(r)=N(r)l(l+1)r2+N(r)rdNdr.

      (10)

      Then, by employing a tortoise coordinate defined as dr=dr/N(r) in Eq. (9), we transform the equation into a Schrödinger-like form:

      [d2dr2+ω2Veff.]Rωl(r)=0.

      (11)

      Various approaches have been employed to explore the wave function and scattering properties [6267]. One frequently used method for QNM analysis, initially suggested by Schutz and Will [68], is the WKB method [22, 23]. Notably, higher-order corrections for this method have been proposed by Konoplya [24, 25]. In our investigation, we aimed to assess the influence of the parameter α associated with a charged AdS black hole featuring PFDM on its QNMs. These QNMs can be computed by solving the wave equation presented in Eq. (11) using the sixth-order WKB approximation method. Within this approach, the QNMs are determined based on boundary conditions that mandate outgoing waves extending to infinity and incoming waves reaching the event horizon exclusively. The expression for the QNMs is given by

      i(ω2V0)2V06j=2Ωj=n+12.

      (12)

      Here, V0 denotes the second derivative of the potential evaluated at its maximum point r0, while Ωj represents correction terms dependent on the effective potential and its derivatives at this maximum. The QNMs are denoted by ω, having a form ω=ωRiωI. Figures 25 show the real and imaginary parts of the QNMs, computed using the sixth-order WKB method, for various values of the parameter α. These figures are categorized based on the multipole number l, as shown in Figs. 2, 3, 4, and 5, corresponding to l=0,1,2, and 3, respectively. The results presented in the tables indicate a trend in which both the real frequency and magnitude of the imaginary components of the QNMs increase with higher values of the α parameter. Our findings suggest that higher α values lead to higher propagating frequencies and faster damping of scalar waves.

      Figure 2.  (color online) Variation of real ωR and imaginary ωR QNMs with α, using M=1, Q=0.2, ˉl=20, l=0, and its related monopole nl.

      Figure 3.  (color online) Variation of real ωR and imaginary ωR QNMs with α, using M=1, Q=0.2, ˉl=20, l=1, and its related monopole nl.

      Figure 4.  (color online) Variation of real ωR and imaginary ωR QNMs with α, using M=1, Q=0.2, ˉl=20, l=2, and its related monopole nl.

      Figure 5.  (color online) Variation of real ωR and imaginary ωR QNMs with α, using M=1, Q=0.2, ˉl=20, l=3, and its related monopole nl.

      Furthermore, the WKB method has demonstrated its efficacy in analyzing wave scattering and absorption cross-sections within black hole spacetimes. The grey-body factor, a critical element in determining the tunneling probability through the potential of a black hole [6971], relies on the probability of either incoming waves being absorbed or outgoing waves reaching infinity. This factor can be determined via the WKB method by imposing appropriate boundary conditions expected near the horizon and at infinity [27]:

      Rωl={eiωr+Reiωrif r(rrh),Teiωrif r+(r),

      (13)

      The reflection coefficient is symbolized as R, while the transmission coefficient is defined as T. When employing the WKB method to investigate the influence of the PFDM parameter α on the scattering problem [24, 27, 72], it can be expressed more precisely as follows:

      |R|2=11+e2iπΥ,

      (14)

      |T|2=11+e+2iπΥ=1|R|2,

      (15)

      where Υ is the phase factor; Υ can be expressed as [73]

      Υ=i(ω2V0)2V06j=2Ωj(Υ).

      (16)

      The coefficients Ωj(Υ) represent functions of the effective potential, its derivatives, and the quantity Υ. With this information, we can calculate the grey-body factor using the following expression:

      |Al|2=1|Rl|2=11+e+2iπΥ,

      (17)

      Figure 6 presents the grey-body factors of the scalar field computed using the WKB method for different values of the dark matter parameter α (l=1). The graphical representation demonstrates that an increase in the α value results in a reduction of the grey-body factors. This decrease means that a smaller portion of the scalar field can traverse through the potential barrier. The computation of the partial absorption cross-section involves utilizing the transmission coefficient [7375], which is defined as follows:

      Figure 6.  (color online) Grey-body factors of the scalar field obtained by the WKB method for M=1, Q=0.2, and ˉl=20. In the upper panel, l=1 and α=0.02,0.04,0.06,0.08; the lower panel shows plots for α=0.02 and l=1,2,3.

      σl=π(2l+1)ω2|Tl(ω)|2,

      (18)

      and the total absorption cross-section is derived by aggregating the partial absorption cross-sections, where l denotes the number of modes and ω denotes the frequency. Hence, it can be expressed as follows:

      σabs=lσl.

      (19)

      Figure 7 depicts the partial absorption cross-sections as a function of the frequency for l=1, showing variations for different values of the dark matter parameter.

      Figure 7.  (color online) Partial absorption cross-section of the scalar field obtained from the WKB method for M=1 and Q=0.2, and ˉl=20. The upper panel presents plots for l=1 and various values of α, whereas the lower panel shows plots for α=0.02 and different l values.

      The plots show that higher values of the dark matter parameter correspond to lower partial absorption cross-sections. Furthermore, we analyzed the behavior of the energy emission rate with respect to the frequency.

      Near black hole horizons, certain particle pairs can be generated by quantum fluctuations. Pairs with positive energy can escape from the black hole via tunneling, leading to the phenomenon known as Hawking radiation. This process is responsible for the eventual evaporation of black holes within a specific timeframe. It is important to mention that, in many scenarios, the final mass of the black hole is not considered zero owing to the existence of a residual mass in the final phase. This residual mass holds significant importance because it ensures that not all information within the black hole is lost.

      According to Ref. [76], the energy emission rates are calculated as

      d2Edωdt=2π2σω3eωT1,

      (20)

      where ω symbolizes the frequency of a photon, and T denotes the Hawking temperature corresponding to the outer event horizon.

      Figure 8 shows the emission rate versus the frequency ω for different values of α. Note that the maximum emission rate can be determined by setting the derivative of the emission rate with respect to the frequency to zero. Additionally, according to the emission rate equation, when frequencies tend towards zero or infinity, the emission rate diminishes. Figure 8 clearly shows that, as the parameter α increases, the peak of the emission rate, representing the maximum energy emission, decreases. This escalation in α also leads to a shift of the peak of the emission rate towards lower frequencies.

      Figure 8.  (color online) Emission rate for M=1 , Q=1, and ˉl=20 for various values of α.

    IV.   SHADOW OF A CHARGED ADS BLACK HOLE WITH PFDM
    • Considering the initial image of a black hole revealed by the Event Horizon Telescope [3], exploring shadows and geodesics becomes profoundly significant for understanding the physical attributes of black holes [7791]. This study aimed to investigate the equation of null geodesics and analyze the influence of the dark matter parameter on photon evolution. Assuming a spherically symmetric spacetime metric as that presented in Eq. (7), the Hamilton-Jacobi action is employed following the methodology proposed in [92]:

      Sτ=12gμνSτμSτν.

      (21)

      In this context, the Jacobi action S and the arbitrary affine parameter τ are utilized. We assume the variables to be separable, allowing us to obtain

      S=12m2τEt+Lϕ+Sr(r)+Sθ(θ).

      (22)

      We express Sr(r) and Sθ(θ) as functions solely dependent on r and θ, respectively. Additionally, we consider a photon with zero mass (m=0) and introduce the constants of motion—energy denoted as E and angular momentum as L—along the trajectory of the photon [92]. We can derive the null geodesic equations using Eqs. (21) and (22) as follows:

      dtdτ=EN(r),drdτ=R(r)r2,dθdτ=±Q(θ)r2,dφdτ=Lcsc2θr2,

      (23)

      where R(r) and Q(θ) are defined as

      R(r)=E2r4(K+L2r2N(r))

      (24)

      Q(θ)=KL2cot2θ,

      (25)

      and K is the Carter constant [93]. In this study, we only consider the equatorial plane (θ=π/2). Now, we concentrate on the radial equation:

      (drdτ)2+Veff.(r)=0.

      (26)

      Here, the effective potential is expressed as

      Veff.(r)=(L2+KN(r)r2E2)

      (27)

      and two impact parameters are introduced

      ξ=LE,η=KE2.

      (28)

      The critical radius of a photon can be determined by applying the unstable condition to the effective potential:

      Veff.=dVeff.dr=0,

      (29)

      which leads to

      2rN(r)N(r)\uparrowvertr=rc2r(2Q2r3+2mr2+αr22rΛ3αln(r|α|)r2)12mr+Q2r213Λr2+αrlnr|α|\uparrowvertr=rc=0.

      (30)

      By applying both Eqs. (27) and (28) to the above expression, we obtain [92]

      ξ2+η=r2cN(rc).

      (31)

      The shadow radius for an AdS black hole can be expressed as

      Rsh=rcN(ro)N(rc)=rc12mro+Q2r2o13Λr2o+αrolnro|α|12mrc+Q2r2c13Λr2c+αrclnrc|α|,

      (32)

      where rcrepresents the radius of the photon sphere obtained from Eq. (30), and ro denotes the distance between the observer and black hole. By employing this formula, we can now establish constraints on the values of Λ and α based on the EHT data. Next, we examine the observational restriction on the parameters Λ and α by analyzing the data acquired from SgrA in connection with the black hole shadow. According to the findings reported in Ref. [6] regarding SgrA, the observational data include the distance of SgrA from Earth, which is determined to be D=8277±33pc, and the mass of the black hole, which is calculated to be MSgrA=4.3±0.013×106M [94, 95].

      Figure 9 shows the shadow shape for different values of α. Note that, when α decreases, the radius of the shadow increases. This indicates that the shadow radius is reduced as the effect of the PFDM parameter increases.

      Figure 9.  (color online) Shadow shape for M=1, Q=1, and ˉl=20 for various values of α.

      Figure 10 shows the shadow shape for different values of ˉl. Note that an increase in ˉl results in a reduction in the radius of the shadow. However, as the value of ˉl increases, the discrepancy in the shadow radius becomes negligible, albeit not entirely absent.

      Figure 10.  (color online) Shadow shape for M=1, Q=0.1 and α=0.01for various values of ˉl.

      The uncertainty at 1σ and 2σ levels is illustrated in Fig. 11, which depicts the variation of the shadow radius with the PFDM parameter α and cosmological constant Λ for SgrA.

      Figure 11.  (color online) The upper panel shows the shadow radius of the black hole as a function of the parameter α for three Λ values: 0.007, 0.005, and 0.003, represented by pink, dark red, and gray lines, respectively. In the lower panel, the shadow radius is plotted as a function of Λ, with α values of 0.01, 0.05, and 0.1, represented by pink, dark red, and gray lines, respectively. The black hole mass is set to M=1. The dark cyan regions indicate α and Λ values inconsistent with observations of stellar dynamics for SgrA. The white and light green shaded areas correspond to the EHT horizon-scale images of Sgr A* at 1σ and 2σ confidence levels, respectively. In both panels, these shaded regions represent α and Λ values consistent with the averaged Keck and VLTI mass-to-distance ratio priors for SgrA.

      We list the upper and lower limits of the acceptable parameter values in Tables 2 and 3.

      α1σ2σ
      UpperLowerUpperLower
      0.050.12
      Λ=0.0072.972.263.341.88
      14.6214.9014.4615.02
      0.030.10
      Λ=0.0053.052.333.411.96
      14.3814.6414.2314.75
      0.020.08
      Λ=0.0033.132.413.492.04
      14.1414.3814.0014.48

      Table 2.  Acceptable values for α can be determined based on three different values of Λ=0.007,0.005, and 0.003, represented in Fig. 11 by pink, dark red, and gray lines, respectively. These values correspond to the shadow radius of the black hole that matches the EHT horizon-scale image of SgrAwithin 1σ and 2σ confidence levels.

      Λ1σ2σ
      UpperLowerUpperLower
      α=0.010.016040.001470.02585
      α=0.050.024760.007130.036650.00031
      α=0.10.034300.013340.048430.00523

      Table 3.  Acceptable values for Λ can be determined based on three different values of α=0.01,0.05, and 0.1, represented in Fig. 11 by pink, dark red, and gray lines, respectively. These values correspond to the shadow radius of the black hole that matches the EHT horizon-scale image of SgrAwithin 1σ and 2σ confidence levels.

      We conclude that, considering various potential ranges of α and Λ at both 1σ and 2σ uncertainty levels, SgrA could be a charged AdS-PFDM black hole within the current precision of astrophysical observations.

    V.   CONCLUSION
    • In summary, this study focused on charged AdS black holes with PFDM as the primary research subject. Our exploration involved probing the Hawking temperature of this type of black holes and analyzing the massless scalar field equation. Additionally, we employed the WKB method to determine the frequencies of the QNMs associated with the black hole while examining its partial absorption cross-section. Using the grey-body factor, we aimed at understanding how the energy emission rate varies according to frequency. Lastly, we aimed at determining the shadow radius of this particular type of black hole. Our findings encompass crucial insights into these aspects of charged AdS black holes with PFDM.

      (i) During the final phase of black hole evaporation, a residual mass is observed, with a maximum remnant mass of 1.17016(Q=20,l=1,α=0.3). The Hawking temperature increases with the dark matter parameter α. Additionally, we provide insights into the critical radius and maximum temperature associated with the phase transition of black holes.

      (ii) Utilizing the WKB method, we calculated the QNM frequency of the black hole. The results reveal that an increase in PFDM parameter α enhances both the real and absolute values of the imaginary part of the QNMs. Higher PFDM values lead to higher frequencies for propagating scalar waves and a faster damping rate.

      (iii) We observed a decrease in the partial absorption cross-section as the PFDM parameter α increases. The energy emission rate diminishes as the frequency approaches zero and infinity, causing the maximum energy emission to decrease with higher parameter values. Consequently, the peak of the emission rate shifts towards lower frequencies in this scenario.

      (iv) Figure 11 shows constraints on α and Λ based on EHT data. The upper panel displays two upper bounds and three lower bounds for both 1σ and 2σ levels. Additionally, the upper panel in Fig. 11 shows changes in the behavior of the shadow radius before and after two extremum points. Specifically, before α=0.68 and after α=9.73, the shadow size decreases with increasing values of α; in the range (0.68,9.73), the shadow size increases with increasing values of α. In the lower panel, both upper and lower bounds are presented for the lines corresponding to α=0.05, 0.1, while for the line corresponding to α=0.01 at the 2σ level, only one upper bound exists without a lower bound. However, at the 1σ level, both upper and lower bounds are provided. This finding may be valuable for astronomical observations of AdS black holes filled with PFDM.

    ACKNOWLEDGMENTS
    • We are grateful to the referees for their insightful comments and suggestions, which have allowed us to improve this paper significantly.

Reference (95)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return