Diffusivity in force-free simulations of global magnetospheres

Abstract

Abstract: Assuming that the numerical diffusivity triggered by violations of the force-free electrodynamics constraints is a proxy for the physical resistivity, we examine its impact on the overall dynamics of force-free aligned pulsar magnetospheres endowed with an equatorial current sheet. We assess the constraint violations as a diffusivity source. The effects of modifications on electric fields used to restore force-free conditions are not confined to the equatorial current sheet, but modify the magnetospheric dynamics on timescales shorter than the pulsar rotational period. These corrections propagate especially via a channel that was unexplored, namely, changes induced to the electric charge density, $\rho$. We quantify the global consequences of diffusivity by comparing different techniques to model $\rho$. By default, we combine a conservative $\rho$-evolution with hyperbolic/parabolic cleaning of inaccuracies in the Maxwell equations. As an alternative, we enforce a constrained evolution where $\rho$ is directly computed as the electric field divergence. The conservative approach reduces the Poynting flux dissipated in the equatorial current sheet by an order of magnitude, along with an increase of the pulsar luminosity driven by a shift of the Y-point location. The luminosity changes according to $L_{\mathrm{Y}}\propto {\alpha }^{0.11}$, where $\alpha$ is the ratio of diffusion to advection timescales, controlling the amount of (numerical) diffusivity. Our models suggest interpreting the luminosity dependence on the Y-point location as differences in resistivities encountered at the equatorial current sheet. Alternatively, they could be interpreted in terms of the pair formation multiplicity, $\kappa$, smaller diffusion being consistent with $\kappa \gg 1$.

(17 pages, 9 figures, Accepted by MNRAS)

References in corpus (13)

• Ab-initio pulsar magnetosphere: three-dimensional particle-in-cell simulations of oblique pulsars
• Particle acceleration in axisymmetric pulsar current sheets
• Simulations of relativistic collisionless shocks: shock structure and particle acceleration
• Electrodynamics of axisymmetric pulsar magnetosphere with electron-positron discharge: a numerical experiment
• Three-dimensional numerical simulations of the pulsar magentoshere: Preliminary results
• Multi-dimensional Numerical Scheme for Resistive Relativistic MHD
• General Relativistic Force-Free Electrodynamics: A New Code and Applications to Black Hole Magnetospheres
• Accurate Simulations of Binary Black-Hole Mergers in Force-Free Electrodynamics
• Tearing instability in relativistic magnetically dominated plasmas
• Electrodynamics of pulsar magnetospheres
• A Constrained Transport Method for the Solution of the Resistive Relativistic MHD Equations
• The role of reconnection in the pulsar magnetosphere
• A ring-of-fire in the pulsar magnetosphere

Cited by in corpus (2)

• Three-dimensional dynamics of strongly twisted magnetar magnetospheres: Kinking flux tubes and global eruptions
• Electromagnetic precursors to black hole - neutron star gravitational wave events: Flares and reconnection-powered fast-radio transients from the late inspiral