The vertical diffusive halo size of the Galaxy, $L$, is a key parameter for dark matter indirect searches. It can be better determined thanks to recent AMS-02 data. We set constraints on $L$ from Be/B and $^{10}$Be/Be data, and we performed a consistency check with positron data. We detail the dependence of Be/B and $^{10}$Be/Be on $L$ and forecast on which energy range better data would be helpful for future $L$ improvements. We used USINE v3.5 for the propagation of nuclei, and $e^+$ were calculated with the pinching method of Boudaud et al. (2017). The current AMS-02 Be/B ($\sim3\%$ precision) and ACE-CRIS $^{10}$Be/Be ($\sim 10\%$ precision) data bring similar and consistent constraints on $L$. The AMS-02 Be/B data alone constrain $L=5^{+3}_{-2}$~kpc at a 68\% confidence level (spanning different benchmark transport configurations), a range for which most models do not overproduce positrons. Future experiments need to deliver percent-level accuracy on $^{10}$Be/$^9$Be anywhere below 10 GV to further constrain $L$. Forthcoming AMS-02, HELIX, and PAMELA $^{10}$Be/$^9$Be results will further test and possibly tighten the limits derived here. Elemental ratios involving radioactive species with different lifetimes (e.g. Al/Mg and Cl/Ar) are also awaited to provide complementary and robuster constraints.

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

The AMS-02 experiment has ushered cosmic-ray physics into precision era. In a companion paper, we designed an improved method to calibrate propagation models on B/C data. Here we provide a robust prediction of the $\bar{p}$ flux, accounting for several sources of uncertainties and their correlations. Combined with a correlation matrix for the $\bar{p}$ data, we show that the latter are consistent with a secondary origin. This paper presents key elements relevant to the dark matter search in this channel, notably by pointing out the inherent difficulties in achieving predictions at the percent-level precision.

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

The distribution of dark matter (DM) subhalos in our Galaxy remains disputed, leading to varying $\gamma$-ray and $\nu$ flux predictions from their annihilation or decay. In this work, we study how, in the inner Galaxy, subhalo tidal disruption from the Galactic baryonic potential impacts these signals. Based on state-of-the art modeling of this effect from numerical simulations and semi-analytical results, updated subhalo spatial distributions are derived and included in the CLUMPY code. The latter is used to produce a thousand realizations of the $\gamma$-ray and $\nu$ sky. Compared to predictions based on DM only, we conclude a decrease of the flux of the brightest subhalo by a factor 2 to 7 for annihilating DM and no impact on decaying DM: the discovery prospects or limits subhalos can set on DM candidates are affected by the same factor. This study also provides probability density functions for the distance, mass, and angular distribution of the brightest subhalo, among which the mass may hint at its nature: it is most likely a dwarf spheroidal galaxy in the case of strong tidal effects from the baryonic potential, whereas it is lighter and possibly a dark halo for DM only or less pronounced tidal effects.

This article aims at establishing new benchmark scenarios for Galactic cosmic-ray propagation in the GV-TV rigidity range, based on fits to the AMS-02 B/C data with the USINE v3.5 propagation code. We employ a new fitting procedure, cautiously taking into account data systematic error correlations in different rigidity bins and considering Solar modulation potential and leading nuclear cross-section as nuisance parameters. We delineate specific low, intermediate, and high-rigidity ranges that can be related to both features in the data and peculiar microphysics mechanisms resulting in spectral breaks. We single out a scenario which yields excellent fits to the data and includes all the presumably relevant complexity, the BIG model. This model has two limiting regimes: (i) the SLIM model, a minimal diffusion-only setup, and (ii) the QUAINT model, a convection-reacceleration model where transport is tuned by non-relativistic effects. All models lead to robust predictions in the high-energy regime ($\gtrsim10$GV), i.e. independent of the propagation scenario: at $1\sigma$, the diffusion slope $\delta$ is $[0.43-0.53]$, whereas $K_{10}$, the diffusion coefficient at 10GV, is $[0.26-0.36]$kpc$^2$Myr$^{-1}$; we confirm the robustness of the high-energy break, with a typical value $\Delta_h\sim 0.2$. We also find a hint for a similar (reversed) feature at low rigidity around the B/C peak ($\sim 4$GV) which might be related to some effective damping scale in the magnetic turbulence.

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

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

Using cosmic-ray boron to carbon ratio (B/C) data recently released by the AMS-02 experiment, we find tantalizing indications ({\it decisive evidence}, in Bayesian terms) in favor of a diffusive origin for the broken power-law spectra found in protons ($p$) and helium nuclei (He). The result is robust with respect to currently estimated uncertainties in the cross sections, and in the presence of a small component of primary boron, expected because of spallation at the acceleration site. Reduced errors at high energy as well as further cosmic ray nuclei data (as absolute spectra of C, N, O, Li, Be) may definitively confirm this scenario.

MeV dark matter (DM) particles annihilating or decaying to electron-positron pairs cannot, in principle, be observed via local cosmic-ray (CR) measurements because of the shielding solar magnetic field. In this letter, we take advantage of spacecraft Voyager 1's capacity for detecting interstellar CRs since it crossed the heliopause in 2012. This opens up a new avenue to probe DM in the sub-GeV energy/mass range that we exploit here for the first time. From a complete description of the transport of electrons and positrons at low energy, we derive predictions for both the secondary astrophysical background and the pair production mechanisms relevant to DM annihilation or decay down to the MeV mass range. Interestingly, we show that reacceleration may push positrons up to energies larger than the DM particle mass. We combine the constraints from the Voyager and AMS-02 data to get novel limits covering a very extended DM particle mass range, from MeV to TeV. In the MeV mass range, our limits reach annihilation cross sections of order $\langle \sigma v\rangle \sim 10^{-28}{\rm cm^3/s}$. An interesting aspect is that these limits barely depend on the details of cosmic-ray propagation in the weak reacceleration case, a configuration which seems to be favored by the most recent boron-to-carbon ($B/C$) data. Though extracted from a completely different and new probe, these bounds have a strength similar to those obtained with the cosmic microwave background --- they are even more stringent for $p$-wave annihilation.

The interaction properties of cold dark matter (CDM) particle candidates, such as those of weakly interacting massive particles (WIMPs), generically lead to the structuring of dark matter on scales much smaller than typical galaxies, potentially down to $\sim 10^{-10}M_\odot$. This clustering translates into a very large population of subhalos in galaxies and affects the predictions for direct and indirect dark matter searches (gamma rays and antimatter cosmic rays). In this paper, we elaborate on previous analytic works to model the Galactic subhalo population, while consistently with current observational dynamical constraints on the Milky Way. In particular, we propose a self-consistent method to account for tidal effects induced by both dark matter and baryons. Our model does not strongly rely on cosmological simulations as they can hardly be fully matched to the real Milky Way, but for setting the initial subhalo mass fraction. Still, it allows to recover the main qualitative features of simulated systems. It can further be easily adapted to any change in the dynamical constraints, and be used to make predictions or derive constraints on dark matter candidates from indirect or direct searches. We compute the annihilation boost factor, including the subhalo-halo cross-product. We confirm that tidal effects induced by the baryonic components of the Galaxy play a very important role, resulting in a local average subhalo mass density $\lesssim 1\%$ of the total local dark matter mass density, while selecting in the most concentrated objects and leading to interesting features in the overall annihilation profile in the case of a sharp subhalo mass function. Values of global annihilation boost factors range from $\sim 2$ to $\sim 20$, while the local annihilation rate is about half as much boosted.

Galactic Cosmic-ray (CR) transport parameters are usually constrained by the boron-to-carbon ratio. This procedure is generically plagued with degeneracies between the diffusion coefficient and the vertical extent of the Galactic magnetic halo. The latter is of paramount importance for indirect dark matter (DM) searches, because it fixes the amount of DM annihilation or decay that contributes to the local antimatter CR flux. These degeneracies could be broken by using secondary radioactive species, but the current data still have large error bars, and this method is extremely sensitive to the very local interstellar medium (ISM) properties. Here, we propose to use the low-energy CR positrons in the GeV range as another direct constraint on diffusion models. We show that the PAMELA data disfavor small diffusion halo ($L\lesssim 3$ kpc) and large diffusion slope models, and exclude the minimal ({\em min}) configuration (Maurin et al. 2001, Donato et al. 2004) widely used in the literature to bracket the uncertainties in the DM signal predictions. This is complementary to indirect constraints (diffuse radio and gamma-ray emissions) and has strong impact on DM searches. Indeed this makes the antiproton constraints more robust while enhancing the discovery/exclusion potential of current and future experiments, like AMS-02 and GAPS, especially in the antiproton and antideuteron channels.

Direct detection (DD) of dark matter (DM) candidates in the $\lesssim$10 GeV mass range is very sensitive to the tail of their velocity distribution. The important quantity is the maximum WIMP speed in the observer's rest frame, i.e. in average the sum of the local Galactic escape speed $v_{\rm esc}$ and of the circular velocity of the Sun $v_c$. While the latter has been receiving continuous attention, the former is more difficult to constrain. The RAVE Collaboration has just released a new estimate of $v_{\rm esc}$ (Piffl {\em et al.}, 2014 --- P14) that supersedes the previous one (Smith {\em et al.}, 2007), which is of interest in the perspective of reducing the astrophysical uncertainties in DD. Nevertheless, these new estimates cannot be used blindly as they rely on assumptions in the dark halo modeling which induce tight correlations between the escape speed and other local astrophysical parameters. We make a self-consistent study of the implications of the RAVE results on DD assuming isotropic DM velocity distributions, both Maxwellian and ergodic. Taking as references the experimental sensitivities currently achieved by LUX, CRESST-II, and SuperCDMS, we show that: (i) the exclusion curves associated with the best-fit points of P14 may be more constraining by up to $\sim 40$% with respect to standard limits, because the underlying astrophysical correlations induce a larger local DM density; (ii) the corresponding relative uncertainties inferred in the low WIMP mass region may be moderate, down to 10-15% below 10 GeV. We finally discuss the level of consistency of these results with other independent astrophysical constraints. This analysis is complementary to others based on rotation curves.

While dark matter (DM) is the key ingredient for a successful theory of structure formation, its microscopic nature remains elusive. Indirect detection may provide a powerful test for some strongly motivated DM particle models. Nevertheless, astrophysical backgrounds are usually expected with amplitudes and spectral features similar to the chased signals. On galactic scales, these backgrounds arise from interactions of cosmic rays (CRs) with the interstellar gas, both being difficult to infer and model in detail from observations. Moreover, the associated predictions unavoidably come with theoretical errors, which are known to be significant. We show that a trustworthy guide for such challenging searches can be obtained by exploiting the full information contained in cosmological simulations of galaxies, which now include baryonic gas dynamics and star formation. We further insert CR production and transport from the identified supernova events and fully calculate the CR distribution in a simulated galaxy. We focus on diffuse gamma-rays, and self-consistently calculate both the astrophysical galactic emission and the dark matter signal. We notably show that adiabatic contraction does not necessarily induce large signal-to-noise ratios in galactic centers, and could anyway be traced from the astrophysical background itself. We finally discuss how all this may be used as a generic diagnostic tool for galaxy formation.

We make a detailed analysis of the indirect diffuse gamma-ray signals from dark matter annihilation in the Galaxy. We include the prompt emission, as well as the emission from inverse Compton scattering whenever the annihilation products contain light leptons. We consider both the contribution from the smooth dark matter halo and that from substructures. The main parameters for the latter are the mass function index and the minimal subhalo mass. We use recent results from N-body simulations to set the most reasonable range of parameters, and find that the signal can be boosted by a factor ranging from 2 to 15 towards the Galactic poles, slightly more towards the Galactic anticenter, with an important dependence on the subhalo mass index. This uncertainty is however much less than that of the extragalactic signal studied in the literature. We derive upper bounds on the dark matter annihilation cross section using the isotropic gamma-ray emission measured by Fermi-LAT, for two directions in the sky, the Galactic anticenter and the Galactic pole(s). The former represents the lowest irreducible signal from dark matter annihilation, and the latter is robust as the astrophysical background, dominated by the hadronic contribution, is rather well established in that direction. Finally, we show how the knowledge of the minimal subhalo mass, which formally depends on the dark matter particle interactions with normal matter, can be used to derive the mass function index.

The astronomical dark matter could be made of weakly interacting and massive particles. If so, these species would be abundant inside the Milky Way, where they would continuously annihilate and produce cosmic rays. Those annihilation products are potentially detectable at the Earth, and could provide indirect clues for the presence of dark matter species within the Galaxy. We will review here the various cosmic radiations which the dark matter can produce. We will examine how they propagate throughout the Milky Way and compare the dark matter yields with what pure astrophysical processes are expected to generate. The presence of dark matter substructures might enhance the signals and will be briefly discussed.

Some direct detection experiments have recently collected excess events that could be interpreted as a dark matter (DM) signal, pointing to particles in the $\sim$10 GeV mass range. We show that scenarios in which DM can self-annihilate with significant couplings to quarks are likely excluded by the cosmic-ray (CR) antiproton data, provided the annihilation is S-wave dominated when DM decouples in the early universe. These limits apply to most of supersymmetric candidates, eg in the minimal supersymmetric standard model (MSSM) and in the next-to-MSSM (NMSSM), and more generally to any thermal DM particle with hadronizing annihilation final states.

The CoGeNT experiment, dedicated to direct detection of dark matter, has recently released excess events that could be interpreted as elastic collisions of $\sim$10 GeV dark matter particles, which might simultaneously explain the still mysterious DAMA/LIBRA modulation signals, while in conflict with results from other experiments such as CDMS, XENON-100 and SIMPLE. It was shown that 5-15 GeV singlino-like dark matter candidates arising in singlet extensions of minimal supersymmetric scenarios can fit these data; annihilation then mostly proceeds into light singlet-dominated Higgs (pseudo)scalar fields. We develop an effective Lagrangian approach to confront these models with the existing data on cosmic-ray antiprotons, including the latest PAMELA data. Focusing on a parameter space consistent with the CoGeNT region, we show that the predicted antiproton flux is generically in tension with the data whenever the produced (pseudo)scalars can decay into quarks energetic enough to produce antiprotons, provided the annihilation S-wave is significant at freeze out in the early universe. In this regime, a bound on the singlino annihilation cross section is obtained, $\sigv\lesssim 10^{-26}\,{\rm cm^3/s}$, assuming a dynamically constrained halo density profile with a local value of $\rho_\odot = 0.4\,{\rm GeV/cm^3}$. Finally, we provide indications on how PAMELA or AMS-02 could further constrain or detect those configurations producing antiprotons which are not yet excluded.

The rise in the cosmic-ray positron fraction measured by the PAMELA satellite is likely due to the presence of astrophysical sources of positrons, e.g. pulsars, on the kpc scale around the Earth. Nevertheless, assessing the properties of these sources from the positron data requires a good knowledge of the secondary positron component generated by the interaction of cosmic rays with the interstellar gas. In this paper, we investigate the impact of the spectral hardening in the cosmic-ray proton and helium fluxes recently reported by the ATIC2 and CREAM balloon experiments, on the predictions of the secondary positron flux. We show that the effect is not negligible, leading to an increase of the secondary positron flux by up to $\sim$60% above $\sim$100 GeV. We provide fitting formulae that allow a straightforward utilization of our results, which can help in deriving constraints on one's favorite primary positron source, e.g. pulsars or dark matter.

We analyse predictions of the CR lepton fluxes at the Earth of both secondary and primary origins, evaluate the theoretical uncertainties, and determine their level of consistency with respect to the available data. For propagation, we use a relativistic treatment of the energy losses for which we provide useful parameterizations. We compute the secondary components by improving on the method that we derived earlier for positrons. For primaries, we estimate the contributions from astrophysical sources (supernova remnants and pulsars) by considering all known local objects within 2 kpc and a smooth distribution beyond. We find that the electron flux in the energy range 5-30 GeV is well reproduced by a smooth distant distribution of sources with index $\gamma\sim 2.3-2.4$, while local sources dominate the flux at higher energy. For positrons, local pulsars have an important effect above 5-10 GeV. Uncertainties affecting the source modeling and propagation are degenerate and each translates into about one order of magnitude error in terms of local flux. The spectral shape at high energy is weakly correlated with the spectral indices of local sources, but more strongly with the hierarchy in their distance, age and power. Despite the large theoretical errors that we describe, our global and self-consistent analysis can explain all available data without over-tuning the parameters, and therefore without the need to consider any exotic physics. Though a \emph{standard paradigm} of Galactic CRs is well established, our results show that we can hardly talk about any \emph{standard model} of CR leptons, because of the very large theoretical uncertainties. Our analysis provides details about the impact of these uncertainties, thereby sketching a roadmap for future improvements.

Recent measurements performed with some direct dark matter detection experiments, e.g. CDMS-II and CoGENT (after DAMA/LIBRA), have unveiled a few events compatible with weakly interacting massive particles. The preferred mass range is around 10 GeV, with a quite large spin-independent cross section of $10^{-43}$-$10^{-41}\,{\rm cm^2}$. In this paper, we recall that a light dark matter particle with dominant couplings to quarks should also generate cosmic-ray antiprotons. Taking advantage of recent works constraining the Galactic dark matter mass profile on the one hand and on cosmic-ray propagation on the other hand, we point out that considering a thermal annihilation cross section for such low mass candidates very likely results in an antiproton flux in tension with the current data, which should be taken into account in subsequent studies.

We study the Sunyaev-Zel'dovich effect potentially generated by relativistic electrons injected from dark matter annihilation or decay in the Galaxy, and check whether it could be observed by Planck or the Atacama Large Millimeter Array (ALMA), or even imprint the current CMB data as, e.g., the specific fluctuation excess claimed from an recent reanalysis of the WMAP-5 data. We focus on high-latitude regions to avoid contamination of the Galactic astrophysical electron foreground, and consider the annihilation or decay coming from the smooth dark matter halo as well as from subhalos, further extending our analysis to a generic modeling of spikes arising around intermediate-mass black holes. We show that all these dark Galactic components are unlikely to produce any observable Sunyaev-Zel'dovich effect. For a self-annihilating dark matter particle of 10 GeV with canonical properties, the largest optical depth we find is $\tau_e \lesssim 10^{-7}$ for massive isolated subhalos hosting intermediate-mass black holes. We conclude that dark matter annihilation or decay on the Galactic scale cannot lead to significant Sunyaev-Zel'dovich distortions of the CMB spectrum.

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.

We revisit the prospects for detecting the Sunyaev Zel'dovich (SZ) effect induced by dark matter (DM) annihilation or decay. We show that with standard (or even extreme) assumptions for DM properties, the optical depth associated with relativistic electrons injected from DM annihilation or decay is much smaller than that associated with thermal electrons, when averaged over the angular resolution of current and future experiments. For example, we find: $\tau_{\rm DM} \sim 10^{-9}-10^{-5}$ (depending on the assumptions) for $\mchi = 1$ GeV and a density profile $\rho\propto r^{-1}$ for a template cluster located at 50 Mpc and observed within an angular resolution of $10"$, compared to $\tau_{\rm th}\sim 10^{-3}-10^{-2}$. This, together with a full spectral analysis, enables us to demonstrate that, for a template cluster with generic properties, the SZ effect due to DM annihilation or decay is far below the sensitivity of the Planck satellite. This is at variance with previous claims regarding heavier annihilating DM particles. Should DM be made of lighter particles, the current constraints from 511 keV observations on the annihilation cross section or decay rate still prevent a detectable SZ effect. Finally, we show that spatial diffusion sets a core of a few kpc in the electron distribution, even for very cuspy DM profiles, such that improving the angular resolution of the instrument, e.g. with ALMA, does not necessarily improve the detection potential. We provide useful analytical formulae parameterized in terms of the DM mass, decay rate or annihilation cross section and DM halo features, that allow quick estimates of the SZ effect induced by any given candidate and any DM halo profile.

The recent observation by the PAMELA satellite of a rising positron fraction up to $\sim$ 100 GeV has triggered a considerable amount of putative interpretations in terms of dark matter (DM) annihilation or decay. Here, we make a critical reassessment of such a possibility, recalling the elementary conditions with respect to the standard astrophysical background that would make it likely, showing that they are not fulfilled. Likewise, we argue that, as now well accepted, DM would need somewhat contrived properties to contribute significantly to the observed positron signal, even when including e.g. clumpiness effects. This means that most of natural DM candidates arising in particle physics beyond the standard model are not expected to be observed in the cosmic antimatter spectrum, unfortunately. However, this does not prevent them from remaining excellent DM candidates, this only points towards the crucial need of developing much more complex detection strategies (multimessenger, multiwavelength, multiscale searches).

Secondary positrons are produced by spallation of cosmic rays within the interstellar gas. Measurements have been typically expressed in terms of the positron fraction, which exhibits an increase above 10 GeV. Many scenarios have been proposed to explain this feature, among them some additional primary positrons originating from dark matter annihilation in the Galaxy. The PAMELA satellite has provided high quality data that has enabled high accuracy statistical analyses to be made, showing that the increase in the positron fraction extends up to about 100 GeV. It is therefore of paramount importance to constrain theoretically the expected secondary positron flux to interpret the observations in an accurate way. We find the secondary positron flux to be reproduced well by the available observations, and to have theoretical uncertainties that we quantify to be as large as about one order of magnitude. We also discuss the positron fraction issue and find that our predictions may be consistent with the data taken before PAMELA. For PAMELA data, we find that an excess is probably present after considering uncertainties in the positron flux, although its amplitude depends strongly on the assumptions made in relation to the electron flux. By fitting the current electron data, we show that when considering a soft electron spectrum, the amplitude of the excess might be far lower than usually claimed. We provide fresh insights that may help to explain the positron data with or without new physical model ingredients. PAMELA observations and the forthcoming AMS-02 mission will allow stronger constraints to be aplaced on the cosmic--ray transport parameters, and are likely to reduce drastically the theoretical uncertainties.

The proposal of a galactic population of intermediate mass black holes (IMBHs), forming dark matter (DM) ``mini-spikes'' around them, has received considerable attention in recent years. In fact, leading in some scenarios to large annihilation fluxes in gamma rays, neutrinos and charged cosmic rays, these objects are sometimes quoted as one of the most promising targets for indirect DM searches. In this letter, we apply a detailed statistical analysis to point out that the existing EGRET data already place very stringent limits on those scenarios, making it rather unlikely that any of these objects will be observed with, e.g., the Fermi/GLAST satellite or upcoming Air Cherenkov telescopes. We also demonstrate that prospects for observing signals in neutrinos or charged cosmic rays seem even worse. Finally, we address the question of whether the excess in the cosmic ray positron/electron flux recently reported by PAMELA/ATIC could be due to a nearby DM point source like a DM clump or mini-spike; gamma-ray bounds, as well as the recently released Fermi cosmic ray electron and positron data, again exclude such a possibility for conventional DM candidates, and strongly constrain it for DM purely annihilating into light leptons.

Recent PAMELA and ATIC data seem to indicate an excess in positron cosmic rays above approximately 10 GeV which might be due to galactic Dark Matter particle annihilation. However the background of this signal suffers many uncertainties that make our task difficult in constraining Dark Matter or any other astrophysical explanation for these recent surprising data.

We derive the covariant formalism associated with the relativistic Sunyaev-Zel'dovich effect due to a non-thermal population of high energy electrons in clusters of galaxies. More precisely, we show that the formalism proposed by Wright in 1979, based on an empirical approach (and widely used in the literature) to compute the inverse Compton scattering of a population of relativistic electrons on CMB photons, can actually be re-interpreted as a Boltzmann-like equation, in the single scattering approximation.

While the PAMELA collaboration has recently confirmed the cosmic ray positron excess, it is interesting to review the effects of dark matter (DM) subhalos on the predicted antimatter signals. We recall that, according to general subhalo properties as inferred from theoretical cosmology, and for DM with constant annihilation cross section, the enhancement cannot be $\gtrsim$ 20 for the antimatter yield. This bound is obviously different from that found for gamma-rays. We also recall some predictions for supersymmetric benchmark models observable at the LHC and derived in the cosmological N-body framework, illustrating in the meantime the existing discrepancy between the profiles derived from N-body experiments and the current observations of the Milky Way.

[Abridged]. We take advantage of the galaxy-like 3D dark matter map extracted from the HORIZON Project results to calculate the positron and antiproton fluxes from dark matter annihilation, in a model-independent approach as well as for dark matter particle benchmarks relevant at the LHC scale (from supersymmetric and extra-dimensional theories). Such a study is dedicated to a better estimate of the theoretical uncertainties affecting predictions, while the PAMELA and GLAST satellites are currently taking data which will soon provide better observational constraints. We discuss the predictions of the antiproton and positron fluxes, and of the positron fraction as well, as compared to the current data. We finally discuss the limits of the Nbody framework in describing the dark matter halo of our Galaxy.

Anti-proton and positron Galactic cosmic ray spectra are among the key targets for indirect detection of dark matter. The boost factors, corresponding to an enhancement of the signal and linked to the clumpiness properties of the dark matter distribution, have been taken as high as thousands in the past. The dramatic impact of these boost factors for indirect detection of antiparticles, for instance with the PAMELA satellite or the coming AMS-02 experiment, asks for their detailed calculation. We take into account the results of high resolution N-body dark matter simulations to calculate the most likely energy dependent boost factors linked to the cosmic ray propagation properties, for anti-protons and positrons. Starting from the mass and space distributions of sub-halos, the anti-proton and positron propagators are used to calculate the mean value and the variance of the boost factor for the primary fluxes. We take advantage of the statistical method introduced in Lavalle et al. (2007) and cross-check the results with Monte Carlo computations. By spanning some extreme configurations of sub-halo and propagation properties, we find that the average contribution of the clumps is negligible compared to that of the smooth dark matter component. Sub-halos do not lead to enhancement of the signals, unless they are taken with some extreme (unexpected) properties. This result is independent of the nature of the self-annihilating dark matter candidate considered, and provides precise estimates of the theoretical and the statistical uncertainties of the antimatter flux from dark matter substructures. Spectral distortions can still be expected in antimatter flux measurements, but scenarios invoking large and even mild clumpiness boost factors are strongly disfavoured by our analysis.

Intermediate Mass Black Holes (IMBHs) are candidates to seed the Supermassive Black Holes (SMBHs), and some could still wander in the Galaxy. In the context of annihilating dark matter (DM), they are expected to drive huge annihilation rates, and could therefore significantly enhance the primary cosmic rays (CRs) expected from annihilation of the DM of the Galactic halo. In this proceeding (the original paper is Brun et al. 2007), we briefly explain the method to derive estimates of such exotic contributions to the anti-proton and positron CR spectra, and the associated statistical uncertainties connected to the properties of IMBHs. We find boost factors of order $10^4$ to the exotic fluxes, but associated with very large statistical uncertainties.

The hierarchical scenario of structure formation, in the frame of the $\Lambda$-CDM cosmology, predicts the existence of dark matter (DM) sub-halos down to very small scales, of which the minimal size depends on the microscopic properties of the DM. In the context of annihilating DM, such substructures are expected to enhance the primary cosmic ray (CR) fluxes originating from DM annihilation in the Galaxy. This enhancement has long been invoked to allow predictions of imprints of DM annihilation on the antimatter CR spectra. Taking advantage of the method developed by Lavalle et al (2007b), we (Lavalle et al, 2007a) accurately compute the boost factors for positrons and anti-protons, as well as the associated theoretical and statistical errors. To this aim, we use a compilation of the latest results of cosmological N-body simulations and the theoretical insights found in the literature. We find that sub-halos are not likely to significantly boost the exotic production of antimatter CRs.

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.

The small-scale distribution of dark matter in Galactic halos is poorly known. Several studies suggest that it could be very clumpy, which turns out to be of paramount importance when investigating the annihilation signal from exotic particles (e.g. supersymmetric or Kaluza-Klein). In this paper we focus on the annihilation signal in positrons. We estimate the associated uncertainty, due to the fact that we do not know exactly how the clumps are distributed in the Galactic halo. To this aim, we perform a statistical study based on analytical computations, as well as numerical simulations. In particular, we study the average and variance of the annihilation signal over many Galactic halos having the same statistical properties. We find that the so-called boost factor used by many authors should be handled with care, as i) it depends on energy and ii) it may be different for positrons, antiprotons and gamma rays, a fact which has not received any attention before. As an illustration, we use our results to discuss the positron spectrum measurements by the HEAT experiment.

We have measured the gamma-ray fluxes of the blazars Mrk 421 and Mrk 501 in the energy range between 50 and 350 GeV (1.2 to 8.3 x 10^25 Hz). The detector, called CELESTE, used first 40, then 53 heliostats of the former solar facility "Themis" in the French Pyrenees to collect Cherenkov light generated in atmospheric particle cascades. The signal from Mrk 421 is often strong. We compare its flux with previously published multi-wavelength studies and infer that we are straddling the high energy peak of the spectral energy distribution. The signal from Mrk 501 in 2000 was weak (3.4 sigma). We obtain an upper limit on the flux from 1ES 1426+428 of less than half that of the Crab flux near 100 GeV. The data analysis and understanding of systematic biases have improved compared to previous work, increasing the detector's sensitivity.

If dark matter is made of neutralinos, annihilation of such Majorana particles should produce high energy cosmic rays, especially in galaxy halo high density regions like galaxy centres. M31 (Andromeda) is our nearest neighbour spiral galaxy, and both its high mass and its low distance make it a source of interest for the indirect search for dark matter through gamma-ray detection. The ground based atmospheric Cherenkov telescope CELESTE observed M31 from 2001 to 2003, in the mostly unexplored energy range 50-500 GeV. These observations provide an upper limit on the flux above 50 GeV around $10^{-10}\rm{cm}^{-2}\rm{s}^{-1}$ in the frame of supersymmetric dark matter, and more generally on any gamma emission from M31.

It is widely believed that dark matter exists within galaxies and clusters of galaxies. Under the assumption that this dark matter is composed of the lightest, stable supersymmetric particle, assumed to be the neutralino, the feasibility of its indirect detection via observations of a diffuse gamma-ray signal due to neutralino annihilation within M31 is examined.

It is widely believed that dark matter exists within galaxies and clusters of galaxies. Under the assumption that this dark matter is composed of the lightest, stable supersymmetric particle, assumed to be the neutralino, the feasibility of its indirect detection via observations of a diffuse gamma-ray signal due to neutralino annihilations within M31 is examined. To this end, first the dark matter halo of the close spiral galaxy M31 is modeled from observations, then the resultant gamma-ray flux is estimated within supersymmetric model configurations. We conclude that under favorable conditions such as the rapid accretion of neutralinos on the central black hole in M31 and/or the presence of many clumps inside its halo with $r^{-3/2}$ inner profiles, a neutralino annihilation gamma-ray signal is marginally detectable by the ongoing collaboration CELESTE.

The distant globular cluster Palomar 13 has been found to have a very high mass-to-light ratio and its profile can be well fitted either by a King model with a tail, or with a NFW model. This cluster may be the first case of the many clumps predicted by CDM simulations that would not be disrupted by the galactic halo potential. We make the hypothesis that Pal 13 is made of neutralinos and run the DarkSuspect code to estimate the high-energy photon flux due to the annihilation of neutralinos through various channels in some benchmark scenarios. These low fluxes may be used as targets to be reached in proposals for future ground-based high altitude Cerenkov telescopes.