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.

arXiv comment24 pages, 16 figures. Accepted to MNRAS

References in corpus (42)

• 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
• Hot Accretion Flows Around Black Holes
• First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring
• Efficient Generation of Jets from Magnetically Arrested Accretion on a Rapidly Spinning Black Hole
• General Relativistic Magnetohydrodynamic Simulations of Magnetically Choked Accretion Flows around Black Holes
• HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
• Magnetically Arrested Disk: An Energetically Efficient Accretion Flow
• First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon
• General Relativistic Magnetohydrodynamic Simulations of Jet Formation and Large-Scale Propagation from Black Hole Accretion Systems
• A Measurement of the Electromagnetic Luminosity of a Kerr Black Hole
• 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
• Magnetically Driven Accretion Flows in the Kerr Metric I: Models and Overall Structure
• The Structure and Dynamics of the Sub-parsec Scale Jet in M87 Based on 50 VLBA Observations Over 17 Years at 43 GHz
• The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
• Numerical simulations of super-critical black hole accretion flows in general relativity
• General relativistic magnetohydrodynamical simulations of the jet in M87
• Primitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics
• The Black Hole Accretion Code
• Efficiency of Magnetic to Kinetic Energy Conversion in a Monopole Magnetosphere
• A public code for general relativistic, polarised radiative transfer around spinning black holes
• 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
• Two-temperature, Magnetically Arrested Disc simulations of the jet from the supermassive black hole in M87
• Parabolic Jets from the Spinning Black Hole in M87
• On the Properties of Time-Dependent, Force-Free, Degenerate Electrodynamics
• 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
• Radiative, two-temperature simulations of low luminosity black hole accretion flows in general relativity
• Horizon-Scale Lepton Acceleration in Jets: Explaining the Compact Radio Emission in M87
• 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
• Evolving Nonthermal Electrons in Simulations of Black Hole Accretion
• A new scheme for matching general relativistic ideal magnetohydrodynamics to its force-free limit
• Pair Drizzle around Sub-Eddington Supermassive Black Holes
• GiRaFFE: An Open-Source General Relativistic Force-Free Electrodynamics Code
• 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