Imprints of dark matter on black hole shadows using spherical accretions

We study the possibility of identifying dark matter in the galactic center from the physical properties of the electromagnetic radiation emitted from an optically-thin disk region around a static and spherically symmetric black hole. In particular, we consider two specific models for the optical-thin disk region: a gas at rest and a gas in a radial free fall. Due to the effect of dark matter on the spacetime geometry, we find that the dark matter can increase or decrease the intensity of the electromagnetic flux radiation depending on the dark matter model. To this end, we analyze two simple dark matter models having different mass functions M(r)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {M}}(r)$$\end{document}, with a matter mass M, thickness Δrs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta r_s$$\end{document} along with a dark matter core radius surrounding the black hole. In addition to that, we explore the scenario of a perfect fluid dark matter surrounding the black hole. We show that in order to have significant effect of dark matter on the intensity of the electromagnetic flux radiation, a high energy density of dark matter near the black hole is needed. We also find that the surrounding dark matter distribution plays a key role on the shadow radius and the intensity of the electromagnetic flux radiation, respectively. Finally we have used the relation between the shadow radius and the quasinormal modes (QNMs) to compute the real part of QNM frequencies.


Introduction
Many modern astrophysical observatories point towards the fact that giant elliptical and spiral galaxies contain supermassive black holes (SMBHs) at their galactic centers. For example, observations suggest that at the center of our Milky Way galaxy there is a supermassive black hole with four million solar masses. These supermassive black holes are characterized by huge masses and spin parameter (or angua e-mail: sbhkmr1999@gmail.com b e-mail: kimet.jusufi@unite.edu.mk (corresponding author) lar momenta). According to Einstein's theory of relativity, black holes generically contain a spacetime singularity at their center and an event horizon where the gravity is so strong that nothing even the electromagnetic radiation or light -can escape from it. Due to the exterior spacetime geometry, black hole (BH) can capture light received from nearby stars or accretion disks into bound orbits. In other words, every black hole is characterized by the photon sphere or a collection of light rays orbiting the BH [1]. In particular, the orbit of light is said to be unstable if the photon can fall into the BH or escape to infinity. The most compelling evidence of the existence of black holes is the first shadow images of the SMBH at the center of M87 galaxy Event Horizon Telescope (EHT) collaboration has detected [2,3] and the detection of gravity waves by LIGO [4]. It is interesting that, by monitoring of motion of stars in the Galactic Center and by collecting precise measurements of the orbital motion of such stars we can probe on the nature of black holes and the surrounding spacetime [5][6][7].
It's quite amazing that we can use these precise observations to constrain different physically viable BH solutions and therefore we can test General Relativity and alternative theories of gravity, say by observing small deviations from the Kerr solution. Toward this goal, we can explore the distortion in the shadow images which encodes valuable information about the black hole mass/spin, and also the spacetime geometry around a given black hole solution. Furthermore, the shadow images can be used to test the existence of other exotic objects such as wormholes and naked singularities. In literature, one can find many detailed studies concerning the black hole shadows, including static and spherically symmetric solutions, rotating solutions, and solutions in different theories of gravity . For further studies concerning the effect of black hole geometry on the electromagnetic radiation see Refs. [58][59][60][61].
From the astrophysical point of view, it is well known that dark matter which is assumed to be some form of elementary particle plays a key role in many astrophysical processes, yet is one of the greatest unsolved mysteries. In recent papers, the effect of dark matter on black hole shadow has been investigated [62][63][64][65][66]. In this paper, we aim to study the possibility of identifying dark matter in the galactic center based on the physical properties of the electromagnetic radiation emitted from a optical-thin disk region around the static and spherically symmetric black hole. Furthermore, we shall consider a radiating optically-thin disk of gas at rest and a radiating thin disk of gas in a free fall.
The structure of our paper is laid out as follows: In Sect. 2, firstly we review the dark matter model proposed in Ref. [62] having dark matter with positive and negative energy density around the black hole. In Sects. 3 and 4, we study the shadow images and intensity of the radiation produced by a spherically thin medium described by a gas at rest and infalling gas model, respectively. In Sect. 5, we study a perfect fluid dark matter model surrounding a black hole and the shadow images/intensity of the radiation produced by a gas at rest and an infalling gas model, respectively. In Sect. 6, we explore the connection between the shadow radius and the quasinormal modes (QNMs). Finally, in Sect. 7, we comment on our results.

Model I
Let us start by considering a toy model of a Schwarzschild black hole surrounded by dark matter proposed recently in [62]. One can then assume a piecewise function to impose three domains [62]: where The expression for G(r ) is chosen so that M(r ) and M (r ) are continuous (see Fig. 1). The Schwarzschild metric surrounded by a spherical shell of dark matter described by Eq.
(2) can be written as follows in which the metric function f (r ) reads In the present work, we shall consider the nontrivial case where the dark matter is concentrated near the black hole such that r s = r h , as shown in Fig. 1. The mass function then is given by m + MG(r ) along with r s ≤ r ≤ r s + r s . From the field equations, we can compute the energy-momentum components and Note that we are going to consider two cases of the surrounding dark matter near the black hole; namely, a dark matter with positive energy density, i.e. when M > 0 (normal matter), and a dark matter with negative energy density i.e. M < 0. Previously, the shadow images of metric (3) and rotating counterpart of metric (3) have been studied in Refs. [62][63][64][65].

Model II: black hole in perfect fluid dark matter
In order to have a more realistic model for dark matter we need to take into account not only the energy density but also the relativistic pressure of dark matter. To do so, let us consider a 3 + 1 gravity theory in the presence of a perfectly fluid dark matter (PFDM) with the action [67][68][69] (see also Refs. [41]) From the action, we get the Einstein field equations as where T DM μν corresponds to the energy-momentum tensor of the perfect fluid dark matter with components (T μ ν ) DM = diag(−ρ, P r , P, P). Let us now investigate a static and a spherically symmetric solutions, for that we need to employ the following metric ds 2 = −e ν dt 2 + e λ dr 2 + r 2 (dθ 2 + sin 2 θ dφ 2 ), adding the ansatz, ν = −λ, from the Einstein's field equations we find and a prime denotes the differentiation with respect to r . If we set T θ θ = T φ φ = T t t (1 − ) with being a constant yielding [69]) It is interesting to note that we can find solutions that are satisfied with the condition λ = −ν, i.e. g tt = −g −1 rr , in the limit T t t → T r r . By using the equation of state for the PFDM to find the solution satisfied with the condition λ = −ν, using ν = ln(1 −U ) with U = U (r ) along with the field equations it can be found [41] We can consider two special cases of the above equation.

Case = 3/2
In this particular case it was found the following solution [41,67]  where r Sh = 2m. Then the corresponding black hole spacetime in perfect fluid dark matter is given by [67,68] with where m is the black hole mass and a is a parameter describing the intensity of the dark matter. In the present work, we are going to assume a > 0. Hence we can write the following approximation if we identify a → 2v 2 0 r 0 , under this scaling we obtain in the large limit r f (r ) 1 + 2v 2 (r ) ln r 2v 2 0 r 0 (20) and the tangential velocity is a function of r , v 2 (r ) = v 2 0 r 0 /r . In order to see the physical interpretation of v 0 let us con-sider the isothermal dark matter density profile given by the following relation along with the mass profile of the dark matter galactic halo given by can explain the flat curve of motions of stars in the galaxy. To see this, let us use the last equation to find the tangential velocity v 2 tg (r ) = M DM (r )/r = v 2 0 which describes a region with the nearly constant tangential speed of the stars. To do so, we can work with a static and spherically symmetric spacetime ansatz with pure dark matter in Schwarzschild coordinates the tangential velocity, one can calculate the radial function F(r ) by the following equation [41] v 2 Assuming a spherically symmetric solution, and solving the last equation we find where r 0 is some constant. From the pure dark matter spacetime it was argued that one can obtain the space-time metric of a black hole surrounded by dark matter halo solving the Einstein field equations (8) and using T ν μ DM = diag[−ρ, p r , p, p] for non-zero energy-momentum corresponding to the pure dark matter space-time metric. One way to include the black hole in our metric is by treating the dark matter as part of the general energy-momentum tensor T μν DM . The space-time metric is chosen as (9) yielding a coefficient functions as follows (see for details [41]) Thus, the above space-time has a black hole in the dark matter halo. Finally, the metric coefficient f (r ) reduces to f (r ) 1 + 2v 2 0 ln(r/r 0 ) which explains the flat curve of motions of stars in the galaxy having v 0 10 −3 . In other words, the last metric function describes the spacetime in the regime where the dark matter effect is not very strong and it is a special case of the metric function (18) having constant velocity v 0 . Thus, the metric function (18) describes a strong regime of dark matter effect near the black hole.

Case
The solution of Eq. (15) in this case reads [67] along with the corresponding metric given by with In this paper, we are interested in the metric that can describe the motions of stars in the galaxy. Therefore let us consider the scaling then the metric function (28) gives This metric describes the extreme case of the dark-matter galaxy with the SMBH at its center and show the corresponding metric function −g tt as a function of r . Furthermore we can consider a Taylor series expansion around α far away from the black hole yielding where it was introduced γ α 2v 2 0 . In what follows we shall consider the shadow images using the dark matter models that we elaborated in this section.

Photon sphere and the shadow radius
Here, we are interested in investigating the shadow of black hole solution surrounded by matter. To do so, we start from Hamilton-Jacobi method for null geodesics in the black hole spacetime written as [70] in which S is the Jacobi action and σ is some affine parameter along the geodesics. If we consider a photon along null geodesics in our spherically symmetrical spacetime surrounded by matter, one can show that the Hamiltonian can be written as Due to the spacetime symmetries related to the coordinates t and φ, there are two constants of motion defined p t = −E and p φ = L, where E and L are the energy and the angular momentum of the photon, respectively. Next, the circular and unstable orbits are related to the maximum value of effective potential in terms of the following conditions Without going into details here, one can now show the following equation of motion Let us consider a light ray sent from a static observer located at a position r 0 and transmitted with an angle ϑ with respect to the radial direction. We, therefore, have [70] cot ϑ = √ g rr Finally, the relation for shadow radius of the black hole as observed by a static observer at the position r 0 can be shown as where r ph represents the photon sphere radius and r 0 is the position of the observer located at a far distance from the black hole. We note here that the spacetime is assumed to be asymptotically flat, hence for the location of the observer in the region far away from the black hole it follows that f (r 0 ) = 1. In general dark matter modifies the photon radius and one should be careful in choosing the correct photon sphere, i.e. the one with the smallest impact parameter. Concerning model I, it was shown that there are two solutions for the photon sphere with opposite sign and this was already analyzed by [62]. On the other hand, for the perfect fluid dark matter case one can obtain the solution for the photon sphere only numerically.

Optically thin radiating gas at rest surrounding the black hole
Let us start by considering a general relativistic analysis for a thin radiating gas at rest in the spacetime of black hole surrounded by dark matter. To do so, we shall consider the following [33]: (i) The proper length corresponding to the coordinate interval dr is √ 1/ f (r )dr. (ii) Every unit of energy emitted in the local frame at r corresponds to √ f (r ) units of energy at infinity. (iii) A time interval dt in the local frame corresponds to an interval dt/ √ f (r ) at infinity. It follows we define the power per unit volume P s , and the emission coefficient per unit solid angle j s in terms of the following relation [33]: where in our case we have g tt = f (r ) along with g rr = 1/ f (r ). The net emitted luminosity in this model, measured at infinity, reads which is the same result as in the in the Newtonian problem. However, it is interesting that only a fraction of the radiation emitted at any given r escapes to infinity and can be thus observed. It is known that the solid angle of the escaping rays is equal to 2π(1 + cos θ) for r outside the photon sphere radius, i.e. r > r ph and 2π(1 − cos θ) for r smaller then the photon radius, ie. r < r ph . Using the fact that the critical rays that barely escape have specific angular momentum equal to that of the photon orbit one can show the relation where b is the impact parameter of the photon orbit. In the case of a black hole spacetime surrounded by dark matter for the escaping rays, the net luminosity observed is found to be [33] For model I, using the parameters M = m = 1, r s = 2m, and r s = 10M, we find In this example, the observer is located at the distance r obs. = 10 6 m. Consider now a more realistic situation when the observer is located within the dark matter halo using the parameters m = 1, M = 10 3 m, r s = 2m, and r s = 10 6 M, and say the observer is located at the distance r obs. = 10 3 m.
In that case, we find We see that when the observer is located inside the dark matter halo, the net luminosity decreases. In the perfect fluid dark matter model using (18), the spacetime is asymptotically flat, hence the net luminosity observed at infinity using the parameters m = 1 along with a = 0.1, and a = 0.2 is found to be and respectively. In this case we see that the net luminosity observed at infinity is considerably higher compared to the Schwarzschild vacuum solution. Moreover, this strong effect shows that metric function (18) describes the spacetime near the black hole surrounded by dark matter with relatively high energy density. The increase in luminosity shows that, the gravitational deflection of light rays when the perfect fluid dark matter is present, is considerably smaller as compared to the Schwarzschild vacuum black hole. In other words, in this case, a larger fraction of the emitted radiation at any r escapes from the black hole. As a result, we observe a higher value of the net luminosity. Finally, using (30) we can find the net luminosity observed at the distance r obs = r 0 = 10 6 using the parameters m = 1 along with α = γ = 0.1, given by Compared to the Schwarzschild vacuum solution, we observe a very small increase in the net luminosity. This shows that the gravitational deflection angle of light in this dark matter model is slightly decreased compared to the Schwarzschild vacuum black hole. We can say that the metric function (30) describes the spacetime around black hole with lower concentration of dark matter while the metric function (18) describes the dark matter effect in the galactic center where the effect is stronger. For example, one can use (30) to explain the flat curve for the observed motions of stars in the outer part of galaxies, while we can use the metric function (18) to describe the motion of stars in the galactic center. Note that in all cases, there is a tiny contribution in the intensity plots even for the region r < r ph . In fact, as we are going to see, the intensity even inside the shadow does not go to zero but has a small finite value. Therefore, if we integrate from the horizon to the photon sphere we end up with a very tiny effect in the luminosity which is included in the shadow/intensity plots. Next, to determine the image observed by a distant observer, we can use the observed specific intensity at the observed photon frequency which can be found by integrating the specific emissivity along the photon path where g = ν obs /ν e is the redshift factor, ν e is the photon frequency as measured in the rest-frame of the emitter, j (ν e ) is the emissivity per unit volume in the rest-frame of the emitter, while dl prop = k α u α e dλ is the infinitesimal proper length as where k μ is the four-velocity of the photons, u α e four-velocity of the accreting gas emitting the radiation, u μ obs = (1, 0, 0, 0) represents the four-velocity of the distant observer, and finally λ is the affine parameter along the photon path γ . Here γ in the integral refers that the integral has to be evaluated along the path of the photon (null geodesics). Using the rest spherical accretion model, the redshift factor is evaluated from g = f (r ) 1/2 . For the specific emissivity we assume a simple model in which the emission is monochromatic with emitter's-rest frame frequency ν , and the emission has a 1/r 2 radial profile: where δ is the Dirac delta function. The proper length can be written as If follows that the specific intensity observed by the infinite observer is where dr dφ = ±r 2 is obtained from the equation of motion obtained for the light ray that moves on the equatorial plane.

Optically thin radiating and infalling gas surrounding a black hole
In this section, we consider a more realistic model, namely an optically thin, radiating accretion flow surrounding the object and then use a numerical technique (Backward Raytracing) to find the shadow cast by the radiating flow. The calculation of the intensity map of the emitting region requires some assumption about the radiating processes and emission mechanisms. The observed specific intensity I ν0 at the observed photon frequency ν obs at the point (X, Y ) of the observer's image (usually measured in erg s −1 cm −2 str −1 Hz −1 ) is given by [36] I obs (ν obs , Here we are considering a simplistic case of the accreting gas. We assume that the gas is in radial free fall with a fourvelocity which in a static and spherically symmetric case reduces to where The four-velocity for the photons has been found in the previous section. To ease our further calculations, we find a relation between the radial and time component of the fourvelocity where the sign +(−) is when the photon approaches (goes) away from the massive object. The redshift function g is therefore given by For the specific emissivity we assume a simple model in which the emission is monochromatic with emitter's-rest frame frequency ν , and the emission has a 1/r 2 radial profile: where δ is the Dirac delta function. The proper length can be written as Integrating the intensity over all the observed frequencies, we obtain the observed flux In Fig. 2 we show the shadow images and the corresponding intensities using the rest gas model as seen by a distant observer for the Schwarzschild black hole surrounded by dark matter using the mass function (1). We see that having a high density of dark matter near the black hole affects the images as well as the intensity. For instance, it can be observed from Fig. 2 that for a given set of values of the parameters, the shadow size is decreased or increased in comparison to the Schwarzschild vacuum black hole. This effect in principle can be used as a tool to observe very dense dark matter located around black holes. On the other hand, in Fig. 3 we show shadow images and the intensities using the infalling gas model as seen by a distant observer for the black hole surrounded by dark matter using the mass function (1). Again we cam observe that for a given domain of parameters for the dark matter we can distinguish such a black hole compared to the Schwarzschild vacuum black hole. Finally, in Figs. 4 and 5 we show the shadow images and intensities for the black hole surrounded by perfect fluid dark matter using the mass function (18) and (30), respectively. It's worth noting that, the perfect fluid dark matter model is more realistic model compared to the Model I which is just a toy model to describe surrounding matter around black hole (not necessarily dark matter). From the Figs. 4 and 5 it can be seen that the particular model of dark matter has a remarkable role on the final image. In other words, depending on the particular gas model as well as the dark matter model, the shadow images and the intensities are different from the Schwarzschild vacuum black hole. This means that, in principle, shadow images can be used as a tool to observe dark matter around black holes. It is interesting to see how the contribution of the intensity of static gas model obtained from the region inside the shadow region (part of the radiation that has escaped to infinity) is higher compared to the infalling gas model. The closer it the dark matter to the black hole the stronger is the effect on the shadow images and this has to do with the fact that the photon radius is changed due to the surrounding matter.

Backward Ray tracing
Now, when we have the desired relation for the observed flux, we are ready to run our simulation. The technique that is used for making shadows of compact objects is Relativistic Ray tracing. There are two different methods which one can be used to do ray tracing; forward ray tracing and backward ray tracing. In forward ray tracing, the light emitted from the source, gets lensed and, red/blue-shifted from the object and then reaches the observer's eye. Here we follow the path of the photon from the point of emission/source to the observer's eye. In terms of computing, we replace the eye with an image plane made of pixels, so when a photon hits one of those many pixels, it gets lightened up indicating that the value of the pixel at that point is greater than zero. This process is repeated multiple times till the image plane fills up and is properly adjusted. The limitation of this technique is, we assume that the rays always hit the image plane, whereas in the real scenario the photons have a very small probability of actually hitting the eye. Potentially, we'll have to shoot millions and billions of rays to get only a single photon actually reaching the eye. One can argue that instead of shooting in all random directions, just shoot where the eye is. This may optimise the problem to some extent, but it is not important for the direction to be specific for diffused surfaces, which is a matter of light-matter interaction. To sum up, this increases the computational power, which means the integration of the rays will neither be that faster nor efficient. Therefore we move on to the second method which is 'backward ray tracing'. Here, instead of tracing the ray from the source towards the eye, we trace the rays backwards in time, which means from the eye towards the source.
To implement a ray tracer, the following steps are performed: 1. Multiple parallel beams of light with varying impact parameters are shot from the observer at r 0 = 30M towards the object. 2. Due to a defined critical impact parameter, the rays which have b > b crit get deflected and reach the source, while the ones with b < b crit plunge into the singularity. 3. Rays that get deflected travel from the observer towards the emitter and possess a certain turning point r t p . So from the observer till r t p the light turns out to be red-shifted and from r t p till the emitter, it turns out to be blue-shifted and illuminates the region. The impact parameter is related to the r t p of a photon (ṙ = 0, 4. Rays that travel below the b crit have no turning points in turn going into the singularity and have only a high redshifted shadow, which darkens the region, hence giving out the shadow. 5. Intensity is then calculated by integrating Eq. (60) from the observer to the emission point (r em ).
We created a routine in python following the above steps mentioned for ray-tracing. A simplified version of the code using the Schwarzschild spacetime has been previously implemented in the EinsteinPY library [71, 72] which was subsequently used in [73] as well.

Correspondence between the shadow radius and QNMs in dark matter spacetime
In this final section we proceed to investigate another possible effect in the spacetime of dark matter based on the field perturbations. In particular we shall investigate the correspondence between the black hole shadow in dark matter halo and the real part of QNMs frequencies. The perturbation of the black hole during the ring down phase can be studied in terms of the QNMs where one must impose an outgoing boundary condition at infinity and an ingoing boundary condition at the horizon. One can show that QNMs can be written in terms of the real part and the imaginary part representing the decaying modes Previously, it has been shown that in the eikonal limit, the real part of the QNMs frequencies is mathematically related to the angular velocity of the unstable null geodesic [11][12][13]. Furthermore, the imaginary part of the QNMs frequencies is linked with the Lyapunov exponent that determines the instability time scale of the orbits [11] where l is the multipole number, c is the angular velocity at the unstable null geodesic, and λ denotes the Lyapunov exponent. In recent works, it was argued that the real part of QNMs frequencies and shadow radius (if we also include a sub-leading regime to half of its value) are related as follows (see for details [66,74]) in which R s is the radius of the black hole shadow. This equation is accurate in the eikonal limit having large values of mul-  Table 2 Numerical values of the shadow radius and the real part of QNMs frequencies obtained for the black hole surrounded by dark matter using m = 1, r s = 2m and negative energy density having M = −1 tipole number, i.e., l >> 1. In that case we can further simplify if as ω = R −1 s l. It is thus interesting to see if we can use the above correspondence to calculate the QNMs frequencies once we have the black hole shadow radius. As we explained, this geometric-optics correspondence becomes accurate only for large l, but sometimes works well even for the fundamental modes (compared to the scalar field perturbations). Given the fact that fundamental modes are more important from observational point of view, we present our numerical results of the real part for some fundamental QNMs frequencies presented in Tables 1, 2, 3 and 4. From Table 1 we observe that the shadow radius decreases with the decrease in dark matter density (increasing r s ). From the inverse relation between the shadow radius and the QNMs frequencies, this means that the real part of QNMs should increase with the decrease in r s . In fact, we can verify this fact from the numerical values of QNMs given in Table 1. In Table 2, we consider the case of having dark matter with negative energy density. Quite interestingly, in this case we find that the values of the shadow radius are smaller compared to the Schwarzschild vacuum solution and increase with the decrease of energy density. The QNMs frequencies decrease with the decrease of the Table 4 Numerical values of the shadow radius and the real part of QNMs frequencies obtained for the black hole surrounded by perfect fluid dark matter as a function of α. Here we chose γ = 0.1 and r 0 = 10 6  of energy density. This suggest that the black hole should oscillate more rapidly if we increase the amount of negative energy density of dark matter. A similar results are obtained for the perfect fluid dark matter with the metric function (18) and (30), respectively. In both cases the shadow radius is smaller compared to the Schwarzschild vacuum solution, however in the first case the shadow decreases with a and increases with α in the latter case. The fact that the shadow radius is smaller as compared to the Schwarzschild vacuum black hole is explained from the fact that dark matter with negative energy density is involved in the perfect fluid model, provided a > 0 and α > 0, respectively. For more details we refer the reader to the plots of QNMs frequencies and the shadow radius shown in Figs. 6, 7, 8 and 9. Note that the fundamental modes of QNMs and the shadow radius in a static perfect fluid dark matter without accretion matter with the metric function (18) have been investigated in Ref. [66] using the WKB method. Our results in Table 3, are in good agreement with the scalar field perturbations presented in [66].

Conclusions
We have studied the influence of dark matter on the shadow images using the electromagnetic radiations emitted from spherical accretion medium which was assumed to be an optically-thin region surrounding the black hole. We have considered two spherical accretion models: optically-thin radiating gas at rest and a gas in a radial free fall around the static and spherically symmetric black hole. We have shown that due to the effect of dark matter on the spacetime geometry the intensity of the electromagnetic flux radiation is altered compared to the Schwarzschild vacuum black hole. In particular we have investigated two models for the dark matter distribution: Firstly, the case of total mass function which consists a black hole with mass m, surrounded by the dark matter having mass M and positive energy density along with a thickness r s has been analyzed. Secondly, we explore the case having exotic dark matter around the black hole with negative density energy and a thickness r s . In order to have a more realistic model we need to take into account the relativistic pressure of dark matter, therefore we can consider a perfect fluid dark matter surrounding the black hole. It is found that in order to have significant effect of dark matter on the intensity of the electromagnetic flux radiation, a high energy density of dark matter near the black hole is needed. Among other things, we have shown the surrounding dark matter mass function which describes the distribution of the dark matter plays a key role on the effect of shadow images. Furthermore, we argue that if the surrounding dark matter has positive (negative) energy density, the shadow radius and the intensity of the electromagnetic flux radiation increases (decreases), respectively. Finally, we have also studied the correspondence between the QNMs frequencies and the black hole shadow radius. We have shown that the shadow radius as well as the QNMs frequencies can decrease or increase depending on the concentration of dark matter near the black hole having positive or negative energy density.
We also note that the observed image will not depend on the distance units being used, as long as the distance between the observer and the black hole is kept the same. However, a more delicate issue is a comparison between images of different black holes. In such a case, saying that the horizon mass m is the same on both systems or that the ADM mass M ADM is the same on both systems will not be the equivalent statements in general. This means that, the conclusion that the shadow image is larger in one system than another strongly depends on which mass scale is set to be the same on both systems.
Finally, we should mention here that the gas medium in the present paper is optically thin accreting medium, we plan in the near future to extend our result by including a rotating accretion disk and a rotating black hole spacetime. This is of course outside the scope of this article.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data underlying the results of simulations in this article can be shared on reasonable request to the corresponding author.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .