High-energy Gamma Rays from the Milky Way: Three-dimensional Spatial Models for the Cosmic-Ray and Radiation Field Densities in the Interstellar Medium

Likely Detection of GeV γ-Ray Emission from Pulsar Wind Nebula G32.64+0.53 with Fermi-LAT

In this study, we report the likely GeV γ-ray emissions originating from the pulsar PSR J1849-0001's pulsar wind nebula (PWN) G32.64+0.53. Our analysis covers approximately 14.7 yr of data from the Fermi Large Area Telescope Pass 8. The position of the source and its spectrum matches those in X-ray and TeV energy bands, so we propose that the GeV γ-ray source is indicative of PWN G32.64+0.53. We interpret the broadband spectral energy distribution (SED) using a time-dependent one-zone model, which assumes that the multiband nonthermal emission of the target source can be generated by synchrotron radiation and inverse Compton scattering (ICS) of the electrons/positrons. Our findings demonstrate that the model substantially elucidates the observed SED. These results lend support to the hypothesis that the γ-ray source originates from the PWN G32.64+0.53 powered by PSR J1849-0001. Furthermore, the γ-rays in TeV bands are likely generated by electrons/positrons within the nebula through ICS.

Weighing the Local Interstellar Medium Using Gamma Rays and Dust

Cold gas forms a significant mass fraction of the Milky Way disk, but is its most uncertain baryonic component. The density and distribution of cold gas is of critical importance for Milky Way dynamics, as well as models of stellar and galactic evolution. Previous studies have used correlations between gas and dust to obtain high-resolution measurements of cold gas, but with large normalization uncertainties. We present a novel approach that uses Fermi-LAT γ-ray data to measure the total gas density, achieving a similar precision as previous works, but with independent systematic uncertainties. Notably, our results have sufficient precision to probe the range of results obtained by current world-leading experiments.

Constraints on Heavy Decaying Dark Matter from 570 Days of LHAASO Observations

The kilometer square array (KM2A) of the large high altitude air shower observatory (LHAASO) aims at surveying the northern γ-ray sky at energies above 10 TeV with unprecedented sensitivity. γ-ray observations have long been one of the most powerful tools for dark matter searches, as, e.g., high-energy γ rays could be produced by the decays of heavy dark matter particles. In this Letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340 days of data from 1/2-KM2A and 230 days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after γ/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest γ-ray constraints on the lifetime of heavy dark matter particles with mass between 10^{5} and 10^{9}  GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

Simulations of cosmic ray propagation

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.

Signatures of Recent Cosmic-Ray Acceleration in the High-latitude Gamma-Ray Sky

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.

On a Possible Origin of the Gamma-ray Excess around the Galactic Center

Recent observations of gamma rays with the Fermi Large Area Telescope (LAT) in the direction of the inner galaxy revealed a mysterious excess of GeV. Its intensity is significantly above predictions of the standard model of cosmic rays (CRs) generation and propagation with a peak in the spectrum around a few GeV. Popular interpretations of this excess are that it is due to either spherically distributed annihilating dark matter (DM) or an abnormal population of millisecond pulsars. We suggest an alternative explanation of the excess through the CR interactions with molecular clouds in the Galactic Center (GC) region. We assumed that the excess could be imitated by the emission of molecular clouds with depleted density of CRs with energies below ∼10 GeV inside. A novelty of our work is in detailed elaboration of the depletion mechanism of CRs with the mentioned energies through the “barrier” near the cloud edge formed by the self-excited MHD turbulence. This depletion of CRs inside the clouds may be a reason for the deficit of gamma rays from the Central Molecular Zone (CMZ) at energies below a few GeV. This in turn changes the ratio between various emission components at those energies and may potentially absorb the GeV excess by a simple renormalization of key components.

Gamma Rays as Probes of Cosmic-Ray Propagation and Interactions in Galaxies

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.

Inference of the Local Interstellar Spectra of Cosmic-Ray Nuclei Z ⩽ 28 with the GalProp-HelMod Framework

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.

Deciphering the Local Interstellar Spectra of Secondary Nuclei with the Galprop/Helmod Framework and a Hint for Primary Lithium in Cosmic Rays

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.

Deciphering Residual Emissions: Time-dependent Models for the Nonthermal Interstellar Radiation from the Milky Way

Cosmic rays (CRs) in the Galaxy are an important dynamical component of the interstellar medium (ISM) that interact with the other major components (interstellar gas and magnetic and radiation fields) to produce broadband interstellar emissions that span the electromagnetic spectrum. The standard modeling of CR propagation and production of the associated emissions is based on a steady-state assumption, where the CR source spatial density is described using a smoothly varying function of position that does not evolve with time. While this is a convenient approximation, reality is otherwise, where primary CRs are produced in and about highly localized regions, e.g., supernova remnants, which have finite lifetimes. In this paper, we use the latest version of the galprop CR propagation code to model time-dependent CR injection and propagation through the ISM from a realistic 3D discretized CR source density distribution, together with full 3D models for the other major ISM components, and make predictions of the associated broadband nonthermal emissions. We compare the predictions for the discretized and equivalent steady-state model, finding that the former predicts novel features in the broadband nonthermal emissions that are absent for the steady-state case. Some of the features predicted by the discretized model may be observable in all-sky observations made by WMAP and Planck, the recently launched eROSITA, the Fermi-LAT, and ground-based observations by HESS, HAWC, and the forthcoming CTA. The nonthermal emissions predicted by the discretized model may also provide explanations of puzzling anomalies in high-energy γ-ray data, such as the Fermi-LAT north/south asymmetry and residuals like the so-called "Fermi bubbles."

Fermi-LAT Observations of γ-Ray Emission toward the Outer Halo of M31

The Andromeda galaxy is the closest spiral galaxy to us and has been the subject of numerous studies. It harbors a massive dark matter halo, which may span up to ~600 kpc across and comprises ~90% of the galaxy's total mass. This halo size translates into a large diameter of 42° on the sky, for an M31-Milky Way (MW) distance of 785 kpc, but its presumably low surface brightness makes it challenging to detect with γ-ray telescopes. Using 7.6 yr of Fermi Large Area Telescope (Fermi-LAT) observations, we make a detailed study of the γ-ray emission between 1-100 GeV toward M31's outer halo, with a total field radius of 60° centered at M31, and perform an in-depth analysis of the systematic uncertainties related to the observations. We use the cosmic-ray propagation code GALPROP to construct specialized interstellar emission models to characterize the foreground γ-ray emission from the MW, including a self-consistent determination of the isotropic component. We find evidence for an extended excess that appears to be distinct from the conventional MW foreground, having a total radial extension upward of ~120-200 kpc from the center of M31. We discuss plausible interpretations of the excess emission, but emphasize that uncertainties in the MW foreground-and in particular, modeling of the H i-related components-have not been fully explored and may impact the results.

Cosmic-Ray Propagation in Light of the Recent Observation of Geminga

The High Altitude Water Cherenkov (HAWC) telescope recently observed extended emission around the Geminga and PSR B0656+14 pulsar wind nebulae (PWNe). These observations have been used to estimate cosmic-ray (CR) diffusion coefficients near the PWNe that appear to be more than two orders of magnitude smaller than the diffusion coefficients typically derived for the interstellar medium from measured abundances of secondary species in CRs. Two-zone diffusion models have been proposed as a solution to this discrepancy, where the slower diffusion zone (SDZ) is confined to a small region around the PWN. Such models are shown to successfully reproduce the HAWC observations of the Geminga PWN while retaining consistency with other CR data. It is found that the size of the SDZ influences the predicted positron flux and the spectral shape of the extended γ-ray emission at lower energies that can be observed with the Fermi Large Area Telescope. If the two observed PWNe are not unique, then it is likely that there are similar pockets of slow diffusion around many CR sources elsewhere in the Milky Way. The consequences of such a picture for Galactic CR propagation is explored.

Cosmic ray e± at high energy

There is a commonly expressed opinion in the literature, that cosmic-ray (CR) e+ come from a primary source, which could be dark matter or pulsars. In these proceedings we review some evidence to the contrary: namely, that e+ come from secondary production due to CR nuclei scattering on interstellar matter. We show that recent measurements of the total e± flux at E ≲ 3 TeV are in good agreement with the predicted flux of secondary e±, that would be obtained if radiative energy losses during CR propagation do not play an important role. If the agreement between data and secondary prediction is not accidental, then the requirement of negligible radiative energy losses implies a very short propagation time for high energy CRs: tesc ≲. 105 yr at rigidities R ≳ 3 TV. Such short propagation history may imply that a recent, near-by source dominates the CRs at these energies. We review independent evidence for a transition in CR propagation, based on the spectral hardening of primary and secondary nuclei around R ~ 100 GV. The transition rigidity of the nuclei matches the rigidity at which the e+ flux saturates its secondary upper bound.

Galactic PeVatrons and helping to find them: Effects of galactic absorption on the observed spectra of very high energy γ-ray sources

Identification of the cosmic-ray (CR) "PeVatrons," which are sources capable of accelerating particles to ~1015 eV energies and higher, may lead to resolving the long-standing question of the origin of the spectral feature in the all-particle CR spectrum known as the "knee." Because CRs with these energies are deflected by interstellar magnetic fields identification of individual sources and determination of their spectral characteristics is more likely via very high energy γ-ray emissions, which provide the necessary directional information. However, pair production on the interstellar radiation field (ISRF) and cosmic microwave background (CMB) leads to steepening of the high energy tails of γ-ray spectra, and should be corrected for to enable true properties of the spectrum at the source to be recovered. Employing recently developed three-dimensional ISRF models this paper quantifies the pair-absorption effect on spectra for sources in the Galactic center (GC) direction at 8.5 and 23.5 kpc distances, with the latter corresponding to the far side of the Galactic stellar disc where it is expected that discrimination of spectral features >10 TeV is possible by the forthcoming Cherenkov Telescope Array (CTA). The estimates made suggest spectral cutoffs could be underestimated by factors of a few in the energy range so far sampled by TeV γ-ray telescopes. As an example to illustrate this, the recent HESS measurements of diffuse γ-ray emissions possibly associated with injection of CRs nearby Sgr A* are ISRF corrected, and estimates of the spectral cutoff are reevaluated. It is found that it could be higher by up to a factor of ~2, indicating that these emissions may be consistent with a CR accelerator with a spectral cutoff of at least 1 PeV at the 95% confidence level.

Deciphering the Local Interstellar Spectra of Primary Cosmic-Ray Species with helmod

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).

The Three-dimensional Spatial Distribution of Interstellar Gas in the Milky Way: Implications for Cosmic Rays and High-energy Gamma-ray Emissions

Direct measurements of cosmic ray (CR) species combined with observations of their associated γ-ray emissions can be used to constrain models of CR propagation, trace the structure of the Galaxy, and search for signatures of new physics. The spatial density distribution of interstellar gas is a vital element for all these studies. So far, models have employed the 2D cylindrically symmetric geometry, but their accuracy is well behind that of the available data. In this paper, 3D spatial density models for neutral and molecular hydrogen are constructed based on empirical model fitting to gas line-survey data. The developed density models incorporate spiral arms and account for the warping of the disk, and the increasing gas scale height with radial distance from the Galactic center. They are employed together with the GALPROP CR propagation code to investigate how the new 3D gas models affect calculations of CR propagation and high-energy γ-ray intensity maps. The calculations reveal non-trivial features that are directly related to the new gas models. The best-fit values for propagation model parameters employing 3D gas models are presented and they differ significantly from those derived with the 2D gas density models that have been widely used. The combination of 3D CR and gas density models provide a more realistic basis for the interpretation of non-thermal emissions from the Galaxy.

HelMod in the Works: From Direct Observations to the Local Interstellar Spectrum of Cosmic-Ray Electrons

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.