MRI turbulence in accretion discs at large magnetic Prandtl numbers

TL;DR

This study explores how large magnetic Prandtl numbers influence MRI-driven turbulence in accretion disks, revealing a scaling relationship and a plateau effect at very high Pm, with implications for astrophysical environments.

Contribution

It provides the first detailed analysis of MRI turbulence at large Pm, identifying a Pm-dependent scaling of turbulent transport and a saturation effect at very high Pm.

Findings

Turbulent transport scales with Pm as α ∼ Pm^δ, with δ ≈ 0.5-0.7.

At very large Pm, turbulent energy and stress plateau, losing dependence on Pm.

Scaling of α-Pm is affected by simulation box aspect ratio and background shear.

Abstract

The effect of large magnetic Prandtl number $\text{Pm}$ (the ratio of viscosity to resistivity) on the turbulent transport and energetics of the magnetorotational instability (MRI) is poorly understood, despite the realization of this regime in astrophysical environments as disparate as discs from binary neutron star mergers, the inner regions of low mass X-ray binaries and active galactic nuclei, and the interiors of protoneutron stars. We investigate the MRI dynamo and associated turbulence in the regime $\text{Pm}>1$ by carrying out fully compressible, 3D MHD shearing box simulations using the finite-volume code \textsc{PLUTO}, focusing mostly on the case of Keplerian shear relevant to accretion discs. We find that when the magnetic Reynolds number is kept fixed, the turbulent transport (as parameterized by $\alpha$, the ratio of stress to thermal pressure) scales with the magnetic Prandtl number as $\alpha \sim \text{Pm}^{\delta}$, with $\delta \sim 0.5-0.7$ up to $\text{Pm} \sim 128$. However, this scaling weakens as the magnetic Reynolds number is increased. Importantly, compared to previous studies, we find a new effect at very large $\text{Pm}$ -- the turbulent energy and stress begin to plateau, no longer depending on ${\rm Pm}$. To understand these results we have carried out a detailed analysis of the turbulent dynamics in Fourier space, focusing on the effect of increasing $\text{Pm}$ on the transverse cascade -- a key non-linear process induced by the disc shear flow that is responsible for the sustenance of MRI turbulence. Finally, we find that $\alpha$-$\text{Pm}$ scaling is sensitive to the box vertical-to-radial aspect ratio, as well as to the background shear.