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.

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.

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 use two hydrodynamical simulations of the self-consistent evolution of molecular clouds undergoing global hierarchical collapse, one with and the other without feedback from the photoionising radiation from massive stars, to study the effect of this feedback on the structural and kinematic properties of the gas and the stellar clusters formed in the clouds. Before the formation of massive stars, the evolution of the two simulations is very similar (implying that the feedback from low mass stars does not affect the cloud-scale evolution significantly) and the star-forming region accretes faster than it can convert gas to stars, thus keeping the instantaneous measured star formation efficiency low even in the absence of significant feedback. Afterwards, the ionising feedback first destroys the filamentary supply to star-forming hubs and ultimately removes the gas from it. The ionising feedback also affects the initial kinematics and spatial distribution of the forming stars, because the gas being dispersed continues to form stars for some time and the stars it forms inherit its motion. In the non-feedback simulation, the groups remain highly compact while, in the feedback run, the dispersal of the gas causes each group to expand, and the cluster expansion thus consists of the combined expansion of the groups, rather than a coherent radial expansion away from the main hub. We do find one example of a peripheral star-forming site that appears only in the feedback run, indicating it was triggered by the expanding shell through the collect-and-collapse mechanism. However, this appears to be the exception rather than the rule. The net effect of the photoionising feedback in our simulation is to globally reduce and ultimately terminate star formation, while simultaneously smoothing the spatial distribution of the cluster stars.

We explore the structure of magnetic field lines in and around filaments in simulations of molecular clouds undergoing global, multi-scale gravitational collapse. In these simulations, filaments are not in a static equilibrium, but are long-lived flow structures that accrete gas from their environment and direct it toward clumps embedded in the filament or at the nodes at the conjunction with other filaments. In this context, the magnetic field is dragged by the collapsing gas, so its structure must reflect the flow that generates the filament. Around the filament, the gas is accreted onto it, and the magnetic lines must then be perpendicular to the filament. As the gas density increases, the gas flow changes direction, becoming almost parallel to the filament, and magnetic lines also tend to align with it. At the spine of the filament, however, magnetic lines become perpendicular again since they must connect to lines on the opposite side of the filament, resulting in "U"-shaped magnetic structures, which tend to be stretched by the longitudinal flow along the filament. Magnetic diffusive processes, however, allow the gas to continue to flow. Assuming a stationary state in which the ram pressure of the flow balances the magnetic tension, the curvature of the field lines is determined by the diffusion rate. We derive an expression relating the curvature of the field lines to the diffusive coefficient, which may be used to observationally determine the nature of the diffusive process.

Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of partially ionized plasmas, including recent advances in different astrophysical areas in which partial ionization plays a fundamental role. We outline outstanding observational and theoretical questions and discuss possible directions for future progress.

We discuss the mechanism of cluster formation in a numerical simulation of a molecular cloud (MC) undergoing global hierarchical collapse (GHC). The global nature of the collapse implies that the SFR increases over time. The hierarchical nature of the collapse consists of small-scale collapses within larger-scale ones. The large-scale collapses culminate a few Myr later than the small-scale ones and consist of filamentary flows that accrete onto massive central clumps. The small-scale collapses form clumps that are embedded in the filaments and falling onto the large-scale collapse centers. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps. Thus, the filaments feed both gaseous and stellar material to the massive central clump. This leads to the presence of a few older stars in a region where new protostars are forming, and also to a self-similar structure, in which each unit is composed of smaller-scale sub-units that approach each other and may merge. Because the older stars formed in the filaments share the infall motion of the gas onto the central clump, they tend to have larger velocities and to be distributed over larger areas than the younger stars formed in the central clump. Finally, interpreting the IMF at face-value as a probability distribution implies that massive stars only form once the {\it local} SFR is large enough to sample the IMF up to high masses. In combination with the increase of the SFR, this implies that massive stars tend to appear late in the evolution of the MC, and only in the central massive clumps. We discuss the correspondence of these features with observed properties of young stellar clusters, finding very good qualitative agreement, thus providing support to the scenario of global, hierarchical collapse of MCs, while explaining the origin of the observed cluster structure.

We discuss the mechanism of cluster formation in hierarchically collapsing molecular clouds. Recent evidence, both observational and numerical, suggests that molecular clouds (MCs) may be undergoing global, hierarchical gravitational collapse. The "hierarchical" regime consists of small-scale collapses within larger-scale ones. The latter implies that the star formation rate increases systematically during the early stages of evolution, and occurs via filamentary flows onto "hubs" of higher density, mass, and velocity dispersion, and culminates a few Myr after than the small-scale collapses have started to form stars. In turn, the small-scale collapses occur in clumps embedded in the filaments, and are themselves falling into the larger potential well of the still-ongoing large-scale collapse. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps towards the larger potential trough, so that the filaments feed both gaseous and stellar material to the hubs. This leads to the presence of older stars in a region where new protostars are still forming, to a scale-free or fractal structure of the clusters, in which each unit is composed of smaller-scale ones, and to the eventual merging of the subunits, explaining the observed structural features of open clusters.

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.

As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more detail than anywhere else in the Universe how galaxies acquire fresh gas to fuel their continuing star formation, how they exchange gas with the surrounding intergalactic medium, and turn warm, diffuse gas into molecular clouds and ultimately stars. The $\lambda$21-cm line of atomic hydrogen (HI) is an excellent tracer of these physical processes. With the SKA we will finally have the combination of surface brightness sensitivity, point source sensitivity and angular resolution to transform our understanding of the evolution of gas in the Milky Way, all the way from the halo down to the formation of individual molecular clouds.

We report on the filaments that develop self-consistently in a new numerical simulation of cloud formation by colliding flows. As in previous studies, the forming cloud begins to undergo gravitational collapse because it rapidly acquires a mass much larger than the average Jeans mass. Thus, the collapse soon becomes nearly pressureless, proceeding along its shortest dimension first. This naturally produces filaments in the cloud, and clumps within the filaments. The filaments are not in equilibrium at any time, but instead are long-lived flow features, through which the gas flows from the cloud to the clumps. The filaments are long-lived because they accrete from their environment while simultaneously accreting onto the clumps within them; they are essentially the locus where the flow changes from accreting in two dimensions to accreting in one dimension. Moreover, the clumps also exhibit a hierarchical nature: the gas in a filament flows onto a main, central clump, but other, smaller-scale clumps form along the infalling gas. Correspondingly, the velocity along the filament exhibits a hierarchy of jumps at the locations of the clumps. Two prominent filaments in the simulation have lengths ~15 pc, and masses ~600 Msun above density n ~ 10^3 cm-3 (~2x10^3 Msun at n > 50 cm-3). The density profile exhibits a central flattened core of size ~0.3 pc and an envelope that decays as r^-2.5, in reasonable agreement with observations. Accretion onto the filament reaches a maximum linear density rate of ~30 Msun Myr^-1 pc^-1.

Observations of molecular clouds reveal a complex structure, with gas and dust often arranged in filamentary rather than spherical geometries. The associations of pre- and proto- stellar cores with the filaments suggest a direct link with the process of star formation. Any study of the properties of such filaments requires a representative samples from different enviroments and so an unbiased detection method. We developed such an approach using the Hessian matrix of a surface-brightness distribution to identify filaments and determine their physical and morphological properties. After testing the method on simulated, but realistic filaments, we apply the algorithms to column-density maps computed from Herschel observations of the Galactic Plane obtained by the Hi-GAL project. We identified ~500 filaments, in the longitude range of l=216.5 to l=225.5, with lengths from ~1 pc up to ~30 pc and widths between 0.1 pc and 2.5 pc. Average column densities are between 10^20 cm^-2 and 10^22 cm^-2. Filaments include the majority of dense material with N_H_2 > 6x10^21cm^-2. We find that the pre- and proto-stellar compact sources already identified in the same region are mostly associated with filaments. However, surface densities in excess of the expected critical values for high-mass star formation are only found on the filaments, indicating that these structures are necessary to channel material into the clumps. Furthermore, we analyze the gravitational stability of filaments and discuss their relationship with star formation.

This chapter reviews the nature of turbulence in the Galactic interstellar medium (ISM) and its connections to the star formation (SF) process. The ISM is turbulent, magnetized, self-gravitating, and is subject to heating and cooling processes that control its thermodynamic behavior. The turbulence in the warm and hot ionized components of the ISM appears to be trans- or subsonic, and thus to behave nearly incompressibly. However, the neutral warm and cold components are highly compressible, as a consequence of both thermal instability in the atomic gas and of moderately-to-strongly supersonic motions in the roughly isothermal cold atomic and molecular components. Within this context, we discuss: i) the production and statistical distribution of turbulent density fluctuations in both isothermal and polytropic media; ii) the nature of the clumps produced by thermal instability, noting that, contrary to classical ideas, they in general accrete mass from their environment; iii) the density-magnetic field correlation (or lack thereof) in turbulent density fluctuations, as a consequence of the superposition of the different wave modes in the turbulent flow; iv) the evolution of the mass-to-magnetic flux ratio (MFR) in density fluctuations as they are built up by dynamic compressions; v) the formation of cold, dense clouds aided by thermal instability; vi) the expectation that star-forming molecular clouds are likely to be undergoing global gravitational contraction, rather than being near equilibrium, and vii) the regulation of the star formation rate (SFR) in such gravitationally contracting clouds by stellar feedback which, rather than keeping the clouds from collapsing, evaporates and diperses them while they collapse.

We present a numerical study of the evolution of molecular clouds, from their formation by converging flows in the warm ISM, to their destruction by the ionizing feedback of the massive stars they form. We improve with respect to our previous simulations by including a different stellar-particle formation algorithm, which allows them to have masses corresponding to single stars rather than to small clusters, and with a mass distribution following a near-Salpeter stellar IMF. We also employ a simplified radiative-transfer algorithm that allows the stellar particles to feed back on the medium at a rate that depends on their mass and the local density. Our results are as follows: a) Contrary to the results from our previous study, where all stellar particles injected energy at a rate corresponding to a star of ~ 10 Msun, the dense gas is now completely evacuated from 10-pc regions around the stars within 10-20 Myr, suggesting that this feat is accomplished essentially by the most massive stars. b) At the scale of the whole numerical simulations, the dense gas mass is reduced by up to an order of magnitude, although star formation (SF) never shuts off completely, indicating that the feedback terminates SF locally, but new SF events continue to occur elesewhere in the clouds. c) The SF efficiency (SFE) is maintained globally at the ~ 10% level, although locally, the cloud with largest degree of focusing of its accretion flow reaches SFE ~ 30%. d) The virial parameter of the clouds approaches unity before the stellar feedback begins to dominate the dynamics, becoming much larger once feedback dominates, suggesting that clouds become unbound as a consequence of the stellar feedback. e) The erosion of the filaments that feed the star-forming clumps produces chains of isolated dense blobs reminiscent of those observed in the vicinity of the dark globule B68.

We report on high-sensitivity and high-angular resolution archival Submillimeter Array (SMA) observations of the large ($\sim$15000 AU) putative circumstellar disk associated with the O-type protostar NGC 7538 S. Observations of the continuum resolve this putative circumstellar disk into five compact sources, with sizes $\sim$ 3000 AU and masses $\sim10 M_\odot$. This confirm the results of recent millimeter observations made with CARMA/BIMA towards this object. However, we find that from most of these compact sources eject collimated bipolar outflows, revealed by our silicon monoxide (SiO {\it J}=5-4) observations and confirm that these sources have a (proto)stellar nature. All outflows are perpendicular to the large and rotating dusty structure. We propose therefore that, rather than being a single massive circumstellar disk, NGC 7538 S could be instead a large and massive contracting or rotating filament that is fragmenting at scales of 0.1 to 0.01 pc to form several B-type stars, via the standard process involving outflows and disks. As in recent high spatial resolution studies of dusty filaments, our observations also suggest that thermal pressure does not seem to be sufficient to support the filament, so that either additional support needs to be invoked, or else the filament must be in the process of collapsing. An SPH numerical simulation of the formation of a molecular cloud by converging warm neutral medium flows produces contracting filaments whose dimensions and spacings between the stars forming within them, as well as their column densities, strongly resemble those observed in the filament reported here.

I discuss the role of self-gravity and radiative heating and cooling in shaping the nature of the turbulence in the interstellar medium (ISM) of our galaxy. The heating and cooling cause it to be highly compressible, and, in some regimes of density and temperature, to become thermally unstable, tending to spontaneously segregate into warm/diffuse and cold/dense phases. On the other hand, turbulence is an inherently mixing process, tending to replenish the density and temperature ranges that would be forbidden under thermal processes alone. The turbulence in the ionized ISM appears to be transonic (i.e, with Mach numbers $\Ms \sim 1$), and thus to behave essentially incompressibly. However, in the neutral medium, thermal instability causes the sound speed of the gas to fluctuate by up to factors of $\sim 30$, and thus the flow can be highly supersonic with respect to the dense/cold gas, although numerical simulations suggest that this behavior corresponds more to the ensemble of cold clumps than to the clumps' internal velocity dispersion. Finally, coherent large-scale compressions in the warm neutral medium (induced by, say, the passage of spiral arms or by supernova shock waves) can produce large, dense molecular clouds that are subject to their own self-gravity, and begin to contract gravitationally. Because they are populated by nonlinear density fluctuations, whose local free-fall times are significantly smaller than that of the whole cloud, the fluctuations terminate their collapse earlier, giving rise to a regime of hierarchical gravitational fragmentation, with small-scale collapses occurring within larger-scale ones. Thus, the "turbulence" in molecular clouds may be dominated by a gravitationally contracting component at all scales.

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 investigate the formation and evolution of giant molecular clouds (GMCs) by the collision of convergent warm neutral medium (WNM) streams in the interstellar medium, in the presence of magnetic fields and ambipolar diffusion (AD), focusing on the evolution of the star formation rate (SFR) and efficiency (SFE), as well as of the mass-to-magnetic-flux ratio (M2FR) in the forming clouds. We find that: 1) Clouds formed by supercritical inflow streams proceed directly to collapse, while clouds formed by subcritical streams first contract and then re-expand, oscillating on the scale of tens of Myr. 2) Our suite of simulations with initial magnetic field strength of 2, 3, and 4 $\mu\G$ show that only supercritical or marginal critical streams lead to reasonable star forming rates. 3) The GMC's M2FR is a generally increasing function of time, whose growth rate depends on the details of how mass is added to the GMC from the WNM. 4) The M2FR is a highly fluctuating function of position in the clouds. 5) In our simulations, the SFE approaches stationarity, because mass is added to the GMC at a similar rate at which it converts mass to stars. In such an approximately stationary regime, the SFE provides a proxy of the supercritical mass fraction in the cloud. 6) We observe the occurrence of buoyancy of the low-M2FR regions within the gravitationally-contracting GMCs, so that the latter naturally segregate into a high-density, high-M2FR "core" and a low-density, low-M2FR "envelope", without the intervention of AD. (Abridged)

I review recent numerical and analytical work on the feedback from both low- and high-mass cluster stars into their gasoeus environment. The main conclusions are that i) outflow driving appears capable of maintaing the turbulence in parsec-sized clumps and retarding their collapse from the free-fall rate, although there exist regions within molecular clouds, and even some examples of whole clouds, which are not actively forming stars, yet are just as turbulent, so that a more universal turbulence-driving mechanism is needed; ii) outflow-driven turbulence exhibits specific spectral features that can be tested observationally; iii) feedback plays an important role in reducing the star formation rate; iv) nevertheless, numerical simulations suggest that feedback cannot completely prevent a net contracting motion of clouds and clumps. Therefore, an appealing source for driving the turbulence everywhere in GMCs is the accretion from the environment, at all scales. In this case, feedback's most important role may be to prevent a fraction of the gas nearest to newly formed stars from actually reaching them, thus reducing the star formation efficiency.

I describe the scenario of molecular cloud (MC) evolution that has emerged over the past decade or so. MCs can start out as cold atomic clouds formed by compressive motions in the warm neutral medium (WNM) of galaxies. Such motions can be driven by large-scale instabilities, or by local turbulence. The compressions induce a phase transition to the cold neutral medium (CNM) to form growing cold atomic clouds, which in their early stages may constitute thin CNM sheets. Several dynamical instabilities soon destabilize a cloud, rendering it turbulent. For solar neighborhood conditions, a cloud is coincidentally expected to become molecular, magnetically supercritical, and gravitationally dominated at roughly the same column density, $N \sim 1.5 \times 10^21 \psc \approx 10 \Msun$ pc$^{-2}$. At this point, the cloud begins to contract gravitationally. However, before its global collapse is completed ($\sim 10^7$ yr later), the nonlinear density fluctuations within the cloud, which have shorter local free-fall times, collapse first and begin forming stars, a few Myr after the global contraction started. Large-scale fluctuations of lower mean densities collapse later, so the formation of massive star-forming regions is expected to occur late in the evolution of a large cloud complex, while scattered low-mass regions are expected to form earlier. Eventually, the local star formation episodes are terminated by stellar feedback, which disperses the local dense gas, although more work is necessary to clarify the details and characteristic scales of this process.

We present numerical simulations aimed at exploring the effects of varying the sub-grid physics parameters on the evolution and the properties of the galaxy formed in a low-mass dark matter halo (~7 times 10^10 Msun/h at redshift z=0). The simulations are run within a cosmological setting with a nominal resolution of 218 pc comoving and are stopped at z = 0.43. In all of our simulations, an extended old/intermediate-age stellar halo and a more compact younger stellar disk are formed. We found that a non negligible fraction of the halo stars are formed in situ in a spheroidal distribution. Changes in the sub-grid physics parameters affect significantly and in a complex way the evolution and properties of the galaxy: (i) Lower threshold densities nsf produce larger stellar effective radii Re, less peaked circular velocity curves V_c(R), and greater amounts of low-density and hot gas in the disk mid-plane; (ii) When stellar feedback is modeled by temporarily switching off radiative cooling in the star forming regions, Re increases (by a factor of ~ 2 in our particular model), the circular velocity curve becomes flatter, and a complex multi-phase gaseous disk structure develops; (iii) A more efficient local conversion of gas mass to stars, measured by a stellar particle mass distribution biased toward larger values, increases the strength of the feedback energy injection -driving outflows and inducing burstier SF histories; iv) If feedback is too strong, gas loss by galactic outflows -which are easier to produce in low-mass galaxies- interrupts SF, whose history becomes episodic. The simulations exhibit two important shortcomings: the baryon fractions are higher, and the specific SF rates are much smaller, than observationally inferred values for redshifts ~ 0.4-1.

We numerically investigate the effect of feedback from the ionizing radiation heating from massive stars on the evolution of giant molecular clouds (GMCs) and their star formation efficiency (SFE). We find that the star-forming regions within the GMCs are invariably formed by gravitational contraction. After an initial period of contraction, the collapsing clouds begin forming stars, whose feedback evaporates part of the clouds' mass, opposing the continuing accretion from the infalling gas. The competition of accretion against dense gas consumption by star formation (SF) and evaporation by the feedback, regulates the clouds' mass and energy balance, as well as their SFE. We find that, in the presence of feedback, the clouds attain levels of the SFE that are consistent at all times with observational determinations for regions of comparable SF rates (SFRs). However, we observe that the dense gas mass is larger in general in the presence of feedback, while the total (dense gas + stars) is nearly insensitive to the presence of feedback, suggesting that the total mass is determined by the accretion, while the feedback inhibits mainly the conversion of dense gas to stars. The factor by which the SFE is reduced upon the inclusion of feedback is a decreasing function of the cloud's mass, for clouds of size ~ 10 pc. This naturally explains the larger observed SFEs of massive-star forming regions. We also find that the clouds may attain a pseudo-virialized state, with a value of the virial mass very similar to the actual cloud mass. However, this state differs from true virialization in that the clouds are the center of a large-scale collapse, continuously accreting mass, rather than being equilibrium entities.

We investigate the response of the star formation efficiency (SFE) to the main parameters of simulations of molecular cloud formation by the collision of warm diffuse medium (WNM) cylindrical streams, neglecting stellar feedback and magnetic fields. The parameters we vary are the Mach number of the inflow velocity of the streams, Msinf, the rms Mach number of the initial background turbulence in the WNM, and the total mass contained in the colliding gas streams, Minf. Because the SFE is a function of time, we define two estimators for it, the "absolute" SFE, measured at t = 25 Myr into the simulation's evolution (sfeabs), and the "relative" SFE, measured 5 Myr after the onset of star formation in each simulation (sferel). The latter is close to the "star formation rate per free-fall time" for gas at n = 100 cm^-3. We find that both estimators decrease with increasing Minf, although by no more than a factor of 2 as Msinf increases from 1.25 to 3.5. Increasing levels of background turbulence similarly reduce the SFE, because the turbulence disrupts the coherence of the colliding streams, fragmenting the cloud, and producing small-scale clumps scattered through the numerical box, which have low SFEs. Finally, the SFE is very sensitive to the mass of the inflows, with sferel decreasing from ~0.4 to ~0.04 as the the virial parameter in the colliding streams increases from ~0.15 to ~1.5. This trend is in partial agreement with the prediction by Krumholz & McKee (2005), since the latter lies within the same range as the observed efficiencies, but with a significantly shallower slope. We conclude that the observed variability of the SFE is a highly sensitive function of the parameters of the cloud formation process, and may be the cause of significant scatter in observational determinations.

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.

Abridged: We study the properties of clumps formed in three-dimensional weakly magnetized magneto-hydrodynamic simulations of converging flows in the thermally bistable, warm neutral medium (WNM). We find that: (1) Similarly to the situation in the classical two-phase medium, cold, dense clumps form through dynamically-triggered thermal instability in the compressed layer between the convergent flows, and are often characterised by a sharp density jump at their boundaries though not always. (2) However, the clumps are bounded by phase-transition fronts rather than by contact discontinuities, and thus they grow in size and mass mainly by accretion of WNM material through their boundaries. (3) The clump boundaries generally consist of thin layers of thermally unstable gas, but these layers are often widened by the turbulence, and penetrate deep into the clumps. (4) The clumps are approximately in both ram and thermal pressure balance with their surroundings, a condition which causes their internal Mach numbers to be comparable to the bulk Mach number of the colliding WNM flows. (5) The clumps typically have mean temperatures 20 < T < 50 K, corresponding to the wide range of densities they contain (20 < n < 5000 pcc) under a nearly-isothermal equation of state. (6) The turbulent ram pressure fluctuations of the WNM induce density fluctuations that then serve as seeds for local gravitational collapse within the clumps. (7) The velocity and magnetic fields tend to be aligned with each other within the clumps, although both are significantly fluctuating, suggesting that the velocity tends to stretch and align the magnetic field with it. (8) The typical mean field strength in the clumps is a few times larger than that in the WNM. (9) The magnetic field strength has a mean value of B ~ 6 mu G ...

The interstellar medium (ISM) is subject, on one hand, to heating and cooling processes that tend to segregate it into distinct phases due to thermal instability (TI), and on the other, to turbulence-driving mechanisms that tend to produce strong nonlinear fluctuations in all the thermodynamic variables. In this regime, large-scale turbulent compressions in the stable warm neutral medium (WNM) dominate the clump-formation process rather than the linear developent of TI. Cold clumps formed by this mechanism are often bounded by sharp density and temperature discontinuities, which however are not contact discontinuities as in the classical 2-phase model, but rather "phase transition fronts", across which there is net mass and momentum flux from the WNM into the clumps. The clumps grow mainly by accretion through their boundaries, are in both thermal and ram pressure balance with their surroundings, and are internally turbulent as well, thus also having significant density fluctuations inside. The temperature and density of the cold and warm gas around the phase transition fronts fluctuate with time and location due to fluctuations in the turbulent pressure. Moreover, shock-compressed diffuse unstable gas can remain in the unstable regime for up to a few Myr before it undergoes a phase transition to the cold phase. These processes populate the classically forbidden density and temperature ranges. Since gas at all temperatures appears to be present in bi- or tri-stable turbulence, we conclude that the word "phase" applies only locally, surrounding phase transition sites in the gas. Globally, the word "phase" must relax its meaning to simply denote a certain temperature or density range.

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.

In this paper, we review some of the properties of dense molecular cloud cores. The results presented here rely on three-dimensional numerical simulations of isothermal, magnetized, turbulent, and self-gravitating molecular clouds (MCs) in which dense core form as a consequence of the gravo-turbulent fragmentation of the clouds. In particular we discuss issues related to the mass spectrum of the cores, their lifetimes and their virial balance.

{It has recently been proposed that giant molecular complexes form at the sites where streams of diffuse warm atomic gas collide at transonic velocities.} {We study the global statistics of molecular clouds formed by large scale colliding flows of warm neutral atomic interstellar gas under ideal MHD conditions. The flows deliver material as well as kinetic energy and trigger thermal instability leading eventually to gravitational collapse.} {We perform adaptive mesh refinement MHD simulations which, for the first time in this context, treat self-consistently cooling and self-gravity.} {The clouds formed in the simulations develop a highly inhomogeneous density and temperature structure, with cold dense filaments and clumps condensing from converging flows of warm atomic gas. In the clouds, the column density probability density distribution (PDF) peaks at $\sim 2 \times 10^{21} \psc$ and decays rapidly at higher values; the magnetic intensity correlates weakly with density from $n \sim 0.1$ to $10^4 \pcc$, and then varies roughly as $n^{1/2}$ for higher densities.} {The global statistical properties of such molecular clouds are reasonably consistent with observational determinations. Our numerical simulations suggest that molecular clouds formed by the moderately supersonic collision of warm atomic gas streams.}

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) The formation of molecular clouds (MCs) from the diffuse interstellar gas appears to be a necessary step for star formation, as young stars invariably occur within them. However, the mechanisms controlling the formation of MCs remain controversial. In this contribution, we focus on their formation in compressive flows driven by interstellar turbulence and large-scale gravitational instability. Turbulent compression driven by supernovae appears insufficient to explain the bulk of cloud and star formation. Rather, gravity must be important at all scales, driving the compressive flows that form both clouds and cores. Cooling and thermal instability allow the formation of dense gas out of moderate, transonic compressions in the warm diffuse gas, and drive turbulence into the dense clouds. MCs may be produced by an overshoot beyond the thermal-pressure equilibrium between the cold and warm phases of atomic gas, caused by some combination of the ram pressure of compression and the self-gravity of the compressed gas. In this case, properties of the clouds such as their mass, mass-to-magnetic flux ratio, and total kinetic and gravitational energies are in general time-variable quantities. MCs may never enter a quasi-equilibrium or virial equilibrium state but rather continuously collapse to stars.

We study the instantaneous virial balance of clumps and cores (CCs) in 3D simulations of driven, MHD, isothermal molecular clouds (MCs). The models represent a range of magnetic field strengths in MCs from subcritical to non-magnetic regimes. We identify CCs at different density thresholds, and for each object, we calculate all the terms that enter the Eulerian form of the virial theorem (EVT). A CC is considered gravitationally bound when the gravitational term in the EVT is larger than the amount for the system to be virialized, which is more stringent than the condition that it be large enough to make the total volume energy negative. We also calculate, quantities commonly used in the observations to indicate the state of gravitational boundedness of CCs such as the Jeans number J_c, the mass-to magnetic flux ratio mu_c, and the virial parameter alpha_vir. Our results show that: a) CCs are dynamical out-of-equilibrium structures. b) The surface energies are of the same order than their volume counterparts c) CCs are either in the process of being compressed or dispersed by the velocity field. Yet, not all CCs that have a compressive net kinetic energy are gravitationally bound. d) There is no 1-to-1 correspondence between the state of gravitational boundedness of a CC as described by the virial analysis or as implied by the classical indicators. In general, in the virial analysis, we observe that only the inner regions of the objects are gravitationally bound, whereas J_c, alpha_vir, and mu_c estimates tend to show that they are more gravitationally bound at the lowest threshold levels and more magnetically supercritical. g) We observe, in the non-magnetic simulation, the existence of a bound core with structural and dynamical properties that resemble those of the Bok globule Barnard 68 (B68).

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.

The workshop on Small Ionized and Neutral Structures in the Interstellar Medium featured many contributions on the theory of the objects which are responsible for ``Tiny Scale Atomic Structures'' (TSAS) and ``Extreme Scattering Events'' (ESE). The main demand on theory is accounting for objects that have the high densities and small sizes apparently required by the observations, but also persist over a sufficiently long time to be observable. One extensively-discussed mechanism is compressions by transonic turbulence in the warm interstellar medium, followed by thermal instabilities leading to an even more compressed state. In addressing the requirements for overpressured but persistent objects, workshop participants also discussed fundamental topics in the physics of the interstellar medium, such as the timescale for evaporation of cool dense clouds, the relevance of thermodynamically-defined phases of the ISM, the effect of magnetic fields, statistical effects, and the length and time scales introduced by interstellar processes.

We summarize recent numerical results on the control of the star formation efficiency (SFE), addressing the effects of turbulence and the magnetic field strength. In closed-box numerical simulations, the effect of the turbulent Mach number $\Ms$ depends on whether the turbulence is driven or decaying: In driven regimes, increasing $\Ms$ decreases the SFE, while in decaying regimes the converse is true. The efficiencies in non-magnetic cases for realistic Mach numbers $\Ms \sim 10$ are somewhat too high compared to observed values. Including the magnetic field can bring the SFE down to levels consistent with observations, but the intensity of the magnetic field necessary to accomplish this depends again on whether the turbulence is driven or decaying. In this kind of simulations, a lifetime of the molecular cloud (MC) needs to be assumed. Further progress requires determining the true nature of the turbulence driving and the lifetimes of the clouds. Simulations of MC formation by large-scale compressions in the warm neutral medium (WNM) show that the clouds' initial turbulence is produced by the accumulation process that forms them, and that the turbulence is driven for as long as this process lasts, producing realistic velocity dispersions and also thermal pressures in excess of the mean WNM value. In simulations including self-gravity, but neglecting the magnetic field and stellar energy feedback, the clouds never reach an equilibrium state, but rather evolve secularly, increasing their mass and gravitational energy until they engage in generalized gravitational collapse. However, local collapse events begin midways through this process, and produce enough stellar objetcs to disperse the cloud or at least halt its collapse before the latter is completed.

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 discuss molecular cloud formation by large-scale supersonic compressions in the diffuse warm neutral medium (WNM). Initially, a shocked layer forms, and within it, a thin cold layer. An analytical model and high-resolution 1D simulations predict the thermodynamic conditions in the cold layer. After $\sim 1$ Myr of evolution, the layer has column density $\sim 2.5 \times 10^{19} \psc$, thickness $\sim 0.03$ pc, temperature $\sim 25$ K and pressure $\sim 6650$ K $\pcc$. These conditions are strongly reminiscent of those recently reported by Heiles and coworkers for cold neutral medium sheets. In the 1D simulations, the inflows into the sheets produce line profiles with a central line of width $\sim 0.5 \kms$ and broad wings of width $\sim 1 \kms$. 3D numerical simulations show that the cold layer develops turbulent motions and increases its thickness, until it becomes a fully three-dimensional turbulent cloud. Fully developed turbulence arises on times ranging from $\sim 7.5$ Myr for inflow Mach number $\Mr = 2.4$ to $> 80$ Myr for $\Mr = 1.03$. These numbers should be considered upper limits. The highest-density turbulent gas (HDG, $n > 100 \pcc$) is always overpressured with respect to the mean WNM pressure by factors 1.5--4, even though we do not include self-gravity. The intermediate-density gas (IDG, $10 < n [{\rm cm}^ {-3}] < 100$) has a significant pressure scatter that increases with $\Mr$, so that at $\Mr = 2.4$, a significant fraction of the IDG is at a higher pressure than the HDG. Our results suggest that the turbulence and at least part of the excess pressure in molecular clouds can be generated by the compressive process that forms the clouds themselves, and that thin CNM sheets may be formed transiently by this mechanism, when the compressions are only weakly supersonic.

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.

We present a systematic numerical study of the effect of turbulent velocity fluctuations on the thermal pressure distribution in thermally bistable flows. The simulations employ a random turbulent driving generated in Fourier space rather than star-like heating. The turbulent fluctuations are characterized by their rms Mach number M and the energy injection wavenumber, k_for. Our results are consistent with the picture that as either of these parameters is increased, the local ratio of turbulent crossing time to cooling time decreases, causing transient structures in which the effective behavior is intermediate between the thermal-equilibrium and adiabatic regimes. As a result, the effective polytropic exponent gamma_ef ranges between ~0.2 to ~1.1. The fraction of high-density zones with P>10^4 Kcm^-3 increases from roughly 0.1% at k_for=2 and M=0.5 to roughly 70% for k_for=16 and M=1.25. A preliminary comparison with the pressure measurements of Jenkins (2004) favors our case with M=0.5 and k_for=2. In all cases, the dynamic range of the pressure summed over the entire density range, typically spans 3-4 orders of magnitude. The total pressure histogram widens as the Mach number is increased, and develops near-power-law tails at high (resp.low) pressures when gamma_ef<~ 0.5 (resp. gamma_ef>~ 1), which occurs at k_for=2 (resp.k_for=16) in our simulations. The opposite side of the pressure histogram decays rapidly, in an approx. lognormal form. Our results show that turbulent advection alone can generate large pressure scatters, with power-law high-P tails for large-scale driving, and provide validation for approaches attempting to derive the shape of the pressure histogram through a change of variable from the known form of the density histogram, such as that performed by MacLow et al.(2004).

[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.

Supersonic turbulence plays a dual role in molecular clouds: On one hand, it contributes to the global support of the clouds, while on the other it promotes the formation of small-scale density fluctuations, identifiable with clumps and cores. Within these, the local Jeans length $\Ljc$ is reduced, and collapse ensues if $\Ljc$ becomes smaller than the clump size and the magnetic support is insufficient (i.e., the core is ``magnetically supercritical''); otherwise, the clumps do not collapse and are expected to re-expand and disperse on a few free-fall times. This case may correspond to a fraction of the observed starless cores. The star formation efficiency (SFE, the fraction of the cloud's mass that ends up in collapsed objects) is smaller than unity because the mass contained in collapsing clumps is smaller than the total cloud mass. However, in non-magnetic numerical simulations with realistic Mach numbers and turbulence driving scales, the SFE is still larger than observational estimates. The presence of a magnetic field, even if magnetically supercritical, appears to further reduce the SFE, but by reducing the probability of core formation rather than by delaying the collapse of individual cores, as was formerly thought. Precise quantification of these effects as a function of global cloud parameters is still needed.

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 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.

I present an overview of the hierarchy of structures existing in the interstellar medium (ISM) and the possible mechanisms that cause the fragmentation of one level into the next, with the formation of stars as its last step. Within this framework, I then give an overview of the contributions to this session. Numerical work addresses, at the largest scales, the shaping and formation of structures in the ISM through turbulence driven by stellar energy injection, and the resulting star formation rate as a function of mean density. At the scales of molecular clouds, results comparing observational and numerical data on the density and velocity structure of turbulence-produced cores, as well as their mass spectra, are summarized, together with existing theories of core and star formation controlled by the turbulence. Observationally, an attempt to discriminate between the standard and turbulent models of star formation is presented, finding inconclusive results, but suggesting that both turbulence and the magnetic field are dynamically important in molecular clouds and their cores. Finally, various determinations of the magnetic field strength and geometry are also presented.

We summarize the current status of the turbulent model of star formation in turbulent molecular clouds. In this model, clouds, clumps and cores form a hierarchy of nested density fluctuations caused by the turbulence, and either collapse or re-expand. Cores that collapse can be either internally sub- or super-sonic. The former cannot further fragment, and can possibly be associated with the formation of a single or a few stars. The latter, instead, can undergo turbulent fragmentation during their collapse, and probably give rise to a cluster of bound objects. The star formation efficiency is low because only a small fraction of the density fluctuations proceed to collapse. Those that do not may constitute a class of ``failed'' cores that can be associated with the observed starless cores. ``Synthetic'' observations of cores in numerical simulations of non-magnetic turbulence show that a large fraction of them have subsonic internal velocity dispersions, can be fitted by Bonnor-Ebert column density profiles, and exhibit ``coherence'' (an apparent independence of linewidth with column density near the projected core centers), in agreement with observed properties of molecular cloud cores.

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.

The kinetic energy dissipation rate in the turbulent ISM of disk galaxies is a key ingredient in galaxy evolution models since it determines the effectiveness of large-scale star formation (SF) feedback. Using magneto-hydro-dynamic simulations, we find that the ISM dissipates efficiently the turbulent kinetic energy injected by sources of stellar nature. Thus, the SF process may be self-regulated by an energy balance only at the level of the disk ISM. The use of the self-regulation SF mechanism in galaxy evolutionary models, where disks form inside growing Cold Dark Matter halos, allows to predict the SF history of disk galaxies, including the Milky Way and the solar neighborhood, as well as the contribution of the whole population of disk galaxies to the cosmic SF history. The results are encouraging.

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.

We discuss magnetic pressure and density fluctuations in strongly turbulent isothermal MHD flows in "1+2/3" dimensions. We first consider "simple" nonlinear MHD waves, which allow us show that the slow and fast modes have different asymptotic dependences of the magnetic pressure B^2 vs. rho. For the slow mode, B^2 ~= c_1-c_2 rho, while for the fast mode, B^2 ~= rho^2. We also perform a perturbative analysis to investigate Alfven wave pressure, recovering previous results that B^2 ~= rho^gamma_e, with gamma_e ~= 2, 3/2 and 1/2 at respectively large, moderate and low M_a. This variety of scalings implies that a single polytropic description of magnetic pressure is not possible in general, since the relation between B^2 and rho depends on which mode dominates the density fluctuation production, which in turn depends on the angle between the magnetic field and the direction of wave propagation, and on the Alfvenic Mach number M_a. Typically, at small M_a, the slow mode dominates, and B is ANTIcorrelated with rho. At large M_a, both modes contribute to density fluctuation production, and the magnetic pressure decorrelates from density, exhibiting a large scatter, which however decreases towards higher densities. In this case, the unsystematic behavior of the magnetic pressure causes the density PDF to generally maintain the lognormal shape corresponding to non-magnetic isothermal turbulence, except when the slow mode dominates, in which case the PDF develops an excess at low densities. Our results are consistent with the low values and apparent lack of correlation between the magnetic field strength and density in surveys of the lower-density molecular gas, and also with the recorrelation apparently seen at higher densities, if M_a is relatively large there.

We review recent results on the nonlinear development of thermal instability in the context of the turbulent atomic interstellar medium (ISM). First, we pre- sent a brief summary of the linear theory, remarking that, in the atomic ISM, the wave mode is stable at small scales. Next, we revisit the growth of isolated entropy perturbations in initially unstable gas, as a function of the ratio $\eta$ of the cooling to the dynamical crossing times. Third, we consider the evolution of {\it velocity} perturbations. These correspond to the wave mode, and are stable at moderate amplitudes and small scales, as confirmed numerically. Fourth, we consider the behavior of magnetic pressure in turbulent regimes. We propose that recent findings of a poor B-rho correlation at low rho are due to the different B-rho scalings for the slow and fast modes of non- linear MHD waves. This implies that, in fully turbulent regimes, the magnetic field may not be a very efficient source of pressure, and that polytropic de- scriptions of magnetic pressure are probably not adequate. Finally, we discuss simulations of the ISM (and resolution issues) concerned with the possibility of significant amounts of gas being in the ``lukewarm'' temperature range be- tween the warm and cold stable phases. The mass fraction in this range in- creases, and the phase segregation decreases, as smaller scales are considered. We attribute this to the enhanced stability of moderate, adiabatic-like veloc- ity fluctuations with $\eta \gg 1$, to the recycling of gas from the dense to the diffuse phase by stellar energy injection, and to the magnetic field not being strongly turbulent there, possibly providing additional stability. Final- ly, we suggest that the lukewarm gas can be observationally distinguished through simultaneous determination of two of its thermodynamic variables.

(Abridged) We discuss the nonlinear development of the isobaric mode of thermal instability (TI) in the context of the atomic interstellar medium (ISM), in order to assess the ability of TI to establish a well-segregated multi-phase structure in the turbulent ISM. The key parameter is the ratio of the cooling time to the dynamical crossing time $\eta$. The isobaric density perturbations of sizes $\gtrsim 15$ pc in media with mean density $\sim 1$ cm$^{-3}$ develop inflow motions with Mach number larger than 0.5 and a shock that propagates outwards from the condensation and brings the surrounding medium out of thermal equilibrium. The time for the dynamical transient state to subside ranges from 4 to 30 Myr for initial density perturbations of 20% and sizes 3 to 75 pc. By the time the condensations have formed, a substantial fraction of the mass is still traversing the unstable range and, therefore, clouds formed by TI should be bounded by accreting gas traversing the unstable range, rather than by sharp transitions to the stable warm phase. In the presence of random velocity forcing the condensation process can be suppressed for arbitrarily long times if the forcing causes a moderate rms Mach number ($\gtrsim 0.3$) and extends to small enough scales that $\eta>1$. We suggest that these mechanisms may be at the origin of the relatively large amounts of gas mass in the unstable regime. We remark that in the (stable) warm diffuse medium, $\eta>1$ for velocity perturbations of scales up to several pc. The flow's response to them is thus nearly adiabatic, and relatively weakly compressible, consistent with observations that suggest a nearly Kolmogorov power spectrum in this medium.

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.

I review recent results from numerical simulations on the structure and dynamics of the ISM, and attempt to put together a coherent dynamical scenario. In particular, I discuss results on 1) the spatial distribution of the gas components, showing that reasonable agreement between simulations and observations exists, but noting that in most models the components are simply defined as temperature intervals, because distinct thermodynamic ``phases'' do not arise; 2) some statistical issues of the physical fields, like the dependence of the one-point statistics of the density field on the effective equation of state of the gas, the poor correlation of magnetic strength with density, the energy spectrum in weakly and highly compressible cases, and the one point statistics of the velocity field; 3) the effects of spectroscopic observation on distorting the physical structures and results from synthetic observations of the simulations, and 4) several dynamical and thermodynamical issues, such as the role of the thermal instability in forming and confining clouds, the continuous, rather than abrupt, transition between ``phases'', which in turn may be consequences of the dynamics rather than the agents controlling it, the possibility of short time scales ($\sim $ a few Myr) for molecular cloud formation, and the star-gas connection, mentioning that the models generally exhibit self-propagating star and cluster formation, while the stars may drive the medium- and small-scale gas motions, and that a ``star formation instability'' may induce chaotic behavior of the star formation rate locally. I conclude with a round-up view, and a discussion of the work needed ahead.

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.

We discuss the temperature distribution in a two-dimensional, thermally unstable numerical simulation of the warm and cold gas in the Galactic disk, including the magnetic field, self-gravity, the Coriolis force, stellar energy injection and a realistic cooling function. We find that ~50% of the turbulent gas mass has temperatures in what would be the thermally unstable range if thermal instability were to be considered alone. This appears to be a consequence of there being many other forces at play than just thermal pressure. We also point out that a bimodal temperature pdf is a simple consequence of the form of the interstellar cooling function and is not necessarily a signature of discontinuous phase transitions.

(Abridged) We discuss the probability distribution function (PDF) of column density resulting from density fields with lognormal PDFs, applicable to isothermal gas (e.g., probably molecular clouds). We suggest that a ``decorrelation length'' can be defined as the distance over which the density auto-correlation function has decayed to, for example, 10% of its zero-lag value, so that the density ``events'' along a line of sight can be assumed to be independent over distances larger than this, and the Central Limit Theorem should be applicable. However, using random realizations of lognormal fields, we show that the convergence to a Gaussian is extremely slow in the high- density tail. Thus, the column density PDF is not expected to exhibit a unique functional shape, but to transit instead from a lognormal to a Gaussian form as the ratio $\eta$ of the column length to the decorrelation length increases. Simultaneously, the PDF's variance decreases. For intermediate values of $\eta$, the column density PDF assumes a nearly exponential decay. We then discuss the density power spectrum and the expected value of $\eta$ in actual molecular clouds. Observationally, our results suggest that $\eta$ may be inferred from the shape and width of the column density PDF in optically-thin-line or extinction studies. Our results should also hold for gas with finite-extent power-law underlying density PDFs, which should be characteristic of the diffuse, non-isothermal neutral medium (temperatures ranging from a few hundred to a few thousand degrees). Finally, we note that for $\eta \gtrsim 100$, the dynamic range in column density is small ($\lesssim$ a factor of 10), but this is only an averaging effect, with no implication on the dynamic range of the underlying density distribution.

In this paper we test the results of a recent analytical study by Lazarian and Pogosyan, on the statistics of emissivity in velocity channel maps, in the case of realistic density and velocity fields obtained from numerical simulations of MHD turbulence in the interstellar medium (ISM). To compensate for the lack of well-developed inertial ranges in the simulations due to the limited resolution, we apply a procedure for modifying the spectral slopes of the fields while preserving the spatial structures. We find that the density and velocity are moderately correlated in space and prove that the analytical results by Lazarian and Pogosyan hold in the case when these fields obey the fluid conservation equations. Our results imply that the spectra of velocity and density can be safely recovered from the position-position-velocity (PPV) data cubes available through observations, and confirm that the relative contributions of the velocity and density fluctuations to those of the emissivity depend on the velocity resolution used and on the steepness of the density spectral index. Furthermore, this paper supports previous reports that an interpretation of the features in the PPV data cubes as simple density enhancements (i.e., ``clouds'') can be often erroneous, as we observe that changes in the velocity statistics substantially modify the statistics of emissivity within the velocity data cubes.

We investigate the behavior of the magnetic pressure, $b^2$, in fully turbulent MHD flows in ``1+2/3'' dimensions by means of its effect on the probability density function (PDF) of the density field. We start by reviewing our previous results for general polytropic flows, according to which the value of the polytropic exponent $\gamma$ determines the functional shape of the PDF. A lognormal density PDF appears in the isothermal ($\gamma=1$) case, but a power-law tail at either large or small densities appears for large Mach numbers when $\gamma >1$ and $\gamma < 1$, respectively. In the isothermal magnetic case, the relevant parameter is the field fluctuation amplitude, $\dbb$. A lognormal PDF still appears for small field fluctuations (generally the case for {\it large mean fields}), but a significant low-density excess appears at large fluctuation amplitudes ({\it weak mean fields}), similar to the behavior at $\gamma > 1$ of polytropic flows. We interpret these results in terms of simple nonlinear MHD waves, for which the magnetic pressure behaves linearly with the density in the case of the slow mode, and quadratically in the case of the fast wave. Finally, we discuss some implications of these results, in particular the fact that the effect of the magnetic field in modifying the PDF is strongest when the mean field is weak.

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 review recent results from numerical simulations and related models of MHD turbulence in the interstellar medium (ISM) and in molecular clouds. We discuss the implications of turbulence for the processes of cloud formation and evolution, and the determination of clouds' physical properties. Numerical simulations of the turbulent ISM to date have included magnetic fields, self-gravity, parameterized heating and cooling, modeled star formation and other turbulent inputs. The structures which form reproduce well observed velocity-size scaling properties, while predicting the non-existence of a general density-size scaling law. Criteria for the formation of gravitationally-bound structures by turbulent compression are summarized. For flows with equations of state $P\propto \rho^\gamma$, the statistics of the density field depend on the exponent $\gamma$. Numerical simulations of both forced and decaying MHD compressible turbulence have shown that the decay rate is comparable to the non-magnetic case. For virialized clouds, the turbulent decay time is shorter than the gravitational free-fall time, so wholesale cloud collapse is only prevented by ongoing turbulent inputs and/or a strong mean magnetic field. Finally, perspectives for future work in this field are briefly discussed.

We discuss HD and MHD compressible turbulence as a cloud-forming and cloud-structuring mechanism in the ISM. Results from a numerical model of the turbulent ISM at large scales suggest that the phase-like appearance of the medium, the typical values of the densities and magnetic field strengths in the intercloud medium, as well as Larson's velocity dispersion-size scaling relation in clouds may be understood as consequences of the interstellar turbulence. However, the density-size relation appears to only hold for the densest simulated clouds, there existing a large population of small, low-density clouds, which, on the other hand, are hardest to observe. We then discuss several tests and implications of a fully dynamical picture of interstellar clouds. The results imply that clouds are transient, constantly being formed, distorted and disrupted by the turbulent velocity field, with a fraction of these fluctuations undergoing gravitational collapse. Simulated line profiles and estimated cloud lifetimes are consistent with observational data. In this scenario, we suggest it is quite unlikely that quasi-hydrostatic structures on any scale can form, and that the near pressure balance between clouds and the intercloud medium is an incidental consequence of the density field driven by the turbulence and in the presence of appropriate cooling, rather than a driving or confining mechanism.

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 discuss HD and MHD compressible turbulence as a cloud-forming and cloud-structuring mechanism in the ISM. Results from a numerical model of the turbulent ISM at large scales suggest that the phase-like appearance of the medium, the typical values of the densities and magnetic field strengths in the intercloud medium, as well as the velocity dispersion-size scaling relation in clouds may be understood as consequences of the interstellar turbulence. However, the density-size relation appears to only hold for the densest clouds, suggesting that low-column density clouds, which are hardest to observe, are turbulent transients. We then explore some properties of highly compressible polytropic turbulence, in one and several dimensions, applicable to molecular cloud scales. At low values of the polytropic index $\gamma$, turbulence may induce the gravitational collapse of otherwise linearly stable clouds, except if they are magnetically subcritical. The nature of the density fluctuations in the high Mach-number limit depends on $\gamma$, and in no case resembles that resulting from Burgers turbulence. In the isothermal ($\gamma=1$) case, the dispersion of $\ln (\rho)$ scales like the turbulent Mach number. The latter case is singular with a lognormal density pdf, while power-law tails develop at high (resp. low) densities for $\gamma <1$ (resp. $\gamma >1$).

We discuss the generation and statistics of the density fluctuations in highly compressible polytropic turbulence, based on a simple model and one-dimensional numerical simulations. Observing that density structures tend to form in a hierarchical manner, we assume that density fluctuations follow a random multiplicative process. When the polytropic exponent $\gamma$ is equal to unity, the local Mach number is independent of the density, and our assumption leads us to expect that the probability density function (PDF) of the density field is a lognormal. This isothermal case is found to be singular, with a dispersion $\sigma_s^2$ which scales like the square turbulent Mach number $\tilde M^2$, where $s\equiv \ln \rho$ and $\rho$ is the fluid density. This leads to much higher fluctuations than those due to shock jump relations. Extrapolating the model to the case $\gamma \not =1$, we find that, as the Mach number becomes large, the density PDF is expected to asymptotically approach a power-law regime, at high densities when $\gamma<1$, and at low densities when $\gamma>1$. This effect can be traced back to the fact that the pressure term in the momentum equation varies exponentially with $s$, thus opposing the growth of fluctuations on one side of the PDF, while being negligible on the other side. This also causes the dispersion $\sigma_s^2$ to grow more slowly than $\tilde M^2$ when $\gamma\not=1$. In view of these results, we suggest that Burgers flow is a singular case not approached by the high-$\tilde M$ limit, with a PDF that develops power laws on both sides.

We investigate the form of the one-point probability distribution function (pdf) for the density field of the interstellar medium using numerical simulations that successively reduce the number of physical processes included. Two-dimensional simulations of self-gravitating supersonic MHD and hydrodynamic turbulence, and of decaying Burgers turbulence, produce in all cases filamentary density structures and a power-law density pdf with logarithmic slope around -1.7. This suggests that the functional form of the pdf and the general filamentary morphology are the signature of the nonlinear advection operator. These results do not support previous claims that the pdf is lognormal. A series of 1D simulations of forced supersonic polytropic turbulence is used to resolve the discrepancy. They suggest that the pdf is lognormal only for effective polytropic indices $\gamma=1$ (or nearly lognormal for $\gamma\not=1$ if the Mach number is sufficiently small), while power laws develop at high densities if $\gamma<1$. We evaluate the polytropic index for conditions relevant to the cool interstellar medium using published cooling functions and different heating sources, finding that a lognormal pdf may occur at densities between 10$^3$ and at least 10$^4$ cm$^{-3}$. Several applications are examined. First, we question a recent derivation of the IMF from the density pdf by Padoan, Nordlund & Jones because a) the pdf does not contain spatial information, and b) their derivation produces the most massive stars in the voids of the density distribution. Second, we illustrate how a distribution of ambient densities can alter the predicted form of the size distribution of expanding shells. Finally, a brief comparison is made with the density pdfs found in cosmological simulations.

We present numerical simulations of an isothermal turbulent gas undergoing gravitational collapse, aimed at testing for ``logatropic'' behavior of the form $P_t \sim \log \rho$, where $P_t$ is the ``turbulent pressure'' and $\rho$ is the density. To this end, we monitor the evolution of the turbulent velocity dispersion $\sigma$ as the density increases during the collapse. A logatropic behavior would require that $\sigma \propto \rho^{-1/2}$, a result which, however, is not verified in the simulations. Instead, the velocity dispersion increases with density, implying a polytropic behavior of $P_t$. This behavior is found both in purely hydrodynamic as well as hydromagnetic runs. For purely hydrodynamic and rapidly-collapsing magnetic cases, the velocity dispersion increases roughly as $\sigma \propto \rho^{1/2}$, implying $P_t\sim \rho^2$, where $P_t$ is the turbulent pressure. For slowly-collapsing magnetic cases the behavior is close to $\sigma \propto \rho^{1/4}$, which implies $P_t \sim \rho^{3/2}$. We thus suggest that the logatropic ``equation of state'' may represent only the statistically most probable state of an ensemble of clouds in equilibrium between self-gravity and kinetic support, but does not adequately represent the behavior of the ``turbulent pressure'' within a cloud undergoing a dynamic compression due to gravitational collapse. Finally, we discuss the importance of the underlying physical model for the clouds (in equilibrium vs. dynamic) on the results obtained.

We investigate the properties of highly compressible turbulence and its ability to produce self-gravitating structures. The compressibility is parameterized by an effective polytropic exponent gama-eff. In the limit of small gama-eff, the density jump at shocks is shown to be of the order of e^{M^2}, and the production of vorticity by the nonlinear terms appears to be negligible. In the presence of self-gravity, we suggest that turbulence can produce bound structures for gama-eff < 2(1-1/n), where 'n' is the typical dimensionality of the turbulent compressions. We show, by means of numerical simulations, that, for sufficiently small gama-eff, small-scale turbulent density fluctuations eventually collapse even though the medium is globally stable. This result is preserved in the presence of a magnetic field for supercritical mass-to-flux ratios.

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.

In this course we review the theory of incompressible homogeneous turbulence at an elementary level, and discuss the similarities and differences expected in the compressible case, relevant to the interstellar medium and molecular clouds. We stress that a general definition of turbulence applicable to the compressible case should not rely on the Kolmogorov $k^{-5/3}$ spectrum nor on an energy cascade from large to small scales. Instead, we discuss the various possibilities for the energy spectrum of compressible turbulence, which numerical simulations suggest should be $\sim k^{-2}$, and the nature of the cascades, if at all present. We then discuss issues concerning molecular clouds which are likely to be directly related to turbulence, such as cloud formation, cloud structure, and cloud support against gravity.

The Virial Thoerem is a mathematical expression obtained from the equation of motion for a fluid, which describes the energy budget of particular regions within the flow. This course reviews the basic theory leading to the Virial Theorem, discusses its applicability and limitations, and then summarizes observational results concerning the physical and statistical properties of interstellar clouds which are normally understood in terms of the Virial Theorem, in particular the so-called ``Larson's Relations''. In all cases, the standard notions as well as relevant counterpoint are presented. In particular, the difficulties arising when the medium is fully turbulent are discussed.

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 investigate the properties of highly compressible turbulence, the compressibility arising from a small effective polytropic exponent $\gamma_e$ due to cooling. In the limit of small $\gamma_e$, the density jump at shocks is shown to be of the order of $e^{M^2}$. Without self-gravity, the density structures arising in the moderately compressible case consist mostly of patches separated by shocks and behaving like waves, while in the highly compressible case clearly defined long-lived object-like clouds emerge. When the forcing in the momentum equation is purely compressible, the rotational energy decays monotonically in time, indicating that the vortex-stretching term is not efficient in transferring energy to rotational modes. This property may be at the origin of the low amount of rotation found in interstellar clouds. Vorticity production is found to rely heavily on the presence of additional terms in the equations. In the presence of self-gravity, we suggest that turbulence can produce bound structures for $\gamma_e < 2(1-1/n)$, where $n$ is the typical dimensionality of the turbulent compressions. We support this result by means of numerical simulations in which, for sufficiently small $\gamma_e$, small-scale turbulent density fluctuations eventually collapse even though the medium is globally stable. This result is preserved in the presence of a magnetic field for supercritical mass-to-flux ratios. At larger polytropic exponents, turbulence alone is not capable of producing bound structures, and collapse can only occur when the medium is globally unstable. This mechanism is a plausible candidate for the differentiation between primordial and present-day stellar-cluster formation and for the low efficiency of star formation.

We discuss the role of turbulence in cloud and star formation, as observed in numerical simulations of the interstellar medium. Turbulent compression at the interfaces of colliding gas streams is responsible for the formation of intermediate ($\simlt 100$ pc) and small clouds (a few tens of pc), although the smallest clouds can also form from fragmentation of expanding shells around stellar heating centers. The largest cloud complexes (several hundred pc) seem to form by slow, gravitational instability-driven merging of individual clouds, which can actually be described as a large-scale tendency towards homogenization of the flow due to gravity rather than cloud collisions. These mechanisms operate as well in the presence of a magnetic field and rotation, although slight variations on the compressibility and cloud morphology are present which depend on the strength and topology of the field. In summary, the role of turbulence in the life-cycle of clouds appears to be twofold: small-scale modes contribute to cloud support, while large-scale modes can both form or disrupt clouds.

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.

We present results from two-dimensional numerical simulations of a supersonic turbulent flow in the plane of the galactic disk, incorporating shear, thresholded and discrete star formation (SF), self-gravity, rotation and magnetic fields. A test of the model in the linear regime supports the results of the linear theory of Elmegreen (1991). In the fully nonlinear turbulent regime, while some results of the linear theory persist, new effects also emerge. Some exclusively nonlinear effects are: a) Even though there is no dynamo in 2D, the simulations are able to maintain or increase their net magnetic energy in the presence of a seed uniform azimuthal component. b) A well-defined power-law magnetic spectrum and an inverse magnetic cascade are observed in the simulations, indicating full MHD turbulence. Thus, magnetic field energy is generated in regions of SF and cascades up to the largest scales. c) The field has a slight but noticeable tendency to be aligned with density features. d) The magnetic field prevents HII regions from expanding freely, as in the recent results of Slavin \& Cox (1993). e) A tendency to exhibit {\it less} filamentary structures at stronger values of the uniform component of the magnetic field is present in several magnetic runs. f) For fiducial values of the parameters, the flow in general appears to be in rough equipartition between magnetic and kinetic energy. There is no clear domination of either the magnetic or the inertial forces. g) A median value of the magnetic field strength within clouds is $\sim 12\mu$G, while for the intercloud medium a value of $\sim 3\mu$G is found. Maximum contrasts of up to a factor of $\sim 10$ are observed.