We use mid-infrared spectroscopy of unobscured active galactic nuclei (AGNs) to reveal their native dusty environments. We concentrate on Seyfert 1 galaxies, observing a sample of 31 with the Infrared Spectrograph aboard the Spitzer Space Telescope, and compare them with 21 higher-luminosity quasar counterparts. Silicate dust reprocessing dominates the mid-infrared spectra, and we generally measure the 10 and 18 micron spectral features weakly in emission in these galaxies. The strengths of the two silicate features together are sensitive to the dust distribution. We present numerical radiative transfer calculations that distinguish between clumpy and smooth geometries, which are applicable to any central heating source, including stars as well as AGNs. In the observations, we detect the obscuring ``torus'' of unified AGN schemes, modeling it as compact and clumpy. We also determine that star formation increases with AGN luminosity, although the proportion of the galaxies' bolometric luminosity attributable to stars decreases with AGN luminosity.

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 10\mic\ 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.

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 10\mic silicate feature, in particular the lack of deep absorption features in AGN of any type. The weak 10\mic 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 10\mic emission feature in QSO2.

We analyze the mid-infrared (MIR) spectra of ultraluminous infrared galaxies (ULIRGs) observed with the Spitzer Space Telescope's Infrared Spectrograph. Dust emission dominates the MIR spectra of ULIRGs, and the reprocessed radiation that emerges is independent of the underlying heating spectrum. Instead, the resulting emission depends sensitively on the geometric distribution of the dust, which we diagnose with comparisons of numerical simulations of radiative transfer. Quantifying the silicate emission and absorption features that appear near 10 and 18um requires a reliable determination of the continuum, and we demonstrate that including a measurement of the continuum at intermediate wavelength (between the features) produces accurate results at all optical depths. With high-quality spectra, we successfully use the silicate features to constrain the dust chemistry. The observations of the ULIRGs and local sightlines require dust that has a relatively high 18/10um absorption ratio of the silicate features (around 0.5). Specifically, the cold dust of Ossenkopf et al. (1992) is consistent with the observations, while other dust models are not. We use the silicate feature strengths to identify two families of ULIRGs, in which the dust distributions are fundamentally different. Optical spectral classifications are related to these families. In ULIRGs that harbor an active galactic nucleus, the spectrally broad lines are detected only when the nuclear surroundings are clumpy. In contrast, the sources of lower ionization optical spectra are deeply embedded in smooth distributions of optically thick dust.

The silicate cross section peak near 10um produces emission and absorption features in the spectra of dusty galactic nuclei observed with the Spitzer Space Telescope. Especially in ultraluminous infrared galaxies, the observed absorption feature can be extremely deep, as IRAS 08572+3915 illustrates. A foreground screen of obscuration cannot reproduce this observed feature, even at large optical depth. Instead, the deep absorption requires a nuclear source to be deeply embedded in a smooth distribution of material that is both geometrically and optically thick. In contrast, a clumpy medium can produce only shallow absorption or emission, which are characteristic of optically-identified active galactic nuclei. In general, the geometry of the dusty region and the total optical depth, rather than the grain composition or heating spectrum, determine the silicate feature's observable properties. The apparent optical depth calculated from the ratio of line to continuum emission generally fails to accurately measure the true optical depth. The obscuring geometry, not the nature of the embedded source, also determines the far-IR spectral shape.