Fuzzy dark matter black holes and droplets

We consider the possibility of having Dark Matter (DM) black holes motivated by the Einasto density profile. This generalizes both the noncommutative mini black hole model and allows DM to enter as the matter constituent which makes up the black hole. We show that it is possible to construct a black hole solution for each value of the Einasto index and for different values of the mass parameter, provided that the we work with the energy-momentum tensor of an anisotropic fluid. In particular, we achieve that by first considering the equation of state (EOS) $p_r=-\rho$. It turns out that the corresponding black hole solution exhibits a horizon structure similar to that of a Reissner-Nordstr\"{o}m black hole and the central singularity is replaced by a regular de Sitter core. We also show that if the previous EOS is replaced by a nonlocal one, it is possible to construct a self-gravitating fuzzy DM droplet but also in this case, the radial pressure is negative. Finally, we contemplate scenarios of different dark matter black holes with moderate mass values which could have formed in galaxies. In particular, we probe the possibility whether such black holes could also be the central galactic objects.


I. INTRODUCTION
To examine the nature of Black Holes (BH) [1] more closely, it would be desirable to infer more about their interior structure, be it from the geometric point of view [2] or probing into the question what kind of matter has contributed to its formation [3]. In view of the fact that baryonic and leptonic matter constitutes only four percent of the content of the universe and secondly recalling that the density profiles of the major component, Dark Matter (DM), grows as we approach the galactic center [4], the question whether the central galactic BH has a close connection to DM is a well-posed problem. Of course, the no-hair theorem [5] prevents us from inferring observationally the inside properties of a BH, but modelling the interior structure of a BH (indeed, theoretical models of the interior BH are quite common in literature [10]) could reveal if a connection between DM and BH is feasible and at the same time be possibly a harbinger of new physics regarding both, the BH and the DM. An attempt in this direction has recently been proposed in [7] where the authors model the central galactic object assuming a DM profile fitted to the outer galactic region. Other possible inter-connections between the two important components in the galactic bulge have been examined in [8] where the authors consider the possibility of the growing of a BH by the DM absorption. This again suggests the BH as a seed for galactic structure. If so, it not unreasonable to think of a BH as made purely from DM. Motivated by this scenario we develop a new model of the galactic central object as a fuzzy BH (or droplet) in close analogy to BH/droplets inspired by non-commutative geometry [9] where the Gaussian matter distribution and the de Sitter EOS play an important role. The above mentioned approach to BH physics can be generalized based on the fact that the Gaussian distribution is a special case of what is known as the Einasto profile of DM. This allows us to establish a possible connection between BH and DM, i.e., in constructing a fuzzy BH we follow the steps of a non-commutative BH with a new density profile, the Einasto parameterization. In this paper, we test the possibility of a connection between the galactic BH and DM at the current stage of the Milky Way. However, in principle, we could also entertain the possibility of a smaller primordial fuzzy BH made out of DM and growing with time by absorbing matter and DM. We will leave such a project to future undertakings.
Before the existence of a supermassive black hole at the centre of the Milky Way [23][24][25], known as Sagittarius A * , was widely accepted, there have been several attempts to construct theoretical models replacing the central black hole with other gravitational objects such as: gravastars [26,27], bosonstars [28,29], naked singularities [30][31][32], burning disks [33], quantum cores (Ruffini-Argüelles-Rueda model) [34,35] and gravitationally bound clamps of dark matter relying on the exponential-sphere density profile [7,36,37]. It is worthwhile noticing that our information on the central galactic BH comes mainly from observing star orbits [23,24]. In any model beyond the standard BH, it is obligatory to pay attention to the condition that the effective potential of the model, be it a droplet or a BH, resembles the standard BH effective potential, at least in the region of the bound orbits, i.e., around the local minimum of the potential (see, e.g, [7]). In this respect, we will look for suitable parameters of the Einasto profile to model a fuzzy BH, a DM droplet or simply a BH made out of DM at the center of the galaxy in such a way that the observational data are confirmed by our model. Preferably, we will favour parameters already fitted to observational data, but we will also entertain the possibility that at the center of the galaxy the density profiles exhibits a different behaviour. To reach our goals, we will couple the aforementioned profile with an energy-momentum tensor for an anisotropic fluid and an equation of state of the form p r = −ρ. Such an equation of state is quite common in the physics of BH [10]. In the present work, it leads to several regular BH models such that for each value of the parameters ξ and h entering in the Einasto profile, the mass parameter can be tuned so that a black hole will be present at the centre of a galaxy. This black hole has a horizon structure reminiscent of that we observe in the case of a Reissner-Nordström geometry. Furthermore, the inner region of the black hole does not exhibit a curvature singularity at r = 0 which is instead replaced by a regular de Sitter core. We also compute the Hawking temperature for the Einasto inspired black hole: we discover that the black hole increases its temperature, as the horizon radius shrinks, until the temperature reaches a maximum after which the temperature decreases sharply and vanishes exactly at the radius of the extremal black hole. But we will also demonstrate that a DM droplet would eventually result in an effective potential in accordance with observational data on the orbits. In the second model, we consider a nonlocal equation of state for an anisotropic fluid. We observe that also in this case a negative pressure cannot be avoided. In particular, we construct a self-gravitating fuzzy DM droplet regular at the origin whose effective potential allows bound states for massive particle. Also in this case there is no central singularity at r = 0.
The paper is organized as follows: In Section II, we introduce the Einasto profile and some relevant formulae needed in the sections to follow. In Section III, we investigate DM objects assuming a de Sitter-like EOS. In section IV we use a nonlocal equation of state to show that a fuzzy DM droplet made of an anisotropic fluid allows for stable orbits of massive particle. Finally, we present our conclusions in section V.

II. THE EINASTO DENSITY PROFILE
In his 1969 seminal paper [11], Einasto showed that any realistic model aiming to give a faithful description of a galactic system should be characterized by certain descriptive functions such as the cumulative mass profile, the gravitational potential, the surface mass density etc., all satisfying a given set of constraints. Since such descriptive functions are integrals of the density profile ρ = ρ(r) with r a radial variable, it is natural to think that the most fundamental descriptive function of a galactic model is represented by the density profile itself which should exhibit the following properties it is a smooth and decreasing function that vanishes asymptotically at space-like infinity; 3. certain moments associated to ρ such as the central gravitational potential, the total mass, and the effective radius of the system must be finite; We recall that the Einasto density profile [11], which is generally adopted to describe cold DM halos [12, 16-18, 20, 38, 39] as well as the surface brightness of early-type galaxies and the bulges of spiral galaxies [40][41][42][43][44][45][46][47][48][49], is represented by the function where ξ is the Einasto index, r s the radius of the sphere enclosing half of the total mass, ρ s the mass density at r = r s and d ξ a numerical constant controlling that r s is indeed the half-mass radius. In the context of DM halos, the above density is also rewritten as [39] ρ(r) = ρ −2 exp −2ξ r r −2 where ρ −2 and r −2 are the density and the radius at which the density profile behaves like r −2 , i.e. d ln ρ/d ln r = −2.
If we introduce the central density and the scale length as in [22], it is straightforward to verify that the density profile becomes and by adjusting the triple of parameters {ρ 0 , h, ξ}, it is possible to model a variety of astrophysical objects. For instance, we have 4.54 ξ 8.33 for DM haloes with masses in the range of dwarfs to clusters [12], ξ ∼ 5.88 for galaxy-sized haloes [17,18], ξ ∼ 4.35 for cluster-sized haloes in the Millenium Run [13] and ξ ∼ 3.33 for the most massive haloes for the Millenium Simulation [13,18]. Since formula (5) is equivalent to (1) and (2), there is no loss in generality if we work with the expression of the density profile given by (5). Furthermore, the mass function m and the gravitational potential Φ can be computed by solving the following ODEs obtained from the Newtonian equations of hydrostatic equilibrium, namely where G N denotes Newton's gravitational constant. As in [22], we immediately find that the total mass M associated to the Einasto density profile is where Γ denotes the Gamma function. The above relation allows to express the central density ρ 0 in terms of the total mass, and hence, we can rewrite (5) as Furthermore, a straightforward integration of the first equation in (6) leads to the following cumulative mass distribution where γ denotes the lower incomplete Gamma function [50]. For a detailed analysis of the classical gravitational potential we refer to [22].

III. FUZZY BLACK HOLES
In this section, we show that it is possible to construct black hole solutions from the Einasto density profile. Let us suppose that the mass density of a static, spherically symmetric, smeared gravitational source of total mass M be modeled by the density profile (8) which contains as a special case (ξ = 1/2 and h = √ θ) the Gaussian profile adopted by [53] in the derivation of the noncommutative geometry inspired Schwarzschild black hole. Furthermore, we consider the following ansatz representing a static, spherically symmetric manifold. We want to find the unknown function g 00 appearing in (10) so that the above line element is a solution of the Einstein field equations coupled to the energy-momentum tensor of a static, anisotropic fluid and in the limit r/h → ∞ such a solution goes over into the usual Schwarzschild metric. As in [53] we introduce the energy-momentum tensor of a static, anistropic fluid with density source (8) given by where p r and p ⊥ are the radial and tangential pressures, respectively, and we consider the Einstein field equations for the line element (10). If we proceed as in [67], we can use the conservation equation T µν ;ν = 0 with µ = 1 in the (µ, ν) = (2, 2) equation coming from (12) to obtain the Tolman-Oppenheimer-Volkoff equation, i.e. the general relativistic hydrostatic equilibrium equation given by dp r dr where the mass function is defined as in the case of the Einasto profile. On the other hand, from the equation T 1ν ;ν = 0 we get − dp r dr = 1 2 g 00 dg 00 dr (p r + ρ) and if we require that then (15) can be solved for p ⊥ giving This procedure ensures that the conservation equation for the energy-momentum tensor is identically satisfied. It is not difficult to verify that the tangential pressure vanishes at and takes its maximum value at Moreover, we have p r (0) = p ⊥ (0). Finally, we observe that equation (13) is trivially satisfied once p r and p ⊥ have been chosen as in (16) and (17), respectively. This approach shows that the Einasto matter distribution ρ may describe a self-gravitating droplet of anistropic fluid. If we consider the (µ, ν) = (0, 0) or equivalently, the (µ, ν) = (1, 1) equations in (12) together with (16) and the requirement that the metric goes over into the Minkowski metric asymptotically at infinity, we end up with the line element where the mass function is given by (14). By means of 6.5.3 in [52] it is straightforward to verify that in the limit r/h → ∞ the metric (20) reproduces the classic Schwarzschild metric. Furthermore, if we consider the g 00 component of the metric as a function of r/h and we introduce the scaled mass µ = M/h, it is possible to show that there exists a value of µ, say µ 0 , such that g 00 has a double root at x 0 = r 0 /h. The numerical values of the extremal rescaled mass µ 0 and the degenerate horizon x 0 for different values of ξ have been displayed in Table I. They have been obtained by expanding the lower incomplete Gamma function in the expression for g 00 with the help of 6.5.29 in [52] where the first hundred terms in the expansion have been considered. Moreover, if µ > µ 0 there exist two distinct horizons r 1 and r 2 , and no horizon if 0 < µ < µ 0 . Figure 3 displays the plot of g 00 which exhibits the behaviours described above, i.e. 2 horizons, 1 horizon and no horizon. Moreover, in the extreme and non extreme regimes, that is µ ≥ µ 0 ,  the behaviour of the metric coefficient g 00 for r → 0 can be obtained by means of 6.5.29 in [52] and we find that The result above signalizes that differently as in the Schwarzschild black hole where there is a singularity at r = 0, the central region is represented by a regular de Sitter core. Hence, the Einasto density profile coupled with an energy momentum tensor for an anisotropic fluid cures the problem of the central singularity provided that an EOS for the radial pressure is assumed to be of the form p r = −ρ. Finally, in the case µ < µ 0 there is no event horizon and also no naked singularity because the central region around the origin is described by a de Sitter geometry. We conclude this section by considering the Hawking temperature for this new class of black holes. The black hole temperature can be computed from the formula [53] T H = 1 4π where r H represents the position of the event horizon and the total mass M has been expressed in terms of r H by using the horizon equation r H = 2m(r H ). Note that in the case r H /h ≫ 1 the expression above reproduces the usual result T H = (4πr H ) −1 . The scenario emerging from Fig. 4, where we plotted the temperature T H as a function of r H , is that an Einasto inspired black hole increases its temperature, as the horizon radius shrinks, until T H reaches a maximum after which T H decreases sharply and vanishes exactly at the radius of the extremal black hole, that is at r H = r 0 . Furthermore, in the case of an extreme black hole the Hawking temperature must be identically zero because the metric component g 00 has a double root at r = r 0 . Hence, instead of observing a blow-up behaviour of the BH temperature, we find that the evaporation process leads to a zero temperature extremal black hole whose final configuration is entirely controlled by the Einasto parameter ξ, the scale factor h and the black hole mass. As already pointed out in [54], a final configuration characterized by a finite temperature inhibits any relevant back reaction, i.e a self-interaction of the radiated energy with its source. This implies that our solution is stable versus back reaction and can describe the entire black hole life until the final configuration. The presence of an inner Cauchy horizon may be a source of concern, in the sense that the inner region of our black holes is unstable, however one may proceed as in [54] to show the stability of the Einasto inspired black hole interior. At this point, a remark on nomenclature is in order. If there is no horizon, we call the object a fuzzy droplet. If at least, one horizon develops, we name it a fuzzy BH. We end this section with a remark on negative pressure, singularity theorems and the involved scales. For instance, it is tempting to attribute the negativity of pressure to some quantum effects which would limit the applicability of the scales one is using in a model. First, we should notice that with this subject we touch the forefront of science and different interpretations exist in literature. For instance in [55] the possibility of a phase transition to a negative pressure p = −ρ EOS is considered as existing already in Einstein gravity. The inclusion of quantum theory would only would complete the picture, but is not a necessary ingredient. In such a picture the question of scales does not really arise. In [56] a quantum "Macro-Quantumness" is advocated with the claim that the quantum effects for the macroscopic black holes are more important than suggested by means of a semi-classical reasoning. It seems that one should treat the black hole as macroscopic quantum phenomenon and again no scales limitations arise. In short, the result is that the assumption that a black hole of macroscopic size can always be described classically leads to a contradiction and the classical description must break down on macroscopic scales. This was initially suggested in the seminal work of [56] by means of the so-called quantum N -portrait, which has been further developed in [57][58][59][60][61]. Moreover, in [63] and [64] a long distance origin for Hawking radiation is considered and a "quantum atmosphere" assumed. This again would hint towards macroscopic quantum phenomena for any size black hole.
Inter alia, we would like to add one more possible interpretation. Our precise EOS, p r = −ρ is coined according to the physics with a cosmological constant. In such a case, it is interesting to notice that the cosmological constant can violate the singularity theorems or the assumptions underlying them. By this token, the regular black hole which we constructed in this section is not limited by scales due to the possible quantum origin of the EOS. To substantiate our claim, we give another example which has to do with violation of singularity theorems in cosmology due to the cosmological constant Λ. We refer to [65] where the authors show that for Λ > Λ crit there is no initial singularity. This violates the global cosmological singularity theorems. In particular, no quantum mechanics is involved. It is rather an effect of the cosmological constant. We think a similar mechanism happens in the construction of regular black hole solutions under the de Sitter EOS coined after the physics with the cosmological constant. If so, the physics is again not really restricted by scales dictated by quantum mechanics. We can look at it from yet another point of view. In general, negative pressure is a concept taken seriously in physics [66]. As pointed out in [66] a negative pressure is not forbidden by the laws of thermodynamics. It is considered mostly in liquids and has an underlying mechanism, which is not necessarily attributed to quantum mechanics. Again we can argue that as such it is not a local phenomenon limited to quantum mechanical scales. A quick comparison with regular black hole physics where a negative pressure is used, tells us that we can start with a de Sitter EOS as done here and in [10] or it emerges naturally like in [55] where the effect is attributed to General Relativity. In the next section with will construct yet another model where a self-gravitating droplet emerges with a negative pressure.

A. The effective potential
If we insist that the black hole solution derived in Section III sits at the centre of our galaxy, all observations should be the same. To this purpose, we study the problem whether our model of a diffuse dark matter black hole is able to fit the central galactic black hole in the Milky Way whose mass and Schwarzschild radius are M BH = 4.1 · 10 6 M ⊙ and R BH = 2G N M BH /c 2 = 17.4 R ⊙ = 3.92 · 10 −7 pc, respectively [23,24]. In order to do that, we need to find estimates for the Einasto parameter ξ and the scaling factor h. This is done in two steps. First of all, we impose that the total mass M entering in the line element (10) through the metric coefficient g 00 coincides with M BH . Secondly, we require that the mass function m provides a good approximation for M BH when it is evaluated at the minimum r min of the Schwarzschild effective potential for a massive particle. More precisely, we exploit the freedom to force that In the analysis to follow, it is convenient to rewrite the above condition in the equivalent form where we made use of (14). As we will soon realize, it will turn out that the above condition not only ensures that the Schwarzschild effective potential and the effective potential of our diffused gravitational object share the same minimum but they both also agree in a large neighbourhood of it and asymptotically away (see for instance Fig. 5). The effective potential V ef f for the problem at hand is obtained from the geodesic equation. More precisely, following the same procedure as in [67], we can bring the radial equation into a form of an energy conservation equation, namelẏ where the dot means differentiation with respect to the proper time or an affine parameter, depending whether we consider the case of a massive or a massless particle. By means of equation (25.27) in [67] we immediately find that the effective potential associated to the geometry described by the line element (10) is where m p denotes the mass of a test particle and ℓ is its total angular momentum per unit mass. At this point, it is also useful to recall that the effective potential in the case of the Schwarzschild metric is Let r s = 2M BH . If we rescale the radial variable and the angular momentum per unit mass as r * = r/r s and L = ℓ/r s , the Schwarzschild effective potential in the massive case becomes and the event horizon is now located at r * = 1. Furthermore, it exhibits a minimum and a maximum at provided that L > √ 3. Introducing the same rescaling for (24) and for the effective potential (26) in the massive case yields and ∆γ := 1 Γ(3ξ) γ 3ξ, The above condition is an inequality in the free parameters H and ξ. To show that its solution set is non empty, we will first consider different choices of H so that the corresponding scale factors h = r s H have the same orders of magnitude of the scaling factors appearing in [11] and [39]. For each choice of H we solve the inequality (31) with respect to the parameter ξ. Since the particular value of r * min depends on the rescaled total angular momentum L, the procedure outlined above requires that we also fix L. For instance, in [11] the scaling factor for a DM halo is h E = 2.121 · 10 −9 Kpc. Hence, if we choose H = 10 the corresponding scaling factor in our model is h = 3.92 · 10 −9 Kpc. To find out which values of ξ will satisfy (31), we consider different values of L and r * min . In the case L = 2 and r * min = 6, we find numerically that ∆γ < 0 in the range 10 −6 ≤ ξ ≤ 13 signalizing that the inequality (31) cannot be satisfied. The situation changes if we increase the value of L. If L = 3 with r * min = 9 + 3 √ 6, it turns out that ∆γ < 0 for ξ < 0.32. If L = 5 with r * min = 25 + 5 √ 22, any ξ < 0.80 will do the job while for L = 100 with r * min = 10 4 + 10 2 √ 9997 it is necessary that ξ < 2.73. Does our model predict a fuzzy BH or a fuzzy droplet when H = 10 and ξ is chosen so that (31) is satisfied? To answer this question, we observe that in geometric units r s = 2M BH so that M BH = r s /2 and the rescaled mass parameter µ entering in our model will be fixed according to where in the last step we used the relation h = r s H. Moreover, by means of the rescaling r * = r/r s together with the expansion 6.5.29 in [52] we can rewrite the metric coefficient g 00 according to As we can see from Fig. 6, the equation g 00 (r * ) = 0 does not admit any real root and therefore, this model predicts a fuzzy droplet. More recent estimates of the DM density using the rotation curve of the Milky Way using different Galactic mass models together with certain DM and baryonic density distributions has been obtained by [39]. In particular, [39] studied two baryonic models: the model B1 where a combination of Plummer's and Miyamoto-Nagai's potentials has been assumed and the model B2 which relies on the assumption of two different double exponential profiles and the Hernquist profile. Three different DM spherical halos were tested where one is described by the Einasto profile. In the B1 model coupled to the Einasto profile, the scaling factor is h B1 = 3.89 · 10 −11 Kpc while the B2 model predicts h B2 = 1.426 · 10 −5 Kpc. Let us first consider the B1 model. In this case, we need to fix H = 0.1 so that the scaling factor predicted by our model is of the same order as h B1 . More precisely, we have h = 3.92 · 10 −11 Kpc. Proceeding as before, we find that, in order for (31) to be satisfied, ξ < 1.67 for L = 2, ξ < 1.98 for L = 3, ξ < 2.31 for L = 5 and ξ < 3.98 for L = 100. As it can be seen in Fig. 7, we observe that (30) is well approximated by the Schwarzschild effective potential in a large neighbourhood of the minimum also in the case H = 0.1. To the value H = 0.01 there corresponds a rescaled mass parameter µ = 5. A close inspection of Fig. 8 shows that differently as in the case H = 10 we have a more complex scenario. If ξ < 1.2865, we have a dark matter black hole with two distinct horizons while for ξ = 1.2865 an extreme black hole with radius r e = 2.33 · 10 −7 pc < R BH = 3.92 · 10 −7 pc. Finally, if ξ > 1.2865, there is a diffused dark matter droplet without horizon. The same scenario occurs if we further reduce the value of H. For instance, if H = 0.02, the extreme value of the Einasto parameter discriminating between a fuzzy black hole and a fuzzy droplet is found to be ξ = 1.9093. Regarding the B2 model in [39], it is necessary to fix H = 10 5 . It turns out that it is not possible to find any value of the parameter ξ such that the inequality (31) is satisfied. This means that the matching procedure at the minimum of the effective potential which ensures at r = r min that the mass function approximates M BH according to the condition (23) cannot be applied.

IV. DIFFUSED SELF-GRAVITATING DARK MATTER DROPLETS FROM A NONLOCAL EQUATION OF STATE
In the previous section, we assumed an equation of state for the radial pressure of the form, p r = −ρ, and an anisotropic fluid with an additional tangential pressure because, if we would have insisted on a hydrostatical equilibrium, expressed through the Tolman-Oppenheimer-Volkov (TOV) equation, and an energy-momentum tensor of a perfect fluid, this approach would have led to an over-determined system of equations as the density ρ is already Kpc in [39]. The minimum is located at r * min = 9 + 3 √ 6 ≈ 16. 35 where V ef f,S (r * min ) = −0.0147 and V ef f (r * min ) = −0.0147. Both potentials share the same minimum and V ef f,S is a good approximation of V ef f in a neighbourhood of the minimum and asymptotically away.
FIG. 8: Plot of the metric coefficient g00 given by (33) in the case H = 0.1 and µ = 5 for ξ = 1 (solid line), ξ = 1.2865 (dot line) and ξ = 1.4 (longdashed line). Note that for ξ = 1 the event horizon is at r * = 1 and coincides with that of a Schwarzschild black hole. The extreme black hole corresponds to ξ = 1.2865 and its event horizon is located at r * e = 0.5947 or equivalently at re = 2.33 · 10 −7 pc. assigned a priori and the pressure is fixed by the EOS. This allowed to show that starting with the Einasto density profile, it is possible to construct a fuzzy black hole or a diffused self gravitating droplet depending on the particular value of the rescaled mass parameter µ. In the present section, we offer a further example through a different EOS signalizing that the regularity of a fuzzy black hole or a fuzzy self-gravitating droplet seems to require as a main feature that the radial pressure has to be negative at least on a subset of the positive real line. To this purpose, we need to fix a certain EOS and an energy-momentum tensor. Since the Einasto energy density has a diffused profile, we would expect that any change in the radial pressure should take into account the effects of the variations of the energy density within the entire volume. Hence, it seems reasonable to think that nonlocal effects may play a certain role when we work with such a distribution. For instance, in order to take into account nonlocality, we could imagine as in [71][72][73] that the components of the energy-momentum tensor besides displaying a dependence on the spacetime event it also exhibits a functional dependence by averaging the energy density over the enclosed configuration . Without further ado, let us derive the complete solution of the gravitational field equations for an Einasto inspired anisotropic fluid described by a nonlocal equation of state of the  [11] and central object is modelled in terms of a droplet (see also Fig. 6). The maximum of V ef f,S is located at the radius of the photon sphere r * γ = 3/2 while the event horizon of the Schwarzschild black hole is r * s = 1. Note that the effective potential given by (30) does not exhibit a maximum and therefore, the droplet does not possess a photon sphere. form [71][72][73] Since we are interested in matter configurations at hydrostatic equilibrium, we can introduce an effective size R of the object by the condition p r (R) = 0. The numerical value of R can be found by plotting h 2 p r versus y = (r/h) 1/ξ . To this purpose, it is convenient to introduce the mass parameter µ = M/h so that (34) becomes where in the last step we made use of 6.5.29 in [52]. For a list of numerical values of R * = R/h we refer to Table II. Moreover, in Figure 11 we plot the radial pressure to explicitly demonstrate that it is indeed positive in a region of finite size R but negative outside. This however does not mean that the gravitational object will have a finite radius because the energy density does not vanish in the region r > R. In order to proceed further, we consider also in this case a spherically symmetric static matter distribution represented by the Einasto density profile ρ but differently as in the previous section we now assume the For different choices of the parameter ξ we present some typical values of the radial distance R * = R/h at which the radial pressure vanishes. The case ξ = 7.072 corresponds to the Dark Matter case studied in [11].
The numerical values has been obtained by setting equal to zero the expression in the brackets appearing in (35) and choosing k = 100. Note that R * grows linearly in the parameter ξ.   Table II. following ansatz for the line element If we suppose that the energy-momentum tensor in the region occupied by the matter distribution is again given by (11), the Einstein field equations G µν = −8πT µν together with the conservation equation T µν ;ν = 0 with µ = r give rise to the following coupled system of ODEs dp r dr + p r (r) + ρ(r) A(r) With the help of equations (37) and (38) it is straightforward to check that equations (39) and (40) represent the same equation. Hence, we will restrict our attention to the differential system given by (37), (38) and (40) and in order to avoid to work with an under-determined system (there are three equations for the four unknown functions A, B, p r and p ⊥ ), we also assume an equation of state for matter represented by (34). Integrating equations (37) (38) yields where the mass function m is given by (14), while the tangential pressure p ⊥ is obtained directly from (40) together with (42), that is At this point a remark is in order. First of all, we observe that the metric coefficient B is the same as the metric coefficient g rr appearing in the line element (20). This implies that the same analysis of the zeroes of g rr performed in Section III applies to the present case as well. Moreover, B appears in (43) in the denominator and this will cause the tangential pressure to become singular at the zeroes of B. On the other hand, the line element (36) can be cast into the form which is reminiscent of a dirty black hole metric provided it satisfies the conditions formulated in [74]. Under the assumption of an anisotropic energy-momentum tensor, one of them is that p ⊥ remains finite. The other one requires B to have two zeros. We will leave the examination if such a dirty black hole is viable at all to future projects and circumvent the latter condition by requiring µ < µ 0 . Note that such a condition also guarantees that the function φ is everywhere regular because it prevents the function B entering in (42) from having real roots. Hence, we conclude that the line element (44) describes a fuzzy self-gravitating dark matter droplet. In Figure 12 we plot the tangential pressure to show that it is indeed well-behaved for any value of r provided that µ < µ 0 . To accomplish that, it is convenient to introduce the variable y = (r/h) 1/ξ which allows to rewrite (43) as In Figure 13 we plot the radial and tangential pressures for the same choice of the Einasto parameter and the rescaled mass. Also in the case of a nonlocal equation of state, it turns out that the self-gravitating droplet does not exhibit any singularity at r = 0. To verify that, we first observe that the Kretschmann scalar associated to the line element (44) is given by Due to the presence of the term 1/r 2 in the above expression, it is not clear if the Kretschmann scalar is singularity free. However, if we expand the lower incomplete Gamma function entering in the metric coefficients by means of 6.5.29 in [52] and we let r → 0, we find that showing that there is no curvature singularity at r = 0. Regarding the asymptotic behaviour of the line element (44) we clearly have B → 1 at space-like infinity while by means of 6.5.3 and 6.5.32 in [52] it is possible to show that the function φ exhibits the asymptotic behaviour where exponentially decaying terms have been neglected. This shows that e φ(r) → 1 as r → ∞ and therefore, the manifold described by (44) is Minkowski flat asymptotically away. We conclude this The plot has been obtained by applying 6.5.29 in [52] to expand the lower incomplete Gamma function in (45) where the first hundred terms in the expansion have been considered. The same behaviour can be observed for all other values of ξ considered in Table II. FIG. 13: Plot of the radial and tangential pressures h 2 pr (solid line) and h 2 p ⊥ (dotted line) versus y = (r/h) 1/ξ in the case ξ = 1/2 and µ = 0.7. 6.5.29 in [52] has been used to expand the lower incomplete Gamma function in (45) where the first hundred terms in the expansion have been considered. The minima of the radial and tangential pressures denoted by yr and y ⊥ do not coincide and they are located at yr = 2.6038 and y ⊥ = 2.2335, respectively. section by showing that our droplet allows for bound states of massive particles. By means of (25. 16) in [67] we immediately find that the effective potential for the droplet is given by where ǫ and ℓ have been already defined in the previous section. Since U ef f is nonnegative, a matching procedure with the effective potential for a Schwarzschild BH cannot be achieved in this case. However, if we insist that the total mass M of the droplet coincides with the mass M BH of the black hole at the galactic centre and we introduce the rescalings r * = r/r s , L = ℓ/r s and H = h/r s with r s = 2M BH , we can rewrite φ and B as follows where we made use of 6.5.29 in [52]. From Fig. 14 we see that in the regime of low L the effective potential exhibits a minimum for which bound states for massive particles can form and a maximum corresponding to an unstable orbit. The situation is very different in the case of a massless particle where no bound states are allowed. We conclude with the observation that this model is not appropriate to reproduce the galactic motion of S-stars but it is nevertheless interesting because it indicates the possibility that a mass of dark matter modelled in terms of an anisotropic fluid with an Einasto energy density profile permits stable trajectories for massive particle over a large region as it can be seen in Fig. 14. Such a droplet albeit not suitable to model the central galactic object could possibly be formed in other parts of the galaxy. We conclude this section by estimating the value of the density and pressures at the centre of the droplet and comparing them with that for degenerate matter. More specifically, we consider the case of a neutron star with typical density 10 17 Kg/m 3 and degenerate pressure of the order 10 31 ÷ 10 34 Pa while in our model we take M = 10M ⊙ . By means of (8) and (34), we find that where in the expression for ρ(0) we used the fact that the function 1/(ξΓ(3ξ)) has a global maximum for ξ = 0.1538. We summarized the values of the density and radial pressure at the centre of the droplet in Table III.   TABLE III: Typical values of ρ(0) and pr(0) at the centre of the droplet. The scaling factors hE, hB1 and hB2 are chosen as in [11,39]. M = 10M⊙ ρ(0) Kg/m 3 pr(0) Pa hE = 2.121 · 10 −9 Kpc ≤ 2.17 · 10 9 ≤ 8.01 · 10 23 hB1 = 3.89 · 10 −11 Kpc ≤ 4.36 · 10 12 ≤ 1.30 · 10 29 hB2 = 1.43 · 10 −5 Kpc ≤ 8.85 · 10 −5 ≤ 2.64 · 10 12

V. CONCLUSIONS
Motivated by the DM dominance of the galaxy we established a connection between DM and the BH physics by considering the Einasto density profile. We showed that starting with this profile and coupling it to an anisotropic energy-momentum tensor, it is possible to derive different black hole solutions by making certain choices for the underlying EOS. In the case of an EOS of the form p r = −ρ, we show that it is possible to construct a self-gravitating droplet or a BH depending on the values of the mass parameter µ. If there is a horizon we could call such an object a fuzzy BH [54]. These objects made of DM are in nature different from the DM clumps constructed in [7] albeit the idea to connect the central galactic object to DM is similar. We also discovered that the Einasto inspired black holes increases its Hawking temperature, as the horizon radius shrinks, until the temperature reaches a maximum after which the temperature decreases sharply and vanishes exactly at the radius of the extremal black hole. In both cases, a fuzzy droplet or a fuzzy BH, it is possible to obtain an effective potential which governs the equation of motion in such way that the orbits will be as in the case of a standard galactic BH. If, instead of the previous EOS, we adopt a nonlocal EOS, it is possible to construct a self-gravitating droplet but it seems that a negative radial pressure cannot be avoided. Moreover, an analysis of the effective potential shows the occurrence of bound states as well as the presence of an unstable orbit in the case of massive particle and low values of their angular momenta.  We also notice that the Einasto profile is a generalization of a Gaussian. Black Holes based on such density profile have been considered before. As a result, issues regarding the perturbation [75] or exotic spacetime structure by considering a more general atlas [76] would proceed along similar lines.
Finally, while finishing this manuscript, we found a paper [77] which is relevant to our work because it shares conclusions similar to those obtained here but by a different method. More precisely, [77] showed by a numerical simulation that if we replace the central supermassive BH by an object made of darkinos, this mass of DM would not only produce the same kinematics for S-stars but can also explain the G2 anomaly [78].