High-energy γ-rays of interstellar origin are produced by the interaction of cosmic-ray (CR) particles with the diffuse gas and radiation fields in the Galaxy. The main features of this emission are ...well understood and are reproduced by existing CR propagation models employing 2D galactocentric cylindrically symmetrical geometry. However, the high-quality data from instruments like the Fermi Large Area Telescope reveal significant deviations from the model predictions on few to tens of degrees scales, indicating the need to include the details of the Galactic spiral structure and thus requiring 3D spatial modeling. In this paper, the high-energy interstellar emissions from the Galaxy are calculated using the new release of the GALPROP code employing 3D spatial models for the CR source and interstellar radiation field (ISRF) densities. Three models for the spatial distribution of CR sources are used that are differentiated by their relative proportion of input luminosity attributed to the smooth disk or spiral arms. Two ISRF models are developed based on stellar and dust spatial density distributions taken from the literature that reproduce local near- to far-infrared observations. The interstellar emission models that include arms and bulges for the CR source and ISRF densities provide plausible physical interpretations for features found in the residual maps from high-energy γ-ray data analysis. The 3D models for CR and ISRF densities provide a more realistic basis that can be used for the interpretation of the nonthermal interstellar emissions from the Galaxy.
ABSTRACT Since 2012 August Voyager 1 has been observing the local interstellar energy spectra of Galactic cosmic-ray nuclei down to 3 MeV nuc−1 and electrons down to 2.7 MeV. The H and He spectra ...have the same energy dependence between 3 and 346 MeV nuc−1, with a broad maximum in the 10-50 MeV nuc−1 range and a H/He ratio of 12.2 0.9. The peak H intensity is ∼15 times that observed at 1 AU, and the observed local interstellar gradient of 3-346 MeV H is −0.009 0.055% AU−1, consistent with models having no local interstellar gradient. The energy spectrum of electrons (e− + e+) with 2.7-74 MeV is consistent with E−1.30 0.05 and exceeds the H intensity at energies below ∼50 MeV. Propagation model fits to the observed spectra indicate that the energy density of cosmic-ray nuclei with >3 MeV nuc−1 and electrons with >3 MeV is 0.83-1.02 eV cm−3 and the ionization rate of atomic H is in the range of 1.51-1.64 × 10−17 s−1. This rate is a factor >10 lower than the ionization rate in diffuse interstellar clouds, suggesting significant spatial inhomogeneity in low-energy cosmic rays or the presence of a suprathermal tail on the energy spectrum at much lower energies. The propagation model fits also provide improved estimates of the elemental abundances in the source of Galactic cosmic rays.
ABSTRACT We present the results of the most complete scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the ...MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine-learning package. This is the first study to separate out low-mass isotopes (p, , and He) from the usual light elements (Be, B, C, N, and O). We find that the propagation parameters that best-fit , and He data are significantly different from those that fit light elements, including the B/C and 10Be/9Be secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests that each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best-fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.
Abstract
The past decade has brought impressive advances in the astrophysics of cosmic rays (CRs) and multiwavelength astronomy, thanks to the new instrumentation launched into space and built on the ...ground. Modern technologies employed by those instruments provide measurements with unmatched precision, enabling searches for subtle signatures of dark matter and new physics. Understanding the astrophysical backgrounds to better precision than the observed data is vital in moving to this new territory. A state-of-the-art CR propagation code, called
GalProp,
is designed to address exactly this challenge. Having 25 yr of development behind it, the
GalProp
framework has become a de facto standard in the astrophysics of CRs, diffuse photon emissions (radio to
γ
-rays), and searches for new physics.
GalProp
uses information from astronomy, particle physics, and nuclear physics to predict CRs and their associated emissions self-consistently, providing a unifying modeling framework. The range of its physical validity covers 18 orders of magnitude in energy, from sub-keV to PeV energies for particles and from
μ
eV to PeV energies for photons. The framework and the data sets are public and are extensively used by many experimental collaborations and by thousands of individual researchers worldwide for interpretation of their data and for making predictions. This paper details the latest release of the
GalProp
framework and updated cross sections, further developments of its initially auxiliary data sets for models of the interstellar medium that grew into independent studies of the Galactic structure—distributions of gas, dust, radiation, and magnetic fields—as well as the extension of its modeling capabilities. Example applications included with the distribution illustrating usage of the new features are also described.
Local interstellar spectra (LIS) for protons, helium, and antiprotons are built using the most recent experimental results combined with state-of-the-art models for propagation in the Galaxy and ...heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of cosmic-ray (CR) species at different modulation levels and at both polarities of the solar magnetic field. To do so in a self-consistent way, an iterative procedure was developed, where the GALPROP LIS output is fed into HelMod, providing modulated spectra for specific time periods of selected experiments to compare with the data; the HelMod parameter optimization is performed at this stage and looped back to adjust the LIS using the new GALPROP run. The parameters were tuned with the maximum likelihood procedure using an extensive data set of proton spectra from 1997 to 2015. The proposed LIS accommodate both the low-energy interstellar CR spectra measured by Voyager 1 and the high-energy observations by BESS, Pamela, AMS-01, and AMS-02 made from the balloons and near-Earth payloads; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The found solution is in a good agreement with proton, helium, and antiproton data by AMS-02, BESS, and PAMELA in the whole energy range.
The local interstellar spectrum (LIS) of cosmic-ray (CR) electrons for the energy range 1 MeV to 1 TeV is derived using the most recent experimental results combined with the state-of-the-art models ...for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The optimized HelMod parameters are then used to adjust GALPROP parameters to predict a refined LIS with the procedure repeated subject to a convergence criterion. The parameter optimization uses an extensive data set of proton spectra from 1997 to 2015. The proposed CR electron LIS accommodates both the low-energy interstellar spectra measured by Voyager 1 as well as the high-energy observations by PAMELA and AMS-02 that are made deep in the heliosphere; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The interstellar and heliospheric propagation parameters derived in this study agree well with our earlier results for CR protons, helium nuclei, and anti-protons propagation and LIS obtained in the same framework.
We use the GALPROP code for cosmic-ray (CR) propagation to calculate the broadband luminosity spectrum of the Milky Way related to CR propagation and interactions in the interstellar medium. This ...includes {gamma}-ray emission from the production and subsequent decay of neutral pions ({pi}{sup 0}), bremsstrahlung, and inverse Compton scattering, and synchrotron radiation. The Galaxy is found to be nearly a CR electron calorimeter, but only if {gamma}-ray emitting processes are taken into account. Synchrotron radiation alone accounts for only one-third of the total electron energy losses with {approx}10%-20% of the total synchrotron emission from secondary CR electrons and positrons. The relationship between far-infrared and radio luminosity that we find from our models is consistent with that found for galaxies in general. The results will be useful for understanding the connection between diffuse emissions from radio through {gamma}-rays in 'normal' (non-active galactic nucleus dominated) galaxies as well as for estimating the broadband extragalactic diffuse background from these kinds of galaxies.
Research in many areas of modern physics such as, e.g., indirect searches for dark matter and particle acceleration in supernova remnant shocks rely heavily on studies of cosmic rays (CRs) and ...associated diffuse emissions (radio, microwave, X-rays, Delta *g-rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions. The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, which uses astrophysical information, and nuclear and particle data as inputs to self-consistently predict CRs, Delta *g-rays, synchrotron, and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.
Local interstellar spectra (LIS) of primary cosmic ray (CR) nuclei, such as helium, oxygen, and mostly primary carbon are derived for the rigidity range from 10 MV to ∼200 TV using the most recent ...experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined into a single framework that is used to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. The developed iterative maximum-likelihood method uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The interstellar and heliospheric propagation parameters derived in this study are consistent with our prior analyses using the same methodology for propagation of CR protons, helium, antiprotons, and electrons. The resulting LIS accommodate a variety of measurements made in the local interstellar space (Voyager 1) and deep inside the heliosphere at low (ACE/CRIS, HEAO-3) and high energies (PAMELA, AMS-02).
ABSTRACT It is widely accepted that cosmic rays (CRs) up to at least PeV energies are Galactic in origin. Accelerated particles are injected into the interstellar medium where they propagate to the ...farthest reaches of the Milky Way, including a surrounding halo. The composition of CRs coming to the solar system can be measured directly and has been used to infer the details of CR propagation that are extrapolated to the whole Galaxy. In contrast, indirect methods, such as observations of γ-ray emission from CR interactions with interstellar gas, have been employed to directly probe the CR densities in distant locations throughout the Galactic plane. In this article we use 73 months of data from the Fermi Large Area Telescope in the energy range between 300 MeV and 10 GeV to search for γ-ray emission produced by CR interactions in several high- and intermediate-velocity clouds (IVCs) located at up to ∼7 kpc above the Galactic plane. We achieve the first detection of IVCs in γ rays and set upper limits on the emission from the remaining targets, thereby tracing the distribution of CR nuclei in the halo for the first time. We find that the γ-ray emissivity per H atom decreases with increasing distance from the plane at 97.5% confidence level. This corroborates the notion that CRs at the relevant energies originate in the Galactic disk. The emissivity of the upper intermediate-velocity Arch hints at a 50% decline of CR densities within 2 kpc from the plane. We compare our results to predictions of CR propagation models.