skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The DOI auto-population feature in the Public Access Repository (PAR) will be unavailable from 4:00 PM ET on Tuesday, July 8 until 4:00 PM ET on Wednesday, July 9 due to scheduled maintenance. We apologize for the inconvenience caused.


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
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 General relativistic magnetohydrodynamics (GRMHD) simulations are an indispensable tool in studying accretion onto compact objects. The Event Horizon Telescope (EHT) frequently uses libraries of ideal GRMHD simulations to interpret polarimetric, event-horizon-scale observations of supermassive black holes at the centers of galaxies. In this work, we present a library of 10 nonradiative, ideal GRMHD simulations that were utilized by the EHT Collaboration in their analysis of Sagittarius A*. The parameter survey explores both low (SANE) and high (MAD) magnetization states across five black hole spinsa* = −15/16, −1/2, 0, +1/2, +15/16 where each simulation was run out to 30,000GM/c−3. We find the angular momentum and energy flux in SANE simulations closely matches the thin-disk value, with minor deviations in prograde models due to fluid forces. This leads to spin equilibrium arounda* ∼ 0.94, consistent with previous studies. We study the flow of conserved quantities in our simulations and find mass, angular momentum, and energy transport in SANE accretion flows to be primarily inward and fluid dominated. MAD models produce powerful jets with outflow efficiency >1 fora* = + 0.94, leading to black hole spin-down in prograde cases. We observe outward directed energy and angular momentum fluxes on the horizon, as expected for the Blandford–Znajek mechanism. MAD accretion flows are sub-Keplerian and exhibit greater variability than their SANE counterpart. They are also hotter than SANE disks withinr ≲ 10GM/c−2. This study is accompanied by a public release of simulation data athttp://thz.astro.illinois.edu/. 
    more » « less
  2. 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
  3. 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
  4. UV and optical continuum reverberation mapping is a powerful tool for probing the accretion disk and inner broad-line region. However, recent reverberation mapping campaigns in the X-ray, UV, and optical have found lags consistently longer than those expected from the standard disk reprocessing picture. The largest discrepancy to date was recently reported in Mrk 335, where UV/optical lags are up to 12 times longer than expected. Here, we perform a frequency-resolved time lag analysis of Mrk 335, using Gaussian processes to account for irregular sampling. For the first time, we compare the Fourier frequency-resolved lags directly to those computed using the popular interpolated cross-correlation function method applied to both the original and detrended light curves. We show that the anticipated disk reverberation lags are recovered by the Fourier lags when zeroing in on the short-timescale variability. This suggests that a separate variability component is present on long timescales. If this separate component is modeled as reverberation from another region beyond the accretion disk, we constrain a size scale of roughly 15 lt-days from the central black hole. This is consistent with the size of the broad-line region inferred from Hβreverberation lags. We also find tentative evidence for a soft X-ray lag, which we propose may be due to light travel time delays between the hard X-ray corona and distant photoionized gas that dominates the soft X-ray spectrum below 2 keV. 
    more » « less
  5. ABSTRACT Turbulent high-energy astrophysical systems often feature asymmetric energy injection: for instance, Alfvén waves propagating from an accretion disc into its corona. Such systems are ‘imbalanced’: the energy fluxes parallel and antiparallel to the large-scale magnetic field are unequal. In the past, numerical studies of imbalanced turbulence have focused on the magnetohydrodynamic regime. In this study, we investigate externally driven imbalanced turbulence in a collision-less, ultrarelativistically hot, magnetized pair plasma using 3D particle-in-cell (PIC) simulations. We find that the injected electromagnetic momentum efficiently converts into plasma momentum, resulting in net motion along the background magnetic field with speeds up to a significant fraction of lightspeed. This discovery has important implications for the launching of accretion disc winds. We also find that although particle acceleration in imbalanced turbulence operates on a slower time-scale than in balanced turbulence, it ultimately produces a power-law energy distribution similar to balanced turbulence. Our results have ramifications for black hole accretion disc coronae, winds, and jets. 
    more » « less