2. Equilibrium

The macroscopic fluid equilibrium whose stability properties we wish to assess arises from the equilibrium condition for an electron plasma in a magnetic field,

(2.1)\begin{equation} m_e n_e \boldsymbol{v}_e \boldsymbol{\cdot} \boldsymbol{\nabla} \boldsymbol{v}_e + \boldsymbol{\nabla}( n_e T_e) -en_e \boldsymbol{\nabla} \phi={-}en_e \boldsymbol{v}_e \times \boldsymbol{B},\end{equation}

where $n_e$ is the number density of the electron plasma, $T_e$ is its temperature and $\boldsymbol {v}_e$ is its mean-flow velocity. The non-negligible equilibrium electrostatic potential $\phi$ is a result of the non-neutrality of the single species plasma and is given by Poisson's equation $\varepsilon _0\nabla ^2 \phi = e n_e$. This potential generates an electrostatically repulsive force on the plasma $e n_e \boldsymbol {\nabla } \phi$ which, along with the pressure gradient force $\boldsymbol {\nabla } (n_e T_e)$ and the inertial force $m_e n_e \boldsymbol {v}_e \boldsymbol {\cdot } \boldsymbol {\nabla } \boldsymbol {v}_e$, must be balanced by the Lorentz force of the electron flow in the confining magnetic field $\boldsymbol {B}$, given by the right-hand side of (2.1) (Pedersen & Boozer Reference Pedersen and Boozer2002).

In a stellarator, the confining magnetic field has a toroidal topology and can be represented in contravariant form as $\boldsymbol {B} = \boldsymbol {\nabla } \psi \times \boldsymbol {\nabla } \theta + \iota (\psi ) \boldsymbol {\nabla } \varphi \times \boldsymbol {\nabla } \psi$, where $2\pi\psi$ is the toroidal magnetic flux that labels the surfaces of constant toroidal magnetic flux (flux surfaces), $\theta$ and $\varphi$ are $2{\rm \pi}$-periodic poloidal and toroidal angles, respectively, and $\iota (\psi )$ is the rotational transform which is defined as

(2.2)\begin{equation} \iota(\psi) = \frac{\text{d}\theta}{\text{d}\varphi}. \end{equation}

The stellarator is set apart from other means of confining a single species plasma, like the magnetic dipole or purely toroidal field traps, by its finite rotational transform, which defines the pitch of magnetic field lines on each of the nested flux surfaces $\psi$. If $\iota$ is irrational (or is a sufficiently high-order rational), then the magnetic field lines touch each part of these toroidal surfaces evenly (Helander Reference Helander2014).

To confine an electron plasma in such a magnetic field, the electron density must satisfy $n_e \ll n_B$, where $n_B = \varepsilon _0 B^2/(2m_e)$ is the Brillouin density (Boozer Reference Boozer2005). In this limit, the inertia of the plasma is negligible and the equilibrium condition reduces to

(2.3)\begin{equation} \boldsymbol{\nabla}(n_e T_e) -en_e \boldsymbol{\nabla} \phi={-}en_e \boldsymbol{v}_e \times \boldsymbol{B}.\end{equation}

We now consider the force balance along magnetic field lines by taking the scalar product of (2.3) with $\boldsymbol {B}$. In a stellarator with well defined flux surfaces, this amounts to considering the force balance over the entire surface $\psi$. If it is assumed that the flux surfaces are isotherms with $T_e = T_e(\psi )$Footnote ¹, we find that the quantity $T_e(\psi )\ln n_e - e\phi = H(\psi )$ is constant on each flux surface. Solving this for the electron density yields the Boltzmann relation $n_e = N(\psi ) \exp {(e\phi /T_e(\psi ))}$. Inserting this expression for the density into Poisson's equation yields the nonlinear Poisson–Boltzmann equation,

(2.4)\begin{equation} \nabla^2 \phi = N(\psi) \exp\left(\frac{e\phi}{T_e(\psi)}\right), \end{equation}

the solution to which gives the equilibrium electrostatic potential $\phi (\psi, \theta, \varphi )$ for a given $N(\psi )$ and temperature $T_e(\psi )$. The function $N(\psi )$, which has remained arbitrary until this point, is related to the flux-surface average of the density by $\langle n_e \rangle = N(\psi )\langle \exp ({e\phi }/{T_e(\psi )})\rangle$, with the flux-surface average defined as (Helander Reference Helander2014)

(2.5)\begin{equation} \langle\ldots\rangle = \frac{\displaystyle \int (\ldots) \sqrt{g}\, {\rm d}\theta \, {\rm d}\varphi}{\displaystyle \int\sqrt{g}\, {\rm d}\theta \, {\rm d}\varphi} = \frac{1}{V'(\psi)}\int (\ldots) \sqrt{g}\, {\rm d}\theta \, {\rm d}\varphi, \end{equation}

where the Jacobian is given by

(2.6)\begin{equation} \frac{1}{\sqrt{g}}= (\boldsymbol{\nabla} \psi \times \boldsymbol{\nabla} \theta) \boldsymbol{\cdot} \boldsymbol{\nabla}\varphi, \end{equation}

and $V(\psi )$ denotes the volume enclosed by the flux surface $\psi$.

The nonlinear partial differential equation (2.4) is the central relation for the equilibria of pure electron plasmas in stellarators. It has been demonstrated numerically that the equilibria which satisfy solutions to this equation behave like minimum-energy states, in that the electron plasma density in equilibrium is attracted to the region of most positive potential in the device, thus minimising the total electrostatic energy of the system (Lefrancois et al. Reference Lefrancois, Pedersen, Boozer and Kremer2005). This has also been shown to be the case experimentally, as the electron density has been found to increase at the point on each flux surface which is closest to the grounded outer wall of the device (Himura et al. Reference Himura, Wakabayashi, Yamamoto, Isobe, Okamura, Matsuoka, Sanpei and Masamune2007).

The stability of this equilibrium has been treated in previous work by Boozer (Reference Boozer2004), where similar stability properties were derived. However, in his analysis, the thermal energy of the plasma was not considered, despite the inclusion of a finite electron temperature $T_e$ in the plasma entropy. Without the inclusion of a finite plasma pressure in the variational principle, the equilibrium, which is a balance of a pressure gradient and electrostatic repulsion, cannot be found to be the minimum-energy state. Here, we conclusively show that the equilibria represented by (2.4) in stellarators minimise the physical energy of the system.

3. Stability theorem

We consider a plasma confined on the magnetic surfaces of a stellarator occupying a confinement volume $V$ bounded by a conducting wall $\partial V$, with an outward-pointing radial normal vector $\boldsymbol {\nabla } \psi /|\boldsymbol {\nabla } \psi |$. The total physical energy of the electron plasma is given by

(3.1)\begin{equation} W = \int_V \left(\frac{1}{2}m_e n_e v_e^2 + \frac{p}{\gamma - 1} +\frac{|\boldsymbol{B}|^2}{2\mu_0} + \varepsilon_0\frac{|\boldsymbol{E}|^2}{2}\right)\,\text{d}V, \end{equation}

where $v_e = |\boldsymbol {v}_e|$ is the magnitude of the flow velocity of the electron fluid, $p = n_eT_e$ is the thermal pressure, $\gamma$ is the adiabatic index, and $\boldsymbol {E}$ and $\boldsymbol {B}$ are the electric and magnetic fields, respectively. We can compare the relative sizes of the contributions to the energy by estimating the equilibrium electric field as $E \sim aen_e/\varepsilon _0$, where $a$ is the minor radius. We find that the sizes of the electrostatic and thermal energies scale as

(3.2)\begin{equation} \frac{\varepsilon_0|\boldsymbol{E}|^2}{n_e T_e} \sim \frac{a^2}{\lambda_D^2}. \end{equation}

In order for a collection of electrons to qualify as a plasma, the parameter $a/\lambda _D$ needs to exceed unity, so that the pressure makes a smaller contribution to the energy than the electric field, but we shall nevertheless treat these energies as comparable, i.e. of the same order in our expansion procedure. We can similarly estimate the ratio of the kinetic energy of the plasma flow to the pressure by using the fact that the electron flow is largely due to the $\boldsymbol {E}\times \boldsymbol {B}$-drift, such that $v_e \sim E/B$. This yields a scaling

(3.3)\begin{equation} \frac{m_e n_e v_e^2}{n_e T_e} \sim \frac{1}{2} \frac{a^2}{\lambda_D^2}\frac{\rho_e^2}{\lambda_D^2} =\frac{a^2}{\lambda_D^2}\frac{n_e}{2 n_B}, \end{equation}

where $\rho _e= \sqrt {(2 T_e/m_e)}/\varOmega _e$ is the electron gyroradius with $\varOmega _e = |e|B/m_e$, and $n_B$ denotes the Brillouin density. It has been shown that in order to confine an electron plasma in a stellarator, it is required that the plasma density be far below the Brillouin density, $n_e/n_B \ll 1$ (Boozer Reference Boozer2005). We can see that this requirement is equivalent to assuming that the Debye length is much larger than the electron gyroradius, such that $\rho _e /\lambda _D \ll 1$. As a result of this necessary criterion for confinement, the inertial contribution to the energy is much smaller than the pressure contribution and can be safely neglected.

Furthermore, while the magnetic energy, ${|\boldsymbol {B}|^2}/({2\mu _0})$, may be large, it does not have a meaningful impact on the variation of the energy under a perturbation of the plasma, given that, by the estimate $\delta B \sim a \mu _0 e n_e v_e$, we find

(3.4)\begin{equation} \frac{\delta B}{B} \sim \left(\frac{v_{Te}}{c}\right)^2\frac{a^2}{\lambda_D^2}\frac{n_e}{4n_B}, \end{equation}

where $v_{Te} = \sqrt {2T_e/m_e}$ is the thermal velocity and $c$ is the speed of light in vacuum. Once more the requirement that $n_e /n_B \ll 1$ makes this ratio small unless $a/\lambda _D$ is extremely large, and therefore we can safely assume that $\boldsymbol {B}$ is entirely generated by the external coils of the confinement device and can be treated as constant in our analysis.

Thus, in the limit $n_e/n_B \ll 1$, the only terms of relevance to the analysis of the stability of the system are the thermal and electrostatic energies, allowing us to reduce (3.1) to

(3.5)\begin{equation} W = \int_V \left(\frac{p}{\gamma - 1} + \varepsilon_0\frac{|\boldsymbol{E}|^2}{2}\right)\,\text{d}V. \end{equation}

This energy can be thought of as a functional, $W = W[n_e, T_e, \boldsymbol {E}]$, such that for any given functions $(n_e, T_e, \boldsymbol {E})$, $W$ gives a real number as an output, the energy of the plasma.

The treatment of this variational energy principle here is twofold. Firstly, it will be shown that equilibria of the equations of motion correspond to unique extrema of the energy functional $W$, such that the first-order variation $\delta W$ vanishes at equilibrium. Secondly, the nature of such extrema will be investigated, and it will be determined whether they are maxima, minima or ‘saddle points’ in the energy space.

The nature of the extremum, which is determined by the sign of $\delta W$ to second order, establishes the stability properties of the equilibrium, with positive definiteness indicating a minimum, negative definiteness indicating maximum and an indeterminate sign indicating a saddle point. The stability of an equilibrium, where $\delta W$ vanishes to first order, is determined by the nature of the corresponding extremum, where so-called ‘formal’ stability is implied by a unique maximum or minimum (Holm et al. Reference Holm, Marsden, Ratiu and Weinstein1985). Because $W$ is conserved by an electron plasma under our assumption of $n_e \ll n_B$, an equilibrium which represents a unique maximum or minimum defined by the functions $(n_1, T_1, \boldsymbol {E}_1)$ cannot access any other state $(n_2, T_2, \boldsymbol {E}_2)$ without a change in energy. Of course, in a physical system, some degree of energy dissipation will be present, so in this sense an energy minimum is more desirable as it is additionally resilient in the face of energy dissipation. If the equilibrium is a degenerate extremum or a saddle point in energy space, stability is not guaranteed.

If one naïvely extremises the energy $W$ in an unconstrained way, an uninteresting $n_e = T_e = \boldsymbol {E} = 0$ solution will be found. Thus, an extremum of (3.5) must be found under a set of physical constraints. Our analysis aims to expand upon the treatment of Boozer (Reference Boozer2004), by considering the same physical constraints with the addition that here, the temperature is not only allowed to contribute to the entropy of the plasma, but also its thermal energy. We now consider the following constraints.

1. (i) The electron plasma must obey Gauss’ law of electrostatics, such that $\boldsymbol {E}$ satisfies $\boldsymbol {\nabla } \boldsymbol {\cdot } \boldsymbol {E} = -{en_e}/{\varepsilon _0}$ everywhere in $V$.

2. (ii) The electrostatic potential $\phi$ vanishes on a surface $\partial V$ surrounding the plasma. This surface thus corresponds to a grounded, electrically conducting wall.

3. (iii) The number of electrons on each flux surface is conserved by the magnetic confinement, and thus the electron flow cannot cross flux surfaces. This condition implies that $\langle n_e \rangle = {G}(\psi )$, where $G(\psi )$ is a fixed function of $\psi$. Since the total charge enclosed by the plasma (or any region bounded by a flux surface within the plasma) is conserved, any variation of the electric field $\delta \boldsymbol {E}$ must then satisfy

   (3.6)\begin{equation} \int_{\partial V} \delta \boldsymbol{E} \boldsymbol{\cdot} \boldsymbol{\nabla} \psi \frac{\text{d}S}{| \boldsymbol{\nabla} \psi |} = 0, \end{equation}
   according to Gauss's law.

4. (iv) The average entropy of the plasma on each flux surface

   (3.7)\begin{equation} S(\psi) =\frac{1}{\gamma - 1}\langle n_e \ln{(T_e/n_e^{\gamma - 1})} \rangle, \end{equation}
   must be conserved or increase on each flux surface $\psi$. This constraint is motivated by the existence of magnetic surfaces which allow for rapid heat conduction along surfaces, but poor conduction between any two flux surfaces.

These constraints can be included when seeking an extremum of (3.5) through the method of Lagrange multipliers. This gives the constrained energy,

(3.8)\begin{align} \bar{W} & = \int_{V} \left[\frac{\varepsilon_0}{2}|{\boldsymbol{E}}|^2+ \frac{n_e T_e}{\gamma -1} -\eta(\boldsymbol{r})(\varepsilon_0\boldsymbol{\nabla} \boldsymbol{\cdot} \boldsymbol{E} + en_e)\right.\nonumber\\ & \quad - \left.\lambda(\psi)(n_e- {G}(\psi)) - \mu(\psi)\left( \frac{1}{\gamma - 1}n_e \ln{(T_e/n_e^{\gamma -1})} - {S}(\psi)\right) \right]\,\text{d}V. \end{align}

Here, the Lagrange multipliers $\eta$, $\lambda$ and $\mu$ guarantee that any extremum of the above expression satisfies our constraints. For example varying $\lambda$ gives

(3.9)\begin{equation} \delta_\lambda \bar{W} = \int \text{d}\psi \delta \lambda(\psi)( \langle n_e \rangle - G(\psi)), \end{equation}

such that, at an extremum where $\delta \bar {W}$ vanishes to first order, the constraint $\langle n_e \rangle = G(\psi )$ must be satisfied. Note that we have stipulated that $\lambda = \lambda (\psi )$ and $\mu = \mu (\psi )$ are functions of $\psi$ alone, indicating that the corresponding constraints are only enforced on average on each flux surface, and not at every point in the plasma, as is the case with Gauss’ law, where $\eta = \eta (\boldsymbol {r})$ has been used to ensure that it is satisfied at every point in space. Now an equilibrium is sought for which $\delta \bar {W}$ vanishes to first order under the individual variation of $\boldsymbol {E}$, $T_e$ and $n_e$. By computing $\delta _{\boldsymbol {E}}\bar {W}$ to first order we find

(3.10)\begin{equation} \delta_{\boldsymbol{E}} \bar{W} = \int_V (\varepsilon_0 \boldsymbol{E} \boldsymbol{\cdot} \delta \boldsymbol{E} - \eta(\boldsymbol{r})(\varepsilon_0\boldsymbol{\nabla} \boldsymbol{\cdot} \delta \boldsymbol{E}))\,\text{d}V . \end{equation}

Integrating by parts and equating to zero gives

(3.11)\begin{equation} \int_{V} (\boldsymbol{E} + \boldsymbol{\nabla} \eta)\boldsymbol{\cdot} \delta \boldsymbol{E}\, \text{d}V -\int_{\partial V} \eta(\boldsymbol{r})\delta \boldsymbol{E}\boldsymbol{\cdot}{\boldsymbol{\nabla}\psi} \frac{\text{d}S}{|\boldsymbol{\nabla} \psi|} = 0, \end{equation}

where the divergence theorem has been used to arrive at the surface integral over the boundary of $V$. It is now required that both integrals vanish for all choices of the function $\delta \boldsymbol {E}$. The volume integral will do so if $\boldsymbol {E} = - \boldsymbol {\nabla } \eta$, implying that, at the extremum, the Lagrange multiplier $\eta (\boldsymbol {r})$ is equal to the electrostatic potential $\phi (\boldsymbol {r})$ within an additive constant. The surface integral over $\partial V$ then also vanishes, which implies $\phi = \text {const.}$ on the boundary, in accordance with the second constraint on our list above.

The first-order variation $\delta _{T_e} \bar {W}$ with respect to the temperature is found to be

(3.12)\begin{equation} \delta_{T_e} \bar{W} = \int_V \left[ n_e\delta T_e \left(1 - \frac{\mu(\psi)}{T_e} \right) \right]\,\text{d}V = 0, \end{equation}

and implies that $T_e = T_e(\psi )$ depends on $\psi$ alone, and that $\mu (\psi ) = T_e(\psi )$ at equilibrium.

Finally, the first-order variation with respect to the density $\delta _{n_e} \bar {W}$ then gives

(3.13)\begin{align} \delta_{n_e} \bar{W}& = \int_V \delta n_e\left[ \frac{T_e(\psi)}{\gamma - 1} - e\phi -{\lambda(\psi)} \right.\nonumber\\ & \quad -\left.\frac{1}{\gamma - 1}T_e(\psi)\ln{T_e(\psi)} +T_e(\psi)\ln{n_e} + T_e(\psi) \right] \,\text{d}V = 0. \end{align}

This admits solutions of the form

(3.14)\begin{equation} T_e(\psi)\ln (n_e) -e \phi = H(\psi), \end{equation}

where the function $H(\psi )$ is given by

(3.15)\begin{equation} H(\psi) = {\lambda(\psi)} - \frac{1}{\gamma - 1}T_e(\psi) + \frac{T_e(\psi)}{\gamma -1}\ln(T_e(\psi)) - T_e(\psi). \end{equation}

Here, $H(\psi )$ amounts to some arbitrary flux function, whose arbitrariness is given by the fact that at equilibrium the Lagrange multiplier $\lambda (\psi )$ is unspecified. From this, the equilibrium condition is found to be

(3.16)\begin{equation} n_e = N(\psi)\exp{(e\phi/T_e(\psi))}. \end{equation}

Thus, it is clear that the finite-temperature equilibria of interest, given by (2.4), extremise $\bar {W}$. With each equilibrium being uniquely specified by the choice of $N(\psi )$ and $T_e(\psi )$, the uniqueness of the equilibrium is guaranteed thanks to the uniqueness theorem of the nonlinear Poisson–Boltzmann equation on magnetic surfaces (Durand de Gevigney Reference Durand de Gevigney2011). Now, all that remains is to determine if such an extremum is a maximum, minimum or saddle point in energy space.

In order to facilitate a more direct interpretation of the role of our constraints on the stability of the system, in computing the second variation we abandon the method of Lagrange multipliers in favour of directly enforcing the constraints when necessary. We now consider the perturbed energy on each flux surface $\psi$ of the equilibrium to second order,

(3.17)\begin{equation} \delta U(\psi) = \left\langle \varepsilon_0 \boldsymbol{\nabla} \phi \boldsymbol{\cdot} \boldsymbol{\nabla} \delta \phi + \varepsilon_0\frac{ |\boldsymbol{\nabla} \delta \phi|^2}{2} + \frac{1}{\gamma -1} (n_e \delta T_e + T_e\delta n_e + \delta T_e \delta n_e) \right\rangle, \end{equation}

where we have used $\boldsymbol {E} = - \boldsymbol {\nabla } \phi$. Under such a general perturbation about the equilibrium, the corresponding change in the flux-surface average of the entropy to second order is

(3.18)\begin{align} \delta S(\psi) & = \frac{1}{\gamma -1 }\left\langle ( \ln T_e - (\gamma -1)\ln n_e - (\gamma -1)) \delta n_e + \frac{n_e}{T_e}\delta T_e \right.\nonumber\\ & \quad - \left.\frac{1}{2}\left(\frac{n_e}{T_e^2}\delta T_e^2 + \frac{(\gamma -1 )}{n_e} \delta n_e^2\right) +\frac{1}{T_e}\delta n_e \delta T_e \right\rangle. \end{align}

By enforcing the constraint that $\langle n_e \rangle = G(\psi )$, such that $\langle \delta n_e \rangle = 0$, and exploiting the fact that the equilibrium temperature is of the form $T_e = T_e (\psi )$, the changes in the energy and entropy of each flux surface, respectively, are found to reduce to

(3.19)\begin{equation} \delta U(\psi) = \left\langle \varepsilon_0 \boldsymbol{\nabla} \phi \boldsymbol{\cdot} \boldsymbol{\nabla} \delta \phi + \varepsilon_0 \frac{|\boldsymbol{\nabla} \delta \phi|^2}{2} + \frac{1}{\gamma -1} (n_e \delta T_e +\delta T_e \delta n_e) \right\rangle, \end{equation}

and

(3.20)\begin{align} \delta S(\psi) & = \left\langle - \ln(n_e)\delta n_e + \frac{n_e\delta T_e}{T_e(\gamma -1)}\right.\nonumber\\ & \quad - \left.\frac{1}{2}\left(\frac{n_e\delta T_e^2}{T_e^2(\gamma -1)} + \frac{\delta n_e^2}{n_e}\right) +\frac{\delta n_e \delta T_e}{T_e(\gamma -1)} \right\rangle. \end{align}

We now consider the linear combination

(3.21)\begin{equation} \delta F(\psi) = \delta U(\psi) - {T_e(\psi)}\delta S(\psi), \end{equation}

where the choice of the factor $T_e(\psi )$ is motivated by the optimal Lagrange multiplier for enforcing the entropy constraint identified above. The function $\delta F(\psi )$ can be interpreted as the change in the Helmholtz free energy of the electron plasma on each flux surface. If we multiply $\delta F(\psi )$ by $V'(\psi )$ and integrate over all the flux surfaces to arrive at a volume integral, we obtain

(3.22)\begin{align} \int \text{d} \psi V'(\psi) \delta F(\psi)& = \int \text{d}\psi V'(\psi)\left\langle \varepsilon_0 \boldsymbol{\nabla} \phi \boldsymbol{\cdot} \boldsymbol{\nabla} \delta \phi + \varepsilon_0\frac{ |\boldsymbol{\nabla} \delta \phi|^2}{2} + T_e(\psi)\ln(n_e)\delta n_e \right.\nonumber\\ & \quad +\left.\frac{T_e(\psi)}{\gamma-1}\left(\frac{n_e}{2T_e^2}\delta T_e^2 + \frac{(\gamma -1 )}{2n_e} \delta n_e^2 \right)\right\rangle. \end{align}

The first term in this expression may be rewritten by integrating by parts, noting that $\phi$ vanishes on the boundary, and using Poisson's equation $\varepsilon _0\nabla ^2 \delta \phi = e\delta n_e$. Using the equilibrium condition $n_e = N(\psi ) \exp ({e \phi /T_e(\psi )})$ and once again enforcing $\langle \delta n_e \rangle = 0$, we find

(3.23)\begin{equation} \int_V \delta F \text{d} V = \int_V \left( \varepsilon_0\frac{ |\boldsymbol{\nabla} \delta \phi|^2}{2} + \frac{n_e}{2 (\gamma - 1) T_e}\delta T_e^2 + \frac{T_e}{2n_e} \delta n_e^2 \right) \, \text{d}V ,\end{equation}

which is clearly positive definite. Hence, the change in energy,

(3.24)\begin{equation} \delta W = \int_V ( \delta F + {T_e \delta S} ) \, \text{d}V, \end{equation}

is positive definite if $\delta S \geq 0$. Thus, the electron plasma equilibrium in a stellarator given by (2.4) minimises the energy $W$ under the constraints we have imposed.

This argument can also be considered from the perspective of equilibrium as a state of maximum entropy, where in this case the average energy on each flux surface is held fixed. Both approaches yield the same result for the equilibrium and stability (an example of the entropy maximisation approach can be seen in Malmberg & O'Neil (Reference Malmberg and O'Neil1977)). This can be most readily seen in (3.21), where one could instead consider $\delta \bar{S} = -\delta F/T_e(\psi )$ as the entropy change constrained by holding the flux-surface average of the energy fixed. The volume-integrated constrained entropy change,

(3.25)\begin{equation} \int_V \delta \bar{S} \text{d}V ={-} \int_V \frac{\delta F}{T_e(\psi)} \text{d}V = \int_V \left(\delta S - \frac{\delta U}{T_e(\psi)} \right) \, \text{d}V, \end{equation}

is then clearly negative definite at equilibrium from (3.23) such that the entropy is maximised.