The sensitivity and angular resolution of photometric surveys executed by the Hubble Space Telescope (HST) enable studies of individual star clusters in galaxies out to a few tens of megaparsecs. The fitting of spectral energy distributions (SEDs) of star clusters is essential for measuring their physical properties and studying their evolution. We report on the use of the publicly available Code Investigating GALaxy Emission (cigale) SED fitting package to derive ages, stellar masses, and reddenings for star clusters identified in the Physics at High Angular resolution in Nearby GalaxieS-HST (PHANGS-HST) survey. Using samples of star clusters in the galaxy NGC 3351, we present results of benchmark analyses performed to validate the code and a comparison to SED fitting results from the Legacy Extragalactic Ultraviolet Survey. We consider procedures for the PHANGS-HST SED fitting pipeline, e.g. the choice of single stellar population models, the treatment of nebular emission and dust, and the use of fluxes versus magnitudes for the SED fitting. We report on the properties of clusters in NGC 3351 and find, on average, the clusters residing in the inner star-forming ring of NGC 3351 are young (<10 Myr) and massive (10(5) M-circle dot) while clusters in the stellar bulge are significantly older. Cluster mass function fits yield beta values around -2, consistent with prior results with a tendency to be shallower at the youngest ages. Finally, we explore a Bayesian analysis with additional physically motivated priors for the distribution of ages and masses and analyse the resulting cluster distributions.

We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy. The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data via false source injection and includes an updated segmentation approach to better deal with blended emission. These improved algorithms are implemented in a publicly available PYTHON package, PYCPROPS. We apply these methods to 10 of the nearest galaxies in the PHANGS-ALMA survey, cataloguing CO emission at a common 90 pc resolution and a matched noise level. We measure the properties of 4986 individual clouds identified in these targets. We investigate the scaling relations among cloud properties and the cloud mass distributions in each galaxy. The physical properties of clouds vary among galaxies, both as a function of galactocentric radius and as a function of dynamical environment. Overall, the clouds in our target galaxies are well-described by approximate energy equipartition, although clouds in stellar bars and galaxy centres show elevated line widths and virial parameters. The mass distribution of clouds in spiral arms has a typical mass scale that is 2.5x larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4).

Statistical analysis of velocity fluctuations in the interstellar medium (ISM) of the Milky Way and NGC 4321 show that the motion of molecular gas over scales ranging from 0.1 to 1,000 pc is similar, and consistent with that generated by a combination of gravity and turbulence. ISM structure at one scale is therefore linked to structure at other scales.

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.