We present $\texttt{galkin}$, a novel compilation of kinematic measurements tracing the rotation curve of our Galaxy, together with a tool to treat the data. The compilation is optimised to Galactocentric radii between 3 and 20 kpc and includes the kinematics of gas, stars and masers in a total of 2780 measurements carefully collected from almost four decades of literature. A simple, user-friendly tool is provided to select, treat and retrieve the data of all source references considered. This tool is especially designed to facilitate the use of kinematic data in dynamical studies of the Milky Way with various applications ranging from dark matter constraints to tests of modified gravity.

We introduce a new type of gamma-ray spectral feature, which we denominate gamma-ray triangle. This spectral feature arises in scenarios where dark matter self-annihilates via a chiral interaction into two Dirac fermions, which subsequently decay in flight into another fermion and a photon. The resulting photon spectrum resembles a sharp triangle and can be readily searched for in the gamma-ray sky. Using data from the Fermi-LAT and H.E.S.S. instruments, we find no evidence for such spectral feature and therefore set strong upper bounds on the corresponding annihilation cross section. A concrete realization of a scenario yielding gamma-ray triangles consists of an asymmetric dark matter model where the dark matter particle carries lepton number. We show explicitly that this class of models can lead to intense gamma-ray spectral features, potentially at the reach of upcoming gamma-ray telescopes, opening a new window to explore asymmetric dark matter through indirect searches.

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.

We compute the gamma-ray output of axion-mediated dark matter and derive the corresponding constraints set by recent data. In such scenarios the dark matter candidate is a Dirac fermion that pair-annihilates into axions and/or scalars. Provided that the axion decays (at least partly) into photons, these models naturally give rise to a box-shaped gamma-ray spectrum that may present two distinct phenomenological behaviours: a narrow box, resembling a line at half the dark matter mass, or a wide box, spanning an extensive energy range up to the dark matter mass. Remarkably, we find that in both cases a sizable gamma-ray flux is predicted for a thermal relic without fine-tuning the model parameters nor invoking boost factors. This large output is in line with recent Fermi-LAT observations towards the Galactic centre region and is on the verge of being excluded. We then make use of the Fermi-LAT and H.E.S.S. data to derive robust, model-independent upper limits on the dark matter annihilation cross section for the narrow and wide box scenarios. H.E.S.S. constraints, in particular, turn out to match the ones from Fermi-LAT at hundreds of GeV and extend to multi-TeV masses. Future Cherenkov telescopes will likely probe gamma-ray boxes from thermal dark matter relics in the whole multi-TeV range, a region hardly accessible to direct detection, collider searches and other indirect detection strategies.

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.

Mapping the dark matter distribution across our Galaxy represents a central challenge for the near future as a new generation of space-borne and ground-based astronomical surveys swiftly comes online. Here we present a synopsis of the present status of the field, reviewing briefly the baryonic content and the kinematics of the Milky Way and outlining the methods used to infer the dark matter component. The discussion then proceeds with some of the latest developments based on our own work. In particular, we present a new compilation of kinematic measurements tracing the rotation curve of the Galaxy and an exhaustive array of observation-based baryonic models setting the contribution of stellar bulge, stellar disc and gas to the total gravitational potential. The discrepancy between these two components is then quantified to derive the latest constraints on the dark matter distribution and on modified Newtonian dynamics. We shall end with an overview of future directions to improve our mapping of the dark matter distribution in the Milky Way.

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.

We present the results of a new, non-parametric method to reconstruct the Galactic dark matter profile directly from observations. Using the latest kinematic data to track the total gravitational potential and the observed distribution of stars and gas to set the baryonic component, we infer the dark matter contribution to the circular velocity across the Galaxy. The radial derivative of this dynamical contribution is then estimated to extract the dark matter profile. The innovative feature of our approach is that it makes no assumption on the functional form nor shape of the profile, thus allowing for a clean determination with no theoretical bias. We illustrate the power of the method by constraining the spherical dark matter profile between 2.5 and 25 kpc away from the Galactic centre. The results show that the proposed method, free of widely used assumptions, can already be applied to pinpoint the dark matter distribution in the Milky Way with competitive accuracy, and paves the way for future developments.

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.

Although proposed long ago, the search for internal bremsstrahlung signatures has only recently been made possible by the excellent energy resolution of ground-based and satellite-borne gamma-ray instruments. Here, we investigate thoroughly the current status of internal bremsstrahlung searches in light of the results of direct dark matter searches and in the framework of a minimal mass-degenerate scenario consisting of a Majorana dark matter particle that couples to a fermion and a scalar via a Yukawa coupling. The upper limits on the annihilation cross section set by Fermi-LAT and H.E.S.S. extend uninterrupted from tens of GeV up to tens of TeV and are rather insensitive to the mass degeneracy in the particle physics model. In contrast, direct searches are best in the moderate to low mass splitting regime, where XENON100 limits overshadow Fermi-LAT and H.E.S.S. up to TeV masses if dark matter couples to one of the light quarks. In our minimal scenario we examine carefully the prospects for GAMMA-400, CTA and XENON1T, all planned to come online in the near future, and find that: (a) CTA and XENON1T are fully complementary, with CTA most sensitive to multi-TeV masses and mass splittings around 10%, and XENON1T probing best small mass splittings up to TeV masses; and (b) current constraints from XENON100 already preclude the observation of any spectral feature with GAMMA-400 in spite of its impressive energy resolution, unless dark matter does not couple predominantly to light quarks. Finally, we point out that, unlike for direct searches, the possibility of detecting thermal relics in upcoming internal bremsstrahlung searches requires, depending on the concrete scenario, boost factors larger than 5-10.

The latest XENON100 data severely constrains dark matter elastic scattering off nuclei, leading to impressive upper limits on the spin-independent cross-section. The main goal of this paper is to stress that the same data set has also an excellent \emph{spin-dependent} sensitivity, which is of utmost importance in probing dark matter models. We show in particular that the constraints set by XENON100 on the spin-dependent neutron cross-section are by far the best at present, whereas the corresponding spin-dependent proton limits lag behind other direct detection results. The effect of nuclear uncertainties on the structure functions of xenon isotopes is analysed in detail and found to lessen the robustness of the constraints, especially for spin-dependent proton couplings. Notwithstanding, the spin-dependent neutron prospects for XENON1T and DARWIN are very encouraging. We apply our constraints to well-motivated dark matter models and demonstrate that in both mass-degenerate scenarios and the minimal supersymmetric standard model the spin-dependent neutron limits can actually override the spin-independent limits. This opens the possibility of probing additional unexplored regions of the dark matter parameter space with the next generation of ton-scale direct detection experiments.

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.

Over the last years both cosmic-ray antiproton measurements and direct dark matter searches have proved particularly effective in constraining the nature of dark matter candidates. The present work focusses on these two types of constraints in a minimal framework which features a Majorana fermion as the dark matter particle and a scalar that mediates the coupling to quarks. Considering a wide range of coupling schemes, we derive antiproton and direct detection constraints using the latest data and paying close attention to astrophysical and nuclear uncertainties. Both signals are strongly enhanced in the presence of degenerate dark matter and scalar masses, but we show that the effect is especially dramatic in direct detection. Accordingly, the latest direct detection limits take the lead over antiprotons. We find that antiproton and direct detection data set stringent lower limits on the mass splitting, reaching 19% at a 300 GeV dark matter mass for a unity coupling. Interestingly, these limits are orthogonal to ongoing collider searches at the Large Hadron Collider, making it feasible to close in on degenerate dark matter scenarios within the next years.

The observation of a sharp spectral feature in the gamma-ray sky would be one of the cleanest ways to identify dark matter and pinpoint its properties. Over the years a lot of attention has been paid to two specific features, namely gamma-ray lines and internal bremsstrahlung. Here, we explore a third class of spectral signatures, box-shaped gamma-ray spectra, that naturally arise in dark matter cascade annihilations or decays into intermediate particles that in turn decay into photons. Using Fermi-LAT data, we derive constraints on the dark matter parameter space for both annihilating and decaying dark matter, and show explicitly that our limits are competitive to strategies employing standard spectral features. More importantly, we find robust limits even in the case of non-degenerate dark matter and intermediate particle masses. This result is particularly relevant in constraining dark matter frameworks with gamma-ray data. We conclude by illustrating the power of box-shaped gamma-ray constraints on concrete particle physics scenarios.

We study the synthesis of lithium isotopes in the hot tori formed around stellar mass black holes by accretion of the companion star. We find that sizable amounts of both stable isotopes 6Li and 7Li can be produced, the exact figures varying with the characteristics of the torus and reaching as much as 1e-2 Msun for each isotope. This mass output is enough to contaminate the entire Galaxy at a level comparable with the original, pre-galactic amount of lithium and to overcome other sources such as cosmic-ray spallation or stellar nucleosynthesis.

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.

Direct searches for dark matter have prompted in recent years a great deal of excitement within the astroparticle physics community, but the compatibility between signal claims and null results of different experiments is far from being a settled issue. In this context, we study here the prospects for constraining the dark matter parameter space with the next generation of ton-scale detectors. Using realistic experimental capabilities for a wide range of targets (including fluorine, sodium, argon, germanium, iodine and xenon), the role of target complementarity is analysed in detail while including the impact of astrophysical uncertainties in a self-consistent manner. We show explicitly that a multi-target signal in future direct detection facilities can determine the sign of the ratio of scalar couplings $f_n/f_p$, but not its scale. This implies that the scalar-proton cross-section is left essentially unconstrained if the assumption $f_p\sim f_n$ is relaxed. Instead, we find that both the axial-proton cross-section and the ratio of axial couplings $a_n/a_p$ can be measured with fair accuracy if multi-ton instruments using sodium and iodine will eventually come online. Moreover, it turns out that future direct detection data can easily discriminate between elastic and inelastic scatterings. Finally, we argue that, with weak assumptions regarding the WIMP couplings and the astrophysics, only the dark matter mass and the inelastic parameter (i.e. mass splitting) may be inferred from the recoil spectra -- specifically, we anticipate an accuracy of tens of GeV (tens of keV) in the measurement of the dark matter mass (inelastic parameter).

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.

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.

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.

The Alpha Magnetic Spectrometer (AMS-02), which is scheduled to be deployed onboard the International Space Station later this year, will be capable of measuring the composition and spectra of GeV-TeV cosmic rays with unprecedented precision. In this paper, we study how the projected measurements from AMS-02 of stable secondary-to-primary and unstable ratios (such as boron-to-carbon and beryllium-10-to-beryllium-9) can constrain the models used to describe the propagation of cosmic rays throughout the Milky Way. We find that within the context of fairly simple propagation models, all of the model parameters can be determined with high precision from the projected AMS-02 data. Such measurements are less constraining in more complex scenarios, however, which allow for departures from a power-law form for the diffusion coefficient, for example, or for inhomogeneity or stochasticity in the distribution and chemical abundances of cosmic ray sources.

Alternative cosmologies, based on extensions of General Relativity, predict modified thermal histories in the Early Universe during the pre Big Bang Nucleosynthesis (BBN) era, epoch which is not directly constrained by cosmological observations. When the expansion rate is enhanced with respect to the standard case, thermal relics typically decouple with larger relic abundances. The correct value of the relic abundance is therefore obtained for larger annihilation cross--sections, as compared to standard cosmology. A direct consequence is that indirect detection rates are enhanced. Extending previous analyses of ours, we derive updated astrophysical bounds on the dark matter annihilation cross sections and use them to constrain alternative cosmologies in the pre--BBN era. We also determine the characteristics of these alternative cosmologies in order to provide the correct value of relic abundance for a thermal relic for the (large) annihilation cross--section required to explain the PAMELA results on the positron fraction, therefore providing a "cosmological boost" solution to the dark matter interpretation of the PAMELA data.

(Abridged) We study the predictions of various annihilating Dark Matter (DM) models in order to interpret the origin of non-thermal phenomena in galaxy clusters. We consider three neutralino DM models with light (9 GeV), intermediate (60 GeV) and high (500 GeV) mass. The secondary particles created by neutralino annihilation produce a multi-frequency Spectral Energy Distribution (SED), as well as heating of the intracluster gas, that are tested against the observations available for the Coma cluster. The DM produced SEDs are normalized to the Coma radio halo spectrum. We find that it is not possible to interpret all non-thermal phenomena observed in Coma in terms of DM annihilation. The DM model with 9 GeV mass produces too small power at all frequencies, while the DM model with 500 GeV produces a large excess power at all frequencies. The DM model with 60 GeV and $\tau^{\pm}$ composition is consistent with the HXR and gamma-ray data but fails to reproduce the EUV and soft X-ray data. The DM model with 60 GeV and $b{\bar b}$ composition is always below the observed fluxes. The radio halo spectrum of Coma is well fitted only in the $b{\bar b}$ or light and intermediate mass DM models. The heating produced by DM annihilation in the center of Coma is always larger than the intracluster gas cooling rate for an NFW DM density profile and it is substantially smaller than the cooling rate only for a cored DM density profile in DM model with 9 GeV. We conclude that the possibility of interpreting the origin of non-thermal phenomena in galaxy clusters with DM annihilation models requires a low neutralino mass and a cored DM density profile. If we then consider the multimessenger constraints to the neutralino annihilation cross-section, it turns out that such scenario would also be excluded unless we introduce a substantial boost factor due to the presence of DM substructures.

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.

A new method to reconstruct the 3-dimensional structure of extensive air showers, seen by fluorescence detectors, is proposed. The observation of the shower is done in 2-dimensional pixels, for consecutive time bins. Time corresponds to a third dimension. Assuming that the cosmic ray shower propagates as a plane wave front moving at the speed of light, a complex 3D volume in space can be associated to each measured charge (per pixel and time bin). The 3D description in space allows a simultaneous access to the longitudinal and lateral profiles of each shower. In the case that several eyes observe the same shower, the method gives a straight-forward combination of all the information. This method is in an early phase of development and is not used for the general reconstruction of the Auger data.