Looking for ultralight dark matter near supermassive black holes
Abstract
Measurements of the dynamical environment of supermassive black holes (SMBHs) are becoming abundant and precise. We use such measurements to look for ultralight dark matter (ULDM), which is predicted to form dense cores (“solitons”) in the centre of galactic halos. We search for the gravitational imprint of an ULDM soliton on stellar orbits near Sgr A* and by combining stellar velocity measurements with Event Horizon Telescope imaging of M87*. Finding no positive evidence, we set limits on the soliton mass for different values of the ULDM particle mass . The constraints we derive exclude the solitons predicted by a naive extrapolation of the soliton–halo relation, found in DMonly numerical simulations, for (from Sgr A*) and (from M87*). However, we present theoretical arguments suggesting that an extrapolation of the soliton–halo relation may not be adequate: in some regions of the parameter space, the dynamical effect of the SMBH could cause this extrapolation to overpredict the soliton mass by orders of magnitude.
I Introduction
Supermassive black holes (SMBHs) reside in most galaxies Kormendy:1995er ; Ferrarese:2004qr ; Narayan:2005ie and their properties (mass, spin, and close environment) are under rapidly improving observational scrutiny. Two SMBHs for which very precise data exists are Sgr A* in the dynamical centre of the Milky Way (MW) and M87* in the elliptical galaxy M87. Near Sgr A*, the orbit of the star S2 in the S star cluster Genzel:2003cn ; Genzel2010 ; 2017ApJ…837…30G has been observed along more than a full lap 2009ApJ…707L.114G ; Abuter:2018drb . In the case of M87*, the event horizon telescope (EHT) collaboration has very recently released a breathtaking image of the BH shadow Akiyama:2019eap .
In this paper we show that measurements of the dynamical environment of SMBHs provide an interesting probe of ultralight dark matter (ULDM) Hu:2000ke ; Amendola:2005ad ; Svrcek:2006yi ; Arvanitaki:2009fg ; Marsh:2015xka . ULDM gained wide interest partially because in the window , it could alleviate smallscale puzzles facing the dark matter paradigm Hu:2000ke ; DelPopolo:2016emo ; Hui:2016ltb . This mass range, moreover, defines the absolute lower bound for the possible mass of dark matter. At the centre of galactic halos ULDM is expected to develop cored density profiles Arbey:2001qi ; Lesgourgues:2002hk ; Chavanis:2011zi ; Chavanis:2011zm ; Schive:2014dra ; Schive:2014hza ; Marsh:2015wka ; Calabrese:2016hmp ; Chen:2016unw ; Schwabe:2016rze ; Veltmaat:2016rxo ; Hui:2016ltb ; GonzalesMorales:2016mkl ; Robles:2012uy ; Bernal:2017oih ; Mocz:2017wlg ; Mukaida:2016hwd ; Vicens:2018kdk ; Bar:2018acw ; Eby:2018ufi ; BarOr:2018pxz ; Marsh:2018zyw ; Chavanis:2018pkx ; Emami:2018rxq ; Levkov:2018kau ; Broadhurst:2019fsl ; Hayashi:2019ynr ; Bar:2019bqz , commonly referred to as “solitons”^{1}^{1}1A more appropriate term is oscillatons; but we will stick to solitons in what follows., corresponding to quasistationary minimum energy solutions of the equations of motion. The ULDM soliton could be detected given detailed knowledge of the mass distribution in the inner halo. Such detailed view is provided by SMBH precision measurements: in the case of Sgr A*, any additional mass distribution between the periastron and apoastron of the S2 orbit ( pc) is constrained at the level of few percent 2017ApJ…837…30G . Measurements of stellar motions at larger distances Genzel2010 ( pc) provide more constraints. For M87*, a combination of the EHT measurement with analysis of stellar velocity dispersion at distances of kpc can be translated into the constraint . We will show that these observations probe ULDM at a meaningful level.
The problem of a minimallycoupled massive scalar field in the strong gravity regime around a BH (including the superradiance phenomenon Brito:2015oca ) was investigated in the literature Unruh:1976fm ; Detweiler1980 ; Arvanitaki:2009fg ; Arvanitaki:2010sy ; Cardoso2011 ; Ferreira:2017pth ; Boskovic:2018rub ; Cardoso:2018tly ; Hui2019 ; Benone2019 ; Chen2019 , recently also in the context of M87* Davoudiasl:2019nlo . Other works Hui:2016ltb ; Bar:2018acw ; BarOr:2018pxz considered the interplay between ULDM and black holes on galactic scales within the Newtonian approximation. Our approach focuses on the intermediate case, where on the one hand a Newtonian analysis is applicable but on the other hand, the SMBH dominates the dynamics.
Several other constraints on ULDM have appeared in the literature. The matter power spectrum revealed by Ly forest analyses is in tension with Armengaud:2017nkf ; Irsic:2017yje ; Zhang:2017chj ; Kobayashi:2017jcf ; Leong:2018opi (see also Bozek:2014uqa ; Hlozek:2017zzf ). Rotation curves of lowsurfacebrightness galaxies (LSBs) also disfavour Bar:2018acw ; Bar:2019bqz , if one accepts the soliton cores predicted by numerical simulations Schive:2014dra ; Schive:2014hza ; Veltmaat:2018dfz . Independent evidence from rotation curve data against ULDM cores was reported in Deng:2018jjz . Dynamical heating of the MW disk Church:2018sro and a preliminary analysis of stellar streams Amorisco:2018dcn disfavour . A weaker bound comes from pulsar timing measurements Porayko:2018sfa of scalar metric perturbations induced by ULDM Khmelnitsky:2013lxt , which exclude . Ref. Marsh:2018zyw showed that a dynamical analysis of a central star cluster in EridanusII could potentially probe ULDM up to .
The paper is outlined as follows.
Sec. II sets the stage for our investigation, introducing a few basic properties of the ULDM soliton, explaining how ballpark numbers for the soliton mass motivate us to look for ULDM near SMBHs, and highlighting a few of the complications we will encounter.
In Sec. III we use observations to search for solitons, considering first Sgr A* (Secs. III.1.1 and III.1.2) and then M87* (Sec. III.2.1). Our goal is to examine how different values of the soliton mass, , affect measurements of the SMBH dynamical environment for different assumed values of . The observational constraints that we found for Sgr A* and for M87* are compared to theoretical expectations in Sec. III.1.3 and Sec. III.2.2, respectively.
Theoretical benchmarks for the soliton are explained in Secs. IV and V, with some details postponed to App. A. The basic benchmark we look at, in Sec. IV, comes from the soliton–host halo relation found in the DMonly numerical simulations of Refs. Schive:2014dra ; Schive:2014hza . We find that the soliton–halo relation is tested by the SMBH data in a new range of compared to previous tests. However, the soliton–halo relation involves caveats that prevent us from turning the constraints on into robust exclusion on . First, the solitons we consider must account for the effect of a SMBH, whereas the numerical simulations included only ULDM. Second, the simulations were only run for a limited range of host halo masses and ULDM particle masses, while we explore more massive halos and more massive particles.
In Sec. V we consider the question of dynamical relaxation. Using the relaxation time estimate of Ref. Levkov:2018kau , combined with observations made in Bar:2018acw ; Bar:2019bqz , we show that dynamical relaxation may become a bottleneck for soliton formation for and . At , for example, the soliton mass prescribed by dynamical relaxation could be an order of magnitude lower than that predicted by naive extrapolation of the scaling relation of Schive:2014dra ; Schive:2014hza , even when one ignores the impact of a SMBH.
We summarise our results in Sec. VI.
We leave some details to appendices. In App. A we review the structure of the soliton in the regime where the dynamics is dominated by a SMBH, but where the Newtonian approximation is still valid. A simple approximation for the soliton profile is introduced to facilitate numerical calculations. We calculate the time scale characterising the absorption of a soliton into the SMBH. Our results suggest that over much of the parameter space of interest, Sgr A* and M87* could absorb ULDM too fast to allow for a soliton to be established.
Finally, in App. B we outline the parametric region where nongravitational selfinteractions, motivated by axionlike particle models of ULDM, could affect our results.
Ii Setting the Stage
The density profile of a selfgravitating soliton is cored with characteristic radius and mass related by (see, e.g. Bar:2018acw )^{2}^{2}2We define as the radius at which the soliton mass density decreases by a factor of 2 compared to its value at the origin Schive:2014dra ; Schive:2014hza .
(1) 
A glance at the properties of Sgr A*, with mass dominating the dynamics out to a few pc Genzel2010 , suggests that stellar orbit measurements with % accuracy in the SMBHdominated region could be sensitive to provided that . Interestingly, this ballpark for is consistent with a naive extrapolation (in ) of the results of DMonly numerical simulations Schive:2014dra ; Schive:2014hza .
These estimates look promising, and we will see that measurements of Sgr A* do lead to constraints that test the extrapolation of Schive:2014dra ; Schive:2014hza . However, the presence of the SMBH complicates the situation. First of all, when the SMBH dominates the dynamics, the soliton shape is distorted. The characteristic radius becomes independent of and, instead of Eq. (1), is given by (see App. A.1)
(2) 
The soliton mass is then an independent parameter. Whether the soliton can be probed by stellar orbits, or not, depends on via its relation to the soliton central density, .
The particular complication due to the transition from Eq. (1) to Eq. (2) does not turn out to be a show stopper, but it does illustrate the impact of the SMBH. We will consider a number of other complications, such as the possible impact of the SMBH on the naive large extrapolation (by large, we mean ) of the scaling of Schive:2014dra ; Schive:2014hza . A key caveat suggested by our findings is that in much of the parameter space where SMBH measurements could naively test ULDM, the soliton may actually be consumed by accretion into the SMBH. Understanding what really happens in this case requires simulating the coevolving SMBH+ULDM systems, which is beyond the scope of this paper. Our takehome message after considering these complications will be that while SMBH measurements open up an interesting avenue to search for ULDM, the theoretical uncertainties are still too large to allow for robust exclusion, at least based on the observables that we analysed.
Iii Looking for ULDM near SMBHs
In this section we derive observational constraints on ULDM solitons from stellar orbits near Sgr A* in the MW (Sec. III.1) and from the EHT measurements and stellar dispersion analyses of M87* (Sec. III.2). While the results can be read and understood without referring to the technical details of the soliton’s structure, the underlying calculations employ tools and results that are explained in subsequent sections. In particular, we use the BHdeformed soliton shape calculation of App. A.1, and compare our constraints to theoretical benchmarks which are explained in Sec. IV, Sec. V and App. A.
iii.1 Milky Way
We now discuss the constraints obtained from observations around the SMBH in the MW, specifically from the orbit of the star S2 (Sec. III.1.1) and from observations of a stellar disk (Sec. III.1.2)^{3}^{3}3Position and polarization measurements in near infrared, attributed to flares of Sgr A* gravity2018 , may also constrain the properties of the SMBH.. These constrain the pc and pc regions, respectively. Farther away, at the few pc region and outwards, the stellar mass contribution becomes comparable to the BH mass 2016ApJ…821…44F . This makes the analysis more involved, beyond the scope of the current work. In Sec. III.1.3 we compare the constraints from observations to theoretical expectations.
iii.1.1 The orbit of S2
Precision measurements of the orbits of stars in the S star cluster at the centre of the MW (see, e.g. Genzel:2003cn ; 2017ApJ…837…30G ) are sensitive to the mass distribution near the SMBH. The discriminatory power between an extended mass distribution to an isolated point mass (BH) arises from the eccentricity of the stellar orbit. In Fig. 1 we show a schematic view of the elliptic orbit of the B2type star S2, for which more than a full orbit has been recorded 2009ApJ…707L.114G . An extended mass distribution , defined as the mass within the green filled shell extending between the periastron and apastron of the orbit, can be constrained independently of the central mass inside of the periastron, shown by the internal white region with a black point representing the SMBH. It should be noted that in the presence of a significant extended mass distribution, the orbit of S2 would exhibit strong precession. In that limit, Fig. 1 should be thought of as showing only the osculating orbit of the star.
Our computation follows Ref. Lacroix:2018zmg . There, VLT measurements of the orbit of S2 up to 2016 were used to constrain the distribution of dark mass, which was assumed to exhibit a density spike towards the SMBH. Here, we use the same data to constrain an ULDM soliton.
We use an orbitfitting procedure as developed in Refs. Ghez2008 ; Gillessen:2008qv ; 2009ApJ…707L.114G ; Boehle2016 ; 2017ApJ…837…30G and described in detail in Ref. Lacroix:2018zmg . The procedure reconstructs the evolution of the position and velocity of the star on its orbit as a function of time, and constrains the properties of the gravitational potential by fitting the parameters of the model to the data, consistently combined in the likelihood. The data includes right ascension, declination, and radial velocity of S2 from VLT measurements 2017ApJ…837…30G .^{4}^{4}4Combining the VLT data with the data from the Keck observatory Boehle2016 , using the procedure of Ref. 2009ApJ…707L.114G , only improves the limits on by an factor at the price of having 4 additional parameters to reconcile the coordinate systems for both data sets, which leads to additional degeneracies and longer computing times. Therefore, in Fig. 2 we only account for the VLT data. The 14 parameters of the problem are the mass of the central object, , and its six phasespace coordinates, namely its distance , its position on the sky (right ascension , declination ), and velocity (, , ), as well as the six phasespace coordinates of the star, and the total soliton mass which characterizes the mass profile of the soliton. We assume that the SMBH and the soliton are concentric.
Given an assumed BH mass , ULDM particle mass and soliton mass , we compute using the formulae given in App. A.1. For the parameters of the star we consider the initial conditions on the polar radius and angle in the plane of the orbit, , , and their corresponding derivatives , , as well as the inclination angle with respect to the plane of the orbit and the standard longitude of the ascending node^{5}^{5}5Note that we cannot rely on the 6 standard orbital elements that characterize Keplerian orbits since in the presence of an extended mass the orbit is no longer Keplerian. We note however that and still define the plane of the orbit even for nonKeplerian motion..
Following Lacroix:2018zmg , we use PyMultiNest Buchner2014 , which relies on the MultiNest multimodal nested Monte Carlo sampling code Feroz2009 , to derive the posterior probability distribution of the parameters of the model, in particular . We fix the ULDM mass for a given Monte Carlo run, and scan over by means of independent runs.
We find that including an ULDM soliton does not modify the Bayesian evidence in a statistically significant way. Therefore, we derive upper limits on at the 95% confidence level^{6}^{6}6We use the corner.py Python module ForemanMackey2016 .. The excluded range in the (,) plane is shown by the redshaded region in Fig. 2, marked by “S2”. Total soliton masses down to are excluded at the 95% confidence level for . It should be noted that above , the entire soliton is confined within the pericentre of the orbit of S2, such that the total soliton mass is degenerate with the BH mass.
iii.1.2 A stellar disk
The orbit reconstruction of S2 and other wellmeasured members of the S star cluster can probe an ULDM soliton in the inner pc around Sgr A*. This translates into constraints that are particularly strong for . Additional measurements of stellar kinematics extending to pc provide somewhat less precise estimates of the central mass Genzel2010 , but the larger distance probed by these measurements makes them sensitive to smaller values of .
Ref. Beloborodov2006 analysed the kinematics of a clockwiserotating disk (CWD) of stars spanning distances of pc around Sgr A*. For an 8 kpc distance estimate to Sgr A*, they found . Combined with the results for the S2 orbit reconstruction, we deduce the constraint:
(3) 
This constraint, translated to the plane, is shown by the greenshaded region, marked by “CWD” in Fig. 2.
iii.1.3 Comparison of constraints with theoretical explanations
It is interesting to compare the observational constraints of Fig. 2 to theoretical expectations for the soliton mass. These are described as follows:

The thin solid blue line in Fig. 2 shows the value of predicted by Eq. (6) using for the MW halo 10.1093/mnras/stw2759 ; Watkins2019 . This benchmark ignores possible effects due to the SMBH.

The thin dashed blue line shows a more conservative prediction, obtained from Eq. (7) subject to the assumption that the SMBH formation preceded the soliton formation: the reasoning behind this prediction is summarised in Sec. IV. Note that if, on the other hand, the soliton formed early preceding the SMBH, then the SMBH formation may actually attract more ULDM mass into the soliton; in which case could exceed not only the thin dashed line, but also the solid blue line.

The thick solid black line shows the constraint on the soliton mass, that arises if one assumes that the soliton is dominantly formed in the kinetic regime via dynamical relaxation. This line is computed using Eq. (10) following the reasoning presented in Sec. V. It is equivalent to the central value of the blue band in Fig. 5. The SMBH is ignored in this computation.

In the blueshaded region, occupying approximately half of the plot^{7}^{7}7With the exception of the upperright corner, where due to the degeneracy between and , the latter is compatible with zero. at , a rough estimate of the time scale for absorption of the soliton by the SMBH, computed in Sec. A.2, is shorter than Gyr. In this region, the soliton may be entirely eaten by the BH and precise determination of the dynamics would require simulating the coevolving SMBH and ULDM systems.

Finally, in the region above the thickdotted magenta line, axionlike particle models of ULDM, where initial field misalignment in a cosine potential determines the ULDM relic abundance, predict that nongravitational selfinteractions could modify the soliton solution. Note that the SMBH causes the field to compress, making nonlinearities more important than in the selfgravitating soliton case. We also note, however, that the precise range in where selfinteractions become important depends strongly on the initial misalignment: a % tuning in the initial conditions would shift the dashed magenta line up by a factor of . The details are given in App. B.
iii.2 Messier 87
In this section we derive observational constraints on ULDM solitons by comparing the EHT measurement and stellar kinematics measurements of M87*, the SMBH in M87. In Sec. III.2.1 we present the constraint and in Sec. III.2.2 we discuss its implications.
iii.2.1 EHT BH shadow vs. stellar kinematics in M87
The EHT collaboration has recently reported an image of the shadow of the SMBH M87* Akiyama:2019eap .^{8}^{8}8For a previous discussion see Lu:2014zja and references therein. The BH shadow observed by the EHT translates into a gravitational angular radius Luminet:1979nyg ; Chandrasekhar:1985kt ; Falcke:1999pj ; Takahashi:2004xh
(4) 
where is the distance to the BH. In what follows we will assume a fiducial distance of Mpc, in which case one finds .
The EHT result allows for a new test of the mass distribution in the inner region of the galaxy.^{9}^{9}9Here we consider the constraints on the mass distribution far away from the BH horizon, in the weak field regime. The mass distribution in the strong field regime could, in principle, be tested too Lacroix:2012nz . Ref. Gebhardt:2011yw analysed stellar kinematics, where the most detailed and precise data used in the analysis fell in the range . Combined with the value of used by Gebhardt:2011yw , their result for translates into . Comparing Eq. (4) to the results of Gebhardt:2011yw , we see that an additional mass distribution within , parametrised by , is constrained by:
(5) 
For simplicity, we interpret Eq. (5) to hold for , understanding that additional information could be deduced in a more detailed analysis extending to somewhat larger or smaller .
Given the values of and , we can calculate the mass profile of an ULDM soliton of total mass and use Eq. (5) to derive constraints in the plane. The result of this exercise is shown in Fig. 3. In the greyshaded region, marked “M87 SMBH”, the soliton mass contribution exceeds the upper limit defined by Eq. (5) for . To compute the plot we use the fiducial values and Mpc; varying these fiducial values within the range and does not affect the results appreciably.
The uncertainty in Eq. (5) is dominantly systematic, and there is room to regard it with care. We comment that Ref. Walsh:2013uua considered gasdynamical models at smaller radii compared to those entering the stellar dispersion of Gebhardt:2011yw , and found the result : a factor of two lower than that found by Gebhardt:2011yw . The discrepancy between the gas models of Ref. Walsh:2013uua and the stellar kinematics analysis of Ref. Gebhardt:2011yw could be due to uncertainties in modelling the gas distribution, e.g. the inclination of the gas disc. Prior to the EHT measurement, one might have argued that the factor of two mismatch between the results of Refs. Gebhardt:2011yw and Walsh:2013uua could in principle come from a dark halo contribution. The EHT closes this window of opportunity with a result at that confirms the conclusions of Ref. Gebhardt:2011yw . In our analysis we therefore used the stellar dispersion analysis of Ref. Gebhardt:2011yw . To get a quick estimate of the impact of changing the allowed soliton mass from , shown in Eq. (5), to , one can multiply the horizontal lower boundary of the greyshaded area (equal to in Eq. (5)) by a factor of .
iii.2.2 Comparison of constraints with theoretical expectations
Here we briefly compare the observational constraint with theoretical expectations.
The thin solid blue line in Fig. 3 shows the value of predicted by Eq. (6) using for M87 Wu:2005wi . The thin dashed blue line shows the more conservative prediction that is obtained from Eq. (7), as explained in Sec. IV. In the region above the thickdotted magenta line, nongravitational selfinteractions may be important, as discussed in App. B. In the blueshaded region, at , the time scale for absorption of the soliton by the SMBH, estimated in App. A.2, is shorter than . The solid black line shows the constraint on the soliton mass, which arises if one assumes that the soliton is dominantly formed in the kinetic regime via dynamical relaxation (see Sec. V).
Iv How much mass in the soliton?
The DMonly numerical simulations of Ref. Schive:2014dra ; Schive:2014hza related the soliton mass, , to the mass of the galactic halo in which it occurs, :
(6) 
Ref. Bar:2018acw ; Bar:2019bqz showed that for DMonly halos, Eq. (6) is equivalent to the statement:
(7) 
Namely, the numerical simulations of Schive:2014dra ; Schive:2014hza are finding that the kinetic energy per particle is the same for ULDM particles in the largescale host halo and in the central soliton.
Independent simulations by Ref. Veltmaat:2018dfz showed a result consistent with Eq. (6). The simulations of Ref. Mocz:2017wlg found a different scaling, but as discussed in Bar:2018acw it remains to be seen if the initial conditions employed in that simulation biased the result. To our knowledge, a direct comparison between the soliton growth found in the simulations of Ref. Levkov:2018kau and those of Schive:2014dra ; Schive:2014hza had not yet been made. There is, therefore, room for caution in accepting the soliton–halo relation. We expect that the theoretical situation will become clearer in the near future as different groups test the validity of Eq. (6).
We have used Eq. (6) as an illustrative benchmark with which observational constraints in the plane can be compared. However, it is important to note that the reference to Eqs. (67) entails two significant caveats:
In the rest of this section we discuss the effect of the SMBH on ULDM energetics in the soliton region, with possible implications on the soliton formation. Some technical details are postponed to App. A. The question of dynamical relaxation is considered in Sec. V.
While Eq. (7)—just like Eq. (6)—was not tested in simulations that include an external baryonic contribution to the gravitational potential, it is suggestive to consider it as evidence for kinetic equilibration between the ULDM reservoir in the halo and in the soliton. This kinetic equilibration is unlikely to represent true steadystate equilibrium, but it could correspond to a bottleneck in the soliton formation which slows down once the soliton grows to saturate Eq. (7). Assuming that this is the case, we could use Eq. (7) to estimate the outcome of dynamical heating of the inner region of the halo due to the SMBH^{10}^{10}10Ref. Bar:2019bqz considered the related effect of stellar and gas mass components in low surfacebrightness galaxies., which would affect the value of on the LHS of Eq. (7).
In the limit of SMBH dominance, the soliton specific kinetic energy becomes where (see App. A.1, in particular Eq. (16)). Therefore, if , then Eq. (7) cannot be satisfied. This would imply that the soliton formation is halted by the dynamical heating due to the SMBH.
In Figs. 2 and 3 the soliton mass predicted by Eq. (7), including the effect of the SMBH on the LHS, was shown by the thin dashed blue lines.
The discussion above is relevant if the SMBH preceded the formation of the soliton. In reality, the timing could be reversed: the soliton may form early and precede the SMBH (and even, perhaps, help to seed the SMBH). If this latter ordering is the relevant one, then we see no reason to expect that SMBH formation would quench the soliton (that is, apart from the possibility that the SMBH could eatup the soliton altogether: this is considered in App. A.2). Adiabatic contraction could actually increase the soliton mass above the prediction of Eq. (6). An example illustrating related dynamics was considered in Ref. 2017arXiv171201947C , which simulated the scenario where a population of massive point particles (“stars”) was added to a halo containing an existing soliton, and the system was then allowed to evolve dynamically. The stars in the simulation flowed to the centre of the halo, causing the soliton to absorb additional mass from the largescale ULDM host halo. After a Hubble time, the system attained near steadystate containing a soliton that was more massive than the prediction of Eq. (6).
V The question of the relaxation time
In this section we consider the question of dynamical relaxation, which, for massive halos with , may imply a bottleneck for the soliton mass for . Throughout this section we ignore the dynamical impact of a SMBH on the soliton: tackling the combined problem of soliton formation via dynamical relaxation alongside a simultaneouslyforming SMBH is beyond the scope of this paper. Our analysis could therefore be justified if the soliton forms before the SMBH. Beyond the particular focus of the current work, our analysis in this section may be more generally useful for the understanding of DMonly numerical simulations.
Ref. Hui:2016ltb pointed out that interference patterns in ULDM facilitate dynamical relaxation by acting as quasiparticles of mass , where is the mass density and is the characteristic velocity dispersion (the precise matching to quasiparticles was derived in Ref. BarOr:2018pxz ). The time scale for dynamical relaxation due to twobody gravitational interactions between the quasiparticles is , where is the impact parameter for gravitational collisions and is the Coulomb logarithm. This scaling was verified in Ref. Levkov:2018kau using numerical simulations, which found that the relaxation time in a region of spatial extent is given by
(8) 
with . Within time , a soliton forms in an ensemble of ULDM particles which initially contained no soliton.
Consider an initial NFW halo density profile Navarro:1996gj characterised by the radius parameter , a concentration parameter and a halo virial mass (defined as the mass within ),
(9) 
with . If we set , we can calculate as a function of the mass coordinate . Inverting this relation gives the mass contained in the relaxed region as a function of the time available for relaxation.
In Fig. 4 we show the mass contained in the relaxed region of an NFW halo, given (blue) or (red) to relax. The width of each shaded band shows for different values of the halo concentration parameter and scale parameter , varying in the range and kpc BoylanKolchin:2009an ; Piffl:2014mfa . For comparison, the dashed black line shows the soliton mass expected from extrapolating the numerical simulations of Schive:2014dra ; Schive:2014hza ; Veltmaat:2018dfz to .
The mass could be considered as an upper bound on the mass of the soliton forming in a halo in the kinetic regime Hui:2016ltb . However, is probably a weak upper bound. If we use Eq. (7) to estimate the soliton by equating
(10) 
where is the velocity dispersion at , we obtain a stronger (and potentially more realistic) upper bound. This comes from noting that, for a selfgravitating soliton, Eq. (10) translates into Bar:2018acw . This upper bound is shown in Fig. 5.
It is interesting to consider the parametric scaling of the soliton mass upper bound derived from Eq. (10). For massive halos like MW and M87, with , and for or so, the mass contained in the relaxed region within a Hubble time is much smaller than . The density, enclosed mass and velocity dispersion in the relevant region can then be approximated by , and . Inserting this into Eq. (8) we have^{11}^{11}11We chose the presentation of Eq. (11) to highlight the fact that the halo parameters enter only via the combination .:
(11) 
In the regime of validity of Eq. (11), the mass contained in the relaxed region at fixed scales as , up to a logarithmic correction. Eq. (10) then leads to , up to a logarithmic correction. Noting that , and ignoring the dependence on the concentration parameter , we have . The scaling at is apparent in Fig. 5 when comparing the shaded band (coming from Eq. (10)) to the dashed line (coming from the naive extrapolation of Eq. (6)) which scales as .
For , the relaxed region extends up to and the small expansion of Eq. (11) does not apply. Instead, becomes approximately independent of ( vanishes at ). Because is approximately independent of , it follows that the soliton determined by Eq. (10) is approximately independent of the time available for relaxation. In this regime the upper bound prescribed by Eq. (10) approximately coincides with Eq. (7).
Finally, note that Eq. (8) and the corresponding dynamical relaxation bottleneck apply to the kinetic regime . These considerations may significantly overestimate the actual time scale for soliton formation, if the initial conditions admit small velocity dispersion, as would be the case in cosmological halos that decouple from the Hubble flow before virialisation Levkov:2018kau .
Vi Summary
Measurements of the mass and dynamical environment of SMBH are becoming increasingly precise. Two SMBHs where precision measurements have become available are Sgr A* (, via stellar orbits) and M87* (, via BH shadow imaging and stellar velocity dispersion). The SMBH measurements provide clean probes of the mass distribution in the region where the BH dominates the dynamics.
We study the implications of the SMBH measurements for ULDM. Analytical arguments and numerical simulations predict that ULDM should form dense cores (“solitons”) in the centre of galactic halos. We present a search for the gravitational imprint of an ULDM soliton with mass on the orbit of the star S2. We also consider constraints from the observations of a stellar disk. We find no evidence for and use the data to derive constraints in the () plain, where is the ULDM particle mass (see Fig. 2). We then use stellar velocity dispersion analyses combined with the Event Horizon Telescope (EHT) measurement of M87* to constrain the presence of a soliton, which could manifest itself as excess mass in the velocity dispersion data as compared to the EHT determination of (see Fig. 3).
DMonly numerical simulations suggest a scaling relation Schive:2014dra ; Schive:2014hza that predicts given the host halo mass and the particle mass . The observational constraints we find from SMBH data exclude a naive extrapolation of this relation for (from Sgr A*) and (from M87*).
However, a number of theoretical arguments lead us to expect that the naive extrapolation of the soliton–halo relation is not adequate.
The most significant caveats are the process of soliton absorption by the SMBH, and the question of dynamical relaxation. Both caveats become more pronounced as is increased, and suggest that a naive extrapolation of the scaling relation of Schive:2014dra ; Schive:2014hza to could overpredict the soliton mass by orders of magnitude. Thus, while the SMBH analysis is a potentially interesting discovery tool for ULDM, the theoretical uncertainties are too large to make it a robust exclusion tool.
Note added: while this manuscript was prepared for publication, Ref. Desjacques2019 appeared, discussing stellar dynamics near SMBHs in the presence of ULDM and claiming the strong constraint . This conclusion appears to assume that the soliton–host halo relation of Schive:2014dra ; Schive:2014hza remains valid up to . However, we note that as explained here, the naive extrapolation of the soliton–host halo relation up to cannot, at present, be used to infer robust constraints on ULDM.
Acknowledgements.
We thank Asimina Arvanitaki, Joshua Eby, Reinhard Genzel, Rainer Schödel and Wei Xue for useful discussions, and Vitor Cardoso and Sergey Sibiryakov for comments on the manuscript. We wish we could consult with Tal Alexander about S star orbits. KB is incumbent of the Dewey David Stone and Harry Levine career development chair. TL has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie SkłodowskaCurie grant agreement No. 713366. The work of TL has also been supported by the Spanish Agencia Estatal de Investigación through the grants PGC2018095161BI00, IFT Centro de Excelencia Severo Ochoa SEV20160597, and Red Consolider MultiDark FPA201790566REDC.Appendix A ULDM soliton near a black hole
In this section we consider the ULDM soliton in the presence of a SMBH. In Sec. A.1 we consider the basic features of the soliton solution, reviewing and extending an earlier analysis in Ref. Bar:2018acw (see App. B there) and introducing some useful approximations. In Sec. A.2 we consider the accretion of mass from the soliton into the SMBH.
a.1 Soliton shape
We consider a real, massive, free scalar field satisfying the KleinGordon equation of motion (EoM) and minimally coupled to gravity.^{12}^{12}12Analyses of interacting fields can be found in, e.g. Refs. Chavanis:2011zi ; Chavanis:2011zm ; RindlerDaller:2012vj ; Desjacques:2017fmf . We look for a sphericallysymmetric, quasistationary bound state solution in the nonrelativistic regime. To this end we decompose as
(12) 
where is the gravitational Newton’s constant and is an eigenvalue of the problem. We do not need the explicit value of in what follows.
We rescale the spatial coordinate,
(13) 
assuming spherical symmetry. The EoM for and for the Newtonian gravitational potential , sourced by , is
(14)  
(15) 
where
(16) 
The lowest energy bound state solutions of Eqs. (1415) are parametrised by a single continuous positive parameter that can be chosen as the value of the field at the origin: . Once we solve for and , all other solutions are obtained via
(17)  
(18) 
The mass density of the soliton is
(19) 
The total mass is given by
(20) 
where is the mass in . In the limit , one finds
(21) 
In the opposite limit one finds
(22) 
Solving numerically for is not particularly difficult. However, it is useful to find a semianalytic approximation for the mass density that applies in both limits and . One such useful approximation is given by
(23)  
(24) 
The functions and are found by a numerical fit to the exact solution, and are shown in Fig. 6. At we have and , which reproduces the Coulomb solution, , where the soliton profile is purely controlled by the BH gravity. At we have , independent of , and . This gives a soliton profile that is independent of and reconstructs with good accuracy the selfgravitating exact numerical solution of .
Using Eqs. (2324), the total mass of the soliton reads
(25) 
Using this expression we can obtain a numerical fit of the parameter , expressed as a function of the physical mass ratio . The results of this fit are shown in Fig. 7. For later use we also compute the soliton mass enclosed in the shell between two radii and (corresponding to physical radii ):
Finally, let us check the domain of validity of the Newtonian approximation of the EoM, which we have been using so far. The characteristic momentum associated with the ULDM field in the soliton is given by
(27) 
When the BH dominates the solution we find
(28) 
Similarly, the kinetic energy per unit mass is when the BH dominates the solution. Therefore, in the regime where the BH dominates the dynamics, the characteristic particle velocity associated with the soliton is given by and the Newtonian approximation remains adequate as long as . Referring back to Eq. (16) we learn that for a SMBH with (appropriate for Sgr A*) the Newtonian approximation breaks down for , while for (appropriate for M87*) the approximation breaks down already for . Once we come close to the relativistic region, , phenomenae such as superradiance come into play (see, e.g. Arvanitaki:2010sy ). Since our analysis did not account for these effects, we stop our exploration of the parameter space at and in Figs. 2 and 3, respectively.
a.2 Absorption by the black hole
A SMBH would accrete mass from the soliton and grow at its expense Hui:2016ltb . Ref. Bar:2018acw estimated the absorption rate of the soliton by the SMBH, in the limit that the soliton’s selfgravity dominates over that due to the SMBH over most of the soliton region. Here we extend this estimate to the case where the SMBH gravity dominates throughout the soliton.
The cross section for absorption of a scalar particle with mass and nonrelativistic momentum by a Schwarzschild BH, whose size is much smaller than the Compton wavelength of the particle, was calculated in Ref. Unruh:1976fm :
(29) 
where the parameter was defined in Eq. (16) and . The mass accretion rate flowing from the soliton into the BH is then
(30) 
where is the soliton central density. Using our parametrisation, and combining Eqs. (23,25,27) with Eq. (30), we can write the characteristic time for the SMBH to absorb the soliton
We are interested in the case where the BH dominates the soliton dynamics (the opposite limit was studied in Ref. Bar:2018acw ), wherein is given by Eq. (28), leading to . Thus, we can neglect in the denominator of Eq. (30). In the same limit, the soliton central density is related to the total soliton mass by . Combining these expressions with Eq. (30), we find the characteristic time for the SMBH to absorb the soliton:
This result is independent of the soliton mass.
In the second line of Eq. (A.2) the BH mass is scaled to that of Sgr A*, and in the third line it is scaled to that of M87*. In each line is scaled such that the soliton absorption time is about 10 Gyr. The strong dependence of on implies that for each of the two systems (Sgr A* and M87*), larger than the values shown in Eq. (A.2) would result in the soliton being rapidly consumed by the SMBH.
Throughout this section we have ignored corrections due to BH rotation. However, the results developed in Benone2019 suggest that this may not lead to a large correction: for a Kerr BH with spin parameter , using Eq. (51) in Benone2019 we find
(33) 
Thus, the absorption time scale could be longer by up to a factor of for a maximally rotating BH.
We conclude that much of the parameter space considered in the body of the paper may be strongly affected by BH absorption. This motivates a more careful analysis, that (i) goes beyond the single wavelength approximation of the soliton, adopted here for simplicity using the characteristic wavelength from Eq. (27), and that (ii) solves the BH and the soliton coevolution.
Appendix B Nongravitational selfinteractions
ULDM could be realised by axionlike particles (see, e.g. Svrcek:2006yi ; for an interesting alternative, see Davoudiasl:2017jke ), with the relic abundance set by initial misalignment of the field in a cosine potential. In this case, nongravitational selfinteractions are expected to become significant in some regions of the plane.
Assuming the potential
(34) 
one finds the following selfinteraction term near the minimum of the potential, with . In the misalignment mechanism, the ULDM relic abundance is given by where is the initial value of the field during inflation. Thus, the observed DM abundance is obtained for GeV. Let us define , where apriori one expects , and fix as needed for the DM abundance. Plugging this into the expression for gives
(35) 
Selfinteractions can be estimated to contribute a fraction of the energy density of the soliton solution when
(36) 
where this estimate holds for . Noting that the ULDM mass density is given by , and utilising formulas from App. A.1, one finds that selfinteractions contribute a fraction to the soliton energy density when
(37) 
In the BHdominated case, this reduces to
(38) 
In Figs. 2 and 3 we used a dashed magenta line to depict Eq. (37), setting and . From the equation above it is clear that smaller values of (for example ) corresponding to mild finetuning of the initial conditions, would push the onset of selfinteractions two orders of magnitude up in , making this effect irrelevant for the parameter space analyzed in this paper. Therefore, while it is interesting to contemplate what changes could occur in the soliton properties due to selfinteractions via a generic cosine potential, we leave this investigation to other works.
References
 (1) J. Kormendy and D. Richstone, “Inward bound: The Search for supermassive black holes in galactic nuclei,” Ann. Rev. Astron. Astrophys. 33 (1995) 581.
 (2) L. Ferrarese and H. Ford, “Supermassive black holes in galactic nuclei: Past, present and future research,” Space Sci. Rev. 116 (2005) 523–624, arXiv:astroph/0411247 [astroph].
 (3) R. Narayan, “Black holes in astrophysics,” New J. Phys. 7 (2005) 199, arXiv:grqc/0506078 [grqc].
 (4) R. Genzel et al., “The Stellar cusp around the supermassive black hole in the Galactic Center,” Astrophys. J. 594 (2003) 812–832, arXiv:astroph/0305423 [astroph].
 (5) R. Genzel, F. Eisenhauer, and S. Gillessen, “The Galactic Center Massive Black Hole and Nuclear Star Cluster,” Rev. Mod. Phys. 82 (2010) 3121–3195, arXiv:1006.0064 [astroph.GA].
 (6) S. Gillessen, P. M. Plewa, F. Eisenhauer, R. Sari, I. Waisberg, M. Habibi, O. Pfuhl, E. George, J. Dexter, S. von Fellenberg, T. Ott, and R. Genzel, “An Update on Monitoring Stellar Orbits in the Galactic Center,” ApJ 837 no. 1, (Mar, 2017) 30, arXiv:1611.09144 [astroph.GA].
 (7) S. Gillessen, F. Eisenhauer, T. K. Fritz, H. Bartko, K. DoddsEden, O. Pfuhl, T. Ott, and R. Genzel, “The Orbit of the Star S2 Around SGR A* from Very Large Telescope and Keck Data,” ApJ 707 no. 2, (Dec, 2009) L114–L117, arXiv:0910.3069 [astroph.GA].
 (8) GRAVITY Collaboration, R. Abuter et al., “Detection of the gravitational redshift in the orbit of the star S2 near the Galactic centre massive black hole,” Astron. Astrophys. 615 (2018) L15, arXiv:1807.09409 [astroph.GA].
 (9) Event Horizon Telescope Collaboration, K. Akiyama et al., “First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole,” Astrophys. J. 875 no. 1, (2019) L6.
 (10) W. Hu, R. Barkana, and A. Gruzinov, “Cold and fuzzy dark matter,” Phys. Rev. Lett. 85 (2000) 1158–1161, arXiv:astroph/0003365 [astroph].
 (11) L. Amendola and R. Barbieri, “Dark matter from an ultralight pseudoGoldsoneboson,” Phys. Lett. B642 (2006) 192–196, arXiv:hepph/0509257 [hepph].
 (12) P. Svrcek and E. Witten, “Axions In String Theory,” JHEP 06 (2006) 051, arXiv:hepth/0605206 [hepth].
 (13) A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper, and J. MarchRussell, “String Axiverse,” Phys. Rev. D81 (2010) 123530, arXiv:0905.4720 [hepth].
 (14) D. J. E. Marsh, “Axion Cosmology,” Phys. Rept. 643 (2016) 1–79, arXiv:1510.07633 [astroph.CO].
 (15) A. Del Popolo and M. Le Delliou, “Small scale problems of the CDM model: a short review,” Galaxies 5 no. 1, (2017) 17, arXiv:1606.07790 [astroph.CO].
 (16) L. Hui, J. P. Ostriker, S. Tremaine, and E. Witten, “Ultralight scalars as cosmological dark matter,” Phys. Rev. D95 no. 4, (2017) 043541, arXiv:1610.08297 [astroph.CO].
 (17) A. Arbey, J. Lesgourgues, and P. Salati, “Quintessential haloes around galaxies,” Phys. Rev. D64 (2001) 123528, arXiv:astroph/0105564 [astroph].
 (18) J. Lesgourgues, A. Arbey, and P. Salati, “A light scalar field at the origin of galaxy rotation curves,” New Astron. Rev. 46 (2002) 791–799.
 (19) P.H. Chavanis, “Massradius relation of Newtonian selfgravitating BoseEinstein condensates with shortrange interactions: I. Analytical results,” Phys. Rev. D84 (2011) 043531, arXiv:1103.2050 [astroph.CO].
 (20) P. H. Chavanis and L. Delfini, “Massradius relation of Newtonian selfgravitating BoseEinstein condensates with shortrange interactions: II. Numerical results,” Phys. Rev. D84 (2011) 043532, arXiv:1103.2054 [astroph.CO].
 (21) H.Y. Schive, T. Chiueh, and T. Broadhurst, “Cosmic Structure as the Quantum Interference of a Coherent Dark Wave,” Nature Phys. 10 (2014) 496–499, arXiv:1406.6586 [astroph.GA].
 (22) H.Y. Schive, M.H. Liao, T.P. Woo, S.K. Wong, T. Chiueh, T. Broadhurst, and W. Y. P. Hwang, “Understanding the CoreHalo Relation of Quantum Wave Dark Matter from 3D Simulations,” Phys. Rev. Lett. 113 no. 26, (2014) 261302, arXiv:1407.7762 [astroph.GA].
 (23) D. J. E. Marsh and A.R. Pop, “Axion dark matter, solitons and the cuspcore problem,” Mon. Not. Roy. Astron. Soc. 451 no. 3, (2015) 2479–2492, arXiv:1502.03456 [astroph.CO].
 (24) E. Calabrese and D. N. Spergel, “UltraLight Dark Matter in UltraFaint Dwarf Galaxies,” Mon. Not. Roy. Astron. Soc. 460 no. 4, (2016) 4397–4402, arXiv:1603.07321 [astroph.CO].
 (25) S.R. Chen, H.Y. Schive, and T. Chiueh, “Jeans Analysis for Dwarf Spheroidal Galaxies in Wave Dark Matter,” Mon. Not. Roy. Astron. Soc. 468 no. 2, (2017) 1338–1348, arXiv:1606.09030 [astroph.GA].
 (26) B. Schwabe, J. C. Niemeyer, and J. F. Engels, “Simulations of solitonic core mergers in ultralight axion dark matter cosmologies,” Phys. Rev. D94 no. 4, (2016) 043513, arXiv:1606.05151 [astroph.CO].
 (27) J. Veltmaat and J. C. Niemeyer, “Cosmological particleincell simulations with ultralight axion dark matter,” Phys. Rev. D94 no. 12, (2016) 123523, arXiv:1608.00802 [astroph.CO].
 (28) A. X. GonzálesMorales, D. J. E. Marsh, J. Peñarrubia, and L. Ureña López, “Unbiased constraints on ultralight axion mass from dwarf spheroidal galaxies,” Mon. Not. Roy. Astron. Soc. 472 (2017) 1346, arXiv:1609.05856 [astroph.CO].
 (29) V. H. Robles and T. Matos, “Flat Central Density Profile and Constant DM Surface Density in Galaxies from Scalar Field Dark Matter,” Mon. Not. Roy. Astron. Soc. 422 (2012) 282–289, arXiv:1201.3032 [astroph.CO].
 (30) T. Bernal, L. M. FernándezHernández, T. Matos, and M. A. Rodrï¿œguezMeza, “Rotation Curves of HighResolution LSB and SPARC Galaxies in Wave (Fuzzy) and Multistate (Ultralight Boson) Scalar Field Dark Matter,” arXiv:1701.00912 [astroph.GA].
 (31) P. Mocz, M. Vogelsberger, V. H. Robles, J. Zavala, M. BoylanKolchin, A. Fialkov, and L. Hernquist, “Galaxy formation with BECDM: I. Turbulence and relaxation of idealized haloes,” Mon. Not. Roy. Astron. Soc. 471 no. 4, (2017) 4559–4570, arXiv:1705.05845 [astroph.CO].
 (32) K. Mukaida, M. Takimoto, and M. Yamada, “On Longevity of Iball/Oscillon,” JHEP 03 (2017) 122, arXiv:1612.07750 [hepph].
 (33) J. Vicens, J. Salvado, and J. MiraldaEscudé, “Bosonic dark matter halos: excited states and relaxation in the potential of the ground state,” arXiv:1802.10513 [astroph.CO].
 (34) N. Bar, D. Blas, K. Blum, and S. Sibiryakov, “Galactic rotation curves versus ultralight dark matter: Implications of the solitonhost halo relation,” Phys. Rev. D98 no. 8, (2018) 083027, arXiv:1805.00122 [astroph.CO].
 (35) J. Eby, K. Mukaida, M. Takimoto, L. C. R. Wijewardhana, and M. Yamada, “Classical Nonrelativistic Effective Field Theory and the Role of Gravitational Interactions,” arXiv:1807.09795 [hepph].
 (36) B. BarOr, J.B. Fouvry, and S. Tremaine, “Relaxation in a Fuzzy Dark Matter Halo,” Astrophys. J. 871 (2019) 28, arXiv:1809.07673 [astroph.GA].
 (37) D. J. E. Marsh and J. C. Niemeyer, “Strong Constraints on Fuzzy Dark Matter from Ultrafaint Dwarf Galaxy Eridanus II,” arXiv:1810.08543 [astroph.CO].
 (38) P.H. Chavanis, “A predictive model of BEC dark matter halos with a solitonic core and an isothermal atmosphere,” arXiv:1810.08948 [grqc].
 (39) R. Emami, T. Broadhurst, G. Smoot, T. Chiueh, and L. H. Nhan, “A Soliton Solution for the Central Dark Masses in Globular Clusters and Implications for the Axiverse,” arXiv:1806.04518 [astroph.CO].
 (40) D. G. Levkov, A. G. Panin, and I. I. Tkachev, “Gravitational BoseEinstein condensation in the kinetic regime,” Phys. Rev. Lett. 121 no. 15, (2018) 151301, arXiv:1804.05857 [astroph.CO].
 (41) T. Broadhurst, I. de Martino, H. N. Luu, G. F. Smoot, and S. H. H. Tye, “Ghostly Galaxies as Solitons of BoseEinstein Dark Matter,” arXiv:1902.10488 [astroph.CO].
 (42) K. Hayashi and I. Obata, “Nonsphericity of ultralight axion dark matter halos in the Galactic dwarf spheroidal galaxies,” arXiv:1902.03054 [astroph.CO].
 (43) N. Bar, K. Blum, J. Eby, and R. Sato, “Ultralight dark matter in disk galaxies,” arXiv:1903.03402 [astroph.CO].
 (44) R. Brito, V. Cardoso, and P. Pani, “Superradiance,” Lect. Notes Phys. 906 (2015) pp.1–237, arXiv:1501.06570 [grqc].
 (45) W. G. Unruh, “Absorption CrossSection of Small Black Holes,” Phys. Rev. D14 (1976) 3251–3259.
 (46) S. L. Detweiler, “KLEINGORDON EQUATION AND ROTATING BLACK HOLES,” Phys. Rev. D22 (1980) 2323–2326.
 (47) A. Arvanitaki and S. Dubovsky, “Exploring the String Axiverse with Precision Black Hole Physics,” Phys. Rev. D83 (2011) 044026, arXiv:1004.3558 [hepth].
 (48) V. Cardoso, S. Chakrabarti, P. Pani, E. Berti, and L. Gualtieri, “Floating and sinking: The Imprint of massive scalars around rotating black holes,” Phys. Rev. Lett. 107 (2011) 241101, arXiv:1109.6021 [grqc].
 (49) M. C. Ferreira, C. F. B. Macedo, and V. Cardoso, “Orbital fingerprints of ultralight scalar fields around black holes,” Phys. Rev. D96 no. 8, (2017) 083017, arXiv:1710.00830 [grqc].
 (50) M. Boškovi?, F. Duque, M. C. Ferreira, F. S. Miguel, and V. Cardoso, “Motion in timeperiodic backgrounds with applications to ultralight dark matter haloes at galactic centers,” Phys. Rev. D98 (2018) 024037, arXiv:1806.07331 [grqc].
 (51) V. Cardoso, . J. C. Dias, G. S. Hartnett, M. Middleton, P. Pani, and J. E. Santos, “Constraining the mass of dark photons and axionlike particles through blackhole superradiance,” JCAP 1803 no. 03, (2018) 043, arXiv:1801.01420 [grqc].
 (52) L. Hui, D. Kabat, X. Li, L. Santoni, and S. S. C. Wong, “Black Hole Hair from Scalar Dark Matter,” arXiv:1904.12803 [grqc].
 (53) C. L. Benone and L. C. B. Crispino, “Massive and charged scalar field in KerrNewman spacetime: Absorption and superradiance,” Phys. Rev. D99 no. 4, (2019) 044009, arXiv:1901.05592 [grqc].
 (54) Y. Chen, J. Shu, X. Xue, Q. Yuan, and Y. Zhao, “Probing Axions with Event Horizon Telescope Polarimetric Measurements,” arXiv:1905.02213 [hepph].
 (55) H. Davoudiasl and P. B. Denton, “Ultra Light Boson Dark Matter and Event Horizon Telescope Observations of M87*,” arXiv:1904.09242 [astroph.CO].
 (56) E. Armengaud, N. PalanqueDelabrouille, D. J. E. Marsh, J. Baur, and C. Yï¿œche, “Constraining the mass of light bosonic dark matter using SDSS Lyman forest,” Mon. Not. Roy. Astron. Soc. 471 no. 4, (2017) 4606–4614, arXiv:1703.09126 [astroph.CO].
 (57) V. Iršič, M. Viel, M. G. Haehnelt, J. S. Bolton, and G. D. Becker, “First constraints on fuzzy dark matter from Lyman forest data and hydrodynamical simulations,” Phys. Rev. Lett. 119 no. 3, (2017) 031302, arXiv:1703.04683 [astroph.CO].
 (58) J. Zhang, J.L. Kuo, H. Liu, Y.L. S. Tsai, K. Cheung, and M.C. Chu, “Is Fuzzy Dark Matter in tension with Lymanalpha forest?,” arXiv:1708.04389 [astroph.CO].
 (59) T. Kobayashi, R. Murgia, A. De Simone, V. Iršič, and M. Viel, “Lymanalpha Constraints on Ultralight Scalar Dark Matter: Implications for the Early and Late Universe,” arXiv:1708.00015 [astroph.CO].
 (60) K.H. Leong, H.Y. Schive, U.H. Zhang, and T. Chiueh, “Testing extremeaxion wave dark matter using the BOSS LymanAlpha forest data,” Mon. Not. Roy. Astron. Soc. 484 (2019) 4273, arXiv:1810.05930 [astroph.CO].
 (61) B. Bozek, D. J. E. Marsh, J. Silk, and R. F. G. Wyse, “Galaxy UVluminosity function and reionization constraints on axion dark matter,” Mon. Not. Roy. Astron. Soc. 450 no. 1, (2015) 209–222, arXiv:1409.3544 [astroph.CO].
 (62) R. Hlozek, D. J. E. Marsh, and D. Grin, “Using the Full Power of the Cosmic Microwave Background to Probe Axion Dark Matter,” Mon. Not. Roy. Astron. Soc. 476 no. 3, (2018) 3063–3085, arXiv:1708.05681 [astroph.CO].
 (63) J. Veltmaat, J. C. Niemeyer, and B. Schwabe, “Formation and structure of ultralight bosonic dark matter halos,” Phys. Rev. D98 no. 4, (2018) 043509, arXiv:1804.09647 [astroph.CO].
 (64) H. Deng, M. P. Hertzberg, M. H. Namjoo, and A. Masoumi, “Can Light Dark Matter Solve the CoreCusp Problem?,” arXiv:1804.05921 [astroph.CO].
 (65) B. V. Church, J. P. Ostriker, and P. Mocz, “Heating of Milky Way disc Stars by Dark Matter Fluctuations in Cold Dark Matter and Fuzzy Dark Matter Paradigms,” arXiv:1809.04744 [astroph.GA].
 (66) N. C. Amorisco and A. Loeb, “First constraints on Fuzzy Dark Matter from the dynamics of stellar streams in the Milky Way,” arXiv:1808.00464 [astroph.GA].
 (67) N. K. Porayko et al., “Parkes Pulsar Timing Array constraints on ultralight scalarfield dark matter,” Phys. Rev. D98 no. 10, (2018) 102002, arXiv:1810.03227 [astroph.CO].
 (68) A. Khmelnitsky and V. Rubakov, “Pulsar timing signal from ultralight scalar dark matter,” JCAP 1402 (2014) 019, arXiv:1309.5888 [astroph.CO].
 (69) GRAVITY Collaboration, R. Abuter et al., “Detection of orbital motions near the last stable circular orbit of the massive black hole SgrA*,” AAP 618 (2018) L10. https://doi.org/10.1051/00046361/201834294.
 (70) T. K. Fritz, S. Chatzopoulos, O. Gerhard, S. Gillessen, R. Genzel, O. Pfuhl, S. Tacchella, F. Eisenhauer, and T. Ott, “The Nuclear Cluster of the Milky Way: Total Mass and Luminosity,” ApJ 821 no. 1, (Apr, 2016) 44, arXiv:1406.7568 [astroph.GA].
 (71) T. Lacroix, “Dynamical constraints on a dark matter spike at the Galactic Centre from stellar orbits,” Astron. Astrophys. 619 (2018) A46, arXiv:1801.01308 [astroph.GA].
 (72) A. M. Ghez et al., “Measuring Distance and Properties of the Milky Way’s Central Supermassive Black Hole with Stellar Orbits,” Astrophys. J. 689 (2008) 1044–1062, arXiv:0808.2870 [astroph].
 (73) S. Gillessen, F. Eisenhauer, S. Trippe, T. Alexander, R. Genzel, F. Martins, and T. Ott, “Monitoring stellar orbits around the Massive Black Hole in the Galactic Center,” Astrophys. J. 692 (2009) 1075–1109, arXiv:0810.4674 [astroph].
 (74) A. Boehle, A. M. Ghez, R. Schödel, L. Meyer, S. Yelda, S. Albers, G. D. Martinez, E. E. Becklin, T. Do, J. R. Lu, K. Matthews, M. R. Morris, B. Sitarski, and G. Witzel, “An Improved Distance and Mass Estimate for Sgr A* from a Multistar Orbit Analysis,” ApJ 830 (Oct., 2016) 17, arXiv:1607.05726.
 (75) J. Buchner, A. Georgakakis, K. Nandra, L. Hsu, C. Rangel, M. Brightman, A. Merloni, M. Salvato, J. Donley, and D. Kocevski, “Xray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue,” Astron. Astroph. 564 (Apr, 2014) A125, arXiv:1402.0004 [astroph.HE].
 (76) F. Feroz, M. P. Hobson, and M. Bridges, “MULTINEST: an efficient and robust Bayesian inference tool for cosmology and particle physics,” MNRAS 398 (Oct., 2009) 1601–1614, arXiv:0809.3437.
 (77) D. ForemanMackey, “corner.py: Scatterplot matrices in Python,” The Journal of Open Source Software 1 (June, 2016) .
 (78) A. M. Beloborodov, Y. Levin, F. Eisenhauer, R. Genzel, T. Paumard, S. Gillessen, and T. Ott, “Clockwise stellar disk and the dark mass in the galactic center,” ApJ 648 (Sept., 2006) 405–410, astroph/0601273. http://adsabs.harvard.edu/abs/2006ApJ...648..405B.
 (79) P. J. McMillan, “The mass distribution and gravitational potential of the Milky Way,” Monthly Notices of the Royal Astronomical Society 465 no. 1, (10, 2016) 76–94. https://doi.org/10.1093/mnras/stw2759.
 (80) L. L. Watkins, R. P. van der Marel, S. T. Sohn, and N. W. Evans, “Evidence for an Intermediatemass Milky Way from Gaia DR2 Halo Globular Cluster Motions,” ApJ 873 (Mar., 2019) 118, arXiv:1804.11348.
 (81) R.S. Lu, A. E. Broderick, F. Baron, J. D. Monnier, V. L. Fish, S. S. Doeleman, and V. Pankratius, “Imaging the Supermassive Black Hole Shadow and Jet Base of M87 with the Event Horizon Telescope,” Astrophys. J. 788 (2014) 120, arXiv:1404.7095 [astroph.IM].
 (82) J. P. Luminet, “Image of a spherical black hole with thin accretion disk,” Astron. Astrophys. 75 (1979) 228–235.
 (83) S. Chandrasekhar, “The mathematical theory of black holes,” in Oxford, UK: Clarendon (1992) 646 p., OXFORD, UK: CLARENDON (1985) 646 P. 1985.
 (84) H. Falcke, F. Melia, and E. Agol, “Viewing the shadow of the black hole at the galactic center,” Astrophys. J. 528 (2000) L13, arXiv:astroph/9912263 [astroph].
 (85) R. Takahashi, “Shapes and positions of black hole shadows in accretion disks and spin parameters of black holes,” J. Korean Phys. Soc. 45 (2004) S1808–S1812, arXiv:astroph/0405099 [astroph]. [Astrophys. J.611,996(2004)].
 (86) T. Lacroix and J. Silk, “Constraining the distribution of dark matter at the Galactic Centre using the highresolution Event Horizon Telescope,” Astron. Astrophys. 554 (2013) A36, arXiv:1211.4861 [astroph.GA].
 (87) K. Gebhardt, J. Adams, D. Richstone, T. R. Lauer, S. M. Faber, K. Gultekin, J. Murphy, and S. Tremaine, “The BlackHole Mass in M87 from Gemini/NIFS Adaptive Optics Observations,” Astrophys. J. 729 (2011) 119, arXiv:1101.1954 [astroph.CO].
 (88) J. L. Walsh, A. J. Barth, L. C. Ho, and M. Sarzi, “The M87 Black Hole Mass from Gasdynamical Models of Space Telescope Imaging Spectrograph Observations,” Astrophys. J. 770 (2013) 86, arXiv:1304.7273 [astroph.CO].
 (89) X.A. Wu and S. Tremaine, “Deriving the mass distribution of m87 from globular clusters,” Astrophys. J. 643 (2006) 210–221, arXiv:astroph/0508463 [astroph].
 (90) J. H. H. Chan, H.Y. Schive, T.P. Woo, and T. Chiueh, “How do stars affect DM halos?,” ArXiv eprints (Dec., 2017) , arXiv:1712.01947.
 (91) J. F. Navarro, C. S. Frenk, and S. D. M. White, “A Universal density profile from hierarchical clustering,” Astrophys. J. 490 (1997) 493–508, arXiv:astroph/9611107 [astroph].
 (92) M. BoylanKolchin, V. Springel, S. D. M. White, and A. Jenkins, “There’s no place like home? Statistics of Milky Waymass dark matter halos,” Mon. Not. Roy. Astron. Soc. 406 (2010) 896, arXiv:0911.4484 [astroph.CO].
 (93) T. Piffl et al., “Constraining the Galaxy’s dark halo with RAVE stars,” Mon. Not. Roy. Astron. Soc. 445 no. 3, (2014) 3133–3151, arXiv:1406.4130 [astroph.GA].
 (94) V. Desjacques and A. Nusser, “Axion core  halo mass and the black hole–halo mass relation: constraints on a few parsec scales,” arXiv:1905.03450 [astroph.CO].
 (95) T. RindlerDaller, T. RindlerDaller, P. R. Shapiro, and P. R. Shapiro, “Finding New Signature Effects on Galactic Dynamics to Constrain BoseEinsteinCondensed Cold Dark Matter,” in 6th International Meeting on Gravitation and Cosmology Guadalajara, Jalisco, Mexico, May 2125, 2012, pp. 163–182. 2012. arXiv:1209.1835 [astroph.CO]. https://inspirehep.net/record/1184893/files/arXiv:1209.1835.pdf.
 (96) V. Desjacques, A. Kehagias, and A. Riotto, “Impact of ultralight axion selfinteractions on the large scale structure of the Universe,” Phys. Rev. D97 no. 2, (2018) 023529, arXiv:1709.07946 [astroph.CO].
 (97) H. Davoudiasl and C. W. Murphy, “Fuzzy Dark Matter from Infrared Confining Dynamics,” Phys. Rev. Lett. 118 no. 14, (2017) 141801, arXiv:1701.01136 [hepph].