A Markov Chain Monte Carlo technique to sample transport and source parameters of Galactic cosmic rays: II. Results for the diffusion model combining B/C and radioactive nuclei

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.Comment: 20 pages, 14 figures, minor language corrections to match the A&A accepted versio

Nuclear Cosmic Rays propagation in the Atmosphere

The transport of the nuclear cosmic ray flux in the atmosphere is studied and the atmospheric corrections to be applied to the measurements are calculated. The contribution of the calculated corrections to the accuracy of the experimental results are discussed and evaluated over the kinetic energy range 10-10$^{3}$ GeV/n. The Boron (B) and Carbon (C) elements system is used as a test case. It is shown that the required corrections become largely dominant at the highest energies investigated. The results are discussed.Comment: Proc. of 30th International Cosmic Ray Conference, Merida, Mexico; 4 page

Simulation of neutrino and charged particle production and propagation in the atmosphere

A precise evaluation of the secondary particle production and propagation in the atmosphere is very important for the atmospheric neutrino oscillation studies. The issue is addressed with the extension of a previously developed full 3-Dimensional Monte-Carlo simulation of particle generation and transport in the atmosphere, to compute the flux of secondary protons, muons and neutrinos. Recent balloon borne experiments have performed a set of accurate flux measurements for different particle species at different altitudes in the atmosphere, which can be used to test the calculations for the atmospheric neutrino production, and constrain the underlying hadronic models. The simulation results are reported and compared with the latest flux measurements. It is shown that the level of precision reached by these experiments could be used to constrain the nuclear models used in the simulation. The implication of these results for the atmospheric neutrino flux calculation are discussed.Comment: 11 pages, 9 figure

Neutron monitors and muon detectors for solar modulation studies: Interstellar flux, yield function, and assessment of critical parameters in count rate calculations

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 $\varphi$ (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 $\varphi$ 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 $\mu$ detectors) to a variation of the solar modulation $\varphi$ is provided.Comment: 28 pages, 16 figures, 9 tables, match accepted version in AdSR (minor corrections, Dorman (1974,2004,2009) reference textbooks added

Neutron monitors and muon detectors for solar modulation studies: 2. ϕ\phi time series

The level of solar modulation at different times (related to the solar activity) is a central question of solar and galactic cosmic-ray physics. In the first paper of this series, we have established a correspondence between the uncertainties on ground-based detectors count rates and the parameter $\phi$ (modulation level in the force-field approximation) reconstructed from these count rates. In this second paper, we detail a procedure to obtain a reference $\phi$ time series from neutron monitor data. We show that we can have an unbiased and accurate $\phi$ reconstruction ($\Delta\phi/\phi\simeq 10\%$). We also discuss the potential of Bonner spheres spectrometers and muon detectors to provide $\phi$ time series. Two by-products of this calculation are updated $\phi$ values for the cosmic-ray database and a web interface to retrieve and plot $\phi$ from the 50's to today (\url{http://lpsc.in2p3.fr/crdb}).Comment: 15 pages, 5 figures, 2 tables. AdSR, in press. Web interface to get modulation parameter phi(t): new tab in http://lpsc.in2p3.fr/crd

An improved gamma-ray limit on the density of PBHs

Gamma-rays are, with antiprotons, a very efficient way to derive upper limits on the density of evaporating black holes. They have been successfully used in the last decades to severely constrain the amount of Primordial Black Holes (PBHs) in our Universe. This article suggests a little refinement, based on the expected background, to improve this limit by a factor of three. The resulting value is : Omega_PBH < 3.3E-9.Comment: To appear in the 28th ICRC conference proceeding

Origin of the high energy proton component below the geomagnetic cutoff in near earth orbit

The high flux proton component observed by AMS below the geomagnetic cutoff can be well accounted for by assuming these particles to be secondaries originating from the interaction of Cosmic Ray protons with the atmosphere. Simulation results are reporte