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 an analysis of , [O iii]₈₈, [N ii]₁₂₂, and far-infrared (FIR) fine-structure line observations obtained with Herschel/PACS, for ∼240 local luminous infrared galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey. We find pronounced declines (“deficits”) of line-to-FIR continuum emission for [N ii]₁₂₂, , and as a function of FIR color and infrared luminosity surface density, . The median electron density of the ionized gas in LIRGs, based on the [N ii]₁₂₂/[N ii]₂₀₅ ratio, is = 41 cm⁻³. We find that the dispersion in the deficit of LIRGs is attributed to a varying fractional contribution of photodissociation regions (PDRs) to the observed emission, f( ) = / , which increases from ∼60% to ∼95% in the warmest LIRGs. The / ratio is tightly correlated with the PDR gas kinetic temperature in sources where is not optically thick or self-absorbed. For each galaxy, we derive the average PDR hydrogen density, , and intensity of the interstellar radiation field, G, in units of and find G/ ratios of ∼0.1–50 cm³, with ULIRGs populating the upper end of the distribution. There is a relation between G/ and , showing a critical break at ≃ 5 × 10¹⁰ L _⊙ kpc⁻². Below , G/ remains constant, ≃0.32 cm³, and variations in are driven by the number density of star-forming regions within a galaxy, with no change in their PDR properties. Above , G/ increases rapidly with , signaling a departure from the typical PDR conditions found in normal star-forming galaxies toward more intense/harder radiation fields and compact geometries typical of starbursting sources.

We use the GALPROP code for cosmic-ray (CR) propagation to calculate the broad-band luminosity spectrum of the Milky Way related to CR propagation and interactions in the interstellar medium. This includes gamma-ray emission from the production and subsequent decay of neutral pions, bremsstrahlung, and inverse Compton scattering, and synchrotron radiation. The Galaxy is found to be nearly a CR electron calorimeter, but {\it only} if gamma ray emitting processes are taken into account. Synchrotron radiation alone accounts for only one third of the total electron energy losses with ~10-20% of the total synchrotron emission from secondary CR electrons and positrons. The relationship between far-infrared and radio luminosity that we find from our models is consistent with that found for galaxies in general. The results will be useful for understanding the connection between diffuse emissions from radio through gamma rays in ``normal'' (non-AGN dominated) galaxies, as well as for estimating the broad-band extragalactic diffuse background from these kinds of galaxies.

Lyman-break galaxy (LBG) candidates at z ≳ 10 are rapidly being identified in James Webb Space Telescope (JWST)/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z ≲ 7) may also mimic the near-infrared (near-IR) colors of z > 10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z ≈ 5.1. We also present a tentative 2.6σ SCUBA-2 detection at 850 μm around a recently identified z ≈ 16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z ∼ 5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z ∼ 4–6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra-high redshift LBG candidates from JWST observations.

We study the spectral energy distribution (SED) of the radio continuum (RC) emission from the Key Insight in Nearby Galaxies Emitting in Radio (KINGFISHER) sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4, 4.8, 8.4, and 10.5 GHz combined with archive data allow us, for the first time, to determine the mid-RC (1–10 GHz, MRC) bolometric luminosities and further present calibration relations versus the monochromatic radio luminosities. The 1–10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo technique leading to measurements for the nonthermal spectral index ( ) and the thermal fraction ( ) with mean values of for the total spectral index) and = (10 ± 9)% at 1.4 GHz. The MRC luminosity changes over ∼3 orders of magnitude in the sample, MRC . The thermal emission is responsible for ∼23% of the MRC on average. We also compare the extinction-corrected diagnostics of the star-formation rate (SFR) with the thermal and nonthermal radio tracers and derive the first star-formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing SFR surface density, indicating the effect of the star-formation feedback on the cosmic-ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with SFR, due to the amplification of the magnetic fields in star-forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star-forming galaxies are detected.