Compton-thick AGN

• Compton-thick AGN are active galaxies with obscured supermassive black holes, featuring hydrogen column densities beyond 1.5×10^24 cm⁻².
• They are identified through X-ray spectral features such as a strong Fe Kα line, flattened continuum, and prominent reflection signatures.
• Multiwavelength analyses reveal complex torus geometries and host galaxy interactions, crucial for understanding AGN evolution and SMBH growth.

Compton-thick Active Galactic Nucleus (AGN) refers to an active galaxy in which the central engine—consisting of an accreting supermassive black hole—is obscured by a line-of-sight hydrogen column density exceeding $$N_{\rm H} \gtrsim 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$, so that the optical depth to Thomson scattering ($\tau_T \geq 1$) dominates over photoelectric absorption. In such systems, the direct X-ray continuum from the AGN is heavily attenuated, and most of the observed hard X-ray emission arises via reflection, scattering, and fluorescence from surrounding material. Compton-thick AGN are a fundamental, though observationally challenging, component of AGN population studies; their demographics, evolutionary history, and role in black hole-host galaxy co-evolution are central themes in extragalactic astrophysics (Gilli et al., 2011, Lanzuisi et al., 2015, Brightman et al., 2015).

1. Physical Definition and Observational Diagnostics

Compton-thick AGN are defined by a line-of-sight absorbing column density $$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$, where $\sigma_T$ is the Thomson cross-section. With such extreme obscuration:

• Direct AGN continuum at $E \lesssim 10$ keV is almost entirely extinguished.
• X-ray spectra are dominated by reflection and reprocessing features:
  • Pronounced Fe K$\alpha$ emission line at 6.4 keV, with equivalent widths often exceeding $1$–$2$ keV (Boorman et al., 2016).
  • Flattened continuum ($\Gamma \lesssim 1$), due to dominance of reflected emission below 10 keV (Gandhi et al., 2014).
• Hardness Ratio (HR): High HR values at high redshift signal heavy obscuration ($HR = 0.23 \pm 0.24$ for $\tau_T \geq 1$0 is strongly suggestive of Compton-thick absorption) (Gilli et al., 2011).
• Absorption Turnover: At energies above 10–20 keV, the spectrum shows a Compton hump and absorption rollover (Annuar et al., 2015, Brightman et al., 2015).

High-energy X-ray observatories (NuSTAR, Chandra, XMM-Newton, Swift-BAT) in conjunction with physically motivated torus models (e.g., MYTorus, borus02, Brightman & Nandra torus) are essential for constraining $\tau_T \geq 1$1 and distinguishing genuine Compton-thick AGN from less obscured, reflection-dominated sources.

2. Multiwavelength Identification and Classification Strategies

Given the suppression of direct X-ray emission, robust identification of Compton-thick AGN requires a multi-pronged approach:

• X-ray Spectral Fitting: Direct measurement of $\tau_T \geq 1$2 from broadband X-ray (0.3–80 keV) spectra using torus and clumpy models (Zhao et al., 2018, Zhao et al., 2018, Lanzuisi et al., 2014). For $\tau_T \geq 1$3, even hard X-rays are significantly diminished, leaving only reflected components.
• Luminosity Ratios: Low observed $\tau_T \geq 1$4/[NeV] and $\tau_T \geq 1$5 ratios, far below thresholds for unobscured AGN, confirm heavy obscuration (Lanzuisi et al., 2015).
• IR and Optical SED decomposition: A dominant torus component (e.g., $\tau_T \geq 1$6 of 5–40 μm emission), and spectral energy distributions requiring heavily obscured AGN to reproduce mid-IR data and [O IV]/12 μm ratios (Lanzuisi et al., 2015, Yamada et al., 2020).
• Optical emission-line diagnostics: Presence of high-ionization lines ([NeV], [FeVII]) with broad, blueshifted components indicate active AGN and outflow, even when broad lines are hidden (Lanzuisi et al., 2015).
• Morphological analysis (HST, Gemini): CT AGN hosts frequently show high merger/disturbance fraction, suggesting gas inflows enhance central obscuration (Lanzuisi et al., 2014).

Notably, many bona fide Compton-thick AGN are undetected or mis-classified in optical surveys due to host-dominated spectra and extinction (Annuar et al., 10 Jun 2025). Mid-IR diagnostics (e.g., 24 μm excess, deep Si 9.7 μm absorption) provide complementary, though not unique, selection (Georgantopoulos, 2012, Shi et al., 2013).

3. Covering Factor, Torus Geometry, and Evolutionary Scenarios

Across the population, Compton-thick AGN exhibit a broad range of torus covering factors ($\tau_T \geq 1$7), from $\tau_T \geq 1$8 to $\tau_T \geq 1$9 (Brightman et al., 2015, Zhao et al., 2018, Zhao et al., 2018, Yamada et al., 2020). There is a strong anti-correlation between covering factor and intrinsic X-ray luminosity:

$$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$0

(Brightman et al., 2015). High-luminosity AGN tend to have low covering factors ("disk-like" or clumpy tori), while lower-luminosity AGN are more heavily buried.

Clumpy/inhomogeneous torus models are favored by observations demonstrating significant differences between the instantaneous line-of-sight $$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$1 versus global/average values for individual AGN (Silver et al., 2022, Zhao et al., 2018). Rapid accretors ($$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$2) are more likely to be CT, while those with low accretion rates are generally unobscured (Annuar et al., 10 Jun 2025, Draper et al., 2010).

CT AGN populations are not a simple orientation-driven extension of type 2 AGN. Instead, population synthesis models requiring composite accretion scenarios—CT AGN arising both in rapid Eddington-limited growth phases (frequently merger fueled), and in weakly accreting, local SMBHs enshrouded by ambient molecular clouds—are most consistent with observed space densities, redshift distributions, and cosmic X-ray background (CXB) constraints (Draper et al., 2010, Shi et al., 2013).

4. Host Galaxy Properties and AGN-Galaxy Co-evolution

Compton-thick AGN are often embedded in massive, star-forming galaxies, with frequent incidence of barred structures, disturbed morphologies, or active merger signatures (Lanzuisi et al., 2015, Lanzuisi et al., 2014, Gandhi et al., 2014, Yamada et al., 2020). At high redshift ($$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$3), CT AGN co-exist with massive star formation rates ($$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$4), indicating coeval assembly of black hole and stellar mass (Gilli et al., 2011).

Key findings include:

• Black hole mass ($$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$5) distributions for CT AGN tend to be lower than those for unobscured AGN—by up to 1.5 dex at low luminosity ($$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$6).
• Eddington ratios are systematically higher ($$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$7–$$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$8) for CT AGN than for unobscured AGN at similar luminosity/redshift (Lanzuisi et al., 2015, Lanzuisi et al., 2014).
• Star formation rates are often on the main sequence, rather than in extreme starburst mode (Lanzuisi et al., 2014, Annuar et al., 10 Jun 2025).
• Host galaxy stellar masses are $$N_{\rm H} \gtrsim \sigma_T^{-1} \simeq 1.5 \times 10^{24}\ \mathrm{cm}^{-2}$$9 dex lower for low-luminosity CT AGN, while SFR distributions are similar to more powerful AGN hosts (Annuar et al., 10 Jun 2025).

Outflows traced by optical forbidden lines ([NeV], [FeVII]) and extended kpc-scale hard X-ray emission in ionization cones further substantiate energizing feedback from CT AGN into their host ISM (Lanzuisi et al., 2015, Silva et al., 2021, Jones et al., 2021). AGN feedback is plausibly implicated in quenching nuclear star formation and shaping galaxy evolution, supported by kinematic and stellar population analyses revealing aged, metal-rich nuclear regions in feedback-active hosts (Silva et al., 2021).

5. Demographics, Cosmic X-ray Background, and Evolutionary Impact

Compton-thick AGN are required by population synthesis models to reproduce the observed CXB spectrum, particularly the 20–30 keV peak (Georgantopoulos, 2012, Shi et al., 2013, Draper et al., 2010). Model predictions indicate:

• CT AGN comprise $\sigma_T$0–$\sigma_T$1 of all AGN at bright fluxes ($\sigma_T$2), but the fraction rises rapidly toward fainter flux levels (Shi et al., 2013).
• Space density peaks at $\sigma_T$3 few $\sigma_T$4 from $\sigma_T$5 to $\sigma_T$6, with luminous CT AGN density increasing strongly toward $\sigma_T$7–$\sigma_T$8 (Shi et al., 2013).
• CT AGN account for $\sigma_T$9 of cumulative SMBH mass accreted and contribute $E \lesssim 10$0 of the CXB at 20 keV (Shi et al., 2013).
• Recent IR-selected, volume-limited local surveys ($E \lesssim 10$1 Mpc) directly measure CT fractions $E \lesssim 10$2, higher than flux-limited hard X-ray samples, and access intrinsically faint AGN ($E \lesssim 10$3) (Annuar et al., 10 Jun 2025).

At low luminosities ($E \lesssim 10$4), the CT fraction ($E \lesssim 10$5) is similar to higher luminosity regimes, indicating a persistent population missed by previous X-ray/optical selections (Annuar et al., 10 Jun 2025).

6. Limitations, Outstanding Issues, and Future Directions

Detection of CT AGN is hampered by strong bias against even hard X-ray selection due to extreme absorption for $E \lesssim 10$6 (Georgantopoulos, 2012, Shi et al., 2013). IR selection, while effective, remains incomplete due to host contamination and overlap with starbursts (Gandhi et al., 2014). Deep, broadband ($E \lesssim 10$7 keV) X-ray observations and multiwavelength selection criteria are essential for robust CT AGN census (Lanzuisi et al., 2015, Jones et al., 2021).

Physically motivated torus modeling is critical; slab and reflection-only models (e.g., pexrav) are unsuited for CT AGN structure determination (Brightman et al., 2015, Boorman et al., 2016). Current torus models suffer from uncertainties related to geometry (clumpiness, covering factor, inclination), elemental abundances (particularly Fe), and spatial scales (Boorman et al., 2016, Zhao et al., 2018).

The inhomogeneous, clumpy nature of AGN tori—revealed by variance in $E \lesssim 10$8 and extended hard X-ray emission—points to complex radiative transfer effects and interplay between AGN and the host galaxy ISM (Zhao et al., 2018, Jones et al., 2021). Large-area, deeper surveys (e.g., COSMOS-Legacy, Chandra Deep Fields) and future high-resolution missions will refine the population census, the evolutionary history, and ultimately the role of CT AGN in SMBH and galaxy growth.

---

Compton-thick AGN thus represent a pivotal, yet challenging, demographic in the cosmic census of SMBH growth. Their detailed study illuminates the nexus of obscuration physics, feedback, multiwavelength selection, and galaxy evolution across cosmic time.