In this paper we examine the factors that shape the distribution of molecular gas surface densities on the 150 pc scale across 67 morphologically diverse star-forming galaxies in the PHANGS-ALMA CO (2-1) survey. Dividing each galaxy into radial bins, we measure molecular gas surface density contrasts, defined here as the ratio between a fixed high percentile of the CO distribution and a fixed reference level in each bin. This reference level captures the level of the faint CO floor that extends between bright filamentary features, while the intensity level of the higher percentile probes the structures visually associated with bright, dense interstellar medium features like spiral arms, bars, and filaments. We compare these contrasts to matched percentile-based measurements of the 3.6 mu m emission measured using Spitzer/IRAC imaging, which trace the underlying stellar mass density. We find that the logarithms of CO contrasts on 150 pc scales are 3-4 times larger than, and positively correlated with, the logarithms of 3.6 mu m contrasts probing smooth nonaxisymmetric stellar bar and spiral structures. The correlation appears steeper than linear, consistent with the compression of gas as it flows supersonically in response to large-scale stellar structures, even in the presence of weak or flocculent spiral arms. Stellar dynamical features appear to play an important role in setting the cloud-scale gas density in our galaxies, with gas self-gravity perhaps playing a weaker role in setting the 150 pc scale distribution of gas densities.

We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO(2-1) emission from nearby galaxies. The observations combine ALMA's main 12 m array, the 7 m array, and total power observations, and use mosaics of dozens to hundreds of individual pointings. We describe the processing of the u - v data, imaging and deconvolution, linear mosaicking, combining interferometer and total power data, noise estimation, masking, data product creation, and quality assurance. Our pipeline has a general design and can also be applied to Very Large Array and ALMA observations of other spectral lines and continuum emission. We highlight our recipe for deconvolution of complex spectral line observations, which combines multiscale clean, single-scale clean, and automatic mask generation in a way that appears robust and effective. We also emphasize our two-track approach to masking and data product creation. We construct one set of "broadly masked" data products, which have high completeness but significant contamination by noise, and another set of "strictly masked" data products, which have high confidence but exclude faint, low signal-to-noise emission. Our quality assurance tests, supported by simulations, demonstrate that 12 m+7 m deconvolved data recover a total flux that is significantly closer to the total power flux than the 7 m deconvolved data alone. In the appendices, we measure the stability of the ALMA total power calibration in PHANGS-ALMA and test the performance of popular short-spacing correction algorithms.

In this work, we measure the Ly alpha escape fraction of 935 [O III]-emitting galaxies between 1.9 < z < 2.35 by comparing stacked spectra from the Hubble Space Telescope/WFC3's near-IR grism to corresponding stacks from the Hobby-Eberly Telescope Dark Energy Experiment's Internal Data Release 2. By measuring the stacks' H beta to Ly alpha ratios, we determine the Ly alpha escape fraction as a function of stellar mass, star-formation rate, internal reddening, size, and [O III]/H beta ratio. We show that the escape fraction of Ly alpha correlates with a number of parameters, such as galaxy size, star-formation rate, and nebular excitation. However, we also demonstrate that most of these relations are indirect, and that the primary variables controlling the escape of Ly alpha are likely to be stellar mass and internal extinction. Overall, the escape of Ly alpha declines from greater than or similar to 16% in galaxies with logM/M-circle dot less than or similar to 9 to less than or similar to 1% for systems with logM/M-circle dot greater than or similar to 10, with the sample's mean escape fraction being 6.0(-0.5%)(+0.6%).

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 (Sigma (stellar) versus Sigma (SFR)), resolved Kennicutt-Schmidt relation (Sigma (mol. gas) versus Sigma (SFR)), and resolved molecular gas main sequence (Sigma (stellar) versus Sigma (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 Sigma (mol. gas+stellar) versus Sigma (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 Sigma (SFR) with respect to Sigma (stellar).

We present PHANGS-ALMA, the first survey to map CO J = 2 -> 1 line emission at similar to 1 '' similar to 100 pc spatial resolution from a representative sample of 90 nearby (d less than or similar to 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.

PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within similar to 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, HII 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 similar to 4 to similar 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.

Central molecular outflows in spiral galaxies are assumed to modulate their host galaxy's star formation rate (SFR) by removing gas from the inner region of the galaxy. Outflows consisting of different gas phases appear to be a common feature in local galaxies, yet, little is known about the frequency of molecular outflows in main sequence galaxies in the nearby universe. We develop a rigorous set of selection criteria, which allow the reliable identification of outflows in large samples of galaxies. Our criteria make use of central spectra, position-velocity diagrams and velocity-integrated intensity maps (line-wing maps). We use this method on high-angular resolution CO (2-1) observations from the PHANGS-ALMA survey, which provides observations of the molecular gas for a homogeneous sample of 90 nearby main sequence galaxies at a resolution of similar to 100 pc. We find correlations between the assigned outflow confidence and stellar mass or global SFR. We determine the frequency of central molecular outflows to be 25 +/- 2% considering all outflow candidates, or 20 +/- 2% for secure outflows only. Our resulting outflow candidate sample of 16-20 galaxies shows an overall enhanced fraction of active galactic nuclei (AGN) (50%) and bars (89%) compared to the full sample (galaxies with AGN: 24%, with bar: 61%). We extend the trend between mass outflow rates and SFR known for high outflow rates down to lower values (log(10) (M) over dot(out) [M-circle dot yr(-1)] < 0). Mass loading factors are of order unity, indicating that these outflows are not efficient in quenching the SFR in main sequence galaxies.

There is a severe tension between the observed star formation rate (SFR)-stellar mass (M,) relations reported by different authors at z = En addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations that tend to predict a factor of 2-4 smaller SFRs at a fixed M. We examine the evolution of the SFR-M, relation of z = 1-4 galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the Evolution and Assembly of GaLaxies and their Environments simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR-M* relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 jtm, MIPS- 24, 70, and 160 um or SPIRE- 250, 350, and 500 um) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by-0.5 dex at z 4 but overpredict the measurements by 0.3 dex at z <^> 1. Relations relying on dust corrected rest-frame ultraviolet luminosities, are flatter since they overpredict/underpredict SFRs for low/high star-forming objects and yield deviations from the intrinsic relation from 0.10 to-0.13 dex at z 4. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs.

We present the results of a proof-of-concept experiment that demonstrates that deep learning can successfully be used for production-scale classification of compact star clusters detected in Hubble Space Telescope (HST) ultraviolet-optical imaging of nearby spiral galaxies (D less than or similar to 20 Mpc) in the Physics at High Angular Resolution in Nearby GalaxieS (PHANGS)-HST survey. Given the relatively small nature of existing, human-labelled star cluster samples, we transfer the knowledge of state-of-the-art neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We perform a series of experiments to determine the dependence of classification performance on neural network architecture (ResNet18 and VGG19-BN), training data sets curated by either a single expert or three astronomers, and the size of the images used for training. We find that the overall classification accuracies are not significantly affected by these choices. The networks are used to classify star cluster candidates in the PHANGS-HST galaxy NGC 1559, which was not included in the training samples. The resulting prediction accuracies are 70 per cent, 40 per cent, 40-50 per cent, and 50-70 per cent for class 1, 2, 3 star clusters, and class 4 non-clusters, respectively. This performance is competitive with consistency achieved in previously published human and automated quantitative classification of star cluster candidate samples (70-80 per cent, 40-50 per cent, 40-50 per cent, and 60-70 per cent). The methods introduced herein lay the foundations to automate classification for star clusters at scale, and exhibit the need to prepare a standardized data set of human-labelled star cluster classifications, agreed upon by a full range of experts in the field, to further improve the performance of the networks introduced in this study.

We study the stellar populations and assembly of the nearby spiral galaxy NGC 2903's bulge, bar, and outer disc using the VIRUS-P Exploration of Nearby Galaxies IFS survey. We observe NGC 2903 with a spatial resolution of 185 pc using the Mitchell Spectrograph's 4.25 arcsec fibres at the 2.7 Harlan J. Smith telescope. Bulge-bar-disc decomposition on the Two Micron All-Sky Survey (2MASS) K-s-band image of NGC 2903 shows that it has similar to 6 per cent, 6 per cent, and 88 per cent, of its stellar mass in the bulge, bar, and outer disc, respectively, and its bulge has a low Sersic index of similar to 0.27, suggestive of a discy bulge. We perform stellar population synthesis and find that the outer disc has 46 per cent of its mass in stars >5 Gyr, 48 per cent in stars between 1 and 5 Gyr, and <10 per cent in younger stars. Its stellar bar has 65 per cent of its mass in ages 1-5 Gyr and has metallicities similar to the outer disc, suggestive of the evolutionary picture where the bar forms from disc material. Its bulge is mainly composed of old high-metallicity stars though it also has a small fraction of young stars. We find enhanced metallicity in the spiral arms and central region, tracing areas of high star formation as seen in the H alpha map. These results are consistent with the idea that galaxies of low bulge-to-total mass ratio and low bulge Sersic index like NGC 2903 has not had a recent major merger event, but has instead grown mostly through minor mergers and secular processes.