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.     

Equilibration Processes in the Warm-Hot Intergalactic Medium

The Warm-Hot Intergalactic Medium (WHIM) is thought to contribute about 40–50% to the baryonic budget at the present evolution stage of the universe. The observed large scale structure is likely to be due to gravitational growth of density fluctuations in the post-inflation era. The evolving cosmic web is governed by non-linear gravitational growth of the initially weak density fluctuations in the dark energy dominated cosmology. Non-linear structure formation, accretion and merging processes, star forming and AGN activity produce gas shocks in the WHIM. Shock waves are converting a fraction of the gravitation power to thermal and non-thermal emission of baryonic/leptonic matter. They provide the most likely way to power the luminous matter in the WHIM. The plasma shocks in the WHIM are expected to be collisionless. Collisionless shocks produce a highly non-equilibrium state with anisotropic temperatures and a large differences in ion and electron temperatures. We discuss the ion and electron heating by the collisionless shocks and then review the plasma processes responsible for the Coulomb equilibration and collisional ionisation equilibrium of oxygen ions in the WHIM. MHD-turbulence produced by the strong collisionless shocks could provide a sizeable non-thermal contribution to the observed Doppler parameter of the UV line spectra of the WHIM.     

Non-Thermal Processes in Cosmological Simulations

Non-thermal components are key ingredients for understanding clusters of galaxies. In the hierarchical model of structure formation, shocks and large-scale turbulence are unavoidable in the cluster formation processes. Understanding the amplification and evolution of the magnetic field in galaxy clusters is necessary for modelling both the heat transport and the dissipative processes in the hot intra-cluster plasma. The acceleration, transport and interactions of non-thermal energetic particles are essential for modelling the observed emissions. Therefore, the inclusion of the non-thermal components will be mandatory for simulating accurately the global dynamical processes in clusters. In this review, we summarise the results obtained with the simulations of the formation of galaxy clusters which address the issues of shocks, magnetic field, cosmic ray particles and turbulence.     

The multifuid solar wind termination shock and its influence on the threedimensional plasma structure upstream and downstream

The solar wind termination shock is described as a multi-fluid phenomenon taking into account the magnetohydrodynamic self-interaction of a multispecies plasma consisting of solar wind ions, pick-up ions and shock-generated anomalous cosmic ray particles. The spatial diffusion of these high energy particles relative to the resulting, pressure-modified solar wind flow structure is described by a coupled system of differential equations describing mass-, momentum-, and energy-flow continuities for all plasma components. The energy loss due to escape of energetic particles (MeV) from the precursor into the inner heliosphere is taken into account. We determine the integrated properties of the anomalous cosmic ray gas and the low-energy solar wind. Also the variation of the compression ratio of the shock structure is quantitatively determined and is related to the pick-up ion energization efficiency and to the mean energy of the downstream anomalous cosmic ray particles. The variation of the resulting shock structure and of the solar wind sheath plasma extent beyond the shock is discussed with respect to its consequences for the LISM neutral gas filtration and the threedimensional shape of the heliosphere.     

Microphysics of Quasi-parallel Shocks in Collisionless Plasmas

Shocks in collisionless plasmas require dissipation mechanisms which couple fields and particles at scales much less than the conventional collisional mean free path. For quasi-parallel geometries, where the upstream magnetic field makes a small angle to the shock normal direction, wave-particle coupling produces a broad transition zone with large amplitude, nonlinear magnetic pulsations playing an important role. At high Mach numbers, ion reflection and acceleration are dominant processes which control the structure and dissipation at the shock. Accelerated particles produce a precursor, or foreshock, characterized by low frequency magnetic waves which are convected by the plasma flow into the shock transition zone. The interplay between energetic particles, waves, ion reflection and acceleration leads to a complicated interdependent system. This review discusses the spacecraft observations which have motivated the current view of the high Mach number quasi-parallel shock, and the theories and simulation studies which have led to a better understanding of the microphysics on which the quasi-parallel shock depends.     

The Dynamic Quasiperpendicular Shock: Cluster Discoveries

The physics of collisionless shocks is a very broad topic which has been studied for more than five decades. However, there are a number of important issues which remain unresolved. The energy repartition amongst particle populations in quasiperpendicular shocks is a multi-scale process related to the spatial and temporal structure of the electromagnetic fields within the shock layer. The most important processes take place in the close vicinity of the major magnetic transition or ramp region. The distribution of electromagnetic fields in this region determines the characteristics of ion reflection and thus defines the conditions for ion heating and energy dissipation for supercritical shocks and also the region where an important part of electron heating takes place. In other words, the ramp region determines the main characteristics of energy repartition. All these processes are crucially dependent upon the characteristic spatial scales of the ramp and foot region provided that the shock is stationary. The process of shock formation consists of the steepening of a large amplitude nonlinear wave. At some point in its evolution the steepening is arrested by processes occurring within the shock transition. From the earliest studies of collisionless shocks these processes were identified as nonlinearity, dissipation, and dispersion. Their relative role determines the scales of electric and magnetic fields, and so control the characteristics of processes such as ion reflection, electron heating and particle acceleration. The determination of the scales of the electric and magnetic field is one of the key issues in the physics of collisionless shocks. Moreover, it is well known that under certain conditions shocks manifest a nonstationary dynamic behaviour called reformation. It was suggested that the transition from stationary to nonstationary quasiperiodic dynamics is related to gradients, e.g. scales of the ramp region and its associated whistler waves that form a precursor wave train. This implies that the ramp region should be considered as the source of these waves. All these questions have been studied making use observations from the Cluster satellites. The Cluster project continues to provide a unique viewpoint from which to study the scales of shocks. During its lifetime the inter-satellite distance between the Cluster satellites has varied from 100 km to 10000 km allowing scientists to use the data best adapted for the given scientific objective. The purpose of this review is to address a subset of unresolved problems in collisionless shock physics from experimental point of view making use multi-point observations onboard Cluster satellites. The problems we address are determination of scales of fields and of a scale of electron heating, identification of energy source of precursor wave train, an estimate of the role of anomalous resistivity in energy dissipation process by means of measuring short scale wave fields, and direct observation of reformation process during one single shock front crossing.     

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.     

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.     

The plasma physics of shock acceleration

The notion that plasma shocks in astrophysical settings can and do accelerate charged particles to high energies is not a new one. However, in recent years considerable progress has been achieved in understanding the role particle acceleration plays both in astrophysics and in the shock process itself. In this paper we briefly review the history and theory of shock acceleration, paying particular attention to theories of parallel shocks which include the backreaction of accelerated particles on the shock structure. We discuss in detail the work that computer simulations, both plasma and Monte Carlo, are playing in revealing how thermal ions interact with shocks and how particle acceleration appears to be an inevitable and necessary part of the basic plasma physics that governs collisionless shocks. We briefly describe some of the outstanding problems that still confront theorists and observers in this field.     

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.     

Shocks: Commonalities in Solar-Terrestrial Chains

Shocks are found throughout the heliosphere, wherever supersonic (or super-magnetosonic) flows encounter obstacles or other, slowly moving, media. Although some of the physical parameters are in different regimes, all shocks heat and decelerate the media incident upon them. Most shocks must propagate in a collisionless plasma, thereby adding importance to the particle interactions with the electromagnetic fields, and enabling some particles to be accelerated to high energies. This paper explores the commonalities, and differences, in shocks throughout the heliosphere, and concentrates on the role of shock microstructure in effecting the shock transition and in governing the resulting energy partition amongst the constituent species. Shocks play a significant role in the solar-terrestrial chain.     

Nonthermal Radiation Mechanisms

In this paper we review the possible radiation mechanisms for the observed non-thermal emission in clusters of galaxies, with a primary focus on the radio and hard X-ray emission. We show that the difficulty with the non-thermal, non-relativistic Bremsstrahlung model for the hard X-ray emission, first pointed out by Petrosian (Astrophys. J. 557, 560, 2001) using a cold target approximation, is somewhat alleviated when one treats the problem more exactly by including the fact that the background plasma particle energies are on average a factor of 10 below the energy of the non-thermal particles. This increases the lifetime of the non-thermal particles, and as a result decreases the extreme energy requirement, but at most by a factor of three. We then review the synchrotron and so-called inverse Compton emission by relativistic electrons, which when compared with observations can constrain the value of the magnetic field and energy of relativistic electrons. This model requires a low value of the magnetic field which is far from the equipartition value. We briefly review the possibilities of gamma-ray emission and prospects for GLAST observations. We also present a toy model of the non-thermal electron spectra that are produced by the acceleration mechanisms discussed in an accompanying paper Petrosian and Bykov (Space Sci. Rev., 2008, this issue, Chap. 11).     

Particle Acceleration Mechanisms

In this paper we review the possible mechanisms for production of non-thermal electrons which are responsible for the observed non-thermal radiation in clusters of galaxies. Our primary focus is on non-thermal Bremsstrahlung and inverse Compton scattering, that produce hard X-ray emission. We first give a brief review of acceleration mechanisms and point out that in most astrophysical situations, and in particular for the intracluster medium, shocks, turbulence and plasma waves play a crucial role. We also outline how the effects of the turbulence can be accounted for. Using a generic model for turbulence and acceleration, we then consider two scenarios for production of non-thermal radiation. The first is motivated by the possibility that hard X-ray emission is due to non-thermal Bremsstrahlung by nonrelativistic particles and attempts to produce non-thermal tails by accelerating the electrons from the background plasma with an initial Maxwellian distribution. For acceleration rates smaller than the Coulomb energy loss rate, the effect of energising the plasma is to primarily heat the plasma with little sign of a distinct non-thermal tail. Such tails are discernible only for acceleration rates comparable or larger than the Coulomb loss rate. However, these tails are accompanied by significant heating and they are present for a short time of <10⁶ years, which is also the time that the tail will be thermalised. A longer period of acceleration at such rates will result in a runaway situation with most particles being accelerated to very high energies. These more exact treatments confirm the difficulty with this model, first pointed out by Petrosian (Astrophys. J. 557:560, 2001). Such non-thermal tails, even if possible, can only explain the hard X-ray but not the radio emission which needs GeV or higher energy electrons. For these and for production of hard X-rays by the inverse Compton model, we need the second scenario where there is injection and subsequent acceleration of relativistic electrons. It is shown that a steady state situation, for example arising from secondary electrons produced from cosmic ray proton scattering by background protons, will most likely lead to flatter than required electron spectra or it requires a short escape time of the electrons from the cluster. An episodic injection of relativistic electrons, presumably from galaxies or AGN, and/or episodic generation of turbulence and shocks by mergers can result in an electron spectrum consistent with observations but for only a short period of less than one billion years.     

Plasma kinetics in the solar wind

An account is given of the observations and theoretical ideas concerning the role of kinetic processes in the solar wind. This includes, first of all, the measurements on distribution functions of plasma electrons and protons, the relation of the observed non-thermal electron features with the concept of an exospheric expansion of the solar corona, and the connection of non-thermal proton distributions with bulk flow inhomogeneities of the wind. A discussion is given of the present understanding of the connection between observed features of the particle distributions and anomalous values of some plasma transport coefficients, which in turn determine the actual values of macroscopic plasma parameters.A further topic of the review is that of possible kinetic processes occurring within small scale structures in the solar wind, like collisionless shocks, various types of discontinuities and D-sheets.     

Hybrid simulation codes with application to shocks and upstream waves

Hybrid codes in which part of the plasma is represented as particles and the rest as a fluid are discussed. In the past few years such codes with particle ions and massless, fluid electrons have been applied to space plasmas, especially to collisionless shocks. All of these simulation codes are one-dimensional and similar in structure, except for how the field equations are solved. We describe in detail the various approaches that are used (resistive Ohm's law, predictor-corrector, Hamiltonian) and compare results from the various codes with examples taken from collisionless shocks and low frequency wave phenomena upstream of shocks.     

Thermodynamical Properties of the ICM from Hydrodynamical Simulations

Modern hydrodynamical simulations offer nowadays a powerful means to trace the evolution of the X-ray properties of the intra-cluster medium (ICM) during the cosmological history of the hierarchical build up of galaxy clusters. In this paper we review the current status of these simulations and how their predictions fare in reproducing the most recent X-ray observations of clusters. After briefly discussing the shortcomings of the self-similar model, based on assuming that gravity only drives the evolution of the ICM, we discuss how the processes of gas cooling and non-gravitational heating are expected to bring model predictions into better agreement with observational data. We then present results from the hydrodynamical simulations, performed by different groups, and how they compare with observational data. As terms of comparison, we use X-ray scaling relations between mass, luminosity, temperature and pressure, as well as the profiles of temperature and entropy. The results of this comparison can be summarised as follows: (a) simulations, which include gas cooling, star formation and supernova feedback, are generally successful in reproducing the X-ray properties of the ICM outside the core regions; (b) simulations generally fail in reproducing the observed “cool core” structure, in that they have serious difficulties in regulating overcooling, thereby producing steep negative central temperature profiles. This discrepancy calls for the need of introducing other physical processes, such as energy feedback from active galactic nuclei, which should compensate the radiative losses of the gas with high density, low entropy and short cooling time, which is observed to reside in the innermost regions of galaxy clusters.     

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.     

The interplanetary medium and its interaction with the Earth's magnetosphere

This paper reviews the principal results of direct measurements of the plasma and magnetic field by spacecraft close to the Earth (within the heliocentric distance range 0.7–1.5 AU). The paper gives an interpretation of the results for periods of decrease, minimum and increase of the solar activity. The following problems are discussed: the interplanetary plasma (chemical composition, density, solar wind flow speed, temperature, temporal and spatial variation of these parameters), the interplanetary magnetic field (intensity, direction, fluctuations and its origin), some derived parameters characterizing the physical condition of the interplanetary medium; the quasi-stationary sector structure and its connection with solar and terrestrial phenomena; the magnetohydrodynamic discontinuities in the interplanetary medium (tangential discontinuities and collisionless shock waves); the solar magnetoplasma interaction with the geomagnetic field (the collisionless bow shock wave, the magnetosheath, the magnetopause, the Earth's magnetic tail, the internal magnetosphere characteristics), the connection between the geomagnetic activity and the interplanetary medium and magnetosphere parameters; peculiarities in behaviour of the interplanetary medium and magnetosphere during geomagnetic storms; energetic aspects of the geomagnetic storms.     

Selected Problems in Collisionless-Shock Physics

The physics of collisionless shocks is a very broad topic, which has been well studied for many decades. However, there are a number of important issues which remain unresolved. Moreover, there have been new findings, which cast doubt on well-established ideas. The purpose of this review is to address a subset of unresolved problems in collisionless shock physics from a theoretical and/or numerical modeling point of view. The topics which are addressed are: the nonstationarity of the shock front, the heating and dynamics of electrons through the shock layer, particle diffusion in turbulent electric and magnetic fields, particle acceleration, and the interaction of pickup ions with collisionless shocks.     

Wave induced energetic particle generation and space plasma modeling

An increasing number of high-resolution spacecraft observations provide access to details of energetic electron and ion velocity-space distribution structures. Since resonant wave-particle interaction processes depend considerably on the distribution function details, space plasma modeling is of particular interest for studies of a variety of plasma environments as planetary magnetospheres, the interplanetary medium or solar flares. After summarizing the most popular particle acceleration processes we focus on wave-powered energization mechanisms induced by Landau interaction and demonstrate from a time-evolutionary scenario that power-law distributions, highly favored by observations in recent years, are generated resonantly by an Alfvén wave spectrum and possibly saturate. This process is further stimulated in non-uniform magnetic field configurations where multiple wave packets at different phase velocities provide the energy source for a continuous acceleration process. Moreover, in this conjunction we demonstrate that in particular κ-distributions are a consequence of a generalized entropy concept, favored by nonextensive statistics, which provides the missing link for power-law plasma models from fundamental physics. With regard to in situ space observations examples are provided illuminating that for non-thermal plasma characteristics the particular structure of the velocity-space distribution dominates as regulating mechanism for the wave-particle interaction process over effects related to changes in space plasma parameters. This revised version was published online in August 2006 with corrections to the Cover Date.