On the 10-micron silicate feature in Active Galactic Nuclei

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)Comment: 11 pages, 9 figures, 4 tables, accepted for publication in ApJ; minor revision, added reference in Section

Nonlocal radiative coupling in non monotonic stellar winds

There is strong observational evidence of shocks and clumping in radiation-driven stellar winds from hot, luminous stars. The resulting non monotonic velocity law allows for radiative coupling between distant locations, which is so far not accounted for in hydrodynamic wind simulations. In the present paper, we determine the Sobolev source function and radiative line force in the presence of radiative coupling in spherically symmetric flows, extending the geometry-free formalism of Rybicki and Hummer (1978) to the case of three-point coupling, which can result from, e.g., corotating interaction regions, wind shocks, or mass overloading. For a simple model of an overloaded wind, we find that, surprisingly, the flow decelerates at all radii above a certain height when nonlocal radiative coupling is accounted for. We discuss whether radiation-driven winds might in general not be able to re-accelerate after a non monotonicity has occurred in the velocity law.Comment: accepted by A&A, 8 pages, 4 figure

A Tale of Three Galaxies: Deciphering the Infrared Emission of the Spectroscopically Anomalous Galaxies IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948

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

Modeling the Infrared Reverberation Response of the Circumnuclear Dusty Torus in AGN: The Effects of Cloud Orientation and Anisotropic Illumination

The obscuring circumnuclear torus of dusty molecular gas is one of the major components of active galactic nuclei (AGN). The torus can be studied by analyzing the time response of its infrared (IR) dust emission to variations in the AGN continuum luminosity, a technique known as reverberation mapping. The IR response is the convolution of the AGN ultraviolet/optical light curve with a transfer function that contains information about the size, geometry, and structure of the torus. Here, we describe a new computer model that simulates the reverberation response of a clumpy torus. Given an input optical light curve, the code computes the emission of a 3D ensemble of dust clouds as a function of time at selected IR wavelengths, taking into account light travel delays. We present simulated dust emission responses at 3.6, 4.5, and 30 $\mu$m that explore the effects of various geometrical and structural properties, dust cloud orientation, and anisotropy of the illuminating radiation field. We also briefly explore the effects of cloud shadowing (clouds are shielded from the AGN continuum source). Example synthetic light curves have also been generated, using the observed optical light curve of the Seyfert 1 galaxy NGC 6418 as the input. The torus response is strongly wavelength-dependent, due to the gradient in cloud surface temperature within the torus, and because the cloud emission is strongly anisotropic at shorter wavelengths. Anisotropic illumination of the torus also significantly modifies the torus response, reducing the lag between the IR and optical variations.Comment: 17 pages, 14 figures, published in the Astrophysical Journal (2017 July 1

A Tale of Three Galaxies: A "Clumpy" View of the Spectroscopically Anomalous Galaxies IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948

We investigate the dust properties in three spectroscopically anomalous galaxies (IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948). Their Spitzer/IRS spectra are characterized by a steep ~5-8 micron emission continuum, strong emission bands from polycyclic aromatic hydrocarbon (PAH) molecules, and prominent 10 micron silicate emission. The steep ~5-8 micron continuum and strong PAH emission features suggest the presence of starbursts, while the silicate emission is indicative of significant heating from AGNs. The simultaneous detection of these two observational properties has rarely been reported on galactic scale. We employ the PAHFIT software to estimate their starlight contributions, and the CLUMPY model for the components contributed by the AGN tori. We find that the CLUMPY model is generally successful in explaining the overall dust infrared emission, although it appears to emit too flat at the ~5-8 micron continuum to be consistent with that observed in IRAS F10398+1455 and IRAS F21013-0739. The flat ~5-8 micron continuum calculated from the CLUMPY model could arise from the adopted specific silicate opacity of Ossenkopf et al. (1992) which exceeds that of the Draine & Lee (1984) "astronomical silicate" by a factor up to 2 in the ~5-8 micron wavelength range. Future models with a variety of dust species incorporated in the CLUMPY radiation transfer regime are needed for a thorough understanding of the dust properties of these spectroscopically anomalous galaxies.Comment: 23 pages, 4 figures, 2 tables; accepted for publication in Planetary and Space Science, special issue on "Cosmic Dust

On the Unification of Active Galactic Nuclei

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.Comment: ApJ Letters, to be published. This is the final, journal version; minor editing revisions from original on

New ATCA, ALMA and VISIR observations of the candidate LBV SK-67266 (S61): the nebular mass from modelling 3D density distributions

We present new observations of the nebula around the Magellanic candidate Luminous Blue Variable S61. These comprise high-resolution data acquired with the Australia Telescope Compact Array (ATCA), the Atacama Large Millimetre/Submillimetre Array (ALMA), and VISIR at the Very Large Telescope (VLT). The nebula was detected only in the radio, up to 17 GHz. The 17 GHz ATCA map, with 0.8 arcsec resolution, allowed a morphological comparison with the H$\alpha$ Hubble Space Telescope image. The radio nebula resembles a spherical shell, as in the optical. The spectral index map indicates that the radio emission is due to free-free transitions in the ionised, optically thin gas, but there are hints of inhomogeneities. We present our new public code RHOCUBE to model 3D density distributions, and determine via Bayesian inference the nebula's geometric parameters. We applied the code to model the electron density distribution in the S61 nebula. We found that different distributions fit the data, but all of them converge to the same ionised mass, ~0.1 $\rm M\odot$, which is an order of magnitude smaller than previous estimates. We show how the nebula models can be used to derive the mass-loss history with high-temporal resolution. The nebula was probably formed through stellar winds, rather than eruptions. From the ALMA and VISIR non-detections, plus the derived extinction map, we deduce that the infrared emission observed by space telescopes must arise from extended, diffuse dust within the ionised region.Comment: 17 pages, 9 figures. Authors list corrected. In press in MNRAS. RHOCUBE code available online ( https://github.com/rnikutta/rhocube

Subaru Spectroscopy and Spectral Modeling of Cygnus A

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