 Incompressible magnetic displacement energy identity

For a static ideal magnetohydrodynamics equilibrium in fixed gravity, let $\nabla \cdot \xi =0$ and $\Pi =p+B^2/(2{\mu }_0)$. Frozen flux gives $b=(B\cdot \nabla )\xi -(\xi \cdot \nabla )B$ and continuity gives ${\rho }^{\prime }=-\xi \cdot \nabla \rho$. Multiplying the linear momentum equation by the conjugate displacement and integrating by parts gives the displayed identity when surface terms vanish. The two magnetic cross terms cancel; the remaining magnetic equilibrium derivative combines with the mass density-gradient force to give the two Hessian matrices. Real symmetric Hessian matrices and the nonnegative magnetic tension term make ${\sigma }^2$ real.