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⊙)< 11.0) nearby (D ≲ 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1′ by 1′ each) and a total of nearly 15 × 106 spectra, covering ∼1.5 × 106 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 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.5× larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4).

The ground state of weakly bound dimers and trimers with a radius extending well into the classically forbidden region is explored, with the goal to test the predicted universality of quantum halo states. The focus of the study is molecules consisting of T↓, D↓, ^{3}He, ^{4}He, and alkali atoms, where the interaction between particles is much better known than in the case of nuclei, which are traditional examples of quantum halos. The study of realistic systems is supplemented by model calculations in order to analyze how low-energy properties depend on the interaction potential. The use of variational and diffusion Monte Carlo methods enabled a very precise calculation of both the size and binding energy of the trimers. In the quantum halo regime, and for large values of scaled binding energies, all clusters follow almost the same universal line. As the scaled binding energy decreases, Borromean states separate from tango trimers.

The PHANGS collaboration has been building a reference data set for the multiscale, multiphase study of star formation and the interstellar medium (ISM) in nearby galaxies. With the successful launch and commissioning of JWST, we can now obtain high-resolution infrared imaging to probe the youngest stellar populations and dust emission on the scales of star clusters and molecular clouds (∼5–50 pc). In Cycle 1, PHANGS is conducting an eight-band imaging survey from 2 to 21 μm of 19 nearby spiral galaxies. Optical integral field spectroscopy, CO(2–1) mapping, and UV-optical imaging for all 19 galaxies have been obtained through large programs with ALMA, VLT-MUSE, and Hubble. PHANGS–JWST enables a full inventory of star formation, accurate measurement of the mass and age of star clusters, identification of the youngest embedded stellar populations, and characterization of the physical state of small dust grains. When combined with Hubble catalogs of ∼10,000 star clusters, MUSE spectroscopic mapping of ∼20,000 H ii regions, and ∼12,000 ALMA-identified molecular clouds, it becomes possible to measure the timescales and efficiencies of the earliest phases of star formation and feedback, build an empirical model of the dependence of small dust grain properties on local ISM conditions, and test our understanding of how dust-reprocessed starlight traces star formation activity, all across a diversity of galactic environments. Here we describe the PHANGS–JWST Treasury survey, present the remarkable imaging obtained in the first few months of science operations, and provide context for the initial results presented in the first series of PHANGS–JWST publications.