No Arabic abstract
We derive the properties of dusty tori in Active Galactic Nuclei (AGN) from the comparison of observed Spectral Energy Distributions (SEDs) of SDSS quasars and a precomputed grid of torus models. The observed SEDs comprise SDSS photometry, 2MASS J, H, and K data, whenever available and mid-Infrared (MIR) data from the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. The adopted model is that of Fritz et al., 2006. The fit is performed by standard chi^2 minimisation, the model however can be multi-component comprising a stellar and a starburst components, whenever necessary. Models with low equatorial optical depth, tau_9.7, were allowed as well as ``traditional models with tau_9.7 > 1.0, corresponding to A_V > 22 and the results were compared. Fits using high optical depth tori models only produced dust more compactly distributed than in the configuration where all tau_9.7 models were permitted. Tori with decreasing dust density with the distance from the centre were favoured while there was no clear preference for models with or without angular variation of the dust density. The computed outer radii of the tori are of some tens of parsecs large but can reach, in a few cases, a few hundreds of parsecs. The mass of dust, M_Dust, and infrared luminosity, L_IR, integrated in the wavelength range between 1 and 1000 micron, do not show significant variations with redshift, once the observational biases are taken into account. Objects with 70 micron detections, representing 25% of the sample, are studied separately and the starburst contribution (whenever present) to the IR luminosity can reach, in the most extreme but very few cases, 80%.
According to unified schemes of Active Galactic Nuclei (AGN), the central engine is surrounded by dusty, optically thick clouds in a toroidal structure. We have recently developed a formalism that for the first time takes proper account of the clumpy nature of the AGN torus. We now provide a detailed report of our findings in a two-paper series. Here we present our general formalism for radiative transfer in clumpy media and construct its building blocks for the AGN problem -- the source functions of individual dusty clouds heated by the AGN radiation field. We show that a fundamental difference from smooth density distributions is that in a clumpy medium, a large range of dust temperatures coexist at the same distance from the radiation central source. This distinct property explains the low dust temperatures found close to the nucleus of NGC1068 in 10 mic interferometric observations. We find that irrespective of the overall geometry, a clumpy dust distribution shows only moderate variation in its spectral energy distribution, and the 10mic absorption feature is never deep. Furthermore, the X-ray attenuating column density is widely scattered around the column density that characterizes the IR emission. All of these properties are characteristic of AGN observations. The assembly of clouds into AGN tori and comparison with observations is presented in the companion paper.
We investigated the gravitational microlensing of active galactic nucleus dusty tori in the case of lensed quasars in the infrared domain. The dusty torus is modeled as a clumpy two-phase medium. To obtain spectral energy distributions and images of tori at different wavelengths, we used the 3D Monte Carlo radiative transfer code SKIRT. A ray-shooting technique has been used to calculate microlensing magnification maps. We simulated microlensing by the stars in the lens galaxy for different configurations of the lensed system and different values of the torus parameters, in order to estimate (a) amplitudes and timescales of high magnification events, and (b) the influence of geometrical and physical properties of dusty tori on light curves in the infrared domain. We found that, despite their large size, dusty tori could be significantly affected by microlensing in some cases, especially in the near-infrared domain (rest-frame). The very long time-scales of such events, in the range from several decades to hundreds of years, are limiting the practical use of this method to study the properties of dusty tori. However, our results indicate that, when studying flux ratios between the images in different wavebands of lensed quasars, one should not disregard the possibility that the near and mid-infrared flux ratios could be under the influence of microlensing.
Warm gas and dust surround the innermost regions of active galactic nuclei (AGN). They provide the material for accretion onto the super-massive black hole and they are held responsible for the orientation-dependent obscuration of the central engine. The AGN-heated dust distributions turn out to be very compact with sizes on scales of about a parsec in the mid-infrared. Only infrared interferometry currently provides the necessary angular resolution to directly study the physical properties of this dust. Size estimates for the dust distributions derived from interferometric observations can be used to construct a size--luminosity relation for the dust distributions. The large scatter about this relation suggests significant differences between the dust tori in the individual galaxies, even for nuclei of the same class of objects and with similar luminosities. This questions the simple picture of the same dusty doughnut in all AGN. The Circinus galaxy is the closest Seyfert 2 galaxy. Because its mid-infrared emission is well resolved interferometrically, it is a prime target for detailed studies of its nuclear dust distribution. An extensive new interferometric data set was obtained for this galaxy. It shows that the dust emission comes from a very dense, disk-like structure which is surrounded by a geometrically thick, similarly warm dust distribution as well as significant amounts of warm dust within the ionisation cone.
From extensive radiative transfer calculations we find that clumpy torus models with No about 5--15 dusty clouds along radial equatorial rays successfully explain AGN infrared observations. The dust has standard Galactic composition, with individual cloud optical depth tV about 30--100 at visual. The models naturally explain the observed behavior of the 10mic silicate feature, in particular the lack of deep absorption features in AGN of any type. The weak 10mic emission feature tentatively detected in type 2 QSO can be reproduced if in these sources No drops to about 2 or tV exceeds about 100. The clouds angular distribution must have a soft-edge, e.g., Gaussian profile, the radial distribution should decrease as $1/r$ or $1/r^2$. Compact tori can explain all observations, in agreement with the recent interferometric evidence that the ratio of the torus outer to inner radius is perhaps as small as about 5--10. Clumpy torus models can produce nearly isotropic IR emission together with highly anisotropic obscuration, as required by observations. In contrast with strict variants of unification schemes where the viewing-angle uniquely determines the classification of an AGN into type 1 or 2, clumpiness implies that it is only a probabilistic effect; a source can display type 1 properties even from directions close to the equatorial plane. The fraction of obscured sources depends not only on the torus angular thickness but also on the cloud number No. The observed decrease of this fraction at increasing luminosity can be explained with a decrease of either torus angular thickness or cloud number, but only the latter option explains also the possible emergence of a 10mic emission feature in QSO2.
We explore the correlation between morphological types and mid-infrared (MIR) properties of an optically flux-limited sample of 154 galaxies from the Forth Data Release (DR4) of Sloan Digital Sky Survey (SDSS), cross-correlated with Spitzer SWIRE (Spitzer Wide-Area InfraRed Extragalactic Survey) fields of ELAIS-N1, ELAIS-N2 and Lockman Hole. Aperture photometry is performed on the SDSS and Spitzer images to obtain optical and MIR properties. The morphological classifications are given based on both visual inspection and bulge-disk decomposition on SDSS g- and r-band images. The average bulge-to-total ratio (B/T) is a smooth function over different morphological types. Both the 8um(dust) and 24um(dust) luminosities and their relative luminosity ratios to 3.6um (MIR dust-to-star ratios) present obvious correlations with both the Hubble T-type and B/T. The early-type galaxies notably differ from the late-types in the MIR properties, especially in the MIR dust-to-star ratios. It is suggested that the MIR dust-to-star ratio is an effective way to separate the early-type galaxies from the late-type ones. Based on the tight correlation between the stellar mass and the 3.6um luminosity, we have derived a formula to calculate the stellar mass from the latter. We have also investigated the MIR properties of both edge-on galaxies and barred galaxies in our sample. Since they present similar MIR properties to the other sample galaxies, they do not influence the MIR properties obtained for the entire sample.