Inhomogeneous accretion discs and the soft states of black hole X-ray binaries

Observations of black hole binaries (BHBs) have established a rich phenomenology of X-ray states. The soft states range from the low variability, accretion disc dominated thermal state (TD) to the higher variability, non-thermal steep power law state (SPL). The disc component in all states is typically modeled with standard thin disc accretion theory. However, this theory is inconsistent with optical/UV spectral, variability, and gravitational microlensing observations of active galactic nuclei (AGNs), the supermassive analogs of BHBs. An inhomogeneous disc (ID) model with large (~0.4 dex) temperature fluctuations in each radial annulus can qualitatively explain all of these AGN observations. The inhomogeneity may be a consequence of instabilities in radiation dominated discs, and therefore may be present in BHBs as well. We show that ID models can explain many features of the TD and SPL states of BHBs. The observed relationships between spectral hardness, disc fraction, and rms variability amplitude in BHBs are reproduced with temperature fluctuations similar to those inferred in AGNs, suggesting a unified picture of luminous accretion discs across orders of magnitude in black hole mass. This picture can be tested with spectral fitting of ID models, X-ray polarization observations, and radiation MHD simulations. If BHB accretion discs are indeed inhomogeneous, only the most disc dominated states (disc fraction > 0.95) can be used to robustly infer black hole spin using current continuum fitting methods.

💡 Research Summary

The paper tackles a long‑standing inconsistency between the standard thin‑disc paradigm and a variety of observations of luminous accretion flows, both in active galactic nuclei (AGN) and in stellar‑mass black‑hole X‑ray binaries (BHBs). While the thin‑disc model successfully predicts the overall shape of the thermal component in the soft states of BHBs, it fails to reproduce the broad optical/UV spectra, high variability, and microlensing size measurements seen in AGN. The authors propose an “inhomogeneous disc” (ID) model in which each radial annulus of the disc possesses a log‑normal temperature distribution. The mean temperature follows the classic T(r)∝r⁻³⁄⁴ law, but the standard deviation of the logarithmic temperature, σ_T, is large—about 0.4 dex—reflecting strong, stochastic fluctuations that may arise from radiation‑pressure‑driven thermal instabilities or MRI‑driven turbulence in radiation‑dominated zones.

First, the authors demonstrate that the ID model reproduces the key AGN observables: a broader-than‑expected UV bump, enhanced short‑timescale variability, and the larger effective sizes inferred from gravitational microlensing. They then apply the same statistical temperature prescription to BHBs, constructing synthetic X‑ray spectra by Monte‑Carlo sampling of the temperature field at each radius and integrating the resulting black‑body emission over the disc surface. By varying σ_T they generate a family of spectra that span the observed phenomenology from the low‑variability, disc‑dominated Thermal Dominant (TD) state to the high‑variability, steep‑power‑law (SPL) state.

The results are strikingly consistent with observations. When σ_T≈0.1 dex, the disc contributes >95 % of the total 2–20 keV flux, the rms variability is ≤2 %, and the spectrum is soft, matching the classic TD state. Increasing σ_T to 0.3–0.4 dex reduces the disc fraction to 70–80 %, raises the rms to 5–10 %, and hardens the spectrum, reproducing the SPL state’s characteristic steep power‑law tail. Moreover, the model naturally reproduces the empirical correlation between spectral hardness, disc fraction, and rms variability that has been documented across many BHB outbursts. The same σ_T≈0.4 dex that works for AGN also explains the BHB data, suggesting that radiation‑pressure‑dominated discs develop comparable temperature fluctuations over a vast range of black‑hole masses (10 M⊙ to 10⁹ M⊙).

A crucial implication concerns black‑hole spin measurements using the continuum‑fitting method. This technique assumes a perfectly smooth, thin disc whose inner edge coincides with the innermost stable circular orbit. The ID model shows that even modest temperature inhomogeneities (σ_T≈0.05–0.1 dex) can artificially harden the high‑energy tail of the spectrum, leading to systematic overestimates of spin unless the disc fraction exceeds ~0.95. Consequently, only the most disc‑dominated observations should be employed for reliable spin determinations.

The authors outline three avenues for testing the ID hypothesis. First, X‑ray polarimetry can probe the anisotropic emission patterns expected from a patchy temperature field, as fluctuations should imprint characteristic signatures on the polarization degree and angle. Second, state‑of‑the‑art radiation‑magnetohydrodynamic (RMHD) simulations can directly measure the statistical properties of temperature fluctuations (σ_T, spatial correlation length) in radiation‑pressure‑dominated discs, providing a first‑principles check on the phenomenological model. Third, systematic spectral fitting of BHB data with ID‑based templates can assess whether the improved fits justify the added complexity over standard thin‑disc plus power‑law models.

In summary, the paper presents a unified, mass‑independent framework for luminous accretion discs that simultaneously accounts for the soft‑state phenomenology of BHBs and the puzzling spectral and variability properties of AGN. By attributing the observed diversity to stochastic temperature fluctuations rather than to separate physical components, the ID model challenges conventional interpretations of the SPL state and calls for a reassessment of spin measurements that rely on the assumption of a perfectly homogeneous disc. Future polarimetric observations and high‑resolution RMHD simulations will be decisive in confirming or refuting this compelling picture.