It remains a major challenge to derive a theory of cloud-scale ( 100 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 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities Σ H 2 8 M pc −2 , the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at Σ H 2 8 M 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 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 HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100−300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

Using spectral-line observations of HNCO, N 2 H + , and HNC, we investigate the kinematics of dense gas in the central ∼ 250 pc of the Galaxy. We present scouse (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line-fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line profiles, scouse accurately describes the {l, b, v LSR } distribution of Central Molecular Zone (CMZ) gas, which is asymmetric about Sgr A* in both position and velocity. Velocity dispersions range from 2.6 km s −1 < σ < 53.1 km s −1 . A median dispersion of 9.8 km s −1 , translates to a Mach number, M 3D 28. The gas is distributed throughout several "streams", with projected lengths ∼ 100 − 250 pc. We link the streams to individual clouds and sub-regions, including Sgr C, the 20 and 50 km s −1 clouds, the dust ridge, and Sgr B2. Shell-like emission features can be explained by the projection of independent molecular clouds in Sgr C and the newly identified conical profile of Sgr B2 in {l, b, v LSR } space. These features have previously invoked supernova-driven shells and cloud-cloud collisions as explanations. We instead caution against structure identification in velocity-integrated emission maps. Three geometries describing the 3-D structure of the CMZ are investigated: i) two spiral arms; ii) a closed elliptical orbit; iii) an open stream. While two spiral arms and an open stream qualitatively reproduce the gas distribution, the most recent parameterisation of the closed elliptical orbit does not. Finally, we discuss how proper motion measurements of masers can distinguish between these geometries, and suggest that this effort should be focused on the 20 km s −1 and 50 km s −1 clouds and Sgr C.

We present PHANGS–ALMA, the first survey to map CO J = 2 → 1 line emission at ∼1″ ∼100 pc spatial resolution from a representative sample of 90 nearby (d ≲ 20 Mpc) galaxies that lie on or near the z = 0 “main sequence” of star-forming galaxies. CO line emission traces the bulk distribution of molecular gas, which is the cold, star-forming phase of the interstellar medium. At the resolution achieved by PHANGS–ALMA, each beam reaches the size of a typical individual giant molecular cloud, so that these data can be used to measure the demographics, life cycle, and physical state of molecular clouds across the population of galaxies where the majority of stars form at z = 0. This paper describes the scientific motivation and background for the survey, sample selection, global properties of the targets, Atacama Large Millimeter/submillimeter Array (ALMA) observations, and characteristics of the delivered data and derived data products. As the ALMA sample serves as the parent sample for parallel surveys with MUSE on the Very Large Telescope, the Hubble Space Telescope, AstroSat, the Very Large Array, and other facilities, we include a detailed discussion of the sample selection. We detail the estimation of galaxy mass, size, star formation rate, CO luminosity, and other properties, compare estimates using different systems and provide best-estimate integrated measurements for each target. We also report the design and execution of the ALMA observations, which combine a Cycle 5 Large Program, a series of smaller programs, and archival observations. Finally, we present the first 1″ resolution atlas of CO emission from nearby galaxies and describe the properties and contents of the first PHANGS–ALMA public data release.

The inner few hundred parsecs of the Milky Way harbours gas densities, pressures, velocity dispersions, an interstellar radiation field and a cosmic ray ionisation rate orders of magnitude higher than the disc; akin to the environment found in star-forming galaxies at high-redshift. Previous studies have shown that this region is forming stars at a rate per unit mass of dense gas which is at least an order of magnitude lower than in the disc, potentially violating theoretical predictions. We show that all observational star formation rate diagnostics -both direct counting of young stellar objects and integrated light measurements -are in agreement within a factor two, hence the low star formation rate is not the result of the systematic uncertainties that affect any one method. As these methods trace the star formation over different timescales, from 0.1 − 5 Myr, we conclude that the star formation rate has been constant to within a factor of a few within this time period. We investigate the progression of star formation within gravitationally bound clouds on ∼ parsec scales and find 1 − 4 per cent of the cloud masses are converted into stars per free-fall time, consistent with a subset of the considered "volumetric" star formation models. However, discriminating between these models is obstructed by the current uncertainties on the input observables and, most importantly and urgently, by their dependence on ill-constrained free parameters. The lack of empirical constraints on these parameters therefore represents a key challenge in the further verification or falsification of current star formation theories.

Infrared Dark Clouds (IRDCs) are unique laboratories to study the initial conditions of high-mass star and star cluster formation. We present high-sensitivity and high-angular resolution IRAM PdBI observations of N 2 H + (1 − 0) towards IRDC G035.39-00.33. It is found that G035.39-00.33 is a highly complex environment, consisting of several mildly supersonic filaments (σ NT /c s ∼ 1.5), separated in velocity by < 1 km s −1 . Where multiple spectral components are evident, moment analysis overestimates the non-thermal contribution to the linewidth by a factor ∼ 2. Large-scale velocity gradients evident in previous single-dish maps may be explained by the presence of substructure now evident in the interferometric maps. Whilst global velocity gradients are small (< 0.7 km s −1 pc −1 ), there is evidence for dynamic processes on local scales (∼ 1.5-2.5 km s −1 pc −1 ). Systematic trends in velocity gradient are observed towards several continuum peaks. This suggests that the kinematics are influenced by dense (and in some cases, starless) cores. These trends are interpreted as either infalling material, with accretion rates ∼ (7 ± 4)×10 −5 M yr −1 , or expanding shells with momentum ∼ 24 ±12 M km s −1 . These observations highlight the importance of high-sensitivity and high-spectral resolution data in disentangling the complex kinematic and physical structure of massive star forming regions.

We present the PHANGS-MUSE survey, a programme that uses the MUSE integral field spectrograph at the ESO VLT to map 19 massive (9.4 <  log(M⋆/M⊙)< 11.0) nearby (D ≲ 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1′ by 1′ each) and a total of nearly 15 × 106 spectra, covering ∼1.5 × 106 independent spectra. PHANGS-MUSE provides the first integral field spectrograph view of star formation across different local environments (including galaxy centres, bars, and spiral arms) in external galaxies at a median resolution of 50 pc, better than the mean inter-cloud distance in the ionised interstellar medium. This ‘cloud-scale’ resolution allows detailed demographics and characterisations of H II regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics and provides constraints on the star-formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2–1) and multi-band HST observations, therefore allowing us to probe the key stages of the star-formation process from molecular clouds to H II regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction, and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre.

Using the PHANGS-ALMA CO(2-1) survey, we characterize molecular gas properties on ∼100 pc scales across 102,778 independent sightlines in 70 nearby galaxies. This yields the best synthetic view of molecular gas properties on cloud scales across the local star-forming galaxy population obtained to date. Consistent with previous studies, we observe a wide range of molecular gas surface densities (3.4 dex), velocity dispersions (1.7 dex), and turbulent pressures (6.5 dex) across the galaxies in our sample. Under simplifying assumptions about subresolution gas structure, the inferred virial parameters suggest that the kinetic energy of the molecular gas typically exceeds its self-gravitational binding energy at ∼100 pc scales by a modest factor (1.3 on average). We find that the cloud-scale surface density, velocity dispersion, and turbulent pressure (1) increase toward the inner parts of galaxies, (2) are exceptionally high in the centers of barred galaxies (where the gas also appears less gravitationally bound), and (3) are moderately higher in spiral arms than in inter-arm regions. The galaxy-wide averages of these gas properties also correlate with the integrated stellar mass, star formation rate, and offset from the star-forming main sequence of the host galaxies. These correlations persist even when we exclude regions with extraordinary gas properties in galaxy centers, which contribute significantly to the inter-galaxy variations. Our results provide key empirical constraints on the physical link between molecular cloud populations and their galactic environment.

We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the cprops algorithms presented by Rosolowsky & Leroy. The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data via false source injection and includes an updated segmentation approach to better deal with blended emission. These improved algorithms are implemented in a publicly available Python package, pycprops. We apply these methods to 10 of the nearest galaxies in the PHANGS-ALMA survey, cataloguing CO emission at a common 90 pc resolution and a matched noise level. We measure the properties of 4986 individual clouds identified in these targets. We investigate the scaling relations among cloud properties and the cloud mass distributions in each galaxy. The physical properties of clouds vary among galaxies, both as a function of galactocentric radius and as a function of dynamical environment. Overall, the clouds in our target galaxies are well-described by approximate energy equipartition, although clouds in stellar bars and galaxy centres show elevated line widths and virial parameters. The mass distribution of clouds in spiral arms has a typical mass scale that is 2.5× larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4).