ABSTRACT We study the young star cluster populations in 23 dwarf and irregular galaxies observed by the Hubble Space Telescope (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 (ΣSFR). A strength of this analysis is the availability of SFRs measured from temporally resolved star formation histories that provide the means to match cluster and host galaxy properties on several time-scales (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 ΣSFR properties. We find good agreement with a well-established relationship ($M_{V}^{\mathrm{ brightest}}$–SFR), but find no significant correlations between ΣSFR and the slopes of the cluster luminosity function, mass function, nor the age distribution. We also find no significant trend between the fraction of stars in bound clusters at different age ranges (Γ1–10, Γ10–100, and Γ1–100) and ΣSFR of the host galaxy. Our data show a decrease in Γ 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 Σ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 ΣSFR range probed.

We present results from a joint ALMA-HST study of the nearby spiral galaxy NGC 628. We combine the HST LEGUS database of over 1000 stellar clusters in NGC 628 with ALMA Cycle 4 millimeter/submillimeter observations of the cold dust continuum that span ~15 square kpc including the nuclear region and western portions of the galaxy's disk. The resolution -- 1.1" or approximately 50 pc at the distance of NGC 628 -- allows us to constrain the spatial variations in the slope of the millimeter dust continuum as a function of the ages and masses of the nearby stellar clusters. Our results indicate an excess of dust emission in the millimeter assuming a typical cold dust model for a normal star-forming galaxy, but little correlation of the dust continuum slope with stellar cluster age or mass. For the depth and spatial coverage of these observations, we cannot substantiate the millimeter/submillimeter excess arising from the processing of dust grains by the local interstellar radiation field. We detect a bright unknown source in NGC 628 in ALMA bands 4 and 7 with no counterparts at other wavelengths from ancillary data. We speculate this is possibly a dust obscured supernova.

ABSTRACT The morphology of H ii regions around young star clusters provides insight into the time-scales and physical processes that clear a cluster’s natal gas. We study ∼700 young clusters (≤10 Myr) 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 α morphology (concentrated, partially exposed, no-emission) and whether they have neighbouring clusters (which could affect the clearing time-scales). 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) colours, we find (1) the median ages indicate a progression from concentrated (∼3 Myr), to partially exposed (∼4 Myr), to no H α emission (>5 Myr), consistent with the expected temporal evolution of H ii regions and previous results. However, (2) similarities in the age distributions for clusters with concentrated and partially exposed H α morphologies imply a short time-scale for gas clearing (≲1 Myr). Also, (3) our cluster sample’s median mass is ∼1000 M⊙, and a significant fraction ($\sim 20{{\ \rm per\ cent}}$) 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 α and those without H α emission appear to be more similar than expected (∼0.18 versus ∼0.14, respectively), but when accounting for stochastic effects, clusters without H α emission are less reddened. To mitigate stochastic effects, we experiment with synthesizing more massive clusters by stacking fluxes of clusters within each H α morphological class. Composite isolated clusters also reveal a colour and age progression for H α morphological classes, consistent with analysis of the individual clusters.

We detected 10 compact galaxy groups (CGs) at z = 0 in the semianalytic galaxy catalog of Guo et al. for the milli-Millennium Cosmological Simulation (sCGs in mGuo2010a). We aimed to identify potential canonical pathways for compact group evolution and thus illuminate the history of observed nearby CGs. By constructing merger trees for z = 0 sCG galaxies, we studied the cosmological evolution of key properties and compared them with z = 0 Hickson CGs (HCGs). We found that, once sCG galaxies come within 1 (0.5) Mpc of their most massive galaxy, they remain within that distance until z = 0, suggesting sCG "birth redshifts." At z = 0 stellar masses of sCG most massive galaxies are within 1010 ≲ M*/M⊙ ≲ 1011. In several cases, especially in the two four- and five-member systems, the amount of cold gas mass anticorrelates with stellar mass, which in turn correlates with hot gas mass. We define the angular difference between group members' 3D velocity vectors, Δθvel, and note that many of the groups are long-lived because their small values of Δθvel indicate a significant parallel component. For triplets in particular, Δθvel values range between 20° and 40° so that galaxies are coming together along roughly parallel paths, and pairwise separations do not show large pronounced changes after close encounters. The best agreement between sCG and HCG physical properties is for M* galaxy values, but HCG values are higher overall, including for star formation rates (SFRs). Unlike HCGs, due to a tail at low SFR and M* and a lack of M* ≳ 1011 M⊙ galaxies, only a few sCG galaxies are on the star-forming main sequence.

We derive the recent star formation histories of 23 active dwarf galaxies using HST observations from the Legacy Extragalactic UV Survey (LEGUS). We apply a color-magnitude diagram fitting technique using two independent sets of stellar models, PARSEC-COLIBRI and MIST. Despite the non-negligible recent activity, none of the 23 star forming dwarfs show enhancements in the last 100 Myr larger than three times the 100-Myr-average. The unweighted mean of the individual SFHs in the last 100 Myr is also consistent with a rather constant activity, irrespective of the atomic gas fraction. We confirm previous results that for dwarf galaxies the CMD-based average star formation rates (SFRs) are generally higher than the FUV-based SFR. For half of the sample, the 60-Myr-average CMD-based SFR is more than two times the FUV SFR. In contrast, we find remarkable agreement between the 10-Myr-average CMD-based SFR and the H${\alpha}$-based SFR. Finally, using core helium burning stars of intermediate mass we study the pattern of star formation spatial progression over the past 60 Myr, and speculate on the possible triggers and connections of the star formation activity with the environment in which these galaxies live. Approximately half of our galaxies show spatial progression of star formation in the last 60 Myr, and/or very recent diffuse and off-center activity compared to RGB stars.

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.

We present an investigation of the relationship between giant molecular cloud (GMC) properties and the associated stellar clusters in the nearby flocculent galaxy NGC 7793. We combine the star cluster catalogue from the HST LEGUS (Legacy ExtraGalactic UV Survey) programme with the 15 pc resolution ALMA CO(2–1) observations. We find a strong spatial correlation between young star clusters and GMCs such that all clusters still associated with a GMC are younger than 11 Myr and display a median age of 2 Myr. The age distribution increases gradually as the cluster–GMC distance increases, with star clusters that are spatially unassociated with molecular gas exhibiting a median age of 7 Myr. Thus, star clusters are able to emerge from their natal clouds long before the time-scale required for clouds to disperse. To investigate if the hierarchy observed in the stellar components is inherited from the GMCs, we quantify the amount of clustering in the spatial distributions of the components and find that the star clusters have a fractal dimension slope of −0.35 ± 0.03, significantly more clustered than the molecular cloud hierarchy with slope of −0.18 ± 0.04 over the range 40–800 pc. We find, however, that the spatial clustering becomes comparable in strength for GMCs and star clusters with slopes of −0.44 ± 0.03 and −0.45 ± 0.06, respectively, when we compare massive (>10^5 M_⊙) GMCs to massive and young star clusters. This shows that massive star clusters trace the same hierarchy as their parent GMCs, under the assumption that the star formation efficiency is a few per cent.

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 ExtragalacticUVSurvey (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) M-circle dot. 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 a 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.

We present a study of the dust-to-gas ratios in five nearby galaxies: NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broadband WFC3/UVIS UV and optical images from the Treasury program Legacy ExtraGalactic UV Survey (LEGUS) combined with archival HST/Advanced Camera for Surveys data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques and correlate these maps with neutral H i and CO gas maps from the literature, including the H i Nearby Galaxy Survey and the HERA CO-Line Extragalactic Survey. We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power-law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H2 is not included. This implies that underestimation of N H 2 in low-metallicity dwarfs from a too-low CO-to-H2 conversion factor XCO could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool for NGC 7793 and find systematically lower extinctions from SED fitting as compared to isochrone matching.

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in 5 bands from the near-UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near-UV color–magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color–magnitude diagrams to identify stars more massive than 14 M⊙, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the timescale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems.

We investigate the relation between gas and star formation in sub-galactic regions, 360 pc to 1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millimeter Telescope, with archival infrared images in the range 8-500 micron from the Spitzer Space Telescope and the Herschel Space Observatory. We extend the dynamic range of our mm (and dust) maps at the faint end, using a correlation between the far-infrared/millimeter colors F(70)/F(1100) [and F(160)/F(1100)] and the mid-infrared color F(8)/F(24) that we establish for the first time for this and other galaxies. Supplementing our data with maps of the extinction-corrected star formation rate (SFR) surface density, we measure both the SFR-molecular gas and the SFR-total gas relations in NGC4449. We find that the SFR-molecular gas relation is described by a power law with exponent that decreases from 1.5 to 1.2 for increasing region size, while the exponent of the SFR-total gas relation remains constant with value 1.5 independent of region size. We attribute the molecular law behavior to the increasingly better sampling of the molecular cloud mass function at larger region sizes; conversely, the total gas law behavior likely results from the balance between the atomic and molecular gas phases achieved in regions of active star formation. Our results indicate a non-linear relation between SFR and gas surface density in NGC4449, similar to what is observed for galaxy samples.

We use Hubble Space Telescope observations from the Legacy Extragalactic UV Survey to reconstruct the recent star formation histories (SFHs) of three actively star-forming dwarf galaxies, NGC 4449, Holmberg II, and NGC 1705, from their UV color–magnitude diagrams (CMDs). We apply a CMD fitting technique using two independent sets of stellar isochrones, PARSEC-COLIBRI and MIST, to assess the uncertainties related to stellar evolution modeling. Irrespective of the adopted stellar models, all three dwarfs are found to have had almost constant star formation rates (SFRs) in the last 100–200 Myr, with modest enhancements (a factor of ∼2) above the 100 Myr averaged SFR. Significant differences among the three dwarfs are found in terms of the overall SFR, the timing of the most recent peak, and the SFR/area. The initial mass function of NGC 1705 and Holmberg II is consistent with a Salpeter slope down to ≈5 M⊙, whereas it is slightly flatter, s = −2.0, in NGC 4449. The SFHs derived with the two different sets of stellar models are consistent with each other, except for some quantitative details, attributable to their input assumptions. They also share the drawback that all synthetic diagrams predict a clear separation in color between the upper main-sequence and helium-burning stars, which is not apparent in the data. Since neither differential reddening, which is significant in NGC 4449, nor unresolved binaries appear to be sufficient to fill the gap, we suggest this calls for a revision of both sets of stellar evolutionary tracks.

We report the large effort that is producing comprehensive high-level young star cluster (YSC) catalogs for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC 628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC 628 is consistent with a power-law distribution of slopes similar to-2 and a truncation of a few times 10(5) M-circle dot. After their formation, YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer time frame, confirming their being potentially bound systems. Associations disappear on timescales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find massindependent cluster disruption in the inner region of NGC 628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass (<= 10(4) M-circle dot) clusters, suggesting that a massdependent component is necessary to fully describe the YSC disruption process in NGC 628.

We present a study of the hierarchical clustering of the young stellar clusters in six local (3–15 Mpc) star-forming galaxies using Hubble Space Telescope broadband WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey). We identified 3685 likely clusters and associations, each visually classified by their morphology, and we use the angular two-point correlation function to study the clustering of these stellar systems. We find that the spatial distribution of the young clusters and associations are clustered with respect to each other, forming large, unbound hierarchical star-forming complexes that are in general very young. The strength of the clustering decreases with increasing age of the star clusters and stellar associations, becoming more homogeneously distributed after ∼40–60 Myr and on scales larger than a few hundred parsecs. In all galaxies, the associations exhibit a global behavior that is distinct and more strongly correlated from compact clusters. Thus, populations of clusters are more evolved than associations in terms of their spatial distribution, traveling significantly from their birth site within a few tens of Myr, whereas associations show evidence of disruption occurring very quickly after their formation. The clustering of the stellar systems resembles that of a turbulent interstellar medium that drives the star formation process, correlating the components in unbound star-forming complexes in a hierarchical manner, dispersing shortly after formation, suggestive of a single, continuous mode of star formation across all galaxies.

Interactions between pairs of isolated dwarf galaxies provide a critical window into low-mass hierarchical, gas-dominated galaxy assembly and the build-up of stellar mass in low-metallicity systems. We present the first Very Large Telescope/Multi Unit Spectroscopic Explorer (VLT/MUSE) optical integral field unit (IFU) observations of the interacting dwarf pair dm1647+21 selected from the TiNy Titans survey. The H α emission is widespread and corresponds to a total unobscured star formation rate (SFR) of 0.44 M ⊙ yr −1 , which is 2.7 times higher than the SFR inferred from Sloan Digital Sky Survey (SDSS) data. The implied specific SFR (sSFR) for the system is elevated by more than an order of magnitude above non-interacting dwarfs in the same mass range. This increase is dominated by the lower-mass galaxy, which has a sSFR enhancement of >50. Examining the spatially resolved maps of classic optical line diagnostics, we find that the interstellar medium (ISM) excitation can be fully explained by star formation. The velocity field of the ionized gas is not consistent with simple rotation. Dynamical simulations indicate that the irregular velocity field and the stellar structure is consistent with the identification of this system as an ongoing interaction between two dwarf galaxies. The widespread, clumpy enhancements in the star formation in this system point to important differences in the effect of mergers on dwarf galaxies, compared to massive galaxies; rather than the funneling of gas to the nucleus and giving rise to a nuclear starburst, starbursts in low-mass galaxy mergers may be triggered by large-scale ISM compression, and thus may be more distributed.

We present the implementation of a Bayesian formalism within the Stochastically Lighting Up Galaxies (SLUG) stellar population synthesis code, which is designed to investigate variations in the initial mass function (IMF) of star clusters. By comparing observed cluster photometry to large libraries of clusters simulated with a continuously varying IMF, our formalism yields the posterior probability distribution function (PDF) of the cluster mass, age and extinction, jointly with the parameters describing the IMF. We apply this formalism to a sample of star clusters from the nearby galaxy NGC 628, for which broad-band photometry in five filters is available as part of the Legacy ExtraGalactic UV Survey (LEGUS). After allowing the upper-end slope of the IMF (a3) to vary, we recover PDFs for the mass, age and extinction that are broadly consistent with what is found when assuming an invariant Kroupa IMF. However, the posterior PDF for a3 is very broad due to a strong degeneracy with the cluster mass, and it is found to be sensitive to the choice of priors, particularly on the cluster mass. We find only a modest improvement in the constraining power of a3 when adding Ha photometry from the companion Ha-LEGUS survey. Conversely, Ha photometry significantly improves the age determination, reducing the frequency of multi-modal PDFs. With the aid of mock clusters, we quantify the degeneracy between physical parameters, showing how constraints on the cluster mass that are independent of photometry can be used to pin down the IMF properties of star clusters.

We investigate the properties and evolution of star particles in two simulations of isolated spiral galaxies, and two galaxies from cosmological simulations. Unlike previous numerical work, where typically each star particle represents one 'cluster', for the isolated galaxies we are able to model features we term 'clusters' with groups of particles. We compute the spatial distribution of stars with different ages, and cluster mass distributions, comparing our findings with observations including the recent LEGUS survey. We find that spiral structure tends to be present in older (100s Myr) stars and clusters in the simulations compared to the observations. This likely reflects differences in the numbers of stars or clusters, the strength of spiral arms, and whether the clusters are allowed to evolve. Where we model clusters with multiple particles, we are able to study their evolution. The evolution of simulated clusters tends to follow that of their natal gas clouds. Massive, dense, long-lived clouds host massive clusters, whilst short-lived clouds host smaller clusters which readily disperse. Most clusters appear to disperse fairly quickly, in basic agreement with observational findings. We note that embedded clusters may be less inclined to disperse in simulations in a galactic environment with continuous accretion of gas on to the clouds than isolated clouds and correspondingly, massive young clusters which are no longer associated with gas tend not to occur in the simulations. Caveats of our models include that the cluster densities are lower than realistic clusters, and the simplistic implementation of stellar feedback.

IC342 is a nearby, late-type spiral galaxy with a young nuclear star cluster surrounded by several giant molecular clouds. The IC342 nuclear region is similar to the Milky Way and therefore provides an interesting comparison. We explore star formation in the nucleus using radio recombination line (RRL) and continuum emission at 5, 6.7, 33, and 35 GHz with the JVLA. These radio tracers are largely unaffected by dust and therefore sensitive to all of the thermal emission from the ionized gas produced by early-type stars. We resolve two components in the RRL and continuum emission within the nuclear region that lie east and west of the central star cluster. These components are associated both spatially and kinematically with two giant molecular clouds. We model these regions in two ways: a simple model consisting of uniform gas radiating in spontaneous emission, or as a collection of many compact HII regions in non-LTE. The multiple HII region model provides a better fit to the data and predicts many dense (ne ~ 10^4-10^5 cm-3), compact (< 0.1 pc) HII regions. For the whole nuclear region as defined by RRL emission, we estimate a hydrogen ionizing rate of NL ~ 2 x 10^{52} s^{-1}, corresponding to equivalent ~ 2000 O6 stars and a star formation rate of ~ 0.15 Msun/year. We detect radio continuum emission west of the southern molecular mini spiral arm, consistent with trailing spiral arms.