skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Toward General Relativistic Magnetohydrodynamics Simulations in Stationary Nonvacuum Spacetimes
Abstract Accretion of magnetized gas on compact astrophysical objects such as black holes (BHs) has been successfully modeled using general relativistic magnetohydrodynamic (GRMHD) simulations. These simulations have largely been performed in the Kerr metric, which describes the spacetime of a vacuum and stationary spinning BH in general relativity (GR). The simulations have revealed important clues to the physics of accretion flows and jets near the BH event horizon and have been used to interpret recent Event Horizon Telescope images of the supermassive BHs M87* and Sgr A*. The GRMHD simulations require the spacetime metric to be given in horizon-penetrating coordinates such that all metric coefficients are regular at the event horizon. Only a few metrics, notably the Kerr metric and its electrically charged spinning analog, the Kerr–Newman metric, are currently available in such coordinates. We report here horizon-penetrating forms of a large class of stationary, axisymmetric, spinning metrics. These can be used to carry out GRMHD simulations of accretion on spinning, nonvacuum BHs and non-BHs within GR, as well as accretion on spinning objects described by non-GR metric theories of gravity.  more » « less
Award ID(s):
1743747
PAR ID:
10468230
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
956
Issue:
1
ISSN:
2041-8205
Format(s):
Medium: X Size: Article No. L11
Size(s):
Article No. L11
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract General relativity, as a diffeomorphism-invariant theory, allows the description of physical phenomena in a wide variety of coordinate systems. In the presence of boundaries, such as event horizons and null infinity, time coordinates must be carefully adapted to the global causal structure of spacetime to ensure a computationally efficient description. Horizon-penetrating time is used to describe the dynamics of infalling matter and radiation across the event horizon, while hyperboloidal time is used to study the propagation of radiation toward the idealized observer at null infinity. In this paper, we explore the historical and mathematical connection between horizon-penetrating and hyperboloidal time coordinates, arguing that both classes of coordinates are simply regular choices of time across null horizons. We review the height-function formalism in stationary spacetimes, providing examples that may be useful in computations, such as source-adapted foliations or Fefferman–Graham–Bondi coordinates near null infinity. We discuss bridges connecting the boundaries of spacetime through a time hypersurface across null horizons, including the event horizon, null infinity, and the cosmological horizon. 
    more » « less
  2. In the past few years, the Event Horizon Telescope (EHT) has provided the first-ever event horizon-scale images of the supermassive black holes (BHs) M87* and Sagittarius A* (Sgr A*). The next-generation EHT project is an extension of the EHT array that promises larger angular resolution and higher sensitivity to the dim, extended flux around the central ring-like structure, possibly connecting the accretion flow and the jet. The ngEHT Analysis Challenges aim to understand the science extractability from synthetic images and movies to inform the ngEHT array design and analysis algorithm development. In this work, we compare the accretion flow structure and dynamics in numerical fluid simulations that specifically target M87* and Sgr A*, and were used to construct the source models in the challenge set. We consider (1) a steady-state axisymmetric radiatively inefficient accretion flow model with a time-dependent shearing hotspot, (2) two time-dependent single fluid general relativistic magnetohydrodynamic (GRMHD) simulations from the H-AMR code, (3) a two-temperature GRMHD simulation from the BHAC code, and (4) a two-temperature radiative GRMHD simulation from the KORAL code. We find that the different models exhibit remarkably similar temporal and spatial properties, except for the electron temperature, since radiative losses substantially cool down electrons near the BH and the jet sheath, signaling the importance of radiative cooling even for slowly accreting BHs such as M87*. We restrict ourselves to standard torus accretion flows, and leave larger explorations of alternate accretion models to future work. 
    more » « less
  3. Abstract We investigate the Teukolsky equation in horizon-penetrating coordinates to study the behavior of perturbation waves crossing the outer horizon. For this purpose, we use the null ingoing/outgoing Eddington–Finkelstein coordinates. The first derivative of the radial equation is a Fuchsian differential equation with an additional regular singularity to the ones the radial one has. The radial functions satisfy the physical boundary conditions without imposing any regularity conditions. We also observe that the Hertz-Weyl scalar equations preserve their angular and radial signatures in these coordinates. Using the angular equation, we construct the metric perturbation for a circularly orbiting perturber around a black hole in Kerr spacetime in a horizon-penetrating setting. Furthermore, we completed the missing metric pieces due to the massMand angular momentumJperturbations. We also provide an explicit formula for the metric perturbation as a function of the radial part, its derivative, and the angular part of the solution to the Teukolsky equation. Finally, we discuss the importance of the extra singularity in the radial derivative for the convergence of the metric expansion. 
    more » « less
  4. The Event Horizon Telescope (EHT) is taking the first images of black holes resolved at horizon scales to measure their shadows and probe accretion physics. A promising avenue for testing the hypothesis that astrophysical black holes are described by the Kerr solution to Einstein’s equations is to compare the size and shape of the shadow a black hole casts on the surrounding emission to the predictions of the Kerr metric. We develop here an efficient parametric framework to perform this test. We carry out ray-tracing simulations for several parametrized non-Kerr metrics to create a large data set of non-Kerr shadows that probe the allowed parameter space for the free parameters of each metric. We then perform principal components analysis (PCA) on this set of shadows and show that only a small number of components are needed to accurately reconstruct all shadows within the set. We further show that the amplitude of the PCA components are smoothly related to the free parameters in the metrics, and therefore, that these PCA components can be fit to EHT observations in order to place constraints on the free parameters of these metrics that will help quantify any potential deviations from the Kerr solution. 
    more » « less
  5. Abstract Relativistic jets from a Kerr black hole (BH) following the core collapse of a massive star (“collapsar”) is a leading model for gamma-ray bursts (GRBs). However, the two key ingredients for a Blandford–Znajek-powered jet—rapid rotation and a strong magnetic field—seem mutually exclusive. Strong fields in the progenitor star’s core transport angular momentum outward more quickly, slowing down the core before collapse. Through innovative multidisciplinary modeling, we first use MESA stellar evolution models followed to core collapse to explicitly show that the small length scale of the instabilities—likely responsible for angular momentum transport in the core (e.g., Tayler–Spruit)—results in a lownetmagnetic flux fed to the BH horizon, far too small to power GRB jets. Instead, we propose a novel scenario in which collapsar BHs acquire their magnetic “hair” from their progenitor proto–neutron star (PNS), which is likely highly magnetized from an internal dynamo. We evaluate the conditions for the BH accretion disk to pin the PNS magnetosphere to its horizon immediately after the collapse. Our results show that the PNS spin-down energy released before collapse matches the kinetic energy of Type Ic-BL supernovae, while the nascent BH’s spin and magnetic flux produce jets consistent with observed GRB characteristics. We map our MESA models to 3D general-relativistic magnetohydrodynamic simulations and confirm that accretion disks confine the strong magnetic flux initiated near a rotating BH, enabling the launch of successful GRB jets, whereas a slower-spinning BH or one without a disk fails to do so. 
    more » « less