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.

Dusty star-forming galaxies at high redshift (1 < z < 3) represent the most intense star-forming regions in the universe. Key aspects to these processes are the gas heating and cooling mechanisms, and although it is well known that these galaxies are gas-rich, little is known about the gas excitation conditions. Only a few detailed radiative transfer studies have been carried out owing to a lack of multiple line detections per galaxy. Here we examine these processes in a sample of 24 strongly lensed star-forming galaxies identified by the Planck satellite (LPs) at z ∼ 1.1–3.5. We analyze 162 CO rotational transitions (ranging from J up = 1 to 12) and 37 atomic carbon fine-structure lines ([C i]) in order to characterize the physical conditions of the gas in the sample of LPs. We simultaneously fit the CO and [C i] lines and the dust continuum emission, using two different non-LTE, radiative transfer models. The first model represents a two-component gas density, while the second assumes a turbulence-driven lognormal gas density distribution. These LPs are among the most gas-rich, IR-luminous galaxies ever observed (μ L L IR ( 8 − 1000 μ m ) ∼ 10 13 − 14.6 L ⊙; 〈 μ L M ISM 〉  = (2.7 ± 1.2) × 1012 M ⊙, with μ L ∼ 10–30 the average lens magnification factor). Our results suggest that the turbulent interstellar medium present in the LPs can be well characterized by a high turbulent velocity dispersion ( 〈 ΔV turb 〉  ∼ 100 km s−1) and ratios of gas kinetic temperature to dust temperature 〈 T kin/T d 〉  ∼ 2.5, sustained on scales larger than a few kiloparsecs. We speculate that the average surface density of the molecular gas mass and IR luminosity, Σ M ISM  ∼ 103–4 M ⊙ pc−2 and Σ L IR  ∼ 1011–12 L ⊙ kpc−2, arise from both stellar mechanical feedback and a steady momentum injection from the accretion of intergalactic gas.

We present EMPIRE, an IRAM 30 m large program that mapped λ = 3–4 mm dense gas tracers at ∼1–2 kpc resolution across the whole star-forming disk of nine nearby massive spiral galaxies. We describe the EMPIRE observing and reduction strategies and show new whole-galaxy maps of HCN(1−0), HCO+(1−0), HNC(1−0), and CO(1−0). We explore how the HCN-to-CO and IR-to-HCN ratios, observational proxies for the dense gas fraction and dense gas star formation efficiency, depend on host galaxy and local environment. We find that the fraction of dense gas correlates with stellar surface density, gas surface density, molecular-to-atomic gas ratio, and dynamical equilibrium pressure. In EMPIRE, the star formation rate per unit dense gas is anticorrelated with these same environmental parameters. Thus, although dense gas appears abundant in the central regions of many spiral galaxies, this gas appears relatively inefficient at forming stars. These results qualitatively agree with previous work on nearby galaxies and the Milky Way’s Central Molecular Zone. To first order, EMPIRE demonstrates that the conditions in a galaxy disk set the gas density distribution and that the dense gas traced by HCN shows an environment-dependent relation to star formation. However, our results also show significant (±0.2 dex) galaxy-to-galaxy variations. We suggest that gas structure below the scale of our observations and dynamical effects likely also play an important role.

We present kinematic orientations and high-resolution (150 pc) rotation curves for 67 main-sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of nonaxisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines of sight contained in the PHANGS CO velocity fields mapped at 1″ resolution. The high-resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log( ( )   M M )=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.