M. Juvela

The Turbulent Shock Origin of Proto-Stellar Cores

The fragmentation of molecular clouds (MC) into proto-stellar cores is a central aspect of the process of star formation. Because of the turbulent nature of supersonic motions in MCs, it has been suggested that dense structures such as filaments and clumps are formed by shocks in a turbulent flow. In this work we present strong evidence in favor of the turbulent origin of the fragmentation of MCs. The most generic result of turbulent fragmentation is that dense postshock gas traces a gas component with a smaller velocity dispersion than lower density gas, since shocks correspond to regions of converging flows, where the kinetic energy of the turbulent motion is dissipated. Using synthetic maps of spectra of molecular transitions, computed from the results of numerical simulations of supersonic turbulence, we show that the dependence of velocity dispersion on gas density generates an observable relation between the rms velocity centroid and the integrated intensity (column density), σ(V0)-I, which is indeed found in the observational data. The comparison between the theoretical model (maps of synthetic 13CO spectra) with 13CO maps from the Perseus, Rosette, and Taurus MC complexes shows excellent agreement in the σ(V0)-I relation. The σ(V0)-I relation of different observational maps with the same total rms velocity are remarkably similar, which is a strong indication of their origin from a very general property of the fluid equations, such as the turbulent fragmentation process.

Dust temperature tracing the ISRF intensity in the Galaxy

New observations withHerschel allow accurate measurement of the equilibrium temperature of large dust grains heated by the interstellar radiation field (ISRF), which is critical in deriving dust column density and masses. We present temperature maps derived from the Herschel SPIRE and PACS data in two fields along the Galactic plane, obtained as part of the Hi-GAL survey during the Herschel science demonstration phase (SDP). We analyze the distribution of the dust temperature spatially, as well as along the two lines-of-sight (LOS) through the Galaxy. The zero-level offsets in the Herschel maps were established by comparison with the IRAS and Planck data at comparable wavelengths. We derive maps of the dust temperature and optical depth by adjusting a detailed model for dust emission at each pixel. The dust temperature maps show variations in the ISRF intensity and reveal the intricate mixture of the warm dust heated by massive stars and the cold filamentary structures of embedded molecular clouds. The dust optical depth at 250 μm is well correlated with the gas column density, but with a significantly higher dust emissivity than in the solar neighborhood. We correlate the optical depth with 3-D cubes of the dust extinction to investigate variations in the ISRF strength and dust abundance along the line of sight through the spiral structure of the Galaxy. We show that the warmest dust along the LOS is located in the spiral arms of the Galaxy, and we quantify their respective IR contribution.

Luminous Infrared Galaxies with the Submillimeter Array. I. Survey Overview and the Central Gas to Dust Ratio

We present new data obtained with the Submillimeter Array for a sample of 14 nearby luminous and ultraluminous infrared galaxies. The galaxies were selected to have distances D_L 11.4. The galaxies were observed with spatial resolutions of order 1 kpc in the CO J = 3–2, CO J = 2–1,^(13)CO J = 2–1, and HCO+ J = 4–3 lines as well as the continuum at 880 μm and 1.3 mm. We have combined our CO and continuum data to measure an average gas-to-dust mass ratio of 120 ± 28 (rms deviation 109) in the central regions of these galaxies, very similar to the value of 150 determined for the Milky Way. This similarity is interesting given the more intense heating from the starburst and possibly accretion activity in the luminous infrared galaxies compared to the Milky Way. We find that the peak H_2 surface density correlates with the far-infrared luminosity, which suggests that galaxies with higher gas surface densities inside the central kiloparsec have a higher star formation rate. The lack of a significant correlation between total H_2 mass and far-infrared luminosity in our sample suggests that the increased star formation rate is due to the increased availability of molecular gas as fuel for star formation in the central regions. In contrast to previous analyses by other authors, we do not find a significant correlation between central gas surface density and the star formation efficiency, as traced by the ratio of far-infrared luminosity to nuclear gas mass. Our data show that it is the star formation rate, not the star formation efficiency, that increases with increasing central gas surface density in these galaxies.

Numerical methods for non-LTE line radiative transfer: Performance and convergence characteristics

Comparison is made between a number of independent computer programs for radiative transfer in molecular rota- tional lines. The test models are spherically symmetric circumstellar envelopes with a given density and temperature profile. The first two test models have a simple power law density distribution, constant temperature and a fictive 2-level molecule, while the other two test models consist of an inside-out collapsing envelope observed in rotational transitions of HCO + .F or the 2-level molecule test problems all codes agree well to within 0.2%, comparable to the accuracy of the individual codes, for low optical depth and up to 2% for high optical depths (τ = 4800). The problem of the collapsing cloud in HCO + has a larger spread in results, ranging up to 12% for the J = 4 population. The spread is largest at the radius where the transition from collisional to radiative excitation occurs. The resulting line profiles for the HCO + J = 4-3 transition agree to within 10%, i.e., within the calibration accuracy of most current telescopes. The comparison project and the results described in this paper provide a benchmark for future code development, and give an indication of the typical accuracy of present day calculations of molecular line transfer.

Galactic cold cores - III. General cloud properties

Context. In the project galactic cold cores we are carrying out Herschel photometric observations of cold regions of the interstellar clouds as previously identified with the Planck satellite. The aim of the project is to derive the physical properties of the population of cold clumps and to study its connection to ongoing and future star formation. Aims. We examine the cloud structure around the Planck detections in 71 fields observed with the Herschel SPIRE instrument by the summer of 2011. We wish to determine the general physical characteristics of the fields and to examine the morphology of the clouds where the cold high column density clumps are found. Methods. Using the Herschel SPIRE data, we derive colour temperature and column density maps of the fields. Together with ancillary data, we examine the infrared spectral energy distributions of the main clumps. The clouds are categorised according to their large scale morphology. With the help of recently released WISE satellite data, we look for signs of enhanced mid-infrared scattering (”coreshine”), an indication of growth of the dust grains, and have a first look at the star formation activity associated with the cold clumps. Results. The mapped clouds have distances ranging from similar to 100 pc to several kiloparsecs and cover a range of sizes and masses from cores of less than 10 M-circle dot to clouds with masses in excess of 10 000 M-circle dot. Most fields contain some filamentary structures and in about half of the cases a filament or a few filaments dominate the morphology. In one case out of ten, the clouds show a cometary shape or have sharp boundaries indicative of compression by an external force. The width of the filaments is typically similar to 0.2-0.3 pc. However, there is significant variation from 0.1 pc to 1 pc and the estimates are sensitive to the methods used and the very definition of a filament. Enhanced mid-infrared scattering, coreshine, was detected in four clouds with six additional tentative detections. The cloud LDN183 is included in our sample and remains the best example of this phenomenon. About half of the fields are associated with active star formation as indicated by the presence of mid-infrared point sources. The mid-infrared sources often coincide with structures whose sub-millimetre spectra are still dominated by the cold dust.

Theoretical Models of Polarized Dust Emission from Protostellar Cores

We model the polarized thermal dust emission from protostellar cores that are assembled by supersonic turbulent flows in molecular clouds. Self-gravitating cores are selected from a three-dimensional simulation of supersonic and super-Alfvenic magnetohydrodynamic (MHD) turbulence. The polarization 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 AV,0 = 3 mag, consistent with theoretical expectations for grain alignment mechanisms. Instead of using a specific set of grain properties, we adopt a maximum degree of polarization Pmax = 15%. Results are therefore sensitive mainly to the topology of the magnetic field (model A) and to the gas distribution that determines the distribution of AV (model B). Furthermore, the radiative transfer in the MHD model is 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 Pmax) are typical, despite the super-Alfvenic nature of the turbulence. (2) A steep decrease of P with increasing 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 AV,0 (model B). (3) The same behavior is hard to reproduce if grains are aligned independently of AV (model A). (4) The Chandrasekhar-Fermi formula, corrected by a factor f ≈ 0.4, provides an approximate estimate of the average magnetic field 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 field inside the cores at visual extinction larger than AV,0 ≈ 3 mag. The use of such maps to constrain models of 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.

Efficient Monte Carlo methods for continuum radiative transfer

We discuss the efficiency of Monte Carlo methods in solving continuum radiative transfer problems. The sampling of the radiation field and convergence of dust temperature calculations in the case of optically thick clouds are both studied. For spherically symmetric clouds we find that the computational cost of Monte Carlo simulations can be reduced, in some cases by orders of magnitude, with simple importance weighting schemes. This is particularly true for models consisting of cells of different sizes, for which the run times would otherwise be determined by the size of the smallest cell. The use of importance weighting is extended to scattered photons, which is found to be useful in calculations of scattered flux and could be important for three-dimensional models when observed intensity is needed only for one general direction of observations. Convergence of dust temperature calculations is studied for models with optical depths

Supersonic Turbulence in the Perseus Molecular Cloud

\tau_{\rm V}=10{-}10^4

The Power Spectrum of Turbulence in NGC 1333: Outflows or Large-Scale Driving?

. We examine acceleration methods where radiative interactions inside a cell or between neighbouring cells are treated explicitly. In optically thick clouds with strong self-coupling between dust temperatures, the run times can be reduced by more than one order of magnitude. Use of a reference field was also examined. It eliminates the need for repeating simulation of constant sources (e.g., background radiation) after the first iteration and it was found to significantly reduce sampling errors. We finally discuss the applicability of our methods to three-dimensional models.

Synthetic Molecular Clouds from Supersonic MHD and Non-LTE Radiative Transfer Calculations

We compare the statistical properties of J=1-0 13CO spectra observed in the Perseus Molecular Cloud with synthetic J=1-0 13CO spectra, computed solving the non-LTE radiative transfer problem for a model cloud obtained as solutions of the three dimensional magneto-hydrodynamic (MHD) equations. The model cloud is a randomly forced super-Alfvenic and highly super-sonic turbulent isothermal flow. The purpose of the present work is to test if idealized turbulent flows, without self-gravity, stellar radiation, stellar outflows, or any other effect of star formation, are inconsistent or not with statistical properties of star forming molecular clouds. We present several statistical results that demonstrate remarkable similarity between real data and the synthetic cloud. Statistical properties of molecular clouds like Perseus are appropriately described by random super-sonic and super-Alfvenic MHD flows. Although the description of gravity and stellar radiation are essential to understand the formation of single protostars and the effects of star formation in the cloud dynamics, the overall description of the cloud and of the initial conditions for star formation can apparently be provided on intermediate scales without accounting for gravity, stellar radiation, and a detailed modeling of stellar outflows. We also show that the relation between equivalent line width and integrated antenna temperature indicates the presence of a relatively strong magnetic field in the core B1, in agreement with Zeeman splitting measurements.