Search

We describe a large simulation of the stars to be observed by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The simulation is based on the TRILEGAL code, which resorts to large databases of stellar evolutionary tracks, synthetic spectra, and pulsation models, added to simple prescriptions for the stellar density and star formation histories of the main structures of the Galaxy, to generate mock stellar samples through a population synthesis approach. The main bodies of the Magellanic Clouds are also included. A complete simulation is provided for single stars, down to the $r=27.5$ mag depth of the co-added wide-fast-deep survey images. A second simulation is provided for a fraction of the binaries, including the interacting ones, as derived with the BinaPSE module of TRILEGAL. We illustrate the main properties and numbers derived from these simulations, including: comparisons with real star counts; the expected numbers of Cepheids, long-period variables and eclipsing binaries; the crowding limits as a function of seeing and filter; the star-to-galaxy ratios, etc. Complete catalogs are accessible through the NOIRLab Astro Data Lab, while the stellar density maps are incorporated in the LSST metrics analysis framework (MAF)., Comment: Accepted for publication in the LSST focused ApJS issue

The upcoming photometric surveys, such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) will monitor unprecedented number of active galactic nuclei (AGN) in a decade long campaign. Motivated by the science goals of LSST, which includes the harnessing of broadband light curves of AGN for photometric reverberation mapping (PhotoRM), we implement the existing formalism to estimate the lagged response of the emission line flux to the continuum variability using only mutli-band photometric light curves. We test the PhotoRM method on a set of 19 artificial light curves simulated using a stochastic model based on the Damped Random Walk process. These light curves are sampled using different observing strategies, including the two proposed by the LSST, in order to compare the accuracy of time-lag retrieval based on different observing cadences. Additionally, we apply the same procedure for time-lag retrieval to the observed photometric light curves of NGC 4395, and compare our results to the existing literature.

Recent infrared interferometric observations revealed sub-parsec scale dust distributions around active galactic nuclei (AGNs). Using images of CLUMPY torus models and NGC 1068 as an example, we demonstrate that the near- and mid-infrared nuclear emission of some nearby AGNs will be resolvable in direct imaging with the next generation of 30~m telescopes, potentially breaking degeneracies from previous studies that used integrated spectral energy distributions of unresolved AGN tori. To that effect we model wavelength-dependent point spread functions from the pupil images of various telescopes: James Webb Space Telescope, Keck, Giant Magellan Telescope, Thirty Meter Telescope, and Extremely Large Telescope. We take into account detector pixel scales and noise, and apply deconvolution techniques for image recovery. We also model 2D maps of the 10-micron silicate feature strength, $S_{10}$, of NGC 1068 and compare with observations. When the torus is resolved, we find $S_{10}$ variations across the image. However, to reproduce the $S_{10}$ measurements of an unresolved torus a dusty screen of $A_V > 9$ mag is required. We also fit the first resolved image of the K-band emission in NGC 1068 recently published by the GRAVITY collaboration, deriving likely model parameters of the underlying dust distribution. We find that both 1) an elongated structure suggestive of a highly inclined emission ring, and 2) a geometrically thin but optically thick flared disk where the emission arises from a narrow strip of hot cloud surface layers on the far inner side of the torus funnel, can explain the observations., 21 pages, 11 figures, accepted for publication in ApJ

AbstarctWe introduce Hypercat, a large set of 2-d AGN torus images computed with the state-of-the-art clumpy radiative transfer code Clumpy. The images are provided as a 9-dimensional hypercube, in addition to a smaller hypercube of corresponding projected dust distribution maps. Hypercat also comprises a software suite for easy use of the hypercubes, quantification of image morphology, and simulation of synthetic observations with single-dish telescopes, interferometers, and Integral Field Units. We apply Hypercat to NGC 1068 and find that it can be spatially resolved in Near- and Mid-IR, for the first time with single-dish apertures, on the upcoming generation of 25–40m class telescopes. We also find that clumpy AGN torus models within a range of the parameter space can explain on scales of several parsec the recently reported polar elongation of MIR emission in several sources, while not upending basic assumptions about AGN unification.

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will detect an unprecedentedly large sample of actively accreting supermassive black holes with typical accretion disk (AD) sizes of a few light days. This brings us to face challenges in the reverberation mapping (RM) measurement of AD sizes in active galactic nuclei (AGNs) using interband continuum delays. We examine the effect of LSST cadence strategies on AD RM using our metric AGNTimeLagMetric. It accounts for redshift, cadence, the magnitude limit, and magnitude corrections for dust extinction. Running our metric on different LSST cadence strategies, we produce an atlas of the performance estimations for LSST photometric RM measurements. We provide an upper limit on the estimated number of quasars for which the AD time lag can be computed within 01000 sources in each Deep Drilling field (DDF, 10 sq. deg) in any filter, with the redshift distribution of these sources peaking at z~1. We find the LSST observation strategies with a good cadence (~ 5 days) and a long cumulative season (~9 yr), as proposed for LSST DDF, are favored for the AD size measurement. We create synthetic LSST light curves for the most suitable DDF cadences and determine RM time lags to demonstrate the impact of the best cadences based on the proposed metric., Comment: accepted for publication in the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

Near- and mid-infrared interferometers have resolved the dusty parsec-scale obscurer (torus) around nearby active galactic nuclei (AGNs). With the arrival of extremely large single-aperture telescopes, the emission morphology will soon be resolvable unambiguously, without modeling directly the underlying brightness distribution probed by interferometers today. Simulations must instead deliver the projected 2D brightness distribution as a result of radiative transfer through a 3D distribution of dusty matter around the AGN. We employ such physically motivated 3D dust distributions in tori around AGNs to compute 2D images of the emergent thermal emission using CLUMPY, a dust radiative transfer code for clumpy media. We demonstrate that CLUMPY models can exhibit morphologies with significant polar elongation in the mid-infrared (i.e. the emission extends perpendicular to the dust distribution) on scales of several parsecs, in line with observations in several nearby AGNs. We characterize the emission and cloud distribution morphologies. The observed emission from near- to mid-infrared wavelengths generally does not trace the bulk of the cloud distribution. The elongation of the emission is sensitive to the torus opening angle or scale height. For cloud distributions with a flat radial profile, polar extended emission is realized only at wavelengths shorter than ~18 micron, and shorter than ~5 micron for steep profiles. We make the full results available through HYPERCAT, a large hypercube of resolved AGN torus brightness maps computed with CLUMPY. HYPERCAT also comprises software to process and analyze such large data cubes and provides tools to simulate observations with various current and future telescopes., 28 pages, 18 figures, accepted for publication in ApJ

There is growing evidence for physical influence between supermassive black holes and their host galaxies. We present a case study of nearby galaxy NGC 7582, for which we find evidence that galactic substructure plays an important role in affecting the collimation of ionized outflows as well as contributing to the heavy active galactic nucleus (AGN) obscuration. This result contrasts with a simple, small-scale AGN torus model, according to which AGN wind collimation may take place inside the torus itself, at subparsec scale. Using 3D spectroscopy with the MUSE instrument, we probe the kinematics of the stellar and ionized gas components as well as the ionization state of the gas from a combination of emission line ratios. We report for the first time a kinematically distinct core (KDC) in NGC 7582, on a scale of ~600pc. This KDC coincides spatially with dust lanes and starbursting complexes previously observed. We interpret it as a circumnuclear ring of stars and dusty, gas-rich material. We obtain a clear view of the outflowing cones over kpc scales, and demonstrate that they are predominantly photoionized by the central engine. We detect the back cone (behind the galaxy), and confirm previous results of a large nuclear obscuration of both the stellar continuum and HII regions. While we tentatively associate the presence of the KDC to a large-scale bar and/or a minor galaxy merger, we stress the importance of gaining a better understanding of the role of galaxy substructure in controlling the fueling, feedback and obscuration of AGN., Comment: 31 pages, 16 figures (including two animations), accepted for publication in the Astrophysical Journal. For a preview of the animations, see https://noirlab.edu/science/videos/rotating-gas-disk-of-the-ngc-7582-galaxy for Figure 15, and see https://noirlab.edu/science/videos/black-hole-wind-from-the-ngc-7582-galaxy for Figure 16

We have conducted an extensive X-ray spectral variability study of a sample of 20 Compton-thin type II galaxies using broad band spectra from XMM-Newton, Chandra, and Suzaku. The aim is to study the variability of the neutral intrinsic X-ray obscuration along the line of sight and investigate the properties and location of the dominant component of the X-ray-obscuring gas. The observations are sensitive to absorption columns of $N_{\rm H} \sim 10^{20.5-24} {\rm cm^{-2}}$ of fully- and partially-covering neutral and/or lowly-ionized gas on timescales spanning days to well over a decade. We detected variability in the column density of the full-covering absorber in 7/20 sources, on timescales of months-years, indicating a component of compact-scale X-ray-obscuring gas lying along the line of sight of each of these objects. Our results imply that torus models incorporating clouds or overdense regions should account for line of sight column densities as low as $\sim$ a few $\times 10^{21}$ cm$^{-2}$. However, 13/20 sources yielded no detection of significant variability in the full-covering obscurer, with upper limits to ${\Delta}N_{\rm H}$ spanning $10^{21-23}$ cm$^{-2}$. The dominant absorbing media in these systems could be distant, such as kpc-scale dusty structures associated with the host galaxy, or a homogeneous medium along the line of sight. Thus, we find that overall, strong variability in full-covering obscurers is not highly prevalent in Compton-thin type IIs, at least for our sample, in contrast to previous results in the literature. Finally, 11/20 sources required a partial-covering, obscuring component in all or some of their observations, consistent with clumpy near-Compton-thick compact-scale gas., Comment: Accepted for publication in ApJ

The size and structure of the dusty circumnuclear torus in active galactic nuclei (AGN) can be investigated by analyzing the temporal response of the torus's infrared (IR) dust emission to variations in the AGN ultraviolet/optical luminosity. This method, reverberation mapping, is applicable over a wide redshift range, but the IR response is sensitive to several poorly constrained variables relating to the dust distribution and its illumination, complicating the interpretation of measured reverberation lags. We have used an enhanced version of our torus reverberation mapping code (TORMAC) to conduct a comprehensive exploration of the torus response functions at selected wavelengths, for the standard interstellar medium grain composition. The shapes of the response functions vary widely over the parameter range covered by our models, with the largest variations occurring at shorter wavelengths ($\leq 4.5\,\mu$m). The reverberation lag, quantified as the response-weighted delay (RWD), is most affected by the radial depth of the torus, the steepness of the radial cloud distribution, the degree of anisotropy of the AGN radiation field, and the volume filling factor. Nevertheless, we find that the RWD provides a reasonably robust estimate, to within a factor of $\sim 3$, of the luminosity-weighted torus radius, confirming the basic assumption underlying reverberation mapping. However, overall, the models predict radii at $2.2\,\mu$m that are a typically factor of $\sim 2$ larger than those derived from K-band reverberation mapping. This is likely an indication that the innermost region of the torus is populated by clouds dominated by large graphite grains., Comment: 38 pages, 36 figures (13 in appendices), published in the Astrophysical Journal

The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows, is not thoroughly studied. Here we report the first detection of polarized thermal emission by means of magnetically aligned dust grains in the dusty torus of NGC 1068 using ALMA Cycle 4 polarimetric dust continuum observations ($0.07"$, $4.2$ pc; 348.5 GHz, $860$ $\mu$m). The polarized torus has an asymmetric variation across the equatorial axis with a peak polarization of $3.7\pm0.5$\% and position angle of $109\pm2^{\circ}$ (B-vector) at $\sim8$ pc east from the core. We compute synthetic polarimetric observations of magnetically aligned dust grains assuming a toroidal magnetic field and homogeneous grain alignment. We conclude that the measured 860 $\mu$m continuum polarization arises from magnetically aligned dust grains in an optically thin region of the torus. The asymmetric polarization across the equatorial axis of the torus arises from 1) an inhomogeneous optical depth, and 2) a variation of the velocity dispersion, i.e. variation of the magnetic field turbulence at sub-pc scales, from the eastern to the western region of the torus. These observations and modeling constrain the torus properties beyond spectral energy distribution results. This study strongly supports that magnetic fields up to a few pc contribute to the accretion flow onto the active nuclei., Comment: 19 pages, 11 figures (Accepted for Publication to ApJ)

We announce the second data release (DR2) of the NOIRLab Source Catalog (NSC), using 412,116 public images from CTIO-4m+DECam, the KPNO-4m+Mosaic3 and the Bok-2.3m+90Prime. NSC DR2 contains over 3.9 billion unique objects, 68 billion individual source measurements, covers $\approx$35,000 square degrees of the sky, has depths of $\approx$23rd magnitude in most broadband filters with $\approx$1-2% photometric precision, and astrometric accuracy of $\approx$7 mas. Approximately 1.9 billion objects within $\approx$30,000 square degrees of sky have photometry in three or more bands. There are several improvements over NSC DR1. DR2 includes 156,662 (61%) more exposures extending over 2 more years than in DR1. The southern photometric zeropoints in $griz$ are more accurate by using the Skymapper DR1 and ATLAS-Ref2 catalogs, and improved extinction corrections were used for high-extinction regions. In addition, the astrometric accuracy is improved by taking advantage of Gaia DR2 proper motions when calibrating the WCS of individual images. This improves the NSC proper motions to $\sim$2.5 mas/yr (precision) and $\sim$0.2 mas/yr (accuracy). The combination of sources into unique objects is performed using a DBSCAN algorithm and mean parameters per object (such as mean magnitudes, proper motion, etc.) are calculated more robustly with outlier rejection. Finally, eight multi-band photometric variability indices are calculated for each object and variable objects are flagged (23 million objects). NSC DR2 will be useful for exploring solar system objects, stellar streams, dwarf satellite galaxies, QSOs, variable stars, high-proper motion stars, and transients. Several examples of these science use cases are presented. The NSC DR2 catalog is publicly available via the NOIRLab's Astro Data Lab science platform., 20 pages, 16 figures, submitted to AJ

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $��$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $��$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($, 13 pages, 7 figures, Accepted for publication by ApJ

We present 1.4 pc resolution observations of 256 GHz nuclear radio continuum and HCN ($J=3 \to 2$) in the molecular torus of NGC 1068. The integrated radio continuum emission has a flat spectrum consistent with free-free emission and resolves into an X-shaped structure resembling an edge-brightened bicone. HCN is detected in absorption against the continuum, and the absorption spectrum shows a pronounced blue wing that suggests a high-velocity molecular outflow with speeds reaching 450 km/s. Analysis of the off-nucleus emission line kinematics and morphology reveals two nested, rotating disk components. The inner disk, inside $r\sim 1.2$ pc, has kinematics consistent with the nearly edge-on, geometrically thin water megamaser disk in Keplerian rotation around a central mass of $1.66\times 10^7\,\mbox{M}_\odot$. The outer disk, which extends to $\sim 7$~pc radius, counter-rotates relative to the inner disk. The rotation curve of the outer disk is consistent with rotation around the same central mass as the megamaser disk but in the opposite sense. The morphology of the molecular gas is asymmetric around the nuclear continuum source. We speculate that the outer disk formed from more recently introduced molecular gas falling out of the host galaxy or from a captured dwarf satellite galaxy. In NGC 1068, we find direct evidence that the molecular torus consists of counter-rotating and misaligned disks on parsec scales., 8 pages, 4 figures; accepted for publication in The Astrophysical Journal Letters

Most of the sky has been imaged with NOAO's telescopes from both hemispheres. While the large majority of these data were obtained for PI-led projects and almost all of the images are publicly available, only a small fraction have been released to the community via well-calibrated and easily accessible catalogs. We are remedying this by creating a catalog of sources from most of the public data taken on the CTIO-4m+DECam and the KPNO-4m+Mosaic3. This catalog, called the NOAO Source Catalog (NSC), contains over 2.9 billion unique objects, 34 billion individual source measurements, covers ~30,000 square degrees of the sky, has depths of ~23rd magnitude in most broadband filters with ~1-2% photometric precision, and astrometric accuracy of ~7 mas. In addition, ~2 billion objects and ~21,000 square degrees of sky have photometry in three or more bands. The NSC will be useful for exploring stellar streams, dwarf satellite galaxies, QSOs, high-proper motion stars, variable stars and other transients. The NSC catalog is publicly available via the NOAO Data Lab service., 18 pages, 15 figures, accepted for publication in AJ. Revisions include more analysis and figures on the achieved data quality. The 34 billion row measurement table is now available through the NOAO Data Lab, see https://datalab.noao.edu

We have detected in ALMA observations CO J = 6 - 5 emission from the nucleus of the Seyfert galaxy NGC 1068. The low-velocity (up to +/- 70 km/s relative to systemic) CO emission resolves into a 12x7 pc structure, roughly aligned with the nuclear radio source. Higher-velocity emission (up to +/- 400 km/s) is consistent with a bipolar outflow in a direction nearly perpendicular (roughly 80 degrees) to the nuclear disk. The position-velocity diagram shows that in addition to the outflow, the velocity field may also contain rotation about the disk axis. These observations provide compelling evidence in support of the disk-wind scenario for the AGN obscuring torus., FONDECYT (Grant ID: 3140436), Science and Technology Facilities Council, This is the author accepted manuscript. The final version is available from Institute of Physics Publishing via http://dx.doi.org/10.3847/2041-8205/829/1/L7

We utilize a Bayesian approach to fit the observed mid-IR-to-submm/mm spectral energy distributions (SEDs) of 22 WISE-selected and submm-detected, hyperluminous hot dust-obscured galaxies. By adopting the Torus+GB model, we decompose the observed IR SEDs of Hot DOGs into torus and cold dust components. The main results are: 1) Hot DOGs in our submm-detected sample are hyperluminous, with torus emission dominating the IR energy output. However, cold dust emission is non-negligible, averagely contributing ~24% of total IR luminosity. 2) Compared to QSO and starburst SED templates, the median SED of Hot DOGs shows the highest luminosity ratio between mid-IR and submm at rest-frame, while it is very similar to that of QSOs at 10-50um suggesting that the heating sources of Hot DOGs should be buried AGNs. 3) Hot DOGs have both high dust temperatures ~73K and IR luminosity of cold dust. The T-L relation of Hot DOGs suggests that the increase in IR luminosity for Hot DOGs is mostly due to the increase of the dust temperature, rather than dust mass. Hot DOGs have lower dust masses than those of submillimeter galaxies (SMGs) and QSOs within the similar redshift range. Both high IR luminosity of cold dust and relatively low dust mass in Hot DOGs can be expected by their relatively high dust temperatures. 4) Hot DOGs have high dust covering factors, which deviate the previously proposed trend of the dust covering factor decreasing with increasing bolometric luminosity. Finally, we can reproduce the observed properties in Hot DOGs by employing a physical model of galaxy evolution. The result suggests that Hot DOGs may lie at or close to peaks of both star formation and black hole growth histories, and represent a transit phase during the evolution of massive galaxies, transforming from the dusty starburst dominated phase to the optically bright QSO phase. (abridged), Comment: 16 pages, 12 figures. Accepted by ApJ

The \textit{Spitzer}/Infrared Spectrograph spectra of three spectroscopically anomalous galaxies (IRAS~F10398+1455, IRAS~F21013-0739 and SDSS~J0808+3948) are modeled in terms of a mixture of warm and cold silicate dust, and warm and cold carbon dust. Their unique infrared (IR) emission spectra are characterized by a steep $\simali$5--8$\mum$ emission continuum, strong emission bands from polycyclic aromatic hydrocarbon (PAH) molecules, and prominent silicate emission. The steep $\simali$5--8$\mum$ emission continuum and strong PAH emission features suggest the dominance of starbursts, while the silicate emission is indicative of significant heating from active galactic nuclei (AGNs). With warm and cold silicate dust of various compositions ("astronomical silicate," amorphous olivine, or amorphous pyroxene) combined with warm and cold carbon dust (amorphous carbon, or graphite), we are able to closely reproduce the observed IR emission of these %spectroscopically anomalous galaxies. We find that the dust temperature is the primary cause in regulating the steep $\sim$5--8$\mum$ continuum and silicate emission, insensitive to the exact silicate or carbon dust mineralogy and grain size $a$ as long as $a\simlt1\mum$. More specifically, the temperature of the $\simali$5--8$\mum$ continuum emitter (which is essentially carbon dust) of these galaxies is $\sim$250--400$\K$, much lower than that of typical quasars which is $\sim$640$\K$. Moreover, it appears that larger dust grains are preferred in quasars. The lower dust temperature and smaller grain sizes inferred for these three galaxies compared with that of quasars could be due to the fact that they may harbor a young/weak AGN which is not maturely developed yet., Comment: 31 pages, 14 figures, accepted for publication in ApJ

A possible source of $\gamma$-ray photons observed from the jets of blazars is inverse Compton scattering by relativistic electrons of infrared seed photons from a hot, dusty torus in the nucleus. We use observations from the Spitzer Space Telescope to search for signatures of such dust in the infrared spectra of four $\gamma$-ray bright blazars, the quasars 4C 21.35, CTA102, and PKS 1510$-$089, and the BL Lacertae object ON231. The spectral energy distribution (SED) of 4C 21.35 contains a prominent infrared excess indicative of dust emission. After subtracting a non-thermal component with a power-law spectrum, we fit a dust model to the residual SED. The model consists of a blackbody with temperature $\sim1200$ K, plus a much weaker optically thin component at $\sim660$ K. The total luminosity of the thermal dust emission is $7.9\pm0.2 \times 10^{45}$ erg s$^{-1}$. If the dust lies in an equatorial torus, the density of IR photons from the torus is sufficient to explain the $\gamma$-ray flux from 4C 21.35 as long as the scattering occurs within a few parsecs of the central engine. We also report a tentative detection of dust in the quasar CTA102, in which the luminosity of the infrared excess is $7 \pm 2 \times 10^{45}$ erg s$^{-1}$. However, in CTA102 the far-IR spectra are too noisy to detect the $10 \mu$m silicate feature. Upper limits to the luminosity from thermal emission from dust in PKS 1510-089, and ON231, are, $2.3\times10^{45}$, and $6.6\times10^{43}$ erg s$^{-1}$, respectively. These upper limits do not rule out the possibility of inverse Compton up-scattering of IR photons to $\gamma$-ray energies in these two sources. The estimated covering factor of the hot dust in 4C 21.35, 22%, is similar to that of non-blazar quasars; however, 4C 21.35 is deficient in cooler dust., Comment: 23 Pages, 5 Figures, 2 Tables, 1 Machine Readable Table. Accepted to ApJ

The 10-micron silicate feature observed with Spitzer in active galactic nuclei (AGN) reveals some puzzling behavior. It (1) has been detected in emission in type 2 sources, (2) shows broad, flat-topped emission peaks shifted toward long wavelengths in several type 1 sources, and (3) is not seen in deep absorption in any source observed so far. We solve all three puzzles with our clumpy dust radiative transfer formalism. (1) We present the spectral energy distribution (SED) of SST1721+6012, the first type 2 quasar observed to show a clear 10-mic silicate feature in emission. We constructed a large database of clumpy torus models and performed extensive fitting of the observed SED, constraining several of the torus parameters. We find that the source bolometric luminosity is ~3*10^12 L_sun. Our modeling suggests that, 11 pages, 9 figures, 4 tables, accepted for publication in ApJ; minor revision, added reference in Section 1

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., Comment: To appear in ApJ September 20, 2008

We present high angular resolution ($\sim$0.5$^\prime$$^\prime$) MIR spectra of the powerful radio galaxy, Cygnus A, obtained with the Subaru telescope. The overall shape of the spectra agree with previous high angular resolution MIR observations, as well as previous Spitzer spectra. Our spectra, both on and off nucleus, show a deep silicate absorption feature. The absorption feature can be modeled with a blackbody obscured by cold dust or a clumpy torus. The deep silicate feature is best fit by a simple model of a screened blackbody, suggesting foreground absorption plays a significant, if not dominant role, in shaping the spectrum of Cygnus A. This foreground absorption prevents a clear view of the central engine and surrounding torus, making it difficult to quantify the extent the torus attributes to the obscuration of the central engine, but does not eliminate the need for a torus in Cygnus A.

We present Spitzer infrared spectra and ultra-violet to mid-infrared spectral energy distributions (SEDs) of 25 luminous type 1 quasars at z \sim 2. In general, the spectra show a bump peaking around 3 {\mu}m, and the 10 {\mu}m silicate emission feature. The 3 {\mu}m emission is identified with hot dust emission at its sublimation temperature. We explore two approaches to modeling the SED: (i) using the Clumpy model SED from Nenkova et al. (2008a), and (ii) the Clumpy model SED, and an additional blackbody component to represent the 3 {\mu}m emission. In the first case, a parameter search of \sim 1.25 million Clumpy models shows: (i) if we ignore the UV-to-near-IR SED, models fit the 2-8 {\mu}m region well, but not the 10 {\mu}m feature; (ii) if we include the UV-to-near-IR SED in the fit, models do not fit the 2-8 {\mu}m region. The observed 10 {\mu}m features are broader and shallower than those in the best-fit models in the first approach. In the second case, the shape of the 10 {\mu}m feature is better reproduced by the Clumpy models. The additional blackbody contribution in the 2-8 {\mu}m range allows Clumpy models dominated by cooler temperatures (T < 800K) to better fit the 8-12{\mu}m SED. A centrally concentrated distribution of a small number of torus clouds is required in the first case, while in the second case the clouds are more spread out radially. The temperature of the blackbody component is ~ 1200 K as expected for graphite grains.

The Wide-field Infrared Survey Explorer (WISE) has scanned the entire sky with unprecedented sensitivity in four infrared bands, at 3.4, 4.6, 12, and 22 micron. The WISE Point Source Catalog contains more than 560 million objects, among them hundreds of thousands of galaxies with Active Nuclei (AGN). While type 1 AGN, owing to their bright and unobscured nature, are easy to detect and constitute a rather complete and unbiased sample, their type 2 counterparts, postulated by AGN unification, are not as straightforward to identify. Matching the WISE catalog with known QSOs in the Sloan Digital Sky Survey we confirm previous identification of the type 1 locus in the WISE color space. Using a very large database of the popular CLUMPY torus models, we find the colors of the putative type 2 counterparts, and also, for the first time, predict their number vs. flux relation that can be expected to be observed in any given WISE color range. This will allow us to put statistically very significant constraints on the torus parameters. Our results are a successful test of the AGN unification scheme., 4 pages, 1 figure, presented at the IAU Symposium #304