The interactions of cosmic rays with the solar atmosphere produce secondary particle which can reach the Earth. In this work we present a comprehensive calculation of the yields of secondary particles as gamma-rays, electrons, positrons, neutrons and neutrinos performed with the FLUKA code. We also estimate the intensity at the Sun and the fluxes at the Earth of these secondary particles by folding their yields with the intensities of cosmic rays impinging on the solar surface. The results are sensitive on the assumptions on the magnetic field nearby the Sun and to the cosmic-ray transport in the magnetic field in the inner solar system.
The measured fluxes of secondary particles produced by the interactions of Cosmic Rays (CRs) with the astronomical environment play a crucial role in understanding the physics of CR transport. In this work we present a comprehensive calculation of the secondary hadron, lepton, gamma-ray and neutrino yields produced by the inelastic interactions between several species of stable or long-lived cosmic rays projectiles (p, D, T, 3He, 4He, 6Li, 7Li, 9Be, 10Be, 10B, 11B, 12C, 13C, 14C, 14N, 15N, 16O, 17O, 18O, 20Ne, 24Mg and 28Si) and different target gas nuclei (p, 4He, 12C, 14N, 16O, 20Ne, 24Mg, 28Si and 40Ar). The yields are calculated using FLUKA, a simulation package designed to compute the energy distributions of secondary products with large accuracy in a wide energy range. The present results provide, for the first time, a complete and self-consistent set of all the relevant inclusive cross sections regarding the whole spectrum of secondary products in nuclear collisions. We cover, for the projectiles, a kinetic energy range extending from $0.1 GeV/n$ up to $100 TeV/n$ in the lab frame. In order to show the importance of our results for multi-messenger studies about the physics of CR propagation, we evaluate the propagated spectra of Galactic secondary nuclei, leptons, and gamma rays produced by the interactions of CRs with the insterstellar gas, exploiting the numerical codes DRAGON and GammaSky. We show that, adopting our cross section database, we are able to provide a good fit of a complete sample of CR observables, including: leptonic and hadronic spectra measured at Earth, the local interstellar spectra measured by Voyager, and the gamma-ray emissivities from Fermi-LAT collaboration. We also show a set of gamma-ray and neutrino full-sky maps and spectra.
FLUKA is a general purpose Monte Carlo transport and interaction code used for fundamental physics and for a wide range of applications. These include Cosmic Ray Physics (muons, neutrinos, EAS, underground physics), both for basic research and applied studies in space and atmospheric flight dosimetry and radiation damage. A review of the hadronic models available in FLUKA and relevant for the description of cosmic ray air showers is presented in this paper. Recent updates concerning these models are discussed. The FLUKA capabilities in the simulation of the formation and propagation of EM and hadronic showers in the Earths atmosphere are shown.
A self-consistent model of a one-dimensional cosmic-ray (CR) halo around the Galactic disk is formulated with the restriction to a minimum number of free parameters. It is demonstrated that the turbulent cascade of MHD waves does not necessarily play an essential role in the halo formation. Instead, an increase of the Alfven velocity with distance to the disk leads to an efficient generic mechanism of the turbulent redshift, enhancing CR scattering by the self-generated MHD waves. As a result, the calculated size of the CR halo at lower energies is determined by the halo sheath, an energy-dependent region around the disk beyond which the CR escape becomes purely advective. At sufficiently high energies, the halo size is set by the characteristic thickness of the ionized gas distribution. The calculated Galactic spectrum of protons shows a remarkable agreement with observations, reproducing the position of spectral break at ~ 0.6 TeV and the spectral shape up to ~ 10 TeV.
Recent results obtained with gamma-ray satellites have established supernova remnants as accelerators of GeV hadronic cosmic rays. In such processes, CRs accelerated in SNR shocks interact with particles from gas clouds in their surrounding. In particular, the rich medium in which core-collapse SNRs explode provides a large target density to boost hadronic gamma-rays. SNR G39.2-0.3 is one of the brightest SNR in infrared wavelengths, and its broad multiwavelength coverage allows detailed modelling of its radiation from radio to high energies. We reanalyzed the Fermi-LAT data on this region and compare it with new radio observations from the MWISP survey. The modelling of the spectral energy distribution from radio to GeV energies favors a hadronic origin of the gamma-ray emission and constrains the SNR magnetic field to be at least ~100 uG. Despite the large magnetic field, the present acceleration of protons seems to be limited to ~10 GeV, which points to a drastic slow down of the shock velocity due to the dense wall traced by the CO observations, surrounding the remnant. Further investigation of the gamma-ray spectral shape points to a dynamically old remnant subjected to severe escape of CRs and a decrease of acceleration efficiency. The low-energy peak of the gamma-ray spectrum also suggests that that the composition of accelerated particles might be enriched by heavy nuclei which is certainly expected for a core-collapse SNR. Alternatively, the contribution of the compressed pre-existing Galactic cosmic rays is discussed, which is, however, found to not likely be the dominant process for gamma-ray production.