The Effects of Complex Accretion Disk Geometry on Broadened Iron K$α$ Lines

Abstract: X-rays are emitted from the corona above the orbiting matter of the accretion disk and travel either directly to us or illuminate the disk. This illumination of the inner disk is enhanced by gravitational light bending, which focuses the rays towards the black hole and therefore towards the inner radii of the disk. These rays that hit the inner radii are reflected back to us, and we observe them in the X-ray reflection spectrum. In this work, we create novel general relativistic ray tracing simulations to investigate the effects of altering the geometry of the accretion disks of black holes on the most dominant part of the reflection spectrum, the iron K$α$ line. Work demonstrating the effect of disk geometry on the iron line has been performed, though many previous analyses have assumed a simplistic system, consisting of a point-source corona with a flat and infinitesimally thin accretion disk. We extend these models to more realistic accretion disk approximations. These include a constant aspect ratio disk, a radiation-pressure-dominated Shakura-Sunyaev disk, an expanded inner disk that has a non-negligible scale height in its inner regions due to radiation pressure, as well as various warped disks. Using measurement uncertainties from XRISM, we find that non-negligible thickness in accretion disks underestimates the black hole spin, corona height, and inclination angle if fitted with a flat disk model. The warped disk model could not be fit with the flat disk approximation.

• The paper demonstrates that complex disk geometries significantly alter iron Kα line profiles, leading to underestimations of black hole spin, corona height, and disk inclination.
• Using full ray tracing in Kerr spacetime, the study quantifies modifications in emissivity profiles and the impacts of self-shadowing, Doppler shifts, and gravitational redshift.
• The findings underscore the need for geometry-sensitive spectral models in high accretion rate regimes to mitigate systematic biases in parameter inference.

Effects of Complex Accretion Disk Geometry on Broadened Iron K$α$ Lines

Introduction

The geometry of accretion disks around black holes critically impacts the interpretation of relativistically broadened iron K$α$ emission lines in X-ray reflection spectra. This paper presents comprehensive general relativistic ray tracing simulations to explore how deviations from the standard infinitesimally thin, flat disk model—commonly used in X-ray spectral fitting—systematically influence iron line profiles, black hole spin estimates, corona height, and inclination angle. The investigations include constant aspect ratio disks, radiation-pressure-dominated Shakura-Sunyaev disks, geometries with expanded or compressed inner regions, and strongly warped disk morphologies. Special attention is given to the implications for parameter inference with future high-resolution instruments such as XRISM.

Methodology: General Relativistic Ray Tracing

The simulations extend the CudaKerr code to account for arbitrary disk geometries, implementing full photon ray tracing in the Kerr spacetime. Both observer-to-disk and corona-to-disk photon trajectories are tracked to compute line transfer functions and photoionization-dependent emissivity maps. The primary illumination is provided by a lamppost (point-source) corona, allowing isolation of geometric effects apart from coronal morphology.

For each geometry, the photon stopping condition is altered to match the disk surface, e.g., by enforcing an inclination-dependent vertical cutoff for thick disks or misalignment at a specified warp radius for warped disks. The approach allows for precise calculations of the gravitational redshift, Doppler shifts, light bending, and self-shadowing effects manifest in the emergent line shapes.

Constant Aspect Ratio and Compressed Disk Geometries

The study first examines the impact of simple geometric thickness, parametrized by a constant scale height-to-radius ratio $\frac{h}{\rho}$. Increasing $\frac{h}{\rho}$ leads to pronounced modifications in the emissivity profiles and line shapes. Disk self-shadowing reduces the flux of the most redshifted photons, especially for configurations with high corona height and large aspect ratio. Enhanced illumination at intermediate radii emerges as the disk intercepts more coronal photons away from the ISCO, altering the classic double-peaked line structure.

Figure 1: Cross-section of the constant aspect ratio accretion disk geometry viewed edge-on ($\theta = \frac{\pi}{2}$).

Figure 2: Emissivity profiles for constant aspect ratio disks showing the dependence on disk thickness and coronal height.

These effects are further modulated in compressed inner disks, where a flat inner region transitions to a thickened structure at a prescribed break radius. The break produces a local enhancement in emissivity and spectral intensity at specific energies, reflecting direct geometric causality in the illumination.

Figure 4: Emissivity profiles for compressed (delayed wedged) accretion disks demonstrate preferential illumination at the transition radius.

Shakura-Sunyaev and Expanded Inner Disks

The simulations extend to the physically motivated Shakura-Sunyaev disk solution, with scale height set by the mass accretion rate and radiative efficiency. For sub-Eddington accretion rates ($\dot{M}/\dot{M}_{\mathrm{Edd}} \lesssim 0.3$), deviations from a flat disk are modest and the flat model remains adequate for fitting. However, at higher rates and in the super-Eddington regime, the scale height becomes significant, and self-obscuration and funneling of coronal photons produce marked changes in line profiles.

Figure 6: Cross-section of the Shakura-Sunyaev accretion disk geometry viewed edge-on; inner radiative-pressure-dominated region is thicker.

Figure 8: Emissivity profiles for the Shakura-Sunyaev accretion disk with varying Eddington ratios. Thick disks at high accretion rates substantially increase illumination of the inner disk.

In super-Eddington disks ($\dot{M}/\dot{M}_{\mathrm{Edd}} \gtrsim 1$), the majority of reflection arises from the thick, inner disk. The result is a strong suppression of blue wing photons and an enhancement in the highly redshifted component, especially at moderate-to-large inclinations—an observational signature directly linked to disk thickness.

Figure 3: Line profiles from a super-Eddington ($\dot{M}/\dot{M}_{\mathrm{Edd}}=17$) Shakura-Sunyaev disk geometry reveal extreme relativistic broadening and suppression of the blue wing.

Warped Disk Geometries

Warped accretion disks, comprising a flat inner disk abruptly tilting into a misaligned outer region, introduce axial asymmetry and highly anisotropic illumination. The azimuthal angle of observation strongly modulates shadowing effects and energy shifts of the observed photons. At specific observer azimuths, the inner disk can fully obscure portions of the outer warped disk, dramatically changing the shape of the K$α$ profile.

Figure 11: Cross-section of the warped accretion disk geometry, highlighting the misalignment $α$ and break radius $α$0.

Figure 5: Schematic of azimuthal dependence and self-shadowing in the warped disk geometry; blue and red shaded regions show areas affected by disk self-obscuration.

Warp angles ($α$1), break radius $α$2, and observer azimuth $α$3 all produce measurable effects, from suppression of one of the Doppler peaks to net broadening and high amplitude centroid shifts. Unlike axisymmetric geometries, warped disks create a diversity of line shapes not reproducible by summing standard thin disk profiles.

Figure 7: Comparison of the flat disk line profile to the family of warped disk profiles for $α$4, $α$5, demonstrating the uniqueness of signatures arising from disk warping.

Biases in Parameter Estimation

Fitting synthetic spectra generated from thick or warped disk models with standard thin disk templates almost invariably leads to underestimation of black hole spin, corona height, and disk inclination. For instance, aspect ratios of $α$6 induce spin underestimations by up to $α$7, corona height by several gravitational radii, and inclination by over $α$8 compared to the actual values. Warped disk profiles could not be fitted at all within the flat disk paradigm at statistically acceptable levels.

Notably, for moderate disk thickness or sub-Eddington rates, the systematic errors remain within or just above typical uncertainties in current observational programs. In the limiting cases of high accretion rates or severely warped disks, the bias exceeds the threshold for secure astrophysical inference, underscoring the imperative for geometry-aware spectral models.

Implications and Future Directions

These findings have significant practical implications for X-ray observatory analysis: the flat, razor-thin disk model is only robust at low accretion rates and in the absence of notable disk warping or thickness. In AGN and X-ray binaries approaching or exceeding the Eddington limit, or exhibiting disk warps (e.g., from Lense-Thirring precession or supernova kicks), geometric complexity is no longer a tolerable subtlety but a dominant source of systematic bias.

Theory-wise, the results call for extended, geometry-inclusive ray tracing and radiative transfer frameworks, particularly as instruments with capabilities comparable to XRISM and Athena come online. These tools must account for self-obscuration, azimuthal asymmetry, and the complex interplay between dynamical and radiative processes in real disks. Furthermore, the study illustrates the need for incorporating wind-driven outflows and returning radiation, which are predicted by recent GRMHD simulations, to achieve a physically self-consistent spectral model for high-$α$9, thick disk regimes.

Conclusion

This investigation demonstrates that the geometry of black hole accretion disks fundamentally alters the observability and interpretation of relativistically broadened iron K$\frac{h}{\rho}$0 lines. The use of general relativistic ray tracing across a broad suite of physically motivated disk models reveals that parameter estimates derived from the canonical thin disk approximation can be significantly biased for thick or warped disks, and that the distinctive spectral signatures of warped geometries cannot be captured within current standard frameworks. These findings have substantial implications for black hole spin measurements, coronal geodesy, and disk structure inference from X-ray reflection spectroscopy in both AGN and X-ray binaries. Progress in both theory and observation will require the adoption of more sophisticated, geometry-sensitive models in the analysis of high-resolution X-ray spectra.