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

Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and size scales. During the last decade a variety of technological and computing advances have transformed our understanding of star formation through the use of multi-wavelength observations, large scale observational surveys, and multi-physics multi-dimensional numerical simulations. Additionally, the use of synthetic observations of simulations have provided a useful tool to interpret observational data and evaluate the importance of various physical processes on different scales in star formation. Here, we review these recent advancements in both high- ($M \gtrsim 8 \, M_{\rm \odot}$) and low-mass star formation.

We study the evolution of the virial parameter, $\alpha$, and the star formation activity of a star forming region in a numerical simulation and two observational samples likely to be at different evolutionary stages, the Pipe and the G14.225 clouds. We consider a numerical simulation of turbulence in the warm atomic gas, in which clouds form by a combination of thermal instability and turbulent compressions. In this simulation we select a region with physical properties similar to those of the Pipe in both its global physical properties and the dense cores it contains, and show that this region evolves to become similar to the G14.225 cloud and its substructure in the course of $\sim$1.5 Myr. Within this region, we follow the evolution of $\alpha$ and the star formation activity of the cores over this time interval, during which they grow in mass by accretion and develop further substructure of higher density. We find that the individual cores evolve by exhibiting first a mild decrease in $\alpha$ followed by a rapid increase when star formation begins. We suggest this evolution is due to an early loss of external compressive kinetic energy followed by an increase of the gravitationally driven motions. Nevertheless, collectively, the ensamble of clumps and cores reproduces the recently observed trend of decreasing $\alpha$ for objects of higher mass. We also find that the star formation rate and star formation efficiency increase monotonically as the region evolves. Thus, we propose that the energy balance, $\alpha$, and the star formation activity determine the evolutionary stage of a star forming region, and that a star forming region like the Pipe may evolve into a G14.225-like stage over the course of a few megayears if it is allowed to accrete from its environment. This evolution is consistent with the recently proposed scenario of hierarchical gravitational collapse.

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

In previous contributions, we have presented an analytical model describing the evolution of molecular clouds (MCs) undergoing hierarchical gravitational contraction. The cloud's evolution is characterized by an initial increase in its mass, density, and star formation rate (SFR) and efficiency (SFE) as it contracts, followed by a decrease of these quantities as newly formed massive stars begin to disrupt the cloud. The main parameter of the model is the maximum mass reached by the cloud during its evolution. Thus, specifying the instantaneous mass and some other variable completely determines the cloud's evolutionary stage. We apply the model to interpret the observed scatter in SFEs of the cloud sample compiled by Lada et al.\ as an evolutionary effect so that, although clouds such as California and Orion A have similar masses, they are in very different evolutionary stages, causing their very different observed SFRs and SFEs. The model predicts that the California cloud will eventually reach a significantly larger total mass than the Orion A cloud. Next, we apply the model to derive estimated ages of the clouds since the time when approximately 25\% of their mass had become molecular. We find ages from $\sim 1.5$ to 27 Myr, with the most inactive clouds being the youngest. Further predictions of the model are that clouds with very low SFEs should have massive atomic envelopes constituting the majority of their gravitational mass, and that low-mass clouds ($M \sim 10^3$-$10^4 \, M_\odot$) end their lives with a mini-burst of star formation, reaching SFRs $\sim 300$-$500\, M_\odot$ Myr$^{-1}$. By this time, they have contracted to become compact ($\sim 1$ pc) massive star-forming clumps, in general embedded within larger GMCs.

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

We revisit the interpretation of blue-excess molecular lines from dense collapsing cores, considering recent numerical results that suggest prestellar core collapse occurs from the outside-in, and not inside-out. We thus create synthetic molecular-line observations of simulated collapsing, spherically-symmetric, density fluctuations of low initial amplitude, embedded in a uniform, globally gravitationally unstable background, without a turbulent component. The collapsing core develops a flattened, Bonnor-Ebert-like density profile, but with an outside-in radial velocity profile, where the peak infall speeds are at large radii, in lower density gas, with cloud-to-core accretion, and no hydrostatic outer envelope. Using optically thick HCO$^+$ J=1-0 and 3-2 rotational lines, we consider several "typical" beamwidths and use a simple line-fitting model to infer infall speeds from the synthetic profiles similarly to what is done in standard line modeling. We find that the inferred infall speeds are ~ 25-30% of the actual peak infall speed as the largest speeds are downweighted by the low density gas in which they occur, due to the outside-in nature of the actual radial collapse profile. Also, with the N$_2$H$^+$ $J_{F_1F}=1_{01}-0_{12}$ hyperfine line, we investigate the change in the asymmetry parameter, $\delta v\equiv(V_{thick}-V_{thin})/\Delta v_{thin}$, during the collapse, finding good agreement with observed values. Finally, the HCO$^+$ J=3-2 lines exhibit extreme $T_b$/$T_r$-ratios like those for evolved cores, for larger beams late in the collapse. Our results suggest that standard dynamical infall reproduces several observed features, but that low-mass core infall speeds are generally undervalued, often interpreted as being subsonic although the actual speeds are supersonic, due to incorrectly assuming an inside-out infall radial velocity profile with a static outer envelope.

Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale.

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

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

We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetised and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse small--sized (~10 pc) regions on timescales of less than 1 Myr. Efficient radiative cooling of the supernova remnant as well as strong compression of the surrounding gas result in non-persistent energy and momentum input from the supernovae. However, if the time between subsequent supernovae is short and they are clustered, large hot bubbles form that disperse larger regions of the parental cloud. On longer timescales, supernova feedback increases the amount of gas with moderate temperatures (T~300-3000 K). Despite its inability to disrupt molecular clouds, supernova feedback leaves a strong imprint on the star formation process. We find an overall reduction of the star formation efficiency by a factor of 2 and of the star formation rate by roughly factors of 2-4.

We investigate the Hierarchical Gravitational Fragmentation scenario through numerical simulations of the prestellar stages of the collapse of a marginally gravitationally unstable isothermal sphere immersed in a strongly gravitationally unstable, uniform background medium. The core developes a Bonnor-Ebert (BE)-like density profile, while at the time of singularity (the protostar) formation the envelope approaches a singular-isothermal-sphere (SIS)-like $r^-2$ density profile. However, these structures are never hydrostatic. In this case, the central flat region is characterized by an infall speed, while the envelope is characterized by a uniform speed. This implies that the hydrostatic SIS initial condition leading to Shu's classical inside-out solution is not expected to occur, and therefore neither should the inside-out solution. Instead, the solution collapses from the outside-in, naturally explaining the observation of extended infall velocities. The core, defined by the radius at which it merges with the background, has a time-variable mass, and evolves along the locus of the ensemble of observed prestellar cores in a plot of $M/M_{BE}$ vs. $M$, where $M$ is the core's mass and $M_{BE}$ is the critical Bonnor-Ebert mass, spanning the range from the "stable" to the "unstable" regimes, even though it is collapsing at all times. We conclude that the presence of an unstable background allows a core to evolve dynamically from the time when it first appears, even when it resembles a pressure-confined, stable BE-sphere. The core can be thought of as a ram-pressure confined BE-sphere, with an increasing mass due to the accretion from the unstable background.

Motivated by recent numerical simulations of molecular cloud (MC)evolution, in which the clouds engage in global gravitational contraction, and local collapse events culminate significantly earlier than the global collapse, we investigate the growth of density perturbations embedded in a collapsing background, to which we refer as an Inverse Hubble Flow (IHF). We use the standard procedure for the growth of perturbations in a universe that first expands (the usual Hubble Flow) and then recollapses (the IHF). We find that linear density perturbations immersed in an IHF grow faster than perturbations evolving in a static background (the standard Jeans analysis). A fundamental distinction between the two regimes is that, in the Jeans case, the time $\tau_\mathrm{nl}$ for a density fluctuation to become nonlinear increases without limit as its initial value approaches zero, while in the IHF case $\tau_\mathrm{nl} \le \tau_\mathrm{ff}$ always, where $\tau_\mathrm{ff}$ is the free-fall time of the background density. We suggest that this effect, although moderate, implies that small-scale density fluctuations embedded in globally-collapsing clouds must collapse earlier than their parent cloud, regardless of whether the initial amplitude of the fluctuations is moderate or strongly nonlinear, thus allowing the classical mechanism of Hoyle fragmentation to operate in multi-Jeans-mass MCs. More fundamentally, our results show that, contrary to the standard paradigm that fluctuations of all scales grow at the same rate in the linear regime, the hierarchical nesting of the fluctuations of different scales does affect their growth even in the linear stage.

We discuss the dependence of various properties of the star formation rate (SFR) and efficiency (SFE) in molecular clouds (MCs) on the maximum mass reached by the clouds, based on a previously-published model for MC and SFR evolution in which the clouds were assumed to be undergoing global collapse, and the SFR was controlled by ioniztion feedback. Because the model neglects various other processes, the results presented are upper limits. We find that clouds with $\Mmax \lesssim 10^4 \Msun$ end their lives with a mini-burst, at which the SFR reaches a peak of $\sim 10^4 ~\Msun \Myr^{-1}$, although its time average is only $\SFRavg \sim \hbox{ a few} \times 10^2 \Msun \Myr^{-1}$. The corresponding efficiencies are $\SFEmax \lesssim $60\%$ and $\SFEavg \lesssim $1\%$. For more massive clouds ($\Mmax \gtrsim 10^5 ~ \Msun$), the SFR first increases and then remains roughly constant for $\sim 10^7$ yr, because the clouds are influenced by the stellar feedback since earlier in their evolution. We find that $\SFRavg$ and $\SFEavg$ are well represented by the fits $\SFRavg \approx 100 (1+\Mmax/2 \times 10^5 ~ \Msun)^{2} ~ \Msun \Myr^{-1}$ and $\SFEavg \approx 0.024 (\Mmax/10^5 ~ \Msun)^{0.28}$, respectively. The massive model clouds follow the SFR-dense gas mass relation obtained by Gao \& Solomon for infrared galaxies, extrapolated down to MCs scales. Low-mass clouds fall above this relation, in agreement with recent observations. An integration of the model-predicted $\SFRavg$ over a Galactic GMC mass spectrum yields a realistic value for the Galactic SFR. Our results reinforce the suggestion that star-forming GMCs may be in global collapse, and still have low net SFRs and SFEs, due to the evaporation of most of the cloud material by the feedback from massive stars.

The cycling of material from the interstellar medium (ISM) into stars and the return of stellar ejecta into the ISM is the engine that drives the "galactic ecology" in normal spirals, a cornerstone in the formation and evolution of galaxies through cosmic time. Major observational and theoretical challenges need to be addressed in determining the processes responsible for converting the low-density ISM into dense molecular clouds, forming dense filaments and clumps, fragmenting them into stars, OB associations and bound clusters, and characterizing the feedback that limits the rate and efficiency of star formation. This formidable task can be now effectively attacked thanks to the combination of new global-scale surveys of the Milky Way Galactic Plane from infrared to radio wavelengths, offering the possibility of bridging the gap between local and extragalactic star formation studies. The Herschel, Spitzer and WISE mid to far infrared continuum surveys, complemented by analogue surveys from ground-based facilities in the millimetre and radio wavelengths, enables us to measure the Galactic distribution and physical properties of dust on all scales and in all components of the ISM from diffuse clouds to filamentary complexes and tens of thousands of dense clumps. A complementary suite of spectroscopic surveys in various atomic and molecular tracers is providing the chemical fingerprinting of dense clumps and filaments, as well as essential kinematic information to derive distances and thus transform panoramic data into a 3D representation. The latest results emerging from these Galaxy-scale surveys are reviewed. New insights into cloud formation and evolution, filaments and their relationship to channeling gas onto gravitationally-bound clumps, the properties of these clumps, density thresholds for gravitational collapse, and star and cluster formation rates are discussed.

(Abridged). We present numerical simulations of isothermal, MHD, supersonic turbulence, designed to test various hypotheses frequently assumed in star formation(SF) theories. We consider three simulations, each with a different combination of physical size, rms sonic Mach number, and Jeans parameter, but chosen as to give the same value of the virial parameter and to conform with Larson's scaling relations. As in the non-magnetic case: we find no simultaneously subsonic and super-Jeans structures in our MHD simulations. We find that the fraction of small-scale super-Jeans structures increases when self gravity is turned on, and that the production of gravitationally unstable dense cores by turbulence alone is very low. This implies that self-gravity is in general necessary not only to induce the collapse of Jeans-unstable cores, but also to form them. We find that denser regions tend to have more negative values of the velocity divergence, implying a net inwards flow towards the regions' centers. We compare the results from our simulations with the predictions from the recent SF theories by Krumholz & McKee, Padoan & Nordlund, and Hennebelle & Chabrier, using the expressions recently provided by Federrath & Klessen. We find that none of these theories reproduces the dependence of the SFEff with Ms observed in our simulations in the MHD case. The SFEff predicted by the theories ranges between half and one order of magnitude larger than what we observe in the simulations in both the HD and the MHD cases. We conclude that the type of flow used in simulations like the ones presented here and assumed in recent SF theories, may not correctly represent the flow within actual clouds, and that theories that assume it does may be missing a fundamental aspect of the flow. We suggest that a more realistic regime may be that of hierarchical gravitational collapse.

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

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

In this contribution, we test our previously published one-dimensional PDR model for deriving total hydrogen volume densities from HI column density measurements in extragalactic regions by applying it to the Taurus molecular cloud, where its predictions can be compared to available data. Also, we make the first direct detailed comparison of our model to CO(1-0) and far-infrared emission. Using an incident UV flux G0 of 4.25 ({\chi} = 5) throughout the main body of the cloud, we derive total hydrogen volume densities of \approx 430 cm-3, consistent with the extensive literature available on Taurus. The distribution of the volume densities shows a log-normal shape with a hint of a power-law shape on the high density end. We convert our volume densities to H2 column densities assuming a cloud depth of 5 parsec and compare these column densities to observed CO emission. We find a slope equivalent to a CO conversion factor relation that is on the low end of reported values for this factor in the literature (0.9 x 1020 cm-2 (K km s-1)-1), although this value is directly proportional to our assumed value of G0 as well as the cloud depth. We seem to under-predict the total hydrogen gas as compared to 100 {\mu}m dust emission, which we speculate may be caused by a higher actual G0 incident on the Taurus cloud than is generally assumed.

We investigate the collapse of non-spherical substructures, such as sheets and filaments, which are ubiquitous in molecular clouds. Such non-spherical substructures collapse homologously in their interiors but are influenced by an edge effect that causes their edges to be preferentially accelerated. We analytically compute the homologous collapse timescales of the interiors of uniform-density, self-gravitating filaments and find that the homologous collapse timescale scales linearly with the aspect ratio. The characteristic timescale for an edge driven collapse mode in a filament, however, is shown to have a square root dependence on the aspect ratio. For both filaments and circular sheets, we find that selective edge acceleration becomes more important with increasing aspect ratio. In general, we find that lower dimensional objects and objects with larger aspect ratios have longer collapse timescales. We show that estimates for star formation rates, based upon gas densities, can be overestimated by an order of magnitude if the geometry of a cloud is not taken into account.

We present an idealized, semi-empirical model for the evolution of gravitationally contracting molecular clouds (MCs) and their star formation rate (SFR) and efficiency (SFE). The model assumes that the instantaneous SFR is given by the mass above a certain density threshold divided by its free-fall time. The instantaneous number of massive stars is computed assuming a Kroupa IMF. These stars feed back on the cloud through ionizing radiation, eroding it. The main controlling parameter of the evolution turns out to be the maximum cloud mass, $\Mmax$. This allows us to compare various properties of the model clouds against their observational counterparts. A giant molecular cloud (GMC) model ($\Mmax \sim 10^5 \Msun$) adheres very well to the evolutionary scenario recently inferred by Kawamura et al. (2009) for GMCs in the Large Magellanic Cloud. A model cloud with $\Mmax \approx 2000 \Msun$ evolves in the Kennicutt-Schmidt diagram first passing through the locus of typical low- to-intermediate mass star-forming clouds, and then moving towards the locus of high-mass star-forming ones over the course of $\sim 10$ Myr. Also, the stellar age histograms for this cloud a few Myr before its destruction agree very well with those observed in the $\rho$-Oph stellar association, whose parent cloud has a similar mass, and imply that the SFR of the clouds increases with time. Our model thus agrees well with various observed properties of star-forming MCs, suggesting that the scenario of gravitationally collapsing MCs, with their SFR regulated by stellar feedback, is entirely feasible and in agreement with key observed properties of molecular clouds.

Molecular clouds often exhibit filamentary or sheet-like shapes. We compute the free-fall time ($\tff$) for finite, uniform, self-gravitating circular sheets and filamentary clouds of small but finite thickness, so that their volume density $\rho$ can still be defined. We find that, for thin sheets, the free-fall time is larger than that of a uniform sphere with the same volume density by a factor proportional to $\sqrt{A}$, where the aspect ratio $A$ is given by $A=R/h$, $R$ being the sheet's radius and $h$ is its thickness. For filamentary clouds, the aspect ratio is defined as $A=L/\calR$, where $L$ is the filament's half length and $\calR$ is its (small) radius, and the modification factor is a more complicated, although in the limit of large $A$ it again reduces to nearly $\sqrt{A}$. We propose that our result for filamentary shapes naturally explains the ubiquitous configuration of clumps fed by filaments observed in the densest structures of molecular clouds. Also, the longer free-fall times for non-spherical geometries in general may contribute towards partially alleviating the "star-formation conundrum", namely, that the star formation rate in the Galaxy appears to be proceeding in a timescale much larger than the total molecular mass in the Galaxy divided by its typical free-fall time. If molecular clouds are in general formed by thin sheets and long filaments, then their relevant free-fall time may have been systematically underestimated, possibly by factors of up to one order of magnitude.

We present the results of radiation-magnetohydrodynamic simulations of the expansion of HII regions and surrounding photodissociation regions in turbulent, magnetised, molecular clouds on scales of up to 4 parsecs, including the effects of ionising and non-ionising ultraviolet radiation and x rays from young star clusters. We find that HII region expansion reduces the disordered component of the B field, imposing a large-scale order on the field around its border. The field in the neutral gas lies along the ionisation front, while the field in the ionised gas tends to be perpendicular to this. The highest pressure compressed neutral/molecular gas is driven towards approximate equipartition between thermal/magnetic/turbulent energy densities, whereas lower pressure neutral/molecular gas divides into quiescent, magnetically dominated regions, and, on the other hand, turbulent, demagnetised regions. The ionised gas shows approximate thermal/turbulent equipartition, but with magnetic energy densities 1 to 3 orders of magnitude lower. A high velocity dispersion (approx 8 km/s) is maintained in the ionised gas throughout our simulations, despite the mean expansion velocity being significantly lower. The B field does not significantly brake the HII region expansion on the length and timescales accessible to our simulations, but it does tend to suppress the small-scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the HII region. However, the relative luminosity of ionising and non-ionising radiation has a much larger influence than the presence or absence of the B field. When the radiation field is relatively soft (as in the case of a lower mass cluster, with earliest spectral type of B0.5), then fragmentation is less vigorous and a thick, relatively smooth PDR forms. Movies available at http://youtube.com/user/divBequals0

Interferometric observations of the W33A massive star-formation region, performed with the Submillimeter Array (SMA) and the Very Large Array (VLA) at resolutions from 5 arcsec (0.1 pc) to 0.5 arcsec (0.01 pc) are presented. Our three main findings are: (1) parsec-scale, filamentary structures of cold molecular gas are detected. Two filaments at different velocities intersect in the zone where the star formation is occurring. This is consistent with triggering of the star-formation activity by the convergence of such filaments, as predicted by numerical simulations of star formation initiated by converging flows. (2) The two dusty cores (MM1 and MM2) at the intersection of the filaments are found to be at different evolutionary stages, and each of them is resolved into multiple condensations. MM1 and MM2 have markedly different temperatures, continuum spectral indices, molecular-line spectra, and masses of both stars and gas. (3) The dynamics of the "hot-core" MM1 indicates the presence of a rotating disk in its center (MM1-Main) around a faint free-free source. The stellar mass is estimated to be approximately 10 Msun. A massive molecular outflow is observed along the rotation axis of the disk.

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

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

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

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

We discuss the dissipation of turbulent kinetic energy Ek in the global ISM by means of 2-D, MHD, non-isothermal simulations in the presence of model radiative heating and cooling. We argue that dissipation in 2D is representative of that in three dimensions as long as it is dominated by shocks rather than by a turbulent cascade. Energy is injected at a few isolated sites in space, over relatively small scales, and over short time periods. This leads to the coexistence of forced and decaying regimes in the same flow. We find that the ISM-like flow dissipates its turbulent energy rapidly. In simulations with forcing, the input parameters are the radius l_f of the forcing region, the total kinetic energy e_k each source deposits into the flow, and the rate of formation of those regions, sfr_OB. The global dissipation time t_d depends mainly on l_f. In terms of measurable properties of the ISM, t_d >= Sigma_g u_rms^2/(e_k sfr_OB), where Sigma_g is the average gas surface density and u_rms is the rms velocity dispersion. For the solar neighborhood, t_d >= 1.5x10^7 yr. The global dissipation time is consistently smaller than the crossing time of the largest energy-containing scales. In decaying simulations, Ek decreases with time as t^-n, where n~0.8-0.9. This suggests a decay with distance d as Ek\propto d^{-2n/(2-n)} in the mixed forced+decaying case. If applicable to the vertical direction, our results support models of galaxy evolution in which stellar energy injection provides significant support for the gas disk thickness, but not models of galaxy formation in which this energy injection is supposed to reheat an intra-halo medium at distances of up to 10-20 times the optical galaxy size, as the dissipation occurs on distances comparable to the disk height.

We study turbulent dissipation in the ISM and explore some implications for galaxy formation and evolution using 2D MHD numerical simulations of compressible fluids. The turbulent kinetic energy E_k is injected by stellar sources formed self-consistently in the simulation. In the ISM-like fluid, regimes of both forced and decaying turbulence coexist. In the active turbulent regions (forced regime), E_k is dissipated locally and efficiently. In the decaying regime (far from input sources), E_k(t) decays ~(1+t)^-0.8. The residual turbulent motions may propagate distances of the order of the observed disk height, suggesting that turbulence may be the responsible of vertical support and star formation self-regulation at the disk level, but not at the level of the whole cosmological halo, as would be required in some models of galaxy formation.

We investigate numerically the role of thermal instability (TI) as a generator of density structures in the interstellar medium (ISM), both by itself and in the context of a globally turbulent medium. Simulations of the instability alone show that the condenstion process which forms a dense phase (``clouds'') is highly dynamical, and that the boundaries of the clouds are accretion shocks, rather than static density discontinuities. The density histograms (PDFs) of these runs exhibit either bimodal shapes or a single peak at low densities plus a slope change at high densities. Final static situations may be established, but the equilibrium is very fragile: small density fluctuations in the warm phase require large variations in the density of the cold phase, probably inducing shocks into the clouds. This result suggests that such configurations are highly unlikely. Simulations including turbulent forcing show that large- scale forcing is incapable of erasing the signature of the TI in the density PDFs, but small-scale, stellar-like forcing causes erasure of the signature of the instability. However, these simulations do not reach stationary regimes, TI driving an ever-increasing star formation rate. Simulations including magnetic fields, self-gravity and the Coriolis force show no significant difference between the PDFs of stable and unstable cases, and reach stationary regimes, suggesting that the combination of the stellar forcing and the extra effective pressure provided by the magnetic field and the Coriolis force overwhelm TI as a density-structure generator in the ISM. We emphasize that a multi-modal temperature PDF is not necessarily an indication of a multi-phase medium, which must contain clearly distinct thermal equilibrium phases.

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

I review recent results derived from numerical simulations of the turbulent interstellar medium (ISM), in particular concerning the nature and formation of turbulent clouds, methods for comparing the structure in simulations and observations, and the effects of projection of three-dimensional structures onto two dimensions. Clouds formed as turbulent density fluctuations are probably not confined by thermal pressure, but rather their morphology may be determined by the large-scale velocity field. Also, they may have shorter lifetimes than normally believed, as the large-scale turbulent modes have larger associated velocities than the clouds' internal velocity dispersions. Structural characterization algorithms have started to distinguish the best fitting simulations to a particular observation, and have opened several new questions, such as the nature of the observed line width-size relation and of the relation between the structures seen in channel maps and the true spatial distribution of the density and velocity fields. The velocity field apparently dominates the morphology seen in intensity channel maps, at least in cases when the density field exhibits power spectra steep enough. Furthermore, the selection of scattered fluid parcels along the line of sight (LOS) by their LOS-velocity inherent to the construction of spectroscopic data may introduce spurious small-scale structure in high spectral resolution channel maps.