skip to main content


Title: Looking for the underlying cause of black hole X-ray variability in GRMHD simulations
ABSTRACT Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a lognormal distribution of flux and a linear rms–flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrinsic to accretion. This variability has been modelled as inward propagating fluctuations of mass accretion rate, although the physical process driving the fluctuations remains puzzling. In this work, we analyse five exceptionally long-duration general relativistic magnetohydrodynamic (GRMHD) simulations of optically thin, geometrically thick, black hole accretion flows to look for hints of propagating fluctuations in the simulation data. We find that the accretion profiles from these simulations do show evidence for inward propagating fluctuations below the viscous frequency by featuring strong radial coherence and positive time lags when comparing smaller to larger radii, although these time lags are generally shorter than the viscous time-scale and are frequency-independent. Our simulations also support the notion that the fluctuations in $\dot{M}$ build up in a multiplicative manner, as the simulations exhibit linear rms–mass flux relations, as well as lognormal distributions of their mass fluxes. When combining the mass fluxes from the simulations with an assumed emissivity profile, we additionally find broad agreement with observed power spectra and time lags, including a recovery of the frequency dependency of the time lags.  more » « less
Award ID(s):
1743747 1907850
NSF-PAR ID:
10188492
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
496
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3808 to 3828
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    The nature and geometry of the accretion flow in the low/hard state of black hole binaries is currently controversial. While most properties are generally explained in the truncated disc/hot inner flow model, the detection of a broad residual around the iron line argues for strong relativistic effects from an untruncated disc. Since spectral fitting alone is somewhat degenerate, we combine it with the additional information in the fast X-ray variability and perform a full spectral-timing analysis for NICER and NuSTAR data on a bright low/hard state of MAXI J1820+070. We model the variability with propagating mass accretion rate fluctuations by combining two separate current insights: that the hot flow is spectrally inhomogeneous, and that there is a discontinuous jump in viscous time-scale between the hot flow and variable disc. Our model naturally gives the double-humped shape of the power spectra, and the increasing high-frequency variability with energy in the second hump. Including reflection and reprocessing from a disc truncated at a few tens of gravitational radii quantitatively reproduces the switch in the lag-frequency spectra, from hard lagging soft at low frequencies (propagation through the variable flow) to the soft lagging hard at the high frequencies (reverberation from the hard X-ray continuum illuminating the disc). The viscous time-scale of the hot flow is derived from the model, and we show how this can be used to observationally test ideas about the origin of the jet.

     
    more » « less
  2. Abstract

    For testing different electron temperature (Te) prescriptions in general relativistic magnetohydrodynamics (GRMHD) simulations through observations, we propose to utilize linear polarization (LP) and circular polarization (CP) images. We calculate the polarization images based on a semi-magnetically arrested disk GRMHD model for variousTeparameters, bearing M87 in mind. We find an LP–CP separation in the images of the low-Tedisk cases at 230GHz; namely, the LP flux mainly originates from downstream of the jet, and the CP flux comes from the counter-side jet, while the total intensity is maximum at the jet base. This can be understood as follows: although the LP flux is generated through synchrotron emission widely around the black hole, most of the LP flux from the jet base does not reach the observer, since it undergoes Faraday rotation (Te2) when passing through the outer cold disk and is thus depolarized. Hence, only the LP flux from the downstream (not passing the cold dense plasmas) can survive. Meanwhile, the CP flux is generated from the LP flux by Faraday conversion ( ∝Te) in the inner hot region. Stronger CP flux is thus observed from the counter-side jet. Moreover, the LP–CP separation is more enhanced at a lower frequency, such as 86 GHz, but is rather weak at 43 GHz, since the media in the latter case is optically thick for synchrotron self-absorption so that all of the fluxes should come from the photosphere. The same is true for cases with higher mass accretion rates and/or larger inclination angles.

     
    more » « less
  3. Abstract

    We report the detection of a long-timescale negative lag, where the blue bands lag the red bands, in the nearby Seyfert 1 galaxy Fairall 9, with two independent methods. Active Galactic Nuclei (AGNs) light curves show variability over a wide range of timescales. By measuring time lags between different wavelengths, the otherwise inaccessible structure and kinematics of the accretion disk can be studied. One common approach, reverberation mapping, quantifies the continuum and line lags moving outward through the disk at the light-travel time, revealing the size and temperature profile of the disk. Inspired by numerical simulations, we expect longer lags to exist in AGN light curves that travel inward on longer timescales, tracing the accretion process itself. By analyzing AGN light curves in both temporal and frequency space, we report the detection of long-timescale lags (∼−70 days) in Fairall 9 that propagate in the opposite direction to the reverberation lag. The short continuum lag (<10 days) is also detected and is consistent with reverberation lags reported in the literature. When fitting the longer lag as a function of frequency with a model motivated by the thin disk model, we find that the disk scale height likely increases outward in the disk. This detection raises the exciting prospect of mapping accretion disk structures across a wide range of AGN parameters.

     
    more » « less
  4. ABSTRACT There exist hitherto unexplained fluctuations in the cosmic infrared background on arcminute scales and larger. These have been shown to cross-correlate with the cosmic X-ray background, leading several authors to attribute the excess to a high-redshift growing black hole population. In order to investigate potential sources that could explain this excess, in this paper, we develop a new framework to compute the power spectrum of undetected sources that do not have constant flux as a function of halo mass. In this formulation, we combine a semi-analytic model for black hole growth and their simulated spectra from hydrodynamical simulations. Revisiting the possible contribution of a high-redshift black hole population, we find that too much black hole growth is required at early epochs for z > 6 accretion to explain these fluctuations. Examining a population of accreting black holes at more moderate redshifts, z ∼ 2–3, we find that such models produce a poor fit to the observed fluctuations while simultaneously overproducing the local black hole mass density. Additionally, we rule out the hypothesis of a missing Galactic foreground of warm dust that produces coherent fluctuations in the X-ray via reflection of Galactic X-ray binary emission. Although we firmly rule out accreting massive black holes as the source of these missing fluctuations, additional studies will be required to determine their origin. 
    more » « less
  5. ABSTRACT

    We present the results of nine simulations of radiatively inefficient magnetically arrested discs (MADs) across different values of the black hole spin parameter a*: −0.9, −0.7, −0.5, −0.3, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $t \gtrsim 100\, 000\, GM/c^3$ to ensure disc inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and consequently also jet power, of MAD discs depends strongly on the black hole spin, confirming previous results. Prograde discs saturate at a much higher relative magnetic flux and have more powerful jets than their retrograde counterparts. MADs with spinning black holes naturally launch jets with generalized parabolic profiles whose widths vary as a power of distance from the black hole. For distances up to 100GM/c2, the power-law index is k ≈ 0.27–0.42. There is a strong correlation between the disc–jet geometry and the dimensionless magnetic flux, resulting in prograde systems displaying thinner equatorial accretion flows near the black hole and wider jets, compared to retrograde systems. Prograde and retrograde MADs also exhibit different trends in disc variability: accretion rate variability increases with increasing spin for a* > 0 and remains almost constant for a* ≲ 0, while magnetic flux variability shows the opposite trend. Jets in the MAD state remove more angular momentum from black holes than is accreted, effectively spinning down the black hole. If powerful jets from MAD systems in Nature are persistent, this loss of angular momentum will notably reduce the black hole spin over cosmic time.

     
    more » « less