skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Large-scale poloidal magnetic field dynamo leads to powerful jets in GRMHD simulations of black hole accretion with toroidal field
Title: Large-scale poloidal magnetic field dynamo leads to powerful jets in GRMHD simulations of black hole accretion with toroidal field
ABSTRACT Accreting black holes (BHs) launch relativistic collimated jets, across many decades in luminosity and mass, suggesting the jet launching mechanism is universal, robust, and scale-free. Theoretical models and general relativistic magnetohydrodynamic (GRMHD) simulations indicate that the key jet-making ingredient is large-scale poloidal magnetic flux. However, its origin is uncertain, and it is unknown if it can be generated in situ or dragged inward from the ambient medium. Here, we use the GPU-accelerated GRMHD code h-amr to study global 3D BH accretion at unusually high resolutions more typical of local shearing box simulations. We demonstrate that turbulence in a radially extended accretion disc can generate large-scale poloidal magnetic flux in situ, even when starting from a purely toroidal magnetic field. The flux accumulates around the BH till it becomes dynamically important, leads to a magnetically arrested disc (MAD), and launches relativistic jets that are more powerful than the accretion flow. The jet power exceeds that of previous GRMHD toroidal field simulations by a factor of 10 000. The jets do not show significant kink or pinch instabilities, accelerate to γ ∼ 10 over three decades in distance, and follow a collimation profile similar to the observed M87 jet.  more » « less
Award ID(s):
1815304
NSF-PAR ID:
10189310
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
494
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3656 to 3662
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

     
    more » « less
  2. ; ; ; ; ( , The Astrophysical Journal)
    Abstract

    Models invoking magnetic reconnection as the particle acceleration mechanism within relativistic jets often adopt a gradual energy dissipation profile within the jet. However, such a profile has yet to be reproduced in first-principles simulations. Here we perform a suite of 3D general relativistic magnetohydrodynamic simulations of post–neutron star merger disks with an initially purely toroidal magnetic field. We explore the variations in both the microphysics (e.g., nuclear recombination, neutrino emission) and system parameters (e.g, disk mass). In all of our simulations, we find the formation of magnetically striped jets. The stripes result from the reversals in the poloidal magnetic flux polarity generated in the accretion disk. The simulations display large variations in the distributions of stripe duration,τ, and power, 〈PΦ〉. We find that more massive disks produce more powerful stripes, the most powerful of which reaches 〈PΦ〉 ∼ 1049erg s−1atτ∼ 20 ms. The power and variability that result from the magnetic reconnection of the stripes agree with those inferred in short-duration gamma-ray bursts. We find that the dissipation profile of the cumulative energy is roughly a power law in both radial distance,z, andτ, with a slope in the range of ∼1.7–3; more massive disks display larger slopes.

     
    more » « less
  3. ; ; ; ; ; ; ; ; ( , The Astrophysical Journal Letters)
    Abstract X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up ≲10% of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole (SMBH), SMBH merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the SMBH sphere of influence of radius R B to the SMBH of gravitational radius R g at the largest scale separation, R B / R g = 10 3 , to date. Initially, our axisymmetric system of constant-density hot gas contains a weak vertical magnetic field and rotates in the equatorial plane of a rapidly spinning SMBH. We seed the gas with small-scale 2% level pressure perturbations. Infalling gas forms an accretion disk, and the SMBH launches relativistically magnetized collimated jets reaching well outside R B . Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the SMBH, becomes dynamically important, and leads to a magnetically arrested disk state. The SMBH accretes at 2% of the Bondi rate ( M ̇ ≃ 2.4 × 10 − 3 M ⊙ yr − 1 for M87*) and launches twin jets at η = 150% efficiency. These jets are powerful enough ( P jets ≃ 2 × 10 44 erg s −1 ) to escape along the SMBH spin axis and end the short-lived intermittent jet state, whose transient nature can account for the rarity of XRGs. 
    more » « less
  4. ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    Wind-fed models offer a unique way to form predictive models of the accretion flow surrounding Sagittarius A*. We present 3D wind-fed magnetohydrodynamic (MHD) and general relativistic magnetohydrodynamic (GRMHD) simulations spanning the entire dynamic range of accretion from parsec scales to the event horizon. We expand on previous work by including non-zero black hole spin and dynamically evolved electron thermodynamics. Initial conditions for these simulations are generated from simulations of the observed Wolf–Rayet stellar winds in the Galactic Centre. The resulting flow tends to be highly magnetized (β ≈ 2) with an ∼r−1 density profile independent of the strength of magnetic fields in the winds. Our simulations reach the magnetically arrested disc (MAD) state for some, but not all cases. In tilted flows, standard and normal evolution (SANE) jets tend to align with the angular momentum of the gas at large scales, even if that direction is perpendicular to the black hole spin axis. Conversely, MAD jets tend to align with the black hole spin axis. The gas angular momentum shows similar behaviour: SANE flows tend to only partially align while MAD flows tend to fully align. With a limited number of dynamical free parameters, our models can produce accretion rates, 230 GHz flux, and unresolved linear polarization fractions roughly consistent with observations for several choices of electron heating fraction. Absent another source of large-scale magnetic field, winds with a higher degree of magnetization (e.g. where the magnetic pressure is 1/100 of the ram pressure in the winds) may be required to get a sufficiently large rotation measure with consistent sign.

     
    more » « less
  5. ; ; ; ; ; ( , The Astrophysical Journal)
    Abstract Astrophysical jets, launched from the immediate vicinity of accreting black holes, carry away large amounts of power in a form of bulk kinetic energy of jet particles and electromagnetic flux. Here we consider a simple analytical model for relativistic jets at larger distances from their launching sites, assuming a cylindrical axisymmetric geometry with a radial velocity shear, and purely toroidal magnetic field. We argue that as long as the jet plasma is in magnetohydrostatic equilibrium, such outflows tend to be particle dominated, i.e., the ratio of the electromagnetic to particle energy flux, integrated over the jet cross-sectional area, is typically below unity, σ < 1. At the same time, for particular magnetic and radial velocity profiles, magnetic pressure may still dominate over particle pressure for certain ranges of the jet radius, i.e., the local jet plasma parameter β pl < 1, and this may be relevant in the context of particle acceleration and production of high-energy emission in such systems. The jet magnetization parameter can be elevated up to the modest values of σ ≲  ( 10 ) only in the case of extreme gradients or discontinuities in the gaseous pressure, and a significantly suppressed velocity shear. Such configurations, which consist of a narrow, unmagnetized jet spine surrounded by an extended, force-free layer, may require an additional poloidal field component to stabilize them against current-driven oscillations, but even this will not substantially elevate their σ parameter. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email