Dynamical constraints on a dark matter spike at the Galactic centre from stellar orbits
Laboratoire Univers & Particules de Montpellier (LUPM), CNRS & Université de Montpellier (UMR5299), Place Eugène Bataillon, 34095 Montpellier Cedex 05, France
email: thomas.lacroix@umontpellier.fr
Received:
16
January
2018
Accepted:
17
August
2018
In this work I use astrometric and spectroscopic data on the S2 star at the Galactic centre (GC) up to 2016 to derive specific constraints on the size of a dark matter (DM) spike around the central supermassive black hole Sgr A*. These limits are the best direct constraints on a DM spike at the GC for nonannihilating dark matter and exclude a spike with radius greater than a few tens of parsecs for cuspy outer halos and a few hundred parsecs for cored outer halos.
Key words: dark matter / Galaxy: center / celestial mechanics / black hole physics
© ESO 2018
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Dark matter profiles in the central regions of galaxies are poorly constrained at present and are the objects of intense debate. While observations seem to favour flat (cored) profiles, numerical simulations favour steeper profiles (cusps), leading to the cusp/core controversy (e.g. de Blok 2010 for a review). At subparsec scales, the dark matter (DM) distribution is even less constrained and can be significantly affected by the central supermassive black hole (SMBH). In particular, if the SMBH grows adiabatically, i.e. on a much longer timescale than the dynamical timescale, the DM density is expected to be significantly enhanced (by up to 10 orders of magnitude at the very centre) in a region corresponding to the sphere of influence of the black hole (BH), typically at parsec scales for the Milky Way. This leads to a very sharp morphological feature referred to as a DM spike, corresponding to a DM profile going as r^{−γsp}, with γ_{sp} typically between 2.25 and 2.5, depending on the slope of the initial DM halo (Gondolo & Silk 1999). DM spikes are of particular interest in the context of indirect DM searches since they lead to very strong signatures of DM annihilation and allow us to probe weakly annihilating DM particles (Gondolo & Silk 1999; Regis & Ullio 2008; Lacroix et al. 2014, 2015, 2017; Fields et al. 2014; Shapiro & Shelton 2016).
There is, however, considerable uncertainty on the formation and survival of DM spikes. In particular, the assumption of adiabaticity may not be verified in general. For instance, dynamical processes such as mergers can lead to weaker cusps (Merritt et al. 2002). However, binary scouring only occurs above parsec scales, while we are interested in the DM profile much closer in when studying the orbits of S stars, as discussed in the following. Moreover, the Milky Way (MW) is unlikely to have suffered such mergers in its recent past, as evidenced by the quiet history of the thick disk since the only major merger which occurred about 12 Gyr ago and is likely to have led to the formation of the bulge and the SMBH (Wyse 2001). A weaker cusp is also formed if the BH does not grow exactly at the centre of the DM halo (within ∼50 pc; Nakano & Makino 1999; Ullio et al. 2001) or if the BH growth cannot be considered adiabatic (Ullio et al. 2001), but the actual impact of these effects on the MW is unclear. Moreover, dynamical heating in the central stellar core would also soften a spike (Gnedin & Primack 2004). Another concern is that the nonobservation of a stellar spike (recent results point to a softer stellar cusp than previously thought, with slope ∼1.15; Schödel et al. 2018) would rule out the existence of a DM spike. However, if the BH grows (for example by gas accretion) mostly before the nuclear star cluster forms in the spike region, then the DM and stellar profiles are decoupled. Additionally, the nuclear star cluster in the most accepted view is formed by merging globular clusters, as in Antonini et al. (2015). This leads to different profiles for the DM and stellar distributions. Therefore, stars and DM essentially decouple, and the absence of a stellar spike in observations does not preclude the existence of a DM spike. On the other hand, additional dynamical processes can have the opposite effect of regenerating a spike, for example enhanced accretion of DM to counteract the depopulation of chaotic orbits in triaxial halos (Merritt & Poon 2004) or gravothermal collapse for selfinteracting DM (Ostriker 2000).
As a result the unclear status of the inner DM profile of galaxies as discussed above calls for direct probes. In particular, there is still no definitive evidence either in favour of or against such a high concentration of DM either in the MW or in any other galaxy. This is due in particular to the small size of the regions involved. Probing such regions requires high angular resolution and astrometric precision to characterize the gravitational potential. However, the inner region of the MW offers a unique window on the DM distribution at the Galactic centre (GC), thanks to the monitoring of the orbits of the S stars within ∼1 arcsec of the central BH. In particular, since it is the closest star to the BH observed so far, the S2 star has been extensively studied through monitoring campaigns based on observations conducted with the Very Large Telescope (VLT; Schödel et al. 2002; Gillessen et al. 2009b, a, 2017; GRAVITY Collaboration 2018) and the Keck observatory (Ghez et al. 2005, 2008; Boehle et al. 2016)^{1}. These series of observations have led to the reconstruction of the orbit of the star over roughly one and a half periods. In addition to tight constraints on the mass of the central SMBH, M_{BH}, and its distance from Earth, R_{0}, these two groups have shown that only a small fraction (typically 1–2%) of the mass of the SMBH can be in the form of an extended distribution (Ghez et al. 2008; Gillessen et al. 2009b, a, 2017; Boehle et al. 2016). Other constraints have been obtained on an extended component by studying the corresponding reconstructed mass profile (Hall & Gondolo 2006) or the pericentre shift of S2 (Zakharov et al. 2007; Iorio 2013).
Here I go a step further and I use astrometric and spectroscopic measurements of the orbit of S2 up to 2016 to set specific constraints on the DM distribution in the inner Galaxy. I present the first direct dynamical constraints from stellar orbits on the size of a DM spike, inside a DM halo constrained by larger scale kinematic data at kpc scales, for example from maser observations. This is especially interesting for nonannihilating or very weakly annihilating DM which is not expected to have significant observational signatures other than gravitational.
In Sect. 2 I describe the model along with the orbitfitting procedure, before presenting my results in Sect. 3. Finally, I conclude in Sect. 4.
2. Model and orbitfitting procedure
2.1. Calibration: the pointmass case
I rely on textbook results of standard mechanics in a central potential (e.g. Bate et al. 1971). I first recall the parameters of the problem in the BHonly case, which has an analytic solution, before moving on to the more general case of an extended mass distribution. The BHonly case serves as calibration for the orbitfitting procedure.
The orbitfitting procedure consists in reconstructing the time evolution of the position and velocity of the star on its orbit to determine the properties of the gravitational potential by fitting the parameters of the model to the data. In the case of one star orbiting a central point mass, the 13 parameters of the problem are the mass of the central object, here denoted M_{BH}, and its six phasespace coordinates, namely its distance R_{0}, its position on the sky (α_{BH}, δ_{BH}), and velocity (v_{α, BH}, v_{δ, BH}, v_{r, BH}), as well as the six phasespace coordinates of the star. However, the orbit of the star is more readily characterized analytically in terms of the six standard orbital elements: the semimajor axis a of the orbit, the eccentricity e, the time of pericentre passage t_{P}, and three angles, namely the inclination I of the orbital plane with respect to the plane of the sky, the longitude of the ascending node Ω, and the angle ω between the directions of the ascending node and the pericentre.
Although the motion of Sgr A* with respect to the local standard of rest (LSR), defined as the circular velocity at the radius of the Sun, is expected to be very small (Reid & Brunthaler 2004; Plewa et al. 2015), its position (α_{BH}, δ_{BH}) on the plane of the sky at a reference time t_{ref} and its velocity (v_{α, BH}, v_{δ, BH}, v_{r, BH}) relative to the LSR are unknown a priori and can be constrained through the orbitfitting procedure. In practice, the motion of the BH is accounted for through a linear term in the angular position of the star as a function of time. The reference time is taken to be 2009 yr (Gillessen et al. 2017) for the data set up to 2016, and 2005.4 yr for the data set up to 2009 (Gillessen et al. 2009a).
In this work, I used the data from the NTT/VLT and Keck observatories compiled in Boehle et al. (2016), Gillessen et al. (2017). In Gillessen et al. (2009a) the authors presented a robust method to consistently combine the two independent data sets for which the astrometric data feature a clear offset due to slight differences in the definition of the coordinate systems. More specifically, to account for the discrepancy between the two data sets, they introduced an offset in angular position (Δα, Δδ) and velocity (Δv_{α}, Δv_{δ}) on the plane of the sky to shift the Keck data back onto the VLT data. This was done by fitting the model with these 4 parameters in addition to the 13 parameters described before. In practice, this is achieved by shifting the observed right ascensions and declinations of the star measured with the Keck observatory by the quantities Δα + Δv_{α}(t − t_{ref}) and Δδ + Δv_{δ}(t − t_{ref}), respectively. I repeated this procedure here for the combined data set.
Throughout this work, I derived the posterior probability density function of model parameters using PyMultiNest (Buchner et al. 2014), which relies on the MultiNest code (Feroz et al. 2009) based on the multimodal nested sampling Monte Carlo technique (Feroz & Hobson 2008). Multimodal nested sampling is particularly suitable for studying highdimensional parameter spaces with possible degeneracies between parameters. The likelihood combines the data on right ascension, declination, and radial velocity of S2. I used uniform priors for all parameters except the position and velocity of the BH, for which I took Gaussian priors based on the results from Plewa et al. (2015): (α_{BH},δ_{BH}) = (0,0) ± (0.2,0.2) mas at t_{ref} = 2009 yr and (υ_{α,BH},υ_{δ,BH}) = (0,0) ± (0.1,0.1) mas yr^{−1}. I recovered the bestfit parameters and errors from Gillessen et al. (2017), as illustrated by the marginalized posterior distributions, for the BHonly case and the full VLT data set up to 2016 (see Fig. A.1). I also recovered the bestfit model for the combined VLT + Keck data set, using the prescriptions of Gillessen et al. (2009a) for the priors on Δα, Δδ, Δv_{α}, and Δv_{δ}. This served as a consistency check of the analysis chain, which was then applied to the study of the impact of a DM spike on the orbit of S2.
2.2. Extended mass
For the extended DM mass distribution, I consider two scenarios: the general case of a nonannihilating cold dark matter (CDM) candidate, and the more specific case of selfannihilating DM, applicable to candidates like weakly interacting massive particles (WIMPs). For nonannihilating DM the spike goes way inside the orbit of S2, down to the close vicinity of the SMBH (Sadeghian et al. 2013):
where R_{sp} is the radial extension of the spike, and the halo profile is assumed to be given by a generalized Navarro–Frenk–White (NFW) profile characterized by a slope index γ,
where r_{s} is the scale radius, and the scale density ρ_{s} is related to the local density ρ_{⊙} via
More specifically, the idea is to consider the various DM halos corresponding to the dynamically constrained Milky Way mass models from the analysis of McMillan (2017), and determine the maximum size of a DM spike inside that halo that does not cause a significant departure from the bestfitting BHonly orbit. The associated values of the local density, scale radius, and R_{0} from the analysis of McMillan (2017) are summarized in Table B.1 ^{2}.
For selfannihilating DM, the inner region of the DM spike is depleted since the DM density is so high that DM particles annihilate more efficiently. This results in a plateau of density ρ_{sat} = m_{DM}/(⟨σv⟩t_{BH}), where m_{DM} is the mass of the DM candidate, ⟨σv⟩ is the velocityaveraged annihilation cross section, and t_{BH} is the age of the central SMBH, which I take conservatively to be ∼10^{10} yr. The saturation plateau extends to a radius
The cases of nonannihilating and selfannihilating DM are both illustrated in Fig. 1. Shown are the density profiles for regular NFWlike halos and halos with a spike in the central region for nonannihilating DM and a selfannihilating 1 TeV DM candidate with three values of ⟨σv⟩ (see figure for details). The corresponding mass profiles are shown in the right panel of Fig. 1 with the same line styles. The profiles are illustrated with γ = 1, which gives a spike slope γ_{sp} = 7/3, and a spike radius R_{sp} ∼ 100 pc which corresponds to the 99.7% upper limit I obtain from deviations of the BHonly orbit, as discussed in Sect. 3. For the selfannihilating case, the values of the cross section are chosen to illustrate the point at which the annihilation plateau becomes as big as the characteristic size of the orbit, given by the semimajor axis constrained to be of the order of 5 mpc by the orbitfitting procedure. For m_{DM} ∼ 1 TeV and ⟨σv⟩ ∼ 10^{−26} cm^{3} s^{−1}, the mass enclosed inside the orbit is significantly reduced with respect to the case of nonannihilating or very weakly annihilating DM, down to values much smaller than a few percent of the BH mass, making deviations from the BHonly orbit undetectable^{3}.
Fig. 1. Left panel: density profiles for a generalized NFW halo (dashed lines) with γ = 1 and the same halo with a spike in the central region (solid lines) for a nonannihilating DM candidate (purple) and a selfannihilating 1 TeV DM candidate with ⟨σv⟩ ∼ 10^{−30} cm^{3} s^{−1}, ⟨σv⟩ ∼ 10^{−27} cm^{3} s^{−1}, and ⟨σv⟩ ∼ 10^{−26} cm^{3} s^{−1} (orange, blue, and red, respectively). Right panel: corresponding mass profiles, with the same line styles as in the left panel. The profiles are shown for illustration purposes for a spike radius R_{sp} ∼ 100 pc corresponding to the 99.7% upper limit from deviations of the BHonly orbit using the VLT data (see Sect. 3). The horizontal solid and dotted lines represent the combined 2009 and 2017 constraints, respectively. The vertical dotdashed line marks the characteristic size of the orbit of S2. 

Open with DEXTER 
In the absence of a spike, the DM halo has a negligible impact on the orbit of S2, as illustrated in the right panel of Fig. 1, due to a much smaller mass enclosed in the orbit. Moreover, for completeness I have also considered the effect of a realistic stellar profile ρ_{star} ∝ r^{−γstar}, with γ_{star} ∼ 1.15 and ρ_{star}(1 pc)=1.5 × 10^{5} M_{⊙} pc^{−3} (Schödel et al. 2018). However, the corresponding mass enclosed in the orbit is about three orders of magnitude below the critical mass needed to have an impact on the orbit. The same conclusion applies to the stellar bulge profile from McMillan (2017).
In the general case of an extended mass distribution around the central point mass, the orbit model is no longer analytic and one must rely on numerical tools to solve the equations of motion. First the polar radius r(t) is determined with Newton’s second law in the Galilean frame of the BH,
where L ≡ r^{2}θ̇ is the angular momentum modulus, G is the gravitational constant, Φ_{ext} is the potential created by the extended mass, and Φ_{S} accounts for the effect of Schwarzschild precession induced by the BH^{4}. The initial conditions r_{0} ≡ r(t_{0}) and ṙ_{0} ≡ ṙ(t_{0}) need to be specified, where t_{0} is chosen as the first epoch in the data, namely t_{0} = 1992.224 yr. For given values of r_{0} and ṙ_{0}, I use the odeint Python routine to solve for r(t). Once r(t) is known, θ(t) is obtained via
where θ_{0} ≡ θ(t_{0}) and , with θ̇_{0} ≡ θ̇(t_{0}). Orbital elements no longer characterize the orbit but only an osculating orbit in the general case of an extended mass, so they cannot be used to parametrize the problem. The free parameters for the star are now the initial conditions r_{0}, ṙ_{0}, θ̇_{0}, as well as I and Ω, which still characterize the plane of the orbit of the star. The parameters of the BH do not change^{5}.
3. Results
In this section, I present constraints on the size of a DM spike as a function of the slope of the DM halo obtained with the multimodal nested sampling analysis implemented in PyMultinest. It should be noted that the nested sampling procedure does find a nonzero bestfit value for the spike radius R_{sp} for all values of γ. However, the very mild increase in Bayesian evidence when adding a DM spike, which remains smaller than ΔlnZ ≈ 3, is insufficient to claim any preference for the BH + spike model (Kass & Raftery 1995). As a result, the data are consistent with the BHonly model. Nevertheless, it is still possible to exclude large values of the spike radius that would lead to a large DM mass inside the orbit, and thus to large deviations of the orbit of S2.
The resulting 95% and 99.7% confidence contours in the γR_{sp} plane are shown in Fig. 2 for the case of nonannihilating DM, for the dynamically constrained halo profiles from McMillan (2017), and using the combined VLT + Keck data sets up to 2009 (left panel), the VLTonly 2016 data set up to 2016 (middle panel), and the combined VLT + Keck data up to 2016 (right panel). The contours are computed for the values of halo slope γ for which McMillan (2017) derived constraints from various data sets including maser observations, namely γ = 0, 0.25, 0.5, 0.75, 1, 1.25, 1.5. I interpolated the results to obtain a smooth limit^{6}. Using the combined VLT+Keck data set up to 2016, I exclude at the 99.7% confidence level a DM spike with a spatial extension larger than 90 pc for an outer halo with γ = 1, and larger than 6 pc for an outer halo with γ = 1.5. For the combined 2009 and the VLTonly 2016 data sets, the limits are about a factor of 2 weaker. The analysis of McMillan (2017) seems to favour cuspy halos (γ ∼ 1), whereas a recent study on the dynamics of the Galactic bar favours a DM halo with slope γ < 0.6 (Portail et al. 2017). For such cored halos, our constraints are weaker, with a maximum spike radius of a few hundred pc. Nevertheless, these limits are the first direct constraints on a DM spike at the GC, valid for nonannihilating DM. This is especially interesting because it is applicable to any CDM candidate with no significant annihilation cross section.
Fig. 2. Median of the marginalized posterior distribution of the spike radius R_{sp} (black dashed), and associated 95% (red shaded) and 99.7% (light red shaded) confidence contours as a function of the slope γ of the corresponding dynamically constrained outer halos, using the combined VLT + Keck data up to 2009 (left panel), VLTonly data up to 2016 (middle panel), and VLT + Keck data up to 2016 (right panel). The upper right white regions correspond to excluded values of the spike radius. The lower left white regions are outside the 99.7% bestfit contours, but cannot robustly be interpreted in terms of exclusion of the associated values of the spike radius since the Bayesian evidence does not favour the presence of a spike over the BHonly model. The contours are derived for values of the slope of the outer halo dynamically constrained by McMillan (2017), and linearly interpolated to get smooth curves. The prediction from Gondolo & Silk (1999) is also shown as a benchmark model (cyan solid). 

Open with DEXTER 
It should be noted that the limit already improves by about a factor of 2 when going from the combined 2009 data set to the combined 2016 one. Significant additional improvements can thus be expected on the constraints with the more recent S2 data.
Constraints on an extended mass component can also be expressed more generally in terms of the total extended mass M_{ext} inside the characteristic size of the orbit. At a 99.7% confidence level, for the combined data up to 2009, I find M_{ext} ≲ 2 × 10^{5} M_{⊙}, corresponding to ∼4%M_{BH}, while for the complete data set up to 2016, I obtain M_{ext} ≲ 4 − 5 × 10^{4} M_{⊙}, corresponding to ∼1%M_{BH}. These results are consistent with the upper limits on a general extended mass component from Gillessen et al. (2009a, 2017). These limits are illustrated in the right panel of Fig. 1 by the horizontal dotted and solid black lines.
For selfannihilating DM, the constraints derived here are valid for ⟨σv⟩ < 10^{−30} cm^{3} s^{−1} for a benchmark particle mass m_{DM} = 1 TeV. As shown in the right panel of Fig. 1, for the same candidate mass, the mass enclosed in a sphere of radius of the characteristic size of the orbit of S2–typically 5 mpc–is decreased by about a factor of 2 for ⟨σv⟩ ∼ 10^{−27} cm^{3} s^{−1} and by about a factor of 20 for ⟨σv⟩ ∼ 10^{−26} cm^{3} s^{−1} with respect to the mass in the absence of annihilation. As a result, the upper limits on R_{sp} are weakened by about a factor of 10 for ⟨σv⟩ ∼ 10^{−27} cm^{3} s^{−1}, while no constraints can be set for ⟨σv⟩ ∼ 10^{−26} cm^{3} s^{−1}. More generally, when the radius R_{sat} of the saturation plateau due to selfannihilations (see Eq. (3)) reaches the size of the orbit of S2, the DM mass enclosed inside the orbit becomes too small to induce significant deviations from the BHonly orbit.
For illustration, the prescription from Gondolo & Silk (1999) is also shown in Fig. 2 (cyan solid line). This corresponds to a spike radius defined by
with α_{γ} ≈ 0.293γ^{4/9}, which only differs from ∼0.1 for γ ≪ 1. It should be noted that these predictions are indicative and can be significantly affected by the various dynamical processes discussed in Sect. 1. For very cuspy halos (γ ≈ 1.5), the combined 2016 data already exclude the prediction from Gondolo & Silk (1999) at a 95% confidence level.
4. Conclusion and outlook
In this work, I have used an orbitfitting procedure similar to those developed in Gillessen et al. (2017), Boehle et al. (2016) to derive specific constraints on the size of a DM spike for given outer DM halos dynamically constrained by larger scale observations. These limits are the best direct constraints on a DM spike at the GC for nonannihilating DM and exclude a spike with radius greater than a few tens of pc for cuspy outer halos and a few hundred pc for cored outer halos.
The addition of the 2017–2018 data from the VLT, which has monitored the pericentre passage of S2 (GRAVITY Collaboration 2018), will make these constraints significantly more stringent, especially thanks to the impressive capabilities of the imaging NACO instrument, the SINFONI spectrometer, and the exquisite astrometric precision of the GRAVITY instrument. However, I postpone the study of the subsequent constraints on the DM profile at the GC to a future work since additional subtleties related to the relativistic effects that have been detected using the new data (GRAVITY Collaboration 2018) may appear, and this warrants a dedicated study. The problem is also complicated further by having to model two pericentre passages when accounting for the entire data set since 1992, which increases the computingtime.
Additional improvements on the data could lead to even stronger constraints on the very inner DM profile at the GC. Firstly, in principle, S stars located further out than S2 would be more suited to probe the extended DM distribution for which the mass increases with radius. However, this comes at the price of longer periods, so that unlike S2, no additional stars have been monitored for about 1.5 periods. As a result, our constraints do not improve when including other stars further out such as S1 or S13 for which no significant precession is detectable yet. However, the situation will change when complete orbits are recorded for these stars. In addition, even more accurate astrometric and spectroscopic data will be instrumental to further improve upon these constraints. In particular, a 30 m extremely large telescope (ELT) should be able to probe an extended mass component as low as a few 10^{3} M_{⊙}, i.e. about one order of magnitude better than the current sensitivity, after only 10 years of observation (Weinberg et al. 2005). Moreover, an ELT would be able to break the degeneracy between relativistic effects and precession from an extended mass component. This would translate into sensitivity to DM spikes as small as a few pc even for cored outer halos, and even smaller for steeper halos.
The bestfit values of R_{0} from McMillan (2017) are consistent with the values obtained with the orbitfitting procedure. I do not keep R_{0} in the generalized NFW profile as a free parameter, but use it to fix the normalization of the halo profile in a way that is consistent with McMillan (2017). R_{0} is only kept free in the position and velocity of the star.
These considerations can be extended to other values of m_{DM} via Eq. (3).
Although the current data on S2 are not yet sensitive to relativistic effects (Gillessen et al. 2017), I include this postNewtonian precession effect for completeness since it is partly degenerate with the precession caused by the extended mass and as such can mildly affect the limits set on the DM profile. Other relativistic effects such as gravitational redshift essentially affect radial velocities, which are much less tightly constrained by observations than the position of the star, as discussed in Gillessen et al. (2009b).
Accounting for the uncertainty on the local DM density ρ_{⊙} from McMillan (2017) only leads to a 4% variation in the limits on R_{sp}.
Acknowledgments
I thank Benjamin Jaillant, Julien Lavalle, Vivian Poulin, and Martin Stref for fruitful discussions on this topic, and the anonymous referee for the very useful comments and suggestions. My work is supported by CNRSIN2P3. I also acknowledge support from the European Union’s Horizon 2020 research and innovation program under the Marie SkłodowskaCurie grant agreements No 690575 and No 674896, in addition to recurrent institutional funding by CNRSIN2P3 and the University of Montpellier.
References
 Antonini, F., Barausse, E., & Silk, J. 2015, ApJ, 812, 72 [NASA ADS] [CrossRef] [Google Scholar]
 Bate, R. R., Mueller, D. D., & White, J. E. 1971, Fundamentals of Astrodynamics (New (York, USA: Dover Publications)) [Google Scholar]
 Boehle, A., Ghez, A. M., Schödel, R., et al. 2016, ApJ, 830, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Buchner, J., Georgakakis, A., Nandra, K., et al. 2014, A&A, 564, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 de Blok, W. J. G. 2010, Adv. Astron., 2010, 789293 [NASA ADS] [CrossRef] [Google Scholar]
 Feroz, F., & Hobson, M. P. 2008, MNRAS, 384, 449 [NASA ADS] [CrossRef] [Google Scholar]
 Feroz, F., Hobson, M. P., & Bridges, M. 2009, MNRAS, 398, 1601 [NASA ADS] [CrossRef] [Google Scholar]
 Fields, B. D., Shapiro, S. L., & Shelton, J. 2014, Phys. Rev. Lett., 113, 151302 [NASA ADS] [CrossRef] [Google Scholar]
 ForemanMackey, D. 2016, J. Open Source Softw., 1 [NASA ADS] [CrossRef] [Google Scholar]
 Ghez, A. M., Salim, S., Hornstein, S. D., et al. 2005, ApJ, 620, 744 [NASA ADS] [CrossRef] [Google Scholar]
 Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, ApJ, 689, 1044 [NASA ADS] [CrossRef] [Google Scholar]
 Gillessen, S., Eisenhauer, F., Fritz, T. K., et al. 2009a, ApJ, 707, L114 [NASA ADS] [CrossRef] [Google Scholar]
 Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009b, ApJ, 692, 1075 [NASA ADS] [CrossRef] [Google Scholar]
 Gillessen, S., Plewa, P. M., Eisenhauer, F., et al. 2017, ApJ, 837, 30 [NASA ADS] [CrossRef] [Google Scholar]
 Gnedin, O. Y., & Primack, J. R. 2004, Phys. Rev. Lett., 93, 061302 [NASA ADS] [CrossRef] [Google Scholar]
 Gondolo, P., & Silk, J. 1999, Phys. Rev. Lett., 83, 1719 [NASA ADS] [CrossRef] [Google Scholar]
 GRAVITY Collaboration(Abuter, R., et al.) 2018, A&A, 615, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hall, J., & Gondolo, P. 2006, Phys. Rev. D, 74, 063511 [NASA ADS] [CrossRef] [Google Scholar]
 Iorio, L. 2013, Galaxies, 1, 6 [NASA ADS] [CrossRef] [Google Scholar]
 Kass, R. E., & Raftery, A. E. 1995, J. Am. Stat. Assoc., 90, 773 [CrossRef] [MathSciNet] [Google Scholar]
 Lacroix, T., Bœhm, C., & Silk, J. 2014, Phys. Rev. D, 89, 063534 [NASA ADS] [CrossRef] [Google Scholar]
 Lacroix, T., Bœhm, C., & Silk, J. 2015, Phys. Rev. D, 92, 043510 [NASA ADS] [CrossRef] [Google Scholar]
 Lacroix, T., Karami, M., Broderick, A. E., Silk, J., & Bœhm, C. 2017, Phys. Rev. D, 96, 063008 [NASA ADS] [CrossRef] [Google Scholar]
 McMillan, P. J. 2017, MNRAS, 465, 76 [NASA ADS] [CrossRef] [Google Scholar]
 Merritt, D., & Poon, M. Y. 2004, ApJ, 606, 788 [NASA ADS] [CrossRef] [Google Scholar]
 Merritt, D., Milosavljević, M., Verde, L., & Jimenez, R. 2002, Phys. Rev. Lett., 88, 191301 [NASA ADS] [CrossRef] [Google Scholar]
 Nakano, T., & Makino, J. 1999, ApJ, 525, L77 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ostriker, J. P. 2000, Phys. Rev. Lett., 84, 5258 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Plewa, P. M., Gillessen, S., Eisenhauer, F., et al. 2015, MNRAS, 453, 3234 [NASA ADS] [CrossRef] [Google Scholar]
 Portail, M., Gerhard, O., Wegg, C., & Ness, M. 2017, MNRAS, 465, 1621 [NASA ADS] [CrossRef] [Google Scholar]
 Regis, M., & Ullio, P. 2008, Phys. Rev. D, 78, 043505 [NASA ADS] [CrossRef] [Google Scholar]
 Reid, M. J., & Brunthaler, A. 2004, ApJ, 616, 872 [NASA ADS] [CrossRef] [Google Scholar]
 Sadeghian, L., Ferrer, F., & Will, C. M. 2013, Phys. Rev. D, 88, 063522 [NASA ADS] [CrossRef] [Google Scholar]
 Schödel, R., Ott, T., Genzel, R., et al. 2002, Nature, 419, 694 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Schödel, R., GallegoCano, E., Dong, H., et al. 2018, A&A, 609, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shapiro, S. L., & Shelton, J. 2016, Phys. Rev. D, 93, 123510 [NASA ADS] [CrossRef] [Google Scholar]
 Ullio, P., Zhao, H., & Kamionkowski, M. 2001, Phys. Rev. D, 64, 043504 [NASA ADS] [CrossRef] [Google Scholar]
 Weinberg, N. N., Milosavljević, M., & Ghez, A. M. 2005, ApJ, 622, 878 [NASA ADS] [CrossRef] [Google Scholar]
 Wyse, R. F. G. 2001, ASP Conf. Ser., 230, 71 [NASA ADS] [Google Scholar]
 Zakharov, A. F., Nucita, A. A., de Paolis, F., & Ingrosso, G. 2007, Phys. Rev. D, 76, 062001 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Bestfit parameters
Here for completeness I provide the posterior probability distributions obtained for the BHonly model and the BH+spike model with an outer halo of slope γ = 0.25.
Fig. A.1. Marginalized posterior probability density functions for the 13 parameters of the BHonly model, using the entire VLT data set up to 2016. This scatterplot matrix was produced using the corner.py Python module (ForemanMackey 2016). 

Open with DEXTER 
Fig. A.2. Marginalized posterior probability density functions for the 14 parameters of the BH+spike model (for fixed halo slope γ = 0.25), using the entire VLT data set up to 2016. This scatterplot matrix was produced using the corner.py Python module (ForemanMackey 2016). 

Open with DEXTER 
Appendix B: Constrained DM halo models used in this work
Parameters of the generalized NFW profiles for the Milky Way mass models constrained by the analysis of McMillan (2017).
All Tables
Parameters of the generalized NFW profiles for the Milky Way mass models constrained by the analysis of McMillan (2017).
All Figures
Fig. 1. Left panel: density profiles for a generalized NFW halo (dashed lines) with γ = 1 and the same halo with a spike in the central region (solid lines) for a nonannihilating DM candidate (purple) and a selfannihilating 1 TeV DM candidate with ⟨σv⟩ ∼ 10^{−30} cm^{3} s^{−1}, ⟨σv⟩ ∼ 10^{−27} cm^{3} s^{−1}, and ⟨σv⟩ ∼ 10^{−26} cm^{3} s^{−1} (orange, blue, and red, respectively). Right panel: corresponding mass profiles, with the same line styles as in the left panel. The profiles are shown for illustration purposes for a spike radius R_{sp} ∼ 100 pc corresponding to the 99.7% upper limit from deviations of the BHonly orbit using the VLT data (see Sect. 3). The horizontal solid and dotted lines represent the combined 2009 and 2017 constraints, respectively. The vertical dotdashed line marks the characteristic size of the orbit of S2. 

Open with DEXTER  
In the text 
Fig. 2. Median of the marginalized posterior distribution of the spike radius R_{sp} (black dashed), and associated 95% (red shaded) and 99.7% (light red shaded) confidence contours as a function of the slope γ of the corresponding dynamically constrained outer halos, using the combined VLT + Keck data up to 2009 (left panel), VLTonly data up to 2016 (middle panel), and VLT + Keck data up to 2016 (right panel). The upper right white regions correspond to excluded values of the spike radius. The lower left white regions are outside the 99.7% bestfit contours, but cannot robustly be interpreted in terms of exclusion of the associated values of the spike radius since the Bayesian evidence does not favour the presence of a spike over the BHonly model. The contours are derived for values of the slope of the outer halo dynamically constrained by McMillan (2017), and linearly interpolated to get smooth curves. The prediction from Gondolo & Silk (1999) is also shown as a benchmark model (cyan solid). 

Open with DEXTER  
In the text 
Fig. A.1. Marginalized posterior probability density functions for the 13 parameters of the BHonly model, using the entire VLT data set up to 2016. This scatterplot matrix was produced using the corner.py Python module (ForemanMackey 2016). 

Open with DEXTER  
In the text 
Fig. A.2. Marginalized posterior probability density functions for the 14 parameters of the BH+spike model (for fixed halo slope γ = 0.25), using the entire VLT data set up to 2016. This scatterplot matrix was produced using the corner.py Python module (ForemanMackey 2016). 

Open with DEXTER  
In the text 