Spontaneous wandering of the magnetic axis in pulsars: 3D magneto-thermal simulations and the imprint on braking indices
Abstract
Braking index measurements of hundreds of pulsars, despite observational caveats, show a clear trend: n<3 at early ages, n>3 at middle ages, and increasingly large, fluctuating values at larger characteristic ages. Crustal magnetic fields in isolated neutron stars evolve through Ohmic dissipation and non-linear Hall redistribution across scales, with additional early-time contributions from the chiral magnetic effect and post-burial re-emergence. Initial conditions are usually assumed dipole-dominated, whereas more complex and realistic fields, with energy spread more evenly over a broad range of scales, have not been explored in this context. We assess how the internal dynamics of initially non-trivial fields change the dipolar field strength and obliquity, and how their combined contribution to the braking index n, together with that of the alignment, compares with the observed values. To this end we perform long-term 3D magneto-thermal simulations, for configurations ranging from the dipole-dominated case to tangled, small-scale-dominated fields. Small-scale dominated configurations with weak dipolar components show rich dynamics: the dipolar mode is continuously fed by the small scales and may grow or decay depending on the initial spectrum, driving n<3 or n>3, while the obliquity varies on timescales as short as cal O(kyr). These spontaneous changes often dominate the alignment torque and can account for the large values of n−3 of either sign observed at characteristic ages τc≳105 yr, unlike the dipole-dominated models. However, short-term magnetospheric fluctuations or superfluid-driven variations of the torque may still be needed to explain the measurements and the structured timing residuals quantitatively. Future work comparing these models with the observed timing properties could fine-tune the most promising initial configurations.
Related
- Magnetic conditions governing the transition to penumbral magneto-convection
- On Photon Propagation in Strong Magnetic Fields: Time Delays from Pulsars
- Interior Magnetic Fields in Magnetars and Radio Pulsars
- Radiation Magneto-hydrodynamic Simulations of MRI with Zero-Net-Vertical-Flux: Necessity of Resolving the Thermal Scale For Strongly Magnetized Initial Conditions
- Simulations of a Conducting Sphere Moving through Magnetized Plasma II: Effect of Internal Magnetic Fields and Enhancement of Alfvénic Drag
BibTeX
@misc{dehman2026spontaneous,
title = {{Spontaneous wandering of the magnetic axis in pulsars: 3D magneto-thermal simulations and the imprint on braking indices}},
author = {Clara Dehman and Daniele Viganò},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2609.17436},
eprint = {2609.17436},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.17436},
}
RIS
TY - UNPB
TI - Spontaneous wandering of the magnetic axis in pulsars: 3D magneto-thermal simulations and the imprint on braking indices
AU - Clara Dehman
AU - Daniele Viganò
PY - 2026
DA - 2026/10/28
PB - arXiv
DO - 10.48550/arxiv.2609.17436
UR - https://arxiv.org/abs/2609.17436
AB - Braking index measurements of hundreds of pulsars, despite observational caveats, show a clear trend: n<3 at early ages, n>3 at middle ages, and increasingly large, fluctuating values at larger characteristic ages. Crustal magnetic fields in isolated neutron stars evolve through Ohmic dissipation and non-linear Hall redistribution across scales, with additional early-time contributions from the chiral magnetic effect and post-burial re-emergence. Initial conditions are usually assumed dipole-dominated, whereas more complex and realistic fields, with energy spread more evenly over a broad range of scales, have not been explored in this context. We assess how the internal dynamics of initially non-trivial fields change the dipolar field strength and obliquity, and how their combined contribution to the braking index n, together with that of the alignment, compares with the observed values. To this end we perform long-term 3D magneto-thermal simulations, for configurations ranging from the dipole-dominated case to tangled, small-scale-dominated fields. Small-scale dominated configurations with weak dipolar components show rich dynamics: the dipolar mode is continuously fed by the small scales and may grow or decay depending on the initial spectrum, driving n<3 or n>3, while the obliquity varies on timescales as short as cal O(kyr). These spontaneous changes often dominate the alignment torque and can account for the large values of n−3 of either sign observed at characteristic ages τ_c≳10^5 yr, unlike the dipole-dominated models. However, short-term magnetospheric fluctuations or superfluid-driven variations of the torque may still be needed to explain the measurements and the structured timing residuals quantitatively. Future work comparing these models with the observed timing properties could fine-tune the most promising initial configurations.
ER -