  Motivated by explosive releases of energy in fusion, space and astrophysical
plasmas, we consider the nonlinear convective stability of stratified
magnetohydrodynamic (MHD) equilibria in 2D. We demonstrate that, unlike the
Schwarzschild criterion in hydrodynamics (``entropy must increase upwards for
convective stability''), the so-called modified Schwarzschild criterion for 2D
MHD (or for any kind of fluid dynamics with more than one source of pressure)
guarantees only linear stability. As a result, in 2D MHD (unlike in
hydrodynamics) there exist metastable equilibria that are unstable to nonlinear
perturbations despite being stable to linear ones. We show that the
minimum-energy configurations attainable by these atmospheres via non-diffusive
reorganization can be determined by solving a combinatorial optimization
problem. We find inter alia that these minimum-energy states are usually 2D,
even when the original metastable equilibrium was 1D. We demonstrate with
direct numerical simulations that these 2D states are fairly accurate
predictors of the final state reached by laminar relaxation of metastable
equilibria at small Reynolds number. To describe relaxation at large Reynolds
number, we construct a statistical mechanical theory based on the maximization
of Boltzmann's mixing entropy that is analogous to the Lynden-Bell statistical
mechanics of self-gravitating systems and collisionless plasmas, and to the
Robert-Sommeria-Miller (RSM) theory of 2D vortex turbulence. The minimum-energy
states described above are, we show, the low-temperature limit of this theory.
We demonstrate that the predictions of the statistical mechanics are in
reasonable agreement with direct numerical simulations.
