Ballistic and diffusive components in the dynamic spectra of ultrahigh energy cosmic rays from nearby transient sources

Ultrahigh energy cosmic rays (UHECRs, E > 10¹⁸ eV) from extragalactic sources deviate in the galactic and intergalactic magnetic fields, which explains the diffusive character of their propagation, the isotropization of their total flux, and the absence of UHECR clusters associated with individual sources. Extremely high energy cosmic rays (E > 10^19.7 eV) are scattered mainly in localized magnetized structures, such as galaxy clusters, filaments, etc., with a mean free path of tens of megaparsecs; therefore, in the case of nearby transient sources, a substantial contribution to the observed flux is expected from unscattered and weakly scattered particles, which may be a decisive factor in the identification of these sources. We propose a method for calculating the time evolution of the UHECR energy spectra based on analytical solutions of the transport equation with the explicit determination of the contributions from scattered and unscattered particles. As examples, we consider the cases of transient activity of the nearest active galactic nucleus, Centaurus A, and the acceleration of UHECRs by a young millisecond pulsar.     

The possibility of investigating ultra-high-energy cosmic-ray sources using data on the extragalactic diffuse gamma-ray emission

We provide our estimates of the intensity of the gamma-ray emission with an energy near 0.1 TeV generated in intergalactic space in the interactions of cosmic rays with background emissions. We assume that the cosmic-ray sources are pointlike and that these are active galactic nuclei. The following possible types of sources are considered: remote and powerful ones, at redshifts up to z = 1.1, with a monoenergetic particle spectrum, E = 10²¹ eV; the same objects, but with a power-law particle spectrum; and nearby sources at redshifts 0 < z ≤ 0.0092, i.e., at distances no larger than 50 Mpc also with a power-law particle spectrum. The contribution of cosmic rays to the extragalactic diffuse gammaray background at an energy of 0.1 TeVhas been found to depend on the type of sources or, more specifically, the contribution ranges from f ? 10^?4 to f ≈ 0.1, depending on the source model. We conclude that the data on the extragalactic background gamma-ray emission can be used to determine the characteristics of extragalactic cosmic-ray sources, i.e., their distances and the pattern of the particle energy spectrum.     

Contribution of Cosmic Rays from Sources with a Monoenergetic Proton Spectrum to the Extragalactic Diffuse Gamma-Ray Emission

The extragalactic sources of ultra-high-energy (E > 4 × 10¹⁹ eV) cosmic rays that make a small contribution to the flux of particles recorded by ground-based arrays are discussed. We show that cosmic rays from such sources can produce a noticeable diffuse gamma-ray flux in intergalactic space compared to the the data obtained with Fermi LAT (onboard the Fermi space observatory). A possible type of active galactic nuclei (AGNs) in which cosmi-ray protons can be accelerated to energies 10²¹ eV is considered as an illustration of such sources. We conclude that ultra-high-energy cosmic rays from the AGNs being discussed can contribute significantly to the extragalactic diffuse gamma-ray emission. In addition, a constraint on the fraction of the AGNs under consideration relative to the BL Lac objects and radio galaxies has been obtained from a comparison with the Fermi LAT data.     

Formation of the spectrum of ultra-high-energy cosmic rays: Effect of extragalactic magnetic fields and expected contribution of individual sources

The effect of the extragalactic magnetic field on the propagation of ultra-high-energy cosmic rays (UHECRs) is investigated. We use the infrared galaxy catalog IRAS PSCz to reconstruct the magnetic field distribution in the Local Universe. The magnetic field induction is considered as a power function of the galactic infrared luminosity density: B = Kρ^β. In contrast to some earlier studies in which the exponent β = 2/3 corresponded to the freezing-in condition, the parameters K and β are estimated from the field inductions normalized by the expected maximum inductions (strong field) and minimum inductions (weak field) in galaxy clusters and voids, respectively. Maps of angular deflections of UHECRs are presented for these magnetic field models. We found that the protons with energies E > 4 × 10¹⁹ eV are not significantly deflected from their sources in a sphere with a radius of 100 Mpc only in the case of the weak magnetic field model (the deflections are comparable to the errors of modern detectors). The effect of the extragalactic magnetic field on the UHECR spectrum is investigated, with Virgo A and Arp 299 taken as potential sources.     

Composition of Cosmic Rays with <Emphasis Type="Italic">E</Emphasis><Superscript>0</Superscript> ⩾ 10<Superscript>17</Superscript> eV Based on Data from the Ground Scintillation Detectors of the Yakutsk EAS Array

The lateral distribution of cascade particles in extensive air showers from cosmic rays with energy E⁰ ? 10¹⁷ eV has been studied at the Yakutsk array by the ground scintillation detectors over the period of continuous observations 1977–2017. The experimental data are compared with those computed with various models for the development of extensive air showers from the CORSIKA software package. The best agreement between the theory and experiment is observed for the QGSjet-01-d model. In the energy range (1?20)× 10¹⁷ eV there is a change in the cosmic-raymass composition from 〈lnA〉 ≈ 2.5 to the proton one.     

Extragalactic cosmic rays and the galactic magnetic field

The origin and behavior of cosmic rays in the Galaxy depends crucially upon whether the galactic magnetic field has a closed topology, as does the field of Earth, or whether a major fraction of the lines of force connect into extragalactic space. If the latter, then cosmic rays could be of extragalactic origin, or they could be of galactic origin, detained in the Galaxy by the scattering offered by hydromagnetic waves, etc. If, on the other hand, the field is largely closed, then cosmic rays cannot be of extragalactic origin (at least below 10¹⁶ eV). They must be of galactic origin and escape because their collective pressure inflates the galactic field and they push their way out.This paper examines the structure of a galactic field that opens initially into intergalactic space and, with the inclusion of turbulent diffusion, finds no possibility for maintaining a significant magnetic connection with an extragalactic field. Unless some mechanism can be found, we are forced to the conclusion that the field is closed, that cosmic rays are of galactic origin, and that cosmic rays escape from the Galaxy only by pushing their way out.     

Pulsars and the origin of cosmic rays

The capabilities and limitations of pulsars as sources of cosmic rays are reviewed in the light of experimental observations. Pulsars can supply the cosmic ray power if they have rotational velocities in excess of 700 rad s^?1 at birth. Though this is theoretically possible, there is no experimental proof for the same. Pulsars can accelerate particles to the highest energies of 10²⁰ eV, but in general, the spectra on simple considerations, turn out to be flatter than the observed cosmic ray spectrum. At the highest energies, absorption processes due to fragmentation and photodisintegration dominate for heavy nuclei. The existence of a steady flux of cosmic rays of energy greater than 10¹⁷ eV demands acceleration of particles to last over fifty years, the time interval between supernovae outbursts, whereas the expected period of activity is less than a few years. Finally, the problem of anisotropy with relevance to pulsars as sources and the possibility of observing pulsar accelerated particles from galactic clusters is considered.     

Gamma-ray bursters as sources of cosmic rays

From the little we know of the physical conditions in γ-ray bursters, it seems that they are potentially effective in the acceleration of high-energy cosmic rays (CRs), especially if the bursters are at cosmological distances. We find that, with the observed statistics and fluxes of γ-ray bursts, cosmological bursters may be an important source of cosmic rays in two regions of the observed spectrum: (1) At the very-high-energy end (E > 10¹⁹ eV), where CRs must be of extragalactic origin. (2) Around and above the spectral feature that has been described as a bump and/or a knee, which occurs around 10¹⁵ eV. The occasional bursters that occur inside the Galaxy — about once in a few hundred thousand years if burst emission is isotropic; more often, if it is beamed — could maintain the density of galactic cosmic rays at the observed level in this range. These two energy ranges might correspond to two typical CR energy scales characteristic of bursters: one pertinent to CR acceleration due to interaction of a magnetized-fireball front with an ambient medium; the other to acceleration in the fireball itself (e.g. shock acceleration).     

Cosmic neutrinos from the sources of galactic and extragalactic cosmic rays

Although kilometer-scale neutrino detectors such as IceCube are discovery instruments, their conceptual design is very much anchored to the observational fact that Nature produces protons and photons with energies in excess of 10²⁰ eV and 10¹³ eV, respectively. The puzzle of where and how Nature accelerates the highest energy cosmic particles is unresolved almost a century after their discovery. From energetics considerations we anticipate on the order of 10–100 neutrino events per kilometer squared per year pointing back at the source(s) of both galactic and extragalactic cosmic rays. In this context, we discuss the results of the AMANDA and IceCube neutrino telescopes which will deliver a kilometer-square-year of data over the next 3 years.     

Extragalactic UHE proton spectrum and prediction for iron-nuclei flux at 10–10 GeV

We investigate the problem of transition from galactic cosmic rays to extragalactic ultra-high energy cosmic rays. Using the model for extragalactic ultra-high energy cosmic rays and observed all-particle cosmic ray spectrum, we calculate the galactic spectrum of iron nuclei in the energy range 10⁸–10⁹ GeV. The flux and spectrum predicted at lower energies agree well with the KASCADE data. The transition from galactic to extragalactic cosmic rays is distinctly seen in spectra of protons and iron nuclei, when they are measured separately. The shape of the predicted iron spectrum agrees with the Hall diffusion.     

Multi messenger astronomy and CTA: TeV cosmic rays and electrons

Cosmic rays are a sample of solar, galactic and extragalactic matter. Their origin and properties are one of the most intriguing question in modern astrophysics. The most energetic events and active objects in the Universe: supernovae explosion, pulsars, relativistic jets, active galactic nuclei, have been proposed as sources of cosmic rays although unambiguous evidences have still to be found. Electrons, while comprising ∼1% of the cosmic radiation, have unique features providing important information regarding the origin and propagation of cosmic rays in the Galaxy that is not accessible from the study of the cosmic-ray nuclear components due to their differing energy-loss processes. In this paper we will analyse, discussing the experimental uncertainties and challenges, the most recent measurements on cosmic-ray nuclei and, in particular, electrons with energies from tens of GeV into the TeV region.     

Centaurus A as a source of extragalactic cosmic rays with arrival energies well beyond the GZK cutoff

The ultra-high energy cosmic rays recently detected by several air shower experiments could have an extragalactic origin. In this case, the nearest active galaxy Centaurus A might be the source of the most energetic particles ever detected on Earth. We have used recent radio observations in order to estimate the arrival energy of the protons accelerated by strong shock fronts in the outer parts of this southern radio source. We expect detections coresponding to particles with energies up to 2.2 × 10²¹ eV and an arrival direction of (l ≈ 310°, b ≈ 20°) in galactic coordinates. The future Southern Hemisphere Pierre Auger Observatory might provide a decisive test for extragalactic models of the origin of the ultra-high energy cosmic rays.     

Cosmic ray acceleration in hypernova remnants expanding in wind bubbles of Wolf-Rayet progenitor stars

The determination of the origin of cosmic rays with observed energies in excess of 10¹⁷ eV that exceed the expected energies of cosmic rays accelerated by supernova remnants in the galaxy is a pressing problem in modern astrophysics. Hypernova remnants are one of the possible galactic sources of cosmic rays with energies of up to 10¹⁹ eV. Hypernovae constitute a class of extremely powerful supernova explosions, whose supposed progenitors are massive Wolf-Rayet stars. We analyze the special aspects of acceleration of cosmic rays in hypernova remnants that expand in wind bubbles of Wolf-Rayet progenitor stars. We show that these cosmic rays may attain maximum energies of 10¹⁸ eV even with a relatively conservative choice of acceleration parameters and account for tens of percent of the total cosmic ray flux observed in the vicinity of the earth in the energy range of 10¹⁶–10¹⁸ eV if the galactic hypernova explosion rate in the modern epoch reaches ? _S ～ 10^?4 year^?1.     

Phenomenological explanation for the shape of the spectrum of cosmic rays with energies <Emphasis Type="Italic">E</Emphasis> > 10 GeV

Assuming that the energy gain by cosmic-ray (CR) particles is a stochastic process with stationary increments, we derive expressions for the shape of their energy spectrum up to energies E ～ 10¹⁸ eV. In the ultrarelativistic case under study, the energy is proportional to the momentum, whose time derivative is the force. According to the Fermi mechanism, a particle accelerates when it passes through a system of shock waves produced by supernova explosions. Since these random forces act on time scales much shorter than the particle lifetime, we assume them to be delta-correlated in time. In this case, due to the linear energy-momentum relationship, the mean square of the energy (increments) is proportional to the differential scale τ(E) ～ Eτ(≥E), where τ (≥E) is the cumulative time it takes for a particle to gain an energy ≥E. The probability of finding a particle with energy ≥E somewhere in the system is inversely proportional to the time it takes to gain the energy E. To estimate an upper limit for the space number density of CR particles, we use estimates of the CR volume energy density, a quantity known for our Galaxy. It is taken to be constant in the range 10 GeV ≤ E ≤ 3 × 10⁶ GeV, where the index of the energy spectrum was found to be ?8/3 ≈ ?2.67 against its empirical value of ?2.7. In the range 3 × 10⁶ GeV ≤ E < 10⁹ GeV, the upper limit for the volume energy density is estimated by using the results from the previous range to be ?28/9 ≈ ?3.11 against its empirical value of ?3.1. The numerical coefficients in the suggested shapes of the spectrum can be determined by comparison with observational data. Thus, the CR energy spectrumis the result of the random walks of ultrarelativistic particles in energy/momentum space caused by the Fermi mechanism.     

Ultra-high energy cosmic rays in the galactic corona

On the basis of recent new information on regular and chaotic magnetic fields in coronae of spiral galaxies, we discuss propagation of ultra-high energy cosmic rays of energies exceeding 10¹⁷ eV in the galactic corona. It is shown that the expected regular magnetic field is able to confine to the corona protons of energies up to 3×10¹⁹ eV. Chaotic magnetic fields of the corona play an important role in dynamics of cosmic-ray protons of energy up to 7×10¹⁸ eV.     

Ground level enhancements of cosmic rays in solar cycle 24

Using data from ground-based observations of cosmic rays (CRs) on the worldwide network of stations and spacecraft, we have investigated the proton spectra and the CR anisotropy during the ground level enhancements of CRs on May 17, 2012 (GLE71) and January 6, 2014 (GLE72) occurred in solar cycle 24 by the spectrographic global survey method. We provide the CR rigidity spectra and the relative changes in the intensity of CRs with a rigidity of 2 GV in the solar–ecliptic geocentric coordinate system in specific periods of these events. We show that the proton acceleration during GLE71 and GLE72 occurred up to rigidities R ~ 2.3?2.5 GV, while the differential rigidity spectra of solar CRs are described neither by a power nor by an exponential function of particle rigidity. At the times of the events considered the Earth was in a loop-like structure of the interplanetary magnetic field.     

A possible universal origin of hadronic cosmic rays from ultrarelativistic ejecta of bipolar supernovae

Based on the “cannonball model” for gamma-ray bursts of Dar and De Rújula, it is proposed that masses of baryonic plasma (“cannonballs”), ejected in bipolar supernova explosions in our Galaxy are the sources of hadronic Galactic cosmic rays (CRs) at all energies. The propagation of the cannonballs in the Galactic disk and halo is studied. Two mechanisms for the acceleration of the observed CRs are proposed. The first is based on ultrarelativistic shocks in the interstellar medium and could accelerate the bulk of CRs up to the “knee” energy of 4×10¹⁵ eV. The second operates with second-order Fermi acceleration within the cannonball. If the total initial energy of the ejected plasmoids in a SN explosion is 10⁵³ erg or higher, this second mechanism may explain the CR spectrum above the knee up to the highest observed energies. It is shown that together with plausible assumptions about CR propagation in the Galactic confinement volume, the observed spectral indices of the CR spectrum can be theoretically understood to first order. The model allows a natural understanding of various basic CR observations like the absence of the Greisen–Zatsepin cutoff, the anisotropy of arrival directions as function of energy and the small Galactocentric gradient of the CR density.     

Contribution of nuclei accelerated by gamma-ray pulsars to cosmic rays in the Galaxy

We consider the galactic population of gamma-ray pulsars as possible sources of cosmic rays at and just above the “knee” in the observed cosmic ray spectrum at 10¹⁵–10¹⁶ eV. We suggest that iron nuclei may be accelerated in the outer gaps of pulsars, and then suffer partial photo-disintegration in the non-thermal radiation fields of the outer gaps. As a result, protons, neutrons, and surviving heavier nuclei are injected into the expanding supernova remnant. We compute the spectra of nuclei escaping from supernova remnants into the interstellar medium, taking into account the observed population of radio pulsars.

Our calculations, which include a realistic model for acceleration and propagation of nuclei in pulsar magnetospheres and supernova remnants, predict that heavy nuclei accelerated directly by gamma-ray pulsars could contribute about 20% of the observed cosmic rays in the knee region. Such a contribution of heavy nuclei to the cosmic ray spectrum at the knee can significantly increase the average value of lnA with increasing energy as is suggested by recent observations.     

Galactic PeV neutrinos

The IceCube experiment has detected two neutrinos with energies between 1 and 10 PeV. They might have originated from Galactic or extragalactic sources of cosmic rays. In the present work we consider hadronic interactions of the diffuse very high energy cosmic rays with the interstellar matter within our Galaxy to explain the PeV neutrino events detected in IceCube. We also expect PeV gamma ray events along with the PeV neutrino events if the observed PeV neutrinos were produced within our Galaxy in hadronic interactions. PeV gamma rays are unlikely to reach us from sources outside our Galaxy due to pair production with cosmic background radiation fields. We suggest that in future with simultaneous detections of PeV gamma rays and neutrinos it would be possible to distinguish between Galactic and extragalactic origins of very high energy neutrinos.