🔭 astrophysics

Unification Model of Active Galactic Nuclei by Photoionization Equilibrium Calculation Based on Radiative Hydrodynamic Simulations

This study utilizes radiative hydrodynamic simulations and photoionization equilibrium calculations to demonstrate that while the Compton-thick covering factor aligns with X-ray observations, the theoretical Compton-thin covering factor matches observations only at high Eddington ratios, suggesting that additional Compton-thin gas extending beyond 10 pc is necessary to explain discrepancies at lower ratios.

Original authors: Atsushi Tanimoto, Keiichi Wada, Yuki Kudoh, Nozomu Kawakatu, Mariko Nomura, Hirokazu Odaka

Published 2026-02-13

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Atsushi Tanimoto, Keiichi Wada, Yuki Kudoh, Nozomu Kawakatu, Mariko Nomura, Hirokazu Odaka

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a Supermassive Black Hole (SMBH) at the center of a galaxy as a giant, blinding spotlight shining out into the darkness. This light is so intense that it can vaporize anything too close to it. Surrounding this spotlight is a swirling cloud of gas and dust, shaped somewhat like a donut or a thick tire. This "donut" is what astronomers call the torus.

The big mystery this paper tries to solve is: How much of the spotlight is hidden by the donut?

If you look at the galaxy from the side, the donut blocks the light, making the galaxy look "obscured" (like looking at a lighthouse through a thick fog). If you look from the top or bottom, you see the light clearly. The "covering factor" is just a fancy way of saying: What percentage of the sky around the black hole is blocked by this dusty donut?

The Problem: The "Blowing Away" Theory vs. Reality

For a long time, scientists thought that if the black hole was shining very brightly (a high "Eddington ratio"), the radiation pressure would be like a giant cosmic hairdryer, blowing the dusty gas away.

Low brightness: The hairdryer is off; the dusty donut stays thick and covers most of the light (80–90% coverage).
High brightness: The hairdryer is on full blast; it blows the dust away, leaving a small hole so you can see the light (30–40% coverage).

However, recent observations showed something weird. Even when the black hole wasn't super bright, the "donut" didn't seem to be as thick as the hairdryer theory predicted. Something else must be happening.

The New Experiment: The "X-Ray Fog"

The authors of this paper decided to run a new simulation. Instead of just looking at the gas as a solid cloud, they treated it like fog that can change its state based on the light hitting it.

The Simulation: They used a supercomputer to simulate how gas moves around a black hole for four different brightness levels.
The Twist (Photoionization): They realized that X-rays from the black hole don't just push gas away; they strip the electrons off atoms.
- Think of the gas atoms as balloons. When the X-ray "sun" hits them, it pops the electron balloons, turning the gas into a hot, transparent plasma (H II).
- Crucial Point: You cannot see this "popped" gas in X-rays because it's invisible to the detectors. Only the "un-popped" neutral gas (H I) can block the light.

What They Found

When they ran the numbers, they discovered a surprising truth about the "donut":

The Inner Core is Always Invisible: Inside a very small distance (about 0.1 light-years) from the black hole, the light is so strong that it instantly "pops" all the electron balloons. The gas becomes transparent plasma. No matter how bright the black hole is, this inner zone is always invisible to X-ray telescopes.
The "Donut" is Actually a Shell: The part of the gas that actually blocks the X-rays (the neutral gas) exists further out.
- The Result: They found that about 30% of the sky is always blocked by this neutral gas, regardless of how bright the black hole is.
- The Match: This 30% number perfectly matches the "Compton-thick" (very heavily obscured) galaxies we see in real life.

The Missing Piece: The "Outer Fence"

Here is the puzzle: Real X-ray observations show that in dimmer galaxies, the covering factor is actually higher (around 80%), not just 30%.

The authors realized their computer simulation was too small. They only looked at the gas within 0.1 light-years.

The Analogy: Imagine you are trying to count how many people are wearing raincoats in a city. You only look at the street right outside your house (0.1 light-years). You see 30% of people wearing raincoats. But if you walk 10 blocks away (10 parsecs), you see a whole new crowd of people in raincoats that you missed.
The Conclusion: To explain why dimmer galaxies look so obscured, there must be a layer of neutral gas farther out (about 10 to 30 light-years away).
- In bright galaxies, the "cosmic hairdryer" (radiation) reaches all the way out there and pops the balloons, clearing the view.
- In dimmer galaxies, the hairdryer is too weak to reach that far, so the "raincoats" (neutral gas) stay intact, blocking the light and making the covering factor look high.

The Takeaway

This paper changes the story of how we see active black holes:

It's not just about wind: The reason we see different amounts of blocked light isn't just because radiation blows dust away. It's mostly because radiation turns the gas invisible (ionizes it).
The "Donut" has two layers: There is a small, inner layer that is always invisible (30% coverage), and a larger, outer layer that only exists in dimmer galaxies.
Unified Model: The "Unified Model" of black holes is still valid, but it needs to include this "photoionization" effect. The black hole doesn't just blow the dust away; it turns the dust into a ghost that X-ray telescopes can't see.

In short: The black hole doesn't just push the fog away; it evaporates it into invisible steam. The amount of fog we see depends on how far out the steam can travel before the light gets too weak to evaporate it.

1. Problem Statement

The Unified Model of Active Galactic Nuclei (AGN) posits that a torus of gas and dust surrounds a supermassive black hole (SMBH), creating an obscuring structure. A critical parameter in this model is the covering factor (the fraction of the central source obscured by this torus), which determines the ratio of unobscured to obscured AGNs.

The Discrepancy: Previous observational studies (e.g., Ricci et al. 2017b) suggested that the covering factor depends heavily on the Eddington ratio ( $\lambda_{Edd}$ ). Specifically, they found that at low Eddington ratios ( $\lambda_{Edd} < 10^{-2}$ ), the Compton-thin covering factor ( $C_{22}$ , gas with column density $10^{22} < N_H < 10^{24}$ cm $^{-2}$ ) is high (80–90%), while at high ratios, it drops to 30–40%. This was attributed to radiation pressure blowing away dusty gas.
The Theoretical Gap: Previous radiative hydrodynamic simulations (e.g., Kudoh et al. 2024) reproduced the torus structure but failed to distinguish between neutral hydrogen (H I) and ionized hydrogen (H II). Since H II does not absorb X-rays via photoelectric absorption, calculating the covering factor based on total gas density leads to inconsistencies when compared to X-ray observations. Furthermore, previous models showed a mismatch where theoretical $C_{22}$ was too high at high Eddington ratios compared to observations.

Objective: To resolve these discrepancies by calculating the covering factors specifically for neutral hydrogen (H I) using photoionization equilibrium calculations, thereby distinguishing between gas that is truly X-ray absorbing and gas that is ionized and transparent.

2. Methodology

The authors employed a two-step hybrid approach combining 2D hydrodynamic simulations with 1D photoionization calculations:

Radiative Hydrodynamic Simulations (RHD):
- Code: CANS+ (Coordinated Astronomical Numerical Software Plus).
- Setup: Two-dimensional axisymmetric simulations of dusty gas around a $10^7 M_\odot$ SMBH.
- Parameters: Simulations were run for four logarithmic Eddington ratios: $\log R_{Edd} = -3, -2, -1, 0$ .
- Physics: Included gravity, radiation pressure (anisotropic on dust), X-ray heating, and viscous forces. The computational domain covered $10^{-4}$ pc to 2 pc.
- Output: Time-averaged 2D density distributions ( $n_H$ ) and velocity fields.
Photoionization Equilibrium Calculations:
- Code: XSTAR.
- Process: The 2D RHD density maps were divided into 9,000 cells (radial and angular). XSTAR was run for each cell to determine the ionization state (H I vs. H II fractions) based on the local radiation field.
- Assumptions:
  - The ionization parameter $\xi = L/n_H r^2$ was calculated using the spectral energy distribution (SED) from the RHD simulations.
  - The disk component was treated as isotropic (dominated by the corona component for ionization).
  - Calculations focused on the inner region (0.0–0.1 pc) where ionization effects are most critical, though the model confirmed results were stable up to 2 pc.
- Definition of Covering Factors:
  - $C_{22}$ : Fraction of solid angle with $N_H > 10^{22}$ cm $^{-2}$ .
  - $C_{24}$ : Fraction of solid angle with $N_H > 10^{24}$ cm $^{-2}$ .
  - Calculations were performed for Total Gas (H I + H II) and Neutral Gas (H I only).

3. Key Results

A. Ionization Structure and H I Distribution

Ionization Parameter ( $\xi$ ): In the inner region ( $< 0.1$ pc), gas with polar angles $\theta < 70^\circ$ is highly ionized ( $\log \xi > 2$ ), meaning hydrogen is almost entirely H II.
Neutral Gas Survival: Neutral hydrogen (H I) survives only in regions where $\theta > 70^\circ$ (near the equatorial plane) and where $\log \xi < -2$ .
Independence from Eddington Ratio: The spatial distribution of H I is remarkably stable across all tested Eddington ratios. Even at the lowest ratio ( $\log R_{Edd} = -3$ ), the gas within 0.1 pc is photoionized.

B. Covering Factors ( $C_{22}$ and $C_{24}$ )

Total Gas (H I + H II): The covering factor is $\approx 70\%$ and is independent of the Eddington ratio.
Neutral Gas (H I Only):
- Compton-thin ( $C_{22}$ ): Approximately 30%, independent of the Eddington ratio.
- Compton-thick ( $C_{24}$ ): Approximately 30%, independent of the Eddington ratio.
- Note: In the simulated region ( $< 0.1$ pc), all Compton-thin gas is photoionized, so $C_{24}$ (H I) equals $C_{22}$ (H I).

C. Comparison with Observations

Compton-thick ( $C_{24}$ ): The theoretical value of ~30% matches X-ray observational data (Ricci et al. 2017b) perfectly across all Eddington ratios.
Compton-thin ( $C_{22}$ ):
- High Eddington Ratio ( $\lambda_{Edd} > 10^{-1}$ ): Theoretical $C_{22}$ (~30%) matches observations.
- Low Eddington Ratio ( $\lambda_{Edd} < 10^{-2}$ ): Observations show high $C_{22}$ (80–90%), but the theory predicts only ~30%.

4. Discussion and Proposed Mechanism

The authors propose that the discrepancy at low Eddington ratios is due to the spatial extent of the absorber, not the physics within the inner 0.1 pc.

Photoionization vs. Radiation Pressure: The study concludes that within 0.1 pc, the covering factor is determined by photoionization (X-rays stripping electrons) rather than radiation pressure blowing dust away. The hydrogen column density in this region is dominated by dust-free, hot gas, which remains gravitationally bound and dense regardless of the Eddington ratio.
The Missing Component: To explain the high $C_{22}$ $C_{22}$ observed at low Eddington ratios, the authors argue that Compton-thin neutral gas must exist beyond the current computational domain.
- Using ionization balance equations, they estimate that for Compton-thin H I to survive at $\lambda_{Edd} = 10^{-2}$ , it must exist at distances $r \gtrsim 30$ pc (assuming a filling factor of $10^{-2}$ ).
- At high Eddington ratios, even this distant gas is likely photoionized, reducing the covering factor to the intrinsic ~30% seen in the inner torus.

5. Significance and Conclusion

This paper refines the AGN Unification Model by integrating detailed photoionization physics into hydrodynamic simulations.

Resolution of the "Low Eddington Ratio" Problem: It demonstrates that the high covering factors observed in low-luminosity AGNs cannot be explained by the inner torus alone. Instead, they require a large-scale distribution of neutral gas (tens of parsecs away) that is shielded from the central X-ray source.
Validation of the Radiation-Driven Fountain: The consistency of the Compton-thick covering factor ( $C_{24} \approx 30\%$ ) with observations supports the radiation-driven fountain model as the mechanism forming the inner torus.
Revised Unification Scheme: The authors propose a new schematic (Figure 8) where:
- Inner Region ( $10^{-3} - 10^{-1}$ pc): Contains Compton-thick neutral gas ( $C_{24} \approx 30\%$ ) formed by the fountain model.
- Outer Region ( $\sim 30$ pc): Contains Compton-thin neutral gas only when the Eddington ratio is low ( $< 10^{-2}$ ). At high ratios, this outer gas is ionized, leading to a lower observed $C_{22}$ .

Final Takeaway: The covering factor of AGNs is not solely a function of radiation pressure blowing away dust; it is critically governed by the photoionization equilibrium of the gas. The observed dependence of the covering factor on the Eddington ratio is likely a result of the ionization state of gas at large radii ( $>10$ pc), rather than the structure of the inner torus itself.