The persistence of magneto-rotational turbulence in gravitationally turbulent accretion disks

Physics Department of Bayreuth, Universitätsstraße 30, Bayreuth, Germany
e-mail: lucas.loehnert@uni-bayreuth.de

Received: 27 July 2022
Accepted: 24 July 2023

Abstract

Aims. Our main goal is to probe the persistence of turbulence originating from the magneto-rotational instability (MRI) in gravito-turbulent disks. This state is referred to here as GI-MRI coexistence, with GI standing for gravitational instability. We test the influence of GI strength, controlled by the cooling law, and the impact of Ohmic resistivity.

Methods. Our starting point was three-dimensional, ideal, magnetohydrodynamic (MHD) simulations of gravitational turbulence in the local shearing-box approximation using the code Athena. We introduced a zero-net-flux magnetic seed field in a GI-turbulent state and investigated the nonlinear evolution. The GI strength was varied by modifying the cooling parameters. We tested the cooling times τ_cΩ₀ = 10, τ_cΩ₀ = 20, and τ_cΩ₀ = 10, with additional background heating. For some resistive cases, ideal-MHD simulations, which had already developed GI-MRI coexistence, were restarted with a finite Ohmic resistivity enabled at the moment of restart.

Results. It appears that there are two possible saturated dynamo states in the ideal-MHD regime: a state of GI-MRI coexistence (for low GI activity) and a strong-GI dynamo. The cases with lower GI activity eventually develop a clearly visible butterfly pattern. For the case with the highest GI activity (τ_cΩ₀ = 10, no heating), a clearly visible butterfly pattern is absent, though more chaotic field reversals are observed above (and below) the mid-plane. We were also able to reproduce the results of previous simulations. With Ohmic resistivity, the simulation outcome can be substantially different. There exists a critical magnetic Reynolds number, ⟨Rm⟩  ∼  500, below which the ideal-MHD outcome is replaced by a new dynamo state. For larger Reynolds numbers, one recovers turbulent states that are more reminiscent of the ideal-MHD states, and especially the strong-GI case. This new state leads to oscillations, which are caused by a significant heat production due to the resistive dissipation of magnetic energy. The additional heat periodically quenches GI, and the quenching events correspond to maxima of the Toomre value, Q.