[Paper Review] The dusty heart of nearby active galaxies -- II. From clumpy torus models to physical properties of dust around AGN

S. Hoenig, Makoto Kishimoto|arXiv (Cornell University)|Sep 24, 2009

Galaxies: Formation, Evolution, PhenomenaPhysics and Astronomy1 references125 citations

TL;DR

This paper presents an upgraded 3D clumpy torus radiative transfer model (CAT3D) that improves the treatment of diffuse radiation in dusty tori around AGN, enabling better translation of SED and interferometric data into physical dust properties. It finds that standard ISM dust with $ au_V \sim 50$ best matches observed IR SEDs, and demonstrates that mid-IR spectral index and interferometric visibility features can constrain the radial dust distribution power-law index $a$, especially in type 1 AGN.

ABSTRACT

The dusty environments (= "dust tori'') of AGN are now in reach of observations. Following our paper I on ground-based mid-IR spectro-photometry (Hönig et al. 2010), we present an upgrade to our radiative transfer model of 3-dimensional clumpy dust tori. The upgrade with respect to Hönig et al. (2006) concerns an improved handling of the diffuse radiation field in the torus which is approximated by a statistical approach. The models are presented as tools to translate classical and interferometric observations into characteristic properties of the dust distribution. We compare model SEDs for different chemical and grain-size compositions of the dust and find that clouds with standard ISM dust and optical depth tau_V~50 appear in overall agreement with observed IR SEDs. By studying parameter dependencies, it is shown that type 1 AGN SEDs, in particular the mid-IR spectral index, can be used to constrain the radial dust cloud distribution power-law index 'a', while other parameters are more difficult to assess using SEDs only. Interferometry adds important additional information for modeling when interpreted simultaneously with the SED. Although type 2 AGN can, in principle, be used to constrain model parameters as well, obscuration effects make the analysis more ambiguous. We propose a simple, interferometry-based method to distinguish between "compact'' and "extended'' radial dust distributions without detailed modeling of the data and introduce a way to easily determine individual or sample average model parameters using the observed optical depth in the silicate feature and the mid-IR spectral index.

Motivation & Objective

• To improve radiative transfer modeling of clumpy dust tori by enhancing the treatment of the diffuse radiation field using a statistical approach.
• To enable direct translation of classical SED and interferometric observations into physical dust properties such as composition, optical depth, and radial distribution.
• To assess the constraints on dust parameters using SEDs alone versus combined SED and interferometry, particularly distinguishing type 1 and type 2 AGN.
• To develop a simple, interferometry-based method to classify dust distributions as 'compact' or 'extended' without full modeling.
• To provide a practical tool (CAT3D) for researchers to derive model parameters from observed silicate feature optical depth and mid-IR spectral index.

Proposed method

• Adopts a 3D clumpy torus geometry with randomly distributed dust clouds, using a statistical approximation for the diffuse radiation field to improve computational efficiency and realism.
• Employs a four-component power-law AGN continuum (from UV to X-ray) as the illuminating source, with spectral shape $F_\nu \propto \nu^{2}, \nu^{1/3}, \nu^{-1}, \/nu^{-2}$ in different frequency bands.
• Calculates dust absorption and re-emission using dust composition-dependent absorption efficiencies $Q_{\text{abs}}(\nu)$, with standard ISM dust (silicate + graphite) as primary test case.
• Simulates SEDs and mid-IR interferometric visibility curves (including spectral features) for various radial distributions characterized by power-law index $a$.
• Compares model outputs with observed IR SEDs and interferometric data from AGN such as NGC 1068 and NGC 4151 to assess parameter degeneracies and constraints.
• Introduces a method to estimate average model parameters using observed silicate feature optical depth and mid-IR spectral index, bypassing full modeling.

Experimental results

Research questions

• RQ1 How does an improved treatment of the diffuse radiation field in clumpy torus models affect the prediction of IR SEDs and visibility curves?
• RQ2 Which dust composition and optical depth best reproduce observed IR SEDs of nearby AGN?
• RQ3 To what extent can the mid-IR spectral index constrain the radial dust cloud distribution power-law index $a$ in type 1 AGN?
• RQ4 How does interferometry improve the constraint on dust torus parameters compared to SEDs alone?
• RQ5 Can a simple, interferometry-based method distinguish between compact and extended radial dust distributions without full radiative transfer modeling?

Key findings

• Standard ISM dust with a total visual optical depth $\tau_V \sim 50$ provides the best overall agreement with observed IR SEDs of nearby AGN.
• The mid-IR spectral index in type 1 AGN is a strong constraint on the radial dust cloud distribution power-law index $a$, with degeneracies reduced when combined with interferometric data.
• Silicate features in visibility curves show complex behavior: they appear in emission or absorption depending on $a$, and are nearly canceled at $a \sim -1$, but their strength is degenerate with multiple parameters.
• Interferometric data add critical information beyond SEDs, especially for constraining geometry and radial distribution, particularly in type 1 AGN.
• The accretion disk contributes only 5–50% to near-IR flux in type 1 AGN, with marginal impact on visibility amplitudes (≤0.1 variation), justifying its omission in the main model.
• A practical method is proposed to estimate average model parameters using observed silicate feature optical depth and mid-IR spectral index, enabling quick parameter estimation without full modeling.