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 just 4 Myr, independent of the lifetime of its host molecular cloud, which may vary from 1 to 90 Myr, with a substantial portion of all star formation in the galaxy occurring in relatively long-lived clouds. The rapid ejection of gas from around young massive stars by early stellar feedback is responsible for this 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 x 10^4 Msol/Myr 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, that often survive up to several tens of Myr. 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.
Past studies have long emphasised the key role played by galactic stellar bars in the context of disc secular evolution, via the redistribution of gas and stars, the triggering of star formation, and ...the formation of prominent structures such as rings and central mass concentrations. However, the exact physical processes acting on those structures, as well as the timescales associated with the building and consumption of central gas reservoirs are still not well understood. We are building a suite of hydro-dynamical RAMSES simulations of isolated, low-redshift galaxies that mimic the properties of the PHANGS sample. The initial conditions of the models reproduce the observed stellar mass, disc scale length, or gas fraction, and this paper presents a first subset of these models. Most of our simulated galaxies develop a prominent bar structure, which itself triggers central gas fuelling and the building of an over-density with a typical scale of 100−1000 pc. We confirm that if the host galaxy features an ellipsoidal component, the formation of the bar and gas fuelling are delayed. We show that most of our simulations follow a common time evolution, when accounting for mass scaling and the bar formation time. In our simulations, the stellar mass of 10 10 M ⊙ seems to mark a change in the phases describing the time evolution of the bar and its impact on the interstellar medium. In massive discs ( M ⋆ ≥ 10 10 M ⊙ ), we observe the formation of a central gas reservoir with star formation mostly occurring within a restricted starburst region, leading to a gas depletion phase. Lower-mass systems ( M ⋆ < 10 10 M ⊙ ) do not exhibit such a depletion phase, and show a more homogeneous spread of star-forming regions along the bar structure, and do not appear to host inner bar-driven discs or rings. Our results seem to be supported by observations, and we briefly discuss how this new suite of simulations can help our understanding of the secular evolution of main sequence disc galaxies.
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, and are associated with large-scale atomic gas flows, driven primarily by (1) expanding bubbles of hot, ionised gas caused by supernova explosions and (2) 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 x 10^6 Msun. 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.
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 kpc-scale observations) gives an over-estimate 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 Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do ...not form stars efficiently. We study here the ETG NGC 524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially resolved observations of the 12CO(2-1)-emitting cold molecular gas from the Atacama Large Millimeter/sub-millimeter Array (ALMA) and of the warm ionized-gas emission lines from SITELLE at the Canada–France–Hawaii Telescope. Although constrained by the resolution of the ALMA observations (≈37 pc), we identify only 52 GMCs with radii ranging from 30 to 140 pc, a low mean molecular gas mass surface density 〈Σgas〉 ≈ 125 M⊙ pc−2 and a high mean virial parameter 〈αobs, vir〉 ≈ 5.3. We measure spatially resolved molecular gas depletion times (τdep ≡ 1/SFE) with a spatial resolution of ≈100 pc within a galactocentric distance of 1.5 kpc. The global depletion time is ≈2.0 Gyr but τdep increases towards the galaxy centre, with a maximum τdep, max ≈ 5.2 Gyr. However, no pure H ii region is identified in NGC 524 using ionized-gas emission-line ratio diagnostics, so the τdep inferred are in fact lower limits. Measuring the GMC properties and dynamical states, we conclude that shear is the dominant mechanism shaping the molecular gas properties and regulating SF in NGC 524. This is supported by analogous analyses of the GMCs in a simulated ETG similar to NGC 524.
Abstract
We combine JWST observations with Atacama Large Millimeter/submillimeter Array CO and Very Large Telescope MUSE H
α
data to examine off-spiral arm star formation in the face-on, grand-design ...spiral galaxy NGC 628. We focus on the northern spiral arm, around a galactocentric radius of 3–4 kpc, and study two spurs. These form an interesting contrast, as one is CO-rich and one CO-poor, and they have a maximum azimuthal offset in MIRI 21
μ
m and MUSE H
α
of around 40° (CO-rich) and 55° (CO-poor) from the spiral arm. The star formation rate is higher in the regions of the spurs near spiral arms, but the star formation efficiency appears relatively constant. Given the spiral pattern speed and rotation curve of this galaxy and assuming material exiting the arms undergoes purely circular motion, these offsets would be reached in 100–150 Myr, significantly longer than the 21
μ
m and H
α
star formation timescales (both < 10 Myr). The invariance of the star formation efficiency in the spurs versus the spiral arms indicates massive star formation is not only triggered in spiral arms, and cannot simply occur in the arms and then drift away from the wave pattern. These early JWST results show that in situ star formation likely occurs in the spurs, and that the observed young stars are not simply the “leftovers” of stellar birth in the spiral arms. The excellent physical resolution and sensitivity that JWST can attain in nearby galaxies will well resolve individual star-forming regions and help us to better understand the earliest phases of star formation.
We present a novel, physically-motivated sub-grid model for HII region feedback within the moving mesh code Arepo, accounting for both the radiation pressure-driven and thermal expansion of the ...ionised gas surrounding young stellar clusters. We apply this framework to isolated disc galaxy simulations with mass resolutions between \(10^3~{\rm M}_\odot\) and \(10^5~{\rm M}_\odot\) 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 HII regions with overlapping ionisation front radii. The Str\"{o}mgren radii of the grouped HII 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 \(10^3~{\rm M}_\odot\), and decreases its turbulent velocity dispersion by \(\sim 10~{\rm kms}^{-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 < 5.6 \times 10^4\)) is reduced from \(\sim 18\) Myr to \(<10\) Myr, indicating that HII region feedback is effective in destroying these clouds. Conversely, the lifetimes of intermediate-mass clouds (\(5.6 \times 10^4 < M/{\rm M}_\odot < 5 \times 10^5\)) are elongated by \(\sim 7\) Myr, likely due to a reduction in supernova clustering. The derived cloud lifetimes span the range from \(10\)-\(40\) Myr, in agreement with observations. All results are independent of whether the momentum is injected from a 'spherical' or a 'blister-type HII region.
The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the EMP-Pathfinder suite of cosmological zoom-in Milky ...Way-mass simulations. These simulations contain a sub-grid 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 multi-phase nature of the interstellar medium (ISM), which enables studies of how the presence of a cold, dense ISM affects 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. Firstly, we find that tidal shock-driven disruption caused by the graininess of the cold ISM produces old (\(\tau>10~\)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, 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 stellar cluster populations represent an important diagnostic tool for constraining baryonic physics models in upcoming galaxy formation simulations that also include a description of the cold ISM.
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 obeys a scaling relation of the form \(\tau_{\rm life} \propto \ell^{-0.3}\) with the cloud size \(\ell\), 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.
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.