Twisted magnetar magnetospheres: a class of semi-analytical force-free non-rotating solutions

Magnetospheric twists, that is magnetospheres with a toroidal component, are under scrutiny due to the key role the twist is believed to play in the behaviour of neutron stars. Notably, its dissipation is believed to power magnetar activity, and is an important element of the evolution of these stars. We exhibit a new class of twisted axisymmetric force-free magnetospheric solutions. We solve the Grad–Shafranov equation by introducing an ansatz akin to a multipolar expansion. We obtain a hierarchical system of ordinary differential equations where lower-order multipoles source the higher-order ones. We show that analytical approximations can be obtained, and that in general solutions can be numerically computed using standard ordinary-differential-equation solvers. We obtain a class of solutions with a great flexibility in initial conditions, and show that a subset of these asymptotically tends to vacuum. The twist is not confined to a subset of field lines. The solutions are symmetric about the equator, with a toroidal component that can be reversed. This symmetry is supported by an equatorial current sheet. We provide a first-order approximation of a particular solution that consists of the superposition of a vacuum dipole and a toroidal magnetic field sourced by the dipole, where the toroidal component decays as $1/r^4$ , where $r$ is the radial coordinate. As an example of strongly multipolar solution, we also exhibit cases with an additional octupole component.