1
|
Abstract
We review numerical methods for simulations of cosmic ray (CR) propagation on galactic and larger scales. We present the development of algorithms designed for phenomenological and self-consistent models of CR propagation in kinetic description based on numerical solutions of the Fokker–Planck equation. The phenomenological models assume a stationary structure of the galactic interstellar medium and incorporate diffusion of particles in physical and momentum space together with advection, spallation, production of secondaries and various radiation mechanisms. The self-consistent propagation models of CRs include the dynamical coupling of the CR population to the thermal plasma. The CR transport equation is discretized and solved numerically together with the set of MHD equations in various approaches treating the CR population as a separate relativistic fluid within the two-fluid approach or as a spectrally resolved population of particles evolving in physical and momentum space. The relevant processes incorporated in self-consistent models include advection, diffusion and streaming propagation as well as adiabatic compression and several radiative loss mechanisms. We discuss, applications of the numerical models for the interpretation of CR data collected by various instruments. We present example models of astrophysical processes influencing galactic evolution such as galactic winds, the amplification of large-scale magnetic fields and instabilities of the interstellar medium.
Collapse
|
2
|
Jóhannesson G, Porter TA. Signatures of Recent Cosmic-Ray Acceleration in the High-latitude Gamma-Ray Sky. THE ASTROPHYSICAL JOURNAL 2021; 917:30. [PMID: 34776518 PMCID: PMC8587633 DOI: 10.3847/1538-4357/ac01c9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Cosmic-ray (CR) sources temporarily enhance the relativistic particle density in their vicinity over the background distribution accumulated from the Galaxy-wide past injection activity and propagation. If individual sources are close enough to the solar system, their localized enhancements may present as features in the measured spectra of the CRs and in the associated secondary electromagnetic emissions. Large-scale loop-like structures visible in the radio sky are possible signatures of such nearby CR sources. If so, these loops may also have counterparts in the high-latitude γ-ray sky. Using ~10 yr of data from the Fermi Large Area Telescope, applying Bayesian analysis including Gaussian Processes, we search for extended enhanced emission associated with putative nearby CR sources in the energy range from 1 GeV to 1 TeV for the sky region |b| > 30°. We carefully control the systematic uncertainty due to imperfect knowledge of the interstellar gas distribution. Radio Loop IV is identified for the first time as a γ-ray emitter, and we also find significant emission from Loop I. Strong evidence is found for asymmetric features about the Galactic l = 0° meridian that may be associated with parts of the so-called "Fermi Bubbles," and some evidence is also found for γ-ray emission from other radio loops. Implications for the CRs producing the features and possible locations of the sources of the emissions are discussed.
Collapse
Affiliation(s)
- Guđlaugur Jóhannesson
- Science Institute, University of Iceland, IS-107 Reykjavik, Iceland
- Nordita, KTH Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden
| | - Troy A Porter
- W.W. Hansen Experimental Physics Laboratory and Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
3
|
Abstract
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.
Collapse
|
4
|
Boschini MJ, Torre SD, Gervasi M, Grandi D, Jóhannesson G, La Vacca G, Masi N, Moskalenko IV, Pensotti S, Porter TA, Quadrani L, Rancoita PG, Rozza D, Tacconi M. Inference of the Local Interstellar Spectra of Cosmic-Ray Nuclei Z ⩽ 28 with the GalProp-HelMod Framework. THE ASTROPHYSICAL JOURNAL. SUPPLEMENT SERIES 2020; 250:27. [PMID: 34711999 PMCID: PMC8549769 DOI: 10.3847/1538-4365/aba901] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Composition and spectra of Galactic cosmic rays (CRs) are vital for studies of high-energy processes in a variety of environments and on different scales, for interpretation of γ-ray and microwave observations, for disentangling possible signatures of new phenomena, and for understanding of our local Galactic neighborhood. Since its launch, AMS-02 has delivered outstanding-quality measurements of the spectra of p ¯ , e ±, and nuclei: 1H-8O, 10Ne, 12Mg, 14Si. These measurements resulted in a number of breakthroughs; however, spectra of heavier nuclei and especially low-abundance nuclei are not expected until later in the mission. Meanwhile, a comparison of published AMS-02 results with earlier data from HEAO-3-C2 indicates that HEAO-3-C2 data may be affected by undocumented systematic errors. Utilizing such data to compensate for the lack of AMS-02 measurements could result in significant errors. In this paper we show that a fraction of HEAO-3-C2 data match available AMS-02 measurements quite well and can be used together with Voyager 1 and ACE-CRIS data to make predictions for the local interstellar spectra (LIS) of nuclei that are not yet released by AMS-02. We are also updating our already-published LIS to provide a complete set from 1H-28Ni in the energy range from 1 MeV nucleon-1 to ~100-500 TeV nucleon-1, thus covering 8-9 orders of magnitude in energy. Our calculations employ the GalProp-HelMod framework, which has proved to be a reliable tool in deriving the LIS of CR p ¯ , e -, and nuclei 1H-8O.
Collapse
Affiliation(s)
- M J Boschini
- INFN, Milano-Bicocca, Milano, Italy
- CINECA, Segrate, Milano, Italy
| | | | - M Gervasi
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - D Grandi
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - G Jóhannesson
- Science Institute, University of Iceland, Dunhaga 3, IS-107 Reykjavik, Iceland
- NORDITA, Roslagstullsbacken 23, 106 91 Stockholm, Sweden
| | - G La Vacca
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - N Masi
- INFN, Bologna, Italy
- Physics Department, University of Bologna, Bologna, Italy
| | - I V Moskalenko
- Hansen Experimental Physics Laboratory, Stanford University, Stanford, CA 94305, USA
- Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA
| | - S Pensotti
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - T A Porter
- Hansen Experimental Physics Laboratory, Stanford University, Stanford, CA 94305, USA
- Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA
| | - L Quadrani
- INFN, Bologna, Italy
- Physics Department, University of Bologna, Bologna, Italy
| | | | - D Rozza
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - M Tacconi
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| |
Collapse
|
5
|
Boschini MJ, Della Torre S, Gervasi M, Grandi D, Jøhannesson G, La Vacca G, Masi N, Moskalenko IV, Pensotti S, Porter TA, Quadrani L, Rancoita PG, Rozza D, Tacconi M. Deciphering the Local Interstellar Spectra of Secondary Nuclei with the Galprop/Helmod Framework and a Hint for Primary Lithium in Cosmic Rays. THE ASTROPHYSICAL JOURNAL 2020; 889:167. [PMID: 34646048 PMCID: PMC8506946 DOI: 10.3847/1538-4357/ab64f1] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Local interstellar spectra (LIS) of secondary cosmic-ray (CR) nuclei, lithium, beryllium, boron, and partially secondary nitrogen, are derived in the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. The lithium spectrum appears somewhat flatter at high energies compared to other secondary species, which may imply a primary lithium component. 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 the model-data comparison. The proposed LIS accommodates the low-energy interstellar spectra measured by Voyager 1, the High Energy Astrophysics Observatory-3 (HEAO-3), and the Cosmic Ray Isotope Spectrometer on board of the Advanced Composition Explorer (ACE/CRIS), as well as the high-energy observations by the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA), Alpha Magnetic Spectrometer-02 (AMS-02), and earlier experiments that are made deep in the heliosphere. The interstellar and heliospheric propagation parameters derived in this study are consistent with our earlier results for propagation of CR protons, helium, carbon, oxygen, antiprotons, and electrons.
Collapse
Affiliation(s)
- M J Boschini
- INFN, Milano-Bicocca, Milano, Italy
- CINECA, Segrate, Milano, Italy
| | | | - M Gervasi
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - D Grandi
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - G Jøhannesson
- Science Institute, University of Iceland, Dunhaga 3, IS-107, Reykjavik, Iceland
- NORDITA, Roslagstullsbacken 23, SE-106 91, Stockholm, Sweden
| | - G La Vacca
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | | | - I V Moskalenko
- Hansen Experimental Physics Laboratory, Stanford University, Stanford, CA 94305, USA
- Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA
| | - S Pensotti
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - T A Porter
- Hansen Experimental Physics Laboratory, Stanford University, Stanford, CA 94305, USA
- Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA
| | - L Quadrani
- INFN, Bologna, Italy
- Physics Department, University of Bologna, Bologna, Italy
| | | | - D Rozza
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| | - M Tacconi
- INFN, Milano-Bicocca, Milano, Italy
- Physics Department, University of Milano-Bicocca, Milano, Italy
| |
Collapse
|