A critical review of the standard paradigm for the origin of Galactic cosmic rays is presented. Recent measurements of local and far-away cosmic rays reveal unexpected behaviours, which challenge the commonly accepted scenario. These recent findings are discussed, together with long-standing open issues. Despite the progress made thanks to ever-improving observational techniques and theoretical investigations, at present our understanding of the origin and of the behaviour of cosmic rays remains incomplete. We believe it is still unclear whether a modification of the standard paradigm, or rather a radical change of the paradigm itself is needed in order to interpret all the available data on cosmic rays within a self-consistent scenario.

Until fairly recently, it was widely accepted that local cosmic ray spectra were largely featureless power laws, containing limited information on their acceleration and transport. This viewpoint is currently being revised in the light of evidence for a variety of spectral breaks in the fluxes of cosmic ray nuclei. Here, we focus on cosmic ray electrons and positrons which at the highest energies must be of local origin due to strong radiative losses. We consider a pure diffusion model for their Galactic transport and determine its free parameters by fitting data in a wide energy range: measurements of the interstellar spectrum by Voyager at MeV energies, radio synchrotron data (sensitive to GeV electrons and positrons) and local observations by AMS up to ~ 1 TeV. For the first time, we also model the time-dependent fluxes of cosmic ray electrons and positrons at GeV energies recently presented by AMS, treating solar modulation in a simple extension of the widely used force-field approximation. We are able to reproduce all the available measurements to date. Our model of the interstellar spectrum of cosmic ray electrons and positrons requires the presence of a number of spectral breaks, both in the source spectra and the diffusion coefficients. While we remain agnostic as to the origin of these spectral breaks, their presence will inform future models of the microphysics of cosmic ray acceleration and transport.

Understanding the isotopic composition of cosmic rays (CRs) observed near Earth represents a milestone towards the identification of their origin. Local fluxes contain all the known stable and long-lived isotopes, reflecting the complex history of primaries and secondaries as they traverse the interstellar medium. For that reason, a numerical code which aims at describing the CR transport in the Galaxy must unavoidably rely on accurate modelling of the production of secondary particles. In this work we provide a detailed description of the nuclear cross sections and decay network as implemented in the forthcoming release of the galactic propagation code DRAGON2. We present the secondary production models implemented in the code and we apply the different prescriptions to compute quantities of interest to interpret local CR fluxes (e.g., nuclear fragmentation timescales, secondary and tertiary source terms). In particular, we develop a nuclear secondary production model aimed at accurately computing the light secondary fluxes (namely: Li, Be, B) above 1 GeV/n. This result is achieved by fitting existing empirical or semi-empirical formalisms to a large sample of measurements in the energy range 100 MeV/n to 100 GeV/n and by considering the contribution of the most relevant decaying isotopes up to iron. Concerning secondary antiparticles (positrons and antiprotons), we describe a collection of models taken from the literature, and provide a detailed quantitative comparison.

A crucial process in Galactic cosmic-ray (CR) transport is the spatial diffusion due to the interaction with the interstellar turbulent magnetic field. Usually, CR diffusion is assumed to be uniform and isotropic all across the Galaxy. However, this picture is clearly inaccurate: Several data-driven and theoretical arguments, as well as dedicated numerical simulations, show that diffusion exhibits highly anisotropic properties with respect to the direction of a background (ordered) magnetic field (i.e., parallel or perpendicular to it). In this paper we focus on a recently discovered anomaly in the hadronic CR spectrum inferred by the Fermi-LAT gamma-ray data at different positions in the Galaxy, i.e. the progressive hardening of the proton slope at low Galactocentric radii. We propose the idea that this feature can be interpreted as a signature of anisotropic diffusion in the complex Galactic magnetic field: In particular, the harder slope in the inner Galaxy is due, in our scenario, to the parallel diffusive escape along the poloidal component of the large-scale, regular, magnetic field. We implement this idea in a numerical framework, based on the DRAGON code, and perform detailed numerical tests on the accuracy of our setup. We discuss how the effect proposed depends on the relevant free parameters involved. Based on low-energy extrapolation of the few focused numerical simulations aimed at determining the scalings of the anisotropic diffusion coefficients, we finally present a set of plausible models that reproduce the behavior of the CR proton slopes inferred by gamma-ray data.

Before being detected at Earth, charged cosmic rays propagate across the Solar System and undergo interactions with the turbulent solar wind and with the heliospheric magnetic field. As a result, they are subject to a series of processes that include diffusion, convection, energy losses and drifts, which significantly affect the shape and the intensity of the cosmic-ray fluxes at low energies. Here we illustrate how all these mechanisms can be realistically modelled with HelioProp, our public tool designed to treat cosmic-ray transport through the heliosphere in a charge-dependent way. We present a detailed description of the features of the code and we illustrate in a quantitative way the effects that the propagation in the heliosphere can have on the different cosmic-ray species with a particular emphasis on the antiparticle channels relevant for dark matter indirect detection.

The search for antideuterons in cosmic rays has been proposed as a promising channel for dark matter indirect detection, especially for dark matter particles with a low or intermediate mass. With the current operational phase of the AMS-02 experiment and the ongoing development of a future dedicated experiment, the General Antiparticle Spectrometer (GAPS), there are exciting prospects for a dark matter detection in the near future. In this paper we develop a detailed and complete re-analysis of the cosmic-ray antideuteron signal, by discussing the main relevant issues related to antideuteron production and propagation through the interstellar medium and the heliosphere. In particular, we first critically revisit the coalescence mechanism for antideuteron production in dark matter annihilation processes. Then, since antideuteron searches have their best prospects of detection at low kinetic energies where the effect of the solar wind and magnetic field are most relevant, we address the impact of solar modulation modeling on the antideuteron flux at the Earth by developing a full numerical 4D solution of cosmic rays transport in the heliosphere. We finally use these improved predictions to provide updated estimates of the reaching capabilities for AMS-02 and GAPS, compatible with the current constraints imposed by the antiprotons measurements of PAMELA. After the antiproton bound is applied, prospects of detection of up to about 15 events in GAPS LDB+ and AMS-02 missions are found, depending on the dark matter mass, annihilation rate and production channel from one side, and on the coalescence process, galactic and solar transport parameters on the other.

We present version 2 of the DRAGON code designed for computing realistic predictions of the CR densities in the Galaxy. The code numerically solves the interstellar CR transport equation (including inhomogeneous and anisotropic diffusion, either in space and momentum, advective transport and energy losses), under realistic conditions. The new version includes an updated numerical solver and several models for the astrophysical ingredients involved in the transport equation. Improvements in the accuracy of the numerical solution are proved against analytical solutions and in reference diffusion scenarios. The novel features implemented in the code allow to simulate the diverse scenarios proposed to reproduce the most recent measurements of local and diffuse CR fluxes, going beyond the limitations of the homogeneous galactic transport paradigm. To this end, several applications using DRAGON2 are presented as well. The new version facilitates the users to include their own physical models by means of a modular C++ structure.

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.

Antibaryons are produced in our Galaxy in collisions of high energy cosmic rays with the interstellar medium and in old supernova remnants, and possibly, in exotic sources such as primordial black hole evaporation or dark matter annihilations and decays. The search for signals from exotic sources in antiproton data is hampered by large backgrounds from spallation which, within theoretical errors, can solely account for the current data. Due to the higher energy threshold for antideuteron production, which translates into a suppression of the low energy flux from spallations, antideuteron searches have been proposed as a probe for exotic sources. We perform in this paper a comprehensive analysis of the antideuteron fluxes at the Earth expected from known and hypothetical sources in our Galaxy, and we calculate their maximal values consistent with current antiproton data from AMS-02. We find that supernova remnants generate a negligible flux, whereas primordial black hole evaporation and dark matter annihilations or decays may dominate the total flux at low energies. On the other hand, we find that the detection of cosmic antideuterons would require, for the scenarios studied in this paper and assuming optimistic values of the coalescence momentum and solar modulation, an increase of the experimental sensitivity compared to ongoing and planned instruments by at least a factor of 2. Finally, we briefly comment on the prospects for antihelium-3 detection.

Galactic Dark Matter (DM) annihilations can produce cosmic-ray anti-nuclei via the nuclear coalescence of the anti-protons and anti-neutrons originated directly from the annihilation process. Since anti-deuterons have been shown to offer a distinctive DM signal, with potentially good prospects of detection in large portions of the DM-particle parameter space, we explore here the production of heavier anti-nuclei, specifically anti-helium. Even more than for anti-deuterons, the DM-produced anti-He flux can be mostly prominent over the astrophysical anti-He background at low kinetic energies, typically below 3-5 GeV/n. However, the larger number of anti-nucleons involved in the formation process makes the anti-He flux extremely small. We therefore explore, for a few DM benchmark cases, whether the yield is sufficient to allow for anti-He detection in current-generation experiments, such as AMS-02. We account for the uncertainties due to the propagation in the Galaxy and to the uncertain details of the coalescence process, and we consider the constraints already imposed by anti-proton searches. We find that only for very optimistic configurations it might be possible to achieve detection with current generation detectors. We estimate that, in more realistic configurations, an increase in experimental sensitivity at low-kinetic energies of about a factor of 500-1000 would allow to start probing DM through the rare cosmic anti-He production.

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.

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.

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.

We present here a quantitative analysis of the recent AMS-02 data with the purpose of investigating the interplay between astrophysical sources and Dark Matter in their interpretation. First, we show that AMS-02 leptonic measurements are in a remarkably good agreement with the hypothesis that all electrons and positrons are the outcome of primary or secondary astrophysical processes. Then, we add Dark Matter to the picture, in order to establish which are the informations on its annihilation cross section (or lifetime) that can be inferred by fitting AMS-02 data within a scenario in which Dark Matter and astrophysical sources jointly contribute to the different leptonic observables. In particular, by performing a Markov Chain Monte Carlo sampling of the parameters space of the theory, we attempt at characterizing the significance of a possible Dark Matter contribution to the observed data and we derive robust upper limits on the Dark Matter annihilation/decay rate.

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).

Cosmic-ray antiprotons represent an important channel for dark matter indirect-detection studies. Current measurements of the antiproton flux at the top of the atmosphere and theoretical determinations of the secondary antiproton production in the Galaxy are in good agreement, with no manifest deviation which could point to an exotic contribution in this channel. Therefore, antiprotons can be used as a powerful tool for constraining particle dark matter properties. By using the spectrum of PAMELA data from 50 MV to 180 GV in rigidity, we derive bounds on the dark matter annihilation cross section (or decay rate, for decaying dark matter) for the whole spectrum of dark matter annihilation (decay) channels and under different hypotheses of cosmic-rays transport in the Galaxy and in the heliosphere. For typical models of galactic propagation, the constraints are significantly strong, setting a lower bound on the dark matter mass of a "thermal" relic at about 50-90 GeV for hadronic annihilation channels. These bounds are enhanced to about 150 GeV on the dark matter mass, when large cosmic-rays confinement volumes in the Galaxy are considered, and are reduced to 4-5 GeV for annihilation to light quarks (no bound for heavy-quark production) when the confinement volume is small. Bounds for dark matter lighter than few tens of GeV are due to the low energy part of the PAMELA spectrum, an energy region where solar modulation is relevant: to this aim, we have implemented a detailed solution of the transport equation in the heliosphere, which allowed us not only to extend bounds to light dark matter, but also to determine the uncertainty on the constraints arising from solar modulation modeling. Finally, we estimate the impact of soon-to-come AMS-02 data on the antiproton constraints.

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 search for an excess of antideuterons in the cosmic rays flux has been proposed as a very promising channel for dark matter indirect detection, especially for WIMPs with a low or intermediate mass. With the development of the AMS experiment and the proposal of a future dedicated experiment, i.e. the General Antiparticle Spectrometer (GAPS), there are exciting possibilities for a dark matter detection in the near future. We give an overview on the principal issues related both to the antideuterons production in dark matter annihilation reactions and to their propagation through the interstellar medium and the heliosphere, with a particular focus on the impact of various solar modulation models on the flux at Earth. Lastly, we provide an updated calculation of the reaching capabilities for current and future experiments compatible with the constraints on the dark matter annihilation cross section imposed by the antiproton measurements of PAMELA.