We present a complete analysis of all wave modes in a cold pair plasma, significantly extending standard textbook treatments. Instead of identifying the maximal number of two propagating waves at fixed frequency $\unicode[STIX]{x1D714}$, we introduce a unique labelling of all 5 mode pairs described by the general dispersion relation $\unicode[STIX]{x1D714}(k)$, starting from their natural ordering at small wavenumber $k$. There, the 5 pairs start off as Alfvén (A), fast magnetosonic (F), modified electrostatic (M) and electromagnetic O and X branches, and each $\unicode[STIX]{x1D714}(k)$ branch smoothly connects to large wavenumber resonances or limits. For cold pair plasmas, these 5 branches show avoided crossings, which become true crossings at exactly parallel or perpendicular orientation. Only for those orientations, we find a changed connectivity between small and large wavenumber behaviour. Analysing phase and group diagrams for all 5 wave modes, distinctly different from the Clemmow–Mullaly–Allis representation, reveals the true anisotropy of the A, M and O branches.

In a magnetized, collisionless plasma, the magnetic moment of the constituent particles is an adiabatic invariant. An increase in the magnetic-field strength in such a plasma thus leads to an increase in the thermal pressure perpendicular to the field lines. Above a $\unicode[STIX]{x1D6FD}$-dependent threshold (where $\unicode[STIX]{x1D6FD}$ is the ratio of thermal to magnetic pressure), this pressure anisotropy drives the mirror instability, producing strong distortions in the field lines on ion-Larmor scales. The impact of this instability on magnetic reconnection is investigated using a simple analytical model for the formation of a current sheet (CS) and the associated production of pressure anisotropy. The difficulty in maintaining an isotropic, Maxwellian particle distribution during the formation and subsequent thinning of a CS in a collisionless plasma, coupled with the low threshold for the mirror instability in a high-$\unicode[STIX]{x1D6FD}$ plasma, imply that the geometry of reconnecting magnetic fields can differ radically from the standard Harris-sheet profile often used in simulations of collisionless reconnection. As a result, depending on the rate of CS formation and the initial CS thickness, tearing modes whose growth rates and wavenumbers are boosted by this difference may disrupt the mirror-infested CS before standard tearing modes can develop. A quantitative theory is developed to illustrate this process, which may find application in the tearing-mediated disruption of kinetic magnetorotational ‘channel’ modes.

The formation of zonal flows from inhomogeneous drift-wave (DW) turbulence is often described using statistical theories derived within the quasilinear approximation. However, this approximation neglects wave–wave collisions. Hence, some important effects such as the Batchelor–Kraichnan inverse-energy cascade are not captured within this approach. Here we derive a wave kinetic equation that includes a DW collision operator in the presence of zonal flows. Our derivation makes use of the Weyl calculus, the quasinormal statistical closure and the geometrical-optics approximation. The obtained model conserves both the total enstrophy and energy of the system. The derived DW collision operator breaks down at the Rayleigh–Kuo threshold. This threshold is missed by homogeneous-turbulence theory but expected from a full-wave quasilinear analysis. In the future, this theory might help better understand the interactions between drift waves and zonal flows, including the validity domain of the quasilinear approximation that is commonly used in the literature.

Energy chirp compensation of the electron bunch (e-bunch) in a laser wakefield accelerator, which is caused by the phase space rotation in the gradient wakefield, has been applied in many schemes for low energy spread e-bunch generation. We report the experimental observation of energy chirp compensation of the e-bunch in a nonlinear laser wakefield accelerator with a negligible beam loading effect. By adjusting the acceleration length using a wedge-roof block, the chirp compensation of the accelerated e-bunch was observed via an electron spectrometer. Apart from this, some significant parameters for the compensation process, such as the longitudinal dispersion and wakefield slope at the bunch position, were also estimated. A detailed comparison between experiment and simulation shows good agreement of the wakefield and bunch parameters. These results give a clear demonstration of the longitudinal characteristics of the wakefield in a plasma and the bunch dynamics, which are important for better control of a compact laser wakefield accelerator.

Quasisymmetric stellarators are appealing intellectually and as fusion reactor candidates since the guiding-centre particle trajectories and neoclassical transport are isomorphic to those in a tokamak, implying good confinement. Previously, quasisymmetric magnetic fields have been identified by applying black-box optimization algorithms to minimize symmetry-breaking Fourier modes of the field strength $B$. Here, instead, we directly construct magnetic fields in cylindrical coordinates that are quasisymmetric to leading order in the distance from the magnetic axis, without using optimization. The method involves solution of a one-dimensional nonlinear ordinary differential equation, originally derived by Garren & Boozer (Phys. Fluids B, vol. 3, 1991, p. 2805). We demonstrate the usefulness and accuracy of this optimization-free approach by providing the results of this construction as input to the codes VMEC and BOOZ_XFORM, confirming the purity and scaling of the magnetic spectrum. The space of magnetic fields that are quasisymmetric to this order is parameterized by the magnetic axis shape along with three other real numbers, one of which reflects the on-axis toroidal current density, and another one of which is zero for stellarator symmetry. The method here could be used to generate good initial conditions for conventional optimization, and its speed enables exhaustive searches of parameter space.

The concept of informativeness of nonlinear plasma physics scenarios is explained. Natural ideas of developing highly informative models of plasma kinetics are spelled out. They are applied to develop a formula that governs the drift of long Langmuir waves in spatial positions and wave vectors in a magnetized plasma due to the plasma inhomogeneity. Together with previous findings (Erofeev, Phys. Plasmas, vol. 22, 2015, 092302), the formula evidences the need for an intelligent generalization of the notion of wave energy density from usual homogeneous plasmas to inhomogeneous ones.

The injection of cryogenic pellets into a magnetically confined plasma is shown to be accompanied by a considerable transfer of thermal energy from the electrons in the background plasma to the ions. The resulting ion heating can be significant, particularly in plasmas with disparate electron and ion temperatures, and can affect the energy balance of the plasma. In recent Wendelstein 7-X experiments, this mechanism can account for a substantial fraction of the ion heating power during pellet injection.

Acoustic solitons obtained through a reductive perturbation scheme are normally governed by a Korteweg–de Vries (KdV) equation. In multispecies plasmas at critical compositions the coefficient of the quadratic nonlinearity vanishes. Extending the analytic treatment then leads to a modified KdV (mKdV) equation, which is characterized by a cubic nonlinearity and is even in the electrostatic potential. The mKdV equation admits solitons having opposite electrostatic polarities, in contrast to KdV solitons which can only be of one polarity at a time. A Hirota formalism has been used to derive the two-soliton solution. That solution covers not only the interaction of same-polarity solitons but also the collision of compressive and rarefactive solitons. For the visualization of the solutions, the focus is on the details of the interaction region. A novel and detailed discussion is included of typical electric field signatures that are often observed in ionospheric and magnetospheric plasmas. It is argued that these signatures can be attributed to solitons and their interactions. As such, they have received little attention.

We study the resistive evolution of a localized self-organizing magnetohydrodynamic equilibrium. In this configuration the magnetic forces are balanced by a pressure force caused by a toroidal depression in the pressure. Equilibrium is attained when this low-pressure region prevents further expansion into the higher-pressure external plasma. We find that, for the parameters investigated, the resistive evolution of the structures follows a universal pattern when rescaled to resistive time. The finite resistivity causes both a decrease in the magnetic field strength and a finite slip of the plasma fluid against the static equilibrium. This slip is caused by a Pfirsch–Schlüter-type diffusion, similar to what is seen in tokamak equilibria. The net effect is that the configuration remains in magnetostatic equilibrium whilst it slowly grows in size. The rotational transform of the structure becomes nearly constant throughout the entire structure, and decreases according to a power law. In simulations this equilibrium is observed when highly tangled field lines relax in a high-pressure (relative to the magnetic field strength) environment, a situation that occurs when the twisted field of a coronal loop is ejected into the interplanetary solar wind. In this paper we relate this localized magnetohydrodynamic equilibrium to magnetic clouds in the solar wind.

In strictly collisionless electrostatic shocks, the ion distribution function can develop discontinuities along phase-space separatrices, due to partial reflection of the ion population. In this paper, we depart from the strictly collisionless regime and present a semi-analytical model for weakly collisional kinetic shocks. The model is used to study the effect of small but finite collisionalities on electrostatic shocks, and they are found to smooth out these discontinuities into growing boundary layers. More importantly, ions diffuse into and accumulate in the previously empty regions of phase space, and, by upsetting the charge balance, lead to growing downstream oscillations of the electrostatic potential. We find that the collisional age of the shock is the more relevant measure of the collisional effects than the collisionality, where the former can become significant during the lifetime of the shock, even for weak collisionalities.

We have combined global magnetohydrodynamic (MHD) simulations of the solar wind and magnetosphere interaction with an implicit particle-in-cell simulation (PIC) and used this approach to model magnetic reconnection at both the dayside magnetopause and in the magnetotail plasma sheet. In this approach, we first model the magnetospheric configuration driven by the solar wind using the MHD simulation. At a time of interest (usually when a thin current sheet has formed in the MHD simulation), we load a large particle-in-cell simulation with plasma and fields based on the MHD state. We use the MHD results to set the boundary conditions on the PIC simulation. The coupling between the two models is one way – the PIC results do not change the MHD results. In these calculations, we use the UCLA global MHD code and the iPic3D implicit particle-in-cell code. In this paper we describe this technique in detail. As an example of this approach, we present PIC results on reconnection in the magnetotail during a magnetospheric substorm.

Dust vortices with a void at the centre are reported in this paper. The role of the spatial variation of the plasma potential in the rotation of dust particles is studied in a parallel plate glow discharge plasma. Probe measurements reveal the existence of a local potential minimum in the region of formation of the dust vortex. The minimum in the potential well attracts positively charged ions, while it repels the negatively charged dust particles. Dust rotation is caused by the interplay of the two oppositely directed ion drag and Coulomb forces. The balance between these two forces is found to play a major role in the radial confinement of the dust particles above the cathode surface. Evolution of the dust vortex is studied by increasing the discharge current from 15 to 20 mA. The local minimum of the potential profile is found to coincide with the location of the dust vortex for both values of discharge currents. Additionally, it is found that the size of the dust vortex as well as the void at the centre increases with the discharge current.

The plasma model used in a recent paper by Kamalam et al. (J. Plasma Phys., vol. 84, 2018, 905840406) assumes a Boltzmann description for two hot ion species, in the presence of two adiabatic (fluid) electron species, for the study of obliquely propagating acoustic-type nonlinear solitary waves with respect to a static magnetic field. We argue that the assumption of Boltzmann distributions for the hot ions is incorrect, thus invalidating their conclusions, in particular about the possible occurrence of supersolitons in magnetized plasmas.

The $\unicode[STIX]{x1D6FF}f$-PIC method is widely used for electrostatic particle-in-cell (PIC) simulations. Its basic idea is the ansatz $f=f_{0}+\unicode[STIX]{x1D6FF}f$ ($\unicode[STIX]{x1D6FF}f$-ansatz) where the particle distribution function $f$ is split into a usually time-independent background $f_{0}$ and a time-dependent perturbation $\unicode[STIX]{x1D6FF}f$. In principle, it can also be used for electromagnetic gyrokinetic PIC simulations, but the required number of markers can be so large that PIC simulations become impractical. The reason is a decreasing efficiency of the $\unicode[STIX]{x1D6FF}f$-ansatz for the so-called ‘Hamiltonian formulation’ using $p_{\Vert }$ as a dynamic variable. As a result, the density and current moment of the distribution function develop large statistical errors. To overcome this obstacle we propose to solve the potential equations in the symplectic formulation using $v_{\Vert }$ as a dynamic variable. The distribution function itself is still evolved in the Hamiltonian formulation which is better suited for the numerical integration of the parallel dynamics. The contributions from the full Jacobian of phase space, a finite velocity sphere of the simulation domain and a shifted Maxwellian as a background are considered. Special care has been taken at the discretisation level to make damped magnetohydrodynamics (MHD) mode simulations within a realistic gyrokinetic model feasible. This includes devices like e.g. large tokamaks with a small aspect ratio.

Global magnetohydrodynamic (MHD) instabilities are investigated in a computationally tractable two-dimensional model of the solar tachocline. The model’s differential rotation yields stability in the absence of a magnetic field, but if a magnetic field is present, a joint instability is observed. We analyse the nonlinear development of the instability via fully nonlinear direct numerical simulation, the generalized quasi-linear approximation (GQL) and direct statistical simulation (DSS) based upon low-order expansion in equal-time cumulants. As the magnetic diffusivity is decreased, the nonlinear development of the instability becomes more complicated until eventually a set of parameters is identified that produces a previously unidentified long-term cycle in which energy is transformed from kinetic energy to magnetic energy and back. We find that the periodic transitions, which mimic some aspects of solar variability – for example, the quasiperiodic seasonal exchange of energy between toroidal field and waves or eddies – are unable to be reproduced when eddy-scattering processes are excluded from the model.

We propose that pressure anisotropy causes weakly collisional turbulent plasmas to self-organize so as to resist changes in magnetic-field strength. We term this effect ‘magneto-immutability’ by analogy with incompressibility (resistance to changes in pressure). The effect is important when the pressure anisotropy becomes comparable to the magnetic pressure, suggesting that in collisionless, weakly magnetized (high-$\unicode[STIX]{x1D6FD}$) plasmas its dynamical relevance is similar to that of incompressibility. Simulations of magnetized turbulence using the weakly collisional Braginskii model show that magneto-immutable turbulence is surprisingly similar, in most statistical measures, to critically balanced magnetohydrodynamic turbulence. However, in order to minimize magnetic-field variation, the flow direction becomes more constrained than in magnetohydrodynamics, and the turbulence is more strongly dominated by magnetic energy (a non-zero ‘residual energy’). These effects represent key differences between pressure-anisotropic and fluid turbulence, and should be observable in the $\unicode[STIX]{x1D6FD}\gtrsim 1$ turbulent solar wind.

We carry out expansions of non-symmetric toroidal ideal magnetohydrodynamic (MHD) equilibria with nested flux surfaces about a periodic cylinder, in which physical quantities are periodic of period $2\unicode[STIX]{x03C0}$ in the cylindrical angle $\unicode[STIX]{x1D703}$ and $z$. The cross-section of a flux surface at a constant toroidal angle is assumed to be approximately circular, and data are given on the cylindrical flux surface $r=1$. Furthermore, we assume that the magnetic field lines are closed on the lowest-order flux surface, and the magnetic shear is relatively small. We extend earlier work in a flat torus by Weitzner (Phys. Plasmas, vol. 23, 2016, 062512) and demonstrate that a power series expansion can be carried out to all orders using magnetic flux as an expansion parameter. The cylindrical metric introduces certain new features to the expansions compared to the flat torus. However, the basic methodology of dealing with resonance singularities remains the same. The results, even though lacking convergence proofs, once again support the possibility of smooth, low-shear non-symmetric toroidal MHD equilibria.

An iterative formulation of an arbitrarily short-wavelength solver for global gyrokinetic simulations is suggested. The solver is verified against solutions of the dispersion relation. It can be used to treat the nonlinear polarisation density which is important at the plasma edge. In the linear case, the solver is shown to be computationally efficient.

An introduction to the physical interpretation of the Coulomb logarithm is given with particular emphasis on the quantum-mechanical corrections that are required at high temperatures. Excerpts from the literature are used to emphasize the historical understanding of the topic, which emerged more than a half-century ago. Several misinterpretations are noted. Quantum-mechanical effects are related to diffraction by scales of the order of the Debye screening length; they are not due to quantum uncertainty related to the much smaller distance of closest approach.

This tutorial covers validation of gyrokinetic turbulent-transport models via comparison of measured turbulence with realistic simulations of fusion plasmas. It presents a brief history of validation of gyrokinetic simulations, the principal challenges encountered, a limited survey of common turbulence diagnostics used on tokamaks and stellarators, an overview of the fundamentals of synthetic diagnostic models and a discussion of frontiers in turbulent-transport model validation.