Resolution analysis of magnetically arrested disc simulations

Abstract

ABSTRACT Polarization measurements by the Event Horizon Telescope from M87* and Sgr A* suggest that there is a dynamically strong, ordered magnetic field, typical of what is expected of a magnetically arrested accretion disc (MAD). In such discs, the strong poloidal magnetic field can suppress the accretion flow and cause episodic flux eruptions. Recent work shows that general relativistic magnetohydrodynamic (GRMHD) MAD simulations feature dynamics of turbulence and mixing instabilities that are becoming resolved at higher resolutions. We perform a convergence study of MAD states exceeding the status quo by an order of magnitude in resolution. We use existing 3D simulations performed with the H-AMR code, up to a resolution of 5376 × 2304 × 2304 in a logarithmic spherical-polar grid. We find consistent time-averaged disc properties across all resolutions. However, higher resolutions reveal signs of inward angular momentum transport attributed to turbulent convection, particularly evident when mixing instabilities occur at the surfaces of flux tubes during flux eruptions. Additionally, we see wave-like features in the jet sheath, which become more prominent at higher resolutions, that may induce mixing between the jet and disc. At higher resolutions, we observe the sheath to be thinner, resulting in increased temperature, reduced magnetization, and greater variability. Those differences could affect the dissipation of energy, which would eventually result in distinct observable radiative emission from high-resolution simulations. With higher resolutions, we can delve into crucial questions about horizon-scale physics and its impact on the dynamics and emission properties of larger-scale jets.