A dark matter overdensity around a black hole may significantly alter the dynamics of the black hole's merger with another compact object. We consider here intermediate mass-ratio inspirals of stellar-mass compact objects with intermediate-mass black holes "dressed" with dark matter. We first demonstrate that previous estimates based on a fixed dark-matter dress are unphysical for a range of binaries and dark-matter distributions by showing that the total energy dissipated by the compact object through dynamical friction, as it inspirals through the dense dark matter environment towards the black hole, is larger than the gravitational binding energy of the dark-matter dress itself. We then introduce a new formalism that allows us to self-consistently follow the evolution of the dark-matter dress due to its gravitational interaction with the binary. We show that the dephasing of the gravitational waveform induced by dark matter is smaller than previously thought, but is still potentially detectable with the LISA space interferometer. The gravitational waves from such binaries could provide powerful diagnostics of the particle nature of dark matter.

Observational constraints on gamma rays produced by the annihilation of weakly interacting massive particles around primordial black holes (PBHs) imply that these two classes of Dark Matter candidates cannot coexist. We show here that the successful detection of one or more PBHs by radio searches (with the Square Kilometer Array) and gravitational waves searches (with LIGO/Virgo and the upcoming Einstein Telescope) would set extraordinarily stringent constraints on virtually all weak-scale extensions of the Standard Model with stable relics, including those predicting a WIMP abundance much smaller than that of Dark Matter. Upcoming PBHs searches have in particular the potential to rule out almost the entire parameter space of popular theories such as the minimal supersymmetric standard model and scalar singlet Dark Matter.

We analyze the distribution of stars along the GD-1 stream with a combination of data from the ${\it Gaia}$ satellite and the Pan-STARRS survey, and we show that the population of subhalos predicted by the cold dark matter paradigm are necessary and sufficient to explain the perturbations observed in the linear density of stars. This allows us to set novel constraints on alternative dark matter scenarios that predict a suppression of the subhalo mass function on scales smaller than the mass of dwarf galaxies. A combined analysis of the density perturbations in the GD-1 and Pal 5 streams leads to a $95\%$ lower limit on the mass of warm dark matter thermal relics $m_{\rm WDM}>4.6$ keV; adding dwarf satellite counts strengthens this to $m_{\rm WDM}>6.3$ keV.

New data from the $\textit{Gaia}$ satellite, when combined with accurate photometry from the Pan-STARRS survey, allow us to accurately estimate the properties of the GD-1 stream. Here, we analyze the stellar density perturbations in the GD-1 stream and show that they cannot be due to known baryonic structures like giant molecular clouds, globular clusters, or the Milky Way's bar or spiral arms. A joint analysis of the GD-1 and Pal 5 streams instead requires a population of dark substructures with masses $\approx 10^{7}$ to $10^9 \ M_{\rm{\odot}}$. We infer a total abundance of dark subhalos normalised to standard cold dark matter $n_{\rm sub}/n_{\rm sub, CDM} = 0.4 ^{+0.3}_{-0.2}$ ($68 \%$), which corresponds to a mass fraction contained in the subhalos $f_{\rm{sub}} = 0.14 ^{+0.11}_{-0.07} \%$, compatible with the predictions of hydrodynamical simulation of cold dark matter with baryons.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of encouraging a strong community effort at the interface between these two exciting fields of research.

If primordial black holes of $\mathcal{O}(1-100) M_{\odot}$ constitute a significant portion of the dark matter in the Universe, they should be very abundant in our Galaxy. We present here a detailed analysis of the radio and X-ray emission that these objects are expected to produce due to the accretion of gas from the interstellar medium. With respect to previous studies, we relax the assumption of a monochromatic mass function, and introduce an improved treatment of the physics of gas accretion onto isolated, moving compact objects, based on a set of state-of-the-art numerical simulations. By comparing our predictions with known radio and X-ray sources in the Galactic center region, we show that the maximum relic density of primordial black holes in the mass range of interest is $\sim 10^{-3}$ smaller than that of dark matter. The new upper bound is two orders of magnitude stronger with respect to previous results, based on a conservative phenomenological treatment of the accretion physics. We also provide a comprehensive critical discussion on the reliability of this bound, and on possible future developments in the field. We argue in particular that future multi-wavelength searches will soon start to probe the galactic population of astrophysical black holes.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

The formation of astrophysical and primordial black holes influences the distribution of dark matter surrounding them. Black holes are thus expected to carry a dark matter `dress' whose properties depend on their formation mechanism and on the properties of the environment. Here we carry out a numerical and analytical study of the merger of dressed black holes, and show that the distribution of dark matter around them dramatically affects the dynamical evolution of the binaries. Although the final impact on the merger rate of primordial black holes is rather small with respect to the case of `naked' black holes, we argue that our analysis places the calculation of this rate on more solid ground, with LIGO-Virgo observations ruling out a dark matter fraction of $10^{-3}$ for primordial black holes of 100 solar masses, and it paves the way to more detailed analyses of environmental effects induced by dark matter on the gravitational wave emission of binary black holes.

There is a growing sense of `crisis' in the dark matter community, due to the absence of evidence for the most popular candidates such as weakly interacting massive particles, axions, and sterile neutrinos, despite the enormous effort that has gone into searching for these particles. Here, we discuss what we have learned about the nature of dark matter from past experiments, and the implications for planned dark matter searches in the next decade. We argue that diversifying the experimental effort, incorporating astronomical surveys and gravitational wave observations, is our best hope to make progress on the dark matter problem.

A key prediction of the standard cosmological model -- which relies on the assumption that dark matter is cold, i.e. non-relativistic at the epoch of structure formation -- is the existence of a large number of dark matter substructures on sub-galactic scales. This assumption can be tested by studying the perturbations induced by dark matter substructures on cold stellar streams. Here, we study the prospects for discriminating cold from warm dark matter by generating mock data for upcoming astronomical surveys such as the Large Synoptic Survey Telescope (LSST), and reconstructing the properties of the dark matter particle from the perturbations induced on the stellar density profile of a stream. We discuss the statistical and systematic uncertainties, and show that the method should allow to set stringent constraints on the mass of thermal dark matter relics, and possibly to yield an actual measurement of the dark matter particle mass if it is in the $\mathcal{O}(1)$ keV range.

One of the most promising strategies to identify the nature of dark matter consists in the search for new particles at accelerators and with so-called direct detection experiments. Working within the framework of simplified models, and making use of machine learning tools to speed up statistical inference, we address the question of what we can learn about dark matter from a detection at the LHC and a forthcoming direct detection experiment. We show that with a combination of accelerator and direct detection data, it is possible to identify newly discovered particles as dark matter, by reconstructing their relic density assuming they are weakly interacting massive particles (WIMPs) thermally produced in the early Universe, and demonstrating that it is consistent with the measured dark matter abundance. An inconsistency between these two quantities would instead point either towards additional physics in the dark sector, or towards a non-standard cosmology, with a thermal history substantially different from that of the standard cosmological model.

We derive the local dark matter density by applying the integrated Jeans equation method from Silverwood et al. (2016) to SDSS-SEGUE G-dwarf data processed and presented by B\"udenbender et al. (2015). We use the MultiNest Bayesian nested sampling software to fit a model for the baryon distribution, dark matter and tracer stars, including a model for the 'tilt term' that couples the vertical and radial motions, to the data. The $\alpha$-young population from B\"udenbender et al. (2015) yields the most reliable result of $\rho_{\rm DM} = 0.46^{+0.07}_{-0.09}\, {{\rm GeV\, cm}^{-3}} = 0.012^{+0.001}_{-0.002}\, {{\rm M}_\odot \, {\rm pc}^{-3}}$. Our analyses yield inconsistent results for the $\alpha$-young and $\alpha$-old data, pointing to problems in the tilt term and its modelling, the data itself, the assumption of a flat rotation curve, or the effects of disequilibria.

In recent years, realistic hydrodynamical simulations of galaxies like the Milky Way have become available, enabling a reliable estimate of the dark matter density and velocity distribution in the Solar neighborhood. We review here the status of hydrodynamical simulations and their implications for the interpretation of direct dark matter searches. We focus in particular on: the criteria to identify Milky Way-like galaxies; the impact of baryonic physics on the dark matter velocity distribution; the possible presence of substructures like clumps, streams, or dark disks; and on the implications for the direct detection of dark matter with standard and non-standard interactions.

We model the accretion of gas onto a population of massive primordial black holes in the Milky Way, and compare the predicted radio and X-ray emission with observational data. We show that under conservative assumptions on the accretion process, the possibility that ${\cal O}(10) \, M_\odot$ primordial black holes can account for all of the dark matter in the Milky Way is excluded at $5\sigma$ by a comparison with a VLA radio catalog at $1.4$ GHz, and at $\simeq 40\sigma$ by a comparison with a Chandra X-ray catalog ($0.5 - 8$ keV). We argue that this method can be used to identify such a population of primordial black holes with more sensitive future radio and X-ray surveys.

The history of the dark matter problem can be traced back to at least the 1930s, but it was not until the early 1970s that the issue of 'missing matter' was widely recognized as problematic. In the latter period, previously separate issues involving missing mass were brought together in a single anomaly. We argue that reference to a straightforward 'accumulation of evidence' alone is inadequate to comprehend this episode. Rather, the rise of cosmological research, the accompanying renewed interest in the theory of relativity and changes in the manpower division of astronomy in the 1960s are key to understanding how dark matter came to matter. At the same time, this story may also enlighten us on the methodological dimensions of past practices of physics and cosmology.

It was recently proposed that the effects usually attributed to particle dark matter on galaxy scales are due to the displacement of dark energy by baryonic matter, a paradigm known as emergent gravity. This formalism leads to predictions similar to Modified Newtonian Dynamics (MOND) in spherical symmetry, but not quite identical. In particular, it leads to a well defined transition between the Newtonian and the modified gravitational regimes, a transition depending on both the Newtonian acceleration and its first derivative with respect to radius. Under the hypothesis of the applicability of this transition to aspherical systems, we investigate whether it can reproduce observed galaxy rotation curves. We conclude that the formula leads to marginally acceptable fits with strikingly low best-fit distances, low stellar mass-to-light ratios, and a low Hubble constant. In particular, some unobserved wiggles are produced in rotation curves because of the dependence of the transition on the derivative of the Newtonian acceleration, leading, even in the most favorable case, to systematically less good fits than MOND. Then, applying the predicted transition from emergent gravity in a regime where it should a priori be applicable, i.e. in spherical symmetry and outside of the bulk of matter, we show that the predictions for the secular advances of Solar System planets' perihelia are discrepant with the data by seven orders of magnitude, ruling out the present emergent gravity weak-field formula with high confidence.

We estimate the sensitivity of the upcoming CTA gamma-ray telescope to DM annihilation at the Galactic centre, improving on previous analyses in a number of significant ways. First, we perform a detailed analyses of all backgrounds, including diffuse astrophysical emission for the first time in a study of this type. Second, we present a statistical framework for including systematic errors and estimate the consequent degradation in sensitivity. These errors may come from e.g. event reconstruction, Monte Carlo determination of the effective area or uncertainty in atmospheric conditions. Third, we show that performing the analysis on a set of suitably optimised regions of interest makes it possible to partially compensate for the degradation in sensitivity caused by systematics and diffuse emission. To probe dark matter with the canonical thermal annihilation cross-section, CTA systematics like non-uniform variations in acceptance over a single field of view must be kept below the 0.3% level, unless the dark matter density rises more steeply in the centre of the Galaxy than predicted by a typical Navarro-Frenk-White or Einasto profile. For a contracted $r^{-1.3}$ profile, and systematics at the 1% level, CTA can probe annihilation to $b\bar{b}$ at the canonical thermal level for dark matter masses between 100 GeV and 10 TeV.

Collider, direct and indirect searches for dark matter have typically little or no sensitivity to weakly interacting massive particles (WIMPs) with masses above a few TeV. This rather unexplored regime can however be probed through the search for distinctive gamma-ray spectral features produced by the annihilation of WIMPs at very high energies. Here we present a dedicated search for gamma-ray boxes -- sharp spectral features that cannot be mimicked by astrophysical sources -- with the upcoming Cherenkov Telescope Array (CTA). Using realistic projections for the instrument performance and detailed background modelling, a profile likelihood analysis is implemented to derive the expected upper limits and sensitivity reach after 100 h of observations towards a $2^\circ\times2^\circ$ region around the Galactic centre. Our results show that CTA will be able to probe gamma-ray boxes down to annihilation cross sections of $10^{-27}-10^{-26}\,\text{cm}^3\text{/s}$ up to tens of TeV. We also identify a number of concrete particle physics models providing thermal dark matter candidates that can be used as target benchmarks in future search campaigns. This constitutes a golden opportunity for CTA to either discover or rule out multi-TeV thermal dark matter in a corner of parameter space where all other experimental efforts are basically insensitive.

Although dark matter is a central element of modern cosmology, the history of how it became accepted as part of the dominant paradigm is often ignored or condensed into a brief anecdotical account focused around the work of a few pioneering scientists. The aim of this review is to provide the reader with a broader historical perspective on the observational discoveries and the theoretical arguments that led the scientific community to adopt dark matter as an essential part of the standard cosmological model.

We study the implications of galaxy formation on dark matter direct detection using high resolution hydrodynamic simulations of Milky Way-like galaxies simulated within the EAGLE and APOSTLE projects. We identify Milky Way analogues that satisfy observational constraints on the Milky Way rotation curve and total stellar mass. We then extract the dark matter density and velocity distribution in the Solar neighbourhood for this set of Milky Way analogues, and use them to analyse the results of current direct detection experiments. For most Milky Way analogues, the event rates in direct detection experiments obtained from the best fit Maxwellian distribution (with peak speed of 223 - 289 km/s) are similar to those obtained directly from the simulations. As a consequence, the allowed regions and exclusion limits set by direct detection experiments in the dark matter mass and spin-independent cross section plane shift by a few GeV compared to the Standard Halo Model, at low dark matter masses. For each dark matter mass, the halo-to-halo variation of the local dark matter density results in an overall shift of the allowed regions and exclusion limits for the cross section. However, the compatibility of the possible hints for a dark matter signal from DAMA and CDMS-Si and null results from LUX and SuperCDMS is not improved.

We present a new method for determining the local dark matter density using kinematic data for a population of tracer stars. The Jeans equation in the $z$-direction is integrated to yield an equation that gives the velocity dispersion as a function of the total mass density, tracer density, and the tilt term that describes the coupling of vertical and radial motions. We then fit a dark matter mass profile to tracer density and velocity dispersion data to derive credible regions on the vertical dark matter density profile. Our method avoids numerical differentiation, leading to lower numerical noise, and is able to deal with the tilt term while remaining one dimensional. In this study we present the method and perform initial tests on idealised mock data. We also demonstrate the importance of dealing with the tilt term for tracers that sample $\gtrsim 1$ kpc above the disc plane. If ignored, this results in a systematic underestimation of the dark matter density.

We present global fits of an effective field theory description of real, and complex scalar dark matter candidates. We simultaneously take into account all possible dimension 6 operators consisting of dark matter bilinears and gauge invariant combinations of quark and gluon fields. We derive constraints on the free model parameters for both the real (five parameters) and complex (seven) scalar dark matter models obtained by combining Planck data on the cosmic microwave background, direct detection limits from LUX, and indirect detection limits from the Fermi Large Area Telescope. We find that for real scalars indirect dark matter searches disfavour a dark matter particle mass below 100 GeV. For the complex scalar dark matter particle current data have a limited impact due to the presence of operators that lead to p-wave annihilation, and also do not contribute to the spin-independent scattering cross- section. Although current data are not informative enough to strongly constrain the theory parameter space, we demonstrate the power of our formalism to reconstruct the theoretical parameters compatible with an actual dark matter detection, by assuming that the excess of gamma rays observed by the Fermi Large Area Telescope towards the Galactic centre is entirely due to dark matter annihilations. Please note that the excess can very well be due to astrophysical sources such as millisecond pulsars. We find that scalar dark matter interacting via effective field theory operators can in principle explain the Galactic centre excess, but that such interpretation is in strong tension with the non-detection of gamma rays from dwarf galaxies in the real scalar case. In the complex scalar case there is enough freedom to relieve the tension.

If dark matter consists of weakly interacting massive particles (WIMPs), dark matter subhalos in the Milky Way could be detectable as gamma-ray point sources due to WIMP annihilation. In this work, we perform an updated study of the detectability of dark matter subhalos as gamma-ray sources with the Fermi Large Area Telescope (Fermi LAT). We use the results of the Via Lactea II simulation, scaled to the Planck 2015 cosmological parameters, to predict the local dark matter subhalo distribution. Under optimistic assumptions for the WIMP parameters --- a 40 GeV particle annihilating to $b\bar{b}$ with a thermal cross-section, as required to explain the Galactic center GeV excess --- we predict that at most $\sim 10$ subhalos might be present in the third Fermi LAT source catalog (3FGL). This is a smaller number than has been predicted by prior studies, and we discuss the origin of this difference. We also compare our predictions for the detectability of subhalos with the number of subhalo candidate sources in 3FGL, and derive upper limits on the WIMP annihilation cross-section as a function of the particle mass. If a dark matter interpretation could be excluded for all 3FGL sources, our constraints would be competitive with those found by indirect searches using other targets, such as known Milky Way satellite galaxies.

We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.

An accurate knowledge of the dark matter distribution in the Milky Way is of crucial importance for galaxy formation studies and current searches for particle dark matter. In this paper we set new dynamical constraints on the Galactic dark matter profile by comparing the observed rotation curve, updated with a comprehensive compilation of kinematic tracers, with that inferred from a wide range of observation-based morphologies of the bulge, disc and gas. The generalised Navarro-Frenk-White (NFW) and Einasto dark matter profiles are fitted to the data in order to determine the favoured ranges of local density, slope and scale radius. For a representative baryonic model, a typical local circular velocity of 230 km/s and a distance of the Sun to the Galactic centre of 8 kpc, we find a local dark matter density of 0.420+0.021-0.018 (2 sigma) +- 0.025 GeV/cm^3 (0.420+0.019-0.021 (2 sigma) +- 0.026 GeV/cm^3) for NFW (Einasto), where the second error is an estimate of the systematic due to baryonic modelling. Apart from the Galactic parameters, the main sources of uncertainty inside and outside the solar circle are baryonic modelling and rotation curve measurements, respectively. Upcoming astronomical observations are expected to reduce all these uncertainties substantially over the coming years.

We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies (APOSTLE project) carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed in the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White (NFW) profile between $r\approx1.5$ kpc and $r\approx10$ kpc. At smaller radii, $r\lesssim1.5$ kpc, the halos develop a flatter than NFW slope. This unexpected feature may be specific to our particular choice of subgrid physics model but nevertheless the dark matter density profiles agree within 30% as the mass resolution is increased by a factor 150. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.97$ within $r=1$ kpc) and there is no offset larger than 45 pc between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $\gamma$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.

Recent studies of M31, the Galactic centre, and galaxy clusters have made tentative detections of an X-ray line at ~3.5 keV that could be produced by decaying dark matter. We use high resolution simulations of the Aquarius project to predict the likely amplitude of the X-ray decay flux observed in the GC relative to that observed in M31, and also of the GC relative to other parts of the Milky Way halo and to dwarf spheroidal galaxies. We show that the reported detections from M31 and Andromeda are compatible with each other, and with upper limits arising from high galactic latitude observations, and imply a decay time {\tau} ~10^28 seconds. We argue that this interpretation can be tested with deep observations of dwarf spheroidal galaxies: in 95 per cent of our mock observations, a 1.3 Msec pointed observation of Draco with XMM-Newton will enable us to discover or rule out at the 3{\sigma} level an X-ray feature from dark matter decay at 3.5 keV, for decay times {\tau} < 0.8x10^28 sec.

Modified Newtonian dynamics (MOND) is an empirical theory originally proposed to explain the rotation curves of spiral galaxies by modifying the gravitational acceleration, rather than by invoking dark matter. Here,we set constraints on MOND using an up-to-date compilation of kinematic tracers of the Milky Way and a comprehensive collection of morphologies of the baryonic component in the Galaxy. In particular, we find that the so-called "standard" interpolating function cannot explain at the same time the rotation curve of the Milky Way and that of external galaxies for any of the baryonic models studied, while the so-called "simple" interpolating function can for a subset of models. Upcoming astronomical observations will refine our knowledge on the morphology of baryons and will ultimately confirm or rule out the validity of MOND in the Milky Way. We also present constraints on MOND-like theories without making any assumptions on the interpolating function.

In a brief note posted recently, the authors of arXiv:1503.07813 raised some concerns on our arxiv:1502.03821, recently published in Nature Physics. We thank them for the interest in our work and respond here to their criticisms.

The ubiquitous presence of dark matter in the universe is today a central tenet in modern cosmology and astrophysics. Ranging from the smallest galaxies to the observable universe, the evidence for dark matter is compelling in dwarfs, spiral galaxies, galaxy clusters as well as at cosmological scales. However, it has been historically difficult to pin down the dark matter contribution to the total mass density in the Milky Way, particularly in the innermost regions of the Galaxy and in the solar neighbourhood. Here we present an up-to-date compilation of Milky Way rotation curve measurements, and compare it with state-of-the-art baryonic mass distribution models. We show that current data strongly disfavour baryons as the sole contribution to the galactic mass budget, even inside the solar circle. Our findings demonstrate the existence of dark matter in the inner Galaxy while making no assumptions on its distribution. We anticipate that this result will compel new model-independent constraints on the dark matter local density and profile, thus reducing uncertainties on direct and indirect dark matter searches, and will shed new light on the structure and evolution of the Galaxy.

The formation of black holes inevitably affects the distribution of dark and baryonic matter in their vicinity, leading to an enhancement of the dark matter density, called spike, and if dark matter is made of WIMPs, to a strong enhancement of the dark matter annihilation rate. Spikes at the center of galaxies like the Milky Way are efficiently disrupted by baryonic processes, but mini-spikes can form and survive undisturbed at the center of dwarf spheroidal galaxies. We show that Fermi LAT satellite data allow to set very stringent limits on the existence of mini-spikes in dwarf galaxies: for thermal WIMPs with mass between 100 GeV and 1 TeV, we obtain a maximum black hole mass between 100 and 1000 $M_\odot$, ruling out black holes masses extrapolated from the M-{\sigma} relationship in a large region of the parameter space. We also performed Monte Carlo simulations of merger histories of black holes in dwarf spheroidals in a scenario where black holes form from the direct collapse of primordial gas in early halos, and found that this specific formation scenario is incompatible at the 84% CL with dark matter being in the form of thermal WIMPs.

We present statistically convergent profile likelihood maps obtained via global fits of a phenomenological Minimal Supersymmetric Standard Model with 15 free parameters (the MSSM-15), based on over 250M points. We derive constraints on the model parameters from direct detection limits on dark matter, the Planck relic density measurement and data from accelerator searches. We provide a detailed analysis of the rich phenomenology of this model, and determine the SUSY mass spectrum and dark matter properties that are preferred by current experimental constraints. We evaluate the impact of the measurement of the anomalous magnetic moment of the muon ($g-2$) on our results, and provide an analysis of scenarios in which the lightest neutralino is a subdominant component of the dark matter. The MSSM-15 parameters are relatively weakly constrained by current data sets, with the exception of the parameters related to dark matter phenomenology ($M_1$, $M_2$, $\mu$), which are restricted to the sub-TeV regime, mainly due to the relic density constraint. The mass of the lightest neutralino is found to be < 1.5 TeV at 99% C.L., but can extend up to 3 TeV when excluding the $g - 2$ constraint from the analysis. Low-mass bino-like neutralinos are strongly favoured, with spin-independent scattering cross-sections extending to very small values, $\sim 10^{-20}$ pb. ATLAS SUSY null searches strongly impact on this mass range, and thus rule out a region of parameter space that is outside the reach of any current or future direct detection experiment. The best-fit point obtained after inclusion of all data corresponds to a squark mass of 2.3 TeV, a gluino mass of 2.1 TeV and a 130 GeV neutralino with a spin-independent cross-section of $2.4 \times 10^{-10}$ pb, which is within the reach of future multi-ton scale direct detection experiments and of the upcoming LHC run at increased centre-of-mass energy.

We investigate the prospects for reconstructing the mass, spin-independent and spin-dependent cross-sections of Dark Matter particles with a combination of future direct detection experiments such as XENON1T, and the IceCube neutrino telescope in the 86-string configuration including the DeepCore array. We quantify the degree of complementarity between the two experiments by adopting realistic values for their exposure, energy threshold and resolution. Starting from benchmark models arising from a supersymmetric model with 25 free parameters, we show that despite the stringent constraints set by the run with 79 strings, IceCube can help break the degeneracies in the Dark Matter cross-section parameter space, even in the unfortunate case where it fails to discover high energy neutrinos from the Sun. We also discuss how the reconstruction of the Dark Matter particle parameters from the combined data sets is affected by the uncertainties associated with the nuclear structure of the target material in case of spin-dependent scattering and those associated with astrophysical quantities such as the Dark Matter density and velocity distribution.

We consider the propagation and energy losses of protons and anti-protons produced by dark matter annihilation at redshifts 100<z<~2000. In the case of dark matter annihilations into quarks, gluons and weak gauge bosons, protons and anti-protons carry about 20% of the energy injected into e^\pm and \gamma's, but their interactions are normally neglected when deriving cosmic microwave background bounds from altered recombination histories. Here, we follow numerically the energy-loss history of typical protons/antiprotons in the cosmological medium. We show that about half of their energy is channeled into photons and e^\pm, and we present a simple prescription to estimate the corresponding strengthening of the cosmic microwave background bounds on the dark matter annihilation cross section.

The nature of dark matter is still unknown and one of the most fundamental scientific mysteries. Although successfully describing large scales, the standard cold dark matter model (CDM) exhibits possible shortcomings on galactic and sub-galactic scales. It is exactly at these highly non-linear scales where strong astrophysical constraints can be set on the nature of the dark matter particle. While observations of the Lyman-$\alpha$ forest probe the matter power spectrum in the mildly non-linear regime, satellite galaxies of the Milky Way provide an excellent laboratory as a test of the underlying cosmology on much smaller scales. Here we present results from a set of high resolution simulations of a Milky Way sized dark matter halo in eight distinct cosmologies: CDM, warm dark matter (WDM) with a particle mass of 2 keV and six different cold plus warm dark matter (C+WDM) models, varying the fraction, $f_{\rm wdm}$, and the mass, $m_{\rm wdm}$, of the warm component. We used three different observational tests based on Milky Way satellite observations: the total satellite abundance, their radial distribution and their mass profile. We show that the requirement of simultaneously satisfying all three constraints sets very strong limits on the nature of dark matter. This shows the power of a multi-dimensional small scale approach in ruling out models which would be still allowed by large scale observations.

A possible solution to the small scale problems of the cold dark matter (CDM) scenario is that the dark matter consists of two components, a cold and a warm one. We perform a set of high resolution simulations of the Milky Way halo varying the mass of the WDM particle ($m_{\rm WDM}$) and the cosmic dark matter mass fraction in the WDM component ($\bar{f}_{\rm W}$). The scaling ansatz introduced in combined analysis of LHC and astroparticle searches postulates that the relative contribution of each dark matter component is the same locally as on average in the Universe (e.g. $f_{\rm W,\odot} = \bar{f}_{\rm W}$). Here we find however, that the normalised local WDM fraction ($f_{\rm W,\odot}$ / $\bar{f}_{\rm W}$) depends strongly on $m_{\rm WDM}$ for $m_{\rm WDM} <$ 1 keV. Using the scaling ansatz can therefore introduce significant errors into the interpretation of dark matter searches. To correct this issue a simple formula that fits the local dark matter densities of each component is provided.

We explore systematic biases in the identification of dark matter in future direct detection experiments and compare the reconstructed dark matter properties when assuming a self-consistent dark matter distribution function and the standard Maxwellian velocity distribution. We find that the systematic bias on the dark matter mass and cross-section determination arising from wrong assumptions for its distribution function is of order ~1\sigma. A much larger systematic bias can arise if wrong assumptions are made on the underlying Milky Way mass model. However, in both cases the bias is substantially mitigated by marginalizing over galactic model parameters. We additionally show that the velocity distribution can be reconstructed in an unbiased manner for typical dark matter parameters. Our results highlight both the robustness of the dark matter mass and cross-section determination using the standard Maxwellian velocity distribution and the importance of accounting for astrophysical uncertainties in a statistically consistent fashion.

We present global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM) and the Non-Universal Higgs Model (NUHM), including the most recent CMS constraint on the Higgs boson mass, 5.8/fb integrated luminosity null Supersymmetry searches by ATLAS, the new LHCb measurement of the Bs to mu+mu- branching ratio and the 7-year WMAP dark matter relic abundance determination. We include the latest dark matter constraints from the XENON100 experiment, marginalising over astrophysical and particle physics uncertainties. We present Bayesian posterior and profile likelihood maps of the highest resolution available today, obtained from up to 350M points. We find that the new constraint on the Higgs boson mass has a dramatic impact, ruling out large regions of previously favoured cMSSM and NUHM parameter space. In the cMSSM, light sparticles and predominantly gaugino-like dark matter with a mass of a few hundred GeV are favoured. The NUHM exhibits a strong preference for heavier sparticle masses and a Higgsino-like neutralino with a mass of 1 TeV. The future ton-scale XENON1T direct detection experiment will probe large portions of the currently favoured cMSSM and NUHM parameter space. The LHC operating at 14 TeV collision energy will explore the favoured regions in the cMSSM, while most of the regions favoured in the NUHM will remain inaccessible. Our best-fit points achieve a satisfactory quality-of-fit, with p-values ranging from 0.21 to 0.35, so that none of the two models studied can be presently excluded at any meaningful significance level.

We study the prospects for studying line features in gamma-ray spectra with upcoming gamma-ray experiments, such as HESS-II, the Cherenkov Telescope Array (CTA), and the GAMMA-400 satellite. As an example we use the narrow feature at 130 GeV seen in public data from the Fermi-LAT satellite. We found that all three experiments should be able to confidently confirm or rule out the presence of this 130 GeV feature. If it is real, it should be confirmed with a confidence level higher than 5 sigma. Assuming it to be a spectral signature of dark matter origin, GAMMA-400, thanks to a projected energy resolution of about 1.5% at 100 GeV, should also be able to resolve both the \gamma\gamma-line and a corresponding Z\gamma- or H\gamma-feature, if the corresponding branching ratio is comparable to that into two photons. It will also allow to distinguish between a gamma-ray line and the similar feature resulting from internal bremsstrahlung photons.

We present new global fits of the cMSSM, including LHC 1/fb integrated luminosity SUSY exclusion limits, recent LHC 5/fb constraints on the mass of the Higgs boson and XENON100 direct detection data. Our analysis fully takes into account astrophysical and hadronic uncertainties that enter the analysis when translating direct detection limits into constraints on the cMSSM parameter space. We provide results for both a Bayesian and a Frequentist statistical analysis. We find that LHC 2011 constraints in combination with XENON100 data can rule out a significant portion of the cMSSM parameter space. Our results further emphasise the complementarity of collider experiments and direct detection searches in constraining extensions of Standard Model physics. The LHC 2011 exclusion limit strongly impacts on low-mass regions of cMSSM parameter space, such as the stau co-annihilation region, while direct detection data can rule out regions of high SUSY masses, such as the Focus-Point region, which is unreachable for the LHC in the near future. We show that, in addition to XENON100 data, the experimental constraint on the anomalous magnetic moment of the muon plays a dominant role in disfavouring large scalar and gaugino masses. We find that, should the LHC 2011 excess hinting towards a Higgs boson at 126 GeV be confirmed, currently favoured regions of the cMSSM parameter space will be robustly ruled out from both a Bayesian and a profile likelihood statistical perspective.

We discuss irreducible statistical limitations of future ton-scale dark matter direct detection experiments. We focus in particular on the coverage of confidence intervals, which quantifies the reliability of the statistical method used to reconstruct the dark matter parameters, and the bias of the reconstructed parameters. We study 36 benchmark dark matter models within the reach of upcoming ton-scale experiments. We find that approximate confidence intervals from a profile-likelihood analysis exactly cover or over-cover the true values of the WIMP parameters, and are hence conservative. We evaluate the probability that unavoidable statistical fluctuations in the data might lead to a biased reconstruction of the dark matter parameters, or large uncertainties on the reconstructed parameter values. We show that this probability can be surprisingly large, even for benchmark models leading to a large event rate of order a hundred counts. We find that combining data sets from two different targets leads to improved coverage properties, as well as a substantial reduction of statistical bias and uncertainty on the dark matter parameters.

We present updated global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM), including the most recent constraints from the ATLAS and CMS detectors at the LHC, as well as the most recent results of the XENON100 experiment. Our robust analysis takes into account both astrophysical and hadronic uncertainties that enter in the calculation of the rate of WIMP-induced recoils in direct detection experiment. We study the consequences for neutralino Dark Matter, and show that current direct detection data already allow to robustly rule out the so-called Focus Point region, therefore demonstrating the importance of particle astrophysics experiments in constraining extensions of the Standard Model of Particle Physics. We also observe an increased compatibility between results obtained from a Bayesian and a Frequentist statistical perspective. We find that upcoming ton-scale direct detection experiments will probe essentially the entire currently favoured region (at the 99% level), almost independently of the statistical approach used. Prospects for indirect detection of the cMSSM are further reduced.

The unfortunate case where the Large Hadron Collider (LHC) fails to discover physics Beyond the Standard Model (BSM) is sometimes referred to as the "Nightmare scenario" of particle physics. We study the consequences of this hypothetical scenario for Dark Matter (DM), in the framework of the constrained Minimal Supersymmetric Standard Model (cMSSM). We evaluate the surviving regions of the cMSSM parameter space after null searches at the LHC, using several different LHC configurations, and study the consequences for DM searches with ton-scale direct detectors and the IceCube neutrino telescope. We demonstrate that ton-scale direct detection experiments will be able to conclusively probe the cMSSM parameter space that would survive null searches at the LHC with 100fb$^{-1}$ of integrated luminosity at 14TeV. We also demonstrate that IceCube (80 strings plus DeepCore) will be able to probe as much as 17% of the currently favoured parameter space after 5 years of observation.

We show that current microlensing and dynamical observations of the Galaxy permit to set interesting constraints on the Dark Matter local density and profile slope towards the galactic centre. Assuming state-of-the-art models for the distribution of baryons in the Galaxy, we find that the most commonly discussed Dark Matter profiles (viz. Navarro-Frenk-White and Einasto) are consistent with microlensing and dynamical observations, while extreme adiabatically compressed profiles are robustly ruled out. When a baryonic model that also includes a description of the gas is adopted, our analysis provides a determination of the local Dark Matter density, \rho_0=0.20-0.56 GeV/cm^3 at 1\sigma, that is found to be compatible with estimates in the literature based on different techniques.

Even if Supersymmetric particles are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the Dark Matter (DM) in the Universe using LHC data alone. We study the complementarity of LHC and DM indirect searches, working out explicitly the reconstruction of the DM properties for a specific benchmark model in the coannihilation region of a 24-parameters supersymmetric model. Combining mock high-luminosity LHC data with present-day null searches for gamma-rays from dwarf galaxies with the Fermi LAT, we show that current Fermi LAT limits already have the capability of ruling out a spurious Wino-like solution that would survive using LHC data only, thus leading to the correct identification of the cosmological solution. We also demonstrate that upcoming Planck constraints on the reionization history will have a similar constraining power, and discuss the impact of a possible detection of gamma-rays from DM annihilation in Draco with a CTA-like experiment. Our results indicate that indirect searches can be strongly complementary to the LHC in identifying the DM particles, even when astrophysical uncertainties are taken into account.

We investigate the reconstruction capabilities of Dark Matter mass and spin-independent cross-section from future ton-scale direct detection experiments using germanium, xenon or argon as targets. Adopting realistic values for the exposure, energy threshold and resolution of Dark Matter experiments which will come online within 5 to 10 years, the degree of complementarity between different targets is quantified. We investigate how the uncertainty in the astrophysical parameters controlling the local Dark Matter density and velocity distribution affects the reconstruction. For a 50 GeV WIMP, astrophysical uncertainties degrade the accuracy in the mass reconstruction by up to a factor of $\sim 4$ for xenon and germanium, compared to the case when astrophysical quantities are fixed. However, combination of argon, germanium and xenon data increases the constraining power by a factor of $\sim 2$ compared to germanium or xenon alone. We show that future direct detection experiments can achieve self-calibration of some astrophysical parameters, and they will be able to constrain the WIMP mass with only very weak external astrophysical constraints.

The injection of secondary particles produced by Dark Matter (DM) annihilation at redshift 100<z<1000 affects the process of recombination, leaving an imprint on Cosmic Microwave Background (CMB) anisotropies. Here we provide a new assessment of the constraints set by CMB data on the mass and self-annihilation cross-section of DM particles. Our new analysis includes the most recent WMAP (7-year) and ACT data, as well as an improved treatment of the time-dependent coupling between the DM annihilation energy with the thermal gas. We show in particular that the improved measurement of the polarization signal places already stringent constraints on light DM particles, ruling out 'thermal' WIMPs with mass less then about 10 GeV.

We study the prospects for discriminating between the dark matter (DM) and pulsar origin of the PAMELA positron excess with the Alpha Magnetic Spectrometer AMS-02. We simulate the response of AMS-02 to positrons (and electrons) originating from DM annihilations, and determine the pulsar parameters (spin-down luminosity, distance and characteristic age) that produce a satisfactory fit to the mock AMS-02 data. It turns out that it is always possible to mimic a DM signal with pulsars. Although the fit in some cases requires values of spin-down luminosity and characteristic age different from those of known pulsars in the ATNF and Fermi-LAT catalogues, these catalogues are known to be incomplete, and therefore the pulsar interpretation can hardly be ruled out. We also show that if the positron excess is due to a single pulsar, it is always possible to find a DM candidate that provides a good fit to the mock AMS-02 data. The discrimination between the two scenarios will thus require a better knowledge of the underlying sources, or complementary data.

In theories with Universal Extra-Dimensions (UED), the gamma_1 particle, first excited state of the hypercharge gauge boson, provides an excellent Dark Matter (DM) candidate. Here we use a modified version of the SuperBayeS code to perform a Bayesian analysis of the minimal UED scenario, in order to assess its detectability at accelerators and with DM experiments. We derive in particular the most probable range of mass and scattering cross sections off nucleons, keeping into account cosmological and electroweak precision constraints. The consequences for the detectability of the gamma_1 with direct and indirect experiments are dramatic. The spin-independent cross section probability distribution peaks at ~ 10^{-11} pb, i.e. below the sensitivity of ton-scale experiments. The spin-dependent cross-section drives the predicted neutrino flux from the center of the Sun below the reach of present and upcoming experiments. The only strategy that remains open appears to be direct detection with ton-scale experiments sensitive to spin-dependent cross-sections. On the other hand, the LHC with 1 1/fb of data should be able to probe the current best-fit UED parameters.

We know that dark matter constitutes 85% of all the matter in the Universe, but we do not know of what it is made. Amongst the many Dark Matter candidates proposed, WIMPs (weakly interacting massive particles) occupy a special place, as they arise naturally from well motivated extensions of the standard model of particle physics. With the advent of the Large Hadron Collider at CERN, and a new generation of astroparticle experiments, the moment of truth has come for WIMPs: either we will discover them in the next five to ten years, or we will witness the inevitable decline of WIMP paradigm.

We study the prospects for detecting the annihilation products of Dark Matter [DM] in the framework of the two highest-resolution numerical simulations currently available, i.e. {\it Via Lactea II} and {\it Aquarius}. We propose a strategy to determine the shape and size of the region around the Galactic center that maximizes the probability of observing a DM signal, and we show that although the predicted flux can differ by a factor of 10 for a given DM candidate in the two simulation setups, the search strategy remains actually unchanged, since it relies on the angular profile of the annihilation flux, not on its normalization. We present mock gamma-ray maps that keep into account the diffuse emission produced by unresolved halos in the Galaxy, and we estimate that in an optimistic DM scenario a few individual clumps can be resolved above the background with the Fermi-LAT. Finally we calculate the energy-dependent boost factors for positrons and antiprotons, and show that they are always of $\cal O$(1), and therefore they cannot lead to the large enhancements of the antimatter fluxes required to explain the recent PAMELA, ATIC, Fermi and HESS data.

Dark matter (DM) is currently searched for with a variety of detection strategies. Accelerator searches are particularly promising, but even if Weakly Interacting Massive Particles (WIMPs) are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the DM in the Universe. We show that a significantly better reconstruction of the DM properties can be obtained with a combined analysis of LHC and direct detection (DD) data, by making a simple Ansatz on the WIMP local density, i.e. by assuming that the local density scales with the cosmological relic abundance. We demonstrate this method in an explicit example in the context of a 24-parameter supersymmetric model, with a neutralino LSP in the stau co-annihilation region. Our results show that future ton-scale DD experiments will allow to break degeneracies in the SUSY parameter space and achieve a significantly better reconstruction of the neutralino composition and its relic density than with LHC data alone.

Indirect Dark Matter searches are based on the observation of secondary particles produced by the annihilation or decay of Dark Matter. Among them, gamma-rays are perhaps the most promising messengers, as they do not suffer deflection or absorption on Galactic scales, so their observation would directly reveal the position and the energy spectrum of the emitting source. Here, we study the detailed gamma-ray energy spectrum of Kaluza--Klein Dark Matter in a theory with 5 Universal Extra Dimensions. We focus in particular on the two body annihilation of Dark Matter particles into a photon and another particle, which produces monochromatic photons, resulting in a line in the energy spectrum of gamma rays. Previous calculations in the context of the five dimensional UED model have computed the line signal from annihilations into \gamma \gamma, but we extend these results to include \gamma Z and \gamma H final states. We find that these spectral lines are subdominant compared to the predicted \gamma \gamma signal, but they would be important as follow-up signals in the event of the observation of the \gamma \gamma line, in order to distinguish the 5d UED model from other theoretical scenarios.

We study the effect of dark matter (DM) particles in the Sun, focusing in particular on the possible reduction of the solar neutrinos flux due to the energy carried away by DM particles from the innermost regions of the Sun, and to the consequent reduction of the temperature of the solar core. We find that in the very low-mass range between 4 and 10 GeV, recently advocated to explain the findings of the DAMA and CoGent experiments, the effects on neutrino fluxes are detectable only for DM models with very small, or vanishing, self-annihilation cross section, such as the so-called asymmetric DM models, and we study the combination of DM masses and Spin Dependent cross sections which can be excluded with current solar neutrino data. Finally, we revisit the recent claim that DM models with large self-interacting cross sections can lead to a modification of the position of the convective zone, alleviating or solving the solar composition problem. We show that when the `geometric' upper limit on the capture rate is correctly taken into account, the effects of DM are reduced by orders of magnitude, and the position of the convective zone remains unchanged.

A convincing identification of dark matter (DM) particles can probably be achieved only through a combined analysis of different detections strategies, which provides an effective way of removing degeneracies in the parameter space of DM models. In practice, however, this program is made complicated by the fact that different strategies depend on different physical quantities, or on the same quantities but in a different way, making the treatment of systematic errors rather tricky. We discuss here the uncertainties on the recoil rate in direct detection experiments and on the muon rate induced by neutrinos from dark matter annihilations in the Sun, and we show that, contrarily to the local DM density or overall cross section scale, irreducible astrophysical uncertainties affect the two rates in a different fashion, therefore limiting our ability to reconstruct the parameters of the dark matter particle. By varying within their respective errors astrophysical parameters such as the escape velocity and the velocity dispersion of dark matter particles, we show that the uncertainty on the relative strength of the neutrino and direct-detection signal is as large as a factor of two for typical values of the parameters, but can be even larger in some circumstances.

A precise determination of the local dark matter density and an accurate control over the corresponding uncertainties are of paramount importance for Dark Matter (DM) searches. Using very recent high-resolution numerical simulations of a Milky Way like object, we study the systematic uncertainties that affect the determination of the local dark matter density based on dynamical measurements in the Galaxy. In particular, extracting from the simulation with baryons the orientation of the Galactic stellar disk with respect to the DM distribution, we study the DM density for an observer located at $\sim$8 kpc from the Galactic center {\it on the stellar disk}, $\rho_0$. This quantity is found to be always larger than the average density in a spherical shell of same radius $\bar{\rho}_0$, which is the quantity inferred from dynamical measurements in the Galaxy, and to vary in the range $\rho_0/\bar{\rho}_0=1.01-1.41$. This suggests that the actual dark matter density in the solar neighbourhood is on average 21\% larger than the value inferred from most dynamical measurements, and that the associated systematic errors are larger than the statistical errors recently discussed in the literature.

Recent measurements of cosmic ray electrons and positrons by PAMELA, ATIC, Fermi and HESS have revealed interesting excesses and features in the GeV-TeV range. Many possible explanations have been suggested, invoking one or more nearby primary sources such as pulsars and supernova remnants, or dark matter. Based on the output of the TANGO in PARIS --Testing Astroparticle with the New GeV/TeV Observations in Positrons And electRons : Identifying the Sources-- workshop held in Paris in May 2009, we review here the latest experimental results and we discuss some virtues and drawbacks of the many theoretical interpretations proposed so far.

The rise in the energy spectrum of the positron ratio, observed by the PAMELA satellite above 10 GeV, and other cosmic ray measurements, have been interpreted as a possible signature of Dark Matter annihilation in the Galaxy. However, the large number of free parameters, and the large astrophysical uncertainties, make it difficult to do conclusive statements about the viability of this scenario. Here, we perform a multi-wavelength, multi-messenger analysis, that combines in a consistent way the constraints arising from different astrophysical observations. We show that if standard assumptions are made for the distribution of Dark Matter (we build models on the recent Via Lactea II and Aquarius simulations) and the propagation of cosmic rays, current DM models cannot explain the observed positron flux without exceeding the observed fluxes of antiprotons or gamma-ray and radio photons. To visualize the multi-messenger constraints, we introduce "star plots", a graphical method that allows to show in the same plot theoretical predictions and observational constraints for different messengers and wavelengths.

The possibility of explaining the positron and electron excess recently found by the PAMELA and ATIC collaborations in terms of dark matter (DM) annihilation has attracted considerable attention. Models surviving bounds from, e.g, antiproton production generally fall into two classes, where either DM annihilates directly with a large branching fraction into light leptons, or, as in the recent models of Arkani-Hamed et al., and of Nomura and Thaler, the annihilation gives low-mass (pseudo)scalars or vectors $\phi$ which then decay into $\mu^+\mu^-$ or $e^+e^-$. While the constraints on the first kind of models have recently been treated by several authors, we study here specifically models of the second type which rely on an efficient Sommerfeld enhancement in order to obtain the necessary boost in the annihilation cross section. We compute the photon flux generated by QED radiative corrections to the decay of $\phi$ and show that this indeed gives a rather spectacular broad peak in $E^2d\sigma/dE$, that for these extreme values of the cross section violate gamma-ray observations of the Galactic center for DM density profiles steeper than that of Navarro, Frenk and White. The most stringent constraint comes from the comparison of the predicted synchrotron radiation in the central part of the Galaxy with radio observations of Sgr A*. For the most commonly adopted DM profiles, the models that provide a good fit to the PAMELA and ATIC data are ruled out, unless there are physical processes that boost the local anti-matter fluxes more than one order of magnitude, while not affecting the gamma-ray or radio fluxes.

The formation and evolution of Black Holes inevitably affects the distribution of dark and baryonic matter in the neighborhood of the Black Hole. These effects may be particularly relevant around Supermassive and Intermediate Mass Black Holes (IMBHs), the formation of which can lead to large Dark Matter overdensities, called {\em spikes} and {\em mini-spikes} respectively. Despite being larger and more dense, spikes evolve at the very centers of galactic halos, in regions where numerous dynamical effects tend to destroy them. Mini-spikes may be more likely to survive, and they have been proposed as worthwhile targets for indirect Dark Matter searches. We review here the formation scenarios and the prospects for detection of mini-spikes, and we present new estimates for the abundances of mini-spikes to illustrate the sensitivity of such predictions to cosmological parameters and uncertainties regarding the astrophysics of Black Hole formation at high redshift. We also connect the IMBHs scenario to the recent measurements of cosmic-ray electron and positron spectra by the PAMELA, ATIC, H.E.S.S., and Fermi collaborations.

The injection of secondary particles produced by Dark Matter (DM) annihilation around redshift 1000 would inevitably affect the process of recombination, leaving an imprint on Cosmic Microwave Background (CMB) anisotropies and polarization. We show that the most recent CMB measurements provided by the WMAP satellite mission place interesting constraints on DM self-annihilation rates, especially for models that exhibit a large Sommerfeld enhancement of the annihilation cross-section, as recently proposed to fit the PAMELA and ATIC results. Furthermore, we argue that upcoming CMB experiments such as Planck, will improve the constraints by at least one order of magnitude, thus providing a sensitive probe of the properties of DM particles.

The spectrum of photons arising from WIMP annihilation carries a detailed imprint of the structure of the dark sector. In particular, loop-level annihilations into a photon and another boson can in principle lead to a series of lines (a WIMP forest) at energies up to the WIMP mass. A specific model which illustrates this feature nicely is a theory of two universal extra dimensions compactified on a chiral square. Aside from the continuum emission, which is a generic prediction of most dark matter candidates, we find a "forest" of prominent annihilation lines that, after convolution with the angular resolution of current experiments, leads to a distinctive (2-bump plus continuum) spectrum, which may be visible in the near future with the Fermi Gamma-Ray Space Telescope (formerly known as GLAST).

PAMELA and ATIC recently reported an excess in e+e- cosmic rays. We show that if it is due to Dark Matter annihilations, the associated gamma-ray flux and the synchrotron emission produced by e+e- in the galactic magnetic field violate HESS and radio observations of the galactic center and HESS observations of dwarf Spheroidals, unless the DM density profile is significantly less steep than the benchmark NFW and Einasto profiles.

We study the flux and the angular power spectrum of gamma-rays produced by Dark Matter (DM) annihilations in the Milky Way (MW) and in extra-galactic halos. The annihilation signal receives contributions from: a) the smooth MW halo, b) resolved and unresolved substructures in the MW, c) external DM halos at all redshifts, including d) their substructures. Adopting a self-consistent description of local and extra-galactic substructures, we show that the annihilation flux from substructures in the MW dominates over all the other components for angles larger than O(1) degrees from the Galactic Center, unless an extreme prescription is adopted for the substructures concentration. We also compute the angular power spectrum of gamma-ray anisotropies and find that, for an optimistic choice of the particle physics parameters, an interesting signature of DM annihilations could soon be discovered by the Fermi LAT satellite at low multipoles, l<100, where the dominant contribution comes from MW substructures with mass M>10^4 solar masses. For the substructures models we have adopted, we find that the contribution of extra-galactic annihilations is instead negligible at all scales.

This is a report on the findings of the dark matter science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses the prospects for dark matter detection with future gamma-ray experiments, and the complementarity of gamma-ray measurements with other indirect, direct or accelerator-based searches. We conclude that any comprehensive search for dark matter should include gamma-ray observations, both to identify the dark matter particle (through the charac- teristics of the gamma-ray spectrum) and to measure the distribution of dark matter in galactic halos.

Dark matter (DM) annihilation could in principle contribute to the diffuse cosmic gamma-ray back- ground (CGB). While with standard assumptions for cosmological and particle physics parameters this contribution is expected to be rather small, a number of processes could boost it, including a larger-than-expected DM annihilation cross-section, or the occurance of DM substructures such as DM mini-spikes around intermediate-mass black holes. We show that angular correlations of the CGB provide a tool to disentangle the signal induced by DM annihilation in mini-spikes from a conventional astrophysical component. Treating blazars as a known background, we study the prospects for detecting DM annihilations with the Fermi Gamma-Ray Space Telescope for different choices of DM mass and annihilation channels.

We investigate the sensitivity of the Gamma-ray Large Area Space Telescope (GLAST) to indirectly detect weakly interacting massive particles (WIMPs) through the $\gamma$-ray signal that their pair annihilation produces. WIMPs are among the favorite candidates to explain the compelling evidence that about 80% of the mass in the Universe is non-baryonic dark matter (DM). They are serendipitously motivated by various extensions of the standard model of particle physics such as Supersymmetry and Universal Extra Dimensions (UED). With its unprecedented sensitivity and its very large energy range (20 MeV to more than 300 GeV) the main instrument on board the GLAST satellite, the Large Area Telescope (LAT), will open a new window of discovery. As our estimates show, the LAT will be able to detect an indirect DM signature for a large class of WIMP models given a cuspy profile for the DM distribution. Using the current state of the art Monte Carlo and event reconstruction software developed within the LAT collaboration, we present preliminary sensitivity studies for several possible sources inside and outside the Galaxy. We also discuss the potential of the LAT to detect UED via the electron/positron channel. Diffuse background modeling and other background issues that will be important in setting limits or seeing a signal are presented.

The H.E.S.S. array of Cherenkov telescopes has performed, from 2004 to 2007, a survey of the inner Galactic plane at photon energies above 100 GeV. About 400 hours of data have been accumulated in the region between -30 and +60 degrees in Galactic longitude, and between -3 and +3 degrees in Galactic latitude. Assuming that dark matter is composed of Weakly Interacting Massive Particles, we calculate here the H.E.S.S. sensitivity map for dark matter annihilations, and derive the first experimental constraints on the ''mini-spikes'' scenario, in which a gamma-ray signal arises from dark matter annihilation around Intermediate Mass Black Holes. The data exclude the proposed scenario at a 90% confidence level for dark matter particles with velocity-weighted annihilation cross section sigma v above 10^28 cm3s^-1 and mass between 800 GeV and 10 TeV.

We study the impact of the capture and annihilation of Weakly Interacting Massive Particles (WIMPs) on the evolution of Pop III stars. With a suitable modification of the Geneva stellar evolution code, we study the evolution of 20 and 200 M$_\odot$ stars in Dark Matter Haloes with densities between 10$^{8}$ and $10^{11}$ GeV/cm$^3$ during the core H-burning phase, and, for selected cases, until the end of the core He-burning phase. We find that for WIMP densities higher than 5.3 $10^{10}(\sigma^{SD}_p/10^{-38} {cm}^2)^{-1}$ GeV cm$^{-3}$ the core H-burning lifetime of $20 M_{\odot}$ and $200 M_{\odot}$ stars exceeds the age of the Universe, and stars are sustained only by WIMP annihilations. We determine the observational properties of these `frozen` objects and show that they can be searched for in the local Universe thanks to their anomalous mass-radius relation, which should allow unambiguous discrimination from normal stars.

Galaxy formation requires a process that continually heats gas and quenches star formation in order to reproduce the observed shape of the luminosity function of bright galaxies. To accomplish this, current models invoke heating from supernovae, and energy injection from active galactic nuclei. However, observations of radio-loud active galactic nuclei suggest that their feedback is likely to not be as efficient as required, signaling the need for additional heating processes. We propose the self-annihilation of weakly interacting massive particles that constitute dark matter as a steady source of heating. In this paper, we explore the circumstances under which this process may provide the required energy input. To do so, dark matter annihilations are incorporated into a galaxy formation model within the Millennium cosmological simulation. Energy input from self-annihilation can compensate for all the required gas cooling and reproduce the observed galaxy luminosity function only for what appear to be extreme values of the relevant key parameters. The key parameters are: the slope of the inner density profile of dark matter haloes and the outer spike radius. The inner density profile needs to be steepened to slopes of -1.5 or more and the outer spike radius needs to extend to a few tens of parsecs on galaxy scales and a kpc or so on cluster scales. If neutralinos or any thermal relic WIMP with s-wave annihilation constitute dark matter, their self-annihilation is inevitable and could provide enough power to modulate galaxy formation. Energy from self-annihilating WIMPs could be yet another piece of the feedback puzzle along with supernovae and active galactic nuclei.

An extraordinarily rich zoo of non-baryonic Dark Matter candidates has been proposed over the last three decades. Here we present a 10-point test that a new particle has to pass, in order to be considered a viable DM candidate: I.) Does it match the appropriate relic density? II.) Is it {\it cold}? III.) Is it neutral? IV.) Is it consistent with BBN? V.) Does it leave stellar evolution unchanged? VI.) Is it compatible with constraints on self-interactions? VII.) Is it consistent with {\it direct} DM searches? VIII.) Is it compatible with gamma-ray constraints? IX.) Is it compatible with other astrophysical bounds? X.) Can it be probed experimentally?

Upcoming $\gamma$-ray satellites will search for Dark Matter annihilations in Milky Way substructures (or 'clumps'). The prospects for detecting these objects strongly depend on the assumptions made on the distribution of Dark Matter in substructures, and on the distribution of substructures in the Milky Way halo. By adopting simplified, yet rather extreme, prescriptions for these quantities, we compute the number of sources that can be detected with upcoming experiments such as GLAST, and show that, for the most optimistic particle physics setup ($m_\chi=40$ GeV and annihilation cross section $\sigma v = 3 \times 10^{-26}$ cm$^3$ s$^{-1}$), the result ranges from zero to $\sim$ hundred sources, all with mass above $10^{5}M\odot$. However, for a fiducial DM candidate with mass $m_\chi=100$ GeV and $\sigma v = 10^{-26}$ cm$^3$ s$^{-1}$, at most a handful of large mass substructures can be detected at $5 \sigma$, with a 1-year exposure time, by a GLAST-like experiment. Scenarios where micro-clumps (i.e. clumps with mass as small as $10^{-6}M\odot$) can be detected are severely constrained by the diffuse $\gamma$-ray background detected by EGRET.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

We review the consequences of the growth and evolution of Black Holes on the distribution of stars and Dark Matter (DM) around them. We focus in particular on Supermassive and Intermediate Mass Black Holes, and discuss under what circumstances they can lead to significant overdensities in the surrounding distribution of DM, thus effectively acting as DM annihilation boosters.

We discuss the consequences of the accretion of dark matter (DM) particles on compact stars such as white dwarfs and neutron stars. We show that in large regions of the DM parameter space, these objects are sensitive probes of the presence of DM and can be used to set constraints both on the DM density and on the physical properties of DM particles.

We review the current status of accelerator, direct and indirect Dark Matter (DM) searches, focusing on the complementarity of different techniques and on the prospects for discovery. After taking a census of present and upcoming DM-related experiments, we review the motivations to go beyond an "accelerator-only" approach, and highlight the benefits of multidisciplinarity in the quest for DM.

We point out a new and largely model-independent constraint on the dark matter scattering cross section with nucleons, applying when this quantity is larger than for typical weakly interacting dark matter candidates. When the dark matter capture rate in Earth is efficient, the rate of energy deposition by dark matter self-annihilation products would grossly exceed the measured heat flow of Earth. This improves the spin-independent cross section constraints by many orders of magnitude, and closes the window between astrophysical constraints (at very large cross sections) and underground detector constraints (at small cross sections). In the applicable mass range, from about 1 to about 10^{10} GeV, the scattering cross section of dark matter with nucleons is then bounded from above by the latter constraints, and hence must be truly weak, as usually assumed.

We study the prospects for detecting gamma-rays from decaying Dark Matter (DM), focusing in particular on gravitino DM in R-parity breaking vacua. Given the substantially different angular distribution of the predicted gamma-ray signal with respect to the case of annihilating DM, and the relatively poor (of order 0.1$^\circ$) angular resolution of gamma-ray detectors, the best strategy for detection is in this case to look for an exotic contribution to the gamma-ray flux at high galactic latitudes, where the decaying DM contribution would resemble an astrophysical extra-galactic component, similar to the one inferred by EGRET observations. Upcoming experiments such as GLAST and AMS-02 may identify this exotic contribution and discriminate it from astrophysical sources, or place significant constraints on the mass and lifetime of DM particles.

We study the impact of astrophysical processes on the gamma-ray background produced by the annihilation of dark matter particles in cosmological halos, with particular attention to the consequences of the formation of supermassive black holes. In scenarios where these objects form adiabatically from the accretion of matter on small seeds, dark matter is first compressed into very dense ``spikes'', then its density progressively decreases due to annihilations and scattering off of stellar cusps. With respect to previous analyses, based on non-evolving halos, the predicted annihilation signal is higher and significantly distorted at low energies, reflecting the large contribution to the total flux from unevolved spikes at high redshifts. The peculiar spectral feature arising from the specific redshift distribution of the signal, would discriminate the proposed scenario from more conventional astrophysical explanations. We discuss how this affects the prospects for detection and demonstrate that the gamma-ray background from DM annihilations might be detectable even in absence of a signal from the Galactic center.

The annihilation of dark matter (DM) in the Galaxy could produce specific imprints on the spectra of antimatter species in Galactic cosmic rays, which could be detected by upcoming experiments such as PAMELA and AMS02. Recent studies show that the presence of substructures can enhance the annihilation signal by a "boost factor" that not only depends on energy, but that is intrinsically a statistical property of the distribution of DM substructures inside the Milky Way. We investigate a scenario in which substructures consist of $\sim 100$ "mini-spikes" around intermediate-mass black holes. Focusing on primary positrons and antiprotons, we find large boost factors, up to a few thousand, that exhibit a large variance at high energy in the case of positrons and at low energy in the case of antiprotons. As a consequence, an estimate of the DM particle mass based on the observed cut-off in the positron spectrum could lead to a substantial underestimate of its actual value.

We study the prospects for detecting Weakly Interacting Massive Particles (WIMPs), in a number of phenomenological scenarios, with a detector composed of a target simultaneously sensitive to both spin-dependent and spin-independent couplings, as is the case of COUPP (Chicagoland Observatory for Underground Particle Physics). First, we show that sensitivity to both couplings optimizes chances of initial WIMP detection. Second, we demonstrate that in case of detection, comparison of the signal on two complementary targets, such as in COUPP CF3I and C4F10 bubble chambers, allows a significantly more precise determination of the dark matter axial and scalar couplings. This strategy would provide crucial information on the nature of the WIMPs, and possibly allow discrimination between neutralino and Kaluza-Klein dark matter.

The existence of a population of wandering Intermediate Mass Black Holes (IMBHs) is a generic prediction of scenarios that seek to explain the formation of Supermassive Black Holes in terms of growth from massive seeds. The growth of IMBHs may lead to the formation of DM overdensities called "mini-spikes", recently proposed as ideal targets for indirect DM searches. Current ground-based gamma-ray experiments, however, cannot search for these objects due to their limited field of view, and it might be challenging to discriminate mini-spikes in the Milky Way from the many astrophysical sources that GLAST is expected to observe. We show here that gamma-ray experiments can effectively search for IMBHs in the nearby Andromeda galaxy (also known as M31), where mini-spikes would appear as a distribution of point-sources, isotropically distributed in a \thickapprox 3^{\circ} circle around the galactic center. For a neutralino-like DM candidate with a mass m_{\chi}=150 GeV, up to 20 sources would be detected with GLAST (at 5\sigma, in 2 months). With Air Cherenkov Telescopes such as MAGIC and VERITAS, up to 10 sources might be detected, provided that the mass of neutralino is in the TeV range or above.

We study the prospects for detecting gamma-rays from point sources of Dark Matter annihilation with the space satellite GLAST. We simulate the instrument response to the gamma-ray spectrum arising from the annihilation of common Dark Matter candidates, and derive full-sky sensitivity maps for the {\it detection} of point sources and for the {\it identification} of the Dark Matter (as opposed to astrophysical) origin of the gamma-ray emission. These maps represent a powerful tool to assess the detectability of point sources, i.e. sources with angular size smaller than the angular resolution of GLAST, ~ 0.1 degrees, in {\it any} DM scenario. As an example, we apply the obtained results to the so-called 'mini-spikes' scenario, where the annihilation signal originates from large Dark Matter overdensities around Intermediate Mass Black Holes. We find that if these objects exist in the Galaxy, not only GLAST should be able to detect them over a timescale as short as 2 months, but in many cases it should be possible to determine with good accuracy the mass of the annihilating Dark Matter particles, while null searches would place stringent constraints on this scenario.

We show that the presence of a rho~1/r^{3/2} dark matter overdensity can be robustly predicted at the center of any galaxy old enough to have grown a power-law density cusp in the stars via the Bahcall-Wolf mechanism. Using both Fokker-Planck and direct N-body integrations, we demonstrate collisional generation of these dark matter "crests" (Collisionally REgenerated STtructures) even in the extreme case that the density of both stars and dark matter were previously lowered by slingshot ejection from a binary supermassive black hole. The time scale for collisional growth of the crest is approximately the two-body relaxation time as defined by the stars, which is < 10 Gyr at the centers of stellar spheroids with luminosities comparable to that of the Milky Way bulge or less. The presence of crests can robustly be predicted in such galaxies, unlike the steeper enhancements, called "spikes," produced by the adiabatic growth of black holes. We discuss special cases where the prospects for detecting dark matter annihilations from the centers of galaxy haloes are significantly affected by the formation of crests.

We review the status of indirect Dark Matter searches, focusing in particular on the connection with gamma-ray Astrophysics, and on the prospects for detection with the upcoming space telescope GLAST and Air Cherenkov Telescopes such as CANGAROO, HESS, MAGIC and VERITAS. After a brief introduction where we review the fundamental motivations for indirect searches, we tackle the question of whether it is possible to obtain strong enough evidence from astrophysical observations, to claim discovery of Dark Matter particles. To this purpose, we discuss some recent conflicting claims that have generated some confusion in the field, and present new strategies that may provide the long-awaited smoking-gun for Dark Matter.

Current strategies of indirect Dark Matter detection with neutrino telescopes are based on the search for high-energy neutrinos from the Solar core or from the center of the Earth. Here, we propose a new strategy based on the detection of neutrinos from Dark Matter annihilations in 'mini-spikes' around Intermediate Mass Black Holes. Neutrino fluxes, in this case, depend on the annihilation cross-section of Dark Matter particles, whereas solar and terrestrial fluxes are sensitive to the scattering cross-section off nucleons, a circumstance that makes the proposed search complementary to the existing ones. We discuss the prospects for detection with upcoming under-water and under-ice experiments such as ANTARES and IceCube, and show that several, up to many, sources could be detected with both experiments. A kilometer scale telescope in the Mediterranean appears to be ideally suited for the proposed search.

A recent observation of the 511 keV electron-positron annihilation line from the Galactic bulge has prompted a debate on the origin of the galactic positrons responsible for this emission. Assuming equilibrium between annihilation and injection of positrons in the Galaxy, we investigate the possibility that positrons were produced by past gamma ray bursts (GRBs). We compare the positron annihilation rate inferred from the observed 511 keV line with the predicted injection rate due to electron-positron pairs leaking out of individual GRBs and those pairs produced by GRB photons in the dense molecular clouds at the Galactic center. We conclude that the proposed scenario is marginally consistent with observations, and can reproduce the observed normalization of the annihilation line only for higher-than-average values of the GRB rate in the Galaxy and/or the molecular clouds optical depth for pair production.

We study the statistical properties of magnification perturbations by substructures in strong lensed systems using linear perturbation theory and an analytical substructure model including tidal truncation and a continuous substructure mass spectrum. We demonstrate that magnification perturbations are dominated by perturbers found within a tidal radius of an image, and that sizable magnification perturbations may arise from small, coherent contributions from several substructures within the lens halo. We find that the root-mean-square (rms) fluctuation of the magnification perturbation is 10% to 20% and both the average and rms perturbations are sensitive to the mass spectrum and density profile of the perturbers. Interestingly, we find that relative to a smooth model of the same mass, the average magnification in clumpy models is lower (higher) than that in smooth models for positive (negative) parity images. This is opposite from what is observed if one assumes that the image magnification predicted by the best-fit smooth model of a lens is a good proxy for what the observed magnification would have been if substructures were absent. While it is possible for this discrepancy to be resolved via nonlinear perturbers, we argue that a more likely explanation is that the assumption that the best-fit lens model is a good proxy for the magnification in the absence of substructure is not correct. We conclude that a better theoretical understanding of the predicted statistical properties of magnification perturbations by CDM substructure is needed in order to affirm that CDM substructures have been unambiguously detected.

The prospects of indirect detection of dark matter at the galactic center depend sensitively on the mass profile within the inner parsec. We calculate the distribution of dark matter on sub-parsec scales by integrating the time-dependent Fokker-Planck equation, including the effects of self-annihilations, scattering of dark matter particles by stars, and capture in the supermassive black hole. We consider a variety of initial dark matter distributions, including models with very high densities ("spikes") near the black hole, and models with "adiabatic compression" of the baryons. The annihilation signal after 10 Gyr is found to be substantially reduced from its initial value, but in dark matter models with an initial spike, order-of-magnitude enhancements can persist compared with the rate in spike-free models, with important implications for indirect dark matter searches with GLAST and Air Cherenkov Telescopes like HESS and CANGAROO.

We study the prospects for detecting gamma-rays from Dark Matter (DM) annihilations in enhancements of the DM density (mini-spikes) around intermediate-mass black holes with masses in the range $10^2 \lsim M / \msun \lsim 10^6$. Focusing on two different IMBH formation scenarios, we show that, for typical values of mass and cross section of common DM candidates, mini-spikes, produced by the adiabatic growth of DM around pregalactic IMBHs, would be bright sources of gamma-rays, which could be easily detected with large field-of-view gamma-ray experiments such as GLAST, and further studied with smaller field-of-view, larger-area experiments like Air Cherenkov Telescopes CANGAROO, HESS, MAGIC and VERITAS. The detection of many gamma-ray sources not associated with a luminous component of the Local Group, and with identical cut-offs in their energy spectra at the mass of the DM particle, would provide a potential smoking-gun signature of DM annihilations and shed new light on the nature of intermediate and supermassive Black Holes.

Non-baryonic, or "dark," matter is believed to be a major component of the total mass budget of the universe. We review the candidates for particle dark matter and discuss the prospects for direct detection (via interaction of dark matter particles with laboratory detectors) and indirect detection (via observations of the products of dark matter self-annihilations), focusing in particular on the Galactic center, which is among the most promising targets for indirect detection studies. The gravitational potential at the Galactic center is dominated by stars and by the supermassive black hole, and the dark matter distribution is expected to evolve on sub-parsec scales due to interaction with these components. We discuss the dominant interaction mechanisms and show how they can be used to rule out certain extreme models for the dark matter distribution, thus increasing the information that can be gleaned from indirect detection searches.

The Galactic positrons, as observed by their annihilation gamma-ray line at 0.511 MeV, are difficult to account for with astrophysical sources. It has been proposed that they are produced instead by dark matter annihilation or decay in the inner Galactic halo. To avoid other constraints, these processes are required to occur "invisibly," such that the eventual positron annihilation is the only detectable signal. However, electromagnetic radiative corrections to these processes inevitably produce real gamma rays (``internal bremsstrahlung"); this emission violates COMPTEL and EGRET constraints unless the dark matter mass is less than about 20 MeV.

In this review article, we discuss the current status of particle dark matter, including experimental evidence and theoretical motivations. We discuss a wide array of candidates for particle dark matter, but focus on neutralinos in models of supersymmetry and Kaluza-Klein dark matter in models of universal extra dimensions. We devote much of our attention to direct and indirect detection techniques, the constraints placed by these experiments and the reach of future experimental efforts.

Direct and indirect dark matter searches can potentially lead to the discovery of Supersymmetry. In the specific context of heterotic orbifold compactification scenarios, it is even possible to get some clues on the SUSY breaking mechanism of the more fundamental underlying theory. In this paper, we investigate the prospects for indirect detection of neutralino dark matter via gamma-ray (continuum and line) and synchrotron radiation from the Galactic center, in the context of AMSB scenarios in an effective heterotic framework, which constitutes a consistent SUSY model between high energy theory and low energy phenomenology.

An intriguing feature of extra dimensions is the possible production of Kaluza--Klein gravitons by nucleon-nucleon bremsstrahlung in the course of core collapse of massive stars, as pointed out by Hannestad and Raffelt. In this event Kaluza--Klein gravitons are copiously produced and a significant fraction of them remains trapped around the newly born neutron stars. They slowly decay into 2 gamma rays, making neutron stars gamma-ray sources. In this letter, we strengthen considerably the limits on the radius of compactification of extra-dimensions for small number $n$ of them, or alternatively the fundamental scale of quantum gravity, considering the gamma-ray emission of the whole population of neutron stars sitting in the Galactic bulge, instead of the closest member of this category. For $n=1$ the constraint on the compactification radius is $R<400\mu$m, overlapping with the distance ($180\mu$m) at which Newtons law is directly measured. In addition, for n=1 and n=2, the fundamental energy scale of quantum gravity is far beyond the collider technology. These results imply that if $n\lesssim4$ and if strong gravity is around a TeV, the compactification topology is to be more complex than that of a torus.

We discuss the prospects for detection of high energy neutrinos from dark matter annihilation at the Galactic centre. Despite the large uncertainties associated with our poor knowledge of the distribution of dark matter in the innermost regions of the Galaxy, we determine an upper limit on the neutrino flux by requiring that the associated gamma-ray emission does not exceed the observed flux. We conclude that if dark matter is made of neutralinos, a neutrino flux from dark matter annihilations at the GC will not be observable by Antares. Conversely, the positive detection of such a flux would either require an alternative explanation, in terms of astrophysical processes, or the adoption of other \DM candidates, disfavouring the case for neutralinos.

One of the models recently proposed to explain the origin of the ultra high energy cosmic rays assumes that these particles may be accelerated by the electromotive force around presently inactive quasar remnants. We study predictions for large and small scale ultra high energy cosmic ray arrival direction anisotropies in a scenario where the particles are injected with a mono-energetic spectrum by a discrete distribution of such sources. We find that known quasar remnants are typically distributed too anisotropically to explain the isotropic ultra high energy cosmic ray flux except in the unrealistic case where extragalactic magnetic fields of ~0.1 microGauss extend over many Mpc.

The annihilation signal of particle dark matter can be strongly enhanced in over-dense regions such as close to the Galactic centre. We summarize some of our recent results on fluxes of gamma-rays, neutrinos and radio waves under different assumptions for the largely uncertain dark matter profile which close to the Galactic centre is strongly influenced by the super-massive black hole. We apply this to two particle dark matter scenarios, namely the case of neutralinos in supersymmetric scenarios and the lowest Kaluza-Klein excitation of the hyper-charge gauge boson in scenarios where the Standard Model fields propagate in one microscopic extra dimension.

We investigate prospects for indirect detection of Kaluza--Klein dark matter, focusing on the annihilation radiation of the first Kaluza--Klein excitation of the Hypercharge gauge boson $B^{(1)}$ in the Galactic halo, in particular we estimate neutrino, gamma-ray and synchrotron fluxes. Comparing the predicted fluxes with observational data we are able to constrain the $B^{(1)}$ mass (and therefore the compactification scale). The constraints depend on the specific model adopted for the dark matter density profile. For a NFW profile the analysis of synchrotron radiation puts a lower bound on the $B^{(1)}$ mass of the order of $\simeq 300$ GeV.

We study annihilation radiation of neutralinos in the Galactic centre, assuming the existence of a 'spike' in the dark matter density profile, due to adiabatic accretion onto the massive black hole lying at the Galactic centre. Under this assumption we find that it is possible to reproduce the observed SgrA* emission at radio and gamma-ray frequencies in a consistent scenario with a magnetic field close to the equipartition strength and with values of $\gamma$ (the density profile power law index) around 0.1.

We extend existing work on the propagation of ultra-high energy cosmic rays in extragalactic magnetic fields to a possible component of heavy nuclei, taking into account photodisintegration, pion production, and creation of e+e- pairs. We focus on the influence of the magnetic field on the spectrum and chemical composition of observed ultra-high energy cosmic rays. We apply our simulations to the scenarios proposed by Anchordoqui et al, in which Iron nuclei are accelerated in nearby starburst galaxies, and show that it is in marginal agreement with the data. We also show that it is highly unlikely to detect He nuclei from M87 at the highest energies observed $\sim3 10^{20} $eV as required for the scenario of Ahn et al. in which the highest energy cosmic rays originate from M87 and are deflected in a Parker spiral Galactic magnetic field.

If the Milky Way is populated by WIMPs as predicted by cosmological models of the large-scale structure of the universe and as motivated by SUSY, the capture of high-mass WIMPs by the Sun would affect the temperature, density and chemical composition of the solar core. We use the sound speed and the density profiles inferred from the helioseismic instruments on the Solar and Heliospheric Observatory (SOHO) to discuss the effect of WIMP accretion and annihilation on the evolution of the Sun. The WIMP transport of energy inside the Sun is not negligible for WIMPs with a mass smaller than 60 GeV and annihilating WIMPs with <sigma_a v>~ 10^{-27}cm^3/sec. WIMP-accreting models with WIMP masses smaller than 30 GeV are in conflict with the most recent seismic data. We combine our new constraints with the analysis of predicted neutrino fluxes from annihilating WIMPs in the solar core. Working in the framework of the Minimal Supersymmetric Standard Model and considering the neutralino as the best dark matter particle candidate, we find that supersymmetric models, consistent with solar seismic data and with recent measurements of dark matter relic density, lead to a measured muon flux on Earth in the range of 1 to 10^4 km^{-2} yr^{-1}, for neutralino masses between 30 and 400 GeV.

We briefly review the Thomas-Fermi statistical model of atoms in the classical non-relativistic formulation and in the generalised finite-nucleus relativistic formulation. We then discuss the classical generalisation of the model to finite temperatures in the non-relativistic approximation and present a new relativistic model at finite temperatures, investigating how to recover the existing theory in the limit of low temperatures. This work is intended to be a propedeutical study for the evaluation of equilibrium configurations of relativistic ``hot'' white dwarfs.

Annihilations of weakly interacting dark matter particles provide an important signature for the possibility of indirect detection of dark matter in galaxy halos. These self-annihilations can be greatly enhanced in the vicinity of a massive black hole. We show that the massive black hole present at the centre of our galaxy accretes dark matter particles, creating a region of very high particle density. Consequently the annihilation rate is considerably increased, with a large number of $e^+e^-$ pairs being produced either directly or by successive decays of mesons. We evaluate the synchrotron emission (and self-absorption) associated with the propagation of these particles through the galactic magnetic field, and are able to constrain the allowed values of masses and cross sections of dark matter particles.

A massive black hole is present at the centre of our galaxy and inevitably accretes dark matter particles, creating a region of very high particle density. The annihilation rate is enhanced with a large number of e+e- pairs produced either directly or by successive decays of mesons. We evaluate the synchrotron emission (and self-absorption) associated with the propagation of these particles through the galactic magnetic field and give constraints on the values of mass and cross section of the dark matter particles.

The Feynman-Metropolis-Teller treatment for compressed atoms is here reconsidered in the framework of the relativistic generalised Fermi-Thomas model, obtained by Ruffini et al. Physical properties of a zero temperature plasma is thus investigated and the resulting equation of state, which keeps into account quantum, relativistic and electromagnetic effects, is applied to the study of equilibrium configurations of relativistic White Dwarfs. It is shown that numerical evaluation of such configuration leads, for the same central density $\rho_c$, to smaller values of radius R and of mass M than in the classical works of Chandrasekhar and Salpeter, the deviations being most marked at the lowest densities (up to 30% from the Chandrasekhar model and 10% from the Salpeter one for $\rho_c \sim 10^6g/cm^3$, corresponding to $M \sim 0.2 M_{\odot}$). At high densities we considered the occurrence of inverse beta decays, whose effect is to introduce gravitational instability of the configurations. We consequently find the maximum mass of White Dwarfs, which, for an Oxygen and an Iron WD, is respectively $1.365 M_{\odot}$ and $1.063 M_{\odot}$.