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.

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.

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.

We use integral field spectroscopy from the PHANGS–MUSE survey, which resolves the ionised interstellar medium structure at ∼50 pc resolution in 19 nearby spiral galaxies, to study the origin of the diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically defined H II regions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines in the DIG. A simple model for the leakage and propagation of ionising radiation from H II regions is able to reproduce the observed distribution of Hα in the DIG. This model infers a typical mean free path for the ionising radiation of 1.9 kpc for photons propagating within the disc plane. Leaking radiation from H II regions also explains the observed decrease in line ratios of low-ionisation species ([S II]/Hα, [N II]/Hα, and [O I]/Hα) with increasing Hα surface brightness (ΣHα). Emission from hot low-mass evolved stars, however, is required to explain: (1) the enhanced low-ionisation line ratios observed in the central regions of some of the galaxies in our sample; (2) the observed trends of a flat or decreasing [O III]/Hβ with ΣHα; and (3) the offset of some DIG regions from the typical locus of H II regions in the Baldwin–Phillips–Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LI[N]ERs). Hot low-mass evolved stars make a small contribution to the energy budget of the DIG (2% of the galaxy-integrated Hα emission), but their harder spectra make them fundamental contributors to [O III] emission. The DIG might result from a superposition of two components, an energetically dominant contribution from young stars and a more diffuse background of harder ionising photons from old stars. This unified framework bridges observations of the Milky Way DIG with LI(N)ER-like emission observed in nearby galaxy bulges.

Every star-forming galaxy has a halo of metal-enriched gas that extends out to at least 100 kiloparsecs, as revealed by the absorption lines that this gas imprints on the spectra of background quasars. However, quasars are sparse and typically probe only one narrow beam of emission through the intervening galaxy. Close quasar pairs and gravitationally lensed quasars have been used to circumvent this inherently one-dimensional technique, but these objects are rare and the structure of the circumgalactic medium remains poorly constrained. As a result, our understanding of the physical processes that drive the recycling of baryons across the lifetime of a galaxy is limited. Here we report integral-field (tomographic) spectroscopy of an extended background source-a bright, giant gravitational arc. We can thus coherently map the spatial and kinematic distribution of Mg ɪɪ absorption-a standard tracer of enriched gas-in an intervening galaxy system at redshift 0.98 (around 8 billion years ago). Our gravitational-arc tomography unveils a clumpy medium in which the absorption strength decreases with increasing distance from the galaxy system, in good agreement with results for quasars. Furthermore, we find strong evidence that the gas is not distributed isotropically. Interestingly, we detect little kinematic variation over a projected area of approximately 600 square kiloparsecs, with all line-of-sight velocities confined to within a few tens of kilometres per second of each other. These results suggest that the detected absorption originates from entrained recycled material, rather than in a galactic outflow.