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.     

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.     

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.     

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.     

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.     

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.     

Radio emission from quasi-parallel shock waves in the corona

Shock waves in the solar corona manifest themselves in type II bursts in dynamic radio spectra. Recently, short large amplitude magnetic structures (SLAMS) have been observed in the vicinity of the quasi-parallel region of Earth's bow shock as an example of a collisionless shock wave in space plasmas. SLAMS are able to accelerate electrons to high energies by shock drift acceleration. Assuming that SLAMS also appear in the vicinity of super-critical, quasi-parallel shocks in the corona, electrons can also be accelerated at quasi-parallel shocks and, subsequently, generate radio waves manifesting in solar type II radio bursts.     

Microphysics in Astrophysical Plasmas

Although macroscale features dominate astrophysical images and energetics, the physics is controlled through microscale transport processes (conduction, diffusion) that mediate the flow of mass, momentum, energy, and charge. These microphysical processes manifest themselves in key (all) boundary layers and also operate within the body of the plasma. Crucially, most plasmas of interest are rarefied to the extent that classical particle collision length- and time-scales are long. Collective plasma kinetic phenomena then serve to scatter or otherwise modify the particle distribution functions and in so-doing govern the transport at the microscale level. Thus collisionless plasmas are capable of supporting thin shocks, current sheets which may be prone to magnetic reconnection, and the dissipation of turbulence cascades at kinetic scales. This paper lays the foundation for the accompanying collection that explores the current state of knowledge in this subject. The richness of plasma kinetic phenomena brings with it a rich diversity of microphysics that does not always, if ever, simply mimic classical collision-dominated transport. This can couple the macro- and microscale physics in profound ways, and in ways which thus depend on the astrophysical context.     

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.     

Quasi-parallel collisionless shocks

The magnetic field and plasma data from the ISEE 1, 2, and 3 spacecraft have greatly increased our knowledge of the quasi-parallel collisionless shock in space. Hybrid-code simulations have provided us with valuable insights into the physics of the quasi-parallel shock. Unfortunately, theoretical understanding of the nonlinear physics of the quasi-parallel shock is still in a qualitative stage of development. Generation of large-amplitude whistler waves and hydromagnetic waves observed in the quasi-parallel shock has been discussed either in terms of linear instabilities or qualitative nonlinear arguments. It appears that the ion reflection, ion heating, and leakage of the shock-heated downstream ions at the quasi-parallel shock can all be explained in terms of nonadiabatic scatterings of ions by the large-amplitude whistler-magnetosonic waves with frequencies near the ion gyrofrequency and wavelength near the ion inertial length. The nonadiabatic scattering is defined by the non-conservation of the magnetic moment. Future study of the quasi-parallel shock should focus on developing quantitative theoretical models for the nonlinear physical processes fundamental to the quasi-parallel shock.     

Optical and UV Diagnostics of Supernova Remnant Shocks

The relatively faint optical and UV emission from non-radiative shock waves provides diagnostics for processes related to cosmic ray acceleration in collisionless shocks. Emission line profiles and intensities can be used to determine the efficiencies of electron-ion and ion-ion thermal equilibration, which influence the population of fast particles injected into the acceleration process. It is found that T _e/T _p declines with shock speed and that T _i is roughly proportional to mass in fast shocks. Important information about cosmic ray precursors may be available, but the interpretation is still somewhat ambiguous. The compression ratios in shocks which efficiently accelerate cosmic rays are predicted to be substantially larger than the factor of 4 expected for a strong shock in a = 5/3 perfect gas, and some limits may be available from observations.     

Points of controversy in the study of magnetic storms

In the present review we concentrate on the most important and so far unresolved problems concerning morphology and physics of magnetic storms, namely, 1) which currents are responsible for the Dst variation; 2) which geophysical factors control the Dst index and what the underlying physics is; 3) why the latitude of the auroral electrojets and energetic particle precipitation region decreases during storms; 4) if there actually is a storm-to-substorm relation. We present major experimental results and theoretical views on each of the above problems.     

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.     

Particle code simulations with injected particles

As problems we are interested in become more complex, we often find our simulations stretching the limits of available computer resources. For example, an interesting problem is simulation of dissipation processes in sub-critical collisionless shocks. To simulate this system our simulation box must contain the shock and its upstream and downstream regions over the entire length of a run. If the shock moves with any appreciable speed the box must then be considerably larger than the shock thickness making it hard to resolve the shock front itself with a reasonable number of grid points. A solution to this problem is to run the simulation in the frame of reference of the shock. Particles are injected upstream of the shock and leave the simulation box downstream. With the shock stationary in the simulation box, we only need to contain enough of the up and downstream regions for the fields, etc., to settle down and separate the shock from the box boundaries. In this tutorial we consider some basic algorithms used in a practical particle injection code, such as the two dimensional WAVE code used at Los Alamos. We will try to present these ideas in a simple format general enough to be easily included in any particle code. Topics covered are:

1. Smoothly Injecting Particles.
2. Generating the Distribution Functions.
3. Time Dependent Injection Density.
4. Boundary Conditions on Fields and Particles.

(Flux and Charge Conservation)     

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.     

Ion Acceleration at the Earth’s Bow Shock

The Earth’s bow shock is the most studied example of a collisionless shock in the solar system. It is also widely used to model or predict the behaviour at other astrophysical shock systems. Spacecraft observations, theoretical modelling and numerical simulations have led to a detailed understanding of the bow shock structure, the spatial organization of the components making up the shock interaction system, as well as fundamental shock processes such as particle heating and acceleration. In this paper we review the observations of accelerated ions at and upstream of the terrestrial bow shock and discuss the models and theories used to explain them. We describe the global morphology of the quasi-perpendicular and quasi-parallel shock regions and the foreshock. The acceleration processes for field-aligned beams and diffuse ion distribution types are discussed with connection to foreshock morphology and shock structure. The different possible mechanisms for extracting solar wind ions into the acceleration processes are also described. Despite several decades of study, there still remain some unsolved problems concerning ion acceleration at the bow shock, and we summarize these challenges.     

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.     

Theory and Simulation of Reconnection

Reconnection is a major commonality of solar and magnetospheric physics. It was conjectured by Giovanelli in 1946 to explain particle acceleration in solar flares near magnetic neutral points. Since than it has been broadly applied in space physics including magnetospheric physics. In a special way this is due to Harry Petschek, who in 1994 published his ground breaking solution for a 2D magnetized plasma flow in regions containing singularities of vanishing magnetic field. Petschek’s reconnection theory was questioned in endless disputes and arguments, but his work stimulated the further investigation of this phenomenon like no other. However, there are questions left open. We consider two of them – “anomalous” resistivity in collisionless space plasma and the nature of reconnection in three dimensions. The CLUSTER and SOHO missions address these two aspects of reconnection in a complementary way -- the resistivity problem in situ in the magnetosphere and the 3D aspect by remote sensing of the Sun. We demonstrate that the search for answers to both questions leads beyond the applicability of analytical theories and that appropriate numerical approaches are necessary to investigate the essentially nonlinear and nonlocal processes involved. Necessary are both micro-physical, kinetic Vlasov-equation based methods of investigation as well as large scale (MHD) simulations to obtain the geometry and topology of the acting fields and flows.     

Simulations of pickup-ion acceleration at quasi-perpendicular shocks

We present results from hybrid simulations (kinetic ion/fluid electron) of the interaction of interstellar pickup ions with collisionless shocks. Since cross-field transport is unphysically suppressed in the one-dimensional geometry used here, an ad hoc scattering algorithm is used to model this effect. This is a necessary step to accelerate the pickup ions from their initial low energies at quasi-perpendicular shocks to the high energies which are often observed associated with traveling interplanetary shocks by Ulysses.