Monte Carlo radiative transfer in embedded prestellar cores

We implement a Monte Carlo radiative transfer method, that uses a large number of monochromatic luminosity packets to represent the radiation transported through a system. These packets are injected into the system and interact stochastically with it. We test our code against various benchmark calculations and determine the number of packets required to obtain accurate results under different circumstances. We then use this method to study cores that are directly exposed to the interstellar radiation field (non-embedded cores). Our code predicts temperature and intensity profiles inside these cores which are in good agreement with previous studies using different radiative transfer methods. We also explore a large number of models of cores that are embedded in the centre of a molecular cloud. We study cores with different density profiles embedded in molecular clouds with various optical extinctions and we calculate temperature profiles, SEDs and intensity profiles. Our study indicates that the temperature profiles in embedded cores are less steep than those in non-embedded cores. Deeply embedded cores (ambient cloud with visual extinction larger than 15-25) are almost isothermal at around 7-8 K. The temperature inside cores surrounded by an ambient cloud of even moderate thickness ( $A_{\rm V}\sim 5$) is less than 12 K, which is lower than previous studies have assumed. Thus, previous mass calculations of embedded cores (for example in the $\rho$ Ophiuchi protocluster), based on mm continuum observations, may underestimate core masses by up to a factor of 2. Our study shows that the best wavelength region to observe embedded cores is between 400 and 500 $\mu {\rm m}$, where the core is quite distinct from the background. We also predict that very sensitive observations (~ 1-3 MJy ${\rm\, sr^{-1}}$) at 170-200 $\mu {\rm m}$ can be used to estimate how deeply a core is embedded in its parent molecular cloud. The upcoming Herschel mission (ESA, 2007) will, in principle, be able to detect these features and test our models.

[1]  C. Curry,et al.  SFChem 2002: Chemistry as a Diagnostic of Star Formation , 2003 .

[2]  N. Santos,et al.  The Origins of Stars and Planets: The VLT View , 2002 .

[3]  J. Alves,et al.  The Origins of Stars and Planets: The VLT View: Proceedings of the ESO Workshop Held in Garching, Germany, 24–27 April 2001 , 2002 .

[4]  D. Johnstone,et al.  Large Area Mapping at 850 Microns. III. Analysis of the Clump Distribution in the Orion B Molecular Cloud , 2001 .

[5]  P. Andre',et al.  The initial conditions of isolated star formation – V. ISOPHOT imaging and the temperature and energy balance of pre stellar cores , 2001, astro-ph/0109173.

[6]  C. Walmsley,et al.  The dust temperature distribution in prestellar cores , 2001, astro-ph/0107055.

[7]  P. Andre',et al.  A SCUBA survey of the NGC 2068/2071 protoclusters , 2001 .

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

[9]  J. Alves,et al.  Internal structure of a cold dark molecular cloud inferred from the extinction of background starlight , 2001, Nature.

[10]  D. Johnstone,et al.  Large-Area Mapping at 850 Microns. II. Analysis of the Clump Distribution in the ρ Ophiuchi Molecular Cloud , 2000 .

[11]  L. Mundy,et al.  Tracing the Mass during Low-Mass Star Formation. II. Modeling the Submillimeter Emission from Preprotostellar Cores , 2000, astro-ph/0006183.

[12]  F. Motte,et al.  The initial conditions of isolated star formation — III. Millimetre continuum mapping of pre-stellar cores , 1999 .

[13]  A. Whitworth,et al.  Gas cooling by dust during dynamical fragmentation , 1998, astro-ph/9804338.

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

[15]  P. Andre',et al.  A submillimetre continuum survey of pre-protostellar cores , 1994 .

[16]  R. Cutri,et al.  The First Symposium on the Infrared Cirrus and Diffuse Interstellar Clouds : conference held 7-9 April 1993, Tucson, Arizona , 1994 .

[17]  B. Shustov Protostars and Planets II , 1987 .

[18]  C. Beichman,et al.  Candidate solar-type protostars in nearby molecular cloud cores , 1986 .

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

[20]  R. White,et al.  Bipolar reflection nebulae : Monte Carlo simulations. , 1984 .

[21]  P. Myers,et al.  Dense cores in dark clouds. II. NH3 observations and star formation. , 1983 .

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

[23]  W. B. Bonnor,et al.  Jeans' Formula for Gravitational Instability , 1957 .

[24]  W. Bonnar,et al.  Boyle's Law and gravitational instability , 1956 .

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