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 boron to carbon ratio (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 nonrelativistic effects. All models lead to robust predictions in the high-energy regime (≳10 GV), i.e., independent of the propagation scenario: at 1σ, the diffusion slope δ is 0.43–0.53, whereas K10, the diffusion coefficient at 10 GV, is 0.26–0.36 kpc2 Myr−1; we confirm the robustness of the high-energy break, with a typical value Δh∼0.2. We also find a hint for a similar (reversed) feature at low rigidity around the B/C peak (∼4 GV) which might be related to some effective damping scale in the magnetic turbulence.
Context. Top-of-atmosphere (TOA) cosmic-ray (CR) fluxes from satellites and balloon-borne experiments are snapshots of the solar activity imprinted on the interstellar (IS) fluxes. Given a series of ...snapshots, the unknown IS flux shape and the level of modulation (for each snapshot) can be recovered. Aims. We wish (i) to provide the most accurate determination of the IS H and He fluxes from TOA data alone; (ii) to obtain the associated modulation levels (and uncertainties) while fully accounting for the correlations with the IS flux uncertainties; and (iii) to inspect whether the minimal force-field approximation is sufficient to explain all the data at hand. Methods. Using H and He TOA measurements, including the recent high-precision AMS, BESS-Polar, and PAMELA data, we performed a non-parametric fit of the IS fluxes JISH,~He and modulation level φi for each data-taking period. We relied on a Markov chain Monte Carlo (MCMC) engine to extract the probability density function and correlations (hence the credible intervals) of the sought parameters. Results. Although H and He are the most abundant and best measured CR species, several datasets had to be excluded from the analysis because of inconsistencies with other measurements. From the subset of data passing our consistency cut, we provide ready-to-use best-fit and credible intervals for the H and He IS fluxes from MeV/n to PeV/n energy (with a relative precision in the range 2−10% at 1σ). Given the strong correlation between JIS and φi parameters, the uncertainties on JIS translate into Δφ ≈ ± 30 MV (at 1σ) for all experiments. We also find that the presence of 3He in He data biases φ towards higher φ values by ~30 MV. The force-field approximation, despite its limitation, gives an excellent (χ2/d.o.f. = 1.02) description of the recent high-precision TOA H and He fluxes. Conclusions. The analysis must be extended to different charge species and more realistic modulation models. It would benefit from the AMS-02 unique capability of providing frequent high-precision snapshots of the TOA fluxes over a full solar cycle.
Primary cosmic-ray elemental spectra have been measured with the balloon-borne Cosmic Ray Energetics And Mass (CREAM) experiment since 2004. The third CREAM payload (CREAM-III) flew for 29 days ...during the 2007-2008 Antarctic season. Energies of incident particles above 1 TeV are measured with a calorimeter. Individual elements are clearly separated with a charge resolution of ∼0.12 e (in charge units) and ∼0.14 e for protons and helium nuclei, respectively, using two layers of silicon charge detectors. The measured proton and helium energy spectra at the top of the atmosphere are harder than other existing measurements at a few tens of GeV. The relative abundance of protons to helium nuclei is 9.53 0.03 for the range of 1 TeV/n to 63 TeV/n. This ratio is considerably smaller than other measurements at a few tens of GeV/n. The spectra become softer above ∼20 TeV. However, our statistical uncertainties are large at these energies and more data are needed.
On-going measurements of the cosmic radiation (nuclear, electronic, and gamma-ray) are shedding new light on cosmic-ray physics. A comprehensive picture of these data relies on an accurate ...determination of the transport and source parameters of propagation models. A Markov Chain Monte Carlo is used to obtain these parameters in a diffusion model. From the measurement of the B/C ratio and radioactive cosmic-ray clocks, we calculate their probability density functions, with a special emphasis on the halo size L of the Galaxy and the local underdense bubble of size r_h. The analysis relies on the USINE code for propagation and on a Markov Chain Monte Carlo technique (Putze et al. 2009, paper I of this series) for the parameter determination. As found in previous studies, the B/C best-fit model favours diffusion/convection/reacceleration (Model III) over diffusion/reacceleration (Model II). A combined fit on B/C and the isotopic ratios (10Be/9Be, 26Al/27Al, 36Cl/Cl) leads to L ~ 8 kpc and r_h ~ 120 pc for the best-fit Model III. This value for r_h is consistent with direct measurements of the local interstallar medium. For Model II, L ~ 4 kpc and r_h is consistent with zero. We showed the potential and usefulness of the Markov Chain Monte Carlo technique in the analysis of cosmic-ray measurements in diffusion models. The size of the diffusive halo depends crucially on the value of the diffusion slope delta, and also on the presence/absence of the local underdensity damping effect on radioactive nuclei. More precise data from on-going experiments are expected to clarify this issue.
Context.
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.
Aims.
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.
Methods.
We used
USINE V3.5
for the propagation of nuclei, and e
+
were calculated with the pinching method.
Results.
The current AMS-02 Be/B (∼3% precision) and ACE-CRIS
10
Be/Be (∼10% precision) data bring similar and consistent constraints on
L
. The AMS-02 Be/B data alone constrain
L
= 5
−2
+3
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
.
Conclusions.
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.
Context. The secondary-to-primary B/C ratio is widely used for studying Galactic cosmic-ray propagation processes. The 2H/4He and 3He/4He ratios probe a different Z/A regime, which provides a test ...for the “universality” of propagation. Aims. We revisit the constraints on diffusion-model parameters set by the quartet (1H, 2H, 3He, 4He), using the most recent data as well as updated formulae for the inelastic and production cross-sections. Methods. Our analysis relies on the USINE propagation package and a Markov Chain Monte Carlo technique to estimate the probability density functions of the parameters. Simulated data were also used to validate analysis strategies. Results. The fragmentation of CNO cosmic rays (resp. NeMgSiFe) on the interstellar medium during their propagation contributes to 20% (resp. 20%) of the 2H and 15% (resp. 10%) of the 3He flux at high energy. The C to Fe elements are also responsible for up to 10% of the 4He flux measured at 1 GeV/n. The analysis of 3He/4He (and to a lesser extent 2H/4He) data shows that the transport parameters are consistent with those from the B/C analysis: the diffusion model with δ ~ 0.7 (diffusion slope), Vc ~ 20 km s-1 (galactic wind), Va ~ 40 km s-1 (reacceleration) is favoured, but the combination δ ~ 0.2, Vc ~ 0, and Va ~ 80 km s-1 is a close second. The confidence intervals on the parameters show that the constraints set by the quartet data can compete with those derived from the B/C data. These constraints are tighter when adding the 3He (or 2H) flux measurements, and the tightest when the He flux is added as well. For the latter, the analysis of simulated and real data shows an increased sensitivity to biases. Using the secondary-to-primary ratio along with a loose prior on the source parameters is recommended to obtain the most robust constraints on the transport parameters. Conclusions. Light nuclei should be systematically considered in the analysis of transport parameters. They provide independent constraints that can compete with those obtained from the B/C analysis.
Particles count rates at given Earth location and altitude result from the convolution of (i) the interstellar (IS) cosmic-ray fluxes outside the solar cavity, (ii) the time-dependent modulation of ...IS into Top-of-Atmosphere (TOA) fluxes, (iii) the rigidity cut-off (or geomagnetic transmission function) and grammage at the counter location, (iv) the atmosphere response to incoming TOA cosmic rays (shower development), and (v) the counter response to the various particles/energies in the shower. Count rates from neutron monitors or muon counters are therefore a proxy to solar activity. In this paper, we review all ingredients, discuss how their uncertainties impact count rate calculations, and how they translate into variation/uncertainties on the level of solar modulation ϕ (in the simple Force-Field approximation). The main uncertainty for neutron monitors is related to the yield function. However, many other effects have a significant impact, at the 5–10% level on ϕ values. We find no clear ranking of the dominant effects, as some depend on the station position and/or the weather and/or the season. An abacus to translate any variation of count rates (for neutron and μ detectors) to a variation of the solar modulation ϕ is provided.
Abstract
The Cosmic Ray Energetics And Mass for the International Space Station (ISS-CREAM) experiment successfully recorded data for 539 days from 2017 August to 2019 February. We report the energy ...spectrum of cosmic-ray protons from the ISS-CREAM experiment at energies from 1.60 × 10
3
to 6.55 × 10
5
GeV. The measured spectrum deviates from a single power law. A smoothly broken power-law fit to the data, including statistical and systematic uncertainties, shows the spectral index change at 9.0 × 10
3
GeV from 2.57 ± 0.03 to 2.82 ± 0.02 with a significance of greater than 3
σ
. This bump-like structure is consistent with a spectral softening recently reported by the balloon-borne CREAM, DAMPE, and NUCLEON, but ISS-CREAM extends measurements to higher energies.
A Bonner sphere spectrometer extended to high energies (HERMEIS) was set up at the summit of the Pic du Midi de Bigorre in the French Pyrenees (altitude: +2,885 m; geomagnetic cutoff: 5.6 GV) in May ...2011. The spectral fluence rate distribution of the cosmic ray induced neutrons was continuously measured over a broad energy range from meV up to several GeV and with a 1 h time resolution. While the Sun's activity was increasing and reaching its 24th maximum in the 11 year solar cycle, some Forbush decreases were observed in the atmospheric secondary radiation at mountain altitude. We investigated the evolution of the cascade fluence rate (i.e., neutrons with energy greater than 20 MeV) during the March 2012 events with a series of strong coronal mass ejections hitting the Earth's magnetosphere. The amplitude of the greatest Forbush decrease peaked at 10%. Then, a simulation work based on the GEANT4 toolkit was carried out to quantify the solar modulation induced on the galactic cosmic ray transportation during these events. We performed calculations of extensive air showers generated by monoenergetic primaries (Hydrogen and Helium nuclei) for several zenith incidences. Hence, a complete database was built and validated. We derived an analytical model to estimate the atmospheric neutron spectrum at the Pic du Midi according to primary spectra which only depend on the solar modulation potential (force field approximation). We compared the solar modulation potentials obtained in March 2012 with the ones derived by the neutron monitor yield method. Finally, a satisfying agreement was found.
Key Points
A neutron spectrometer is operating at the Pic du Midi de Bigorre
Extensive air shower simulations were performed with GEANT4 toolkit
Solar modulation potential is derived with two methods