- The paper introduces a general relativistic framework that self-consistently models energy exchange between a Novikov–Thorne disk and a static slab corona in XRBs.
- It employs an upgraded MONK code to simulate radiative transfer in Kerr spacetime, capturing the effects of Compton scattering and disk reflection.
- The findings highlight that realistic coronae require alternative geometries or energy transport mechanisms to reproduce the hardest observed X-ray spectra.
General Relativistic Disk–Corona Interaction: Static Slab Corona in Global Energy Balance
Introduction
This work presents a comprehensive general relativistic treatment of the energy exchange between an accretion disk and an extended, static slab corona, focusing on the global energy equilibrium of the system. The study is motivated by the need to reconcile observed hard X-ray spectra in X-ray binaries (XRBs) with physically consistent models of disk–corona interaction, particularly in the hard spectral state. The authors extend the Monte Carlo radiative transfer code MONK to self-consistently couple disk and corona emission, including feedback via illumination and reflection, within the Kerr metric. The approach enables the determination of equilibrium solutions for a given set of disk and corona parameters, constrained by the total available accretion power.
The upgraded MONK code models the disk as a Novikov–Thorne thin disk, parameterized by black hole mass, spin, accretion rate, and radial extent. The corona is implemented as a static, co-rotating slab with specified electron temperature, optical depth, and vertical thickness. The code tracks photon geodesics in Kerr spacetime, simulates Compton scattering in the corona, and models disk reflection using Chandrasekhar’s formalism for a semi-infinite electron atmosphere. An albedo parameter controls the fraction of incident flux reflected versus absorbed by the disk.
The equilibrium solution is achieved via iterative adjustment of the fraction $\alpha$ of accretion power dissipated in the disk, with $(1-\alpha)$ heating the corona. The iteration proceeds until the total radiative luminosity matches the available accretion power, $L_{\rm tot} = L_{\rm acc}$, and both $L_{\rm tot}$ and $\alpha$ converge to within 1%. The code accounts for feedback: absorbed illumination increases the disk temperature, modifying the seed photon distribution for Comptonization, which in turn affects coronal output.
Results: Spectral Properties and Physical Constraints
Equilibrium Spectra and Component Fractions
The equilibrium solution yields a self-consistent partitioning of luminosity between disk, corona, and reflection components. For typical hard state parameters (e.g., $T_e = 120$ keV, albedo = 0.5, $M = 10 M_\odot$, $a = 0$), the lowest achievable photon index for the Comptonized spectrum is $\Gamma \approx 1.7$–$1.8$. Increasing the coronal optical depth hardens the spectrum, but equilibrium cannot be maintained for $\tau \gtrsim 0.16$; beyond this, the feedback loop destabilizes, and no physical solution exists. The disk fraction in the total luminosity decreases with increasing optical depth, while the coronal and reflection fractions increase, with reflection saturating at high $\tau$ due to photon trapping.
Dependence on Spin, Temperature, and Albedo
Higher black hole spin and higher coronal temperature both yield harder spectra at fixed optical depth, primarily by increasing disk illumination and reducing the required $\alpha$ for equilibrium. However, these effects do not significantly lower the minimum photon index. In contrast, increasing the disk albedo (i.e., ionization) allows for much harder spectra and lower $\Gamma_{\rm min}$, as the reprocessed soft photon flux is suppressed. For albedo = 1, the system can reach $\Gamma \sim 1.0$ at high optical depth, but this scenario is physically extreme and inconsistent with observed reflection features.
Anisotropy and Escaping Fractions
The study quantifies the escaping fraction $f_e$ of disk photons and the anisotropic fraction $f_a$ of coronal luminosity illuminating the disk. $f_e$ is systematically lower than the naive $\exp(-\tau)$ expectation due to geometric and angular effects. $f_a$ increases with optical depth at low $\tau$ due to backward scattering, then decreases as multiple scatterings isotropize the emission. These diagnostics are essential for understanding the angular distribution of energy and the feedback loop.
Local vs. Global Energy Balance
The model imposes global energy equilibrium but does not enforce local balance at each radius. Analysis of the radial profiles of coronal cooling and available heating power reveals that, for a uniform slab corona, energy must flow from the inner to the outer regions to maintain local equilibrium. This suggests that realistic coronae should have radially decreasing temperature or require energy transport mechanisms, such as magnetic reconnection, to redistribute power.
Implications and Limitations
Constraints on Coronal Geometry
The results demonstrate that a static slab corona fully covering the disk cannot reproduce the hardest observed spectra in XRBs (i.e., $\Gamma < 1.7$) for plausible disk albedo and coronal temperatures. This finding is robust and consistent with previous studies in Newtonian and local equilibrium frameworks. Harder spectra require either extreme disk ionization (which suppresses reflection features) or alternative geometries, such as truncated disks or outflowing coronae.
Reflection Modeling and Atomic Features
The current reflection treatment neglects atomic transitions, limiting the ability to fit observed reflection spectra. Incorporating external reflection tables (e.g., XILLVER, reflionx) is nontrivial due to their assumptions about illumination, lack of intrinsic disk dissipation, and computational challenges in achieving convergence. Future work should address these limitations to enable direct comparison with observational data.
Model Simplifications and Future Directions
The model assumes a homogeneous corona and uniform energy partitioning, which are idealizations. Realistic coronae likely exhibit radial gradients in temperature and optical depth, and the disk ionization should decrease with radius. Improvements should include energy-dependent albedo, hydrodynamic disk modeling, and exploration of alternative geometries. The upgraded MONK code provides a foundation for such studies, including polarization diagnostics.
Conclusion
This study establishes a general relativistic, self-consistent framework for modeling disk–corona interaction in accreting black hole systems, with a focus on static slab coronae. The main findings are:
- Static slab coronae cannot explain the hardest XRB spectra for realistic disk albedo and coronal temperatures.
- Higher spin and temperature harden the spectrum but do not significantly lower $\Gamma_{\rm min}$.
- Increased disk albedo allows for harder spectra but at the expense of reflection features.
- Local energy balance requires either radial temperature gradients or energy transport within the corona.
- The framework enables constraints on coronal geometry using spectral and polarization data.
Further development is needed to incorporate atomic reflection features and to model more complex geometries, which are essential for interpreting the full range of observed XRB states.