Three-dimensional radiative transfer models of clumpy tori in Seyfert galaxies

Context. Tori of Active Galactic Nuclei (AGN) are made up of a mixture of hot and cold gas, as well as dust. In order to protect the dust grains from destruction by the surrounding hot gas as well as by the energetic (UV/optical) radiation from the accretion disk, the dust is often assumed to be distributed in clouds. Aims. A new three-dimensional model of AGN dust tori is extensively investigated. The torus is modelled as a wedge-shaped disk within which dusty clouds are randomly distributed throughout the volume, by taking the dust density distribution of the corresponding continuous model into account. We especially concentrate on the differences between clumpy and continuous models in terms of the temperature distributions, the surface brightness distributions and interferometric visibilities, as well as spectral energy distributions. Methods. Radiative transfer calculations with the help of the three-dimensional Monte Carlo radiative transfer code MC3D are used in order to simulate spectral energy distributions as well as surface brightness distributions at various wavelengths. In a second step, interferometric visibilities for various inclination as well as position angles and baselines are calculated, which can be used to directly compare our models to interferometric observations with the MIDI instrument. Results. We find that the radial temperature distributions of clumpy models possess significantly enhanced scatter compared to the continuous cases. Even at large distances, clouds can be heated directly by the central accretion disk. The existence of the silicate 10 µm-feature in absorption or in emission depends sensitively on the distribution, the size and optical depth of clouds in the innermost part of the dust distribution. With this explanation, failure and success of previous modelling efforts of clumpy tori can be understood. The main reason for this outcome are shadowing effects of clouds within the central region. We underline this result with the help of several parameter variations. After adapting the parameters of our clumpy standard model to the circumstances of the Seyfert 2 Circinus galaxy, it can qualitatively explain recent mid-infrared interferometric observations performed with MIDI, as well as high resolution spectral data.

[1]  K. Nordsieck,et al.  The Size distribution of interstellar grains , 1977 .

[2]  J. Miller,et al.  Evidence for a highly polarized continuum in the nucleus of NGC 1068 , 1983 .

[3]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[4]  Julian H. Krolik,et al.  Molecular tori in Seyfert galaxies - Feeding the monster and hiding it , 1988 .

[5]  William H. Press,et al.  Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing , 1989 .

[6]  Julian H. Krolik,et al.  Infrared spectra of obscuring dust tori around active galactic nuclei. I - Calculational method and basic trends , 1992 .

[7]  E. Pier,et al.  Radiation-pressure-supported obscuring tori around active galactic nuclei , 1992 .

[8]  William H. Press,et al.  Numerical recipes in FORTRAN (2nd ed.): the art of scientific computing , 1992 .

[9]  Robert Antonucci,et al.  Unified models for active galactic nuclei and quasars , 1993 .

[10]  A. Königl,et al.  DISK-DRIVEN HYDROMAGNETIC WINDS AS A KEY INGREDIENT OF ACTIVE GALACTIC NUCLEI UNIFICATION SCHEMES , 1994 .

[11]  L. Danese,et al.  Thick tori around active galactic nuclei: a comparison of model predictions with observations of the infrared continuum and silicate features , 1994 .

[12]  P. Padovani,et al.  UNIFIED SCHEMES FOR RADIO-LOUD ACTIVE GALACTIC NUCLEI , 1995, astro-ph/9506063.

[13]  George H. Rieke,et al.  Low-luminosity and obscured seyfert nuclei in nearby galaxies , 1995 .

[14]  Sebastian Wolf,et al.  Accelerated self-consistent radiative transfer based on the Monte Carlo method , 2000 .

[15]  P. Drummond,et al.  Computer Physics Communications , 2001 .

[16]  A. Marconi,et al.  Dust in active nuclei. I. evidence for anomalous properties , 2000, astro-ph/0010009.

[17]  M. Elvis,et al.  Ubiquitous Variability of X-Ray-absorbing Column Densities in Seyfert 2 Galaxies , 2001, astro-ph/0107510.

[18]  E. Oliva,et al.  Dust in active nuclei - II. Powder or gravel? , 2000, astro-ph/0010066.

[19]  Z. Ivezic,et al.  Dust Emission from Active Galactic Nuclei , 2002 .

[20]  G. Granato,et al.  Revisiting the location and environment of the central engine in NGC 1068 , 2003, astro-ph/0304293.

[21]  R. Norris,et al.  A Warped Accretion Disk and Wide-Angle Outflow in the Inner Parsec of the Circinus Galaxy , 2003, astro-ph/0302533.

[22]  C. Dullemond,et al.  New radiative transfer models for obscuring tori in active galaxies , 2003, astro-ph/0303496.

[23]  Sebastiano Ligori,et al.  Ten-micron instrument MIDI: getting ready for observations on the VLTI , 2003, SPIE Astronomical Telescopes + Instrumentation.

[24]  S. Wolf Efficient Radiative Transfer in Dust Grain Mixtures , 2002, astro-ph/0209281.

[25]  K. Meisenheimer,et al.  Unveiling the Central Parsec Region of an Active Galactic Nucleus: The Circinus Nucleus in the Near-Infrared with the Very Large Telescope , 2004 .

[26]  The 2D continuum radiative transfer problem - Benchmark results for disk configurations , 2004, astro-ph/0402357.

[27]  W. Duschl,et al.  The dynamical state of a thick cloudy torus around an AGN , 2004 .

[28]  Supermassive black hole demography: the match between the local and accreted mass functions , 2004, astro-ph/0405585.

[29]  M. Schoeller,et al.  The central dusty torus in the active nucleus of NGC 1068 , 2004, Nature.

[30]  Bonn,et al.  From the Circumnuclear Disk in the Galactic Center to thick, obscuring tori of AGNs , 2003, astro-ph/0310329.

[31]  Klaus Meisenheimer,et al.  Towards a physical model of dust tori in Active Galactic Nuclei , 2005 .

[32]  Heidelberg,et al.  Towards a physical model of dust tori in Active Galactic Nuclei Radiative transfer calculations for a hydrostatic torus model , 2005, astro-ph/0504105.

[33]  B. Brandl,et al.  Mid-Infrared Spectra of Classical AGNs Observed with the Spitzer Space Telescope , 2005, astro-ph/0507423.

[34]  R. Siebenmorgen,et al.  Discovery of 10 μm silicate emission in quasars : Evidence of the agn unification scheme , 2005, astro-ph/0504263.

[35]  A. Tielens,et al.  The Detection of Silicate Emission from Quasars at 10 and 18 Microns , 2005, astro-ph/0504423.

[36]  C. Dullemond,et al.  Clumpy tori around active galactic nuclei , 2005, astro-ph/0501570.

[37]  D. Maoz,et al.  Accepted for publication in the Astrophysical Journal Letters Silicate emissions in active galaxies- From LINERs to QSOs , 2005 .

[38]  G. Weigelt,et al.  Radiative transfer modeling of three-dimensional clumpy AGN tori and its application to NGC 1068 , 2006, astro-ph/0602494.

[39]  Paul S. Smith,et al.  9.7 μm Silicate Features in Active Galactic Nuclei: New Insights into Unification Models , 2006, astro-ph/0608645.

[40]  Observatoire de la Cote d'Azur,et al.  Resolving the complex structure of the dust torus in the active nucleus of the Circinus galaxy , 2007, 0709.0209.

[41]  Julian H. Krolik,et al.  AGN obscuring tori supported by infrared radiation pressure , 2007 .