Reducing theoretical uncertainties in Galactic dark matter (DM) searches is an important challenge as several experiments are now delving into the parameter space relevant to popular (particle or not) candidates. Since many DM signal predictions rely on the knowledge of the DM velocity distribution---direct searches, capture by stars, $p$-wave-suppressed or Sommerfeld-enhanced annihilation rate, microlensing of primordial black holes, etc.---it is necessary to assess the accuracy of our current theoretical handle. Beyond Maxwellian approximations or ad-hoc extrapolations of fits on cosmological simulations, approaches have been proposed to self-consistently derive the DM phase-space distribution only from the detailed mass content of the Galaxy and some symmetry assumptions (e.g. the Eddington inversion and its anisotropic extensions). Although theoretically sound, these methods are still based on simplifying assumptions and their relevance to real galaxies can be questioned. In this paper, we use zoomed-in cosmological simulations to quantify the associated uncertainties. Assuming isotropy, we predict the speed distribution and its moments from the DM and baryonic content measured in simulations, and compare them with the true ones. We reach a predictivity down to ~ 10% for some observables, significantly better than some Maxwellian models. This moderate theoretical error is particularly encouraging at a time when stellar surveys like the Gaia mission should allow us to improve constraints on Galactic mass models.

Dark matter subhalos, predicted in large numbers in the cold dark matter scenario, should have an impact on particle dark matter searches. Recent results show that tidal disruption of these objects in computer simulations is over-efficient due to numerical artifacts and resolution effects. Accounting for these results, we re-estimate the subhalo abundance in the Milky Way using semi-analytical techniques. In particular, we show that the boost factor for gamma rays and cosmic-ray antiprotons is increased by roughly a factor of two.

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 $m$. The constraints we derive exclude the solitons predicted by a naive extrapolation of the soliton-halo relation, found in DM-only numerical simulations, for $2\times10^{-20}~{\rm eV}\lesssim m\lesssim8\times10^{-19}~{\rm eV}$ (from Sgr A*) and $m\lesssim4\times10^{-22}~{\rm eV}$ (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 over-predict the soliton mass by orders of magnitude.

Analysing 9 years of Fermi-LAT observations, we recently studied the spectral properties of the prominent globular cluster 47 Tuc (Brown et al. 2018). In particular, we investigated several models to explain the observed gamma-ray emission, ranging from millisecond pulsars (MSP) to Dark Matter (DM), with the motivation for the latter model driven by recent evidence that 47 Tuc harbours an intermediate-mass black hole (IMBH). This investigation found evidence that the observed gamma-ray emission from 47 Tuc is due to two source populations of MSPs and DM. In Bartels \& Edwards (2018), the authors comment that this evidence is an artifact of the MSP spectra used in (Brown et al. 2019). Here we reply to this comment and argue that the authors of Bartels \& Edwards (2018) (i) do not give due consideration to a very important implication of their result and (ii) there is tension between our MSP fit and their MSP fit when taking uncertainties into consideration. As such, we still conclude there is evidence for a DM component which motivates a deeper radio study of the prominent globular cluster 47 Tuc.

Dark matter (DM) is the most abundant material in the Universe, but has so far been detected only via its gravitational effects. Several theories suggest that pairs of DM particles can annihilate into a flash of light at gamma-ray wavelengths. While gamma-ray emission has been observed from environments where DM is expected to accumulate, such as the centre of our Galaxy, other high energy sources can create a contaminating astrophysical gamma-ray background, thus making DM detection difficult. In principle, dwarf galaxies around the Milky Way are a better place to look -- they contain a greater fraction of DM with no astrophysical gamma-ray background -- but they are too distant for gamma-rays to have been seen. A range of observational evidence suggests that Omega Centauri (omega Cen or NGC 5139), usually classified as the Milky Way's largest globular cluster, is really the core of a captured and stripped dwarf galaxy. Importantly, Omega Cen is ten times closer to us than known dwarfs. Here we show that not only does Omega Cen contain DM with density as high as compact dwarf galaxies, but also that it emits gamma-rays with an energy spectrum matching that expected from the annihilation of DM particles with mass 31$\pm$4 GeV (68\% confidence limit). No astrophysical sources have been found that would otherwise explain Omega Cen's gamma-ray emission, despite deep multi-wavelength searches. We anticipate our results to be the starting point for even deeper radio observations of Omega Cen. If multi-wavelength searches continue to find no astrophysical explanations, this pristine, nearby clump of DM will become the best place to study DM interactions through forces other than gravity.

We recently proposed a method to constrain $s$-wave annihilating MeV dark matter from a combination of the Voyager 1 and the AMS-02 data on cosmic-ray electrons and positrons. Voyager 1 actually provides an unprecedented probe of dark matter annihilation to cosmic rays down to $\sim 10$ MeV in an energy range where the signal is mostly immune to uncertainties in cosmic-ray propagation. In this article, we derive for the first time new constraints on $p$-wave annihilation down to the MeV mass range using cosmic-ray data. To proceed, we derive a self-consistent velocity distribution for the dark matter across the Milky Way by means of the Eddington inversion technique and its extension to anisotropic systems. As inputs, we consider state-of-the-art Galactic mass models including baryons and constrained on recent kinematic data. They allow for both a cored or a cuspy halo. We then calculate the flux of cosmic-ray electrons and positrons induced by $p$-wave annihilating dark matter and obtain very stringent limits in the MeV mass range, robustly excluding cross sections greater than $\sim 10^{-22}{\rm cm^3/s}$ (including theoretical uncertainties), about 5 orders of magnitude better than current CMB constraints. This limit assumes that dark matter annihilation is the sole source of cosmic rays and could therefore be made even more stringent when reliable models of astrophysical backgrounds are included.

Irrespective of the dark matter (DM) candidate, several potentially observable signatures derive from the velocity distribution of DM in halos, in particular in the Milky Way (MW) halo. Examples include direct searches for weakly-interacting massive particles (WIMPs), $p$-wave suppressed or Sommerfeld-enhanced annihilation signals, microlensing events of primordial black holes (PBHs), {\em etc}. Most current predictions are based on the Maxwellian approximation which is not only theoretically inconsistent in bounded systems, but also not supported by cosmological simulations. A more consistent method sometimes used in calculations for direct WIMP searches relies on the so-called Eddington inversion method, which relates the DM phase-space distribution function (DF) to its mass density profile and the total gravitational potential of the system. Originally built upon the isotropy assumption, this method can be extended to anisotropic systems. We investigate these inversion methods in the context of Galactic DM searches, motivated by the fact that the MW is a strongly constrained system, and should be even more so with the ongoing Gaia survey. We still draw conclusions that apply to the general case. In particular, we illustrate how neglecting the radial boundary of the DM halo leads to theoretical inconsistencies. We also show that several realistic configurations of the DM halo and the MW baryonic content entail ill-defined DFs, significantly restricting the configuration space over which these inversion methods can apply. We propose consistent solutions to these issues. Finally, we compute several observables inferred from constrained Galactic mass models relevant to DM searches (WIMPs or PBHs), {\em e.g.} moments and inverse moments of the DM speed and relative speed distributions.

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 non-annihilating DM 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.

47 Tuc was the first globular cluster observed to be $\gamma$-ray bright, with the $\gamma$-rays being attributed to a population of unresolved millisecond pulsars (MSPs). Recent kinematic data, combined with detailed simulations, appears to be consistent with the presence of an intermediate mass black hole (IMBH) at the centre of 47 Tuc. Building upon this, we analyse 9 years of \textit{Fermi}-LAT observations to study the spectral properties of 47 Tuc with unprecedented accuracy and sensitivity. This 9-year $\gamma$-ray spectrum shows that 47 Tuc's $\gamma$-ray flux cannot be explained by MSPs alone, due to a systematic discrepancy between the predicted and observed flux. Rather, we find a significant preference (TS $=40$) for describing 47 Tuc's spectrum with a two source population model, consisting of an ensemble of MSPs and annihilating dark matter (DM) with an enhanced density around the IMBH, when compared to an MSP-only explanation. The best-fit DM mass of 34 GeV is essentially the same as the best-fit DM explanation for the Galactic centre "excess" when assuming DM annihilation into $b\bar{b}$ quarks. Our work constitutes the first possible evidence of dark matter within a globular cluster.

Could there be a large population of intermediate-mass black holes (IMBHs) formed in the early universe? Whether primordial or formed in Population III, these are likely to be very subdominant compared to the dark matter density, but could seed early dwarf galaxy/globular cluster and supermassive black hole formation. Via survival of dark matter density spikes, we show here that a centrally concentrated relic population of IMBHs, along with ambient dark matter, could account for the Fermi gamma-ray "excess" in the Galactic center because of dark matter particle annihilations.

We demonstrate the unprecedented capabilities of the Event Horizon Telescope (EHT) to image the innermost dark matter profile in the vicinity of the supermassive black hole at the center of the M87 radio galaxy. We present the first model of the synchrotron emission induced by dark matter annihilations from a spiky profile in the close vicinity of a supermassive black hole, accounting for strong gravitational lensing effects. Our results show that the EHT should readily resolve dark matter spikes if present. Moreover, the photon ring surrounding the silhouette of the black hole is clearly visible in the spike emission, which introduces observable small-scale structure into the signal. We find that the dark matter-induced emission provides an adequate fit to the existing EHT data, implying that in addition to the jet, a dark matter spike may account for a sizable portion of the millimeter emission from the innermost (subparsec) region of M87. Regardless, our results show that the EHT can probe very weakly annihilating dark matter. Current EHT observations already constrain very small cross sections, typically down to a few 10^{-31} cm^{3} s^{-1} for a 10 GeV candidate, close to characteristic values for p-wave-suppressed annihilation. Future EHT observations will further improve constraints on the DM scenario.

We report the discovery of a statistically significant hardening in the Fermi-LAT $\gamma$-ray spectrum of Centaurus A's core, with the spectral index hardening from $\Gamma_{1}=2.73 \pm 0.02$ to $\Gamma_{1}=2.29 \pm 0.07$ at a break energy of ($2.6 \pm 0.3$) GeV. Using a likelihood analysis, we find no evidence for flux variability in Cen A's core lightcurve above or below the spectral break when considering the entire 8 year period. Interestingly, however, the first $\sim3.5$ years of the low energy lightcurve shows evidence of flux variability at the $\sim3.5 \sigma$ confidence level. To understand the origin of this spectral break, we assume that the low energy component below the break feature originates from leptons in Centaurus A's radio jet and we investigate the possibility that the high energy component above the spectral break is due to an additional source of very high energy particles near the core of Cen A. We show for the first time that the observed $\gamma$-ray spectrum of an Active Galactic Nucleus is compatible with either a very large localized enhancement (referred to as a spike) in the dark matter halo profile or a population of millisecond pulsars. Our work constitutes the first robust indication that new $\gamma$-ray production mechanisms can explain the emission from active galaxies and could provide tantalizing first evidence for the clustering of heavy dark matter particles around black holes.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

The H.E.S.S. collaboration has reported a high-energy spherically symmetric diffuse gamma-ray emission in the inner 50 pc of the Milky Way, up to ~ 50 TeV. Here we propose a leptonic model which provides an alternative to the hadronic scenario presented by the H.E.S.S. collaboration, and connects the newly reported TeV emission to the Fermi-LAT Galactic center GeV excess. Our model relies on a combination of inverse Compton emission from a population of millisecond pulsars---which can account for the GeV excess---and a supermassive black hole-induced spike of heavy (~ 60 TeV) dark matter particles annihilating into electrons with a sub-thermal cross-section. With an up-to-date interstellar radiation field, as well as a standard magnetic field and diffusion set-up, our model accounts for the spectral morphology of the detected emission. Moreover, we show that the dark matter induced emission reproduces the spatial morphology of the H.E.S.S. signal above ~ 10 TeV, while we obtain a slightly more extended component from pulsars at lower energies, which could be used as a prediction for future H.E.S.S. observations.

Excess GeV gamma rays from the Galactic Center (GC) have been measured with the Fermi Large Area Telescope (LAT). The presence of the GC excess (GCE) appears to be robust with respect to changes in the diffuse galactic background modeling. The three main proposals for the GCE are an unresolved population of millisecond pulsars (MSPs), outbursts of cosmic rays from the GC region, and self-annihilating dark matter (DM). The injection of secondary electrons and positrons into the interstellar medium (ISM) by an unresolved population of MSPs or DM annihilations can lead to observable gamma-ray emission via inverse Compton scattering or bremsstrahlung. Here we investigate how to determine whether secondaries are important in a model for the GCE. We develop a method of testing model fit which accounts for the different spatial morphologies of the secondary emission. We examine several models which give secondary emission and illustrate a case where a broadband analysis is not sufficient to determine the need for secondary emission.

Using the spectral energy distribution of M87, a nearby radio galaxy in the Virgo cluster, and assuming a supermassive black hole induced spike in the dark matter halo profile, we exclude any dark matter candidate with a velocity-independent (s-wave) annihilation cross-section of the order of sigma v ~ 10^{-26} cm^3/s and a mass up to O(100) TeV. These limits supersede all previous constraints on thermal, s-wave, annihilating dark matter candidates by orders of magnitude, and rule out the entire canonical mass range. We remark in addition that, under the assumption of a spike, dark matter particles with a mass of a few TeV and an annihilation cross-section of ~ 10^{-27} cm^3/s could explain the TeV gamma-ray emission observed in M87. A central dark matter spike is plausibly present around the supermassive black hole at the center of M87, for various, although not all, formation scenarios, and would have profound implications for our understanding of the dark matter microphysics.

An excess of gamma rays at GeV energies has been detected in the Fermi-LAT data. This signal comes from a narrow region around the Galactic Center and has been interpreted as possible evidence for light (30 GeV) dark matter particles. Focussing on the prompt gamma-ray emission, previous works found that the best fit to the data corresponds to annihilations proceeding into b quarks, with a dark matter profile going as r^{-1.2}. We show that this is not the only possible annihilation set-up. More specifically, we show how including the contributions to the gamma-ray spectrum from inverse Compton scattering and bremsstrahlung from electrons produced in dark matter annihilations, and undergoing diffusion through the Galactic magnetic field, significantly affects the spectrum for leptonic final states. This drastically changes the interpretation of the excess in terms of dark matter.

An excess of gamma rays at GeV energies has been pointed out in the Fermi-LAT data. This signal comes from a narrow region centred around the Galactic center and has been interpreted as possible evidence for light dark matter particles annihilating either into a mixture of leptons-antileptons and b-bbar or into b-bbar only. Focusing on the prompt gamma-ray emission, previous works found that the best fit to the data corresponds to annihilations proceeding predominantly into b-bbar. However, here we show that omitting the photon emission originating from primary and secondary electrons produced in dark matter annihilations, and undergoing diffusion through the Galactic magnetic field, can actually lead to the wrong conclusion. Accounting for this emission, we find that not only are annihilations of ~ 10 GeV particles into a purely leptonic final state allowed, but the democratic scenario actually provides a better fit to the spectrum of the excess than the pure b-bbar channel. We conclude our work with a discussion on constraints on these leptophilic scenarios based on the AMS data and the morphology of the excess.

The morphology and characteristics of the so-called GeV gamma-ray excess detected in the Milky Way lead us to speculate about a possible common origin with the 511 keV line mapped by the SPI experiment about ten years ago. In the previous version of our paper, we assumed 30 GeV dark matter particles annihilating into $b \bar{b}$ and obtained both a morphology and a 511 keV flux (\phi_{511 keV} ~ 10^{-3} ph/cm^2/s) in agreement with SPI observation. However our estimates assumed a negligible number density of electrons in the bulge which lead to an artificial increase in the flux (mostly due to negligible Coulomb losses in this configuration). Assuming a number density greater than $n_e > 10^{-3} cm^{-3}$, we now obtain a flux of 511 keV photons that is smaller than \phi_{511 keV} ~ 10^{-6} ph/cm^2/s and is essentially in agreement with the 511 keV flux that one can infer from the total number of positrons injected by dark matter annihilations into $b \bar{b}$. We thus conclude that -- even if 30 GeV dark matter particles were to exist-- it is impossible to establish a connexion between the two types of signals, even though they are located within the same 10 deg region in the galactic centre.

The dark matter halo profile in the inner Galaxy is very uncertain. Yet its radial dependence toward the Galactic Center is of crucial importance for the determination of the gamma-ray and radio fluxes originating from dark matter annihilations. Here we use synchrotron emission to probe the dark matter energy distribution in the inner Galaxy. We first solve the problem of the cosmic ray diffusion on very small scales, typically smaller than 10^{-3} pc, by using a Green's function approach and use this technique to quantify the effect of a spiky profile (rho(r) ~ r^{-7/3}) on the morphology and intensity of the synchrotron emission expected from dark matter. We illustrate our results using 10 and 800 GeV candidate weakly interacting dark matter particles annihilating directly into e+ e-. Our most critical assumptions are that the dark matter is heavier than a few GeV and directly produces a reasonable amount of electrons and positrons in the Galaxy. We conclude that dark matter indirect detection techniques (including the Planck experiment) could be used to shed light on the dark matter halo profile on scales that lie beyond the capability of any current numerical simulations.

We investigate constraints on the distribution of dark matter in the neighbourhood of the Galactic centre that may eventually be attained with the high-resolution Event Horizon Telescope (EHT). The shadow of a black hole in vacuum is used to generate a toy model describing how dark matter affects the size of the shadow of the supermassive black hole located at the Galactic centre. Observations by the EHT may constrain the properties of the dark matter distribution in a possible density spike around the black hole. Current uncertainties due to both the resolution of the telescope and the analysis of stellar orbits prevent one from discerning the effect of dark matter on the measured size of the shadow. The change in the size of the shadow induced by dark matter can be seen as an additional uncertainty in any test of general relativity that relies on using the angular size of the shadow to estimate the Schwarzschild radius of the black hole.

In this paper we put forward a proposal to use Near Earth Objects as radiation shield for deep space exploration. In principle these objects can provide also a spacious habitat for the astronauts and their supplies on their journeys. We undertake also a detailed assessment of this proposal for a mission from Earth to Mars.

In 2011, the XENON100 experiment has set unprecedented constraints on dark matter-nucleon interactions, excluding dark matter candidates with masses down to 6 GeV if the corresponding cross section is larger than 10^{-39} cm^2. The dependence of the exclusion limit in terms of the scintillation efficiency (Leff) has been debated at length. To overcome possible criticisms XENON100 performed an analysis in which Leff was considered as a nuisance parameter and its uncertainties were profiled out by using a Gaussian likelihood in which the mean value corresponds to the best fit Leff value smoothly extrapolated to zero below 3 keVnr. Although such a method seems fairly robust, it does not account for more extreme types of extrapolation nor does it enable to anticipate on how much the exclusion limit would vary if new data were to support a flat behaviour for Leff below 3 keVnr, for example. Yet, such a question is crucial for light dark matter models which are close to the published XENON100 limit. To answer this issue, we use a maximum Likelihood ratio analysis, as done by the XENON100 collaboration, but do not consider Leff as a nuisance parameter. Instead, Leff is obtained directly from the fits to the data. This enables us to define frequentist confidence intervals by marginalising over Leff.