Understanding the spatial distribution of metals within galaxies allows us to study the processes of chemical enrichment and mixing in the interstellar medium. In this work, we map the 2D distribution of metals using a Gaussian Process Regression (GPR) for 19 star-forming galaxies observed with the Very Large Telescope/Multi Unit Spectroscopic Explorer (VLF-MUSE) as a part of the PRANGS-MUSE survey. We find that 12 of our 19 galaxies show significant 2D metallicity variation. Those without significant variations typically have fewer metallicity measurements, indicating this is due to the dearth of H II regions in these galaxies, rather than a lack of higher-order variation. After subtracting a linear radial gradient, we see no enrichment in the spiral arms versus the disc. We measure the 50 per cent correlation scale from the two-point correlation function of these radially subtracted maps, finding it to typically be an order of magnitude smaller than the fitted GPR kernel scale length. We study the dependence of the two-point correlation scale length with a number of global galaxy properties. We find no relationship between the 50 per cent correlation scale and the overall gas turbulence, in tension with existing theoretical models. We also find more actively star-forming galaxies, and earlier type galaxies have a larger 50 per cent correlation scale. The size and stellar mass surface density do not appear to correlate with the 50 per cent correlation scale, indicating that perhaps the evolutionary state of the galaxy and its current star formation activity is the strongest indicator of the homogeneity of the metal distribution.

It is a major open question which physical processes stop gas accretion on to giant molecular clouds (GMCs) and limit the efficiency at which gas is converted into stars. While feedback from supernova explosions has been the popular feedback mechanism included in simulations of galaxy formation and evolution, 'early' feedback mechanisms such as stellar winds, photoionization, and radiation pressure are expected to play an important role in dispersing the gas after the onset of star formation. These feedback processes typically take place on small scales (similar to 10-100 pc) and their effects have therefore been difficult to constrain in environments other than the Milky Way. We apply a novel statistical method to similar to 1 arcsec resolution maps of CO and H a across a sample of nine nearby galaxies, to measure the time over which GMCs are dispersed by feedback from young, high-mass stars, as a function of the galactic environment. We find that GMCs are typically dispersed within similar to 3 Myr on average after the emergence of unembedded high-mass stars, with variations within galaxies associated with morphological features rather than radial trends. Comparison with analytical predictions demonstrates that, independently of the environment, early feedback mechanisms (particularly photoionization and stellar winds) play a crucial role in dispersing GMCs and limiting their star formation efficiency in nearby galaxies. Finally, we show that the efficiency at which the energy injected by these early feedback mechanisms couples with the parent GMC is relatively low (a few tens of per cent), such that the vast majority of momentum and energy emitted by the young stellar populations escapes the parent GMC.

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-circle dot)< 11.0) nearby (D less than or similar to 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1 ' by 1 ' each) and a total of nearly 15 x 10(6) spectra, covering similar to 1.5 x 10(6) 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 show the results of a study using the spectral synthesis technique study for the full MaNGA sample showing their chemical enrichment history (ChEH) as well as the evolution of the stellar mass-metallicity relation (MZR) over cosmic time. We find that the more massive galaxies became enriched first and the lower-mass galaxies did so later, producing a change in the MZR that becomes shallower in time. Separating the sample into morphology and star-forming status bins, some particularly interesting results appear: The mass dependence of the MZR becomes less relevant for later morphological types, to the extent that it inverts for Sd/Irr galaxies, suggesting that morphology is at least as important a factor as mass in the chemical evolution. The MZR for the full sample shows a flattening at the high-mass end and another in the low-mass range, but the former only appears for retired galaxies, while the latter only appears for star-forming galaxies. We also find that the average metallicity gradient is currently negative for all mass bins, but for low-mass galaxies, it was inverted at some point in the past, before which all galaxies had a positive gradient. We also compare how diverse the ChEHs are in the different bins we considered, as well as what primarily drives the diversity: By how much galaxies become enriched, or how quickly they do so.

We present a rich, multiwavelength, multiscale database built around the PHANGS-ALMA CO (2 - 1) survey and ancillary data. We use this database to present the distributions of molecular cloud populations and subgalactic environments in 80 PHANGS galaxies, to characterize the relationship between population-averaged cloud properties and host galaxy properties, and to assess key timescales relevant to molecular cloud evolution and star formation. We show that PHANGS probes a wide range of kpc-scale gas, stellar, and star formation rate (SFR) surface densities, as well as orbital velocities and shear. The population-averaged cloud properties in each aperture correlate strongly with both local environmental properties and host galaxy global properties. Leveraging a variable selection analysis, we find that the kpc-scale surface densities of molecular gas and SFR tend to possess the most predictive power for the population-averaged cloud properties. Once their variations are controlled for, galaxy global properties contain little additional information, which implies that the apparent galaxy-to-galaxy variations in cloud populations are likely mediated by kpc-scale environmental conditions. We further estimate a suite of important timescales from our multiwavelength measurements. The cloud-scale freefall time and turbulence crossing time are similar to 5-20 Myr, comparable to previous cloud lifetime estimates. The timescales for orbital motion, shearing, and cloud-cloud collisions are longer, similar to 100 Myr. The molecular gas depletion time is 1-3 Gyr and shows weak to no correlations with the other timescales in our data. We publish our measurements online, and expect them to have broad utility to future studies of molecular clouds and star formation.

Aims. There exists some consensus that the stellar mass surface density (Sigma(star)) and molecular gas mass surface density (Sigma(mol)) are the main quantities responsible for locally setting the star formation rate. This regulation is inferred from locally resolved scaling relations between these two quantities and the star formation rate surface density (Sigma(SFR)), which have been extensively studied in a wide variety of works. However, the universality of these relations is debated. Here, we probe the interplay between these three quantities across different galactic environments at a spatial resolution of 150 pc.