skip to main content

Title: Magnetically modified spherical accretion in GRMHD: reconnection-driven convection and jet propagation
ABSTRACT We present 3D general relativistic magnetohydrodynamic simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider several angles between the magnetic field direction and the black hole spin. In the resulting flows, the mid-plane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ∼10–200 per cent occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ∼30○ with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40○–80○ with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or more » M87. « less
 ;  ;  ;  
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Sponsoring Org:
National Science Foundation
More Like this

    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,more »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.

    « less

    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.more »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.

    « less

    We present the results of nine simulations of radiatively inefficient magnetically arrested discs (MADs) across different values of the black hole spin parameter a*: −0.9, −0.7, −0.5, −0.3, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $t \gtrsim 100\, 000\, GM/c^3$ to ensure disc inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and consequently also jet power, of MAD discs depends strongly on the black hole spin, confirming previous results. Prograde discs saturate at a much higher relative magnetic flux and have more powerful jets than their retrograde counterparts. MADs with spinning black holes naturally launch jets with generalized parabolic profiles whose widths vary as a power of distance from the black hole. For distances up to 100GM/c2, the power-law index is k ≈ 0.27–0.42. There is a strong correlation between the disc–jet geometry and the dimensionless magnetic flux, resulting in prograde systems displaying thinner equatorial accretion flows near the black hole and wider jets, compared to retrograde systems. Prograde and retrograde MADs also exhibit different trends in disc variability: accretion rate variability increases with increasing spin for a* > 0 and remains almost constant for a* ≲ 0,more »while magnetic flux variability shows the opposite trend. Jets in the MAD state remove more angular momentum from black holes than is accreted, effectively spinning down the black hole. If powerful jets from MAD systems in Nature are persistent, this loss of angular momentum will notably reduce the black hole spin over cosmic time.

    « less

    Magnetically arrested accretion discs (MADs) around black holes (BHs) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here, we propose instead that it is mainly a dynamically important toroidal field – the result of dynamo action triggered by the significant but still relatively weak vertical field – that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotational instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (standard and normal evolution; strongly magnetized but not ‘arrested’) discs, and suggest amore »set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field.

    « less
  5. Abstract

    Using a combination of general-relativistic magnetohydrodynamics simulations and ray tracing of synchrotron emission, we study the effect of modest (24°) misalignment between the black hole spin and plasma angular momentum, focusing on the variability of total flux, image centroids, and image sizes. We consider both millimeter and infrared (IR) observables motivated by Sagittarius A* (Sgr A*), though our results apply more generally to optically thin flows. For most quantities, tilted accretion is more variable, primarily due to a significantly hotter and densercoronalregion well off the disk midplane. We find (1) a 150% increase in millimeter light-curve variability when adding tilt to the flow; (2) the tilted image centroid in the millimeter shifts on a scale of 3.7μas over 28 hr (5000 gravitational times) for some electron temperature models; (3) tilted disk image diameters in the millimeter can be 10% larger (52 versus 47μas) than those of aligned disks at certain viewing angles; (4) the tilted models produce significant IR flux, similar to that seen in Sgr A*, with comparable or even greater variability than observed; and (5) for some electron models, the tilted IR centroid moves by more than 50μas over several hours, in a similar fashion to themore »centroid motion detected by the GRAVITY interferometer.

    « less