Magnetic fields are believed to play a crucial role in the process of star formation. We compare high-angular resolution observations of the submillimeter polarized emission of NGC 1333 IRAS 4A, tracing the magnetic field around a low-mass protostar, with models of the collapse of magnetized molecular cloud cores. Assuming a uniform dust alignment efficiency, we computed the Stokes parameters and synthetic polarization maps from the model density and magnetic field distribution by integrations along the line-of-sight and convolution with the interferometric response. The synthetic maps are in good agreement with the data. The best-fitting models were obtained for a protostellar mass of 0.8 solar masses, of age 9e4 yr, formed in a cloud with an initial mass-to-flux ratio ~2 times the critical value. The magnetic field morphology in NGC 1333 IRAS 4A is consistent with the standard theoretical scenario for the formation of solar-type stars, where well-ordered, large-scale, rather than turbulent, magnetic fields control the evolution and collapse of the molecular cloud cores from which stars form.

We compute polarization maps for molecular cloud cores modeled as magnetized singular isothermal toroids, under the assumption that the emitting dust grains are aspherical and aligned with the large-scale magnetic field. We show that, depending on the inclination of the toroid with the line-of-sight, the bending of the magnetic field lines resulting from the need to counteract the inward pull of gravity naturally produces a depolarization effect toward the centre of the map. We compute the decrease of polarization degree with increasing intensity for different viewing angles and frequencies, and we show that an outward increasing temperature gradient, as expected in starless cores heated by the external radiation field, enhances the decrease of polarization. We compare our results with recent observations, and we conclude that this geometrical effect, together with other mechanisms of depolarization, may significantly contribute to the decrease of polarization degrees with intensity observed in the majority of molecular cloud cores. Finally, we consider the dependence of the polarization degree on the dust temperature gradient predicted for externally heated clouds, and we briefly comment on the limits of the Chandrasekhar-Fermi formula to estimate the magnetic field strength in molecular cloud cores.

We present the results of a three-dimensional Monte Carlo radiative transfer code for starless molecular cloud cores heated by an external isotropic or non-isotropic interstellar radiation field. The code computes the dust temperature distribution inside model clouds with specified but arbitrary density profiles. In particular we examine in detail spherical (Bonnor-Ebert) clouds, axisymmetric and non-axisymmetric toroids, and clouds heated by an external stellar source in addition to the general interstellar field. For these configurations, the code also computes maps of the emergent intensity at different wavelengths and arbitrary viewing angle, that can be compared directly with continuum maps of prestellar cores. In the approximation where the dust temperature is independent of interactions with the gas and where the gas is heated both by collisions with dust grains and ionization by cosmic rays, the temperature distribution of the gas is also calculated. For cloud models with parameters typical of dense cores, the results show that the dust temperature decreases monotonically from a maximum value near the cloud's edge (14-15 K) to a minimum value at the cloud's center (6-7 K). Conversely, the gas temperature varies in a similar range, but, due to efficient dust-gas coupling in the inner regions and inefficient cosmic-ray heating in the outer regions, the gradient is non-monotonic and the gas temperature reaches a maximum value at intermediate radii. The emission computed for these models (at 350 micron and 1.3 mm) shows that deviations from spherical symmetry in the density and/or temperature distributions are generally reduced in the simulated intensity maps (even without beam convolution), especially at the longer wavelengths.

We have observed the dark cloud Barnard 68 with SCUBA at 850 um and with SIMBA at 1.2 mm. The submillimetre and millimetre dust emission correlate well with the extinction map of Alves, Lada and Lada (2001).The A_V/850um correlation is clearly not linear and suggests lower temperatures for the dust in the inner core of the cloud. Assuming a model for the temperature gradient, we derive the cloud-averaged dust emissivities (normalised to the V-Band extinction efficiency) at 850 um and 1.2 mm. We find k_850um/k_V = 4.0 +/- 1.0 x 10^-5 and k_1.2mm/k_V = 9.0 +/- 3.0 x 10^-6. These values are compared with other determinations in this wavelength regime and with expectations for models of diffuse dust and grain growth in dense clouds.

We have considered the thermal equilibrium in pre-protostellar cores in the approximation where the dust temperature is independent of interactions with the gas and where the gas is heated both by collisions with dust grains and ionization by cosmic rays. We have then used these results to study the stability of cores in the limit where thermal pressure dominates over magnetic field and turbulence. We find that for cores with characteristics similar to those observed, the gas and dust temperatures are coupled in the core interior. As a consequence, the gas temperature like the dust temperature decreases towards the center of these objects. The density structure computed taking into account such deviations from isothermality are not greatly different from that expected for an isothermal Bonnor-Ebert sphere. It is impossible in the framework of these models to have a stable equilibrium core with mass above about 5 solar masses and column density compatible with observed values. We conclude from this that observed high mass cores are either supported by magnetic field or turbulence or are already in a state of collapse. Lower mass cores on the other hand have stable states and we conclude that the much studied object B68 may be in a state of stable equilibrium if the internal gas temperature is computed in self-consistent fashion. Finally we note that in molecular clouds such as Ophiuchus and Orion with high radiation fields and pressures, gas and dust temperatures are expected to be well coupled and hence one expects temperatures to be relatively high as compared to low pressure clouds like Taurus.