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 c
ontrol 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.
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 a
nd 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.
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 $muG$ show that only supercritical or marginal critical streams lead to reasonable star forming rates. 3) The GMCs 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 clum
ps 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, feedbacks 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 l
arge-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 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 feedba
ck. 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.
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 non
linear 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 o
btain 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.
