Radiation GRMHD Simulations of the Hard State of Black Hole X-ray Binaries and the Collapse of a Hot Accretion Flow

We present global radiation GRMHD simulations of strongly magnetized accretion onto a spinning, stellar mass black hole at sub-Eddington rates. Using a frequency-dependent Monte Carlo procedure for Compton scattering, we self-consistently evolve a two-temperature description of the ion-electron fluid and its radiation field. For an Eddington ratio $L/L_{rm Edd} gtrsim 10^{-3}$, the emergent spectrum forms an apparent power law shape from thermal Comptonization up to a cutoff at $simeq 100$ keV, characteristic of that seen in the hard spectral states of black hole X-ray binary systems. At these luminosities, the radiative efficiency is high ($approx 24%$) and results in a denser midplane region where magnetic fields are dynamically important. For $L/L_{rm Edd} sim 10^{-2}$, our hot accretion flow appears to undergo thermal runaway and collapse. Our simulations demonstrate that hot accretion flows can be radiatively efficient and provide an estimate of their maximum luminosity.

Radiation GRMHD simulations of M87: funnel properties and prospects for gap acceleration

We use the public code ebhlight to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion onto the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several sub-grid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) gamma-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $rm sim 10^{41} , erg , s^{-1}$ for MAD models and $rm sim 2times 10^{40} , erg , s^{-1}$ for SANE models, which are comparable to measurements of M87s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behavior. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.

Sgr A* X-ray flares from non-thermal particle acceleration in a magnetically arrested disc

Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of 10-100 for ~60 minutes. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from PIC simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies {gamma_e} <~ 1e6 and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simpified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard (<~ 100 keV) X-ray energies.

QPOs in compact binaries from small scale eruptions in an inner magnetized disk

Dwarf novae (DNe) and low mass X-ray binaries (LMXBs) are compact binaries showing variability on time scales from years to less than seconds. Here, we focus on explaining part of the rapid fluctuations in DNe, following the framework of recent studies on the monthly eruptions of DNe that use a hybrid disk composed of an outer standard disk and an inner magnetized disk. We show that the ionization instability, that is responsible for the monthly eruptions of DNe, is also able to operate in the inner magnetized disk. Given the low density and the fast accretion time scale of the inner magnetized disk, the ionization instability generates small, rapid heating and cooling fronts propagating back and forth in the inner disk. This leads to quasi-periodic oscillations (QPOs) with a period of the order of $1000$ s. A strong prediction of our model is that these QPOs can only develop in quiescence or at the beginning/end of an outburst. We propose that these rapid fluctuations might explain a subclass of already observed QPOs in DNe as well as a, still to observe, subclass of QPOs in LMXBs. We also extrapolate to the possibility that the radiation pressure instability might be related to Type B QPOs in LMXBs.

Infrared interferometry to spatially and spectrally resolve jets in X-ray binaries

Infrared interferometry is a new frontier for precision ground based observing, with new instrumentation achieving milliarcsecond (mas) spatial resolutions for faint sources, along with astrometry on the order of 10 microarcseconds. This technique has already led to breakthroughs in the observations of the supermassive black hole at the Galactic centre and its orbiting stars, AGN, and exo-planets, and can be employed for studying X-ray binaries (XRBs), microquasars in particular. Beyond constraining the orbital parameters of the system using the centroid wobble and spatially resolving jet discrete ejections on mas scales, we also propose a novel method to discern between the various components contributing to the infrared bands: accretion disk, jets and companion star. We demonstrate that the GRAVITY instrument on the Very Large Telescope Interferometer (VLTI) should be able to detect a centroid shift in a number of sources, opening a new avenue of exploration for the myriad of transients expected to be discovered in the coming decade of radio all-sky surveys. We also present the first proof-of-concept GRAVITY observation of a low-mass X-ray binary transient, MAXI J1820+070, to search for extended jets on mas scales. We place the tightest constraints yet via direct imaging on the size of the infrared emitting region of the compact jet in a hard state XRB.

The Effects of Tilt on the Images of Black Hole Accretion Flows

We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets.

Studying Black Holes on Horizon Scales with VLBI Ground Arrays

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

The Sloan Digital Sky Survey Reverberation Mapping Project: Accretion and Broad Emission Line Physics from a Hypervariable Quasar

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical changing look quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

ALMA Polarimetry of Sgr A*: Probing the Accretion Flow from the Event Horizon to the Bondi Radius

Millimeter polarimetry of Sgr A* probes the linearly polarized emission region on a scale of $sim 10$ Schwarzschild radii ($R_S$) as well as the dense, magnetized accretion flow on scales out to the Bondi radius ($sim 10^5 R_S$) through Faraday rotation. We present here multi-epoch ALMA Band 6 (230 GHz) polarimetry of Sgr A*. The results confirm a mean rotation measure, ${rm RM} approx -5 times 10^5 {rm rad m^{-2}}$, consistent with measurements over the past 20 years and support an interpretation of the RM as originating from a radiatively inefficient accretion flow (RIAF) with $dot{M} approx 10^{-8} { rm M_{odot} y^{-1} }$. Variability is observed for the first time in the RM on time scales that range from hours to months. The long-term variations may be the result of changes in the line of sight properties in a turbulent accretion flow. Short-term variations in the apparent RM are not necessarily the result of Faraday rotation and may be the result of complex emission and propagatation effects close to the black hole, some of which have been predicted in numerical modeling. We also confirm the detection of circular polarization at a mean value of $-1.1 pm 0.2 %$. It is variable in amplitude on time scales from hours to months but the handedness remains unchanged from that observed in past centimeter- and millimeter-wavelength detections. These results provide critical constraints for the analysis and interpretation of Event Horizon Telescope data of Sgr A*, M87, and similar sources.

Extreme AGN variability: evidence of magnetically elevated accretion?

Rapid, large amplitude variability at optical to X-ray wavelengths is now seen in an increasing number of Seyfert galaxies and luminous quasars. The variations imply a global change in accretion power, but are too rapid to be communicated by inflow through a standard thin accretion disc. Such discs are long known to have difficulty explaining the observed optical/UV emission from active galactic nuclei. Here we show that alternative models developed to explain these observations have larger scale heights and shorter inflow times. Accretion discs supported by magnetic pressure in particular are geometrically thick at all luminosities, with inflow times as short as the observed few year timescales in extreme variability events to date. Future time-resolved, multi-wavelength observations can distinguish between inflow through a geometrically thick disc as proposed here, and alternative scenarios of extreme reprocessing of a central source or instability-driven limit cycles.