No Arabic abstract
The Picard code for the numerical solution of the Galactic cosmic ray propagation problem allows for high-resolution models that acknowledge the 3D structure of our Galaxy. Picard was used to determine diffuse gamma-ray emission of the Galaxy over the energy range from 100 MeV to 100 TeV. We discuss the impact of a cosmic-ray source distribution aligned with the Galactic spiral arms for a range of such spiral-arm models. As expected, the impact on the gamma-ray emission is most distinct in the inverse-Compton channel, where imprints of the spiral arms are visible and yield predictions that are no longer symmetric to the rotational axis of the Milkyway. We will illustrate these differences by a direct comparison to results from previous axially symmetric Galactic propagation models: we find differences in the gamma-ray flux both on global scales and on local scales related to the spiral arm tangents. We compare gamma-ray flux and spectra at on-arm vs. off-arm projections and characterize the differences to axially symmetric models.
The Tibet ASgamma experiment just reported their measurement of sub-PeV diffuse gamma ray emission from the Galactic disk, with the highest energy up to 957 TeV. These gamma-rays are most likely the hadronic origin by cosmic ray interaction with interstellar gas in the Galaxy. This measurement provides direct evidence to the hypothesis that the Galactic cosmic rays can be accelerated beyond PeV energies. In this work, we try to explain the sub-PeV diffuse gamma-ray spectrum within cosmic rays diffusive propagation model. We find there is a tension between the sub-PeV diffuse gamma rays and the local cosmic ray spectrum. To describe the sub-PeV diffuse gamma-ray flux, it generally requires larger local cosmic-ray flux than measurement in the knee region. We further calculate the PeV neutrino flux from the cosmic ray propagation model. Even all of these sub-PeV diffuse gamma rays originate from the propagation, the Galactic neutrinos only account for less than ~15% of observed flux, most of which are still from extragalactic sources.
The propagation of particles accelerated at supernova remnant shocks and escaping the parent remnants is likely to proceed in a strongly non-linear regime, due to the efficient self-generation of Alfven waves excited through streaming instability near the sources. Depending on the amount of neutral hydrogen present in the regions around the sites of supernova explosions, cosmic rays may accumulate an appreciable grammage in the same regions and get self-confined for non-negligible times, which in turn results in an enhanced rate of production of secondaries. Here we calculate the contribution to the diffuse gamma-ray background due to the overlap along lines of sight of several of these extended halos as due to pion production induced by self-confined cosmic rays. We find that if the density of neutrals is low, the halos can account for a substantial fraction of the diffuse emission observed by Fermi-LAT, depending on the orientation of the line of sight with respect to the direction of the Galactic centre.
More than 90% of the Galactic gas-related gamma-ray emissivity above 1 GeV is attributed to the decay of neutral pions formed in collisions between cosmic rays and interstellar matter, with lepton-induced processes becoming increasingly important below 1 GeV. Given the high-quality measurements of the gamma-ray emissivity of local interstellar gas between ~50 MeV and ~4 GeV obtained with the Large Area Telescope on board the Fermi space observatory, it is timely to re-investigate this topic in detail, including the hadronic production mechanisms. The emissivity spectrum will allow the interstellar cosmic-ray spectrum to be determined reliably, providing a reference for origin and propagation studies as well as input to solar modulation models. A method for such an analysis and illustrative results are presented.
Most of the diffuse Galactic GeV gamma-ray emission is produced via collisions of cosmic ray (CR) protons with ISM protons. As such the observed spectra of the gamma-rays and the CRs should be strongly linked. Recent observations of Fermi-LAT exhibit a hardening of the gamma-ray spectrum at around a hundred GeV, between the Sagittarius and Carina tangents, and a further hardening at a few degrees above and below the Galactic plane. However, standard CR propagation models that assume a time independent source distribution and a location independent diffusion cannot give rise to a spatially dependent CR (and hence gamma-ray) spectral slopes. Here we consider a dynamic spiral arm model in which the distribution of CR sources is concentrated in the (dynamic) spiral arms, and we study the effects of this model on the $pi^0$-decay produced gamma-ray spectra. Within this model, near the Galactic arms the observed gamma-ray spectral slope is not trivially related to the CR injection spectrum and energy dependence of the diffusion coefficient. We find unique signatures that agree with the Fermi-LAT observations. This model also provides a physical explanation for the difference between the local CR spectral slope and the CR slope inferred from the average gamma-ray spectrum.
Continuum gamma-ray emission produced by interactions of cosmic rays with interstellar matter and radiation fields is a probe of non-thermal particle populations in galaxies. After decades of continuous improvements in experimental techniques and an ever-increasing sky and energy coverage, gamma-ray observations reveal in unprecedented detail the properties of galactic cosmic rays. A variety of scales and environments are now accessible to us, from the local interstellar medium near the Sun and the vicinity of cosmic-ray accelerators, out to the Milky Way at large and beyond, with a growing number of gamma-ray emitting star-forming galaxies. Gamma-ray observations have been pushing forward our understanding of the life cycle of cosmic rays in galaxies and, combined with advances in related domains, they have been challenging standard assumptions in the field and have spurred new developments in modelling approaches and data analysis methods. We provide a review of the status of the subject and discuss perspectives on future progress.