Low-mass stars like our Sun begin their evolution within cold (10 K) and dense (similar to 10(5) cm(-3)) cores of gas and dust. The physical structure of starless cores is best probed by thermal emission of dust grains. We present a high-resolution dust continuum study of the starless cores in the B10 region of the Taurus Molecular Cloud. New observations at 1.2 and 2.0 mm (12 and 18 arcsec resolution) with the NIKA2 instrument on the IRAM 30m have probed the inner regions of 14 low-mass starless cores. We perform sophisticated 3D radiative transfer modelling for each of these cores through the radiative transfer framework pandora, which utilizes RADMC-3D. Model best-fits constrain each cores' central density, density slope, aspect ratio, opacity, and interstellar radiation field strength. These 'typical' cores in B10 span central densities from 5 x 10(4) to 1 x 10(6) cm(-3), with a mean value of 2.6 x 10(5) cm(-3). We find the dust opacity laws assumed in the 3D modelling, as well as the estimates from Herschel, have dust emissivity indices, beta's, on the lower end of the distribution constrained directly from the NIKA2 maps, which averages to beta = 2.01 +/- 0.48. From our 3D density structures and archival NH3 data, we perform a self-consistent virial analysis to assess each core's stability. Ignoring magnetic field contributions, we find nine out of the 14 cores (64 per cent) are either in virial equilibrium or are bound by gravity and external pressure. To push the bounded cores back to equilibrium, an effective magnetic field difference of only similar to 15 mu G is needed.

We develop a method to identify and determine the physical properties of stellar associations using Hubble Space Telescope (HST) NUV-U-B-V-I imaging of nearby galaxies from the PHANGS-HST survey. We apply a watershed algorithm to density maps constructed from point source catalogues Gaussian smoothed to multiple physical scales from 8 to 64 pc. We develop our method on two galaxies that span the distance range in the PHANGS-HST sample: NGC 3351 (10 Mpc), NGC 1566 (18 Mpc). We test our algorithm with different parameters such as the choice of detection band for the point source catalogue (NUV or V), source density image filtering methods, and absolute magnitude limits. We characterise the properties of the resulting multi-scale associations, including sizes, number of tracer stars, number of associations, photometry, as well as ages, masses, and reddening from Spectral Energy Distribution fitting. Our method successfully identifies structures that occupy loci in the UBVI colour-colour diagram consistent with previously published catalogues of clusters and associations. The median ages of the associations increases from log(age/yr) = 6.6 to log(age/yr) = 6.9 as the spatial scale increases from 8 pc to 64 pc for both galaxies. We find that the youngest stellar associations, with ages < 3 Myr, indeed closely trace H ii regions in Hα imaging, and that older associations are increasingly anti-correlated with the Hα emission. Owing to our new method, the PHANGS-HST multi-scale associations provide a far more complete census of recent star formation activity than found with previous cluster and compact association catalogues. The method presented here will be applied to the full sample of 38 PHANGS-HST galaxies.

A long-standing problem when deriving the physical properties of stellar populations is the degeneracy between age, reddening, and metallicity. When a single metallicity is used for all star clusters in a galaxy, this degeneracy can result in $`$catastrophic$'$ errors for old globular clusters. Typically, approximately 10 - 20 % of all clusters detected in spiral galaxies can have ages that are incorrect by a factor of ten or more. In this paper we present a pilot study for four galaxies (NGC 628, NGC 1433, NGC 1365, and NGC 3351) from the PHANGS-HST survey. We describe methods to correct the age-dating for old globular clusters, by first identifying candidates using their colors, and then reassigning ages and reddening based on a lower metallicity solution. We find that young $`$interlopers$'$ can be identified from their Halpha flux. CO (2-1) intensity or the presence of dust can also be used, but our tests show that they do not work as well. Improvements in the success fraction are possible at the $\sim$ 15 % level (reducing the fraction of catastrophic age-estimates from between 13 - 21 % to 3 - 8 %). A large fraction of the incorrectly age-dated globular clusters are systematically given ages around 100 Myr, polluting the younger populations as well. Incorrectly age-dated globular clusters significantly impact the observed cluster age distribution in NGC 628, which affects the physical interpretation of cluster disruption in this galaxy. For NGC 1365, we also demonstrate how to fix a second major age-dating problem, where very dusty young clusters with E(B-V) $>$ 1.5 mag are assigned old, globular-cluster like ages. Finally, we note the discovery of a dense population of $\sim$ 300 Myr clusters around the central region of NGC 1365. and discuss how this results naturally from the dynamics in a barred galaxy.

In the hierarchical view of star formation, giant molecular clouds (GMCs) undergo fragmentation to form small-scale structures made up of stars and star clusters. Here we study the connection between young star clusters and cold gas across a range of extragalactic environments by combining the high resolution (1 '') PHANGS-ALMA catalogue of GMCs with the star cluster catalogues from PHANGS-HST. The star clusters are spatially matched with the GMCs across a sample of 11 nearby star-forming galaxies with a range of galactic environments (centres, bars, spiral arms, etc.). We find that after 4 - 6 Myr the star clusters are no longer associated with any gas clouds. Additionally, we measure the autocorrelation of the star clusters and GMCs as well as their cross-correlation to quantify the fractal nature of hierarchical star formation. Young (<= 10 Myr) star clusters are more strongly autocorrelated on kpc and smaller spatial scales than the > 10 Myr stellar populations, indicating that the hierarchical structure dissolves over time.

The PHANGS program is building the first data set to enable the multiphase, multiscale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with the Atacama Large Millimeter/submillimeter Array (ALMA), MUSE on the Very Large Telescope, and the Hubble Space Telescope (HST), with which we have obtained CO(2–1) imaging, optical spectroscopic mapping, and high-resolution UV–optical imaging, respectively. Here, we present PHANGS-HST, which has obtained NUV– U – B – V – I imaging of the disks of 38 spiral galaxies at distances of 4–23 Mpc, and parallel V - and I -band imaging of their halos, to provide a census of tens of thousands of compact star clusters and multiscale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ∼100,000 star clusters, associations, H ii regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey and provide an overview of the HST data processing pipeline and first results. We highlight new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. The PHANGS high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with the James Webb Space Telescope, for which a PHANGS Cycle 1 Treasury program to obtain eight-band 2–21 μ m imaging has been approved.