Shadows of Sgr A$^*$ Black Hole Surrounded by Superfluid Dark Matter Halo

In this paper we construct a black hole solution surrounded by superfluid dark matter (BH-SFDM) and baryonic matter, and study their effects on the shadow images of the Sgr A$^{*}$ black hole. To achieve this goal, we have considered two density profiles for the baryonic matter described by the spherical exponential profile and the power law profile including a special case describing a totally dominated dark matter galaxy (TDDMG). Using the present values for the parameters of the superfluid dark matter and baryonic density profiles for the Sgr A$^{*}$ black hole, we find that the effects of the superfluid dark matter and baryonic matter on the size of shadows are almost negligible compared to the Kerr vacuum black hole. In addition, we find that by increasing the baryonic mass the shadow size increases considerably. This result can be linked to the matter distribution in the galaxy, namely the baryonic matter is mostly located in the galactic center and therefore, by increasing the baryonic matter can affect the size of black hole shadow compared to the totally dominated dark matter galaxy where we observe an increase of the angular diameter of the Sgr A$^{*}$ black hole of the magnitude $10^{-5} \mu$arcsec.


I. INTRODUCTION
It is widely believed that the region at the centre of many galaxies contain black holes (BHs). They are astrophysical objects which perform manifestations of extremely strong gravity such as formation of gigantic jets of particles and disruption of neighboring stars. From theoretical perspective, BHs serve as a perfect lab to test gravity in strong field regime. Recently, the Event Horizon Telescope (EHT) Collaboration announced their first results concerning the detection of an event horizon of a supermassive black hole at the center of a neighboring elliptical M87 galaxy [1]. Due to the gravitational lensing effect, the bright structure surrounding the black hole also known as the accretion disk appears distorted. Moreover, the region of accretion disk behind the black hole also gets visible due to the bending of light by black hole. The shadow image captured by EHT is in good agreement with the predictions of the spacetime geometry of black hole described by the Kerr metric. The studies of the shadow images of various black holes, together with the current and future observations, are expected to provide an important approach to understand the geometric structure of black holes or small deviations from Kerr metric in the strong field regime.
On the other hand, it is believed that up to 90% of matter in the host galaxy of a central black hole con- * Electronic address: kimet.jusufi@unite.edu.mk † Electronic address: mjamil@zjut.edu.cn(corresponding author) ‡ Electronic address: zhut05@zjut.edu.cn sists of dark matter. Thus, it is natural to expect that the dark matter halo which surrounds the central black hole could lead to small deviations from Kerr metric near the black hole horizon. With this motivation, the black hole solutions with dark matter halo and their effects on the black hole shadow has been studied in [2][3][4][5][6][7]. In this paper, we consider a scenario of a BH solution surrounded by a halo, which contains both the dark matter and baryonic matter. For dark matter, we assume a superfluid dark matter model recently proposed in [8] and construct the corresponding black hole solution with the halo of dark matter and baryonic matter. As a specific example we are going to elaborate the shadow images of Sgr A * to estimate its angular diameter with the effect of the halo. One aim here is to study the shadows of rotating black holes surrounded by SFDM and baryonic matter. As we mentioned, shadows possess interesting observational signatures of the black hole spacetime in the strong gravity regime and, in the future we expect that shadow observations can impose constrains on different gravity theories. It is interesting to mention that the characteristic map of a shadow image depends on the details of the surrounding environment around the black hole, while the shadow contour is determined only by the spacetime metric itself. In light of this, there have been extensive efforts to investigate shadows cast by different black hole and compact object spacetimes. The shadow of a Schwarzschild black hole was first studied in [9] and later in [10] and the same for Kerr black hole was studied by [11]. Since then various authors have studied shadows of black holes in modified theories of gravity and wormholes geometries .
The plan of the paper is as follows: In Sec. II, we use the superfluid density profile to obtain the radial function of a spherically symmetric spacetime. In Sec. III, we study the effect of baryonic and SFDM in the spacetime metric, using the exponential and power law profile for the baryonic matter. In Sec. IV, we find a spherically symmetric black hole metric surrounded by SFDM. And then in Sec. V, we first construct the spinning black hole surrounded by dark matter halo based the spherically symmetric spacetime we obtained by applying the Newman-Janis method. With this spinning black hole, we then study null geodesics and circular orbits to explore the shadows of Sgr A * black hole. In Sec. VI we present the main conclusions of this paper. We shall use the natural units G = c =h = 1 throughout the paper.

A. Superfluid dark matter
From the view of field theory, one can describe the superfluid by the theory of a spontaneously broken global U(1) symmetry, in a state of finite U(1) charge density. In particular, at low energy regime the relevant degree of freedom is the Goldstone boson for the broken symmetry of the phonon field ψ. On the other hand, the U(1) symmetry acts non-linearly on ψ as a shift symmetry, ψ → ψ + c. In the non-relativistic regime, at finite chemical potential µ, the most general effective theory is described by [8] L T=0 = P(X) , (1) in which X = µ − mΦ +ψ − ( ∇ψ) 2 /2m. Note that m is the particle mass and Φ represents the Newtonian gravitational potential. The main idea put forward in [8] is that the DM superfluid phonons are described by the modified Newtonian dynamics (MOND) type Lagrangian At first sight the fractional power of X may seem strange if (2) represents for example a scalar field, however from the theory of phonons one can argue that exponent is crucial since it determines the superfluid equation of state. On the other hand, to obtain a mediate a force between baryons and the DM phonons one must have the coupling where α is a constant while ρ B gives the baryonic density. Let us point out that at zero temperature, this superfluid theory has three parameters, namely the particle mass m, a self-interaction strength parameter Λ, and finally the coupling constant α between phonons and baryons. Starting from the action composed of the (2) and (3) one can obtain a MOND type force law. In particular the equation of motion for phonons is the following [8] ∇ · withμ ≡ µ − mΦ. Using the limit ( ∇ψ) 2 2mμ and ignoring the curl term one has with a B being the Newtonian acceleration due to baryons only. While the mediated acceleration resulting from (3) is In this way it can be shown that From this result it is easy to see well known MOND acceleration by identifying From the galactic rotation curves for a 0 it is found a MOND 0 1.2 × 10 −8 cm/s 2 .
As we saw, the MOND type law obtained in the superfluid model is not exact value and only applies in the regime ( ∇ψ) 2 2mμ. More importantly, the total acceleration experienced by a test particle is a contribution of a B , a phonon and a DM .

B. Spacetime metric in the presence of superfluid dark matter
The density profile can be given in terms of dimensionless variables Ξ and ξ, defined by [8] where ρ 0 ≡ ρ(0) is the central density. From the Lane-Emden equation: and using the boundary conditions Ξ(0) = 1 and Ξ (0) = 0, the numerical solution yields The last equation determines the size of the condensate A simple analytical form that provides a good fit is The central density is related to the mass of the halo condensate as follows From the numerics we find Ξ (ξ 1 ) −0.5. Substituting (13) in Eq. (15), we can solve for the central density   Remarkably, for m ∼ eV and Λ ∼ meV we obtain DM halos of realistic size. The mass profile of the dark matter galactic halo is given by To solve the last integral we use Eqs. (10)- (14), and by rewriting the mass profile in terms of the new coordinate ξ. In particular we obtain From the last equation one can find the tangential velocity v 2 tg (r) = M DM (r)/r of a test particle moving in the dark halo in spherical symmetric space-time In this section, we derive the space-time geometry for pure dark matter. To do so, let us consider a static and spherically symmetric spacetime ansatz with pure dark matter in Schwarzschild coordinates can be written as follows in which f (r) = g(r) are known as the redshift and shape functions, respectively.

III. EFFECT OF BARYONIC MATTER AND SUPERFLUID DARK MATTER
In a realistic situation, galaxies consist of a baryonic (normal) matter (consisting of stars of mass M star , ionized gas of mass M gas , neutral hydrogen of mass M H I etc.,), the dark matter of mass M DM , which we assume to be in the form of a superfluid. The total mass of the galaxy is therefore M B = M star + M gas + M H I as the total baryonic mass in the galaxy. But, as we saw by introducing the baryonic matter we must take into account a phonon-mediated force which describes the interaction between SFDM and baryonic matter.
To calculate the radial acceleration on a test baryonic particle within the superfluid core this is given by the sum of three contributions [52] a(r) = a B (r) + a DM (r) + a phonon (r), (25) where the first term is just the baryonic Newtonian acceleration a B = M B (r)/r 2 . The second term is the gravitational acceleration from the superfluid core, a DM = M DM (r)/r 2 . The third term is the phonon-mediated acceleration a phonon (r) = αΛ ∇ψ.
The strength of this term is set by α and Λ, or equivalently the critical acceleration which is nothing but the phonon force closely matches the deep-MOND acceleration a 0 given by Eq. (8). Thus, in total, the centrifugal acceleration of the particle of a test particle moving in the dark halo in spherical symmetric space-time is given by Outside the size of the superfluid "core" noted as R, there should be a different phase for DM matter, namely it is surrounded by DM particles in the normal phase, most likely described by a Navarro-Frenk-White (NFW) profile or some other effective type dark matter profile, most probably described by the empirical Burkert profile (see for example, [53][54][55][56][57][58][59]).
However, as was argued in [52] we do not expect a sharp transition, but instead, there should be a transition region around at which the DM is nearly in equilibrium but is incapable of maintaining long range coherence. In the present work we are mostly interested to explore the effect of the SFDM core and, therefore, we are going to neglect the outer DM region in the normal phase given by the NFW profile.

A. Exponential profile
As a toy model for describing the baryonic distribution, we consider a spherical exponential profile of density given by [52] where the characteristic scale L plays the role of a radial scale length and M B is the baryonic mass. In that case the total mass can be written as a sum of masses of dark matter and the disc. The mass function of the stellar thin disk is written as We observe from these profiles that by going in the outside region of core i.e. r >> L, the mass M B (r) reduces to a constant M B . In that way the last term in (28) becomes independent of r and the empirical Tully-Fisher relation which describes the rotational curves of galaxies. In our analyses, we shall simplify the calculations by considering a fixed baryonic mass M B in the last term in (28). Let us introduce the quantity and use the tangent velocity while assuming a spherically symmetric solution resulting with valid for r ≤ R. In the present work, we shall focues on the HSB type galaxyes, such as the Milky Way galaxy. We can use M B = 5 − 7 × 10 10 M and L = 2 − 3 kpc [60]. Outside the superfluid core DM exists in a normal phase probably described by the NFW profile or Burkert profile.

B. Power law profile
As a second example, we assume the baryonic matter to be concentrated into an inner core of radius r c , and that its mass profile M B (r) can be described by the simple relation [61] in the present paper we shall be interested in the case β = 1, for high surface brightness galaxies (HSB). Making use of the tangent velocity and assuming a spherically symmetric solution for the baryonic matter contribution in the case of HSB galaxies we find Note that the above analyses holds only inside the superfluid core with r ≤ R.

C. Totally dominated dark matter galaxies
There is one particular solution of interest which describes a totally dominated dark matter galaxies (TD-DMG) or sometimes known as dark galaxies. To do this, we neglect the baryonic mass M BM → 0, and therefore a 0 → 0, yielding valid for r ≤ R. This result is obtained directly from Eqs. (32) and (34) in the limit M B = 0. In other words this solution reduces to (23), as expected. Note that the integration constant C can be absorbed in the time coordinate using dt 2 → C dt 2 .

IV. BLACK HOLE METRIC SURROUNDED BY BARYONIC MATTER AND SUPERFLUID DARK MATTER
Let us now consider a more interesting scenario by adding a black hole present in the SFDM halo. To achieve this aim, we shall use the method introduced in [6]. Namely we need to compute the spacetime metric by assuming the Einstein field equations given by [6] where the corresponding energy-momentum tensors T ν µ = diag[−ρ, p r , p, p], where ρ = ρ DM + ρ BM , encodes the total contribution coming from the surrounded dark matter and baryonic matter, respectively. The spacetime metric including black hole is thus given by [6] where dΩ 2 = dθ 2 + sin 2 θdφ 2 is the metric of a unit sphere. Moreover we have introduced the following quantities In terms of these relations, the space-time metric describing a black hole in dark matter halo becomes [6] In the absence of the dark matter halo, i.e., f (r) = g(r) = 1, the indefinite integral reduces to a constant [6] yielding F 1 (r) = −2M/r. In that way one can obtain the black hole space-time metric with a surrounding matter Thus the general black hole solution surrounded by SFDM with F (r) = G(r) are given by (43) As a special case we obtain a black hole in a totally dominated dark matter galaxy given by valid for r ≤ R. Note that M is the black hole mass. In the limit M = 0 our solution reduces to (23), as expected.

V. SHADOW OF SGR A * BLACK HOLE SURROUNDED BY SUPERFLUID DARK MATTER
In this section we generalize the static and the spherical symmetric black hole solution to a spinning black hole surrounded by dark matter halo applying the Newman-Janis method and adapting the approach introduced in [62,63]. Applying the Newman-Janis method we obtain the spinning black hole metric surrounded by dark matter halo as follows (see Appendix A) where we have introduced and identified F 1 = F (r) exp and F 2 = F (r) power , respectively. As a special case we find the Kerr BH limit, In addition we can explore the shape of the ergoregion of our black hole metric (45), that is we can plot the shape of the ergoregion, say in the xz-plane. On the other hand the horizons of the black hole are found by solving ∆ 1,2 = 0, t.e., while the inner and outer ergo-surfaces are obtained from solving g tt = 0, i.e., r 2 F (r) 1,2 + a 2 cos 2 θ = 0.
In Fig. 2 we plot ∆ 1,2 as a function of r for different values of a. In the special case for some critical value a = a E (blue line) the two horizons coincide, in other words we have an extremal black hole with degenerate horizons. Beyond this critical value, a > a E , there is no event horizon and the solution corresponds to a naked singularity.

A. Geodesic equations
In order to find the contour of a black hole shadow of the rotating spacetime (45), first we need to find the null geodesic equations using the Hamilton-Jacobi method given by where τ is the affine parameter, and S is the Jacobi action. Due to the spectime symmetries there are two conserved quantities, namely the conserved energy E = −p t and the conserved angular momentum L = p φ , respectively.
To find the separable solution of Eq. (49), we need to express the action in the following form in which µ gives the mass of the test particle. Of course, this gives µ = 0 in the case of the photon. That being said, it is straightforward to obtain the following equa- tions of motions from the Hamilton-Jacobi equation where R(r) and Θ(θ) are given by and K is the separation constant known as the Carter constant. We observe that the dashed blue curve corresponding to Sgr A * with the above parameters is almost indistinguishable from the black dashed curve describing the Kerr vacuum BH. The solid blue curve corresponds to the case of increasing the baryonic mass by a factor of 10 2 and decrease of core radius by a factor of 10 3 .

B. Circular Orbits
Due to the strong gravity near the black hole, it is thus expected that the photons emitted near a black hole will eventually fall into the black hole or eventually scatter away from it. In this way, the photons captured by the black hole will form a dark region defining the contour of the shadow. To elaborate the presence of unstable circular orbits around the black hole we need to study the radial geodesic by introducing the effective potential V eff which can be written as follows where , gives the exponential and the power law respectively. At this point, it is convenient to introduce two parameters ξ and η, defined as For the effective potential then we obtain the following relation where we have replaced V 1,2 eff /E 2 by V 1,2 eff . For more details in Figure (3) we plot the variation of the effective potential associated with the radial motion of photons. The circular photon orbits exists when at some constant r = r cir. the conditions Combining all these equations it is possible to show that ξ 1,2 = (r 2 + a 2 )(rF 1,2 (r) + 2F 1,2 (r)) − 4(r 2 F 1,2 (r) + a 2 ) a(rF 1,2 (r) + 2F 1,2 (r)) One can recover the Kerr vacuum case by letting a 0 = M B = ρ 0 = 0, yielding To obtain the shadow images of our black hole in the presence of dark matter we assume that the observer is located at the position with coordinates (r o , θ o ), where r o and θ o represents the angular coordinate on observer's sky. Furthermore we need to introduce two celestial coordinates, α and β, for the observer by using the following relations [29] where (p (t) , p (r) , p (θ) , p (φ) ) are the tetrad components of the photon momentum with respect to locally nonrotating reference frame. The observer bases e µ (ν) can be expanded as a form in the coordinate bases (see, [64]) where the quantities E = −p t and p φ = L are conserved due to the associated Killing vectors. If we use ξ = L/E, η = K/E 2 and p θ = ± Θ(θ) we can rewrite these coordinates in terms our parameters ξ and η, as follows [43] and Finally to simplify the problem further we are going to consider that our observer is located in the equatorial plane (θ = π/2) and very large but finite r o = D = 8.3 kpc, we find Note that f 1,2 are given by Eqs. (32) and (34), respectively. We see that due to the presence of SFDM our solution is non-asymptotically flat. In the special case when SFDM is absent we recover f 1,2 → 1, hence the above relations reduces to the asymptotically flat case. Note that for the observer located at the position with coordinates (r o , θ o ) we have neglected the effect of rotation while ξ 1,2 and η 1,2 given by Eq. (61) and are evaluated at the circular photon orbits r = r cir. For more details in Figs. 5 and 6 we plot the shape of Sgr A * black hole shadows. Using the parameter values specifying the SFDM and baryonic matter we obtain almost no effect on the shadow images. In fact one can observe that the dashed curve (exponential profile) and dashed red curve (power law profile) are almost indistinguishable from the black dashed curve (Kerr vacuum BH). However we can observe that by increasing the baryonic mass the shadow size increases considerably (solid blue/red curve). The angular radius of the Sgr A * black hole shadow can be estimated using the observable R s as θ s = R s M/D, where M is the black hole mass where we have introduced where we have used the approximation cos( π r 2 R ) 1 since our solution is valid for r ≤ R. From the circular photon orbit conditions one can show [41] 2 − rF (r) Now if we expand around X in the above function we obtain By solving this equation and considering only the leading order terms one can determine the radius of the photon sphere r ps yielding On the other hand, one also can show the relation [41] The shadow radius R s can be expressed in terms of the celestial coordinates (α, β) as follows Note that in the last equation we have used Eqs. (69) and (75). Using the values: m = 0.6 eV, Λ = 0.2 meV and ρ 0 = 9.2 × 10 −8 eV 4 , in units of the Sgr A * black hole we find R s = 5.196158313M. Using this result, we can find that the angular diameter increases by δθ s = 3 × 10 −5 µas compared to the Schwarzschild vacuum. This result is consistent with the result reported in Ref. [4] where authors estimated that the dark matter halo could influence the shadow of Sgr A * at a level of order of magnitude of 10 −3 µas and 10 −5 µas, respectively. They have used the so called Cold Dark Matter and Scalar Field Dark Matter models to study the apparent shapes of the shadow. A similar result has been obtained for the dark matter effect on the M87 black hole (see, [5]) using the Burkert dark matter profile.

VI. CONCLUSION
In this paper we have obtained a rotating black hole solution surrounded by superfluid dark matter along with baryonic matter. To achieve this purpose, in the present work, we considered the superfluid dark matter model and used two specific profiles to describe the baryonic matter distribution, namely the spherical exponential profile and the power law profile followed by the special case of a totally dominated dark matter case. Using the current values for the baryonic mass, central density of the superfluid dark matter, halo radius, as well as the radial scale length for the baryonic matter in our galaxy, we found that the shadow size of Sgr A * black hole remains almost unchanged compared to the Kerr vacuum BH. This result is consistent with a recent work reported in [5] and also [4]. For entirely dominated dark matter galaxies, we find that the angular diameter increases by 10 −5 µarcsec. This result shows that it is very difficult to constrain the dark matter parameters using the shadow images. Such tiny effects on the angular diameter are out of reach for the existing space technology and it remains an open question if future astronomical observations can potentially detect such effects. That being said, the expected value for the angular diameter in the case of Sgr A * is of the order of θ s 53 µas as predicted in GR.
As an interesting observation, we show that an increase of baryonic mass followed by a decrease of the radial scale length can increase the shadow size considerably. This can be explained by the fact that the baryonic matter is mostly located in the interior of the galaxy, on the other hand, dark matter is mostly located in the outer region of the galaxy. In other words, for the totally dominated dark matter galaxies we observe almost no effect on black hole shadows, but a more precise measurement of the Sgr A * shadow radius can play a significant role in determining the baryonic mass in our galaxy.