The dusty MOCASSIN: fully self-consistent 3D photoionization and dust radiative transfer models

We present the first 3D Monte Carlo (MC) photoionization code to include a fully self-consistent treatment of dust radiative transfer (RT) within a photoionized region. This is the latest development (version 2.0) of the gas-only photoionization code MOCASSIN and employs a stochastic approach to the transport of radiation, allowing both the primary and secondary components of the radiation field to be treated self-consistently, whilst accounting for the scattering of radiation by dust grains mixed with the gas, as well as the absorption and emission of radiation by both the gas and the dust components. An escape probability method is implemented for the transfer of resonance lines that may be absorbed by the grains, thus contributing to their energy balance. The energetics of the co-existing dust and gas components must also take into account the effects of dust‐gas collisions and photoelectric emission from the dust grains, which are dependent on the grain charge. These are included in our code using the average grain potential approximation scheme. A set of rigorous benchmark tests have been carried out for dust-only spherically symmetric geometries and 2D disc configurations. The results of MOCASSIN are found to be in agreement with those obtained by well-established dust-only RT codes that employ various approaches to the solution of the RT problem. A model of the dust and of the photoionized gas components of the planetary nebula NGC 3918 is also presented as a means of testing the correct functioning of the RT procedures in a case where both gas and dust opacities are present. The two components are coupled via the heating of dust grains by the absorption of both UV continuum photons and resonance line photons emitted by the gas. The MOCASSIN results show agreement with those of a 1D dust and gas model of this nebula published previously, showing the reliability of the new code, which can be applied to a variety of astrophysical environments. Ke yw ords: radiative transfer ‐ dust, extinction ‐ H II regions ‐ planetary nebulae: general.

[1]  J. Baldwin,et al.  Physical conditions in the Orion Nebula and an assessment of its helium abundance , 1991 .

[2]  J. Krolik,et al.  The formation of emission lines in quasars and Seyfert nuclei , 1981 .

[3]  D. Harper,et al.  Observations of cool dust in planetary nebulae. , 1980 .

[4]  J. Walsh,et al.  Planetary Nebulae Beyond the Milky Way , 2006 .

[5]  G. Plass,et al.  Electromagnetic scattering from absorbing spheres. , 1967, Applied optics.

[6]  MOCASSIN: a fully three-dimensional Monte Carlo photoionization code , 2002, astro-ph/0209378.

[7]  D. G. Hummer,et al.  Energy loss by resonance line photons in an absorbing medium , 1980 .

[8]  Antonella Natta Michael R. Meyer Steven V.W. Beckwith.  Silicate Emission in T Tauri Stars: Evidence for Disk Atmospheres? , 1999, astro-ph/9911490.

[9]  Ž. Ivezić,et al.  Erratum: Self-similarity and scaling behaviour of infrared emission from radiatively heated dust — I. Theory , 1997 .

[10]  Benchmark problems for dust radiative transfer , 1997 .

[11]  H. M. Lee,et al.  Optical properties of interstellar graphite and silicate grains , 1984 .

[12]  L. Lucy,et al.  Multiline Transfer and the Dynamics of Stellar Winds , 1985 .

[13]  P. B. Fellgett,et al.  Vistas in astronomy , 1957 .

[14]  L. C. Henyey,et al.  Diffuse radiation in the Galaxy , 1940 .

[15]  R. Corradi,et al.  Jets, Knots, and Tails in Planetary Nebulae: NGC 3918, K1-2, and Wray 17-1 , 1999, astro-ph/9912466.

[16]  B. Draine,et al.  Collisional charging of interstellar grains , 1987 .

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

[18]  Thermal infrared emission by dust in the planetary nebula NGC 3918: a model analysis of IRAS observations , 1988 .

[19]  B. Draine,et al.  Temperature fluctuations in interstellar grains. I. Computational method and sublimation of small grains , 1989 .

[20]  M. Hanner.  Infrared Observations of Comets Halley and Wilson and Properties of the Grains , 1988 .

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

[22]  France.,et al.  Three-dimensional photoionization modelling of the planetary nebula NGC 3918 , 2002, astro-ph/0209417.

[23]  J. Weingartner,et al.  Photoelectric Emission from Interstellar Dust: Grain Charging and Gas Heating , 1999, astro-ph/9907251.

[24]  G. Ferland,et al.  CLOUDY 90: Numerical Simulation of Plasmas and Their Spectra , 1998 .

[25]  M. Barlow,et al.  INFRARED PHOTOMETRY OF SOUTHERN PLANETARY-NEBULAE AND EMISSION-LINE OBJECTS , 1980 .

[26]  G. Ferland,et al.  Grain size distributions and photoelectric heating in ionized media , 2004, astro-ph/0402381.

[27]  M. Barlow,et al.  The dual dust chemistries of planetary nebulae with [WCL] central stars , 2002 .

[28]  Ari Laor,et al.  Spectroscopic constraints on the properties of dust in active galactic nuclei , 1993 .

[29]  D. Axon,et al.  Polarization Profiles of Scattered Emission Lines. III. Effects of Multiple Scattering and Non-Rayleigh Phase Functions , 1995 .

[30]  K. Wood,et al.  Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer , 2001, astro-ph/0103249.

[31]  J. Harrington,et al.  A Grain-heated, Dusty Planetary Nebula in M22 , 1991 .