Toby S. Wood

Energy Conservation and Gravity Waves in Sound-proof Treatments of Stellar Interiors. II. Lagrangian Constrained Analysis

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.

Three Dimensional Simulation of the Magnetic Stress in a Neutron Star Crust

We present the first fully self-consistent three dimensional model of a neutron stars magnetic field, generated by electric currents in the stars 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} G can still be subject to starquakes more than 10^{5} yr after their formation.

Magnetic field evolution in magnetar crusts through three-dimensional simulations

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. Current models of magnetars require extremely strong magnetic fields to explain their observed quiescent and bursting emission, implying that the field strength within the star’s outer crust is orders of magnitude larger than the dipole component inferred from spin-down measurements. This presents a serious challenge to theories of magnetic field generation in a proto-neutron star. Here, we present detailed modeling of the evolution of the magnetic field in the crust of a neutron star through 3D simulations. We find that, in the plausible scenario of equipartition of energy between global-scale poloidal and toroidal magnetic components, magnetic instabilities transfer energy to nonaxisymmetric, kilometer-sized magnetic features, in which the local field strength can greatly exceed that of the global-scale field. These intense small-scale magnetic features can induce high-energy bursts through local crust yielding, and the localized enhancement of Ohmic heating can power the star’s persistent emission. Thus, the observed diversity in magnetar behavior can be explained with mixed poloidal−toroidal fields of comparable energies.

Dynamics of the solar tachocline – III. Numerical solutions of the Gough and McIntyre model

We present the first numerical simulations of the solar interior to exhibit a tachocline consistent with the Gough and McIntyre (1998) model. We find nonlinear, axisymmetric, steady-state numerical solutions in which: (1) a large-scale primordial field is confined within the radiation zone by downwelling meridional flows that are gyroscopically pumped in the convection zone (2) the radiation zone is in almost-uniform rotation, with a rotation rate consistent with observations (3) the bulk of the tachocline is magnetic free, in thermal-wind balance and in thermal equilibrium and (4) the interaction between the field and the flows takes place within a very thin magnetic boundary layer, the tachopause, located at the bottom of the tachocline. We show that the thickness of the tachocline scales with the amplitude of the meridional flows exactly as predicted by Gough and McIntyre. We also determine the parameter conditions under which such solutions can be obtained, and provide a simple explanation for the failure of previous numerical attempts at reproducing the Gough and McIntyre model. Finally, we discuss the implications of our findings for future numerical models of the solar interior, and for future observations of the Sun and other stars.

Polar confinement of the Sun's interior magnetic field by laminar magnetostrophic flow

The global-scale interior magnetic field B i needed to account for the Suns observed differential rotation can be effective only if confined below the convection zone in all latitudes including, most critically, the polar caps. Axisymmetric solutions are obtained to the nonlinear magnetohydrodynamic equations showing that such polar confinement can be brought about by a very weak downwelling flow U ~ 10 −5 cms −1 over each pole. Such downwelling is consistent with the helioseismic evidence. All three components of the magnetic field B decay exponentially with altitude across a thin, laminar ‘magnetic confinement layer’ located at the bottom of the tachocline and permeated by spiralling field lines. With realistic parameter values, the thickness of the confinement layer ~10 −3 of the Suns radius, the thickness scale being the magnetic advection–diffusion scale δ = η/ U where the magnetic (ohmic) diffusivity η ≈ 4.1 × 10 2 cm 2 s −1 . Alongside baroclinic effects and stable thermal stratification, the solutions take into account the stable compositional stratification of the helium settling layer, if present as in todays Sun, and the small diffusivity of helium through hydrogen, χ ≈ 0.9 × 10 1 cm 2 s −1 . The small value of χ relative to η produces a double boundary-layer structure in which a ‘helium sublayer‘ of smaller vertical scale (χ/η) 1/2 δ is sandwiched between the top of the helium settling layer and the rest of the confinement layer. Solutions are obtained using both semi-analytical and purely numerical, finite-difference techniques. The confinement-layer flows are magnetostrophic to excellent approximation. More precisely, the principal force balances are between Lorentz, Coriolis, pressure-gradient and buoyancy forces, with relative accelerations negligible to excellent approximation. Viscous forces are also negligible, even in the helium sublayer where shears are greatest. This is despite the kinematic viscosity being somewhat greater than χ. We discuss how the confinement layers s at each pole might fit into a global dynamical picture of the solar tachocline. That picture, in turn, suggests a new insight into the early Sun and into the longstanding enigma of solar lithium depletion.

Density-shear instability in electron magneto-hydrodynamics

We discuss a novel instability in inertia-less electron magneto-hydrodynamics (EMHD), which arises from a combination of electron velocity shear and electron density gradients. The unstable modes have a lengthscale longer than the transverse density scale, and a growth-rate of the order of the inverse Hall timescale. We suggest that this density-shear instability may be of importance in magnetic reconnection regions on scales smaller than the ion skin depth, and in neutron star crusts. We demonstrate that the so-called Hall drift instability, previously argued to be relevant in neutron star crusts, is a resistive tearing instability rather than an instability of the Hall term itself. We argue that the density-shear instability is of greater significance in neutron stars than the tearing instability, because it generally has a faster growth-rate and is less sensitive to geometry and boundary conditions. We prove that, for uniform electron density, EMHD is “at least as stable” as regular, incompressible MHD, in the sense that any field configuration that is stable in MHD is also stable in EMHD. We present a connection between the density-shear instability in EMHD and the magneto-buoyancy instability in anelastic MHD.We discuss a novel instability in inertia-less electron magneto-hydrodynamics (EMHD), which arises from a combination of electron velocity shear and electron density gradients. The unstable modes have a lengthscale longer than the transverse density scale, and a growth-rate of the order of the inverse Hall timescale. We suggest that this density-shear instability may be of importance in magnetic reconnection regions on scales smaller than the ion skin depth, and in neutron star crusts. We demonstrate that the so-called Hall drift instability, previously argued to be relevant in neutron star crusts, is a resistive tearing instability rather than an instability of the Hall term itself. We argue that the density-shear instability is of greater significance in neutron stars than the tearing instability, because it generally has a faster growthrate and is less sensitive to geometry and boundary conditions. We prove that, for uniform electron density, EMHD is “at least as stable” as regular, incompressible MHD, in the sense that any field configuration that is stable in MHD is also stable in EMHD. We present a connection between the density-shear instability in EMHD and the magneto-buoyancy instability in anelastic MHD.

The Sun's meridional circulation and interior magnetic field

To date, no self-consistent numerical simulation of the solar interior has succeeded in reproducing the observed thinness of the solar tachocline and the persistence of uniform rotation beneath it. Although it is known that the uniform rotation can be explained by the presence of a global-scale confined magnetic field, numerical simulations have thus far failed to produce any solution where such a field remains confined against outward diffusion. We argue that the problem lies in the choice of parameters for which these numerical simulations have been performed. We construct a simple analytical magnetohydrodynamic model of the solar interior and identify several distinct parameter regimes. For realistic solar parameter values, our results are in broad agreement with the tachocline model of Gough & McIntyre. In this regime, meridional flows driven at the base of the convection zone are of sufficient amplitude to hold back the interior magnetic field against diffusion. For the parameter values used in existing numerical simulations, on the other hand, we find that meridional flows are significantly weaker and, we argue, unable to confine the interior field. We propose a method for selecting parameter values in future numerical models.

TRANSPORT BY MERIDIONAL CIRCULATIONS IN SOLAR-TYPE STARS

Transport by meridional flows has significant consequences for stellar evolution, but is difficult to capture in global-scale numerical simulations because of the wide range of timescales involved. Stellar evolution models therefore usually adopt parameterizations for such transport based on idealized laminar or mean-field models. Unfortunately, recent attempts to model this transport in global simulations have produced results that are not consistent with any of these idealized models. In an effort to explain the discrepancies between global simulations and idealized models, here we use three-dimensional local Cartesian simulations of compressible convection to study the efficiency of transport by meridional flows below a convection zone in several parameter regimes of relevance to the Sun and solar-type stars. In these local simulations we are able to establish the correct ordering of dynamical timescales, although the separation of the timescales remains unrealistic. We find that, even though the generation of internal waves by convective overshoot produces a high degree of time dependence in the meridional flow field, the mean flow has the qualitative behavior predicted by laminar, balanced models. In particular, we observe a progressive deepening, or burrowing, of the mean circulation if the local Eddington-Sweet timescale is shorter than the viscous diffusion timescale. Such burrowing is a robust prediction of laminar models in this parameter regime, but has never been observed in any previous numerical simulation. We argue that previous simulations therefore underestimate the transport by meridional flows.

The solar tachocline: a self-consistent model of magnetic confinement

This work was supported by a Research Studentship from the Science and Technology Facilities Council

Cosmic Rays in Intermittent Magnetic Fields

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 relying on Gaussian random magnetic fields. Using test particle simulations, we investigate 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.