Computer modeling of test particle acceleration at oblique shocks

We review the basic techniques and results of numerical codes used to model the acceleration of charged particles at oblique, fast-mode, collisionless shocks. The emphasis is upon models in which accelerated particles (ions) are treated as test particles, and particle dynamics is calculated by numerically integrating along exact phase-space orbits. We first review the case where ions are sufficiently energetic so that the shock can be approximated by a planar discontinuity, and where the electromagnetic fields on both sides of the shock are defined at the outset of each computer run. When the fields are uniform and static, particles are accelerated by the scatter-free drift acceleration process at a single shock encounter. We review the characteristics of scatter-free drift acceleration by considering how an incident particle distribution is modified by interacting with a shock. Next we discuss drift acceleration when magnetic fluctuations are introduced on both sides of the shock, and compare these results with those obtained under scatter-free conditions. We describe the modeling of multiple shock encounters, discuss specific applications, and compare the model predictions with theory. Finally, we review some recent numerical simulations that illustrate the importance of shock structure to both the ion injection process and to the acceleration of ions to high energies at quasi-perpendicular shocks.     

Cosmological Shock Waves

Large-scale structure formation, accretion and merging processes, AGN activity produce cosmological gas shocks. The shocks convert a fraction of the energy of gravitationally accelerated flows to internal energy of the gas. Being the main gas-heating agent, cosmological shocks could amplify magnetic fields and accelerate energetic particles via the multi-fluid plasma relaxation processes. We first discuss the basic properties of standard single-fluid shocks. Cosmological plasma shocks are expected to be collisionless. We then review the plasma processes responsible for the microscopic structure of collisionless shocks. A tiny fraction of the particles crossing the shock is injected into the non-thermal energetic component that could get a substantial part of the ram pressure power dissipated at the shock. The energetic particles penetrate deep into the shock upstream producing an extended shock precursor. Scaling relations for postshock ion temperature and entropy as functions of shock velocity in strong collisionless multi-fluid shocks are discussed. We show that the multi-fluid nature of collisionless shocks results in excessive gas compression, energetic particle acceleration, precursor gas heating, magnetic field amplification and non-thermal emission. Multi-fluid shocks provide a reduced gas entropy production and could also modify the observable thermodynamic scaling relations for clusters of galaxies.     

Collisionless Shocks in Partly Ionized Plasma with Cosmic Rays: Microphysics of Non-thermal Components

In this review we discuss some observational aspects and theoretical models of astrophysical collisionless shocks in partly ionized plasma with the presence of non-thermal components. A specific feature of fast strong collisionless shocks is their ability to accelerate energetic particles that can modify the shock upstream flow and form the shock precursors. We discuss the effects of energetic particle acceleration and associated magnetic field amplification and decay in the extended shock precursors on the line and continuum multi-wavelength emission spectra of the shocks. Both Balmer-type and radiative astrophysical shocks are discussed in connection to supernova remnants interacting with partially neutral clouds. Quantitative models described in the review predict a number of observable line-like emission features that can be used to reveal the physical state of the matter in the shock precursors and the character of nonthermal processes in the shocks. Implications of recent progress of gamma-ray observations of supernova remnants in molecular clouds are highlighted.     

Shock Acceleration of Ions in the Heliosphere

Energetic particles constitute an important component of the heliospheric plasma environment. They range from solar energetic particles in the inner heliosphere to the anomalous cosmic rays accelerated at the interface of the heliosphere with the local interstellar medium. Although stochastic acceleration by fluctuating electric fields and processes associated with magnetic reconnection may account for some of the particle populations, the majority are accelerated by the variety of shock waves present in the solar wind. This review focuses on “gradual” solar energetic particle (SEP) events including their energetic storm particle (ESP) phase, which is observed if and when an associated shock wave passes Earth. Gradual SEP events are the intense long-duration events responsible for most space weather disturbances of Earth’s magnetosphere and upper atmosphere. The major characteristics of gradual SEP events are first described including their association with shocks and coronal mass ejections (CMEs), their ion composition, and their energy spectra. In the context of acceleration mechanisms in general, the acceleration mechanism responsible for SEP events, diffusive shock acceleration, is then described in some detail including its predictions for a planar stationary shock, shock modification by the energetic particles, and wave excitation by the accelerating ions. Finally, some complexities of shock acceleration are addressed, which affect the predictive ability of the theory. These include the role of temporal and spatial variations, the distinction between the plasma and wave compression ratios at the shock, the injection of thermal plasma at the shock into the process of shock acceleration, and the nonlinear evolution of ion-excited waves in the vicinity of the shock.     

Electron Acceleration in the Heliosphere

We review the evidence for electron acceleration in the heliosphere putting emphasis on the acceleration processes. There are essentially four classes of such processes: shock acceleration, reconnection, wave particle interaction, and direct acceleration by electric fields. We believe that only shock and electric field acceleration can in principle accelerate electrons to very high energies. The shocks known in the heliosphere are coronal shocks, traveling interplanetary shocks, CME shocks related to solar type II radio bursts, planetary bow shocks, and the termination shock of the heliosphere. Even in shocks the acceleration of electrons requires the action of wave particle resonances of which beam driven whistlers are the most probable. Other mechanisms of acceleration make use of current driven instabilities which lead to electron and ion hole formation. In reconnection acceleration is in the current sheet itself where the particles perform Speiser orbits. Otherwise, acceleration takes place in the slow shocks which are generated in the reconnection process and emanate from the diffusion region in the Petschek reconnection model and its variants. Electric field acceleration is found in the auroral zones of the planetary magnetospheres and may also exist on the sun and other stars including neutron stars. The electric potentials are caused by field aligned currents and are concentrated in narrow double layers which physically are phase space holes in the ion and electron distributions. Many of them add up to a large scale electric field in which the electrons may be impulsively accelerated to high energies and heated to large temperatures.     

Fe/O Ratios in Interplanetary Shock Accelerated Particles

It is widely accepted that diffusive shock acceleration is an important process in the heliosphere, in particular in producing the energetic particles associated with interplanetary shocks driven by coronal mass ejections. In its simplest formulation shock acceleration is expected to accelerate ions with higher mass to charge ratios less efficiently than those with lower mass to charge. Thus it is anticipated that the Fe/O ratio in shock-accelerated ion populations will decrease with increasing energy above some energy. We examine the circumstances of five interplanetary shocks that have been reported to have associated populations in which Fe/O increases with increasing energy. In each event, the situation is complex, with particle contributions from other sources in addition to the shock. Furthermore, we show that the Fe/O ratio in shock-accelerated ions can decrease even when the shock is traveling through an Fe-rich ambient ion population. Thus, although shock acceleration of an Fe-rich suprathermal population has been proposed to explain large Fe-rich solar particle events, we find no support for this proposal in these observations.     

CME Theory and Models

This chapter provides an overview of current efforts in the theory and modeling of CMEs. Five key areas are discussed: (1) CME initiation; (2) CME evolution and propagation; (3) the structure of interplanetary CMEs derived from flux rope modeling; (4) CME shock formation in the inner corona; and (5) particle acceleration and transport at CME driven shocks. In the section on CME initiation three contemporary models are highlighted. Two of these focus on how energy stored in the coronal magnetic field can be released violently to drive CMEs. The third model assumes that CMEs can be directly driven by currents from below the photosphere. CMEs evolve considerably as they expand from the magnetically dominated lower corona into the advectively dominated solar wind. The section on evolution and propagation presents two approaches to the problem. One is primarily analytical and focuses on the key physical processes involved. The other is primarily numerical and illustrates the complexity of possible interactions between the CME and the ambient medium. The section on flux rope fitting reviews the accuracy and reliability of various methods. The section on shock formation considers the effect of the rapid decrease in the magnetic field and plasma density with height. Finally, in the section on particle acceleration and transport, some recent developments in the theory of diffusive particle acceleration at CME shocks are discussed. These include efforts to combine self-consistently the process of particle acceleration in the vicinity of the shock with the subsequent escape and transport of particles to distant regions.     

Stochastic Acceleration by Turbulence

The subject of this paper is stochastic acceleration by plasma turbulence, a process akin to the original model proposed by Fermi. We review the relative merits of different acceleration models, in particular the so called first order Fermi acceleration by shocks and second order Fermi by stochastic processes, and point out that plasma waves or turbulence play an important role in all mechanisms of acceleration. Thus, stochastic acceleration by turbulence is active in most situations. We also show that it is the most efficient mechanism of acceleration of relatively cool non relativistic thermal background magnetized plasma particles. In addition, it can preferentially accelerate electrons relative to protons as is needed in many astrophysical radiating sources, where usually there are no indications of presence of shocks. We also point out that a hybrid acceleration mechanism consisting of initial acceleration by turbulence of background particles followed by a second stage acceleration by a shock has many attractive features. It is demonstrated that the above scenarios can account for many signatures of the accelerated electrons, protons and other ions, in particular ³He and ⁴He, seen directly as Solar Energetic Particles and through the radiation they produce in solar flares.     

Magnetic Fields, Relativistic Particles, and Shock Waves in Cluster Outskirts

It is only now, with low-frequency radio telescopes, long exposures with high-resolution X-ray satellites and γ-ray telescopes, that we are beginning to learn about the physics in the periphery of galaxy clusters. In the coming years, Sunyaev-Zel’dovich telescopes are going to deliver further great insights into the plasma physics of these special regions in the Universe. The last years have already shown tremendous progress with detections of shocks, estimates of magnetic field strengths and constraints on the particle acceleration efficiency. X-ray observations have revealed shock fronts in cluster outskirts which have allowed inferences about the microphysical structure of shocks fronts in such extreme environments. The best indications for magnetic fields and relativistic particles in cluster outskirts come from observations of so-called radio relics, which are megaparsec-sized regions of radio emission from the edges of galaxy clusters. As these are difficult to detect due to their low surface brightness, only few of these objects are known. But they have provided unprecedented evidence for the acceleration of relativistic particles at shock fronts and the existence of μG strength fields as far out as the virial radius of clusters. In this review we summarise the observational and theoretical state of our knowledge of magnetic fields, relativistic particles and shocks in cluster outskirts.     

Particle Acceleration at Interplanetary Shocks

This paper briefly reviews proton acceleration at interplanetary shocks. This is key to describing the acceleration of heavy ions at interplanetary shocks because wave excitation—and hence particle scattering—at oblique shocks is controlled by the protons and not the heavy ions. Heavy ions behave as test particles, and their acceleration characteristics are controlled by the properties of proton-excited turbulence. As a result, the resonance condition for heavy ions introduces distinctly different signatures in abundance, spectra, and intensity profiles, depending on ion mass and charge. Self-consistent models of heavy-ion acceleration and the resulting fractionation are discussed. This includes discussion of the injection problem and the acceleration characteristics of quasi-parallel and quasi-perpendicular shocks.     

Perpendicular Shock Reformation and Ion Acceleration

Kinetic simulations of supercritical, quasi-perpendicular shocks yield time-varying solutions that cyclically reform on proton spatio-temporal scales. Whether a shock solution is stationary or reforming depends upon the plasma parameters which, for SNR shocks and the heliospheric termination shock, are ill defined but believed to be within this time-dependent regime. We first review the time-dependent solutions and the acceleration processes of the ions for a proton–electron plasma. We then present recent results for a three-component plasma: background protons, electrons and a second ion population appropriate for SNR (heavy ions) or the termination shock (pickup protons). This ion acceleration generates a suprathermal “injection” population – a seed population for subsequent acceleration at the shock, which may in turn generate ions at cosmic ray energies.     

Acceleration of particles by shocks in a cosmic plasma

The theory and observational evidence pertaining to particle acceleration by shock waves in astrophysical objects and in space are systematized. Recent works showing observational and theoretical aspects of the problem dealing with shocks in turbulent media are emphasized. The acceleration of particles by shocks in turbulent media is observed in interplanetary space. This acceleration mechanism is of particular interest from the point of view of the origin of cosmic rays, providing the degree form of the spectrum. The index of the spectrum is close to the observable one for galactic cosmic rays. It depends slightly on specific conditions in the acceleration region. Electron and nucleus acceleration in supernova remnants and in radiogalaxies is discussed, and theory and observational data are compared. The theory of particle acceleration by supersonic turbulence is outlined.     

Stochastic Acceleration of Energetic Particles in the Heliosphere

In this paper we assess possible roles of stochastic acceleration by random electric field and plasma motion in the production and transport of energetic particles in the heliosphere. Stochastic acceleration can occur in the presence of multiple small-scale magnetohydrodynamic waves propagating in different directions. Usually, this type of stochastic acceleration is closely related to particle pitch angle scattering or parallel diffusion. Given the values of the parallel diffusion coefficient inferred from the observations of cosmic ray modulation or other energetic particle phenomena in the heliosphere, stochastic acceleration by small-scale waves is much slower than acceleration by shock waves and it is also much slower than adiabatic cooling by the expansion of the solar wind; thus it is considered as inefficient for producing heliospheric energetic particles or for the modulation of cosmic rays. Another type of stochastic acceleration occurs when particles go through random compressions or expansions due to large-scale plasma motion. This acceleration mechanism could be very fast when the correlation time of the fluctuations in plasma compression is short compared to the diffusion time. Particle acceleration by an ensemble of small shock waves or intermittent long wavelength compressible turbulence belongs to this category. It tends to establish an asymptotic p ^?3 universal distribution function quickly if there is no or little large-scale adiabatic cooling. Such a particle distribution will contain an infinite amount of pressure. Back reaction from the pressure is expected to modify the amplitude of plasma waves to an equilibrium state. At that point, the pressure of accelerated particles must remain finite and the accelerated particles could approach a p ^?5 distribution function.     

Shock acceleration of energetic particles in the heliosphere

The theory of shock acceleration of energetic particles is briefly discussed and reviewed with an emphasis on clarifying the apparent distinction between the V × B and Fermi mechanisms. Attention is restricted to those situations in which the energetic particles do not themselves influence the given shock structure. In particular, application of the theory to the acceleration of energetic particles in corotating interaction regions (CIR) in the solar wind is presented. Here particles are accelerated at the forward and reverse shocks which bound the CIR by being compressed between the shock fronts and magnetic irregularities upstream from the shocks, or by being compressed between upstream irregularities and those downstream from the shocks. Particles also suffer adiabatic deceleration in the expanding solar wind, an effect not included in previous shock models for acceleration in CIRs. The model is able to account for the observed exponential spectra at Earth, the observed behavior of the spectra with radial distance, the observed radial gradients in the intensity, and the observed differences in the intensity and spectra at the forward and reverse shocks.Calculations and resulting energy spectra are also presented for shock acceleration of energetic particles in large solar flare events. Based on the simplifying assumption that the shock evolves as a spherically symmetric Sedov blast wave, the calculation yields the time-integrated spectrum of particles initially injected at the shock which eventually escape ahead of the shock into interplanetary space. The spectra are similar to those observed at Earth. Finally further applications are suggested.An invited paper presented at STIP Workshop on Shock Waves in the Solar Corona and Interplanetary Space, 15–19 June, 1980, Smolenice, Czechoslovakia.     

Particle Acceleration in the Heliosphere: Implications for Astrophysics

There has been a remarkable discovery concerning particles that are accelerated in the solar wind. At low energies, in the region where the particles are being accelerated, the spectrum of the accelerated particles is always the same: when expressed as a distribution function, the spectrum is a power law in particle speed with a spectral index of ?5, and a rollover at higher particle speeds that can often be described as exponential. This common spectral shape cannot be accounted for by any conventional acceleration mechanism, such as diffusive shock acceleration or traditional stochastic acceleration. It has thus been necessary to invent a new acceleration mechanism to account for these observations, a pump mechanism in which particles are pumped up in energy through a series of adiabatic compressions and expansions. The conditions under which the pump acceleration is the dominant acceleration mechanism are quite general and are likely to occur in other astrophysical plasmas. In this paper, the most compelling observations of the ?5 spectra are reviewed; the governing equation of the pump acceleration mechanism is derived in detail; the pump acceleration mechanism is applied to acceleration at shocks; and, as an illustration of the potential applicability of the pump acceleration mechanism to other astrophysical plasmas, the pump mechanism is applied to the acceleration of galactic cosmic rays in the interstellar medium.     

Ion acceleration at shocks in interplanetary space: A brief review of recent observations

The mechanism by which ions are accelerated near the Earth's bow shock and near shocks propagating outward from the Sun in response to solar activity appears to be essentially the same. For both types of shock the solar wind thermal distribution acts as a seed population. Leaked magnetospheric ions and resident flare ions are additional seed populations for the bow shock and outward propagating shocks respectively. The acceleration of solar wind ions at these shocks begins with either the reflection of ions off the shock or leakage of shocked plasma back through the shock. Interaction with a disruption wave field self-generated by these backstreaming ions is responsible for the remainder of the acceleration at the bow shock. Both the disruption wave field and the ambient interplanetary wave field play important roles in accelerating ions at outward propagating shocks, but on different time scales. The geometry of the shock and the duration of field line connection to the shock play decisive roles in determining what is observed.     

Particle Acceleration and Nonthermal Phenomena in Superbubbles

Models of nonthermal particle acceleration in the vicinity of active star forming regions are reviewed. We discuss a collective effect of both stellar winds of massive stars and core collapsed supernovae as particle acceleration agents. Collective supernova explosions with great energy release in the form of multiple interacting shock waves inside the superbubbles are argued as a favourable site of nonthermal particle acceleration. The acceleration mechanism provides efficient creation of a nonthermal nuclei population with a hard low-energy spectrum, containing a substantial part of the kinetic energy released by the winds of young massive stars and supernovae. We discuss a model of temporal evolution of particle distribution function accounting for the nonlinear effect of the reaction of the accelerated particles on the shock turbulence inside the superbubble. The model illustrates that both the low-energy metal-rich nonthermal component and the standard galactic cosmic rays could be efficiently produced by superbubbles at different evolution stages.     

Particle acceleration in solar flares and some related phenomena

Observational features concerning solar energetic particles are compactly reviewed with some emphasis on the spectra and time histories. Velocity dependent characteristics in the energy spectra are pointed out, and compared to the results of the interplanetary shocks. A shock drift acceleration is introduced in order to interpret the observational features, especially a very fast acceleration to MeV energies within an order of second. There is a strong evidence of the shock drift acceleration in the interplanetary shocks. When some conditions are satisfied in the corona, only one or several encounters of particles with a near perpendicular shock accelerates protons to gamma-ray emitting energies (> 10 MeV). Pre-acceleration is inevitable for any kind of acceleration mechanisms in solar flares. To fulfill the requirements from the abundance ratios between various species of accelerated ions, pre-acceleration to the same velocities before the injection into a main acceleration process turns out to be absolutely necessary.     

Particle acceleration at the Sun and in the heliosphere

Energetic particles are accelerated in rich profusion at sites throughout the heliosphere. They come from solar flares in the low corona, from shock waves driven outward by coronal mass ejections (CMEs), from planetary magnetospheres and bow shocks. They come from corotating interaction regions (CIRs) produced by high-speed streams in the solar wind, and from the heliospheric termination shock at the outer edge of the heliospheric cavity. We sample many populations near Earth, but can distinguish them readily by their element and isotope abundances, ionization states, energy spectra, angular distributions and time behavior. Remote spacecraft have probed the spatial distributions of the particles and examined new sources in situ. Most acceleration sources can be ‘seen’ only by direct observation of the particles; few photons are produced at these sites. Wave-particle interactions are an essential feature in acceleration sources and, for shock acceleration, new evidence of energetic-proton-generated waves has come from abundance variations and from local cross-field scattering. Element abundances often tell us the physics of the source plasma itself, prior to acceleration. By comparing different populations, we learn more about the sources, and about the physics of acceleration and transport, than we can possibly learn from one source alone. This revised version was published online in August 2006 with corrections to the Cover Date.     

Heliospheric Lessons for Galactic Cosmic-ray Acceleration

The heliosphere is bathed in the supersonic solar wind, which generally creates shocks at any obstacles it encounters: magnetic structures such as coronal mass ejections and planetary magnetospheres, or fast-slow stream interactions such as corotating interaction regions (CIRs) or the termination shock. Each of these shock structures has an associated energetic particle population whose spectra and composition contain clues to the acceleration process and the sources of the particles. Over the past several years, the solar wind composition has been systematically studied, and the long-standing gap between high energy (>1 MeV amu^–1) and the plasma ion populations has been closed by instruments capable of measuring the suprathermal ion composition. In CIRs, where it has been possible to observe all the relevant populations, it turns out that the suprathermal ion population near 1.8–2.5 times the solar wind speed is the seed population that gets accelerated, not the bulk particles near the solar wind peak. These new results are of interest to the problem of Galactic Cosmic-Ray (GCR) Acceleration, since the injection and acceleration of GCRs to modest energies is likely to share many features with processes we can observe in detail in the heliosphere.