Hybrid GRMHD and Force-Free Simulations of Black Hole Accretion

Abstract

We present a new approach for stably evolving general relativistic magnetohydrodynamic (GRMHD) simulations in regions where the magnetization $\sigma =b^2/\rho c^2$ becomes large. GRMHD codes typically struggle to evolve plasma above $\sigma \approx 100$ in simulations of black hole accretion. To ensure stability, GRMHD codes will inject mass density artificially to the simulation as necessary to keep the magnetization below a ceiling value ${\sigma }_{\max }$. We propose an alternative approach where the simulation transitions to solving the equations of general relativistic force-free electrodynamics (GRFFE) above a magnetization ${\sigma }_{\mathrm{trans}}$. We augment the GRFFE equations in the highly magnetized region with approximate equations to evolve the decoupled field-parallel velocity, plasma energy density, and plasma mass density. Our hybrid scheme is explicit and easily added to the framework of standard-volume GRMHD codes. We present a variety of tests of our method, implemented in the GRMHD code KORAL, and we show first results from a 3D hybrid GRMHD+GRFFE simulation of a magnetically arrested disc (MAD) around a spinning black hole. Our hybrid MAD simulation closely matches the average properties of a standard GRMHD MAD simulation with the same initial conditions in low magnetization regions, but it achieves a magnetization $\sigma \approx 10^6$ in the evacuated jet funnel. We present simulated horizon-scale images of both simulations at 230 GHz with the black hole mass and accretion rate matched to M87*. Images from the hybrid simulation are less affected by the choice of magnetization cutoff ${\sigma }_{\mathrm{cut}}$ imposed in radiative transfer than images from the standard GRMHD simulation.

24 pages, 16 figures. Accepted to MNRAS

References in corpus (18)

• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole
• First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole
• First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring
• First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon
• First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole
• First Sagittarius A* Event Horizon Telescope Results. III: Imaging of the Galactic Center Supermassive Black Hole
• Observing the Inner Shadow of a Black Hole: A Direct View of the Event Horizon
• Global simulations of axisymmetric radiative black hole accretion disks in general relativity with a sub-grid magnetic dynamo
• Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes
• A ring-like accretion structure in M87 connecting its black hole and jet
• The Disc-Jet Symbiosis Emerges: Modeling the Emission of Sagittarius A* with Electron Thermodynamics
• First M87 Event Horizon Telescope Results. IX. Detection of Near-horizon Circular Polarization
• Kinematics of the M87 jet in the collimation zone: gradual acceleration and velocity stratification
• General-Relativistic Simulations of Four States of Accretion onto Millisecond Pulsars
• Pair Drizzle around Sub-Eddington Supermassive Black Holes
• Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration
• Impacts of nonthermal emission on the images of a black hole shadow and extended jets in two-temperature GRMHD simulations

Cited by in corpus (7)

• Survey of Radiative, Two-Temperature Magnetically Arrested Simulations of the Black Hole M87* I: Turbulent Electron Heating
• Signatures of Black Hole Spin and Plasma Acceleration in Jet Polarimetry
• Discovery of Limb Brightening in the Parsec-scale Jet of NGC 315 through Global Very Long Baseline Interferometry Observations and Its Implications for Jet Models
• GR-RMHD Simulations of Super-Eddington Accretion Flows onto a Neutron Star with Dipole and Quadrupole Magnetic Fields
• Non-thermal Synchrotron Emission and Polarization Signatures during Black Hole Flux Eruptions
• Jet Archaeology and Forecasting: Image Variability and Magnetic Field Configuration
• Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts