skip to main content


Title: On large-scale dynamos with stable stratification and the application to stellar radiative zones
Abstract Our understanding of large-scale magnetic fields in stellar radiative zones remains fragmented and incomplete. Such magnetic fields, which must be produced by some form of dynamo mechanism, are thought to dominate angular-momentum transport, making them crucial to stellar evolution. A major difficulty is the effect of stable stratification, which generally suppresses dynamo action. We explore the effects of stable stratification on mean-field dynamo theory with a particular focus on a non-helical large-scale dynamo (LSD) mechanism known as the magnetic shear-current effect. We find that the mechanism is robust to increasing stable stratification as long as the original requirements for its operation are met: a source of shear and non-helical magnetic fluctuations (e.g. from a small-scale dynamo). Both are plausibly sourced in the presence of differential rotation. Our idealized direct numerical simulations, supported by mean-field theory, demonstrate the generation of near equipartition large-scale toroidal fields. Additionally, a scan over magnetic Reynolds number shows no change in the growth or saturation of the LSD, providing good numerical evidence of a dynamo mechanism resilient to catastrophic quenching, which has been an issue for helical dynamos. These properties – the absence of catastrophic quenching and robustness to stable stratification – make the mechanism a plausible candidate for generating in situ large-scale magnetic fields in stellar radiative zones.  more » « less
Award ID(s):
1804048
NSF-PAR ID:
10402209
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
517
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
526 to 542
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Magnetic fields are now widely recognized as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Using imposed magnetic fields cannot reliably model the interstellar medium's (ISM) dynamical structure nor phase interactions. Dynamos must be modeled. ISM models exist of turbulent magnetic fields using small-scale dynamo (SSD). Others model the large-scale dynamo (LSD) organizing magnetic fields at the scale of the disk or spiral arms. Separately, neither can fully describe the galactic magnetic field dynamics nor topology. We model the LSD and SSD together at a sufficient resolution to use the low explicit Lagrangian resistivity required. The galactic SSD saturates within 20 Myr. We show that the SSD is quite insensitive to the presence of an LSD and is even stronger in the presence of a large-scale shear flow. The LSD grows more slowly in the presence of SSD, saturating after 5 Gyr versus 1–2 Gyr in studies where the SSD is weak or absent. The LSD primarily grows in warm gas in the galactic midplane. Saturation of the LSD occurs due toα-quenching near the midplane as the growing mean-field produces a magneticαthat opposes the kineticα. The magnetic energy in our models of the LSD shows a slightly sublinear response to increasing resolution, indicating that we are converging toward the physical solution at 1 pc resolution. Clustering supernovae in OB associations increases the growth rates for both the SSD and the LSD, compared to a horizontally uniform supernova distribution.

     
    more » « less
  2. Abstract

    Fingering convection (also known as thermohaline convection) is a process that drives the vertical transport of chemical elements in regions of stellar radiative zones where the mean molecular weight increases with radius. Recently, Harrington & Garaud used three-dimensional direct numerical simulations (DNS) to show that a vertical magnetic field can dramatically enhance the rate of chemical mixing by fingering convection. Furthermore, they proposed a so-called “parasitic saturation” theory to model this process. Here, we test their model over a broad range of parameter space, using a suite of DNS of magnetized fingering convection, varying the magnetic Prandtl number, magnetic field strength, and composition gradient. We find that the rate of chemical mixing measured in the simulations is not always predicted accurately by their existing model, in particular when the magnetic diffusivity is large. We then present an extension of the Harrington & Garaud model which resolves this issue. When applied to stellar parameters, it recovers the results of Harrington & Garaud except in the limit where fingering convection becomes marginally stable, where the new model is preferred. We discuss the implications of our findings for stellar structure and evolution.

     
    more » « less
  3. ABSTRACT

    Convection in massive main-sequence stars generates large-scale magnetic fields in their cores that persists as they evolve up the red giant branch. The remnants of these fields may take the form of the Prendergast magnetic field, a combination of poloidal and toroidal field components that are expected to stabilize each other. Previous analytic and numerical calculations did not find any evidence for instability of the Prendergast field over short time-scales. In this paper, we present numerical simulations which show a long time-scale, linear instability of this magnetic field. We find the instability to be robust to changes in boundary conditions and it is not stabilized by strong stable stratification. The instability is a resistive instability, and the growth rate has a power-law dependence on the resistivity, in which the growth rate decreases as the resistivity decreases. We estimate the growth rate of the instability in stars by extrapolating this power law to stellar values of the resistivity. The instability is sufficiently rapid to destabilize the magnetic field on time-scales shorter than the stellar evolution time-scale, indicating that the Prendergast field is not a good model to use in studies of magnetic fields in stars.

     
    more » « less
  4. Constraining dynamo theories of magnetic field origin by observation is indispensable but challenging, in part because the basic quantities measured by observers and predicted by modelers are different. We clarify these differences and sketch out ways to bridge the divide. Based on archival and previously unpublished data, we then compile various important properties of galactic magnetic fields for nearby spiral galaxies. We consistently compute strengths of total, ordered, and regular fields, pitch angles of ordered and regular fields, and we summarize the present knowledge on azimuthal modes, field parities, and the properties of non-axisymmetric spiral features called magnetic arms. We review related aspects of dynamo theory, with a focus on mean-field models and their predictions for large-scale magnetic fields in galactic discs and halos. Furthermore, we measure the velocity dispersion of H i gas in arm and inter-arm regions in three galaxies, M 51, M 74, and NGC 6946, since spiral modulation of the root-mean-square turbulent speed has been proposed as a driver of non-axisymmetry in large-scale dynamos. We find no evidence for such a modulation and place upper limits on its strength, helping to narrow down the list of mechanisms to explain magnetic arms. Successes and remaining challenges of dynamo models with respect to explaining observations are briefly summarized, and possible strategies are suggested. With new instruments like the Square Kilometre Array (SKA), large data sets of magnetic and non-magnetic properties from thousands of galaxies will become available, to be compared with theory. 
    more » « less
  5. Visco-resistive magnetohydrodynamic turbulence, driven by a two-dimensional unstable shear layer that is maintained by an imposed body force, is examined by decomposing it into dissipationless linear eigenmodes of the initial profiles. The down-gradient momentum flux, as expected, originates from the large-scale instability. However, continual up-gradient momentum transport by large-scale linearly stable but nonlinearly excited eigenmodes is identified, and found to nearly cancel the down-gradient transport by unstable modes. The stable modes effectuate this by depleting the large-scale turbulent fluctuations via energy transfer to the mean flow. This establishes a physical mechanism underlying the long-known observation that coherent vortices formed from nonlinear saturation of the instability reduce turbulent transport and fluctuations, as such vortices are composed of both the stable and unstable modes, which are nearly equal in their amplitudes. The impact of magnetic fields on the nonlinearly excited stable modes is then quantified. Even when imposing a strong magnetic field that almost completely suppresses the instability, the up-gradient transport by the stable modes is at least two-thirds of the down-gradient transport by the unstable modes, whereas for weaker fields, this fraction reaches up to 98% . These effects are persistent with variations in magnetic Prandtl number and forcing strength. Finally, continuum modes are shown to be energetically less important, but essential for capturing the magnetic fluctuations and Maxwell stress. A simple analytical scaling law is derived for their saturated turbulent amplitudes. It predicts the fall-off rate as the inverse of the Fourier wavenumber, a property which is confirmed in numerical simulations. 
    more » « less