A brief critique is presented of some different classes of magnetohydrodynamic equilibrium solutions based on their continuity properties and whether the magnetic field is integrable or not. A generalized energy functional is introduced that is comprised of alternating ideal regions, with nested flux surfaces with an irrational rotational transform, and Taylor-relaxed regions, possibly with magnetic islands and chaos. The equilibrium states have globally continuous magnetic fields, and may be constructed for arbitrary three-dimensional plasma boundaries and appropriately prescribed pressure and rotational-transform profiles.

A one-dimensional electrostatic filtered Vlasov–Poisson simulation study is discussed. The transition from persisting to arrested Landau damping that is produced by increasing the strength of a sinusoidal perturbation on a background Vlasov–Poisson equilibrium is explored. Emphasis is placed on observed features of the electron phase-space distribution when the perturbation strength is near the transition value. A single ubiquitous waveform is found perturbing the space-averaged phase-space distribution at almost any time in all of the simulations; the sole exception is the saturation stage that can occur at the end of the arrested damping scenario. This waveform contains relatively strong, very narrow structures in velocity bracketing $\pm v_{\text{res}}$ – the velocities at which electrons must move to traverse the dominant field mode wavelength in one of its oscillation periods – and propagating with $\pm v_{\text{res}}$ respectively. Local streams of electrons are found in these structures crossing the resonant velocities from low speed to high speed during Landau damping and from high speed to low speed during Landau growth. At the arrest time, when the field strength is briefly constant, these streams vanish. It is conjectured that the expected transfer of energy between electrons and field during Landau growth or damping has been visualized for the first time. No evidence is found in the phase-space distribution to support recent well-established discoveries of a second-order phase transition in the electric field evolution. While trapping is known to play a role for larger perturbation strengths, it is shown that trapping plays no role at any time in any of the simulations near the transition perturbation strength.

The energy budget of a collisionless plasma subject to electrostatic fluctuations is considered, and the excess of thermal energy over the minimum accessible to it under various constraints that limit the possible forms of plasma motion is calculated. This excess measures how much thermal energy is ‘available’ for conversion into plasma instabilities, and therefore constitutes a nonlinear measure of plasma stability. A distribution function with zero available energy defines a ‘ground state’ in the sense that its energy cannot decrease by any linear or nonlinear plasma motion. In a Vlasov plasma with small density and temperature fluctuations, the available energy is proportional to the mean square of these quantities, and exceeds the corresponding energy in ideal or resistive magnetohydrodynamics. If the first or second adiabatic invariant is conserved, ground states generally have inhomogeneous density and temperature. Magnetically confined plasmas are usually not in any ground state, but certain types of stellarator plasmas are so with respect to fluctuations that conserve both these adiabatic invariants, making the plasma linearly and nonlinearly stable to such fluctuations. Similar stability properties can also be enjoyed by plasmas confined by a dipole magnetic field.