Engineering Papers⌕ Search

Engineering topics

Gammie, Charles

Publications and source records attributed to Gammie, Charles.

Probing plasma physics with spectral index maps of accreting black holes on event horizon scales

The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. Here, in this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code IPOLE to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically show that the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximized near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disc (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma.

(magnetohydrodynamics) MHD↗

Two-temperature GRRMHD Simulations of M87

Here, we present axisymmetric two-temperature general relativistic radiation magnetohydrodynamic simulations of the inner region of the accretion flow onto the supermassive black hole M87. We address uncertainties from previous modeling efforts through inclusion of models for (1) self-consistent dissipative and Coulomb electron heating (2) radiation transport (3) frequency-dependent synchrotron emission, self-absorption, and Compton scattering. We adopt a distance D = 16.7 Mpc, an observer angle θ = 20°, and consider black hole masses M/M ⊙ = (3.3 x 10 9 , 6.2 x 10 9 ) and spins a star = (0.5, 0.9375) in a four-simulation suite. For each (M, a star ), we identify the accretion rate that recovers the 230 GHz flux from very long baseline interferometry measurements. We report on disk thermodynamics at these accretion rates (M · /M · Edd ~ 10 -5 ). The disk remains geometrically thick; cooling does not lead to a thin disk component. While electron heating is dominated by Coulomb rather than dissipation for r ≳ 10GM/c 2 , the accretion disk remains two-temperature. Radiative cooling of electrons is not negligible, especially for r ≲ 10GM/c 2 . The Compton y parameter is of order unity. We then compare derived and observed or inferred spectra, millimeter images, and jet powers. Simulations with M/M ⊙ = 3.3 × 10 9 are in conflict with observations. These simulations produce millimeter images that are too small, while the low-spin simulation also overproduces X-rays. For M/M ⊙ = 6.2 x 10 9 , both simulations agree with constraints on radio/IR/X-ray fluxes and millimeter image sizes. Simulation jet power is a factor 10 2 –10 3 below inferred values, a possible consequence of the modest net magnetic flux in our models.

79 ASTRONOMY AND ASTROPHYSICS↗

The Radiative Efficiency and Spectra of Slowly Accreting Black Holes from Two-temperature GRRMHD Simulations

In this paper, we present axisymmetric numerical simulations of radiatively inefficient accretion flows onto black holes combining general relativity, magnetohydrodynamics, self-consistent electron thermodynamics, and frequency-dependent radiation transport. We investigate a range of accretion rates up to ${10}^{-5}\,{\dot{M}}_{\mathrm{Edd}}$ onto a ${10}^{8}\,{M}_{\odot }$ black hole with spin ${a}_{\star }=0.5$. We report on averaged flow thermodynamics as a function of accretion rate. We present the spectra of outgoing radiation and find that it varies strongly with accretion rate, from synchrotron-dominated in the radio at low $\dot{M}$ to inverse-Compton-dominated at our highest $\dot{M}$. In contrast to canonical analytic models, we find that by $\dot{M}\approx {10}^{-5}\,{\dot{M}}_{\mathrm{Edd}}$, the flow approaches $\sim 1 \% $ radiative efficiency, with much of the radiation due to inverse-Compton scattering off Coulomb-heated electrons far from the black hole. Finally, these results have broad implications for modeling of accreting black holes across a large fraction of the accretion rates realized in observed systems.

79 ASTRONOMY AND ASTROPHYSICS↗

THE LONG-TERM EVOLUTION OF PHOTOEVAPORATING PROTOPLANETARY DISKS

We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting for {approx}< 20% of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that the initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the M-dot-R{sub wall} plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of the order of 10 AU.

79 ASTRONOMY AND ASTROPHYSICS↗

VARIABLE ACCRETION OUTBURSTS IN PROTOSTELLAR EVOLUTION

We extend the one-dimensional, two-zone models of long-term protostellar disk evolution with infall of Zhu et al. to consider the potential effects of a finite viscosity in regions where the ionization is too low for the magnetorotational instability (MRI) to operate (the {sup d}ead zone{sup )}. We find that the presence of a small but finite dead zone viscosity, as suggested by simulations of stratified disks with MRI-active outer layers, can trigger inside-out bursts of accretion, starting at or near the inner edge of the disk, instead of the previously found outside-in bursts with zero dead zone viscosity, which originate at a few AU in radius. These inside-out bursts of accretion bear a qualitative resemblance to the outburst behavior of one FU Ori object, V1515 Cyg, in contrast to the outside-in burst models, which more closely resemble the accretion events in FU Ori and V1057 Cyg. Our results suggest that the type and frequency of outbursts are potentially a probe of transport efficiency in the dead zone. Simulations must treat the inner disk regions, R {approx}< 0.5 AU, to show the detailed time evolution of accretion outbursts in general and to observe the inside-out bursts in particular.

79 ASTRONOMY AND ASTROPHYSICS↗