Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract

Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces and a continuous rotational transform. These currents have two components: a surface current (Dirac $\hat{I}´$-function in flux surface labeling) that prevents the formation of magnetic islands and an algebraically divergent Pfirsch--Schlüter current density when a pressure gradient is present across the rational surface. At flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch--Schlüter current density scaling as $J\sim 1/\hat{I}”{\hat{I}}^1$, where $\hat{I}”{\hat{I}}^1$ is the difference of the rotational transform relative to the rational surface. If the distance $s$ between flux surfaces is proportional to $\hat{I}”{\hat{I}}^1$, the scaling relation $J\sim 1/\hat{I}”{\hat{I}}^1\sim 1/s$ will lead to a paradox that the Pfirsch--Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm\textendash Kulsrud\textendash Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch--Schlüter current density but also the diamagnetic current density are divergent as $\sim 1/\hat{I}”{\hat{I}}^1$. However, due to the formation of a Dirac $\hat{I}´$-function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with $s\sim (\hat{I}”{\hat{I}}^1{)}^2$. Consequently, the singular current density $J\sim 1/\sqrt{s}$, making the total current finite, thus resolving the paradox.

Minor corrections of typos