Modeling the black holes surrounded by a dark matter halo in the galactic center of M87

In this paper, the structure of a dark matter halo can be well described by the mass model of M87 and the Einasto profile for the cold dark matter model, i.e., $\rho_{\text{eina}} (r)=\rho_\text{e} \exp ( -2 \alpha ^{-1} ((r/r_\text{e})^\alpha -1 ) )$ (Wang et al. in Nature 585:39-42, 2020). Under these conditions, we construct a solution of a static spherically symmetric black hole in a dark matter halo. Then, using the Newman-janis algorithm, we extend this static solution to the case of rotation, and obtain a solution for the Kerr-like black hole. We prove that this solution of the Kerr-like black hole is indeed a solution to the Einstein field equations. Finally, taking M87 as an example, we study and analyze some physical properties of this Kerr-like black hole, and then compare them with the Kerr black hole. Particularly, from the perspective of the black hole shadow and the fact that the Kerr-like black hole and the Kerr black hole is distinguishable, we give the upper limit of the shape parameter of the Einasto density profile, that is approximately $\alpha<0.22$, which may provide a new method to further improve and perfect the density profile of dark matter model. These research results for the black hole in a dark matter halo may indirectly provide an effective method for detecting the existence of dark matter.


Introduction
Black hole is a fascinating and mysterious celestial body predicted by General Relativity.General Relativity predicts that a sufficiently compact mass can warp spacetime, forming the gravitational field.In other words, when the mass of the celestial body is large enough, the celestial body will make the gravitational field so strong that it can prevent light from escaping, and eventually a black hole may be emerged.The famous physicists K. Schwarzschild and R. Kerr respectively solved the Einstein field equations of the vacuum, and then they obtained the well-known the exact solution of the Schwarzschild black hole [1] and Kerr black hole [2].In 2019, the Event Horizon Telescope (EHT) released the first image of a black hole, that is a supermassive black hole at the center of M87 in the Virgo elliptical galaxy [3].The photo of the black hole is composed of two parts, that is the outer circular bright area and the inner dark area.The bright region near the black hole represents the accretion disk, which appears distorted due to gravitational lensing.And the side of the detector away from the accretion disk will be relatively darker, forming a dark region.The images of the black hole shadow can help us to understand the spacetime geometry of the event horizon for the black hole and its speed of rotation.From these Refs.[4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20], we can learn more about the study of the black hole shadow.In 2022, the EHT Collaboration released the first photo of Sgr A, that is the supermassive black hole in the center of the Milky Way [21].Therefore, it is generally accepted for us that black hole exists in our universe.
On the other hand, more and more observational data indicate the existence of dark matter, such as the rotation curves of spiral galaxies, the cosmic microwave background radiation and gravitational lensing, etc.Based on these observational data, astronomers have proposed a series of dark matter models to study dark matter, such as cold dark matter (CDM) model [22,23], warm dark matter (WDM) model [24], scalar field dark matter (SFDM) model [25] and self-interacting dark matter (SIDM) model [26].In these dark matter models, the density profile of dark matter can always be well described by relatively few variable parameters, such as the density parameter ρ and the characteristic radius r.So far, various density profiles of dark matter models have been proposed, and we list some of them as follows.F. Navarro et al. gave the NFW profile for studying the dark matter halo structure in the CDM model [22].Then, T. Matos et al. gave the spacetime geometry of the pure dark matter halo for this NFW profile [27].L. Urena-Lopez et al. gave the density profile of scalar field dark matter model by solving Einstein Klein-Gordon equation coupled with the scalar field [25].A. Burkert gave the density of the URC profile from the dark matter halos of seven dwarf spiral galaxies [28].The accuracy of the density profile of dark matter may directly determine the difficulty of finding dark matter.Here, we main focus on the density profile proposed by J. Einasto for the CDM model [29].J. Wang et al. showed in their recent research that the NFW profile and the Einasto profile can well describe the structure of dark matter halo in the entire range of the mass [30].More density profiles and models of dark matter can be found in Refs.[31,32].These studies do make positive contributions to the indirect search for dark matter.
Based on these interesting physical backgrounds above, a large number of the individ-uals are beginning to study the interaction between black holes and dark matter.If there is dark matter near the black hole, due to the strong gravity of the black hole itself, the dark matter will be distributed in a certain structure near the black hole and disappear at a certain distance away from the black hole.Using Newtonian approximations, P. Gondolo et al. first obtained the density profile of dark matter near black hole at the center of the Milky Way [33].L. Sadeghian et al. used accurate Schwarzschild geometry to perform relativistic calculations, and corrected the density profile of P. Gondolo's dark matter [34].Xu et al. combined the NFW profile in the dark matter halo structure and the spike structure, and then they obtained the spacetime geometry of the black hole immersed in dark matter [35][36][37].S. Nampalliwar et al. tested the impacts of dark matter spike on the black hole at the center of the Milky Way [38].Ning Dai et al. analyzed the possibility of the existence of a Schwarzschild black hole immersed in a dark matter halo with Hernquist type density distribution from the perspective of gravitational waves [39].Combining the previously mentioned photos of the black holes at the center of the M87 and Sgr A, the authors began to study these black holes in a dark matter halo of M87 and Sgr A and we list some below.K. Jusuf et al. obtained the spacetime geometry of the URC profile dark matter in the center of the M87 galaxy, and studied the impacts of dark matter parameters on the shadow of the black hole [40].H. Davoudiasl et al. used the Event Horizon Telescope to measure the ultralight bosons in M87 and then study the distribution of dark matter on the large-scale structure [41]. A. K. Saha et al. used the Event Horizon Telescope to detect the ultralight axions of Sgr A, and thus gave the mass limit of ultralight bosons [42].V.  [46].At the same time, a recent article published in Nature once again showed that the NFW profile and the Einasto profile of dark matter can well describe the astronomically predicted the structure of dark matter halo [30].For the studies of the Einasto profiles and black holes, we list some below.V. Cardoso et al. studied and explored black holes immersed in the Einasto profile from geodesic equations and axial gravitational perturbations [47].F. Rahaman et al. obtained a solution of a transmissible wormhole in the Einasto profile from Einstein's field equations [48].If it is said that there is dark matter near the black hole, then this dark matter must be able to change the spacetime geometry of the black hole.Therefore, it is meaningful to study the impacts of dark matter halo on the spacetime geometry of the black hole.Therefore, in this work, we will study and explore the impacts of the Einasto profile of dark matter on the black hole in the center of M87.In addition, over the years, the gravitational waves produced by the merger of binary black holes have also been gradually detected by LIGO/Virgo [49][50][51][52].The source of gravitational waves is usually due to perturbations in the black hole or additional sources such as dark matter.Different gravitational waves may be generated according to different black holes, and we can use this to judge the type of black hole or the background space-time.V. Vishveshwara et al. indicated that a black hole in the perturbed state can emit gravitational waves [53].And this gravitational wave is usually dominated by the excited oscillation mode of complex frequency, that is quasinormal mode (QNM) [54].As a characteristic sound, the QNM from the black hole can provide us a new method to identify the black hole in our universe.The recent studies about quasinormal modes of the black hole in the dark matter can refer to these Refs.[55][56][57][58][59].In the future, with the increase of LIGO/Virgo observation data and the photos of the black hole shadow, it is possible to detect such black holes in a dark matter halo, which also provides an effective method for indirect detection of the existence of dark matter.Furthermore, from the geometry of the black hole shadow, we can define observables and establish a connection between them and black hole parameters [18,[60][61][62].In other words, if we observe a shadow image of a black hole in our universe, we can deduce the relevant parameters of the black hole according to the model, even the metric of the black hole.S. Vagnozzi et al. studied the impact of the observables of black hole shadows on some general black holes in M87 and Sgr A, which has received widespread attention [63][64][65][66][67].They show that the observables of black hole shadows can, to a certain extent, be used to test the black hole parameters or the parameters of physical models.These results may be expected to promote the development of physics.On the other hand, for some black holes coupled with other matter, such as dark matter or dark energy [68], these observables of the black hole shadow may help to examine the density profile in a dark matter halo, or the models of dark matter and dark energy.In this work, we are considering black holes coupled with dark matter, and we will try to study these dark matter profiles from the perspective of black hole shadows, which may provide a new method for verifying and improving the density profile of dark matter and even the dark matter model.This is also one of the research objectives of this work.This paper is organized as follows.In section 2, we introduce a much-concerned Einasto profile that can well describe the cold dark matter model, and then we give the metric for this pure dark matter spacetime.In section 3, starting from Einstein field equations, we introduce a spherically symmetric black hole into this pure dark matter spacetime, and then we give the metric of the Schwarzschild-like black hole in a dark matter halo.In section 4, using the Newman-Janis algorithm, we extend the solution of this Schwarzschild-like black hole to the case of the Kerr-like black hole.Besides, in Appendix A, we prove that this solution of the Kerr-like black hole is a solution of Einstein field equations.In section 5, we study and analyze some basic physical properties of this Kerr-like black hole in a dark matter halo, and then we compare these results with the Kerr black hole.Here, we list the following physical properties: 1) The horizons of the Kerr-like black hole; 2) The ergosphere of the Kerr-like black hole; 3) The motion equation of neutral particle near the Kerr-like black hole; 4) The geodesic equations of lightlike particle and photon regions; 5) The shadow of the Kerr-like black hole in a dark matter halo.6) The shadow observables of the Kerr-like black hole in a dark matter halo.Finally, section 6 is our conclusion and discussion.The relevant parameters of the galaxy M87 are ρ e = 6.9 × 10 6 M ⊙ /kpc 3 , r e = 91.2kpcand the mass of the black hole at the center of the M87 is 6.5 × 10 9 M ⊙ [3].In this paper, we mainly use the black hole units of G = c = M SBH = 1, the radius of the Schwarzschild black hole (SBH) is given by r SBH = 2GM SBH /c 2 .

Density profile of dark matter and its spacetime geometry
Recently, one density profile firstly proposed by J. Einasto has received widespread attention [29].More studies on this density profile of dark matter can be found in these Refs.[69][70][71][72][73]. Particularly, J. Wang et al. use this density profile to fit the observed data of the halo masses (over 20 orders of magnitude), and they find that the observed data can be well described by this density profile of the simple two-parameter [30].This density profile is also known as the Einasto profile, and it reads where, r e is the characteristic radius and its logarithmic slope is d lnρ/d lnr = −2, ρ e is the halo density and α is a shape parameter1 .Based on the density profile (2.1), the mass distribution of the dark matter halo is given by ) where, Ei[n, z] is the exponential integral function Ei n (z) and Γ[z] is the gamma function Γ(z).According to Newtonian theory, for a test particle in the equatorial plane of spherical symmetric spacetime, its tangential velocity can be determined by the mass distribution of dark matter [76].Therefore, the tangential velocity V can be defined as (2.3) On the other hand, for the pure dark matter spacetime, its metric can always be written as a static spherical symmetric form where, f (r) is the redshift function and g(r) is the shape function.Here, we mainly consider the simpler case in the spherical symmetric spacetime, that is, f (r) = g(r).The relationship between the redshift function f (r) and the tangential velocity V suggested in Ref. [27]: According to the definitions (2.3) and (2.5), we can get an analytical solution of the metric coefficient f (r) in a dark matter halo, and it finds that where, Γ[a, z 0 , z 1 ] is the generalized incomplete gamma function as Γ[a, z 0 , z 1 ] = Γ(a, z 0 ) − Γ(a, z 1 ), and then the incomplete gamma function Γ[a, z] can be defined by an integral expression as Γ[a, z] = ∞ z t a−1 e −t dt.As an example, we consider the galaxy of M87 and its relevant parameter are ρ e = 6.9 × 10 6 M ⊙ /kpc 3 , r e = 91.2kpc.The mass of the black hole in the center of M87 is M BH = 6.5 × 10 9 M .In this case, the dark matter parameters can be converted in the black hole units (BHU) by the following formula: r e (BHU) = r e /(2GM BH /c 2 ) × r BH and ρ e (BHU) = ρ e /(M BH /(4/3π(2GM BH /c 2 ) 3 )) × ρ BH .For the f (r) in function (2.6), it is not difficult to find that The Eqs. (2.7) show that when the density of dark matter is 0 or the distance is infinity away from the black hole, the dark matter is absent.Therefore, the spacetime return to Minkowski flat spacetime.

Schwarzschild-like black hole in a dark matter halo
In this section, we extend our results (2.6) to the case of black hole in a dark matter halo.We will strictly follow the method recorded in the Refs.[35,36].First of all, we need to solve the Einstein field equation of pure dark matter spacetime, which has the following form where, R µν is Ricci tensor, g µν is the metric of pure dark matter spacetime and R is Ricci scalar.The energy momentum tensor of this spacetime can be defined as T ν µ = g να T µα = diag[−ρ, p r , p, p], and g να is the inverse metric of the spacetime.Taking the metric (2.4) into Eq.(3.1), we can obtain the following equations When the black hole is in a dark matter halo, then the energy momentum tensor can be written as T ν µ = T ν µ (BH) + T ν µ (DM-halo).For Schwarzschild black hole, the energy momentum tensor is 0, that is, T ν µ (BH) = 0.That is to say, we only need to consider the energy-momentum tensor T ν µ (DM-halo) of the dark matter halo.Therefore, based on the above analysis, we can assume the new metric of black hole in a dark matter halo that where, f (r) and g(r) are the coefficients of the pure dark matter halo.Then Eq. (3.1) can be written as Taking the metric (3.3) into Eq.(3.4) and then comparing with Eq. (3.2), we can obtain Eq. (3.5) is the energy-momentum tensors that can be described by two different metrics.
The first equation of Eq. (3.5) is only related to function G1.The second equation is related to functions G1, F 1 and it can be simplified as The general solution to the first equation of (3.5) is G1(r) = C/r.C is an undetermined coefficient, which needs to be determined by boundary condition.Inspired by the Eq.(2.7), here, we use the Schwarzschild black hole as the boundary condition.Under this boundary condition, the new metric G(r) = g(r) + G1(r) of the black hole in a dark matter halo will degenerate into the Schwarzschild black hole at the boundary (G(r) = 1 − 2M/r).At this point, a particular solution that satisfies the boundary conditions is g(r)=1, and then G1(r) = −2M/r.In the end, we determined that the constant C is C = −2M .Therefore, we can obtain the analytical expressions of functions F 1, G1, which are as follows where, M BH is the mass of the black hole.Here, using the case of f (r) = g(r) we assumed before in Eq. (3.5), it is not difficult to find that F 1(r) = G1(r) = −2M BH /r.Therefore, the new metric for the Schwarzschild-like black hole in a dark matter halo is given by where, The terms f (r) and g(r) represent the factor terms for considering pure dark matter halo in the metric (2.6).Finally, the metric coefficients for the Schwarzschild-like black hole in a dark matter halo are as follows If the dark matter is absent, that is, ρ e = 0, and then the metric (3.9) degenerates into the case of the Schwarzschild black hole.

Kerr-like black hole in a dark matter halo
In the previous section, we derived the metric of Schwarzschild-like black hole in a dark matter halo from the Einasto profile.Next, we will extend this Schwarzschild-like black hole metric to the Kerr-like metric.The method we used is Newman-Janis (N-J) algorithm, which is one of the commonly used and efficient methods [36,40,[77][78][79]. Strictly following this method, we firstly transform the metric (3.8) from Boyer-Lindquist (BL) coordinates (t, r, θ, ϕ) to Eddington-Finkelstein (EF) coordinates (u, r, θ, ϕ) and then one can obtain where, F (r), G(r), H(r) = r 2 are the metric functions.In the null tetrad, the inverse metric of Schwarzschild-like black hole can be composed of basis vectors l µ , n µ , m µ , mµ , and it is given by where, these basis vectors l µ , n µ , m µ , mµ are defined as follows where, the basic vectors l and n are real and basic vectors m and m are mutually complex conjugate [80].The relationship between these basis vectors satisfies the conditions of normalization, orthogonality and isotropy.According to the Newman-Janis algorithm, the spacetime coordinates between different observers satisfy the complex transformation.Therefore, we can write the new coordinates as follows where, a is rotation parameter, x ′µ and x µ are the EF coordinate components of the case of rotation and non-rotation, respectively.Under this condition, these new metric functions can be rewritten as F (r) → F(r, a, θ), G(r) → G(r, a, θ), H(r) → Σ(r, a, θ).Besides, in the case of rotation, a new null tetrad Z ′µ can be transformed from the null tetrad of nonrotation through the relation Z ′µ = (∂x ′µ /∂x µ )Z µ [40,81].Therefore, we can obtain this rotation null tetrad with the new metric functions F(r, a, θ), G(r, a, θ) and Σ(r, a, θ) similar to the definitions (4.4), where, m ′µ means the complex conjugate of m ′µ , Σ 1 = r − ia cos θ and Σ 2 = r + ia cos θ.
Taking the results of (4.6) into the equation of new contravariant metric, and then we can obtain the contravariant non-zero components of this new metric which are as follows, where, Σ = Σ 1 Σ 2 = r 2 + a 2 cos 2 θ.Based on these inverse metric and then taking the inverse of (4.8), the new metric of the Kerr-like black hole in EF coordinates becomes (4.9) Finally, we only need to convert this EF coordinates into the coordinates we are familiar with, that is, BL coordinates.Here, we follow the method from the Ref. [78], which introduces more physical parameters and symmetry properties, simplifying the complex procedure of the Newman-Janis algorithm.This method is further studied and applied in these Refs.[15,38,40,82,83].Firstly, we introduce the following transformations, where, K(r) = H(r) G(r)/F (r) and and then taking (4.10), (4.11) into (4.9),we can obtain the final form of metric under BLC coordinates ) where, G, H, K are the function of r we introduced before and Σ = r 2 + a 2 cos 2 θ.With the case F (r) = G(r), H(r) = r 2 we introduced previous section, we find that K(r) = H(r) = r 2 .Finally, the metric (4.12) can be written as Kerr-like form, where, It is easy to find that in (4.13), if the dark matter is absent, (i.e., ρ e = 0, and then g(r) = 1)), the metric of black hole in a dark matter halo degenerates into the Kerr metric (4.15) 5 Some properties of black holes in a dark matter halo Through the introduction of the previous sections, we have strictly given the metrics of the Schwarzschild-like black hole and the Kerr-like black hole in a dark matter halo.The difference between them and Schwarzschild black hole and Kerr black hole is that there is an additional dark matter term, that is f (r), g(r) in (3.8) and (4.13).In this section, we will introduce the physical properties of this special black hole from several aspects such as horizons of black hole, shape of the ergosphere, motion equation of the neutral particle and so on.Finally, we compare these properties with the Kerr black hole, so that readers can understand the black hole in a dark matter halo more intuitively.

Horizons of black hole in a dark matter halo
The event horizon of a black hole is the outermost boundary of the black hole.In black hole physics, the event horizon of a black hole is a type of zero-surface with spacetime symmetry.General relativity shows that the Schwarzschild black hole has only one event horizon, which is determined entirely by the mass of the black hole.And the Kerr black hole has two event horizons, namely the event horizon and the Cauchy horizon, which are completely determined by the mass and angular momentum of the black hole.
For the Kerr-like black hole in a dark matter halo, its horizons can be found by solving ∆(r) = 0 in the formula (4.14), that is In Figure 1, we present the functional image of ∆(r) as a function of the independent variate r both in dark matter spacetime (left panel) and Kerr spacetime (middle panel).The number of these intersections with the x-axis in the image represents the number of roots of the equation ∆(r) = 0, which in black hole physics represents the horizons of the black hole.Our results show that black hole has the same number of horizons as Kerr black hole in a dark matter halo.The outer horizon decreases with the increase of the rotation parameter a, while the inner horizon increases with the increase of a until the inner and outer event horizons are equal, that is, the extreme black hole.In order to distinguish the difference between the black hole in a dark matter halo and the Kerr spacetime, we define their difference δ∆(r) = ∆ Kerr (r) − ∆ DM (r), and then their differences as a function of rotation parameter a are given in Figure 1 (right panel) at a fixed r.Our results show that the ∆(r) of the Kerr spacetime is greater than that of the dark matter halo, that is, ∆ Kerr > ∆ DM .In the right panel of Figure 1, the maximum difference between them is approximately 7 × 10 −11 , and the difference between them will increase with the increasing of variable parameter r.For example, when the rotation parameter a = 0.8, the inner and outer horizons of the Kerr-like black hole in a dark matter halo are r inner = 0.40000000000000002220, r outer = 1.6000000000000000888, respectively, while the inner and outer horizons of the Kerr black hole are r inner = 0.399999999999999996669, r outer = 1.5999999999999998668.This case shows that the inner and outer horizons of the Kerr-like black hole in a dark matter halo are slight larger than those of the Kerr black hole.

Shape of ergosphere of black hole in a dark matter halo
Next, let us study the shape of the ergosphere of this Kerr-like black hole in a dark matter halo.Please note that the ergosphere is only the outer region of the Kerr-like black hole.Firstly, in Figure 2, we give the cross-section diagram of Kerr-like black hole in the dark matter halo in the xz-plane.In the black hole physics, the ergosphere is composed of two parts, one is the area surrounded by the outer infinite redshift surface and the outer event horizon, and the other is the area surrounded by the inner infinite redshift surface and the inner event horizon.The central area in the middle is a singular ring.Then, the inner and outer horizons and the infinite redshift surface of the Kerr-like black hole can be obtained by solving the Eqs.∆(r) = 0 and g tt = 0, respectively, that is Eq.(5.1) and ) From these figures in Figure 2, for the Kerr-like black hole in a dark matter halo, we find that the outer horizon and outer infinite redshift surface decrease with the increase of rotation parameter a, while the inner horizon and inner infinite redshift surface increase with the increase of rotation parameter a.In addition, according to the calculation of the horizons in Figure 1, dark matter may also increase the range of the ergosphere for the Kerr-like black hole.

Motion equations of neutral particle near black hole in a dark matter halo
In this subsection, we are interested in the neutral particle motion near the Kerr-like black hole in a dark matter halo, which can be solved by geodesic equations.Firstly, we rewrite the metric (4.13) in the Kerr-like form where, The geodesic equation of the neutral particle in a dark matter halo satisfy the Lagrangian equation [10], The shape of ergosphere of the Kerr-like black hole in a dark matter halo with different rotation parameter a in the xz-plane.The main calculation parameters are M = 1, α = 0.16, ρ e = 6.9 × 10 6 M /kpc 3 , r e = 91.2kpc.We have converted these main calculation parameters to the black hole units before plotting.
where, g µν is the metric component of (5.3) and ẋµ denotes the first-order partial derivative to the affine parameter of the coordinate x µ .The 4-momentum p µ of the neutral particle from (5.3) is p µ = ∂L /∂ ẋµ = g µν ẋν , and then where, E and L are the energy and angular momentum of the neutral particle, respectively.In addition, p t and p ϕ are a conserved quantity because the Lagrangian is not related to the coordinates t and ϕ.Therefore, the Kerr-like black hole in a dark matter halo has the properties of steady state and axisymmetric.Through the first two equations in (5.6), we can reverse solve ṫ and φ, that is (5.7) Based on the momentum and Lagrangian above, the Hamiltonian of the neutral particle motion can be written as Taking (5.3), (5.5) and (5.6) into account, (5.8) can be written as where, m is the mass of the neutral particle, and m = 0 corresponds to the motion equation of photons.The results (5.7), (5.8), (5.9) are very important, because they can be used to study the properties of particle near black hole.Eq. (5.9) is a second order partial differential equation.Usually, a separation constant can be introduced to reduce it to a radial equation and an angular equation.Where, the radial equation describes the motion of the neutral particle near the black hole.As an example, we study the neutral particle whose orbits lie in the equatorial plane (i.e., θ = π/2, m = 0), and its motion equation is given by (5.10) In particular, when the dark matter is absent (ρ e = 0), that is g(r) = 1 in R(r), the radial equation will return to the case of the Kerr vacuum.

The geodesic equations of lightlike particle and photon regions
In the previous subsection, we calculated the motion equation for the neutral particle using the momentum and Lagrangian based on the metric (5.3).Here, we rewrite this metric as where, (5.12) In this subsection, we will focus on the motion of light particles and their photon regions.Therefore, we need to know the geodesic equations describing the lightlike particles.The lightlike geodesic equation can be obtained by combining the momentum and the Hamilton-Jacobi equation where, S is Hamiltonian-Jacobi action, σ is an affine parameter, x µ is the coordinate component and g µν is the inverse metric of (5.11).For a test particle moving along the geodesic of the Kerr-like black hole, the Hamiltonian-Jacobi action S can always be written in the separated form, where, m is the mass of the particle (here the photon takes m = 0).Taking (5.11) and (5.14) into account, the Eq.(5.13) can be divided into the following two equations, that is and where, S ′ r (r) denotes ∂S r (r)/∂r, S ′ θ (θ) denotes ∂S θ (θ)/∂θ and K is a constant obtained by separating variables, which is called Carter constant [84].Besides, we note that ∂S/∂x µ = p µ , and then it is not difficult to obtain the geodetic equations of the lightlike particles where, the dot (•) in the Eq. ( 5.17) denotes the derivative to the affine parameter σ, that is d/dσ.
Next, let's study and analyze the motion of photons in circular orbits.The motion of a photon usually satisfies the radial equation, that is the first equation of Eq. (5.17), and it can be rewritten as where, the effective potential is where, the function ∆ can be obtained in (5.12).For a circular photon orbit, it satisfies the following condition where, r c is the orbital radius of photon sphere and then we can obtain where, ∆ ′ (r) represents the derivative of ∆ with respect to r. Taking (5.22) into the third equation of (5.18) can be obtained In particular, if V ′′ eff (r c ) > 0, the radial equation corresponding to the lightlike geodesics is unstable, but V ′′ eff (r c ) < 0 is stable.The unstable orbit of the photon will determine the shape of shadow of the black hole.The range of the photon regions r c can be obtained by the equation Θ ≥ 0 with case of θ = π/2, and then taking (5.22) into the last equation of (5.18) and then we can obtain the following form From Eq. (5.24), once the parameters of the black hole are determined, the range of photon region r c may be r cmin ≤ r c ≤ r cmax .Here, r cmin , r cmax are the minimum and maximum positions of the photon regions, respectively.It can be seen that there are maximum and minimum values of the photon sphere.

The shadow of the Kerr-like black hole in a dark matter halo
In the previous subsection, we calculated the geodesic equations for the motion of lightlike particle and the photon regions in a dark matter halo.Next, we will study the black hole shadow of the Kerr-like black hole in a dark matter halo based on these geodesic equations, and compare them with that of the Kerr black hole.For an asymptotically flat spacetime, we can easily calculate the shadow of the black hole in celestial coordinates, and the formula of the corner radius of the shadow has the same form as that of the shadow corner radius of the black hole in Minkowski spacetime [60], that is where, (r 0 , θ 0 ) is the position of the observer, and the motion of the light particle is described by dϕ/dr and dθ/dr from Eqs. (5.17).Here, we choose to observe the shape of shadow of the Kerr-like black hole on the equator plane, that is, θ 0 = π/2.Observed in the equatorial plane, the shadow of Kerr-like black hole is the outgoing light perpendicular to the rotation axis of the black hole received by an observer at a certain distance.For more details on the introduction of black hole shadow, we can refer to this brief review in Ref. [14].This viewing angle can provide some interesting and important information, because from this angle, we can see the impacts of the rotation parameter a of the black hole on the shadow of the black hole.Then the above equation can be written as (5.26) Taking (4.14) and (5.22) into account, we can obtain the Celestial coordinates in a dark matter halo, that is (5.27) where, g ′ (r) denotes dg(r)/dr and g(r) = f (r) can be obtained in (2.6) in this work.In the Figure 3, we give the shadow of the Kerr-like black hole and the Kerr black hole with the different rotation parameter a for observers located at equatorial plane.At the same time, we also give the shadow of the non-rotating black hole, that is, the Schwarzschild black hole, and the critical value of its shadow is r sch = 3 √ 3M BH [14,85].Our results show that with the increasing of the rotation parameter a, the degree of distortion of the black hole shadow gradually increases.Besides, in Figure 4, we also study the impact of different dark matter parameters on the black hole shadow, and compare them with the Kerr black hole (ρ e = 0).We find that the shape of the black hole shadow in a dark matter halo is similar to the Kerr black hole, which also may show that the dark matter has little impact on the shape of the black hole shadow.On the other hand, the results in Figure 1 show that the horizons of a black hole in a dark matter halo is slightly larger than the GR black hole, but the difference between them is very small, which is a challenge for astronomical observation.But what may be certain is that the shadow of a dark matter black hole is slightly larger than that of a GR black hole.Therefore, when plotting Figures 3 and 4, in order to ensure the integrity of the shadows, we did not show these small differences as the right panel of Figure 1 like.In the next subsection of discussion, we can quantitatively analyze these differences between them from another perspective.

The shadow observables of the Kerr-like black hole in a dark matter halo
In the previous subsection, we introduced the shadow shape of a black hole in a dark matter halo.These black hole shadows have the potential to provide important information about dark matter black holes.Therefore, in this subsection, we start from the observables of black hole shadows and study their impacts on the parameters of the Kerr-like black hole.In other words, we can conversely determine the parameters of the Kerr-like black hole through the shadow shape of the Kerr-like black hole obtained, such as a spin parameter of the Kerr-like black hole.Based on these observables, we may be able to determine or verify the black holes in the universe.
Here, as an example, in Figure 5, we give a coordinate diagram of the Kerr-like black hole shadow (a = 0.999) in a dark matter halo.The point C corresponds to an unstable circular orbit when the observer looks from the equatorial plane.Therefore, in Figure 5, we mainly use the radius of this approximate circle (the black dotted line) to determine the radius R s of the shadow of the black hole.Besides, we also define the distortion parameter d s = α l − α ′ l , which represents the deformation degree of the left side of the shadow.These two parameters R s and d s are the important observables of the black hole shadow.In .The observables R s , δ s and K s separately as a function of the rotation parameter a of the Kerr spacetime (left panels) and the Kerr-like black hole (middle panels) in the galactic center of M87.The right panels are the differences of these observables between Kerr black hole and Kerr-like black hole.The main calculation parameters are M = 1, α = 0.16, ρ e = 6.9 × 10 6 M /kpc 3 , r e = 91.2kpc and the step size of the rotation parameter a is 0.03 from 0 to 1.We have converted these main calculation parameters to the black hole units before plotting. .The observables R s , δ s and K s separately as a function of the shape parameter α of the Kerr spacetime (left panels) and the Kerr-like black hole (middle panels) in the galactic center of M87.The right panels are the differences of these observables between Kerr black hole and Kerr-like black hole.The main calculation parameters are M = 1, ρ e = 6.9 × 10 6 M /kpc 3 , r e = 91.2kpc and the step size of the shape parameter α is 0.005 from 0.10 to 0.30.We have converted these main calculation parameters to the black hole units before plotting.
addition, the authors also studied two other observables, the so-called the area of shadow and the oblateness [15,61].Next, we will study these observables strictly following the procedures in Refs.[15,18,[60][61][62].For the convenience of discussion, we introduce the coordinates of some important points in Figure 5.They are the top point (α t , α t ), the right point (α r , 0), the bottom point (α b , α b ), the left point (α l , 0) and the left point (α ′ l , 0) of the approximate circle.C is the center of the approximate circle with a radius of R s .The center C is on the x-axis, and its coordinate is (α r − R s , 0).For the shadow of a Kerr-like black hole, the observables can be defined as the shadow radius R s and the deformation δ s respectively.According to the geometry of shadows, these observations R s and δ s can be written as and Besides, a new observable K s , the ratio of ∆β and ∆α, could be used to fit the data of the M87.This ratio K s has the following form in shadow geometry where, α r = 2R s + α ′ l .The definitions (5.28), (5.29) and (5.30) are of great significance that they establish the relationship between the observables and the parameters of a Kerr-like black hole.Therefore, we can determine some important parameters of a black hole based on the shape of its shadow.These coordinates in above equations can be calculated with the help of (5.27).Next, we first study the impacts of the observables on the rotation parameter of this Kerr-like black hole in the M87, and compared with the case of the Kerr black hole.In Figure 6, we present these three observables separately as a function of the rotation parameter.Our results show that these observables increase with the increasing of the rotation parameter a both in the Kerr spacetime and the Kerr-like black hole.However, the difference of these observables between these two black holes is not obvious.Therefore, we define three observable differences to determine the size relationship between them, that is Particularly, it can be seen from the observable difference ∆R s that the shadow radius of a Kerr-like black hole is slightly larger than that of a Kerr black hole and the difference ∆R s increases with the increasing of the rotation parameter..The results of the black hole shadow radius ∆R s are consistent with the results in subsection 5.1, which shows that the event horizons of a Kerr-like black hole in a dark matter halo are slightly larger than the Kerr black hole.Similarly, the size between the Kerr-like black hole and the Kerr black hole in the observable differences ∆δ s and ∆K s are the same as the ∆R s and the values of the differences are negative.These results indicate that these two observables ∆δ s and ∆K s of the Kerr-like black hole is slightly larger than that of the Kerr black hole.For example, when the rotation parameter is a = 0.8, the observables of these Kerr-like black holes in the M87 are R s = 5.1979167414914803103, δ s = 0.35114055434480332020 and K s = 1.2125479519509213481.While the Kerr black hole with the rotation parameter a = 0.8, these observables are R s = 5.1979167414914775730, δ s = 0.35114055434480303609 and K s = 1.2125479519509211361. Besides, these two observable differences ∆δ s and ∆K s decrease with the increasing of the rotation parameter, and then increase with the increasing of the rotation parameter in the nearly extremal case.Therefore, the observables of a black hole in a dark matter halo are different from those of the Kerr black hole.In other words, with the development of observation technology, if we observe these observables in the galaxy center of M87, it is possible to provide a reliable basis for the existence of the Kerr-like black hole in a dark matter halo.
On the other hand, these observables of the Kerr-like black hole may help to examine the density profile in a dark matter halo.In the density profile of dark matter, the parameter α is the shape parameter of the halo.Therefore, here we study the impacts of observables on the shape parameter α in the galactic center of M87.In Figure 7, we present these observables as a function of the shape parameter α both in the Kerr spacetime and the Kerr-like black holes in a dark matter halo.In the left panels, since the Kerr black hole is independent of the parameter α, these observables do not change as α increases.In the middle panels of Figure 7, we present these observables of the Kerr-like black hole.We find that these observable changes are small with shape parameter, but they all increase with the increasing of the rotation parameter.Therefore, in order distinguish Kerr black hole and Kerr-like black hole, we define these differences (∆R s , ∆δ s , ∆K s ) as discussed above.In the right panels of Figure 7, we present the differences between these observables.We find that, for a particular rotation parameter, these observable differences are trend to zero as the shape parameter increases.This shows that as the shape parameter increases, it becomes difficult to distinguish the Kerr black holes from the Kerr-like black holes.However, the black hole in a dark matter halo must be different from the GR black hole.Therefore, there is a theoretical upper limit to the shape parameter from the perspective of the black hole shadow observables.Based on the numerical results of the black hole shadow observable differences in the right panels of Figure 7, the upper limit we given is approximately α < 0.22, and this is equivalent to verifying the rationality of the authors' hypothesis in [30] based on the density profile of dark matter from the perspective of the black hole shadow observables.Besides, based on these differences, we find that these observables for Kerrlike black holes are always slightly larger than for Kerr black holes under the same shape and rotation parameter.As the rotation parameter increases, these differences ∆δ s and ∆K s gradually increase, but the difference ∆R s is not obvious.Similarly, this method can also be used to constrain and test other parameters of the dark matter model, such as the upper limit of the density of dark matter.Since the basis of our research is based on the mass model of M87, the halo density and characteristic radius of dark matter are uniquely determined.Therefore, we only studied the impacts of shape parameters to black hole shadow observables in a dark matter halo.Overall, it is possible that the observables from the black hole shadow provide a new method to examine the density parameter of dark matter.With the increase of observational data, these observables of the black hole shadows may further enrich and improve dark matter density profiles and dark matter models.

Conclusions and Discussions
In this paper, using the Einasto profile in a dark matter halo, we obtained the metric of the Schwarzschild-like black hole located in the galactic center of M87.And then using the Newman-Janis algorithm, we extended our solution to the case of the Kerr-like black hole, and this solution was satisfied with the Einstein field equations.Besides, we also studied and analyzed the basic physical properties of these black holes in the galactic center of M87 and then compared some results of them with the Kerr black hole.We found that dark matter may have a positive impact on the physical properties of black hole.Our conclusions are as follows: (1) When the black hole is in a dark matter halo, our results show that the dark matter may effectively increase the event horizons of the black hole and expand the range of the ergosphere of the Kerr-like black hole.
(2) In a dark matter halo, for the area of ergosphere of the Kerr-like black hole, that is the outer horizon and outer infinite redshift surface decrease with the increase of rotation parameter a, while the inner horizon and inner infinite redshift surface increase with the increase of rotation parameter a.
(3) In a dark matter halo, the motion equation for neutral particles near the black hole is given by (5.9).In particular, for photons in the equatorial plane, the motion equation of the photons is given by Eq. (5.10).where, R ′ (r) denotes dR(r)/dr.These results indicate that (4.13) is indeed a solution to Einstein field equations.

Figure 1 .
Figure1.The functional image of ∆ as a function of independent variate r in dark matter spacetime (left panel), Kerr spacetime (middle panel), respectively.The function image of ∆ for different rotation parameters a are given in the first two panels.The right panel is their comparison δ∆ as a function of rotation parameter a at the fixed r.The main calculation parameters are M = 1, ρ e = 6.9 × 10 6 M /kpc3 , r e = 91.2kpc.We have converted these main calculation parameters to the black hole units before plotting.

6 aFigure 3 .
Figure 3.The shape of shadow of the Kerr-like black hole and the Kerr black hole both with the different rotation parameter a for observers located at equatorial plane.The black circle in the center of the coordinate axis represents the shadow of the Schwarzschild black hole.The main calculation parameters are M = 1, α = 0.16, ρ e = 6.9 × 10 6 M /kpc 3 , r e = 91.2kpc.We have converted these main calculation parameters to the black hole units before plotting.

) 2 Figure 4 .
Figure 4.The shape of shadow of the Kerr-like black hole with the different dark matter densities ρ e (M ⊙ /kpc 3 ) for observers located at equatorial plane.The black circle in the center of the coordinate axis represents the shadow of the Schwarzschild black hole.The main calculation parameters are M = 1, α = 0.16, r e = 91.2kpc.We have converted these main calculation parameters to the black hole units before plotting.

Figure 5 .Figure 6
Figure 5.The coordinate diagram of the Kerr-like black hole shadow (a = 0.999) in a dark matter halo.Here, we define the radius R s , the distortion parameter δ s of the shadow in the Kerrlike black hole and the coordinates of some specific points.The main calculation parameters are M = 1,α = 0.16, ρ e = 6.9×10 6 M /kpc 3 , r e = 91.2kpc.We have converted these main calculation parameters to the black hole units before plotting.

Figure 7
Figure 7.The observables R s , δ s and K s separately as a function of the shape parameter α of the Kerr spacetime (left panels) and the Kerr-like black hole (middle panels) in the galactic center of M87.The right panels are the differences of these observables between Kerr black hole and Kerr-like black hole.The main calculation parameters are M = 1, ρ e = 6.9 × 10 6 M /kpc 3 , r e = 91.2kpc and the step size of the shape parameter α is 0.005 from 0.10 to 0.30.We have converted these main calculation parameters to the black hole units before plotting.