In this chapter, we review some historical understanding and recent advances on the Initial Mass Function (IMF) and the Core Mass Function (CMF), both in terms of observations and theories. We focus mostly on star formation in clustered environment since this is suggested by observations to be the dominant mode of star formation. The statistical properties and the fragmentation behaviour of turbulent gas is discussed, and we also discuss the formation of binaries and small multiple systems.

Molecular clouds are a fundamental ingredient of galaxies: they are the channels that transform the diffuse gas into stars. The detailed process of how they do it is not completely understood. We review the current knowledge of molecular clouds and their substructure from scales $\sim~$1~kpc down to the filament and core scale. We first review the mechanisms of cloud formation from the warm diffuse interstellar medium down to the cold and dense molecular clouds, the process of molecule formation and the role of the thermal and gravitational instabilities. We also discuss the main physical mechanisms through which clouds gather their mass, and note that all of them may have a role at various stages of the process. In order to understand the dynamics of clouds we then give a critical review of the widely used virial theorem, and its relation to the measurable properties of molecular clouds. Since these properties are the tools we have for understanding the dynamical state of clouds, we critically analyse them. We finally discuss the ubiquitous filamentary structure of molecular clouds and its connection to prestellar cores and star formation.

Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work.

Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

We present a unified description of the scenario of Global Hierarchical Collapse and fragmentation (GHC) in molecular clouds (MCs), owing to the continuous decrease of the average Jeans mass in the contracting cloud. GHC constitutes a regime of collapses within collapses, in which small-scale collapses begin at later times, but occur on shorter timescales than large-scale ones. The difference in timescales allows for most of the clouds' mass to be dispersed by feedback from the first massive stars, maintaining the global star formation rate low. All scales accrete from their parent structures. The main features of GHC are: star-forming MCs are in an essentially pressureless regime, which produces filaments that accrete onto clumps and cores ("hubs"). The filaments constitute the collapse flow from cloud to hub scales and may approach a quasi-stationary state; the molecular and dense mass fractions of the clouds increase over time; the first (low-mass) stars appear several Myr after global contraction began; more massive stars appear after a few Myr in massive hubs resulting from the collapse of larger scales; the minimum fragment mass may extend well into the brown-dwarf regime; Bondi-Hoyle-Lyttleton accretion occurs at the protostellar and core scales, accounting for a near-Salpeter IMF; the extreme anisotropy of the filamentary network explains the difficulty in detecting large-scale infall signatures; the balance between inertial and gravitationally-driven motions in clumps evolves during the contraction; prestellar cores adopt Bonnor-Ebert-like profiles, but are contracting ever since early times when they may appear to be unbound and to require pressure confinement; stellar clusters develop radial age and mass segregation gradients. Finally, we discuss the incompatibility between supersonic turbulence and the observed scalings in the molecular hierarchy.

Different studies have reported a power-law mass-size relation $M \propto R^q$ for ensembles of molecular clouds. In the case of nearby clouds, the index of the power-law $q$ is close to 2. However, for clouds spread all over the Galaxy, indexes larger than 2 are reported. We show that indexes larger than 2 could be the result of line-of-sight superposition of emission that does not belong to the cloud itself. We found that a random factor of gas contamination, between 0.001\%\ and 10\%\ of the line-of-sight, allows to reproduce the mass-size relation with $q \sim 2.2-2.3$ observed in Galactic CO surveys. Furthermore, for dense cores within a single cloud, or molecular clouds within a single galaxy, we argue that, even in these cases, there is observational and theoretical evidence that some degree of superposition may be occurring. However, additional effects may be present in each case, and are briefly discussed. We also argue that defining the fractal dimension of clouds via the mass-size relation is not adequate, since the mass is not {necessarily} a proxy to the area, and the size reported in $M-R$ relations is typically obtained from the square root of the area, rather than from an estimation of the size independent from the area. Finally, we argue that the statistical analysis of finding clouds satisfying the Larson's relations does not mean that each individual cloud is in virial equilibrium.

Recent analyses of Gaia data have provided direct evidence that most young stellar clusters are in a state of expansion, with velocities of the order of ~0.5 km/s. Traditionally, expanding young clusters have been pictured as entities that became unbound due to the lack of gravitational binding once the gas from the parental cloud that formed the cluster has been expelled by the stellar radiation of the massive stars in the cluster. In the present contribution, we used radiation-magnetohydrodynamic numerical simulations of molecular cloud formation and evolution to understand how stellar clusters form and disperse. We found that the ionising feedback from the newborn massive stars expels the gas from the collapse centre, flipping-up the gravitational potential as a consequence of the mass removal from the inside-out. Since neither the parental clouds nor the formed shells are distributed symmetrically around the HII region, net forces pulling out the stars are present, accelerating them towards the edges of the cavity. We call this mechanism ``gravitational feedback", in which the gravity from the expelled gas appears to be the crucial mechanism producing unbound clusters that expand away from their formation centre in an accelerated way in young stellar clusters. This mechanism naturally explains the "Hubble flow-like" expansion observed in several young clusters.

High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2, harbors one of the closest such systems, making it an excellent target for case studies. We investigate the morphology, stability and dynamical properties of the hub-filament system on basis of 13CO and C18O observations obtained with the IRAM-30m telescope and H2 column density maps derived from Herschel dust emission observations. We identified the filamentary network and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100 Msun/pc and show signs of fragmentation. The secondary filaments have line masses of 12-60 Msun/pc and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally super-critical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly trans-critical regardless of whether non-thermal motions are included or not. From the main filaments, we estimate a mass accretion rate of 10(-4)-10(-3) Msun/pc into the hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4x10(-4) Msun/pc. The main filaments extend into the hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas. We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~2.5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the HII region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern.

We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventually stops the accretion (through filaments) toward the massive star-forming regions. The new over-pressured H\,{\normalsize II} regions push away the dense gas, thus disrupting the more massive collapse centres. Also, because of the complex density structure in the cloud, the H\,{\normalsize II} regions expand in a hybrid manner: they virtually do not expand toward the densest regions (cores), while they expand according to the classical analytical result towards the rest of the cloud, and in an accelerated way, as a blister region, towards the diffuse medium. Thus, the ionised regions grow anisotropically, and the ionising stars generally appear off-centre of the regions. Finally, we find that the hypotheses assumed in standard H\,{\normalsize II}-region expansion models (fully embedded region, blister-type, or expansion in a density gradient) apply simultaneously in different parts of our simulated H\,{\normalsize II} regions, producing a net expansion law ($R \propto t^\alpha$, with $\alpha$ in the range of 0.93-1.47 and a mean value of $1.2 \pm 0.17$) that differs from any of those of the standard models.

Numerical simulations of star formation have found that a power-law mass function can develop at high masses. In a previous paper, we employed isothermal simulations which created large numbers of sinks over a large range in masses to show that the power law exponent of the mass function, $dN/d\log M \propto M^{\Gamma}$, asymptotically and accurately approaches $\Gamma = -1.$ Simple analytic models show that such a power law can develop if the mass accretion rate $\dot{M} \propto M^2$, as in Bondi-Hoyle accretion; however, the sink mass accretion rates in the simulations show significant departures from this relation. In this paper we show that the expected accretion rate dependence is more closely realized provided the gravitating mass is taken to be the sum of the sink mass and the mass in the near environment. This reconciles the observed mass functions with the accretion rate dependencies, and demonstrates that power-law upper mass functions are essentially the result of gravitational focusing, a mechanism present in, for example, the competitive accretion model.

We present \emph{Herschel} PACS observations of 8 Classical T Tauri Stars in the $\sim 7-10$ Myr old OB1a and the $\sim 4-5$ Myr old OB1b Orion sub-asscociations. Detailed modeling of the broadband spectral energy distributions, particularly the strong silicate emission at 10 $\mu$m, shows that these objects are (pre)transitional disks with some amount of small optically thin dust inside their cavities, ranging from $\sim 4$ AU to $\sim 90$ AU in size. We analyzed \emph{Spitzer} IRS spectra for two objects in the sample: CVSO-107 and CVSO-109. The IRS spectrum of CVSO-107 indicates the presence of crystalline material inside its gap while the silicate feature of CVSO-109 is characterized by a pristine profile produced by amorphous silicates; the mechanisms creating the optically thin dust seem to depend on disk local conditions. Using millimeter photometry we estimated dust disk masses for CVSO-107 and CVSO-109 lower than the minimum mass of solids needed to form the planets in our Solar System, which suggests that giant planet formation should be over in these disks. We speculate that the presence and maintenance of optically thick material in the inner regions of these pre-transitional disks might point to low-mass planet formation.

We study the dynamical state of cores by using a simple analytical model, a sample of observational massive cores, and numerical simulations of collapsing massive cores. From the model, we find that, if cores are formed from turbulent compressions, they evolve from small to large column densities, increasing their velocity dispersion as they collapse, while they tend to equipartition between gravity and kinetic energy. From the observational sample, we find that: (a) cores with substantially different column densities in the sample do not follow a Larson-like linewidth-size relation. Instead, cores with higher column densities tend to be located in the upper-left corner of the Larson velocity dispersion-size diagram, a result predicted previously (Ballesteros-Paredes et al. 2011a). (b) The data exhibit cores with overvirial values. Finally, in the simulations we reproduce the behavior depicted by the model and observational sample: cores evolve towards larger velocity dispersions and smaller sizes as they collapse and increase their column density. However, collapsing cores appear to approach overvirial states within a free-fall time. The cause of this apparent excess of kinetic energy is an underestimation of the gravitational energy, due to the assumption that the gravitational energy is given by the energy of an isolated sphere of constant column density. This excess disappears when the gravitational energy is correctly calculated from the actual spatial mass distribution, where inhomogeneities, as well as the potential due to the mass outside of the core, also contribute to the gravitational energy. We conclude that the observed energy budget of cores in surveys is consistent with their non-thermal motions being driven by their self-gravity and in a hierarchical and chaotic collapse.

The origin of the observed morphological and kinematic substructure of young star forming regions is a matter of debate. We offer a new analysis of data from simulations of globally gravitationally collapsing clouds of progenitor gas to answer questions about sub-structured star formation in the context of cold collapse. As a specific example, we compare our models to recent radial velocity survey data from the IN-SYNC survey of Orion and new observations of dense gas kinematics, and offer possible interpretations of kinematic and morphological signatures in the region. In the context of our model, we find the frequently-observed hub-filament morphology of the gas naturally arises during gravitational evolution, as well as the dynamically-distinct kinematic substructure of stars. We emphasize that the global and not just the local gravitational potential plays an important role in determining the dynamics of both clusters and filaments.

We analyse the morphology and kinematics of dense filamentary structures produced in a numerical simulation of a star--forming cloud of $1.4 \times 10^4 \, \Msun$ evolving under their own self--gravity in a magnetized medium. This study is motivated by recent observations of velocity--coherent substructures ("fibres") in star-forming filaments. We find such "fibres" ubiquitously in our simulated filament. We found that a fibre in one projection is not necessarily a fibre in another projection, and thus, caution should be taken into account when considering them as real objects. We found that only the densest parts of the filament ($\sim$30\% of the densest volume, which contains $\sim$70\% of the mass) belong to fibres in 2 projections. Moreover, it is quite common that they are formed by separated density enhancements superposed along the line of sight. Observations of fibres can yield insight into the level of turbulent substructure driven by gravity, but care should be taken in interpreting the results given the problem of line of sight superposition. We also studied the morphology and kinematics of the 3D skeleton (spine), finding that subfilaments accrete structured material mainly along the magnetic field lines, which are preferentially perpendicular to the skeleton. The magnetic field is at the same time dragged by the velocity field due to the gravitational collapse.

Herschel has revolutionized our ability to measure column densities (N$_{\rm H}$) and temperatures (T) of molecular clouds thanks to its far infrared multiwavelength coverage. However, the lack of a well defined background intensity level in the Herschel data limits the accuracy of the N$_{\rm H}$ and T maps. We provide a method that corrects the missing Herschel background intensity levels using the Planck model for foreground Galactic thermal dust emission. We present a Fourier method that combines the publicly available Planck model on large angular scales with the Herschel images on smaller angular scales. We apply our method to two regions spanning a range of Galactic environments: Perseus and the Galactic plane region around $l = 11\deg$ (HiGal--11). We post-process the combined dust continuum emission images to generate column density and temperature maps. We compare these to previously adopted constant--offset corrections. We find significant differences ($\gtrsim$20\%) over significant ($\sim$15\%) areas of the maps, at low column densities ($N_{\rm H}\lesssim10^{22}$\,cm$^{-2}$) and relatively high temperatures ($T\gtrsim20$\,K). We also apply our method to synthetic observations of a simulated molecular cloud to validate our method. Our method successfully corrects the Herschel images, including both the constant--offset intensity level and the scale-dependent background variations measured by Planck. Our method improves the previous constant--offset corrections, which did not account for variations in the background emission levels.

We present Submillimeter Array (SMA) observations at 345 GHz towards the intermediate/high-mass cluster-forming region NGC 6334 V. From the dust emission we spatially resolve three dense condensations, the brightest one presenting the typical chemistry of a hot core. The magnetic field (derived from the dust polarized emission) shows a bimodal converging pattern towards the hot core. The molecular emission traces two filamentary structures at two different velocities, separated by 2 km/s, converging to the hot core and following the magnetic field distribution. We compare the velocity field and the magnetic field derived from the SMA observations with MHD simulations of star-forming regions dominated by gravity. This comparison allows us to show how the gas falls in from the larger-scale extended dense core (~0.1 pc) of NGC 6334 V towards the higher-density hot core region (~0.02 pc) through two distinctive converging flows dragging the magnetic field, whose strength seems to have been overcome by gravity.

We analyze the physical properties and energy balance of density enhancements in two SPH simulations of the formation, evolution, and collapse of giant molecular clouds. In the simulations, no feedback is included, so all motions are due either to the initial, decaying turbulence, or to gravitational contraction. We define clumps as connected regions above a series of density thresholds. The resultingfull set of clumps follows the generalized energy-equipartition relation $\sigma_{v}/R^{1/2} \propto \Sigma^{1/2}$, where $\sigma_{v}$ is the velocity dispersion, $R$ is the "radius", and $\Sigma$ is the column density. We interpret this as a natural consequence of gravitational contraction at all scales, rather than virial equilibrium. Nevertheless, clumps with low $\Sigma$ tend to show a large scatter around equipartition. In more than half of the cases, this scatter is dominated by external turbulent compressions that {\it assemble} the clumps, rather than by small-scale random motions that would disperse them. The other half does actually disperse. Moreover, clump sub-samples selected by means of different criteria exhibit different scalings. Sub-samples with narrow $\Sigma$ ranges follow Larson-like relations, although characterized by their respective value of $\Sigma$. Finally, we find that: i) clumps lying in filaments tend to appear sub-virial; ii) high-density cores ($n \ge 10^5$ cm$^3$) that exhibit moderate kinetic energy excesses often contain sink ("stellar") particles, and the excess disappears when the stellar mass is taken into account in the energy balance; iii) cores with kinetic energy excess but no stellar particles are truly in a state of dispersal.

We present new Herschel PACS observations of 32 T Tauri stars in the young ($\sim$3 Myr) $\sigma$ Ori cluster. Most of our objects are K & M stars with large excesses at 24 $\mu$m. We used irradiated accretion disk models of D'Alessio et al. (2006) to compare their spectral energy distributions with our observational data. We arrive at the following six conclusions. (i) The observed disks are consistent with irradiated accretion disks systems. (ii) Most of our objects (60%) can be explained by significant dust depletion from the upper disk layers. (iii) Similarly, 61% of our objects can be modeled with large disk sizes ($\rm R_{\rm d} \geq$ 100 AU). (iv) The masses of our disks range between 0.03 to 39 $\rm M_{Jup}$, where 35% of our objects have disk masses lower than 1 Jupiter. Although these are lower limits, high mass ($>$ 0.05 M$_{\odot}$) disks, which are present e.g, in Taurus, are missing. (v) By assuming a uniform distribution of objects around the brightest stars at the center of the cluster, we found that 80% of our disks are exposed to external FUV radiation of $300 \leq G_{0} \leq 1000$, which can be strong enough to photoevaporate the outer edges of the closer disks. (vi) Within 0.6 pc from $\sigma$ Ori we found forbidden emission lines of [NII] in the spectrum of one of our large disk (SO662), but no emission in any of our small ones. This suggests that this object may be an example of a photoevaporating disk.

Sub-virial gravitational collapse is one mechanism by which star clusters may form. Here we investigate whether this mechanism can be inferred from observations of young clusters. To address this question, we have computed SPH simulations of the initial formation and evolution of a dynamically young star cluster through cold (sub-virial) collapse, starting with an ellipsoidal, turbulently seeded distribution of gas, and forming sink particles representing (proto)stars. While the initial density distributions of the clouds do not have large initial mass concentrations, gravitational focusing due to the global morphology leads to cluster formation. We use the resulting structures to extract observable morphological and kinematic signatures for the case of sub-virial collapse. We find that the signatures of the initial conditions can be erased rapidly as the gas and stars collapse, suggesting that kinematic observations need to be made either early in cluster formation and/or at larger scales, away from the growing cluster core. Our results emphasize that a dynamically young system is inherently evolving on short timescales, so that it can be highly misleading to use current-epoch conditions to study aspects such as star formation rates as a function of local density. Our simulations serve as a starting point for further studies of collapse including other factors such as magnetic fields and stellar feedback.

We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected when the gravitationally unstable mass is calculated including various prescriptions for `turbulent support'. On the other hand, the best correlation is found for the case of pure thermal Jeans fragmentation, for which we infer a core formation efficiency around 13 per cent, consistent with previous works. We conclude that the dominant factor determining the fragmentation level of star-forming massive dense cores at 0.1 pc scale seems to be thermal Jeans fragmentation.

We report on a series of numerical simulations of gas clouds with self-gravity forming sink particles, adopting an isothermal equation of state to isolate the effects of gravity from thermal physics on the resulting sink mass distributions. Simulations starting with supersonic velocity fluctuations develop sink mass functions with a high-mass power-law tail $dN/d\log M \propto M^{\Gamma}$, $\Gamma = -1 \pm 0.1$, independent of the initial Mach number of the velocity field. Similar results but with weaker statistical significance hold for a simulation starting with initial density fluctuations. This mass function power-law dependence agrees with the asymptotic limit found by Zinnecker assuming Bondi-Hoyle-Littleton (BHL) accretion, even though the mass accretion rates of individual sinks show significant departures from the predicted $\mdot \propto M^2$ behavior. While BHL accretion is not strictly applicable due to the complexity of the environment, we argue that the final mass functions are the result of a {\em relative} $M^2$ dependence resulting from gravitationally-focused accretion. Our simulations may show the power-law mass function particularly clearly compared with others because our adoption of an isothermal equation of state limits the effects of thermal physics in producing a broad initial fragmentation spectrum; $\Gamma \rightarrow -1$ is an asymptotic limit found only when sink masses grow well beyond their initial values. While we have purposely eliminated many additional physical processes (radiative transfer, feedback) which can affect the stellar mass function, our results emphasize the importance of gravitational focusing for massive star formation.

We present a study of 15 new brown dwarfs belonging to the $\sim7$ Myr old 25 Orionis group and Orion OB1a sub-association with spectral types between M6 and M9 and estimated masses between $\sim0.07$M$_\odot$ and $\sim0.01$ M$_\odot$. By comparing them through a Bayesian method with low mass stars ($0.8\lesssim$ M/M$_\odot\lesssim0.1$) from previous works in the 25 Orionis group, we found statistically significant differences in the number fraction of classical T Tauri stars, weak T Tauri stars, class II, evolved discs and purely photospheric emitters at both sides of the sub-stellar mass limit. Particularly we found a fraction of $3.9^{+2.4}_{-1.6}~\%$ low mass stars classified as CTTS and class II or evolved discs, against a fraction of $33.3^{+10.8}_{-9.8}~\%$ in the sub-stellar mass domain. Our results support the suggested scenario in which the dissipation of discs is less efficient for decreasing mass of the central object.

We present synthetic Hi and CO observations of a simulation of decaying turbulence in the thermally bistable neutral medium. We first present the simulation, with clouds initially consisting of clustered clumps. Self-gravity causes these clump clusters to form more homogeneous dense clouds. We apply a simple radiative transfer algorithm, and defining every cell with <Av> > 1 as molecular. We then produce maps of Hi, CO-free molecular gas, and CO, and investigate the following aspects: i) The spatial distribution of the warm, cold, and molecular gas, finding the well-known layered structure, with molecular gas surrounded by cold Hi, surrounded by warm Hi. ii) The velocity of the various components, with atomic gas generally flowing towards the molecular gas, and that this motion is reflected in the frequently observed bimodal shape of the Hi profiles. This conclusion is tentative, because we do not include feedback. iii) The production of Hi self-absorption (HISA) profiles, and the correlation of HISA with molecular gas. We test the suggestion of using the second derivative of the brightness temperature Hi profile to trace HISA and molecular gas, finding limitations. On a scale of ~parsecs, some agreement is obtained between this technique and actual HISA, as well as a correlation between HISA and N(mol). It quickly deteriorates towards sub-parsec scales. iv) The N-PDFs of the actual Hi gas and those recovered from the Hi line profiles, with the latter having a cutoff at column densities where the gas becomes optically thick, thus missing the contribution from the HISA-producing gas. We find that the power-law tail typical of gravitational contraction is only observed in the molecular gas, and that, before the power-law tail develops in the total gas density PDF, no CO is yet present, reinforcing the notion that gravitational contraction is needed to produce this component. (abridged)

We present an analysis of the dust continuum emission at 870 um in order to investigate the mass distribution of clumps within infrared dark clouds (IRDCs). We map six IRDCs with the Large APEX BOlometer CAmera (LABOCA) at APEX, reaching an rms noise level of 28-44 mJy/beam. The dust continuum emission coming from these IRDCs was decomposed by using two automated algorithms, Gaussclumps and Clumpfind. Moreover, we carried out single-pointing observations of the N_2H^+ (3-2) line toward selected positions to obtain kinematic information. The mapped IRDCs are located in the range of kinematic distances of 2.7-3.2 kpc. We identify 510 and 352 sources with Gaussclumps and Clumpfind, respectively, and estimate masses and other physical properties assuming a uniform dust temperature. The mass ranges are 6-2692 Msun (Gaussclumps) and 7-4254 Msun (Clumpfind) and the ranges in effective radius are around 0.10-0.74 pc (Gaussclumps) and 0.16-0.99 pc (Clumpfind). The mass distribution, independent of the decomposition method used, is fitted by a power law, dN/dM propto M^alpha, with an index of -1.60 +/- 0.06, consistent with the CO mass distribution and other high-mass star-forming regions.

Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field.

We study the effect that the dark matter background (DMB) has on the gravitational energy content and, in general, on the star formation efficiency of a molecular cloud (MC). We first analyze the effect that a dark matter halo, described by the Navarro et al. (1996) density profile, has on the energy budget of a spherical, homogeneous, cloud located at different distances from the halo center. We found that MCs located in the innermost regions of a massive galaxy can feel a contraction force greater than their self-gravity due to the incorporation of the potential of the galaxy's dark matter halo. We also calculated analytically the gravitational perturbation that a MC produces over a uniform DMB (uniform at the scales of a MC) and how this perturbation will affect the evolution of the MC itself. The study shows that the star formation in a MC will be considerably enhanced if the cloud is located in a dense and low velocity dark matter environment. We confirm our results by measuring the star formation efficiency in numerical simulations of the formation and evolution of MCs within different DMBs. Our study indicates that there are situations where the dark matter's gravitational contribution to the evolution of the molecular clouds should not be neglected.

Most young stars in nearby molecular clouds have estimated ages of 1-2 Myr, suggesting that star formation is rapid. However, small numbers of stars in these regions with inferred ages of >= 5-10 Myr have been cited to argue that star formation is instead a slow, quasi-static process. When considering these alternative pictures it is important to recognize that the age spread in a given star-forming cloud is necessarily an upper limit to the timescales of local collapse, as not all spatially-distinct regions will start contracting at precisely the same instant. Moreover, star-forming clouds may dynamically evolve on timescales of a few Myr; in particular, global gravitational contraction will tend to yield increasing star formation rates with time due to generally increasing local gas densities. We show that two different numerical simulations of dynamic, flow-driven molecular cloud formation and evolution 1) predict age spreads for the main stellar population roughly consistent with observations, and 2) raise the possibility of forming small numbers of stars early in cloud evolution, before global contraction concentrates the gas and the bulk of the stellar population is produced. In general, the existence of a small number of older stars among a generally much-younger population is consistent with the picture of dynamic star formation, and may even provide clues to the time evolution of star-forming clouds.

It has been recently shown that molecular clouds do not exhibit a unique shape for the column density probability distribution function (Npdf). Instead, clouds without star formation seem to possess a lognormal distribution, while clouds with active star formation develope a power-law tail at high column densities. The lognormal behavior of the Npdf has been interpreted in terms of turbulent motions dominating the dynamics of the clouds, while the power-law behavior occurs when the cloud is dominated by gravity. In the present contribution we use thermally bi-stable numerical simulations of cloud formation and evolution to show that, indeed, these two regimes can be understood in terms of the formation and evolution of molecular clouds: a very narrow lognormal regime appears when the cloud is being assembled. However, as the global gravitational contraction occurs, the initial density fluctuations are enhanced, resulting, first, in a wider lognormal Npdf, and later, in a power-law Npdf. We thus suggest that the observed Npdf of molecular clouds are a manifestation of their global gravitationally contracting state. We also show that, contrary to recent suggestions, the exact value of the power-law slope is not unique, as it depends on the projection in which the cloud is being observed.

We discuss the nature of the velocity dispersion vs. size relation for molecular clouds. In particular, we add to previous observational results showing that the velocity dispersions in molecular clouds and cores are not purely functions of spatial scale but involve surface gas densities as well. We emphasize that hydrodynamic turbulence is required to produce the first condensations in the progenitor medium. However, as the cloud is forming, it also becomes bound, and gravitational accelerations dominate the motions. Energy conservation in this case implies $|E_g| \sim E_k$, in agreement with observational data, and providing an interpretation for two recent observational results: the scatter in the $\delta v-R$ plane, and the dependence of the velocity dispersion on the surface density ${\delta v^2/ R} \propto \Sigma$. We argue that the observational data are consistent with molecular clouds in a state of hierarchical gravitational collapse, i.e., developing local centers of collapse throughout the whole cloud while the cloud itself is collapsing, and making equilibrium unnecessary at all stages prior to the formation of actual stars. Finally, we discuss how this mechanism need not be in conflict with the observed star formation rate.

We present optical and IR integrated colours and SBF magnitudes, computed from stellar population synthesis models that include emission from the dusty envelopes surrounding TP-AGB stars undergoing mass-loss. We explore the effects of varying the mass-loss rate by one order of magnitude around the fiducial value, modifying accordingly both the stellar parameters and the output spectra of the TP-AGB stars plus their dusty envelopes. The models are single burst, and range in age from a few Myr to 14 Gyr, and in metallicity between $Z$ = 0.0001 and $Z$ = 0.07; they combine new calculations for the evolution of stars in the TP-AGB phase, with star plus envelope SEDs produced with the radiative transfer code DUSTY. We compare these models to optical and near-IR data of single AGB stars and Magellanic star clusters. This comparison validates the current understanding of the role of mass-loss in determining stellar parameters and spectra in the TP-AGB. However, neither broad-band colours nor SBF measurements in the optical or the near-IR can discern global changes in the mass-loss rate of a stellar population. We predict that mid-IR SBF measurements can pick out such changes, and actually resolve whether a relation between metallicity and mass-loss exists.

We investigate the properties of "star forming regions" in a previously published numerical simulation of molecular cloud formation out of compressive motions in the warm neutral atomic interstellar medium, neglecting magnetic fields and stellar feedback. In this simulation, the velocity dispersions at all scales are caused primarily by infall motions rather than by random turbulence. We study the properties (density, total gas+stars mass, stellar mass, velocity dispersion, and star formation rate) of the cloud hosting the first local, isolated "star formation" event in the simulation and compare them with those of the cloud formed by a later central, global collapse event. We suggest that the small-scale, isolated collapse may be representative of low- to intermediate-mass star-forming regions, while the large-scale, massive one may be representative of massive star forming regions. We also find that the statistical distributions of physical properties of the dense cores in the region of massive collapse compare very well with those from a recent survey of the massive star forming region in the Cygnus X molecular cloud. The star formation efficiency per free-fall time (SFE_ff) of the high-mass SF clump is low, ~0.04. This occurs because the clump is accreting mass at a high rate, not because its specific SFR (SSFR) is low. This implies that a low value of the SFE_ff does not necessarily imply a low SSFR, but may rather indicate a large gas accretion rate. We suggest that a globally low SSFR at the GMC level can be attained even if local star forming sites have much larger values of the SSFR if star formation is a spatially intermittent process, so that most of the mass in a GMC is not participating of the SF process at any given time.

Recent models of molecular cloud formation and evolution suggest that such clouds are dynamic and generally exhibit gravitational collapse. We present a simple analytic model of global collapse onto a filament and compare this with our numerical simulations of the flow-driven formation of an isolated molecular cloud to illustrate the supersonic motions and infall ram pressures expected in models of gravity-driven cloud evolution. We apply our results to observations of the Pipe Nebula, an especially suitable object for our purposes as its low star formation activity implies insignifcant perturbations from stellar feedback. We show that our collapsing cloud model can explain the magnitude of the velocity dispersions seen in the $^{13}$CO filamentary structure by Onishi et al. and the ram pressures required by Lada et al. to confine the lower-mass cores in the Pipe nebula. We further conjecture that higher-resolution simulations will show small velocity dispersions in the densest core gas, as observed, but which are infall motions and not supporting turbulence. Our results point out the inevitability of ram pressures as boundary conditions for molecular cloud filaments, and the possibility that especially lower-mass cores still can be accreting mass at significant rates, as suggested by observations.

Only a few molecular clouds in the Solar Neighborhood exhibit the formation of only low-mass stars. Traditionally, these clouds have been assumed to be supported against more vigorous collapse by magnetic fields. The existence of strong magnetic fields in molecular clouds, however, poses serious problems for the formation of stars and of the clouds themselves. In this {\em Letter}, we review the three-dimensional structure and kinematics of Taurus --the archetype of a region forming only low-mass stars-- as well as its orientation within the Milky way. We conclude that the particularly low star-formation efficiency in Taurus may naturally be explained by tidal forces from the Galaxy, with no need for magnetic regulation or stellar feedback.

(Abridged) The gravitational term for clouds and cores entering in the virial theorem is usually assumed to be equal to the gravitational energy, since the contribution to the gravitational force from the mass distribution outside the volume of integration is assumed to be negligible. Such approximation may not be valid in the presence of an important external net potential. In the present work we analyze the effect of an external gravitational field on the gravitational budget of a density structure. Our cases under analysis are (a) a giant molecular cloud (GMC) with different aspect ratios embedded within a galactic net potential, and (b) a molecular cloud core embedded within the gravitational potential of its parent molecular cloud. We find that for roundish GMCs, the tidal tearing due to the shear in the plane of the galaxy is compensated by the tidal compression in the z direction. The influence of the external effective potential on the total gravitational budget of these clouds is relatively small, although not necessarily negligible. However, for more filamentary GMCs, the external effective potential can be dominant and can even overwhelm self-gravity, regardless of whether its main effect on the cloud is to disrupt it or compress it. This may explain the presence of some GMCs with few or no signs of massive star formation, such as the Taurus or the Maddalena's clouds. In the case of dense cores embedded in their parent molecular cloud, we found that the gravitational content due to the external field may be more important than the gravitational energy of the cores themselves. This effect works in the same direction as the gravitational energy, i.e., favoring the collapse of cores. We speculate on the implications of these results for star formation models.

We present three numerical simulations of randomly driven, isothermal, non-magnetic, self-gravitating turbulence with different rms Mach numbers Ms and physical sizes L, but approximately the same value of the virial parameter, alpha approx 1.2. We obtain the following results: a) We test the hypothesis that the collapsing centers originate from locally Jeans-unstable ("super-Jeans"), subsonic fragments; we find no such structures. b) We find that the fraction of small-scale super-Jeans structures is larger in the presence of self-gravity. c) The velocity divergence of subregions of the simulations exhibits a negative correlation with their mean density. d) The density probability density function (PDF) deviates from a lognormal in the presence of self-gravity. e) Turbulence alone in the large-scale simulation does not produce regions with the same size and mean density as those of the small-scale simulation. Items (b)-(e) suggest that self-gravity is not only involved in causing the collapse of Jeans-unstable density fluctuations produced by the turbulence, but also in their {it formation}. We also measure the star formation rate per free-fall time, as a function of Ms for the three runs, and compare with the predictions of recent semi-analytical models. We find marginal agreement to within the uncertainties of the measurements. However, the hypotheses of those models neglect the net negative divergence of dense regions we find in our simulations. We conclude that a) part of the observed velocity dispersion in clumps must arise from clump-scale inwards motions, and b) analytical models of clump and star formation need to take into account this dynamical connection with the external flow and the fact that, in the presence of self-gravity, the density PDF may deviate from a lognormal.

We present preliminary numerical evidence that the physical conditions in high-mass star forming regions can arise from global gravitational infall, with the velocity dispersions being caused primarily by infall motions rather than random turbulence. To this end, we study the clumps and cores appearing in the region of central collapse in a numerical simulation of the formation, evolution, and subsequent collapse of a dense cloud out of a transonic compression in the diffuse atomic ISM. The clumps have sizes $\sim 1$ pc, masses of several hundred $M_\odot$, and three-dimensional velocity dispersions $\sim 3$ km s$^{-1}$, in agreement with typical observed values for such structures. The clumps break down into massive cores of sizes $\sim 0.1$ pc, densities $\sim 10^5$, masses 2-300 $M_\odot$, with distributions of these quantities that peak at the same values as the massive core sample in a recent survey of the Cygnus X molecular cloud complex. Although preliminary, these results suggest that high-mass star forming clumps may be in a state of global gravitational collapse rather than in equilibrium supported by strong turbulence.

(Abridged) We present numerical hydrodynamical simulations of the formation, evolution and gravitational collapse of isothermal molecular cloud cores. A compressive wave is set up in a constant sub-Jeans density distribution of radius r = 1 pc. As the wave travels through the simulation grid, a shock-bounded spherical shell is formed. The inner shock of this shell reaches and bounces off the center, leaving behind a central core with an initially almost uniform density distribution, surrounded by an envelope consisting of the material in the shock-bounded shell, with a power-law density profile that at late times approaches a logarithmic slope of -2 even in non-collapsing cases. The resulting density structure resembles a quiescent core of radius < 0.1 pc, with a Bonnor-Ebert-like (BE-like) profile, although it has significant dynamical differences: it is initially non-self-gravitating and confined by the ram pressure of the infalling material, and consequently, growing continuously in mass and size. With the appropriate parameters, the core mass eventually reaches an effective Jeans mass, at which time the core begins to collapse. Thus, there is necessarily a time delay between the appearance of the core and the onset of its collapse, but this is not due to the dissipation of its internal turbulence as it is often believed. These results suggest that pre-stellar cores may approximate Bonnor-Ebert structures which are however of variable mass and may or may not experience gravitational collapse, in qualitative agreement with the large observed frequency of cores with BE-like profiles.

We present measurements of the mean dense core lifetimes in numerical simulations of magnetically supercritical, turbulent, isothermal molecular clouds, in order to compare with observational determinations. "Prestellar" lifetimes (given as a function of the mean density within the cores, which in turn is determined by the density threshold n_thr used to define them) are consistent with observationally reported values, ranging from a few to several free-fall times. We also present estimates of the fraction of cores in the "prestellar", "stellar'', and "failed" (those cores that redisperse back into the environment) stages as a function of n_thr. The number ratios are measured indirectly in the simulations due to their resolution limitations. Our approach contains one free parameter, the lifetime of a protostellar object t_yso (Class 0 + Class I stages), which is outside the realm of the simulations. Assuming a value t_yso = 0.46 Myr, we obtain number ratios of starless to stellar cores ranging from 4-5 at n_thr = 1.5 x 10^4 cm^-3 to 1 at n_thr = 1.2 x 10^5 cm^-3, again in good agreement with observational determinations. We also find that the mass in the failed cores is comparable to that in stellar cores at n_thr = 1.5 x 10^4 cm^-3, but becomes negligible at n_thr = 1.2 x 10^5 cm^-3, in agreement with recent observational suggestions that at the latter densities the cores are in general gravitationally dominated. We conclude by noting that the timescale for core contraction and collapse is virtually the same in the subcritical, ambipolar diffusion-mediated model of star formation, in the model of star formation in turbulent supercritical clouds, and in a model intermediate between the previous two, for currently accepted values of the clouds' magnetic criticality.

We have mapped the starless core TMC-1C in a variety of molecular lines with the IRAM 30m telescope. High density tracers show clear signs of self-absorption and sub-sonic infall asymmetries are present in N2H+ (1-0) and DCO+ (2-1) lines. The inward velocity profile in N2H+ (1-0) is extended over a region of about 7,000 AU in radius around the dust continuum peak, which is the most extended ``infalling'' region observed in a starless core with this tracer. The kinetic temperature (~12 K) measured from C17O and C18O suggests that their emission comes from a shell outside the colder interior traced by the mm continuum dust. The C18O (2-1) excitation temperature drops from 12 K to ~10 K away from the center. This is consistent with a volume density drop of the gas traced by the C18O lines, from ~4x10^4 cm^-3 towards the dust peak to ~6x10^3 cm^-3 at a projected distance from the dust peak of 80" (or 11,000 AU). The column density implied by the gas and dust show similar N2H+ and CO depletion factors (f_D < 6). This can be explained with a simple scenario in which: (i) the TMC-1C core is embedded in a relatively dense environment (H2 ~10^4 cm^-3), where CO is mostly in the gas phase and the N2H+ abundance had time to reach equilibrium values; (ii) the surrounding material (rich in CO and N2H+) is accreting onto the dense core nucleus; (iii) TMC-1C is older than 3x10^5 yr, to account for the observed abundance of N2H+ across the core (~10^-10 w.r.t. H2); and (iv) the core nucleus is either much younger (~10^4 yr) or ``undepleted'' material from the surrounding envelope has fallen towards it in the past 10,000 yr.

We study the formation of giant dense cloud complexes and of stars within them by means of SPH numerical simulations of the mildly supersonic collision of gas streams (``inflows'') in the warm neutral medium (WNM). The resulting compressions cause cooling and turbulence generation in the gas, forming a cloud that then becomes self-gravitating and undergoes global collapse. Simultaneously, the turbulent, nonlinear density fluctuations induce fast, local collapse events. The simulations show that: a) The clouds are not in a state of equilibrium. Instead, they undergo secular evolution. Initially, their mass and gravitational energy |Eg| increase steadily, while the turbulent energy Ek reaches a plateau. b) When |Eg| becomes comparable to Ek, global collapse begins, causing a simultaneous increase in both that maintains a near-equipartition condition |Eg| ~ 2 Ek. c) Longer inflow durations delay the onset of global and local collapse, by maintaining a higher turbulent velocity dispersion in the cloud over longer times. d) The star formation rate is large from the beginning, without any period of slow and accelerating star formation. e) The column densities of the local star-forming clumps are very similar to reported values of the column density required for molecule formation, suggesting that locally molecular gas and star formation occur nearly simultaneously. The MC formation mechanism discussed here naturally explains the apparent ``virialized'' state of MCs and the ubiquitous presence of HI halos around them. Within their assumptions, our simulations support the scenario of rapid star formation after MCs are formed, although long (>~ 15 Myr) accumulation periods do occur during which the clouds build up their gravitational energy, and which are expected to be spent in the atomic phase.

We discuss problems with some observational estimates indicating long protostellar core lifetimes and large stellar age spreads in molecular clouds. We also point out some additional observational constraints which suggest that protostellar cores do not have long lifetimes before collapsing. For external galaxies, we argue that the widths of spiral arms does not imply a long star-formation process, since the formation of massive stars will disrupt molecular clouds, move material around, compress it in other regions which produce new star-forming clouds. Thus, it seems unavoidable that this cyclical process will result in an extended period of enhanced star formation, which does not represent the survival time of any individual molecular cloud. We argue that the rapid star formation indicated observationally is also easier to understand theoretically than the traditional scenario of slow quasi-static contraction with ambipolar diffusion.

It has been paid little or no attention to the implications that turbulent fragmentation has on the validity of at least six common assumptions on the Virial Theorem (VT), which are: (i) the only role of turbulent motions within a cloud is to provide support against collapse, (ii) the surface terms are negligible compared to the volumetric ones, (iii) the gravitational term is a binding source for the clouds, (iv) the sign of the second-time derivative of the moment of inertia determines whether the cloud is contracting or expanding, (v) interstellar clouds are in Virial Equilibrium (VE), and (vi) Larson's (1981) relations are the observational proof that clouds are in VE. Interstellar clouds cannot fulfill these assumptions, however, because turbulent fragmentation will induce flux of mass, moment and energy between the clouds and their environment, and will favor local collapse while may disrupt the clouds within a dynamical timescale. It is argued that, although the observational and numerical evidence suggests that interstellar clouds are not in VE, the so-called ``Virial Mass'' estimations, which actually should be called ``energy-equipartition mass'' estimations, are good order-of magnitude estimations of the actual mass of the clouds just because observational surveys will tend to detect interstellar clouds appearing to be close to energy equipartition. However, since clouds are actually out of VE, as suggested by asymmetrical line profiles, they should be transient entities. These results are compatible with observationally-based estimations for rapid star formation. , and call into question the models for the star formation efficiency based on clouds being in VE.

We review the properties of turbulent molecular clouds (MCs), focusing on the physical processes that influence star formation (SF). MC formation appears to occur during large-scale compression of the diffuse ISM driven by supernovae, magnetorotational instability, or gravitational instability in galactic disks of stars and gas. The compressions generate turbulence that can accelerate molecule production and produce the observed morphology. We then review the properties of MC turbulence, including density enhancements observed as clumps and cores, magnetic field structure, driving scales, the relation to observed scaling relations, and the interaction with gas thermodynamics. We argue that MC cores are dynamical, not quasistatic, objects with relatively short lifetimes not exceeding a few megayears. We review their morphology, magnetic fields, density and velocity profiles, and virial budget. Next, we discuss how MC turbulence controls SF. On global scales turbulence prevents monolithic collapse of the clouds; on small scales it promotes local collapse. We discuss its effects on the SF efficiency, and critically examine the possible relation between the clump mass distribution and the initial mass function, and then turn to the redistribution of angular momentum during collapse and how it determines the multiplicity of stellar systems. Finally, we discuss the importance of dynamical interactions between protostars in dense clusters, and the effect of the ionization and winds from those protostars on the surrounding cloud. We conclude that the interaction of self-gravity and turbulence controls MC formation and behavior, as well as the core and star formation processes within them.

We investigate the core mass distribution (CMD) resulting from numerical models of turbulent fragmentation of molecular clouds. In particular we study its dependence on the sonic root-mean-square Mach number $\Ms$. We analyze simulations with $\Ms$ ranging from 1 to 15 to show that, as $\Ms$ increases, the number of cores increases as well while their average mass decreases. This stems from the fact that high-Mach number flows produce many and strong shocks on intermediate to small spatial scales, leading to a highly-fragmented density structure. We also show that the CMD from purely turbulent fragmentation does not follow a single power-law, but it may be described by a function that changes continuously its shape, probably more similar to a log-normal function. The CMD in supersonic turbulent flows does not have a universal slope, and as consequence, cast some doubt on attempts to directly relate the CMD to a universal Initial Mass Function.

We present measurements of the star formation efficiency (SFE) in 3D numerical simulations of driven turbulence in supercritical, ideal-MHD, and non-magnetic regimes, characterized by their mean normalized mass-to-flux ratio $\mu$, all with 64 Jeans masses and similar rms Mach numbers ($\sim 10$). In most cases, the moderately supercritical runs with $\mu = 2.8$ have significantly lower SFEs than the non-magnetic cases, being comparable to observational estimates for whole molecular clouds ($\lesssim$ 5% over 4 Myr). Also, as the mean field is increased, the number of collapsed objects decreases, and the median mass of the collapsed objects increases. However, the largest collapsed-object masses systematically occur in the weak-field case $\mu = 8.8$. The high-density tails of the density histograms in the simulations are depressed as the mean magnetic field strength is increased. This suggests that the smaller numbers and larger masses of the collapsed objects in the magnetic cases may be due to a greater scarcity and lower mean densities (implying larger Jeans masses) of the collapse candidates. In this scenario, the effect of a weak field is to reduce the probability of a core reaching its thermal Jeans mass, even if it is supercritical. We thus suggest that the SFE may be monotonically reduced as the field strength increases from zero to subcritical values, rather than there being a discontinuous transition between the sub- and supercritical regimes, and that a crucial question to address is whether the turbulence in molecular clouds is driven or decaying, with current observational and theoretical evidence favoring (albeit inconclusively) the driven regime.

[abridged] We investigate the velocity structure of protostellar cores that result from non-magnetic numerical models of the gravoturbulent fragmentation of molecular cloud material. A large fraction of the cores analyzed are ``quiescent'', and more than half are identified as ``coherent''. The fact that dynamically evolving cores in highly supersonic turbulent flows can be quiescent may be understood because cores lie at the stagnation points of convergent turbulent flows, where compression is at a maximum, and relative velocity differences are at a minimum. The coherence may be due to an observational effect related to the length and concentration of the material contributing to the line. The velocity dispersion of the our cores often has its local maximum at small offsets from the column density maximum, suggesting that the core is the dense region behind a shock. Such a configuration is often found in observations of molecular cloud cores, and argues in favor of the gravoturbulent scenario of stellar birth as it is not expected in star-formation models based on magnetic mediation. Cores with collapsed objects tend to be near equipartition between their gravitational and kinetic energies, while cores without collapsed objects tend to be gravitationally unbound, suggesting that gravitational collapse occurs immediately after gravity becomes dominant. Finally, cores in simulations driven at large scales are more frequently quiescent and coherent, and have more realistic ratios of $M_{\rm vir}/M$, supporting the notion that molecular cloud turbulence is driven at large scales.

We discuss the lifetimes and evolution of clumps and cores formed as turbulent density fluctuations in nearly isothermal molecular clouds. In the non-magnetic case, clumps are unlikely to reach a hydrostatic state, and instead are expected to either proceed directly to collapse, or else ``rebound'' towards the mean pressure and density of the parent cloud. Rebounding clumps are delayed in their re-expansion by their self-gravity. From a simple virial calculation, we find re-expansion times of a few free-fall times. In the magnetic case, we present a series of driven-turbulence, ideal-MHD isothermal numerical simulations in which we follow the evolution of clumps and cores in relation to the magnetic criticality of their ``parent clouds'' (the numerical boxes). In subcritical boxes, magnetostatic clumps do not form. A few moderately-gravitationally bound clumps form which however are dispersed by the turbulence in < 1.3 Myr. An estimate of the ambipolar diffusion (AD) time scale t_AD in these cores gives t_AD > 1.3 Myr, only slightly longer than the dynamical times. In supercritical boxes, some cores become locally supercritical and collapse in typical times ~ 1 Myr. We also observe longer-lived supercritical cores that however do not collapse because they are smaller than the local Jeans length. Fewer clumps and cores form in these simulations than in their non-magnetic counterpart. Our results suggest that a) A fraction of the cores may not form stars, and may correspond to some of the observed starless cores. b) Cores may be out-of-equilibrium structures, rather than quasi-magnetostatic ones. c) The magnetic field may help reduce the star formation efficiency by reducing the probability of core formation, rather than by significantly delaying the collapse of individual cores.

We examine the formation of molecular gas behind shocks in atomic gas using a chemical/dynamical model, particular emphasis is given to constraints the chemistry places on the dynamical evolution. The most important result of this study is to stress the importance of shielding the molecular gas from the destructive effects of UV radiation. For shock ram pressures comparable to or exceeding typical local ISM pressures, self-shielding controls the formation time of H2 but CO formation requires shielding of the interstellar radiation field by dust grains. We find that the molecular hydrogen fractional abundance can become significant well before CO forms. The timescale for (CO) molecular cloud formation is not set by H2 formation, but rather by the timescale for accumulating a sufficient column density or extinction, A_V > 0.7. The local ratio of atomic to molecular gas (4:1), coupled with short estimates for cloud lifetimes (3-5 Myr), suggests that the timescales for accumulating molecular clouds from atomic material typically must be no longer than about 12-20 Myr. Based on the shielding requirement, this implies that the typical product of pre-shock density and velocity must be n*v > 20 cm^-3 km s^-1. Based on these results we find that flow-driven formation of molecular clouds in the local interstellar medium can occur sufficiently rapidly to account for observations. We also provide detailed predictions of atomic and molecular emission and absorption that track molecular cloud formation, with a view toward helping to verify cloud formation by shock waves. Finally, we provide an analytic solution for time-dependent H2 formation which may be of use in numerical hydrodynamic calculations.

Stars form by gravoturbulent fragmentation of interstellar gas clouds. The supersonic turbulence ubiquitously observed in Galactic molecular gas generates strong density fluctuations with gravity taking over in the densest and most massive regions. Collapse sets in to build up stars and star clusters. Turbulence plays a dual role. On global scales it provides support, while at the same time it can promote local collapse. Stellar birth is thus intimately linked to the dynamical behavior of parental gas cloud, which determines when and where protostellar cores form, and how they contract and grow in mass via accretion from the surrounding cloud material to build up stars. Slow, inefficient, isolated star formation is a hallmark of turbulent support, whereas fast, efficient, clustered star formation occurs in its absence. The fact that Galactic molecular clouds are highly filamentary can be explained by a combination of compressional flows and shear. The dynamical evolution of nascent star clusters is very complex. This strongly influences the stellar mass spectrum. The equation of state (EOS) plays a pivotal role in the fragmentation process. Under typical cloud conditions, massive stars form as part of dense clusters. However, for gas with effective polytropic index greater than unity star formation becomes biased towards isolated massive stars, which may be of relevance for understanding Pop III stars.

We discuss the lifetimes and evolution of dense cores formed as turbulent density fluctuations in magnetized, isothermal molecular clouds. We consider numerical simulations in which we measure the cores' magnetic criticality and Jeans stability in relation to the magnetic criticality of their ``parent clouds'' (the numerical boxes). In subcritical boxes, dense cores do not form, and collapse does not occur. In supercritical boxes, some cores collapse, being part of larger clumps that are supercritical from the start, and whose central, densest regions (the cores) are initially subcritical, but rapidly become supercritical, presumably by accretion along field lines. Numerical artifacts are ruled out. The time scales for cores to go from subcritical to supercritical and then collapse are a few times the core free-fall time, $\tfc$. Our results suggest that cores are out-of-equilibrium, transient structures, rather than quasi-magnetostatic configurations.

We discuss star formation in the turbulent interstellar medium. We argue that morphological appearance and dynamical evolution of the gas is primarily determined by supersonic turbulence, and that stars form via a process we call gravoturbulent fragmentation. Turbulence that is dominated by large-scale shocks or is free to decay leads to an efficient, clustered, and synchronized mode of star formation. On the other hand, when turbulence carries most of its energy on very small scales star formation is inefficient and biased towards single objects. The fact that Galactic molecular clouds are highly filamentary can be explained by a combination of compressional flows and shear. Some filaments may accumulate sufficient mass and density to become gravitationally unstable and form stars. This is observed in the Taurus molecular cloud. Timescales and spatial distribution of protostars are well explained by the linear theory of gravitational fragmentation of filaments. The dynamical evolution, especially at late times, and the final mass distribution strongly depend on the global properties of the turbulence. In dense embedded clusters mutual protostellar interactions and competition for the available mass reservoir lead to considerable stochastic variations between the mass growth histories of individual stars.

We review the main results from recent numerical simulations of turbulent fragmentation and star formation. Specifically, we discuss the observed scaling relationships, the ``quiescent'' (subsonic) nature of many star-forming cores, their energy balance, their synthesized polarized dust emission, the ages of stars associated with the molecular gas from which they have formed, the mass spectra of clumps, and the density and column density probability distribution function of the gas. We then give a critical discussion on recent attempts to explain and/or predict the star formation efficiency and the stellar initial mass function from the statistical nature of turbulent fields. Finally, it appears that turbulent fragmentation alone cannot account for the final stages of fragmentation: although the turbulent velocity field is able to produce filaments, the spatial distribution of cores in such filaments is better explained in terms of gravitational fragmentation.

Under the assumptions that molecular clouds are nearly spatially and temporally isothermal and that the density peaks (``cores'') within them are formed by turbulent fluctuations, we argue that cores cannot reach a hydrostatic (or magneto-static) state as a consequence of their formation process. In the non-magnetic case, stabilization requires a Bonnor-Ebert truncation at a finite radius, which is not feasible for a single-temperature flow, unless it amounts to a shock, which is clearly a dynamical feature. Instead, in this case, cores must be dynamical entities that can either be pushed into collapse, or else ``rebound'' towards the mean pressure and density as the parent cloud. Nevertheless, rebounding cores are delayed in their re-expansion by their own self-gravity. We give a crude estimate for the re-expansion time as a function of the closeness of the final compression state to the threshold of instability, finding typical values $\sim 1$ Myr, i.e., of the order of a few free-fall times. Our results support the notion that not all cores observed in molecular clouds need to be on route to forming stars, but that instead a class of ``failed cores'' should exist, which must eventually re-expand and disperse, and which can be identified with observed starless cores. In the magnetic case, recent observational and theoretical work suggests that all cores are critical or supercritical, and are thus qualitatively equivalent to the non-magnetic case. Our results support the notion that the entire star formation process is dynamical, with no intermediate hydrostatic stages.

(Abrigded) We discuss the column density profiles of "cores" in 3D SPH numerical simulations of turbulent molecular clouds. The SPH scheme allows us to perform a high spatial resolution analysis of the density maxima (cores) at scales between ~0.003 and 0.3pc. We analyze simulations in three different physical conditions: large scale driving, small scale driving, and random Gaussian initial conditions without driving; each one at two different timesteps: just before self-gravity is turned-on, and when gravity has been operating such that 5% of the total mass in the box has been accretted into cores. For this dataset, we perform Bonnor-Ebert fits to the column density profiles of cores. We find that 65% of the cores can be matched to Bonnor-Ebert profiles, and of these, 47% correspond to stable equilibrium configurations with xi_max < 6.5, even though the cores analyzed in the simulations are not in equilibrium, but instead are dynamically evolving. We also find in some cases substantial superposition effects when we analyze the projection of the density structures. As a consequence, different projections of the same core may give very different values of the BE fits. Finally, we briefly discuss recent results claiming that Bok globule B68 is in hydrostatic equilibrium, stressing that they imply that this core is unstable by a wide margin. We conclude that fitting BE profiles to observed cores is not an unambiguous test of hydrostatic equilibrium, and that fit-estimated parameters of the BE sphere may differ significantly from the actual values in the cores.

We compile a holistic scenario for molecular cloud (MC) evolution and control of the star formation efficiency (SFE), and present a first set of numerical tests of it. A {\it lossy} compressible cascade can generate density fluctuations and further turbulence at small scales from large-scale motions, implying that the turbulence in MCs may originate from the compressions that form them. Below a {\it sonic} scale $\ls$, turbulence cannot induce any further subfragmentation, nor be a dominant support agent against gravity. Since progressively smaller density peaks contain progressively smaller fractions of the mass, we expect the SFE to decrease with decreasing $\ls$, at least when the cloud is globally supported by turbulence. Our numerical experiments confirm this prediction. We also find that the collapsed mass fraction in the simulations always saturates below 100% efficiency. This may be due to the decreased mean density of the leftover interclump medium, which in real clouds (not confined to a box) should then be more easily dispersed, marking the ``death'' of the cloud. We identify two different functional dependences (``modes'') of the SFE on $\ls$, which roughly correspond to globally supported and unsupported cases. Globally supported runs with most of the turbulent energy at the largest scales have similar SFEs to those of unsupported runs, providing numerical evidence of the dual role of turbulence, whereby large-scale turbulent modes induce collapse at smaller scales. We tentatively suggest that these modes may correspond to the clustered and isolated modes of star formation, although here they are seen to form part of a continuum rather than being separate modes. Finally, we compare with previous proposals that the relevant parameter is the energy injection scale.

Molecular clouds (MC) are the densest and coldest component of the interstellar gas, and the sites of star formation. They are also turbulent and fractal and their masses and sizes span several orders of magnitude. It is also generally believed that they are close to Virial equilibrium (VE). Since this statement has been questioned by a number of authors, with important implications on molecular clouds' lifetimes, we will review this subject within the context of a turbulent ISM. In this framework, there is significant numerical evidence that MCs are not in VE, that there is a strong exchange of mass, momentum and energy between clouds and their surrounding medium, and that it is difficult (if not impossible) to form quasistatic cores inside MCs, suggesting that they must be transient, short-lived phenomena. Thus, their formation and disruption must be primarily dynamical, and probably not due to just a single mechanism, but rather to the combination of several processes. This picture seems consistent with recent estimates of ages of stars in the solar neighborhood.

We suggest that molecular cloud (MC) turbulence is a consequence of the very process of MC formation by collisions of larger-scale flows in the diffuse atomic gas, which generate turbulence in the accumulated gas through bending- mode instabilities. Turbulence is thus maintained for as long as the accumulation process lasts ($\sim$ several Myr). Assuming that supersonic turbulence in MCs has the double role of preventing global collapse while promoting the formation of smaller-scale structures by turbulent compression (``turbulent fragmentation''), we then note the following properties: a) Turbulent fragmentation necessarily deposits progressively smaller fractions of the total mass in regions of progressively smaller sizes, because the smaller structures are subsets of the larger ones. b) The turbulent spectrum implies that smaller scales have smaller velocity differences. Therefore, below some scale, denoted $l_{\rm eq}$, the turbulent motions become subsonic. This is an energy {\it distribution} phenomenon, not a dissipative one. On this basis, we propose that the star formation efficiency (SFE) is determined by the fraction of the total mass that is deposited in clumps with masses larger than $M_{\rm J}(l_{\rm eq})$, the Jeans mass at scale $l_{\rm eq}$, because subsonic turbulence cannot promote any further subfragmentation. In this scenario, the SFE should be a monotonically increasing function of the sonic and turbulent equality scale, $l_{\rm eq}$. We present preliminary numerical tests supporting this prediction, and thus the suggestion that (one of the) relevant parameter(s) is $l_{\rm eq}$, and compare with previous proposals that the relevant parameter is the energy injection scale.

(Abridged) We use the statistical tool known as the ``Spectral Correlation Function" [SCF] to intercompare simulations and observations of the atomic interstellar medium. The simulations considered mimic three distinct sets of physical conditions. One of them (run "ISM") is intended to represent a mixture of cool and warm atomic gas, and includes self-gravity and magnetic fields. For each simulation, H I spectral-line maps are synthesized and intercompared, both with each other, and with observations, using the SCF. We find that, when thermal broadening is large in comparison with fine-scale turbulent velocity structure, it masks sub-thermal velocity sub-structure in the synthesized spectra. The H I observations we use here for comparison are of the North Celestial Pole (NCP) Loop. None of the simulations match the NCP Loop data very well. The most realistic sets of line profiles and SCF statistics comes from artifically rescaling the velocity axis of run ISM. Without rescaling, almost all velocity structure is smeared out by thermal broadening. However, if the velocity axis is expanded by a factor of 6, the SCF distributions of run ISM an the NCP Loop match up fairly well. This means that the ratio of thermal to turbulent pressure in run ISM is much too large as it stands, and that the simulation is deficient in turbulent energy. This is a consequence of run ISM not including the effects of supernovae. We conclude that the SCF is a useful tool for understanding and fine-tuning simulations of interstellar gas, and in particular that realistic simulations of the atomic ISM need to include the effects of energetic stellar winds (e.g. supernovae) in order for the ratio of thermal-to-turbulent pressure to give spectra representative of the observed interstellar medium in our Galaxy.

We examine how well the physical properties of clumps in turbulent molecular clouds can be determined by measurements of observed clump structures. We compare simulated observations of three-dimensional numerical models of isothermal, magnetized, supersonic turbulence to the actual physical structure of the models. We determine how changing the parameters of turbulence changes the structure of the simulations. Stronger driving produces greater density fluctuations, and longer wavelength driving produces larger structures. Magnetic fields have a less pronounced effect on structure, and one that is not monotonic with field strength. Aligned structures are seen only with low-density tracers, and when the intensity of the field is large. Comparing different regions with the same tracers (or the same region with different tracers) can give information about physical conditions: different density tracers can help determine the size of the density fluctuations and thus the strength of the driving. Nevertheless, velocity superposition of multiple physical clumps can fully obscure the physical properties of those clumps, and short wavelength driving worsens this effect. Regarding Larson's relationships, we confirm previous claims that the mean density-size relationship is an artifact of the observations; and the velocity dispersion-size relationship, is consistent with observations. Regarding the mass spectra, we show that, when we look for clumps with high enough resolution, they converge to a log-normal function, rather than the power-law obtained in the literature.

We discuss the virial balance of all members of a cloud ensemble in numerical simulations of self-gravitating MHD turbulence. We first discuss the choice of reference frame for evaluating the terms entering the virial theorem (VT), concluding that the balance of each cloud should be measured in its own reference frame. We then report preliminary results suggesting that a) the clouds are far from virial equilibrium, with the ``geometric'' (time derivative) terms dominating the VT. b) The surface terms in the VT are as important as the volume ones, and tend to decrease the action of the latter. c) This implies that gravitational binding should be considered including the surface terms in the overall balance.

Observations of both star-forming regions and young, gas-free stellar associations indicate that most nearby molecular clouds form stars only over a short time span before dispersal; large-scale flows in the diffuse interstellar medium have the potential for forming clouds sufficiently rapidly, and for producing stellar populations with ages much less than the lateral crossing times of their host molecular clouds. We identify four important factors for understanding rapid star formation and short cloud lifetimes. First, much of the accumulation and dispersal of clouds near the solar circle might occur in the atomic phase; only the high-density portion of a cloud's lifecycle is spent in the molecular phase, thus helping to limit molecular cloud ``lifetimes''. Second, once a cloud achieves a high enough column density to form $\h2$ and CO, gravitational forces become larger than typical interstellar pressure forces; thus star formation can follow rapidly upon molecular gas formation and turbulent dissipation in limited areas of each cloud complex. Third, typical magnetic fields are not strong enough to prevent rapid cloud formation and gravitational collapse. Fourth, rapid dispersal of gas by newly-formed stars, and reduction of shielding by a small expansion of the cloud after the first events of star formation, might limit the length of the star formation epoch and the lifetime of a cloud in its molecular state. This picture emphasizes the importance of large-scale boundary conditions for understanding molecular cloud formation, and implies that star formation is a highly dynamic, rather than quasi-static, process.

We study the effects of projection of three-dimensional (3D) data onto the plane of the sky by means of numerical simulations of turbulence in the interstellar medium including the magnetic field, parameterized cooling and diffuse and stellar heating, self-gravity and rotation. We compare the physical-space density and velocity distributions with their representation in position-position-velocity (PPV) space (``channel maps''), noting that the latter can be interpreted in two ways: either as maps of the column density's spatial distribution (at a given line-of-sight (LOS) velocity), or as maps of the spatial distribution of a given value of the LOS velocity (weighted by density). This ambivalence appears related to the fact that the spatial and PPV representations of the data give significantly different views. First, the morphology in the channel maps more closely resembles that of the spatial distribution of the LOS velocity component than that of the density field, as measured by pixel-to-pixel correlations between images. Second, the channel maps contain more small-scale structure than 3D slices of the density and velocity fields, a fact evident both in subjective appearance and in the power spectra of the images. This effect may be due to a pseudo-random sampling (along the LOS) of the gas contributing to the structure in a channel map: the positions sampled along the LOS (chosen by their LOS velocity) may vary significantly from one position in the channel map to the next.

We suggest that molecular clouds can be formed on short time scales by compressions from large scale streams in the interstellar medium (ISM). In particular, we argue that the Taurus-Auriga complex, with filaments of 10-20 pc $\times$ 2-5 pc, most have been formed by H I flows in $\lesssim 3$Myr, explaining the absence of post-T Tauri stars in the region with ages $\gtrsim 3$ Myr. Observations in the 21 cm line of the H I `halos' around the Taurus molecular gas show many features (broad asymmetric profiles, velocity shifts of H I relative to $^{12}$CO) predicted by our MHD numerical simulations, in which large-scale H I streams collide to produce dense filamentary structures. This rapid evolution is possible because the H I flows producing and disrupting the cloud have much higher velocities (5-10 kms) than present in the molecular gas resulting from the colliding flows. The simulations suggest that such flows can occur from the global ISM turbulence without requiring a single triggering event such as a SN explosion.

We examine the idea that diffuse and giant molecular clouds and their substructure form as density fluctuations induced by large scale interstellar turbulence. We do this by investigating the topology of various fields in realistic simulations of the ISM. We find that a) the velocity field is continuous across threshold-defined cloud boundaries; b) such cloud boundaries are rather arbitrary, with no correspondence to any actual physical boundary, such as a density discontinuity; c) abrupt velocity jumps are coincident with the density maxima; d) the volume and surface kinetic terms in the Eulerian Virial Theorem for a cloud ensemble are comparable in general; e) the magnetic field exhibits bends and reversals highly correlated with similar density features. These results suggest that clouds are formed by colliding gas streams. Within this framework, we argue that thermal pressure equilibrium is irrelevant for cloud confinement in a turbulent medium, since inertial motions can still distort or disrupt a cloud. Turbulent pressure confinement appears self-defeating, because turbulence contains large-scale motions which necessarily distort cloud boundaries. Density-weighted velocity histograms show similar FWHMs and similar multi-component structure to those of observational line profiles, though the histogram features do not correspond to isolated "clumps", but rather to extended regions throughout a cloud. We argue that the results presented here may be also applicable to small scales with larger densities (molecular clouds and cores) and suggest that quasi- hydrostatic configurations cannot be produced from turbulent fluctuations unless the thermodynamic behavior of the flow becomes nearly adiabatic at late stages of collapse. We expect this to occur only at protostellar densities.

We present results from a virial analysis of fully nonlinear two-dimensional (2D) simulations of the ISM. We discuss the Eulerian Virial Theorem in 2D, and describe preliminary results on the virial budget of clouds in the simulations. The clouds are far from a static equilibrium, and the Virial Theorem is dominated by the time-derivative terms, indicating the importance of flux through the cloud boundaries and mass redistributions. A trend towards greater importance of the gravitational term at larger scales is observed, although a few small clouds are strongly self-gravitating. The magnetic and kinetic terms scale linearly with each other.

We present results from a statistical study of clouds in two-dimensional numerical simulations of the interstellar medium. The clouds in the simulations exhibit a differential mass spectrum $dN(M)/dM \sim M^{-1.44 \pm 0.1}$ and a velocity dispersion-size relation $\Delta v \sim R^{0.41 \pm 0.08}$. However, the clouds do {\it not} exhibit a clear density-size relation. At a given mean density, clouds span a range of sizes from the smallest resolved scales up to a maximum given by a Larson-type relation $R_{\rm max} \sim \rho^\alpha$, with $\alpha = -0.81 \pm .15$, although numerical effects cannot be ruled out as responsible for the latter correlation. Clouds span a range of column densities $N$ of two orders of magnitude, supporting the suggestion that the observational density-size relation may be an artifact of survey limitations. In this case, the $\Delta v$--$R$ relation can be interpreted as a direct consequence of a $k^{-2}$ turbulent spectrum, characteristic of a field of shocks, verified in the simulations, rather than of virial equilibrium of clouds with a $\rho \propto R^{-1}$ law. However, we also discuss the possibility that the clouds are in balance between self-gravity and turbulence, but with a scatter of at least a factor of 10 in the $Delta v$--$R$ relation, and of 100 in the density-size relation, according to the equilibrium relation $\Delta v \sim (NR)^{1/2}$. Finally, we compare these results with observational data. We propose a simple model suggesting that recent results finding nearly constant column densities for dark IRAS clouds may be an artifact of a temperature gradient within the clouds induced by external radiative heating. As a consequence, we emphasize that IRAS surface brightness maps are not appropriate for measuring column densities.

We present preliminary results on the energy budgets of clouds in two-dimensional numerical simulations of the interstellar medium. Using an automated cloud-identification algorithm, we calculate the gravitational, internal, kinetic and magnetic energies of the clouds. We find that, within a dispersion of roughly one order of magnitude, the gravitational energy in the clouds is balanced by the remaining energies. Furthermore, within the same dispersion, there appears to be equipartition between the kinetic and magnetic energies.