Theoretical Models of Polarized Dust Emission from Protostellar Cores

Abstract

We model the polarized thermal dust emission from protostellar cores that are assembled by super-sonic turbulent Ñows in molecular clouds. Self-gravitating cores are selected from a three-dimensional simulation of supersonic and magnetohydrodynamic (MHD) turbulence. The polarization super-Alfve nic is computed in two ways. In model A it is assumed that dust properties and grain alignment efficiency are uniform ; in model B it is assumed that grains are not aligned at visual extinction larger than A V,0 \ 3 mag, consistent with theoretical expectations for grain alignment mechanisms. Instead of using a spe-ciÐc set of grain properties, we adopt a maximum degree of polarization Results are there-P max \ 15%. fore sensitive mainly to the topology of the magnetic Ðeld (model A) and to the gas distribution that determines the distribution of (model B). Furthermore, the radiative transfer in the MHD model is A V solved with a non-LTE Monte Carlo method, to compute spectral maps of the J \ 1È0 transition of CS. The CS spectral maps are used to estimate the turbulent velocity, as in the observations. The main results of this work are the following : (1) Values of P between 1% and 10% (up to almost are P max) typical, despite the nature of the turbulence. (2) A steep decrease of P with increasing super-Alfve nic values of the submillimeter dust continuum intensity I is always found in self-gravitating cores selected from the MHD simulations if grains are not aligned above a certain value of visual extinction A V,0 (model B). (3) The same behavior is hard to reproduce if grains are aligned independently of (model A V A). (4) The Chandrasekhar-Fermi formula, corrected by a factor f B 0.4, provides an approximate esti-mate of the average magnetic Ðeld strength in the cores. Submillimeter dust continuum polarization maps of quiescent protostellar cores and Bok globules have recently been obtained. They always show a decrease in P with increasing value of I consistent with the predictions of our model B. We therefore conclude that submillimeter polarization maps of quiescent cores do not map the magnetic Ðeld inside the cores at visual extinction larger than mag. The use of such maps to constrain models of A V,0 B 3 protostellar core formation and evolution is questionable. This conclusion suggests that there is no inconsistency between the results from optical and near-IR polarized absorption of background stars and the observed polarization of submillimeter dust continuum from quiescent cores. In both cases, grains at large visual extinction appear to be virtually unaligned. Subject headings : ISM : clouds È ISM : kinematics and dynamics È polarization È radio continuum : ISM È submillimeter È turbulence On-line material : color Ðgures 1. INTRODUCTION The topology of the magnetic Ðeld inside and around protostellar cores can be predicted with models for their formation and evolution. Models of subcritical cores con-tracting under the e ect of gravity and ambipolar drift, for example, predict a rather uniform magnetic Ðeld roughly perpendicular to the core major axis. Such models can be tested if the magnetic Ðeld topology is constrained obser-vationally. Polarization maps of submillimeter thermal dust emis-sion have recently been obtained for a number of protostel-1 ppadoan=cfa.harvard.edu.