The spatial distribution of metals reflects, and can be used to constrain, the processes of chemical enrichment and mixing. Using PHANGS-MUSE optical integral field spectroscopy, we measure the gas-phase oxygen abundances (metallicities) across 7138 HII regions in a sample of eight nearby disc galaxies. In Paper I, we measure and report linear radial gradients in the metallicities of each galaxy, and qualitatively searched for azimuthal abundance variations. Here, we examine the 2D variation in abundances once the radial gradient is subtracted, Delta(O/H), in order to quantify the homogeneity of the metal distribution and to measure the mixing scale over which HII region metallicities are correlated. We observe low (0.03-0.05 dex) scatter in Delta(O/H) globally in all galaxies, with significantly lower (0.02-0.03 dex) scatter on small (<600 pc) spatial scales. This is consistent with the measurement uncertainties, and implies the 2D metallicity distribution is highly correlated on scales of less than or similar to 600 pc. We compute the two-point correlation function for metals in the disc in order to quantify the scale lengths associated with the observed homogeneity. This mixing scale is observed to correlate better with the local gas velocity dispersion (of both cold and ionized gas) than with the star formation rate. Selecting only HII regions with enhanced abundances relative to a linear radial gradient, we do not observe increased homogeneity on small scales. This suggests that the observed homogeneity is driven by the mixing introducing material from large scales rather than by pollution from recent and on-going star formation.

We measured gas-phase metallicity, ionisation parameter, and dust extinction for a representative sample of 1795 local star-forming galaxies using integral field spectroscopy from the SDSS-IV MaNGA survey. We self-consistently derive these quantities by comparing observed line fluxes with photoionisation models using a Bayesian framework. We also present the first comprehensive study of the [SIII]lambda lambda 9069,9532 nebular lines, which have long been predicted to be ideal tracers of the ionisation parameter. However, we find that current photoionisation model predictions substantially over-predict the intensity of the [SIII] lines, while broadly reproducing other observed optical line ratios. We discuss how to nonetheless make use of the information provided by the [SIII] lines by setting a prior on the ionisation parameter. Following this approach, we derive spatially resolved maps and radial profiles of metallicity and ionisation parameter. The metallicity radial profiles derived are comparable with previous works, with metallicity declining toward the outer parts and showing a flattening in the central regions. This is in agreement with infall models of galaxy formation, which predict that spiral discs build up through accretion of material, leading to an inside-out growth. On the other hand, ionisation parameter radial profiles are flat for low-mass galaxies, while their slope becomes positive as galaxy mass increases. However, the ionisation parameter maps we obtain are clumpy, especially for low-mass galaxies. The ionisation parameter is tightly correlated with the equivalent width of H alpha [EW(H alpha)], following a nearly universal relation, which we attribute to the change of the spectral shape of ionising sources due to ageing of HII regions. We derive a positive correlation between ionisation parameter and metallicity at fixed EW(H alpha), in disagreement with previous theoretical work that predict an anti-correlation.

The Sloan Digital Sky Survey V (SDSS-V) is an all-sky spectroscopic survey of >6 million objects, designed to decode the history of the Milky Way, reveal the inner workings of stars, investigate the origin of solar systems, and track the growth of supermassive black holes across the Universe. The Local Volume Mapper (LVM) is one of three surveys that form SDSS-V. LVM will employ a coordinated system of four telescopes feeding three fiber spectrographs at Las Campanas Observatory in Chile. The goal is to map approximately 2500 square degrees of the Galactic plane over the wavelength range 360-980 nm with R-4000 spectral resolution. These observations will reveal for the first time how distinct gaseous environments within the Galaxy interact with each other and with the stellar population, producing the large-scale interstellar medium that we observe.

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%).