1 Introduction
In a previous publication (Helander Reference Helander2017), hereafter referred to as I, a quantity called the ‘available energy’ was proposed as a measure of nonlinear plasma stability. The available energy is defined as the difference between the total plasma energy and the lowest value it could possibly attain under the various constraints that limit the motion of the plasma. It thus measures how much energy can, in principle, be released from the plasma in the form of linear and nonlinear instabilities. It is analogous to a quantity called the ‘available potential energy’ in meteorology, which is defined as the difference between the potential energy of the atmosphere and the minimum attainable by any adiabatic redistribution of mass (Lorenz Reference Lorenz1955). This analogy is made mathematically explicit in the appendix A.
In the present paper, the concept of available energy is explored further. Its simplest version was considered many years ago by Gardner (Reference Gardner1963), who identified the lowest-energy state of a Vlasov plasma under the sole constraint that its motion should satisfy Liouville’s theorem. As pointed out in I, the concept of available energy becomes much more interesting and relevant if one also accounts for the fact that adiabatic invariants are usually conserved.
Indeed, most important types of instabilities and turbulence affecting magnetically confined plasmas are characterised by frequencies below the ion gyrofrequency 
$\unicode[STIX]{x1D6FA}_{i}$, implying that the magnetic moment,
$$\begin{eqnarray}\unicode[STIX]{x1D707}_{a}=\frac{m_{a}v_{\bot }^{2}}{2B},\end{eqnarray}$$ of each particle is conserved. Here 
$m_{a}$ denotes the mass of the species 
$a$ in question and 
$v_{\bot }$ the speed perpendicular to the magnetic field 
$\boldsymbol{B}=B\boldsymbol{b}$. If, moreover, the frequency is below the electron bounce/transit frequency 
$\unicode[STIX]{x1D714}_{be}$, the parallel adiabatic invariant
$$\begin{eqnarray}J_{a}=\int m_{a}v_{\Vert }\,\text{d}l\end{eqnarray}$$ is conserved for the electrons (
$a=e$), where 
$l$ denotes the arc length along the magnetic field and
$$\begin{eqnarray}m_{a}v_{\Vert }=\sqrt{2m_{a}(H_{a}-\unicode[STIX]{x1D707}_{a}B-e_{a}\unicode[STIX]{x1D6F7})},\end{eqnarray}$$the parallel momentum. The energy
$$\begin{eqnarray}H_{a}=\frac{m_{a}v^{2}}{2}+e_{a}\unicode[STIX]{x1D6F7},\end{eqnarray}$$ (where 
$\unicode[STIX]{x1D6F7}$ denotes the electrostatic potential) and 
$\unicode[STIX]{x1D707}_{a}$ are kept constant in the integration, which is taken between two consecutive bounce points, defined by 
$v_{\Vert }=0$.
For drift-wave instabilities and turbulence (e.g. ion-temperature-gradient and trapped-electron modes), the characteristic frequency 
$\unicode[STIX]{x1D714}$ usually satisfies the ordering
$$\begin{eqnarray}\unicode[STIX]{x1D714}_{bi}\sim \unicode[STIX]{x1D714}\ll \unicode[STIX]{x1D714}_{be}\lesssim \unicode[STIX]{x1D6FA}_{i}\ll \unicode[STIX]{x1D6FA}_{e},\end{eqnarray}$$ so that 
$\unicode[STIX]{x1D707}_{e}$, 
$\unicode[STIX]{x1D707}_{i}$ and 
$J_{e}$ are conserved, but not 
$J_{i}$. (The subscript 
$i$ refers to the ions, which will be assumed to be singly charged.) Different constraints thus pertain to the electron and ion distribution functions, and their lowest-energy states will therefore be different. However, the plasma must remain quasineutral, which imposes an additional constraint on the evolution of the plasma that was not accounted for in I.
The aim of the present paper is to discuss how this and other relevant constraints affect the available energy of magnetically confined plasmas. The energy available to fluctuations subject to constraints such as (1.5) is lower than that available to less constrained fluctuations, and our goal is to develop a general method for incorporating various types of constraints. We begin in § 2 by briefly recapitulating the basic results of I, deriving them in a particularly expeditious way. This derivation naturally lends itself to the inclusion of additional constraints, and in the following several sections we show how these can be incorporated. In § 3, the available energy of a single plasma species with conserved 
$\unicode[STIX]{x1D707}$ and 
$J$ is calculated; in §§ 4 and 5 the constraint of a fixed density profile is explored; and in § 6 ground states of a plasma satisfying the ordering (1.5) are investigated, where different constraints apply to the ions and electrons.
2 Minimum-energy states
As in I, let us first consider the distribution function 
$f(\boldsymbol{x},t)$ of a single particle species governed by the Vlasov equation or any other kinetic equation satisfying the Liouville theorem,
$$\begin{eqnarray}\displaystyle & \displaystyle \frac{\unicode[STIX]{x2202}f}{\unicode[STIX]{x2202}t}+{\dot{x}}_{k}\frac{\unicode[STIX]{x2202}f}{\unicode[STIX]{x2202}x_{k}}=0, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}x_{k}}({\dot{x}}_{k}\sqrt{g})=0, & \displaystyle\end{eqnarray}$$ where the phase-space coordinates 
$\boldsymbol{x}$ are arbitrary except that the Jacobian 
$\sqrt{g}$ is assumed to be independent of time. Summation over repeated indices is understood. Liouville’s theorem implies that the flux of particles in phase space is incompressible, so that, for any positive number 
$\unicode[STIX]{x1D719}$, the volume of the set of points satisfying 
$f(\boldsymbol{x},t)>\unicode[STIX]{x1D719}$ is constant, i.e.
$$\begin{eqnarray}\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}t}H[f(\boldsymbol{x},t),\unicode[STIX]{x1D719}]=0,\end{eqnarray}$$ where the functional 
$H$ is defined by
$$\begin{eqnarray}H[f(\boldsymbol{x},t),\unicode[STIX]{x1D719}]=\int \unicode[STIX]{x1D6E9}[f(\boldsymbol{x},t)-\unicode[STIX]{x1D719}]\sqrt{g}\,\text{d}\boldsymbol{x}\end{eqnarray}$$ and 
$\unicode[STIX]{x1D6E9}$ denotes the Heaviside step function.
Now, let 
$\unicode[STIX]{x1D716}(\boldsymbol{x})$ denote the particle energy as a function of the phase-space coordinates 
$\boldsymbol{x}$. (Typically 
$\unicode[STIX]{x1D716}=mv^{2}/2$, where 
$v$ denotes speed.) The total energy associated with the distribution function 
$f$ is then
$$\begin{eqnarray}E[f]=\int \unicode[STIX]{x1D716}f\sqrt{g}\,\text{d}\boldsymbol{x},\end{eqnarray}$$ and we may ask for the minimum energy that can be attained starting from some initial state 
$f(\boldsymbol{x},0)$. We denote the energy-minimising distribution function by 
$f_{0}(\boldsymbol{x})$ and note that, because of the constraint (2.3), we require
$$\begin{eqnarray}H[f_{0}(\boldsymbol{x}),\unicode[STIX]{x1D719}]=H[f(\boldsymbol{x},0),\unicode[STIX]{x1D719}]\equiv H_{0}(\unicode[STIX]{x1D719}).\end{eqnarray}$$ We can account for this constraint by introducing a continuous set of Lagrange multipliers 
$\unicode[STIX]{x1D706}(\unicode[STIX]{x1D719})$ and minimising the functional
$$\begin{eqnarray}W[f_{0},\unicode[STIX]{x1D706}]=E[f_{0}]+\int _{0}^{\infty }\unicode[STIX]{x1D706}(\unicode[STIX]{x1D719})\,\text{d}\unicode[STIX]{x1D719}\left[\int \unicode[STIX]{x1D6E9}[f_{0}(\boldsymbol{x})-\unicode[STIX]{x1D719}]\sqrt{g}\,\text{d}\boldsymbol{x}-H_{0}(\unicode[STIX]{x1D719})\right].\end{eqnarray}$$ The minimisation with respect to 
$\unicode[STIX]{x1D706}$ then ensures that the constraint (2.6) is satisfied, whereas the variation with respect to 
$f_{0}$ gives
$$\begin{eqnarray}\unicode[STIX]{x1D6FF}W=\int \unicode[STIX]{x1D6FF}f_{0}\left[\unicode[STIX]{x1D716}+\unicode[STIX]{x1D706}(f_{0})\right]\sqrt{g}\,\text{d}\boldsymbol{x}=0,\end{eqnarray}$$ implying that 
$\unicode[STIX]{x1D706}(f_{0})=-\unicode[STIX]{x1D716}(\boldsymbol{x})$. Inverting this relation, we come to the conclusion that the lowest-energy state can only depend on the phase-space coordinates through 
$\unicode[STIX]{x1D716}(\boldsymbol{x})$, i.e. it must be possible to write 
$f_{0}$ as
$$\begin{eqnarray}f_{0}(\boldsymbol{x})=F[\unicode[STIX]{x1D716}(\boldsymbol{x})],\end{eqnarray}$$ for some function 
$F$.
As originally shown by Gardner (Reference Gardner1963), this function should be monotonically decreasing. Otherwise, the distribution function does not correspond to a minimum-energy state, because the energy can be lowered by interchanging two neighbouring phase-space elements of equal volume 
$\text{d}\unicode[STIX]{x1D70F}=\sqrt{g}\,\text{d}\boldsymbol{x}$ but different energies. To see this, let us denote the distribution function before the interchange by 
$f_{0}$ and afterwards by 
$f_{1}$. Then
$$\begin{eqnarray}\displaystyle & \displaystyle f_{1}(\boldsymbol{x}_{b})=f_{0}(\boldsymbol{x}_{a}), & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle f_{1}(\boldsymbol{x}_{a})=f_{0}(\boldsymbol{x}_{b}), & \displaystyle\end{eqnarray}$$ where points in the two respective phase-space elements are denoted by 
$\boldsymbol{x}_{a,b}$. As a result of the interchange, the energy changes by the amount
$$\begin{eqnarray}\unicode[STIX]{x0394}E=\int \unicode[STIX]{x1D716}(f_{1}-f_{0})\sqrt{g}\,\text{d}\boldsymbol{x}=-[\unicode[STIX]{x1D716}(\boldsymbol{x}_{a})-\unicode[STIX]{x1D716}(\boldsymbol{x}_{b})][f_{0}(\boldsymbol{x}_{a})-f_{0}(\boldsymbol{x}_{b})]\,\text{d}\unicode[STIX]{x1D70F},\end{eqnarray}$$and becomes equal to
$$\begin{eqnarray}\unicode[STIX]{x0394}E=-[\unicode[STIX]{x1D716}(\boldsymbol{x}_{a})-\unicode[STIX]{x1D716}(\boldsymbol{x}_{b})]^{2}F_{0}^{\prime }(\unicode[STIX]{x1D716})\end{eqnarray}$$ if the distribution function is of the form (2.9). Hence, it is clear that 
$F^{\prime }(\unicode[STIX]{x1D716})$ must be negative (or zero) everywhere in a lowest-energy state, which we called the ground state in I.
Since 
$F$ is a monotonically decreasing function we have
$$\begin{eqnarray}\unicode[STIX]{x1D6E9}\left[F(\unicode[STIX]{x1D716}(\boldsymbol{x}))-\unicode[STIX]{x1D719}\right]=\unicode[STIX]{x1D6E9}\left[w-\unicode[STIX]{x1D716}(\boldsymbol{x})\right],\end{eqnarray}$$ where 
$F(w)=\unicode[STIX]{x1D719}$. It follows that this function is determined by the condition
$$\begin{eqnarray}H_{0}[F(w)]=\unicode[STIX]{x1D6FA}(w),\end{eqnarray}$$ for each 
$w\in [0,\infty )$, where
$$\begin{eqnarray}\unicode[STIX]{x1D6FA}(w)=\int \unicode[STIX]{x1D6E9}\left[w-\unicode[STIX]{x1D716}(\boldsymbol{x})\right]\sqrt{g}\,\text{d}\boldsymbol{x},\end{eqnarray}$$ measures the volume of the subset of phase space in which 
$\unicode[STIX]{x1D716}(\boldsymbol{x})<w$. Equation (2.15) states that this volume must equal the volume in which the distribution function exceeds 
$F(y)$, as follows from Liouville’s theorem.
Equation (2.15) determines the ground states of a plasma whose evolution is only constrained by Liouville’s theorem. As demonstrated in I and discussed further below, it is not difficult to account for more constraints by using additional Lagrange multipliers.
The available energy is defined as the difference in energy between the initial state and the ground state,
$$\begin{eqnarray}A=\int (f-f_{0})\unicode[STIX]{x1D716}\sqrt{g}\,\text{d}\boldsymbol{x}.\end{eqnarray}$$ It evidently constitutes an upper bound on the energy that can be extracted from the plasma by linear and nonlinear instabilities, but this upper bound is in general unattainable. For instance, it does not account for the fact that the fluctuations created by the instability in question must also be contained in 
$f$. Moreover, because coherent waves tend to diffuse particles or because internal wave turbulence interacts incoherently with particles, particles tend to diffuse in phase space (looking at realistic granularity) rather than to conserve the phase-space density. Gardner restacking therefore overestimates the available free energy (Fisch & Rax Reference Fisch and Rax1993; Hay, Schiff & Fisch Reference Hay, Schiff and Fisch2015).
3 Conservation of 
$\unicode[STIX]{x1D707}$ and $J$
As already remarked in the Introduction and discussed at length in I, constraints implied by the adiabatic invariants 
$\unicode[STIX]{x1D707}$ and 
$J$ can place severe limits on the available energy. We are particularly interested in the ordering (1.5), according to which both these invariants are conserved for trapped electrons.Footnote 1 In the present section, we calculate the ground state of any such particle species for which 
$\unicode[STIX]{x1D707}$ and 
$J$ are conserved.
3.1 The role of omnigeneity
For a plasma with magnetic surfaces, it is useful to write the magnetic field as 
$\boldsymbol{B}=\unicode[STIX]{x1D735}\unicode[STIX]{x1D713}\times \unicode[STIX]{x1D735}\unicode[STIX]{x1D6FC}$, where 
$\unicode[STIX]{x1D713}$ labels the magnetic surfaces and 
$\unicode[STIX]{x1D6FC}\in [0,2\unicode[STIX]{x03C0})$ the different field lines thereon. The distribution function is then usefully considered as a function of the phase-space coordinates 
$(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)$, and we shall take it to be initially Maxwellian,
$$\begin{eqnarray}f(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)=f_{M}[\unicode[STIX]{x1D716}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J),\unicode[STIX]{x1D713}]=n(\unicode[STIX]{x1D713})\left(\frac{m}{2\unicode[STIX]{x03C0}T(\unicode[STIX]{x1D713})}\right)^{3/2}\exp (-\unicode[STIX]{x1D716}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)/T(\unicode[STIX]{x1D713})),\end{eqnarray}$$ with a density 
$n(\unicode[STIX]{x1D713})$ and temperature 
$T(\unicode[STIX]{x1D713})$ that are constant on each flux surface 
$\unicode[STIX]{x1D713}$.Footnote 2 We shall assume that the distribution function evolves on a time scale longer than the bounce time for trapped particles, so that 
$\unicode[STIX]{x1D707}$ and 
$J$ remain constant for each such particle and the distribution function is independent of the position 
$l$ along 
$\boldsymbol{B}$. Our main task is to find the ground state as described by some distribution function 
$f_{0}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)$ of the trapped particles. As shown in I, ground states can be written as functions of energy and the conserved quantities alone, i.e. it must be possible to write 
$f_{0}$ as
$$\begin{eqnarray}f_{0}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)=F[\unicode[STIX]{x1D716}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J),\unicode[STIX]{x1D707},J]\end{eqnarray}$$ for some function 
$F$.
Before proceeding to find this function, it is useful ask the question as to whether a Maxwellian flux function (3.1) can represent a ground state. This can evidently only be the case if
$$\begin{eqnarray}\left(\frac{\unicode[STIX]{x2202}f}{\unicode[STIX]{x2202}\unicode[STIX]{x1D6FC}}\right)_{\unicode[STIX]{x1D716},\unicode[STIX]{x1D707},J}=0,\end{eqnarray}$$but for a Maxwellian the left-hand side is equal to
$$\begin{eqnarray}\left(\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D6FC}}f_{M}(\unicode[STIX]{x1D716},\unicode[STIX]{x1D713})\right)_{\unicode[STIX]{x1D716},\unicode[STIX]{x1D707},J}=\frac{\unicode[STIX]{x2202}f_{M}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}\left(\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D6FC}}\right)_{\unicode[STIX]{x1D716},\unicode[STIX]{x1D707},J}=-\frac{\unicode[STIX]{x2202}f_{M}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}\left(\frac{\unicode[STIX]{x2202}J}{\unicode[STIX]{x2202}\unicode[STIX]{x1D6FC}}\right)_{\unicode[STIX]{x1D716},\unicode[STIX]{x1D707},\unicode[STIX]{x1D713}}\left/\left(\frac{\unicode[STIX]{x2202}J}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}\right)_{\unicode[STIX]{x1D716},\unicode[STIX]{x1D707},\unicode[STIX]{x1D6FC}}\right.\!.\end{eqnarray}$$ We thus conclude that (3.1) can only represent a ground state if either the density and temperature are constant, so that 
$\unicode[STIX]{x2202}f_{M}/\unicode[STIX]{x2202}\unicode[STIX]{x1D713}=0$, or the magnetic field is omnigenous, 
$\unicode[STIX]{x2202}J/\unicode[STIX]{x2202}\unicode[STIX]{x1D6FC}=0$. By definition, in an omnigenous field, 
$J$ can be written as a function of energy, magnetic moment and the flux-surface label 
$\unicode[STIX]{x1D713}$, without any dependence on 
$\unicode[STIX]{x1D6FC}$, i.e. we can write 
$J=J(\unicode[STIX]{x1D716},\unicode[STIX]{x1D707},\unicode[STIX]{x1D713})$. Since trapped-particle orbits precess on surfaces of constant 
$J$, omnigenous fields have the property that collisionless particle trajectories are radially well confined, making the neoclassical transport small (Hall & McNamara Reference Hall and McNamara1975; Cary & Shasharina Reference Cary and Shasharina1997; Helander Reference Helander2014b). It is interesting to note that there is thus a connection between low neoclassical transport – an important goal of stellarator optimisation – and small available energy, which could be beneficial for turbulence reduction.
Of course, omnigeneity is merely a necessary, but not sufficient, condition for a ground state of the form (3.1). The further condition that 
$\unicode[STIX]{x2202}F/\unicode[STIX]{x2202}\unicode[STIX]{x1D716}\leqslant 0$ was shown in I to be equivalent to
$$\begin{eqnarray}\frac{\unicode[STIX]{x1D714}_{\ast }^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}}\leqslant 1,\end{eqnarray}$$where we have written
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D714}_{\ast }=\frac{T}{q}\frac{\text{d}\ln n}{\text{d}\unicode[STIX]{x1D713}}, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D714}_{\ast }^{T}=\unicode[STIX]{x1D714}_{\ast }\left[1+\unicode[STIX]{x1D702}\left(\frac{\unicode[STIX]{x1D716}}{T}-\frac{3}{2}\right)\right], & \displaystyle\end{eqnarray}$$ and 
$\unicode[STIX]{x1D702}=(\text{d}\ln T/\text{d}\unicode[STIX]{x1D713})/(\text{d}\ln n/\text{d}\unicode[STIX]{x1D713})$ as in I but we now denote the charge by 
$q$. We have also introduced the bounce-averaged precession frequency for trapped particles,
$$\begin{eqnarray}\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}=\frac{1}{\unicode[STIX]{x1D70F}_{b}}\int _{l_{1}}^{l_{2}}(\boldsymbol{v}_{d}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D6FC})\frac{\text{d}l}{v_{\Vert }}=\frac{1}{q}\left(\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D716}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}\right)_{\unicode[STIX]{x1D707},J,\unicode[STIX]{x1D6FC}},\end{eqnarray}$$ where 
$\boldsymbol{v}_{d}$ denotes the drift velocity and
$$\begin{eqnarray}\unicode[STIX]{x1D70F}_{b}=\int _{l_{1}}^{l_{2}}\frac{\text{d}l}{|v_{\Vert }|}=\left(\frac{\unicode[STIX]{x2202}J}{\unicode[STIX]{x2202}\unicode[STIX]{x1D716}}\right)_{\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707}}\end{eqnarray}$$ the bounce time between two consecutive bounce points, which have been denoted by 
$l_{1,2}$. The precession frequency can also be written as
$$\begin{eqnarray}\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}=-\frac{1}{q}\left(\frac{\unicode[STIX]{x2202}J}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}\right)_{\unicode[STIX]{x1D716},\unicode[STIX]{x1D707},\unicode[STIX]{x1D6FC}}\left/\left(\frac{\unicode[STIX]{x2202}J}{\unicode[STIX]{x2202}\unicode[STIX]{x1D716}}\right)_{\unicode[STIX]{x1D707},\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC}}\right.\!,\end{eqnarray}$$ and is thus seen to be independent of 
$\unicode[STIX]{x1D6FC}$ in an omnigenous field, where 
$\unicode[STIX]{x2202}J/\unicode[STIX]{x2202}\unicode[STIX]{x1D6FC}=0$.
3.2 Equations for the ground state
We now turn to the calculation of the available energy when the distribution function (3.1) is not a ground state, i.e. when the condition (3.5) is not satisfied. For simplicity, we only consider omnigenous magnetic fields, such as that of a tokamak or well-optimised stellarator.
When certain quantities 
$\boldsymbol{y}$ such as 
$\unicode[STIX]{x1D707}$ and 
$J$ are conserved, it is useful to employ these as (some of the) phase-space coordinates. We thus write the latter as 
$\boldsymbol{x}=(\boldsymbol{y},\boldsymbol{z})$ and recall from I that the minimum-energy state is given by a simple generalisation of (2.15),
$$\begin{eqnarray}H[F(w,\boldsymbol{y}),\boldsymbol{y}]=\unicode[STIX]{x1D6FA}(w,\boldsymbol{y}),\end{eqnarray}$$where
$$\begin{eqnarray}\displaystyle & \displaystyle H(\unicode[STIX]{x1D719},\boldsymbol{y})=\int \unicode[STIX]{x1D6E9}[f(\boldsymbol{y},\boldsymbol{z})-\unicode[STIX]{x1D719}]\sqrt{g}\,\text{d}\boldsymbol{z}, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D6FA}(w,\boldsymbol{y})=\int \unicode[STIX]{x1D6E9}[w-\unicode[STIX]{x1D716}(\boldsymbol{y},\boldsymbol{z})]\sqrt{g}\,\text{d}\boldsymbol{z}. & \displaystyle\end{eqnarray}$$ In the present case, we have two adiabatic invariants for the electrons, 
$\boldsymbol{y}=(\unicode[STIX]{x1D707},J)$, and the phase-space volume element is (see I)
$$\begin{eqnarray}\text{d}\unicode[STIX]{x1D6E4}=\text{d}\boldsymbol{r}\text{d}\boldsymbol{v}=\frac{4\unicode[STIX]{x03C0}\text{d}l}{m^{2}|v_{\Vert }|\unicode[STIX]{x1D70F}_{b}}\,\text{d}\unicode[STIX]{x1D707}\,\text{d}J\,\text{d}\unicode[STIX]{x1D713}\,\text{d}\unicode[STIX]{x1D6FC}.\end{eqnarray}$$We thus have the relations
$$\begin{eqnarray}\displaystyle & \displaystyle H(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},J)=\frac{4\unicode[STIX]{x03C0}}{m^{2}}\oint \text{d}\unicode[STIX]{x1D6FC}\int \unicode[STIX]{x1D6E9}\left(f-\unicode[STIX]{x1D719}\right)\text{d}\unicode[STIX]{x1D713}, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D6FA}(w,\unicode[STIX]{x1D707},J)=\frac{4\unicode[STIX]{x03C0}}{m^{2}}\oint \text{d}\unicode[STIX]{x1D6FC}\int \unicode[STIX]{x1D6E9}\left[w-\unicode[STIX]{x1D716}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)\right]\text{d}\unicode[STIX]{x1D713}, & \displaystyle\end{eqnarray}$$ which together with (3.11) determine the ground state. Here and in the following, integrals over 
$\unicode[STIX]{x1D6FC}$ are to be taken over all field lines on the flux surface, 
$0\leqslant \unicode[STIX]{x1D6FC}\leqslant 2\unicode[STIX]{x03C0}$, and we have used the fact that 
$f$ does not vary much along the field line, i.e. 
$f$ is nearly independent of 
$l$ at fixed 
$(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)$.
Rather than solving (3.11) directly, we consider its derivative (Dodin & Fisch Reference Dodin and Fisch2005),
$$\begin{eqnarray}\left(\frac{\unicode[STIX]{x2202}F}{\unicode[STIX]{x2202}w}\right)_{\unicode[STIX]{x1D707},J}=\frac{\unicode[STIX]{x1D6FA}_{1}(w,\unicode[STIX]{x1D707},J)}{H_{1}\left[F(w,\unicode[STIX]{x1D707},J),\unicode[STIX]{x1D707},J\right]},\end{eqnarray}$$ where the subscript 
$1$ denotes a derivative with respect to the first argument of the function in question. Thus
$$\begin{eqnarray}\unicode[STIX]{x1D6FA}_{1}(w,\unicode[STIX]{x1D707},J)=\oint \text{d}\unicode[STIX]{x1D6FC}\int \unicode[STIX]{x1D6FF}\left[w-\unicode[STIX]{x1D716}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)\right]\frac{4\unicode[STIX]{x03C0}}{m^{2}}\,\text{d}\unicode[STIX]{x1D713}=\frac{4\unicode[STIX]{x03C0}}{m^{2}}\oint \frac{\text{d}\unicode[STIX]{x1D6FC}}{|q\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}|},\end{eqnarray}$$where we have used (3.8) and
$$\begin{eqnarray}\int \unicode[STIX]{x1D6FF}\left[w-\unicode[STIX]{x1D716}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)\right]\text{d}\unicode[STIX]{x1D713}=\frac{1}{|\unicode[STIX]{x2202}\unicode[STIX]{x1D716}/\unicode[STIX]{x2202}\unicode[STIX]{x1D713}|}.\end{eqnarray}$$For the denominator of (3.17) we need
$$\begin{eqnarray}H_{1}\left(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},J\right)=-\oint \text{d}\unicode[STIX]{x1D6FC}\int \unicode[STIX]{x1D6FF}[f_{M}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)-\unicode[STIX]{x1D719}]\frac{4\unicode[STIX]{x03C0}}{m^{2}}\,\text{d}\unicode[STIX]{x1D713}=-\frac{4\unicode[STIX]{x03C0}}{m^{2}}\oint \left|\left(\frac{\unicode[STIX]{x2202}f_{M}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}\right)_{\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J}\right|^{-1}\text{d}\unicode[STIX]{x1D6FC},\end{eqnarray}$$which is equal to
$$\begin{eqnarray}H_{1}\left(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},J\right)=\oint \frac{4\unicode[STIX]{x03C0}T\,\text{d}\unicode[STIX]{x1D6FC}}{m^{2}\left|q\left(\unicode[STIX]{x1D714}_{\ast }^{T}-\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}\right)\right|f_{M}},\end{eqnarray}$$where we have used
$$\begin{eqnarray}\left(\frac{\unicode[STIX]{x2202}f_{M}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}\right)_{\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J}=\frac{qf_{M}}{T}(\unicode[STIX]{x1D714}_{\ast }^{T}-\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}).\end{eqnarray}$$Our (3.17) for the ground state thus becomes
$$\begin{eqnarray}\left(\frac{\unicode[STIX]{x2202}F}{\unicode[STIX]{x2202}w}\right)_{\unicode[STIX]{x1D707},J}=-\frac{f_{M}}{T}\left|\frac{\unicode[STIX]{x1D714}_{\ast }^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}}-1\right|.\end{eqnarray}$$3.3 Available energy
We can solve this equation in the neighbourhood of a flux surface, 
$\unicode[STIX]{x1D713}_{0}-\unicode[STIX]{x0394}\unicode[STIX]{x1D713}<\unicode[STIX]{x1D713}<\unicode[STIX]{x1D713}_{0}+\unicode[STIX]{x0394}\unicode[STIX]{x1D713}$, by expanding in 
$\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0}$,
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D716}(\unicode[STIX]{x1D707},J,\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC})=\unicode[STIX]{x1D716}_{0}(\unicode[STIX]{x1D707},J)+\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D716}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})+\cdots \,, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle f_{M}(\unicode[STIX]{x1D707},J,\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC})=f_{M0}+\frac{\unicode[STIX]{x2202}f_{M}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})+\cdots \,, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle F(\unicode[STIX]{x1D707},J,\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC})=F_{0}+\frac{\unicode[STIX]{x2202}F}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})+\cdots \,, & \displaystyle\end{eqnarray}$$ where all derivatives and all quantities with a subscript 
$0$ are to be evaluated at 
$\unicode[STIX]{x1D713}=\unicode[STIX]{x1D713}_{0}$, and we recall (3.8) and
$$\begin{eqnarray}\frac{\unicode[STIX]{x2202}F}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}=-\frac{qf_{M0}}{T}\left|\unicode[STIX]{x1D714}_{\ast }^{T}-\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}\right|.\end{eqnarray}$$ Because the volume element in phase space is given by (3.14) and the total number of particles having any given values of 
$\unicode[STIX]{x1D707}$ and 
$J$ must be the same in the initial state and the ground state, we have
$$\begin{eqnarray}\int _{\unicode[STIX]{x1D713}_{0}-\unicode[STIX]{x0394}\unicode[STIX]{x1D713}}^{\unicode[STIX]{x1D713}_{0}+\unicode[STIX]{x0394}\unicode[STIX]{x1D713}}(f_{M}-F)\,\text{d}\unicode[STIX]{x1D713}=\int _{\unicode[STIX]{x1D713}_{0}-\unicode[STIX]{x0394}\unicode[STIX]{x1D713}}^{\unicode[STIX]{x1D713}_{0}+\unicode[STIX]{x0394}\unicode[STIX]{x1D713}}\left[f_{M0}-F_{0}+\frac{1}{2}\left(\frac{\unicode[STIX]{x2202}^{2}f_{M}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}^{2}}-\frac{\unicode[STIX]{x2202}^{2}F}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}^{2}}\right)(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})^{2}\right]\text{d}\unicode[STIX]{x1D713}=0.\end{eqnarray}$$The available energy can therefore to leading order be calculated from
$$\begin{eqnarray}A=\int \unicode[STIX]{x1D716}(f_{M}-F)\,\text{d}\unicode[STIX]{x1D6E4}\end{eqnarray}$$as follows:
$$\begin{eqnarray}\displaystyle A & \simeq & \displaystyle \int [\unicode[STIX]{x1D716}_{0}+q\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})]\left[f_{M0}-F_{0}+\frac{qf_{M0}}{T}\left(\unicode[STIX]{x1D714}_{\ast }^{T}-\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}+\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}\left|\frac{\unicode[STIX]{x1D714}_{\ast }^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}}-1\right|\right)(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})\right.\nonumber\\ \displaystyle & & \displaystyle \left.\hphantom{\int \left[\unicode[STIX]{x1D716}_{0}+q\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})\right]~}+\frac{1}{2}\left(\frac{\unicode[STIX]{x2202}^{2}f_{M}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}^{2}}-\frac{\unicode[STIX]{x2202}^{2}F}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}^{2}}\right)(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})^{2}\right]\text{d}\unicode[STIX]{x1D6E4}\nonumber\\ \displaystyle & = & \displaystyle \frac{q^{2}}{T}\int f_{M0}\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}^{2}\left(\frac{\unicode[STIX]{x1D714}_{\ast }^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}}-1+\left|\frac{\unicode[STIX]{x1D714}_{\ast }^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}}-1\right|\right)(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})^{2}\,\text{d}\unicode[STIX]{x1D6E4},\end{eqnarray}$$and the result becomes
$$\begin{eqnarray}A=\frac{16\unicode[STIX]{x03C0}q^{2}(\unicode[STIX]{x0394}\unicode[STIX]{x1D713})^{3}}{3m^{2}T}\int f_{M0}\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}^{2}r\left(\frac{\unicode[STIX]{x1D714}_{\ast }^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}}-1\right)\,\text{d}\unicode[STIX]{x1D707}\,\text{d}J\,\text{d}\unicode[STIX]{x1D6FC},\end{eqnarray}$$ where 
$r(x)=x\unicode[STIX]{x1D6E9}(x)$. It can also be written as
$$\begin{eqnarray}A=\frac{4q^{2}(\unicode[STIX]{x0394}\unicode[STIX]{x1D713})^{3}}{3T}U^{\prime }(\unicode[STIX]{x1D713})\int _{\text{trapped}}f_{M0}\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}^{2}r\left(\frac{\unicode[STIX]{x1D714}_{\ast }^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}}-1\right)\text{d}\boldsymbol{v},\end{eqnarray}$$ where 
$U(\unicode[STIX]{x1D713})$ denotes the volume enclosed by the flux surface 
$\unicode[STIX]{x1D713}$.
This expression for the maximum energy that can be extracted from electrons by fluctuations satisfying the ordering (1.5) is one of the principal results of this paper. In order to explore how it scales with plasma parameters and magnetic-field geometry, we note that
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}\sim \frac{T}{m\unicode[STIX]{x1D6FA}aR}, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D714}_{\ast }^{T}\sim \frac{T}{eaL_{\bot }B}, & \displaystyle\end{eqnarray}$$ where 
$a$ denotes the minor radius of the flux surface in question, 
$R$ the radius of curvature of the magnetic field lines and 
$L_{\bot }$ the length scale of the pressure profile, which is assumed to be shorter than 
$R$. The volume enclosed by the flux surface 
$\unicode[STIX]{x1D713}$ is
$$\begin{eqnarray}U(\unicode[STIX]{x1D713})\sim \frac{R\unicode[STIX]{x1D713}}{B},\end{eqnarray}$$ and that of the region between 
$\unicode[STIX]{x1D713}-\unicode[STIX]{x0394}\unicode[STIX]{x1D713}$ and 
$\unicode[STIX]{x1D713}+\unicode[STIX]{x0394}\unicode[STIX]{x1D713}$
$$\begin{eqnarray}V\sim aR\unicode[STIX]{x0394}r,\end{eqnarray}$$ where 
$\unicode[STIX]{x0394}\unicode[STIX]{x1D713}\sim aB\unicode[STIX]{x0394}r$. Within factors of order unity, the available energy (3.32) is thus given by
$$\begin{eqnarray}\frac{A}{V}\sim f_{t}nT\frac{(\unicode[STIX]{x0394}r)^{2}}{RL_{\bot }},\end{eqnarray}$$ where 
$f_{t}$ denotes the fraction of electrons that are trapped in regions of average unfavourable curvature.
The result (3.32) for the available energy of electrons conserving 
$\unicode[STIX]{x1D707}$ and 
$J$ is always smaller than the corresponding energy without these constraints. The latter was calculated in I and was found to be equal to
$$\begin{eqnarray}\frac{A}{V}=nT\left\langle \frac{1}{2}\left(\frac{\unicode[STIX]{x1D6FF}n}{n}\right)^{2}+\frac{3}{4}\left(\frac{\unicode[STIX]{x1D6FF}T}{T}\right)^{2}\right\rangle \sim nT\frac{(\unicode[STIX]{x0394}r)^{2}}{L_{\bot }^{2}}\end{eqnarray}$$ for a Maxwellian plasma with small density and temperature variations, 
$\unicode[STIX]{x1D6FF}n\ll n$ and 
$\unicode[STIX]{x1D6FF}T\ll T$, respectively.
4 Fixed density profile
In many numerical simulations of plasma instabilities and turbulence, the electrons are assumed to be ‘adiabatic’, meaning that the electron density fluctuations 
$\unicode[STIX]{x1D6FF}n_{e}$ are proportional to those of the electrostatic potential 
$\unicode[STIX]{x1D6FF}\unicode[STIX]{x1D6F7}$ according to the linearised Boltzmann relation
$$\begin{eqnarray}\frac{\unicode[STIX]{x1D6FF}n_{e}}{n_{e}}=\frac{e\unicode[STIX]{x1D6FF}\unicode[STIX]{x1D6F7}}{T_{e}},\end{eqnarray}$$ where 
$-e$ denotes the electron charge and 
$T_{e}$ the electron temperature. (In configurations with closed magnetic field lines or magnetic flux surfaces, an appropriate average of 
$\unicode[STIX]{x1D6FF}\unicode[STIX]{x1D6F7}$ should be subtracted from the right-hand side.) Because the 
$\boldsymbol{E}\times \boldsymbol{B}$-drift is then out of phase with the density fluctuations, there is no electron particle flux across the magnetic field and the density profile must remain constant.
In a plasma, as opposed to a numerical simulation, the assumption of ‘adiabatic’ (or, better, Boltzmann-distributed) electrons is valid if the ordering (1.5) holds and, additionally, no (or sufficiently few) electrons are trapped in magnetic wells. In this ordering, the transit frequency of the circulating electrons substantially exceeds the fluctuation frequency. These electrons thus travel many times around the torus on the time scale of one fluctuation, and therefore experience little transport. Mathematically, this conclusion follows, for instance, from a familiar argument for the confinement of circulating particles in a stellarator in the face of cross-field drifts (Helander Reference Helander2014b). The drift velocity, including the 
$\boldsymbol{E}\times \boldsymbol{B}$ drift, is
$$\begin{eqnarray}\boldsymbol{v}_{da}=\frac{v_{\Vert }}{B}\unicode[STIX]{x1D735}\times \left(\frac{v_{\Vert }\boldsymbol{B}}{\unicode[STIX]{x1D6FA}_{a}}\right),\end{eqnarray}$$ where the derivative is taken at constant 
$\unicode[STIX]{x1D707}_{a}$ and 
$H_{a}$. To first order in the smallness of the drift, the net displacement across the flux surface is thus
$$\begin{eqnarray}\unicode[STIX]{x0394}\unicode[STIX]{x1D713}=\int (\boldsymbol{v}_{da}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D713})\frac{\text{d}l}{v_{\Vert }}=\int \unicode[STIX]{x1D735}\boldsymbol{\cdot }\left(\frac{v_{\Vert }\boldsymbol{B}\times \unicode[STIX]{x1D735}\unicode[STIX]{x1D713}}{\unicode[STIX]{x1D6FA}_{a}}\right)\frac{\text{d}l}{B},\end{eqnarray}$$ where the integral is taken along the field line. If the latter wraps many times around the torus, the integral becomes proportional to the flux-surface average of 
$\unicode[STIX]{x1D735}\boldsymbol{\cdot }[(v_{\Vert }/\unicode[STIX]{x1D6FA}_{a})\boldsymbol{B}\times \unicode[STIX]{x1D735}\unicode[STIX]{x1D713}]$, which vanishes since the vector within the divergence is perpendicular to 
$\unicode[STIX]{x1D735}\unicode[STIX]{x1D713}$. An exception occurs close to surfaces on which the rotational transform (winding number of the magnetic field lines) is a rational number. In such regions, it takes many turns around the torus before a field line has come close to most points on the magnetic surface, and the electrons therefore do not respond ‘adiabatically’ despite the ordering (1.5) (Hallatschek & Dorland Reference Hallatschek and Dorland2005; Dominski et al. Reference Dominski, Brunner, Görler, Jenko, Told and Villard2015).
4.1 Ground state
Suppose we ask for the minimum-energy state of the ions in a plasma with fixed density profile. It will in general not be given by (2.9), because the latter does not account for the constraint that the density profile cannot change. Equation (2.9) predicts a state of constant density throughout the plasma, which is not attainable unless the initial state had constant density. Instead, the ground state must be determined by minimising the energy subject to both the constraint (2.6) and the requirement that the density profile is given,
$$\begin{eqnarray}\int f_{0}(\boldsymbol{r},\boldsymbol{v})\,\text{d}\boldsymbol{v}=n(\boldsymbol{r})=\text{fixed}.\end{eqnarray}$$ The latter can be accounted for by an additional set of Lagrange multipliers 
$\unicode[STIX]{x1D705}(\boldsymbol{r})$, so that instead of (2.7) we seek to minimise
$$\begin{eqnarray}U[f_{0},\unicode[STIX]{x1D705},\unicode[STIX]{x1D706}]=W[f_{0},\unicode[STIX]{x1D706}]+\int \unicode[STIX]{x1D705}(\boldsymbol{r})\,\text{d}\boldsymbol{r}\left[\int f_{0}(\boldsymbol{r},\boldsymbol{v})\,\text{d}\boldsymbol{v}-n(\boldsymbol{r})\right].\end{eqnarray}$$ The variation with respect to 
$\unicode[STIX]{x1D705}(\boldsymbol{r})$ now ensures that the density profile is equal to 
$n(\boldsymbol{r})$, and the variation with respect to 
$f_{0}$ now gives
$$\begin{eqnarray}\unicode[STIX]{x1D6FF}U=\int \unicode[STIX]{x1D6FF}f_{0}[\unicode[STIX]{x1D716}+\unicode[STIX]{x1D706}(f_{0})+\unicode[STIX]{x1D705}(\boldsymbol{r})]\sqrt{g}\,\text{d}\boldsymbol{x}=0,\end{eqnarray}$$ so that 
$\unicode[STIX]{x1D706}(f_{0})=-\unicode[STIX]{x1D716}-\unicode[STIX]{x1D705}(\boldsymbol{r})$. The ground-state distribution function must therefore be a function of the form
$$\begin{eqnarray}f_{0}(\boldsymbol{r},\boldsymbol{v})=F[\unicode[STIX]{x1D716}(\boldsymbol{r},\boldsymbol{v})+\unicode[STIX]{x1D705}(\boldsymbol{r})],\end{eqnarray}$$rather than (2.9), and the density constraint becomes
$$\begin{eqnarray}\int F[\unicode[STIX]{x1D716}(\boldsymbol{r},\boldsymbol{v})+\unicode[STIX]{x1D705}(\boldsymbol{r})]\,\text{d}\boldsymbol{v}=n(\boldsymbol{r}).\end{eqnarray}$$ If we take 
$\unicode[STIX]{x1D716}=mv^{2}/2$, we can conclude that a particular example of a ground state is an isothermal Maxwellian,
$$\begin{eqnarray}f=n(\boldsymbol{r})\left(\frac{m}{2\unicode[STIX]{x03C0}T}\right)^{3/2}\exp \left(-\frac{mv^{2}}{2T}\right),\end{eqnarray}$$ where 
$\unicode[STIX]{x1D735}T=0$. It follows that such a distribution function is linearly and nonlinearly stable to fluctuations that do not disturb the density profile. This conclusion is not surprising but has perhaps not been proved before. Moreover, it follows that any Maxwellian that is not isothermal is not in a ground state, even if it happens to be linearly stable. For instance, a plasma with 
$0<\unicode[STIX]{x1D702}=|\unicode[STIX]{x1D735}\ln T|/|\unicode[STIX]{x1D735}\ln n|<2/3$ is linearly stable to ion-temperature-gradient modes with adiabatic electrons (Hazeltine & Waelbroeck Reference Hazeltine and Waelbroeck1998) but has non-zero available energy and could thus potentially support subcritical turbulence (Newton, Cowley & Loureiro Reference Newton, Cowley and Loureiro2010; Plunk Reference Plunk2013; van Wyk et al. Reference van Wyk, Highcock, Schekochihin, Roach, Field and Dorland2016; McMillan, Pringle & Teaca Reference McMillan, Pringle and Teaca2018).
4.2 Available energy
In order to calculate the available energy of such a plasma, we write the initial distribution function as
$$\begin{eqnarray}f(\boldsymbol{x},0)=M(\boldsymbol{r})\exp \left(-\frac{mv^{2}}{2T(\boldsymbol{r})}\right),\end{eqnarray}$$ and explore the consequences of Liouville’s theorem (2.6). Since 
$f(\boldsymbol{x},0)>\unicode[STIX]{x1D719}$ if and only if
$$\begin{eqnarray}v^{2}<\frac{2T}{m}\ln \frac{M}{\unicode[STIX]{x1D719}},\end{eqnarray}$$which defines a sphere in velocity space of volume
$$\begin{eqnarray}\frac{4\unicode[STIX]{x03C0}}{3}\left(\frac{2T}{m}\ln \frac{M}{\unicode[STIX]{x1D719}}\right)^{3/2},\end{eqnarray}$$we have
$$\begin{eqnarray}H[f(\boldsymbol{x},0),\unicode[STIX]{x1D719}]=\frac{4\unicode[STIX]{x03C0}V}{3}\left\langle \left(\frac{2T}{m}\ln \frac{M}{\unicode[STIX]{x1D719}}\right)^{3/2}\unicode[STIX]{x1D6E9}(M-\unicode[STIX]{x1D719})\right\rangle\end{eqnarray}$$and
$$\begin{eqnarray}H[f_{0}(\boldsymbol{x}),\unicode[STIX]{x1D719}]=\int \unicode[STIX]{x1D6E9}\left[F\left(\unicode[STIX]{x1D716}(\boldsymbol{v})+\unicode[STIX]{x1D705}(\boldsymbol{r})\right)-\unicode[STIX]{x1D719}\right]\text{d}\boldsymbol{x}.\end{eqnarray}$$ In these equations, 
$\boldsymbol{x}=(\boldsymbol{r},\boldsymbol{v})$, 
$V$ denotes the plasma volume and angular brackets a volume average. If we denote the inverse of 
$F$ by 
$w$, so that 
$F[w(\unicode[STIX]{x1D719})]=\unicode[STIX]{x1D719}$, then 
$F(\unicode[STIX]{x1D716}(\boldsymbol{v})+\unicode[STIX]{x1D705}(\boldsymbol{r}))=\unicode[STIX]{x1D719}$ when
$$\begin{eqnarray}v=\sqrt{\frac{2}{m}\left[w(\unicode[STIX]{x1D719})-\unicode[STIX]{x1D705}(\boldsymbol{r})\right]},\end{eqnarray}$$ and we conclude that the phase-space volume in which 
$f_{0}(\boldsymbol{x})$ exceeds 
$\unicode[STIX]{x1D719}$ is equal to
$$\begin{eqnarray}H[f_{0}(\boldsymbol{x}),\unicode[STIX]{x1D719}]=\frac{4\unicode[STIX]{x03C0}V}{3}\left\langle \left[\frac{2}{m}(w(\unicode[STIX]{x1D719})-\unicode[STIX]{x1D705}(\boldsymbol{r}))\right]_{\unicode[STIX]{x1D6E9}}^{3/2}\right\rangle ,\end{eqnarray}$$where we have introduced the notation
$$\begin{eqnarray}\left[y\right]_{\unicode[STIX]{x1D6E9}}^{3/2}=y^{3/2}\unicode[STIX]{x1D6E9}(y).\end{eqnarray}$$Equation (2.6) thus becomes
$$\begin{eqnarray}\left\langle \left[w(\unicode[STIX]{x1D719})-\unicode[STIX]{x1D705}(\boldsymbol{r})\right]_{\unicode[STIX]{x1D6E9}}^{3/2}\right\rangle =\left\langle \left[T(\boldsymbol{r})\ln \frac{M(\boldsymbol{r})}{\unicode[STIX]{x1D719}}\right]_{\unicode[STIX]{x1D6E9}}^{3/2}\right\rangle ,\end{eqnarray}$$ and constitutes an integral equation for the function 
$w(\unicode[STIX]{x1D719})$ – or, alternatively, for its inverse 
$F(w)$. (The corresponding equation (3.14) in I is obtained by setting 
$\unicode[STIX]{x1D705}(\boldsymbol{r})=0$.) We thus have two nonlinear integral equations, (4.8) and (4.18), for the two unknown functions 
$\unicode[STIX]{x1D705}$ and 
$F$.
We cannot construct the general solution explicitly, but if the density and temperature of the initial state only vary slightly, it is possible to make progress. We then write
$$\begin{eqnarray}\displaystyle & \displaystyle n(\boldsymbol{r})=\left\langle n\right\rangle [1+\unicode[STIX]{x1D708}(\boldsymbol{r})], & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle T(\boldsymbol{r})=\left\langle T\right\rangle [1+\unicode[STIX]{x1D70F}(\boldsymbol{r})], & \displaystyle\end{eqnarray}$$ with 
$\unicode[STIX]{x1D708}\sim \unicode[STIX]{x1D70F}\ll 1$, so that the initial energy density is
$$\begin{eqnarray}E_{M}={\textstyle \frac{3}{2}}\left\langle nT\right\rangle ={\textstyle \frac{3}{2}}\left\langle n\right\rangle \left\langle T\right\rangle \left\langle 1+\unicode[STIX]{x1D708}\unicode[STIX]{x1D70F}\right\rangle ,\end{eqnarray}$$and seek an approximate solution of the form
$$\begin{eqnarray}F(w)=\bar{M}\exp (-w/\left\langle T\right\rangle +h(w)),\end{eqnarray}$$ where 
$\bar{M}=\left\langle n\right\rangle (m/2\unicode[STIX]{x03C0}\left\langle T\right\rangle )^{3/2}$ and 
$h\ll 1$. Before trying to solve (4.18), it is helpful to consider the density constraint (4.8), which becomes
$$\begin{eqnarray}\left\langle n\right\rangle [1+\unicode[STIX]{x1D708}(\boldsymbol{r})]=\bar{M}\exp (-\unicode[STIX]{x1D705}(\boldsymbol{r})/\left\langle T\right\rangle )\int \exp \left[-\frac{\unicode[STIX]{x1D716}(\boldsymbol{v})}{\left\langle T\right\rangle }+h(\unicode[STIX]{x1D716}(\boldsymbol{v})+\unicode[STIX]{x1D705}(\boldsymbol{r}))\right]\text{d}\boldsymbol{v}.\end{eqnarray}$$ Anticipating that 
$h$ will turn out to be of order 
$O(\unicode[STIX]{x1D708}^{2})$ and 
$\unicode[STIX]{x1D705}$ of order 
$O(\unicode[STIX]{x1D708})$, we can approximate the right-hand side by
$$\begin{eqnarray}\bar{M}\exp (-\unicode[STIX]{x1D705}(\boldsymbol{r})/\left\langle T\right\rangle )\int \exp \left[-\frac{\unicode[STIX]{x1D716}(\boldsymbol{v})}{\left\langle T\right\rangle }+h(\unicode[STIX]{x1D716}(\boldsymbol{v}))\right]\,\text{d}\boldsymbol{v},\end{eqnarray}$$ where only 
$\unicode[STIX]{x1D705}(\boldsymbol{r})$ depends on 
$\boldsymbol{r}$. It follows that
$$\begin{eqnarray}\frac{\unicode[STIX]{x1D705}}{\left\langle T\right\rangle }=-\unicode[STIX]{x1D708}+\frac{\unicode[STIX]{x1D708}^{2}}{2}+O(\unicode[STIX]{x1D708}^{3}).\end{eqnarray}$$The ground-state equation (4.18) now becomes, to the requisite accuracy,
$$\begin{eqnarray}\left\langle \left(x+\unicode[STIX]{x1D708}-\frac{\unicode[STIX]{x1D708}^{2}}{2}\right)^{3/2}\right\rangle =\left\langle (1+\unicode[STIX]{x1D70F})^{3/2}\left(x+\ln \frac{1+\unicode[STIX]{x1D708}}{(1+\unicode[STIX]{x1D70F})^{3/2}}-h\right)^{3/2}\right\rangle ,\end{eqnarray}$$ where we have written 
$x=w/\left\langle T\right\rangle$. Expanding both sides to second order in 
$\unicode[STIX]{x1D708}\sim \unicode[STIX]{x1D70F}\ll 1$ and solving for 
$h$ gives, after some algebra,
$$\begin{eqnarray}h(w)=\frac{\left\langle \unicode[STIX]{x1D70F}^{2}\right\rangle }{4}x+\frac{3}{2}\left\langle \unicode[STIX]{x1D708}\unicode[STIX]{x1D70F}-\unicode[STIX]{x1D70F}^{2}\right\rangle +\frac{3}{4x}\left\langle \frac{3\unicode[STIX]{x1D70F}^{2}}{4}-\unicode[STIX]{x1D708}\unicode[STIX]{x1D70F}\right\rangle ,\end{eqnarray}$$ confirming our earlier assumption that 
$h=O(\unicode[STIX]{x1D708}^{2})$. The ground-state distribution function is thus given by
$$\begin{eqnarray}F(\unicode[STIX]{x1D716}+\unicode[STIX]{x1D705})=\bar{M}\exp (-\unicode[STIX]{x1D716}/\left\langle T\right\rangle )\left[1+\unicode[STIX]{x1D708}+h(\unicode[STIX]{x1D716})\right],\end{eqnarray}$$and the associated energy per unit volume
$$\begin{eqnarray}E_{0}=\left\langle \int \unicode[STIX]{x1D716}F(\unicode[STIX]{x1D716}+\unicode[STIX]{x1D705})\,\text{d}\boldsymbol{v}\right\rangle =\frac{3}{2}\left\langle n\right\rangle \left\langle T\right\rangle \left\langle \unicode[STIX]{x1D708}\unicode[STIX]{x1D70F}-\frac{\unicode[STIX]{x1D70F}^{2}}{2}\right\rangle .\end{eqnarray}$$The available energy is thus
$$\begin{eqnarray}\frac{A}{V}=E_{M}-E_{0}=E_{0}\left\langle \frac{\unicode[STIX]{x1D70F}^{2}}{2}\right\rangle .\end{eqnarray}$$In contrast, if the density profile is free to evolve, the available energy becomes
$$\begin{eqnarray}\frac{A}{V}=E_{0}\left\langle \frac{\unicode[STIX]{x1D708}^{2}}{3}+\frac{\unicode[STIX]{x1D70F}^{2}}{2}\right\rangle ,\end{eqnarray}$$as originally found in I without explicitly calculating the distribution function. Thus, as one might naively expect, in a plasma with fixed density profile only the temperature variations contribute to the free energy and their contribution is the same as if the density profile were free to evolve. Lest this conclusion be thought to be self-evident, it should be noted that it does not actually hold unless the density and temperature variations are small. This fact will be illustrated by a simple example in the next section.
5 Mixing of two different plasmas
We now demonstrate that the amount of available energy from temperature variations within the plasma in general depends on the density profile even if the latter is held fixed. To this end, consider a Maxwellian plasma occupying a volume consisting of two subvolumes, 
$V=V_{1}+V_{2}$, where the density and temperature are constant in each subvolume, so that 
$n(\boldsymbol{r})=n_{j}$ and 
$T(\boldsymbol{r})=T_{j}$ for 
$\boldsymbol{r}\in V_{j}$, 
$j=1,2$. Moreover, we assume that 
$T_{2}/T_{1}\ll 1$ but 
$n_{1}V_{1}\sim n_{2}V_{2}$ and anticipate that the ground state is one where the hot plasma in 
$V_{1}$ and the cold plasma in 
$V_{2}$ have mixed in such a way that the total energy is reduced. The initial plasma energy is
$$\begin{eqnarray}E=\frac{3}{2}(n_{1}T_{1}V_{1}+n_{2}T_{2}V_{2})\simeq \frac{3n_{1}T_{1}V_{1}}{2},\end{eqnarray}$$and will be compared with the minimum attainable energy assuming either that the density is free to evolve or that it is fixed in each of the two subvolumes.
5.1 Unconstrained density profile
If the density profile is free to evolve, the equation (2.15) for the ground state becomes
$$\begin{eqnarray}\unicode[STIX]{x1D716}^{3/2}=\frac{V_{1}}{V}\left[T_{1}\ln \frac{M_{1}}{F(\unicode[STIX]{x1D716})}\right]_{\unicode[STIX]{x1D6E9}}^{3/2}+\frac{V_{2}}{V}\left[T_{2}\ln \frac{M_{2}}{F(\unicode[STIX]{x1D716})}\right]_{\unicode[STIX]{x1D6E9}}^{3/2},\end{eqnarray}$$ where 
$M_{j}=n_{j}(m/2\unicode[STIX]{x03C0}T_{j})^{3/2}$. This equation follows immediately from (3.14) in I and is similar to (4.18) with 
$\unicode[STIX]{x1D705}=0$ in the present paper. In the limit 
$T_{2}/T_{1}\rightarrow 0$, the first term on the right dominates in most of velocity space, where 
$\unicode[STIX]{x1D716}\sim T_{1}\gg T_{2}$, whereas the second term dominates for energies comparable to 
$T_{2}$. The solution is thus
$$\begin{eqnarray}F(\unicode[STIX]{x1D716})=\left\{\begin{array}{@{}ll@{}}M_{1}\exp (-\unicode[STIX]{x1D716}/\tilde{T}_{1}), & \unicode[STIX]{x1D716}\sim T_{1},\\ M_{2}\exp (-\unicode[STIX]{x1D716}/\tilde{T}_{2}), & \unicode[STIX]{x1D716}\sim T_{2},\end{array}\right.\end{eqnarray}$$ where 
$\tilde{T}_{j}=T_{j}(V_{j}/V)^{2/3}$.Footnote 3 This solution describes a bi-Maxwellian distribution such that the density of each component is constant throughout the entire plasma volume 
$V$. The energy associated with this distribution function is
$$\begin{eqnarray}E_{0}=V\int \unicode[STIX]{x1D716}F\,\text{d}\boldsymbol{v}\simeq \frac{3n_{1}T_{1}V}{2}\left(\frac{V_{1}}{V}\right)^{5/3},\end{eqnarray}$$and the available energy thus becomes
$$\begin{eqnarray}A=E-E_{0}=E\left[1-\left(\frac{V_{1}}{V}\right)^{2/3}\right].\end{eqnarray}$$ This is also the energy that is released if the hot plasma expands adiabatically (with pressure 
$p\propto V^{-5/3}$) from 
$V_{1}$ into the greater volume 
$V$.
Note that the result (5.5) does not lead to any definite prediction of what will actually happen to the plasma, only that it is possible that thermal energy can be converted into kinetic energy. It should come as no surprise that, even in the collisionless case under consideration, the amount of convertible thermal energy is bounded by classical thermodynamics.
5.2 Fixed density profile
If the density is constrained to remain equal to its initial value in each of the two subvolumes, the distribution function of the ground state no longer depends on energy alone but is a function of the form (4.7), where the function 
$\unicode[STIX]{x1D705}(\boldsymbol{r})$ enforcing the density constraint can be anticipated to be constant in each of the two subvolumes, 
$\unicode[STIX]{x1D705}(\boldsymbol{r})=\unicode[STIX]{x1D705}_{j}$ for 
$\boldsymbol{r}\in V_{j}$. The ground-state equation (4.18) then becomes
$$\begin{eqnarray}V_{1}\left[w-\unicode[STIX]{x1D705}_{1}\right]_{\unicode[STIX]{x1D6E9}}^{3/2}+V_{2}\left[w-\unicode[STIX]{x1D705}_{2}\right]_{\unicode[STIX]{x1D6E9}}^{3/2}=V_{1}\left[T_{1}\ln \frac{M_{1}}{\unicode[STIX]{x1D719}}\right]_{\unicode[STIX]{x1D6E9}}^{3/2}+V_{2}\left[T_{2}\ln \frac{M_{2}}{\unicode[STIX]{x1D719}}\right]_{\unicode[STIX]{x1D6E9}}^{3/2}.\end{eqnarray}$$ It does not seem possible to write down an explicit solution to this algebraic equation for 
$w(\unicode[STIX]{x1D719})$ except in the special case that the density vanishes in one of the two regions, 
$M_{2}\rightarrow 0$, say. The second term on the right then vanishes and 
$\unicode[STIX]{x1D705}_{2}\rightarrow \infty$, so that the solution becomes
$$\begin{eqnarray}F[\unicode[STIX]{x1D716}+\unicode[STIX]{x1D705}(\boldsymbol{r})]=\left\{\begin{array}{@{}ll@{}}M_{1}\exp (-\unicode[STIX]{x1D716}/T_{1}), & \boldsymbol{r}\in V_{1}\\ 0, & \boldsymbol{r}\in V_{2}.\end{array}\right.\end{eqnarray}$$ As one would expect from the fact that the density cannot change, the particles cannot move from one subvolume to the other, the distribution function remains equal to the initial condition, and the available energy thus vanishes, 
$A=0$. We conclude that the contribution to the available energy from temperature variations within the plasma in general depends on the density profile and on whether this profile is free to evolve or is constrained to remain constant.
6 Two species with different constraints
We finally turn to the case discussed in the Introduction of a plasma where the magnetic moment is conserved for both the ions and the electrons, but the second adiabatic invariant only for the latter. The central question is whether it is possible to confine a plasma in a Maxwellian ground state for both species with respect to perturbations that conserve 
$\unicode[STIX]{x1D707}_{i}$, 
$\unicode[STIX]{x1D707}_{e}$ and 
$J_{e}$, and, additionally, satisfy the requirement of quasineutrality. Following the method developed in §§ 2–4, possible ground states can be found by minimising the energy functional augmented by terms with Lagrange multipliers ensuring that the appropriate constraints are satisfied.
Some care is required because different constraints pertain to the trapped and passing electrons. As discussed at the beginning of § 4, the latter experience little transport when the ordering (1.5) is adopted. Since 
$\unicode[STIX]{x1D713}$ is thus conserved for each passing electron, it is advantageous to use the phase-space coordinates 
$(\unicode[STIX]{x1D707},\unicode[STIX]{x1D716},\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l)$ for these particles. Because of rapid streaming along the field lines, the distribution function 
$f_{ec}$ of circulating electrons is approximately independent of 
$\unicode[STIX]{x1D6FC}$ and 
$l$, and the local density of these particles becomes
$$\begin{eqnarray}n_{ec}=\frac{4\unicode[STIX]{x03C0}B}{m_{e}^{2}}\int _{\unicode[STIX]{x1D707}B<\unicode[STIX]{x1D716}}f_{ec}(\unicode[STIX]{x1D707},\unicode[STIX]{x1D716},\unicode[STIX]{x1D713})\frac{\text{d}\unicode[STIX]{x1D707}\,\text{d}\unicode[STIX]{x1D716}}{|v_{\Vert }|},\end{eqnarray}$$ if, for simplicity, the numbers of co- and counter-moving particles are taken to be equal. A similar expression holds for the ion density if we use analogous phase-space coordinates for the ions. For trapped electrons, however, it is more appropriate to use the coordinates 
$(\unicode[STIX]{x1D707},J,\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l)$, so that the volume element in phase space is given by (3.14) and the local density becomes
$$\begin{eqnarray}\displaystyle & & \displaystyle \hspace{-10.00002pt}n_{et}(\boldsymbol{r}_{0})=\int f_{et}\unicode[STIX]{x1D6FF}(\boldsymbol{r}-\boldsymbol{r}_{0})\,\text{d}\unicode[STIX]{x1D6E4}\nonumber\\ \displaystyle & & \displaystyle =\frac{4\unicode[STIX]{x03C0}}{m_{e}^{2}}\int f_{et}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707},J)\left|\frac{\unicode[STIX]{x2202}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l)}{\unicode[STIX]{x2202}\boldsymbol{r}}\right|\unicode[STIX]{x1D6FF}(\unicode[STIX]{x1D713}-\unicode[STIX]{x1D713}_{0})\unicode[STIX]{x1D6FF}(\unicode[STIX]{x1D6FC}-\unicode[STIX]{x1D6FC}_{0})\unicode[STIX]{x1D6FF}(l-l_{0})\frac{\text{d}l}{\unicode[STIX]{x1D70F}_{b}|v_{\Vert }|}\,\text{d}\unicode[STIX]{x1D707}\,\text{d}J\,\text{d}\unicode[STIX]{x1D713}\,\text{d}\unicode[STIX]{x1D6FC},\nonumber\\ \displaystyle & & \displaystyle\end{eqnarray}$$ where subscripts 
$0$ refer to the point 
$\boldsymbol{r}_{0}$ and 
$f_{et}$ denotes the distribution function of trapped electrons. Since the Jacobian is
$$\begin{eqnarray}\left|\frac{\unicode[STIX]{x2202}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l)}{\unicode[STIX]{x2202}\boldsymbol{r}}\right|=B,\end{eqnarray}$$the density of trapped electrons becomes
$$\begin{eqnarray}n_{et}(\boldsymbol{r}_{0})=\frac{4\unicode[STIX]{x03C0}B_{0}}{m_{e}^{2}}\int f_{et}(\unicode[STIX]{x1D713}_{0},\unicode[STIX]{x1D6FC}_{0},\unicode[STIX]{x1D707},J)\left|v_{\Vert }\unicode[STIX]{x1D70F}_{b}\right|_{0}^{-1}\,\text{d}\unicode[STIX]{x1D707}\,\text{d}J.\end{eqnarray}$$ The ambipolarity condition is thus 
$Q[f_{i},f_{ec},f_{et},\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l]=0$, where
$$\begin{eqnarray}Q[f_{i},f_{ec},f_{et},\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l]=\frac{4\unicode[STIX]{x03C0}B}{m_{i}^{2}}\int _{T\cup C}f_{i}\frac{\text{d}\unicode[STIX]{x1D707}\,\text{d}\unicode[STIX]{x1D716}}{|v_{\Vert }|}-\frac{4\unicode[STIX]{x03C0}B}{m_{e}^{2}}\left(\int _{C}f_{ec}\frac{\text{d}\unicode[STIX]{x1D707}\,\text{d}\unicode[STIX]{x1D716}}{|v_{\Vert }|}+\int _{T}f_{et}\frac{\text{d}\unicode[STIX]{x1D707}\,\text{d}J}{\left|v_{\Vert }\unicode[STIX]{x1D70F}_{b}\right|}\right).\end{eqnarray}$$ Here, 
$C$ refers to the circulating region of velocity space,
$$\begin{eqnarray}0\leqslant \unicode[STIX]{x1D707}<\frac{\unicode[STIX]{x1D716}}{B_{\max }},\end{eqnarray}$$
$T$ to the trapped region,
$$\begin{eqnarray}\frac{\unicode[STIX]{x1D716}}{B_{\max }}\leqslant \unicode[STIX]{x1D707}\leqslant \frac{\unicode[STIX]{x1D716}}{B},\end{eqnarray}$$ and 
$B_{\max }(\unicode[STIX]{x1D713})$ denotes the maximum magnetic-field strength on the surface labelled by 
$\unicode[STIX]{x1D713}$. Similarly, the energy functional is
$$\begin{eqnarray}\displaystyle E[f_{i},f_{ec},f_{et}] & = & \displaystyle 4\unicode[STIX]{x03C0}\int \text{d}\unicode[STIX]{x1D713}\,\text{d}\unicode[STIX]{x1D6FC}\,\text{d}l\left(\int _{T\cup C}\unicode[STIX]{x1D716}f_{i}\frac{\text{d}\unicode[STIX]{x1D707}\,\text{d}\unicode[STIX]{x1D716}}{m_{i}^{2}|v_{\Vert }|}\right.\nonumber\\ \displaystyle & & \displaystyle \left.+\int _{C}\unicode[STIX]{x1D716}f_{ec}\frac{\text{d}\unicode[STIX]{x1D707}\,\text{d}\unicode[STIX]{x1D716}}{m_{e}^{2}|v_{\Vert }|}+\int _{T}\unicode[STIX]{x1D716}f_{et}\frac{\text{d}\unicode[STIX]{x1D707}\,\text{d}J}{m_{e}^{2}\left|v_{\Vert }\unicode[STIX]{x1D70F}_{b}\right|}\right),\end{eqnarray}$$and possible ground states can be found be seeking the minimum of the functional
$$\begin{eqnarray}\displaystyle & & \displaystyle W_{1}[f_{e},f_{i},\unicode[STIX]{x1D706}_{e},\unicode[STIX]{x1D706}_{i},\unicode[STIX]{x1D708}]=E[f_{i},f_{ec},f_{et}]+\int \unicode[STIX]{x1D708}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l)Q[f_{i},f_{ec},f_{et},\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l]\,\text{d}\unicode[STIX]{x1D713}\,\text{d}\unicode[STIX]{x1D6FC}\,\text{d}l\nonumber\\ \displaystyle & & \displaystyle \quad +\,\frac{4\unicode[STIX]{x03C0}}{m_{i}^{2}}\int _{T\cup C}\unicode[STIX]{x1D706}_{i}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707})\,\text{d}\unicode[STIX]{x1D719}\,\text{d}\unicode[STIX]{x1D707}\int [\unicode[STIX]{x1D6E9}(f_{i}-\unicode[STIX]{x1D719})-H_{i}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707})]\frac{\text{d}\unicode[STIX]{x1D716}}{|v_{\Vert }|}\,\text{d}\unicode[STIX]{x1D713}\,\text{d}\unicode[STIX]{x1D6FC}\,\text{d}l\nonumber\\ \displaystyle & & \displaystyle \quad +\,\frac{4\unicode[STIX]{x03C0}}{m_{e}^{2}}\int _{C}\unicode[STIX]{x1D706}_{ec}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},\unicode[STIX]{x1D713})\,\text{d}\unicode[STIX]{x1D719}\,\text{d}\unicode[STIX]{x1D707}\,\text{d}\unicode[STIX]{x1D713}\int [\unicode[STIX]{x1D6E9}(f_{ec}-\unicode[STIX]{x1D719})-H_{ec}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},\unicode[STIX]{x1D713})]\frac{\text{d}\unicode[STIX]{x1D716}}{|v_{\Vert }|}\,\text{d}\unicode[STIX]{x1D6FC}\,\text{d}l\nonumber\\ \displaystyle & & \displaystyle \quad +\,\frac{4\unicode[STIX]{x03C0}}{m_{e}^{2}}\int _{T}\unicode[STIX]{x1D706}_{et}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},J)\,\text{d}\unicode[STIX]{x1D719}\,\text{d}\unicode[STIX]{x1D707}\,\text{d}J\int [\unicode[STIX]{x1D6E9}(f_{et}-\unicode[STIX]{x1D719})-H_{et}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},J)]\text{d}\unicode[STIX]{x1D713}\,\text{d}\unicode[STIX]{x1D6FC}.\end{eqnarray}$$ In this expression, the first term on the right represents the total energy of the ions and electrons, and the next term, which contains the Lagrange multiplier 
$\unicode[STIX]{x1D708}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC},l)$, enforces quasineutrality on each field line 
$(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC})$. The term containing 
$\unicode[STIX]{x1D706}_{i}(\unicode[STIX]{x1D719})$ ensures the conservation of 
$\unicode[STIX]{x1D707}$ for all ions, that containing 
$\unicode[STIX]{x1D706}_{ec}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},\unicode[STIX]{x1D713})$ makes sure that 
$\unicode[STIX]{x1D707}$ and 
$\unicode[STIX]{x1D713}$ are conserved for circulating electrons and the term with 
$\unicode[STIX]{x1D706}_{et}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707},J)$ enforces conservation of 
$\unicode[STIX]{x1D707}$ and 
$J$ for trapped electrons.
The variation with respect to 
$f_{i}$ gives
$$\begin{eqnarray}\int \unicode[STIX]{x1D6FF}f_{i}\left[\unicode[STIX]{x1D716}+\unicode[STIX]{x1D708}+\int \unicode[STIX]{x1D706}_{i}(\unicode[STIX]{x1D719},\unicode[STIX]{x1D707})\unicode[STIX]{x1D6FF}(f_{i}-\unicode[STIX]{x1D719})\,\text{d}\unicode[STIX]{x1D719}\right]\frac{\text{d}\unicode[STIX]{x1D707}_{i}\,\text{d}\unicode[STIX]{x1D716}}{|v_{\Vert }|}\,\text{d}\unicode[STIX]{x1D713}\,\text{d}\unicode[STIX]{x1D6FC}\,\text{d}l=0,\end{eqnarray}$$which implies
$$\begin{eqnarray}\unicode[STIX]{x1D716}+\unicode[STIX]{x1D708}+\unicode[STIX]{x1D706}_{i}(f_{i},\unicode[STIX]{x1D707})=0.\end{eqnarray}$$ It must therefore be possible to write 
$f_{i}$ as a function of the form
$$\begin{eqnarray}f_{i}=F_{i}[\unicode[STIX]{x1D716}+\unicode[STIX]{x1D708}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC}),\unicode[STIX]{x1D707}],\end{eqnarray}$$ where we have noted that rapid streaming along 
$\boldsymbol{B}$ implies that the distribution function must be independent of 
$l$, and so must therefore the function 
$\unicode[STIX]{x1D708}$. Similarly, the variation of 
$W_{1}$ with respect to the electron distribution function shows that the latter can be written as
$$\begin{eqnarray}\displaystyle & \displaystyle f_{ec}=F_{ec}[\unicode[STIX]{x1D716}-\unicode[STIX]{x1D708}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC}),\unicode[STIX]{x1D707},\unicode[STIX]{x1D713}], & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle f_{et}=F_{et}[\unicode[STIX]{x1D716}-\unicode[STIX]{x1D708}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC}),\unicode[STIX]{x1D707},J]. & \displaystyle\end{eqnarray}$$6.1 Maxwellian ground states
Any possible ground states in plasmas with conservation properties derived from (1.5) must be of the form indicated by (6.12)–(6.14), but the most important ones are of course those where the distribution functions of both species are Maxwellian. Equation (6.12) implies that this is only possible if the ions are isothermal, so that
$$\begin{eqnarray}f_{i}=n(\unicode[STIX]{x1D713})\left(\frac{m_{i}}{2\unicode[STIX]{x03C0}T_{i}}\right)^{3/2}\exp (-\unicode[STIX]{x1D716}/T_{i}),\end{eqnarray}$$ where 
$T_{i}$ is constant and we have taken 
$\unicode[STIX]{x1D708}(\unicode[STIX]{x1D713},\unicode[STIX]{x1D6FC})=-T_{i}\ln n(\unicode[STIX]{x1D713})$. For the electrons, however, the temperature may vary across magnetic flux surfaces,
$$\begin{eqnarray}f_{e}=n(\unicode[STIX]{x1D713})\left(\frac{m_{e}}{2\unicode[STIX]{x03C0}T_{e}(\unicode[STIX]{x1D713})}\right)^{3/2}\exp (-\unicode[STIX]{x1D716}/T_{e}(\unicode[STIX]{x1D713})).\end{eqnarray}$$ In order that 
$f_{e}$ be of the form (6.13) and (6.14), we write 
$\unicode[STIX]{x1D705}=\unicode[STIX]{x1D716}-\unicode[STIX]{x1D708}(\unicode[STIX]{x1D713})$ and define
$$\begin{eqnarray}\displaystyle & \hspace{-5.0pt}\displaystyle F_{ec}(\unicode[STIX]{x1D705},\unicode[STIX]{x1D707},\unicode[STIX]{x1D713})=N(\unicode[STIX]{x1D713})\left(\frac{m_{e}}{2\unicode[STIX]{x03C0}T_{e}(\unicode[STIX]{x1D713})}\right)^{3/2}\exp (-\unicode[STIX]{x1D705}/T_{e}(\unicode[STIX]{x1D713})),\qquad & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \hspace{-5.0pt}\displaystyle F_{et}(\unicode[STIX]{x1D705},\unicode[STIX]{x1D707},J)=N[\unicode[STIX]{x1D713}(\unicode[STIX]{x1D705},\unicode[STIX]{x1D707},J)]\left(\frac{m_{e}}{2\unicode[STIX]{x03C0}T_{e}[\unicode[STIX]{x1D713}(\unicode[STIX]{x1D705},\unicode[STIX]{x1D707},J)]}\right)^{3/2}\exp (-\unicode[STIX]{x1D705}/T_{e}[\unicode[STIX]{x1D713}(\unicode[STIX]{x1D705},\unicode[STIX]{x1D707},J)]),\qquad & \displaystyle\end{eqnarray}$$ where the flux function 
$N(\unicode[STIX]{x1D713})$ is related to the density by
$$\begin{eqnarray}n(\unicode[STIX]{x1D713})=N(\unicode[STIX]{x1D713})\exp (\unicode[STIX]{x1D708}(\unicode[STIX]{x1D713})/T_{e}(\unicode[STIX]{x1D713})).\end{eqnarray}$$ This prescription suggests that it could be possible to confine a plasma in a Maxwellian ground state for both species with respect to perturbations that satisfy the ordering (1.5) if the ions are isothermal. However, we must remember that (6.12)–(6.14) only provide necessary, and not sufficient, conditions on any ground states. In the case of a single species without conserved quantities, the corresponding condition is given by (2.9), but a ground state must also satisfy 
$\text{d}F/\text{d}\unicode[STIX]{x1D716}\leqslant 0$ everywhere. To find the corresponding condition for the present two-species case, let us again consider what happens to the energy of the system when the particles in two neighbouring and equally large volume elements in phase space are interchanged, but we now allow such an interchange to occur both in the population of ions and in that of trapped electrons.
As in § 2, the energy released is equal to
$$\begin{eqnarray}\displaystyle & \displaystyle \text{d}E_{i}=-\text{d}\unicode[STIX]{x1D716}_{i}\,\text{d}f_{i}\,\text{d}\unicode[STIX]{x1D6E4}_{i}, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \text{d}E_{e}=-\text{d}\unicode[STIX]{x1D716}_{e}\,\text{d}f_{e}\,\text{d}\unicode[STIX]{x1D6E4}_{e}, & \displaystyle\end{eqnarray}$$ where 
$\text{d}\unicode[STIX]{x1D6E4}_{a}$ denotes the phase-space volume of the elements, 
$\text{d}\unicode[STIX]{x1D716}_{a}$ their difference in energy and 
$\text{d}f_{a}$ the difference in the distribution functions. In order that the interchanges be ambipolar, we must require
$$\begin{eqnarray}\text{d}\unicode[STIX]{x1D713}_{i}\,\text{d}f_{i}\,\text{d}\unicode[STIX]{x1D6E4}_{i}=\text{d}\unicode[STIX]{x1D713}_{e}\,\text{d}f_{e}\,\text{d}\unicode[STIX]{x1D6E4}_{e},\end{eqnarray}$$and conclude that the total amount of energy released is
$$\begin{eqnarray}\text{d}E=\text{d}E_{i}+\text{d}E_{e}=-\text{d}\unicode[STIX]{x1D716}_{e}\text{d}f_{e}\text{d}\unicode[STIX]{x1D6E4}_{e}\left(1+\frac{\text{d}\unicode[STIX]{x1D716}_{i}\,\text{d}\unicode[STIX]{x1D713}_{e}}{\text{d}\unicode[STIX]{x1D716}_{e}\,\text{d}\unicode[STIX]{x1D713}_{i}}\right).\end{eqnarray}$$ For the trapped electrons, 
$\text{d}\unicode[STIX]{x1D713}_{e}$ and 
$\text{d}\unicode[STIX]{x1D716}_{e}$ are related by (3.8),
$$\begin{eqnarray}\text{d}\unicode[STIX]{x1D716}_{e}=-e\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}\,\text{d}\unicode[STIX]{x1D713}_{e},\end{eqnarray}$$ which has to do with the conservation of 
$J_{e}$: it is in general impossible to change the energy of a trapped electron without also moving it radially. We also have
$$\begin{eqnarray}\text{d}f_{e}=\left(\frac{\unicode[STIX]{x2202}f_{et}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D716}}\right)_{\unicode[STIX]{x1D707},J}\,\text{d}\unicode[STIX]{x1D716}_{e}=\frac{f_{et}}{T_{e}}\left(\frac{\unicode[STIX]{x1D714}_{\ast e}^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}}-1\right)\text{d}\unicode[STIX]{x1D716}_{e},\end{eqnarray}$$and thus
$$\begin{eqnarray}\text{d}E=-\frac{f_{et}}{T_{e}}\left(\frac{\unicode[STIX]{x1D714}_{\ast e}^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}}-1\right)\left(\text{d}\unicode[STIX]{x1D716}_{e}\right)^{2}\left(1-\frac{1}{e\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}}\frac{\text{d}\unicode[STIX]{x1D716}_{i}}{\text{d}\unicode[STIX]{x1D713}_{i}}\right).\end{eqnarray}$$ The first term in the last bracket describes the energy change 
$\text{d}E_{e}$ in the electron channel and second one that in the ions. As found in I, the former is always positive if the condition (3.5) is satisfied,
$$\begin{eqnarray}\frac{\unicode[STIX]{x1D714}_{\ast e}^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}}-1<0,\end{eqnarray}$$ but the total energy may nevertheless decrease. If more energy is released by the ions than taken up the electrons, so that 
$\text{d}E<0$, the system could not have been in a ground state before the interchange. This is evidently the case if 
$\text{d}\unicode[STIX]{x1D716}_{i}$ and 
$\text{d}\unicode[STIX]{x1D713}_{i}$ are chosen in such a way that the condition
$$\begin{eqnarray}\frac{1}{e\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}}\frac{\text{d}\unicode[STIX]{x1D716}_{i}}{\text{d}\unicode[STIX]{x1D713}_{i}}>1\end{eqnarray}$$is satisfied. Whether this is possible could, in principle, depend on the density gradient and the magnetic geometry. Since
$$\begin{eqnarray}\text{d}f_{i}=\left(\frac{\unicode[STIX]{x2202}f_{i}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D716}}\right)_{\unicode[STIX]{x1D713}}\,\text{d}\unicode[STIX]{x1D716}_{i}+\left(\frac{\unicode[STIX]{x2202}f_{i}}{\unicode[STIX]{x2202}\unicode[STIX]{x1D713}}\right)_{\unicode[STIX]{x1D716}}\,\text{d}\unicode[STIX]{x1D713}_{i}=\frac{f_{i}}{T_{i}}\left(e\unicode[STIX]{x1D714}_{\ast i}\,\text{d}\unicode[STIX]{x1D713}_{i}-\text{d}\unicode[STIX]{x1D716}_{i}\right),\end{eqnarray}$$the ambipolarity condition (6.22) can be written as
$$\begin{eqnarray}\frac{f_{i}}{T_{i}}\left(e\unicode[STIX]{x1D714}_{\ast i}\,\text{d}\unicode[STIX]{x1D713}_{i}-\text{d}\unicode[STIX]{x1D716}_{i}\right)\text{d}\unicode[STIX]{x1D713}_{i}\,\text{d}\unicode[STIX]{x1D6E4}_{i}=\frac{ef_{e}}{T_{e}}\left(\frac{\unicode[STIX]{x1D714}_{\ast e}^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}}-1\right)\text{d}\unicode[STIX]{x1D716}_{e}\,\text{d}\unicode[STIX]{x1D713}_{e}\,\text{d}\unicode[STIX]{x1D6E4}_{e},\end{eqnarray}$$and implies
$$\begin{eqnarray}\frac{1}{e\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}}\frac{\text{d}\unicode[STIX]{x1D716}_{i}}{\text{d}\unicode[STIX]{x1D713}_{i}}=\frac{\unicode[STIX]{x1D714}_{\ast i}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}}+\left(\frac{\unicode[STIX]{x1D714}_{\ast e}^{T}}{\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}}-1\right)\left(\frac{\text{d}\unicode[STIX]{x1D713}_{e}}{\text{d}\unicode[STIX]{x1D713}_{i}}\right)^{2}\frac{T_{i}f_{e}\,\text{d}\unicode[STIX]{x1D6E4}_{e}}{T_{e}f_{i}\,\text{d}\unicode[STIX]{x1D6E4}_{i}}.\end{eqnarray}$$ The question is whether this expression can be made to exceed unity by any choice of the coordinates and volumes of the phase-space volume elements that are to be interchanged. If so, the plasma is not in a ground state according to (6.28). Since 
$\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}e}$ is proportional to energy, the two terms on the right that have this quantity in the denominator can be made arbitrarily large by choosing the energy to be small.Footnote 4 And since 
$\unicode[STIX]{x1D714}_{\ast e}$ and 
$\unicode[STIX]{x1D714}_{\ast i}$ have opposite signs, it is clear that (6.31) will always exceed unity for some choice of parameters. We thus conclude that non-trivial ground states do not exist in the present situation. In other words, it does not seem possible to confine a two-component plasma in a Maxwellian minimum-energy state with respect to fluctuations satisfying (1.5), unless the density and temperature gradients vanish. If the criterion (3.5) is satisfied for the electrons, their energy cannot be lowered, but it is apparently always possible to redistribute ions in such a way that the total energy decreases and quasineutrality is maintained.
Finally, let us put this conclusion in the context of microinstabilities in stellarators and dipoles. As noted in I, the criterion (3.5) for the ground state of a single species can be met in two different ways. It holds if either the average magnetic curvature is favourable for all trapped particles and the temperature gradient is moderate,
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D714}_{\ast }\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}<0, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle 0<\unicode[STIX]{x1D702}<{\textstyle \frac{2}{3}}, & \displaystyle\end{eqnarray}$$or, alternatively, if the curvature is unfavourable and the temperature gradient large enough, i.e.
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D714}_{\ast }\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}>0, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D702}>{\textstyle \frac{2}{3}}. & \displaystyle\end{eqnarray}$$ The first one of these cases holds in a maximum-
$J$ stellarator and implies that ordinary density-gradient-driven trapped-electron modes are linearly stable. Proll et al. (Reference Proll, Helander, Connor and Plunk2012) and Helander, Proll & Plunk (Reference Helander, Proll and Plunk2013) concluded that any collisionless instability must then draw energy from the ions rather than the electrons, and Plunk, Connor & Helander (Reference Plunk, Connor and Helander2017) showed that there is indeed a remnant ion-driven trapped-electron mode if 
$\unicode[STIX]{x1D702}_{i}>0$. From our considerations of the available energy, this is not surprising, since it is apparently impossible to achieve a non-trivial ground state for both species simultaneously in the ordering (1.5). On the other hand, if the electrons are in a ground state, they act to stabilise any ion-driven instability and one would expect that the available energy is relatively small. Gyrokinetic simulations of turbulence in approximately maximum-
$J$ stellarator configurations do indeed suggest that the transport under these conditions is up to an order of magnitude smaller than in a typical tokamak (Helander et al. Reference Helander, Bird, Jenko, Kleiber, Plunk, Proll, Riemann and Xanthopoulos2015).
As discussed in I, the second case (6.34) can apply in the field of a magnetic dipole, which is then nonlinearly stable with respect to instabilities that conserve 
$\unicode[STIX]{x1D707}$ and 
$J$ for all species. This is of particular importance for electron–positron plasmas, where identical constraints apply to both species (Helander Reference Helander2014a). In electron–ion dipole plasmas, however, no such statement can be made concerning instabilities that do not conserve 
$J_{i}$, but low-frequency modes can be remarkably stable (Hasegawa, Chen & Mauel Reference Hasegawa, Chen and Mauel1990; Kesner & Hastie Reference Kesner and Hastie2002; Simakov, Hastie & Catto Reference Simakov, Hastie and Catto2002).
7 Conclusions
The available energy of a plasma – the amount of thermal energy that can be converted into linear and nonlinear fluctuations – depends on what constraints limit the possible forms of plasma motion. In the present paper, this concept has been explored in the context of magnetically confined plasmas in local thermodynamic equilibrium. Specifically, an explicit formula (3.32) has been derived for the energy available from a species for which the first and second adiabatic invariants, 
$\unicode[STIX]{x1D707}$ and 
$J$, are conserved. This formula shows that the available energy is equal to a weighted average of 
$\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}^{2}(\unicode[STIX]{x1D714}_{\ast }^{T}/\unicode[STIX]{x1D714}_{\unicode[STIX]{x1D6FC}}-1)$ over the subset of phase space where this quantity is positive.
Frequently, different constraints apply to the electrons and the ions, but their combined motion must be such that quasineutrality is preserved. This condition adds another constraint and thus lowers the available energy. It also changes the functional form of the distribution function for the ‘ground state’ of lowest accessible energy. In a plasma with fixed density profile, this distribution function no longer depends on energy alone, but on energy plus a function of position ensuring ambipolarity, see (4.7). Since the density profile cannot relax, there is no available energy associated with the density gradient, and in a Vlasov plasma with small fluctuations this energy is given by (4.30) rather than (4.31), which applies in a plasma where the density is free to evolve. Thus, small density variations do not affect the available energy of a plasma with fixed density profile. However, if the density variations are not small, then they can affect the amount of energy that can be extracted from temperature variations in the plasma, as demonstrated in § 5.
Finally, we have noted that, for some of the most important forms of turbulence in magnetically confined fusion plasmas, the magnetic moment 
$\unicode[STIX]{x1D707}$ is conserved for the ions whereas both 
$\unicode[STIX]{x1D707}$ and the parallel invariant 
$J$ are conserved for trapped electrons, while 
$\unicode[STIX]{x1D707}$ and 
$\unicode[STIX]{x1D713}$ are conserved for circulating electrons. This is, for instance, usually the case for ion-temperature-gradient-driven modes and trapped-electron modes, which are believed to cause much of the turbulence and transport observed in fusion experiments. In the collisionless, electrostatic limit, Proll et al. (Reference Proll, Helander, Connor and Plunk2012) and Helander et al. (Reference Helander, Proll and Plunk2013) derived the criterion (3.5) for the stability of trapped-electron modes taking energy from the electrons in the plasma. This condition came out of lengthy gyrokinetic calculations but can be simply understood in terms of the available energy. When this criterion is satisfied, any instability satisfying the orderings must derive its energy from the ions in the plasma, and Plunk et al. (Reference Plunk, Connor and Helander2017) identified a particular such instability. The considerations of available energy given in § 6 suggest that it is indeed always possible to lower the total energy of a Maxwellian plasma with non-zero density or temperature gradients.
Before closing, we note that all our results are, strictly speaking, only valid in zero-
$\unicode[STIX]{x1D6FD}$, collisionless plasmas. Collisions violate Liouville’s theorem and can enable the plasma to access lower-energy states than otherwise possible. Similarly, allowing for magnetic fluctuations can fundamentally change the stability properties of a plasma. However, experience from gyrokinetic simulations suggests that a small amount of collisions or plasma pressure does not substantially affect the turbulence. Insofar that the properties of the latter are correlated with the available energy, one can thus hope that this remains true beyond the strict realm of applicability of the calculations presented here.
Appendix A. Available potential energy of the atmosphere
In this appendix A, we establish the relation between our formalism and the results of Lorenz (Reference Lorenz1955). Consider an atmosphere consisting of an ideal gas satisfying the equations
$$\begin{eqnarray}\displaystyle & \displaystyle \frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}{\unicode[STIX]{x2202}t}+\unicode[STIX]{x1D735}\boldsymbol{\cdot }(\unicode[STIX]{x1D70C}\boldsymbol{v})=0, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \frac{\unicode[STIX]{x2202}s}{\unicode[STIX]{x2202}t}+\boldsymbol{v}\boldsymbol{\cdot }\unicode[STIX]{x1D735}s=0, & \displaystyle\end{eqnarray}$$ where 
$\unicode[STIX]{x1D70C}$ denotes the density, 
$\boldsymbol{v}$ the flow velocity and 
$s=p/\unicode[STIX]{x1D70C}^{\unicode[STIX]{x1D6FE}}$ the specific entropy. Here, 
$p$ is the pressure and 
$\unicode[STIX]{x1D6FE}$ the adiabatic index (ratio of specific heats). The first of these equations expresses conservation of mass and the second that of entropy. Any function 
$f(s)$ of 
$s$ alone then satisfies the conservation law
$$\begin{eqnarray}\frac{\text{d}}{\text{d}t}\int \unicode[STIX]{x1D70C}f(s)\,\text{d}V=0,\end{eqnarray}$$ where the volume integral is taken over the entire atmosphere and the normal component of 
$\unicode[STIX]{x1D70C}\boldsymbol{v}$ is assumed to vanish on the boundary. Specifically, we may choose 
$f$ to be a Heaviside function and conclude that the quantity
$$\begin{eqnarray}H(\unicode[STIX]{x1D70E})=\int \unicode[STIX]{x1D70C}\unicode[STIX]{x1D6E9}(s-\unicode[STIX]{x1D70E})\,\text{d}V\end{eqnarray}$$ does not change with time for any value of 
$\unicode[STIX]{x1D70E}\in [0,\infty )$.
To find the ground state, we follow the recipe from § 2 and minimise the energy functional
$$\begin{eqnarray}E=\int \left(\frac{p}{\unicode[STIX]{x1D6FE}-1}+\unicode[STIX]{x1D70C}\unicode[STIX]{x1D719}\right)\text{d}V,\end{eqnarray}$$ where 
$\unicode[STIX]{x1D719}=gz$ denotes the gravitational potential, subject to the constraint that 
$H(\unicode[STIX]{x1D70E})$ is fixed. We are thus led to consider the functional
$$\begin{eqnarray}W[\unicode[STIX]{x1D70C},s,\unicode[STIX]{x1D706}]=\int \left(u+\unicode[STIX]{x1D70C}\unicode[STIX]{x1D719}\right)\text{d}V+\int _{0}^{\infty }\unicode[STIX]{x1D706}(\unicode[STIX]{x1D70E})\,\text{d}\unicode[STIX]{x1D70E}\left(\int \unicode[STIX]{x1D70C}\unicode[STIX]{x1D6E9}(s-\unicode[STIX]{x1D70E})\,\text{d}V-H(\unicode[STIX]{x1D70E})\right),\end{eqnarray}$$where
$$\begin{eqnarray}u(\unicode[STIX]{x1D70C},s)=\frac{s\unicode[STIX]{x1D70C}^{\unicode[STIX]{x1D6FE}}}{\unicode[STIX]{x1D6FE}-1}\end{eqnarray}$$ and the Lagrange multipliers 
$\unicode[STIX]{x1D706}(\unicode[STIX]{x1D70E})$ ensure that 
$H(\unicode[STIX]{x1D70E})$ is indeed conserved for each 
$\unicode[STIX]{x1D70E}\geqslant 0$. The variations of 
$W$ with respect to 
$\unicode[STIX]{x1D70C}$ and 
$s$ give
$$\begin{eqnarray}\displaystyle & \displaystyle \unicode[STIX]{x1D719}+\frac{\unicode[STIX]{x2202}u}{\unicode[STIX]{x2202}\unicode[STIX]{x1D70C}}+\int _{0}^{\infty }\unicode[STIX]{x1D706}(\unicode[STIX]{x1D70E})\unicode[STIX]{x1D6E9}(s-\unicode[STIX]{x1D70E})\,\text{d}\unicode[STIX]{x1D70E}=0, & \displaystyle\end{eqnarray}$$
$$\begin{eqnarray}\displaystyle & \displaystyle \frac{\unicode[STIX]{x2202}u}{\unicode[STIX]{x2202}s}+\unicode[STIX]{x1D706}(s)\unicode[STIX]{x1D70C}=0, & \displaystyle\end{eqnarray}$$respectively. The second of these equations implies
$$\begin{eqnarray}\unicode[STIX]{x1D706}(s)=-\frac{\unicode[STIX]{x1D70C}^{\unicode[STIX]{x1D6FE}-1}}{\unicode[STIX]{x1D6FE}-1},\end{eqnarray}$$which, when substituted into the gradient of the first equation,
$$\begin{eqnarray}\unicode[STIX]{x1D735}\unicode[STIX]{x1D719}+\frac{\unicode[STIX]{x1D6FE}}{\unicode[STIX]{x1D6FE}-1}\unicode[STIX]{x1D735}(s\unicode[STIX]{x1D70C}^{\unicode[STIX]{x1D6FE}-1})+\unicode[STIX]{x1D706}(s)\unicode[STIX]{x1D735}s=0,\end{eqnarray}$$gives
$$\begin{eqnarray}\unicode[STIX]{x1D70C}\unicode[STIX]{x1D735}\unicode[STIX]{x1D719}+\unicode[STIX]{x1D735}p=0.\end{eqnarray}$$ The lowest-energy state is thus one in which the gas is in mechanical equilibrium. Moreover, according to (A 10), in this equilibrium 
$p$ is a function of 
$s$ alone, so that the pressure varies with height in the same way at every location on Earth. In other words, as noted by Lorenz, the minimum total energy which can result from adiabatic rearrangement of air occurs when the pressure is everywhere equal to its horizontal average.
Open access
