Cosmic Rays; High Energy Astrophysics; Acceleration of Particles; Magnetic Turbulence Ultra High Energy Cosmic Rays (UHECRs) hit the Earth’s atmosphere with energies exceeding eV. This is the same energy as carried by a tennis ball moving at 100 km/h, but concentrated on a sub-atomic particle. UHECRs are so rare (the flux of particles with eV is /km/century) that only a few such particles have been detected over the past 50 years.
Recently, the HiRes and Auger experiments have reported the discovery of a high-energy cut-off in the UHECR spectrum, and Auger has found an apparent clustering of the highest energy events towards nearby active galactic nuclei. Consensus is building that the highest energy particles are accelerated within the radio-bright lobes of these objects, but it remains unclear how this actually happens, and whether the cut-off is due to propagation effects or reflects an intrinsically physical limitation of the acceleration process.
The low event statistics presently allows for many different plausible models; nevertheless observations are beginning to impose strong constraints on them. These observations have also motivated suggestions that new physics may be implicated. We present a review of the key theoretical and observational issues related to the processes of propagation and acceleration of UHECRs and proposed solutions.
On the UHECRs acceleration]On the acceleration of
Ultra-High-Energy Cosmic Rays
F. Fraschetti]Federico Fraschetti
Hess (1912) discovered cosmic rays by measuring ionisation at increasing altitudes in a balloon flight. Auger et al. (1938) proved the existence of extended air showers of secondary particles caused by the interaction of primary particles, having energies above eV (eV J), with the Earth’s atmosphere. Auger simultaneously observed the arrival of secondary particles in ground detectors set many meters apart. After almost 90 years of research, the origin of cosmic rays and the physical mechanism by which they are accelerated has remained an open question.
The Universe is filled with the remnant radiation emitted by the Big Bang, the Cosmic Microwave Background (CMB) radiation. Since we do not observe degradation of Ultra High Energy Cosmic Rays (UHECRs) due to the interaction with CMB, their sources must be located within 75 Mpc ( Mpc pc and pc cm).
The supernova remnant, i.e. the material ejected in the interstellar space by the explosion of a massive star, is believed to be a strong candidate for the source of medium-energy cosmic rays, although direct evidence for this is still lacking. Moreover, the processes of acceleration of particles in the supernova remnant cannot overcome the threshold of eV since higher energy particles cannot be confined in such objects (Lagage & Cesarsky 1983). Therefore the sources of UHECRs must be located elsewhere, very likely in nearby extragalactic sources.
This review is intended to cover theoretical and observational aspects of the acceleration mechanisms which lead to Ultra High Energy (UHE), without entering the mathematical details. Inevitably not all works have been properly discussed or mentioned; for a more complete review see for instance Torres & Anchordoqui (2004) and Bhattacharjee & Sigl (2000).
2 Review of observations of UHECRs
The cosmic ray spectrum extends over 13 orders of magnitude in energy and more than 30 orders of magnitude in observed flux (see figure 1a). The cosmic ray flux at energy eV is estimated to be about particle/kmcentury. Qualitatively, the current experiments agree that the overall cosmic rays spectrum can be divided into the following four regions: 1) below eV it is flat and the particles are decelerated by the magnetized plasma ejected by the Sun’s surface; 2) between eV and eV (the so-called “knee”) it can be fitted by a unique power law with index and is likely to be dominated by galactic supernova remnants; 3) between eV and eV (the so-called “ankle”) the galactic contribution of primary cosmic rays heavier than protons, such as He or Fe, may be boosted with acceleration in other environments, such as pulsar winds, namely wind of relativistic particles emitted by a rapidly rotating and pulsating star; 4) at energies eV, in the UHE region, some experiments have provided indication for a slope around , and the steepening of the galactic component would progressively be drowned by a flatter extra-galactic flux. In the knee region, KASCADE data have confirmed (2008) an overlap of heavier elements, such as Fe. A similarly non negligible contribution of heavy nuclei has been found in UHE region by The Auger Collaboration (2008).
At an energy greater than eV, the cosmic rays impinging on the Earth’s atmosphere produce extended showers; from a primary cosmic ray, the shower spreads into a number of detectable secondary particles over a diameter of 100 meters. At eV this diameter has increased to 5 km (see figure 1 b). From the particle densities measured in the detectors, the energy of the extended shower, and therefore of the primary cosmic ray, can be inferred. The main difficulty in the analysis of such observations is that the energy of the collision of primaries with the Earth’s atmosphere (in the center of mass frame) is well beyond the limits of the particle accelerators, even the Large Hadron Collider (LHC) at CERN (acronym standing in French for Conseil Européen pour la Recherche Nucléaire). Hence the reconstruction of the showers relies on phenomenological extrapolation of models of particle interactions, and uncertainties on the nature of the primaries grow with the energy. In the near future, new measurements at the LHC will perhaps help the development of models of these interactions.
In the last twenty years or so, observations have been made by the following experiments: Haverah Park in UK, Akeno-AGASA in Japan, covering 100 km, Yakutsk in north-east Siberia, High Resolution Fly’s Eye, in Dugway, Utah, and Auger South in western Argentina, spread over an area of 3000 km, similar in size to Rhode Island.
Haverah Park and Akeno measure the density of the secondary particles in the shower at ground level at a distance of 0.6 km from the core. The Yakutsk experiment uses this density technique but calibrated by an array of wide angle photomultiplier detectors that observe the Cerenkov light, namely the radiation emitted by particles propagating through a medium at a velocity greater than the velocity of light in that medium. HiRes measures the energy from the fluorescence of the showers in the atmosphere, by observing the flash of nitrogen molecules which de-excite after the passage of the shower. The Auger South Observatory (see figure 1, b) uses both these complementary techniques to measure the extended showers on the ground with 1600 detectors spaced 1.5 km apart and in the atmosphere using four fluorescence telescope arrays. The latest observations by The Auger Collaboration (2007) have shown that the sources of cosmic rays with energy eV are correlated with the matter distribution within 75 Mpc, specifically Active Galactic Nuclei (AGN). This correlation does not identify unequivocally the nearby AGNs as the sources, but rather indicates that any astrophysical class of objects with a local spatial distribution sufficiently similar to that of AGNs is a possible accelerator of UHECRs.
The GZK cut-off (Greisen 1966, Zatsepin & Kuz’min 1966) is the upper limit on the energy that a cosmic ray proton can retain while propagating through the CMB. This should be evident in the cosmic ray energy spectrum if the sources of UHECRs are homogeneously distributed through the Universe. Recent observations by The Auger Collaboration (2008b) The HiRes Collaboration (2008)HiRes (2008) and Auger (2008) have reported the discovery of a cut-off at UHE in the cosmic ray spectrum. Whether the spectral steepening is a manifestation of the long-sought GZK cut-off or an intrisic limitation of the acceleration process, such as the geometry or the nature of the magnetic turbulence of the source, still needs to be understood. Moreover, difficulties in the definition of the GZK cut-off lie in the model dependence of the shape of the cosmic ray spectrum beyond the ankle: as an example, a local overdensity of UHECR sources, assumed to be isotropically distributed within a distance of a few tens of Mpc, drastically changes the spectrum in the GZK cut-off region (Berezinski et al. 2006). A precise definition of the GZK cut-off is also affected by the fact that the details of the spectral shape in the UHE range depend on the nature of the primaries, as the cut-off of heavy nuclei such as Fe is expected at lower energies than that of protons.
3 Propagation of high energy particles
Beyond the ankle ( eV), the energy of a microscopic particle becomes macroscopic, i.e. eV is the same amount of energy as carried by a tennis ball moving at 100 km . Thus the existence of such high-energy particles, in itself, represents a major theoretical challenge and an explanation of the mechanism of acceleration based on first principles needs to be found.
The length scale characterizing the motion of a particle of energy , electric charge in a magnetic field is the gyration radius , or gyroradius: ; a greater gyroradius implies a less curved trajectory. If the UHE is attained through an acceleration process, a general estimate of the maximal energy can be obtained by requiring the gyroradius of the particle to be no larger than the linear size of the accelerator:
if is measured in T, in eV, then is measured in kpc. Equation (1) is also called the Hillas criterion (Hillas (1984)). It neglects the finite lifetime of the acceleration region and energy losses due to interactions with the environment such as synchrotron radiation in the magnetic field and production of secondary particles. A more general expression for the maximal energy is
with , where and are respectively the electric and magnetic fields. The formula (2) can be derived from dimensional arguments, representing the electromotive force associated with the time variation of magnetic field. Charged particles are accelerated by the electric field; however cosmic plasmas do not allow high voltages to be maintained. With regard to cosmic rays at eV, there are no known astrophysical environments where voltages high enough can be generated. However, time variation of a strong magnetic field can generate electromotive forces which through a strong induced electric field may lead to extremely high energies.
Even if only qualitative, equation (2) provides an interesting criterion to identify possible sources of UHECR by simply looking at the largest values of the product . It is remarkable (see figure 2) that the possible accelerators range from Neutron stars ( km), namely remnants that can result from the gravitational collapse of a massive star in a supernova explosion, up to galaxy clusters ( Mpc).
We should keep in mind that at a given energy the gyroradius is larger for smaller charge of the particle. Therefore, if UHECRs are mostly protons, due to their large gyration radius, they are not deviated significantly by magnetic fields, so they should point back to their sources within an angle that depends on the intensity of the intergalactic magnetic field. For heavy nuclei, the effect of the magnetic field becomes more important.
Theoretical models can be divided in two classes: the so-called ”bottom-up”, in which a specific process of acceleration in a particular astrophysical object leads to UHEs (see the review by Torres & Anchordoqui 2004); and the so-called ”top-down”, where UHE particles are produced by the decay of super-heavy dark matter particles or collision among cosmic strings or topological defects (see the review by Bhattacharjee & Sigl 2000).
4 Astrophysical models
In this section we give a short presentation of some of the most promising astrophysical models of acceleration to UHE. Owing to the low statistics involved it is difficult to favour one model over another: many of these models are equally plausible. Only upcoming observations will allow us to discriminate between them. These models may be divided in two classes according to the acceleration process, either magnetohydrodynamical (), MHD, or stochastic ().
4.1 Non-stochastic models
UHE can be attained through the drift of charged particles in magnetized, cylindrically collimated and ultrarelativistic flows of plasma emitted by active galaxies (Lyutikov & Ouyed 2007). If the velocity of the flow increases towards the axis of the jet, there will be a radial induced electric field which can accelerate particles. The acceleration efficiency is inversely proportional to the charge, thus protons are expected to dominate over heavy nuclei in UHECR composition. In order for the particles to reach an energy of eV, a minimal luminosity W is required; therefore only very high power radio-galaxies are the favourite accelerator candidates. This mechanism needs a seed population of pre-accelerated particles at energy eV at the base of the centre of the galaxy or even outside the jet.
Recently it has been proposed that supermassive black holes (with mass , where kg is the mass of the Sun) might also be accelerators up to UHE (Neronov et al. 2007). If the axis of magnetic field is misaligned with the rotation axis of the black hole, charged particles, initially spiralling into the black hole along the equatorial plane, penetrate into polar regions of the black hole and are ejected to infinity with very high energies. It is a necessary condition in this scenario that the magnetic field is almost aligned to the rotation axis of the black hole.
It has also been speculated (Boldt & Ghosh 1999) that presently radio-quiet quasars, containing black holes with mass in the centre of galaxies located within a distance of Mpc, could be the accelerators to UHE. An external magnetic field, generated by the external currents flowing in the surrounding disc, threads the black hole; thus, an electromotive force is induced. The loss of energy of the accelerated particles would allow only a limited number of particles to reach energies eV, but this model is not yet falsified by observations.
Young neutron stars with surface magnetic field in the range T could accelerate iron nuclei through relativistic hydrodynamical winds in a magnetic field up to energy eV (Blasi et al. 2000). The iron nuclei can pass through the surrounding supernova remnant without suffering significant deceleration and thus be observed. However, the discovery of the high-energy cut-off disfavours this model which predicts observations above the cut-off.
4.2 Fermi second-order mechanism
The original idea of explaining the acceleration of high energy cosmic rays by stochastic collisions is due to Fermi (1949) who imagined that a particle can be accelerated through random collisions with the irregularities of the interstellar magnetic field. The energy source is therefore the time variation of the intergalactic magnetic field which, through the induced electric field, accelerates particles. This is the only possible origin of acceleration because a constant magnetic field, with any possible orientation, only affects the trajectory but does not do any work on the particle. On the other hand, in intergalactic space any background electric field vanishes, because the ionized gas in the interstellar medium is highly conductive and therefore globally neutral.
As shown by Fermi, the average energy gain is given by , where is the mean non-relativistic velocity of the magnetic irregularities in units of the velocity of light in vacuum, m s. It is remarkable that the energy gain increases as , depends only on the velocity of the clouds, and not on the magnetic field strength and/or the electric charge of the particle.
Now we consider a number of particles continuously injected in the magnetically turbulent region described previously. The second-order Fermi mechanism yields a spectral distribution in the form of a power law: , where dE is the number of particles with energy between and dE. It is remarkable that, if the spectrum has a power law dependence such as in a certain energy range, the total power emitted also has the form of a power law in the corresponding energy range. This property is at the basis of the widespread use of the power law distribution of ultrarelativistic electrons in many branches of high energy astrophysics.
However, a few issues affect the original model of Fermi. First of all, the random velocity of the magnetic clouds is of order of m s; thus, the average gain of energy is very small, that is . Since the average distance between clouds can be estimated to be of order of pc, with one collision per year the energy of a particle increases only of a factor 3 in one billion years. We know currently that the average age of cosmic rays is Myr, therefore this mechanism in its original form is too slow to attain UHEs. Second, the value of is not unequivocally predicted by the theory, since it has a complicated dependence on the propagation properties of the particles. The third issue concerns the fact that the particle injected at low energy can lose its energy in coulombian interactions; in order to avoid this, a lower limit on the initial energy was already inferred by Fermi: 200 MeV for protons, but higher for heavy nuclei as Si or Fe.
In a modern version of Fermi’s mechanism, the clouds of turbulent magnetic field can be considered an idealization of the clumps of matter homogeneously distributed in the galactic space moving at Alfven wave velocity , given by , where is the strength of background magnetic field and is the matter density of background plasma. For mercury at room temperature [m s][T], while the sound speed is m s. In astrophysical environments, can, however, be much larger owing to the low density: in the Sun’s photosphere the density is of order of kg m, so that [m s][T]. Since in this scenario the scatterings are much more frequent than in the original Fermi model, the acceleration is much more efficient.
4.3 Acceleration at shocks
More popular and efficient than Fermi second-order is the acceleration at the non-relativistic shock preceding the expansion of matter flowing at speeds larger than the speed of sound in that medium, e.g. the remnant of supernovae. This is the so-called first-order Fermi process (Blandford & Eichler 1987). In this case, the average energy gain is , where is roughly the shock velocity. The basic idea is that the second-order process is not very efficient because the energy gain results as the small positive balance between the gains in the head-on collisions with the clumps of matter and the losses in the overtaking collisions. By contrast, in environments such as supernova shocks the collision of a particle with a magnetically turbulent plasma, extended both upstream or downstream, is always head-on.
The first-order mechanism is more efficient than the second-order one as long as the velocity of the shock is larger than . At late times in the supernova expansion, when the kinetic energy has been mainly dissipated in the expansion or sapped by the escape of high energy particles, the shock velocity becomes comparable to and the second-order effect can even dominate over the first-order one. Moreover, if the shock wave is relativistic, the first-order Fermi acceleration seems to be unable to attain ultra-high energies (see Ostrowski, in press, and references therein).
4.4 Radio lobes of AGN
As confirmed by the latest Auger observations (2007), AGN and their extensions are possible candidate sites of cosmic ray acceleration. The stochastic mechanism proposed by Fermi deserves to be explored in more detail. An attempt in this direction has been recently proposed by Fraschetti & Melia (2008). The three-dimensional motion of individual particles is followed during their propagation within a turbulent magnetic field. The turbulence is generated by the prescription of Giacalone & Jokipii (1994), modified to allow for temporal variations.
In the treatment by Fraschetti & Melia (2008), the remaining unknowns are the partition of energy between turbulent and background fields, and the turbulent energy distribution, though this may reasonably be assumed to be isotropic (Kolmogorov).
Furthermore, this acceleration mechanism is sustained over 10 orders of magnitude in particle energy, up to eV and the UHECRs therefore emerge naturally from the physical conditions thought to be prevalent within AGN radio lobes. The theoretically predicted spectrum is in accordance with the recent observations if only the closest sources are considered (Gelmini et al. 2007).
A new cut-off consisting of the size of the acceleration region is introduced in addition to or substitution for the GZK cut-off. This result reopens the question of whether the high-energy cut-off is due to interaction with CMB, or it represents the maximum energy permitted by the acceleration process itself.
5 Non-astrophysical models
The difficulty of finding a physical mechanism of acceleration based on simple physical principles and the distance limit posed by the interaction with the CMB have lead to the development of other non-astrophysical models for the origin of UHECRs. In these models, the acceleration problem is circumvented by assuming the existence of new particles or new interactions beyond the standard model of elementary particles. Within the standard model the only particle whose path is not limited by the interaction with the CMB is the neutrino, which may have been produced as secondary particle by a primary proton accelerated to UHE in an AGN environment at a distance larger than Mpc. The neutrino-nucleon interaction in Earth’s atmosphere may be enhanced by possible new physics, but this is unlikely to be able to simulate hadronic interactions (Anchordoqui et al. 2003).
By contrast, in the so-called “top-down” models (Bhattacharjee & Sigl 2000), the UHE particles are the decay product of some currently unknown supermassive “X” particle with mass eV, and energy up to . Therefore no acceleration mechanism needs to be searched for. The X particles can be produced in two ways: if they have a short lifetime, they have to be produced continuously by certain topological defects left over from a cosmological phase transition that occurred in the early Universe and their sources are evenly distributed in the Universe. On the other hand, if they have been produced only in the early Universe, they would be distributed like the local dark matter in the halo of our Galaxy and we would observe the UHECRs coming from their slow decays.
These non-astrophysical solutions to the UHECRs problem have the advantage of connecting with new ideas beyond the standard model of particle physics, such as grand unification and supersymmetry, as well as affording a possible probe of conditions in the early Universe. However, by their very nature these are very speculative ideas (Sarkar 2004).
In such “top-down” models the arrival directions of UHECRs are not expected to cluster on small angular regions of the sky, as in fact the recent Auger observations indicate. Moreover, The Auger Collaboration (2008) has recently excluded all such models as possible sources of UHECRs at energy eV by demonstrating that the primaries are not photons as expected in such models.
6 New perspectives
The importance of elucidating the mechanism responsible for the acceleration of UHECRs is confirmed by the growing number of dedicated experiments which will join Auger South: Auger North, with an exposure times larger than the Auger South’s; the JEM-EUSO mission, looking downward from the International Space Station to the dark side of the Earth; the proposed orbiting wide angle light collector (OWL), etc. The forthcoming data on UHECRs from Auger along with the other observatories will significantly boost the statistics and drive the UHECR source identification. We should be able to investigate the nature of the spectral steepening and address the question of whether the cut-off in the cosmic ray distribution is indeed due to photo-meson interaction or the result of limitations in the acceleration process itself. The primary particles have been shown not to be photons, thus favouring the bottom-up scenario with respect to top-down models. Presently the inferences concerning the composition of UHECR rely on hadronic interaction models that extrapolate the current particle physics phenomenology up to energies well below the LHC region. It is quite likely that further investigations on the nature of UHECRs will have consequences for our understanding of the physics of fundamental interactions.
Acknowledgements.The author thanks F. Melia for useful and stimulating discussions and S. Sarkar for the revision of the manuscript. The work was supported by CNES (French Space Agency) and was carried out at CEA/Saclay and LUTh.
- Anchordoqui, L., Paul, T., Reucroft, S. & Swain, J., 2003 Ultrahigh Energy Cosmic Rays: The state of the art before the Auger Observatory Int. J. Mod. Phys. A18, 2229–2366.
- Auger, P., Maze, R. & Grivet-Meyer, T. 1938 Grandes gerbes cosmiques atmospériques contenant des corpuscules ultra-pénétrantes Comptes rendus Acad. Sci. 206, 1721–1722.
- Berezinsky, V, Gazizov, A. & Grigorieva, S., 2006 On astrophysical solution to ultrahigh energy cosmic rays Phys. Rev. D 74 043005, 1–35.
- Bhattacharjee, P. & Sigl, G. 2000 Origin and propagation of extremely high-energy cosmic rays Phys. Rep. 327, 109–247.
- Blandford, R. D. & Eichler, D. 1987 Particle acceleration at astrophysical shocks: A theory of cosmic ray origin Astroph. J. 154 1–75.
- Blasi, P., Epstein, R.I. & Olinto, A. V. 2000 Ultra-high-energy cosmic rays from young neutron star winds Astroph. J. 533 L123–L126.
- Boldt, E. & Ghosh, P. 1999 Cosmic rays from remnants of quasars? Mon. Not. R. Astron. Soc. 307, 491–494.
- Fermi, E. 1949 On the Origin of the Cosmic Radiation. Phys. Rev. 75, 1169–1174.
- Fraschetti, F. & Melia F., 2008 Ultra-High-Energy Cosmic Rays from the Radio Lobes of AGNs, MNRAS, accepted.
- Gelmini, G., Kalashev, O. & Semikoz, D. V. 2007 GZK Photons in the Minimal Ultra High Energy Cosmic Rays Model (available at http://arxiv.org/abs/astro-ph/0702464).
- Giacalone, J. & Jokipii, J. R., 1994 Charged-particle motion in multidimensional magnetic-field turbulence Astrophys. J. 430, L137–L140.
- Greisen, K. 1966 End to the Cosmic-Ray Spectrum? Phys. Rev. Lett. 16, 748–750.
- Hess, V. F. 1912 Über Beobachtungen der durchdringenden Strahlung bei sieben Freiballonfahrten. Phys. Zeitschr. 13, 1084–1091.
- Hillas, A. M. 1984 The Origin of Ultra-High-Energy Cosmic Rays Ann. Rev. Astron. Astrophys. 22, 425–444.
- Lagage, P. O. & Cesarsky C. J. 1983 The maximum energy of cosmic rays accelerated by supernova shocks, Astron. Astrophys. 125, 249–257.
- Lyutikov, M. & Ouyed, R. 2007 Inductive acceleration of UHECRs in sheared relativistic jets Astropart. Phys. 27, 473–489
- Neronov A. YU., Semikoz D.V. & Tkachev I.I., 2007 Ultra-High Energy Cosmic Ray production in the polar cap regions of black hole magnetospheres (available at http://arxiv.org/abs/0712.1737).
- Ostrowski, M. 2008 Cosmic ray acceleration at relativistic shocks, shear layers,… New Astron. Rev. in press, (available at http://arXiv.org/abs/astro-ph/0801.1339).
- The Auger Collaboration, 2007 Correlation of the Highest-Energy Cosmic Rays with Nearby Extragalactic Objects Science 318, 939–943.
- The Auger Collaboration, 2008 Upper limit on the cosmic-ray photon flux above eV using the surface detector of the Pierre Auger Observatory Astropart. Phys. 29, 243–256.
- The Auger Collaboration, 2008 Observation of the suppression of the flux of cosmic rays above 4x eV (available at http://arxiv.org/pdf/0806.4302).
- The HiRes collaboration 2008 First Observation of the Greisen-Zatsepin-Kuzmin Suppression Phys. Rev. Lett. 100, 101101, 1–5.
- The KASCADE collaboration 2008 Applying Shower Development Universality to KASCADE Data (available at http://arxiv.org/pdf/0804.4274v1)
- Sarkar, S. 2004 New physics from ultrahigh energy cosmic rays Acta Phys.Polon. B35 351–364.
- Torres, D.F. & Anchordoqui, L.A. 2004 Astrophysical origins of ultrahigh energy cosmic rays Rep. Prog. Phys. 67, 1663-1730.
- Zatsepin, G. T. & Kuz’min, V. A. 1966 Upper Limit of the Spectrum of Cosmic Rays Soviet Journal of Experimental and Theoretical Physics Letters 4, 78–80.