Magnetic fields play an important role in astrophysical accretion discs and in the interstellar and intergalactic medium. They drive jets, suppress fragmentation in star-forming clouds and can have a significant impact on the accretion rate of stars. However, the exact amplification mechanisms of cosmic magnetic fields remain relatively poorly understood. Here, I start by reviewing recent advances in the numerical and theoretical modelling of the turbulent dynamo, which may explain the origin of galactic and intergalactic magnetic fields. While dynamo action was previously investigated in great detail for incompressible plasmas, I here place particular emphasis on highly compressible astrophysical plasmas, which are characterised by strong density fluctuations and shocks, such as the interstellar medium. I find that dynamo action works not only in subsonic plasmas, but also in highly supersonic, compressible plasmas, as well as for low and high magnetic Prandtl numbers. I further present new numerical simulations from which I determine the growth of the turbulent (un-ordered) magnetic field component ( $B_{turb}$ ) in the presence of weak and strong guide fields ( $B_{0}$ ). I vary $B_{0}$ over five orders of magnitude and find that the dependence of $B_{turb}$ on $B_{0}$ is relatively weak, and can be explained with a simple theoretical model in which the turbulence provides the energy to amplify $B_{turb}$ . Finally, I discuss some important implications of magnetic fields for the structure of accretion discs, the launching of jets and the star-formation rate of interstellar clouds.

The existence of a potential energy functional in the zero-Larmor-radius collisionless plasma theory of Kruskal & Oberman (Phys. Fluids, vol. 1, 1958 p. 275), Rosenbluth & Rostoker (Phys. Fluids, vol. 2, 1959, p. 23) allows us to derive easily sufficient conditions for linear stability. However, this kinetic magnetohydrodynamics (KMHD) theory does not have a self-adjointness property, making it difficult to derive necessary conditions. In particular, the standard methods to prove that an instability follows if some trial perturbation makes the incremental potential energy negative, which rely on the self-adjointness of the force operator or on the existence of a complete basis of normal modes, are not applicable to KMHD. This paper investigates KMHD linear stability criteria based on the time evolution of initial-value solutions, without recourse to the classic bounds or comparison theorems of Kruskal–Oberman and Rosenbluth–Rostoker for the KMHD potential energy. The adopted approach does not solve the kinetic equations by integration along characteristics and does not require that the particle orbits be periodic or nearly periodic. Most importantly, the investigation of a necessary condition for stability does not require the self-adjointness of the force operator or the existence of a complete basis of normal modes. It is thereby shown that stability in isothermal ideal-MHD is a sufficient condition for stability in KMHD and that, with a proviso on the long-time behaviour of oscillations about stable equilibria, stability in the double-adiabatic fluid theory, including the variation of the parallel fluid displacement, would be a necessary condition for stability in KMHD.

The self-organisation of turbulence into regular zonal flows can be fruitfully investigated with quasi-linear methods and statistical descriptions. A wave-kinetic equation that assumes asymptotically large-scale zonal flows leads to ultraviolet divergence. From an exact description of quasi-linear dynamics emerges two better geometrical optics approximations. These involve not only the mean flow shear but also the second and third derivative of the mean flow. One approximation takes the form of a new wave-kinetic equation, but is only valid when the zonal flow is quasi-static and wave action is conserved.

Dust acoustic (or dust density) waves have been observed in many laboratory dusty plasmas. These low-frequency waves involve the dynamics of highly charged and massive dust grains, and can be excited by the flow of ions relative to dust. In this paper, we consider the nonlinear development of the dust acoustic instability, excited by thermal ion flow, in a collisional plasma containing dust with high kinetic temperature (warm dust). It is shown that under certain conditions there may be a long-wavelength secondary instability in the nonlinear stage as dust gets heated by the waves. The characteristics of the nonlinear development are considered as a function of the relative charge density of the dust. Application to possible experimental parameters is discussed.

The dynamics of passage and capture into resonance of a distribution of particles driven by a chirped frequency perturbation is discussed. The resonant capture in this case involves crossing of the separatrix by individual particles and, therefore, the adiabatic theorem cannot be used in studying this problem no matter how slow the variation of the driving frequency is. It is shown that, if instead of analysing complicated single particle dynamics in passage through resonance, one considers the slow evolution of a whole distribution of initial conditions in phase space, the adiabaticity and phase space incompressibility arguments yield a solution to the resonant passage problem. This approach is illustrated in the case of an ensemble of electrons driven by a chirped frequency wave passing through Cherenkov resonances with the velocity distribution of electrons.

The stability of a collisionless, magnetised plasma to local convective disturbances is examined, with a focus on kinetic and finite-Larmor-radius effects. Specific application is made to the outskirts of galaxy clusters, which contain hot and tenuous plasma whose temperature increases in the direction of gravity. At long wavelengths (the ‘drift-kinetic’ limit), we obtain the kinetic version of the magnetothermal instability (MTI) and its Alfvénic counterpart (Alfvénic MTI), which were previously discovered and analysed using a magnetofluid (i.e. Braginskii) description. At sub-ion-Larmor scales, we discover an overstability driven by the electron-temperature gradient of kinetic-Alfvén drift waves – the electron MTI (eMTI) – whose growth rate is even larger than the standard MTI. At intermediate scales, we find that ion finite-Larmor-radius effects tend to stabilise the plasma. We discuss the physical interpretation of these instabilities in detail, and compare them both with previous work on magnetised convection in a collisional plasma and with temperature-gradient-driven drift-wave instabilities well known to the magnetic-confinement-fusion community. The implications of having both fluid and kinetic scales simultaneously driven unstable by the same temperature gradient are briefly discussed.

Rotating plasmas have the potential to offer unique capabilities for isotope and mass separation. Among the various electric and magnetic field configurations offering mass separation capabilities, rotating plasmas produced through static or oscillating fields are shown to be a leading candidate for tackling the unsolved problem of large-scale plasma separation. The successful development and deployment of industrial-scale plasma separation technologies could, among many other applications, provide an innovative path towards advanced sustainable nuclear energy. In this context, the potential and versatility of plasma rotation induced by rotating magnetic fields is uncovered and analysed. Analytical stability diagrams are derived from rotating ion orbits as a function of ion mass. Based on these findings, the basic principles of a rotating field plasma separator are then introduced. In light of these results, challenges associated with this original separation process are underlined, and the main directions for a future research program aimed at this important unsolved problem of applied plasma physics are identified.

In recent papers, we show how accretion disks around massive black holes could act as dynamos producing magnetic jets similar to the jets that create spheromaks in the laboratory. In this paper, we discuss how these magnetic astrophysical jets might naturally produce runaway ion beams accelerated to $10^{20}$  eV or more, finally ejected as ultra high energy cosmic rays (UHECRs) long regarded as one of the mysteries of astrophysics. The acceleration is mainly due to the drift cyclotron loss cone kinetic instability known from plasma research. Experiments and simulations are suggested to verify the acceleration process.

The Universe is permeated by hot, turbulent, magnetized plasmas. Turbulent plasma is a major constituent of active galactic nuclei, supernova remnants, the intergalactic and interstellar medium, the solar corona, the solar wind and the Earth’s magnetosphere, just to mention a few examples. Energy dissipation of turbulent fluctuations plays a key role in plasma heating and energization, yet we still do not understand the underlying physical mechanisms involved. THOR is a mission designed to answer the questions of how turbulent plasma is heated and particles accelerated, how the dissipated energy is partitioned and how dissipation operates in different regimes of turbulence. THOR is a single-spacecraft mission with an orbit tuned to maximize data return from regions in near-Earth space – magnetosheath, shock, foreshock and pristine solar wind – featuring different kinds of turbulence. Here we summarize the THOR proposal submitted on 15 January 2015 to the ‘Call for a Medium-size mission opportunity in ESAs Science Programme for a launch in 2025 (M4)’. THOR has been selected by European Space Agency (ESA) for the study phase.

Neutron stars are fascinating astrophysical objects immersed in strong gravitational and electromagnetic fields, at the edge of our current theories. These stars manifest themselves mostly as pulsars, emitting a timely very stable and regular electromagnetic signal. Even though discovered almost fifty years ago, they still remain mysterious compact stellar objects. In this review, we summarize the most fundamental theoretical aspects of neutron star magnetospheres and winds. The main competing models explaining their radiative properties like multi-wavelength pulse shapes and spectra and the underlying physical processes such as pair creation and radiation mechanisms are scrutinized. A global but still rather qualitative picture slowly emerges thanks to recent advances in numerical simulations on the largest scales. However considerations about pulsar magnetospheres remain speculative. For instance, the exact composition of the magnetospheric plasma is not yet known. Is it solely filled with a mixture of $e^{\pm }$  leptons or does it contain a non-negligible fraction of protons and/or ions? Is it almost entirely filled or mostly empty except for some small anecdotal plasma filled regions? Answers to these questions will strongly direct the description of the magnetosphere to seemingly contradictory results leading sometimes to inconsistencies. Nevertheless, accounts are given as to the latest developments in the theory of pulsar magnetospheres and winds, the existence of a possible electrosphere and physical insight obtained from related observational signatures of multi-wavelength pulsed emission.

Two in situ experimental methods are presented in which dust particles are used to determine the extent of the sheath and gain information about the time-averaged electric force profile within a radio frequency (RF) plasma sheath. These methods are advantageous because they are not only simple and quick to carry out, but they also can be performed using standard dusty plasma experimental equipment. In the first method, dust particles are tracked as they fall through the plasma towards the lower electrode. These trajectories are then used to determine the electric force on the particle as a function of height as well as the extent of the sheath. In the second method, dust particle levitation height is measured across a wide range of RF voltages. Similarities were observed between the two experiments, but in order to understand the underlying physics behind these observations, the same conditions were replicated using a self-consistent fluid model. Through comparison of the fluid model and experimental results, it is shown that the particles exhibiting a levitation height that is independent of RF voltage indicate the sheath edge – the boundary between the quasineutral bulk plasma and the sheath. Therefore, both of these simple and inexpensive, yet effective, methods can be applied across a wide range of experimental parameters in any ground-based RF plasma chamber to gain useful information regarding the sheath, which is needed for interpretation of dusty plasma experiments.

This article presents a mini-tutorial aimed at a wide readership not familiar with the field of solar plasma physics. The exposition is centred around the issue of excess/free magnetic energy stored in the solar corona. A general consideration is followed with a particular example of coronal magnetic arcade, where free magnetic energy builds up by photospheric convective flows. In the context of solar physics the major task is to explain how this free energy can be released quickly enough to match what is observed in coronal explosive events such as solar flares. Therefore, in the last section of the paper we discuss briefly a possible role of magnetic reconnection in these processes. This is done in quite simple qualitative physical terms, so that an interested reader can follow it up in more detail with help of the provided references.

Successful phenomenological models of pulsar wind nebulae assume efficient dissipation of the Poynting flux of the magnetized electron–positron wind as well as efficient acceleration of the pairs in the vicinity of the termination shock, but how this is realized is not yet well understood. This paper suggests that the corrugation of the termination shock, at the onset of nonlinearity, may lead towards the desired phenomenology. Nonlinear corrugation of the termination shock would convert a fraction of order unity of the incoming ordered magnetic field into downstream turbulence, slowing down the flow to sub-relativistic velocities. The dissipation of turbulence would further preheat the pair population on short length scales, close to equipartition with the magnetic field, thereby reducing the initial high magnetization to values of order unity. Furthermore, it is speculated that the turbulence generated by the corrugation pattern may sustain a relativistic Fermi process, accelerating particles close to the radiation reaction limit, as observed in the Crab nebula. The required corrugation could be induced by the fast magnetosonic modes of downstream nebular turbulence; but it could also be produced by upstream turbulence, either carried by the wind or seeded in the precursor by the accelerated particles themselves.

Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been carried out using the hybrid method, wherein the ions are treated as particles and the electron response is modelled as a massless fluid. We confirm the results of earlier two-dimensional (2-D) simulations showing both field-parallel aligned propagating fluctuations and fluctuations carried by the reflected-gyrating ions. In addition, it is shown that, for 2-D simulations of the shock coplanarity plane, the presence of short-wavelength fluctuations in all magnetic components is associated with the ion Weibel instability driven at the upstream edge of the foot by the reflected-gyrating ions. In 3-D simulations we show for the first time that the dominant microstructure is due to a coupling between field-parallel propagating fluctuations in the ramp and the motion of the reflected ions. This results in a pattern of fluctuations counter-propagating across the surface of the shock at an angle inclined to the magnetic field direction, due to a combination of field-parallel motion at the Alfvén speed of the ramp and motion in the sense of gyration of the reflected ions.

Enhancement of the quality of laser wake-field accelerated (LWFA) electron beams implies the improvement and controllability of the properties of the wake waves generated by ultra-short pulse lasers in underdense plasmas. In this work we present a compendium of useful formulas giving relations between the laser and plasma target parameters allowing one to obtain basic dependences, e.g. the energy scaling of the electrons accelerated by the wake field excited in inhomogeneous media including multi-stage LWFA accelerators. Consideration of the effects of using the chirped laser pulse driver allows us to find the regimes where the chirp enhances the wake field amplitude. We present an analysis of the three-dimensional effects on the electron beam loading and on the unlimited LWFA acceleration in inhomogeneous plasmas. Using the conditions of electron trapping to the wake-field acceleration phase we analyse the multi-equal stage and multiuneven stage LWFA configurations. In the first configuration the energy of fast electrons is a linear function of the number of stages, and in the second case, the accelerated electron energy grows exponentially with the number of stages. The results of the two-dimensional particle-in-cell simulations presented here show the high quality electron acceleration in the triple stage injection–acceleration configuration.

The most massive baryonic component of galaxy clusters is the ‘intracluster medium’ (ICM), a diffuse, hot, weakly magnetized plasma that is most easily observed in the X-ray band. Despite being observed for decades, the macroscopic transport properties of the ICM are still not well constrained. A path to determine macroscopic ICM properties opened up with the discovery of ‘cold fronts’. These were observed as sharp discontinuities in surface brightness and temperature in the ICM, with the property that the denser side of the discontinuity is the colder one. The high spatial resolution of the Chandra X-ray Observatory revealed two puzzles about cold fronts. First, they should be subject to Kelvin–Helmholtz instabilities, yet in many cases they appear relatively smooth and undisturbed. Second, the width of the interface between the two gas phases is typically narrower than the mean free path of the particles in the plasma, indicating negligible thermal conduction. It was thus realized that these special characteristics of cold fronts may be used to probe the properties of the cluster plasma. In this review, we will discuss the recent simulations of cold fronts in galaxy clusters, focusing on those which have attempted to use these features to constrain ICM physics. In particular, we will examine the effects of magnetic fields, viscosity, and thermal conductivity on the stability properties and long-term evolution of cold fronts. We conclude with a discussion on what important questions remain unanswered, and the future role of simulations and the next generation of X-ray observatories.

The large number of ${\it\gamma}$ -ray pulsars discovered by the Fermi Gamma-Ray Space Telescope since its launch in 2008 dwarfs the handful that were previously known. The variety of observed light curves makes possible a tomography of both the ensemble-averaged field structure and the high-energy emission regions of a pulsar magnetosphere. Fitting the ${\it\gamma}$ -ray pulsar light curves with model magnetospheres and emission models has revealed that most of the high-energy emission, and the particles acceleration, takes place near or beyond the light cylinder, near the current sheet. As pulsar magnetosphere models become more sophisticated, it is possible to probe magnetic field structure and emission that are self-consistently determined. Light curve modelling will continue to be a powerful tool for constraining the pulsar magnetosphere physics.

These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is ‘frozen in’ to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as their stability and their response to tidal forcing.

Particle image velocimetry is a fluid measurement technique that has been used for more than 20 years to characterize the particle transport and thermal state of dusty plasma systems. This manuscript provides an overview of this diagnostic technique, highlighting the strengths and limitations that are associated with its use. Additionally, the variations of this technique that have been applied in the study of dusty plasma systems will be discussed, along with a small selection of measurements that can be made with the technique. Potential future directions for this diagnostic tool within the dusty plasma community will also be discussed.

Counter-propagating and suitably polarized light (laser) beams can provide conditions for pair production. Here, we consider in more detail the following two situations: (i) in the homogeneity regions of anti-nodes of linearly polarized ultra-high intensity laser beams, the Schwinger process is dynamically assisted by a second high-frequency field, e.g. by an XFEL beam; and (ii) a high-energy probe photon beam colliding with a superposition of co-propagating intense laser and XFEL beams gives rise to the laser-assisted Breit–Wheeler process. The prospects of such bi-frequent field constellations with respect to the feasibility of conversion of light into matter are discussed.