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

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

The AMS-02 experiment measured several secondary-to-primary ratios enabling a detailed study of Galactic cosmic-ray transport. We constrain previously derived benchmark scenarios (based on AMS-02 B/C data only) using other secondary-to-primary ratios, to test the universality of transport and the presence of a low-rigidity diffusion break. We use the 1D thin disc/thick halo propagation model of USINE and a $\chi^2$ minimisation accounting for a covariance matrix of errors (AMS-02 systematics) and nuisance parameters (cross-sections and solar modulation uncertainties). The combined analysis of AMS-02 Li/C, Be/C, and B/C strengthens the case for a diffusion slope of $\delta=0.50\pm 0.03$ with a low-rigidity break or upturn of the diffusion coefficient at GV rigidities. Our simple model can successfully reproduce all considered data (Li/C, Be/C, B/C, N/O, and 3He/4He), although several issues remain: (i) the quantitative agreement depends on the assumptions made on the not well constrained correlation lengths of AMS-02 data systematics; (ii) combined analyses are very sensitive to production cross sections, and we find post-fit values differing by $\sim5-15\%$ from their most likely values (roughly within currently estimated nuclear uncertainties); (iii) two very distinct regions of the parameter space remain viable, either with reacceleration and convection, or with purely diffusive transport. To take full benefit of combined analyses of AMS-02 data, better nuclear data and a better handle on energy correlations in the data systematic are required. AMS-02 data on heavier species are eagerly awaited to further explore cosmic-ray propagation scenarios.

AMS-02 on the International Space Station has been releasing data of unprecedented accuracy. This poses new challenges for their interpretation. We refine the methodology to get a statistically sound determination of the cosmic-ray propagation parameters. We inspect the numerical precision of the model calculation, nuclear cross-section uncertainties, and energy correlations in data systematic errors. We used the 1D diffusion model in USINE. Our $\chi^2$ analysis includes a covariance matrix of errors for AMS-02 systematics and nuisance parameters to account for cross-section uncertainties. Mock data were used to validate some of our choices. We show that any mis-modelling of nuclear cross-section values or the energy correlation length of the covariance matrix of errors biases the analysis. It also makes good models ($\chi^2_{\rm min}/{\rm dof}\approx1$) appear as excluded ($\chi^2_{\rm min}/{\rm dof}\gg1$). We provide a framework to mitigate these effects (AMS-02 data are interpreted in a companion paper). New production cross-section data and the publication by the AMS-02 collaboration of a covariance matrix of errors for each data set would be an important step towards an unbiased view of cosmic-ray propagation in the Galaxy.

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

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

Primordial Black Holes (PBHs) with a mass $M \lesssim {10^{17}}$g are expected to inject sub-GeV electrons and positrons in the Galaxy via Hawking radiation. These cosmic rays are shielded by the solar magnetic field for Earth-bound detectors, but not for Voyager-1, which is now beyond the heliopause. We use its data to constrain the fraction of PBHs to the dark matter in the Galaxy, finding that PBHs with $M<10^{16}$g cannot contribute more than 0.1% (or less for a lognormal mass distribution). Our limits are based on local galactic measurements and are thus complementary to those derived from cosmological observations.

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

We study the implications of dark matter searches, together with collider constraints, on the phenomenological MSSM with neutralino dark matter and focus on the consequences of the related uncertainties in some detail. We consider, inter alia, the latest results from AMS-02, Fermi-LAT and XENON1T. In particular, we examine the impact of the choice of the dark matter halo profile, as well as the propagation model for cosmic rays, for dark matter indirect detection and show that the constraints on the MSSM differ by one to two orders of magnitude depending on the astrophysical hypotheses. On the other hand, our limited knowledge of the local relic density in the vicinity of the Earth and the velocity of Earth in the dark matter halo leads to a factor 3 in the exclusion limits obtained by direct detection experiments. We identified the astrophysical models leading to the most conservative and the most stringent constraints and for each case studied the complementarities with the latest LHC measurements and limits from Higgs, SUSY and monojet searches. We show that combining all data from dark matter searches and colliders, a large fraction of our supersymmetric sample could be probed. Whereas the direct detection constraints are rather robust under the astrophysical assumptions, the uncertainties related to indirect detection can have an important impact on the number of the excluded points.

It is well known that dark matter density measurements, indirect and direct detection experiments, importantly complement the LHC in setting strong constraints on new physics scenarios. Yet, dark matter searches are subject to limitations which need to be considered for realistic analyses. For illustration, we explore the parameter space of the phenomenological MSSM and discuss the interplay of the constraints from dark matter searches and the LHC, and analyse the impact of the astrophysical uncertainties in some detail.

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

Two years ago, the AMS collaboration released the most precise measurement of the cosmic ray positron flux. It confirms that pure secondary predictions fall below the data above 10 GeV, suggesting the presence of a primary component, e.g. annihilations of WIMPs dark matter. Most analyses have focused on the high-energy part of the spectrum, disregarding the GeV energy region where cosmic ray transport is harder to model and solar modulation comes into play. Given the high quality of AMS measurements, it is timely to re-examine the positron anomaly over the entire energy range, taking into account transport processes so far neglected, e.g. convection or diffusive re-acceleration. We devise a new semi-analytical method to take into account transport processes so far neglected, but important below a few GeV. It is based on the pinching of inverse Compton and synchrotron energy losses inside the Galactic disc. It allows to carry out extensive scans over the cosmic ray propagation parameters, which we strongly constrain by requiring that the secondary component does not overshoot the AMS measurements. Only models with large diffusion coefficients survive this test. The positron flux is a powerful and independent probe of cosmic ray propagation, complementary to the boron-to-carbon ratio. We then scan over WIMP mass to fit the annihilation cross section and branching ratios, exploring both direct annihilations into standard model particles or through light mediators. In the former case, the best fit yields a $p$-value of 0.4% for a mass of 264 GeV, a value that does not allow to reproduce the highest energy data points. Worse quality fits are found in the latter case. The interpretation of the positron excess in terms of single DM species annihilations is strongly disfavored. This conclusion is based solely on the positron data, and no other observation needs to be invoked.

Antiprotons are regarded as a powerful probe for Dark Matter (DM) indirect detection and indeed current data from \PAMELA\ have been shown to lead to stringent constraints. However, in order to exploit their constraining/discovery power properly, great attention must be put into effects (linked to their propagation in the Galaxy) which may be perceived as subleading but actually prove to be quite relevant. We revisit the computation of the astrophysical background and of the DM antiproton fluxes fully including the effects of: diffusive reacceleration, energy losses including tertiary component and solar modulation (in a force field approximation). Using the updated proton and helium fluxes just released by the \AMS\ experiment we reevaluate the secondary astrophysical antiproton to proton ratio and its uncertainties, and compare it with the ratio preliminarly reported by \AMS. We find no unambiguous evidence for a significant excess with respect to expectations. Yet, some preference for a flatter energy dependence of the diffusion coefficient (with respect to the {\sc Med} benchmark often used in the literature) starts to emerge. Finally, we provide a first assessment of the room left for exotic components such as Galactic Dark Matter annihilation, deriving new stringent constraints.

The positron fraction in cosmic rays was found to be steadily increasing in function of energy, above $\sim$10 GeV. This behaviour contradicts standard astrophysical mechanisms, in which positrons are secondary particles, produced in the interactions of primary cosmic rays during the propagation in the interstellar medium. The observed anomaly in the positron fraction triggered a lot of excitement, as it could be interpreted as an indirect signature of the presence of dark matter species in the Galaxy. Alternatively, it could be produced by nearby astrophysical sources, such as pulsars. Both hypotheses are probed in this work in light of the latest AMS-02 positron fraction measurements. The transport of primary and secondary positrons in the Galaxy is described using a semi-analytic two-zone model. MicrOMEGAs is used to model the positron flux generated by dark matter species. We provide mass and annihilating cross section that best fit AMS-02 data for each single annihilating channel as well as for combinations of channels. We find that the mass of the favoured dark matter candidates is always larger than 500 GeV. The description of the positron fraction from astrophysical sources is based on the pulsar observations included in the ATNF catalogue. The region of the distance-to-age plane that best fits the positron fraction for a single source is determined and a list of five pulsars from the ATNF catalogue is given. Those results are obtained with the cosmic ray transport parameters that best fit the B/C ratio. Uncertainties in the propagation parameters turn out to be very significant.

Using the updated proton and helium fluxes just released by the AMS-02 experiment we reevaluate the secondary astrophysical antiproton to proton ratio and its uncertainties, and compare it with the ratio preliminarly reported by AMS-02. We find no unambiguous evidence for a significant excess with respect to expectations. Yet, some preference for a flatter energy dependence of the diffusion coefficient starts to emerge. Also, we provide a first assessment of the room left for exotic components such as Galactic Dark Matter annihilation or decay, deriving new stringent constraints.

Antiprotons are regarded as a powerful probe for Dark Matter (DM) indirect detection and indeed current data from PAMELA have been shown to lead to stringent constraints. However, in order to exploit their constraining/discovery power properly and especially in anticipation of the exquisite accuracy of upcoming data from AMS, great attention must be put into effects (linked to their propagation in the Galaxy) which may be perceived as subleasing but actually prove to be quite relevant. We revisit the computation of the astrophysical background and of the DM antiproton fluxes fully including the effects of: diffusive reacceleration, energy losses including tertiary component and solar modulation (in a force field approximation). We show that their inclusion can somewhat modify the current bounds, even at large DM masses, and that a wrong interpretation of the data may arise if they are not taken into account. The numerical results for the astrophysical background are provided in terms of fit functions; the results for Dark Matter are incorporated in the new release of the PPPC4DMID.

The positron fraction in cosmic rays was found to be a steadily increasing in function of energy, above $\sim$ 10 GeV. This behaviour contradicts standard astrophysical mechanisms, in which positrons are secondary particles, produced in the interactions of primary cosmic rays during the propagation in the interstellar medium. The observed anomaly in the positron fraction triggered a lot of excitement, as it could be interpreted as an indirect signature of the presence of dark matter species in the Galaxy. Alternatively, it could be produced by nearby astrophysical sources, such as pulsars. Both hypotheses are probed in this work in light of the latest AMS-02 positron fraction measurements. The transport of the primary and secondary positrons in the Galaxy is described using a semi-analytic two-zone model. MicrOMEGAs is used to model the positron flux generated by dark matter species. The description of the positron fraction from astrophysical sources is based on the pulsar observations included in the ATNF catalogue. We find that the mass of the favoured dark matter candidates is always larger than 500 GeV. The only dark matter species that fulfils the numerous gamma ray and cosmic microwave background bounds is a particle annihilating into four leptons through a light scalar or vector mediator, with a mixture of tau (75%) and electron (25%) channels, and a mass between 0.5 and 1 TeV. The positron anomaly can also be explained by a single astrophysical source and a list of five pulsars from the ATNF catalogue is given. Those results are obtained with the cosmic ray transport parameters that best fit the B/C ratio. Uncertainties in the propagation parameters turn out to be very significant. In the WIMP annihilation cross section to mass plane for instance, they overshadow the error contours derived from the positron data.