The global contribution of unresolved gamma-ray point sources to the extragalactic gamma-ray background has been recently measured down to gamma-ray fluxes lower than those reached with standard source detection techniques, and by employing the statistical properties of the observed gamma-ray counts. We investigate and exploit the complementarity of the information brought by the one-point statistics of photon counts (using more than 10 years of Fermi-LAT data) and by the recent measurement of the angular power spectrum of the unresolved gamma-ray background (based on 8 years of Fermi-LAT data). We determine, under the assumption that the source-count distribution of the brightest unresolved objects is dominated by blazars, their gamma-ray luminosity function and spectral energy distribution down to fluxes almost two orders of magnitude smaller than the threshold for detecting resolved sources. The different approaches provide consistent predictions for the gamma-ray luminosity function of blazars, and they show a significant complementarity.

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

The hypothesis that pulsar wind nebulae (PWNe) can significantly contribute to the excess of the positron ($e^+$) cosmic-ray flux has been consolidated after the observation of a $\gamma$-ray emission at TeV energies of a few degree size around Geminga and Monogem PWNe, and at GeV energies for Geminga at a much larger extension. The $\gamma$-ray halos around these PWNe are interpreted as due to electrons ($e^-$) and $e^+$ accelerated and escaped by their PWNe, and inverse Compton scattering low-energy photons of the interstellar radiation fields. The extension of these halos suggests that the diffusion around these PWNe is suppressed by two orders of magnitude with respect to the average in the Galaxy. We implement a two-zone diffusion model for the propagation of $e^+$ accelerated by the Galactic population of PWNe. We consider pulsars from the ATNF catalog and build up simulations of the PWN Galactic population. In both scenarios, we find that within a two-zone diffusion model, the total contribution from PWNe and secondary $e^+$ is at the level of AMS-02 data, for an efficiency of conversion of the pulsar spin down energy in $e^\pm$ of $\eta\sim0.1$. For the simulated PWNe, a $1\sigma$ uncertainty band is determined, which is of at least one order of magnitude from 10 GeV up to few TeV. The hint for a decreasing $e^+$ flux at TeV energies is found, even if it is strongly connected to the chosen value of the radius of the low diffusion bubble around each source.

The HAWC Collaboration has discovered a $\gamma$-ray emission extended about 2 degrees around the Geminga and Monogem pulsar wind nebulae (PWNe) at $\gamma$-ray energies $E_\gamma >5$ TeV. We analyze, for the first time, almost 10 years of $\gamma$-ray data obtained with the Fermi Large Area Telescope at $E_\gamma >$ 8 GeV in the direction of Geminga and Monogem. Since these two pulsars are close the Galactic plane we run our analysis with 10 different interstellar emission models (IEMs) to study the systematics due to the modeling of this component. We detect a $\gamma$-ray halo around Geminga with a significance in the range $7.8-11.8\sigma$ depending on the IEM considered. This measurement is compatible with $e^+$ and $e^-$ emitted by the PWN, which inverse-Compton scatter (ICS) with photon fields located within a distance of about 100 pc from the pulsar, where the diffusion coefficient is estimated to be around $1.1 \times 10^{27}$ cm$^2$/s at 100 GeV. We include in our analysis the proper motion of the Geminga pulsar which is relevant for $\gamma$ rays produced for ICS in the Fermi-LAT energy range. We find that an efficiency of about $1\%$ for the conversion of the spin-down energy of the pulsar into $e^+$ and $e^-$ is required to be consistent with $\gamma$-ray data from Fermi-LAT and HAWC. The inferred contribution of Geminga to the $e^+$ flux is at most $20\%$ at the highest energy AMS-02 data. Our results are compatible with the interpretation that the cumulative emission from Galactic pulsars explains the positron excess.

Recently, a few-degrees extended $\gamma$-ray halo in the direction of Geminga pulsar has been detected by HAWC, Milagro and Fermi-LAT. These observations can be interpreted with positrons ($e^+$) and electrons ($e^-$) accelerated by Geminga pulsar wind nebula (PWN), released in a Galactic environment with a low diffusion coefficient ($D_0$), and inverse Compton scattering (ICS) with the interstellar radiation fields. We inspect here how the morphology of the ICS $\gamma$-ray flux depends on the energy, the pulsar age and distance, and the strength and extension of the low-diffusion bubble. In particular we show that $\gamma$-ray experiments with a peak of sensitivity at TeV energies are the most promising ones to detect ICS halos. We perform a study of the sensitivity of HAWC, HESS and the future CTA experiment finding that, with pulsar efficiencies, for the conversion of spin-down luminosity into $e^{\pm}$, of the order of few %, the first two experiments should have already detected a few tens of ICS halos while the latter will increase the number of detections by a factor of 4. We then consider a sample of sources associated to PWNe and detected in the HESS Galactic plane survey and in the second HAWC catalog. We use the information available in these catalogs for the $\gamma$-ray spatial morphology and flux of these sources to inspect the value of $D_0$ around them and the $e^{\pm}$ injection spectrum. All sources are detected as extended with a $\gamma$-ray emission extended about $15-80$ pc. Assuming that most of the $e^{\pm}$ accelerated by these sources have been released in the interstellar medium, the diffusion coefficient is $2-30 \cdot 10^{26}$ cm$^2$/s at 1 TeV, i.e. two orders of magnitude smaller than the value considered to be the average in the Galaxy. These observations imply that Galactic PWNe have low-diffusion bubbles with a size of at least 80 pc.

The data on the inclusive flux of cosmic positrons and electrons ($e^++e^{-}$) have been recently collected from GeV to tens of TeV energies by several experiments with unprecedented precision. In addition, the Fermi-LAT Collaboration has provided a new energy spectrum for the upper bounds on the $e^++e^{-}$ dipole anisotropy. This observable can bring information on the emission from local Galactic sources, notably measured with high precision at radio frequencies. We develop a framework in which $e^-$ and $e^+$ measured at Earth from GeV up to tens of TeV energies have a composite origin. A dedicated analysis is deserved to Vela YZ and Cygnus Loop Supernova Remnants (SNRs), for which we consider two different models for the injection of $e^-$. We investigate the consistency of these models using the three physical observables: the radio flux from Vela YZ and Cygnus Loop at all the available frequencies, the $e^++e^-$ flux from five experiments from the GeV to tens of TeV energy, the $e^++e^-$ dipole anisotropy upper limits from 50 GeV to about 1 TeV. We find that the radio flux for these nearby SNRs strongly constraints the properties of the injection electron spectrum, partially compatible with the looser constraints derived from the $e^+ + e^-$ flux data. We also perform a multi-wavelength multi-messenger analysis by fitting simultaneously the radio flux on Vela YZ and Cygnus Loop and the $e^+ + e^-$ flux, and checking the outputs against the $e^+ + e^-$ dipole anisotropy data. Remarkably, we find a model which is compatible with all the $e^++e^-$ flux data, the radio data for Vela YZ and Cygnus Loop, and with the anisotropy upper bounds. We show the severe constraints imposed by the most recent data on the $e^+ + e^-$ dipole anisotropy.

The origin of cosmic-ray positrons detected with an energy above 10 GeV is one of the most intriguing mysteries in Astroparticle Physics. Different interpretations have been invoked to solve this puzzle such as pulsar wind nebulae (PWNe), supernovae remnants or interactions of dark matter particles. The HAWC and Milagro experiments have measured an extended emission of photons above 10 TeV in the direction of Geminga and Monogem PWNe which can be interpreted by photons emitted through inverse Compton scattering by positrons and electrons accelerated by these sources. These HAWC and Milagro detections are extremely important to estimate the PWN contribution to the positron excess. Positrons and electrons emitted by PWNe and injected in the interstellar medium can produce photons with an energy between radio and X-ray through synchrotron radiation with the Galactic magnetic field. In this white paper we show that Astrophysics Probe mission concepts, such as AMEGO (the All-sky Medium Energy Gamma-ray Observatory) and AdEPT (The Advanced Energetic Pair Telescope), are suitable to detect the synchrotron halos from Monogem and Geminga PWNe and we will report the list of the most promising PWNe that AMEGO could be able to detect.

Gamma-ray searches for new physics such as dark matter are often driven by investigating the composition of the extragalactic gamma-ray background (EGB). Classic approaches to EGB decomposition manifest in resolving individual point sources and dissecting the intensity spectrum of the remaining unresolved component. Furthermore, statistical methods have recently been proven to outperform the sensitivity of classic source detection algorithms in finding point-source populations in the unresolved flux regime. In this article, we employ the 1-point photon count statistics of eight years of Fermi-LAT data to resolve the population of extragalactic point sources and to decompose the diffuse isotropic background contribution for Galactic latitudes |b|>30 deg. We use three adjacent energy bins between 1 and 10 GeV. For the first time, we extend the analysis to incorporate a potential contribution from annihilating dark matter smoothly distributed in the Galaxy. We investigate the sensitivity reach of 1-point statistics for constraining the thermally-averaged self-annihilation cross section <\sigma v> of dark matter, using different template models for the Galactic foreground emission. Given the official Fermi-LAT interstellar emission model, we set upper bounds on the DM self-annihilation cross section <\sigma v> that are comparable with the constraints obtained by other indirect detection methods, in particular by the stacking analysis of several dwarf spheroidal galaxies.

Theoretical predictions for the cosmic antiproton spectrum currently fall short of the corresponding experimental level of accuracy. Among the main sources of uncertainty are the antiproton production cross sections in cosmic ray inelastic interactions. We analyse existing data on antiproton production in $pp$ scattering, including for the first time the measurements performed by the NA49 Collaboration. We compute the antiproton spectrum finding that in the energy range where data are available (antiproton energies of about 4-550 GeV) different approaches lead to almost equivalent results, with an uncertainty of 10-20\%. Extrapolations outside this region lead to different estimates, with the uncertainties reaching the 50\% level around $1$ TeV, degrading the diagnostic power of the antiproton channel at those energies. We also comment on the uncertainties in the antiproton production source term coming from nuclei heavier than protons and from neutrons produced in $pp$ scatterings, and point out the need for dedicated experimental campaigns for all processes involving antiproton production in collisions of light nuclei.

After the breakthrough from the satellite-borne PAMELA detector, the flux of cosmic-ray (CR) antiprotons has been provided with unprecedented accuracy by AMS-02 on the International Space Station. Its data spans an energy range from below 1 GeV up to 400 GeV and most of the data points contain errors below the amazing level of 5%. The bulk of the antiproton flux is expected to be produced by the scatterings of CR protons and helium off interstellar hydrogen and helium atoms at rest. The modeling of these interactions, which requires the relevant production cross sections, induces an uncertainty in the determination of the antiproton source term that can even exceed the uncertainties in the CR $\bar{p}$ data itself. The aim of the present analysis is to determine the uncertainty required for $p+p\rightarrow \bar{p} + X$ cross section measurements such that the induced uncertainties on the $\bar{p}$ flux are at the same level. Our results are discussed both in the center-of-mass reference frame, suitable for collider experiments, and in the laboratory frame, as occurring in the Galaxy. We find that cross section data should be collected with accuracy better that few percent with proton beams from 10 GeV to 6 TeV and a pseudorapidity $\eta$ ranging from 2 to almost 8 or, alternatively, with $p_T$ from 0.04 to 2 GeV and $x_R$ from 0.02 to 0.7. Similar considerations hold for the $p$He production channel. The present collection of data is far from these requirements. Nevertheless, they could, in principle, be reached by fixed target experiments with beam energies in the reach of CERN accelerators.

The cosmic-ray flux of antiprotons is measured with high precision by the space-borne particle spectrometers AMS-02.Its interpretation requires a correct description of the dominant production process for antiprotons in our Galaxy, namely, the interaction of cosmic-ray proton and helium with the interstellar medium. In the light of new cross section measurements by the NA61 experiment of $p + p \rightarrow \bar{p} + X$ and the first ever measurement of $p + \mathrm{He} \rightarrow \bar{p} + X$ by the LHCb experiment, we update the parametrization of proton-proton and proton-nucleon cross sections.We find that the LHCb $p$He data constrain a shape for the cross section at high energies and show for the first time how well the rescaling from the $pp$ channel applies to a helium target. By using $pp$, $p$He and $p$C data we estimate the uncertainty on the Lorentz invariant cross section for $p + \mathrm{He} \rightarrow \bar{p} + X$. We use these new cross sections to compute the source term for all the production channels, considering also nuclei heavier than He both in cosmic rays and the interstellar medium. The uncertainties on the total source term are at the level of $\pm20$% and slightly increase below antiproton energies of 5 GeV. This uncertainty is dominated by the $p+p \rightarrow \bar{p} + X$ cross section, which translates into all channels since we derive them using the $pp$ cross sections. The cross sections to calculate the source spectra from all relevant cosmic-ray isotopes are provided in the Supplemental Material. We finally quantify the necessity of new data on antiproton production cross sections, and pin down the kinematic parameter space which should be covered by future data.

Cosmic rays are an important tool to study dark matter annihilation in our Galaxy. Recently, a possible hint for dark matter annihilation was found in the antiproton spectrum measured by AMS-02, even though the result might be affected by theoretical uncertainties. A complementary way to test its dark matter interpretation would be the observation of low-energy antinuclei in cosmic rays. We determine the chances to observe antideuterons with GAPS and AMS-02, and the implications for the ongoing AMS-02 antihelium searches. We find that the corresponding antideuteron signal is within the GAPS and AMS-02 detection potential. If, more conservatively, the putative signal was considered as an upper limit on DM annihilation, our results would indicate the highest possible fluxes for antideuterons and antihelium compatible with current antiproton data.

The Fermi Large Area Telescope (LAT) Collaboration has recently released the Third Catalog of Hard Fermi-LAT Sources (3FHL), which contains 1556 sources detected above 10 GeV with seven years of Pass 8 data. We investigate the source count distribution of 3FHL sources at Galactic latitudes $|b|>20^{\circ}$, where the sources are mostly blazars. We use two complementary techniques: 1) a source-detection efficiency correction method and 2) an analysis of pixel photon count statistics with the 1-point probability distribution function (1pPDF). With the first method, using realistic Monte Carlo simulations of the $\gamma$-ray sky, we calculate the efficiency of the LAT to detect point sources. This enables us to find the intrinsic source count distribution at photon fluxes down to $7.5\times10^{-12}$ ph cm$^{-2}$s$^{-1}$. With this method we detect a flux break at $(3.5\pm0.4) \times 10^{-11} \mathrm{ph}\,\mathrm{cm}^{-2}\mathrm{s}^{-1}$ with a significance of at least $5.4 \sigma$. The power-law indexes of the source count distribution above and below the break are $2.09\pm0.04$ and $1.07\pm0.27$, respectively. This result is confirmed with the 1pPDF method, which has a sensitivity reach of $\sim10^{-11}$ ph cm$^{-2}$s$^{-1}$. Integrating the derived source count distribution above the sensitivity of our analysis, we find that $(42\pm8)\%$ of the extragalactic $\gamma$-ray background originates from blazars.

A new estimation of the isotropic diffuse gamma-ray background (IGRB) observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) has been presented for 50 months of data, in the energy range 100 MeV-820 GeV and for different modelings of the Galactic foreground. We attempt here the interpretation of the Fermi-LAT IGRB data in terms of the gamma-ray unresolved emission from different extragalactic populations. We find very good fits to the experimental IGRB, obtained with theoretical predictions for the emission from active galactic nuclei and star-forming galaxies. In addition, we probe a possible emission coming from the annihilation of weakly interacting dark matter (DM) particles in the halo of our Galaxy. We set stringent limits on its annihilation cross section into gamma rays, which are about the thermal relic value for a wide range of DM masses. We also identify regions in the DM mass and annihilation cross section parameter space which can significantly improve the fit to the IGRB data. Our analysis is conducted within the different IGRB data sets obtained from different models for the Galactic emission, which is shown to add a significant ambiguity on the IGRB interpretation.

The flux of positrons and electrons ($e^+ + e^-$) has been measured by the $Fermi$ Large Area Telescope (LAT) in the energy range between 7~GeV and 2~TeV. We discuss a number of interpretations of Pass 8 $Fermi$-LAT $e^+ + e^-$ spectrum, combining electron and positron emission from supernova remnants (SNRs) and pulsar wind nebulae (PWNe), or produced by the collision of cosmic rays with the interstellar medium. We find that the $Fermi$-LAT spectrum is compatible with the sum of electrons from a smooth SNR population, positrons from cataloged PWNe, and a secondary component. If we include in our analysis constraints from AMS-02 positron spectrum, we obtain a slightly worse fit to the $e^+ + e^-$ $Fermi$-LAT spectrum, depending on the propagation model. As an additional scenario, we replace the smooth SNR component within 0.7 kpc with the { individual sources} found in Green's catalog of Galactic SNRs. We find that separate consideration of far and near sources helps to reproduce the $e^+ + e^-$ $Fermi$-LAT spectrum. However, we show that the fit degrades when the radio constraints on the positron emission from Vela SNR (which is the main contributor at high energies) are taken into account. We find that a break in the power-law injection spectrum at about 100 GeV can also reproduce the measured $e^+ + e^-$ spectrum} and, among the cosmic-ray propagation models that we consider, no reasonable break of the power-law dependence of the diffusion coefficient can modify the electron flux enough to reproduce the observed shape.

Photon counts statistics have recently been proven to provide a sensitive observable for characterizing gamma-ray source populations and for measuring the composition of the gamma-ray sky. In this work, we generalize the use of the standard 1-point probability distribution function (1pPDF) to decompose the high-latitude gamma-ray emission observed with Fermi-LAT into: (i) point-source contributions, (ii) the Galactic foreground contribution, and (iii) a diffuse isotropic background contribution. We analyze gamma-ray data in five adjacent energy bands between 1 and 171 GeV. We measure the source-count distribution dN/dS as a function of energy, and demonstrate that our results extend current measurements from source catalogs to the regime of so far undetected sources. Our method improves the sensitivity for resolving point-source populations by about one order of magnitude in flux. The dN/dS distribution as a function of flux is found to be compatible with a broken power law. We derive upper limits on further possible breaks as well as the angular power of unresolved sources. We discuss the composition of the gamma-ray sky and capabilities of the 1pPDF method.

High energy cosmic ray electrons and positrons probe the local properties of our Galaxy. In fact, electromagnetic energy losses limit the typical propagation scale of GeV-TeV electrons and positrons to a few kpc. In the diffusion model, nearby and dominant sources may produce an observable dipole anisotropy in the cosmic ray fluxes. We present a detailed study on the role of anisotropies from nearby sources in the interpretation of the observed GeV-TeV cosmic ray electron and positrons fluxes. We compute predictions for the anisotropies from known astrophysical sources as supernova remnants and pulsar wind nebulae of the ATNF catalog. Our results are compared with current anisotropy upper limits from the Fermi- LAT, AMS-02 and PAMELA experiments.

The cosmic electrons and positrons have been measured with unprecedented statistics up to several hundreds GeV, thus permitting to explore the role that close single sources can have in shaping the flux at different energies. In the present analysis, we consider electrons and positrons in cosmic rays to be produced by spallations of hadron fluxes with the interstellar medium, by a smooth Supernova Remnant (SNR) population, by all the ATNF catalog pulsars, and by few discrete, local SNRs. We test several source models on the $e^++e^-$ and $e^+$ AMS-02 flux data. For the configurations compatible with the data, we compute the dipole anisotropy in $e^++e^-$, $e^+$, $e^+/e^-$ from single sources. Our study includes a dedicated analysis to the Vela SNR. We show that Fermi-LAT present data on dipole anisotropy of $e^++e^-$ start to explore some of the models for the Vela SNR selected by AMS-02 flux data. We also investigate how the observed anisotropy could result from a combination of local sources. Our analysis shows that the search of anisotropy in the lepton fluxes up to TeV energies can be an interesting tool for the inspection of properties of close SNRs, complementary to the high precision flux data.

Numerical simulations of structure formation have recorded a remarkable progress in the recent years, in particular due to the inclusion of baryonic physics evolving with the dark matter component. We generate Monte Carlo realizations of the dark matter sub-halo population based on the results of the recent hydrodynamical simulation suite of Milky Way-sized galaxies. We then simulate the gamma-ray sky for both the setup of the 3FGL and 2FHL Fermi Large Area Telescope (LAT) catalogs, including the contribution from the annihilation of dark matter in the sub-halos. We find that the flux sensitivity threshold strongly depends on the particle dark matter mass, and more mildly also on its annihilation channel and the observation latitude. The results differ for the 3FGL and 2FHL catalogs, given their different energy thresholds. We also predict that the number of dark matter sub-halos among the unassociated sources is very small. A null number of detectable sub-halos in the Fermi-LAT 3FGL catalog would imply upper limits on the dark matter annihilation cross section into $b\bar{b}$ of $2 \cdot 10^{-26}$ ($5 \cdot 10^{-25}$) cm$^3$/s with $M_{\rm DM}$= 50 (1000) GeV. We find less than one extended sub-halo in the Fermi-LAT 3FGL catalog. As a matter of fact, the differences in the spatial and mass distribution of sub-halos between hydrodynamic and dark matter-only runs do not have significant impact on the gamma-ray dark matter phenomenology.

Analogously to globular clusters, the dense stellar environment of the Galactic center has been proposed to host a large population of as-yet undetected millisecond pulsars (MSPs). Recently, this hypothesis found support in the analysis of gamma rays from the inner Galaxy seen by the Large Area Telescope (LAT) aboard the Fermi satellite, which revealed a possible excess of diffuse GeV photons in the inner 15 deg about the Galactic center (Fermi GeV excess). The excess can be interpreted as the collective emission of thousands of MSPs in the Galactic bulge, with a spherical distribution that strongly peaks towards the Galactic center. In order to fully establish the MSP interpretation, it is essential to find corroborating evidence in multi-wavelength searches, most notably through the detection of radio pulsation from individual bulge MSPs. Based on globular cluster observations and the gamma-ray emission from the inner Galaxy, we investigate the prospects for detecting MSPs in the Galactic bulge. While previous pulsar surveys failed to identify this population, we demonstrate that, in the upcoming years, new large-area surveys with focus on regions a few degrees north or south of the Galactic center should lead to the detection of dozens of bulge MSPs. Additionally, we show that, in the near future, deep targeted searches of unassociated Fermi sources should be able to detect the first few MSPs in the bulge. The prospects for these deep searches are enhanced by a tentative gamma-ray/radio correlation that we infer from high-latitude gamma-ray MSPs. Such detections would constitute the first clear discoveries of field MSPs in the Galactic bulge, with far-reaching implications for gamma-ray observations, the formation history of the central Milky Way and strategy optimization for future radio observations.

The source-count distribution as a function of their flux, dN/dS, is one of the main quantities characterizing gamma-ray source populations. We employ statistical properties of the Fermi Large Area Telescope (LAT) photon counts map to measure the composition of the extragalactic gamma-ray sky at high latitudes (|b|>30 deg) between 1 GeV and 10 GeV. We present a new method, generalizing the use of standard pixel-count statistics, to decompose the total observed gamma-ray emission into (a) point-source contributions, (b) the Galactic foreground contribution, and (c) a truly diffuse isotropic background contribution. Using the 6-year Fermi-LAT data set (P7REP), we show that the dN/dS distribution in the regime of so far undetected point sources can be consistently described with a power law of index between 1.9 and 2.0. We measure dN/dS down to an integral flux of ~2x10^{-11} cm^{-2}s^{-1}, improving beyond the 3FGL catalog detection limit by about one order of magnitude. The overall dN/dS distribution is consistent with a broken power law, with a break at 2.1^{+1.0}_{-1.3}x10^{-8} cm^{-2}s^{-1}. The power-law index n_1=3.1^{+0.7}_{-0.5} for bright sources above the break hardens to n_2=1.97+-0.03 for fainter sources below the break. A possible second break of the dN/dS distribution is constrained to be at fluxes below 6.4x10^{-11} cm^{-2}s^{-1} at 95% confidence level. The high-latitude gamma-ray sky between 1 GeV and 10 GeV is shown to be composed of ~25% point sources, ~69.3% diffuse Galactic foreground emission, and ~6% isotropic diffuse background.

Statistical properties of photon count maps have recently been proven as a new tool to study the composition of the gamma-ray sky with high precision. We employ the 1-point probability distribution function of 6 years of Fermi-LAT data to measure the source-count distribution dN/dS and the diffuse components of the high-latitude gamma-ray sky as a function of energy. To that aim, we analyze the gamma-ray emission in five adjacent energy bands between 1 GeV and 171 GeV. It is demonstrated that the source-count distribution as a function of flux is compatible with a broken power law up to energies of ~50 GeV. The index below the break is between 1.95 and 2.0. For higher energies, a simple power-law fits the data, with an index of 2.2^{+0.7}_{-0.3} in the energy band between 50 GeV and 171 GeV. Upper limits on further possible breaks as well as the angular power of unresolved sources are derived. We find that point-source populations probed by this method can explain 83^{+7}_{-13}% (81^{+52}_{-19}%) of the extragalactic gamma-ray background between 1.04 GeV and 1.99 GeV (50 GeV and 171 GeV). The method has excellent capabilities for constraining the gamma-ray luminosity function and the spectra of unresolved blazars.

We perform a detailed quantitative analysis of the recent AMS-02 electron and positron data. We investigate the interplay between the emission from primary astrophysical sources, namely Supernova Remnants and Pulsar Wind Nebulae, and the contribution from a dark matter annihilation or decay signal. Our aim is to assess the information that can be derived on dark matter properties when both dark matter and primary astrophysical sources are assumed to jointly contribute to the leptonic observables measured by the AMS-02 experiment. We investigate both the possibility to set robust constraints on the dark matter annihilation/decay rate and the possibility to look for dark matter signals within realistic models that take into account the full complexity of the astrophysical background. Our results show that AMS-02 data enable to probe efficiently vast regions of the dark matter parameter space and, in some cases, to set constraints on the dark matter annihilation/decay rate that are comparable or even stronger than the ones derived from other indirect detection channels.

Blazars represent the most abundant class of high-energy extragalactic $\gamma$-ray sources. The subset of blazars known as BL Lac objects is on average closer to Earth and characterized by harder spectra at high energy than the whole sample. The fraction of BL Lacs that is too dim to be detected and resolved by current $\gamma$-ray telescopes is therefore expected to contribute to the high-energy isotropic diffuse $\gamma$-ray background (IGRB). The IGRB has been recently measured over a wide energy range by the Large Area Telescope (LAT) on board the Gamma-ray Space Telescope ({\it Fermi}). We present a new prediction of the diffuse $\gamma$-ray flux due to the unresolved BL Lac blazar population. The model is built upon the spectral energy distribution and the luminosity function derived from the fraction of BL Lacs detected (and spectrally characterized) in the $\gamma$-ray energy range. We focus our attention on the ${\cal O}(100)$ GeV energy range, predicting the emission up to the TeV scale and taking into account the absorption on the extragalactic background light. In order to better shape the BL Lac spectral energy distribution, we combine the {\it Fermi}-LAT data with Imaging Atmospheric Cerenkov Telescopes measurements of the most energetic sources. Our analysis is carried on separately for low- and intermediate-synchrotron-peaked BL Lacs on one hand, and high-synchrotron-peaked BL Lacs on the other one: we find in fact statistically different features for the two. The diffuse emission from the sum of both BL Lac classes increases from about 10$\%$ of the measured IGRB at 100 MeV to $\sim$100$\%$ of the data level at 100 GeV. At energies greater than 100 GeV, our predictions naturally explain the IGRB data, accommodating their softening with increasing energy. Uncertainties are estimated to be within of a factor of two of the best-fit flux up to 500 GeV.

Recent years have seen increased theoretical and experimental effort towards the first-ever detection of cosmic-ray antideuterons, in particular as an indirect signature of dark matter annihilation or decay. In contrast to indirect dark matter searches using positrons, antiprotons, or gamma-rays, which suffer from relatively high and uncertain astrophysical backgrounds, searches with antideuterons benefit from very suppressed conventional backgrounds, offering a potential breakthrough in unexplored phase space for dark matter. This article is based on the first dedicated cosmic-ray antideuteron workshop, which was held at UCLA in June 2014. It reviews broad classes of dark matter candidates that result in detectable cosmic-ray antideuteron fluxes, as well as the status and prospects of current experimental searches. The coalescence model of antideuteron production and the influence of antideuteron measurements at particle colliders are discussed. This is followed by a review of the modeling of antideuteron propagation through the magnetic fields, plasma currents, and molecular material of our Galaxy, the solar system, the Earth's geomagnetic field, and the atmosphere. Finally, the three ongoing or planned experiments that are sensitive to cosmic-ray antideuterons, BESS, AMS-02, and GAPS, are detailed. As cosmic-ray antideuteron detection is a rare event search, multiple experiments with orthogonal techniques and backgrounds are essential. Many theoretical and experimental groups have contributed to these studies over the last decade, this review aims to provide the first coherent discussion of the relevant dark matter theories that antideuterons probe, the challenges to predictions and interpretations of antideuteron signals, and the experimental efforts toward cosmic antideuteron detection.

The precise measurement of cosmic-ray antiparticles serves as important means for identifying the nature of dark matter. Recent years showed that identifying the nature of dark matter with cosmic-ray positrons and higher energy antiprotons is difficult, and has lead to a significantly increased interest in cosmic-ray antideuteron searches. Antideuterons may also be generated in dark matter annihilations or decays, offering a potential breakthrough in unexplored phase space for dark matter. Low-energy antideuterons are an important approach because the flux from dark matter interactions exceeds the background flux by more than two orders of magnitude in the low-energy range for a wide variety of models. This review is based on the "dbar14 - dedicated cosmic-ray antideuteron workshop", which brought together theorists and experimentalists in the field to discuss the current status, perspectives, and challenges for cosmic-ray antideuteron searches and discusses the motivation for antideuteron searches, the theoretical and experimental uncertainties of antideuteron production and propagation in our Galaxy, as well as give an experimental cosmic-ray antideuteron search status update. This report is a condensed summary of the article "Review of the theoretical and experimental status of dark matter identification with cosmic-ray antideuteron" (arXiv:1505.07785).

Recent data on Galactic cosmic-ray (CR) leptons and hadrons gave rise to two exciting problems: on the lepton side, the origin of the rise of the CR positron fraction e+/(e- + e+) at ~10 - 300 GeV of energy; on the hadron side, the nature of the spectral hardening observed in CR protons and nuclei at ~TeV energies. The lepton anomaly indicates the existence of a nearby e+/- source. It has been proposed that high-energy positrons can be produced inside nearby supernova remnants (SNRs) via interactions of CR hadrons with the ambient medium. A distinctive prediction of this mechanism is a high-energy rise of the boron-to-carbon ratio, which has not been observed. It also requires old SNRs at work (with ineffective magnetic field amplification and slow shock speed), that cannot account for the CR hadronic spectra observed up to the knee energies (~5 PeV). We propose a new picture where, in addition to such a nearby CR accelerator, the high-energy spectrum of CR hadrons is provided by the large-scale population of SNRs, on average younger, that can efficiently accelerate CRs up to the knee. Under this scenario, the spectral hardening of CR hadrons can be naturally interpreted as the transition between the two components. As we will show, our two-component model breaks the connection between the positron fraction and the boron-to-carbon ratio, which is now predicted to decrease with energy in accordance with the data. Forthcoming data from AMS will be crucial for testing this model.

Gamma rays from young pulsars and milli-second pulsars are expected to contribute to the diffuse gamma-ray emission measured by the {\it Fermi} Large Area Telescope (LAT) at high latitudes. We derive the contribution of the pulsars undetected counterpart by using information from radio to gamma rays and we show that they explain only a small fraction of the isotropic diffuse gamma-ray background.

Radio-loud active galactic nuclei (AGN) are expected to contribute substantially to both the intensity and anisotropy of the isotropic gamma-ray background (IGRB). In turn, the measured properties of the IGRB can be used to constrain the characteristics of proposed contributing source classes. We consider individual subclasses of radio-loud AGN, including low-, intermediate-, and high-synchrotron-peaked BL Lacertae objects, flat-spectrum radio quasars, and misaligned AGN. Using updated models of the gamma-ray luminosity functions of these populations, we evaluate the energy-dependent contribution of each source class to the intensity and anisotropy of the IGRB. We find that collectively radio-loud AGN can account for the entirety of the IGRB intensity and anisotropy as measured by the Fermi Large Area Telescope (LAT). Misaligned AGN provide the bulk of the measured intensity but a negligible contribution to the anisotropy, while high-synchrotron-peaked BL Lacertae objects provide the dominant contribution to the anisotropy. In anticipation of upcoming measurements with the Fermi-LAT and the forthcoming Cherenkov Telescope Array, we predict the anisotropy in the broader energy range that will be accessible to future observations.

Millisecond Pulsars are second most abundant source population discovered by the Fermi-LAT. They might contribute non-negligibly to the diffuse emission measured at high latitudes by Fermi-LAT, the IDGRB. Gamma-ray sources also contribute to the anisotropy of the IDGRB measured on small scales by Fermi-LAT. We aim to assess the contribution of the unresolved counterpart of the detected MSPs population to the IDGRB and the maximal fraction of the measured anisotropy produced by this source class. We model the MSPs spatial distribution in the Galaxy and the gamma-ray emission parameters by considering radio and gamma-ray observational constraints. By simulating a large number of MSPs populations, we compute the average diffuse emission and the anisotropy 1-sigma upper limit. The emission from unresolved MSPs at 2 GeV, where the peak of the spectrum is located, is at most 0.9% of the measured IDGRB above 10 degrees in latitude. The 1-sigma upper limit on the angular power for unresolved MSP sources turns out to be about a factor of 60 smaller than Fermi-LAT measurements above 30 degrees. Our results indicate that this galactic source class represents a negligible contributor to the high-latitude gamma-ray sky and confirm that most of the intensity and geometrical properties of the measured diffuse emission are imputable to other extragalactic source classes. Nevertheless, given the MSP distribution, we expect them to contribute significantly to the gamma-ray diffuse emission at low latitudes. Since, along the galactic disk, the population of young Pulsars overcomes in number the one of MSPs, we compute the gamma-ray emission from the whole population of unresolved Pulsars in two low-latitude regions: the inner Galaxy and the galactic center.

The annihilation of dark matter particles in the halo of galaxies may end up into gamma rays, which travel almost unperturbed till to their detection at Earth. This annihilation signal can exhibit an anisotropic behavior quantified by the angular power spectrum, whose properties strongly depend on the dark matter distribution and its clumpiness. We use high resolution pure dark matter N-body simulations to quantify the contribution of different components (main halo and satellites) to the global signal as a function of the analytical profile adopted to describe the numerical results. We find that the smooth main halo dominates the angular power spectrum of the gamma-ray signal up to quite large multipoles, where the sub-haloes anisotropy signal starts to emerge, but the transition multipole strongly depends on the assumed radial profile. The extrapolation down to radii not resolved by current numerical simulations can affect both the normalization and the shape of the gamma-ray angular power spectrum. For the sub-haloes described by an asymptotically cored dark matter distribution, the angular power spectrum shows an overall smaller normalization and a flattening at high multipoles. Our results show the criticality of the dark matter density profile shape in gamma-ray anisotropy searches, and evaluate quantitatively the intrinsic errors occurring when extrapolating the dark matter radial profiles down to spatial scales not yet explored by numerical simulations.

We perform a combined analysis of the recent AMS-02 data on electrons, positrons, electrons plus positrons and positron fraction, in a self-consistent framework where we realize a theoretical modeling of all the astrophysical components that can contribute to the observed fluxes in the whole energy range. The primary electron contribution is modeled through the sum of an average flux from distant sources and the fluxes from the local supernova remnants in the Green catalog. The secondary electron and positron fluxes originate from interactions on the interstellar medium of primary cosmic rays, for which we derive a novel determination by using AMS-02 proton and helium data. Primary positrons and electrons from pulsar wind nebulae in the ATNF catalog are included and studied in terms of their most significant (while loosely known) properties and under different assumptions (average contribution from the whole catalog, single dominant pulsar, a few dominant pulsars). We obtain a remarkable agreement between our various modeling and the AMS-02 data for all types of analysis, demonstrating that the whole AMS-02 leptonic data admit a self-consistent interpretation in terms of astrophysical contributions.

The nature of the Isotropic $\gamma$-ray Background (IGRB) measured by the Large Area Telescope (LAT) on the Fermi $\gamma$-ray space Telescope ({\it Fermi}) remains partially unexplained. Non-negligible contributions may originate from extragalactic populations of unresolved sources such as blazars, star-forming galaxies or galactic milli-second pulsars. A recent prediction of the diffuse $\gamma$-ray emission from Active Galactic Nuclei (AGN) with a large viewing angle with respect to the line-of-sight (l.o.s.) has demonstrated that this faint but numerous population is also expected to contribute significantly to the total IGRB intensity. A more exotic contribution to the IGRB invokes the pair annihilation of dark matter (DM) weakly interacting massive particles (WIMPs) into $\gamma$-rays. In this work, we evaluate the room left for galactic DM at high latitudes ($>10^\circ$) by including photons from both prompt emission and inverse Compton scattering, emphasizing the impact of the newly discovered contribution from misaligned AGN (MAGN) for such an analysis. Summing up all significant galactic and extragalactic components of the IGRB, we find that an improved understanding of the associated astrophysical uncertainties is still mandatory to put stringent bounds on thermally produced DM. On the other hand, we also demonstrate that the IGRB has the potential to be one of the most competitive {\it future} ways to test the DM WIMP hypothesis, once the present uncertainties are even slightly reduced. In fact, if MAGN contribute even at 90% of the maximal level consistent with our current understanding, thermally produced WIMPs would be severely constrained as DM candidates for masses up to several TeV.

Active galactic nuclei (AGN) with jets seen at small viewing angles are the most luminous and abundant objects in the $\gamma$-ray sky. AGN with jets misaligned along the line-of-sight appear fainter in the sky, but are more numerous than the brighter blazars. We calculate the diffuse $\gamma$-ray emission due to the population of misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the {\it Fermi} Gamma-ray Space Telescope ({\it Fermi}). A correlation between the $\gamma$-ray luminosity and the radio-core luminosity is established and demonstrated to be physical by statistical tests, as well as compatible with upper limits based on {\it Fermi}-LAT data for a large sample of radio-loud MAGN. We constrain the derived $\gamma$-ray luminosity function by means of the source count distribution of the radio galaxies (RGs) detected by the {\it Fermi}-LAT. We finally calculate the diffuse $\gamma$-ray flux due to the whole MAGN population. Our results demonstrate that the MAGN can contribute from 10% up to nearly the entire measured Isotropic Gamma-Ray Background (IGRB). We evaluate a theoretical uncertainty on the flux of almost an order of magnitude.

We calculate the diffuse $\gamma$-ray emission due to the population of misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the {\it Fermi} Gamma-ray Space Telescope ({\it Fermi}). A correlation between the $\gamma$-ray luminosity and the radio-core luminosity is established and demonstrated to be physical by statistical tests, as well as compatible with upper limits based on {\it Fermi}-LAT data for a large sample of radio-loud MAGN. We constrain the derived $\gamma$-ray luminosity function by means of the source count distribution of the MAGN detected by the {\it Fermi}-LAT. We finally estimate the diffuse $\gamma$-ray flux due to the whole MAGN population which ranges from 10% up to nearly the entire measured Isotropic Gamma-Ray Background (IGRB). We evaluate also the room left to galactic DM at high latitudes ($>10^\circ$), by taking into account the results on the MAGN together with the other significant galactic and extragalactic $\gamma$-rays emitting sources.

The secondary-to-primary B/C ratio is widely used to study the cosmic ray (CR) propagation processes in the Galaxy. It is usually assumed that secondary nuclei such as Li-Be-B are entirely generated by collisions of heavier CR nuclei with the interstellar medium (ISM). We study the CR propagation under a scenario where secondary nuclei can also be produced or accelerated from galactic sources. We consider the processes of hadronic interactions inside supernova remnants (SNRs) and re-acceleration of background CRs in strong shocks. Thus, we investigate their impact in the propagation parameter determination within present and future data. The spectra of Li-Be-B nuclei emitted from SNRs are harder than those due to CR collisions with the ISM. The secondary-to-primary ratios flatten significantly at ~TeV/n energies, both from spallation and re-acceleration in the sources. The two mechanisms are complementary to each other and depend on the properties of the local ISM around the expanding remnants. The secondary production in SNRs is significant for dense background media, n ~1 cm^-3, while the amount of re-accelerated CRs is relevant for SNRs expanding into rarefied media, n ~0.1 cm-3. Due to these effects, the the diffusion parameter 'delta' may be misunderstood by a factor of ~5-15%. Our estimations indicate that an experiment of the AMS-02 caliber can constrain the key propagation parameters while breaking the source-transport degeneracy, for a wide class of B/C-consistent models. Given the precision of the data expected from on-going experiments, the SNR production/acceleration of secondary nuclei should be considered, if any, to prevent a possible mis-determination of the CR transport parameters.

It is well established that phenomenological two-zone diffusion models of the galactic halo can very well reproduce cosmic-ray nuclear data and the observed antiproton flux. Here, we consider lepton propagation in such models and compute the expected galactic population of electrons, as well as the diffuse synchrotron emission that results from their interaction with galactic magnetic fields. We find models in agreement not only with cosmic ray data but also with radio surveys at essentially all frequencies. Requiring such a globally consistent description strongly disfavors very large ($L\gtrsim 15$ kpc) and, even stronger, small ($L\lesssim 1$ kpc) effective diffusive halo sizes. This has profound implications for, e.g., indirect dark matter searches.

The spectrum of an isotropic extragalactic $\gamma$-ray background (EGB) has been measured by the Fermi-LAT telescope at high latitudes. Two new models for the EGB are derived from the subtraction of unresolved point sources and extragalactic diffuse processes, which could explain from 30% to 70% of the Fermi-LAT EGB. Within the hypothesis that the two residual EGBs are entirely due to the annihilation of dark matter (DM) particles in the Galactic halo, we obtain $conservative$ upper limits on their annihilation cross section \sigmav. Severe bounds on a possible Sommerfeld enhancement of the annihilation cross section are set as well. Finally, would {\sigmav} be inversely proportional to the WIMP velocity, very severe limits are derived for the velocity-independent part of the annihilation cross section.

The spectrum of an isotropic extragalactic {\gamma}-ray background (EGB) has been measured by the Fermi-LAT telescope at high latitudes. Two new models for the EGB are derived from the subraction of unresolved point sources and extragalactic diffuse processes, which could explain from 30% to 70% of the Fermi-LAT EGB. Within the hypothesis that the two residual EGBs are entirely due to the annihilation of dark matter (DM) particles in the Galactic halo, we obtain stringent upper limits on their annihilation cross section. Severe bounds on a possible Sommerfeld enhancement of the annihilation cross section are set as well. Finally, we consider models for DM annihilation depending on the inverse of the velocity and associate the EGBs to photons arising from the annihilation of DM in primordial halos. For DM velocities in protohalos of v ~ 10^-8 c, the annihilation cross section is constrained down to 10^-33 cm^3 /s for DM masses below 100 GeV. Given our choices for the EGB and the minimal DM modelling, the derived upper bounds are claimed to be conservative.

A review is given of the main properties of the charged component of galactic cosmic rays, particles detected at Earth with an energy spanning from tens of MeV up to about 10^19 eV. After a short introduction to the topic and an historical overview, the properties of cosmic rays are discussed with respect to different energy ranges. The origin and the propagation of nuclei in the Galaxy are dealt with from a theoretical point of view. The mechanisms leading to the acceleration of nuclei by supernova remnants and to their subsequent diffusion through the inhomogeneities of the galactic magnetic field are discussed and some clue is given on the predictions and observations of fluxes of antimatter, both from astrophysical sources and from dark matter annihilation in the galactic halo. The experimental techniques and instrumentations employed for the detection of cosmic rays at Earth are described. Direct methods are viable up to 10^14 eV, by means of experiments flown on balloons or satellites, while above that energy, due to their very low flux, cosmic rays can be studied only indirectly by exploiting the particle cascades they produce in the atmosphere. The possible physical interpretation of the peculiar features observed in the energy spectrum of galactic cosmic rays, and in particular the so-called "knee" at about 4x10^15 eV, are discussed. A section is devoted to the region between about 10^18 and 10^19 eV, which is believed to host the transition between galactic and extragalactic cosmic rays. The conclusion gives some perspectives on the cosmic ray astrophysics field. Thanks to a wealth of different experiments, this research area is living a very flourishing era. The activity is exciting both from the theoretical and the instrumental sides, and its interconnection with astronomy, astrophysics and particle physics experiences non-stop growth.

Recent data from CREAM seem to confirm early suggestions that primary cosmic ray (CR) spectra at few TeV/nucleon are harder than in the 10-100 GeV range. Also, helium and heavier nuclei spectra appear systematically harder than the proton fluxes at corresponding energies. We note here that if the measurements reflect intrinsic features in the interstellar fluxes (as opposed to local effects) appreciable modifications are expected in the sub-TeV range for the secondary yields, such as antiprotons and diffuse gamma-rays. Presently, the ignorance on the origin of the features represents a systematic error in the extraction of astrophysical parameters as well as for background estimates for indirect dark matter searches. We find that the spectral modifications are appreciable above 100 GeV, and can be responsible for ~30% effects for antiprotons at energies close to 1 TeV or for gamma's at energies close to 300 GeV, compared to currently considered predictions based on simple extrapolation of input fluxes from low energy data. Alternatively, if the feature originates from local sources, uncorrelated spectral changes might show up in antiproton and high-energy gamma-rays, with the latter ones likely dependent from the line-of-sight.

The propagated fluxes of proton, helium, and heavier primary cosmic-ray species (up to Fe) are a means to indirectly access the source spectrum of cosmic rays. We check the compatibility of the primary fluxes with the transport parameters derived from the B/C analysis, but also if they bring further constraints. Proton data are well described in the simplest model defined by a power-law source spectrum and plain diffusion. They can also be accommodated by models with, e.g., convection and/or reacceleration. There is no need for breaks in the source spectral indices below $\sim 1$ TeV/n. Fits on the primary fluxes alone do not provide physical constraints on the transport parameters. If we let free the source spectrum $dQ/dE = q \beta^{\eta_S} {\cal R}^{-\alpha}$ and fix the diffusion coefficient $K(R)= K_0\beta^{\eta_T} {\cal R}^{\delta}$ such as to reproduce the B/C ratio, the MCMC analysis constrains the source spectral index $\alpha$ to be in the range $2.2-2.5$ for all primary species up to Fe, regardless of the value of the diffusion slope $\delta$. The $\eta_S$ low-energy shape of the source spectrum is degenerate with the low-energy shape $\eta_T$ of the diffusion coefficient: we find $\eta_S-\eta_T\approx 0$ for p and He data, but $\eta_S-\eta_T\approx 1$ for C to Fe primary species. This is consistent with the toy-model calculation in which the shape of the p/He and C/O to Fe/O data is reproduced if $\eta_S-\eta_T\approx 0-1$ (no need for different slopes $\alpha$). When plotted as a function of the kinetic energy per nucleon, the low-energy p/He ratio is shaped mostly by the modulation effect, whereas primary/O ratios are mostly shaped by their destruction rate.

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.

The CDMS Collaboration has presented its results for the final exposure of the CDMS II experiment and reports that two candidate events for dark matter would survive after application of the various discrimination and subtraction procedures inherent in their analysis. We show that a population of relic neutralinos, which was already proved to fit the DAMA/LIBRA data on the annual modulation effect, could naturally explain the two candidate CDMS II events, if these are actually due to a dark matter signal.

The indirect detection of particle dark matter (DM) is based on the search for anomalous components in cosmic rays (CRs) due to the annihilation of DM pairs in the galactic halo, on the top of the standard astrophysical production. These additional exotic components are potentially detectable at Earth as spectral distortions for the various cosmic radiations: $\chi + \chi \to q \bar{q}, W^+ W^-, ... \to \bar{p}, \bar{D}, e^+ \gamma and \nu's $. Detection of the DM annihilation products has motivated the spectacular development of several new experimental techniques. They range from detectors on ballons or in space for the study of antimatter and gamma-rays, to large area cosmic-ray and gamma-ray detecors on the ground to neutrino telescopes underground for the study of the neutrino component. In the following, we will discuss in detail the antimatter component of DM indirect searches, namely antiprotons, antideuterons, and positrons.

The predictions of our previous analyses about possible low-mass (lower than 50 GeV) relic neutralinos are discussed in the light of the most recent results from WIMP direct detection experiments. It is proved that these light neutralinos are quite compatible with the new annual-modulation data of the DAMA Collaboration; our theoretical predictions are also compared with the upper bounds of the CDMS and EDELWEISS Collaborations.

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.

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 confirm and extend the recent finding that the central surface density r_0*rho_0 galaxy dark matter halos, where r_0 and rho_0 are the halo core radius and central density, is nearly constant and independent of galaxy luminosity. Based on the co-added rotation curves of about 1000 spiral galaxies, mass models of individual dwarf irregular and spiral galaxies of late and early types with high-quality rotation curves and, galaxy-galaxy weak lensing signals from a sample of spiral and elliptical galaxies, we find that log(r_0*rho_0) = 2.15 +- 0.2, in units of log(Msol/pc^2). We also show that the observed kinematics of Local Group dwarf spheroidal galaxies are consistent with this value. Our results are obtained for galactic systems spanning over 14 magnitudes, belonging to different Hubble Types, and whose mass profiles have been determined by several independent methods. In the same objects, the approximate constancy of rho_0*r_0 is in sharp contrast to the systematical variations, by several orders of magnitude, of galaxy properties, including rho_0 and central stellar surface density.

A new calculation of the $\bar{p}/p$ ratio in cosmic rays is compared to the recent PAMELA data. The good match up to 100 GeV allows to set constraints on exotic contributions from thermal WIMP dark matter candidates. We derive stringent limits on possible enhancements of the WIMP \pbar flux: a $m_{\rm WIMP}$=100 GeV (1 TeV) signal cannot be increased by more than a factor 6 (40) without overrunning PAMELA data. Annihilation through the $W^+W^-$ channel is also inspected and cross-checked with $e^+/(e^-+e^+)$ data. This scenario is strongly disfavored as it fails to simultaneously reproduce positron and antiproton measurements.

The most recent results from direct searches for dark matter particles in the galactic halo are examined in terms of an effective Minimal Supersymmetric extension of the Standard Model at the electroweak scale without gaugino masses unification. We show that the annual modulation effect at 8.2 $\sigma$ C.L. recently presented by the DAMA Collaboration, as the result of a combined analysis of the DAMA/NaI and the DAMA/LIBRA experiments for a total exposure of 0.82 ton yr, fits remarkably well with what expected for relic neutralinos for a wide variety of WIMP distribution functions. Bounds derivable from other measurements of direct searches for dark matter particles are analyzed. We stress the role played by the uncertainties affecting the neutralino--quark couplings arising from the involved hadronic quantities. We also examine how present data on cosmic antiprotons can help in constraining the neutralino configurations selected by the DAMA effect, in connection with the values of the astrophysical parameters. Perspectives for measurement of antideuterons possibly produced in the galactic halo by self--annihilation of neutralinos belonging to the DAMA configurations are examined. Finally, we discuss how findings at LHC would impact on these issues.

Antideuterons are among the most promising galactic cosmic ray-related targets for dark matter indirect detection. Currently only upper limits exist on the flux, but the development of new experiments, such as GAPS and AMS-02, provides exciting perspectives for a positive measurement in the near future. In this Paper, we present a novel and updated calculation of both the secondary and primary antideuteron fluxes. We employ a two-zone diffusion model which successfully reproduces cosmic-ray nuclear data and the observed antiproton flux. We review the nuclear and astrophysical uncertainties and provide an up to date secondary (i.e. background) antideuteron flux. The primary (i.e. signal) contribution is calculated for generic WIMPs annihilating in the galactic halo: we explicitly consider and quantify the various sources of uncertainty in the theoretical evaluations. Propagation uncertainties, as is the case of antiprotons, are sizeable. Nevertheless, antideuterons offer an exciting target for indirect dark matter detection for low and intermediate mass WIMP dark matter. We then show the reaching capabilities of the future experiments for neutralino dark matter in a variety of supersymmetric models.

Two recent papers reconsider the possibility that the excess of microwave emission from a region within $\sim 20^0$ of the galactic center (the {\it WMAP haze}), measured by WMAP, can be due to the synchrotron emission originated by neutralino self-annihilation; on the basis of this possible occurrence, also upper bounds on the neutralino self-annihilation cross--section are suggested. In the present note, we show that in the common case of thermal WIMPs in a standard cosmological model, when the rescaling of the galactic WIMP density is duly taken into account for subdominant WIMPs, the upper bound applicable generically to {\it any} signal due to self-conjugate WIMPs is more stringent than the ones obtained from analysis of the WMAP haze. We also argue that an experimental upper bound, which can compete with our generic upper limit, can rather be derived from measurements of cosmic antiproton fluxes, for some values of the parameters of the astrophysical propagation model. Finally, we comment on the possible impact of our generic upper bound on the interpretation of the WMAP haze in terms of thermal neutralinos in a standard cosmological scheme.

Indirect detection signals from dark matter annihilation are studied in the positron channel. We discuss in detail the positron propagation inside the galactic medium: we present novel solutions of the diffusion and propagation equations and we focus on the determination of the astrophysical uncertainties which affect the positron dark matter signal. We find dark matter scenarios and propagation models that nicely fit existing data on the positron fraction. Finally, we present predictions both on the positron fraction and on the flux for already running or planned space experiments, concluding that they have the potential to discriminate a possible signal from the background and, in some cases, to distinguish among different astrophysical propagation models.

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.

The DAMA Collaboration has recently analyzed its data of the extensive WIMP direct search (DAMA/NaI) which detected an annual modulation, by taking into account the channelling effect which occurs when an ion traverses a detector with a crystalline structure. Among possible implications, this Collaboration has considered the case of a coherent WIMP-nucleus interaction and then derived the form of the annual modulation region in the plane of the WIMP-nucleon cross section versus the WIMP mass, using a specific modelling for the channelling effect. In the present paper we first show that light neutralinos fit the annual modulation region also when channelling is taken into account. To discuss the connection with indirect signals consisting in galactic antimatter, in our analysis we pick up a specific galactic model, the cored isothermal-sphere. In this scheme we determine the sets of supersymmetric models selected by the annual modulation regions and then prove that these sets are compatible with the available data on galactic antiprotons. We comment on implications when other galactic distribution functions are employed. Finally, we show that future measurements on galactic antiprotons and antideuterons will be able to shed further light on the populations of light neutralinos singled out by the annual modulation data.

Indirect detection signals from dark matter annihilation are studied in the positron channel. We discuss in detail the positron propagation inside the galactic medium: we present novel solutions of the diffusion and propagation equations and we focus on the determination of the astrophysical uncertainties which affect the positron dark matter signal. We show that, especially in the low energy tail of the positron spectra at Earth, the uncertainty is sizeable and we quantify the effect. Comparison of our predictions with current available and foreseen experimental data are derived.

We calculate the antideuteron flux expected from dark matter annihilation in the galactic halo. The propagation is treated in a full 2-D propagation model consistent with the results obtained from the propagation of B/C and other galactic species. We discuss the potentials of this indirect dark matter detection means, with special emphasis on the possible sources of uncertainties affecting future measurements

The possibility to use $\gamma$--ray data from the Galactic Center (GC) to constrain the cosmological evolution of the Universe in a phase prior to primordial nucleosyntesis, namely around the time of cold dark matter (CDM) decoupling, is analyzed. The basic idea is that in a modified cosmological scenario, where the Hubble expansion rate is enhanced with respect to the standard case, the CDM decoupling is anticipated and the relic abundance of a given dark matter (DM) candidate enhanced. This implies that the present amount of CDM in the Universe may be explained by a Weakly Interacting Massive Particle (WIMP) which possesses annihilation cross section (much) larger than in standard cosmology. This enhanced annihilation implies larger fluxes of indirect detection signals of CDM. We show that the HESS measurements can set bounds for WIMPs heavier than a few hundreds of GeV, depending on the actual DM halo profile. These results are complementary to those obtained in a previous analysis based on cosmic antiprotons. For a Moore DM profile, $\gamma$--ray data limit the maximal Hubble rate enhancement to be below a factor of 100. Moreover, a WIMP heavier than 1 TeV is not compatible with a cosmological scenario with enhanced expansion rate prior to Big Bang Nucleosynthesis (BBN). Less steep DM profiles provide less stringent bounds, depending of the cosmological scenario.

New upper bounds on direct detection rates have recently been presented by a number of experimental collaborations working on searches for WIMPs. In this paper we analyze how the constraints on relic neutralinos which can be derived from these results is affected by the uncertainties in the distribution function of WIMPs in the halo. Various different categories of velocity distribution functions are considered, and the ensuing implications for supersymmetric configurations derived. We conservatively conclude that current experimental data do not constrain neutralinos of small mass (below 50 GeV).

We analyze how the measurements of the low-energy spectrum of cosmic antiprotons can provide information on relic neutralinos. The analysis is focused on the light neutralinos which emerge in supersymmetric schemes where gaugino-mass unification is not assumed. We determine which ranges of the astrophysical parameters already imply stringent constraints on the supersymmetric configurations and those ranges which make the antiproton flux sensitive to the primary component generated by the neutralino self-annihilation. Our results are derived from some general properties of the antiproton flux proved to be valid for a generic cold dark matter candidate.

We discuss the potentials of several experimental configurations dedicated to direct measurements of charged cosmic ray (CR) nuclei at energies $\gsim$ 100 GeV/n. Within a two-zone propagation model for stable CRs, we calculate light primary and secondary nuclei fluxes for different diffusion and acceleration schemes. We show that the new detectors exploiting the long and ultra long duration balloon flights could determine the diffusion coefficient power index $\delta$ through the measurement of the boron-to-carbon ratio with an uncertainty of about 10-15 %, if systematic errors are low enough. Only space-based or satellite detectors will be able to determine $\delta$ with very high accuracy even in the case of important systematic errors, thanks to the higher energy reach and the less severe limitations in the exposure. We show that no uncertainties other than those on $\delta$ affect the determination of the acceleration slope $\alpha$, so that measures of light primary nuclei, such as the carbon one, performed with the same experiments, will provide valuable information on the acceleration mechanisms.

Extra dimensions offer new ways to address long-standing problems in beyond the standard model particle physics. In some classes of extra-dimensional models, the lightest Kaluza-Klein particle is a viable dark matter candidate. In this work, we study indirect detection of Kaluza-Klein dark matter via its annihilation into antiprotons. We use a sophisticated galactic cosmic ray diffusion model whose parameters are fully constrained by an extensive set of experimental data. We discuss how fluxes of cosmic antiprotons can be used to exclude low Kaluza-Klein masses.

We investigate in detail the mass distribution obtained by means of high resolution rotation curves of 25 galaxies of different morphological types. The dark matter contribution to the circular rotation velocity is well-described by resorting to a dark component whose density shows an inner core, i.e. a central constant density region. We find a very strong correlation between the core radius size $R_C$ and the stellar exponential scale length $R_D$: $R_C \simeq 13 (\frac {R_D} {5 {\rm kpc}})^{1.05} {\rm kpc} $, and between $R_C$ and the galaxy dynamical mass at this distance, $M_{dyn}(R_C)$. These relationships would not be expected if the core radii were the product of an incorrect decomposition procedure, or the biased result of wrong or misunderstood observational data. The very strong correlation between the dark and luminous scale lengths found here seems to hold also for different Hubble types and opens new scenarios for the nature of the dark matter in galaxies.

We examine indirect signals produced by neutralino self-annihilations, in the galactic halo or inside celestial bodies, in the frame of an effective MSSM model without gaugino-mass unification at a grand unification scale. We compare our theoretical predictions with current experimental data of gamma-rays and antiprotons in space and of upgoing muons at neutrino telescopes. Results are presented for a wide range of the neutralino mass, though our discussions are focused on light neutralinos. We find that only the antiproton signal is potentially able to set constraints on very low-mass neutralinos, below 20 GeV. The gamma-ray signal, both from the galactic center and from high galactic latitudes, requires significantly steep profiles or substantial clumpiness in order to reach detectable levels. The up-going muon signal is largely below experimental sensitivities for the neutrino flux coming from the Sun; for the flux from the Earth an improvement of about one order of magnitude in experimental sensitivities (with a low energy threshold) can make accessible neutralino masses close to O, Si and Mg nuclei masses, for which resonant capture is operative.

We calculate the antiproton flux due to relic neutralino annihilations, in a two-dimensional diffusion model compatible with stable and radioactive cosmic ray nuclei. We find that the uncertainty in the primary flux induced by the propagation parameters alone is about two orders of magnitude at low energies, and it is mainly determined by the lack of knowledge on the thickness of the diffusive halo. On the contrary, different dark matter density profiles do not significantly alter the flux: a NFW distribution produces fluxes which are at most 20% higher than an isothermal sphere. The most conservative choice for propagation parameters and dark matter distribution normalization, together with current data on cosmic antiprotons, cannot lead to any definitive constraint on the supersymmetric parameter space, neither in a low-energy effective MSSM, or in a minimal SUGRA scheme. However, if the best choice for propagation parameters - corresponding to a diffusive halo of L=4 kpc - is adopted, some supersymmetric configurations with the neutralino mass of about 100 GeV should be considered as excluded. An enhancement flux factor - due for instance to a clumpy dark halo or to a higher local dark matter density - would imply a more severe cut on the supersymmetric parameters.

In the frame of an effective MSSM model without gaugino-mass unification at a grand unification (GUT) scale, we set a lower bound on the neutralino mass based on the new WMAP data on $\Omega_{CDM}$. Our lower bound, $m_{\chi} \gsim 6$ GeV, leaves much room for relic neutralinos significantly lighter than those commonly considered ($m_{\chi} \gsim 50$ GeV). We prove that these light neutralinos may produce measurable effects in WIMP direct detection experiments of low energy threshold and large exposure.

We present the results on the source spectrum function for primary nuclei in galactic cosmic rays, where two distinct energy dependences are used for the source spectra. We discuss the evolution of the goodness of fit to B/C data with the propagation parameters and also show that the results are not much affected by a different choice for the diffusion scheme. We apply the constraints on the diffusion scheme as derived from stable nuclei to calculate the propagation of beta-radioactive isotopes. The diffusion model is refined to properly take into account the effect of a local bubble in the interstellar medium. Our calculations are compared to existing data, which prefer the local bubble description instead of the homogeneus scheme.

We use the constraints on the diffusion parameters as obtained with stable nuclei to calculate the cosmic antiproton flux from annihilating relic neutralinos. We discuss the relevance of each characteristic parameter, describing our two dimension diffusion model, on the flux of antiprotons produced in the dark halo of our Galaxy. We estimate a two orders of magnitude uncertainty on the flux due to the unknowledge of the propagation parameters. A conservative and systematic evaluation of the flux in the supersymmetric parameter space is done in order to exclude configurations providing a total (secondary plus primary) flux in excess with observations. We also study the effect on the flux induced by modifications in the distribution of cold dark matter in the Galaxy.

In most cosmological models, primordial black holes (PBHs) should have formed in the early Universe. Their Hawking evaporation into particles could eventually lead to the formation of antideuterium nuclei. This paper is devoted to a first computation of this antideuteron flux. The production of these antinuclei is studied with a simple coalescence scheme, and their propagation in the Galaxy is treated with a well-constrained diffusion model. We compare the resulting primary flux to the secondary background, due to the spallation of protons on the interstellar matter. Antideuterons are shown to be a very sensitive probe for primordial black holes in our Galaxy. The next generation of experiments should allow investigators to significantly improve the current upper limit, nor even provide the first evidence of the existence of evaporating black holes.

Cosmic Ray nuclei in the energy range 100 MeV/nuc - 100 GeV/nuc provide crucial information about the physical properties of the Galaxy. They can also be used to answer questions related to astroparticle physics. This paper reviews the results obtained in this direction, with a strong bias towards the work done by the authors at {\sc lapth}, {\sc isn} and {\sc iap}. The propagation of these nuclei is studied quantitatively in the framework of a semi-analytical two-zone diffusion model taking into account the effect of galactic wind, diffuse reacceleration and energy losses. The parameters of this model are severely constrained by an analysis of the observed B/C ratio. These constraints are then used to study other species such as radioactive species and light antinuclei. Finally, we focus on the astroparticle subject and we study the flux of antiprotons and antideuterons that might be due to neutralino annihilations or primordial black hole evaporation. The question of the spatial origin of all these species is also addressed.

In a previous study (Maurin et al., 2001), we explored the set of parameters describing diffusive propagation of cosmic rays (galactic convection, reacceleration, halo thickness, spectral index and normalization of the diffusion coefficient), and we identified those giving a good fit to the measured B/C ratio. This study is now extended to take into account a sixth free parameter, namely the spectral index of sources. We use an updated version of our code where the reacceleration term comes from standard minimal reacceleration models. The goal of this paper is to present a general view of the evolution of the goodness of fit to B/C data with the propagation parameters. In particular, we find that, unlike the well accepted picture, and in accordance with our previous study, a Kolmogorov-like power spectrum for diffusion is strongly disfavored. Rather, the $\chi^2$ analysis points towards $\delta\gtrsim 0.7$ along with source spectra index $\lesssim 2.0$. Two distinct energy dependences are used for the source spectra: the usual power-law in rigidity and a law modified at low energy, the second choice being only slightly preferred. We also show that the results are not much affected by a different choice for the diffusion scheme. Finally, we compare our findings to recent works, using other propagation models. This study will be further refined in a companion paper, focusing on the fluxes of cosmic ray nuclei.

Primordial black holes (PBHs) have motivated many studies since it was shown that they should evaporate and produce all kinds of particles. Recent experimental measurements of cosmic rays with great accuracy, theoretical investigations on the possible formation mechanisms and detailed evaporation processes have revived the interest in such astrophysical objects. This article aims at using the latest developments on antiproton propagation models together with new data from BESS, CAPRICE and AMS experiments to constrain the local amount of PBH dark matter. Depending on the diffusion halo parameters and on the details of emission mechanism, we derive an average upper limit of the order of rho < 1.7E-33 g cm^-3.

We discuss the implications of a new determination of the pion-nucleon sigma term for the evaluation of the neutralino-nucleon cross-section, and, in turn, for detection rates for relic neutralinos in WIMP direct experiments and in some of the indirect searches. It is shown that the new range for the pion-nucleon sigma term, taken at its face value, favours values of the neutralino-nucleon cross-sections which are sizeably larger than some of the current estimates. Implications for neutralino cosmological properties are derived.

In the present paper, we extend the analysis of Maurin et al. (2001) and Donato et al. (2001) to the $\beta$-radioactive nuclei $^{10}$Be, $^{26}$Al and $^{36}$Cl. These species are be shown to be particularly sensitive to the properties of the local interstellar medium (LISM). As studies of the LISM suggest that we live in an underdense bubble of extent $r_{hole} \sim 50 - 200 \unit{pc}$, this local feature must be taken into account. We present a modified version of our diffusion model which describes the underdensity as a hole in the galactic disc. It is found that the presence of the bubble leads to a decrease in the radioactive fluxes which can be approximated by a simple factor $\exp(-r_{hole}/l_{rad})$ where $l_{rad}=\sqrt{K \gamma \tau_0}$ is the typical distance travelled by a radioactive nucleus before it decays. We find that each of the radioactive nuclei independently point towards a bubble radius $\lesssim 100 \unit{pc}$. If these nuclei are considered simultaneously, only models with a bubble radius $r_{hole} \sim 60 - 100 \unit{pc}$ are marginally consistent with data. In particular, the standard case $r_{hole}=0 \unit{pc}$ is disfavoured. Our main concern is about the consistency of the currently available data, especially $^{26}$Al/$^{27}$Al.

Cosmic ray antiprotons provide an important probe for the study of the galactic Dark Matter, as they could be produced by exotic sources. On the other hand, antiprotons are anyway produced by standard nuclear reactions of cosmic ray nuclei on interstellar matter. This process is responsible for a background flux that must be carefully determined to estimate the detectability of an hypothetical exotic signal. Estimates of this background suffer from potential uncertainties of various origins. The propagation of cosmic antiprotons depends on several physical characteristics of the Galaxy which are poorly known. Antiprotons are created from cosmic protons and helium nuclei whose fluxes were not measured with great accuracy until very recently. Calculations of antiproton fluxes make use of nuclear physics models with their own shortcomings and uncertainties. The goal of this paper is to give a new evaluation of the cosmic antiproton flux along with the associated uncertainties. The propagation parameters were tightly constrained in Maurin et al. 2001 by an analysis of cosmic ray nuclei data in the framework of a two-zone diffusion model and we apply these parameters to the propagation of antiprotons. We use the recently published data on proton and helion fluxes, and we find that this particular source of uncertainty has become negligible. The Monte Carlo program DTUNUC was used to carefully examine nuclear reactions. We find that all the cosmic antiproton fluxes naturally coming out of the calculation are fully compatible with experimental data. Uncertainties in this flux have been strongly reduced. Those related to propagation are less than 25%. All other possible sources of uncertainty have also been studied.

After a short introduction on particle candidates for dark matter within possible extensions of the standard model, we concentrate on Weakly Interacting Massive Particles, and on one of their most interesting physical realizations: the neutralino. We analyze how detectability of relic neutralinos by direct and indirect means is related to their local and cosmological densities; we use simple general arguments to discusss different scenarios where relic neutralinos make up the dominant bulk of dark matter or only a small fraction of it. Our general arguments are further corroborated by specific numerical results. We show to which extent the present experiments of direct searches for WIMPs, when interpreted in terms of relic neutralinos, probe interesting regions of the supersymmetric parameter space. Our analysis is performed in a number of different supersymmetric schemes.

We discuss to which extent the present experiments of direct search for WIMPs, when interpreted in terms of relic neutralinos, probe interesting regions of the supersymmetric parameter space, which are also being progressively explored at accelerators. Our analysis is performed in a number of different supersymmetric schemes. We derive the relevant neutralino cosmological properties, locally and on the average in the universe. We prove that part of the susy configurations probed by current WIMP experiments entail relic neutralinos of cosmological interest. The main astrophysical and particle physics uncertainties, relevant for a proper comparison between theory and experimental data, are stressed and taken into account.

Cosmic ray nuclei fluxes are expected to be measured with high precision in the near future. For instance, high quality data on the antiproton component could give important clues about the nature of the astronomical dark matter. A very good understanding of the different aspects of cosmic ray propagation is therefore necessary. In this paper, we use cosmic ray nuclei data to give constraints on the diffusion parameters. Propagation is studied with semi-analytical solutions of a diffusion model, and we give new analytical solutions for radioactively produced species. Our model includes convection and reacceleration as well as the standard energy losses. We perform a $\chi^2$ analysis over B/C data for a large number of configurations obtained by varying the relevant parameters of the diffusion model. A very good agreement with B/C data arises for a number of configurations, all of which are compatible with sub-Fe/Fe data. Different source spectra $Q(E)$ and diffusion coefficients $K(E)$ have been tried, but for both parameters only one form gives a good fit. Another important result is that models without convection or without reacceleration are excluded. We find that the various parameters, i.e. the diffusion coefficient normalisation $K_0$ and spectral index $\delta$, the halo thickness $L$, the Alfv\'en velocity $V_a$, and the convection velocity $V_c$ are strongly correlated. We obtain limits on the spectral index $\delta$ of thediffusion coefficient, and in particular we exclude a Kolmogorov spectrum ($\delta = 1/3$).

We propose low-energy antideuterons in cosmic rays as a new possible signature for indirect detection of supersymmetric dark matter. Since the energy spectrum of the antiproton secondary component is still spoilt by considerable theoretical uncertainties, looking for low-energy antideuterons seems a plausible alternative. We apply our calculation to the AMS experiment, when mounted on the International Spatial Station. If a few low-energy antideuterons will be discovered by AMS, this should be seriously taken as a clue for the existence of relic, massive neutralinos in the dark halo of our Galaxy.

We analyze the annual-modulation effect, measured by the DAMA Collaboration with the new implementation of a further two-years running, in the context of a possible interpretation in terms of relic neutralinos. We impose over the set of supersymmetric configurations, selected by the annual-modulation data, the constraints derived from WIMP indirect measurements, and discuss the features of the ensuing relic neutralinos. We critically discuss the sources of the main theoretical uncertainties in the analysis of event rates for direct and indirect WIMP searches.

We discuss the effect induced on the neutralino-nucleon cross-section by the present uncertainties in the values of the quark masses and of the quark scalar densities in the nucleon. We examine the implications of this aspect on the determination of the neutralino cosmological properties, as derived from measurements of WIMP direct detection. We show that, within current theoretical uncertainties, the DAMA annual modulation data are compatible with a neutralino as a major dark matter component, to an extent which is even larger than the one previously derived. We also comment on implications of the mentioned uncertainties for experiments of indirect dark matter detection.

Once the energy spectrum of the secondary component is well understood, measurements of the antiproton cosmic-ray flux at the Earth will be a powerful way to indirectly probe for the existence of supersymmetric relics in the galactic halo. Unfortunately, it is still spoilt by considerable theoretical uncertainties. As shown in this work, searches for low-energy antideuterons appear in the mean time as a plausible alternative, worth being explored. Above a few GeV/n, a dozen spallation antideuterons should be collected by the future AMS experiment on board ISSA. For energies less than about 3 GeV/n, the antideuteron spallation component becomes negligible and may be supplanted by a potential supersymmetric signal. If a few low-energy antideuterons are discovered, this should be seriously taken as a clue for the existence of massive neutralinos in the Milky Way.

The original annual-modulation region, singled out by the DAMA/NaI experiment for direct detection of WIMPs, is extended by taking into account the uncertainties in the galactic astrophysical velocities. Also the effect due to a possible bulk rotation for the dark matter halo is considered. We find that the range for the WIMP mass becomes 30 GeV < m_chi < 130 GeV at 1-sigma C.L. with a further extension in the upper bound, when a possible bulk rotation of the dark matter halo is taken into account. We show that the DAMA results, when interpreted in the framework of the Minimal Supersymmetric extension of the Standard Model, are consistent with a relic neutralino as a dominant component of cold dark matter (on the average in our universe and in our galactic halo). It is also discussed the discovery potential for the relevant supersymmetric configurations at accelerators of present generation.

We report the main results of a paper where recent data of low-energy cosmic-ray $\bar p$ spectrum have been analyzed in terms of newly calculated fluxes for secondary antiprotons and for a possible contribution of an exotic signal due to neutralino annihilation in the galactic halo. We also present the results of a paper in which we have proved that a sizeable fraction of the supersymmetric neutralino configurations, singled out by the DAMA/NaI data on a possible annual modulation effect in WIMP direct search, may provide detectable cosmic-ray antiproton signals.

In previous papers we showed that the data of the DAMA/NaI experiment for WIMP direct detection, which indicate a possible annual modulation effect, are widely compatible with an interpretation in terms of a relic neutralino as the major component of dark matter in the Universe. In the present note we discuss the detectability of the relevant supersymmetric neutralino configurations by two of the most promising methods of indirect search for relic particles: measurement of cosmic-ray antiprotons and measurement of neutrino fluxes from Earth and Sun.

Recent results of the DAMA/NaI experiment for WIMP direct detection point to a possible annual modulation effect in the detection rate. We show that these results, when interpreted in terms of a relic neutralino, are compatible with supergravity models. Together with the universal SUGRA scheme, we also consider SUGRA models where the unification condition in the Higgs mass parameters at GUT scale is relaxed.

We demonstrate that the further indication of a possible annual modulation effect, singled out by the DAMA/NaI experiment for WIMP direct detection, is widely compatible with an interpretation in terms of a relic neutralino as the major component of dark matter in the Universe. We discuss the supersymmetric features of this neutralino in the Minimal Supersymmetric extension of the Standard Model (MSSM) and their implications for searches at accelerators.

We analyze the data of low-energy cosmic-ray antiproton spectrum, recently published by the BESS Collaboration, in terms of newly calculated fluxes for secondary antiprotons and for a possible contribution of an exotic signal due to neutralino annihilation in the galactic halo. We single out the relevant supersymmetric configurations and discuss their explorability with experiments of direct search for particle dark matter and at accelerators. We discuss how future measurements with the Alpha Magnetic Spectrometer (AMS) on the Shuttle flight may disentangle the possible neutralino-induced contribution from the secondary one.

The effects of a possible rotation of the galactic dark halo on the calculation of the direct detection rates for particle dark matter are analyzed, with special attention to the extraction of the upper limits on the WIMP--nucleon scalar cross section from the experimental data. We employ a model of dark halo rotation which describes the maximal possible effects. For WIMP masses above 50 GeV, the upper limit exclusion plot is modified by less than a factor of two when rotation is included. For lighter masses the effect can be stronger, suggesting the necessity to develop specific models of halo rotation in order to provide more accurate conclusions.

In this note we present an extension of our previous analysis of some preliminary experimental results of the DAMA/NaI Collaboration which might be indicative of a yearly modulation effect. Here we present a direct way for obtaining from the experimental data the relevant cosmological implications for relic neutralinos. We find that some of the configurations singled out by the DAMA/NaI results would have cosmological properties compatible with a neutralino as a dominant component of cold dark matter (on the average in the Universe and in our galactic halo).

It is shown that the new data of the DAMA/NaI detector, combined with recent calculations of nuclear structure functions, provide the most stringent upper bounds on the WIMP-nucleon cross-section both for spin-dependent and for scalar interactions. It is proved that, for the first time, data from an experiment of direct search for WIMPs allow the exploration of a region of the supersymmetric parameter space not yet fully investigated at accelerators.