Significance: The observed diversity of magnetars indicates that their magnetic topology is more complicated than a simple dipole. Current models of their radiative emission, based on axially symmetric simulations, require the presence of a concealed toroidal magnetic field having up to 100 times more energy than the observed dipole component, but the physical origin of such a field is unclear. Our fully 3D simulations of the crustal magnetic field demonstrate that magnetic instabilities operate under a range of plausible conditions and generate small-scale field structures that are an order of magnitude stronger than the large-scale field. The Maxwell stress and Ohmic heating from these structures can explain magnetar bursts and surface hotspots, using comparable poloidal and toroidal magnetic fields.

We present the first fully self-consistent three dimensional model of a neutron star’s magnetic field, generated by electric currents in the star’s crust via the Hall effect. We find that the global-scale field converges to a dipolar Hall-attractor state, as seen in recent axisymmetric models, but that small-scale features in the magnetic field survive even on much longer time scales. These small-scale features propagate toward the dipole equator, where the crustal electric currents organize themselves into a strong equatorial jet. By calculating the distribution of magnetic stresses in the crust, we predict that neutron stars with fields stronger than $10^{14}\mathrm{G}$ can still be subject to starquakes more than $10^5\mathrm{yr}$ after their formation.

The propagation of cosmic rays in turbulent magnetic fields is a diffusive process driven by the scattering of the charged particles by random magnetic fluctuations. Such fields are usually highly intermittent, consisting of intense magnetic filaments and ribbons surrounded by weaker, unstructured fluctuations. Studies of cosmic-ray propagation have largely overlooked intermittency, instead adopting Gaussian random magnetic fields. Using test particle simulations, we calculate cosmic-ray diffusivity in intermittent, dynamo-generated magnetic fields. The results are compared with those obtained from non-intermittent magnetic fields having identical power spectra. The presence of magnetic intermittency significantly enhances cosmic-ray diffusion over a wide range of particle energies. We demonstrate that the results can be interpreted in terms of a correlated random walk.

The speed of sound greatly exceeds typical flow velocities in many stellar and planetary interiors. To follow the slow evolution of subsonic motions, various sound-proof models attempt to remove fast acoustic waves while retaining stratified convection and buoyancy dynamics. In astrophysics, anelastic models typically receive the most attention in the class of sound-filtered stratified models. Generally, anelastic models remain valid in nearly adiabatically stratified regions like stellar convection zones, but may break down in strongly sub-adiabatic, stably stratified layers common in stellar radiative zones. However, studying stellar rotation, circulation, and dynamos requires understanding the complex coupling between convection and radiative zones, and this requires robust equations valid in both regimes. Here we extend the analysis of equation sets begun in Brown et al., which studied anelastic models, to two types of pseudo-incompressible models. This class of models has received attention in atmospheric applications, and more recently in studies of white-dwarf supernova progenitors. We demonstrate that one model conserves energy but the other does not. We use Lagrangian variational methods to extend the energy conserving model to a general equation of state, and dub the resulting equation set the generalized pseudo-incompressible (GPI) model. We show that the GPI equations suitably capture low-frequency phenomena in both convection and radiative zones in stars and other stratified systems, and we provide recommendations for converting low-Mach number codes to this equation set.