An official website of the United States government
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 M87.
more »
« less
The Role of Strong Magnetic Fields in Stabilizing Highly Luminous Thin Disks
Abstract We present and analyze a set of three-dimensional, global, general relativistic radiation magnetohydrodynamic simulations of thin, radiation-pressure-dominated accretion disks surrounding a nonrotating, stellar-mass black hole. The simulations are initialized using the Shakura–Sunyaev model with a mass accretion rate of M ̇ = 3 L Edd / c 2 (corresponding to L = 0.17 L Edd ). Our previous work demonstrated that such disks are thermally unstable when accretion is driven by an α -viscosity. In the present work, we test the hypothesis that strong magnetic fields can both drive accretion through magnetorotational instability and restore stability to such disks. We test four initial magnetic field configurations: (1) a zero-net-flux case with a single, radially extended set of magnetic field loops (dipole), (2) a zero-net-flux case with two radially extended sets of magnetic field loops of opposite polarity stacked vertically (quadrupole), (3) a zero-net-flux case with multiple radially concentric rings of alternating polarity (multiloop), and (4) a net-flux, vertical magnetic field configuration (vertical). In all cases, the fields are initially weak, with a gas-to-magnetic pressure ratio ≳100. Based on the results of these simulations, we find that the dipole and multiloop configurations remain thermally unstable like their α -viscosity counterpart, in our case collapsing vertically on the local thermal timescale and never fully recovering. The vertical case, on the other hand, stabilizes and remains so for the duration of our tests (many thermal timescales). The quadrupole case is intermediate, showing signs of both stability and instability. The key stabilizing factor is the ability of specific field configurations to build up and sustain strong, P mag ≳ 0.5 P tot , toroidal fields near the midplane of the disk. We discuss the reasons why certain configurations are able to do this effectively and others are not. We then compare our stable simulations to the standard Shakura–Sunyaev disk.
more »
« less
- Award ID(s):
- 1907850
- PAR ID:
- 10450380
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 939
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 31
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Previously we demonstrated that the magnetorotational instability (MRI) grows vigorously in eccentric disks, much as it does in circular disks, and we investigated the nonlinear development of the eccentric MRI without vertical gravity. Here we explore how vertical gravity influences the magnetohydrodynamic (MHD) turbulence stirred by the eccentric MRI. Similar to eccentric disks without vertical gravity, the ratio of Maxwell stress to pressure, or the Shakura–Sunyaevαparameter, remains ∼10−2, and the local sign flip in the Maxwell stress persists. Vertical gravity also introduces two new effects. Strong vertical compression near pericenter amplifies reconnection and dissipation, weakening the magnetic field. Angular momentum transport by MHD stresses broadens the mass distribution over eccentricity at much faster rates than without vertical gravity; as a result, spatial distributions of mass and eccentricity can be substantially modified in just ∼5 to 10 orbits. MHD stresses in the eccentric debris of tidal disruption events may power emission ≳1 yr after disruption.more » « less
-
null (Ed.)ABSTRACT We present a systematic shearing-box investigation of magnetorotational instability (MRI)-driven turbulence in a weakly collisional plasma by including the effects of an anisotropic pressure stress, i.e. anisotropic (Braginskii) viscosity. We constrain the pressure anisotropy (Δp) to lie within the stability bounds that would be otherwise imposed by kinetic microinstabilities. We explore a broad region of parameter space by considering different Reynolds numbers and magnetic-field configurations, including net vertical flux, net toroidal-vertical flux, and zero net flux. Remarkably, we find that the level of turbulence and angular-momentum transport are not greatly affected by large anisotropic viscosities: the Maxwell and Reynolds stresses do not differ much from the MHD result. Angular-momentum transport in Braginskii MHD still depends strongly on isotropic dissipation, e.g. the isotropic magnetic Prandtl number, even when the anisotropic viscosity is orders of magnitude larger than the isotropic diffusivities. Braginskii viscosity nevertheless changes the flow structure, rearranging the turbulence to largely counter the parallel rate of strain from the background shear. We also show that the volume-averaged pressure anisotropy and anisotropic viscous transport decrease with increasing isotropic Reynolds number (Re); e.g. in simulations with net vertical field, the ratio of anisotropic to Maxwell stress (αA/αM) decreases from ∼0.5 to ∼0.1 as we move from Re ∼ 103 to Re ∼ 104, while 〈4$$\pi$$Δp/B2〉 → 0. Anisotropic transport may thus become negligible at high Re. Anisotropic viscosity nevertheless becomes the dominant source of heating at large Re, accounting for $${\gtrsim } 50 {{\ \rm per\ cent}}$$ of the plasma heating. We conclude by briefly discussing the implications of our results for radiatively inefficient accretion flows on to black holes.more » « less
-
null (Ed.)ABSTRACT FU Ori is the prototype of FU Orionis systems that are outbursting protoplanetary discs. Magnetic fields in FU Ori’s accretion discs have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori’s inner accretion disc. We find that (1) when the disc is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z ∼ R, similar to the ‘surface accretion’ mechanism in previous locally isothermal MHD simulations. (2) A moderate disc wind is launched in the vertical field simulations with a terminal speed of ∼300–500 km s−1 and a mass-loss rate of 1–10 per cent the disc accretion rate, which is consistent with observations. Disc wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disc photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disc wind, the disc photosphere temperature is lower than that predicted by the thin disc theory, and the previously inferred disc accretion rate may be lower than the real accretion rate by a factor of ∼2–3.more » « less
-
Abstract The matter in an accretion disk must lose angular momentum when moving radially inwards but how this works has long been a mystery. By calculating the trajectories of individual colliding neutrals, ions, and electrons in a weakly ionized 2D plasma containing gravitational and magnetic fields, we numerically simulate accretion disk dynamics at the particle level. As predicted by Lagrangian mechanics, the fundamental conserved global quantity is the total canonical angular momentum, not the ordinary angular momentum. When the Kepler angular velocity and the magnetic field have opposite polarity, collisions between neutrals and charged particles cause: (i) ions to move radially inwards, (ii) electrons to move radially outwards, (iii) neutrals to lose ordinary angular momentum, and (iv) charged particles to gain canonical angular momentum. Neutrals thus spiral inward due to their decrease of ordinary angular momentum while the accumulation of ions at small radius and accumulation of electrons at large radius produces a radially outward electric field. In 3D, this radial electric field would drive an out-of-plane poloidal current that produces the magnetic forces that drive bidirectional astrophysical jets. Because this neutral angular momentum loss depends only on neutrals colliding with charged particles, it should be ubiquitous. Quantitative scaling of the model using plausible disk density, temperature, and magnetic field strength gives an accretion rate of 3 × 10−8solar mass per year, which is in good agreement with observed accretion rates.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/1538-4357/ac938b
