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.     

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.     

Clusters of Galaxies: Setting the Stage

Clusters of galaxies are self-gravitating systems of mass ∼10¹⁴–10¹⁵ h ⁻¹ M_⊙ and size ∼1–3h ⁻¹ Mpc. Their mass budget consists of dark matter (∼80%, on average), hot diffuse intracluster plasma (≲20%) and a small fraction of stars, dust, and cold gas, mostly locked in galaxies. In most clusters, scaling relations between their properties, like mass, galaxy velocity dispersion, X-ray luminosity and temperature, testify that the cluster components are in approximate dynamical equilibrium within the cluster gravitational potential well. However, spatially inhomogeneous thermal and non-thermal emission of the intracluster medium (ICM), observed in some clusters in the X-ray and radio bands, and the kinematic and morphological segregation of galaxies are a signature of non-gravitational processes, ongoing cluster merging and interactions. Both the fraction of clusters with these features, and the correlation between the dynamical and morphological properties of irregular clusters and the surrounding large-scale structure increase with redshift. In the current bottom-up scenario for the formation of cosmic structure, where tiny fluctuations of the otherwise homogeneous primordial density field are amplified by gravity, clusters are the most massive nodes of the filamentary large-scale structure of the cosmic web and form by anisotropic and episodic accretion of mass, in agreement with most of the observational evidence. In this model of the universe dominated by cold dark matter, at the present time most baryons are expected to be in a diffuse component rather than in stars and galaxies; moreover, ∼50% of this diffuse component has temperature ∼0.01–1 keV and permeates the filamentary distribution of the dark matter. The temperature of this Warm-Hot Intergalactic Medium (WHIM) increases with the local density and its search in the outer regions of clusters and lower density regions has been the quest of much recent observational effort. Over the last thirty years, an impressive coherent picture of the formation and evolution of cosmic structures has emerged from the intense interplay between observations, theory and numerical experiments. Future efforts will continue to test whether this picture keeps being valid, needs corrections or suffers dramatic failures in its predictive power.     

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.     

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.     

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.     

Linking Primordial to Solar and Galactic Composition

Evolution and composition of baryonic matter is influenced by the evolution of other forms of matter and energy in the universe. At the time of primordial nucleosynthesis the universal expansion and thus the decrease of the density and temperature of baryonic matter were controlled by leptons and photons. Non-baryonic dark matter initiated the formation of clusters and galaxies, and to this day, dark matter largely determines the dynamics and geometries of these baryonic structures and indirectly influences their chemical evolution. Chemical analyses and isotopic abundance measurements in the solar system established the composition in the protosolar cloud (PSC). The abundances of nuclear species in the PSC led to the discovery of the magic numbers and the nuclear shell model, and they allowed the identification of nucleosynthetic sites and processes. To this day, we know the abundances of the ∼300 stable and long-lived nuclides infinitely better in the PSC than in any other sample of matter in the universe. Thus, we know the exact composition of a Galactic sample of intermediate age, allowing us to check on theories of Galactic evolution before and after the formation of the solar system. This paper specifically discusses the nucleosynthesis in the early universe and the Galactic evolution during the last 5 Gyr.     

Numerical Simulations of the Warm-Hot Intergalactic Medium

In this paper we review the current predictions of numerical simulations for the origin and observability of the warm hot intergalactic medium (WHIM), the diffuse gas that contains up to 50 per cent of the baryons at z∼0. During structure formation, gravitational accretion shocks emerging from collapsing regions gradually heat the intergalactic medium (IGM) to temperatures in the range T∼10⁵–10⁷ K. The WHIM is predicted to radiate most of its energy in the ultraviolet (UV) and X-ray bands and to contribute a significant fraction of the soft X-ray background emission. While O vi and C iv absorption systems arising in the cooler fraction of the WHIM with T∼10⁵–10^5.5 K are seen in FUSE and Hubble Space Telescope observations, models agree that current X-ray telescopes such as Chandra and XMM-Newton do not have enough sensitivity to detect the hotter WHIM. However, future missions such as Constellation-X and XEUS might be able to detect both emission lines and absorption systems from highly ionised atoms such as O vii, O viii and Fe xvii.     

Two-fluid 2.5D MHD Simulations of the Fast Solar Wind in Coronal Holes and the Relation to UVCS Observations

Recent SOHO/UVCS observations indicate that the perpendicular proton and ion temperatures are much larger than electron temperatures. In the present study we simulate numerically the solar wind flow in a coronal hole with the two-fluid approach. We investigate the effects of electron and proton temperatures on the solar wind acceleration by nonlinear waves. In the model the nonlinear waves are generated by Alfvén waves with frequencies in the 10^-3 Hz range, driven at the base of the coronal hole. The resulting electron and proton flow profile exhibits density and velocity fluctuations. The fluctuations may steepen into shocks as they propagate away from the sun. We calculate the effective proton temperature by combining the thermal and wave velocity of the protons, and find qualitative agreement with the proton kinetic temperature increase with height deduced from the UVCS Ly-α observations by Kohl et al. (1998). This revised version was published online in June 2006 with corrections to the Cover Date.     

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.     

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.     

SUMER Observations of Coronal-Hole Temperatures

Observations of emission lines in the vacuum-ultraviolet spectral range with calibrated instrumentation provide crucial information on the prevailing plasma temperatures in the solar atmosphere. Coronal-hole temperatures measured by the SUMER spectrometer on SOHO will be presented in this contribution. Electron temperatures can be estimated from the formation temperatures of the observed emission lines. Line-ratio and emission-measure analyses, however, offer higher accuracies. Typical electron temperatures at altitudes of H<200?Mm in coronal holes are below 1?MK in bright structures—the coronal plumes—with higher values in darker areas—the inter-plume regions. Line-width measurements yield effective ion temperatures, which are much higher than the electron temperatures. Observations of line profiles emitted from species with different masses allow a separation of the effective temperatures into ion temperatures and unresolved non-thermal motions along the line of sight.     

Non-Linear Growth of Trapped Particle Modes in Linearly Stable, Current-Carrying Plasmas – A Fundamental Process in Plasma Turbulence and Anomalous Transport

The two-stream instability as a fundamental process in a current-carrying plasma is reconsidered. Its well-established linear version, based on kinetic Landau theory, predicts a threshold for the drift velocity between both species below which the plasma should be stable. We report on simulations which, however, show that a plasma as a non-linearly responding medium can be destabilized well below this threshold. Responsible for this unexpected behaviour are coherent, electrostatic, trapped particle structures such as phase space vortices or holes which can grow non-linearly out of thermal noise receiving their energy from the net imbalance of loss of electron kinetic energy and gain of ion kinetic energy. The birth of predominantly zero-energy holes is shown numerically being associated with initial, non-topological fluctuations. The latter are not subject to Landau damping, as they lie outside the realm of linear wave theory. For a pair plasma a typical scenario is presented, which encompasses several regimes such as non-linear growth of multiple holes, saturation and fully developed structural turbulence as well as an asymptotic approach to a new collisionless equilibrium. During the transient, structural state the plasma transport appears to be highly anomalous.     

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.     

Electron-Ion Temperature Equilibration in Collisionless Shocks: The Supernova Remnant-Solar Wind Connection

Collisionless shocks are loosely defined as shocks where the transition between pre-and post-shock states happens on a length scale much shorter than the collisional mean free path. In the absence of collision to enforce thermal equilibrium post-shock, electrons and ions need not have the same temperatures. While the acceleration of electrons for injection into shock acceleration processes to produce cosmic rays has received considerable attention, the related problem of the shock heating of quasi-thermal electrons has been relatively neglected. In this paper we review the state of our knowledge of electron heating in astrophysical shocks, mainly associated with supernova remnants (SNRs), shocks in the solar wind associated with the terrestrial and Saturnian bowshocks, and galaxy cluster shocks. The solar wind and SNR samples indicate that the ratio of electron temperature, (T _e ) to ion temperature (T _p ) declining with increasing shock speed or Alfvén Mach number. We discuss the extent to which such behavior can be understood on the basis of waves generated by cosmic rays in a shock precursor, which then subsequently damp by heating electrons, and speculate that a similar explanation may work for both solar wind and SNR shocks.     

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.     

Coronal Holes and the High-Speed Solar Wind

Coronal holes are the lowest density plasma components of the Sun's outer atmosphere, and are associated with rapidly expanding magnetic fields and the acceleration of the high-speed solar wind. Spectroscopic and polarimetric observations of the extended corona, coupled with interplanetary particle and radio sounding measurements going back several decades, have put strong constraints on possible explanations for how the plasma in coronal holes receives its extreme kinetic properties. The Ultraviolet Coronagraph Spectrometer (UVCS) aboard the Solar and Heliospheric Observatory (SOHO) spacecraft has revealed surprisingly large temperatures, outflow speeds, and velocity distribution anisotropies for positive ions in coronal holes. We review recent observations, modeling techniques, and proposed heating and acceleration processes for protons, electrons, and heavy ions. We emphasize that an understanding of the acceleration region of the wind (in the nearly collisionless extended corona) is indispensable for building a complete picture of the physics of coronal holes.     

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.     

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.