Jouni Kainulainen

Probing the evolution of molecular cloud structure: From quiescence to birth

Context. Probability distribution of densities is a fundamental measure of molecular cloud structure, containing information on how the material arranges itself in molecular clouds. Aims. We derive the probability density functions (PDFs) of column density for a complete sample of prominent molecular cloud complexes closer than d < 200 pc. For comparison, additional complexes at d ≈ 250−700 pc are included in the study. Methods. We derive near-infrared dust extinction maps for 23 molecular cloud complexes, using the nicest colour excess mapping technique and data from the 2MASS archive. The extinction maps are then used to examine the column density PDFs in the clouds. Results. The column density PDFs of most molecular clouds are well-fitted by log-normal functions at low column densities (0. 5m ag< AV < 3− 5m ag, or−0.5 < lnAV /AV < 1). But at higher column densities prominent power-law-like wings are common. In particular, we identify a trend among the PDFs: active star-forming clouds always have prominent non-log-normal wings. In contrast, clouds without active star formation resemble log-normals over the whole observed column density range or show only low excess of higher column densities. This trend is also reflected in the cumulative forms of the PDFs, showing that the fraction of high column density material is significantly larger in star-forming clouds. These observations agree with an evolutionary trend where turbulent motions are the main cloud-shaping mechanism for quiescent clouds, but the density enhancements induced by them quickly become dominated by gravity (and other mechanisms), which is in turn strongly reflected by the shape of the column density PDFs. The dominant role of the turbulence is restricted to the very early stages of molecular cloud evolution, comparable to the onset of active star formation in the clouds.

Probing the evolution of molecular cloud structure - II. From chaos to confinement

We present an analysis of the large-scale molecular cloud structure and of the stability of clumpy structures in nearby molecular clouds. In our recent work, we identified a structural transition in molecular clouds by studying the probability distributions of their gas column densities. In this paper, we further examine the nature of this transition. The transition takes place at the visual extinction of A tail = 2−4 mag, or equivalently, at Σ tail ≈ 40−80 Mpc −2 . The clumps identified above this limit have wide ranges of masses and sizes, but a remarkably constant mean volume density of n ≈ 10 3 cm −3 .T his is 5−10 times higher than the density of the medium surrounding the clumps. By examining the stability of the clumps, we show that they are gravitationally unbound entities, and that the external pressure from the parental molecular cloud is a significant source of confining pressure for them. Then, the structural transition at A tail may be linked to a transition between this population and the surrounding medium. The star-formation rates in the clouds correlate strongly with the total mass in the clumps, i.e., with the mass above A tail , and drops abruptly below that threshold. These results imply that the formation of pressure-confined clumps introduces a prerequisite for star formation. Furthermore, they give a physically motivated explanation for the recently reported relation between the star-formation rates and the amount of dense material in molecular clouds. Likewise, they give rise to a natural threshold for star formation at A tail .

DENSE CORES IN THE PIPE NEBULA: AN IMPROVED CORE MASS FUNCTION

In this paper we derive an improved core mass function (CMF) for the Pipe Nebula from a detailed comparison between measurements of visual extinction and molecular-line emission. We have compiled a refined sample of 201 dense cores toward the Pipe Nebula using a 2-dimensional threshold identification algorithm informed by recent simulations of dense core populations. Measurements of radial velocities using complimentary C18O (1-0) observations enable us to cull out from this sample those 43 extinction peaks that are either not associated with dense gas or are not physically associated with the Pipe Nebula. Moreover, we use the derived C18O, central velocities to differentiate between single cores with internal structure and blends of two or more physically distinct cores, superposed along the same line-of-sight. We then are able to produce a more robust dense core sample for future follow-up studies and a more reliable CMF than was possible previously. We confirm earlier indications that the CMF for the Pipe Nebula departs from a single power-law like form with a break or knee at M ~ 2.7 +/- 1.3 Msun. Moreover, we also confirm that the CMF exhibits a similar shape to the stellar IMF, but is scaled to higher masses by a factor of ~4.5. We interpret this difference in scaling to be a measure of the star formation efficiency (22 +/- 8%). This supports earlier suggestions that the stellar IMF may originate more or less directly from the CMF.

The Earliest Phases of Star formation (EPoS) observed with Herschel : the dust temperature and density distributions of B68

Context. Isolated starless cores within molecular clouds can be used as a testbed to investigate the conditions prior to the onset of fragmentation and gravitational proto-stellar collapse. Aims. We aim to determine the distribution of the dust temperature and the density of the starless core B68. Methods. In the framework of the Herschel guaranteed-time key programme “The Earliest Phases of Star formation” (EPoS), we have imaged B68 between 100 and 500 μm. Ancillary data at (sub)millimetre wavelengths, spectral line maps of the 12 CO (2–1), and 13 CO (2–1) transitions, as well as an NIR extinction map were added to the analysis. We employed a ray-tracing algorithm to derive the 2D mid-plane dust temperature and volume density distribution without suffering from the line-of-sight averaging effects of simple SED fitting procedures. Additional 3D radiative transfer calculations were employed to investigate the connection between the external irradiation and the peculiar crescent-shaped morphology found in the FIR maps. Results. For the first time, we spatially resolve the dust temperature and density distribution of B68, convolved to a beam size of 36. �� 4. We find a temperature gradient dropping from (16.7 +1.3 −1.0 ) K at the edge to (8.2

Dust-temperature of an isolated star-forming cloud: Herschel observations of the Bok globule CB244

We present Herschel observations of the isolated, low-mass star-forming Bok globule CB244. It contains two cold sources, a low-mass Class 0 protostar and a starless core, which is likely to be prestellar in nature, separated by 90 �� (∼18 000 AU). The Herschel data sample the peak of the Planck spectrum for these sources, and are therefore ideal for dust-temperature and column density modeling. With these data and a near-IR extinction map, the MIPS 70 μm mosaic, the SCUBA 850 μm map, and the IRAM 1.3 mm map, we model the dust-temperature and column density of CB 244 and present the first measured dust-temperature map of an entire starforming molecular cloud. We find that the column-averaged dust-temperature near the protostar is ∼17.7 K, while for the starless core it is ∼10.6 K, and that the effect of external heating causes the cloud dust-temperature to rise to ∼17 K where the hydrogen column density drops below 10 21 cm −2 . The total hydrogen mass of CB 244 (assuming a distance of 200 pc) is 15 ± 5 M� . The mass of the

Giant molecular filaments in the Milky Way

Throughout the Milky Way, molecular clouds typically appear filamentary, and mounting evidence indicates that this morphology plays an important role in star formation. What is not known is to what extent the dense filaments most closely associated with star formation are connected to the surrounding diffuse clouds up to arbitrarily large scales. How are these cradles of star formation linked to the Milky Way’s spiral structure? Using archival Galactic plane survey data, we have used multiple datasets in search of large-scale, velocity-coherent filaments in the Galactic plane. In this paper, we present our methods employed to identify coherent filamentary structures first in extinction and confirmed using Galactic Ring Survey data. We present a sample of seven giant molecular filaments (GMFs) that have lengths on the order of ~100 pc, total masses of 104–105 M⊙, and exhibit velocity coherence over their full length. The GMFs we study appear to be inter-arm clouds and may be the Milky Way analogs to spurs observed in nearby spiral galaxies. We find that between 2 and 12% of the total mass (above ~1020 cm-2) is “dense” (above 1022 cm-2), where filaments near spiral arms in the Galactic midplane tend to have higher dense gas mass fractions than those further from the arms.

Unfolding the laws of star formation: The density distribution of molecular clouds

Mapping Stardust A galaxys structure throughout time depends largely on its ability to convert the raw material of molecular clouds into stars. One of the most influential properties in determining star formation rates is the distribution of densities among individual molecular clouds, which can be described by a probability density function of volume densities. Kainulainen et al. (p. 183) devised a method to quantify these distributions from empirical dust extinction maps of nearby clouds. The threshold for star formation in these observationally based calculations was significantly lower than theoretical predictions. Models of star formation are better constrained through an empirical assessment of galactic density structures. The formation of stars shapes the structure and evolution of entire galaxies. The rate and efficiency of this process are affected substantially by the density structure of the individual molecular clouds in which stars form. The most fundamental measure of this structure is the probability density function of volume densities (ρ-PDF), which determines the star formation rates predicted with analytical models. This function has remained unconstrained by observations. We have developed an approach to quantify ρ-PDFs and establish their relation to star formation. The ρ-PDFs instigate a density threshold of star formation and allow us to quantify the star formation efficiency above it. The ρ-PDFs provide new constraints for star formation theories and correctly predict several key properties of the star-forming interstellar medium.

The star formation in the L1615/L1616 Cometary Cloud

The present work aims at performing a comprehensive census and characterization of the pre-main-sequence (PMS) population in the cometary cloud L1615/L1616, in order to assess the significance of the triggered star formation scenario and investigate the impact of massive stars on its star formation history and mass spectrum. Our study is based on UBVRCIC and JHKs photometry, as well as optical multiobject spectroscopy. We performed a physical parameterization of the young stellar population in L1615/L1616. We identified 25 new T Tauri stars mainly projected on the dense head of the cometary cloud, almost doubling the current number of known members. We studied the spatial distribution of the cloud members as a function of the age and H? emission. The star formation efficiency (SFE) in the cloud is ~7%-8%, as expected for molecular clouds in the vicinity of OB associations. The slope of the initial mass function (IMF), in the mass range -->0.1 M? ? M ? 5.5 M?, is consistent with that of other T and OB associations, providing further support of a universal IMF down to the hydrogen-burning limit, regardless of environmental conditions. The cometary appearance, as well as the high SFE, can be explained in terms of triggered star formation induced by the strong UV radiation from OB stars or supernova shock waves. The age spread and both the spatial and age distribution of the PMS objects provide strong evidence of sequential, multiple events and possibly still ongoing star formation activity in the cloud.

Star Formation in the Taurus Filament L 1495: From Dense Cores to Stars

We present a study of dense structures in the L1495 filament in the Taurus Molecular Cloud and examine its star-forming properties. In particular we construct a dust extinction map of the filament using deep near-infrared observations, exposing its small-scale structure in unprecedented detail. The filament shows highly fragmented substructures and a high mass-per-length value of Mline = 17Mpc −1 , reflecting star-forming potential in all parts of it. However, a part of the filament, namely B211, is remarkably devoid of young stellar objects. We argue that in this region the initial filament collapse and fragmentation is still taking place and star formation is yet to occur. In the star-forming part of the filament, we identify 39 cores with masses from 0.4...10Mand preferred separations in agreement with the local Jeans length. Most of these cores exceed the Bonnor-Ebert critical mass, and are therefore likely to collapse and form stars. The Dense Core Mass Function follows a power law with exponent = 1 .2 ± 0.2, a form commonly observed in star-forming regions. Subject headings: ISM: clouds — (ISM:) dust, extinction — ISM: individual (L1495) — ISM: structure — stars: formation

The Optical-Infrared Extinction Curve and its Variation in the Milky Way

The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity.