Dusty winds in active galactic nuclei: reconciling observations with models

TL;DR

This paper introduces a new radiative transfer model for AGN IR emission that accounts for polar IR emission observed in many AGN, reconciling previous models with recent spatially resolved observations.

Contribution

The paper presents the CAT3D-WIND model, incorporating a disk and wind structure to better match IR observations of AGN, especially the polar emission features.

Findings

Disk+wind models match observed IR spectral features.

Reproduce the 3-5μm bump in type 1 AGN.

Explain IR interferometry polar emission features.

Abstract

This letter presents a revised radiative transfer model for the infrared (IR) emission of active galactic nuclei (AGN). While current models assume that the IR is emitted from a dusty torus in the equatorial plane of the AGN, spatially resolved observations indicate that the majority of the IR emission from 100 pc in many AGN originates from the polar region, contradicting classical torus models. The new model CAT3D-WIND builds upon the suggestion that the dusty gas around the AGN consists of an inflowing disk and an outflowing wind. Here, it is demonstrated that (1) such disk+wind models cover overall a similar parameter range of observed spectral features in the IR as classical clumpy torus models, e.g. the silicate feature strengths and mid-IR spectral slopes, (2) they reproduce the 3-5{\mu}m bump observed in many type 1 AGN unlike torus models, and (3) they are able to explain polar emission features seen in IR interferometry, even for type 1 AGN at relatively low inclination, as demonstrated for NGC3783. These characteristics make it possible to reconcile radiative transfer models with observations and provide further evidence of a two-component parsec-scaled dusty medium around AGN: the disk gives rise to the 3-5{\mu}m near-IR component, while the wind produces the mid-IR emission. The model SEDs will be made available for download.

Tables

Table 1. Table 1: Range of model parameters used to create the parameter grid.

disk parameters
$a$	$[-0.5,-1.0,-1.5,-2.0,-2.5,-3.0]$
$N_0$	$[5,7.5,10]$
$h$	$[0.1,0.2,0.3,0.4,0.5]$
wind parameters
$a_w$	$[-0.5,-1.0,-1.5,-2.0,-2.5]$
${\theta }_w$	$[45°,30°]$
${\sigma }_{\theta }$	$[7.5,10,15]$
global parameters
$f_{wd}$	$[0.15,0.3,0.45,0.6,0.75]$
	additionally $[1,1.25,1.5,1.75]$ for ${\theta }_w=30°$
$R_{\mathrm{out}}$ (1,2)	$500$ (large graphite) / $151$ (ISM)
$R_{\mathrm{cl}}$ (1)	$0.035\cdot r$
${\tau }_V$	$50$
line of sight
inclination	$[0°,15°,30°,45°,60°,75°,90°]$

Equations

$$\alpha =-\frac{\log \left(F_{\nu }({\lambda }_j)/F_{\nu }({\lambda }_i)\right)}{\log \left({\lambda }_j/{\lambda }_i\right)}$$

Full text

galaxies: active — infrared: galaxies — radiative transfer — galaxies: individual (NGC3783)

††journal: ApJL††software: CAT3D (Hönig & Kishimoto, 2010)

1 Introduction

It was held that the infrared (IR) emission from active galactic nuclei (AGN) originates from a ”dusty torus” around the central black hole and accretion disk (Antonucci, 1993; Netzer, 2016). This geometrically thick structure would explain the relative strength of the IR emission as well as the angle-dependent obscuration of the accretion disk and broad line region, leading to the well-established type 1/type 2 dichotomy of unobscured and obscured objects. Radiative transfer models based on such a toroidal geometry have proven successful in reproducing the IR spectral energy distribution (SED) of AGN (e.g. Fritz et al., 2006; Ramos Almeida et al., 2009; Alonso Herrero et al., 2011). However, given the milliarcsecond sizes of the IR emission in most AGN, the true geometric distribution of the emission was never scrutinised up until recently.

Advances in IR interferometry made it possible to resolve the near- and mid-IR emission in about three dozen AGN (for a recent overview, see Burtscher et al., 2016). A growing number of nearby Seyfert galaxies have a sufficient interferometric measurements to determine the wavelength- and angular-dependent size of the emission. It has been shown that contrary to the classical torus picture, the bulk of the mid-IR emission originates from the polar region of the AGN rather than the equatorial plane (e.g. Raban et al., 2009; Hönig et al., 2012, 2013; Tristram et al., 2014; López-Gonzaga et al., 2016). Detailed studies of the type 1 AGN NGC3783 and the type 2 sources NGC 1068 and the Circinus galaxy imply that the IR emission emerges from at least two distinct parsec-scale structures: (1) a compact, geometrically thin disk in the equatorial region of the AGN, and (2) an extended, elongated polar structure, which is cospatial with the outflow region of the AGN on larger scales (Asmus et al., 2016). The disk seems to be contributing a relatively larger fraction to shorter wavelengths, while the polar region dominates the overall IR emission energetically, which requires a significant covering fraction. Individual aspects of this behaviour may be explained with radiative transfer effects within a torus; however, the combination of extension, direction and emitted radiation strongly disfavours current torus models, either smooth or clumpy.

One of the hypotheses put forward to explain this structure leans on the idea that matter accreting onto the AGN is subject to radiative and hydrodynamic pressure producing strong outflows (e.g. Königl & Kartje, 1994; Elvis, 2000; Roth et al., 2012). In the specific picture outlined in Hönig et al. (2012), the disk component represents the inflowing dusty gas. Near the sublimation radius, the radiation pressure from both the AGN and the dust itself causes material to lift up and form a hollow outflow cone. Both outflow and disk contribute to the obscuration. This paper presents the radiative transfer model CAT3D-WIND of a clumpy disk and outflow with the geometry and dust composition built upon the hypothesis outlined in Hönig et al. (2013). It will be shown that such a configuration reproduces key observational characteristics of AGN, including the overall IR SEDs, the type 1/type 2 dichotomy, the observed disk plus polar emission structure, and the $3-5\,\micron$ bump seen in many unobscured AGN.

2 CAT3D-WIND: disk + wind model description

The CAT3D-WIND model is based on the parametrisation of a clumpy dust torus used in the underlying CAT3D model (for details see Hönig & Kishimoto, 2010). It distributes dust clumps according to a radial power law $\propto r^{a}$, where $a$ is the power-law index and $r$ denotes the distance from the AGN in units of the sublimation radius $r_{\mathrm{sub}}$. The sublimation radius forms the inner boundary of the dust distribution where the dust becomes too hot to survive. The vertical distribution is characterised by a Gaussian distribution $\propto\exp(-z^{2}/(2(hr)^{2}))$ with dimensionless scale height $h$ and vertical distance from the mid-plane $z$ in units of $r_{\mathrm{sub}}$. The model space is limited from $r_{\mathrm{sub}}$ as the inner boundary to the outer radius $R_{\mathrm{out}}$. Please note that $R_{\mathrm{out}}$ should not be considered a free model parameter, but rather a minimum boundary that has to be chosen to not cause an artificial brightness cut-off in the wavelength range of interest (for a detailed discussion of this issue, see Hönig & Kishimoto, 2010, Sect. 4.1.4). Obscuration is defined by the average number of dust clouds $N_{0}$ along an equatorial line-of-sight. The clouds are characterised by their optical depth $\tau_{\mathrm{cl}}$ and radius $R_{\mathrm{cl}}$, but the size of the clouds only controls the total number of clouds and not the typical model SED (we choose a power law of the form $R_{\mathrm{cl}}=R_{\mathrm{cl;0}}\cdot r$ for numerical reasons; for further details see Hönig & Kishimoto, 2010). Moreover, as long as the cloud is optically thick in the IR, which most clumpy models assume, the dominating emission from the hot side of the clouds is insensitive to the exact choice of $\tau_{\mathrm{cl}}$. In summary, these 6 parameters fully define typical clumpy torus models.

Classical torus models do not account for the observed predominance of polar emission in the mid-IR emission of AGN. Thus, a second component is added to the standard model in the form of a polar outflow. This structure is modelled as a hollow cone and characterised by three parameters: (1) the radial distribution of dust clouds in the wind $a_{w}$, (2) the half-opening angle of the wind $\theta_{w}$, and (3) its angular width $\sigma_{\theta}$. Finally, a wind-to-disk ratio $f_{wd}$ defines the ratio between number of clouds along the cone and $N_{0}$. Note that converting $f_{wd}$ to a mass ratio would require knowledge of $\tau_{\mathrm{cl}}$ and $R_{\mathrm{cl}}$. However, these are not constrained by the SEDs. For each set of model parameters, the emission as seen from inclinations of $0\arcdeg-90\arcdeg$ is simulated, in steps of $15\arcdeg$. An illustration of the dust distribution and typical view of a disk+wind is shown in Fig. 1.

The new model makes use of a more physical description of dust sublimation near the AGN than the commonly assumed single sublimation radius. It accounts for the facts that for any given density, the silicate sublimation temperature is lower than the graphite sublimation temperature (e.g. Phinney, 1989) and that small grains are hotter for a fixed distance from the AGN. When the dust heats above 1200 K, silicates are removed from the dust composition, leaving only graphites that can heat up to 1900 K. In the intermediate stage from 1200 K leading up to 1900 K, the smallest graphite grains are removed, so that at the innermost radius where dust can survive only grains with a size between 0.075$\,\micron$ and 1$\,\micron$ are present. In this framework, it is possible to reconcile the observed small near-IR reverberation mapping sizes and near-IR interferometry sizes with dust sublimation physics (e.g. Kishimoto et al., 2007).

One specific aspect of the disk+wind hypothesis put forward in Hönig et al. (2013) is that the wind is launched near the sublimation zone of the dusty disk. Thus, the chemical composition of the dust in the wind is expected to be very similar to the one seen at about the sublimation radius. Accordingly, the model dust clouds in the polar region are devoid of silicates and small grains. This will also help near edge-on views of the disk+wind structure to display silicate absorption features at $\sim$10 $\micron$ from the disk despite the direct view to hotter dust in the outflow region. Wind-only models with standard silicate and graphite mixtures predominantly show silicate emission features (e.g. Keating et al., 2012; Gallagher et al., 2015). This is confirmed with the presented models when using normal ISM dust for the wind clouds. Such models show exclusively silicate emission features for the range of observed SED mid-IR slopes of both type 1 and type 2 AGN (see Sect. 3.1).

3 Results and discussion

The parameter space of the model has been explored in the range of parameters as listed in Table 1. A total of 132,300 model SEDs have been simulated. It should be noted that the range of $f_{wd}$ is larger for $\theta_{w}=30\arcdeg$ than for $45\arcdeg$. In general, the covered surface for a cone with half-opening angle closer to the polar axis is smaller than for a cone with a larger half-opening angle. To counter this effect, the effectively filled area in the cone for the models with smaller half-opening angle has been increased.

3.1 SED shapes and silicate features

To obtain a swift overview of the typical SED shapes of all the models, each model SED has been broken down into three parameters: (1) the spectral index $\alpha_{\mathrm{MIR}}$ of the mid-IR emission, expressed by the flux ratio between $14\,\micron$ and $8\,\micron$, (2) the spectral index $\alpha_{\mathrm{IR}}$ from the near-IR to the mid-IR, expressed by the flux ratio between $6\,\micron$ and $3\,\micron$, and (3) the silicate feature optical depth $\tau_{\mathrm{Si}}$. Here, the spectral index is defined as

[TABLE]

where $F_{\nu}(\lambda)$ refers to the flux at wavelength $\lambda$ and $\lambda_{j}>\lambda_{i}$. Fig. 2 shows an overview of the range of SED shapes obtained by the CAT3D-WIND model based on the parameter grid in Table 1. In addition, a corresponding grid for the classical clumpy torus based on CAT3D is shown for comparison by setting $f_{wd}=0$. Greater spectral indices relate to bluer SEDs. For $\tau_{\mathrm{Si}}$, blue colours indicate silicate emission features and red colours silicate absorption features. The range of $\alpha_{\mathrm{MIR}}$ has been studied in the past for a variety of AGN classes (e.g. Hao et al., 2007; Hönig et al., 2010; Alonso Herrero et al., 2016; García-González et al., 2015). In general, Seyfert galaxies and quasars fall in the range of $-0.5>\alpha_{\mathrm{MIR}}>-3$, with type 1 AGN showing slightly bluer SEDs than type 2 AGN. This range of observed parameters is well covered by both model grids, illustrating the general viability of the disk+wind models as well as the classical torus models.

One interesting aspect of the covered model space is the relation of $\alpha_{\mathrm{MIR}}$ and $\alpha_{\mathrm{IR}}$. Fig. 2 shows a black dashed line that marks equal spectral index in the near- and mid-IR. Generally, classic clumpy torus models cover the parameter space above this line, where the near-IR rises relatively steeply towards the mid-IR (see Hönig et al., 2010, Figs. 5 & 6). Then, the SEDs flatten in the mid-IR and turn over longward of $\sim 20-30\,\micron$. Such models struggle to reproduce the $3-5\,\micron$ bump that is seen in several type 1 AGN, including quasars (e.g. Edelson & Malkan, 1986; Kishimoto et al., 2011; Mor & Netzer, 2012). Those SEDs are characterised by a first, local emission peak in the near-IR regime of the SED before they become steeper again towards the global $20-30\,\micron$ peak. SEDs with a $3-5\,\micron$ bump are found below the dashed line where the near-to-mid-IR spectral index $\alpha_{\mathrm{IR}}$ is bluer than the mid-IR spectral index $\alpha_{\mathrm{IR}}$. Indeed, the disk+wind models of CAT3D-WIND naturally produce such SEDs for a significant number of parameter combinations. Most notably, the $3-5\,\micron$ bump emerge from models where the disk has a relatively steep dust cloud distribution of $-2>a>-3$ while the distribution of clouds in the wind is much shallower with $-0.5>a>-1.5$. In these cases, the near-IR emission is dominated by the hot dust in the disk and the wind primarily contributes at longer wavelengths, which is the scenario implied by interferometry observations of NGC3783 (Hönig et al., 2013).

The silicate emission features in most of the CAT3D-WIND model SEDs are relatively shallow while the absorption features can be deeper. Qualitatively, this is consistent with observational constraints for both space-based and ground-based AGN samples (e.g. Hao et al., 2007; Hönig et al., 2010; Alonso Herrero et al., 2016). The distribution of observed feature strengths peak in the region of about $0.1\lesssim\tau_{\mathrm{Si}}\lesssim-0.4$ for quasars, $0.3\lesssim\tau_{\mathrm{Si}}\lesssim-0.3$ for Seyfert 1 galaxies, and $1\lesssim\tau_{\mathrm{Si}}\lesssim-0.1$ for Seyfert 2 AGN. The model space of the presented parameter grid covers these values.

In summary, the model SEDs of the disk+wind model cover the observed SED shapes and silicate feature depths quite well, very similar to classical clumpy torus models. In addition and beyond the capabilities of the torus models, the disk+wind model is able to reproduce the observed $3-5\,\micron$ bump in type 1 AGN from a compact disk, while preserving the mid-IR bump produced by the wind. This makes it superior to classical clumpy torus models on the criterium of reproducing SEDs alone.

3.2 Polar elongations in type 1 AGN: NGC3783 as a case study

Type 1 AGN NGC3783 displays a pronounced $3-5\,\micron$ bump (see Fig. 3). Decomposing the SED of NGC3783 similarly to the CAT3D-WIND model SEDs results in $\alpha_{\mathrm{MIR}}=-2.05\pm 0.20$, $\alpha_{\mathrm{IR}}=-1.0\pm 0.2$, and $\tau_{\mathrm{Si}}=-0.07\pm 0.06$. This puts NGC3783 in a range that is only covered by the disk+wind models, not the classical clumpy torus models (see grey contours in Fig. 2 left & right). In Fig. 3 we plot the range of model SEDs that fall within the uncertainty interval of these parameters as a grey band. The majority of these models are characterised by a compact, relatively flat disk ($a\sim-2.0\ldots-3$; $h\lesssim 0.3$) and an extended wind ($a_{w}\gtrsim-1$). Inclinations range from $0\arcdeg\leq i\leq 60\arcdeg$ indicating a moderately inclined AGN, as can be expected for a type 1. The SED-only models are inconclusive with regards to the half-opening angle of the outflow.

The major motivation for the disk+wind model stems from IR interferometry of nearby AGN. Therefore, viability of the new models rests on the ability to not only reproduce the SEDs but also the interferometry measurements. In general, it can be expected that a type 2 AGN seen edge on will indeed show elongation towards the polar region in these models. On the other hand, radiative transfer effects might also produce some polar emission in obscured AGN in classical torus models (e.g. Schartmann et al., 2005, 2008). The more challenging objects are moderately-inclined type 1 AGN where classical models would produce an elongated, disk-like structure with orientation along the disk plane.

NGC3783 can be considered a prototypical inclined type 1 AGN with dominating polar IR emission. It has been extensively covered with near- and mid-IR interferometry. Hönig et al. (2013) report simple modelling of 41 VLTI MIDI observations in the mid-IR and 6 VLTI AMBER measurements in the near-IR (AMBER data originally from Weigelt et al., 2012). The IR SED shows the characteristic $3-5\,\micron$ bump with an estimated hot-dust covering factor of about 30%. The near-to-mid-IR photometric and interferometric data is best reproduced by two IR emitting components, characterised by two perpendicular elongated structures, with the polar-oriented dominating in the mid-IR. The mid-IR emission is very elongated with a major-to-minor axis ratio of $\sim$3:1, which is challenging to achieve in classical torus models for an unobscured type 1 AGN at moderate inclinations.

The best SED model fits (see grey contour in Fig. 2, right) were selected to simulate multi-wavelength images, which is computationally expensive. Interferometric visibilities were obtained for the model images at $12\,\micron$ and compared to the corresponding interferometric visibilities of NGC3783 (see Hönig et al., 2013, Fig. 7). It is immediately found that the strong elongation seen in NGC3783 disfavour low inclinations at $0\arcdeg-15\arcdeg$, as expected. The strongest elongations are seen either at higher inclination, where the cone is seen from the side, or when the line-of-sight into the cone is close to its edge. Stronger polar elongations are found for $\theta_{w}=30\arcdeg$ than for $\theta_{w}=45\arcdeg$.

Fig. 4 shows model visibilities for several position angles for a representative model with $a=-3$, $a_{w}=-0.5$, $N_{0}=10$, $\theta_{w}=38\arcdeg$, $\sigma_{\theta}=10\arcdeg$ and $f_{wd}=1.2$ seen under an inclination of $i=30\arcdeg$. This orientation and cone shape implied by this model are in between the suggestions of Müller-Sánchez et al. (2011) and Fischer et al. (2013), who infer $i=60\arcdeg/\theta_{w}=30\arcdeg$ and $i=15\arcdeg/\theta_{w}=45\arcdeg$, respectively. While Fischer et al. (2013) argue that kinematics disfavour the high-inclination solution, the strong elongation seen in IR interferometry might indicate a higher inclination. On the other hand, it is shown here that a configuration where the line-of-sight is close to the cone edge produces strong elongations as well.

3.3 Anisotropy of the mid-IR emission

An optically thick dusty torus causes significant anisotropy of the emerging IR emission depending on inclination: More IR radiation escapes towards low inclinations than towards highly inclined line-of-sights. However, when taking the intrinsic X-ray emission as an isotropic tracer for the AGN luminosity, it has been shown that the mid-IR emission is surprisingly isotropic (e.g. Lutz et al., 2004; Gandhi et al., 2015; Asmus et al., 2015), unlike expected from smooth torus models. While this small scatter in the mid-IR-to-X-ray luminosity correlation can be partly explained in the framework of clumpy torus models, the range of torus parameters leading to low anisotropy is very restricted, which would imply a physical mechanism that restricts those parameters (e.g. Hönig et al., 2011).

In the disk+wind picture, the hollow cone is visible to the observer at any inclination, naturally leading to low anisotropy of the mid-IR emission. Quantitatively, the peak of the luminosity-normalised flux distributions at $12\,\micron$ shifts by about 0.4 dex from $0\arcdeg$ to $90\arcdeg$. When considering typical inclinations of $30\arcdeg$ for type 1 AGN and $75\arcdeg$ for type 2s, the anisotropy is $0.24^{+0.11}_{-0.16}$ dex. This is consistent with estimates from the relation between X-ray and mid-IR luminosities in AGN, finding $<0.3$ dex (Asmus et al., 2015) and samples isotropically selected in the radio at $\sim 0.22$ dex (Hönig et al., 2011).

4 Summary and Conclusions

In this letter, a new radiative transfer model CAT3D-WIND is presented to reconcile spatially resolved observations of AGN in the IR with model predictions. The model consists of a geometrically thin disk of optically-thick dust clumps and a hollow cone of putatively outflowing dust clouds. It is phenomenologically motivated by recent results in IR interferometry that hint towards such a two-component structure. In this picture, dust clouds are accreted in the plane of the disk. As dust partially sublimates near the sublimation radius (leading to a dearth of silicates), some dust clouds are lifted up by radiation pressure and flow out into the polar region. Based on a parametrisation of this picture, the following has been found:

• •

The new disk+wind models cover a similar parameter space in terms of mid-IR spectral slope and strength of the silicate feature as classical clumpy torus models.

• •

It is found that the new models are able to explain the $3-5\,\micron$ SED bump seen in many type 1 AGN. This local emission peak appears from the disk when it has a more compact distribution of dust clouds than the wind.

• •

Many of the disk+wind models do indeed show polar elongation for inclinations typically associated with both type 1 and type 2 AGN. Indeed, a model has been found to simultaneously reproduce the IR SED and IR interferometry of the type 1 AGN NGC3783.

• •

Exposure of the cone to the observer at all inclinations naturally results in a low degree of anisotropy of the mid-IR emission with respect to the inclination or viewing angle. This is consistent with the observed low scatter of the mid-IR/X-ray luminosity correlation.

The presented models are radiative transfer models of an empirical hypothesis. However, the physical mechanisms to drive the wind are not clear yet. It is possible that a combination of radiation pressure on the dust and neutral gas, both from the AGN radiation and from the IR-emitting medium itself, will play a role. Furthermore, hydrodynamic pressure close to the sublimation radius may help in launching such a wind. Work on radiative-hydrodynamical simulations will be needed to find a physical foundation for the observed polar structure (e.g. Dorotnisyn et al., 2011; Wada, 2012, 2015; Dorotnisyn et al., 2016; Chan & Krolik, 2016). The SEDs presented in this paper will be made available for download at http://cat3d.sungrazer.org.

We thank the anonymous referee for helpful comments that improved the manuscript. SFH acknowledges support for this work from the European Research Council Horizon 2020 grant DUST-IN-THE-WIND (677117). Part of the numerical computations were done on the Sciama High Performance Compute (HPC) cluster which is supported by the ICG, SEPNet and the University of Portsmouth. Special thanks go to D. Asmus, P. Gandhi, C. Ricci, M. Stalevski, and K. Tristram for helpful discussions.