Monte Carlo radiative transfer in protoplanetary disks

Aims.We present a new continuum 3D radiative transfer code, MCFOST, based on a Monte-Carlo method. MCFOST can be used to calculate (i) monochromatic images in scattered light and/or thermal emission; (ii) polarisation maps; (iii) interferometric visibilities; (iv) spectral energy distributions; and (v) dust temperature distributions of protoplanetary disks. Methods: .Several improvements to the standard Monte Carlo method are implemented in MCFOST to increase efficiency and reduce convergence time, including wavelength distribution adjustments, mean intensity calculations, and an adaptive sampling of the radiation field. The reliability and efficiency of the code are tested against a previously-defined benchmark, using a 2D disk configuration. No significant difference (no more than 10% and usually much less) is found between the temperatures and SEDs calculated by MCFOST and by other codes included in the benchmark. Results: . A study of the lowest disk mass detectable by Spitzer, around young stars, is presented and the colours of "representative" parametric disks compared to recent IRAC and MIPS Spitzer colours of solar-like young stars located in nearby star-forming regions.

[1]  A. Ghez,et al.  A Multiwavelength Scattered Light Analysis of the Dust Grain Population in the GG Tauri Circumbinary Ring , 2004, astro-ph/0401560.

[2]  M. Tamura,et al.  Polarimetry of young stellar objects — I. Linear polarization of GSS 30 , 1996 .

[3]  G. Fazio,et al.  IRAC Observations of Taurus Pre-Main-Sequence Stars , 2005, astro-ph/0505323.

[4]  F. Allard,et al.  Evolutionary Models for Very Low-Mass Stars and Brown Dwarfs with Dusty Atmospheres , 2000 .

[5]  Radiative equilibrium in Monte Carlo radiative transfer using frequency distribution adjustment , 2005, astro-ph/0504037.

[6]  John E. Krist,et al.  An Edge-on Circumstellar Disk in the Young Binary System HK Tauri , 1998 .

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

[8]  F. Allard,et al.  Evolutionary models for low-mass stars and brown dwarfs: uncertainties and limits at very young ages , 2002 .

[9]  J. Mathis,et al.  Composite interstellar grains , 1989 .

[10]  Lynne Hillenbrand,et al.  On the Evolutionary Status of Class I Stars and Herbig-Haro Energy Sources in Taurus-Auriga , 2004 .

[11]  L. Hartmann,et al.  Pre-Main-Sequence Evolution in the Taurus-Auriga Molecular Cloud , 1995 .

[12]  T. Fusco,et al.  Adaptive optics imaging of the MBM 12 association: Seven binaries and an edge-on disk in a quadruple system , 2002, astro-ph/0209290.

[13]  Sao,et al.  Effects of Dust Growth and Settling in T Tauri Disks , 2005, astro-ph/0511564.

[14]  A. Ghez,et al.  NICMOS Images of the GG Tauri Circumbinary Disk , 2002, astro-ph/0204465.

[15]  Radiative transfer in protoplanetary disks , 2006, astro-ph/0604067.

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

[17]  John E. Krist,et al.  WFPC2 Images of a Face-on Disk Surrounding TW Hydrae , 2000 .

[18]  D. Padgett,et al.  Hubble Space Telescope WFPC2 Imaging of the Disk and Jet of HV Tauri C , 2003 .

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

[20]  E. al.,et al.  From molecular cores to planet-forming disks: An SIRTF legacy program , 2003, astro-ph/0305127.