No Arabic abstract
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 $1sigma$, 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_hsim 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.
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.
We study the propagation and injection models of cosmic rays using the latest measurements of the Boron-to-Carbon ratio and fluxes of protons, Helium, Carbon, and Oxygen nuclei by the Alpha Magnetic Spectrometer and the Advanced Composition Explorer at top of the Earth, and the Voyager spacecraft outside the heliosphere. The ACE data during the same time interval of the AMS-02 data are extracted to minimize the complexity of the solar modulation effect. We find that the cosmic ray nucleus data favor a modified version of the diffusion-reacceleration scenario of the propagation. The diffusion coefficient is, however, required to increase moderately with decreasing rigidity at low energies, which has interesting implications on the particle and plasma interaction in the Milky Way. We further find that the low rigidity ($<$ a few GV) injection spectra are different for different compositions. The injection spectra are softer for lighter nuclei. These results are expected to be helpful in understanding the acceleration process of cosmic rays.
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.50pm 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.
In this work, we considered 2 schemes (a high-rigidity break in primary source injections and a high-rigidity break in diffusion coefficient) to reproduce the newly released AMS-02 nuclei spectra (He, C, N, O, Li, Be, and B) when the rigidity larger than 50 GV. The fitting results show that current data set favors a high-rigidity break at $sim 325 mathrm{GV}$ in diffusion coefficient rather than a break at $sim 365 mathrm{GV}$ in primary source injections. Meanwhile, the fitted values of the factors to rescale the cosmic-ray (CR) flux of secondary species/components after propagation show us that the secondary flux are underestimated in current propagation model. It implies that we might locate in a slow diffusion zone, in which the CRs propagate with a small value of diffusion coefficient compared with the averaged value in the galaxy. Another hint from the fitting results show that extra secondary CR nuclei injection may be needed in current data set. All these new hints should be paid more attention in future research.
The latest AMS-02 data on cosmic ray electrons show a break in the energy spectrum around 40 GeV, with a change in the slope of about 0.1. We perform a combined fit to the newest AMS-02 positron and electron flux data using a model which includes production of pairs from pulsar wind nebulae (PWNe), electrons from supernova remnants (SNRs) and both species from spallation of hadronic cosmic rays with interstellar medium atoms. We demonstrate that the change of slope in the AMS-02 electron data is well explained by the interplay between the flux contributions from SNRs and from PWNe. In fact, the relative contribution to the data of these two populations changes by a factor of about 13 from 10 to 1000 GeV. The effect of the energy losses alone, when the inverse Compton scattering is properly computed within a fully numerical treatment of the Klein-Nishina cross section, cannot explain the break in the $e^-$ flux data, as recently proposed in the literature.