2 The general KMHD system and its conservation laws

The KMHD theory describes a collisionless, quasineutral plasma with a single ion species (of unit charge for simplicity), in the strong magnetization limit of zero Larmor radii where particles are perfectly tied to the magnetic field lines. The kinetic pressures of each species are formally taken as comparable to the magnetic pressure (so the skin depths are of the order of the Larmor radii) and the macroscopic flow velocities are formally taken as comparable to the sound speed. Thus, the Hall and other two-fluid effects become negligible in the leading order as the Larmor radii tend to zero, in which limit the two species share a common, single-fluid macroscopic velocity $\boldsymbol{u}$ . The quasineutrality condition implies that they share also a common particle density $n$ . The result is a hybrid model, whereby a system of single-fluid, dissipation-free, hydromagnetic equations is closed with a pressure tensor which is the sum of ion and electron contributions determined kinetically:

(2.1) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}\boldsymbol{B}}{\unicode[STIX]{x2202}t}=\unicode[STIX]{x1D735}\times (\boldsymbol{u}\times \boldsymbol{B})\end{eqnarray}$$

(2.2) $$\begin{eqnarray}\boldsymbol{j}=\unicode[STIX]{x1D735}\times \boldsymbol{B}\end{eqnarray}$$

(2.3) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}n}{\unicode[STIX]{x2202}t}+\unicode[STIX]{x1D735}\boldsymbol{\cdot }(n\boldsymbol{u})=0\end{eqnarray}$$

(2.4) $$\begin{eqnarray}\unicode[STIX]{x1D70C}\left[\frac{\unicode[STIX]{x2202}\boldsymbol{u}}{\unicode[STIX]{x2202}t}+(\boldsymbol{u}\boldsymbol{\cdot }\unicode[STIX]{x1D735})\boldsymbol{u}\right]-\boldsymbol{j}\times \boldsymbol{B}+\mathop{\sum }_{s=i,e}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D617}_{s}=0,\end{eqnarray}$$

where $\unicode[STIX]{x1D70C}=(m_{i}+m_{e})n$ denotes the mass density and the ion and electron pressure tensors have the zero-Larmor-radius collisionless form (Chew et al. Reference Chew, Goldberger and Low1956)

(2.5) $$\begin{eqnarray}\unicode[STIX]{x1D617}_{s}=p_{s\bot }\unicode[STIX]{x1D644}+(p_{s\Vert }-p_{s\bot })\boldsymbol{b}\boldsymbol{b},\end{eqnarray}$$

with $\boldsymbol{b}=\boldsymbol{B}/B$ denoting the magnetic unit vector and $\unicode[STIX]{x1D644}$ denoting the identity tensor. The parallel and perpendicular pressures are the corresponding moments of the distribution functions of each species,

(2.6a,b ) $$\begin{eqnarray}p_{s\Vert }=m_{s}\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }^{2}f_{s},\quad p_{s\bot }=\frac{m_{s}}{2}\int \text{d}^{3}\boldsymbol{w}~w_{\bot }^{2}~f_{s}\end{eqnarray}$$

and the phase-space velocity variable $\boldsymbol{w}$ is the random velocity in the reference frame of the macroscopic flow, $\boldsymbol{w}=\boldsymbol{v}-\boldsymbol{u}$ . The distribution functions are independent of the gyrophase angle and satisfy a zero-Larmor-radius drift-kinetic equation. It is most advantageous to refer this drift-kinetic equation to the shifted velocity variable $\boldsymbol{w}$ so that its solution, $f_{s}=f_{s}(w_{\Vert },w_{\bot },\boldsymbol{x},t)$ , can be used directly in the evaluation of the pressure moments (2.6), without having to subtract any components of the macroscopic velocity. The transformation to this reference frame of the macroscopic flow allows also to eliminate exactly the electric field from the drift-kinetic equation, which becomes (Ramos Reference Ramos2008)

(2.7) $$\begin{eqnarray}\displaystyle & & \displaystyle \frac{\unicode[STIX]{x2202}f_{s}}{\unicode[STIX]{x2202}t}+\left(\boldsymbol{u}+w_{\Vert }\boldsymbol{b}\right)\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}f_{s}}{\unicode[STIX]{x2202}\boldsymbol{x}}\Bigg|_{w_{\Vert },w_{\bot }}+\frac{w_{\bot }}{2}\left[(\boldsymbol{b}\boldsymbol{b}-\unicode[STIX]{x1D644})\boldsymbol{ : }(\unicode[STIX]{x1D735}\boldsymbol{u})+w_{\Vert }\boldsymbol{b}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\ln B\right]\frac{\unicode[STIX]{x2202}f_{s}}{\unicode[STIX]{x2202}w_{\bot }}\nonumber\\ \displaystyle & & \displaystyle \quad +\left[\frac{\boldsymbol{b}\boldsymbol{\cdot }\left(\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D617}_{s}\right)}{m_{s}n}-w_{\Vert }(\boldsymbol{b}\boldsymbol{b})\boldsymbol{ : }(\unicode[STIX]{x1D735}\boldsymbol{u})-\frac{w_{\bot }^{2}}{2}\boldsymbol{b}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\ln B\right]\frac{\unicode[STIX]{x2202}f_{s}}{\unicode[STIX]{x2202}w_{\Vert }}=0.\end{eqnarray}$$

The closed system of (2.1)–(2.7) (which has been written without any explicit reference to the electric field but has taken it properly into account) is the formulation of the KMHD theoretical model that will be investigated in this work. This general, nonlinear system has several exact conservation properties. The continuity equation (2.3) is in the form of conservation law for the number of particles and the equation of motion (2.4) combined with (2.3) can be expressed as the conservation law for momentum,

(2.8) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}(\unicode[STIX]{x1D70C}\boldsymbol{u})}{\unicode[STIX]{x2202}t}+\unicode[STIX]{x1D735}\boldsymbol{\cdot }\left(\unicode[STIX]{x1D70C}\boldsymbol{u}\boldsymbol{u}+\frac{1}{2}B^{2}\unicode[STIX]{x1D644}-\boldsymbol{B}\boldsymbol{B}+\mathop{\sum }_{s=i,e}\unicode[STIX]{x1D617}_{s}\right)=0.\end{eqnarray}$$

The drift-kinetic equation (2.7) upholds these two conservation laws. For consistency with particle conservation and quasineutrality, the ‘kinetically defined’ particle densities $n_{s}^{kin}\equiv \int \text{d}^{3}\boldsymbol{w}~f_{s}$ should always be the same for both species and equal to the solution $n$ of the fluid continuity equation (2.3). Since the momentum conservation equation was subtracted from the kinetic equation when transforming it to the macroscopic flow frame and eliminating the electric field, the $w_{\Vert }$ moment of (2.7) is not redundant with the parallel component of (2.4). Instead, the $w_{\Vert }$ moment of $f_{s}$ should by definition be always zero. Now, calling $c_{s\Vert }\equiv \int \text{d}^{3}\boldsymbol{w}~w_{\Vert }f_{s}$ , the $1$ and $w_{\Vert }$ moments of the drift-kinetic equation (2.7) combined with the continuity equation (2.3) yield

(2.9) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}(n_{s}^{kin}-n)}{\unicode[STIX]{x2202}t}+\unicode[STIX]{x1D735}\boldsymbol{\cdot }\left[(n_{s}^{kin}-n)\boldsymbol{u}+c_{s\Vert }\boldsymbol{b}\right]=0,\end{eqnarray}$$

(2.10) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}c_{s\Vert }}{\unicode[STIX]{x2202}t}+\unicode[STIX]{x1D735}\boldsymbol{\cdot }(c_{s\Vert }\boldsymbol{u})+c_{s\Vert }(\boldsymbol{b}\boldsymbol{b})\boldsymbol{ : }(\unicode[STIX]{x1D735}\boldsymbol{u})-\frac{(n_{s}^{kin}-n)}{m_{s}n}~\boldsymbol{b}\boldsymbol{\cdot }(\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D617}_{s})=0.\end{eqnarray}$$

This is a homogeneous system for $(n_{s}^{kin}-n)$ and $c_{s\Vert }$ , which guarantees that the required constraints $n_{s}^{kin}-n=0$ and $c_{s\Vert }=0$ are satisfied at all times provided they are satisfied by the initial condition. Finally, the KMHD system has an energy conservation law. Combining equation (2.1) dotted with $\boldsymbol{B}$ , equation (2.4) dotted with $\boldsymbol{u}$ , equation (2.3) and the $w^{2}$ moment of equation (2.7), one gets

(2.11) $$\begin{eqnarray}\displaystyle & & \displaystyle \hspace{-10.00002pt}\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}t}\left[\frac{\unicode[STIX]{x1D70C}u^{2}}{2}+\frac{B^{2}}{2}+\frac{1}{2}\mathop{\sum }_{s=i,e}\left(p_{s\Vert }+2p_{s\bot }\right)\right]\nonumber\\ \displaystyle & & \displaystyle \quad +\,\unicode[STIX]{x1D735}\boldsymbol{\cdot }\left\{\frac{\unicode[STIX]{x1D70C}u^{2}}{2}\boldsymbol{u}-(\boldsymbol{u}\times \boldsymbol{B})\times \boldsymbol{B}+\mathop{\sum }_{s=i,e}\left[\frac{1}{2}\left(p_{s\Vert }+2p_{s\bot }\right)\boldsymbol{u}+\unicode[STIX]{x1D617}_{s}\boldsymbol{\cdot }\boldsymbol{u}+q_{s\Vert }\boldsymbol{b}\right]\right\}=0,\nonumber\\ \displaystyle & & \displaystyle\end{eqnarray}$$

where $q_{s\Vert }=m_{s}/2\int \text{d}^{3}\boldsymbol{w}w_{\Vert }w^{2}f_{s}$ is the parallel heat flux. In order to streamline the exposition, this paper will assume that the plasma domain is bounded by a rigid and perfectly conducting wall where the normal components of $\boldsymbol{u}$ and $\boldsymbol{B}$ vanish. As indicated by the normal-mode analysis of Ramos (Reference Ramos2015b ), the results can be generalized to the case where the plasma domain is surrounded by a vacuum region and the plasma–vacuum interface is a closed magnetic surface where the density drops continuously to zero, but the details of this generalization will be skipped here. Thus, for the purposes of the present paper, the global version of the energy conservation law is

(2.12) $$\begin{eqnarray}\frac{\text{d}}{\text{d}t}\int \text{d}^{3}\boldsymbol{x}\left[\frac{\unicode[STIX]{x1D70C}u^{2}}{2}+\frac{B^{2}}{2}+\frac{1}{2}\mathop{\sum }_{s=i,e}(p_{s\Vert }+2p_{s\bot })\right]=0\end{eqnarray}$$

and the space integral $\int \text{d}^{3}\boldsymbol{x}$ extends over the plasma domain.

3 Linearized KMHD system and specification of initial conditions

The linear stability problem that this work will address concerns small-amplitude perturbations about a static equilibrium, and the subscripts $0$ and $1$ will be used to denote respectively equilibrium and first-order perturbation variables. The equilibrium distribution functions will be assumed here to be isotropic Maxwellians,

(3.1) $$\begin{eqnarray}f_{Ms0}=\left(\frac{m_{s}}{2\unicode[STIX]{x1D70B}}\right)^{3/2}\frac{n_{0}}{T_{s0}^{3/2}}~\exp \left(-\frac{m_{s}w^{2}}{2T_{s0}}\right),\end{eqnarray}$$

with density and temperatures that are uniform in the direction of the equilibrium magnetic field,

(3.2a,b ) $$\begin{eqnarray}\boldsymbol{B}_{0}\boldsymbol{\cdot }\unicode[STIX]{x1D735}n_{0}=0,\quad \boldsymbol{B}_{0}\boldsymbol{\cdot }\unicode[STIX]{x1D735}T_{s0}=0\end{eqnarray}$$

and satisfy the force balance equation

(3.3) $$\begin{eqnarray}\boldsymbol{j}_{0}\times \boldsymbol{B}_{0}=(\unicode[STIX]{x1D735}\times \boldsymbol{B}_{0})\times \boldsymbol{B}_{0}=\unicode[STIX]{x1D735}[n_{0}(T_{i0}+T_{e0})]\end{eqnarray}$$

so that, with a vanishing macroscopic flow, $\boldsymbol{u}_{0}=0$ , this constitutes a time-independent solution of the KMHD system (2.1)–(2.7). The Maxwellian equilibrium distribution function is a preferred choice on physical grounds, but the analysis is equally applicable to other equilibrium distribution functions that are constant along the magnetic field, isotropic in velocity space and monotonically decreasing with $w^{2}$ .

The first-order perturbation satisfies the linearized KMHD system

(3.4) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}\boldsymbol{B}_{1}}{\unicode[STIX]{x2202}t}=\unicode[STIX]{x1D735}\times (\boldsymbol{u}_{1}\times \boldsymbol{B}_{0}),\end{eqnarray}$$

(3.5) $$\begin{eqnarray}\boldsymbol{j}_{1}=\unicode[STIX]{x1D735}\times \boldsymbol{B}_{1},\end{eqnarray}$$

(3.6) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}n_{1}}{\unicode[STIX]{x2202}t}+\unicode[STIX]{x1D735}\boldsymbol{\cdot }(n_{0}\boldsymbol{u}_{1})=0,\end{eqnarray}$$

(3.7) $$\begin{eqnarray}\unicode[STIX]{x1D70C}_{0}\frac{\unicode[STIX]{x2202}\boldsymbol{u}_{1}}{\unicode[STIX]{x2202}t}-\boldsymbol{j}_{0}\times \boldsymbol{B}_{1}-\boldsymbol{j}_{1}\times \boldsymbol{B}_{0}+\mathop{\sum }_{s=i,e}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D617}_{s1}=0,\end{eqnarray}$$

(3.8) $$\begin{eqnarray}\unicode[STIX]{x1D617}_{s1}=p_{s\bot 1}\unicode[STIX]{x1D644}+(p_{s\Vert 1}-p_{s\bot 1})\boldsymbol{b}_{0}\boldsymbol{b}_{0},\end{eqnarray}$$

(3.9a,b ) $$\begin{eqnarray}p_{s\Vert 1}=m_{s}\int \text{d}^{3}\boldsymbol{w}w_{\Vert }^{2}f_{s1},\quad p_{s\bot 1}=\frac{m_{s}}{2}\int \text{d}^{3}\boldsymbol{w}~w_{\bot }^{2}f_{s1},\end{eqnarray}$$

(3.10) $$\begin{eqnarray}\displaystyle & & \displaystyle \frac{\unicode[STIX]{x2202}f_{s1}}{\unicode[STIX]{x2202}t}+w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}f_{s1}}{\unicode[STIX]{x2202}\boldsymbol{x}}\Bigg|_{w_{\Vert },w_{\bot }}+~\frac{w_{\bot }}{2}(\boldsymbol{b}_{0}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\ln B_{0})\left(w_{\Vert }\frac{\unicode[STIX]{x2202}f_{s1}}{\unicode[STIX]{x2202}w_{\bot }}-w_{\bot }\frac{\unicode[STIX]{x2202}f_{s1}}{\unicode[STIX]{x2202}w_{\Vert }}\right)\nonumber\\ \displaystyle & & \displaystyle \qquad +\left(\boldsymbol{u}_{1}+w_{\Vert }\boldsymbol{b}_{1}\right)\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}f_{Ms0}}{\unicode[STIX]{x2202}\boldsymbol{x}}\Bigg|_{w_{\Vert },w_{\bot }}+~\frac{w_{\bot }}{2}[(\boldsymbol{b}_{0}\boldsymbol{b}_{0}-\unicode[STIX]{x1D644})\boldsymbol{ : }(\unicode[STIX]{x1D735}\boldsymbol{u}_{1})]\frac{\unicode[STIX]{x2202}f_{Ms0}}{\unicode[STIX]{x2202}w_{\bot }}\nonumber\\ \displaystyle & & \displaystyle \qquad +\left[\frac{\boldsymbol{b}_{0}\boldsymbol{\cdot }(\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D617}_{s1})+\boldsymbol{b}_{1}\boldsymbol{\cdot }(\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D617}_{s0})}{m_{s}n_{0}}-w_{\Vert }(\boldsymbol{b}_{0}\boldsymbol{b}_{0})\boldsymbol{ : }(\unicode[STIX]{x1D735}\boldsymbol{u}_{1})\right]\frac{\unicode[STIX]{x2202}f_{Ms0}}{\unicode[STIX]{x2202}w_{\Vert }}=0.\end{eqnarray}$$

The initial-value solution of this system requires specification of the initial conditions $\boldsymbol{B}_{1}(0)$ , $n_{1}(0)$ , $\boldsymbol{u}_{1}(0)$ and $f_{s1}(0)$ , subject to the constraints $\unicode[STIX]{x1D735}\boldsymbol{\cdot }\boldsymbol{B}_{1}(0)=0$ , $\int \text{d}^{3}\boldsymbol{w}~f_{s1}(0)=n_{1}(0)$ and $\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }f_{s1}(0)=0$ . The zero-divergence constraint on $\boldsymbol{B}_{1}(0)$ leaves it with only two independent components, so the initial conditions for $\boldsymbol{B}_{1}$ and $n_{1}$ can be implemented in terms of a three-component vector $\unicode[STIX]{x1D743}(0)=\unicode[STIX]{x1D709}_{\Vert }(0)\boldsymbol{b}_{0}+\unicode[STIX]{x1D743}_{\bot }(0)$ , through $\boldsymbol{B}_{1}(0)=\unicode[STIX]{x1D735}\times [\unicode[STIX]{x1D743}_{\bot }(0)\times \boldsymbol{B}_{0}]$ and $n_{1}(0)=-\unicode[STIX]{x1D743}_{\bot }(0)\boldsymbol{\cdot }\unicode[STIX]{x1D735}n_{0}-n_{0}~\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743}(0)$ . Now, defining the time-dependent displacement vector

(3.11) $$\begin{eqnarray}\unicode[STIX]{x1D743}(t)=\unicode[STIX]{x1D743}(0)+\int _{0}^{t}\,\text{d}t^{\prime }\,\boldsymbol{u}_{1}(t^{\prime }),\end{eqnarray}$$

equations (3.4) and (3.6) can be integrated to get the relationships

(3.12) $$\begin{eqnarray}\boldsymbol{B}_{1}=\unicode[STIX]{x1D735}\times (\unicode[STIX]{x1D743}_{\bot }\times \boldsymbol{B}_{0}),\end{eqnarray}$$

(3.13) $$\begin{eqnarray}n_{1}=-\unicode[STIX]{x1D743}_{\bot }\boldsymbol{\cdot }\unicode[STIX]{x1D735}n_{0}-n_{0}~\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743},\end{eqnarray}$$

valid at all times, which allow us to eliminate $\boldsymbol{B}_{1}$ , $\boldsymbol{j}_{1}$ and $n_{1}$ in favour of $\unicode[STIX]{x1D743}$ in the dynamical system. The dynamical displacement vector $\unicode[STIX]{x1D743}$ also serves to define the convection-subtracted (Lagrangian) perturbations of the distribution functions,

(3.14) $$\begin{eqnarray}\hat{f}_{s}=f_{s1}+\unicode[STIX]{x1D743}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}f_{Ms0}}{\unicode[STIX]{x2202}\boldsymbol{x}},\end{eqnarray}$$

so the pressure tensors become

(3.15) $$\begin{eqnarray}\unicode[STIX]{x1D617}_{s1}=-[\unicode[STIX]{x1D743}_{\bot }\boldsymbol{\cdot }\unicode[STIX]{x1D735}(n_{0}T_{s0})]\unicode[STIX]{x1D644}+\hat{\unicode[STIX]{x1D617}}_{s}\end{eqnarray}$$

with

(3.16) $$\begin{eqnarray}\hat{\unicode[STIX]{x1D617}}_{s}=\left(\frac{m_{s}}{2}\int \text{d}^{3}\boldsymbol{w}\,w_{\bot }^{2}~\hat{f}_{s}\right)\unicode[STIX]{x1D644}+\left[m_{s}\int \text{d}^{3}\boldsymbol{w}\left(w_{\Vert }^{2}-\frac{w_{\bot }^{2}}{2}\right)\hat{f}_{s}\right]\boldsymbol{b}_{0}\boldsymbol{b}_{0},\end{eqnarray}$$

and the parallel and perpendicular components of $-\unicode[STIX]{x1D735}\boldsymbol{\cdot }\hat{\unicode[STIX]{x1D617}}_{s}$ are the following functionals of $\hat{f}_{s}$ :

(3.17) $$\begin{eqnarray}\hat{F}_{s\Vert }[\hat{f}_{s}]=-\boldsymbol{b}_{0}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\left[m_{s}\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }^{2}\hat{f}_{s}\right]+\left[m_{s}\int \text{d}^{3}\boldsymbol{w}\,\left(w_{\Vert }^{2}-\frac{w_{\bot }^{2}}{2}\right)\hat{f}_{s}\right](\boldsymbol{b}_{0}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\ln B_{0})\end{eqnarray}$$

and

(3.18) $$\begin{eqnarray}\hat{\boldsymbol{F}}_{s\bot }[\hat{f}_{s}]=-\unicode[STIX]{x1D735}_{\bot }\left[\frac{m_{s}}{2}\int \text{d}^{3}\boldsymbol{w}~w_{\bot }^{2}\hat{f}_{s}\right]-\left[m_{s}\int \text{d}^{3}\boldsymbol{w}\,\left(w_{\Vert }^{2}-\frac{w_{\bot }^{2}}{2}\right)\hat{f}_{s}\right]\unicode[STIX]{x1D73F}_{0},\end{eqnarray}$$

where $\unicode[STIX]{x1D73F}_{0}=(\boldsymbol{b}_{0}\boldsymbol{\cdot }\unicode[STIX]{x1D735})\boldsymbol{b}_{0}$ is the equilibrium magnetic curvature. Then, recalling that the force operator in ‘perpendicular ideal MHD’ is the functional of $\unicode[STIX]{x1D743}_{\bot }$

(3.19) $$\begin{eqnarray}\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]=(\boldsymbol{j}_{0}\times \boldsymbol{B}_{1})_{\bot }+\boldsymbol{j}_{1}\times \boldsymbol{B}_{0}+\unicode[STIX]{x1D735}_{\bot }[\unicode[STIX]{x1D743}_{\bot }\boldsymbol{\cdot }\unicode[STIX]{x1D735}(n_{0}T_{i0}+n_{0}T_{e0})],\end{eqnarray}$$

with $\boldsymbol{B}_{1}$ and $\boldsymbol{j}_{1}$ given respectively by (3.12) and (3.5), and that, as a consequence of (3.3) and (3.12) the parallel component of $(\boldsymbol{j}_{0}\times \boldsymbol{B}_{1})+\unicode[STIX]{x1D735}[\unicode[STIX]{x1D743}_{\bot }\boldsymbol{\cdot }\unicode[STIX]{x1D735}(n_{0}T_{i0}+n_{0}T_{e0})]$ vanishes, the linearized momentum equation (3.7) can be written as

(3.20) $$\begin{eqnarray}\unicode[STIX]{x1D70C}_{0}\frac{\unicode[STIX]{x2202}\boldsymbol{u}_{1}}{\unicode[STIX]{x2202}t}=\unicode[STIX]{x1D70C}_{0}\frac{\unicode[STIX]{x2202}^{2}\unicode[STIX]{x1D743}}{\unicode[STIX]{x2202}t^{2}}=\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]+\mathop{\sum }_{s=i,e}\left(\hat{F}_{s\Vert }[\hat{f}_{s}]\boldsymbol{b}_{0}+\hat{\boldsymbol{F}}_{s\bot }[\hat{f}_{s}]\right).\end{eqnarray}$$

The transformation of the linearized drift-kinetic equation (3.10) to the variable $\hat{f}_{s}$ yields

(3.21) $$\begin{eqnarray}\displaystyle & & \displaystyle \frac{\unicode[STIX]{x2202}\hat{f}_{s}}{\unicode[STIX]{x2202}t}+w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\hat{f}_{s}}{\unicode[STIX]{x2202}\boldsymbol{x}}\Bigg|_{w_{\Vert },w_{\bot }}+\frac{w_{\bot }}{2}(\boldsymbol{b}_{0}\boldsymbol{\cdot }\unicode[STIX]{x1D735}\ln B_{0})\left(w_{\Vert }\frac{\unicode[STIX]{x2202}\hat{f}_{s}}{\unicode[STIX]{x2202}w_{\bot }}-w_{\bot }\frac{\unicode[STIX]{x2202}\hat{f}_{s}}{\unicode[STIX]{x2202}w_{\Vert }}\right)\nonumber\\ \displaystyle & & \displaystyle \quad +\left[\frac{w_{\Vert }\hat{F}_{s\Vert }}{m_{s}n_{0}}+\frac{w_{\bot }^{2}}{2}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\boldsymbol{u}_{1}+\left(w_{\Vert }^{2}-\frac{w_{\bot }^{2}}{2}\right)(\boldsymbol{b}_{0}\boldsymbol{b}_{0})\boldsymbol{ : }(\unicode[STIX]{x1D735}\boldsymbol{u}_{1})\right]\frac{m_{s}f_{Ms0}}{T_{s0}}=0,\quad\end{eqnarray}$$

subject to the initial-value constraints $\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}(0)=-n_{0}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743}(0)$ and $\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }\hat{f}_{s}(0)=0$ . Now, equations (3.17)–(3.21) constitute the sought-after formulation of the linearized KMHD system that uses $(\unicode[STIX]{x1D743},\hat{f}_{s})$ as the primary dynamical variables.

It is convenient to express $\hat{f}_{s}$ as the sum $\hat{f}_{s}^{even}+\hat{f}_{s}^{odd}$ of its even and odd parts with respect to $w_{\Vert }$ and split equation (3.21) into a system for $\hat{f}_{s}^{even}$ and $\hat{f}_{s}^{odd}$ . Also, it is useful to change the phase-space coordinates from $(w_{\Vert },w_{\bot })$ to the coordinates $w=(w_{\Vert }^{2}+w_{\bot }^{2})^{1/2}$ and $\unicode[STIX]{x1D706}=w_{\bot }^{2}/(w^{2}B_{0})$ related to the kinetic energy and the magnetic moment. Denoting the phase-space parallel derivative at constant $w$ and $\unicode[STIX]{x1D706}$ simply by $\boldsymbol{b}_{0}\boldsymbol{\cdot }\unicode[STIX]{x2202}/\unicode[STIX]{x2202}\boldsymbol{x}$ and defining the $Q$ operator

(3.22) $$\begin{eqnarray}Q[\boldsymbol{\unicode[STIX]{x1D702}}]\equiv \frac{1}{w^{2}}\left[\frac{w_{\bot }^{2}}{2}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\boldsymbol{\unicode[STIX]{x1D702}}+\left(w_{\Vert }^{2}-\frac{w_{\bot }^{2}}{2}\right)(\boldsymbol{b}_{0}\boldsymbol{b}_{0})\boldsymbol{ : }(\unicode[STIX]{x1D735}\boldsymbol{\unicode[STIX]{x1D702}})\right],\end{eqnarray}$$

one obtains

(3.23) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{even}}{\unicode[STIX]{x2202}t}+w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{odd}}{\unicode[STIX]{x2202}\boldsymbol{x}}+Q[\boldsymbol{u}_{1}]~\frac{m_{s}w^{2}f_{Ms0}}{T_{s0}}=0,\end{eqnarray}$$

(3.24) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{odd}}{\unicode[STIX]{x2202}t}+w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{even}}{\unicode[STIX]{x2202}\boldsymbol{x}}+w_{\Vert }\hat{F}_{s\Vert }~\frac{f_{Ms0}}{n_{0}T_{s0}}=0.\end{eqnarray}$$

Several useful identities involving the parallel derivative $\boldsymbol{b}_{0}\boldsymbol{\cdot }\unicode[STIX]{x2202}/\unicode[STIX]{x2202}\boldsymbol{x}$ at constant $w$ and $\unicode[STIX]{x1D706}$ are listed in appendix A. By virtue of (A 4), the functional $\hat{F}_{s\Vert }[\hat{f}_{s}]$ (3.17) reduces to

(3.25) $$\begin{eqnarray}\hat{F}_{s\Vert }[\hat{f}_{s}]=-m_{s}\int \text{d}^{3}\boldsymbol{w}w_{\Vert }^{2}\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{even}}{\unicode[STIX]{x2202}\boldsymbol{x}}.\end{eqnarray}$$

Then, the $w_{\Vert }$ moment of (3.24) yields

(3.26) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}t}\left(\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }\hat{f}_{s}^{odd}\right)=0\end{eqnarray}$$

and, using (A 3), the $1$ moment of (3.23) yields

(3.27) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}t}\left[\left(\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}^{even}\right)+n_{0}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743}\right]+\unicode[STIX]{x1D735}\boldsymbol{\cdot }\left[\left(\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }\hat{f}_{s}^{odd}\right)\boldsymbol{b}_{0}\right]=0.\end{eqnarray}$$

Thus, as expected from the general nonlinear result (2.9) and (2.10), the linearized drift-kinetic equation (3.21) also preserves the constraints $\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}=-n_{0}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743}$ and $\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }\hat{f}_{s}=0$ at all times, once they are imposed on the initial condition.

4 The potential energy functional

The first term of the nonlinearly conserved energy (2.12) is the kinetic energy of the macroscopic flow, $K=\int \text{d}^{3}\boldsymbol{x}~\unicode[STIX]{x1D70C}u^{2}/2$ . One may call the remainder a ‘potential energy $W$ ’ and state the energy conservation law as $\text{d}(K+W)/\text{d}t=0$ . Within the context of small-amplitude perturbations about a static equilibrium, it can be shown that such $W$ agrees with the physical definition of potential energy as the integral along a dynamic trajectory of the force times the infinitesimal displacement, the result of which is a functional of the final state independent of the details of the trajectory leading to it. For the considered small-amplitude perturbations about a static equilibrium, the kinetic energy is quadratic in the perturbation. Therefore, retaining quantities to the second perturbative order, one has

(4.1) $$\begin{eqnarray}K=K_{2}=\int \text{d}^{3}\boldsymbol{x}~\frac{\unicode[STIX]{x1D70C}_{0}u_{1}^{2}}{2}=\int \text{d}^{3}\boldsymbol{x}~\frac{\unicode[STIX]{x1D70C}_{0}}{2}~\Big|\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D743}}{\unicode[STIX]{x2202}t}\Big|^{2}\end{eqnarray}$$

and $W=W_{0}+W_{1}+W_{2}$ , where

(4.2) $$\begin{eqnarray}\displaystyle & \displaystyle W_{0}=\int \text{d}^{3}\boldsymbol{x}\left(\frac{B_{0}^{2}}{2}+\frac{3}{2}\mathop{\sum }_{s=i,e}n_{0}T_{s0}\right), & \displaystyle\end{eqnarray}$$

(4.3) $$\begin{eqnarray}\displaystyle & \displaystyle W_{1}=\int \text{d}^{3}\boldsymbol{x}\left[\boldsymbol{B}_{0}\boldsymbol{\cdot }\boldsymbol{B}_{1}+\frac{1}{2}\mathop{\sum }_{s=i,e}(p_{s\Vert 1}+2p_{s\bot 1})\right], & \displaystyle\end{eqnarray}$$

(4.4) $$\begin{eqnarray}\displaystyle & \displaystyle W_{2}=\int \text{d}^{3}\boldsymbol{x}\left[\frac{B_{1}^{2}}{2}+\boldsymbol{B}_{0}\boldsymbol{\cdot }\boldsymbol{B}_{2}+\frac{1}{2}\mathop{\sum }_{s=i,e}(p_{s\Vert 2}+2p_{s\bot 2})\right]. & \displaystyle\end{eqnarray}$$

Since $K$ is a second-order quantity, energy conservation implies that $W_{1}$ must be independent of time like $W_{0}$ . This is also a consequence of the equilibrium condition (3.3), as shown in appendix B. Moreover, instead of using (4.4) that involves the second-order perturbation of the magnetic field and the pressures, $W_{2}$ can be inferred from the energy conservation law. Thus, up to a constant $C$ independent of time, one can write $W(t)=W_{0}+W_{1}+W_{2}(t)=C+\unicode[STIX]{x1D6FF}W(t)$ , with

(4.5) $$\begin{eqnarray}\frac{\text{d}\unicode[STIX]{x1D6FF}W}{\text{d}t}=-\frac{\text{d}K}{\text{d}t}=-\int \text{d}^{3}\boldsymbol{x}~\unicode[STIX]{x1D70C}_{0}\boldsymbol{u}_{1}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\boldsymbol{u}_{1}}{\unicode[STIX]{x2202}t}.\end{eqnarray}$$

The linearized momentum equation (3.20) implies that $\unicode[STIX]{x1D6FF}W$ is the sum of a fluid term plus a kinetic term, $\unicode[STIX]{x1D6FF}W=\unicode[STIX]{x1D6FF}W_{\bot }^{F}+\unicode[STIX]{x1D6FF}W^{K}$ , given by

(4.6) $$\begin{eqnarray}\frac{\text{d}\unicode[STIX]{x1D6FF}\,W_{\bot }^{F}}{\text{d}t}=-\int \text{d}^{3}\boldsymbol{x}\,\boldsymbol{u}_{1\bot }\boldsymbol{\cdot }\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]\end{eqnarray}$$

and

(4.7) $$\begin{eqnarray}\frac{\text{d}\unicode[STIX]{x1D6FF}\,W^{K}}{\text{d}t}=-\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\,\left(u_{1\Vert }\hat{F}_{s\Vert }[\hat{f}_{s}]+\boldsymbol{u}_{1\bot }\boldsymbol{\cdot }\hat{F}_{s\bot }[\hat{f}_{s}]\right).\end{eqnarray}$$

The ‘perpendicular ideal-MHD’ force operator $\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]$ is known to be self-adjoint (Bernstein et al. Reference Bernstein, Frieman, Kruskal and Kurlsrud1957, Reference Bernstein, Frieman, Kruskal and Kurlsrud1958), therefore

(4.8) $$\begin{eqnarray}\int \text{d}^{3}\boldsymbol{x}\boldsymbol{u}_{1\bot }\boldsymbol{\cdot }\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]=\int \text{d}^{3}\boldsymbol{x}~\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D743}_{\bot }}{\unicode[STIX]{x2202}t}\boldsymbol{\cdot }\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]=\frac{\text{d}}{\text{d}t}\left(\frac{1}{2}\int \text{d}^{3}\boldsymbol{x}\,\unicode[STIX]{x1D743}_{\bot }\boldsymbol{\cdot }\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]\right)\end{eqnarray}$$

and (4.6) can be integrated to express $\unicode[STIX]{x1D6FF}W_{\bot }^{F}$ as the functional of $\unicode[STIX]{x1D743}_{\bot }$

(4.9) $$\begin{eqnarray}\unicode[STIX]{x1D6FF}W_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]=-\frac{1}{2}\int \text{d}^{3}\boldsymbol{x}\,\unicode[STIX]{x1D743}_{\bot }\boldsymbol{\cdot }\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }].\end{eqnarray}$$

It will be shown next that, for physical trajectories that are solutions of the linearized KMHD system, equation (4.7) can also be integrated to get an expression of $\unicode[STIX]{x1D6FF}W^{K}$ as a functional of $\hat{f}_{s}$ only. Substituting (3.25) and (3.18) for $\hat{F}_{s\Vert }$ and $\hat{F}_{s\bot }$ , integrating by parts and using the identity (A 3), one obtains

(4.10) $$\begin{eqnarray}\displaystyle & & \displaystyle \hspace{-10.00002pt}\frac{\text{d}\unicode[STIX]{x1D6FF}\,W^{K}}{\text{d}t}=\mathop{\sum }_{s=i,e}m_{s}\int \text{d}^{3}\boldsymbol{x}~\unicode[STIX]{x1D735}\boldsymbol{\cdot }\left[\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}^{even}\left(w_{\Vert }^{2}u_{1\Vert }\boldsymbol{b}_{0}+\frac{w_{\bot }^{2}}{2}\boldsymbol{u}_{1\bot }\right)\right]\nonumber\\ \displaystyle & & \displaystyle \hspace{-10.00002pt}\quad -\mathop{\sum }_{s=i,e}m_{s}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,\hat{f}_{s}^{even}\left[w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}(w_{\Vert }u_{1\Vert })}{\unicode[STIX]{x2202}\boldsymbol{x}}+\frac{w_{\bot }^{2}}{2}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\boldsymbol{u}_{1\bot }-\left(w_{\Vert }^{2}-\frac{w_{\bot }^{2}}{2}\right)\boldsymbol{u}_{1\bot }\boldsymbol{\cdot }\unicode[STIX]{x1D73F}_{0}\right].\nonumber\\ \displaystyle & & \displaystyle\end{eqnarray}$$

The first term of this expression vanishes by virtue of Green’s divergence theorem and the ideal wall boundary conditions, and the second term can be rewritten using the identity (A 5) and the definition (3.22). The result is

(4.11) $$\begin{eqnarray}\frac{\text{d}\unicode[STIX]{x1D6FF}\,W^{K}}{\text{d}t}=-\mathop{\sum }_{s=i,e}m_{s}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,\hat{f}_{s}^{even}w^{2}Q[\boldsymbol{u}_{1}]\end{eqnarray}$$

and, substituting the drift-kinetic evolution equation for $\hat{f}_{s}^{even}$ (3.23),

(4.12) $$\begin{eqnarray}\frac{\text{d}\unicode[STIX]{x1D6FF}\,W^{K}}{\text{d}t}=\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\frac{T_{s0}}{f_{Ms0}}\hat{f}_{s}^{even}\left[\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{even}}{\unicode[STIX]{x2202}t}+w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{odd}}{\unicode[STIX]{x2202}\boldsymbol{x}}\right].\end{eqnarray}$$

Further integration by parts with a vanishing boundary term yields

(4.13) $$\begin{eqnarray}\frac{\text{d}\unicode[STIX]{x1D6FF}\,W^{K}}{\text{d}t}=\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\frac{T_{s0}}{f_{Ms0}}\left[\frac{1}{2}\frac{\unicode[STIX]{x2202}(\hat{f}_{s}^{even})^{2}}{\unicode[STIX]{x2202}t}-w_{\Vert }\hat{f}_{s}^{odd}~\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{even}}{\unicode[STIX]{x2202}\boldsymbol{x}}\right]\end{eqnarray}$$

and, substituting the drift-kinetic equation for $\hat{f}_{s}^{odd}$  (3.24),

(4.14) $$\begin{eqnarray}\frac{\text{d}\unicode[STIX]{x1D6FF}\,W^{K}}{\text{d}t}=\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\left\{\frac{T_{s0}}{2f_{Ms0}}\left[\frac{\unicode[STIX]{x2202}(\hat{f}_{s}^{even})^{2}}{\unicode[STIX]{x2202}t}+\frac{\unicode[STIX]{x2202}(\hat{f}_{s}^{odd})^{2}}{\unicode[STIX]{x2202}t}\right]+\frac{\hat{F}_{s\Vert }w_{\Vert }\hat{f}_{s}^{odd}}{n_{0}}\right\}.\end{eqnarray}$$

Finally, recalling $\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }\hat{f}_{s}^{odd}=0$ , this becomes

(4.15) $$\begin{eqnarray}\frac{\text{d}\unicode[STIX]{x1D6FF}\,W^{K}}{\text{d}t}=\frac{\text{d}}{\text{d}t}\left\{\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\frac{T_{s0}}{2f_{Ms0}}\left[(\hat{f}_{s}^{even})^{2}+(\hat{f}_{s}^{odd})^{2}\right]\right\},\end{eqnarray}$$

which can be integrated to arrive at the desired result

(4.16) $$\begin{eqnarray}\unicode[STIX]{x1D6FF}W^{K}[\hat{f}_{s}]=\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\frac{T_{s0}}{2f_{Ms0}}\left[(\hat{f}_{s}^{even})^{2}+(\hat{f}_{s}^{odd})^{2}\right]=\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\frac{T_{s0}}{2f_{Ms0}}\hat{f}_{s}^{2}.\end{eqnarray}$$

In summary, adding its fluid and kinetic terms, the KMHD incremental potential energy functional is

(4.17) $$\begin{eqnarray}\unicode[STIX]{x1D6FF}W[\unicode[STIX]{x1D743},\hat{f}_{s}]=\unicode[STIX]{x1D6FF}W_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]+\unicode[STIX]{x1D6FF}W^{K}[\hat{f}_{s}]=-\frac{1}{2}\int \text{d}^{3}\boldsymbol{x}\,\unicode[STIX]{x1D743}_{\bot }\boldsymbol{\cdot }\boldsymbol{F}_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }]+\frac{1}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,\frac{T_{s0}}{f_{Ms0}}\hat{f}_{s}^{2},\end{eqnarray}$$

which is equivalent to the expression advanced in (1.1) since, for the equilibrium Maxwellian distribution functions, $\unicode[STIX]{x2202}f_{Ms0}/\unicode[STIX]{x2202}\unicode[STIX]{x1D700}=-f_{Ms0}/T_{s0}$ . The second-order energy conservation law becomes

(4.18) $$\begin{eqnarray}\frac{\text{d}}{\text{d}t}(K[\unicode[STIX]{x2202}\unicode[STIX]{x1D743}/\unicode[STIX]{x2202}t]+\unicode[STIX]{x1D6FF}W[\unicode[STIX]{x1D743},\hat{f}_{s}])=0,\end{eqnarray}$$

with the $K$ and $\unicode[STIX]{x1D6FF}W$ functionals given respectively by (4.1) and (4.17). The relationship between this energy conservation law and the quadratic form for normal modes derived in Ramos (Reference Ramos2015b ) is discussed in appendix C.

The criterion that a positive definite $\unicode[STIX]{x1D6FF}W[\unicode[STIX]{x1D743},\hat{f}_{s}]$ is sufficient for KMHD stability is an immediate consequence of the second-order energy conservation law because then, for any dynamical evolution,

(4.19) $$\begin{eqnarray}K(t)=K(0)+\unicode[STIX]{x1D6FF}W(0)-\unicode[STIX]{x1D6FF}W(t)\leqslant K(0)+\unicode[STIX]{x1D6FF}W(0)\end{eqnarray}$$

and $K(t)$ would be a bounded function of time. The Cauchy–Schwarz inequality and the particle conservation constraint yield

(4.20) $$\begin{eqnarray}\int \text{d}^{3}\boldsymbol{w}~\frac{T_{s0}}{f_{Ms0}}\hat{f}_{s}^{2}\geqslant \left(\int \text{d}^{3}\boldsymbol{w}~\frac{f_{Ms0}}{T_{s0}}\right)^{-1}\left(\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}\right)^{2}=n_{0}T_{s0}(\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743})^{2},\end{eqnarray}$$

which implies the lower bound (1.2) for the KMHD incremental potential energy, hence the result that stability in isothermal ideal MHD is a sufficient condition for stability in KMHD.

5 The ‘force-times-displacement’ functional

The search for a necessary condition for KMHD stability will involve the ‘force-times-displacement’ functional $U$ defined in (1.6). From the linearized momentum equation (3.20) and the expression for the ‘perpendicular ideal-MHD’ potential energy (4.9), this functional is $U=\unicode[STIX]{x1D6FF}W_{\bot }^{F}+U^{K}$ with

(5.1) $$\begin{eqnarray}U^{K}=-\frac{1}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}(\unicode[STIX]{x1D709}_{\Vert }\hat{F}_{s\Vert }[\hat{f}_{s}]+\unicode[STIX]{x1D743}_{\bot }\boldsymbol{\cdot }\hat{F}_{s\bot }[\hat{f}_{s}]).\end{eqnarray}$$

Unlike the fluid part of $U$ that equals $\unicode[STIX]{x1D6FF}W_{\bot }^{F}$ , and as a consequence of the lack of self-adjointness in the kinetic part of the force operator, $U^{K}$ is in general different from $\unicode[STIX]{x1D6FF}W^{K}$ and an explicit evaluation of $U^{K}$ is necessary. This can be carried out along the lines of the evaluation of $\unicode[STIX]{x1D6FF}W^{K}$ in the previous section. In complete analogy with the procedure followed to arrive at (4.11), one obtains

(5.2) $$\begin{eqnarray}U^{K}=-\frac{1}{2}\mathop{\sum }_{s=i,e}m_{s}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,\hat{f}_{s}^{even}w^{2}Q[\unicode[STIX]{x1D743}].\end{eqnarray}$$

The next step is to substitute the drift-kinetic evolution equation for $\hat{f}_{s}^{even}$ . However, what is needed here is not (3.23) but its integrated form with respect to time,

(5.3) $$\begin{eqnarray}\hat{f}_{s}^{even}-\hat{f}_{s}^{even}(0)+w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\boldsymbol{x}}\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })+Q[\unicode[STIX]{x1D743}-\unicode[STIX]{x1D743}(0)]~\frac{m_{s}w^{2}f_{Ms0}}{T_{s0}}=0.\end{eqnarray}$$

After substituting this time-integrated form, (5.2) becomes

(5.4) $$\begin{eqnarray}\displaystyle U^{K} & = & \displaystyle \frac{1}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,\frac{T_{s0}}{f_{Ms0}}~\hat{f}_{s}^{even}\left[\hat{f}_{s}^{even}-\hat{f}_{s}^{even}(0)+w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}\boldsymbol{x}}\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })\right]\nonumber\\ \displaystyle & & \displaystyle -\,\frac{1}{2}\mathop{\sum }_{s=i,e}m_{s}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}^{even}w^{2}Q[\unicode[STIX]{x1D743}(0)]\end{eqnarray}$$

and, after integration by parts with a vanishing boundary term,

(5.5) $$\begin{eqnarray}\displaystyle U^{K} & = & \displaystyle \frac{1}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}~\frac{T_{s0}}{f_{Ms0}}\left[(\hat{f}_{s}^{even})^{2}-\left(w_{\Vert }\boldsymbol{b}_{0}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{even}}{\unicode[STIX]{x2202}\boldsymbol{x}}\right)~\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })\right]\nonumber\\ \displaystyle & & \displaystyle -\,\frac{1}{2}\mathop{\sum }_{s=i,e}m_{s}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}^{even}\left\{\frac{T_{s0}}{m_{s}f_{Ms0}}\hat{f}_{s}^{even}(0)+w^{2}Q[\unicode[STIX]{x1D743}(0)]\right\}.\end{eqnarray}$$

Finally, substituting the drift-kinetic equation for $\hat{f}_{s}^{odd}$ (3.24) and using again $\int \text{d}^{3}\boldsymbol{w}\,w_{\Vert }\hat{f}_{s}^{odd}=0$ ,

(5.6) $$\begin{eqnarray}\displaystyle U^{K} & = & \displaystyle \frac{1}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\frac{T_{s0}}{f_{Ms0}}\left[(\hat{f}_{s}^{even})^{2}+\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{odd}}{\unicode[STIX]{x2202}t}\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })\right]\nonumber\\ \displaystyle & & \displaystyle -~\frac{1}{2}\mathop{\sum }_{s=i,e}m_{s}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}^{even}\left\{\frac{T_{s0}}{m_{s}f_{Ms0}}\hat{f}_{s}^{even}(0)+w^{2}Q[\unicode[STIX]{x1D743}(0)]\right\},\end{eqnarray}$$

hence

(5.7) $$\begin{eqnarray}\displaystyle U^{K} & = & \displaystyle \unicode[STIX]{x1D6FF}W^{K}-\frac{1}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,\frac{T_{s0}}{f_{Ms0}}\left[(\hat{f}_{s}^{odd})^{2}-\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{odd}}{\unicode[STIX]{x2202}t}\,\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })\right]\nonumber\\ \displaystyle & & \displaystyle -\,\frac{1}{2}\mathop{\sum }_{s=i,e}m_{s}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}^{even}\left\{\frac{T_{s0}}{m_{s}f_{Ms0}}\hat{f}_{s}^{even}(0)+w^{2}Q[\unicode[STIX]{x1D743}(0)]\right\}.\end{eqnarray}$$

This result provides an explicit formula for the difference $\unicode[STIX]{x1D6FF}W-U$ which, due to the lack of KMHD self-adjointness, is not zero in general, as anticipated.

A special class of perturbations make the difference $\unicode[STIX]{x1D6FF}W-U$ particularly simple, albeit not zero. They are those for which the initial value of the distribution functions is

(5.8) $$\begin{eqnarray}\hat{f}_{s}(0)=\hat{f}_{s}^{even}(0)=-Q[\unicode[STIX]{x1D743}(0)]\frac{m_{s}w^{2}f_{Ms0}}{T_{s0}},\end{eqnarray}$$

that satisfies the constraints $\int \text{d}^{3}\boldsymbol{w}~\hat{f}_{s}(0)=-n_{0}\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743}(0)$ and $\int \text{d}^{3}\boldsymbol{w}~w_{\Vert }\hat{f}_{s}(0)=0$ , and makes the last term of (5.7) vanish. For this special class, the ‘force-times-displacement’ functional reduces to $U=\unicode[STIX]{x1D6FF}W-R/2$ , with

(5.9) $$\begin{eqnarray}R=\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}~\frac{T_{s0}}{f_{Ms0}}\left[(\hat{f}_{s}^{odd})^{2}-\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{odd}}{\unicode[STIX]{x2202}t}~\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })\right],\end{eqnarray}$$

and the initial value of the kinetic term of the potential energy functional is

(5.10) $$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D6FF}W^{K}[\hat{f}_{s}(0)] & = & \displaystyle \frac{1}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,Q[\unicode[STIX]{x1D743}(0)]^{2}\frac{m_{s}^{2}w^{4}f_{Ms0}}{T_{s0}}\nonumber\\ \displaystyle & = & \displaystyle \frac{1}{6}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}~n_{0}T_{s0}\left\{5[\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743}(0)]^{2}+[\unicode[STIX]{x1D735}\boldsymbol{\cdot }\unicode[STIX]{x1D743}(0)-3(\boldsymbol{b}_{0}\boldsymbol{b}_{0})\boldsymbol{ : }(\unicode[STIX]{x1D735}\unicode[STIX]{x1D743}(0))]^{2}\right\}.\nonumber\\ \displaystyle & & \displaystyle\end{eqnarray}$$

Thus, for this class of perturbations, the initial value of the KMHD incremental potential energy is equal to that of the double-adiabatic fluid theory (Chew et al. Reference Chew, Goldberger and Low1956; Bernstein et al. Reference Bernstein, Frieman, Kruskal and Kurlsrud1958),

(5.11) $$\begin{eqnarray}\unicode[STIX]{x1D6FF}W[\unicode[STIX]{x1D743}(0),\hat{f}_{s}(0)]=\unicode[STIX]{x1D6FF}W_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }(0)]+\unicode[STIX]{x1D6FF}W^{K}[\hat{f}_{s}(0)]=\unicode[STIX]{x1D6FF}W^{DA}[\unicode[STIX]{x1D743}(0)],\end{eqnarray}$$

where $\unicode[STIX]{x1D6FF}W^{DA}[\unicode[STIX]{x1D743}(0)]$ includes the contribution of the parallel displacement $\unicode[STIX]{x1D709}_{\Vert }(0)$ . This choice of initial condition will be used next in the investigation of a necessary condition for KMHD stability.

6 Necessary condition for KMHD stability

The energy principle provides a necessary and sufficient condition for stability in the double-adiabatic fluid theory (Bernstein et al. Reference Bernstein, Frieman, Kruskal and Kurlsrud1958; Laval et al. Reference Laval, Mercier and Pellat1965). This means that, if an equilibrium is unstable in the double-adiabatic theory, a trial fluid displacement $\unicode[STIX]{x1D743}^{tr}$ exists such that $\unicode[STIX]{x1D6FF}W^{DA}[\unicode[STIX]{x1D743}^{tr}]<0$ . Such a trial fluid displacement may have a non-zero parallel component. Consider then the following initial condition for the KMHD system:

(6.1a,b ) $$\begin{eqnarray}\unicode[STIX]{x1D743}(0)=\unicode[STIX]{x1D743}^{tr},\quad \boldsymbol{u}_{1}(0)=\unicode[STIX]{x2202}\unicode[STIX]{x1D743}(0)/\unicode[STIX]{x2202}t=0,\end{eqnarray}$$

(6.2) $$\begin{eqnarray}\hat{f}_{s}(0)=\hat{f}_{s}^{even}(0)=-Q[\unicode[STIX]{x1D743}^{tr}]\frac{m_{s}w^{2}f_{Ms0}}{T_{s0}}.\end{eqnarray}$$

This perturbation belongs to the special class discussed in the previous section, therefore the initial value of the KMHD incremental potential energy is

(6.3) $$\begin{eqnarray}\unicode[STIX]{x1D6FF}W(0)=\unicode[STIX]{x1D6FF}W_{\bot }^{F}[\unicode[STIX]{x1D743}_{\bot }(0)]+\unicode[STIX]{x1D6FF}W^{K}[\hat{f}_{s}(0)]=\unicode[STIX]{x1D6FF}W^{DA}[\unicode[STIX]{x1D743}^{tr}]<0\end{eqnarray}$$

and the initial value of the kinetic energy is $K(0)=0$ . Defining the time-dependent fluid displacement norm,

(6.4) $$\begin{eqnarray}N(t)\equiv \frac{1}{2}\int \text{d}^{3}\boldsymbol{x}\,\unicode[STIX]{x1D70C}_{0}|\unicode[STIX]{x1D743}(t)|^{2},\end{eqnarray}$$

its time derivatives are

(6.5) $$\begin{eqnarray}\displaystyle & \displaystyle \frac{\text{d}N\,(t)}{\text{d}t}=\int \text{d}^{3}\boldsymbol{x}\,\unicode[STIX]{x1D70C}_{0}~\unicode[STIX]{x1D743}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D743}}{\unicode[STIX]{x2202}t}, & \displaystyle\end{eqnarray}$$

(6.6) $$\begin{eqnarray}\displaystyle & \displaystyle \frac{\text{d}^{2}N(t)}{\text{d}t^{2}}=\int \text{d}^{3}\boldsymbol{x}\,\unicode[STIX]{x1D70C}_{0}\left[\Big|\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D743}}{\unicode[STIX]{x2202}t}\Big|^{2}+\unicode[STIX]{x1D743}\boldsymbol{\cdot }\frac{\unicode[STIX]{x2202}^{2}\unicode[STIX]{x1D743}}{\unicode[STIX]{x2202}t^{2}}\right]=2K(t)-2U(t). & \displaystyle\end{eqnarray}$$

The considered perturbation belongs to the class for which $U(t)=\unicode[STIX]{x1D6FF}W(t)-R(t)/2$ and energy conservation implies $\unicode[STIX]{x1D6FF}W(t)+K(t)=\unicode[STIX]{x1D6FF}W(0)+K(0)=\unicode[STIX]{x1D6FF}W(0)$ , hence

(6.7) $$\begin{eqnarray}\frac{\text{d}^{2}N(t)}{\text{d}t^{2}}=4K(t)-2\unicode[STIX]{x1D6FF}W(0)+R(t)\end{eqnarray}$$

with $R(t)$ as defined in (5.9). The initial condition $\unicode[STIX]{x2202}\unicode[STIX]{x1D743}(0)/\unicode[STIX]{x2202}t=0$ implies $\text{d}N\,(0)/\text{d}t=0$ . Then, after time integration,

(6.8) $$\begin{eqnarray}N(t)=N(0)-\unicode[STIX]{x1D6FF}W(0)~t^{2}+N_{R}(t)+4\int _{0}^{t}\,\text{d}t^{\prime }\int _{0}^{t^{\prime }}\,\text{d}t^{\prime \prime }\,K(t^{\prime \prime }),\end{eqnarray}$$

with

(6.9) $$\begin{eqnarray}N_{R}(t)=\int _{0}^{t}\,\text{d}t^{\prime }\int _{0}^{t^{\prime }}\,\text{d}t^{\prime \prime }\,R(t^{\prime \prime }).\end{eqnarray}$$

A general lower bound for $N_{R}(t)$ can be derived. Calling

(6.10) $$\begin{eqnarray}g_{s}\equiv (\hat{f}_{s}^{odd})^{2}-\frac{\unicode[STIX]{x2202}\hat{f}_{s}^{odd}}{\unicode[STIX]{x2202}t}\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime }),\end{eqnarray}$$

one has

(6.11) $$\begin{eqnarray}N_{R}(t)=\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}~\frac{T_{s0}}{f_{Ms0}}\int _{0}^{t}\,\text{d}t^{\prime }\int _{0}^{t^{\prime }}\,\text{d}t^{\prime \prime }\,g_{s}(t^{\prime \prime }).\end{eqnarray}$$

Simple integration by parts yields

(6.12) $$\begin{eqnarray}\displaystyle \int _{0}^{t}\,\text{d}t^{\prime }\int _{0}^{t^{\prime }}\,\text{d}t^{\prime \prime }\,g_{s}(t^{\prime \prime }) & = & \displaystyle 2\int _{0}^{t}\,\text{d}t^{\prime }\int _{0}^{t^{\prime }}\,\text{d}t^{\prime \prime }\,[\hat{f}_{s}^{odd}(t^{\prime \prime })]^{2}-\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })\int _{0}^{t^{\prime }}\,\text{d}t^{\prime \prime }\,\hat{f}_{s}^{odd}(t^{\prime \prime })\nonumber\\ \displaystyle & = & \displaystyle 2\int _{0}^{t}\,\text{d}t^{\prime }\,(t-t^{\prime })[\hat{f}_{s}^{odd}(t^{\prime })]^{2}-\frac{1}{2}\left[\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })\right]^{2}\end{eqnarray}$$

and, by virtue of the Cauchy–Schwarz inequality

(6.13) $$\begin{eqnarray}0\leqslant \left[\int _{0}^{t}\,\text{d}t^{\prime }\,\hat{f}_{s}^{odd}(t^{\prime })\right]^{2}~\leqslant t\int _{0}^{t}\,\text{d}t^{\prime }\,[\hat{f}_{s}^{odd}(t^{\prime })]^{2},\end{eqnarray}$$

one obtains

(6.14) $$\begin{eqnarray}\int _{0}^{t}\,\text{d}t^{\prime }\int _{0}^{t^{\prime }}\,\text{d}t^{\prime \prime }\,g_{s}(t^{\prime \prime })~\geqslant \frac{1}{2}\int _{0}^{t}\,\text{d}t^{\prime }\,(3t-4t^{\prime })[\hat{f}_{s}^{odd}(t^{\prime })]^{2}.\end{eqnarray}$$

This yields the lower bound for $N_{R}(t)$

(6.15) $$\begin{eqnarray}N_{R}(t)\geqslant \frac{t^{2}}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,\frac{T_{s0}}{f_{Ms0}}\unicode[STIX]{x1D6F7}_{s}(t),\end{eqnarray}$$

where

(6.16) $$\begin{eqnarray}\unicode[STIX]{x1D6F7}_{s}(t)=\int _{0}^{1}\,\text{d}\unicode[STIX]{x1D708}\,(3-4\unicode[STIX]{x1D708})[\hat{f}_{s}^{odd}(\unicode[STIX]{x1D708}t)]^{2}.\end{eqnarray}$$

Since $K(t)\geqslant 0$ , this implies also the lower bound for $N(t)$

(6.17) $$\begin{eqnarray}N(t)\geqslant N(0)-\unicode[STIX]{x1D6FF}W(0)t^{2}+\frac{t^{2}}{2}\mathop{\sum }_{s=i,e}\int \text{d}^{3}\boldsymbol{x}\int \text{d}^{3}\boldsymbol{w}\,\frac{T_{s0}}{f_{Ms0}}\unicode[STIX]{x1D6F7}_{s}(t).\end{eqnarray}$$

The considered perturbation has $\unicode[STIX]{x1D6FF}W(0)=\unicode[STIX]{x1D6FF}W^{DA}(0)<0$ , so $-\unicode[STIX]{x1D6FF}W(0)t^{2}$ grows in time without bound. If the equilibrium under consideration is KMHD stable, the norm $N(t)$ must be a bounded function of time and (6.17) implies that this can only happen if $\unicode[STIX]{x1D6F7}_{s}(t)$ becomes negative and not negligible as $t\rightarrow \infty$ . Therefore, if it were possible to prove that for any KMHD-stable equilibrium the long-time behaviour of $\hat{f}_{s}^{odd}(t)$ is such that, as $t\rightarrow \infty$ , $\unicode[STIX]{x1D6F7}_{s}(t)$ either stays non-negative or tends to zero, the double-adiabatic unstable equilibrium under consideration must be KMHD unstable and that would complete the proof that stability according to the double-adiabatic fluid theory is a necessary condition for KMHD stability. The KMHD theory is time reversible and dissipationless, but includes the possible Landau damping of initial-value solutions as a consequence of the mixing of the phases of different spectral components. It would then seem reasonable to expect that, for perturbations about a KMHD-stable equilibrium, the long-time behaviour of $[\hat{f}_{s}^{odd}(t)]^{2}$ is such that it either tends to a time-independent limit or oscillates about some time-independent average, with deviations of finite amplitude evenly distributed above and below the average. If one characterizes such a behaviour by the condition that the limit

(6.18) $$\begin{eqnarray}(\hat{f}_{\infty s}^{odd})^{2}=\lim _{t\rightarrow \infty }\frac{1}{t}\int _{0}^{t}\,\text{d}t^{\prime }\,[\hat{f}_{s}^{odd}(t^{\prime })]^{2}\end{eqnarray}$$

exists and $0\leqslant [\hat{f}_{s}^{odd}(t\rightarrow \infty )]^{2}\leqslant 2(\hat{f}_{\infty s}^{odd})^{2}$ , then the lemma given in appendix D guarantees that, as $t\rightarrow \infty$ , $\unicode[STIX]{x1D6F7}_{s}(t)$ either stays non-negative or tends to zero. Actually, appendix D shows that this is guaranteed by the weaker condition that $\hat{f}_{s}^{odd}(t)$ be bounded and, for sufficiently large values of $t$ ,

(6.19) $$\begin{eqnarray}[\hat{f}_{s}^{odd}(t)]^{2}\leqslant \frac{3}{t}\int _{0}^{t}\,\text{d}t^{\prime }\,[\hat{f}_{s}^{odd}(t^{\prime })]^{2}\end{eqnarray}$$

regardless of the existence of the limit (6.18). The validity of this condition means that linear KMHD would not allow the possibility that, after a sufficiently long time has elapsed, $[\hat{f}_{s}^{odd}(t)]^{2}$ underwent short bursts during which it would reach large values relative to the prevalent, time-averaged ones. Such bursty time variation, reminiscent of the so-called ‘rogue wave’ phenomenon, is precluded by the condition (6.19). In conclusion, if the long-time behaviour of oscillations about a KMHD-stable equilibrium satisfies the condition (6.19), it has been demonstrated that stability in the double-adiabatic fluid theory would be a necessary condition for stability in KMHD.

The form of the initial condition for the distribution functions (6.2) used here is related to that of the zero-frequency normal modes, although it is not the same. The distribution functions of the zero-frequency normal mode solution of Rosenbluth & Rostoker (Reference Rosenbluth and Rostoker1959) are given in terms of the corresponding fluid displacement eigenvector $\unicode[STIX]{x1D743}$ by

(6.20) $$\begin{eqnarray}\hat{f}_{s}^{\unicode[STIX]{x1D714}=0}=-\langle Q[\unicode[STIX]{x1D743}]\rangle _{\unicode[STIX]{x1D70F}}\frac{m_{s}w^{2}f_{Ms0}}{T_{s0}},\end{eqnarray}$$

where $\langle \cdots \rangle _{\unicode[STIX]{x1D70F}}$ is the transit-bounce average

(6.21) $$\begin{eqnarray}\langle Q\rangle _{\unicode[STIX]{x1D70F}}=\frac{\oint \text{d}l\,\mathop{w}_{\Vert }^{-1}Q}{\oint \text{d}l\,\mathop{w}_{\Vert }^{-1}}.\end{eqnarray}$$

These contour integrals are taken along one period of the particle phase-space trajectories that follow the magnetic field lines at constant $w$ and $\unicode[STIX]{x1D706}$ , under the assumption that such orbits are periodic (trapped particles or passing particles on closed magnetic loops) or nearly periodic (passing particles on ergodic magnetic lines, with the argument that such lines approximately close on themselves). As a consequence of the identity (A 5), the contribution of $\unicode[STIX]{x1D709}_{\Vert }$ to $\langle Q[\unicode[STIX]{x1D743}]\rangle _{\unicode[STIX]{x1D70F}}$ vanishes, so the zero-frequency normal-mode distribution functions depend only on the perpendicular displacement, $\langle Q[\unicode[STIX]{x1D743}]\rangle _{\unicode[STIX]{x1D70F}}=\langle Q[\unicode[STIX]{x1D743}_{\bot }]\rangle _{\unicode[STIX]{x1D70F}}$ . The classic comparison theorem (1.5) (Kruskal & Oberman Reference Kruskal and Oberman1958; Rosenbluth & Rostoker Reference Rosenbluth and Rostoker1959) applies to a generic perpendicular displacement $\unicode[STIX]{x1D743}_{\bot }$ and can be stated as

(6.22) $$\begin{eqnarray}\unicode[STIX]{x1D6FF}W^{RR}[\unicode[STIX]{x1D743}_{\bot }]=\unicode[STIX]{x1D6FF}W\left[\unicode[STIX]{x1D743}_{\bot },-\langle Q[\unicode[STIX]{x1D743}_{\bot }]\rangle _{\unicode[STIX]{x1D70F}}\frac{m_{s}w^{2}f_{Ms0}}{T_{s0}}\right]\leqslant \unicode[STIX]{x1D6FF}W\left[\unicode[STIX]{x1D743}_{\bot },-Q[\unicode[STIX]{x1D743}_{\bot }]\frac{m_{s}w^{2}f_{Ms0}}{T_{s0}}\right]=\unicode[STIX]{x1D6FF}W^{DA}[\unicode[STIX]{x1D743}_{\bot }]\end{eqnarray}$$

but, in order to prove the inequality, one needs the periodic or nearly periodic orbit assumption. The present study uses the initial distribution functions $\hat{f}_{s}(0)=-Q[\unicode[STIX]{x1D743}(0)]~m_{s}w^{2}T_{s0}^{-1}f_{Ms0}$ with an initial $\unicode[STIX]{x1D743}(0)$ that can have a non-zero parallel component and makes $\unicode[STIX]{x1D6FF}W^{DA}[\unicode[STIX]{x1D743}(0)]$ negative. Therefore, the initial KMHD incremental potential energy is known to be negative irrespective of orbit periodicity considerations and the analysis is equally valid for non-periodic passing particles. The argument of near periodicity on ergodic magnetic lines is not needed here. Besides, the present study looked for a proof of instability by actually demonstrating that the chosen initial condition would result in an unbounded growth of the perturbation, not just that the incremental potential energy would become negative.

7 Summary

This work has investigated the linear stability properties of KMHD from the initial-value point of view. The analysis has been based on the formulation introduced in Ramos (Reference Ramos2015a ,Reference Ramos b ) to study the KMHD normal modes. This formulation uses drift-kinetic equations in the reference frame of the complete macroscopic fluid velocity, eliminates the electric field and guarantees that the quasineutrality condition is satisfied automatically at all times if it is satisfied by the initial condition. Moreover, it does not solve the kinetic equations by integration along characteristics and does not require that the particle orbits be periodic or nearly periodic. Such a formulation is unconventional compared with more traditional approaches, so the article was written as a self-contained presentation and §§ 2 and 3 were devoted to the detailed introduction of the adopted KMHD description, both in its nonlinear and linearized versions. The derivation of the KMHD potential energy functional within this framework was then given in § 4. The main new results are in §§ 5 and 6. Section 5 evaluated explicitly the ‘force-times-displacement’ functional, emphasizing the fact that it is different from the potential energy in KMHD, a result that can be taken as an independent proof that the KMHD force operator is not self-adjoint. Finally, § 6 looked for a proof that stability according to the double-adiabatic fluid theory is a necessary condition for KMHD stability by studying the possibility that, if an isotropic-pressure static equilibrium is double-adiabatic unstable, an appropriately chosen KMHD initial condition would grow in time without bound. This was demonstrated with the proviso that the long-time behaviour of oscillations about KMHD-stable equilibria satisfies the condition (6.19).

The two results regarding KMHD stability that this work has arrived at, namely that stability in isothermal ideal MHD is sufficient for stability in KMHD and that, provided the condition (6.19) holds, stability in the Chew–Goldberger–Low double-adiabatic theory including the variation of the parallel fluid displacement would be necessary for stability in KMHD, do not rely on the classic comparison theorems of Kruskal & Oberman (Reference Kruskal and Oberman1958), Rosenbluth & Rostoker (Reference Rosenbluth and Rostoker1959). They do not require the periodicity or near periodicity of particle orbits and are equally valid for non-periodic passing particles. In addition, the necessary condition analysis does not require the self-adjointness of the force operator or the existence of a complete basis of normal modes, which are properties that have not been proved for KMHD. However, the necessary condition result hinges on the validity of the hypothesis (6.19) on the long-time behaviour of oscillations about KMHD-stable equilibria. This is a reasonable assumption but remains a yet to be proved mathematical conjecture, which underscores the difficulty of finding a completely rigorous proof of a necessary condition for KMHD stability.