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.

Aims. The complexity of star formation at the physical scale of molecular clouds is not yet fully understood. We investigate the mechanisms regulating the formation of stars in different environments within nearby star-forming galaxies from the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) sample. Methods. Integral field spectroscopic data and radio-interferometric observations of 18 galaxies were combined to explore the existence of the resolved star formation main sequence (Σstellar versus ΣSFR), resolved Kennicutt–Schmidt relation (Σmol. gas versus ΣSFR), and resolved molecular gas main sequence (Σstellar versus Σmol. gas), and we derived their slope and scatter at spatial resolutions from 100 pc to 1 kpc (under various assumptions). Results. All three relations were recovered at the highest spatial resolution (100 pc). Furthermore, significant variations in these scaling relations were observed across different galactic environments. The exclusion of non-detections has a systematic impact on the inferred slope as a function of the spatial scale. Finally, the scatter of the Σmol. gas + stellar versus ΣSFR correlation is smaller than that of the resolved star formation main sequence, but higher than that found for the resolved Kennicutt–Schmidt relation. Conclusions. The resolved molecular gas main sequence has the tightest relation at a spatial scale of 100 pc (scatter of 0.34 dex), followed by the resolved Kennicutt–Schmidt relation (0.41 dex) and then the resolved star formation main sequence (0.51 dex). This is consistent with expectations from the timescales involved in the evolutionary cycle of molecular clouds. Surprisingly, the resolved Kennicutt–Schmidt relation shows the least variation across galaxies and environments, suggesting a tight link between molecular gas and subsequent star formation. The scatter of the three relations decreases at lower spatial resolutions, with the resolved Kennicutt–Schmidt relation being the tightest (0.27 dex) at a spatial scale of 1 kpc. Variation in the slope of the resolved star formation main sequence among galaxies is partially due to different detection fractions of ΣSFR with respect to Σstellar.

We use an unprecedented sample of about 23 000 H II regions detected at an average physical resolution of 67 pc in the PHANGS–MUSE sample to study the extragalactic H II region Hα luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z = 0. The mean LF slope, α, in our sample is =1.73 with a σ of 0.15. We find that α decreases with the galaxy’s star formation rate surface density, ΣSFR, and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the H II regions within single galaxies, we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds’ arm–inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in α between samples of H II regions with a high and low ionization parameter, q, and argue that they are driven by the aging of H II regions.

The PHANGS collaboration has been building a reference data set for the multiscale, multiphase study of star formation and the interstellar medium (ISM) in nearby galaxies. With the successful launch and commissioning of JWST, we can now obtain high-resolution infrared imaging to probe the youngest stellar populations and dust emission on the scales of star clusters and molecular clouds (∼5–50 pc). In Cycle 1, PHANGS is conducting an eight-band imaging survey from 2 to 21 μm of 19 nearby spiral galaxies. Optical integral field spectroscopy, CO(2–1) mapping, and UV-optical imaging for all 19 galaxies have been obtained through large programs with ALMA, VLT-MUSE, and Hubble. PHANGS–JWST enables a full inventory of star formation, accurate measurement of the mass and age of star clusters, identification of the youngest embedded stellar populations, and characterization of the physical state of small dust grains. When combined with Hubble catalogs of ∼10,000 star clusters, MUSE spectroscopic mapping of ∼20,000 H ii regions, and ∼12,000 ALMA-identified molecular clouds, it becomes possible to measure the timescales and efficiencies of the earliest phases of star formation and feedback, build an empirical model of the dependence of small dust grain properties on local ISM conditions, and test our understanding of how dust-reprocessed starlight traces star formation activity, all across a diversity of galactic environments. Here we describe the PHANGS–JWST Treasury survey, present the remarkable imaging obtained in the first few months of science operations, and provide context for the initial results presented in the first series of PHANGS–JWST publications.