Local stability of differential rotation in magnetised radiation zones and the solar tachocline

Robert W. Dymott,Adrian J. Barker,Chris A. Jones,Steven M. Tobias
2024-10-09
Abstract:We study local magnetohydrodynamical (MHD) instabilities of differential rotation in magnetised, stably-stratified regions of stars and planets using a Cartesian Boussinesq model. We consider arbitrary latitudes and general shears (with gravity direction misaligned from this by an angle $\phi$), to model radial ($\phi=0$), latitudinal ($\phi=\pm 90^\circ$), and mixed differential rotations, and study both non-diffusive (including magnetorotational, MRI, and Solberg-Høiland instabilities) and diffusive instabilities (including Goldreich-Schubert-Fricke, GSF, and MRI with diffusion). These instabilities could drive turbulent transport and mixing in radiative regions, including the solar tachocline and the cores of red giant stars, but their dynamics are incompletely understood. We revisit linear axisymmetric instabilities with and without diffusion and analyse their properties in the presence of magnetic fields, including deriving stability criteria and computing growth rates, wavevectors and energetics, both analytically and numerically. We present a more comprehensive analysis of axisymmetric local instabilities than prior work, exploring arbitrary differential rotations and diffusive processes. The presence of a magnetic field leads to stability criteria depending upon angular velocity rather than angular momentum gradients. We find MRI operates for much weaker differential rotations than the hydrodynamic GSF instability, and that it typically prefers much larger lengthscales, while the GSF instability is impeded by realistic strength magnetic fields. We anticipate MRI to be more important for turbulent transport in the solar tachocline than the GSF instability when $\phi>0$ in the northern (and vice versa in the southern) hemisphere, though the latter could operate just below the convection zone when MRI is absent for $\phi<0$.
Solar and Stellar Astrophysics,Earth and Planetary Astrophysics
What problem does this paper attempt to address?
This paper aims to study the local magnetohydrodynamic (MHD) instabilities of differential rotation in magnetized and stably stratified regions. Specifically, the authors use the Cartesian Boussinesq model to study these instabilities under arbitrary latitudes and general shear conditions in stars and planets, including the misalignment between the direction of gravity and the latitude angle φ. The research covers non - diffusive types (such as the magnetorotational instability MRI and the Solberg - Høiland instability) and diffusive types (such as the Goldreich - Schubert - Fricke instability GSF and MRI with diffusion). These instabilities may drive turbulent transport and mixing in radiative regions (such as the solar transition layer and the red - giant core), but their dynamics are not fully understood yet. ### Main objectives of the paper: 1. **Understand the linear properties of local instabilities**: Analyze the linear properties of local instabilities of differential rotation in the radiative regions of stars and planets in the presence of a locally uniform magnetic field, especially the axisymmetric linear stability. 2. **Determine the roles of MRI and GSF instabilities**: Study the potential roles of these instabilities in angular - momentum transport, chemical mixing, and dynamo generation. 3. **Analyze instabilities under different parameters**: Explore how the properties of unstable modes depend on the magnetic field strength \( B_0 \) and the magnetic Prandtl number \( P_m=\nu / \eta \) (the ratio of kinematic viscosity \( \nu \) to Ohmic diffusivity \( \eta \)). ### Main methods: - **Local Cartesian model**: Use a local Cartesian coordinate system to study small - scale differential - rotation instabilities. - **Linear stability analysis**: Conduct axisymmetric linear stability analysis, derive stability criteria, and calculate growth rates, wave vectors, and energy. - **Numerical simulation**: Plan to conduct three - dimensional nonlinear numerical simulations to consider the effects of more general magnetic field directions. ### Key findings: - **Operating conditions of MRI**: MRI can operate under differential - rotation conditions much weaker than those of the hydrodynamic GSF instability, and it generally prefers larger scales. - **Influence of GSF instability**: Magnetic fields of realistic strength will impede the GSF instability. - **Predictions for the solar transition layer**: When \( \phi>0 \), MRI is more important than the GSF instability in the solar transition layer in the northern hemisphere; the opposite is true in the southern hemisphere. When \( \phi < 0 \), the GSF instability may operate below the convective zone, while MRI is absent. ### Examples of formulas: - **Square of Alfvén frequency**: \[ \omega_A^2=(B_0\cdot k)^2 = k^2 B_0^2\sin^2\theta_k \] - **Stability criterion**: \[ \Delta = 4a\omega_A^2-4\xi+(a + b)^2=(a + b)^2+16\omega_A^2(\hat{k}\cdot\Omega)^2 \] - **Condition for maximum growth rate**: \[ s^2=\Omega S\sin(\phi)\sin(\Lambda)-B_0^2\omega_A^2\sin^2(\phi) \] Through these studies, the authors hope to better understand the dynamic behavior of differential rotation inside stars, thereby providing a more accurate physical basis for stellar - evolution models.