Abstract
We propose a simple analytic theory for environmentally dependent molecular cloud lifetimes, based on the large-scale (galactic) dynamics of the interstellar medium. Within this theory, the ...cloud lifetime is set by the time-scales for gravitational collapse, galactic shear, spiral arm interactions, epicyclic perturbations, and cloud–cloud collisions. It is dependent on five observable quantities, accessible through measurements of the galactic rotation curve, the gas and stellar surface densities, and the gas and stellar velocity dispersions of the host galaxy. We determine how the relative importance of each dynamical mechanism varies throughout the space of observable galactic properties, and conclude that gravitational collapse and galactic shear play the greatest role in setting the cloud lifetime for the considered range of galaxy properties, while cloud–cloud collisions exert a much lesser influence. All five environmental mechanisms are nevertheless required to obtain a complete picture of cloud evolution. We apply our theory to the galaxies M31, M51, M83, and the Milky Way, and find a strong dependence of the cloud lifetime upon galactocentric radius in each case, with a typical cloud lifetime between 10 and 50 Myr. Our theory is ideally suited for systematic observational tests with the Atacama Large Millimetre/submillimetre array.
ABSTRACT
We study the physical drivers of slow molecular cloud mergers within a simulation of a Milky Way-like galaxy in the moving-mesh code arepo, and determine the influence of these mergers on ...the mass distribution and star formation efficiency of the galactic cloud population. We find that 83 per cent of these mergers occur at a relative velocity below 5 km s−1, and are associated with large-scale atomic gas flows, driven primarily by expanding bubbles of hot, ionized gas caused by supernova explosions and galactic rotation. The major effect of these mergers is to aggregate molecular mass into higher-mass clouds: mergers account for over 50 per cent of the molecular mass contained in clouds of mass M > 2 × 106 M⊙. These high-mass clouds have higher densities, internal velocity dispersions and instantaneous star formation efficiencies than their unmerged, lower mass precursors. As such, the mean instantaneous star formation efficiency in our simulated galaxy, with its merger rate of just 1 per cent of clouds per Myr, is 25 per cent higher than in a similar population of clouds containing no mergers.
ABSTRACT
It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic ...environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially resolved (∼100 pc) CO-to-H α flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically $10\!-\!30\,{\rm Myr}$, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $\Sigma _{\rm H_2}\ge 8\,\rm {M_\odot}\,{{\rm pc}}^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $\Sigma _{\rm H_2}\le 8\,\rm {M_\odot}\,{{\rm pc}}^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H α (75–90 per cent of the cloud lifetime), GMCs disperse within just $1\!-\!5\,{\rm Myr}$ once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4–10 per cent. These results show that galactic star formation is governed by cloud-scale, environmentally dependent, dynamical processes driving rapid evolutionary cycling. GMCs and H ii regions are the fundamental units undergoing these lifecycles, with mean separations of $100\!-\!300\,{{\rm pc}}$ in star-forming discs. Future work should characterize the multiscale physics and mass flows driving these lifecycles.
ABSTRACT
We examine the role of the large-scale galactic-dynamical environment in setting the properties of giant molecular clouds in Milky Way-like galaxies. We perform three high-resolution ...simulations of Milky Way-like discs with the moving-mesh hydrodynamics code arepo, yielding a statistical sample of ${\sim}80\, 000$ giant molecular clouds and ${\sim}55\, 000$ H i clouds. We account for the self-gravity of the gas, momentum, and thermal energy injection from supernovae and H ii regions, mass injection from stellar winds, and the non-equilibrium chemistry of hydrogen, carbon, and oxygen. By varying the external gravitational potential, we probe galactic-dynamical environments spanning an order of magnitude in the orbital angular velocity, gravitational stability, mid-plane pressure, and the gradient of the galactic rotation curve. The simulated molecular clouds are highly overdense (∼100×) and overpressured (∼25×) relative to the ambient interstellar medium. Their gravoturbulent and star-forming properties are decoupled from the dynamics of the galactic mid-plane, so that the kpc-scale star formation rate surface density is related only to the number of molecular clouds per unit area of the galactic mid-plane. Despite this, the clouds display clear, statistically significant correlations of their rotational properties with the rates of galactic shearing and gravitational free-fall. We find that galactic rotation and gravitational instability can influence their elongation, angular momenta, and tangential velocity dispersions. The lower pressures and densities of the H i clouds allow for a greater range of significant dynamical correlations, mirroring the rotational properties of the molecular clouds, while also displaying a coupling of their gravitational and turbulent properties to the galactic-dynamical environment.
ABSTRACT
We present a novel, physically motivated sub-grid model for H ii region feedback within the moving mesh code arepo, accounting for both the radiation pressure-driven and thermal expansion of ...the ionized gas surrounding young stellar clusters. We apply this framework to isolated disc galaxy simulations with mass resolutions between 103 and 105 M⊙ per gas cell. Each simulation accounts for the self-gravity of the gas, the momentum and thermal energy from supernovae, the injection of mass by stellar winds, and the non-equilibrium chemistry of hydrogen, carbon, and oxygen. We reduce the resolution dependence of our model by grouping those H ii regions with overlapping ionization front radii. The Strömgren radii of the grouped H ii regions are at best marginally resolved, so that the injection of purely thermal energy within these radii has no effect on the interstellar medium. By contrast, the injection of momentum increases the fraction of cold and molecular gas by more than 50 per cent at mass resolutions of 103 M⊙, and decreases its turbulent velocity dispersion by ∼10 km s−1. The mass-loading of galactic outflows is decreased by an order of magnitude. The characteristic lifetime of the least-massive molecular clouds ($M/{\rm M}_\odot \lesssim 5.6 \times 10^4$) is reduced from ∼18 to $\lesssim 10$ Myr, indicating that H ii region feedback is effective in destroying these clouds. Conversely, the lifetimes of intermediate-mass clouds ($5.6 \times 10^4 \lesssim M/{\rm M}_\odot \lesssim 5 \times 10^5$) are elongated by ∼7 Myr, likely due to a reduction in supernova clustering. The derived cloud lifetimes span the range from 10 to 40 Myr, in agreement with observations. All results are independent of whether the momentum is injected from a ‘spherical’ or a ‘blister-type’ H ii region.
ABSTRACT
The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the emp-Pathfindersuite of cosmological zoom-in ...Milky Way-mass simulations. These simulations contain a subgrid description for stellar cluster formation and evolution, allowing us to study the simultaneous formation and evolution of stellar clusters alongside their host galaxies across cosmic time. As a key ingredient in these simulations, we include the physics of the multiphase nature of the interstellar medium (ISM), which enables studies of how the presence of a cold, dense ISM affects star cluster formation and evolution. We consider two different star formation prescriptions: a constant star formation efficiency per free-fall time, as well as an environmentally dependent, turbulence-based prescription. We identify two key results drawn from these simulations. First, we find that the tidal shock-driven disruption caused by the graininess of the cold ISM produces old ($\tau \gt 10~\mbox{${\rm Gyr}$}$) stellar cluster populations with properties that are in excellent agreement with the observed populations in the Milky Way and M31. Importantly, the addition of the cold ISM addresses the areas of disagreement found in previous simulations that lacked the cold gas phase. Secondly, we find that the formation of stellar clusters is extremely sensitive to the baryonic physics that govern the properties of the cold, dense gas reservoir in the galaxy. This implies that the demographics of the stellar cluster population represent an important diagnostic tool for constraining baryonic physics models in upcoming galaxy formation simulations that also include a description of the cold ISM.
ABSTRACT
We study the relationship between the scale height of the molecular gas disc and the turbulent velocity dispersion of the molecular interstellar medium within a simulation of a Milky ...Way-like galaxy in the moving-mesh code arepo.
We find that the vertical distribution of molecular gas can be described by a Gaussian function with a uniform scale height of ∼50 pc. We investigate whether this scale height is consistent with a state of hydrostatic balance between gravity and turbulent pressure. We find that the hydrostatic prediction using the total turbulent velocity dispersion (as one would measure from kiloparsec-scale observations) gives an overestimate of the true molecular disc scale height. The hydrostatic prediction using the velocity dispersion between the centroids of discrete giant molecular clouds (cloud–cloud velocity dispersion) leads to more accurate estimates. The velocity dispersion internal to molecular clouds is elevated by the locally enhanced gravitational field. Our results suggest that observations of molecular gas need to reach the scale of individual molecular clouds in order to accurately determine the molecular disc scale height.
ABSTRACT
We study the time evolution of molecular clouds across three Milky Way-like isolated disc galaxy simulations at a temporal resolution of 1 Myr and at a range of spatial resolutions spanning ...two orders of magnitude in spatial scale from ∼10 pc up to ∼1 kpc. The cloud evolution networks generated at the highest spatial resolution contain a cumulative total of ∼80 000 separate molecular clouds in different galactic–dynamical environments. We find that clouds undergo mergers at a rate proportional to the crossing time between their centroids, but that their physical properties are largely insensitive to these interactions. Below the gas–disc scale height, the cloud lifetime τlife obeys a scaling relation of the form τlife∝ℓ−0.3 with the cloud size ℓ, consistent with over-densities that collapse, form stars, and are dispersed by stellar feedback. Above the disc scale height, these self-gravitating regions are no longer resolved, so the scaling relation flattens to a constant value of ∼13 Myr, consistent with the turbulent crossing time of the gas disc, as observed in nearby disc galaxies.
ABSTRACT
We use passive gas tracer particles in an Arepo simulation of a dwarf spiral galaxy to relate the Lagrangian evolution of star-forming gas parcels and their H2 molecules to the evolution of ...their host giant molecular clouds. We find that the median chemical lifetime of H2 is 4 Myr, with an interquartile range between 2 and 9 Myr. This chemical lifetime is independent of the lifetime of the host molecular cloud, which may extend up to 90 Myr, with around 50 per cent of star formation occurring in longer lived clouds (>25 Myr). The rapid ejection of gas from around young massive stars by early stellar feedback is responsible for the short H2 survival time, driving down the density of the surrounding gas, so that its H2 molecules are dissociated by the interstellar radiation field. This ejection of gas from the H2-dominated state is balanced by the constant accretion of new gas from the galactic environment, constituting a ‘competition model’ for molecular cloud evolution. Gas ejection occurs at a rate that is proportional to the molecular cloud mass, so that the cloud lifetime is determined by the accretion rate, which may be as high as 4 × 104 M⊙ Myr−1 in the longest lived clouds. Our findings therefore resolve the conflict between observations of rapid gas ejection around young massive stars and observations of long-lived molecular clouds in galaxies. We show that the fastest-accreting, longest lived, highest mass clouds drive supernova clustering on sub-cloud scales, which in turn is a key driver of galactic outflows.
ABSTRACT
Connecting the gas in H ii regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how H ii regions evolve over time. ...With PHANGS–MUSE, we detect nearly 24 000 H ii regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, and density). We use catalogues of multiscale stellar associations from PHANGS–HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4177 H ii regions that are clearly linked to a single ionizing association. A weak anticorrelation is observed between the association ages and the $\mathrm{H}\, \alpha$ equivalent width $\mathrm{EW}(\mathrm{H}\, \alpha)$, the $\mathrm{H}\, \alpha/\mathrm{FUV}$ flux ratio, and the ionization parameter, log q. As all three are expected to decrease as the stellar population ages, this could indicate that we observe an evolutionary sequence. This interpretation is further supported by correlations between all three properties. Interpreting these as evolutionary tracers, we find younger nebulae to be more attenuated by dust and closer to giant molecular clouds, in line with recent models of feedback-regulated star formation. We also observe strong correlations with the local metallicity variations and all three proposed age tracers, suggestive of star formation preferentially occurring in locations of locally enhanced metallicity. Overall, $\mathrm{EW}(\mathrm{H}\, \alpha)$ and log q show the most consistent trends and appear to be most reliable tracers for the age of an H ii region.