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.

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.

(abridged) We review methods for measuring the sizes, line widths, and luminosities of giant molecular clouds (GMCs) in molecular-line data cubes with low resolution and sensitivity. We find that moment methods are robust and sensitive -- making full use of both position and intensity information -- and we recommend a standard method to measure the position angle, major and minor axis sizes, line width, and luminosity using moment methods. Without corrections for the effects of beam convolution and sensitivity to GMC properties, the resulting properties may be severely biased. This is particularly true for extragalactic observations, where resolution and sensitivity effects often bias measured values by 40% or more. We correct for finite spatial and spectral resolutions with a simple deconvolution and we correct for sensitivity biases by extrapolating properties of a GMC to those we would expect to measure with perfect sensitivity. The resulting method recovers the properties of a GMC to within 10% over a large range of resolutions and sensitivities, provided the clouds are marginally resolved with a peak signal-to-noise ratio greater than 10. We note that interferometers systematically underestimate cloud properties, particularly the flux from a cloud. The degree of bias depends on the sensitivity of the observations and the (u,v) coverage of the observations. In the Appendix to the paper we present a conservative, new decomposition algorithm for identifying GMCs in molecular-line observations. This algorithm treats the data in physical rather than observational units, does not produce spurious clouds in the presence of noise, and is sensitive to a range of morphologies. As a result, the output of this decomposition should be directly comparable among disparate data sets.

We measure the velocity dispersion, σ, and surface density, Σ, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution 45–120 pc, matched to the size of individual giant molecular clouds. Combining 11 galaxies from the PHANGS-ALMA survey with four targets from the literature, we characterize ∼30,000 independent sightlines where CO is detected at good significance. Σ and σ show a strong positive correlation, with the best-fit power-law slope close to the expected value for resolved, self-gravitating clouds. This indicates only a weak variation in the virial parameter α _vir ∝ σ ²/Σ, which is ∼1.5–3.0 for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure P _turb ∝ Σσ ², which spans ∼5 dex across our sample. We find Σ, σ, and P _turb to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kiloparsec of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the σ–Σ relation flattens at Σ ≲ 10 M _⊙ pc⁻², leading to high σ for a given Σ and high apparent α _vir. This echoes results found in the Milky Way and likely originates from a combination of lower beam-filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low-density regime.

PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ∼20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the data set will provide an unprecedented look into the connections between young stars, H ii regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g. brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy haloes in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 10 PHANGS galaxies from ∼4 to ∼15 Mpc, based on the first year of PHANGS-HST observations. Four of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, and NGC 4321), and seven of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release.