The threshold for the resonant destabilisation of ion-temperature-gradient (ITG) driven instabilities that render the modes ubiquitous in both tokamaks and stellarators is investigated. We discover remarkably similar results for both confinement concepts if care is taken in the analysis of the effect of the global shear ${\hat{s}}$. We revisit, analytically and by means of gyrokinetic simulations, accepted tokamak results and discover inadequacies of some aspects of their theoretical interpretation. In particular, for standard tokamak configurations, we find that global shear effects on the critical gradient cannot be attributed to the wave–particle resonance destabilising mechanism of Hahm & Tang (Phys. Plasmas, vol. 1, 1989, pp. 1185–1192), but are consistent with a stabilising contribution predicted by Biglari et al. (Phys. Plasmas, vol. 1, 1989, pp. 109–118). Extensive analytical and numerical investigations show that virtually no previous tokamak theoretical predictions capture the temperature dependence of the mode frequency at marginality, thus leading to incorrect instability thresholds. In the asymptotic limit ${\hat{s}}\unicode[STIX]{x1D704}\ll 1$, where $\unicode[STIX]{x1D704}$ is the rotational transform, and such a threshold should be solely determined by the resonant toroidal branch of the ITG mode, we discover a family of unstable solutions below the previously known threshold of instability. This is true for a tokamak case described by a local ${\hat{s}}-\unicode[STIX]{x1D6FC}$ local equilibrium, and for the stellarator Wendelstein 7-X, where these unstable solutions are present even for configurations with a small trapped-particle population. We conjecture they are of the Floquet type and derive their properties from the Fourier analysis of toroidal drift modes of Connor & Taylor (Phys. Fluids, vol. 30, 1987, pp. 3180–3185), and to Hill’s theory of the motion of the lunar perigee (Acta Math., vol. 8, 1886, pp. 1–36). The temperature dependence of the newly determined threshold is given for both confinement concepts. In the first case, the new temperature-gradient threshold is found to be rather insensitive to the temperature ratio $T_{i}/T_{e}$, at least for $T_{i}/T_{e}\lesssim 1$, and to be a growing function of the density gradient scale for $T_{i}/T_{e}\gtrsim 1$. For Wendelstein 7-X, the new critical temperature gradient is a growing function of the temperature ratio. The importance of these findings for the assessment of turbulence in stellarators and low-shear tokamak configurations is discussed.

Evolving magnetic fields are shown to generically reach a state of fast magnetic reconnection in which magnetic field line connections change and magnetic energy is released at an Alfvénic rate. This occurs even in plasmas with zero resistivity; only the finiteness of the mass of the lightest charged particle, an electron, is required. The speed and prevalence of Alfvénic or fast magnetic reconnection imply that its cause must be contained within the ideal evolution equation for magnetic fields, $\unicode[STIX]{x2202}\boldsymbol{B}/\unicode[STIX]{x2202}t=\unicode[STIX]{x1D735}\times (\boldsymbol{u}\times \boldsymbol{B})$, where $\boldsymbol{u}(\boldsymbol{x},t)$ is the velocity of the magnetic field lines. For a generic $\boldsymbol{u}(\boldsymbol{x},t)$, neighbouring magnetic field lines develop a separation that increases exponentially, as $e^{\unicode[STIX]{x1D70E}(\ell ,t)}$ with $\ell$ the distance along a line. This exponentially enhances the sensitivity of the evolution to non-ideal effects. An analogous effect, the importance of stirring to produce a large-scale flow and enhance mixing, has been recognized by cooks through many millennia, but the importance of the large-scale flow $\boldsymbol{u}$ to reconnection is customarily ignored. In part this is due to the sixty-year focus of recognition theory on two-coordinate models, which eliminate the exponential enhancement that is generic with three coordinates. A simple three-coordinate model is developed, which could be used to address many unanswered questions.

Reconnection outflows have been under intense recent scrutiny, from in situ observations and from simulations. These regions are host to a variety of instabilities and intense energy exchanges, often even superior to the main reconnection site. We report here a number of results drawn from an investigation of simulations. First, the outflows are observed to become unstable to drift instabilities. Second, these instabilities lead to the formation of secondary reconnection sites. Third, the secondary processes are responsible for large energy exchanges and particle energization. Finally, the particle distribution function are modified to become non-Maxwellian and include multiple interpenetrating populations.

It is demonstrated that the azimuthal magnetorotational instability (AMRI) also works with radially increasing rotation rates contrary to the standard magnetorotational instability for axial fields which requires negative shear. The stability against non-axisymmetric perturbations of a conducting Taylor–Couette flow with positive shear under the influence of a toroidal magnetic field is considered if the background field between the cylinders is current free. For small magnetic Prandtl number $Pm\rightarrow 0$ the curves of neutral stability converge in the (Hartmann number,Reynolds number) plane approximating the stability curve obtained in the inductionless limit $Pm=0$. The numerical solutions for $Pm=0$ indicate the existence of a lower limit of the shear rate. For large $Pm$ the curves scale with the magnetic Reynolds number of the outer cylinder but the flow is always stable for magnetic Prandtl number unity as is typical for double-diffusive instabilities. We are particularly interested to know the minimum Hartmann number for neutral stability. For models with resting or almost resting inner cylinder and with perfectly conducting cylinder material the minimum Hartmann number occurs for a radius ratio of $r_{\text{in}}=0.9$. The corresponding critical Reynolds numbers are smaller than $10^{4}$.

We study standing waves on the surface of a tangential discontinuity in an incompressible plasma. The plasma is moving with constant velocity at one side of the discontinuity, while it is at rest at the other side. The moving plasma is ideal and the plasma at rest is viscous. We only consider the long wavelength limit where the viscous Reynolds number is large. A standing wave is a superposition of a forward and a backward wave. When the flow speed is between the critical speed and the Kelvin–Helmholtz threshold the backward wave is a negative energy wave, while the forward wave is always a positive energy wave. We show that viscosity causes the standing wave to grow. Its increment is equal to the difference between the negative energy wave increment and the positive energy wave decrement.

Large-angle Coulomb collisions lead to an avalanching generation of runaway electrons in a plasma. We present the first fully conservative large-angle collision operator, derived from the relativistic Boltzmann operator. The relation to previous models for large-angle collisions is investigated, and their validity assessed. We present a form of the generalized collision operator which is suitable for implementation in a numerical kinetic equation solver, and demonstrate the effect on the runaway-electron growth rate. Finally we consider the reverse avalanche effect, where runaways are slowed down by large-angle collisions, and show that the choice of operator is important if the electric field is close to the avalanche threshold.

In space and astrophysical plasmas, turbulence is responsible for transferring energy from large scales driven by violent events or instabilities, to smaller scales where turbulent energy is ultimately converted into plasma heat by dissipative mechanisms. The nonlinear interaction between counterpropagating Alfvén waves, denoted Alfvén wave collisions, drives this turbulent energy cascade, as recognized by early work with incompressible magnetohydrodynamic (MHD) equations. Recent work employing analytical calculations and nonlinear gyrokinetic simulations of Alfvén wave collisions in an idealized periodic initial state have demonstrated the key properties that strong Alfvén wave collisions mediate effectively the transfer of energy to smaller perpendicular scales and self-consistently generate current sheets. For the more realistic case of the collision between two initially separated Alfvén wavepackets, we use a nonlinear gyrokinetic simulation to show here that these key properties persist: strong Alfvén wavepacket collisions indeed facilitate the perpendicular cascade of energy and give rise to current sheets. Furthermore, the evolution shows that nonlinear interactions occur only while the wavepackets overlap, followed by a clean separation of the wavepackets with straight uniform magnetic fields and the cessation of nonlinear evolution in between collisions, even in the gyrokinetic simulation presented here which resolves dispersive and kinetic effects beyond the reach of the MHD theory.

Wave–particle interaction is a key process in particle diffusion in collisionless plasmas. We look into the interaction of single plasma waves with individual particles and discuss under which circumstances this is a chaotic process, leading to diffusion. We derive the equations of motion for a particle in the fields of a magnetostatic, circularly polarized, monochromatic wave and show that no chaotic particle motion can arise under such circumstances. A novel and exact analytic solution for the equations is presented. Additional plasma waves lead to a breakdown of the analytic solution and chaotic particle trajectories become possible. We demonstrate this effect by considering a linearly polarized, monochromatic wave, which can be seen as the superposition of two circularly polarized waves. Test particle simulations are provided to illustrate and expand our analytical considerations.

Understanding the removal of energy from turbulent fluctuations in a magnetized plasma and the consequent energization of the constituent plasma particles is a major goal of heliophysics and astrophysics. Previous work has shown that nonlinear interactions among counterpropagating Alfvén waves – or Alfvén wave collisions – are the fundamental building block of astrophysical plasma turbulence and naturally generate current sheets in the strongly nonlinear limit. A nonlinear gyrokinetic simulation of a strong Alfvén wave collision is used to examine the damping of the electromagnetic fluctuations and the associated energization of particles that occurs in self-consistently generated current sheets. A simple model explains the flow of energy due to the collisionless damping and the associated particle energization, as well as the subsequent thermalization of the particle energy by collisions. The net particle energization by the parallel electric field is shown to be spatially localized, and the nonlinear evolution is essential in enabling spatial non-uniformity. Using the recently developed field–particle correlation technique, we show that particles resonant with the Alfvén waves in the simulation dominate the energy transfer, demonstrating conclusively that Landau damping plays a key role in the spatially localized damping of the electromagnetic fluctuations and consequent energization of the particles in this strongly nonlinear simulation.

In this paper, weak-turbulence theory is used to investigate the nonlinear evolution of the parametric instability in three-dimensional low-$\unicode[STIX]{x1D6FD}$ plasmas at wavelengths much greater than the ion inertial length under the assumption that slow magnetosonic waves are strongly damped. It is shown analytically that the parametric instability leads to an inverse cascade of Alfvén wave quanta, and several exact solutions to the wave kinetic equations are presented. The main results of the paper concern the parametric decay of Alfvén waves that initially satisfy $e^{+}\gg e^{-}$, where $e^{+}$ and $e^{-}$ are the frequency ($f$) spectra of Alfvén waves propagating in opposite directions along the magnetic field lines. If $e^{+}$ initially has a peak frequency $f_{0}$ (at which $fe^{+}$ is maximized) and an ‘infrared’ scaling $f^{p}$ at smaller $f$ with $-1<p<1$, then $e^{+}$ acquires an $f^{-1}$ scaling throughout a range of frequencies that spreads out in both directions from $f_{0}$. At the same time, $e^{-}$ acquires an $f^{-2}$ scaling within this same frequency range. If the plasma parameters and infrared $e^{+}$ spectrum are chosen to match conditions in the fast solar wind at a heliocentric distance of 0.3 astronomical units (AU), then the nonlinear evolution of the parametric instability leads to an $e^{+}$ spectrum that matches fast-wind measurements from the Helios spacecraft at 0.3 AU, including the observed $f^{-1}$ scaling at $f\gtrsim 3\times 10^{-4}~\text{Hz}$. The results of this paper suggest that the $f^{-1}$ spectrum seen by Helios in the fast solar wind at $f\gtrsim 3\times 10^{-4}~\text{Hz}$ is produced in situ by parametric decay and that the $f^{-1}$ range of $e^{+}$ extends over an increasingly narrow range of frequencies as $r$ decreases below 0.3 AU. This prediction will be tested by measurements from the Parker Solar Probe.

A class of simple kinetic systems is considered, described by the one-dimensional Vlasov–Landau equation with Poisson or Boltzmann electrostatic response and an energy source. Assuming a stochastic electric field, a solvable model is constructed for the phase-space turbulence of the particle distribution. The model is a kinetic analogue of the Kraichnan–Batchelor model of chaotic advection. The solution of the model is found in Fourier–Hermite space and shows that the free-energy flux from low to high Hermite moments is suppressed, with phase mixing cancelled on average by anti-phase-mixing (stochastic plasma echo). This implies that Landau damping is an ineffective route to dissipation (i.e. to thermalisation of electric energy via velocity space). The full Fourier–Hermite spectrum is derived. Its asymptotics are $m^{-3/2}$ at low wavenumbers and high Hermite moments ($m$) and $m^{-1/2}k^{-2}$ at low Hermite moments and high wavenumbers ($k$). These conclusions hold at wavenumbers below a certain cutoff (analogue of Kolmogorov scale), which increases with the amplitude of the stochastic electric field and scales as inverse square of the collision rate. The energy distribution and flows in phase space are a simple and, therefore, useful example of competition between phase mixing and nonlinear dynamics in kinetic turbulence, reminiscent of more realistic but more complicated multi-dimensional systems that have not so far been amenable to complete analytical solution.

First-principles simulations of tokamak turbulence have proven to be of great value in recent decades. We develop a pseudo-spectral velocity formulation of the turbulence equations that smoothly interpolates between the highly efficient but lower resolution three-dimensional (3-D) gyrofluid representation and the conventional but more expensive 5-D gyrokinetic representation. Our formulation is a projection of the nonlinear gyrokinetic equation onto a Laguerre–Hermite velocity-space basis. We discuss issues related to collisions, closures and entropy. While any collision operator can be used in the formulation, we highlight a model operator that has a particularly sparse Laguerre–Hermite representation, while satisfying conservation laws and the H theorem. Free streaming, magnetic drifts and nonlinear phase mixing each give rise to closure problems, which we discuss in relation to the instabilities of interest and to free energy conservation. We show that the model is capable of reproducing gyrokinetic results for linear instabilities and zonal flow dynamics. Thus the final model is appropriate for the study of instabilities, turbulence and transport in a wide range of geometries, including tokamaks and stellarators.

In space and astrophysical plasmas, violent events or instabilities inject energy into turbulent motions at large scales. Nonlinear interactions among the turbulent fluctuations drive a cascade of energy to small perpendicular scales at which the energy is ultimately converted into plasma heat. Previous work with the incompressible magnetohydrodynamic (MHD) equations has shown that this turbulent energy cascade is driven by the nonlinear interaction between counterpropagating Alfvén waves – also known as Alfvén wave collisions. Direct numerical simulations of weakly collisional plasma turbulence enables deeper insight into the nature of the nonlinear interactions underlying the turbulent cascade of energy. In this paper, we directly compare four cases: both periodic and localized Alfvén wave collisions in the weakly and strongly nonlinear limits. Our results reveal that in the more realistic case of localized Alfvén wave collisions (rather than the periodic case), all nonlinearly generated fluctuations are Alfvén waves, which mediates nonlinear energy transfer to smaller perpendicular scales.

Decaying homogeneous and isotropic magnetohydrodynamics (MHD) turbulence is investigated numerically at large Reynolds numbers thanks to the eddy-damped quasi-normal Markovian (EDQNM) approximation. Without any background mean magnetic field, the total energy spectrum $E$ scales as $k^{-3/2}$ in the inertial range as a consequence of the modelling. Moreover, the total energy is shown, both analytically and numerically, to decay at the same rate as kinetic energy in hydrodynamic isotropic turbulence: this differs from a previous prediction, and thus physical arguments are proposed to reconcile both results. Afterwards, the MHD turbulence is made imbalanced by an initial non-zero cross-helicity. A spectral modelling is developed for the velocity–magnetic correlation in a general homogeneous framework, which reveals that cross-helicity can contain subtle anisotropic effects. In the inertial range, as the Reynolds number increases, the slope of the cross-helical spectrum becomes closer to $k^{-5/3}$ than $k^{-2}$ . Furthermore, the Elsässer spectra deviate from $k^{-3/2}$ with cross-helicity at large Reynolds numbers. Regarding the pressure spectrum $E_{P}$ , its kinetic and magnetic parts are found to scale with $k^{-2}$ in the inertial range, whereas the part due to cross-helicity rather scales in $k^{-7/3}$ . Finally, the two $4/3$ rd laws for the total energy and cross-helicity are assessed numerically at large Reynolds numbers.

A drift-kinetic calculation in an axisymmetric tokamak, with super-diamagnetic flows, is presented to elucidate the relation between the radial electric field, $E_{r}$, zonal flows and the rapid precession of the trapped particle (TP) population. It has been shown earlier (Rosenbluth & Hinton, Phys. Rev. Lett., vol. 80(4), 1998, p. 724, hereafter RH) that an initial radial electric field results in geodesic acoustic mode oscillations which subsequently Landau damp, resulting in a much smaller final residual electric field, and accompanying parallel zonal flows. We observe an apparent paradox: the final angular momentum in the RH parallel zonal flow is much smaller than the angular momentum expected from the well-known rapid precession of the trapped particle population in the RH residual electric field. We reconcile this paradox by illuminating the presence of a population of reverse circulating particle flows that, dominantly, are equal and opposite to the rapid TP precession. Mathematically, the calculation is facilitated by transforming to an energy coordinate shifted from conventional by an amount proportional to $E_{r}$. We also discuss the well-known RH coefficient in the context of effective mass and show how the TP precession and the opposite circulating flows contribute to this mass. Finally, we show that in the long wavelength limit, the RH flows and RH coefficient arise as a natural consequence of conservation of toroidal angular momentum and the second adiabatic invariant.

We compute the real and imaginary parts of the electric permittivities and magnetic permeabilities of relativistic electrons from quantum electrodynamics at finite temperatures and densities, for weak fields, neglecting electron–electron interactions. For non-zero temperatures, electromagnetic responses are reduced to one-dimensional integrals computed numerically. For zero temperature, we find analytic expressions for both their real/dispersive and imaginary/absorptive parts. As an application of our results, we obtain the dispersion relation for longitudinal electric plasmons. Present calculations support our recent claim that, at low frequencies and long wavelengths, the system will exhibit simultaneously negative electric and magnetic responses.

The structure of a whistler precursor in a quasi-perpendicular shock is studied within two-fluid approach in one-dimensional case. The complete set of equations is reduced to the KdV equation, if no dissipation is included. With a phenomenological resistive dissipation the structure is described with the KdV–Burgers equation. The shock profile is intrinsically time dependent. For sufficiently strong dissipation, temporal evolution of a steepening profile results in generation of a stationary decaying whistler ahead of the shock front. With the decrease of the dissipation parameter, whistler wave trains begin to detach and propagate toward the upstream and the ramp is weakly time dependent. In the weakly dissipative regime the shock front experiences reformation.

We have explored the thermodynamics of compressed magnetized plasmas in laboratory experiments and we call these studies ‘magnetothermodynamics’. The experiments are carried out in the Swarthmore Spheromak eXperiment device. In this device, a magnetized plasma source is located at one end and at the other end, a closed conducting can is installed. We generate parcels of magnetized plasma and observe their compression against the end wall of the conducting cylinder. The plasma parameters such as plasma density, temperature and magnetic field are measured during compression using HeNe laser interferometry, ion Doppler spectroscopy and a linear ${\dot{B}}$ probe array, respectively. To identify the instances of ion heating during compression, a PV diagram is constructed using measured density, temperature and a proxy for the volume of the magnetized plasma. Different equations of state are analysed to evaluate the adiabatic nature of the compressed plasma. A three-dimensional resistive magnetohydrodynamic code (NIMROD) is employed to simulate the twisted Taylor states and shows stagnation against the end wall of the closed conducting can. The simulation results are consistent to what we observe in our experiments.

Sweet–Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. Phys. Fluids, vol. 6 (4), 1963, pp. 459–484, Rutherford, Phys. Fluids, vol. 16 (11), 1973, pp. 1903–1908, Coppi et al. Fizika Plazmy, vol. 2, 1976, pp. 961–966) is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal magnetohydrodynamics, a coordinate transformation reveals that the flow compresses and stretches perturbations. A multiple scale formulation describes how linear tearing mode theory (Furth et al. Phys. Fluids, vol. 6 (4), 1963, pp. 459–484, Coppi et al. Fizika Plazmy, vol. 2, 1976, pp. 961–966) can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

The gyrokinetic stability of electron–positron plasmas contaminated by an ion (proton) admixture is studied in a slab geometry. The appropriate dispersion relation is derived and solved. Stable K-modes, the universal instability, the ion-temperature-gradient-driven instability, the electron-temperature-gradient-driven instability and the shear Alfvén wave are considered. It is found that the contaminated plasma remains stable if the contamination degree is below some threshold and that the shear Alfvén wave can be present in a contaminated plasma in cases where it is absent without ion contamination.