Sources Of Cosmic Rays Research Articles

Cosmic-ray transport and gamma-ray emission in M31

We study the possibility that an extended cosmic-ray leptonic and/or hadronic halo is at the origin of the large-scale gamma-ray emission detected from the Andromeda Galaxy (M31). We consider a broad ensemble of non-homogeneous diffusion scenarios and of cosmic-ray injection sources. We find that cosmic-ray electrons and protons could be, and very likely are, responsible for part, or all, of the gamma-ray emission from M31, including out to more than 100 kpc from the center of the galaxy. We also simulate possible emission from pulsars in M31, and consider the effect of regions of highly inefficient diffusion around cosmic-ray acceleration sites, as suggested by recent TeV halo observations with Cherenkov telescopes.

Mirror Diffusion of Cosmic Rays in Highly Compressible Turbulence Near Supernova Remnants

Abstract Recent gamma-ray observations have revealed inhomogeneous diffusion of cosmic rays (CRs) in the interstellar medium (ISM). This is expected, as the diffusion of CRs depends on the properties of turbulence, which can vary widely in the multiphase ISM. We focus on the mirror diffusion arising in highly compressible turbulence in molecular clouds (MCs) around supernova remnants (SNRs), where the magnetic mirroring effect results in significant suppression of diffusion of CRs near CR sources. Significant energy loss via proton–proton interactions due to slow diffusion flattens the low-energy CR spectrum, while the high-energy CR spectrum is steepened due to the strong dependence of mirror diffusion on CR energy. The resulting broken power-law spectrum of CRs matches well the gamma-ray spectrum observed from SNR/MC systems, e.g., IC443 and W44.

Multimessenger Astronomy of Transient Point Sources at the Pierre Auger Observatory

One of the key challenges in astroparticle physics is the identification of the sources of cosmic rays at the highest energies (above 1018 eV). In this context, the search for neutral messenger particles in the ultra-high energy (UHE) regime is of high interest. The sources of the gravitational waves (GWs) that can be observed with the current generation of GW detectors provide extreme astrophysical environments that are most likely to be unique in the universe. Another extraordinary source candidate is the anomalous blazar TXS 0506+056 which has been found to be coincident with two periods of enhanced high-energy neutrino flux in 2014/15 and 2017 as reported by the IceCube Collaboration. Due to their distance and transient nature, the capabilities of these sources to produce UHE radiation can only be studied through neutral messengers like photons and neutrinos. The Pierre Auger Observatory near Malargue, Argentina, is the largest air shower experiment for the detection of ultra-high energy cosmic rays. With its surface detector, consisting of a grid of 1660 water-Cherenkov detectors covering an area of 3000 km2, it has a unique exposure to UHE photons and neutrinos and has published first constraints on these particles from GW sources and TXS 0506+056.

Further evidence for an increase of the mean mass of the highest-energy cosmic-rays with energy

Sokolsky and D'Avignon have recently reported an examination of a variety of measurements that relate to the question as to how the mass composition of the highest-energy cosmic-rays evolves with energy. They assert that cosmic rays arriving from the Northern Hemisphere have a different mass composition from those arriving from the Southern Hemisphere, implying a diversity of sources of high-energy cosmic-rays in the two hemispheres. Were this conclusion to be correct, it would have profound implications for theories of cosmic-ray origin and would influence planning of future projects. Their claim thus merits careful scrutiny. In this paper their analysis is examined in detail with the verdict being that evidence for a North/South difference is not proven, a conclusion supported by other data from the Northern and Southern Hemispheres. However, what is of major importance is that the study of Sokolsky and D'Avignon provides long-awaited confirmation of the claim that the mean mass of cosmic rays increases with energy above ∼3 EeV made by the Pierre Auger Collaboration in 2014.

A GeV-TeV particle component and the barrier of cosmic-ray sea in the Central Molecular Zone

Cosmic rays are important probe of a number of fundamental physical problems such as the acceleration of high and very high energy particles in extreme astrophysical environments. The Galactic center is widely anticipated to be an important cosmic-ray source and the observations of some Imaging Atmospheric Cherenkov Telescopes did successfully reveal a component of TeV-PeV cosmic rays in the vicinity of the Galactic center. Here we report the identification of GeV-TeV cosmic rays in the central molecular zone with the γ-ray observations of the Fermi Large Area Telescope, whose spectrum and spatial gradient are consistent with that measured by the Imaging Atmospheric Cherenkov Telescopes but the corresponding cosmic-ray energy density is substantially lower than the so-called cosmic-ray sea component, suggesting the presence of a high energy particle accelerator at the Galactic center and the existence of a barrier that can effectively suppress the penetration of the particles from the cosmic-ray sea to the central molecular zone.

LMC N132D: A mature supernova remnant with a power-law gamma-ray spectrum extending beyond 8 TeV

Context. Supernova remnants (SNRs) are commonly thought to be the dominant sources of Galactic cosmic rays up to the knee of the cosmic-ray spectrum at a few PeV. Imaging Atmospheric Cherenkov Telescopes have revealed young SNRs as very-high-energy (VHE, >100 GeV) gamma-ray sources, but for only a few SNRs the hadronic cosmic-ray origin of their gamma-ray emission is indisputably established. In all these cases, the gamma-ray spectra exhibit a spectral cutoff at energies much below 100 TeV and thus do not reach the PeVatron regime. Aims. The aim of this work was to achieve a firm detection for the oxygen-rich SNR LMC N132D in the VHE gamma-ray domain with an extended set of data, and to clarify the spectral characteristics and the localization of the gamma-ray emission from this exceptionally powerful gamma-ray-emitting SNR. Methods. We analyzed 252 h of High Energy Stereoscopic System (H.E.S.S.) observations towards SNR N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 h of observations to the previously published data set to ensure a > 5σ detection. To broaden the gamma-ray spectral coverage required for modeling the spectral energy distribution, an analysis of Fermi-LAT Pass 8 data was also included. Results. We unambiguously detect N132D at VHE with a significance of 5.7σ. We report the results of a detailed analysis of its spectrum and localization based on the extended H.E.S.S. data set. The joint analysis of the extended H.E.S.S and Fermi-LAT data results in a spectral energy distribution in the energy range from 1.7 GeV to 14.8 TeV, which suggests a high luminosity of N132D at GeV and TeV energies. We set a lower limit on a gamma-ray cutoff energy of 8 TeV with a confidence level of 95%. The new gamma-ray spectrum as well as multiwavelength observations of N132D when compared to physical models suggests a hadronic origin of the VHE gamma-ray emission. Conclusions. SNR N132D is a VHE gamma-ray source that shows a spectrum extending to the VHE domain without a spectral cutoff at a few TeV, unlike the younger oxygen-rich SNR Cassiopeia A. The gamma-ray emission is best explained by a dominant hadronic component formed by diffusive shock acceleration. The gamma-ray properties of N132D may be affected by an interaction with a nearby molecular cloud that partially lies inside the 95% confidence region of the source position.

Cosmic rays, neutrinos, and GeV-TeV gamma rays from starburst galaxy NGC 4945

The detection of high-energy astrophysical neutrinos and ultra-high-energy cosmic rays (UHECRs) provides a new way to explore sources of cosmic rays. One of the highest energy neutrino events detected by IceCube, tagged as IC35, is close to the UHECR anisotropy region detected by Pierre Auger Observatory. The nearby starburst (SB) galaxy, NGC 4945, is close to this anisotropic region and inside the mean angular error of the IC35 event. Considering the hypernovae contribution located in the SB region of NGC 4945, which can accelerate protons up to $\ensuremath{\sim}{10}^{17}\text{ }\text{ }\mathrm{eV}$ and inject them into the interstellar medium, we investigate the origin of this event around this starburst galaxy. We show that the interaction of these protons with the SB region's gas density could explain Fermi-LAT gamma-ray and radio observations if the magnetic field's strength in the SB region is the order of $\ensuremath{\sim}\mathrm{mG}$. Our estimated PeV neutrino events, in ten years, for this source is approximately 0.01 ($4\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}$) if a proton spectral index of 2.4 (2.7) is considered, which would demonstrate that IC35 is not produced in the central region of this SB galaxy. Additionally, we consider the superwind region of NGC 4945 and show that protons can hardly be accelerated in it up to UHECRs.

Stochastic Fluctuations of Low-Energy Cosmic Rays and the Interpretation of Voyager Data.

Data from the Voyager probes have provided us with the first measurement of cosmic ray intensities at MeV energies, an energy range that had previously not been explored. Simple extrapolations of models that fit data at GeV energies, e.g., from AMS-02, however, fail to reproduce the Voyager data in that the predicted intensities are too high. Oftentimes, this discrepancy is addressed by adding a break to the source spectrum or the diffusion coefficient in an adhoc fashion, with a convincing physical explanation yet to be provided. Here, we argue that the discrete nature of cosmic ray sources, which is usually ignored, is instead a more likely explanation. We model the distribution of intensities expected from a statistical model of discrete sources and show that its expectation value is not representative but has a spectral shape different from that for a typical configuration of sources. The Voyager proton and electron data are however compatible with the median of the intensity distribution.

LHAASO and the origin of cosmic rays

The extension of the cosmic ray spectrum beyond 1 PeV (1 PeV = 1015 eV) indicates the existence of particle accelerators operating as PeVatrons in the Milky Way. However, recent gamma-ray observations reveal that supernova remnants (SNRs), considered the major source of galactic cosmic rays, have difficulty in accelerating cosmic rays to PeV. Therefore, searching and identifying PeVatrons is crucial for understanding the origin of cosmic rays. The ultra-high-energy (UHE) gamma-ray observations will provide decisive evidence for such studies. Due to its high sensitivity in the UHE gamma-ray band, Large High Altitude Air Shower Observatory (LHAASO) would play a key role in the research on PeVatron. After running the 1/2 LHAASO-KM2A array, 12 UHE gamma-ray sources were detected. This paper summarizes the observation results and discusses the prospect of the full LHAASO array.

Particles accelerated in Hertzian magnetic dipole field of a pulsar * *Supported by grants from NSFC (11675182, 11690022)

Supernova remnants are supposed to be the most possible sources of cosmic rays. However, alternative sources of cosmic rays, such as an active galactic nucleus, gamma-ray bursts, and pulsars, have not be excluded. In this study, we investigate the possibility of cosmic rays being generated by pulsars. The pulsar is simply described as a rotational magnetic dipole, the so-called Hertzian magnetic dipole, an exact solution of the d'Alembert equations. In the rotational magnetic dipole field, charged particles experience an accelerated electric field with their radiation reaction. The particles, which are initially static out of the light cylinder radius, can be accelerated up to a high energy.

Magnetic diffusion and interaction effects on ultrahigh energy cosmic rays: Protons and nuclei

The flux of ultrahigh energy cosmic rays reaching the Earth is affected by the interactions with the cosmic radiation backgrounds as well as with the magnetic fields that are present along their trajectories. We combine the SimProp cosmic ray propagation code with a routine that allows to account for the average effects of a turbulent magnetic field on the direction of propagation of the particles. We compute in this way the modification of the spectrum which is due to the magnetic horizon effect, both for primary nuclei as well as for the secondary nuclei resulting from the photo-disintegration of the primary ones. We also provide analytic parameterizations of the attenuation effects, as a function of the magnetic field parameters and of the density of cosmic ray sources, which make it possible to obtain the expected spectra in the presence of the magnetic fields from the spectra that would be obtained in the absence of magnetic fields. The discrete nature of the distribution of sources with finite density also affects the spectrum of cosmic rays at the highest energies where the flux is suppressed due to the interactions with the radiation backgrounds, and parameterizations of these effects are obtained.

New limit on high Galactic latitude PeV γ-ray flux from Tibet ASγ data

The Tibet ASγ collaboration has recently reported the detection of γ-rays with energies up to Peta-electronvolt from parts of the Galactic plane. We note that the analysis of γ-ray flux by the Tibet-ASγ experiment also implies an upper bound on the diffuse γ-ray flux from high Galactic latitudes (|b|> 20°) in the energy range between 100 TeV and 1 PeV. This bound is up to an order of magnitude stronger than previously derived bounds from GRAPES3, KASCADE, and CASA-MIA experiments. We discuss the new Tibet-ASγ limit on the high Galactic latitude γ-ray flux in the context of possible mechanisms of multi-messenger (γ-ray and neutrino) emission from nearby cosmic ray sources, dark matter decays, and the large-scale cosmic ray halo of the Milky Way.

The Hunt for Pevatrons: The Case of Supernova Remnants

The search for Galactic pevatrons is now a well-identified key science project of all instruments operating in the very-high-energy domain. Indeed, in this energy range, the detection of gamma rays clearly indicates that efficient particle acceleration is taking place, and observations can thus help identify which astrophysical sources can energize particles up to the ~PeV range, thus being pevatrons. In the search for the origin of Galactic cosmic rays (CRs), the PeV range is an important milestone, since the sources of Galactic CRs are expected to accelerate PeV particles. This is how the central scientific goal that is ’solving the mystery of the origin of CRs’ has often been distorted into ’finding (a) pevatron(s)’. Since supernova remnants (SNRs) are often cited as the most likely candidates for the origin of CRs, ’finding (a) pevatron(s)’ has often become ’confirming that SNRs are pevatrons’. Pleasingly, the first detection(s) of pevatron(s) were not associated to SNRs. Moreover, all clearly detected SNRs have yet revealed to not be pevatrons, and the detection from VHE gamma rays from regions unassociated with SNRs, are reminding us that other astrophysical sites might well be pevatrons. This short review aims at highlighting a few important results on the search for Galactic pevatrons.

Studying Cosmic Ray Sources Using Intergalactic Electromagnetic Cascades

In this paper, intergalactic electromagnetic cascades are used as a probe of cosmic ray sources. This is achieved as follows. In extragalactic space, cosmic rays initiate electromagnetic cascades, in which gamma-ray and neutrino emission arises. We used the joint analysis of cosmic ray data, along with extragalactic gamma-ray and neutrino emissions, to study particle acceleration in the vicinity of supermassive black holes. Particle injection spectra depend on processes of particle acceleration, and here we discuss models with various injection spectra. The computations of the propagation of cosmic rays in space were performed using the publicly available TransportCR code. It was found that a new subclass of sources might exist that does not contribute to the particle flux on Earth, instead to gamma-ray and neutrino emissions arising in electromagnetic cascades. In addition, the upper limit of the relative number of ‘exotic’ supermassive black holes surrounded by a superstrong magnetic field is derived.

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.

Cosmic-Ray Acceleration from Turbulence in Molecular Clouds

Abstract Low-energy cosmic rays, in particular protons with energies below 1 GeV, are significant drivers of the thermochemistry of molecular clouds. However, these cosmic rays are also greatly impacted by energy losses and magnetic field transport effects in molecular gas. Explaining cosmic-ray ionization rates of 10−16 s−1 or greater in dense gas requires either a high external cosmic-ray flux, or local sources of MeV–GeV cosmic-ray protons. We present a new local source of low-energy cosmic rays in molecular clouds: first-order Fermi acceleration of protons in regions undergoing turbulent reconnection in molecular clouds. We show from energetic-based arguments there is sufficient energy within the magnetohydrodynamic turbulent cascade to produce ionization rates compatible with inferred ionization rates in molecular clouds. As turbulent reconnection is a volume-filling process, the proposed mechanism can produce a near-homogeneous distribution of low-energy cosmic rays within molecular clouds.

High-energy neutrino production in clusters of galaxies

ABSTRACT Clusters of galaxies can potentially produce cosmic rays (CRs) up to very high energies via large-scale shocks and turbulent acceleration. Due to their unique magnetic-field configuration, CRs with energy ≤1017 eV can be trapped within these structures over cosmological time-scales, and generate secondary particles, including neutrinos and gamma rays, through interactions with the background gas and photons. In this work, we compute the contribution from clusters of galaxies to the diffuse neutrino background. We employ 3D cosmological magnetohydrodynamical simulations of structure formation to model the turbulent intergalactic medium. We use the distribution of clusters within this cosmological volume to extract the properties of this population, including mass, magnetic field, temperature, and density. We propagate CRs in this environment using multidimensional Monte Carlo simulations across different redshifts (from z ∼ 5 to z = 0), considering all relevant photohadronic, photonuclear, and hadronuclear interaction processes. We find that, for CRs injected with a spectral index α = 1.5–2.7 and cutoff energy Emax = 1016–5 × 1017 eV, clusters contribute to a sizeable fraction to the diffuse flux observed by the IceCube Neutrino Observatory, but most of the contribution comes from clusters with M ≳ 1014 M⊙ and redshift z ≲ 0.3. If we include the cosmological evolution of the CR sources, this flux can be even higher.

Origin of Galactic Sub-PeV Diffuse Gamma-Ray Emission: Constraints from High-energy Neutrino Observations

Abstract Very recently, diffuse gamma-rays with 0.1 PeV < E γ < 1 PeV have been discovered from the Galactic disk by the Tibet air shower array and muon detector array (Tibet AS+MD array). While the measured sub-PeV flux may be compatible with the hadronic origin in the conventional Galactic cosmic-ray propagation model, we find that it is in possible tension with the nondetection of Galactic neutrino emissions by the IceCube neutrino telescope. We further find that the presence of an extra cosmic-ray component of relatively hard spectrum, which is probably related to the Cygnus Cocoon region and other PeV cosmic-ray sources in the Galactic disk, would alleviate the tension. This scenario implies the existence of an extreme accelerator of either protons or electrons beyond PeV in the Cygnus region, and predicts the continuation of the gamma-ray spectrum of Cygnus Cocoon up to 1 PeV with a possible hardening beyond ∼30–100 TeV.

Cosmic-ray production from neutron escape in microquasar jets

Context. The origin of Galactic cosmic rays remains a matter of debate, but supernova remnants are commonly considered to be the main place where high-energy cosmic rays are accelerated. Nevertheless, current models predict cosmic-ray spectra that do not match observations and the efficiency of the acceleration mechanism is still undetermined. On the other hand, the contribution of other kinds of sources to the Galactic cosmic-ray population is still unclear, and merits investigation. Aims. In this work we explore a novel mechanism through which microquasars might produce cosmic rays. In this scenario, microquasar jets generate relativistic neutrons, which escape and decay outside the system; protons and electrons, created when these neutrons decay, escape to the interstellar medium as cosmic rays. Methods. We introduce the relativistic neutron component through a coupling term in the transport equation that governs the jet proton population. We compute the escape rate and decay distribution of these neutrons, and follow the propagation of the decay products until they escape the system and become cosmic rays. We then compute the spectra of these cosmic rays. Results. Neutrons can drain only a small fraction of the jet power as cosmic rays. The most promising scenarios arise in extremely luminous systems (Ljet ∼ 1040 erg s−1), in which the fraction of jet power deposited in cosmic rays can reach ∼0.001. Slow jets (Γ ≲ 2, where Γ is the bulk Lorentz factor) favour neutron production. The resulting cosmic-ray spectrum is similar for protons and electrons, which share the power in the ratio given by neutron decay. The spectrum peaks at roughly half the minimum energy of the relativistic protons in the jet; it is soft (spectral index ∼3) above this energy, and almost flat below. Conclusions. The proposed mechanism produces more energetic cosmic rays from microquasars than those presented by previous works in which the particles escape through the jet terminal shock. Values of spectral index steeper than 2 are possible for cosmic rays in our model and these indeed agree with those required to explain the spectral signatures of Galactic cosmic rays, although only the most extreme microquasars provide power comparable to that of a typical supernova remnant. The mechanism explored in this work may provide stronger and softer cosmic-ray sources in the early Universe, and therefore contribute to the heating and reionisation of the intergalactic medium.

Colliding-Wind Binaries as a Source of TeV Cosmic Rays

In addition to gamma-ray binaries which contain a compact object, high-energy and very high–energy gamma rays have also been detected from colliding-wind binaries. The collision of the winds produces two strong shock fronts, one for each wind, both surrounding a shock region of compressed and heated plasma, where particles are accelerated to very high energies. Magnetic field is also amplified in the shocked region on which the acceleration of particles greatly depends. In this work, we performed full three-dimensional magnetohydrodynamic simulations of colliding winds coupled to a code that evolves the kinematics of passive charged test particles subject to the plasma fluctuations. After the run of a large ensemble of test particles with initial thermal distributions, we show that such shocks produce a nonthermal population (nearly 1% of total particles) of few tens of GeVs up to few TeVs, depending on the initial magnetization level of the stellar winds. We were able to determine the loci of fastest acceleration, in the range of MeV/s to GeV/s, to be related to the turbulent plasma with amplified magnetic field of the shock. These results show that colliding-wind binaries are indeed able to produce a significant population of high-energy particles, in relatively short timescales, compared to the dynamical and diffusion timescales.