Search

The earliest stages of star and cluster formation are hidden within dense cocoons of gas and dust, limiting their detection at optical wavelengths. With the unprecedented infrared capabilities of JWST, we can now observe dust-enshrouded star formation with $\sim$10 pc resolution out to $\sim$20 Mpc. Early findings from PHANGS-JWST suggest that 3.3 $\mu$m polycyclic aromatic hydrocarbon (PAH) emission can identify star clusters in their dust-embedded phases. Here, we extend this analysis to 19 galaxies from the PHANGS-JWST Cycle 1 Treasury Survey, providing the first characterization of compact sources exhibiting 3.3$\mu$m PAH emission across a diverse sample of nearby star-forming galaxies. We establish selection criteria, a median color threshold of F300M-F335M=0.67 at F335M=20, and identify of 1816 sources. These sources are predominantly located in dust lanes, spiral arms, rings, and galaxy centers, with $\sim$87% showing concentration indices similar to optically detected star clusters. Comparison with the PHANGS-HST catalogs suggests that PAH emission fades within $\sim$3 Myr. The H$\alpha$ equivalent width of PAH emitters is 1-2.8 times higher than that of young PHANGS-HST clusters, providing evidence that PAH emitters are on average younger. Analysis of the bright portions of luminosity functions (which should not suffer from incompleteness) shows that young dusty clusters may increase the number of optically visible $\leq$ 3 Myr-old clusters in PHANGS-HST by a factor between $\sim$1.8x-8.5x., Comment: Submitted to ApJ

Observations of young star clusters in a variety of galaxies have been used to constrain basic properties related to star-formation, such as the fraction of stars found in clusters (Gam) and the shape of the cluster mass function. However, the results can depend heavily on the reliability of the cluster age-dating process and other assumptions. One of the biggest challenges for successful age-dating lies in breaking the age-reddening degeneracy, where older, dust-free clusters and young, reddened clusters can have similar broad-band colors. While this degeneracy affects cluster populations in all galaxies, it is particularly challenging in dusty, extreme star-forming environments systems. We study the cluster demographics in the luminous infrared galaxy NGC1614 using Hubble imaging taken in 8 optical-NIR passbands. For age-dating, we adopt a spectral energy distribution fitting process that limits the maximum allowed reddening by region, and includes Ha photometry directly. We find that without these assumptions, essentially all clusters in the dust-free UV-bright arm which should have ages 50-250Myr are incorrectly assigned ages younger than 10Myr. We find this method greatly reduces the number of clusters in the youngest (tau<10Myrs) age bin and shows a fairly uniform distribution of massive clusters, the most massive being few10^7M. A maximum likelihood fit shows that the cluster mass function is well fitted by a power-law with an index -1.8, with no statistically significant high-mass cutoff. We calculate the fraction of stars born in clusters to be Gam1-10=22.4+_5.7%. The fraction of stars in clusters decreases quickly over time, with Gam10-100= 4.5+_1.1% and Gam100-400=1.7+_0.4%, suggesting that clusters dissolve rapidly over the first ~0.5Gyr. The decreasing fraction of stars in clusters is consistent with the declining shape observed for the cluster age distribution.

We present new catalogs of likely globular clusters (GCs) in 17 nearby spiral galaxies studied as part of the PHANGS-HST Treasury Survey. The galaxies were imaged in five broad-band filters from the near-ultraviolet through the $I$ band. PHANGS-HST has produced catalogs of stellar clusters of all ages by selecting extended sources (from multiple concentration index measurements) followed by morphological classification (centrally concentrated and symmetric or asymmetric, multiple peaks, contaminant) by visually examining the V-band image and separately by a machine-learning algorithm which classified larger samples to reach fainter limits. From both cluster catalogs, we select an initial list of candidate GCs to have $B-V \geq 0.5$ and $V-I \geq 0.73$~mag, then remove likely contaminants (including reddened young clusters, background galaxies misclassified by the neural network, and chance superpositions/blends of stars) after a careful visual inspection. We find that $\approx86$ % of the color-selected candidates classified as spherically symmetric, and $\approx68$ of those classified as centrally concentrated but asymmetric are likely to be GCs. The luminosity functions of the GC candidates in 2 of our 17 galaxies, NGC 628 and NGC 3627, are atypical, and continue to rise at least 1~mag fainter than the expected turnover near $M_V \sim -7.4$. These faint candidate GCs have more extended spatial distributions than their bright counterparts, and may reside in the disk rather than the bulge/halo, similar to faint GCs previously discovered in M101. These faint clusters may be somewhat younger since the age-metallicity degeneracy makes it difficult to determine precise cluster ages from integrated colors once they reach $\approx1$~Gyr., Comment: 21 pages, 12 figures

We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which have passed human inspection and morphological classification, and a larger sample of $\sim$100,000 classified by neural network models. Here, we report on the observed properties of these samples, and demonstrate that tremendous insight can be gained from just the observed properties of clusters, even in the absence of their transformation into physical quantities. In particular, we show the utility of the UBVI color-color diagram, and the three principal features revealed by the PHANGS-HST cluster sample: the young cluster locus, the middle-age plume, and the old globular cluster clump. We present an atlas of maps of the 2D spatial distribution of clusters and compact associations in the context of the molecular clouds from PHANGS-ALMA. We explore new ways of understanding this large dataset in a multi-scale context by bringing together once-separate techniques for the characterization of clusters (color-color diagrams and spatial distributions) and their parent galaxies (galaxy morphology and location relative to the galaxy main sequence). A companion paper presents the physical properties: ages, masses, and dust reddenings derived using improved spectral energy distribution (SED) fitting techniques., Comment: 48 pages, 26 figures, article in review at APJ

We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (${\sim}10^5\;\mathrm{M_\odot}$), small radii (${\lesssim}\,5\;\mathrm{pc}$), large escape velocities ($6{-}10\;\mathrm{km/s}$), and short free-fall times ($0.5{-}1\;\mathrm{Myr}$). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed HII region-cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline going from $\sim$1-2 Myr before cluster formation as starless clumps, to $\sim$4-6 Myr after as exposed HII region-cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit an non-uniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows., Comment: 27 pages, 12 figures; ApJ accepted

We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, linewidths, pressures, and virial parameters in the arms of NGC 1313 and center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environments, allowing more clouds with the necessary conditions to arise, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as clouds in the arms of NGC 1313, but have densities and pressures more similar to the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface densities between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions of lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density., Comment: 24 pages, 12 figures, accepted for publication in ApJ

We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times., Comment: 25 pages, 10 figures, accepted for publication in ApJ

Currently available star cluster catalogues from HST imaging of nearby galaxies heavily rely on visual inspection and classification of candidate clusters. The time-consuming nature of this process has limited the production of reliable catalogues and thus also post-observation analysis. To address this problem, deep transfer learning has recently been used to create neural network models which accurately classify star cluster morphologies at production scale for nearby spiral galaxies (D < 20 Mpc). Here, we use HST UV-optical imaging of over 20,000 sources in 23 galaxies from the Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) survey to train and evaluate two new sets of models: i) distance-dependent models, based on cluster candidates binned by galaxy distance (9-12 Mpc, 14-18 Mpc, 18-24 Mpc), and ii) distance-independent models, based on the combined sample of candidates from all galaxies. We find that the overall accuracy of both sets of models is comparable to previous automated star cluster classification studies (~60-80 per cent) and show improvement by a factor of two in classifying asymmetric and multi-peaked clusters from PHANGS-HST. Somewhat surprisingly, while we observe a weak negative correlation between model accuracy and galactic distance, we find that training separate models for the three distance bins does not significantly improve classification accuracy. We also evaluate model accuracy as a function of cluster properties such as brightness, colour, and SED-fit age. Based on the success of these experiments, our models will provide classifications for the full set of PHANGS-HST candidate clusters (N ~ 200,000) for public release., Comment: 16 pages, 10 figures

We investigate the ionised gas morphology, excitation properties, and kinematics in 19 nearby star-forming galaxies from the PHANGS-MUSE survey. We directly compare the kinetic energy of expanding superbubbles and the turbulent motions in the interstellar medium with the mechanical energy deposited by massive stars in the form of winds and supernovae, with the aim to answer whether the stellar feedback is responsible for the observed turbulent motions and to quantify the fraction of mechanical energy retained in the superbubbles. Based on the distribution of the flux and velocity dispersion in the H$\alpha$ line, we select 1484 regions of locally elevated velocity dispersion ($\sigma$(H$\alpha$)>45 km/s), including at least 171 expanding superbubbles. We analyse these regions and relate their properties to those of the young stellar associations and star clusters identified in PHANGS-HST data. We find a good correlation between the kinetic energy of the ionised gas and the total mechanical energy input from supernovae and stellar winds from the stellar associations, with a typical coupling efficiency of 10-20%. The contribution of mechanical energy by the supernovae alone is not sufficient to explain the measured kinetic energy of the ionised gas, which implies that pre-supernova feedback in the form of radiation/thermal pressure and winds is necessary. We find that the gas kinetic energy decreases with metallicity for our sample covering Z=0.5-1.0 Zsun, reflecting the lower impact of stellar feedback. For the sample of superbubbles, we find that about 40% of the young stellar associations are preferentially located in their rims. We also find a slightly higher (by ~15%) fraction of the youngest (<3 Myr) stellar associations in the rims of the superbubbles than in the centres, and the opposite for older associations, which implies possible propagation or triggering of star formation., Comment: Accepted by A&A. The abstract is abridged. 31 pages (including 8 pages in appendix), 20 figures, 2 tables

Connecting the gas in HII regions to the underlying source of the ionizing radiation can help us constrain the physical processes of stellar feedback and how HII regions evolve over time. With PHANGS$\unicode{x2013}$MUSE we detect nearly 24,000 HII regions across 19 galaxies and measure the physical properties of the ionized gas (e.g. metallicity, ionization parameter, density). We use catalogues of multi-scale stellar associations from PHANGS$\unicode{x2013}$HST to obtain constraints on the age of the ionizing sources. We construct a matched catalogue of 4,177 HII regions that are clearly linked to a single ionizing association. A weak anti-correlation is observed between the association ages and the H$\alpha$ equivalent width EW(H$\alpha$), the H$\alpha$/FUV flux ratio and the ionization parameter, log q. As all three are expected to decrease as the stellar population ages, this could indicate that we observe an evolutionary sequence. This interpretation is further supported by correlations between all three properties. Interpreting these as evolutionary tracers, we find younger nebulae to be more attenuated by dust and closer to giant molecular clouds, in line with recent models of feedback-regulated star formation. We also observe strong correlations with the local metallicity variations and all three proposed age tracers, suggestive of star formation preferentially occurring in locations of locally enhanced metallicity. Overall, EW(H$\alpha$) and log q show the most consistent trends and appear to be most reliable tracers for the age of an HII region., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

Star formation and stellar feedback are interlinked processes that redistribute energy and matter throughout galaxies. When young, massive stars form in spatially clustered environments, they create pockets of expanding gas termed superbubbles. As these processes play a critical role in shaping galaxy discs and regulating the baryon cycle, measuring the properties of superbubbles provides important input for galaxy evolution models. With wide coverage and high angular resolution (50-150 pc) of the PHANGS-ALMA $^{12}$CO (2-1) survey, we can now resolve and identify a statistically representative number of superbubbles with molecular gas in nearby galaxies. We identify superbubbles by requiring spatial correspondence between shells in CO with stellar populations identified in PHANGS-HST, and combine the properties of the stellar populations with CO to constrain feedback models and quantify their energetics. We visually identify 325 cavities across 18 PHANGS-ALMA galaxies, 88 of which have clear superbubble signatures (unbroken shells, central clusters, kinematic signatures of expansion). We measure their radii and expansion velocities using CO to dynamically derive their ages and the mechanical power driving the bubbles, which we use to compute the expected properties of the parent stellar populations driving the bubbles. We find consistency between the predicted and derived stellar ages and masses of the stellar populations if we use a supernova blast wave model that injects energy with a coupling efficiency of 10%, whereas continuous models fail to explain stellar ages we measure. Not only does this confirm molecular gas accurately traces superbubble properties, but it also provides key observational constraints for superbubble models. We also find evidence that the bubbles sweep up gas as they expand and speculate that these sites have the potential to host new generations of stars., Comment: 21 pages, 15 figures, 3 tables. Accepted to A&A. Abstract abridged for arXiv

The ultra-luminous infrared galaxy Arp2 20 is a late-stage merger with several tidal structures in the outskirts and two very compact, dusty nuclei that show evidence for extreme star formation and host at least one AGN. New and archival high-resolution images taken by the Hubble Space Telescope provide a state-of-the-art view of the structures, dust, and stellar clusters in Arp 220. We find that ~90% of the Halpha emission arises from a shock-ionized bubble emanating from the AGN in the western nucleus, while the nuclear disks dominate the Pbeta emission. Four very young (~3-6 Myr) but lower mass (< 10^4 Msun) clusters are detected in Halpha within a few arcsec of the nuclei, but produce less than 1% of the line emission. We see little evidence for a population of massive clusters younger than 100Myr anywhere in Arp 220. From the masses and ages of the detected clusters, we find that star formation took place more-or-less continuously starting ~few Gyr ago with a rate between ~3-12 Msun/yr. Approximately 100Myr ago, star formation shut off suddenly everywhere, except in the nuclear disks. A very recent flicker of weak star formation produced the four young, low-mass clusters, while the rest of the galaxy appears to have remained in a post-starburst state. Cluster ages indicate that the tidal structures on the west side of the galaxy are older than those on the east side, but all appear to pre-date the shutoff of star formation. Arp 220 has many of the characteristics expected of a 'Shocked Post-Starburst Galaxy' or SPOG, since most of the system has been in a post-starburst state for the past ~100Myr and the detected Halpha emission arises from shocked rather than photo-ionized gas., Comment: accepted for publication in the Astrophysical Journal

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., Comment: 27 pages, 17 figures (NOTE: 1, 5, 9, 12, 14 are lower resolution than in the journal)

PHANGS-JWST mid-infrared (MIR) imaging of nearby spiral galaxies has revealed ubiquitous filaments of dust emission in intricate detail. We present a pilot study to systematically map the dust filament network (DFN) at multiple scales between 25-400 pc in NGC 628. MIRI images at 7.7, 10, 11.3 and 21$\mu$m of NGC 628 are used to generate maps of the filaments in emission, while PHANGS-HST B-band imaging yields maps of dust attenuation features. We quantify the correspondence between filaments traced by MIR thermal continuum / polycyclic aromatic hydrocarbon (PAH) emission and filaments detected via extinction / scattering of visible light; the fraction of MIR flux contained in the DFN; and the fraction of HII regions, young star clusters and associations within the DFN. We examine the dependence of these quantities with the physical scale at which the DFN is extracted. With our highest resolution DFN maps (25 pc filament width), we find that filaments in emission and attenuation are co-spatial in 40% of sight lines, often exhibiting detailed morphological agreement; that ~30% of the MIR flux is associated with the DFN; and that 75-80% of HII regions and 60% of star clusters younger than 5 Myr are contained within the DFN. However, the DFN at this scale is anti-correlated with looser associations of stars younger than 5 Myr identified using PHANGS-HST near-UV imaging. We discuss the impact of these findings for studies of star formation and the ISM, and the broad range of new investigations enabled with multi-scale maps of the DFN., Comment: 21 pages, 14 figures, accepted for publication as part of PHANGS-JWST ApJL Focus Issue

A primary new capability of JWST is the ability to penetrate the dust in star forming galaxies to identify and study the properties of young star clusters that remain embedded in dust and gas. In this paper we combine new infrared images taken with JWST with our optical HST images of the star-bursting barred (Seyfert2) spiral galaxy NGC 1365. We find that this galaxy has the richest population of massive young clusters of any known galaxy within 30 Mpc, with $\sim$ 30 star clusters that are more massive than 10$^6$ Msolar and younger than 10 Myr. Sixteen of these clusters are newly discovered from our JWST observations. An examination of the optical images reveals that 4 of 30 ($\sim$13$\%$) are so deeply embedded that they cannot be seen in the I band (AV $\gt$ 10 mag), and that 11 of 30 ($\sim$37$\%$) are missing in the HST B band, so age and mass estimates from optical measurements alone are challenging. These numbers suggest that massive clusters in NGC 1365 remain obscured in the visible for $\sim$ 1.3 $\pm$ 0.7 Myr, and are either completely or partially obscured for $\sim$ 3.7 $\pm$ 1.1 Myr. We also use the JWST observations to gain new insights into the triggering of star cluster formation by the collision of gas and dust streamers with gas and dust in the bar. The JWST images reveal previously unknown structures (e.g., bridges and overshoot regions from stars that form in the bar) that help us better understand the orbital dynamics of barred galaxies and associated star-forming rings. Finally, we note that the excellent spatial resolution of the NIRCAM F200W filter provides a better way to separate barely resolved compact clusters from individual stars based on their sizes., Comment: 17 pages, 9 figures, accepted for publications as part of PHANGS-JWST ApJL Focus Issue

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$\alpha$ imaging, and that older associations are increasingly anti-correlated with the H$\alpha$ 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., Comment: Submitted to MNRAS. Referee report received with minor comments, and "request to clarify if the smaller associations are always included in the larger ones and how this may affect the photometric fitting of the larger association if the groups have different ages." Revision in progress

We compare embedded young massive star clusters (YMCs) to (sub-)millimeter line observations tracing the excitation and dissociation of molecular gas in the starburst ring of NGC 1365. This galaxy hosts one of the strongest nuclear starbursts and richest populations of YMCs within 20 Mpc. Here we combine near-/mid-IR PHANGS-JWST imaging with new ALMA multi-J CO (1-0, 2-1 and 4-3) and [CI](1-0) mapping, which we use to trace CO excitation via R42 = I_CO(4-3)/I_CO(2-1) and R21 = I_CO(2-1)/I_CO(1-0) and dissociation via RCICO = I_[CI](1-0)/I_CO(2-1) at 330 pc resolution. We find that the gas flowing into the starburst ring from northeast to southwest appears strongly affected by stellar feedback, showing decreased excitation (lower R42) and increased signatures of dissociation (higher RCICO) in the downstream regions. There, radiative transfer modeling suggests that the molecular gas density decreases and temperature and [CI/CO] abundance ratio increase. We compare R42 and RCICO with local conditions across the regions and find that both correlate with near-IR 2 um emission tracing the YMCs and with both PAH (11.3 um) and dust continuum (21 um) emission. In general, RCICO exhibits ~ 0.1 dex tighter correlations than R42, suggesting CI to be a more sensitive tracer of changing physical conditions in the NGC 1365 starburst than CO (4-3). Our results are consistent with a scenario where gas flows into the two arm regions along the bar, becomes condensed/shocked, forms YMCs, and then these YMCs heat and dissociate the gas., Comment: 19 pages, 7 figures and 2 tables in total (12 pages and 6 figures in main text). Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJL

Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ring-like structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ~0.3" (30pc) resolution new JWST/MIRI imaging with archival ALMA CO(2-1) mapping of the central ~5kpc of the nearby barred spiral galaxy NGC1365, to investigate the physical mechanisms responsible for this extreme star formation. The molecular gas morphology is resolved into two well-known bright bar lanes that surround a smooth dynamically cold gas disk (R_gal ~ 475pc) reminiscent of non-star-forming disks in early type galaxies and likely fed by gas inflow triggered by stellar feedback in the lanes. The lanes host a large number of JWST-identified massive young star clusters. We find some evidence for temporal star formation evolution along the ring. The complex kinematics in the gas lanes reveal strong streaming motions and may be consistent with convergence of gas streamlines expected there. Indeed, the extreme line-widths are found to be the result of inter-`cloud' motion between gas peaks; ScousePy decomposition reveals multiple components with line widths of ~ 19km/s and surface densities of ~ 800M_sun/pc^2, similar to the properties observed throughout the rest of the central molecular gas structure. Tailored hydro-dynamical simulations exhibit many of the observed properties and imply that the observed structures are transient and highly time-variable. From our study of NGC1365, we conclude that it is predominantly the high gas inflow triggered by the bar that is setting the star formation in its CMZ., Comment: 24 pages, 8 figures, accepted for publications as part of PHANGS-JWST ApJL Focus Issue

We study the young star cluster populations in 23 dwarf and irregular galaxies observed by the HST Legacy ExtraGalactic Ultraviolet Survey (LEGUS), and examine relationships between the ensemble properties of the cluster populations and those of their host galaxies: star formation rate (SFR) density ($\Sigma_{\rm SFR}$). A strength of this analysis is the availability of SFRs measured from temporally resolved star formation histories which provide the means to match cluster and host-galaxy properties on several timescales (1-10, 1-100, and 10-100~Myr). Nevertheless, studies of this kind are challenging for dwarf galaxies due to the small numbers of clusters in each system. We mitigate these issues by combining the clusters across different galaxies with similar $\Sigma_{\rm SFR}$ properties. We find good agreement with a well-established relationship ($M_{V}^{brightest}$-SFR), but find no significant correlations between $\Sigma_{\rm SFR}$ and the slopes of the cluster luminosity function, mass function, nor the age distribution. We also find no significant trend between the the fraction of stars in bound clusters at different age ranges ($\Gamma_{1-10}$, $\Gamma_{10-100}$, and $\Gamma_{1-100}$) and $\Sigma_{\rm SFR}$ of the host galaxy. Our data show a decrease in $\Gamma$ over time (from 1-10 to 10-100~Myr) suggesting early cluster dissolution, though the presence of unbound clusters in the youngest time bin makes it difficult to quantify the degree of dissolution. While our data do not exhibit strong correlations between $\Sigma_{\rm SFR}$ and ensemble cluster properties, we cannot rule out that a weak trend might exist given the relatively large uncertainties due to low number statistics and the limited $\Sigma_{\rm SFR}$ range probed., Comment: 30 pages, 15 figures, 7 tables, Accepted for publication in MNRAS

The PHANGS collaboration has been building a reference dataset for the multi-scale, multi-phase study of star formation and the interstellar medium 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 ($\sim$5-50 pc). In Cycle 1, PHANGS is conducting an 8-band imaging survey from 2-21$\mu$m of 19 nearby spiral galaxies. CO(2-1) mapping, optical integral field spectroscopy, 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 $\sim$10,000 star clusters, MUSE spectroscopic mapping of $\sim$20,000 HII regions, and $\sim$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., Comment: Re-submitted after addressing minor comments from referee. To be published as part of PHANGS-JWST ApJL Focus Issue

We use PHANGS-JWST data to identify and classify 1271 compact 21 $\mu$m sources in four nearby galaxies using MIRI F2100W data. We identify sources using a dendrogram-based algorithm, and we measure the background-subtracted flux densities for JWST bands from 2 $\mu$m to 21 $\mu$m. Using the SED in JWST as well as HST bands, plus ALMA and MUSE/VLT observations, we classify the sources by eye. Then we use this classification to define regions in color-color space, and so establish a quantitative framework for classifying sources. We identify 1085 sources as belonging to the ISM of the target galaxies with the remainder being dusty stars or background galaxies. These 21 $\mu$m sources are strongly spatially associated with HII regions ($>92\%$ of sources), while 74$\%$ of sources are coincident with a stellar association defined in the HST data. Using SED fitting, we find that the stellar masses of the 21 $\mu$m sources span a range of 10$^{2}$ to 10$^{4}~M_\odot$ with mass-weighted ages down to 2 Myr. There is a tight correlation between attenuation-corrected H$\alpha$ and 21 $\mu$m luminosity for $L_{\nu,\mathrm{F2100W}}>10^{19}~\mathrm{W~Hz}^{-1}$. Young embedded source candidates selected at 21 $\mu$m are found below this threshold and have $M_\star < 10^{3}~M_\odot$., Comment: 15 pages, 5 figures, Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

We present a high-resolution view of bubbles within The Phantom Galaxy (NGC 628); a nearby (~10Mpc), star-forming (~2Msun/yr), face-on (i~9deg) grand-design spiral galaxy. With new data obtained as part of the PHANGS-JWST treasury program, we perform a detailed case-study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (The Phantom Void; over 1kpc in diameter), and the other being a smaller region that may be the precursor to such a large bubble (The Precursor Phantom Void). When comparing to matched resolution Halpha observations from the Hubble Space Telescope (HST), we see that the ionized gas is brightest in the shells of both bubbles, and is coincident with the youngest (~1Myr) and most massive (~100,000Msun) stellar associations. We also find an older generation (~20Myr) of stellar associations is present within the bubble of The Phantom Void. From our kinematic analysis of the HI, H2 (CO) and HII gas across The Phantom Void, we infer a high expansion speed of around 15 to 50km/s. The large size and high expansion speed of The Phantom Void suggest that the driving mechanism is sustained stellar feedback due to multiple mechanisms, where early feedback first cleared a bubble (as we observe now in The Precursor Phantom Void), and since then SNe have been exploding within the cavity, and have accelerated the shell. Finally, comparison to simulations shows a striking resemblance to our JWST observations, and suggests that such large-scale stellar feedback-driven bubbles should be common within other galaxies., Comment: 12 pages total, 8 figures, and 1 table. Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

The first JWST observations of nearby galaxies have unveiled a rich population of bubbles that trace the stellar feedback mechanisms responsible for their creation. Studying these bubbles therefore allows us to chart the interaction between stellar feedback and the interstellar medium, and the larger galactic flows needed to regulate star formation processes globally. We present the first catalog of bubbles in NGC628, visually identified using MIRI F770W PHANGS-JWST observations, and use them to statistically evaluate bubble characteristics. We classify 1694 structures as bubbles with radii between 6-552 pc. Of these, 31% contain at least one smaller bubble at their edge, indicating that previous generations of star formation have a local impact on where new stars form. On large scales, most bubbles lie near a spiral arm, and their radii increase downstream compared to upstream. Furthermore, bubbles are elongated in a similar direction to the spiral arm ridge-line. These azimuthal trends demonstrate that star formation is intimately connected to the spiral arm passage. Finally, the bubble size distribution follows a power-law of index $p=-2.2\pm0.1$, which is slightly shallower than the theoretical value by 1-3.5$\sigma$ that did not include bubble mergers. The fraction of bubbles identified within the shells of larger bubbles suggests that bubble merging is a common process. Our analysis therefore allows us to quantify the number of star-forming regions that are influenced by an earlier generation, and the role feedback processes have in setting the global star formation rate. With the full PHANGS-JWST sample, we can do this for more galaxies., Comment: 25 pages total, 13 Figures and 1 Table. Accepted for publication in ApJL as part of a PHANGS-JWST First Results Focus issue

We present a comparison of theoretical predictions of dust continuum and polycyclic aromatic hydrocarbon (PAH) emission with new JWST observations in three nearby galaxies: NGC 628, NGC 1365, and NGC 7496. Our analysis focuses on a total of 1063 compact stellar clusters and 2654 stellar associations previously characterized by HST in the three galaxies. We find that the distributions and trends in the observed PAH-focused infrared colors generally agree with theoretical expectations, and that the bulk of the observations is more aligned with models of larger, ionized PAHs. These JWST data usher in a new era of probing interstellar dust and studying how the intense radiation fields near stellar clusters and associations play a role in shaping the physical properties of PAHs., Comment: 8 pages, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

We combine JWST observations with ALMA CO and VLT-MUSE H$\alpha$ data to examine off-spiral arm star formation in the face-on, grand-design spiral galaxy NGC 628. We focus on the northern spiral arm, around a galactocentric radius of 3-4 kpc, and study two spurs. These form an interesting contrast, as one is CO-rich and one CO-poor, and they have a maximum azimuthal offset in MIRI 21$\mu$m and MUSE H$\alpha$ of around 40$^\circ$ (CO-rich) and 55$^\circ$ (CO-poor) from the spiral arm. The star formation rate is higher in the regions of the spurs near to spiral arms, but the star formation efficiency appears relatively constant. Given the spiral pattern speed and rotation curve of this galaxy and assuming material exiting the arms undergoes purely circular motion, these offsets would be reached in 100-150 Myr, significantly longer than the 21$\mu$m and H$\alpha$ star formation timescales (both <10 Myr). The invariance of the star formation efficiency in the spurs versus the spiral arms indicates massive star formation is not only triggered in spiral arms, and cannot simply occur in the arms and then drift away from the wave pattern. These early JWST results show that in-situ star formation likely occurs in the spurs, and that the observed young stars are not simply the `leftovers' of stellar birth in the spiral arms. The excellent physical resolution and sensitivity that JWST can attain in nearby galaxies will well resolve individual star-forming regions and help us to better understand the earliest phases of star formation., Comment: 11 pages, 5 Figures, accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ. Updated to include missing author and published paper references

In the hierarchical view of star formation, giant molecular gas 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 ($\leq$ 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., Comment: 15 pages, 11 figures, 4 tables. Accepted to MNRAS Sept 6 2022

We present Hubble Space Telescope WFC3/UVIS (F275W, F336W) and ACS/WFC optical (F435W, F555W, and F814W) observations of the nearby grand-design spiral galaxy M101 as part of the Legacy Extragalactic UV Survey (LEGUS). Compact sources detected in at least four bands were classified by both human experts and the convolutional neural network StarcNet. Human experts classified the 2,351 brightest sources, retrieving $N_{c} = 965$ star clusters. StarcNet, trained on LEGUS data not including M101, classified all 4,725 sources detected in four bands, retrieving $N_{c} = 2,270$ star clusters. The combined catalog represents the most complete census to date of compact star clusters in M101. We find that for the 2,351 sources with both a visual- and ML-classification StarcNet is able to reproduce the human classifications at high levels of accuracy ($\sim 80-90\%$), which is equivalent to the level of agreement between human classifiers in LEGUS. The derived cluster age distribution implies a disruption rate of $dN/d\tau \propto \tau^{-0.45 \pm 0.14}$ over $10^{7} < \tau < 10^{8.5}$yr for cluster masses $\geq 10^{3.55} M_{\odot}$ for the central region of M101 and $dN/d\tau \propto \tau^{-0.02 \pm 0.15}$ for cluster masses $\geq 10^{3.38} M_{\odot}$ in the northwest region of the galaxy. The trends we recover are weaker than those of other nearby spirals (e.g. M51) and starbursts, consistent with the M101 environment having a lower-density interstellar medium, and providing evidence in favor of environmentally-dependent cluster disruption in the central, southeast, and northwest regions of M101., Comment: 24 pages, 18 figures, 3 tables. Accepted for publication in the Astrophysical Journal

The analysis of star cluster ages in tandem with the morphology of their HII regions can provide insight into the processes that clear a cluster's natal gas, as well as the accuracy of cluster ages and dust reddening derived from Spectral Energy Distribution (SED) fitting. We classify 3757 star clusters in 16 nearby galaxies according to their H$\alpha$ morphology (concentrated, partially exposed, no emission), using Hubble Space Telescope (HST) imaging from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We find: 1) The mean SED ages of clusters with concentrated (1-2 Myr) and partially exposed HII region morphologies (2-3 Myr) indicate a relatively early onset of gas clearing and a short (1-2 Myr) clearing timescale. 2) The reddening of clusters can be overestimated due to the presence of red supergiants, which is a result of stochastic sampling of the IMF in low mass clusters. 3) The age-reddening degeneracy impacts the results of the SED fitting - out of 1408 clusters with M$_*$ $\geq$ 5000 M$_{\odot}$, we find that at least 46 (3%) have SED ages which appear significantly underestimated or overestimated based on H$\alpha$ and their environment, while the total percentage of poor age estimates is expected to be several times larger. 4) Lastly, we examine the dependence of the morphological classifications on spatial resolution. At HST resolution, our conclusions are robust to the distance range spanned by the sample (3-10 Mpc). However, analysis of ground-based H$\alpha$ images shows that compact and partially exposed morphologies frequently cannot be distinguished from each other., Comment: 24 pages, 14 figures, 9 tables; accepted for publication in MNRAS

Using PHANGS-HST NUV-U-B-V-I imaging of 17 nearby spiral galaxies, we study samples of star clusters and stellar associations, visually selected to be bright and relatively isolated, for three purposes: to compute aperture corrections for star cluster photometry, to explore the utility of quantitative morphologies in the analysis of clusters and associations, and to compare to synthetic stellar population models. We provide a technical summary of our procedures to determine aperture corrections, a standard step in the production of star cluster candidate catalogs, and compare to prior work. We also use this specialized sample to launch an analysis into the measurement of star cluster light profiles. We focus on one measure, $M_{20}$ (normalized second order moment of the brightest 20% of pixels), applied previously to study the morphologies of galaxies. We find that $M_{20}$ in combination with UB-VI colors, yields a parameter space where distinct loci are formed by single-peaked symmetric clusters, single-peaked asymmetric clusters, and multi-peaked associations. We discuss the potential applications for using $M_{20}$ to gain insight into the formation and evolution of clusters and associations. Finally, we compare the color distributions of this sample with various synthetic stellar population models. One finding is that the standard procedure of using a single-metallicity SSP track to fit the entire population of clusters in a given galaxy should be revisited, as the oldest globular clusters will be more metal-poor compared to clusters formed recently., Comment: 20 figures, 4 tables, 23 pages. Accepted for publication in MNRAS

When completed, the PHANGS-HST project will provide a census of roughly 50,000 compact star clusters and associations, as well as human morphological classifications for roughly 20,000 of those objects. These large numbers motivated the development of a more objective and repeatable method to help perform source classifications. In this paper we consider the results for five PHANGS-HST galaxies (NGC 628, NGC 1433, NGC 1566, NGC 3351, NGC 3627) using classifications from two convolutional neural network architectures (RESNET and VGG) trained using deep transfer learning techniques. The results are compared to classifications performed by humans. The primary result is that the neural network classifications are comparable in quality to the human classifications with typical agreement around 70 to 80$\%$ for Class 1 clusters (symmetric, centrally concentrated) and 40 to 70$\%$ for Class 2 clusters (asymmetric, centrally concentrated). If Class 1 and 2 are considered together the agreement is 82 $\pm$ 3$\%$. Dependencies on magnitudes, crowding, and background surface brightness are examined. A detailed description of the criteria and methodology used for the human classifications is included along with an examination of systematic differences between PHANGS-HST and LEGUS. The distribution of data points in a colour-colour diagram is used as a 'figure of merit' to further test the relative performances of the different methods. The effects on science results (e.g., determinations of mass and age functions) of using different cluster classification methods are examined and found to be minimal., Comment: 28 pages, 25 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society. Version with full resolution figures found at https://sites.google.com/view/phangs/publications

We present an innovative and widely applicable approach for the detection and classification of stellar clusters, developed for the PHANGS-HST Treasury Program, an $NUV$-to-$I$ band imaging campaign of 38 spiral galaxies. Our pipeline first generates a unified master source list for stars and candidate clusters, to enable a self-consistent inventory of all star formation products. To distinguish cluster candidates from stars, we introduce the Multiple Concentration Index (MCI) parameter, and measure inner and outer MCIs to probe morphology in more detail than with a single, standard concentration index (CI). We improve upon cluster candidate selection, jointly basing our criteria on expectations for MCI derived from synthetic cluster populations and existing cluster catalogues, yielding model and semi-empirical selection regions (respectively). Selection purity (confirmed clusters versus candidates, assessed via human-based classification) is high (up to 70$\%$) for moderately luminous sources in the semi-empirical selection region, and somewhat lower overall (outside the region or fainter). The number of candidates rises steeply with decreasing luminosity, but pipeline-integrated Machine Learning (ML) classification prevents this from being problematic. We quantify the performance of our PHANGS-HST methods in comparison to LEGUS for a sample of four galaxies in common to both surveys, finding overall agreement with $50-75\%$ of human verified star clusters appearing in both catalogues, but also subtle differences attributable to specific choices adopted by each project. The PHANGS-HST ML-classified Class 1 or 2 catalogues reach $\sim1$ magnitude fainter, $\sim2\times$ lower stellar mass, and are $2{-}5\times$ larger in number, than attained in the human classified samples., Comment: 36 pages, 18 figures, updated to v2, accepted to MNRAS

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and 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 five band 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 multi-scale 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, HII 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, highlighting 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. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $\mu$m imaging has been approved., Comment: ApJS in press. Accepted draft includes moderate revision and updates to the originally posted submitted version. 27 pages, 11 figures, 1 table. Data products released at https://archive.stsci.edu/hlsp/phangs-hst/

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 (LEGUS). 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_{\odot}$) 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 analyze the resulting cluster distributions., Comment: Accepted for publication in MNRAS (2021-01-06), 22 pages, 18 figures, 2 tables

We carry out a comparative analysis of the metallicities from the stellar, neutral-gas, and ionized-gas components in the metal-rich spiral galaxy M83. We analyze spectroscopic observations taken with the Hubble Space Telescope (HST), the Large Binocular Telescope (LBT) and the Very Large Telescope (VLT). We detect a clear depletion of the HI gas, as observed from the HI column densities in the nuclear region of this spiral galaxy. We find column densities of log[$N$(HI) cm$^{-2}$] $<$ 20.0 at galactocentric distances of $<$ 0.18 kpc, in contrast to column densities of log[$N$(HI) cm$^{-2}$] $\sim$ 21.0 in the galactic disk, a trend observed in other nearby spiral galaxies. We measure a metallicity gradient of $-$0.03 $\pm$ 0.01 dex kpc$^{-1}$ for the ionized gas, comparable to the metallicity gradient of a local benchmark of 49 nearby star-forming galaxies of $-$0.026 $\pm$ 0.002 dex kpc$^{-1}$. Our co-spatial metallicity comparison of the multi-phase gas and stellar populations shows excellent agreement outside of the nucleus of the galaxy hinting at a scenario where the mixing of newly synthesized metals from the most massive stars in the star clusters takes longer than their lifetimes ($\sim$10 Myr). Finally, our work shows that caution must be taken when studying the metallicity gradient of the neutral-gas component in star-forming galaxies, since this can be strongly biased, as these environments can be dominated by molecular gas. In these regions the typical metallicity tracers can provide inaccurate abundances as they may trace both the neutral- and molecular-gas components., Comment: 29 pages, 14 tables, 9 figures, Accepted for Publication in ApJ

The Clusters, Clumps, Dust, and Gas in Extreme Star-Forming Galaxies (CCDG) survey with the Hubble Space Telescope includes multi-wavelength imaging of 13 galaxies less than 100 Mpc away spanning a range of morphologies and sizes, from Blue Compact Dwarfs (BCDs) to luminous infrared galaxies (LIRGs), all with star formation rates in excess of hundreds of solar masses per year. Images of 7 merging galaxies in the CCDG survey were artificially redshifted to compare with galaxies at z=0.5, 1, and 2. Most redshifted tails have surface brightnesses that would be visible at z=0.5 or 1 but not at z=2 due to cosmological dimming. Giant star clumps are apparent in these galaxies; the 51 measured have similar sizes, masses and colors as clumps in observed high-z systems in UDF, GEMS, GOODS, and CANDELS surveys. These results suggest that some clumpy high-z galaxies without observable tidal features could be the result of mergers. The local clumps also have the same star formation rate per unit area and stellar surface density as clumps observed at intermediate and high redshift, so they provide insight into the substructure of distant clumps. A total of 1596 star clusters brighter than MV = -9 were identified within the boundaries of the local clumps. The cluster magnitude distribution function is a power law with approximately the same slope (approximately -1 for a number-log luminosity plot) for all the galaxies both inside and outside the clumps and independent of clump surface brightness., Comment: 32 pages, 17 figures, 2 tables

We present a new catalog and results for the cluster system of the starburst galaxy NGC 4449 based on multi-band imaging observations taken as part of the LEGUS and Halpha-LEGUS surveys. We improve the spectral energy fitting method used to estimate cluster ages and find that the results, particularly for older clusters, are in better agreement with those from spectroscopy. The inclusion of Halpha measurements, the role of stochasticity for low mass clusters, the assumptions about reddening, and the choices of SSP model and metallicity all have important impacts on the age-dating of clusters. A comparison with ages derived from stellar color-magnitude diagrams for partially resolved clusters shows reasonable agreement, but large scatter in some cases. The fraction of light found in clusters relative to the total light (i.e., T_L) in the U, B, and V filters in 25 different ~kpc-size regions throughout NGC 4449 correlates with both the specific Region Luminosity, R_L, and the dominant age of the underlying stellar population in each region. The observed cluster age distribution is found to decline over time as dN/dt ~ t^g, with g=-0.85+/-0.15, independent of cluster mass, and is consistent with strong, early cluster disruption. The mass functions of the clusters can be described by a power law with dN/dM ~ M^b and b=-1.86+/-0.2, independent of cluster age. The mass and age distributions are quite resilient to differences in age-dating methods. There is tentative evidence for a factor of 2-3 enhancement in both the star and cluster formation rate ~100 - 300 Myr ago, indicating that cluster formation tracks star formation generally. The enhancement is probably associated with an earlier interaction event.

The morphology of HII regions around young star clusters provides insight into the timescales and physical processes that clear a cluster's natal gas. We study ~700 young clusters (<10Myr) in three nearby spiral galaxies (NGC 7793, NGC 4395, and NGC 1313) using Hubble Space Telescope (HST) imaging from LEGUS (Legacy ExtraGalactic Ultraviolet Survey). Clusters are classified by their H$\alpha$ morphology (concentrated, partially exposed, no-emission) and whether they have neighboring clusters (which could affect the clearing timescales). Through visual inspection of the HST images, and analysis of ages, reddenings, and stellar masses from spectral energy distributions fitting, together with the (U-B), (V-I) colors, we find: 1) the median ages indicate a progression from concentrated (~3 Myr), to partially exposed (~4 Myr), to no H$\alpha$ emission (>5Myr), consistent with the expected temporal evolution of HII regions and previous results. However, 2) similarities in the age distributions for clusters with concentrated and partially exposed H$\alpha$ morphologies imply a short timescale for gas clearing (<1Myr). 3) our cluster sample's median mass is ~1000 M, and a significant fraction (~20%) contain one or more bright red sources (presumably supergiants), which can mimic reddening effects. Finally, 4) the median E(B-V) values for clusters with concentrated H$\alpha$ and those without H$\alpha$ emission appear to be more similar than expected (~0.18 vs. ~0.14, respectively), but when accounting for stochastic effects, clusters without H$\alpha$ emission are less reddened. To mitigate stochastic effects, we experiment with synthesizing more massive clusters by stacking fluxes of clusters within each H$\alpha$ morphological class. Composite isolated clusters also reveal a color and age progression for H$\alpha$ morphological classes, consistent with analysis of the individual clusters., Comment: 20 pages, 16 figures, 7 tables; accepted for publication in MNRAS

The Hubble Space Telescope (HST) has obtained multi-epoch observations providing the opportunity for a comprehensive variability search aiming to uncover new variables. We have therefore undertaken the task of creating a catalog of variable sources based on version 3 of the Hubble Source Catalog (HSC), which relies on publicly available images obtained with the WFPC2, ACS, and WFC3 instruments onboard the HST. We adopted magnitude-dependent thresholding in median absolute deviation (a robust measure of light curve scatter) combined with sophisticated preprocessing techniques and visual quality control to identify and validate variable sources observed by Hubble with the same instrument and filter combination five or more times. The Hubble Catalog of Variables (HCV) includes 84,428 candidate variable sources (out of 3.7 million HSC sources that were searched for variability) with $V \leq 27$ mag; for 11,115 of them the variability is detected in more than one filter. The data points in the light curves of the variables in the HCV catalog range from five to 120 points (typically having less than ten points); the time baseline ranges from under a day to over 15 years; while $\sim$8% of all variables have amplitudes in excess of 1 mag. Visual inspection performed on a subset of the candidate variables suggests that at least 80% of the candidate variables that passed our automated quality control are true variable sources rather than spurious detections resulting from blending, residual cosmic rays, and calibration errors. The HCV is the first, homogeneous catalog of variable sources created from archival HST data and currently is the deepest catalog of variables available. The catalog includes variable stars in our Galaxy and nearby galaxies, as well as transients and variable active galactic nuclei. (abbreviated), Comment: 33 pages, 15 figures, 11 tables; accepted in A&A. The HCV is available from the ESA Hubble Science Archive (eHST) (http://archives.esac.esa.int/ehst) at ESAC, and can be easily explored using the HCV Explorer (http://archives.esac.esa.int/hcv-explorer). The HCV is also hosted at STScI in the form of a High Level Science Product (HLSP) in MAST (https://archive.stsci.edu/hlsp/hcv)

We present the results of a proof-of-concept experiment which demonstrates that deep learning can successfully be used for production-scale classification of compact star clusters detected in HST UV-optical imaging of nearby spiral galaxies (D < 20 Mpc) in the PHANGS-HST survey. Given the relatively small nature of existing, human-labelled star cluster samples, we transfer the knowledge of state-of-the-art neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We perform a series of experiments to determine the dependence of classification performance on: neural network architecture (ResNet18 and VGG19-BN); training data sets curated by either a single expert or three astronomers; and the size of the images used for training. We find that the overall classification accuracies are not significantly affected by these choices. The networks are used to classify star cluster candidates in the PHANGS-HST galaxy NGC 1559, which was not included in the training samples. The resulting prediction accuracies are 70%, 40%, 40-50%, 50-70% for class 1, 2, 3 star clusters, and class 4 non-clusters respectively. This performance is competitive with consistency achieved in previously published human and automated quantitative classification of star cluster candidate samples (70-80%, 40-50%, 40-50%, and 60-70%). The methods introduced herein lay the foundations to automate classification for star clusters at scale, and exhibit the need to prepare a standardized dataset of human-labelled star cluster classifications, agreed upon by a full range of experts in the field, to further improve the performance of the networks introduced in this study., Comment: 17 pages, 7 figures, 3 appendices

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses., Comment: 25 pages; 28 figures (Accepted to MNRAS)

Stellar populations are powerful tools for investigating the evolution of extragalactic environments. We present the first UV integrated-light spectroscopic observations for 15 young star clusters in the starburst M83 with a special focus on metallicity measurements. The data were obtained with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope. We analyse the data applying an abundance technique previously used to study an optical set of star clusters. We estimate a central metallicity of [Z] = $+$0.20 $\pm$ 0.15 dex in agreement with those obtained through independent methods, i.e. $J$-band and blue supergiants. We estimate a UV metallicity gradient of $-$0.041 $\pm$ 0.022 dex kpc$^{-1}$ consistent with the optical metallicity gradient of $-$0.040 $\pm$ 0.032 dex kpc$^{-1}$ for $R/R_{25}<0.5$. Combining our stellar metallicities, UV and optical, with those from HII regions (strong-line abundances based on empirical calibrations) we identify two possible breaks in the gradient of M83 at galactocentric distances of $R\sim0.5$ and $1.0\:R_{25}$. If the abundance breaks are genuine, the metallicity gradient of this galaxy follows a steep-shallow-steep trend, a scenario predicted by three-dimensional (3D) numerical simulations of disc galaxies. The first break is located near the corotation radius. This first steep gradient may have originated by recent star formation episodes and a relatively young bar ($<$1 Gyr). In the numerical simulations the shallow gradient is created by the effects of dilution by outflow where low-metallicity material is mixed with enriched gas. And finally, the second break and last steep gradient mark the farthest galactocentric distances where the outward flow has penetrated., Comment: 18 pages, 7 figures, Accepted for Publication in ApJ

We present a study correlating the spatial locations of young star clusters with those of molecular clouds in NGC~5194, in order to investigate the timescale over which clusters separate from their birth clouds. The star cluster catalogues are from the Legacy ExtraGalactic UV Survey (LEGUS) and the molecular clouds from the Plateau de Bure Interefrometer Arcsecond Whirpool Survey (PAWS). We find that younger star clusters are spatially closer to molecular clouds than older star clusters. The median ages for clusters associated with clouds is 4~Myr whereas it is 50~Myr for clusters that are sufficiently separated from a molecular cloud to be considered unassociated. After $\sim$6~Myr, the majority of the star clusters lose association with their molecular gas. Younger star clusters are also preferentially located in stellar spiral arms where they are hierarchically distributed in kpc-size regions for 50-100~Myr before dispersing. The youngest star clusters are more strongly clustered, yielding a two-point correlation function with $\alpha=-0.28\pm0.04$, than the GMCs ($\alpha=-0.09\pm0.03$) within the same PAWS field. However, the clustering strength of the most massive GMCs, supposedly the progenitors of the young clusters for a star formation efficiency of a few percent, is comparable ($\alpha=-0.35\pm0.05$) to that of the clusters. We find a galactocentric-dependence for the coherence of star formation, in which clusters located in the inner region of the galaxy reside in smaller star-forming complexes and display more homogeneous distributions than clusters further from the centre. This result suggests a correlation between the survival of a cluster complex and its environment., Comment: 17 pages, 12 figures. Accepted for publication in MNRAS

We have explored the role environmental factors play in determining characteristics of young stellar objects in nearby dwarf irregular and Blue Compact Dwarf galaxies. Star clusters are characterized by concentrations, masses, and formation rates, OB associations by mass and mass surface density, O stars by their numbers and near-ultraviolet absolute magnitudes, and HII regions by Halpha surface brightnesses. These characteristics are compared to surrounding galactic pressure, stellar mass density, HI surface density, and star formation rate surface density. We find no trend of cluster characteristics with environmental properties, implying that larger scale effects are more important in determining cluster characteristics or that rapid dynamical evolution erases memory of the initial conditions. On the other hand, the most massive OB associations are found at higher pressure and HI surface density, and there is a trend of higher HII region Halpha surface brightness with higher pressure, suggesting that a higher concentration of massive stars and gas are found preferentially in regions of higher pressure. At low pressures we find massive stars but not bound clusters and OB associations. We do not find evidence for an increase of cluster formation efficiency as a function of star formation rate density. However, there is an increase in the ratio of the number of clusters to number of O stars with pressure, perhaps reflecting an increase in clustering properties with star formation rate., Comment: 32 pages, 28 figures, published in 2018, AJ, 156, 21

We estimate the fraction of stars that form in compact clusters (bound and unbound), Gamma_F, in a diverse sample of eight star-forming galaxies, including two irregulars, two dwarf starbursts, two spirals, and two mergers. The average value for our sample is Gamma_F ~ 24 +/- 9%. We also calculate the fraction of stars in clusters that survive to ages between t1 and t2, denoted by Gamma_S(t1,t2), and find Gamma_S(10,100)=4.6 +/- 2.5% and Gamma_S(100,400)=2.4 +/- 1.1 %, significantly lower than Gamma_F for the same galaxies. We do not find any systematic trends in Gamma_F or Gamma_S with the star formation rate (SFR), the SFR per unit area (Sigma_SFR), or the surface density of molecular gas (Sigma_H2) within the host galaxy. Our results are consistent with those found previously from the CMF/SFR statistic (where CMF is the cluster mass function), and with the quasi-universal model in which clusters in different galaxies form and disrupt in similar ways. Our results, however, contradict many previous claims that the fraction of stars in bound clusters increases strongly with Sigma_SFR and Sigma_H_2. We find that the previously reported trends are largely driven by comparisons that mixed Gamma_F ~ Gamma_S(0,10) and Gamma_S(10,100), where Gamma_S(0,10) was systematically used for galaxies with higher Sigma_SFR and Sigma_H2, and Gamma_S(10,100) for galaxies with lower Sigma_SFR and Sigma_H_2., Comment: Published in ApJ, posting a version so everyone has access

We use Atacama Large Millimeter Array CO(3-2) observations in conjunction with optical observations from the Hubble Space Telescope to determine the ratio of stellar to gas mass for regions in the Antennae Galaxies. We adopt the term "instantaneous mass ratio" IMR(t) = M$_{stars}$/(M$_{gas}$ +M$_{stars}$), that is equivalent to the star formation efficiency for an idealized system at t = 0. We use two complementary approaches to determining the IMR(t) based on 1) the enclosed stellar and molecular mass within circular apertures centered on optically-identified clusters, and 2) a tessellation algorithm that defines regions based on CO emission. We find that only a small number of clusters appear to have IMR(0) = SFE > 0.2, which suggests that only a small fraction of these clusters will remain bound. The results suggest that by ages of $10^{6.7}$ years, some clusters will have lost all of their associated molecular gas, and by $10^{7.5}$ years this is true for the majority of clusters. There appears to be slight dependence of the IMR(t) on the CO surface brightness, which could support the idea that dense molecular environments are more likely to form bound clusters. However, the IMR(t) appears to have a strong dependence on extinction, which likely traces the evolutionary state of clusters., Comment: 17 pages, 10 figures, accepted for publication in ApJ

It has recently been established that the properties of young star clusters (YSCs) can vary as a function of the galactic environment in which they are found. We use the cluster catalogue produced by the Legacy Extragalactic UV Survey (LEGUS) collaboration to investigate cluster properties in the spiral galaxy M51. We analyse the cluster population as a function of galactocentric distance and in arm and inter-arm regions. The cluster mass function exhibits a similar shape at all radial bins, described by a power law with a slope close to $-2$ and an exponential truncation around $10^5\ \rm{M}_{\odot}$ . While the mass functions of the YSCs in the spiral arm and inter-arm regions have similar truncation masses, the inter-arm region mass function has a significantly steeper slope than the one in the arm region; a trend that is also observed in the giant molecular cloud mass function and predicted by simulations. The age distribution of clusters is dependent on the region considered, and is consistent with rapid disruption only in dense regions, while little disruption is observed at large galactocentric distances and in the inter-arm region. The fraction of stars forming in clusters does not show radial variations, despite the drop in the $H_2$ surface density measured as function of galactocentric distance. We suggest that the higher disruption rate observed in the inner part of the galaxy is likely at the origin of the observed flat cluster formation efficiency radial profile., Comment: 22 pages, 18 figures, 9 tables (+appendix: 6 pages, 3 figures, 4 tables). Accepted for publication in MNRAS

We present a detailed study of the Magellanic irregular galaxy NGC 4449 based on both archival and new photometric data from the Legacy Extragalactic UV Survey, obtained with the Hubble Space Telescope Advanced Camera for Surveys and Wide Field Camera 3. Thanks to its proximity ($D=3.82\pm 0.27$ Mpc) we reach stars 3 magnitudes fainter than the tip of the red giant branch in the F814W filter. The recovered star formation history spans the whole Hubble time, but due to the age-metallicity degeneracy of the red giant branch stars, it is robust only over the lookback time reached by our photometry, i.e. $\sim 3$ Gyr. The most recent peak of star formation is around 10 Myr ago. The average surface density star formation rate over the whole galaxy lifetime is $0.01$ M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$. From our study it emerges that NGC 4449 has experienced a fairly continuous star formation regime in the last 1 Gyr with peaks and dips whose star formation rates differ only by a factor of a few. The very complex and disturbed morphology of NGC 4449 makes it an interesting galaxy for studies of the relationship between interactions and starbursts, and our detailed and spatially resolved analysis of its star formation history does indeed provide some hints on the connection between these two phenomena in this peculiar dwarf galaxy., Comment: 16 pages, 25 figures. Accepted for publication in the Astrophysical Journal