Nested toroidal flux surfaces in magnetohydrostatics. Generation via soliton theory

It is shown that the classical magnetohydrostatic equations of an infinitely conducting fluid reduce to the integrable potential Heisenberg spin equation subject to a Jacobian condition if the magnitude of the magnetic field is constant along individual magnetic field lines. Any solution of the constrained potential Heisenberg spin equation gives rise to a multiplicity of magnetohydrostatic equilibria which share the magnetic field line geometry. The multiplicity of equilibria is reflected by the local arbitrariness of the total pressure profile. A connection with the classical Da Rios equations is exploited to establish the existence of associated helically and rotationally symmetric equilibria. As an illustration, Palumbo's ‘unique’ toroidal isodynamic equilibrium is retrieved.