The origin of galactic cosmic rays

One century ago Viktor Hess carried out several balloon flights that led him to conclude that the penetrating radiation responsible for the discharge of electroscopes was of extraterrestrial origin. One century from the discovery of this phenomenon seems to be a good time to stop and think about what we have understood about Cosmic Rays. The aim of this review is to illustrate the ideas that have been and are being explored in order to account for the observable quantities related to cosmic rays and to summarize the numerous new pieces of observation that are becoming available. In fact, despite the possible impression that development in this field is somewhat slow, the rate of new discoveries in the last decade or so has been impressive, and mainly driven by beautiful pieces of observation. At the same time scientists in this field have been able to propose new, fascinating ways to investigate particle acceleration inside the sources, making use of multifrequency observations that range from the radio, to the optical, to X-rays and gamma rays. These ideas can now be confronted with data.I will mostly focus on supernova remnants as the most plausible sources of Galactic cosmic rays, and I will review the main aspects of the modern theory of diffusive particle acceleration at supernova remnant shocks, with special attention for the dynamical reaction of accelerated particles on the shock and the phenomenon of magnetic field amplification at the shock. Cosmic-ray escape from the sources is discussed as a necessary step to determine the spectrum of cosmic rays at the Earth. The discussion of these theoretical ideas will always proceed parallel to an account of the data being collected especially in X-ray and gamma-ray astronomy.In the end of this review I will also discuss the phenomenon of cosmic-ray acceleration at shocks propagating in partially ionized media and the implications of this phenomenon in terms of width of the Balmer line emission. This field of research has recently experienced a remarkable growth, in that Hα lines have been found to bear information on the cosmic-ray acceleration efficiency of supernova shocks.     

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.     

Propagation of galactic cosmic rays in the outer heliosphere

The propagation of galactic cosmic rays in heliospheric magnetic fields is studied. An approximate solution to the cosmic ray transport equation has been derived on the basis of a method that takes into account the small value of anisotropy of particle angular distribution. The spatial and energy distributions of the cosmic ray intensity and anisotropy have been investigated, and estimates of cosmic ray energy flux have been carried out.     

The age distribution of galactic cosmic rays

The general solution for the distribution of ages for primary galactic cosmic rays is given for a class of steady-state, bounded models of cosmic-ray diffusion in the Galaxy. Both one-and threedimensional models are considered, with point sources and distributed sources. In all models, the age distribution, is approximately exponential for ages longer than the average age, although for shorter ages the distribution depends on the details of the model. The leaky-box model, with an exponential age distribution, is thus a good approximation to most diffusive models. It is shown that one-dimensional (disk) models are consistent with both age and anisotropy data for galactic cosmic rays regardless of whether production takes place near the galactic center or throughout the disk.     

Radial gradients and anisotropies due to galactic cosmic rays

Numerical calculations have been made of the radial gradients and the anisotropyvector atr=1 AU due to galactic cosmic-ray protons and helium nuclei. The model used assumes transport by convection and anisotropic diffusion, and includes the energy losses due to adiabatic deceleration. The present calculations are for the 1964–65 solar minimum. An important constraint applied ineach case was that the model reproduces the electron modulation known from deductions of the galactic spectrum and observations of the near-Earth spectrum; and also reproduces the near-Earth proton and helium nuclei spectra. The diffusion coefficients have been based upon those deduced from magnetic-field power spectra.The principal aim has been to provide estimates of radial gradients and anisotropies, particularly at kinetic energiesT100 MeV/nucleon, by the complete solution of realistic models. Typical values for protons, obtained with a galactic differential number density (total energy)^–2.5, atT50 MeV are: radial gradient, 25%/AU; radial anisotropy, –0.2%; azimuthal anisotropy, 0.2%. These values change markedly when the galactic spectrum is cut-off or greatly enhanced atT<150 MeV, but the intensity spectrum near Earth remains substantially unchanged.It has been shown that it is possible to obtain negative radial gradients and positive radial anisotropies atT50 MeV for galactic particles and thus to mimic solar sources. The radial gradient for 1964–65 reported by Anderson (1968) and by Krimigis and Venkatesan (1969) are shown to be consistent with the diffusion coefficient deduced from the magnetic-field power spectrum; those reported by O'Gallagher are higher than expected and that for 20T30 MeV protons appears to be inconsistent. More precise data on conditions throughout the solar cavity are required if more definitive gradients and anisotropies are to be determined.     

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.     

The age distribution of galactic cosmic rays

The distribution of ages for secondary and interacting cosmic-ray nuclei is derived for a class of steady-state, bounded models of diffusion in the galaxy. Results are presented in detail for the model in which diffusion is in one dimension, sources are uniformly distributed throughout the scattering region, and particles are observed near the central plane of the galaxy. The leakage-lifetime approximation is shown to be accurate as long as the ratio of the cosmic-ray diffusive lifetime to the mean interaction or decay lifetime is less than about ten.     

Intensity of galactic cosmic rays in the early sun epoch

The process of heliospheric modulation of intensity of galactic cosmic rays is investigated by solving the transport equation. The spatial-energetic distribution of cosmic rays in the present epoch and in the past is analyzed. It is demonstrated that the particle density and the energy density of cosmic rays in the Solar System in the distant past were much lower than the corresponding current values. The cosmic ray intensity modulation in the early heliosphere was especially strong in the case of low-energy particles.     

Energy losses of galactic cosmic rays in the interplanetary medium

Using realistic models of cosmic-ray propagation in interplanetary space we present, for electrons, protons and helium nuclei of a given energy near Earth, calculations of their distribution in energy before entering the solar cavity and their mean energy loss. Interplanetary conditions appropriate for the epochs 1965 and 1969 have been used. Cosmic-ray energies in the range of 20 MeV/nucleon to 1000 MeV/nucleon have been considered.     

Energy losses and modulation of galactic cosmic rays

Numerical solutions of the cosmic-ray transport equation for the interplanetary region and including the effects of diffusion, convection, and energy loss, have been obtained for a galactic differential number-density spectrum of gaussian form, central kinetic energyT ₀, and a width at half height maximum of 0.1T ₀. These solutions have been used to show the redistribution in energy, and the reduction in number density, within the solar cavity. The mean energy loss is shown to be usefully approximated by the force-field potential when /T ₀1/2. The principal finding is that galactic cosmic-ray protons and helium nuclei, with kinetic energies less than about 80 MeV/nucleon, virtually, are completely excluded from near Earth by convective effects. As a result of this exclusion, it is found that it is not possible to demodulate the proton and helium-nuclei spectra, observed in 1964–65, to obtain information about low-energy galactic cosmic rays.     

Pulsed galactic nuclei and the origin of cosmic rays

The possibility that a series of explosions of the galactic nuclei every 5×10⁶ yr can cause a substantial flux of cosmic ray particles at the vicinity of the Earth is investigated. The steady flux of cosmic radiation forces the conclusion that there have been explosions back to 10⁹ yr if this is a dominant source of cosmic rays.     

The intensity variations of solar and galactic cosmic rays with azimuthal angle in the polar region

A balloon-born multidirectional detector is used to measure the intensity variation of galactic and solar cosmic rays with the azimuthal angle, the zenith angle being maintained at 60°. In polar regions, the intensity towards the north is found to be 20% larger than that towards the south. It is shown that this anisotropy does not originate in interplanetary space and is not produced by a magnetospheric source. It is suggested that the effect is due to propagation effects within the magnetosphere.     

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.     

Diffusion of the galactic cosmic rays in the vicinity of the solar system

Cosmic-ray propagation in the vicinity of 1 kpc from the Sun is considered. The data on the 10¹²–10¹⁵ eV particle anisotropy, on 10¹² eV electron spectrum, and on temporal cosmic-ray variations are analyzed. The diffusion coefficientD(10¹²–10¹³ eV)=10²⁹–10³⁰ cm²s^–1 inferred from the analysis coincides with its standard value in the large-halo model withh=15 kpc. The total power of cosmic-ray generation, about 3×10⁴⁹ erg per SN in the proton component and about 10⁴⁸ erg per SN in the electron component, typical of the galactic diffusion model is in agreement with the obtained parameters of local sources.     

Reconstruction of the direction of the true anisotropy of galactic cosmic rays

A method of reconstructing the declination of galactic cosmic ray anisotropy is described, and its results are presented. The method is based on analysis of delay distributions in symmetrically arranged detectors of an air shower array, and it represents a modification of the crossed telescopes method. It is shown that the declination of the true anisotropy vector is close to 60° (i.e., this vector lies approximately within the galactic plane). Because of this, the true degree of anisotropy of galactic cosmic rays is severalfold higher than the first harmonic of intensity in the sidereal time (the quantity measured directly), and it equals about 0.2%.     

Active galactic nuclei flow velocities and the highest energy cosmic rays

The dependence of the maximum energy cosmic rays can reach via diffusive shock acceleration in AGN jets on flow speeds is discussed. It is shown that in highly inclined termination shocks where the speed of the de Hoffman-Teller frame is crucial, a good independent knowledge of the jet speed is required to properly assess the extent of the cosmic ray spectrum.     

On the magnitude and direction of the anisotropy of galactic cosmic rays

The results of measuring the diurnal cosmic-ray intensity variations in the energy range 1–100 TeV are discussed. Whereas the phase of the first harmonic of the sidereal daily wave directly determines the phase (right ascension) of the cosmic-ray anisotropy vector, the amplitude and declination of the true anisotropy cannot be reconstructed directly from the amplitude of the first harmonic. However, they can be determined by invoking data on the zero harmonic. The results of some recent experiments purporting to measure the cosmic-ray anisotropy with a particularly high accuracy are shown to be interpreted erroneously.     

The peculiarities of diffusion of the galactic cosmic rays at low energies

I give an interpretation of a result of Simpsonet al. (1988) on the variation with kinetic energyT _i of the mean pathlengthX _m (T _i ) of the galactic cosmic rays (CGRs) in the range 0.1T _i 10.0 GeV nucl^–1. I argue that the data onX _m (T _i ) may be interpreted in terms of a model of GCR diffusion on the one-dimensional Alfvén-wave turbulence, having a cutoff in the spectrum at frequencies _h , where _h is the proton gyrofrequency. The cutoff results in changing of the character of variation of the GCR diffusion coefficientD(T _i )T ^a in the rangeT _i 1 GeV nucl^–1 towards some more complicated variation at 0.1T _i 1.0 GeV nucl^–1 due to the peculiarities of the pitch-angle scattering at 90⁰.