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.

Due to their extreme density and low temperature, neutron stars (NS) are efficient probes to unveil interactions between standard model and dark matter (DM) particles. From elastic scatterings on NS material, DM can get gravitationally trapped by the star. The cooling of DM through further collisions may lead to the formation of a dense core which could collapse into a black hole, thus destroying the whole NS. From the observation of old NS, such a scenario leads to very stringent constraints on the parameter space of asymmetric DM. In this work we reexamine this possibility in detail. This includes: (a) a new detailed determination of the number of DM particles captured, properly taking into account the fact that neutrons form a highly degenerate Fermi material; (b) the determination of the time evolution of the DM density and energy profiles inside the NS, which allows us to understand how, as a function of time, DM thermalizes with NS material; (c) the determination of the corresponding constraints which hold on the DM-neutron cross section, including for the case where a large fraction of DM particles have not thermalized; (d) the first determination of the stringent constraints which also hold in a similar way on the DM-muon cross section, particularly relevant for leptophilic DM models; and (e) the use of realistic NS equations of state in determining these constraints.

The DAMPE experiment recently published its first results on the lepton ($e^+ + e^-$) cosmic-ray (CRs) flux. These results are of importance since they account for the first direct detection of the lepton break around the energy of 1 TeV and confirm the discoveries of ground-based Cherenkov detectors. Meanwhile they reveal a new high-energy feature in the spectrum which triggered a lot of excitement on the theory side, when interpreted as the typical signature of leptophilic dark-matter annihilation. In this proceeding I mainly focus on the theoretical understanding of the lepton break. Then I quickly review the status of the more speculative line-like DAMPE excess, whose astrophysical (pulsar) or exotic (dark matter) explanation is strongly constrained by multi-messenger astronomy.

The accuracy of the current generation of cosmic-ray (CR) experiments, such as AMS-02, PAMELA, CALET, and ISS-CREAM, is now reaching $\sim$1--3\% in a wide range in energy per nucleon from GeV/n to multi-TeV/n. Their correct interpretation could potentially lead to discoveries of new physics and subtle effects that were unthinkable just a decade ago. However, a major obstacle in doing so is the current uncertainty in the isotopic production cross sections that can be as high as 20--50\% or even larger in some cases. While there is a recently reached consensus in the astrophysics community that new measurements of cross sections are desirable, no attempt to evaluate the importance of particular reaction channels and their required accuracy has been made yet. It is, however, clear that it is a huge work that requires an incremental approach. The goal of this study is to provide the ranking of the isotopic cross sections contributing to the production of the most astrophysically important CR Li, Be, B, C, and N species. In this paper, we (i) rank the reaction channels by their importance for a production of a particular isotope, (ii) provide comparisons plots between the models and data used, and (iii) evaluate a generic beam time necessary to reach a 3\% precision in the production cross-sections pertinent to the AMS-02 experiment. This first roadmap may become a starting point in the planning of new measurement campaigns that could be carried out in several nuclear and/or particle physics facilities around the world. A comprehensive evaluation of other isotopes $Z\leq30$ will be a subject of follow-up studies.

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.

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.

In the new precision era for cosmic ray astrophysics, theoretical predictions cannot content themselves with average trends, but need to correctly take into account intrinsic uncertainties. The space-time discreteness of the cosmic ray sources, together with a substantial ignorance of their precise epochs and locations (with the possible exception of the most recent and close ones) play an important role in this sense. We elaborate a statistical theory to deal with this problem, relating the composite probability P({\Psi}) to obtain a flux {\Psi} at the Earth and the single-source probability p({\psi}) to contribute with a flux {\psi}. The main difficulty arises from the fact that p({\psi}) is a fat tail distribution, characterized by power-law or broken power-law behavior up to very large fluxes for which central limit theorem does not hold, and leading to well-known stable laws as opposed Gaussian. We find that relatively simple recipes provide a satisfactory description of the probability P({\Psi}). We also find that a naive Gaussian fit to simulation results would underestimate the probability of very large fluxes, that is, several times above the average, while overestimating the probability of relatively milder excursions. At large energies, large flux fluctuations are prevented by causal considerations, while at low energies a partial knowledge on the recent and nearby population of sources plays an important role. A few proposal have been discussed in the literature to account for spectral breaks recently reported in cosmic ray data in terms of local contributions. We apply our newly developed theory to assess their probabilities, finding that they are relatively small.

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.

PAMELA and, more recently, AMS-02, are ushering us into a new era of greatly reduced statistical uncertainties in experimental measurements of cosmic-ray fluxes. In particular, new determinations of traditional diagnostic tools such as the boron-to-carbon ratio (B/C) are expected to significantly reduce errors on cosmic-ray diffusion parameters, with important implications for astroparticle physics, ranging from inferring primary source spectra to indirect dark matter searches. It is timely to stress, however, that the conclusions obtained crucially depend on the framework in which the data are interpreted as well as from some nuclear input parameters. We aim at assessing the theoretical uncertainties affecting the outcome, with models as simple as possible while still retaining the key dependencies. We compare different semi-analytical, two-zone model descriptions of cosmic-ray transport in the Galaxy. We test for the effect of a primary source contamination in the boron flux by parametrically altering its flux, as well as for nuclear cross section uncertainties. Our study on preliminary results from AMS-02 suggests that, differently for instance from the leptonic case, realistic modelling of the geometry of the Galaxy and of the source distribution are of minor importance to correctly reproduce B/C data at high energies and thus, to a large extent, for the extraction of diffusion parameters. The Ansatz on the lack of primary injection of boron represents the most serious bias, and requires multi-messenger studies to be addressed. If this uncertainty could be lifted, nuclear uncertainties would still represent a serious concern, which degrade the systematic error on the inferred parameters to the 20% level, or three times the estimated experimental sensitivity. In order to reduce this, a new nuclear cross section measurement campaign is probably required.

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.