Black hole surrounded by the pseudo-isothermal dark matter halo

. The abundance of dark matter in the actual universe motivates us to construct the black hole spacetime enveloped by dark matter. In this paper, we derive a new spherically symmetric black hole surrounded by the pseudo-isothermal dark matter halo, and then explore the effects of the pseudo-isothermal halo profile on a rotating black hole at the M87 galactic center, aiming to achieve a black hole solution that aligns with those found in the real universe. Using the Newman-Janis method, we derive a rotating black hole solution encompassed by the pseudo-isothermal halo, which is consistent with observations of actual black holes that are believed to possess spin. Our investigation focuses on the impact of the pseudo-isothermal halo on the black hole event horizon, time-like and null orbits, as well as the black hole shadow. We find that as the spin parameter a increases, the interval between the inner event horizon and the outer event horizon of the rotating black hole surrounded by the pseudo-isothermal halo in M87 diminishes. This leads to the formation of an extreme black hole. The presence of dark matter, however, has minimal effect on the event horizon. Moreover, in the M87 as the spin parameter a increases, the black hole shadow deviates increasingly from a standard circle, with larger spin parameters causing more pronounced distortion relative to the standard circle. Surprisingly, we observe that the dark matter density has very little influence on the shadow of the black hole surrounded by the pseudo-isothermal halo in the M87. This study contributes to a deeper understanding of black hole structures and the role of dark matter in the universe.


Introduction
In 1915, Einstein put forth the groundbreaking theory of general relativity (GR) in the realm of physics, where he linked the gravitational field to the curvature of spacetime.Within this theory emerged the concept of black holes (BH), extraordinary celestial objects in the universe.Their existence was not observed directly but rather predicted based on the principles of general relativity.In 1916, Schwarzschild [1] became the first to solve the field equations of general relativity, resulting in a static spherically symmetric vacuum solution that described a non-rotating black hole with mass alone.Later, in 1963, Kerr derived a rotating, steadystate axisymmetric vacuum black hole solution [2].However, within classical gravitational theory, a singularity exists inside the black hole [3][4][5][6].
On the other hand, the true nature of dark matter is one of the most profound mysteries in astrophysics and cosmology.The widely accepted Λ CDM model has been remarkably successful in explaining the dynamics of the large-scale Universe, suggesting that about 26.8% of our Universe is composed of dark matter, making up a significant 85% of its total mass [7].In the early universe, the dark matter halo has been suggested as a potential factor in understanding this phenomenon [8,9].Although numerical simulations have been used to analyze the growth of black holes within dark matter halos over time, obtaining an analytical form for black holes surrounded by dark matter halo remains challenging.This difficulty arises from the unknown and uncertain behavior of dark matter particles when they interact with black holes.
The presence of a supermassive black hole at the center of a galaxy can significantly enhance the dark matter density, leading to a phenomenon called the "Spike" [10][11][12].However, the Navarro-Frenk-White (NFW) density profile exhibits a "cusp" problem [13], which contradicts observations showing a rather flat density profile.Different dark matter models, such as scalar field dark matter, modified Newtonian dynamics dark matter, and warm dark matter, do not produce the "cusp" in small scales.It remains uncertain whether the "Spike" and "Cusp" manifest in the galactic center.
These challenges serve as motivation for our study, where we examine black holes (spherically symmetric and rotating) in dark matter halos under stationary conditions.By exploring the results obtained, we can investigate various dynamic processes occurring near the black hole, and the energy density of dark matter under relativistic conditions.This research seeks to shed light on the behavior of dark matter in the vicinity of black holes and provides insights into potential solutions for the observed discrepancies between theoretical predictions and observational data.One notable method for extending Schwarzschild-like black holes to Kerr-like black holes is the Newman-Janis (NJ) algorithm [14][15][16], which has proven to be highly suitable for this purpose.Additionally, reference [17] also contributes to the topic.
In 2019, the Event Horizon Telescope (EHT) released the first image of a black hole, capturing the shadow of the supermassive black hole M87 at the center of the Virgo elliptical galaxy [18].Moreover, in 2022, EHT released an image of Sagittarius A*, the supermassive black hole at the center of the Milky Way [19].The information revealed in these black hole photos has provided valuable insights into the shadow, jet, and accretion processes of black holes.As a result of these groundbreaking discoveries, the phenomenon of black hole shadows has attracted significant attention in the scientific literature.Numerous studies have been conducted to delve into this intriguing aspect .In Ref. [120] Tang and Xu presented a spacetime metric for a rotating short-hairy black hole, and studied the influence of short-hairy on the shadow of black holes.These studies have contributed significantly to our understanding of black hole properties and their impact on the surrounding environment.In our study, we will utilize the black hole shadow to gain valuable insights into the effects of pseudo-isothermal halo (PIH) profile [121] on a rotating black hole geometry.To achieve this, we will primarily rely on the formalism developed by Xu et al. [122,123] to obtain the combined dark matter and black hole geometries.By employing these techniques, we aim to elucidate the impact of pseudo-isothermal halo profile on a rotating black hole, specifically in the context of their observable shadows.This investigation can contribute to a deeper understanding of the interaction between black holes and dark matter and may offer valuable clues for unraveling the enigmatic nature of dark matter through astrophysical observations.
In our research, we aim to explore the shadows of rotating black holes surrounded by dark matter and analyze their unique characteristics.By doing so, we hope to provide valuable insights and directions for the experimental detection of this particular spacetime configuration.This investigation will contribute to our understanding of black hole properties and enhance our ability to interpret the observations of black hole shadows in astrophysical contexts.This work is organized as follows.In Sec. 2, we give a pure dark matter spacetime metric by considering the pseudo-isothermal halo profile, and solve the Einstein field equations to obtain the spherically symmetric BH metric surrounded by the pseudo-isothermal halo.In Sec. 3 we derive the rotating black hole surrounded by the pseudo-isothermal halo profile derived by using the Newman-Janis method.In Sec. 4, we obtain the geodesic equation by solving the Hamilton-Jacobi equation in the rotating black hole.In Sec. 5, we present the shadow images of a rotating black hole and analyze the influence of spin and dark matter parameters on the shadow.Sec.6 provides our main conclusion of this work.

Spherically symmetric BH metric surrounded by the PIH
In this work, we study the pseudo-isothermal halo profile [121].The dark matter density profile reads where ρ 0 denotes the central halo density and r 0 denotes the halo core radius.In order to obtain a spherically symmetric black hole surrounded by dark matter, we first need to calculate the mass distribution of the dark matter halo.According to the above dark matter halo profile, we can obtain the mass profile [124]) In the spherically symmetric spacetime, we can find the tangential velocity of the test particle moving in the dark matter halo determined by the mass distribution of the dark matter halo.Therefore, we obtain the tangential velocity (in units G = 1, c = 1) 3) The spherically symmetric spacetime line elements describing pure dark matter halo can be written as (using the method introduced in Ref. [124] and later used by Xu et al. in Ref. [122]) where the A(r) represents the redshift functions, and B(r) represents the shape functions.The tangential velocity is closely related to the redshift function A(r) of the spacetime metric, which reads It should be noted that in our work, we consider the case of A(r) = B(r).Substituting equation ( 2.3) into the above equation, we can obtain the analytical expression of the redshift function and the shape function in the pure dark matter metric Now, we will solve the Einstein field equations to obtain the spherically symmetric BH metric surrounded by the pseudo-isothermal halo.As proposed by Xu et al. [122], the Einstein field equation for the pure pseudo-isothermal dark matter spacetime needs to be solved, which can be written as where g µν denotes the metric tensor of the pure pseudo-isothermal dark matter spacetime, R µν and R denotes the Ricci tensor and Ricci scalar, respectively.The energymomentum tensor of the pure pseudo-isothermal dark matter spacetime can be written as . Therefore, we can obtain . (2.8) If consider a black hole surrounded by the PIH, the energy-momentum tensor becomes T ν µ = T ν µ (BH) + T ν µ (PIH).According to GR, a Schwarzschild black hole is a vacuum solution, which satisfies the condition T ν µ (BH) = 0.If we consider a Schwarzschild-like BH surrounded by the PIH, we can only consider the energy-momentum tensor of the pure PIH.Therefore, we assume the black hole metric surrounded by the PIH has the following form where A P I (r) and B P I (r) denote the metric coefficients of the pure PIH and the unknown functions X 1 (r) and X 2 (r) are determined by the black hole parameters and PIH parameters.As a consequence, the Einstein field equations can be cast in the form Substituting the new black hole metric (2.9) into the above Einstein field equation, we can obtain Using the Schwarzschild black hole as the boundary condition, we can obtain the analytical solution of the above two differential equations dr − 1 r dr − A P I (r). (2.12) For the pseudo-isothermal halo, if A P I (r) = B P I (r), one can find that X 1 (r) = X 2 (r) = −2M/r.The black hole surrounded by the pseudo-isothermal halo can be hence written as where R(r) = r 2 , and 3 Rotating black hole surrounded by the pseudo-isothermal halo The black hole photos taken by EHT and the detection of gravitational wave events indicate that the black holes in the actual universe are Kerr-like black holes, i.e., they usually have spins.Therefore, we will introduce spin into the spherically symmetric black hole spacetime.To achieve this, we will use the Newman-Janis algorithm, which allows us to introduce the spin into the black hole surrounded by the pseudo-isothermal halo.In the NJ method, in order to obtain the spacetime metric of a rotating black hole surrounded by the pseudo-isothermal halo, it is necessary to transform the spherically symmetric metric from Boyer-Lindquist to Eddington-Finkelstein coordinates.The relationship between the Boyer-Lindquist coordinates (t, r, θ, ϕ) and the advanced null coordinate (u, r, θ, ϕ) is The contravariant metric tensor g µν for a spherically symmetric metric can be given by the four basis vectors (l µ , n µ , m µ , mµ ) where the four basis vectors satisfy the orthogonality conditions and the modulus of the basis vector is one, i.e.
For the metric of a rotating black hole surrounded by the pseudo-isothermal halo, the basis vectors are (3.4) According to the NJ method, the spacetime coordinates for different observers satisfy the following complex transformation Using this transformation, the metric coefficient will include spin, i.e., the metric coefficient becomes a function of (r, θ, a).According to the above complex transformation, the metric functions become: and R(r) → Σ(r, θ, a).After considering this change, the null tetrad become Therefore, according to these null tetrads, the contravariant metric tensor of the rotating black hole surrounded by the pseudo-isothermal halo can be written as Moreover, the covariant non-zero metric tensor components are given by Then, we can obtain the rotating black hole surrounded by the pseudo-isothermal halo in Eddington-Finkelstein coordinates (3.9)One significant challenge in describing black holes is the presence of coordinate singularities, particularly at the event horizon.The Boyer-Lindquist coordinates are chosen to avoid singularities at the event horizon, making them regular there, which is crucial for understanding the physical properties of the black hole, especially when studying the behavior of objects and light near the horizon.Moreover, The Boyer-Lindquist coordinates are manifestly covariant under coordinate transformations, which means that the mathematical expressions for physical quantities maintain the same form regardless of the chosen coordinates.The number of non-zero off-diagonal components of the metric is minimal in Boyer-Lindquist coordinates [125].Therefore, transforming the Eddington-Finkelstein coordinates back to Boyer-Lindquist coordinates is necessary.To this purpose, we use the following coordinate transformation du = dt + T 1 (r)dr, dϕ = dϕ ′ + T 2 (r)dr, (3.10) where the T 1 (r) and T 2 (r) take the form [126] T with and Therefore, the spacetime metric of the rotating black hole surrounded by the pseudo-isothermal halo can be given by (3.14) where Σ and ∆ are defined by For a = 0, the rotating black hole surrounded by the pseudo-isothermal halo degenerates into the non-rotating black hole surrounded by the pseudo-isothermal halo, while for ρ 0 = 0, the rotating black hole surrounded by the pseudo-isothermal halo degenerates into the Kerr black hole.The Schwarzschild black hole metric is recovered for a = 0, ρ 0 = 0.The Kerr black hole in general relativity usually has two event horizons, the outer event horizon and the inner event horizon.For the rotating black hole surrounded by the pseudo-isothermal halo we obtained in the previous section, it is worth studying whether its event horizon is affected by dark matter.For a Kerr-like black hole, the event horizon is obtained by solving the equation ∆(r) = 0.For the rotating black hole surrounded by the pseudo-isothermal halo, the event horizon is obtained by solving the following equation For this equation, we cannot give an analytical solution, so we use numerical analysis to study the event horizon structure.In Fig. 1 and Fig. 2, we give ∆(r) as the function of r.From Fig. 1 and Fig. 2, we can see that the rotating black hole surrounded by the pseudo-isothermal halo still has a inner event horizon and event horizon.From Fig. 1, one can find that as the spin parameter a increases, the interval between the inner event horizon and the outer event horizon becomes smaller and smaller, and the rotating dark matter black hole finally becomes an extreme black hole.In Fig. 2, the green curve represents ρ 0 = 0M ⊙ /kpc 3 , the black curve represents ρ 0 = 10M ⊙ /kpc 3 , the red curve represents ρ 0 = 1000M ⊙ /kpc 3 , and the blue curve represents ρ 0 = 100000M ⊙ /kpc 3 .From Fig. 2, one can see that the existence of dark matter has a tiny effect on the event horizon.Here we study the effect of different dark matter densities on the event horizon.The green curve represents ρ 0 = 0M ⊙ /kpc 3 , the black curve represents ρ 0 = 10M ⊙ /kpc 3 , the red curve represents ρ 0 = 1000M ⊙ /kpc 3 , and the blue curve represents ρ 0 = 100000M ⊙ /kpc 3 .
into the black hole and cannot be detected; ii) photons with impact parameter greater than a certain critical value, and they can be detected by the observer.In this section, we will study the geodesic equation of the rotating spacetime surrounded by the pseudo-isothermal dark matter halo.Our main concern is the Hamilton-Jacobi equation, which is defined as where I denotes the Jacobi action, and λ denotes the geodesic affine parameter.Furthermore, the Hamiltonian H reads For the rotational black hole surrounded by the pseudo-isothermal halo, the Hamilton-Jacobi equation is given by By using the method of the separation of variables, the action I assumes the form where µ is the mass of the particle.E and L are the conserved quantities in the motion of photons along the geodesic.Moreover, S r (r) and S θ (θ) denote the radial and angular functions, respectively.Substituting the action I into Eq.( 4.3), and using the functions {S r (r), S θ (θ)} and the conserved quantities {E, L}, we obtain where R(r) and Θ(θ) are defined as and Q represents the Carter constant [127], which is an integral constant introduced by Carter.The equations (4.5) and (4.6) are complete time-like geodesics and null geodesics equations.When µ = 1 they describe the motion of time-like particles, and when µ = 0 they describe the motion of null particles.For the black hole shadow, we consider null geodesics, i.e. µ = 0. Note that R(r) and Θ(θ) must be greater than 0, i.e. Θ(θ) where the impact parameters are defined as The condition satisfied by the stable motion orbit region of the photon is R ′′ (r) > 0, and the condition satisfied by the unstable motion orbit region of the photon is R ′′ (r) < 0. The critical orbit between the stable orbit and the unstable orbit satisfies the condition where r ph denotes the radius of the photo-sphere.According to the above conditions, the critical impact parameters can be written as 5 The shadow shape of rotating BH surrounded by the pseudo-isothermal halo Through the research in the previous section we obtained the geodesics of photons, and based on these results we can calculate the motion of photons detected in a rotating black hole surrounded by the pseudo-isothermal halo for any observer.To draw shadow images, people usually use celestial coordinates (x and y), which is a two-dimensional coordinate system.
For an observer at the position (r, θ), the celestial coordinates are defined as [128,129] x = lim r→∞ −r 2 sin θ dϕ dr , y = lim r→∞ r 2 dθ dr . (5.1) where r represents the space distance from the observer to the rotating black hole surrounded by the pseudo-isothermal halo, while θ represents the angle between the straight line determined by the observer and the center of the black hole, and the rotation axis of the black hole.The motion of photons is described by dϕ/dt and dθ/dt.For asymptotically flat black hole metrics, the celestial coordinates can be simplified as and y = ± η + a 2 cos 2 θ − ξ 2 cot 2 θ.
(5.3)In our work, we consider that the observer is located on the equatorial plane of the rotating black hole surrounded by the pseudo-isothermal halo, i.e., θ = π/2.Therefore, the celestial coordinates can be written as x = −ξ, and y = ± √ η.Using x and y we can present the shadow image of the rotating black hole surrounded by the pseudo-isothermal halo.It allows us to study the influence of the presence of dark matter on the black hole through the image of the shadow.In Fig. 3 we investigate the effect of the spin parameter a on the shadow of the rotating black hole surrounded by the pseudo-isothermal halo.The left panel of Fig. 3 is the shadow of the Kerr black hole when the dark matter density ρ 0 is zero.The right panel of Fig. 3 is the shadow of rotating black hole surrounded by the pseudo-isothermal dark matter halo, where the parameters are set to M = 1, ρ 0 = 6.9 × 10 6 M ⊙ /kpc 3 , r 0 = 91.2kpc.Moreover, the purple solid line is the shadow of the Schwarzschild black hole, which is a standard circle.From Fig. 3, one can find that as the spin parameter a increases, the shadow of the rotating black hole surrounded by the pseudo-isothermal halo gradually deviates from a standard circle.The left half of the shadow image is gradually distorted, while the right half of the shadow is still a standard semicircle.Furthermore, the shadow of the rotating black hole surrounded by the pseudo-isothermal halo is very similar to that of the Kerr black hole.
In Fig. 4 we study the effect of dark matter density on the black hole shadow.From left to right and from top to bottom in Fig. 4 are different spin parameters a = 0.1, a = 0.6, a = 0.9, and a = 0.999.The blue solid line is the shadow of the Schwarzschild black hole, and ρ 0 = 0 corresponds to the shadow of the kerr black hole.The other parameters are set to M = 1, r 0 = 91.2kpc.From Fig. 4, one can find that the dark matter density has a very little influence on the black hole shadow.

Conclusion
The abundance of dark matter in our universe is substantial, leading us to consider a black hole enveloped by dark matter.Our research involves a link of the spherically symmetric Schwarzschild black hole with the pseudo-isothermal halo profile, resulting in a spherically symmetric black hole metric surrounded by dark matter.Observations of black hole photos taken by EHT and gravitational wave events have shown that black holes likely possess spin.To account for this, we employ the Newman-Janis method to derive a rotating black hole surrounded by the pseudo-isothermal halo, which brings us closer to a black hole solution akin to those in the real universe.
Our investigation delves into the effects of the pseudo-isothermal halo profile on the black hole event horizon, time-like and null orbits, and the shadow of the rotating black hole.We find that as the spin parameter a increases, the interval between the inner event horizon and the outer event horizon of the rotating black hole surrounded by the pseudoisothermal halo in M87 diminishes, yielding eventually the formation of an extreme black hole.Moreover, our research reveals that the presence of dark matter has minimal impact on the event horizon.However, as the spin parameter a increases, the shadow of the rotating black hole surrounded by the pseudo-isothermal halo in M87 deviates increasingly from a standard circle, with larger spin parameters causing more pronounced distortion relative to the standard circle.
Interestingly, the dark matter density appears to have very little influence on the black hole shadow.The shadows we observe suggest that dark matter density does not significantly affect the black hole shadow, even with varying densities.The observed shadow images closely resemble those of the Kerr black hole, even in the presence of dark matter, indicating that dark matter has minimal influence on the shadow.This raises the question of whether gravitational waves emitted by the rotating black hole surrounded by the pseudo-isothermal halo would also be challenging to distinguish from those of the Kerr black hole.A relatively important stage in the gravitational wave signal is the quasinormal mode, and it carries the unique characteristics of the background spacetime.Therefore, it is necessary to further study the quasinormal mode of these black holes surrounded by the pseudo-isothermal halo.We intend to explore the quasinormal mode of the rotating black hole surrounded by the pseudo-isothermal halo in M87, which will be presented in a future work.The insights gained from gravitational waves and shadow observations hold the potential to guide experimental detection efforts.

Figure 1 .
Figure 1.∆(r) as the function of r (see Eq. (3.15b)), with ρ 0 = 6.9×10 6 M ⊙ /kpc 3 , r 0 = 91.2kpc.The three curves in the left panel correspond to the results of a = 0.1, a = 0.3 and a = 0.5, respectively.The three curves in the right panel correspond to the results of a = 0.7, a = 0.9 and a = 1, respectively.

Figure 3 .
Figure 3.The shadow of the black hole as the spin parameter a varies.The left panel is the shadow of the Kerr black hole when the dark matter density ρ 0 is zero.The right panel is the shadow of a rotating black hole surrounded by the pseudo-isothermal dark matter halo, where the parameters are set to M = 1, ρ 0 = 6.9 × 10 6 M ⊙ /kpc 3 , r 0 = 91.2kpc.Moreover, the purple solid line is the shadow of the Schwarzschild black hole.

Figure 4 .
Figure 4.The shadow of the black hole as the different dark matter densities ρ 0 varies.From left to right and from top to bottom are different spin parameters a = 0.1, a = 0.6, a = 0.9, and a = 0.999.The blue solid line is the shadow of the Schwarzschild black hole, and ρ 0 = 0 corresponds to the shadow of the Kerr black hole.The other parameters are set to M = 1, r 0 = 91.2kpc.