skip to main content


Title: Neutral-charged-particle Collisions as the Mechanism for Accretion Disk Angular Momentum Transport
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
Award ID(s):
2105492 1914599
NSF-PAR ID:
10367271
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
930
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 167
Size(s):
["Article No. 167"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important large-scale magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows, known as winds,that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular momentum transport is still poorly understood. To probe this question, we analyze a long-duration (1.2 × 105rg/c) simulation of a rapidly rotating (a= 0.9) black hole feeding from a thick (H/r∼ 0.3), adiabatic, magnetically arrested disk (MAD), whose dynamically important magnetic field regulates mass inflow and drives both uncollimated and collimated outflows (i.e., winds and jets, respectively). By carefully disentangling the various angular momentum transport processes within the system, we demonstrate the novel result that disk winds and disk turbulence both extract roughly equal amounts of angular momentum from the disk. We find cumulative angular momentum and mass accretion outflow rates ofL̇r0.9andṀr0.4, respectively. This result suggests that understanding both turbulent and laminar stresses is key to understanding the evolution of systems where geometrically thick MADs can occur, such as the hard state of X-ray binaries, low-luminosity active galactic nuclei, some tidal disruption events, and possibly gamma-ray bursts.

     
    more » « less
  2. Abstract

    Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius,RB= 103Rg. The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outsideRBand suppress BH accretion to 1.5% of the Bondi rate,ṀB. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% ofṀBreach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2%ṀB. Because the BAD and RAD states tangle up the jets and destroy them well insideRB, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies.

     
    more » « less
  3. Abstract

    Absorption-line measurements of the circumgalactic medium (CGM) display a highly nonuniform distribution of lower ionization state species accompanied by more widespread higher ionization state material. This suggests that the CGM is a dynamic, multiphase medium, such as arises in the presence of turbulence. To better understand this evolution, we perform hydrodynamic and magnetohydrodynamic (MHD) simulations of the CGM surrounding Milky Way–like galaxies. In both cases, the CGM is initially in hydrostatic balance in a 1012Mdark matter gravitational potential, and the simulations include rotation in the inner halo and turbulence that decreases radially. They also track ionizations, recombinations, and species-by-species radiative cooling in the presence of the redshift-zero UV background, employing the MAIHEM nonequilibrium chemistry package. We find that after 9 Gyr of evolution, the presence of a magnetic field leads to an overall hotter CGM, with cool gas in the center where magnetic pressure dominates. While the non-MHD run produces more cold clouds overall, we find similar Siiv/Oviand Nv/Oviratios between the MHD and non-MHD runs, which are both very different from their equilibrium values. The non-MHD halo develops cool, low angular momentum filaments above the central disk, in comparison to the MHD run that has more efficient angular momentum transport, especially for the cold gas, which forms a more ordered and extended disk late into its evolution.

     
    more » « less
  4. Abstract

    The behavior of energetic particles in the Earth's dipole and stretched magnetic‐field models is explored with a focus on the shift of the atmospheric loss cone away from the magnetic‐field direction and on non‐adiabatic behavior occurring when the particle gyroradius becomes comparable with the gradient scale length of the local magnetic field. It is shown that the equatorial loss cone is aberrated away from the magnetic‐field direction (pitch angle of 0°) by the perpendicular drift of a charged particle (which is referred to as the “Mozer transform” described in Mozer (1966,https://doi.org/10.1029/JZ071i011p02701). It is found that the Mozer coordinate transformation in pitch‐angle/gyrophase‐angle space better organizes the behavior of bounce times, mirror altitudes, drift speeds, and first adiabatic invariant. It also describes the loss‐cone shift accurately for a certain range of values of the adiabaticity parameterϵ, defined by the ratio of the particle gyroradius to the local radius of curvature of the magnetic field. For particles with largerϵ, the Mozer transform (evaluated with the standard gradient‐curvature drift) breaks down and the “central trajectory” theory can be used to calculate the angular shift of the loss cone, which now includes an Earthward shift. Whenϵbecomes relatively large, stochastic field line curvature (FLC) scattering occurs: we show the intimate connection between the strength of FLC scattering and the loss‐cone shift. By comparing ion orbits in the dipole and Tsyganenko (Ts89 and Ts04) magnetic fields, we conclude that the loss cone of ring‐current ions is significantly modified during geomagnetically active times.

     
    more » « less
  5. 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