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.

We present a formalism for continuum and line emission from random clumpy media together with its application to problems of current interest, including CO spectral lines from ensembles of clouds and radio emission from HII regions, supernovae and star-forming regions. For line emission we find that the effects of clump opacity on observed line ratios can be indistinguishable from variations of intrinsic line strengths, adding to the difficulties in determining abundances from line observations. Our formalism is applicable to arbitrary distributions of cloud properties, provided the cloud volume filling factor is small; numerical simulations show it to hold up to filling factors of about 10%. We show that irrespective of the complexity of the cloud ensemble, the radiative effect of clumpiness can be parametrized at each frequency by a single multiplicative correction to the overall optical depth; this multiplier is derived from appropriate averaging over individual cloud properties. Our main finding is that cloud shapes have only a negligible effect on radiation propagation in clumpy media; the results of calculations employing point-like clouds are practically indistinguishable from those for finite-size clouds with arbitrary geometrical shapes.

We present MOLPOP-CEP, a universal line transfer code that allows the exact calculation of multi-level line emission from a slab with variable physical conditions for any arbitrary atom or molecule for which atomic data exist. The code includes error control to achieve any desired level of accuracy, providing full confidence in its results. Publicly available, MOLPOP-CEP employs our recently developed Coupled Escape Probability (CEP) technique, whose performance exceeds other exact methods by orders of magnitude. The program also offers the option of an approximate solution with different variants of the familiar escape probability method. As an illustration of the MOLPOP-CEP capabilities we present an exact calculation of the Spectral Line Energy Distribution (SLED) of the CO molecule and compare it with escape probability results. We find that the popular large-velocity gradient (LVG) approximation is unreliable at large CO column densities. Providing a solution of the multi-level line transfer problem at any prescribed level of accuracy, MOLPOP-CEP is removing any doubts about the validity of its final results.

We present models of deeply buried ultraluminous infrared galaxy (ULIRG) spectral energy distributions (SEDs) and use them to construct a three-dimensional diagram for diagnosing the nature of observed ULIRGs. Our goal is to construct a suite of SEDs for a very simple model ULIRG structure, and to explore how well this simple model can (by itself) explain the full range of observed ULIRG properties. We use our diagnostic to analyze archival Spitzer Space Telescope IRS spectra of ULIRGs and find that: (1) In general, our model does provide a comprehensive explanation of the distribution of mid-IR ULIRG properties; (2) >75% (in some cases 100%) of the bolometric luminosities of the most deeply buried ULIRGs must be powered by a dust-enshrouded active galactic nucleus; (3) an unobscured "keyhole" view through <~10% of the obscuring medium surrounding a deeply buried ULIRG is sufficient to make it appear nearly unobscured in the mid-IR; and (4) the observed absence of deeply buried ULIRGs with large PAH equivalent widths is naturally explained by our models showing that deep absorption features are "filled-in" by small quantities of foreground unobscured PAH emission (e.g., from the host galaxy disk) at the level of ~1% the bolometric nuclear luminosity. The modeling and analysis we present will also serve as a powerful tool for interpreting the high angular resolution spectra of high-redshift sources to be obtained with the James Webb Space Telescope.

We present mid-infrared (MIR, 7.5-13.5 $\mu$m) imaging and spectroscopy observations obtained with the CanariCam (CC) instrument on the 10.4m Gran Telescopio CANARIAS for a sample of 20 nearby, MIR bright and X-ray luminous QSOs. We find that for the majority of QSOs the MIR emission is unresolved at angular scales nearly 0.3 arcsec, corresponding to physical scales $<600$ pc. We find that the higher-spatial resolution CC spectra have similar shapes to those obtained with Spitzer/IRS, and hence we can assume that the spectra are not heavily contaminated by extended emission in the host galaxy. We thus take advantage of the higher signal to noise Spitzer/IRS spectra, as a fair representation of the nuclear emission, to decompose it into a combination of active galactic nuclei (AGN), polycyclic aromatic hydrocarbon (PAH) and stellar components. In most cases the AGN is the dominant component, with a median contribution of 85 per cent of the continuum light at MIR (5-15 $\mu$m) within the IRS slit. This IR AGN emission is well reproduced by clumpy torus models. We find evidence for significant differences in the parameters that describe the dusty tori of QSOs when compared with the same parameters of Seyfert 1 and 2 nuclei. In particular, we find a lower number of clouds ($N_{0}<12$), steeper radial distribution of clouds ($q=1.5-3.0$), and clouds that are less optically thick ($\tau_{V}<100$) than in Seyfert 1, which could be attributed to dusty structures that have been partially evaporated and piled up by the higher radiation field in QSOs. We find that the combination of the angular width $\sigma_{torus}$, viewing angle $i$, and number of clouds along the equatorial line $N_{0}$, produces large escape probabilities ($P_{esc} > 2$ per cent) and low geometrical covering factors ($f_{2}<0.6$), as expected for AGN with broad lines in their optical spectra.

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.

The absence of intrinsic broad line emission has been reported in a number of active galactic nuclei (AGN), including some with high Eddington ratios. Such "true type 2 AGN" are inherent to the disk-wind scenario for the broad line region: Broad line emission requires a minimal column density, implying a minimal outflow rate and thus a minimal accretion rate. Here we perform a detailed analysis of the consequences of mass conservation in the process of accretion through a central disk. The resulting constraints on luminosity are consistent with all the cases where claimed detections of true type 2 AGN pass stringent criteria, and predict that intrinsic broad line emission can disappear at luminosities as high as about 4x$10^{46}$ erg s$^{-1}$ and any Eddington ratio, though more detections can be expected at Eddington ratios below about 1%. Our results are applicable to every disk outflow model, whatever its details and whether clumpy or smooth, irrespective of the wind structure and its underlying dynamics. While other factors, such as changes in spectral energy distribution or covering factor, can affect the intensities of broad emission lines, within this scenario they can only produce true type 2 AGN of higher luminosity then those prescribed by mass conservation.

We present an analysis of the nuclear infrared (IR, 1.6 to 18 $\mu$m) emission of the ultraluminous IR galaxy UGC 5101 to derive the properties of its active galactic nucleus (AGN) and its obscuring material. We use new mid-IR high angular resolution ($0.3-0.5$ arcsec) imaging using the Si-2 filter ($\lambda_{C}=8.7\, \mu$m) and $7.5-13$ $\mu$m spectroscopy taken with CanariCam (CC) on the 10.4m Gran Telescopio CANARIAS. We also use archival HST/NICMOS and Subaru/COMICS imaging and Spitzer/IRS spectroscopy. We estimate the near- and mid-IR unresolved nuclear emission by modelling the imaging data with GALFIT. We decompose the Spitzer/IRS and CC spectra using a power-law component, which represents the emission due to dust heated by the AGN, and a starburst component, both affected by foreground extinction. We model the resulting unresolved near- and mid-IR, and the starburst subtracted CC spectrum with the CLUMPY torus models of Nenkova et al. The derived geometrical properties of the torus, including the large covering factor and the high foreground extinction needed to reproduce the deep $9.7\, \mu$m silicate feature, are consistent with the lack of strong AGN signatures in the optical. We derive an AGN bolometric luminosity $L_{bol}\sim1.9\times10^{45}\,$erg s$^{-1}$ that is in good agreement with other estimates in the literature.

We present results from the fitting of infrared (IR) spectral energy distributions of 21 active galactic nuclei (AGN) with clumpy torus models. We compiled high spatial resolution ($\sim 0.3$--$0.7$ arcsec) mid-IR $N$-band spectroscopy, $Q$-band imaging and nuclear near- and mid-IR photometry from the literature. Combining these nuclear near- and mid-IR observations, far-IR photometry and clumpy torus models, enables us to put constraints on the torus properties and geometry. We divide the sample into three types according to the broad line region (BLR) properties; type-1s, type-2s with scattered or hidden broad line region (HBLR) previously observed, and type-2s without any published HBLR signature (NHBLR). Comparing the torus model parameters gives us the first quantitative torus geometrical view for each subgroup. We find that NHBLR AGN have smaller torus opening angles and larger covering factors than those of HBLR AGN. This suggests that the chance to observe scattered (polarized) flux from the BLR in NHBLR could be reduced by the dual effects of (a) less scattering medium due to the reduced scattering volume given the small torus opening angle and (b) the increased torus obscuration between the observer and the scattering region. These effects give a reasonable explanation for the lack of observed HBLR in some type-2 AGN.

We present a multiwavelength study of the OH Megamaser galaxy (OHMG) IRAS16399-0937, based on new HST/ACS F814W and H$\alpha$+[NII] images and archive data from HST, 2MASS, Spitzer, Herschel and the VLA. This system has a double nucleus, whose northern (IRAS16399N) and southern (IRAS16399S) components have a projected separation of $\sim$ 6'' (3.4 kpc) and have previously been identified based on optical spectra as a Low Ionization Nuclear Emission Line Region (LINER) and starburst nucleus, respectively. The nuclei are embedded in a tidally distorted common envelope, in which star formation is mostly heavily obscured. The infrared spectrum is dominated by strong polycyclic aromatic hydrocarbon (PAH), but deep silicate and molecular absorption features are also present, and are strongest in the IRAS16399N nucleus. The 0.435 - 500$\mu$m SED was fitted with a model including stellar, ISM and AGN torus components using our new MCMC code, clumpyDREAM. The results indicate that the IRAS16399N contains an AGN (L$_{bol} \sim 10^{44}$ ergs/s) deeply embedded in a quasi-spherical distribution of optically-thick clumps with a covering fraction $\approx1$. We suggest that these clumps are the source of the OHM emission in IRAS16399-0937. The high torus covering fraction precludes AGN-photoionization as the origin of the LINER spectrum, however, the spectrum is consistent with shocks (v $\sim100-200$ km s$^{-1}$). We infer that the $\sim10^8$ M$_{\odot}$ black-hole in IRAS16399N is accreting at a small fraction ($\sim1$%) of its Eddington rate. The low accretion-rate and modest nuclear SFRs suggest that while the gas-rich major merger forming the IRAS16399-0937 system has triggered widespread star formation, the massive gas inflows expected from merger simulations have not yet fully developed.

Through matches with the Sloan Digital Sky Survey (SDSS) catalogue we identify the location of various families of astronomical objects in WISE colour space. We identify reliable indicators that separate Galactic/local from extragalactic sources and concentrate here on the objects in our Galaxy and its closest satellites. We develop colour and magnitude criteria that are based only on WISE data to select asymptotic giant branch (AGB) stars with circumstellar dust shells, and separate them into O-rich and C-rich classes. With these criteria we produce an all-sky map for the count ratio of the two populations. The map reveals differences between the Galactic disc, the Magellanic Clouds and the Sgr Dwarf Spheroidal galaxy, as well as a radial gradient in the Large Magellanic Cloud (LMC) disc. We find that the C:O number ratio for dusty AGB stars increases with distance from the LMC centre about twice as fast as measured for near-IR selected samples of early AGB stars. Detailed radiative transfer models show that WISE colours are well explained by the emission of centrally heated dusty shells where the dust has standard properties of interstellar medium (ISM) grains. The segregation of different classes of objects in WISE colour space arises from differences in properties of the dust shells: those around young stellar objects have uniform density distributions while in evolved stars they have steep radial profiles.

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.

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.

Apart from viewing-dependent obscuration, intrinsic broad-line emission from active galactic nuclei (AGNs) follows an evolutionary sequence: Type $1 \to 1.2/1.5 \to 1.8/1.9 \to 2$ as the accretion rate onto the central black hole is decreasing. This spectral evolution is controlled, at least in part, by the parameter $L_{\rm bol}/M^{2/3}$, where $L_{\rm bol}$ is the AGN bolometric luminosity and $M$ is the black hole mass. Both this dependence and the double-peaked profiles that emerge along the sequence arise naturally in the disk-wind scenario for the AGN broad-line region.

We report the results of a full-Stokes survey of all four 18 cm OH lines in 77 OH megamasers (OHMs) using the Arecibo Observatory. This is the first survey of OHMs that included observations of the OH satellite lines; only 4 of the 77 OHMs have existing satellite line observations in the literature. In 5 sources, satellite line emission is detected, with 3 of the 5 sources re-detections of previously published sources. The 2 sources with new detections of satellite line emission are IRAS F10173+0829, which was detected at 1720 MHz, and IRAS F15107+0724, for which both the 1612 MHz and 1720 MHz lines were detected. In IRAS F15107+0724, the satellite lines are partially conjugate, as 1720 MHz absorption and 1612 MHz emission have the same structure at some velocities within the source, along with additional broader 1612 MHz emission. This is the first observed example of conjugate satellite lines in an OHM. In the remaining sources, no satellite line emission is observed. The detections and upper limits are generally consistent with models of OHM emission in which all of the 18 cm OH lines have the same excitation temperature. There is no evidence for a significant population of strong satellite line emitters among OHMs.

We present a model in which the 22 GHz H$_2$O masers observed in star-forming regions occur behind shocks propagating in dense regions (preshock density $n_0 \sim 10^6 - 10^8$ cm$^{-3}$). We focus on high-velocity ($v_s > 30$ km s$^{-1}$) dissociative J shocks in which the heat of H$_2$ re-formation maintains a large column of $\sim 300-400$ K gas; at these temperatures the chemistry drives a considerable fraction of the oxygen not in CO to form H$_2$O. The H$_2$O column densities, the hydrogen densities, and the warm temperatures produced by these shocks are sufficiently high to enable powerful maser action. The observed brightness temperatures (generally $\sim 10^{11} - 10^{14}$ K) are the result of coherent velocity regions that have dimensions in the shock plane that are 10 to 100 times the shock thickness of $\sim 10^{13}$ cm. The masers are therefore beamed towards the observer, who typically views the shock "edge-on", or perpendicular to the shock velocity; the brightest masers are then observed with the lowest line of sight velocities with respect to the ambient gas. We present numerical and analytic studies of the dependence of the maser inversion, the resultant brightness temperature, the maser spot size and shape, the isotropic luminosity, and the maser region magnetic field on the shock parameters and the coherence path length; the overall result is that in galactic H$_2$O 22 GHz masers these observed parameters can be produced in J shocks with $n_0\sim 10^6 - 10^8$ cm$^{-3}$ and $v_s \sim 30 -200$ km s$^{-1}$. A number of key observables such as maser shape, brightness temperature, and global isotropic luminosity depend only on the particle flux into the shock, $j=n_0v_s$, rather than on $n_0$ and $v_s$ separately.

An optically and geometrically thick torus obscures the central engine of Active Galactic Nuclei (AGN) from some lines of sight. From a magnetohydrodynamical framework, the torus can be considered to be a particular region of clouds surrounding the central engine where the clouds are dusty and optically thick. In this framework, the magnetic field plays an important role in the creation, morphology and evolution of the torus. If the dust grains within the clouds are assumed to be aligned by paramagnetic alignment, then the ratio of the intrinsic polarisation and visual extinction, P(%)/Av, is a function of the magnetic field strength. To estimate the visual extinction through the torus and constrain the polarisation mechanisms in the nucleus of AGN, we developed a polarisation model to fit both the total and polarised flux in a 1.2" (~263pc) aperture of the type 2 AGN, IC5063. The polarisation model is consistent with the nuclear polarisation observed at K being produced by dichroic absorption from aligned dust grains with a visual extinction through the torus of 48$\pm$2 mag. We estimated the intrinsic polarisation arising from dichroic absorption to be P$_{K}^{dic}$=12.5$\pm$2.7%. We consider the physical conditions and environment of the gas and dust for the torus of IC5063. Then, through paramagnetic alignment, we estimate a magnetic field strength in the range of 12-128mG in the NIR emitting regions of the torus of IC5063. Alternatively, we estimate the magnetic field strength in the plane of the sky using the Chandrasekhar-Fermi method. The minimum magnetic field strength in the plane of the sky is estimated to be 13 and 41 mG depending of the conditions within the torus of IC5063. These techniques afford the chance to make a survey of AGN, to investigate the effects of magnetic field strength on the torus, accretion, and interaction to the host galaxy.

We present observations of the molecular gas in the nuclear environment of three prototypical low luminosity AGN (LLAGN), based on VLT/SINFONI AO-assisted integral-field spectroscopy of H2 1-0 S(1) emission at angular resolutions of ~0.17". On scales of 50-150 pc the spatial distribution and kinematics of the molecular gas are consistent with a rotating thin disk, where the ratio of rotation (V) to dispersion (sigma) exceeds unity. However, in the central 50 pc, the observations reveal a geometrically and optically thick structure of molecular gas (V/sigma<1 and N_H>10^{23} cm^{-2}) that is likely to be associated with the outer extent of any smaller scale obscuring structure. In contrast to Seyfert galaxies, the molecular gas in LLAGN has a V/sigma<1 over an area that is ~9 times smaller and column densities that are in average ~3 times smaller. We interpret these results as evidence for a gradual disappearance of the nuclear obscuring structure. While a disk wind may not be able to maintain a thick rotating structure at these luminosities, inflow of material into the nuclear region could provide sufficient energy to sustain it. In this context, LLAGN may represent the final phase of accretion in current theories of torus evolution. While the inflow rate is considerable during the Seyfert phase, it is slowly decreasing, and the collisional disk is gradually transitioning to become geometrically thin. Furthermore, the nuclear region of these LLAGN is dominated by intermediate-age/old stellar populations (with little or no on-going star formation), consistent with a late stage of evolution.

We present high-resolution mid-infrared (MIR) imaging, nuclear spectral energy distributions (SEDs) and archival Spitzer spectra for 22 low-luminosity active galactic nuclei (LLAGN; Lbol \lesssim 10^42 erg/sec). Infrared (IR) observations may advance our understanding of the accretion flows in LLAGN, the fate of the obscuring torus at low accretion rates, and, perhaps, the star formation histories of these objects. However, while comprehensively studied in higher-luminosity Seyferts and quasars, the nuclear IR properties of LLAGN have not yet been well-determined. We separate the present LLAGN sample into three categories depending on their Eddington ratio and radio emission, finding different IR characteristics for each class. (I) At the low-luminosity, low-Eddington ratio (log Lbol/LEdd < -4.6) end of the sample, we identify "host-dominated" galaxies with strong polycyclic aromatic hydrocarbon bands that may indicate active (circum-)nuclear star formation. (II) Some very radio-loud objects are also present at these low Eddington ratios. The IR emission in these nuclei is dominated by synchrotron radiation, and some are likely to be unobscured type 2 AGN that genuinely lack a broad line region. (III) At higher Eddington ratios, strong, compact nuclear sources are visible in the MIR images. The nuclear SEDs of these galaxies are diverse; some resemble typical Seyfert nuclei, while others lack a well-defined MIR "dust bump". Strong silicate emission is present in many of these objects. We speculate that this, together with high ratios of silicate strength to hydrogen column density, could suggest optically thin dust and low dust-to-gas ratios, in accordance with model predictions that LLAGN do not host a Seyfert-like obscuring torus.

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: 1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; 2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image reconstruction process was "blind" in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with 6 to 10 telescopes, should be well suited to provide images of complex sceneries.

The inevitable spread in properties of the toroidal obscuration of active galactic nuclei (AGNs) invalidates the widespread notion that type 1 and 2 AGNs are intrinsically the same objects, drawn randomly from the distribution of torus covering factors. Instead, AGNs are drawn \emph{preferentially} from this distribution; type 2 are more likely drawn from the distribution higher end, type 1 from its lower end. Type 2 AGNs have a higher IR luminosity, lower narrow-line luminosity and a higher fraction of Compton thick X-ray obscuration than type 1. Meaningful studies of unification statistics cannot be conducted without first determining the intrinsic distribution function of torus covering factors.

Double-peaked line profiles are commonly considered a hallmark of rotating disks, with the distance between the peaks a measure of the rotation velocity. However, double-peaks can arise also from radiative transfer effects in optically thick non-rotating sources. Utilizing exact solutions of the line transfer problem we present a detailed study of line emission from geometrically thin Keplerian disks. We derive the conditions for emergence of kinematic double peaks in optically thin and thick disks, and find that it is generally impossible to disentangle the effects of kinematics and line opacity in observed double-peaked profiles. Unless supplemented by additional information, a double-peaked profile alone is not a reliable indicator of a rotating disk. In certain circumstances, triple and quadruple profiles might be better indicators of rotation in optically thick disks.

We have obtained high-resolution mid-infrared (MIR) imaging, nuclear spectral energy distributions (SEDs) and archival Spitzer spectra for 22 low-luminosity active galactic nuclei (LLAGN; L_bol < 5 x 10^42 erg/s). Infrared (IR) observations may advance our understanding of the accretion flows in LLAGN, the fate of the obscuring torus at low accretion rates, and, perhaps, the star formation histories of these objects. However, while comprehensively studied in higher-luminosity Seyferts and quasars, the nuclear IR properties of LLAGN have not yet been well-determined. In these proceedings we summarise the results for the LLAGN at the relatively high-luminosity, high-Eddington ratio end of the sample. Strong, compact nuclear sources are visible in the MIR images of these objects, with luminosities consistent with or slightly in execss of that predicted by the standard MIR/X-ray relation. Their broadband nuclear SEDs are diverse; some resemble typical Seyfert nuclei, while others possess less of a well-defined MIR ``dust bump''. Strong silicate emission is present in many of these objects. We speculate that this, together with high ratios of silicate strength to hydrogen column density, could suggest optically thin dust and low dust-to-gas ratios, in accordance with model predictions that LLAGN do not host a Seyfert-like obscuring torus.

We used the CLUMPY torus models and a Bayesian approach to fit the infrared spectral energy distributions (SEDs) and ground-based high-angular resolution mid-infrared spectroscopy of 13 nearby Seyfert galaxies. This allowed us to put tight constraints on torus model parameters such as the viewing angle, the radial thickness of the torus Y, the angular size of the cloud distribution sigma_torus, and the average number of clouds along radial equatorial rays N_0. The viewing angle is not the only parameter controlling the classification of a galaxy into a type 1 or a type 2. In principle type 2s could be viewed at any viewing angle as long as there is one cloud along the line of sight. A more relevant quantity for clumpy media is the probability for an AGN photon to escape unabsorbed. In our sample, type 1s have relatively high escape probabilities, while in type 2s, as expected, tend to be low. Our fits also confirmed that the tori of Seyfert galaxies are compact with torus model radii in the range 1-6pc. The scaling of the models to the data also provided the AGN bolometric luminosities, which were found to be in good agreement with estimates from the literature. When we combined our sample of Seyfert galaxies with a sample of PG quasars from the literature to span a range of L_bol(AGN)~10^{43}-10^{47}erg/s, we found plausible evidence of the receding torus. That is, there is a tendency for the torus geometrical covering factor to be lower at high AGN luminosities than at low AGN luminosities. This is because at low AGN luminosities the tori appear to have wider angular sizes and more clouds along radial equatorial rays. We cannot, however rule out the possibility that this is due to contamination by extended dust structures not associated with the dusty torus at low AGN luminosities, since most of these in our sample are hosted in highly inclined galaxies. (Abridged)

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.

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 nature of maser emission means that the apparent angular size of an individual maser spot is determined by the amplification process as well as by the instrinsic size of the emitting cloud. Highly sensitive MERLIN radio interferometry images spatially and spectrally resolve water maser clouds around evolved stars. We measured the properties of clouds around the red supergiant S Per and the AGB stars IK Tau, RT Vir, U Her and U Ori, to test maser beaming theory. Spherical clouds are expected to produce an inverse relationship between maser intensity and apparent size, which would not be seen from cylindrical or slab-like regions. We analysed the maser properties, in order to estimate the saturation state, and investigated the variation of observed spot size with intensity and across the spectral line profiles. Circumstellar masers emanate from discrete clouds from about one to 20 AU in diameter depending on the star. Most of the maser features have negative excitation temperatures close to zero and modest optical depths, showing that they are mainly unsaturated. Around S Per and (at most epochs) RT Vir and IK Tau, the maser component size shrinks with increasing intensity. In contrast, the masers around U Ori and U Her tend to increase in size, with a larger scatter. The water masers from S Per, RT Vir and IK Tau are mainly beamed into spots with an observed angular size much smaller than the emitting clouds and smallest of all at the line peaks. This suggests that the masers are amplification-bounded, emanating from approximately spherical clouds. Many of the masers around U Her and U Ori have apparent sizes which are more similar to the emitting clouds and have less or no dependence on intensity, suggesting that these masers are matter-bounded. This is consistent with an origin in flattened clouds and these two stars have shown other behaviour indicating the presence of shocks.

After recent sensitivity upgrades at the Keck Interferometer (KI), systematic interferometric 2um studies of the innermost dust in nearby Seyfert nuclei are within observational reach. Here, we present the analysis of new interferometric data of NGC 4151, discussed in context of the results from recent dust reverberation, spectro-photometric and interferometric campaigns. The complete data set gives a complex picture, in particular the measured visibilities from now three different nights appear to be rather insensitive to the variation of the nuclear luminosity. KI data alone indicate two scenarios: the K-band emission is either dominated to ~90% by size scales smaller than 30mpc, which falls short of any dust reverberation measurement in NGC 4151 and of theoretical models of circum-nuclear dust distributions. Or contrary, and more likely, the K-band continuum emission is dominated by hot dust (>= 1300K) at linear scales of about 50mpc. The linear size estimate varies by a few tens of percent depending on the exact morphology observed. Our interferometric, deprojected centro-nuclear dust radius estimate of 55+-5mpc is roughly consistent with the earlier published expectations from circum-nuclear, dusty radiative transfer models, and spectro-photometric modeling. However, our data do not support the notion that the dust emission size scale follows the nuclear variability of NGC 4151 as a R_dust \propto L_nuc^0.5 scaling relation. Instead variable nuclear activity, lagging, and variable dust response to illumination changes need to be combined to explain the observations.

The mid-far-infrared spectral energy distributions (SEDs) of 83 active galaxies, mostly Seyfert galaxies, selected from the extended 12 micron sample are presented. The data were collected using all three instruments, IRAC, IRS, and MIPS, aboard the Spitzer Space Telescope. The IRS data were obtained in spectral mapping mode, and the photometric data from IRAC and IRS were extracted from matched, 20 arcsec diameter circular apertures. The MIPS data were obtained in SED mode, providing very low resolution spectroscopy (R ~ 20) between ~ 55 and 90 microns in a larger, 20 by 30 arcsec synthetic aperture. We further present the data from a spectral decomposition of the SEDs, including equivalent widths and fluxes of key emission lines; silicate 10 and 18 micron emission and absorption strengths; IRAC magnitudes; and mid-far infrared spectral indices. Finally, we examine the SEDs averaged within optical classifications of activity. We find that the infrared SEDs of Seyfert 1s and Seyfert 2s with hidden broad line regions (HBLR, as revealed by spectropolarimetry or other technique) are qualitatively similar, except that Seyfert 1s show silicate emission and HBLR Seyfert 2s show silicate absorption. The infrared SEDs of other classes with the 12 micron sample, including Seyfert 1.8-1.9, non-HBLR Seyfert 2 (not yet shown to hide a type 1 nucleus), LINER and HII galaxies, appear to be dominated by star-formation, as evidenced by blue IRAC colors, strong PAH emission, and strong far-infrared continuum emission, measured relative to mid-infrared continuum emission.

We compare observations of AGB stars and predictions of the Elitzur & Ivezic (2001) steady-state radiatively driven dusty wind model. The model results are described by a set of similarity functions of a single independent variable, and imply general scaling relations among the system parameters. We find that the model properly reproduces various correlations among the observed quantities and demonstrate that dust drift through the gas has a major impact on the structure of most winds. From data for nearby oxygen-rich and carbon-rich mass-losing stars we find that (1) the dispersion in grain properties within each group is rather small; (2) both the dust cross-section per gas particle and the dust-to-gas mass ratio are similar for the two samples even though the stellar atmospheres and grain properties are very different; (3) the dust abundance in both outflows is significantly below the Galactic average, indicating that most of the Galactic dust is not stardust - contrary to popular belief, but in support of Draine (2009). Our model results can be easily applied to recent massive data sets, such as the Spitzer SAGE survey of the Large Magellanic Cloud, and incorporated in galaxy evolution models.

We present an analysis of Spitzer IRS spectroscopy of 83 active galaxies from the extended 12 micron sample. We find rank correlations between several tracers of star formation which suggest that (1) the PAH feature is a reliable tracer of star formation, (2) there is a significant contribution to the heating of the cool dust by stars, (3) the H$_2$ emission is also primarily excited by star formation. The 55-90 vs. 20-30 spectral index plot is also a diagnostic of the relative contribution of Starburst to AGN. We see there is a large change in spectral index across the sample. Thus, the contribution to the IR spectrum from the AGN and starburst components can be comparable in magnitude but the relative contribution also varies widely across the sample. We find rank correlations between several AGN tracers. We search for correlations between AGN and Starburst tracers and we conclude that the AGN and Starburst tracers are not correlated. This is consistent with our conclusion that the relative strength of the AGN and Starburst components varies widely across the sample. Thus, there is no simple link between AGN fueling and Black Hole Growth and star formation in these galaxies. The distribution of Sil 10 micron and 18 micron strengths is consistent with the clumpy torus models of Sirocky et al. We find a rank correlation between the [NeV] 14 micron line and the 6.7 micron continuum which may be due to an extended component of hot dust. The Sy 2s with a Hidden Broad Line Region (HBLR) have a higher ratio of AGN to Starburst contribution to the SED than Sy 2s without an HBLR. This may contribute to the detection of the HBLR in polarized light. The Sy 2s with an HBLR are more similar to the Sy 1s than they are to the Sy 2s without an HBLR.

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 <35% of objects with tori sharing characteristics and geometry similar to the best fit would have their central engines obscured. This relatively low obscuration probability can explain the clear appearance of the 10-mic emission feature in SST1721+6012 together with its rarity among other QSO2. (2) We also fitted the SED of PG1211+143, one of the first type 1 QSOs with a 10-mic silicate feature in emission. Among similar sources, this QSO appears to display an unusually broadened feature whose peak is shifted toward longer wavelengths. Although this led to suggestions of non-standard dust chemistry in these sources, our analysis fits such SEDs with standard galactic dust; the apparent peak shifts arise from radiative transfer effects. (3) We find that the distribution of silicate feature strengths among clumpy torus models closely resembles the observed distribution, and the feature never occurs deeply absorbed. (abridged)

We present the first detections of the 3um water ice and 3.4um amorphous hydrocarbon (HAC) absorption features in z~2 ULIRGs. These are based on deep rest-frame 2-8um Spitzer IRS spectra of 11 sources selected for their appreciable silicate absorption. The HAC-to-silicate ratio for our z~2 sources is typically higher by a factor of 2-5 than that observed in the Milky Way. This HAC `excess' suggests compact nuclei with steep temperature gradients as opposed to predominantly host obscuration. Beside the above molecular absorption features, we detect the 3.3um PAH emission feature in one of our sources with three more individual spectra showing evidence for it. Stacking analysis suggests that water ice, hydrocarbons, and PAH are likely present in the bulk of this sample even when not individually detected. The most unexpected result of our study is the lack of clear detections of the 4.67um CO gas absorption feature. Only three of the sources show tentative signs of this feature and at significantly lower levels than has been observed in local ULIRGs. Overall, we find that the closest local analogs to our sources, in terms of 3-4um color, HAC-to-silicate and ice-to-silicate ratios, as well as low PAH equivalent widths are sources dominated by deeply obscured nuclei. Such sources form only a small fraction of ULIRGs locally and are commonly believed to be dominated by buried AGN. Our sample suggests that, in absolute number, such buried AGN are at least an order of magnitude more common at z~2 than today. The presence of PAH suggests that significant levels of star-formation are present even if the obscured AGN typically dominate the power budget.

The disk-wind scenario for the broad-line region (BLR) and toroidal obscuration in active galactic nuclei predicts the disappearance of the BLR at low luminosities. In accordance with the model predictions, data from a nearly complete sample of nearby AGNs show that the BLR disappears at luminosities lower than $5\times\E{39} (M/10^7\Mo)^{2/3}$ erg s$^{-1)$, where $M$ is the black hole mass. The radiative efficiency of accretion onto the black hole is $\la \E{-3}$ for these sources, indicating that their accretion is advection-dominated.

We fitted Spitzer/IRS ~2-35 um spectra of 26 luminous QSOs in attempt to define the main emission components. Our model has three major components: a clumpy torus, dusty narrow line region (NLR) clouds and a blackbody-like dust. The models utilize the clumpy torus of Nenkova et al. (2008) and are the first to allow its consistent check in type-I AGNs. Single torus models and combined torus-NLR models fail to fit the spectra of most sources but three component models adequately fit the spectra of all sources. We present torus inclination, cloud distribution, covering factor and torus mass for all sources and compare them with bolometric luminosity, black hole mass and accretion rate. The torus covering factor and mass are found to be correlated with the bolometric luminosity of the sources. We find that a substantial amount of the ~2-7 um radiation originates from a hot dust component, which likely situated in the innermost part of the torus. The luminosity radiated by this component and its covering factor are comparable to those of the torus. We quantify the emission by the NLR clouds and estimate their distance from the center. The distances are ~700 times larger than the dust sublimation radius and the NLR covering factor is about 0.07. The total covering factor by all components is in good agreement with the known AGN type-I:type-II ratio.

A strong signature of a circumstellar disc around a high-mass protostar has been inferred from high resolution methanol maser observations in NGC7538-IRS1 N. This interpretation has however been challenged with a bipolar outflow proposed as an alternative explanation. We compare the two proposed scenarios for best consistency with the observations. Using a newly developed formalism we model the optical depth of the maser emission at each observed point in the map and LOS velocity for the two scenarios. We find that if the emission is symmetric around a central peak in both space and LOS velocity then it has to arise from an edge-on disc in sufficiently fast differential rotation. Disc models successfully fit ~100 independent measurement points in position-velocity space with 4 free parameters to an overall accuracy of 3-4%. Solutions for Keplerian rotation require a central mass of at least 4 solar masses. Close to best-fitting models are obtained if Keplerian motion is assumed around a central mass equaling ~30 solar masses as inferred from other observations. In contrast we find that classical bipolar outflow models cannot fit the data, although could be applicable in other sources. Our results strongly favour the differentially rotating disc hypothesis to describe the main feature of the 12.2 (and 6.7) GHz methanol maser emission in NGC7538 IRS1 N. Furthermore, for Keplerian rotation around a ~30 solar masses protostar we predict the position and velocity at which tangentially amplified masers should be detected in high dynamic range observations. [abridged]

Toroidal obscuration is a keystone of AGN unification. There is now direct evidence for the torus emission in infrared, and possibly water masers. Here I summarize the torus properties, its possible relation to the immediate molecular environment of the AGN and present some speculations on how it might evolve with the AGN luminosity.

The infrared is a key wavelength regime for probing the dusty, obscured nuclear regions of active galaxies. We present results from an infrared study of 87 nearby Seyfert galaxies using the Spitzer Space Telescope and ground-based telescopes. Combining detailed modelling of the 3 - 100 micron spectral energy distributions with mid-IR spectral diagnostics and near-infrared observations, we find broad support for the unified model of AGNs. The IR emission of Seyfert 1s and 2s is consistent with their having the same type of central engine viewed at a different orientation. The nature of the putative torus is becoming clearer; in particular we present evidence that it is likely a clumpy medium. Mid-infrared correlations between tracers of star formation and AGN ionizing luminosity reveal the starburst-AGN connection implied by the black hole/bulge mass relation, however it is not yet clear if this is due to feedback.

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.

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.

Observations give strong support for the unification scheme of active galactic nuclei. The scheme is premised on toroidal obscuration of the central engine by dusty clouds that are individually very optically thick. These lectures summarize the torus properties, describe the handling and implications of its clumpy nature and present speculations about its dynamic origin.

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.

OH megamasers (OHMs) emit primarily in the main lines at 1667 and 1665 MHz, and differ from their Galactic counterparts due to their immense luminosities, large linewidths and 1667/1665 MHz flux ratios, which are always greater than one. We find that these maser properties result from strong 53 micron radiative pumping combined with line overlap effects caused by turbulent linewidths of about 20 km/s; pumping calculations that do not include line overlap are unreliable. A minimum dust temperature of about 45 K is needed for inversion, and maximum maser efficiency occurs for dust temperatures in the range 80 - 140 K. We find that warmer dust can support inversion at lower IR luminosities, in agreement with observations. Our results are in good agreement with a clumpy model of OHMs, with clouds sizes about 1 pc and OH column densities about 5e16 cm^2, that is able to explain both the diffuse and compact emission observed for OHMs. We suggest that all OH main line masers may be pumped by far-IR radiation, with the major differences between OHMs and Galactic OH masers caused by differences in linewidth produced by line overlap. Small Galactic maser linewidths tend to produce stronger 1665 MHz emission. The large OHM linewidths lead to inverted ground state transitions having approximately the same excitation temperature, producing 1667/1665 MHz flux ratios greater than one and weak satellite line emission. Finally, the small observed ratio of pumping radiation to dense molecular gas, as traced by HCN and HCO$^+$, is a possible reason for the lack of OH megamaser emission in NGC 6240.

The collisionally pumped, ground-state 1720 MHz maser line of OH is widely recognized as a tracer for shocked regions and observed in star forming regions and supernova remnants. Whereas some lines of excited states of OH have been detected and studied in star forming regions, the subject of excited-state OH in supernova remnants -- where high collision rates are to be expected -- is only recently being addressed. Modeling of collisional excitation of OH demonstrates that 1720, 4765 and 6049 MHz masers can occur under similar conditions in regions of shocked gas. In particular, the 6049 and 4765 MHz masers become more significant at increased OH column densities where the 1720 MHz masers begin to be quenched. In supernova remnants, the detection of excited-state OH line maser emission could therefore serve as a probe of regions of higher column densities. Using the Very Large Array, we searched for excited-state OH in the 4.7, 7.8, 8.2 and 23.8 GHz lines in four well studied supernova remnants with strong 1720 MHz maser emission (SgrAEast, W28, W44 and IC443). No detections were made, at typical detection limits of around 10 mJy/beam. The search for the 6 GHz lines were done using Effelsberg since the VLA receivers did not cover those frequencies, and are reported on in an accompanying letter (Fish, Sjouwerman & Pihlstrom 2007). We also cross-correlated the positions of known supernova remnants with the positions of 1612 MHz maser emission obtained from blind surveys. No probable associations were found, perhaps except in the SgrAEast region. The lack of detections of excited-state OH indicates that the OH column densities suffice for 1720 MHz inversion but not for inversion of excited-state transitions, consistent with the expected results for C-type shocks.

We present Spitzer Space Telescope mid-infrared spectra of 12 Seyfert 1.8 and 1.9 galaxies over the 5-38 um region. We compare the spectral characteristics of this sample to those of 58 Seyfert 1 and Seyfert 2 galaxies from the Spitzer archives. An analysis of the spectral shapes, the silicate 10 um feature and the emission line fluxes have enabled us to characterize the mid-IR properties of Seyfert 1.8/1.9s. We find that the equivalent widths of the 10 um silicate feature are generally weak in all Seyfert galaxies, as previously reported by several studies. The few Seyfert galaxies in this sample that show deep 10 um silicate absorption features are highly inclined and/or merging galaxies. It is likely that these absorption features originate primarily in the dusty interstellar medium of the host galaxy rather than in a dusty torus on parsec scales close to the central engine. We find that the equivalent width of the polycyclic aromatic hydrocarbon (PAH) band at 6.2 um correlates strongly with the 20-30 um spectral index. Either of these quantities are good indicators of the amount of starburst contribution to the mid-IR spectra. The spectra of Seyfert 1.8 and 1.9s are dominated by these starburst features, similar to most Seyfert 2s. They show strong PAH bands and a strong red continuum toward 30 um. The strengths of the high-ionization forbidden narrow emission lines [O IV] 25.89 um, [Ne III] 15.56 um and [Ne V] 14.32 um relative to [Ne II] 12.81 um are weaker in Seyfert 1.8/1.9s and Seyfert 2s as compared to Seyfert 1s. The weakness of high-ionization lines in Seyfert 1.8-1.9s is suggestive of intrinsically weak active galactic nuclei (AGN) continua, and/or stronger star formation activity leading to enhanced [Ne II]. We discuss the implications of these observational results in the context of the Unified Model of AGN.

%Context: {Previous studies have indicated that the 372.4 GHz ground transition of ortho-H$_2$D$^+$ might be a powerful probe of Proto-Planetary Disks. The line could be especially suited for study of the disk mid-plane, where the bulk of the mass resides and where planet formation takes place.} %Aims: {Provide detailed theoretical predictions for the line intensity, profile and maps expected for representative disk models.} %Methods: {We determine the physical and chemical structure of the disks from the model developed by Ceccarelli & Dominik (2005). The line emission is computed with the new radiative transfer method developed recently by Elitzur & Asensio Ramos (2006).} %Results: {We present intensity maps convolved with the expected ALMA resolution, which delineate the origin of the H$_2$D$^+$ 372.4 GHz line. In the disk inner regions, the line probes the conditions in the mid-plane out to radial distances of a few tens of AU, where Solar-like planetary systems might form. In the disk outermost regions, the line originates from slightly above the mid-plane. When the disk is spatially resolved, the variation of line profile across the image provides important information about the velocity field. Spectral profiles of the entire disk flux show a double peak shape at most inclination angles.} %Conclusions: {Our study confirms that the 372.4 GHz H$_2$D$^+$ line provides powerful diagnostics of the mid-plane of Proto-Planetary Disks. Current submillimeter telescopes are capable of observing this line, though with some difficulties. The future ALMA interferometer will have the sensitivity to observe and even spatially resolve the H$_2$D$^+$ line emission.}

This review covers selected developments in maser theory since the previous meeting, "Cosmic Masers: From Proto-Stars to Black Holes" (Migenes & Reid 2002). Topics included are time variability of fundamental constants, pumping of OH megamasers and indicators for differentiating disks from bi-directional outflows.

We discuss Subaru and Spitzer Space Telescope imaging and spectroscopy of M87 in the mid-infrared from 5-35 um. These observations allow us to investigate mid-IR emission mechanisms in the core of M87 and to establish that the flaring, variable jet component HST-1 is not a major contributor to the mid-IR flux. The Spitzer data include a high signal-to-noise 15-35 $\mu$m spectrum of the knot A/B complex in the jet, which is consistent with synchrotron emission. However, a synchrotron model cannot account for the observed {\it nuclear} spectrum, even when contributions from the jet, necessary due to the degrading of resolution with wavelength, are included. The Spitzer data show a clear excess in the spectrum of the nucleus at wavelengths longer than 25 um, which we model as thermal emission from cool dust at a characteristic temperature of 55 \pm 10 K, with an IR luminosity \sim 10^{39} {\rm ~erg ~s^{-1}}. Given Spitzer's few-arcsecond angular resolution, the dust seen in the nuclear spectrum could be located anywhere within ~5'' (390 pc) of the nucleus. In any case, the ratio of AGN thermal to bolometric luminosity indicates that M87 does not contain the IR-bright torus that classical unified AGN schemes invoke. However, this result is consistent with theoretical predictions for low-luminosity AGNs

We present diffraction limited, 10um imaging polarimetry data for the central regions of the archetypal Seyfert AGN, NGC1068. The position angle of polarization is consistent with three dominant polarizing mechanisms. We identify three distinct regions of polarization: (a) north of the nucleus, arising from aligned dust in the NLR, (b) south, east and west of the nucleus, consistent with dust being channeled toward the central engine and (c) a central minimum of polarization consistent with a compact (<22pc) torus. These observations provide continuity between the geometrically and optically thick torus and the host galaxy's nuclear environments. These images represent the first published mid-IR polarimetry from an 8-m class telescope and illustrate the potential of such observations.

We present an analysis of the effects of luminosity on the shape of the mid-infrared spectral energy distributions (SEDs) of 234 radio-quiet quasars originally presented by Richards et al. In quasars without evident dust extinction, the spectrally integrated optical and infrared luminosities are linearly correlated over nearly three decades in luminosity. We find a significant (>99.99% confidence) correlation between the 1.8-8.0 micron spectral index and infrared luminosity that indicates an enhancement of the mid-infrared continuum with increasing luminosity. Coupled with strong evidence for spectral curvature in more luminous quasars, we conclude this trend is likely a manifestation of the `near-infrared (3-5 micron) bump' noticed in earlier quasar SED surveys. The strength of this feature is indicative of the contribution of emission from the hottest (>1000 K) dust to the mid-infrared spectrum; higher luminosity quasars tend to show more hot dust emission. Finally, the comparable distribution of bolometric corrections from the monochromatic 3 micron luminosity as well as its lack of sensitivity to dust extinction as compared to the standard bolometric correction from nu*L_{5100A} suggest that the former may be a more robust indicator of bolometric quasar luminosity. The close link between the power in the mid-infrared and optical and the effect of luminosity on the shape of the mid-infrared continuum indicate that considering mid-infrared emission independent of the properties of the quasar itself is inadequate for understanding the parsec-scale quasar environment.

The nucleus of the nearby galaxy, NGC 1097, is known to host a young, compact (r < 9 parsec) nuclear star cluster as well as a low-luminosity active galactic nucleus (AGN). It has been suggested both that the nuclear stellar cluster is associated with a dusty torus, and that low-luminosity AGN like NGC 1097 do not have the torus predicted by the unified model of AGN. To investigate these contradictory possibilities we have acquired Gemini/T-ReCS 11.7 micron and 18.3 micron images of the central few hundred parsecs of this galaxy at < 45 parsec angular resolution, in which the nucleus and spectacular, kiloparsec-scale star-forming ring are detected in both bands. The small-scale mid-infrared (mid-IR) luminosity implies thermal emission from warm dust close to the central engine of this galaxy. Fitting of torus models shows that the observed mid-IR emission cannot be accounted for by dust heated by the central engine. Rather, the principal source heating the dust in this object is the nuclear star cluster itself, suggesting that the dust that we detect is not the torus of AGN unified schemes (although it is also possible that the dusty starburst itself could provide the obscuration invoked by the unified model). Comparison of Spitzer/IRS and Gemini/GNIRS spectra shows that, although polyaromatic hydrocarbon emission (PAH) bands are strong in the immediate circumnuclear region of the galaxy, PAH emission is weak or absent in the central 19 parsecs. The lack of PAH emission can probably be explained largely by destruction/ionization of PAH molecules by hard photons from the nuclear star cluster. If NGC 1097 is typical, PAH emission bands may not be a useful tool with which to find very compact nuclear starbursts even in low-luminosity AGN.

A sample of 196 AGNs and ULIRGs observed by the Infrared Spectrograph (IRS) on Spitzer is analyzed to study the distribution of the strength of the 9.7 micron silicate feature. Average spectra are derived for quasars, Seyfert 1 and Seyfert 2 AGNs, and ULIRGs. We find that quasars are characterized by silicate features in emission and Seyfert 1s equally by emission or weak absorption. Seyfert 2s are dominated by weak silicate absorption, and ULIRGs are characterized by strong silicate absorption (mean apparent optical depth about 1.5). Luminosity distributions show that luminosities at rest frame 5.5 micron are similar for the most luminous quasars and ULIRGs and are almost 10^5 times more luminous than the least luminous AGN in the sample. The distributions of spectral characteristics and luminosities are compared to those of optically faint infrared sources at z~2 being discovered by the IRS, which are also characterized by strong silicate absorption. It is found that local ULIRGs are a similar population, although they have lower luminosities and somewhat stronger absorption compared to the high redshift sources.

Observations give strong support for the unification scheme of active galactic nuclei. Clumpiness of the toroidal obscuration is crucial for explaining the IR observations and has significance consequences for AGN classification: type 1 and type 2 viewing is an angle-dependent probability, not an absolute property. The broad line region (BLR) and the dusty torus are, respectively, the inner and outer segments, across the dust sublimation radius, of a continuous cloud distribution. Continuum X-ray obscuration comes mostly from the inner, BLR clouds. All clouds are embedded in a disk wind, whose mass outflow rate is diminishing as the accretion rate, i.e., AGN luminosity, is decreasing. The torus disappears when $L \la$ \E{42} erg/sec, the BLR at some lower, yet to be determined luminosities.

We present a new diagnostic diagram for mid-infrared spectra of infrared galaxies based on the equivalent width of the 6.2 micron PAH emission feature and the strength of the 9.7 micron silicate feature. Based on the position in this diagram we classify galaxies into 9 classes ranging from continuum-dominated AGN hot dust spectra and PAH-dominated starburst spectra to absorption-dominated spectra of deeply obscured galactic nuclei. We find that galaxies are systematically distributed along two distinct branches: one of AGN and starburst-dominated spectra and one of deeply obscured nuclei and starburst-dominated spectra. The separation into two branches likely reflects a fundamental difference in the dust geometry in the two sets of sources: clumpy versus non-clumpy obscuration. Spectra of ULIRGs are found along the full length of both branches, reflecting the diverse nature of the ULIRG family.

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.

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central engine. The compact sizes (only a few pc) determined in recent high-resolution observations require that the obscuring matter be clumpy and located inside the region where the black-hole gravity dominates over the galactic bulge. This location is in line with the scenario depicting the torus as the region of the clumpy wind coming off the accretion disk in which the clouds are dusty and optically thick. We study here the outflow scenario within the framework of hydromagnetic disk winds, incorporating the cloud properties determined from detailed modeling of the IR emission from clumpy tori. We find that torus clouds were likely detected in recent water maser observations of NGC 3079. In the wind scenario, the AGN main dynamic channel for release of accreted mass seems to be switching at low luminosities from torus outflow to radio jets. The torus disappears when the bolometric luminosity decreases below about \E{42} erg/sec because the accretion onto the central black hole can no longer sustain the required cloud outflow rate. This disappearance seems to have been observed in both LINERs and radio galaxies. With further luminosity decrease, suppression of cloud outflow spreads radially inward from the disk's dusty, molecular region into its atomic, ionized zone, resulting in disappearance of the broad emission line region at lower luminosities, yet to be determined.

New Spitzer imaging observations have revealed the structure around the Mira variable star R Hya to be a one-sided parabolic arc 100 arcsec to the West stretching from North to South. We successfully model R Hya and its surroundings in terms of an interaction of the stellar wind from an asymptotic giant branch (AGB) star with the interstellar medium (ISM) the star moves through. Our three-dimensional hydrodynamic simulation reproduces the structure as a bow shock into the oncoming ISM. We propose this as another explanation of detached shells around such stars which should be considered alongside current theories of internal origin. The simulation predicts the existence of a tail of ram-pressure-stripped AGB material stretching downstream. Indications for such a tail behind R Hya are seen in IRAS maps.

We present the first results of the MIRIAD (MIPS [Multiband Imaging Photometer for Spitzer] Infra-Red Imaging of AGB [asymptotic giant branch] Dustshells) project using the Spitzer Space Telescope. The primary aim of the project is to probe the material distribution in the extended circumstellar envelopes (CSE) of evolved stars and recover the fossil record of their mass loss history. Hence, we must map the whole of the CSEs plus the surrounding sky for background subtraction, while avoiding the central star that is brighter than the detector saturation limit. With our unique mapping strategy, we have achieved better than one MJy/sr sensitivity in three hours of integration and successfully detected a faint (< 5 MJy/sr), extended (~400 arcsec) far-infrared nebula around the AGB star R Hya. Based on the parabolic structure of the nebula, the direction of the space motion of the star with respect to the nebula shape, and the presence of extended H alpha emission co-spatial to the nebula, we suggest that the detected far-IR nebula is due to a bow shock at the interface of the interstellar medium and the AGB wind of this moving star. This is the first detection of the stellar-wind bow-shock interaction for an AGB star and exemplifies the potential of Spitzer as a tool to examine the detailed structure of extended far-IR nebulae around bright central sources. \

NGC 7538 IRS1 is a high-mass (approx. 30 M_sun) protostar with a CO outflow, an associated UCHII region, and a linear methanol maser structure, which might trace a Keplerian-rotating circumstellar disk. The directions of the various associated axes are misaligned with each other. We investigate the near-infrared morphology of the source to clarify the relations among the various axes. K'-band bispectrum speckle interferometry was performed at two 6-meter-class telescopes -- the BTA 6m telescope and the 6.5m MMT. Complementary IRAC images from the Spitzer Space Telescope Archive were used to relate the structures detected with the outflow at larger scales. High-dynamic range images show fan-shaped outflow structure in which we detect 18 stars and several blobs of diffuse emission. We interpret the misalignment of various outflow axes in the context of a disk precession model, including numerical hydrodynamic simulations of the molecular emission. The precession period is approx. 280 years and its half-opening angle is 40 degrees. A possible triggering mechanism is non-coplanar tidal interaction of an (undiscovered) close companion with the circumbinary protostellar disk. Our observations resolve the nearby massive protostar NGC 7538 IRS2 as a close binary with separation of 195 mas. We find indications for shock interaction between the outflow activities in IRS1 and IRS2. Indications of outflow precession have been discovered to date in a number of massive protostars, all with large precession angles 20--45 degrees. This might explain the difference between the outflow widths in low- and high-mass stars and add support to a common collimation mechanism.

We are conducting a large observing program with the Spitzer Space Telescope to determine the mid-to-far infrared spectral energy distributions of a well-defined sample of 87 nearby, 12 micron-selected Seyfert galaxies. In this paper we present the results of IRS low-resolution spectroscopy of a statistically representative subsample of 51 of the galaxies (59%), with an analysis of the continuum shapes and a comparison of the Seyfert types. We find that the spectra clearly divide into groups based on their continuum shapes and spectral features. Some spectral features are clearly related to a starburst contribution to the IR spectrum, while the mechanisms producing observed power-law continuum shapes, attributed to an AGN component, may be dust or non-thermal emission. The infrared spectral types appear to be related to the Seyfert types. Principal component analysis results suggest that the relative contribution of starburst emission may be the dominant cause of variance in the observed spectra. We find that Sy 1's have higher ratios of IR/radio emission than Sy 2's, as predicted by the unified model if the torus is optically thick in the mid-IR. However, smooth-density torus models predict a much larger difference between type 1's and 2's than observed in our sample. Our observations may be consistent with clumpy torus models containing a steep radial distribution of optically thick dense clumps. (Abridged)

We present the results of IRS low-resolution spectroscopy of 51 Seyfert galaxies, part of a large Spitzer observing program to determine the mid-to-far infrared spectral energy distributions of a well-defined sample of 87 nearby, 12 micron-selected Seyferts. We find that the spectra clearly divide into groups based on their continuum shapes and spectral features. The infrared spectral types appear to be related to the Seyfert types. Some features are clearly related to a starburst contribution to the IR spectrum, while the observed power-law continuum shapes, attributed to the AGN, may be dust or non-thermal emission. Principal component analysis results suggest that the relative contribution of starburst emission is the dominant cause of variance in the spectra. We find that the Sy 2's show on average stronger starburst contributions than the Sy 1's.

We present spatially-resolved, near-diffraction-limited 10 micron spectra of the nucleus of the Seyfert 2 galaxy NGC 1068, obtained with Michelle, the mid-IR imager and spectrometer on the 8.1 m Gemini North telescope. The spectra cover the nucleus and the central 6.0" x 0.4" of the ionization cones at a spatial resolution of approximately 0.4" (approx. 30 parsecs). The spectra extracted in 0.4" steps along the slit reveal striking variations in continuum slope, silicate feature profile and depth, and fine structure line fluxes on subarcsecond scales, illustrating in unprecedented detail the complexity of the circumnuclear regions of this galaxy at mid-IR wavelengths. A comparison of photometry in various apertures reveals two distinct components: a compact (radius <15 pc), bright source within the central 0.4" x 0.4" and extended, lower brightness emission. We identify the compact source with the AGN obscuring torus, and the diffuse component with the AGN-heated dust in the ionization cones. While the torus emission dominates the flux observed in the near-IR, the mid-IR flux measured with apertures larger than about 1" is dominated instead by the dust emission from the ionization cones. Many previous attempts to determine the torus spectral energy distribution are thus likely to be significantly affected by contamination from the extended emission. The observed spectrum of the compact source is compared with clumpy torus models, which require most of the mid-IR emitting clouds to be located within a few parsecs of the central engine. We also present a UKIRT/CGS4 5 micron spectrum covering the R(0) -- R(4) lines of the fundamental vibration-rotation band of 12CO. None of these lines was detected, and we discuss these non-detections in terms of the filling factor and composition of the nuclear clouds. (Abridged)

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central engine. Torus sizes of hundreds of parsecs were deduced from early theoretical modeling efforts, but high-resolution IR observations now show that the torus size is no more than a few parsecs. This conflict is resolved when the clumpy nature of the torus is taken into account. The compact torus may be best understood when identified with the dusty, optically thick region of the wind coming off the central accretion disk.

We present a new method, the Coupled Escape Probability (CEP), for exact calculation of line emission from multi-level systems, solving only algebraic equations for the level populations. The CEP formulation of the classical two-level problem is a set of {\em linear equations}, and we uncover an exact analytic expression for the emission from two-level optically thick sources that holds as long as they are in the "effectively thin" regime. In comparative study of a number of standard problems, the CEP method outperformed the leading line transfer methods by substantial margins. The algebraic equations employed by our new method are already incorporated in numerous codes based on the escape probability approximation. All that is required for an exact solution with these existing codes is to augment the expression for the escape probability with simple zone-coupling terms. As an application, we find that standard escape probability calculations generally produce the correct cooling emission by the CII 158 mic line but not by the $^3$P lines of OI.

We examine the ``puffed-up inner disk'' model (Dullemond, Dominik & Natta 2001), proposed for explaining the near-IR excess radiation from Herbig Ae/Be stars. Detailed model computations show that the observed near-IR excess requires more hot dust than is contained in the puffed-up disk rim. The rim can produce the observed near-IR excess only if its dust has perfectly gray opacity, but such dust is in conflict with the observed 10$\mu$m spectral feature. We find that a compact ($\sim$ 10 AU) tenuous ($\tau_V \la 0.4$) dusty halo around the disk inner regions contains enough dust to readily explain the observations. Furthermore, this model also resolves the puzzling relationship noted by Monnier and Millan-Gabet (2002) between luminosity and the interferometric inner radii of disks.

We model the OH megamaser emission from the luminous infrared galaxy IIIZw35 as arising from a narrow rotating starburst ring of radius 22 pc enclosing a mass of 7 x 10^6 Msun. We show how both the compact and apparently diffuse maser emission from this ring can arise from a single phase of unsaturated maser clouds amplifying background radio continuum. The masing clouds are estimated to have a diameter of <0.7 pc and internal velocity dispersion of about 20 km/s. We find that the clouds are neither self-gravitating nor pressure confined but are freely expanding. Their dispersal lifetimes may set the vertical thickness of the ring. For an estimated internal density of 3 x 10^3 cm^-3, cloud masses are of order 24 Msun. The observed spectral features and velocity gradients indicate that the clouds must be outflowing and escaping the nucleus. The cloud mass outflow rate is estimated to be 0.8 Msun/yr, while the star formation rate is about 19 Msun/yr. Associated ionised gas, possibly generated from dissipated clouds, provides free-free absorption along the source axis, explaining the observed East-West asymmetries. We show that the clumpiness of a maser medium can have a dramatic effect on what is observed even in a relatively low gain OH megamaser. Specifically, in IIIZw35 our clumpy maser model naturally explains the large line to continuum ratios, the large 1667MHz:1665MHz line ratios and the wide velocity dispersions seen in the compact maser spots. Other astrophysical masers showing both compact and apparently diffuse emission might be explained by similar clumpy structures.

We present Odin observations of the AGB star W Hya in the ground-state transition of ortho-H2O, 1(10)-1(01), at 557GHz. The line is clearly of circumstellar origin. Radiative transfer modelling of the water lines observed by Odin and ISO results in a mass-loss rate of (2.5 +/- 0.5)E-7 Msol/yr, and a circumstellar H2O abundance of (2.0 +/- 1.0)E-3. The inferred mass-loss rate is consistent with that obtained from modelling the circumstellar CO radio line emission, and also with that obtained from the dust emission modelling combined with a dynamical model for the outflow. The very high water abundance, higher than the cosmic oxygen abundance, can be explained by invoking an injection of excess water from evaporating icy bodies in the system. The required extra mass of water is quite small, on the order of ~ 0.1 earth mass.

We present an edge-on Keplerian disk model to explain the main component of the 12.2 and 6.7 GHz methanol maser emission detected toward NGC7538-IRS1 N. The brightness distribution and spectrum of the line of bright masers are successfully modeled with high amplification of background radio continuum emission along velocity coherent paths through a maser disk. The bend seen in the position-velocity diagram is a characteristic signature of differentially rotating disks. For a central mass of 30 solar masses, suggested by other observations, our model fixes the masing disk to have inner and outer radii of about 270 AU and 750 AU.

We present a model for the OH megamaser emission in the starburst galaxy IIIZw35. The observed diffuse and compact OH maser components in this source are explained by a single phase of unsaturated clumpy gas distributed in a thin ring structure and amplifying background continuum. We emphasize the importance of clumpiness in the OH masing medium, an effect that has not been fully appreciated previously. The model explains why multiple bright spots are seen only at the ring tangents while smoother emission is found elsewhere. Both the observed velocity gradients and the line to continuum ratios around the ring enquire a geometry where most of the seed photons come from a continuum emission which lies outside the OH ring. To explain both the OH and continuum brightness, free-free absorbing gas is required along the ring axis to partially absorb the far side of the ring. It is proposed that the required geometry arises from an inwardly propagating ring of starburst activity.

We discuss the propagation of spectral line and continuum radiation in a clumpy medium and give general expressions for the observed absorption or emission from a cloud population. We show that the affect of the medium clumpiness can usually be characterised by a single number multiplying the mean column opacity. Our result provides a simpler proof and generalisation of the result of Martin et al (1984). The formalism provides a simple way to understand the effects of clumping on molecular line profiles and ratios, for example how clumping effects the interpretation of 13CO(1-0) to 12CO(1-0) line ratios. It can also be used as a propagation operator in physical models of clumpy media where the incident radiation effects the spectral line emissivity. We are working to extend the formalism to the propagation of masers in a clumpy medium, but in this case there are special difficulties because formal expectation values are not characteristic of observations because they are biased by rare events.

Near-IR imaging of the AGB star IRC+10011 (= CIT3) reveals the presence of a bipolar structure within the central ~0.1 arcsec of a spherical dusty wind. We show that the image asymmetries originate from ~1E-4 Msun of swept-up wind material in an elongated cocoon whose expansion is driven by bipolar jets. We perform detailed 2D radiative transfer calculations with the cocoon modeled as two cones extending to ~1,100 AU within an opening angle of ~30deg, imbedded in a wind with the standard r^{-2} density profile. The cocoon expansion started <~200 years ago, while the total lifetime of the circumstellar shell is ~5,500 years. Similar bipolar expansion, at various stages of evolution, has been recently observed in a number of other AGB stars, culminating in jet breakout from the confining spherical wind. The bipolar outflow is triggered at a late stage in the evolution of AGB winds, and IRC+10011 provides its earliest example thus far. These new developments enable us to identify the first instance of symmetry breaking in the evolution from AGB to planetary nebula.

We present an edge-on Keplerian disc model to explain the main component of the 12.2 and 6.7 GHz methanol maser emission detected toward NGC7538-IRS1 N. The brightness distribution and spectrum of the line of bright masers are successfully modeled with high amplification of background radio continuum emission along velocity coherent paths through a maser disc. The bend seen in the position-velocity diagram is a characteristic signature of differentially rotating discs. For a central mass of 30 solar masses, suggested by other observations, our model fixes the masing disc to have inner and outer radii of about 350 AU and 1000 AU.

Bipolar expansion at various stages of evolution has been recently observed in a number of AGB stars. The expansion is driven by bipolar jets that emerge late in the evolution of AGB winds. The wind traps the jets, resulting in an expanding, elongated cocoon. Eventually the jets break-out from the confining spherical wind, as recently observed in W43A. This source displays the most advanced evolutionary stage of jets in AGB winds. The earliest example is IRC+10011, where the asymmetry is revealed in high-resolution near-IR imaging. In this source the jets turned on only ~200 years ago, while the spherical wind is ~4000 years old.

High-mass star formation is not well understood chiefly because examples are deeply embedded, relatively distant, and crowded with sources of emission. Using VLA and VLBA observations of water and SiO maser emission, we have mapped in detail the structure and proper motion of material 20-500 AU from the closest high-mass YSO, radio source-I in the Orion KL region. We observe streams of material driven in a rotating, wide angle, bipolar wind from the surface of an edge-on accretion disk. The example of source-I provides strong evidence that high-mass star formation proceeds via accretion

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central source, giving rise to type 1 line spectrum for pole-on viewing and type 2 characteristics in edge-on sources. Infrared radiation at its different wavelengths is the best probe of the dust distribution, whether the torus orientation is edge-on or pole-on. The observed IR is in broad agreement with the unified scheme but serious problems remained in all early modeling efforts. In spite of a general awareness that the dust must be concentrated in clouds, clumpiness remained the one major ingredient missing in those radiative transfer studies because of the inherent difficulties it presents. We have recently developed the formalism to handle dust clumpiness and our results indicate that its inclusion may resolve the difficulties encountered by the previous theoretical efforts. We show that these problems find a natural explanation if the dust is contained in \about\ 5--10 clouds along radial rays through the torus. The spectral energy distributions (SED) of both type 1 and type 2 sources are properly reproduced from different viewpoints of the same object if the optical depth of each cloud is \ga 40 at visual wavelengths and the clouds' mean free path increases roughly in proportion to radial distance.

We study the IR emission from flared disks with and without additional optically thin halos. Flux calculations of a flared disk in vacuum can be considered a special case of the more general family of models in which the disk is imbedded in an optically thin halo. In the absence of such halo, flux measurements can never rule out its existence because the disk flaring surface defines a mathematically equivalent halo that produces the exact same flux at all IR wavelengths. When a flared disk with height H at its outer radius R is imbedded in a halo whose optical depth at visual wavelengths is \tau\sub{halo}, the system IR flux is dominated by the halo whenever \tau\sub{halo} > (1/4)H/R. Even when its optical depth is much smaller, the halo can still have a significant effect on the disk temperature profile. Imaging is the only way to rule out the existence of a potential halo, and we identify a decisive test that extracts a signature unique to flared disks from imaging observations.

We report on single-dish monitoring and extremely high angular resolution observations of the flaring H2O megamaser in the Seyfert 2 galaxy Mrk 348. The H2O line is redshifted by ~130 km/s with respect to the systemic velocity, is very broad, with a FWHM of 130 km/s, and has no detectable high velocity components within 1500 km/s on either side of the strong line. Monitoring observations made with the Effelsberg 100m telescope show that the maser varies significantly on timescales as short as one day and that the integrated line flux is loosely correlated with the continuum flux. VLBA observations indicate that the maser emission arises entirely from a region less than 0.25 pc in extent, located toward a continuum component thought to be associated with the receding jet. We also report on integrated flux monitoring with the VLA between 1.4 and 43 GHz, and VLBA continuum observations of the milliarcsecond scale jets at 1.7, 8, 15, and 22 GHz. These observations have allowed us to tentatively pinpoint the location of the core, and also show the ejection of a new jet component during the current radio "flare."

Most dusty winds are described by a set of similarity functions of a single independent variable that can be chosen as the overall optical depth at visual. The self-similarity implies general scaling relations among the system parameters, in agreement with observations. Dust drift through the gas has a major impact on the structure of most winds.

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central source, giving rise to Seyfert 1 line spectrum for pole-on viewing and Seyfert 2 characteristics in edge-on sources. Although the observed IR is in broad agreement with this scheme, the behavior of the 10 micron silicate feature and the width of the far-IR emission peak remained serious problems in all previous modeling efforts. We show that these problems find a natural explanation if the dust is contained in about 5-10 clouds along radial rays through the torus. The spectral energy distributions (SED) of both type 1 and type 2 sources are properly reproduced from different viewpoints of the same object if the visual optical depth of each cloud is larger than about 60 and the clouds' mean free path increases roughly in proportion to radial distance.

We address the dusty wind problem, from the point where dust formation has been completed and outward. Given grain properties, both radiative transfer and hydrodynamics components of the problem are fully defined by four additional input parameters. The wind radiative emission and the shape of its velocity profile are both independent of the actual magnitude of the velocity and are determined by just three dimensionless free parameters. Of the three, only one is always significant---for most of phase space the solution is described by a set of similarity functions of a single independent variable, which can be chosen as the overall optical depth at visual \tV. The self-similarity implies general scaling relations among mass loss rate (\Mdot), luminosity ($L$) and terminal velocity (\vf). Systems with different \Mdot, $L$ and \vf but the same combination $\Mdot/L^{3/4}$ necessarily have also the same $\Mdot\vf/L$. ... Eliminating \tV produces $\vf^3 = A \Mdot(1 + B \Mdot^{4/3}/L)^{-1.5}$, where $A$ and $B$ are coefficients that contain the only dependence of this universal correlation on chemical composition. At a given $L$, the maximal velocity of a dusty wind is $v\sub{max} \propto L^{1/4}$ attained at $\Mdot \propto L^{3/4}$, with proportionality coefficients derived from $A$ and $B$.

This review covers some of the developments in basic theory of astronomical masers over the past ten years. Topics included are the effects of three dimensional geometry and polarization, with special emphasis on the differences between maser and non-maser radiation.

We construct a self-consistent model for the wind around W Hya by solving the coupled equations describing the hydrodynamics and dust radiative transfer problems. The model matches simultaneously the observed continuum radiation and wind velocity profile. The water line emission is calculated next using the water abundance as the only free parameter, fitted from the ISO observations of Neufeld et al. (1996) and Barlow et al. (1996). The gas temperature is determined from a thermal balance calculation that includes water as one of its main components. Our model successfully fits all the observed water lines, resolving a major discrepancy between the modeling results of the two observing teams. The mass loss rate is 2.3 x 10^{-6} M_solar yr^{-1}, the water abundance is 1.0 x 10^{-4} and the ortho:para ratio is 1:1.3.

Unbiased analysis shows that IRAS data reliably differentiate between the early and late stages of stellar evolution because objects at these stages clearly segregate in infrared color-color diagrams. Structure in these diagrams is primarily controlled by the density distribution of circumstellar dust. The density profile around older objects is the steepest, declining as $r^{-2}$, while young objects have profiles that vary as $r^{-3/2}$ and flatter. The different density profiles reflect the different dynamics that govern the different environments. Our analysis also shows that high mass star formation is strongly concentrated within \about 5 kpc around the Galactic center, in support of other studies.

DUSTY solves the problem of radiation transport in a dusty environment. The code can handle both spherical and planar geometries. The user specifies the properties of the radiation source and dusty region, and the code calculates the dust temperature distribution and the radiation field in it. The solution method is based on a self-consistent equation for the radiative energy density, including dust scattering, absorption and emission, and does not introduce any approximations. The solution is exact to within the specified numerical accuracy. DUSTY has built in optical properties for the most common types of astronomical dust and comes with a library for many other grains. It supports various analytical forms for the density distribution, and can perform a full dynamical calculation for radiatively driven winds around AGB stars. The spectral energy distribution of the source can be specified analytically as either Planckian or broken power-law. In addition, arbitrary dust optical properties, density distributions and external radiation can be entered in user supplied files. Furthermore, the wavelength grid can be modified to accommodate spectral features. A single DUSTY run can process an unlimited number of models, with each input set producing a run of optical depths, as specified. The user controls the detail level of the output, which can include both spectral and imaging properties as well as other quantities of interest.

IR and mm-wave emission from Herbig Ae/Be stars has produced conflicting conclusions regarding the dust geometry in these objects. We show that the compact dimensions of the mm-wave emitting regions are a decisive indication for disks. But a disk cannot explain the spectral energy distribution (SED) unless it is embedded in an extended envelope that (1) dominates the IR emission and (2) provides additional disk heating on top of the direct stellar radiation. Detailed radiative transfer calculations based on the simplest model for envelope-embedded disks successfully fit the data from UV to mm wavelengths and show that the disks have central holes. This model also resolves naturally some puzzling results of IR imaging.

Recent spectroscopic and VLBI-imaging observations of bright extragalactic water maser sources have revealed that the megamaser emission often originates in thin circumnuclear disks near the centers of active galactic nuclei (AGNs). Using general radiative and kinematic considerations and taking account of the observed flux variability, we argue that the maser emission regions are clumpy, a conclusion that is independent of the detailed mechanism (X-ray heating, shocks, etc.) driving the collisionally pumped masers. We examine scenarios in which the clumps represent discrete gas condensations (i.e., clouds) and do not merely correspond to velocity irregularities in the disk. We show that even two clouds that overlap within the velocity coherence length along the line of sight could account (through self-amplification) for the entire maser flux of a high-velocity ``satellite'' feature in sources like NGC 4258 and NGC 1068, and we suggest that cloud self-amplification likely contributes also to the flux of the background-amplifying ``systemic'' features in these objects. Analogous interpretations have previously been proposed for water maser sources in Galactic star-forming regions. We argue that this picture provides a natural explanation of the time-variability characteristics of extragalactic megamaser sources and of their apparent association with Seyfert 2-like galaxies. We also show that the requisite cloud space densities and internal densities are consistent with the typical values of nuclear (broad emission-line region-type) clouds.

Recent observations show that the OH 1720 MHz maser is a powerful probe of the shocked region where a supernova remnant strikes a molecular cloud. We perform a thorough study of the pumping of this maser and find tight constraints on the physical conditions needed for its production. The presence of the maser implies moderate temperatures (50 -- 125 K) and densities ($\sim 10^5 cm^{-3}$), and OH column densities of order $10^{16} cm^{-2}$. We show that these conditions can exist only if the shocks are of C-type. J-shocks fail by such a wide margin that the presence of this maser could become the most powerful indicator of C-shocks. These conditions also mean that the 1720 MHz maser will be inherently weak compared to the other ground state OH masers. All the model predictions are in good agreement with the observations.

We present a detailed, self-consistent model of radiatively driven stellar outflows which couples the radiative transfer and hydrodynamics equations. The circumstellar envelope, which consists of gas and dust, is described as a two-component fluid to account for relative drifts. Our results agree with both molecular line observations and infrared continuum spectra, and show that steady-state outflows driven by radiation pressure on dust grains adequately describe the surroundings of the majority of cool luminous evolved stars.

Profile comparison of the Stokes parameters $V$ and $I$ is a powerful tool for maser data analysis, providing the first direct methods for unambiguous determination of (1) the maser saturation stage, (2) the amplification optical depth and intrinsic Doppler width of unsaturated masers, and (3) the comparative magnitudes of Zeeman splitting and Doppler linewidth. Circular polarization recently detected in OH 1720 MHz emission from the Galactic center appears to provide the first direct evidence for maser saturation.

Dust infrared emission possesses scaling properties. Overall luminosity is never an input parameter of the radiative transfer problem, spectral shape is the only relevant property of the heating radiation when the inner boundary of the dusty region is controlled by dust sublimation. Similarly, the absolute scales of densities and distances are irrelevant; the geometry enters only through angles, relative thicknesses and aspect ratios, and the actual magnitudes of densities and distances enter only through one independent parameter, the overall optical depth. Dust properties enter only through dimensionless, normalized distributions that describe the spatial variation of density and the wavelength dependence of scattering and absorption efficiencies. Scaling enables a systematic approach to modeling and classification of IR spectra. We develop a new, fully scale-free method for solving radiative transfer, present exact numerical results, and derive approximate analytical solutions for spherical geometry, covering the entire range of parameter space relevant to observations. Scaling implies tight correlations among the SEDs of various members of the same class of sources such as young stellar objects, late-type stars, etc. In particular, all members of the same class occupy common, well defined regions in color-color diagrams. The observational data corroborate the existence of these correlations.

The spectral shape of IR emission from Herbig Ae/Be stars has been invoked as evidence for accretion disks around high-mass protostars. Instead, we present here models based on spherical envelopes with $r^{-1.5}$ dust density profile that successfully explain the observed spectral shapes. The spectral energy distributions (SEDs) of eight primary candidates for protostellar disks are fitted in detail for all wavelengths available, from visual to far IR. The only envelope property adjusted in individual sources is the overall visual optical depth, and it ranges from 0.3 to 3. In each case, our models properly reproduce the data for both IR excess, visual extinction and reddening. The success of our models shows that accretion disks cannot make a significant contribution to the radiation observed in these pre-main sequence stars.

Infrared emission from radiatively heated dust possesses general scaling properties. The spectral shape is independent of overall luminosity when the inner boundary of the dusty region is controlled by dust sublimation; the only relevant property of the heating radiation is its frequency profile. For a given type of dust grains, the emission from red giants and supergiants is essentially controlled by a single parameter --- the overall dust optical depth. This leads to tight correlations among different spectral properties that explain many available observations and enable systematic studies of large data bases.

Infrared imaging properties of dusty winds around late-type stars are investigated in detail, employing a self-consistent model that couples the equations of motion and radiative transfer. Because of general scaling properties, the angular profiles of surface brightness are self-similar. In any given star, the profile shape is determined essentially by overall optical depth at each wavelength and it is self-similarly scaled by the size of the dust condensation zone. We find that the mid-IR is the best wavelength range to measure directly the angular size of this zone, and from {\it IRAS} data we identify the 15 best candidates for such future observations. We also show that the visibility function at short wavelengths ($\la$ 2 \mic) directly determines the scattering optical depth, and produce theoretical visibility curves for various characteristic wavelengths and the entire parameter range relevant to late-type stars. The infrared emission should display time variability because of cyclical changes in overall optical depth, reflecting luminosity-induced movement of the dust condensation point. Calculations of the wavelength dependence of photometric amplitudes and time variability of envelope sizes are in agreement with observations; envelopes are bigger and bluer at maximum light.

Infrared emission from the dust shell around IRC+10216 is analysed in detail, employing a self-consistent model for radiatively driven winds around late-type stars that couples the equations of motion and radiative transfer in the dust. The resulting model provides agreement with the wealth of available data, including the spectral energy distribution in the range 0.5--1000 \mic, and visibility and array observations. Previous conclusions about two dust shells, derived from modelling the data with a few single-temperature components of different radii, are not supported by our results. The extended, continuous temperature and density distributions derived from our model obviate the need for such discrete shells. The IR properties vary with the stellar phase, reflecting changes in both the dust condensation radius $r_1$ and the overall optical depth $\tau$ -- as the luminosity increases from minimum to maximum, $r_1$ increases while $\tau$ decreases. We find that the angular size of the dust condensation zone varies from 0.3 arcsec at minimum light to 0.5 arcsec at maximum. The shortage of flux at short wavelengths encountered in previous studies is resolved by employing a grain size distribution that includes grains larger than \about\ 0.1 \mic, required also for the visibility fits. This distribution is in agreement with the one recently proposed by Jura in a study that probed the outer regions of the envelope. Since our constraints on the size distribution mostly reflect the envelope's inner regions, the agreement of these independent studies is evidence against significant changes in grain sizes through effects like sputtering or grain growth after the initial formation at the dust condensation zone.

We present the analytical solution of a maser shaped like a spherical shell. We determine the general condition on the size of the central cavity at which a sphere becomes a shell maser, and derive the intensity, beaming angle and observed size of both unsaturated and saturated shells.

General solutions of the maser polarization problem are presented for arbitrary absorption coefficients. The results are used to calculate polarization for masers permeated by magnetic fields with arbitrary values of \xB, the ratio of Zeeman splitting to Doppler linewidth, and for anisotropic pumping. The $\xb \to 0$ limit of the magnetic solution reproduces the linear polarization derived in previous studies, which were always conducted at this unphysical limit. While terms of higher order in \xb\ have a negligible effect on the magnitude of $q$, they produce some major new results. In particular, the linear polarization is accompanied by circular polarization, proportional to \xb. Because \xb\ is proportional to the transition wavelength, the circular polarization of SiO masers should decrease with rotation quantum number, as observed. In the absence of theory for $\xb < 1$, previous estimates of magnetic fields from detected maser circular polarization had to rely on conjectures in this case and generally need to be revised downward. The fields in SiO masers are \about\ 2--10 G and were overestimated by a factor of 8. The OH maser regions around supergiants have fields of \about\ 0.1--0.5 mG, which were overestimated by factors of 10--100. The fields were properly estimated for OH/IR masers ($\la$ 0.1 mG) and \H2O masers in star-forming regions (\about\ 15--50 mG). Spurious solutions that required stability analysis for their removal in all previous studies are never reproduced here; in particular, there are no stationary physical solutions for propagation at $\sin^2\theta < \third$, where $\theta$ is the angle from the direction of the magnetic field, so such radiation is unpolarized. These spurious solutions can be identified as the \xb\ = 0 limits of non-physical solutions and they never arise at finite

The dynamical structure and infrared emission of winds around late-type stars are studied in a self-consistent model that couples the equations of motion and radiative transfer. Both the dynamics and IR spectrum of the solution are fully characterized by the flux averaged optical depth. Five types of dust grains are considered: astronomical silicate, SiC, crystalline olivine, graphite and amorphous carbon as well as mixtures. Analysis of infrared signatures provides constraints on the grain chemical composition and indications for the simultaneous existence of silicate and carbon grains. The region of parameter space where radiation pressure can support a given mass-loss rate is identified and a new method to derive mass loss rates from IR data is presented. IRAS LRS classes are associated with optical depth for various grain materials and the regions of color-color diagrams expected to be populated by late-type stars are identified. We show that cirrus emission can severely affect point source measurements to the extent that their listed IRAS long wavelength fluxes are unreliable. After accounting for cirrus contamination, essentially all IRAS point sources located in the relevant regions of the color-color diagrams can be explained as late-type stars.