These notes are based on lectures given by one of us (J.B.T.) at the University of Texas in Austin in 1991. Part I concerns some basic features of plasma confinement by magnetic fields as an introduction to an account of plasma relaxation in Part II. Part III discusses confinement by magnetic mirrors, especially minimum- $B$ systems. It also includes a general discussion of adiabatic invariants and of the principle of maximal ordering in perturbation theory. Part IV is devoted to the analysis of perturbations in toroidal plasmas and the stability of ballooning modes.

The effective potential acting on particles in plasmas being essentially the Debye-shielded Coulomb potential, the particles collisional transport in thermal equilibrium is calculated for all impact parameters b, with a convergent expression reducing to Rutherford scattering for small b. No cutoff at the Debye length scale is needed, and the Coulomb logarithm is only slightly modified.

In honour of the 50th anniversary of the influential review/monograph on plasma turbulence by B. B. Kadomtsev as well as the seminal works of T. H. Dupree and J. Weinstock on resonance-broadening theory, an introductory tutorial is given about some highlights of the statistical–dynamical description of turbulent plasmas and fluids, including the ideas of nonlinear incoherent noise, coherent damping, and self-consistent dielectric response. The statistical closure problem is introduced. Incoherent noise and coherent damping are illustrated with a solvable model of passive advection. Self-consistency introduces turbulent polarization effects that are described by the dielectric function ${\mathcal{D}}$ . Dupree’s method of using ${\mathcal{D}}$ to estimate the saturation level of turbulence is described; then it is explained why a more complete theory that includes nonlinear noise is required. The general theory is best formulated in terms of Dyson equations for the covariance $C$ and an infinitesimal response function $R$ , which subsumes ${\mathcal{D}}$ . An important example is the direct-interaction approximation (DIA). It is shown how to use Novikov’s theorem to develop an $\boldsymbol{x}$ -space approach to the DIA that is complementary to the original $\boldsymbol{k}$ -space approach of Kraichnan. A dielectric function is defined for arbitrary quadratically nonlinear systems, including the Navier–Stokes equation, and an algorithm for determining the form of ${\mathcal{D}}$ in the DIA is sketched. The independent insights of Kadomtsev and Kraichnan about the problem of the DIA with random Galilean invariance are described. The mixing-length formula for drift-wave saturation is discussed in the context of closures that include nonlinear noise (shielded by ${\mathcal{D}}$ ). The role of $R$ in the calculation of the symmetry-breaking (zonostrophic) instability of homogeneous turbulence to the generation of inhomogeneous mean flows is addressed. The second-order cumulant expansion and the stochastic structural stability theory are also discussed in that context. Various historical research threads are mentioned and representative entry points to the literature are given. Some outstanding conceptual issues are enumerated.

Magnetic reconnection and associated heating of ions and electrons in strongly magnetized, weakly collisional plasmas are studied by means of gyrokinetic simulations. It is shown that an appreciable amount of the released magnetic energy is dissipated to yield (irreversible) electron and ion heating via phase mixing. Electron heating is mostly localized to the magnetic island, not the current sheet, and occurs after the dynamical reconnection stage. Ion heating is comparable to electron heating only in high-β plasmas, and results from both parallel and perpendicular phase mixing due to finite Larmor radius (FLR) effects; in space, ion heating is mostly localized to the interior of a secondary island (plasmoid) that arises from the instability of the current sheet.

Many astrophysical disks, such as protoplanetary disks, are in a regime where non-ideal, plasma-specific magnetohydrodynamic (MHD) effects can significantly influence the behaviour of the magnetorotational instability (MRI). The possibility of studying these effects in the plasma Couette experiment (PCX) is discussed. An incompressible, dissipative global stability analysis is developed to include plasma-specific two-fluid effects and neutral collisions, which are inherently absent in analyses of Taylor–Couette flows (TCFs) in liquid metal experiments. It is shown that with boundary driven flows, a ion-neutral collision drag body force significantly affects the azimuthal velocity profile, thus limiting the flows to regime where the MRI is not present. Electrically driven flow (EDF) is proposed as an alternative body force flow drive in which the MRI can destabilize at more easily achievable plasma parameters. Scenarios for reaching MRI relevant parameter space and necessary hardware upgrades are described.

A rotating and magnetized three-dimensional axisymmetric equilibrium for hot plasma confined by a gravitational field is found. The plasma density and current can exhibit strong equatorial plane localization, resulting in disk equilibria with open magnetic field lines. The associated equatorial plane pinching results in magnetic field flaring, implying a strong gravitational squeezing of the plasma carrying ambient magnetic field lines toward the gravitational source. At high plasma pressure, the magnetic field becomes strongly radial outside the disk. The model predicts the rotation frequency bound, the condition for a plasma disk, and the requirement for strong magnetic field flaring.

This work shows the basic foundation of the particle-based representation of low temperature plasma description. In particular, the Monte Carlo Collision (MCC) recipe has been described for the case of electron-atom and ion-atom collisions. The model has been applied to the problem of plasma plume expansion from an electric Hall-effect type thruster. The presence of low energy secondary electrons from electron-atom ionization on the electron energy distribution function (EEDF) have been identified in the first 3 mm from the exit plane where, due to the azimuthal heating the ionization continues to play an important role. In addition, low energy charge-exchange ions from ion-atom electron transfer collisions are evident in the ion energy distribution functions (IEDF) 1 m from the exit plane.

We briefly review some aspects of magnetic turbulence intermittency observed in space plasmas. Deviation of statistical characteristics of a system (e.g. its high statistical momenta) from the Gaussian can manifest itself as domination of rare large intensity peaks often associated with the intermittency in the system's dynamics. Thirty years ago, Zeldovich stressed the importance of the non-Gaussian appearance of the sharp values of vector and scalar physical parameters in random media as a factor of magnetic field amplification in cosmic structures. Magnetic turbulence is governing the behavior of collisionless plasmas in space and especially the physics of shocks and magnetic reconnections. Clear evidence of intermittent magnetic turbulence was found in recent in-situ spacecraft measurements of magnetic fields in the near-Earth and interplanetary plasma environments. We discuss the potentially promising approaches of incorporating the knowledge gained from spacecraft in-situ measurements into modern models describing plasma dynamics and radiation in various astrophysical systems. As an example, we discuss supernova remnants (SNRs) which are known to be the sources of energy, momentum, chemical elements, and high-energy cosmic rays (CRs) in galaxies. Supernova shocks accelerate charged particles to very high energies and may strongly amplify turbulent magnetic fields via instabilities driven by CRs. Relativistic electrons accelerated in SNRs radiate polarized synchrotron emission in a broad range of frequencies spanning from the radio to gamma-rays. We discuss the effects of intermittency of magnetic turbulence on the images of polarized synchrotron X-ray emission of young SNRs and emission spectra of pulsar wind nebula.

JET experiments have compared the efficacy of low- and high-field side ion cyclotron resonance heating (ICRH) as an actuator to deliberately minimise the sawtooth period. It is found that low-field side ICRH with low minority concentration is optimal for sawtooth control for two main reasons. Firstly, low-field side heating means that any toroidal phasing of the ICRH ( $-90^{\circ }$ , $+90^{\circ }$ or dipole) has a destabilising effect on the sawteeth, meaning that dipole phasing can be employed, since this is preferable due to less plasma wall interaction from Resonant Frequency (RF) sheaths. Secondly, the resonance position of the low-field side ICRH does not have to be very accurately placed to achieve sawtooth control, relaxing the requirement for real-time control of the RF frequency. These empirical observations have been confirmed by hybrid kinetic–magnetohydrodynamic modelling, and suggest that the ICRH antenna design for ITER is well positioned to provide a control actuator capable of having a significant effect on the sawtooth behaviour.

Plasmas are ubiquitous in nature, surround our local geospace environment, and permeate the universe. Plasma phenomena in space give rise to energetic particles, the aurora, solar flares and coronal mass ejections, as well as many energetic phenomena in interstellar space. Although plasmas can be studied in laboratory settings, it is often difficult, if not impossible, to replicate the conditions (density, temperature, magnetic and electric fields, etc.) of space. Single-point space missions too numerous to list have described many properties of near-Earth and heliospheric plasmas as measured both in situ and remotely (see http://www.nasa.gov/missions/#.U1mcVmeweRY for a list of NASA-related missions). However, a full description of our plasma environment requires three-dimensional spatial measurements. Cluster is the first, and until data begin flowing from the Magnetospheric Multiscale Mission (MMS), the only mission designed to describe the three-dimensional spatial structure of plasma phenomena in geospace. In this paper, we concentrate on some of the many plasma phenomena that have been studied using data from Cluster. To date, there have been more than 2000 refereed papers published using Cluster data but in this paper we will, of necessity, refer to only a small fraction of the published work. We have focused on a few basic plasma phenomena, but, for example, have not dealt with most of the vast body of work describing dynamical phenomena in Earth's magnetosphere, including the dynamics of current sheets in Earth's magnetotail and the morphology of the dayside high latitude cusp. Several review articles and special publications are available that describe aspects of that research in detail and interested readers are referred to them (see for example, Escoubet et al. 2005Multiscale Coupling of Sun-Earth Processes, p. 459, Keith et al. 2005Sur. Geophys.26, 307–339, Paschmann et al. 2005Outer Magnetospheric Boundaries: Cluster Results, Space Sciences Series of ISSI. Berlin: Springer, Goldstein et al. 2006Adv. Space Res.38, 21–36, Taylor et al. 2010The Cluster Mission: Space Plasma in Three Dimensions, Springer, pp. 309–330 and Escoubet et al. 2013Ann. Geophys.31, 1045–1059).

A simple 1D model is presented for the heating caused by cylindrically-symmetric plasma waves in a non-neutral plasma column due to the addition of a symmetric squeeze potential applied to the center of the column. We study this model by using analytical techniques and by using a numerical grids method solution, and we compare the results of this model to previous work (Ashourvan and Dubin (2014)). squeeze divides the plasma into passing and trapped particles; the latter cannot pass over the squeeze potential. In collisionless theory, enhanced heating is caused by additional bounce harmonics induced by the squeeze in the particle distribution, leading to Landau resonances at energies En for which the bounce frequency ωb(E) and wave frequency ωm satisfy ωm = nωb(En). As a result, heating is substantially higher than the case with no squeeze, even when ωm is much greater than the thermal bounce frequency ωb(T). Adding collisions to the theory creates a boundary layer at the separatrix between trapped and passing particles that further enhances the heating at small ωm/kmvs, where km is the axial wavenumber and vs is the velocity at the separatrix. However, at large ωm/vs, the heating from the separatrix boundary layer is only a small correction to the heating from collisionless resonances in the trapped particle distribution function.

A theoretical framework for low-frequency electromagnetic (drift-)kinetic turbulence in a collisionless, multi-species plasma is presented. The result generalises reduced magnetohydrodynamics (RMHD) and kinetic RMHD (Schekochihin et al., Astrophys. J. Suppl. Ser., vol. 182, 2009, pp. 310–377) to the case where the mean distribution function of the plasma is pressure-anisotropic and different ion species are allowed to drift with respect to each other – a situation routinely encountered in the solar wind and presumably ubiquitous in hot dilute astrophysical plasmas such as the intracluster medium. Two main objectives are achieved. First, in a non-Maxwellian plasma, the relationships between fluctuating fields (e.g. the Alfvén ratio) are order-unity modified compared to the more commonly considered Maxwellian case, and so a quantitative theory is developed to support quantitative measurements now possible in the solar wind. Beyond these order-unity corrections, the main physical feature of low-frequency plasma turbulence survives the generalisation to non-Maxwellian distributions: Alfvénic and compressive fluctuations are energetically decoupled, with the latter passively advected by the former; the Alfvénic cascade is fluid, satisfying RMHD equations (with the Alfvén speed modified by pressure anisotropy and species drifts), whereas the compressive cascade is kinetic and subject to collisionless damping (and for a bi-Maxwellian plasma splits into three independent collisionless cascades). Secondly, the organising principle of this turbulence is elucidated in the form of a conservation law for the appropriately generalised kinetic free energy. It is shown that non-Maxwellian features in the distribution function reduce the rate of phase mixing and the efficacy of magnetic stresses, and that these changes influence the partitioning of free energy amongst the various cascade channels. As the firehose or mirror instability thresholds are approached, the dynamics of the plasma are modified so as to reduce the energetic cost of bending magnetic-field lines or of compressing/rarefying them. Finally, it is shown that this theory can be derived as a long-wavelength limit of non-Maxwellian slab gyrokinetics.

Recent advances in dynamo theory have been made by examining the competition between small- and large-scale dynamos at high magnetic Reynolds number $\mathit{Rm}$ . Small-scale dynamos rely on the presence of chaotic stretching whilst the generation of large-scale fields occurs in flows lacking reflectional symmetry via a systematic electromotive force (EMF). In this paper we discuss how the statistics of the EMF (at high  $\mathit{Rm}$ ) depend on the properties of the multi-scale velocity that is generating it. In particular, we determine that different scales of flow have different contributions to the statistics of the EMF, with smaller scales contributing to the mean without increasing the variance. Moreover, we determine when scales in such a flow act independently in their contribution to the EMF. We further examine the role of large-scale shear in modifying the EMF. We conjecture that the distribution of the EMF, and not simply the mean, largely determines the dominant scale of the magnetic field generated by the flow.

The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility designed to study a range of astrophysically relevant plasma processes as well as novel geometries that mimic astrophysical systems. A multi-cusp magnetic bucket constructed from strong samarium cobalt permanent magnets now confines a $10~\text{m}^{3}$ , fully ionized, magnetic-field-free plasma in a spherical geometry. Plasma parameters of $T_{e}\approx 5$ to $20~\text{eV}$ and $n_{e}\approx 10^{11}$ to $5\times 10^{12}~\text{cm}^{-3}$ provide an ideal testbed for a range of astrophysical experiments, including self-exciting dynamos, collisionless magnetic reconnection, jet stability, stellar winds and more. This article describes the capabilities of WiPAL, along with several experiments, in both operating and planning stages, that illustrate the range of possibilities for future users.

The TORPEX basic plasma physics device at the Center for Plasma Physics Research (CRPP) in Lausanne, Switzerland is described. In TORPEX, simple magnetized toroidal configurations, a paradigm for the tokamak scrape-off layer (SOL), as well as more complex magnetic geometries of direct relevance for fusion are produced. Plasmas of different gases are created and sustained by microwaves in the electron-cyclotron (EC) frequency range. Full diagnostic access allows for a complete characterization of plasma fluctuations and wave fields throughout the entire plasma volume, opening new avenues to validate numerical codes. We detail recent advances in the understanding of basic aspects of plasma turbulence, including its development from linearly unstable electrostatic modes, the formation of filamentary structures, or blobs, and its influence on the transport of energy, plasma bulk and suprathermal ions. We present a methodology for the validation of plasma turbulence codes, which focuses on quantitative assessment of the agreement between numerical simulations and TORPEX experimental data.

We study Landau damping in the 1+1D Vlasov–Poisson system using a Fourier–Hermite spectral representation. We describe the propagation of free energy in Fourier–Hermite phase space using forwards and backwards propagating Hermite modes recently developed for gyrokinetic theory. We derive a free energy equation that relates the change in the electric field to the net Hermite flux out of the zeroth Hermite mode. In linear Landau damping, decay in the electric field corresponds to forward propagating Hermite modes; in nonlinear damping, the initial decay is followed by a growth phase characterized by the generation of backwards propagating Hermite modes by the nonlinear term. The free energy content of the backwards propagating modes increases exponentially until balancing that of the forward propagating modes. Thereafter there is no systematic net Hermite flux, so the electric field cannot decay and the nonlinearity effectively suppresses Landau damping. These simulations are performed using the fully-spectral 5D gyrokinetics code SpectroGK, modified to solve the 1+1D Vlasov–Poisson system. This captures Landau damping via Hou–Li filtering in velocity space. Therefore the code is applicable even in regimes where phase mixing and filamentation are dominant.

Kinetic instabilities in weakly collisional, high beta plasmas are investigated using two-dimensional hybrid expanding box simulations with Coulomb collisions modeled through the Langevin equation (corresponding to the Fokker-Planck one). The expansion drives a parallel or perpendicular temperature anisotropy (depending on the orientation of the ambient magnetic field). For the chosen parameters the Coulomb collisions are important with respect to the driver but are not strong enough to keep the system stable with respect to instabilities driven by the proton temperature anisotropy. In the case of the parallel temperature anisotropy the dominant oblique fire hose instability efficiently reduces the anisotropy in a quasilinear manner. In the case of the perpendicular temperature anisotropy the dominant mirror instability generates coherent compressive structures which scatter protons and reduce the temperature anisotropy. For both the cases the instabilities generate temporarily enough wave energy so that the corresponding (anomalous) transport coefficients dominate over the collisional ones and their properties are similar to those in collisionless plasmas.

A spectacular phenomenon of intermittency, i.e. a progressive growth of higher statistical moments of a physical field excited by an instability in a random medium, attracted the attention of Zeldovich in the last years of his life. At that time, the mathematical aspects underlying the physical description of this phenomenon were still under development and relations between various findings in the field remained obscure. Contemporary results from the theory of the product of independent random matrices (the Furstenberg theory) allowed the elaboration of the phenomenon of intermittency in a systematic way. We consider applications of the Furstenberg theory to some problems in cosmology and dynamo theory.

This issue commemorates an outstanding scientist of the twentieth century, Yakov Borisovich Zeldovich, in connection with the centenary of his birth (8 March 1914), with a collection of reviews and research articles broadly related to large-scale random phenomena in astrophysical plasmas.

A linearised kinetic equation describing electrostatic perturbations of a Maxwellian equilibrium in a weakly collisional plasma forced by a random source is considered. The problem is treated as a kinetic analogue of the Langevin equation and the corresponding fluctuation-dissipation relations are derived. The kinetic fluctuation-dissipation relation reduces to the standard “fluid” one in the regime where the Landau damping rate is small and the system has no real frequency; in this case the simplest possible Landau-fluid closure of the kinetic equation coincides with the standard Langevin equation. Phase mixing of density fluctuations and emergence of fine scales in velocity space is diagnosed as a constant flux of free energy in Hermite space; the fluctuation-dissipation relations for the perturbations of the distribution function are derived, in the form of a universal expression for the Hermite spectrum of the free energy. Finite-collisionality effects are included. This work is aimed at establishing the simplest fluctuation-dissipation relations for a kinetic plasma, clarifying the connection between Landau and Hermite-space formalisms, and setting a benchmark case for a study of phase mixing in turbulent plasmas.