Authors: Mélanie Chevance, J. M. Diederik Kruijssen, Enrique Vazquez-Semadeni, Fumitaka Nakamura, Ralf Klessen, Javier Ballesteros-Paredes, Shu-ichiro Inutsuka, Angela Adamo, Patrick Hennebelle
Categories: Article, Galaxy evolution, Interstellar medium, Molecular clouds, Star formation
Source: Space Science Reviews
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.
Keywords: Star formation, Interstellar medium, Molecular clouds, Galaxy evolution
The processes of star formation and feedback happen at cloud scales (∼100pc) within galaxies, but they also play a critical role in galaxy evolution. The properties (e.g. size, mass, surface density, temperature, pressure) of the clouds in which stars form are directly linked to the large-scale properties and structure of their host galaxies (e.g. Colombo et al. 2014; Sun et al. 2018), which therefore sets the cloud-scale conditions for star formation and the properties and structure of the interstellar medium (ISM) in which stellar feedback occurs. In turn, stellar feedback deposits energy, momentum, mass, and metals in the surrounding ISM through photoionisation, stellar winds and supernovae (SNe), therefore playing a major role in the continuous evolution of the ISM in galaxies. This cycling of matter between gas and stars, generated by star formation and feedback on the cloud scale, therefore affects global galaxy properties, which in turn again influence the star formation and feedback processes. Describing the physical processes of star formation and feedback is critical to understand galaxy evolution through cosmic time. However, the detailed physics regulating these processes and how they depend on the large scale environment (e.g. galaxy morphology and dynamics, gas and star surface densities) remain major uncertainties in simulations of galaxy formation and evolution (e.g. Scannapieco et al. 2012; Haas et al. 2013; Hopkins et al. 2013). This situation is mostly caused by a lack of observational constraints on these processes.
The baryon cycle in the Universe is regulated by the interaction between these small-scale processes of star formation and feedback acting on the scales of giant molecular clouds (GMCs), and the large-scale accretion flows acting on galactic scales (e.g. Tumlinson et al. 2017). The large scale balance between galactic-scale gas inflow and outflow rates set the global ISM properties (e.g. gas content, star formation rate SFR) and cloud properties, therefore regulating the global rate and efficiency of star formation. In turn, star formation and feedback on the cloud scale compete with the global gas inflow to set these large scale equilibria. Characterising the mass and energy flows between the different components of this multi-scale system is therefore a key step toward a comprehensive model for galaxy-scale star formation. To achieve this goal, it is necessary to characterise the evolutionary lifecycle of the transitions between gas and stars, driven by star formation and feedback.
Star formation takes place in the dense cores of GMCs.1 Measuring the assembly time from the diffuse atomic gas to the dense molecular gas and the collapse time of these molecular clouds until they start forming stars, provides constraints on which physical mechanisms drive the star formation process in galaxies. Various mechanisms are likely to trigger cloud collapse and therefore limit their lifetime (such as the gravitational collapse of the ISM, interactions with spiral arms, epicyclic perturbations or cloud-cloud collisions). These mechanisms act on different timescales (e.g. Jeffreson and Kruijssen 2018) and a comparison between observations and theoretical predictions can reveal the dominant mechanism(s) setting the cloud lifetime. Similarly, eventual support by shear or magnetic fields can be expected if long cloud lifetimes are measured. Quantitatively describing the cloud formation and evolution lifecycle as a function of the environment (i.e. galactic properties and structure) is critical to identify the relevant physical mechanisms regulating this evolution.
In a similar way, several stellar feedback processes, such as SN explosions, stellar winds, photoionisation and radiation pressure feedback, are capable of disrupting the parent molecular cloud on different timescales (e.g. Agertz et al. 2013; Dale et al. 2014; Hopkins et al. 2018; Kruijssen et al. 2019b). By doing so, they halt star formation and limit the fraction of the gas cloud effectively converted into stars, therefore regulating the integrated, cloud-scale star formation efficiency (SFE). The timescale on which this destruction happens after the onset of star formation can be observationally measured to bring strong constraints on which feedback mechanism(s) play(s) a major role in limiting the efficiency of the conversion of molecular gas to stars. Measuring this SFE on the cloud scale is key to understanding the observed difference between the long depletion time measured on galaxy scales (∼2Gyr, required for all the gas in a galaxy to be converted into stars at the current SFR; e.g. Bigiel et al. 2008; Leroy et al. 2008; Blanc et al. 2009; Bigiel et al. 2011; Schruba et al. 2011; Leroy et al. 2013), compared to the short dynamical timescale of GMCs (about 2 orders of magnitude smaller; e.g. Zuckerman and Palmer 1974). This can be explained either by an efficient, but slow (compared to the cloud-scale dynamical time) star formation process, or by rapid, but inefficient star formation, requiring multiple cycle of cloud formation and destruction in order to convert a significant fraction of gas into stars. Unveiling the physical processes of star formation and feedback is necessary to break the degeneracy between these two possible scenarios.
Finally, the physical processes that regulate star formation and feedback in galaxies as described above also drive the formation and evolution of star clusters, because the mechanisms regulating cloud assembly and collapse set the initial conditions for star cluster formation. After the onset of star formation in a GMC, the stellar feedback mechanisms induced by young stellar clusters set the timescale for gas removal, as well as the efficiency of the conversion of gas to stars, which influence the properties of the star cluster population. We refer to Adamo et al. (2020) for more details.
In order to obtain a definitive answer to which physical mechanisms drive the matter cycle in galaxies, several critical questions remain to be answered. We outline some of them below, focusing on how characterising the molecular cloud lifecycle enables us to advance our knowledge of the physics driving star formation and feedback in galaxies.
It remains unclear what processes govern cloud formation, evolution and collapse in galaxies, from the diffuse atomic gas, to the dense molecular cores in which stars form. These processes cannot be observed directly as their durations are well in excess of a human lifetime and indirect methods need to be used. Early CO observations of the nearby spiral galaxy M51 revealed the presence of a molecular cloud population between the spiral arms of the galaxy, suggesting that these structures are not transient but rather live for ∼100Myr (e.g. Scoville et al. 1979). This cloud longevity, greatly in excess of the free-fall times of such structures (approximately a few Myr), suggests the existence of support against gravitational collapse by turbulence or magnetic fields (e.g. Fleck 1980; Shu et al. 1987; Krumholz et al. 2006). However, other studies measured much shorter molecular cloud lifetimes (∼ a few tens of Myr, both in the Milky Way and in nearby galaxies; e.g. Bash et al. 1977; Leisawitz et al. 1989; Elmegreen 2000; Hartmann et al. 2001; Engargiola et al. 2003; Kawamura et al. 2009; Meidt et al. 2015), not necessarily requiring the presence of an additional support mechanism. Until recently, the great diversity of methods used to measure the molecular cloud lifetime in various environments left it undecided whether the apparent discrepancies between measurements resulted from different physical mechanisms acting in different environments, or were merely methodological. The ambiguity of definitions (between the lifetime of molecular clouds and the lifetime of actual H2 molecules), the subjective classification of clouds required by some methods, and the difficulty of extragalactic observations at sufficient resolution, have hindered progress on this question for a long time.
The physical mechanisms driving cloud dispersal are also still debated. Are clouds mainly dispersed by stellar feedback or under the effect of dynamical processes (such as galactic shear), or are both mechanisms contributing significantly to the dispersal? In the case where stellar feedback plays the major role in destroying the molecular clouds, are early feedback mechanisms (such as photoionisation or stellar winds) rapidly inhibiting star formation or are the clouds only efficiently dispersed after the first SN explosions (with a delay time of several Myr)? Also in this case, measuring the duration of the successive phases of the cloud lifecycle is a crucial step. Resolving these questions will bring important insight on how the ISM in the neighbourhood of a young stellar region is shaped by stellar feedback, and more generally, how stellar feedback affects galaxy properties.
The answers to all of the above questions are critical to determine the resulting SFE of molecular clouds and how this efficiency depends on the environment such as galactic structure and other galaxy properties (i.e. gas and stellar surface densities, rotation curve, gas pressure, metallicity). The above considerations clearly show that the field is now making steps towards a dynamical (rather than static) view of star formation and feedback in galaxies. To make the link between this small-scale cloud lifecycle and the larger scale galactic environment, the next question will then be to describe the coupling between the cloud-scale processes of star formation and feedback and the galactic baryon cycle in terms of mass flows. Looking forward, we will need to combine the answers to these questions to construct a comprehensive, galaxy-wide description of star formation.
These questions have been notoriously difficult to address, especially due to the lack of observational constraints. The latest state-of-the-art facilities, such as the Atacama Large Millimeter/submillimeter Array (ALMA), the NOrthern Extended Millimeter Array (NOEMA), and the MUSE spectrograph on the Very Large Telescope (VLT), are pushing the sensitivity and resolution limits of observations, making it now possible to collect the required data in a statistically representative sample of galaxies, on the cloud scale. These observations are key to better understand the physical processes governing star formation and feedback in galaxies, as a function of the environment.
In this review, we aim at describing the evolutionary cycle of molecular clouds. Clouds in the ISM are typically described as discrete entities, with boundaries defined based on the chemical state of the gas (atomic or molecular), density distribution, or gravitational potential (unbound or bound). However, it is not clear that these definitions are physically meaningful, due to the hierarchical structure of the ISM (Efremov and Elmegreen 1998) and the fact that chemical and density thresholds vary with galactic environment. In practice, observational studies often use the detection of CO emission as a proxy for molecular clouds (e.g. using a cloud identification algorithm such as CPROPS, Rosolowsky and Leroy 2006) or extinction measurements (e.g. Lombardi et al. 2014). The exact meaning of what a ‘molecular cloud’ represents in this context depends on the tracers and observational technique used, as well as on the galactic environment. By contrast, clouds in numerical simulations are generally defined based on physical considerations (e.g. relying on the presence of a mass overdensity or on boundedness) and direct comparison with observational studies are therefore non-trivial. In the context of the above discussion, we therefore adopt the following definition, based on an evolutionary point of view. The term ‘molecular cloud’ here refers to over-densities of molecular gas within galaxies, with a typical observed scale of a few tens of pc (we note that this value is likely environmentally dependent and may change towards high-redshift or high-pressure environments in particular), and which are not necessarily gravitationally bound. Recent observational work reveals a correspondence between molecular clouds and the units within galaxies that undergo evolutionary lifecycles independently of their neighbours (Kruijssen et al. 2019b; Chevance et al. 2020). This mirrors theoretical results suggesting that molecular clouds represent the largest size scale that can decouple from galactic dynamics, i.e. that can become self-gravitating objects, evolving predominantly due to internal physical processes (Hopkins 2012). As such, molecular clouds can be regarded as the fundamental building blocks defining how star formation proceeds in galaxies.
This review explores some of the open questions listed above, combining both observational and theoretical perspectives. We focus in particular on giving an overview of the state of the art in understanding the molecular cloud lifecycle in galaxies. We first review the observed statistical (instantaneous) properties of the molecular cloud population in galaxies in Sect. 2. We then describe the successive phases of the evolutionary lifecycle of molecular clouds, from their assembly to their destruction in Sect. 5. In Sect. 3, we investigate the processes of star formation at the scale of GMCs. In Sect. 4, we review the various feedback mechanisms susceptible of disrupting parent molecular clouds and the associated SFE. Finally, we conclude and discuss how the cycle of star formation and feedback occurring on the cloud scale participates in the multi-scale baryon cycle regulating galaxy evolution in Sect. 6.
A promising approach to fully characterise and understand GMC formation is to look at the properties of the GMC population as a function of the environment, both within the same galaxy and across a large variety of galaxies. This approach to studying GMCs has been limited to our own Galaxy until recently, because it requires both high spatial and spectral resolution, high sensitivity, and the capability to cover a significant fraction of the galaxy area. All three requirements are impossible to achieve for external galaxies with single dish telescopes and can only be fulfilled by interferometric observatories. In particular, ALMA has been the step change that has enabled us to detect GMCs in a large set of galactic environments, well beyond the limits of the Local Group and in mass ranges comparable to those observed in our own Galaxy (≳104M⊙). Moreover, thanks to ALMA, and with the aid of gravitational lenses, it is now possible to resolve GMCs in typical star-forming galaxies at redshift z∼1, offering a unique view into the onset of star formation in a significantly younger Universe (Swinbank et al. 2015; Dessauges-Zavadsky et al. 2019).
The variety of observations of GMC populations now in hand unlocks the prospect of measuring the link between galactic environment, cloud properties and structure, and star formation, well beyond the solar neighbourhood. In this section we provide a short overview of observational constraints on the resolved cold gas content of galaxies, from z∼1 down to the local Universe. Knowing the average galactic gas conditions, distributions, and reservoirs sets the basis to probe where GMCs form. We then use observable characteristics of molecular clouds, e.g. their sizes (R), masses (M), velocity dispersions (σ), virial parameters (αvir) and cloud surface densities (Σ), to offer a snapshot view of their properties and dynamical states across a variety of galactic environments.
It is well established that ∼90% of the cosmic star formation rate is traced by normal star-forming galaxies (e.g. Rodighiero et al. 2011). These systems lie on a near-linear relationship between stellar mass (M∗) and SFR, referred to as the main-sequence of galaxies (e.g. Elbaz et al. 2007; Noeske et al. 2007). The main sequence is observed to evolve with redshift, with the typical SFR of main sequence galaxies of a given M∗ increasing with z (Whitaker et al. 2014), resulting in the rapid rise of the specific SFR of galaxies from present-time to the peak of the cosmic star formation (z∼2; Madau and Dickinson 2014). This rise has been recently explained by a similar rise of the molecular gas fraction with z (e.g. Tacconi et al. 2013), i.e. galaxies at high z are more gas-rich. These high-z systems are showing ordered disc rotation already in place at z∼1–3. However, contrary to their local counterparts, they are highly turbulent with high average velocity dispersions resulting in marginally stable discs (e.g. Wisnioski et al. 2015) in terms of the Toomre (1964) Q parameter.
In optical and UV rest-frame wavelengths, galaxies at the peak of the cosmic star formation history are characterised by clumpy, irregular morphologies (e.g. Genzel et al. 2011; Elmegreen et al. 2013; Shibuya et al. 2016). Their stellar clumps have median masses of 107M⊙ (e.g. Adamo et al. 2013; Dessauges-Zavadsky et al. 2017), SFRs higher than 0.5M⊙yr−1 (e.g. Livermore et al. 2015) and sizes between 30–300 pc (e.g. Cava et al. 2018). Recent sub-millimeter observations of a z∼1 galaxy with ALMA have enabled the determination of the characteristics of GMCs in such environments (Dessauges-Zavadsky et al. 2019), as presented in Fig. 1. This shows that the bulk of star formation in z∼1–3 main sequence galaxies seems to take place in giant clumps, which look like ‘scaled-up’ versions of local star clusters or H ii regions forming in local spiral galaxies. In addition, the internal GMC velocity dispersion increases for clouds forming in gas with high external pressures and densities, i.e. starbursts (e.g. Wei et al. 2012; Leroy et al. 2018), centres of galaxies (e.g. Oka et al. 2001; Shetty et al. 2012), and high-redshift systems (e.g. Swinbank et al. 2012; Dessauges-Zavadsky et al. 2019).
Fig. 1 Larson scaling relations from Dessauges-Zavadsky et al. (2019). Molecular gas masses (left panel) and internal velocity dispersions (right panel) are shown as a function of the size for a sample of clouds identified in the Cosmic Snake galaxy at z∼1 (pink; Dessauges-Zavadsky et al. 2019), in the SDP81 star forming galaxy z∼3 (blue; Swinbank et al. 2015), in local starburst galaxies (black) and in local quiescent galaxies (grey; see full list of references in Dessauges-Zavadsky et al. 2019). The black lines show fixed molecular gas mass surface densities of 100M⊙pc−2, 1000M⊙pc−2 and 2600M⊙pc−2. High redshift and local starburst galaxies show on average GMCs with larger surface densities, achieving energy equipartition at larger internal velocity dispersions
It is interesting to note that the gas conditions of typical galaxies at z∼1–3, although rare, can also be found in the local Universe. Observations of interacting or merging galaxies in the local Universe show that these galaxies experience an increase in molecular gas fraction and achieve shorter depletion timescales than local main sequence galaxies (see e.g. Kennicutt and Evans 2012, for a review). This results in starburst phases that displace these galaxies with respect to the local main sequence of star-forming galaxies, but overlap with the main sequence at z∼1 (Genzel et al. 2010; Saintonge et al. 2012, see also Fig. 1). Similarly, the central regions of local disc galaxies, such as the Central Molecular Zone of the Milky Way, are analogues of typical star-forming galaxies at high-redshift (e.g. Kruijssen and Longmore 2013). In summary, it is clear that average GMC properties vary with the large-scale galactic environment, but it remains unclear how the change in galactic physical conditions affects the formation of GMCs. This question can only be answered by studying the physical properties of GMCs, at high resolution, over a large range of galactic environments.
Early observational works in the Milky Way suggested a common set of cloud properties described by a size–linewidth relation (σ∝R1/2), approximate virial equilibrium (αvir≡5σ2R/GM=5σ2/πGRΣ≈1, where αvir is called the virial parameter, G is the gravitational constant, and Σ the gas mass surface density), and a roughly constant surface density (e.g. Larson 1981; Solomon et al. 1987). The virial parameter is a dimensionless quantity that expresses the importance of the cloud’s kinetic energy relative to the gravitational potential energy. The above definition of αvir assumes a sphere of constant density with no surface pressure and magnetic support (e.g. Bertoldi and McKee 1992, who also derive αvir for non-spherical clouds). Irrespective of the density profile and geometry, αvir<2 for a gravitationally bound cloud and αvir=1 if the cloud is in virial balance. Observations extending to nearby galaxies show that GMCs often have virial parameters in the range αvir=1.5–3 (e.g. Sun et al. 2018, also see below). Much of the recent discussion in the field has focused on what this virial parameter means physically. Although it indicates some level of energy equipartition between the kinetic and gravitational energy and marginal boundedness (for which αvir=2) for at least some part of the GMC population, the question is whether GMCs are in virial equilibrium (for which αvir=1), collapsing (for which αvir=1–2), unbound and transient (for which αvir>2), or possibly confined by external pressure (potentially allowing all values of αvir). Throughout the discussion, two things should be kept in mind. Firstly, the relatively small differences between these numbers do not always allow observations to discriminate between the scenarios put forward due to the uncertainties associated with the measurement. Secondly, GMCs are hierarchically structured objects, within which some part may be gravitationally bound and collapsing, whereas the GMC is globally unbound and transient. Historically, these two aspects have been sources of confusion and should be carefully considered as the field moves on to resolve this discussion.
Blitz et al. (2007) pointed out that, for most GMCs, the entire clouds exhibit roughly equal potential and kinetic energies, a feature that can be interpreted as marginal gravitational binding, with the gravitational energy being about half that necessary for virial equilibrium. An alternative interpretation has been put forward by Vázquez-Semadeni et al. (2019, see also the discussions in Sect. 3 below and in Sect. 3.2 in Krause et al. 2020) that molecular clouds may not be in equilibrium, but rather regions undergoing global hierarchical collapse. The reason is that gravitational collapse has a similar energy signature (αvir∼2) as virial equilibrium (Ballesteros-Paredes et al. 2011), since the free-fall velocity is 2 times larger than the virial velocity. For generality, and given typical observational uncertainties on the measured αvir, we therefore refer to the condition αvir=1–2 as approximate ‘energy equipartition’, without referring specifically to virial equilibrium or collapse. This interpretation may apply to (at least) some part of GMCs, as evidenced by filamentary accretion flows extending up to scales of several parsecs (e.g. Schneider et al. 2010; Sugitani et al. 2011; Kirk et al. 2013; Peretto et al. 2014; Chen et al. 2019), and by systematic shifts between ^12^CO and ^13^CO lines that indicate converging cloud-scale flows, operating on time scales of ∼30Myr (Barnes et al. 2018). We note that these timescales are consistent with upper end of the GMC lifetime measurements inferred by Chevance et al. (2020, also see Sect. 5 below). Nonetheless, these dense clumps only constitute the “tip of the iceberg” of a potential cloud-scale collapse hierarchy. It is currently unclear whether the collapse of these clumps extends to the GMC at large – especially in galaxies with low gas fractions, GMCs are often observed to have high virial parameters (αvir=2–10, e.g. Sun et al. 2018; Schruba et al. 2019), suggesting that they evolve on a crossing time rather than by gravitational free-fall. Distinguishing between these cases is complicated further by the fact that the crossing and free-fall times are correlated and typically differ by a factor of <2 (Chevance et al. 2020), making it challenging to determine which of these timescales best traces the GMC lifetime.
In the past decade, it has become clear that Larson’s other two scaling relations for GMCs are likely also restricted to the Solar Neighbourhood. In particular, the near constancy of Σ in Galactic GMCs is most likely the result of considering an environment with a single gas pressure, and observations outside of the Solar Neighbourhood have revealed a wide spectrum of GMC surface densities (see e.g. Heyer et al. 2009 for Milky Way clouds and Sun et al. 2018 for extragalactic clouds). Similarly, the size-linewidth relation may simply result from the manifestation of energy equipartition at a roughly constant GMC surface density. Despite being close to energy equipartition, it is not clear what sets the velocity dispersions of GMCs. They might represent (some combination of) turbulent motion, induced by stellar feedback (e.g. Krumholz et al. 2018), shear (e.g. Meidt et al. 2018; Kruijssen et al. 2019a), or the ambient pressure (e.g. Schruba et al. 2019), or they might arise from gravitational collapse (Ballesteros-Paredes et al. 2011; Ibáñez-Mejía et al. 2016).
Irrespective of the ongoing discussion regarding the interpretation of these relations, it is common practice to characterise GMCs by considering their position in the plane spanned by Σ and σ/R1/2, as a probe of their dynamical state and internal gas pressure (e.g. Keto and Myers 1986; Heyer et al. 2009; Ballesteros-Paredes et al. 2011; Leroy et al. 2015, 2016; Vázquez-Semadeni et al. 2019). We show an example at fixed spatial scale R in Fig. 2 (Sun et al. 2018). The data have been obtained by combining the GMC populations from two Local Group galaxies (M31 and M33), 15 galaxies from the PHANGS-ALMA survey (Leroy et al., in prep.), M51 (e.g. Pety et al. 2013; Schinnerer et al. 2013), and the Antennae system (see Sun et al. 2018, for a full description of the different datasets), at a fixed spatial resolution of R=120pc. This extensive compilation of GMCs in local galaxies reveals some of their fundamental properties. First of all, we see that the GMC population spans ∼2 orders of magnitude in velocity dispersion (σ120pc), and ∼4 orders of magnitude in surface density (Σ120pc). These large variations are most likely related to the galactic environment where GMCs form. The most extreme values are found the interacting Antennae galaxies, where the gas is experiencing elevated external pressures induced by the merger (Sun et al. 2018) and is undergoing violent compression, possibly leading to gravitational collapse (e.g. Elmegreen 2018; Vázquez-Semadeni et al. 2018). The GMC population detected in main sequence galaxies at z∼1 (Dessauges-Zavadsky et al. 2019) has physical properties similar to those detected in local starbursts (e.g. the Antennae but also NGC253; see Leroy et al. 2015, 2018). In spite of these large variations, GMCs both in local galaxies and in environments with high densities and pressures show a nearly linear scaling relation between σ/R1/2 and Σ1/2, suggesting that a significant part of the GMC population is close to gravitational boundedness and energy equipartition.
Fig. 2 The top panel shows the relation between the CO line width σ and the gas surface density Σ at a common resolution of 120 pc for the discs of a sample of 15 nearby galaxies. The bottom left panel presents the mass-weighted distribution of the virial parameter αvir and the bottom right panel the distribution of turbulent pressure Pturb for the disc (circles) and centre (star symbols) of all galaxies. The spread in molecular gas dynamical state and internal turbulent pressure is clearly visible within and between galaxies, in particular when comparing normal star forming disc galaxies with a merger system such as the Antennae, or more quiescent galaxies such as M31 and M33. Figure taken from Sun et al. (2018)
GMCs in nearby galaxies show a larger scatter of αvir=1.5–3.0 (see Fig. 2). This is true for all GMC populations in the galaxies studied by Sun et al. (2018), except in M31 and M33, where clouds seem to have on average larger kinetic energy than their gravitational energy. The latter may be interpreted as an indication that GMCs forming in low-pressure environments exhibit turbulent motions that are confined by the ambient pressure (either from the gravitational potential or from an atomic gas layer) rather than by gravity (e.g. Field et al. 2011; Schruba et al. 2019), that they are dispersed by shear-driven motion (e.g. Meidt et al. 2018, 2020), or that they are assembled by external compressive motions (Camacho et al. 2016; Ballesteros-Paredes et al. 2018). These interpretations are not truly independent. For instance, compressive motions may result from large-scale turbulence in the medium or the gravitational potential of the stellar component (e.g. Li et al. 2005b). This underlines that GMCs and the large-scale ISM form a multi-scale system, where the statistical average properties of the medium on large scales couple to the properties of individual GMCs. This may result in a large-scale statistical equilibrium (e.g. Ostriker and Shetty 2011; Krumholz et al. 2018; Sun et al. 2020) even if individual objects are out of equilibrium (see Sects. 3.3–5).
From early studies of the GMC population in our Galaxy and other galaxies, it appears that the GMC mass distribution can be described by a power law function with a slope between −1.6 and −2.0 (Kennicutt and Evans 2012). However, in the past years, evidence has been provided that galaxy dynamics, gas content and distribution might be responsible to cause noticeable changes in the mass distribution of molecular clouds (Hughes et al. 2013a; Colombo et al. 2014; Kruijssen 2014; Hughes et al. 2016). Both Hughes et al. (2013b) and Colombo et al. (2014) show that the GMCs forming in distinct dynamical regions of M51 (e.g. arm or inter-arm, upstream or downstream, disc or molecular ring) have different spectral mass distributions (with different measured normalisation and slope of the power law) and different maximum mass. The GMC mass spectrum also changes from galaxy to galaxy. Hughes et al. (2016) show that the number of clouds detected above a certain mass increases and the upper end of the mass distribution flattens towards increasing galaxy stellar mass and SFR surface density. Finally, Reina-Campos and Kruijssen (2017) find that the upper end of the GMC mass spectrum varies in a way that is expected for the competition between gravitational collapse and dispersal by centrifugal forces and stellar feedback.
It is more challenging to extend the study of the GMC mass function to high-redshift main sequence systems. So far, GMC-like objects have been detected only in a handful of galaxies. Dessauges-Zavadsky et al. (2019) show that GMCs in the Cosmic Snake main sequence galaxy at z∼1 have masses comparable to those detected in local starbursts such as the centre of NGC253 Leroy et al. (2015) and the Antennae (Whitmore et al. 2014). While the number of GMCs is still too small to allow for a robust study of their mass distributions, it has recently been reported that stellar clumps detected at z∼1 have mass distributions compatible with a power-law function with a slope of −2 (Dessauges-Zavadsky and Adamo 2018), similarly to GMCs, clusters and H ii regions forming in local galaxies. Therefore, we expect that, in spite of the extreme conditions of GMCs at high-redshift, their formation is driven by similar physical mechanisms as in local main sequence galaxies.
As reported above, GMC mass functions depend on the large-scale environment in galaxies. It is therefore crucial to combine a sample of observations in a large range of environments in order to be able to understand the origin of these variations. However, studying the mass function of molecular clouds in the Milky Way is difficult, since many of the GMCs reside in the mid-plane of the Galactic disc and overlap along the line of sight. In addition, they probe mostly a single environment. These problems are greatly reduced by looking at nearby face-on galaxies (e.g. Colombo et al. 2014; Hughes et al. 2016; Sun et al. 2018). Based on these observations, theories of molecular cloud formation and evolution should explain the origin of the mass function of GMCs and its dependence with environment. We describe below two different approaches to understand the formation of molecular clouds.
Some studies discuss the possibility of direct formation of GMCs in a top-down manner, via gravitational instabilities of the Galactic disc (e.g. Dobbs et al. 2014). This process of self-gravitating fragmentation is more likely to happen at high gas surface densities environments such as in spiral arms, where it has been studied with linear theory (e.g. Elmegreen 1979; Tubbs 1980; Marochnik et al. 1983; Balbus and Cowie 1985; Balbus 1988). In this scenario, the typical mass of resulting clouds corresponds to the characteristic mass of the gravitational instabilities and tends to be very large (up to several 106M⊙; Wada et al. 2000; Shetty and Ostriker 2008; Tasker and Tan 2009; Dobbs et al. 2011; Hopkins 2012; Kruijssen 2014). However, in our Galaxy, observed clouds are dominated in number by clouds with masses lower than 105M⊙. In this top-down scenario, low-mass clouds may arise when stellar feedback from new-born stars disperses the gravitationally unstable gas reservoir before its collapse has been completed (Reina-Campos and Kruijssen 2017).
In contrast to the above scenario, the “bottom-up” scenario for GMC assembly explains the formation of molecular clouds of various masses starting from the formation of clouds with very small masses (∼100M⊙). The typical density of molecular gas in the Galaxy is on the order of 102cm−3 or larger, i.e. at least two orders of magnitude larger than the average density of warm neutral medium that occupies most of the volume of the Galactic thin disc. In general, the formation of molecular clouds requires a phase transition from the warm neutral medium to the cold neutral medium (e.g. Inutsuka et al. 2005; Glover and Mac Low 2007; Krumholz 2014; Hennebelle and Inutsuka 2019). This phase transition inevitably generates a turbulent velocity that is supersonic relative to the sound speed of the cold neutral medium, yet subsonic relative to the sound speed of warm neutral medium, and hence remains in the system without decay due to shock dissipation (Kritsuk and Norman 2002; Koyama and Inutsuka 2002). This picture is studied in detail by many authors (e.g. Audit and Hennebelle 2005; Vázquez-Semadeni et al. 2006; Hennebelle and Audit 2007; Heitsch et al. 2006). The effect of magnetic fields on this phase transition has been studied by e.g. Inoue and Inutsuka (2008, 2009), Körtgen and Banerjee (2015), Valdivia et al. (2016), van Loo et al. (2007), Mandal et al. (2020).
In general, a single compression of the warm neutral medium by a shock wave can create a dense molecular cloud if the compression lasts for 10–20 Myr, and is oriented (at intermediate angles) along the magnetic field lines (e.g. Hennebelle and Pérault 2000; Vázquez-Semadeni et al. 2011; Fogerty et al. 2016, 2017). This approximate alignment is not unlikely, because the mean magnetic field and the mean motions of the gas in disc galaxies are nearly circular, dominated by the total (stellar and gas) galactic potential (Li et al. 2005a). If the alignment is substantially different (e.g. due to motions resulting from stellar feedback), multiple compressions of the ISM would be needed, and the timescales for molecular cloud formation could be substantially larger (Inutsuka et al. 2015).
The distribution of GMCs as a function of mass spans many orders of magnitudes and star formation occurs in molecular clouds of various masses (although most of the star formation in a galaxy may take place in the most massive GMCs, see e.g. Murray 2011). Thus, the understanding of the overall star formation rates in a region of a galaxy or galaxies as a whole requires the determination of not only the star formation rate in an individual cloud, but also the mass distribution of molecular clouds. While it is technically challenging to accurately determine the mass function of molecular clouds in our Galaxy (due to line-of-sight confusion, limited knowledge of the distances to the clouds, and limited resolution), extragalactic observations enabled by development of recent sub-millimeter observatories, such as ALMA and NOEMA, have revealed the environmental dependency of the cloud mass function (Colombo et al. 2014; Hughes et al. 2016).
Theoretical studies describing cloud properties are needed for a better understanding of the formation and destruction of molecular clouds, as a function of the environment. Until recently, it has been challenging to perform direct numerical simulations of a statistically representative ensemble of molecular clouds, with sufficient resolution to study the small scale physics such as cloud formation and destruction. This has motivated the formulation of simplified, semi-analytical models for the evolution of the GMC mass function, which can be a useful way to understand the origins of the basic properties and environmental variations of GMCs. Earlier attempts in this direction can be found in Kwan (1979), Scoville and Hersh (1979), and Tomisaka (1986) that formulated the so-called coagulation equation for molecular clouds. In these investigations, the growth of clouds is driven by cloud-cloud collision and omit any smooth accretion of molecular clouds. The recent theoretical finding of the long timescale for molecular cloud formation (Inoue and Inutsuka 2009) and the importance of the gradual growth process by accretion of dense H i gas (Inoue and Inutsuka 2012) gives a crucial need for a self-growth term in the coagulation equation (Kobayashi et al. 2017).
According to this recent development, we can adopt coarse graining of short-timescale (∼ a few Myr) events of the growth and destruction of clouds, and describe the long timescale evolution by the continuity equation of the differential number density N of molecular cloud of mass M (Kobayashi et al. 2017, 2018).
where N(dM/dt) denotes the flux of mass function in mass space, td is the cloud disruption timescale, dM/dt describes the growth rate of the molecular cloud, and the last term accounts for the growth due to cloud-cloud collisions.
If the contribution from cloud-cloud collisions is negligible (Kobayashi et al. 2017, 2018), then the mass growth can be approximated by dM/dt=M/tf with the growth timescale tf. In this case, a steady state solution to the above equation is N(M)=M−α, where the slope of the GMC mass function can then be expressed as α=1+tf/td (Inutsuka et al. 2015). In typical conditions of spiral arm regions, we expect the timescale for a massive star to form once the cloud is created to be similar to its formation timescale, i.e. t∗∼tf. The timescale for cloud dispersal by feedback after the cloud has been created must be td=t∗+tfb, where tfb is the ‘feedback timescale’ for cloud destruction after the massive stars have formed. As a result, we have tf≲td, which corresponds to 1<α≲2. For example, if tf=10Myr, then α≈1.7, which agrees with observations (Solomon et al. 1987; Kramer et al. 1998; Heyer et al. 2001; Roman-Duval et al. 2010). However, in quiescent regions with a very limited amount of gaseous material, away from spiral arms in the Galactic disc, tf is expected to be large, breaking the above assumption that t∗∼tf and instead giving t∗<tf. In this case, we can have tf>td and hence α>2. Such steep mass functions are observed in M33 (Gratier et al. 2012) and in M51 (Colombo et al. 2014). The prediction of the above model for these galaxies is that tf>td.
In order to test this prediction, it is necessary to directly measure the GMC lifetime and the feedback timescale and relate these to tf and td. In order to do so, the key question is when GMCs are visible in CO observations, because CO is the most commonly used molecule for tracing molecular gas. In the context of the above discussion, the limiting cases are either that CO is always visible, during the formation and dispersal phases, such that tCO=tf+td, or that it is only visible during the dispersal phase, such that tCO=td (this is the case preferred by Kobayashi et al. 2017). For the sake of this example, we assume that t∗∼tf. If instead tf>t∗, any slope α>2 is possible. The above cases now enable the slope of the GMC mass function predicted by equation (1) to be expressed in terms of the cloud lifetime tCO and the feedback timescale tfb. If we define x≡tfb/tCO, the prediction that α=1+tf/td implies
For the 11 galaxies considered in the discussion of the GMC lifecycle in Sect. 5, the observed range of GMC lifetimes and feedback timescales implies x=0.10–0.27 (Kruijssen et al. 2019b; Chevance et al. 2020; Hygate et al. 2019). This means that the Kobayashi et al. (2017) model predicts GMC mass function slopes of α=1.7–1.9 (if tCO=td) or α=1.6–1.8 (if tCO=tf+td). It is clear that an unambiguous prediction by this model requires a clear definition of the CO-bright phase in the context of the model, as well as providing a quantitative prediction for the GMC formation timescale tf. For the time being, the firm prediction of the model is that GMC mass functions with α<1.6 are excluded in these galaxies.
In this section, we describe the physical conditions under which clouds collapse and form stars, as well as the different possible mechanisms at the origin of the collapse. We then discuss the low observed efficiency of the conversion of gas to stars in galaxies and link this to models of rapid and slow star formation.
All present star formation occurs in the densest regions of molecular clouds. A natural question is thus to determine the characteristic physical conditions under which these regions within molecular clouds collapse. Here we assume that these conditions are set by the competition between self-gravity and thermal pressure, allowing us to define the characteristic Jeans length λJ, and the corresponding spherical2 Jeans mass MJ, above which an isothermal parcel of fluid collapses due to its self-gravity
and
where cs is the sound speed, G the gravitational constant, ρ the density of the gas, and n the gas number density.
The minimum duration of the collapse for such an object under the influence of self-gravity only is given by the free-fall
In reality, the collapse takes longer by a factor of a few, because the thermal pressure gradient delays the collapse, especially at its early stages (e.g. Larson 1969; Galván-Madrid et al. 2007).
At any level of the hierarchy in the ISM, from GMCs with masses of 105–106M⊙ and regions therein with mean densities of n∼102cm−3, down to dense cores with masses of a few M⊙ and densities of n∼105cm−3, the identified structures typically contain many Jeans masses (see e.g. Krause et al. 2020, and references therein). Nonetheless, the simple set of equations provides a timescale estimation for collapse, and thus for star formation within these collapsing regions. If not delayed by other physical mechanisms (see Sect. 3.3), collapse should occur on timescales ranging from tff∼3Myr for the regions within GMCs that have volume-averaged densities of n∼102cm−3, down to tff∼0.1Myr for dense cores of n∼105cm−3. Equations (3) and (4) show that, as collapse proceeds, the density increases and thus the Jeans length and the Jeans mass decrease, inducing fragmentation of the collapsing cloud (Hoyle 1953). Such fragmentation should stop once the gas ceases to behave isothermally. This occurs at large volume densities, when the gas becomes optically thick.
During the early discussions in the literature on Jeans fragmentation, it was argued that fragmentation should not occur, because the largest scales of a homogeneous, gravitationally unstable isothermal medium have the largest growth rates (Tohline 1980). However, during the cloud assembly different non-linear instabilities can develop, producing inhomogeneous clouds (for a discussion of the mechanisms for cloud formation, see e.g. Ballesteros-Paredes et al. 2020, in this volume). The growth timescales for such non-linear density fluctuations are substantially shorter than the timescale for the collapse of the whole cloud, as shown in numerical simulations (see, e.g. Koyama and Inutsuka 2002; Heitsch et al. 2005, 2008a,b; Audit and Hennebelle 2005, 2010; Vázquez-Semadeni et al. 2007, 2010). Nonetheless, Clark and Bonnell (2005) show that turbulent density fluctuations often do not reach high enough densities to become Jeans unstable on their own. Therefore, Vázquez-Semadeni et al. (2019) suggest that the smaller scales initiate their collapse later than the larger scales (see also Krause et al. 2020, in this volume). In this scenario, local collapse starts when the global Jeans mass of the cloud has decreased sufficiently to match the masses of local density fluctuations. At this point, the turbulent density fluctuations become unstable and begin to collapse. Smaller mass scales at a given density therefore start their collapse at a later time than the larger scales, but terminate it earlier, because their free-fall time is significantly shorter than that of the entire cloud.
If all the molecular clouds observed in our Galaxy are self-gravitating and collapse within a free-fall time, we can roughly estimate the free-fall rate of star formation (SFRff) in the Milky Way. The total molecular gas mass in our Galaxy is derived from ^12^CO observations to be about 109M⊙ (Bolatto et al. 2013; Heyer and Dame 2015), and the typical free-fall time of molecular clouds is evaluated to be ∼10Myr. If all the gas is converted to stars within a free-fall time, the free-fall rate of star formation of the Milky Way is SFRff∼100M⊙yr−1. However, the observed rate of star formation, both in the Milky Way and in nearby spiral galaxies are generally estimated to be much lower, about a factor 100 times smaller than the typical free-fall rate (e.g. McKee and Williams 1997; Robitaille and Whitney 2010; Leroy et al. 2012). In other words, the galaxy-wide SFE per free-fall time ϵff∼0.01. The immediate question at hand is therefore to determine whether this discrepancy is due to the fact that not all gas is converted into stars (i.e. the integrated SFE is low), or if the timescale of star formation is much longer than the free-fall time.
We define the instantaneous efficiency of star formation at any moment in time, ϵ(t), by the amount of mass that has gone into stars (M∗) divided by the total mass involved in the collapse, i.e. the mass of the cloud (MGMC) plus the mass of the
Under the assumption that star formation proceeds by gravitational collapse on a free-fall time, it is useful to express the SFR as the SFE per unit free-fall time (ϵff). Typical values of the SFE per free-fall time on the cloud scale are of the order of 0.3–3% for GMCs in the Milky Way and in nearby galaxies (see Fig. 3 and e.g. Krumholz and Tan 2007; Kennicutt and Evans 2012; Leroy et al. 2017; Utomo et al. 2018; Krumholz et al. 2019), although it increases slightly towards dense subregions within GMCs, reaching 3–6% at n∼5×102cm−3 (Evans et al. 2009). In other words, it would take several tens of free-fall times to convert all gas into stars.
Fig. 3 SFE per free-fall time ϵff (median and 16th–84th percentile range) measured in a sample of nearby galaxies, as a function of galaxy stellar mass. The red data points are calculated using a constant value of the CO-to-H
2conversion factor (αCO), while the grey data points are calculated using a metallicity-dependent αCO. Figure taken from Utomo et al. (2018)
Before proceeding, we note that care needs to be taken when defining the SFE (both integrated and per unit time), and especially when comparing observations and simulations, because the measurements can differ significantly due to fundamental differences in the definition (Grudić et al. 2019). In simulations, the integrated SFE and the instantaneous SFE per unit time are trivial to measure, because the total mass involved in the simulation is known from the initial conditions. However, in observed clouds the molecular mass reservoir MGMC(t) evolves – it can grow by the accretion and cooling of atomic gas and it can be depleted not only by star formation, but also by feedback-driven dispersal. As a result, any instantaneous observational measurement of MGMC(t) in a single GMC never encompasses the entire mass reservoir, resulting in a bias towards high SFEs. This limitation strongly affects studies relying on cloud matching (i.e. associating individual clouds to individual young stellar regions, e.g. Ochsendorf et al. 2017; Dessauges-Zavadsky et al. 2019, also see the discussion in Krumholz et al. 2019). Statistical approaches for measuring the SFE per free-fall time (e.g. Leroy et al. 2017; Utomo et al. 2018) and the integrated SFE (e.g. Kruijssen et al. 2019b; Chevance et al. 2020) are able to avoid this problem by comparing the free-fall time or the GMC lifetime to the galaxy-wide gas depletion time. This category of papers finds ϵff∼0.01 and integrated SFEs of 2–10%.
If most GMCs are transient and disperse within a few free-fall times without forming stars, then it may be straightforward to reproduce the observed level of star formation in galaxies. However, if most GMCs contain gravitationally bound regions and form stars at some point during their lifecycle (which is supported by observations, see Sect. 5), this either requires that star formation is slow (relative to the free-fall time) or inefficient (only a small fraction of the gas is actually converted into stars before the cloud is dispersed). These two possibilities are examined further below. Determining which physical mechanisms keep the SFE low in galaxies has remained one of the most challenging questions over the past 50 years. A promising way forward in distinguishing between the above scenarios is by measuring the molecular cloud lifetime, which is discussed in Sect. 5.
Based on the above discussion, two models of star formation can be (1) one in which star formation is slow (i.e. it takes many free-fall times) and efficient (i.e. a large fraction of the GMC is turned into stars), or (2) one in which star formation is rapid (i.e. it takes of the order a free-fall time) and inefficient (i.e. only a few percent of the GMC is turned into stars). The roles played by stellar feedback and magnetic fields are critical in the distinction between these two models.
One hypothesis is that the low SFE and SFR measured in galaxies result from slow star formation, and that the efficient conversion of gas into stars happens on timescales much longer that the free-fall time (e.g. Krumholz and Federrath 2019). In this case, some mechanisms must support the clouds against gravitational collapse (e.g. Zuckerman and Palmer 1974; Klessen et al. 2000; Krumholz and Federrath 2019). As mentioned above, the timescale suggested by equation (5) is a lower-limit to the actual collapse time. In addition to the delay induced by the thermal pressure gradient at early times in the collapse, several other factors may make the collection and contraction time of the gas longer. These include magnetic fields (even if the clouds are magnetically supercritical, i.e. the magnetic energy is less than the binding energy, see e.g. Inoue and Inutsuka 2009; Vázquez-Semadeni et al. 2011; Girichidis et al. 2018), turbulence (e.g. Klessen et al. 2000; Dobbs et al. 2014, and references therein), galactic differential rotation through shear and Coriolis forces (e.g. Dobbs et al. 2014; Meidt et al. 2018, 2020), and the non-spherical (planar or filamentary) shape of the clouds (Toalá et al. 2012; Pon et al. 2012).
Numerical simulations of GMC formation by compressions in the warm atomic gas under Solar Neigbourhood conditions suggest that global contraction begins several Myr after the compressions first occurred, and that local collapse starts several Myr after the whole GMC has engaged in gravitational contraction (e.g. Vázquez-Semadeni et al. 2007, 2011; Heitsch and Hartmann 2008; Heitsch et al. 2008a,b; Carroll-Nellenback et al. 2014). However, for most of the assembly time, the clouds may be in a mostly atomic form, and molecule formation may start almost simultaneously with local collapse and star formation (e.g. Hartmann et al. 2001; Bergin et al. 2004; Clark et al. 2012; Heiner et al. 2015; Vázquez-Semadeni et al. 2018), explaining the observation that, in the Solar Neighborhood (i.e. within 1 kpc from the Sun) most of the molecular clouds exhibit signs of (low-mass) star formation (e.g. Ballesteros-Paredes and Hartmann 2007; Kainulainen et al. 2009).
About 30 years ago, GMCs were thought to be fully supported by a relatively strong magnetic field, i.e. they were thought to be magnetically subcritical. In such a case, gravitational collapse can only be initiated after magnetic flux loss due to magnetic (e.g. ambipolar) diffusion (Shu et al. 1987). This model proposed that the star formation timescale is determined by the magnetic diffusion time (>100tff). This timescale is considerably longer than GMC lifetimes estimated from recent observations (see Sect. 5). However, Nakamura and Li (2008) demonstrated that local turbulent compression accelerates ambipolar diffusion, because the diffusion timescale is proportional to the magnetic force (which is enhanced in the compressed regions). Thus, even in the presence of subcritical magnetic fields, the gravitational collapse of small compact parts created by the turbulent compression may be initiated within a few free-fall times (see also Kudoh and Basu 2011). Nonetheless, the initial magnetic fields should be close to critical in order for the collapse timescale to not be much longer than the observed one. An advantage of this magnetically-supported model is that the cloud and the core envelopes remain magnetically-supported and therefore only some fraction of the gas (i.e. the magnetically supercritical part) can form stars (Nakamura and Li 2005). However, more recent measurements suggest that observed magnetic field strengths are insufficient to support GMCs against collapse (see the review by Crutcher 2012).
If the magnetic support is weaker, star formation is expected to proceed more efficiently and star clusters can be formed. For clustered star formation, numerical simulations show that stellar feedback such as protostellar jets, outflows, and stellar winds can inject supersonic turbulence in molecular clumps (Nakamura and Li 2007; Offner and Arce 2015), and the clumps can be kept near virial equilibrium for several dynamical timescales. Thus, in the context of slow star formation, the relative importance of turbulence and magnetic fields determines the mode of star formation (i.e. clustered or distributed), while thermal stellar feedback maintains the star-forming clumps within GMCs close to virial equilibrium.
Non-thermal motions can play a role in fragmenting the clouds (e.g. Kim and Ryu 2005; Ballesteros-Paredes et al. 2006) and delaying their collapse, therefore lowering the efficiency per free-fall time (e.g. Federrath 2015). However, 3D magnetohydrodynamic turbulent simulations demonstrate that supersonic turbulence decays quickly, in a turbulent crossing time (Stone et al. 1998; Mac Low et al. 1998), which implies that turbulence has to be driven continuously and homogeneously in order significantly delay collapse, either on very small scales or aided by magnetic fields (Klessen et al. 2000; Heitsch et al. 2001; Federrath and Klessen 2012). It is therefore unclear whether realistically driven turbulence can significantly delay collapse (e.g. Hennebelle and Iffrig 2014; Iffrig and Hennebelle 2015).
The second scenario for explaining the low SFE and SFR in molecular clouds is that star formation is fast and inefficient. In this scenario, star formation takes place within of the order a dynamical time (e.g. Elmegreen 2000; Hartmann et al. 2012), with quick disruption of the clouds due to the effect of feedback mechanisms (e.g. photodissociating radiation, stellar winds and SN explosions, see the review by Dale 2015, and references therein). If cloud disruption occurs soon after the onset of star formation within a cloud, then the resulting efficiency of the conversion of gas into stars will be low. This idea is supported by recent observations showing that young stellar clusters and H ii regions with ages in the range 1–10 Myr have already dispersed their parent cloud (Kawamura et al. 2009; Whitmore et al. 2014; Hollyhead et al. 2015; Corbelli et al. 2017; Grasha et al. 2019; Hannon et al. 2019; Kruijssen et al. 2019b; Chevance et al. 2020). The various feedback mechanisms that may be responsible for the rapid dispersal of the natal cloud are discussed in Sect. 4.
In this model of rapid star formation, large-scale flows or turbulent motions, such as converging flows, SNe or bubbles, primarily control cloud-scale star formation (Hennebelle et al. 2008). Accretion of gas can bring significant turbulent motions into a compressed layer where molecular clouds form. The accretion-driven turbulence injected in the compressed layer is transonic in the early stage, but later on, strongly supersonic turbulence appears due to global gravitational contraction. Magnetic fields need to be relatively weak in order to enable star formation to take place on a dynamical time. This scenario assumes that star formation proceeds over a short period of time, but is sufficient to result in the rapid dispersal of the whole parent cloud (see Fig. 4 and Sect. 4). Only a small fraction of the cloud is effectively converted into stars and star formation is therefore globally inefficient. Statistically speaking, rapid star formation is self-regulated, because it is feedback-limited. Each individual region may exhibit strong (and long-lived) excursions from equilibrium, either by collapse or by feedback-driven dispersal once a sufficient integrated SFE has been reached. Due to this dynamic cycle, the GMC population will on average exhibit a low effective SFE per free-fall time, because only a small fraction of each cloud is converted into stars once a sufficient number of massive stars has formed (e.g. Grudić et al. 2018; Rahner et al. 2019).
Fig. 4 Evolution of the gas surface density (orange) and young stars (blue) for a simulated GMC of 3×107M⊙ and an initial radius of 50 pc, from left (initial conditions) to right (after 3.6 Myr). The top row is showing the edge-on view of the disc and the bottom row, the face on view. Star formation starts within a free-fall time (here 1.2 Myr) and stops after a few Myr, when the gas has been completely blown out by feedback. Figure taken from Grudić et al. (2018)
The hierarchical structure of GMCs enables star formation to start before a global free-fall time has been completed. This contributes to rapid GMC dispersal – once the densest regions begin to form massive stars, stellar feedback begins to erode the cloud, and gas dispersal can occur well before a substantial fraction of the mass is converted into stars. If star formation indeed proceeds rapidly and inefficiently, then the fact that only a small amount of mass is locked into dense cores (e.g. Kainulainen et al. 2009) occurs because GMCs are young and have not completed their fragmentation into cores yet, as opposed to fragmentation being slowed down by an opposing force such as magnetic fields.
In the two scenarios contrasted above, i.e. of slow and efficient star formation and of a fast and inefficient star formation, critical roles are taken up by stellar feedback, turbulence, and magnetic fields. Observations of these processes on the cloud scale should be used to distinguish these scenarios. A promising way of doing this is by measuring the GMC lifetime as directly as possible. We discuss this in Sect. 5.
As discussed above, high-mass star and star cluster formation occur in dense clumps within GMCs. Cloud-cloud collisions have been proposed as a possible mechanism for forming such dense clumps. In the last few decades, there is a growing volume of literature claiming observational evidence of cloud-cloud collisions triggering star formation (Furukawa et al. 2009; Duarte-Cabral et al. 2011; Nakamura et al. 2012, 2014; Kang et al. 2010; Fukui et al. 2014). These papers are mainly motivated by the idea that a collision between clouds can rapidly create dense compact parts on the cloud interface of the clouds, so that the conditions that favour the formation of high-mass stars and star clusters are easily achieved.
It is observationally challenging to find clear evidence of ongoing cloud-cloud collisions. Most observational indicators used to identify cloud-cloud collisions are highly degenerate with other mechanisms.
Even if star formation may sometimes be triggered by a cloud-cloud collision, the general role of this hypothetical process in driving the galactic SFR remains uncertain. Recent numerical simulations of GMC evolution and star formation in galaxies show that the average GMC experiences between zero and a few collisions over its lifetime (Tasker and Tan 2009; Dobbs et al. 2015) and the analytic theory of Jeffreson and Kruijssen (2018) shows that the rate of cloud-cloud collisions in most galactic environments is too low to significantly affect the average cloud lifetime. While some studies suggest that cloud-cloud collisions represent at least one of the important mechanisms to trigger high-mass star and star cluster formation in galaxies (e.g. Scoville et al. 1986; Kobayashi et al. 2018), these results depend strongly on the assumed parameters such as the cloud number density, cloud size, and cloud lifetime. High numbers need to be adopted for each of these quantities to make cloud-cloud collisions a viable contributor to the galactic SFR (e.g. Tan 2000).
After accumulating material from the larger-scale ISM, some regions of the clouds eventually gain enough mass and reach sufficiently high densities for gravitational instability to set in. The standard interpretation of observational data is that typically ∼5–10% of the total cloud mass becomes unstable and proceeds to collapse. The rest of the material is supported by turbulent motions in approximate energy equipartition with the clouds self-gravity (see e.g. Blitz et al. 2007).
The onset of star formation in these regions marks a primary transition in the evolution of the cloud (e.g. Krumholz 2014; Klessen and Glover 2016) with stellar feedback processes from massive stars initiating the final cloud dispersal phase. This transition phase may be long, with mass gain and mass loss being approximately equal, because it may take a few Myr until the star-forming clumps have grown to sufficiently high masses and densities to form massive stars (e.g. Vázquez-Semadeni et al. 2017; Krause et al. 2020). Eventually, the energy and momentum input from newly formed star-forming regions begins to dominate and the parent cloud is dispersed by stellar feedback (e.g. Krumholz et al. 2014; Lopez et al. 2014; Rahner et al. 2017, 2019; Grudić et al. 2018; Haid et al. 2018; Kim et al. 2018; Kruijssen et al. 2019b; McLeod et al. 2020, among many others).
For young stellar clusters or associations with M∗≳103M⊙ there are three main forms of ultraviolet radiation, stellar winds, and SNe. Each of these processes provides a source of energy and momentum that acts to oppose the forces of gravity. We refer the reader to Krumholz et al. (2014, 2019) for a summary of stellar feedback processes, and to the reviews by Dale (2015) and Hennebelle (2018) for their treatment in numerical simulations.
Typically, stellar winds are the first form of feedback that becomes noticeable. As the winds from intermediate- to high-mass stars collide and thermalise, they produce very hot (T∼106–108K) bubbles (Dunne et al. 2003; Townsley et al. 2003), which are filled with collisionially ionised gas. While this gas remains hot, its high thermal pressure drives the expansion into the surrounding medium which is swept up and compressed into a thin dense shell (Weaver et al. 1977). Once the gas cools, the winds from the central stellar population push the remainder of the gas from the bubble into the shell. Thereafter, the wind momentum is deposited directly into the shell in the form of ram pressure. Once the first SN explosions occur they also deposit momentum at the inner working surface of the expanding shell.
In intermediate-mass protocluster clouds (up to several times 104M⊙), UV photoionising radiation is often the most important feedback mechanism for regulating the SFR (e.g. Matzner 2002). Numerical simulations of cloud-scale star formation and feedback confirm this picture. The addition of photoionising radiation to simulations including stellar winds and SNe reduces the SN remnant density by up to two orders of magnitude with respect to simulations that do not include photoionisation (e.g. Peters et al. 2017). Because the optical depth of the gas inside a wind bubble is very low in most environments, the radiation from the central stellar population is able to easily reach the dense shell surrounding the bubble (Gupta et al. 2016). Ultraviolet photons with energies E>13.6eV photoionise hydrogen in this shell, resulting in one of two either the entire shell becomes ionised, and the remaining UV photons can leak into the surrounding ISM, or only the inner layers become ionised and the outer layers of the shell remain neutral (e.g. Martínez-González et al. 2014). The photons absorbed in the shell do not only provide heat and potentially change the chemical state of the gas, but they also deposit momentum (Rybicki and Lightman 1986). This results in a force acting at the inner working surface of the shell pushing it radially outwards, away from the central stellar population. If this radiation pressure is sufficiently large, then it can become dynamically significant and can play a major role in driving the evolution of the shell (Tielens 2010; Draine 2011). One of the key factors that determines whether or not radiation pressure becomes significant is the efficiency with which radiation couples to the material making up the shell. This is set by the opacity and column density of the gas and dust absorbing the radiation. When all the ionising photons are absorbed, i.e. when the shell is optically thick to ionising radiation, the system is called “radiation bounded”. In that case, the coupling is efficient and momentum is transferred effectively. By contrast, there are also many H ii regions where the observed shell structure is optically thin to ionising radiation, so that the coupling between radiation and matter is not very effective (Seon 2009; Pellegrini et al. 2012).
Several models have been developed to describe the evolution of wind-blown or radiation-driven bubbles surrounding young massive stellar populations. Calculations in which the dynamics of the shell is dominated by the effect of winds are presented by Weaver et al. (1977), Chevalier and Clegg (1985), Mac Low and McCray (1988), Koo and McKee (1992), Canto et al. (2000), and Silich and Tenorio-Tagle (2013). Models that focus on radiation pressure are introduced by Krumholz and Matzner (2009), Murray et al. (2010), and Kim et al. (2016). It is important to note here that both processes should be considered simultaneously and self-consistently in order to get the right expansion dynamics of the shell, because the various feedback mechanisms combine in a non-linear way (Dale 2015). As pointed out by Rahner et al. (2017), one needs to include winds to obtain the correct density structure of the shell, which is a prerequisite for correctly computing the number of photons absorbed in the shell. The reason is that the recombination rate depends quadratically on the density, implying that knowledge of the column density alone is insufficient for computing the balance between absorption and recombination as photons travel outwards. Instead, one needs a complete description of the density and chemical state of the material as function of radius.
Rahner et al. (2017, 2019) present a simple, yet detailed model that includes all physical processes currently considered to be relevant for GMC dispersal, under the assumption of spherical symmetry. This model includes self-consistent descriptions for stellar winds, SNe, radiation pressure, ionisation and gravity, solving explicitly for the density structure adopted by the gas in response to the action of stellar feedback. It assumes that the internal pressure of the feedback-blown bubble is larger than the external pressure and that the shell surrounding the bubble is in quasi-hydrostatic equilibrium with the forces acting at its inner working surface (based on e.g. Pellegrini et al. 2007, 2011). The model predicts that radiation pressure dominates over winds only for the dispersal of very massive and dense clouds (at and above ∼106M⊙). For less massive or dense systems, stellar winds dominate the force budget. This still holds at low the momentum output by winds is decreased but radiation also couples more weakly with the shell, and therefore winds can still dominate over radiation, assuming that their momentum and energy couples sufficiently to the expanding shell.
At later times (≳4Myr; Leitherer et al. 2014), SNe become the main drivers of shell expansion. However, when integrated over the entire cloud lifetime, the momentum and energy input from SNe does not exceed the contributions from either winds or radiation pressure. Several simulations suggest that the effects of a SN on the parent cloud are relatively limited compared to other feedback mechanisms, especially without the effects of pre-SN feedback (e.g. Geen et al. 2016; Körtgen et al. 2016; Rey-Raposo et al. 2017). However, SNe can have a large impact on cloud dispersal after pre-processing by early feedback mechanisms such as photoionisation and stellar winds (Geen et al. 2016). Recent observations support the idea that early feedback mechanisms in the form of photoionisation and stellar winds as described above, are playing a major role in dispersing the cloud before the first SN explosion. This is discussed in more detail in Sects. 3.3 and 5.
The effects of the feedback mechanisms described above are dominated by massive stars, which are unlikely to form in low-mass GMCs. From stochastic stellar population models, only stellar populations with masses ≳100M⊙ are expected to have at least one SN, and even stellar populations with masses ≲700M⊙ still show ionising luminosities 50% lower than a fully sampled initial mass function (Krumholz et al. 2019). Low-mass clouds are therefore unlikely to be disrupted by these types of feedback. Protostellar outflows are thus the only feedback mechanism that can potentially quench star formation for stellar populations of a few 100M⊙. Low-mass GMCs may also simply disperse under the influence of local dynamics (if they are gravitationally unbound) or galactic shear.
Altogether, analytical calculations and numerical simulations suggest that an integrated SFE of order of a few percent is sufficient to disrupt the parent cloud (Colín et al. 2013; Geen et al. 2016; Grudić et al. 2018; Kim et al. 2018; Rahner et al. 2019). Simulations shows that early feedback mechanisms, such as photoionising radiation and radiation pressure can efficiently destroy the clouds a few Myr after the onset star formation, therefore strongly restricting the fraction of the cloud forming stars. This results in a low star formation efficiency integrated over the cloud lifetime (see Fig. 5), in line with observational estimates in nearby star-forming galaxies (Kruijssen et al. 2019b; Chevance et al. 2020).
Fig. 5 SFE as a function of the initial GMC surface density for simulations with photoionisation only (PH-only, red circles), with radiation pressure only (RP-only, blue squares), and with both photoionisation and radiation pressure (PH+RP, dark grey stars). The size of the symbols reflects the initial GMC mass. The green triangle shows the SFE for a fiducial GMC without feedback. Radiation pressure and photoionisation (dominating in particular at gas surface densities Σ0<300M⊙pc−2) efficiently destroy the clouds, limiting the SFE to a few percent at the GMC surface densities typically observed in nearby disc galaxies (Σ0<100M⊙pc−2, e.g. Sun et al. 2018). Figure taken from Kim et al. (2018)
For efficiencies that fall considerably short of this value, the central star stellar population may be too weak to provide the momentum and energy input needed to fully disrupt the parent cloud. In this case, stellar feedback may have produced an ionised bubble surrounded by a dense shell, but this shell has never managed to sweep up the entire cloud. Once the momentum and energy input from winds, radiation and SNe fades away, the self-gravity of the gas takes over again and leads to recollapse followed by a second phase of star formation. This cycle could potentially occur several times, and it could explain why there are clusters which consist of two or more apparently distinct generations of stars. For example, this is proposed for 30 Doradus in the LMC (Brandl et al. 1996; Rahner et al. 2018). Such significant age spreads are likely restricted to stellar associations (see e.g. Efremov and Elmegreen 1998), because the age spreads observed in stellar clusters are small (e.g. Kudryavtseva et al. 2012; Longmore et al. 2014).
While multiple populations have also been observed in many globular clusters (see Adamo et al. 2020; Krause et al. 2020), the above recollapse scenario cannot explain this observation, because globular clusters do not exhibit the [Fe/H] spread that would be expected for the chemical enrichment by the type II SNe occurring over multiple collapse cycles. Generally speaking, the competition between feedback and gravity in one-dimensional models can have no other outcome than radial expansion or radial (re)collapse. However, the extension to two or three dimensions enables the introduction of shear, turbulence, torques, and external feedback, all of which imply that the focus of any form of recollapse likely deviates from the source of feedback. As a result, it is possible that subsequent generations of star formation would not occupy the same volume as the original one, implying that they do not necessarily affect its age spread or integrated SFE.
The spatial structure of a GMC subjected to stellar feedback can play an important role in setting the impact of stellar feedback on GMC dispersal. A common setup for investigating the effect of feedback on the clouds is to start with a spherical cloud and stir it with a numerical turbulence driver (see e.g. Dale 2015). However, it is not clear that this is a realistic setup, because clouds are likely to form in a sheet-like or filamentary fashion. If GMCs form by turbulent compression in the warm atomic gas, or by falling into the stellar potential of spiral arms, they are expected to form flattened sheets, because one-dimensional compressions are more likely than two- or three-dimensional ones. This anisotropy is likely to grow with time, because thermal and gravitational instabilities tend to grow faster along the shortest dimension of a perturbation, producing sheets and filaments from initially triaxial or ellipsoidal configurations (Field 1965; Lin et al. 1965; Zel’Dovich 1970; Heitsch et al. 2008a).
The gravitational potential of sheet-like or filamentary clouds is less deep than that of a spheroidal cloud filled with the same gas volume density and covering the same spatial extent, and thus these substructured clouds are easier to disperse by feedback. Numerical simulations of intermediate-mass clouds that inherit their initial structure from compressions of the warm diffuse medium show that the clouds are readily destroyed within ∼10Myr after the first massive stars appear, leaving an unembedded stellar population (Colín et al. 2013; Zamora-Avilés et al. 2019). Similar simulations using a spherical cloud with an initial turbulent velocity field show less efficient cloud destruction by feedback (e.g. Dale et al. 2012).
Summarising the preceding sections of this review, delineates three important evolutionary phases constituting the GMC lifecycle.
The evolutionary cycling between these three phases is visualised in Fig. 6. It remains a major open question how the physical mechanisms governing each of these phases may change with the galactic environment.
Fig. 6 Distribution of the gas in the galaxy simulation of Semenov et al. (2018) in the plane spanned by gas number density (n) and velocity dispersion (σtot). The black contours enclose 68%, 95%, and 99% of the total gas mass. In the simulation, star formation is assumed to take place in gas at high density and low velocity dispersion (reflecting the conditions expected in the real Universe), as indicated by the grey dashed line. In the analytical model proposed by Semenov et al. (2017), gas can enter this regime through gravitational collapse, cooling, and turbulent energy dissipation (at a mass flow rate F+, indicated by the blue arrows; this corresponds to phase 1 in the text). Gas below the star formation threshold can potentially for stars (at a rate M˙⋆, indicated by the black arrow; this corresponds to phase 2 in the text). Gas can be ejected from the star-forming regime by stellar feedback and dynamical processes (at a rate F−, indicated by the red arrows; this corresponds to phase 3 in the text). Altogether, these rates characterise the matter flow between the three phases of the GMC lifecycle in galaxies. Figure taken from Semenov et al. (2018)
While the above summary of the key phases in GMC evolution sketches a relatively comprehensive picture of the physical mechanisms that each must be understood in order to describe the molecular cloud lifecycle, the underlying timescales on which these phases proceed are not known a priori. However, as discussed throughout this review, knowledge of these timescales holds the key to identifying several of the dominant physical processes and enables a comprehensive view of the GMC lifecycle. Initial studies of the GMC lifecycle often focused on a single (and often differing) evolutionary phase (such as GMC assembly, low-mass star formation, or dispersal by feedback from massive stars, see e.g. Scoville et al. 1979; Sanders et al. 1985; Elmegreen 2000; Hartmann et al. 2001). In addition, these studies generally adopted highly dissimilar methodological approaches to the problem, leading to greatly differing evolutionary timescales. While they all rely on some form of statistical inference, some previous works rely on object classification and number counts to infer timescales (e.g. Kawamura et al. 2009; Corbelli et al. 2017), whereas others follow GMCs on evolutionary streamlines (e.g. Engargiola et al. 2003; Meidt et al. 2015; Kruijssen et al. 2015), consider stellar age spreads (e.g. Hartmann et al. 2001; Grasha et al. 2018, thus excluding any ‘inert’ phase of GMC evolution), or consider the lifetimes of molecules rather than those of GMCs (e.g. Scoville et al. 1979; Koda et al. 2009). Finally, the empirical constraints on the GMC lifecycle were not only limited by the lack of a single methodological framework, but also by the lack of large data sets enabling a systematic census of the GMC lifecycle as a function of the galactic environment. Thanks to the recent development of novel analysis frameworks (e.g. Kruijssen et al. 2018; Semenov et al. 2018) and the arrival of surveys of the molecular ISM with ALMA, combining a high spatial resolution with a large field of view (e.g. Sun et al. 2018; Schinnerer et al. 2019b, Leroy et al. in prep.), both hurdles have recently been overcome. In conjunction with the recent major progress in numerical simulations of cloud-scale star formation and feedback (e.g. Dale 2015; Walch et al. 2015; Grudić et al. 2018; Haid et al. 2018; Kim and Ostriker 2018; Kim et al. 2018; Semenov et al. 2018), a consistent picture of the GMC lifecycle is now emerging.
Measurements of the timescales governing GMC evolution can answer two main empirical
The first of these questions can help understand why the galaxy-wide gas depletion time is two orders of magnitude longer than the dynamical times of GMCs (see Sects. 1 and 3), i.e. whether GMCs live for many dynamical times and convert a large fraction of their mass into stars, or if they live for one or few dynamical times and reach a low SFE. In addition, it may help address what fraction of GMCs disperses without forming stars (see below). The second of these questions can help understand which feedback mechanisms drive GMC dispersal, e.g. whether early, pre-SN feedback is responsible, or if SNe play an important role in GMC destruction. These different cases are quite straightforward to distinguish observationally (Schruba et al. 2010; Kruijssen and Longmore 2014). If feedback operates slowly and GMCs are long-lived, we expect tracers of molecular gas and massive star formation to be co-spatial on the cloud scale. However, if molecular gas and massive stars represent distinct evolutionary phases of a rapid lifecycle, then they should not be correlated on small scales, but often be observed in isolation.
Empirically, ALMA has enabled a major step towards characterising the GMC lifecycle. Observations can now reach resolutions of 50–100 pc across the nearby galaxy population out to 20 Mpc, both for molecular gas traced by CO and massive star formation traced by Hα or ultraviolet emission. This latter observation is an essential complement to CO data, because it provides an absolute ‘reference timescale’ that enables translating the relative lifetimes of regions bright in CO and star formation rate tracers to absolute timescales (see Haydon et al. 2018). High-resolution observations of gas and star formation in nearby galaxies now show that CO and Hα emission rarely coincide on the cloud scale (Kreckel et al. 2018; Kruijssen et al. 2019b; Schinnerer et al. 2019a). Kruijssen et al. (2019b) and Chevance et al. (2020) used this empirical result to constrain the GMC lifecycle in the nearby flocculent spiral galaxy NGC300 and to nine nearby star-forming spiral galaxies observed as part of the PHANGS-ALMA survey (Leroy et al. in prep.), respectively. The observations exhibit a universal decorrelation of molecular gas and massive stars on GMC scales, implying a rapid evolutionary lifecycle, with short-lived clouds and rapid GMC dispersal by pre-SN feedback. By measuring how the CO-to-Hα flux ratio deviates from the galactic average near regions bright in CO or near those bright in Hα, they obtain a quantitative measurement of the GMC lifetime and the time taken for stellar feedback to drive GMC dispersal in these galaxies (the ‘feedback timescale’) (see Kruijssen et al. 2018 for details).
Figure 7 shows the evolutionary timelines of GMC evolution, star formation, and feedback measured by Kruijssen et al. (2019b, NGC300), Chevance et al. (2020, nine nearby disc galaxies), and Hygate et al. (2019, M33). GMC lifetimes range from 10–30 Myr and exhibit a slight trend of lifetimes increasing with galaxy mass. The GMC lifecycle is characterised by a long inert phase, without any unembedded massive star formation, that encompasses 75–90% of the GMC lifetime. Once unembedded massive stars appear, GMCs are dispersed rapidly, within 1–5 Myr, often due to early, pre-SN feedback (e.g. photoionisation and stellar winds), because SN explosions only occur after a ∼4Myr delay (Leitherer et al. 2014). By measuring the GMC lifetime, it is possible to infer the integrated star formation efficiency per star formation event. This efficiency is otherwise inaccessible, because it is defined as the ratio between the GMC lifetime and the galaxy-wide molecular gas depletion time. The measured GMC lifetimes are much shorter than the molecular gas depletion time (∼1Gyr, e.g. Bigiel et al. 2008; Leroy et al. 2008), implying that GMCs achieve low integrated star formation efficiencies, ranging from 2–10%. The homogeneous census across eleven nearby star-forming galaxies shown in Fig. 7 thus demonstrates that star formation is fast and inefficient – GMCs disperse rapidly, most likely due to early stellar feedback, such that only a small fraction of their mass is converted into stars.
Fig. 7 Evolutionary timeline of the GMC lifecycle from molecular gas to star formation and feedback, for a sample of eleven nearby galaxies. The first phase (in orange) indicates the duration of the ‘inert’ CO phase, without any signs of massive star formation. During the second phase (in maroon), gas and massive stars coexist. The third phase (in blue) represents the isolated young stellar phase, after the gas has dispersed and when only Hα emission is visible. Galaxies are ordered from top to bottom by increasing stellar mass. This diagram is based on Fig. 4 of Chevance et al. (2020), with the addition of NGC300 (from Kruijssen et al. 2019b) and M33 (from Hygate et al. 2019)
A key result of Fig. 7 is that GMC lifetimes are not universal, but vary from galaxy to galaxy. Chevance et al. (2020) and Kruijssen et al. (2019b) show that this variation does not only hold between galaxies, but also within individual galaxies, when characterising the GMC lifecycle in bins of galactocentric radius. The obvious question is what drives this environmental variation. Previous studies had already argued that GMCs are dynamical entities, evolving either on an internal dynamical (i.e. free-fall or crossing) time (e.g. Elmegreen 2000) or on a dynamical time-scale set by galactic dynamical processes (e.g. Dobbs et al. 2014). Examples of galactic dynamical processes that have been proposed to set GMC lifetimes are free-fall collapse of the midplane gas (e.g. Krumholz et al. 2012; Vázquez-Semadeni et al. 2019), shear (e.g. Elmegreen 1987; Dobbs and Pringle 2013), spiral arm passages (e.g. Meidt et al. 2013; Dobbs and Baba 2014), cloud-cloud collisions (e.g. Tan 2000; Takahira et al. 2014), and pericentre passages (or ‘epicyclic perturbations’, e.g. Longmore et al. 2013; Jeffreson et al. 2018). Jeffreson and Kruijssen (2018) derived an analytical model for GMC lifetimes under the influence of galactic dynamics that combines the timescales for the above processes through a harmonic sum (and thus assumes that the corresponding rates can be linearly added or subtracted).
Chevance et al. (2020) compare their observational measurements to the predicted internal and external dynamical times listed above. They find that evidence of two regimes of GMC lifetimes, separated by a critical kpc-scale mean gas surface density. At high surface densities (Σ>8M⊙pc−2), the GMC lifetime best matches the timescale predicted based on galactic dynamics from Jeffreson and Kruijssen (2018), with gravitational free-fall of the midplane ISM and shear being the dominant processes. At low surface densities (Σ<8M⊙pc−2), it best matches the internal dynamical timescale (i.e. the cloud-scale free-fall or crossing time). The physical interpretation of this result is that GMCs in high surface density environments reside in a (mostly) molecular medium, such that the detectable, CO-bright part of the cloud can extend beyond its tidal radius and the visible part of the GMC is sensitive to galactic dynamics. By contrast, GMCs in low surface density environments are ‘island GMCs’ that are decoupled from galactic dynamics, because they are the ‘tip of the iceberg’, surrounded by an extended atomic gas reservoir, and therefore evolve in isolation.
Interestingly, the fact that the observed GMC lifetime largely matches an (internal or external) dynamical time implies that GMCs in nearby galaxies on average do not undergo evolutionary cycles without massive star formation. The reason is that the methodology applied to measure the timescales in Fig. 7 measures the total time spent in a CO-bright phase before a Hα-bright phase emerges. If a GMC undergoes a lifecycle in which it does not form stars, disperses dynamically, forms again, and then does experience massive star formation, the starless cycle is added onto the measured total cloud lifetime. In such a scenario, the measured GMC lifetime would need to span at least three dynamical times (one to form, one to disperse, and one to form again). The observations rule out this possibility for the luminosity-weighted average GMC population, which is mostly biased towards massive (≳105M⊙) GMCs, but do still allow the possibility that low-mass, fainter GMCs experience starless lifecycles.
While previous literature results did not provide as wide a variety of galactic environments or as homogeneous an analysis as in Fig. 7, it is important to compare these recent measurements to previous results. The measurement of ‘short’ GMC lifetimes (i.e. of the order of a dynamical time) is qualitatively consistent with other measurements made during the past decade. By classifying and counting GMCs and H ii regions, Kawamura et al. (2009) found that GMCs in the Large Magellanic Cloud live for 20–30 Myr. Applying the same methodology, Corbelli et al. (2017) find a GMC lifetime of 14 Myr in M33 (without any quoted uncertainties). For the same galaxy, Hygate et al. (2019) find a lifetime of 16–23 Myr, which is qualitatively consistent. Using evolutionary streamlines, Meidt et al. (2015) find a GMC lifetime of 20–30 Myr in M51, again consistent with the measurement of 26–40 Myr obtained for the same galaxy by Chevance et al. (2020). While these measurements achieve broad consistency, the homogeneous application of a single analysis framework to a large sample of galaxies now rules out the possibility that differences between observed GMC lifetimes are caused by differences in methodology, and thus enables environmental trends to be cleanly identified.
The most compelling galaxy for which long (∼100Myr) GMC lifetimes had been reported in the literature has traditionally been M51 (Scoville et al. 1979; Koda et al. 2009). While it is now clear that part of the disagreement with other studies (e.g. Meidt et al. 2015) comes from the fact that the long timescale likely refers to a molecule lifetime rather than a GMC lifetime, there is also an important environmental factor. In the sample studied by Chevance et al. (2020), M51 is the galaxy with the longest GMC lifetime, even though it is only 30 Myr on average – a factor of three shorter than the inferred molecule lifetime. This example demonstrates the importance of both using homogenised methods and obtaining a sample large enough to reveal any environmental dependences.
The evolutionary timelines shown in Fig. 7 suggest rapid GMC dispersal by stellar feedback, on timescales of 1–5 Myr. In many cases, this requires early, pre-SN feedback. Kruijssen et al. (2019b) and Chevance et al. (in prep.) compare these measured feedback timescales to the expectations for various feedback mechanisms (assuming full coupling between the expanding shell and the ambient medium) and find that GMC dispersal is dominated by (predominantly) photoionisation, as well as stellar winds. GMCs in galaxies with the longest feedback timescales (4–5 Myr) may receive the final push towards dispersal from SNe. These results are consistent with those from previous observational studies, which used e.g. the ages of young stellar clusters and their association to nearby GMCs to estimate how quickly the feedback from the young stellar population drives GMC dispersal. For instance, Hollyhead et al. (2015) and Hannon et al. (2019) found that clusters in four nearby galaxies are no longer embedded once they reach ages of ∼4–5Myr. For clusters in M51, Grasha et al. (2019) find a GMC dispersal timescale of ∼6Myr, consistent with the measurement for the same galaxy by Chevance et al. (2020), who find 4.8−1.1+2.1Myr. The rapid dispersal of GMCs after massive star formation naturally explains why the integrated SFE is low. If star formation in GMCs typically accelerates with time (as has been suggested by e.g. Murray 2011), then the restriction of their lifetimes to 1–2 free-fall times also explains why the SFE per free-fall time is low.
The measured feedback timescales can be translated into characteristic velocities for GMC dispersal by dividing the GMC radius by the feedback timescale. This ‘feedback velocity’ is found to be 7–21kms−1 for the galaxies shown in Fig. 7. These numbers are consistent with the expansion velocities of H ii regions, which are directly observed in the Milky Way and the Large Magellanic Cloud using optical spectroscopy (6–30kms−1, see e.g. Bertoldi and McKee 1990; Murray and Rahman 2010; McLeod et al. 2019, 2020). The similarity between directly measured feedback velocities and those inferred from the feedback timescales is encouraging and shows that the measured feedback timescales are plausible.
Numerical models for the GMC lifecycle reveal a similar picture of highly dynamical, feedback-regulated, short GMC lifecycles with low star formation efficiencies (e.g. Dale et al. 2014; Gatto et al. 2017; Rahner et al. 2017, 2019; Grudić et al. 2018; Haid et al. 2018; Kim et al. 2018). In accordance with the interpretation of the observational measurements, these simulations highlight the importance of early, pre-SN feedback from photoionisation and stellar winds, as well as radiation pressure. These early feedback mechanisms are critical for reproducing the observed cloud lifecycle, but it is not clear how important they are for other GMC demographics such as masses, radii, and densities. Fujimoto et al. (2019) present a numerical simulation of an isolated Milky Way-like galaxy, which reproduces the observed GMC demographics and global gas depletion time typically found in nearby star-forming disc galaxies. However, the feedback from young stellar populations in this simulation is insufficient to disperse GMCs, causing CO and Hα emission to be correlated down to the cloud scale, in strong disagreement with observations. The observed cloud-scale decorrelation between tracers of molecular gas and massive star formation (Schruba et al. 2010; Kreckel et al. 2018; Kruijssen et al. 2019b; Schinnerer et al. 2019a; Chevance et al. 2020; Hygate et al. 2019) thus provides a fundamental test of how well numerical simulations reproduce the evolutionary lifecycle of GMCs in the real Universe, because it probes the GMC lifecycle more directly than the demographics of the GMC population (Fujimoto et al. 2019).
A major step made by numerical simulations during the last decade is to model the interplay between GMC-scale physics (such as star formation and feedback) and galaxy-scale processes (such as galactic dynamics). The ongoing growth of the spatial dynamic range spanned by simulations has recently made it possible to follow the galactic processes driving convergence (e.g. gravitational collapse, spiral arms) and dispersal (e.g. shear) down to resolving star formation and feedback in individual GMCs (e.g. Dobbs and Pringle 2013; Kim and Ostriker 2018; Tress et al. 2020, Jeffreson et al. in prep.). This is a major step towards understanding how the GMC lifecycle both drives and responds to galaxy evolution.
Taking together the results discussed above, the field has now reached the point at which the key phases of GMC formation, massive star formation, and feedback can be placed on an evolutionary timeline. Recent observations and simulations have made first steps towards understanding how this timeline may depend on the galactic environment. Across a wide variety of studies, the GMC lifecycle is now found to take place on a (galactic or internal) dynamical time (mostly governed by gravitational free-fall and shear), after which it is truncated by early stellar feedback from massive stars (mostly from photoionisation and stellar winds), resulting in low star formation efficiencies of up to a few percent (both integrated and per unit free-fall time). With large observational surveys and comprehensive numerical simulations that cover a wide parameter space of galactic environments at high spatial resolution, the community is very close to obtaining a systematic census of how the GMC lifecycle changes with the galactic environment, how it connects inflow and outflow processes in the ISM, and how it feeds the galactic baryon cycle.
In this review, we have described the characteristics of the GMC population and of the GMC lifecycle. We have shown that observationally measuring the durations of the successive phases of the evolutionary cycle of molecular clouds and star formation, from cloud assembly to cloud collapse and dispersal allows us to identify the relevant physical mechanisms at play, on the cloud scale in nearby disc galaxies. New theoretical developments combined with recent observations show that molecular clouds can be seen as the building blocks of galaxies. The cycle between molecular clouds and young stellar regions is rapid, driven by dynamics, self-gravity, and early stellar feedback (e.g. photoionisation and stellar winds), which disperses the clouds within a few Myr after the onset of massive star formation. In addition, this cycle is not universal, but the physical mechanisms controlling the different phases of this process likely depend on the environmental conditions. We have shown in particular that cloud lifetime may be set by the galactic dynamical timescale at high kpc-scale gas surface densities, whereas at low kpc-scale gas surface densities, GMCs appear to decouple from the galactic dynamics and their lifetime is regulated by internal dynamical processes.
To comprehensively constrain the relative roles of these mechanisms, and determine quantitatively how they depend on galactic structure and properties, future high-resolution, high-sensitivity, multi-wavelength observations across a large range of environments will be necessary, from the most quiescent (e.g. early-type galaxies), to the most star-forming (e.g. starburst galaxies, galaxy centres), also probing the particular case of high-redshift galaxies. This is only now becoming possible, thanks to the large recent and future observatories (ALMA, VLT, JWST) enabling multi-wavelength surveys of galaxies at the cloud-scale (such as PHANGS, see e.g. Sun et al. 2018; Schinnerer et al. 2019b; Leroy et al. in prep.). These observations will enable building a multi-scale model for star formation and feedback in galaxies, applicable across cosmic time.
Constructing a multi-scale model for star formation and feedback is becoming critical, because galaxy formation and evolution simulations are starting to reach these small (cloud) scale resolutions, even in large cosmological volumes (see e.g. Fig. 1 in Nelson et al. 2019). However, it remains computationally too demanding to treat the actual mechanisms of star-formation and feedback, which happen on the scales of individual stars (at sub-pc resolution), from first principles. Therefore, these simulations need to use sub-grid models for describing how gas is converted into stars and how energy and momentum is deposited by stellar feedback in the surrounding ISM. The new generation of these sub-grid models operates at the cloud scale, and can be informed by the new observational state of the art defined by ALMA, MUSE, and JWST. Additionally, in order to make reliable predictions for the demographics of the observed galaxy population at large, the cloud-scale predictions of simulations also need to be tested against similarly high resolution observations, as a function of the galactic environment. These cloud-scale predictions need to replicate specific observables, and most prominently the observed molecular cloud lifecycle. The recent study by Fujimoto et al. (2019) shows that this is not necessarily the case, even for simulations reproducing other observed, macroscopic (SFR, total gas mass, depletion time) and cloud-scale quantities (cloud sizes, masses, velocity dispersion). Comparing the observed and simulated molecular cloud lifecycles will make major contributions to better constraining the sub-grid physics used in galaxy formation and evolution simulations.
This dynamical vision of star formation and feedback in galaxies can be extended to larger scales. The next challenge is to characterise the physical processes driving the mass flows coupling the small-scale molecular cloud lifecycle to the galactic-scale baryon cycle, as a function of the environment. Eventually, combining all of these different elements will allow us to construct a multi-scale description of star formation across cosmic history.
Open Access funding provided by Projekt DEAL. We thank the staff of the International Space Science Institute (ISSI) for their generous hospitality and for creating a stimulating collaborative environment. M.C. and J.M.D.K. gratefully acknowledge funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through an Emmy Noether Research Group (grant number KR4801/1-1) and the DFG Sachbeihilfe (grant number KR4801/2-1). J.M.D.K. gratefully acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme via the ERC Starting Grant MUSTANG (grant agreement number 714907), and from the DFG via the SFB 881 “The Milky Way System” (subproject B2). R.S.K. acknowledges support from the DFG via the SFB 881 “The Milky Way System” (subprojects B1, B2, and B8) as well as funding from the Heidelberg Cluster of Excellence STRUCTURES in the framework of Germany’s Excellence Strategy (grant EXC-2181/1 – 390900948). J.B.P. acknowledges UNAM-DGAPA-PAPIIT support through grant number IN-111-219. A.A. acknowledges the support of the Swedish Research Council, Vetenskapsrådet, and the Swedish National Space Agency (SNSA).