skip to main content

Title: Global Electron Thermodynamics in Radiatively Inefficient Accretion Flows

In the collisionless plasmas of radiatively inefficient accretion flows, heating and acceleration of ions and electrons are not well understood. Recent studies in the gyrokinetic limit revealed the importance of incorporating both the compressive and Alfvénic cascades when calculating the partition of dissipated energy between the plasma species. In this paper, we use a covariant analytic model of the accretion flow to explore the impact of compressive and Alfvénic heating, Coulomb collisions, compressional heating, and radiative cooling on the radial temperature profiles of ions and electrons. We show that, independent of the partition of heat between the plasma species, even a small fraction of turbulent energy dissipated to the electrons makes their temperature scale with a virial profile and the ion-to-electron temperature ratio smaller than in the case of pure Coulomb heating. In contrast, the presence of compressive cascades makes this ratio larger because compressive turbulent energy is channeled primarily into the ions. We calculate the ion-to-electron temperature in the inner accretion flow for a broad range of plasma properties, mass accretion rates, and black hole spins and show that it ranges between 5 ≲Ti/Te≲ 40. We provide a physically motivated expression for this ratio that can be used to calculate observables from simulations of black hole accretion flows for a wide range of conditions.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 47
["Article No. 47"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The centers of our Galaxy and the nearby Messier 87 are known to contain supermassive black holes, which support accretion flows that radiate across the electromagnetic spectrum. Although the composition of the accreting gas is unknown, it is likely a mix of ionized hydrogen and helium. We use a simple analytic model and a suite of numerical general relativistic magnetohydrodynamic accretion simulations to study how polarimetric images and spectral energy distributions of the source are influenced by the hydrogen/helium content of the accreting matter. We aim to identify general trends rather than make quantitatively precise predictions, since it is not possible to fully explore the parameter space of accretion models. If the ion-to-electron temperature ratio is fixed, then increasing the helium fraction increases the gas temperature; to match the observational flux density constraints, the number density of electrons and magnetic field strengths must therefore decrease. In our numerical simulations, emission shifts from regions of low to high plasmaβ—both altering the morphology of the image and decreasing the variability of the light curve—especially in strongly magnetized models with emission close to the midplane. In polarized images, we find that the model gas composition influences the degree to which linear polarization is (de)scrambled and therefore affects estimates for the resolved linear polarization fraction. We also find that the spectra of helium-composition flows peak at higher frequencies and exhibit higher luminosities. We conclude that gas composition may play an important role in predictive models for black hole accretion.

    more » « less

    We perform 2D particle-in-cell simulations of magnetic reconnection in electron-ion plasmas subject to strong Compton cooling and calculate the X-ray spectra produced by this process. The simulations are performed for trans-relativistic reconnection with magnetization 1 ≤ σ ≤ 3 (defined as the ratio of magnetic tension to plasma rest-mass energy density), which is expected in the coronae of accretion discs around black holes. We find that magnetic dissipation proceeds with inefficient energy exchange between the heated ions and the Compton-cooled electrons. As a result, most electrons are kept at a low temperature in Compton equilibrium with radiation, and so thermal Comptonization cannot reach photon energies $\sim 100\,$ keV observed from accreting black holes. Nevertheless, magnetic reconnection efficiently generates $\sim 100\,$ keV photons because of mildly relativistic bulk motions of the plasmoid chain formed in the reconnection layer. Comptonization by the plasmoid motions dominates the radiative output and controls the peak of the radiation spectrum Epk. We find Epk ∼ 40 keV for σ = 1 and Epk ∼ 100 keV for σ = 3. In addition to the X-ray peak around 100 keV, the simulations show a non-thermal MeV tail emitted by a non-thermal electron population generated near X-points of the reconnection layer. The results are consistent with the typical hard state of accreting black holes. In particular, we find that the spectrum of Cygnus X-1 is well explained by electron-ion reconnection with σ ∼ 3.

    more » « less
  3. Abstract

    We demonstrate using linear theory and particle-in-cell (PIC) simulations that a synchrotron-cooling collisionless plasma acquires pressure anisotropy and, if the plasma beta is sufficiently high, becomes unstable to the firehose instability, in a process that we dub the synchrotron firehose instability (SFHI). The SFHI channels free energy from the pressure anisotropy of the radiating, relativistic electrons (and/or positrons) into small-amplitude, kinetic-scale, magnetic-field fluctuations, which pitch-angle scatter the particles and bring the plasma to a near-thermal state of marginal instability. The PIC simulations reveal a nonlinear cyclic evolution of firehose bursts interspersed by periods of stable cooling. We compare the SFHI for electron–positron and electron–ion plasmas. As a byproduct of the growing electron-firehose magnetic-field fluctuations, magnetized ions gain a pressure anisotropy opposite to that of the electrons. If these ions are relativistically hot, we find that they also experience cooling due to collisionless thermal coupling with the electrons, which we argue is mediated by a secondary ion-cyclotron instability. We suggest that the SFHI may be activated in a number of astrophysical scenarios, such as within ejecta from black hole accretion flows and relativistic jets, where the redistribution of energetic electrons from low to high pitch angles may cause transient bursts of radiation.

    more » « less
  4. Abstract

    The temperatures of the heavy ions (Ti) in the solar corona provide critical information about the heating mechanism of the million-degree corona. However, the measurement ofTiis usually challenging due to the nonthermal motion, instrumental limitations, and optically thin nature of the coronal emissions. We present the measurement ofTiand its dependency on the ion charge-to-mass ratio (Z/A) at the polar coronal hole boundary, only assuming that heavy ions have the same nonthermal velocity. To improve theZ/Acoverage and study the influence of the instrumental broadening, we used a coordinated observation from the EUV Imaging Spectrometer on board the Hinode satellite and the Solar Ultraviolet Measurements of Emitted Radiation (SUMER) on board the Solar and Heliospheric Observatory. We found that theTiof ions withZ/Aless than 0.20 or greater than 0.33 are much higher than the local electron temperature. We ran the Alfvén Wave Solar Model-realtime to investigate the formation of optically thin emissions along the line of sight (LOS). The simulation suggested that plasma bulk motions along the LOS broaden the widths of hot emission lines in the coronal hole (e.g., Fexii, Fexiii). We discussed other factors that might affect theTimeasurement, including the non-Gaussian wings in some bright SUMER lines, which can be fitted by a double-Gaussian or aκdistribution. Our study confirms the preferential heating of heavy ions in coronal holes and provides new constraints on coronal heating models.

    more » « less
  5. Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and$P_{1/3}$multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number$M_A \approx 1.5$), supersonic (Sonic Mach number$M_S \approx 4-5$) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio$L/\delta \approx 100$, Lundquist number$S_L \approx 400$). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate =$\tau _{\text {cool}}^{-1}/\tau _{H}^{-1}\approx 100$); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

    more » « less