 Solution

Use the magnetic energy $E_B=(2{\mu }_0{)}^{-1}{\int }_V|\mathbf{B}{|}^{2}dV$. Vary the Euler potentials while keeping their boundary values fixed. The field variation isLet $\delta \mathbf{A}=\delta \alpha \nabla \beta -\delta \beta \nabla \alpha$. Integration by parts givesbecause $\delta \mathbf{A}=0$ on the boundary. ThereforeThe two interior variations are independent. At a smooth minimum their coefficients must vanish, so the curl is perpendicular to both potential gradients and hence parallel to $\mathbf{B}$ wherever $\mathbf{B}\ne 0$:This is a force-free magnetic field: the Lorentz force density vanishes. Taking the divergence also gives $\mathbf{B}\cdot \nabla \lambda =0$, so the proportionality factor is constant along a field line but need not be the same on different lines.

The fixed-label variation preserves the encoded field-line connection rather than allowing arbitrary reconnection during minimization. The argument proves a necessary condition for a smooth minimum, not that every force-free stationary field is a stable energy minimum. It also does not include gas pressure or gravity in the minimized functional.