Thermal instability in X-ray photoionized media in Active Galactic Nuclei: II. Role of the thermal conduction in warm absorber

A photoionized gas under constant pressure can display a thermal instability, with three or more solutions for possible thermal equilibrium. A unique solution of the structure of the irradiated medium is obtained only if electron conduction is considered. The subject of our study is to estimate how the effect of thermal conduction affects the structure and transmitted spectrum of the warm absorber computed by solving radiative transfer with the code TITAN. We developed a new computational mode for the code TITAN to obtain several solutions for a given external conditions and we test a posteriori which solution is the closest one to the required integral condition based on conduction. We demonstrate that the automatic mode of the code TITAN provides the solution to the radiative transfer which is generally consistent with the estimated exact solution within a few per cent accuracy, with larger errors for some line intensities (up to 20 per cent) for iron lines at intermediate ionization state.

💡 Research Summary

The paper investigates the thermal instability that arises in X‑ray photo‑ionized gas under constant pressure, a situation commonly encountered in the warm absorbers of active galactic nuclei (AGN). In such environments the temperature–pressure (T‑P) curve can develop an S‑shaped branch, producing three or more possible thermal equilibrium solutions (cold, intermediate, hot) for the same pressure. Traditional radiative‑transfer calculations, including those performed with the widely used TITAN code, typically select one of these solutions arbitrarily or enforce numerical stability by suppressing the multi‑solution region, thereby ignoring the physical role of electron thermal conduction.

The authors introduce a new computational mode in TITAN that explicitly searches for all admissible equilibrium branches for a given set of external parameters (incident spectrum, gas density, total pressure). For each branch they compute the thickness of the conductive layer and the associated conductive heat flux using a Spitzer‑type conductivity and the local temperature gradient. The key physical constraint is the integral conduction condition: the conductive heat flux must exactly balance the net radiative heating across the layer. By evaluating this condition a posteriori, the authors identify which of the numerically obtained branches corresponds to the physically correct, conduction‑stabilized solution.

A systematic comparison between the “conduction‑selected” solution and the solution automatically produced by TITAN’s standard mode shows that the latter is generally very close to the exact one. The temperature structure differs by less than a few percent in most of the slab, and the intensities of the majority of diagnostic lines (e.g., O VII, O VIII, Ne IX) agree within about 5 %. However, lines arising from iron ions in the intermediate ionization range (Fe XVII–Fe XXIII) exhibit larger discrepancies, up to ~20 % in intensity. This sensitivity stems from the fact that the conductive layer, although extremely thin (∼10⁸–10⁹ cm), can shift the depth at which these iron lines form, thereby altering the local ionization balance and line optical depth.

The study demonstrates that electron thermal conduction, even when confined to a very narrow transition zone, plays a decisive role in selecting a unique thermal equilibrium for constant‑pressure photo‑ionized gas. Ignoring this effect can lead to significant errors in the predicted spectra, especially for temperature‑sensitive features. Nevertheless, the standard TITAN mode, which does not explicitly enforce the conduction condition, still provides a remarkably accurate approximation for most practical purposes, justifying its continued use in large‑scale AGN modeling.

Beyond the immediate technical findings, the authors discuss broader implications. The presence of a thin conductive layer could, in principle, be inferred from high‑resolution X‑ray spectra obtained with upcoming missions such as XRISM and Athena, where subtle changes in iron line ratios might serve as indirect diagnostics of conduction. Future work is suggested to incorporate more sophisticated physics—anisotropic conduction, magnetic field suppression, and time‑dependent hydrodynamics—to refine the model further. Ultimately, the paper establishes that a self‑consistent treatment of radiative transfer and thermal conduction is essential for a reliable interpretation of warm‑absorber spectra and for constraining the physical conditions in the inner regions of AGN.