Using the star cluster catalogs from the Hubble Space Telescope program Legacy ExtraGalactic UV survey (LEGUS) and 8 $\mu$m images from the IRAC camera on the Spitzer Space Telescope for 5 galaxies within 5 Mpc, we investigate how the 8 $\mu$m dust luminosity correlates with the stellar age on the 30--50 pc scale of star forming regions. We construct a sample of 97 regions centered at local peaks of 8 $\mu$m emission, each containing one or more young star cluster candidates from the LEGUS catalogs. We find a tight anti-correlation with a Pearson correlation coefficient of $r=-0.84\pm0.05$ between the mass-normalized dust-only 8 $\mu$m luminosity and the age of stellar clusters younger than 1 Gyr; the 8 $\mu$m luminosity decreases with increasing age of the stellar population. Simple assumptions on a combination of stellar and dust emission models reproduce the observed trend. We also explore how the scatter of the observed trend depends on assumptions of stellar metallicity, PAH abundance, fraction of stellar light absorbed by dust, and instantaneous versus continuous star formation models. We find that variations in stellar metallicity have little effect on the scatter, while PAH abundance and the fraction of dust-absorbed light bracket the full range of the data. We also find that the trend is better explained by continuous star formation, rather than instantaneous burst models. We ascribe this result to the presence of multiple star clusters with different ages in many of the regions. Upper limits of the dust-only 8 $\mu$m emission as a function of age are provided.

Star clusters are fundamental units of stellar feedback and unique tracers of their host galactic properties. In this review, we will first focus on their constituents, i.e.\ detailed insight into their stellar populations and their surrounding ionised, warm, neutral, and molecular gas. We, then, move beyond the Local Group to review star cluster populations at various evolutionary stages, and in diverse galactic environmental conditions accessible in the local Universe. At high redshift, where conditions for cluster formation and evolution are more extreme, we are only able to observe the integrated light of a handful of objects that we believe will become globular clusters. We therefore discuss how numerical and analytical methods, informed by the observed properties of cluster populations in the local Universe, are used to develop sophisticated simulations potentially capable of disentangling the genetic map of galaxy formation and assembly that is carried by globular cluster populations.

Molecular clouds are a fundamental ingredient of galaxies: they are the channels that transform the diffuse gas into stars. The detailed process of how they do it is not completely understood. We review the current knowledge of molecular clouds and their substructure from scales $\sim~$1~kpc down to the filament and core scale. We first review the mechanisms of cloud formation from the warm diffuse interstellar medium down to the cold and dense molecular clouds, the process of molecule formation and the role of the thermal and gravitational instabilities. We also discuss the main physical mechanisms through which clouds gather their mass, and note that all of them may have a role at various stages of the process. In order to understand the dynamics of clouds we then give a critical review of the widely used virial theorem, and its relation to the measurable properties of molecular clouds. Since these properties are the tools we have for understanding the dynamical state of clouds, we critically analyse them. We finally discuss the ubiquitous filamentary structure of molecular clouds and its connection to prestellar cores and star formation.

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.

Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

Using MUSE AO data, we probe the ISM in the local spiral galaxy NGC 7793 at a spatial resolution of $\sim$ 10 pc. We identify HII regions and compile a catalogue of supernova remnants (SNRs), planetary nebulae (PNe) and Wolf Rayet stars (WR). We compute electron densities and temperatures from the [SII]6716/6731 and [SIII]6312/9069 line ratios. We study the properties of the ionised gas through BPT diagrams combined with gas velocity dispersion. We spectroscopically confirm 2 previously detected WR stars and a SNR and report the discovery of 7 WR, 1 SNR, and 2 PNe. The diffuse ionized gas (DIG) fraction is between $\sim$ 27 and 42% depending on the method used to define the HII region boundaries. In agreement with previous studies, we find that the DIG exhibits enhanced [SII]/H$\alpha$ and [NII]/H$\alpha$ ratios and a median temperature that is $\sim$ 3000 K higher than in HII regions. We also observe an apparent inverse correlation between temperature and H$\alpha$ surface brightness. Overall, the observed [SII]6716/6731 ratio is consistent within 1$\sigma$ with $n_e$ < 30 cm$^{-3}$, with an almost identical distribution for the DIG and HII regions. The velocity dispersion of the ionised gas indicates that the DIG has a higher degree of turbulence than the HII regions. Comparison with photoionisation and shock models reveals that the diffuse component can only partially be explained via shocks and that it is most likely consistent with photons leaking from density bounded HII regions or with radiation from evolved field stars. [abridged]

Only about 19 Galactic and 25 extra-galactic bona-fide Luminous Blue Variables (LBVs) are known to date. This incomplete census prevents our understanding of this crucial phase of massive star evolution which leads to the formation of heavy binary black holes via the classical channel. With large samples of LBVs one could better determine the duration and maximum stellar luminosity which characterize this phase. We search for candidate LBVs (cLBVs) in a new galaxy, NGC 7793. For this purpose, we combine high spatial resolution images from two Hubble Space Telescope (HST) programs with optical spectroscopy from the Multi Unit Spectroscopic Explorer (MUSE). By combining PSF-fitting photometry measured on F547M, F657N, and F814W images, with restrictions on point-like appearance (at HST resolution) and H{\alpha} luminosity, we find 100 potential cLBVs, 36 of which fall in the MUSE fields. Five of the latter 36 sources are promising cLBVs which have M(V) less than or equal to -7 and a combination of: H{\alpha} with a P-Cygni profile; no [O I] 6300 emission; weak or no [O III] 5007 emission; large [N II]/H{\alpha} relative to H II regions; and [S II] 6716 / [S II] 6731 similar to 1. It is not clear if these five cLBVs are isolated from O-type stars, which would favor the binary formation scenario of LBVs. Our study, which approximately covers one fourth of the optical disc of NGC 7793, demonstrates how by combining the above HST surveys with multi-object spectroscopy from 8-m class telescopes, one can efficiently find large samples of cLBVs in nearby galaxies.

The relative average minimum projected separations of star clusters in the Legacy ExtraGalactic UV Survey (LEGUS) and in tidal dwarfs around the interacting galaxy NGC 5291 are determined as a function of cluster mass to look for cluster-cluster mass segregation. Class 2 and 3 LEGUS clusters, which have a more irregular internal structure than the compact and symmetric class 1 clusters, are found to be mass segregated in low mass galaxies, which means that the more massive clusters are systematically bunched together compared to the lower mass clusters. This mass segregation is not present in high-mass galaxies nor for class 1 clusters. We consider possible causes for this segregation including differences in cluster formation and scattering in the shallow gravitational potentials of low mass galaxies.

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.

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.

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.

SIGNALS, the Star formation, Ionized Gas, and Nebular Abundances Legacy Survey, is a large observing program designed to investigate massive star formation and HII regions in a sample of local extended galaxies. The program will use the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope. Over 355 hours (54.7 nights) have been allocated beginning in fall 2018 for eight consecutive semesters. Once completed, SIGNALS will provide a statistically reliable laboratory to investigate massive star formation, including over 50 000 resolved HII regions : the largest, most complete, and homogeneous database of spectroscopically and spatially resolved extragalactic HII regions ever assembled. For each field observed, three datacubes covering the spectral bands of the filters SN1 (363 -386 nm), SN2 (482 - 513 nm), and SN3 (647 - 685 nm) are gathered. The spectral resolution selected for each spectral band is 1000, 1000, and 5000, respectively. As defined, the project sample will facilitate the study of small-scale nebular physics and many other phenomena linked to star formation at a mean spatial resolution of 20 pc. This survey also has considerable legacy value for additional topics including planetary nebulae, diffuse ionized gas, andsupernova remnants. The purpose of this paper is to present a general outlook of the survey, notably the observing strategy, galaxy sample, and science requirements.

We have used the capability of the MUSE instrument to explore the impact of stellar feedback at large scales in Haro 11, a galaxy under extreme starburst condition and one of the first galaxies where Lyman continuum (LyC) has been detected. Using Ha, [OIII] and [OI] emission lines from deep MUSE observations, we have constructed a sequence of velocity-dependent maps of the Ha emission, the state of the ionised gas and a tracer of fast shocks. These allowed us to investigate the ionisation structure of the galaxy in 50 kms^2 bins over a velocity range of -400 to 350 kms. The ionised gas in Haro 11 is assembled by a rich arrangement of structures, such as superbubbles, filaments, arcs and galactic ionised channels, whose appearances change drastically with velocity. The central star forming knots and the star forming dusty arm are the main engines that power the strong mechanical feedback in this galaxy, although with different impact on the ionisation structure. Haro 11 appears to leak LyC radiation in many directions. We found evidence of a kpc-scale fragmented superbubble, that may have cleared galactic-scale channels in the ISM. Additionally, the southwestern hemisphere is highly ionised in all velocities, hinting at a density bound scenario. A compact kpc-scale structure of lowly ionised gas coincides with the diffuse Lya emission and the presence of fast shocks. Finally, we find evidence that a significant fraction of the ionised gas mass may escape the gravitational potential of the galaxy.

We study young star-forming clumps on physical scales of 10-500 pc in the Lyman-Alpha Reference Sample (LARS), a collection of low-redshift (z = 0.03-0.2) UV-selected star-forming galaxies. In each of the 14 galaxies of the sample, we detect clumps for which we derive sizes and magnitudes in 5 UV-optical filters. The final sample includes $\sim$1400 clumps, of which $\sim$600 have magnitude uncertainties below 0.3 in all filters. The UV luminosity function for the total sample of clumps is described by a power-law with slope $\alpha = -2.03^{+0.11}_{-0.13}$. Clumps in the LARS galaxies have on average $\Sigma_{SFR}$ values higher than what observed in HII regions of local galaxies and comparable to typical SFR densities of clumps in z = 1-3 galaxies. We derive the clumpiness as the relative contribution from clumps to the UV emission of each galaxy, and study it as a function of galactic-scale properties, i.e. $\Sigma_{SFR}$ and the ratio between rotational and dispersion velocities of the gas ($v_s/\sigma_0$). We find that in galaxies with higher $\Sigma_{SFR}$ or lower $v_s/\sigma_0$, clumps dominate the UV emission of their host systems. All LARS galaxies with Ly$\alpha$ escape fractions larger than 10% have more than 50% of the UV luminosity from clumps. We tested the robustness of these results against the effect of different physical resolutions. At low resolution, the measured clumpiness appears more elevated than if we could resolve clumps down to single clusters. This effect is small in the redshift range covered by LARS, thus our results are not driven by the physical resolution.

X-ray radiation from accreting compact objects is an important part of stellar feedback. The metal-poor galaxy ESO 338-4 has experienced vigorous starburst during the last 40 Myr and contains some of the most massive super star clusters in the nearby Universe. Given its starburst age and its star-formation rate, ESO 338-4 is one of the most efficient nearby manufactures of neutron stars and black holes, hence providing an excellent laboratory for feedback studies. We compared X-ray images and spectra obtained by XMM-Newton and Chandra telescopes with integral field spectroscopic VLT MUSE observations in the optical to constrain the nature of strong X-ray emitters. X-ray observations uncover three ultraluminous X-ray sources (ULXs) in ESO 338-4. The brightest among them, ESO 338~X-1, has X-ray luminosity in excess of 10^{40} erg/s. We speculate that ESO 338-4 is powered by accretion on an intermediate-mass (~300Msun) black hole. We show that X-ray radiation from ULXs and hot superbubbles strongly contributes to HeII ionization and general stellar feedback in this template starburst galaxy.

We use two unsupervised machine learning algorithms, Shared Nearest Neighbor (SNN) and \textit{EnLink}, as a single approach to identify stellar groupings in the Orion star-forming complex as an application to the 5-dimensional astrometric data from Gaia DR2. The algorithms present two different ways to limit user bias when evaluating the relative weights among the astrometric parameters, automatically determined by the machine and through a standard procedure by monitoring several outcome measures. Both algorithms complement each other and produce similar stellar groups. Because SNN groups have a much smaller spread in proper motions compared to \textit{EnLink}, we use \textit{EnLink}, which requires no input, as a first pass tool for group identification and validation. We then used the SNN algorithm to dissect the Orion star-forming complex. We identify 21 spatially- and kinematically-coherent groups in the Orion complex, nine of which previously unknown. The groups show a wide distribution of distances extending as far as about 150 pc in front of the star-forming Orion clouds, to about 50 pc beyond them where we find, unexpectedly, three groups. Our results expose to view the wealth of sub-structure in the OB association, within and beyond the classical Blaauw Orion OBI sub-groups. A full characterization of the new groups is of the essence as it offers the potential to unveil how star formation proceeds globally in large complexes such as Orion. The data and code that generated the groups is provided in this Letter.

We present a detailed study of the flocculent spiral galaxy NGC 7793, part of the Sculptor group. By analyzing the resolved stellar populations of the galaxy, located at a distance of ~3.7 Mpc, we infer for the first time its radial star formation history (SFH) from Hubble Space Telescope photometry, thanks to both archival and new data from the Legacy ExtraGalactic UV Survey. We determine an average star formation rate (SFR) for the galaxy portion covered by our F555W and F814W data of 0.23 +- 0.02 Msun/yr over the whole Hubble time, corresponding to a total stellar mass of 3.09 +- 0.33 x 10^9 Msun in agreement with previous determinations. Thanks to the new data extending to the F336W band, we are able to analyze the youngest stellar populations with a higher time resolution. Most importantly, we recover the resolved SFH in different radial regions of the galaxy; this shows an indication of a growing trend of the present-to-past SFR ratio, increasing from internal to more external regions, supporting previous findings of the inside-out growth of the galaxy.

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 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.

We describe a novel method for determining the demographics of a population of star clusters, for example distributions of cluster mass and age, from unresolved photometry. This method has a number of desirable properties: it fully exploits all the information available in a data set without any binning, correctly accounts for both measurement error and sample incompleteness, naturally handles heterogenous data (for example fields that have been imaged with different sets of filters or to different depths), marginalises over uncertain extinctions, and returns the full posterior distributions of the parameters describing star cluster demographics. We demonstrate the method using mock star cluster catalogs and show that our method is robust and accurate, and that it can recover the demographics of star cluster populations significantly better than traditional fitting methods. For realistic sample sizes, our method is sufficiently powerful that its accuracy is ultimately limited by the accuracy of the underlying physical models for stellar evolution and interstellar dust, rather than by statistical uncertainties. Our method is implemented as part of the Stochastically Lighting Up Galaxies (SLUG) stellar populations code, and is freely available.

Stellar feedback strongly affects the interstellar medium (ISM) of galaxies. Stellar feedback in the first galaxies likely plays a major role in enabling the escape of LyC photons, which contribute to the re-ionization of the Universe. Nearby starburst galaxies serve as local analogues allowing for a spatially resolved assessment of the feedback processes in these galaxies. We characterize the feedback effects from the star clusters in the local high-redshift analogue ESO 338-IG04 on the ISM and compare the results with the properties of the most massive clusters. We use high quality VLT/MUSE optical integral field data to derive the physical properties of the ISM such as ionization, density, shocks, and perform new fitting of the spectral energy distributions of the brightest clusters in ESO 338-IG04 from HST imaging. ESO 338-IG04 has a large ionized halo which we detect to a distance of 9 kpc. We identify 4 Wolf-Rayet (WR) clusters based on the blue and red WR bump. We follow previously identified ionization cones and find that the ionization of the halo increases with distance. Analysis of the galaxy kinematics shows two complex outflows driven by the numerous young clusters in the galaxy. We find a ring of shocked emission traced by an enhanced [OI]/H$\alpha$ ratio surrounding the starburst and at the end of the outflow. Finally we detect nitrogen enriched gas associated with the outflow, likely caused by the WR stars in the massive star clusters. Photo-ionization dominates the central starburst and sets the ionization structure of the entire halo, resulting in a density bounded halo, facilitating the escape of LyC photons. Outside the central starburst, shocks triggered by an expanding super bubble become important. The shocks at the end of the outflow suggest interaction between the hot outflowing material and the more quiescent halo gas.

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 catalog from the HST LEGUS (Legacy ExtraGalactic UV Survey) program with the 15 parsec 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 timescale 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 \pm 0.03$, significantly more clustered than the molecular cloud hierarchy with slope of $-0.18 \pm 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\pm0.03$ and $-0.45\pm0.06$ respectively, when we compare massive ($>$10$^5$ M$_{\odot}$) 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 percent.

We explore the possibility of using UV spectroscopy in combination with broad-band photometry as diagnostic tools for understanding the shape of the Initial Mass Function (IMF) in unresolved stellar populations. Building on our previous work, we extend the Stochastically Lighting Up Galaxies code (slug) to include a high-resolution UV spectral synthesiser and equivalent width calculation capabilities. We first gain a qualitative understanding of how UV spectral features behave as the parameters that define a star cluster in slug (mass, age, extinction, and IMF slope alpha) are changed. We then exploit Bayesian inference techniques to recover the alpha values for clusters simulated with slug, using mock observations of these clusters comprised of broad-band photometry and equivalent width measurements of a selection of UV spectral features. We find some improvement when compared to attempts using broad-band photometry alone (with the interquartile range of the alpha posterior PDF shrinking by 32%), although we still do not yet fully break the known degeneracy between the cluster mass and alpha. Finally, we make predictions about how effective real observations will be by quantifying our ability to constrain alpha as a function of limiting equivalent width. We find that observations sensitive to a modest equivalent width of 9 A are sufficient to improve the recovery of the IMF slope parameter by 32% (interquartile range of posterior PDF median residuals), moving to 39% when we include all the significant spectral features in the wavelength range 900-3000 A.

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.

Star-forming clumps dominate the rest-frame ultraviolet morphology of galaxies at the peak of cosmic star formation. If turbulence driven fragmentation is the mechanism responsible for their formation, we expect their stellar mass function to follow a power-law of slope close to $-2$. We test this hypothesis performing the first analysis of the stellar mass function of clumps hosted in galaxies at $z\sim 1-3.5$. The clump sample is gathered from the literature with similar detection thresholds and stellar masses determined in a homogeneous way. To overcome the small number statistics per galaxy (each galaxy hosts up to a few tens of clumps only), we combine all high-redshift clumps. The resulting clump mass function follows a power-law of slope $\sim -1.7$ and flattens at masses below $2\times 10^7$ M$_{\odot}$. By means of randomly sampled clump populations, drawn out of a power-law mass function of slope $-2$, we test the effect of combining small clump populations, detection limits of the surveys, and blending on the mass function. Our numerical exercise reproduces all the features observed in the real clump mass function confirming that it is consistent with a power-law of slope $\simeq -2$. This result supports the high-redshift clump formation through fragmentation in a similar fashion as in local galaxies, but under different gas conditions.

One of the main theories for explaining the formation of spiral arms in galaxies is the stationary density wave theory. This theory predicts the existence of an age gradient across the arms. We use the stellar cluster catalogues of the galaxies NGC 1566, M51a, and NGC 628 from the Legacy Extragalactic UV Survey (LEGUS) program. In order to test for the possible existence of an age sequence across the spiral arms, we quantified the azimuthal offset between star clusters of different ages in our target galaxies. We found that NGC 1566, a grand-design spiral galaxy with bisymmetric arms and a strong bar, shows a significant age gradient across the spiral arms that appears to be consistent with the prediction of the stationary density wave theory. In contrast, M51a with its two well-defined spiral arms and a weaker bar does not show an age gradient across the arms. In addition, a comparison with non LEGUS star cluster catalogues for M51a yields similar results. We believe that the spiral structure of M51a is not the result of a stationary density wave with a fixed pattern speed. Instead, tidal interactions could be the dominant mechanism for the formation of spiral arms. We also found no offset in the azimuthal distribution of star clusters with different ages across the weak spiral arms of NGC 628.

We have identified a progenitor candidate in archival Hubble Space Telescope (HST) images for the Type Ic SN 2017ein in NGC 3938, pinpointing the candidate's location via HST Target-of-Opportunity imaging of the SN itself. This would be the first identification of a stellar-like object as a progenitor candidate for any Type Ic supernova to date. We also present observations of SN 2017ein during the first ~49 days since explosion. We find that SN 2017ein most resembles the well-studied Type Ic SN 2007gr. We infer that SN 2017ein experienced a total visual extinction of A_V~1.0--1.9 mag, predominantly because of dust within the host galaxy. Although the distance is not well known, if this object is the progenitor, it was likely of high initial mass, ~47--48 M_sun if a single star, or ~60--80 M_sun if in a binary system. However, we also find that the progenitor candidate could be a very blue and young compact cluster, further implying a very massive (>65 M_sun) progenitor. Furthermore, the actual progenitor might not be associated with the candidate at all and could be far less massive. From the immediate stellar environment, we find possible evidence for three different populations; if the SN progenitor was a member of the youngest population, this would be consistent with an initial mass of ~57 M_sun. After it has faded, the SN should be reobserved at high spatial resolution and sensitivity, to determine whether the candidate is indeed the progenitor.

Minor mergers are important processes contributing significantly to how galaxies evolve across the age of the Universe. Their impact on supermassive black hole growth and star formation is profound. The detailed study of dense molecular gas in galaxies provides an important test of the validity of the relation between star formation rate and HCN luminosity on different galactic scales. We use observations of HCN, HCO+1-0 and CO3-2 to study the dense gas properties in the Medusa merger. We calculate the brightness temperature ratios and use them in conjunction with a non-LTE radiative line transfer model. The HCN and HCO+1-0, and CO3-2 emission do not occupy the same structures as the less dense gas associated with the lower-J CO emission. The only emission from dense gas is detected in a 200pc region within the "Eye of the Medusa". No HCN or HCO+ is detected for the extended starburst. The CO3-2/2-1 brightness temperature ratio inside "the Eye" is ~2.5 - the highest ratio found so far. The line ratios reveal an extreme, fragmented molecular cloud population inside "the Eye" with large temperatures (>300K) and high gas densities (>10^4 cm^-3). "The Eye" is found at an interface between a large-scale minor axis inflow and the Medusa central region. The extreme conditions inside "the Eye" may be the result of the radiative and mechanical feedback from a deeply embedded, young, massive super star cluster, formed due to the gas pile-up at the intersection. Alternatively, shocks from the inflowing gas may be strong enough to shock and fragment the gas. For both scenarios, however, it appears that the HCN and HCO+ dense gas tracers are not probing star formation, but instead a post-starburst and/or shocked ISM that is too hot and fragmented to form new stars. Thus, caution is advised in linking the detection of emission from dense gas tracers to evidence of ongoing or imminent star formation.

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.

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.

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 broad band WFC3/UVIS UV and optical images from the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with archival HST/ACS 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 HI and CO gas maps from literature, including The HI Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey (HERACLES). 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 H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$ 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 (BEAST) for NGC 7793 and find systematically lower extinctions from SED-fitting as compared to isochrone matching.

We use HST observations from the Legacy Extragalactic UV Survey to reconstruct the recent star formation histories (SFHs) of three actively star-forming dwarf galaxies, NGC4449, Holmberg II and NGC1705, 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 modelling. Irrespective of the adopted stellar models, all the 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 $\sim$2) above the 100 Myr-averaged-SFR. Significant differences among the three dwarfs are found in the overall SFR, the timing of the most recent peak and the SFR$/$area. The Initial Mass Function (IMF) of NGC1705 and Holmberg II is consistent with a Salpeter slope down to $\approx$ 5 M$_{\odot}$, whereas it is slightly flatter, s$=-2.0$, in NGC4449. 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 upper main sequence and helium burning stars, which is not apparent in the data. Since differential reddening, significant in NGC4449, or unresolved binaries don't appear to be sufficient to fill the gap, we suggest this calls for a revision of both sets of stellar evolutionary tracks.

Two dwarf irregular galaxies DDO 187 and NGC 3738 exhibit a striking pattern of star formation: intense star formation is taking place in a large region occupying roughly half of the inner part of the optical galaxy. We use data on the HI distribution and kinematics and stellar images and colors to examine the properties of the environment in the high star formation rate (HSF) halves of the galaxies in comparison with the low star formation rate (LSF) halves. We find that the pressure and gas density are higher on the HSF sides by 30-70%. In addition we find in both galaxies that the HI velocity fields exhibit significant deviations from ordered rotation and there are large regions of high velocity dispersion and multiple velocity components in the gas beyond the inner regions of the galaxies. The conditions in the HSF regions are likely the result of large-scale external processes affecting the internal environment of the galaxies and enabling the current star formation there.

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 five 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 Mo, 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 time scale 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.

Recently acquired WFC3 UV (F275W and F336W) imaging mosaics under the Legacy Extragalactic UV Survey (LEGUS) combined with archival ACS data of M51 are used to study the young star cluster (YSC) population of this interacting system. Our newly extracted source catalogue contains 2834 cluster candidates, morphologically classified to be compact and uniform in colour, for which ages, masses and extinction are derived. In this first work we study the main properties of the YSC population of the whole galaxy, considering a mass-limited sample. Both luminosity and mass functions follow a power law shape with slope -2, but at high luminosities and masses a dearth of sources is observed. The analysis of the mass function suggests that it is best fitted by a Schechter function with slope -2 and a truncation mass at $1.00\pm0.12\times10^5\ M_\odot$. Through Monte Carlo simulations we confirm this result and link the shape of the luminosity function to the presence of a truncation in the mass function. A mass limited age function analysis, between 10 and 200 Myr, suggests that the cluster population is undergoing only moderate disruption. We observe little variation in the shape of the mass function at masses above $1\times10^4\ M_\odot$, over this age range. The fraction of star formation happening in the form of bound clusters in M51 is $\sim20\%$ in the age range 10 to 100 Myr and little variation is observed over the whole range from 1 to 200 Myr.

We present an analysis of the positions and ages of young star clusters in eight local galaxies to investigate the connection between the age difference and separation of cluster pairs. We find that star clusters do not form uniformly but instead are distributed such that the age difference increases with the cluster pair separation to the 0.25-0.6 power, and that the maximum size over which star formation is physically correlated ranges from ~200 pc to ~1 kpc. The observed trends between age difference and separation suggest that cluster formation is hierarchical both in space and time: clusters that are close to each other are more similar in age than clusters born further apart. The temporal correlations between stellar aggregates have slopes that are consistent with turbulence acting as the primary driver of star formation. The velocity associated with the maximum size is proportional to the galaxy's shear, suggesting that the galactic environment influences the maximum size of the star-forming structures.

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 ($\alpha_3$) 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 $\alpha_3$ 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 $\alpha_3$ when adding H$\alpha$ photometry from the companion H$\alpha$-LEGUS survey. Conversely, H$\alpha$ 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 present a study of the effective (half-light) radii and other structural properties of a systematically selected sample of young, massive star clusters (YMCs, $\geq$$5\times10^3$ M$_{\odot}$ and $\leq$200 Myr) in two nearby spiral galaxies, NGC 628 and NGC 1313. We use Hubble Space Telescope WFC3/UVIS and archival ACS/WFC data obtained by the Legacy Extragalactic UV Survey (LEGUS), an HST Treasury Program. We measure effective radii with GALFIT, a two-dimensional image-fitting package, and with a new technique to estimate effective radii from the concentration index (CI) of observed clusters. The distribution of effective radii from both techniques spans $\sim$0.5-10 pc and peaks at 2-3 pc for both galaxies. We find slight positive correlations between effective radius and cluster age in both galaxies, but no significant relationship between effective radius and galactocentric distance. Clusters in NGC 1313 display a mild increase in effective radius with cluster mass, but the trend disappears when the sample is divided into age bins. We show that the vast majority of the clusters in both galaxies are much older than their dynamical times, suggesting they are gravitationally bound objects. We find that about half of the clusters in NGC 628 are underfilling their Roche lobes, based on their Jacobi radii. Our results suggest that the young, massive clusters in NGC 628 and NGC 1313 are expanding due to stellar mass loss or two-body relaxation and are not significantly influenced by the tidal fields of their host galaxies.

We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues 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 NGC628 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 NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628.

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 broad band WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey). We have 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.

Scientific illustrations, thanks to the vision of great artists fascinated by astronomical research and astronautics, have provided us with an accurate depiction of the possible views which mankind will one day observed from locations other than our planet. In this talk I will pay homage to some of these geniuses who serve science, and underline the scientific, artistic, political, and social implications deriving from a wise use of space-art.

The work of Adamo et al. (2015) showed that the mass distributions of young massive stellar clusters were truncated above a maximum-mass scale in the nearby galaxy M83 and that this truncation mass varies with galactocentric radius. Here, we present a cloud-based analysis of ALMA CO($1\to 0$) observations of M83 to search for such a truncation mass in the molecular cloud population. We identify a population of 873 molecular clouds in M83 that is largely similar to those found in the Milky Way and Local Group galaxies, though clouds in the centre of the galaxy show high surface densities and enhanced turbulence, as is common for clouds in high-density nuclear environments. Like the young massive clusters, we find a maximum-mass scale for the molecular clouds that decreases radially in the galaxy. We find the most massive young massive cluster tracks the most massive molecular cloud with the cluster mass being $10^{-2}$ times that of the most massive molecular cloud. Outside the nuclear region of M83 ($R_{g}>0.5$ kpc), there is no evidence for changing internal conditions in the population of molecular clouds, with the average internal pressures, densities, and free-fall times remaining constant for the cloud population over the galaxy. This result is consistent with the bound cluster formation efficiency depending only on the large-scale properties of the ISM, rather than the internal conditions of individual clouds.

We investigate how star formation is spatially organized in the grand-design spiral NGC 1566 from deep HST photometry with the Legacy ExtraGalactic UV Survey (LEGUS). Our contour-based clustering analysis reveals 890 distinct stellar conglomerations at various levels of significance. These star-forming complexes are organized in a hierarchical fashion with the larger congregations consisting of smaller structures, which themselves fragment into even smaller and more compact stellar groupings. Their size distribution, covering a wide range in length-scales, shows a power-law as expected from scale-free processes. We explain this shape with a simple "fragmentation and enrichment" model. The hierarchical morphology of the complexes is confirmed by their mass--size relation which can be represented by a power-law with a fractional exponent, analogous to that determined for fractal molecular clouds. The surface stellar density distribution of the complexes shows a log-normal shape similar to that for supersonic non-gravitating turbulent gas. Between 50 and 65 per cent of the recently-formed stars, as well as about 90 per cent of the young star clusters, are found inside the stellar complexes, located along the spiral arms. We find an age-difference between young stars inside the complexes and those in their direct vicinity in the arms of at least 10 Myr. This timescale may relate to the minimum time for stellar evaporation, although we cannot exclude the in situ formation of stars. As expected, star formation preferentially occurs in spiral arms. Our findings reveal turbulent-driven hierarchical star formation along the arms of a grand-design galaxy.

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 Myrs) 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 onto 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.

Interarm star formation contributes significantly to a galaxy's star formation budget, and provides an opportunity to study stellar birthplaces unperturbed by spiral arm dynamics. Using optical integral field spectroscopy of the nearby galaxy NGC 628 with VLT/MUSE, we construct Halpha maps including detailed corrections for dust extinction and stellar absorption to identify 391 HII regions at 35pc resolution over 12 kpc^2. Using tracers sensitive to the underlying gravitational potential, we associate HII regions with either arm (271) or interarm (120) environments. Using our full spectral coverage of each region, we find that most HII region physical properties (luminosity, size, metallicity, ionization parameter) are independent of environment. We calculate the fraction of Halpha luminosity due to the diffuse ionized gas (DIG) background contaminating each HII region, and find the DIG surface brightness to be higher within HII regions compared to the surroundings, and slightly higher within arm HII regions. Use of the temperature sensitive [SII]/Halpha line ratio map instead of the Halpha surface brightness to identify HII region boundaries does not change this result. Using the dust attenuation as a tracer of the gas, we find depletion times consistent with previous work (2 x 10^9 yr) with no differences between the arm and interarm, however this is very sensitive to the DIG correction. Unlike molecular clouds, which can be dynamically affected by the galactic environment, we see fairly consistent HII region properties in both arm and interarm environments. This suggests either a difference in arm star formation and feedback, or a decoupling of dense star forming clumps from the more extended surrounding molecular gas.

We present results of a photometric study into the cluster population of NGC 1566, a nearby grand design spiral galaxy, sampled out to a Galactocentric radius of $\approx 5.5$ kpc. The shape of the mass-limited age distribution shows negligible variation with radial distance from the centre of the galaxy, and demonstrates three separate sections, with a steep beginning, flat middle and steep end. The luminosity function can be approximated by a power law at lower luminosities with evidence of a truncation at higher luminosity. The power law section of the luminosity function of the galaxy is best fitted by an index $\approx -2$, in agreement with other studies, and is found to agree with a model luminosity function, which uses an underlying Schechter mass function. The recovered power law slope of the mass distribution shows a slight steepening as a function of galactocentric distance, but this is within error estimates. It also displays a possible truncation at the high mass end. Additionally, the cluster formation efficiency ($\Gamma$) and the specific U-band luminosity of clusters ($T_L(U)$) are calculated for NGC 1566 and are consistent with values for similar galaxies. A difference in NGC 1566, however, is that the fairly high star formation rate is in contrast with a low $\Sigma_{SFR}$ and $\Gamma$, indicating that $\Gamma$ can only be said to depend strongly on $\Sigma_{SFR}$, not the star formation rate.

We test the predictions of spectral synthesis models based on seven different massive-star prescriptions against Legacy ExtraGalactic UV Survey (LEGUS) observations of eight young massive clusters in two local galaxies, NGC 1566 and NGC 5253, chosen because predictions of all seven models are available at the published galactic metallicities. The high angular resolution, extensive cluster inventory and full near-ultraviolet to near-infrared photometric coverage make the LEGUS dataset excellent for this study. We account for both stellar and nebular emission in the models and try two different prescriptions for attenuation by dust. From Bayesian fits of model libraries to the observations, we find remarkably low dispersion in the median E(B-V) (~0.03 mag), stellar masses (~10^4 M_\odot) and ages (~1 Myr) derived for individual clusters using different models, although maximum discrepancies in these quantities can reach 0.09 mag and factors of 2.8 and 2.5, respectively. This is for ranges in median properties of 0.05-0.54 mag, 1.8-10x10^4 M_\odot and 1.6-40 Myr spanned by the clusters in our sample. In terms of best fit, the observations are slightly better reproduced by models with interacting binaries and least well reproduced by models with single rotating stars. Our study provides a first quantitative estimate of the accuracies and uncertainties of the most recent spectral synthesis models of young stellar populations, demonstrates the good progress of models in fitting high-quality observations, and highlights the needs for a larger cluster sample and more extensive tests of the model parameter space.

While numerical simulations suggest that the strength of the Lyman alpha (Lya) line of star-forming disk galaxies strongly depends on the inclination at which they are observed (i.e. from edge-on to face-on, we expect to see a change from an attenuated Lya line to a strong Lya emission line), recent observations with the Hubble space telescope (HST) have highlighted few low-redshift highly inclined (edge-on) disk galaxies that breaks this trend. We aim to understand how a strong Lya emission line is able to escape from one of those inclined disk galaxies, named Mrk1486 (z=0.0338). For that purpose we used a large set of HST imaging and spectroscopic data to investigate both the ISM structure and the dominant source of Lya radiation inside Mrk1486. Moreover, we used a 3D Monte Carlo Lya radiation transfer code to study the radiative transfer of Lya and UV continuum photons inside a 3D geometry of neutral hydrogen (HI) and dust that models the ISM structure at the galaxy center. The analysis of IFU Halpha spectroscopic data of Mrk1486 indicates the presence of two bipolar galactic winds of HI gas above and bellow the disk plane of Mrk1486. Furthermore, comparing different diagnostic diagrams (such as [OIII]5007/Hbeta versus [OI]6300/Halpha) to photo- and shock-ionization models, we find that the Lya production of Mrk1486 is dominated by photoionization inside the galaxy disk. From this perspective, our numerical simulations succeed in reproducing the strength and spectral shape of the observed Lya line of Mrk1486 by assuming a scenario in which the Lya photons are produced inside the disk, travel along the galactic winds and scatter on cool HI materials toward the observer. As bipolar galactic winds are ubiquitous in star-forming disk galaxies, this mechanism may explain the origin of strong Lya emission lines recently observed from highly inclined galaxies at high-redshift.

We review many of the basic properties of star cluster systems, and focus in particular on how they relate to their host galaxy properties and ambient environment. The cluster mass and luminosity functions are well approximated by power-laws of the form $Ndm \propto M^{\alpha}dm$, with $\alpha\sim-2$ over most of the observable range. However, there is now clear evidence that both become steeper at high masses/luminosities, with the value of the downward turn dependent on environment. The host galaxy properties also appear to affect the cluster formation efficiency ($\Gamma$ - i.e., the fraction of stars that form in bound clusters), with higher star-formation rate density galaxies having higher $\Gamma$ values. Within individual galaxies, there is evidence for $\Gamma$ to vary by a factor of 3-4, likely following the molecular gas surface density, in agreement with recent predictions. Finally, we discuss cluster disruption and its effect on the observed properties of a population, focussing on the age distribution of clusters. We briefly discuss the expectations of theoretical and numerical studies, and also the observed distributions in a number of galaxies. Most observational studies now find agreement with theoretical expectations, namely nearly a constant cluster age distribution for ages up to ~100 Myr (i.e. little disruption), and a drastic steepening above this value caused by a combination of cluster disruption and incompleteness. Rapid cluster disruption for clusters with ages < 100 Myr is ruled out for most galaxies.

Young star clusters (YSCs) appear to be a ubiquitous product of star formation in local galaxies, thus, they can be used to study the star formation process at work in their host galaxies. Moreover, YSCs are intrinsically brighter that single stars, potentially becoming the most important tracers of the recent star formation history in galaxies in the local Universe. In local galaxies, we also witness the presence of a large population of evolved star clusters, commonly called globular clusters (GCs). GCs peak formation history is very close to the redshift (z~2) when the cosmic star formation history reached the maximum. Therefore, GCs are usually associated to extreme star formation episodes in high-redshift galaxies. It is yet not clear whether YSCs and GCs share a similar formation process (same physics under different interstellar medium conditions) and evolution process, and whether the former can be used as progenitor analogs of the latter. In this invited contribution, I review general properties of YSC populations in local galaxies. I will summarise some of the current open questions in the field, with particular emphasis to whether or not galactic environments, where YSCs form, leave imprints on the nested populations. The importance of this rapidly developing field can be crucial in understanding GC formation and possibly the galactic environment condition where this ancient population formed.

We present integral field spectroscopic observations with the Potsdam Multi Aperture Spectrophotometer of all 14 galaxies in the $z\sim 0.1$ Lyman Alpha Reference Sample (LARS). We produce 2D line of sight velocity maps and velocity dispersion maps from the Balmer $\alpha$ (H$\alpha$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$\alpha$ radiation field. We show our kinematic maps spatially registered onto the Hubble Space Telescope H$\alpha$ and Lyman $\alpha$ (Ly$\alpha$) images. Only for individual galaxies a causal connection between spatially resolved H$\alpha$ kinematics and Ly$\alpha$ photometry can be conjectured. However, no general trend can be established for the whole sample. Furthermore, we compute non-parametric global kinematical statistics -- intrinsic velocity dispersion $\sigma_0$, shearing velocity $v_\mathrm{shear}$, and the $v_\mathrm{shear}/\sigma_0$ ratio -- from our kinematic maps. In general LARS galaxies are characterised by high intrinsic velocity dispersions (54\,km\,s$^{-1}$ median) and low shearing velocities (65\,km\,s$^{-1}$ median). $v_\mathrm{shear}/\sigma_0$ values range from 0.5 to 3.2 with an average of 1.5. Noteworthy, five galaxies of the sample are dispersion dominated systems with $v_\mathrm{shear}/\sigma_0 <1$ and are thus kinematically similar to turbulent star forming galaxies seen at high redshift. When linking our kinematical statistics to the global LARS Ly$\alpha$ properties, we find that dispersion dominated systems show higher Ly$\alpha$ equivalent widths and higher Ly$\alpha$ escape fractions than systems with $v_\mathrm{shear}/\sigma_0 > 1$. Our result indicates that turbulence in actively star-forming systems is causally connected to interstellar medium conditions that favour an escape of Ly$\alpha$ radiation.

We present a study of the spatial distribution of the stellar cluster populations in the star forming galaxy NGC 628. Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey), we have identified 1392 potential young (<100 Myr) stellar clusters within the galaxy, identified from a combination of visual inspection and automatic selection. We investigate the clustering of these young stellar clusters and quantify the strength and change of clustering strength with scale using the two-point correlation function. We also investigate how image boundary conditions and dust lanes affect the observed clustering. The distribution of the clusters is well fit by a broken power law with negative exponent $\alpha$. We recover a weighted mean index of $\alpha$ ~ -0.8 for all spatial scales below the break at 3".3 (158 pc at a distance of 9.9 Mpc) and an index of $\alpha$ ~ -0.18 above 158 pc for the accumulation of all cluster types. The strength of the clustering increases with decreasing age and clusters older than 40 Myr lose their clustered structure very rapidly and tend to be randomly distributed in this galaxy whereas the mass of the star cluster has little effect on the clustering strength. This is consistent with results from other studies that the morphological hierarchy in stellar clustering resembles the same hierarchy as the turbulent interstellar medium.

The stellar content of young massive star clusters emit large amounts of Lyman continuum photons and inject momentum into the inter stellar medium (ISM) by the strong stellar winds of the most massive stars in the cluster. When the most massive stars explode as supernovae, large amounts of mechanical energy are injected in the ISM. A detailed study of the ISM around these massive cluster provides insights on the effect of cluster feedback. We present high quality integral field spectroscopy taken with VLT/MUSE of two starburst galaxies: ESO 338-IG04 and Haro 11. Both galaxies contain a significant number of super star clusters. The MUSE data provide us with an unprecedented view of the state and kinematics of the ionized gas in the galaxy allowing us to study the effect of stellar feedback on small and large spatial scales. We present our recent results on studying the ISM state of these two galaxies. The data of both galaxies show that the mechanical and ionization feedback of the super star clusters in the galaxy modify the state and kinematics of the ISM substancially by creating highly ionized bubbles around the cluster, making the central part of the galaxy highly ionized. This shows that the HII regions around the individual clusters are density bounded, allowing the ionizing photons to escape and ionize the ISM further out.

We investigate a novel Bayesian analysis method, based on the Stochastically Lighting Up Galaxies (slug) code, to derive the masses, ages, and extinctions of star clusters from integrated light photometry. Unlike many analysis methods, slug correctly accounts for incomplete IMF sampling, and returns full posterior probability distributions rather than simply probability maxima. We apply our technique to 621 visually-confirmed clusters in two nearby galaxies, NGC 628 and NGC 7793, that are part of the Legacy Extragalactic UV Survey (LEGUS). LEGUS provides Hubble Space Telescope photometry in the NUV, U, B, V, and I bands. We analyze the sensitivity of the derived cluster properties to choices of prior probability distribution, evolutionary tracks, IMF, metallicity, treatment of nebular emission, and extinction curve. We find that slug's results for individual clusters are insensitive to most of these choices, but that the posterior probability distributions we derive are often quite broad, and sometimes multi-peaked and quite sensitive to the choice of priors. In contrast, the properties of the cluster population as a whole are relatively robust against all of these choices. We also compare our results from slug to those derived with a conventional non-stochastic fitting code, Yggdrasil. We show that slug's stochastic models are generally a better fit to the observations than the deterministic ones used by Yggdrasil. However, the overall properties of the cluster populations recovered by both codes are qualitatively similar.

The nearby dwarf starburst galaxy NGC5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the `radio nebula'). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC5253 with wavelength coverage from 1500 Ang to 1.9 micron in 13 filters. These include H-alpha, P-beta, and P-alpha, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the 9 optically brightest clusters (M_V < -8.8) and the two young radio nebula clusters. The clusters have ages ~1-15 Myr and masses ~1x10^4 - 2.5x10^5 M_sun. The clusters' spatial location and ages indicate that star formation has become more concentrated towards the radio nebula over the last ~15 Myr. The most massive cluster is in the radio nebula; with a mass 2.5x10^5 M_sun and an age ~1 Myr, it is 2-4 times less massive and younger than previously estimated. It is within a dust cloud with A_V~50 mag, and shows a clear nearIR excess, likely from hot dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars.

(abridged) Luminous blue compact galaxies are among the most active galaxies in the local universe in terms of their star formation rate per unit mass. They may be seen as the local analogs of higher redshift Lyman Break Galaxies. Studies of their kinematics is key to understanding what triggers their unusually active star formation In this work we investigate the kinematics of stars and ionised gas in Haro11, one of the most luminous blue compact galaxies in the local universe. Previous works have indicated that many such galaxies may be triggered by galaxy mergers. We have employed Fabry-Perot interferometry, long-slit spectroscopy and Integral Field Unit (IFU) spectroscopy to explore the kinematics of Haro11. We target the near infrared Calcium triplet to derive the stellar velocity field and velocity dispersion. Ionised gas is analysed through emission lines from hydrogen, [OIII] , and [SIII]. When spectral resolution and signal to noise allows we investigate the the line profile in detail and identify multiple velocity components when present. We find that to first order, the velocity field and velocity dispersions derived from stars and ionised gas agree. Hence the complexities reveal real dynamical disturbances providing further evidence for a merger in Haro11. Through decomposition of emission lines we find evidence for kinematically distinct components, for instance a tidal arm behind the galaxy. The ionised gas velocity field can be traced to large galactocentric radii, and shows significant velocity dispersion even far out in the halo. We discuss the origin of the line width, and interpreted as virial motions it indicates a mass of ~1E11 M_sun. Morphologically and kinematically Haro11 shows many resemblances with the famous Antennae galaxies, but is much denser which is the likely explanation for the higher star formation efficiency in Haro11.

We present a study of the M83 cluster population, covering the disc of the galaxy between radii of 0.45 and 4.5 kpc. We aim to probe the properties of the cluster population as a function of distance from the galactic centre. We observe a net decline in cluster formation efficiency ($\Gamma$, i.e. the amount of star formation happening in bound clusters) from about 26 % in the inner region to 8 % in the outer part of the galaxy. The recovered $\Gamma$ values within different regions of M83 follow the same $\Gamma$ versus star formation rate density relation observed for entire galaxies. We also probe the initial cluster mass function (ICMF) as a function of galactocentric distance. We observe a significant steepening of the ICMF in the outer regions (from $-1.90\pm0.11$ to $-2.70\pm0.14$) and for the whole galactic cluster population (slope of $-2.18\pm0.07$) of M83. We show that this change of slope reflects a more fundamental change of the 'truncation mass' at the high-mass end of the distribution. This can be modelled as a Schechter function of slope $-2$ with an exponential cut-off mass ($M_{\rm c}$) that decreases significantly from the inner to the outer regions (from 4.00 to $0.25\times 10^5$ M$_\odot$) while the galactic $M_{\rm c}$ is $\approx1.60\times10^5$ M$_\odot$. The trends in \Gamma and ICMF are consistent with the observed radial decrease of the $\Sigma({\rm H}_2)$, hence in gas pressure. As gas pressure declines cluster formation becomes less efficient. We conclude that the host galaxy environment appears to regulate 1) the fraction of stars locked in clusters; 2) the upper mass limit of the ICMF, consistently described by a near-universal slope $-2$ truncated at the high-mass end.

We present a detailed clustering analysis of the young stellar population across the star-forming ring galaxy NGC 6503, based on the deep HST photometry obtained with the Legacy ExtraGalactic UV Survey (LEGUS). We apply a contour-based map analysis technique and identify in the stellar surface density map 244 distinct star-forming structures at various levels of significance. These stellar complexes are found to be organized in a hierarchical fashion with 95% being members of three dominant super-structures located along the star-forming ring. The size distribution of the identified structures and the correlation between their radii and numbers of stellar members show power-law behaviors, as expected from scale-free processes. The self-similar distribution of young stars is further quantified from their autocorrelation function, with a fractal dimension of ~1.7 for length-scales between ~20 pc and 2.5 kpc. The young stellar radial distribution sets the extent of the star-forming ring at radial distances between 1 and 2.5 kpc. About 60% of the young stars belong to the detected stellar structures, while the remaining stars are distributed among the complexes, still inside the ring of the galaxy. The analysis of the time-dependent clustering of young populations shows a significant change from a more clustered to a more distributed behavior in a time-scale of ~60 Myr. The observed hierarchy in stellar clustering is consistent with star formation being regulated by turbulence across the ring. The rotational velocity difference between the edges of the ring suggests shear as the driving mechanism for this process. Our findings reveal the interesting case of an inner ring forming stars in a hierarchical fashion.

We measure the radii and two-dimensional light profiles of a large sample of young, massive star clusters in M83 using archival HST/WFC3 imaging of seven adjacent fields. We use GALFIT to fit the two-dimensional light profiles of the clusters, from which we find effective (half-light) radii, core radii, and slopes of the power-law (EFF) profile ($\eta$). We find lognormal distributions of effective radius and core radius, with medians of $\approx$2.5 pc and $\approx$1.3 pc, respectively. Our results provide strong evidence for a characteristic size of young, massive clusters. The average effective radius and core radius increase somewhat with cluster age. Little to no change in effective radius is observed with increasing galactocentric distance, except perhaps for clusters younger than 100 Myr. We find a shallow correlation between effective radius and mass for the full cluster sample, but a stronger correlation is present for clusters 200-300 Myr in age. Finally, the majority of the clusters are best fit by an EFF model with index $\eta\leq3.0$. There is no strong evidence for change in $\eta$ with cluster age, mass, or galactocentric distance. Our results suggest that clusters emerge from early evolution with similar radii and are not strongly affected by the tidal field of M83. Mass loss due to stellar evolution and/or GMC interactions appear to dominate cluster expansion in the age range we study.

The Ly$\alpha$ line is an important diagnostic for star formation at high redshift, but interpreting its flux and line profile is difficult because of the resonance nature of Ly$\alpha$. Trends between the escape of Ly$\alpha$ photons and dust and properties of the interstellar medium (ISM) have been found, but detailed comparisons between Ly$\alpha$ emission and the properties of the gas in local high-redshift analogs are vital for understanding the relation between Ly$\alpha$ emission and galaxy properties. For the first time, we can directly infer the properties of the ionized gas at the same location and similar spatial scales of the extended Ly$\alpha$ halo around ESO 338-IG04. We obtained VLT/MUSE integral field spectra. We used ionization parameter mapping of the [SII]/[OIII] line ratio and the kinematics of H$\alpha$ to study the ionization state and kinematics of the ISM of ESO338-IG04. The velocity map reveals two outflows. The entire central area of the galaxy is highly ionized by photons leaking from the HII regions around the youngest star clusters. Three highly ionized cones have been identified, of which one is associated with an outflow. We propose a scenario where the outflows are created by mechanical feedback of the older clusters, while the highly ionized gas is caused by the hard ionizing photons emitted by the youngest clusters. A comparison with the Ly$\alpha$ map shows that the (approximately bipolar) asymmetries observed in the Ly$\alpha$ emission are consistent with the base of the outflows detected in H$\alpha$. No clear correlation with the ionization cones is found. The mechanical and ionization feedback of star clusters significantly changes the state of the ISM by creating ionized cones and outflows. The comparison with Ly$\alpha$ suggests that especially the outflows could facilitate the escape of Ly$\alpha$ photons [Abridged].

We present high-resolution far-UV spectroscopy of the 14 galaxies of the Lyman Alpha Reference Sample; a sample of strongly star-forming galaxies at low redshifts ($0.028 < z < 0.18$). We compare the derived properties to global properties derived from multi band imaging and 21 cm HI interferometry and single dish observations, as well as archival optical SDSS spectra. Besides the Lyman $\alpha$ line, the spectra contain a number of metal absorption features allowing us to probe the kinematics of the neutral ISM and evaluate the optical depth and and covering fraction of the neutral medium as a function of line-of-sight velocity. Furthermore, we show how this, in combination with precise determination of systemic velocity and good Ly$\alpha$ spectra, can be used to distinguish a model in which separate clumps together fully cover the background source, from the "picket fence" model named by Heckman et al. (2011). We find that no one single effect dominates in governing Ly$\alpha$ radiative transfer and escape. Ly$\alpha$ escape in our sample coincides with a maximum velocity-binned covering fraction of $\lesssim 0.9$ and bulk outflow velocities of $\gtrsim 50$ km s$^{-1}$, although a number of galaxies show these characteristics and yet little or no Ly$\alpha$ escape. We find that Ly$\alpha$ peak velocities, where available, are not consistent with a strong backscattered component, but rather with a simpler model of an intrinsic emission line overlaid by a blueshifted absorption profile from the outflowing wind. Finally, we find a strong anticorrelation between H$\alpha$ equivalent width and maximum velocity-binned covering factor, and propose a heuristic explanatory model.

The study of young massive clusters can provide key information for the formation of globular clusters, as they are often considered analogues. A currently unanswered question in this field is how long these massive clusters remain embedded in their natal gas, with important implications for the formation of multiple populations that have been used to explain phenomena observed in globular clusters. We present an analysis of ages and masses of the young massive cluster population of M83. Through visual inspection of the clusters, and comparison of their SEDs and position in colour-colour space, the clusters are all exposed (no longer embedded) by < 4 Myr, most likely less, indicating that current proposed age spreads within older clusters are unlikely. We also present several methods of constraining the ages of very young massive clusters. This can often be difficult using SED fitting due to a lack of information to disentangle age-extinction degeneracies and possible inaccurate assumptions in the models used for the fitting. The individual morphology of the Halpha around each cluster has a significant effect on the measured fluxes, which contributes to inaccuracies in the age estimates for clusters younger than 10 Myr using SED fitting. This is due to model uncertainties and aperture effects. Our methods to help constrain ages of young clusters include using the near-infrared and spectral features, such as Wolf-Rayet stars.

The Legacy ExtraGalactic UV Survey (LEGUS) is a Cycle 21 Treasury program on the Hubble Space Telescope, aimed at the investigation of star formation and its relation with galactic environment in nearby galaxies, from the scales of individual stars to those of ~kpc-size clustered structures. Five-band imaging, from the near-ultraviolet to the I-band, with the Wide Field Camera 3, plus parallel optical imaging with the Advanced Camera for Surveys, is being collected for selected pointings of 50 galaxies within the local 12 Mpc. The filters used for the observations with the Wide Field Camera 3 are: F275W(2,704 A), F336W(3,355 A), F438W(4,325 A), F555W(5,308 A), and F814W(8,024 A); the parallel observations with the Advanced Camera for Surveys use the filters: F435W(4,328 A), F606W(5,921 A), and F814W(8,057 A). The multi-band images are yielding accurate recent (<~50 Myr) star formation histories from resolved massive stars and the extinction-corrected ages and masses of star clusters and associations. The extensive inventories of massive stars and clustered systems will be used to investigate the spatial and temporal evolution of star formation within galaxies. This will, in turn, inform theories of galaxy evolution and improve the understanding of the physical underpinning of the gas-star formation relation and the nature of star formation at high redshift. This paper describes the survey, its goals and observational strategy, and the initial science results. Because LEGUS will provide a reference survey and a foundation for future observations with JWST and with ALMA, a large number of data products are planned for delivery to the community.

The Lyman-alpha reference sample (LARS) is a program with the Hubble Space Telescope (HST) that provides a sample of local universe laboratory galaxies in which to study the astrophysics of the visibility and strength of the Lyman-alpha (Lya) line of hydrogen. This article presents an overview of the survey, its selection function and HST imaging observations. The sample was selected from the GALEX+SDSS catalogue at z=0.028-0.19, in order to allow Lya to be captured with combinations of long pass filters in the Solar Blind Channel (SBC) of HST/ACS. In addition, LARS utilises Halpha and Hbeta narrow, and U, B, i broad-band imaging with ACS and WFC3. In order to study galaxies in which large numbers of Lya photons are produced we demanded an Halpha equivalent width > 100{\AA}. The sample of 14 galaxies covers far UV (FUV) luminosities that overlaps with those of high-z Lya emitters and Lyman Break Galaxies (LBGs). We present the reduction steps used to obtain the Lya images, including our LARS eXtraction software (LaXs) which utilises pixel-by-pixel spectral synthesis fitting to subtract the continuum at Lya. We demonstrate that the use of SBC long pass filters dramatically increase the signal to noise compared to the nominal Lya filter in SBC. To exemplify, we present results for Mrk 259. This irregular galaxy shows extended but strongly asymmetric Lya emission. Spectroscopy from HST/COS show a moderate outflow in the neutral ISM and an asymmetric Lya profile. Radiative transfer modeling is able to reproduce the Lya line profile, confirming the presence of an outflow. The integrated photometry provides a Lya luminosity of 1.3E42 erg/s an equivalent width W(Lya) =45{\AA} and a FUV absolute magnitude -19.2. Mrk 259 would hence be detectable in high-z Lya and LBG surveys. The total Lya escape fraction is 12%, similar to that observed at z>4. (shortened)

We present new H I imaging and spectroscopy of the 14 UV-selected star-forming galaxies in the Lyman Alpha Reference Sample (LARS), aimed for a detailed study of the processes governing the production, propagation, and escape of Ly$\alpha$ photons. New H I spectroscopy, obtained with the 100m Green Bank Telescope (GBT), robustly detects the H I spectral line in 11 of the 14 observed LARS galaxies (although the profiles of two of the galaxies are likely confused by other sources within the GBT beam); the three highest redshift galaxies are not detected at our current sensitivity limits. The GBT profiles are used to derive fundamental H I line properties of the LARS galaxies. We also present new pilot H I spectral line imaging of 5 of the LARS galaxies obtained with the Karl G. Jansky Very Large Array (VLA). This imaging localizes the H I gas and provides a measurement of the total H I mass in each galaxy. In one system, LARS 03 (UGC 8335 or Arp 238), VLA observations reveal an enormous tidal structure that extends over 160 kpc from the main interacting systems and that contains $>$10$^9$ M$_{\odot}$ of H I. We compare various H I properties with global Ly$\alpha$ quantities derived from HST measurements. The measurements of the Ly$\alpha$ escape fraction are coupled with the new direct measurements of H I mass and significantly disturbed H I velocities. Our robustly detected sample reveals that both total H I mass and linewidth are tentatively correlated with key Ly$\alpha$ tracers. Further, on global scales, these data support a complex coupling between Ly$\alpha$ propagation and the H I properties of the surrounding medium.

Recent ground based near-IR studies of stellar clusters in nearby galaxies have suggested that young clusters remain embedded for 7-10Myr in their progenitor molecular cloud, in conflict with optical based studies which find that clusters are exposed after 1-3Myr. Here, we investigate the role that spatial resolution plays in this apparent conflict. We use a recent catalogue of young ($<10$~Myr) massive ($>5000$~\msun) clusters in the nearby spiral galaxy, M83, along with Hubble Space Telescope (HST) imaging in the optical and near-IR, and ground based near-IR imaging, to see how the colours (and hence estimated properties such as age and extinction) are affected by the aperture size employed, in order to simulate studies of differing resolution. We find that the near-IR is heavily affected by the resolution, and when aperture sizes $>40$~pc are used, all young/blue clusters move red-ward in colour space, which results in their appearance as heavily extincted clusters. However, this is due to contamination from nearby sources and nebular emission, and is not an extinction effect. Optical colours are much less affected by resolution. Due to the larger affect of contamination in the near-IR, we find that, in some cases, clusters will appear to show near-IR excess when large ($>20$~pc) apertures are used. Our results explain why few young ($<6$~Myr), low extinction ($\av < 1$~mag) clusters have been found in recent ground based near-IR studies of cluster populations, while many such clusters have been found in higher resolution HST based studies. Additionally, resolution effects appear to (at least partially) explain the origin of the near-IR excess that has been found in a number of extragalactic YMCs.

We study the star cluster population of NGC 2997, a giant spiral galaxy located at 9.5 Mpc and targeted by the Snapshot Hubble U-band Cluster Survey (SHUCS). Combining our U-band imaging from SHUCS with archival BVI imaging from HST, we select a high confidence sample of clusters in the circumnuclear ring and disk through a combination of automatic detection procedures and visual inspection. The cluster luminosity functions in all four filters can be approximated by power-laws with indices of $-1.7$ to $-2.3$. Some deviations from pure power-law shape are observed, hinting at the presence of a high-mass truncation in the cluster mass function. However, upon inspection of the cluster mass function, we find it is consistent with a pure power-law of index $-2.2\pm0.2$ despite a slight bend at $\sim$$2.5\times10^{4}$ M$_{\odot}$. No statistically significant truncation is observed. From the cluster age distributions, we find a low rate of disruption ($\zeta\sim-0.1$) in both the disk and circumnuclear ring. Finally, we estimate the cluster formation efficiency ($\Gamma$) over the last 100 Myr in each region, finding $7\pm2$% for the disk, $12\pm4$% for the circumnuclear ring, and $10\pm3$% for the entire UBVI footprint. This study highlights the need for wide-field UBVI coverage of galaxies to study cluster populations in detail, though a small sample of clusters can provide significant insight into the characteristics of the population.

Hierarchical structure in ultraviolet images of 12 late-type LEGUS galaxies is studied by determining the numbers and fluxes of nested regions as a function of size from ~1 to ~200 pc, and the number as a function of flux. Two starburst dwarfs, NGC 1705 and NGC 5253, have steeper number-size and flux-size distributions than the others, indicating high fractions of the projected areas filled with star formation. Nine subregions in 7 galaxies have similarly steep number-size slopes, even when the whole galaxies have shallower slopes. The results suggest that hierarchically structured star-forming regions several hundred parsecs or larger represent common unit structures. Small galaxies dominated by only a few of these units tend to be starbursts. The self-similarity of young stellar structures down to parsec scales suggests that star clusters form in the densest parts of a turbulent medium that also forms loose stellar groupings on larger scales. The presence of super star clusters in two of our starburst dwarfs would follow from the observed structure if cloud and stellar subregions more readily coalesce when self-gravity in the unit cell contributes more to the total gravitational potential.

In order to empirically determine the timescale and environmental dependence of stellar cluster disruption, we have undertaken an analysis of the unprecedented multi-pointing (seven), multi-wavelength (U, B, V, H$\alpha$, and I) Hubble Space Telescope imaging survey of the nearby, face-on spiral galaxy M83. The images are used to locate stellar clusters and stellar associations throughout the galaxy. Estimation of cluster properties (age, mass, and extinction) was done through a comparison of their spectral energy distributions with simple stellar population models. We constructed the largest catalog of stellar clusters and associations in this galaxy to-date, with ~1800 sources with masses above ~5000 M$_{\odot}$ and ages younger than ~300 Myr. In the present letter, we focus on the age distribution of the resulting clusters and associations. In particular, we explicitly test whether the age distributions are related with the ambient environment. Our results are in excellent agreement with previous studies of age distributions in the centre of the galaxy, which gives us confidence to expand out to search for similarities or differences in the other fields which sample different environments. We find that the age distribution of the clusters inside M83 varies strongly as a function of position within the galaxy, indicating a strong correlation with the galactic environment. If the age distributions are approximated as a power-law of the form $dN/dt \propto t^{\zeta}$, we find $\zeta$ values between 0 and -0.62 ( $\zeta$~-0.40 for the whole galaxy), in good agreement with previous results and theoretical predictions.

Embedded clusters like W3 Main are complex and dynamically evolving systems that represent an important phase of the star formation process. We aim at the characterization of the entire stellar content of W3 Main in a statistical sense to identify possible differences in evolutionary phase of the stellar populations and find clues about the formation mechanism of this massive embedded cluster. Methods. Deep JHKs imaging is used to derive the disk fraction, Ks-band luminosity functions and mass functions for several subregions in W3 Main. A two dimensional completeness analysis using artificial star experiments is applied as a crucial ingredient to assess realistic completeness limits for our photometry. We find an overall disk fraction of 7.7 $\pm$ 2.3%, radially varying from 9.4 $\pm$ 3.0 % in the central 1 pc to 5.6 $\pm$ 2.2 % in the outer parts of W3 Main. The mass functions derived for three subregions are consistent with a Kroupa and Chabrier mass function. The mass function of IRSN3 is complete down to 0.14 Msun and shows a break at M $\sim$ 0.5 Msun. We interpret the higher disk fraction in the center as evidence for a younger age of the cluster center. We find that the evolutionary sequence observed in the low-mass stellar population is consistent with the observed age spread among the massive stars. An analysis of the mass function variations does not show evidence for mass segregation. W3 Main is currently still actively forming stars, showing that the ionizing feedback of OB stars is confined to small areas ($\sim$ 0.5 pc). The FUV feedback might be influencing large regions of the cluster as suggested by the low overall disk fraction.

We report upon new results regarding the Lya output of galaxies, derived from the Lyman alpha Reference Sample (LARS), focusing on Hubble Space Telescope imaging. For 14 galaxies we present intensity images in Lya, Halpha, and UV, and maps of Halpha/Hbeta, Lya equivalent width (EW), and Lya/Halpha. We present Lya and UV light profiles and show they are well-fitted by S\'ersic profiles, but Lya profiles show indices systematically lower than those of the UV (n approx 1-2 instead of >~4). This reveals a general lack of the central concentration in Lya that is ubiquitous in the UV. Photometric growth curves increase more slowly for Lya than the FUV, showing that small apertures may underestimate the EW. For most galaxies, however, flux and EW curves flatten by radii ~10 kpc, suggesting that if placed at high-z, only a few of our galaxies would suffer from large flux losses. We compute global properties of the sample in large apertures, and show total luminosities to be independent of all other quantities. Normalized Lya throughput, however, shows significant correlations: escape is found to be higher in galaxies of lower star formation rate, dust content, mass, and several quantities that suggest harder ionizing continuum and lower metallicity. Eight galaxies could be selected as high-z Lya emitters, based upon their luminosity and EW. We discuss the results in the context of high-z Lya and UV samples. A few galaxies have EWs above 50 AA, and one shows f_escLya of 80%; such objects have not previously been reported at low-z.

Recent observations, as well as theoretical studies, have suggested that stellar cluster formation may depend on local and global environmental properties. In particular, the fraction of stars that form within long-lived bound clusters (\Gamma) may depend on environment, with indications that it may be higher in the more extreme environments of high star-formation rate density galaxies. How \Gamma varies has important implications on the use of clusters to determine the star-formation histories of galaxies as well as our understanding of the star-formation process itself. Previous studies have estimated \Gamma over full galaxies, making it difficult to discern the physical cause of the reported variations. Here, we use existing star cluster catalogues and HST-WFC3 V and I images of the grand design, face-on spiral galaxy M83 in order to see if and how \Gamma varies within a single galaxy. We find that \Gamma decreases strongly as a function of galactocentric radius, by a factor of ~5 over the inner ~6kpc, in agreement with recent theoretical predictions and decreasing trends observed in the gas surface density of the galaxy.

We report on new imaging observations of the Lyman alpha emission line (Lya), performed with the Hubble Space Telescope, that comprise the backbone of the Lyman alpha Reference Sample (LARS). We present images of 14 starburst galaxies at redshifts 0.028 < z < 0.18 in continuum-subtracted Lya, Halpha, and the far ultraviolet continuum. We show that Lya is emitted on scales that systematically exceed those of the massive stellar population and recombination nebulae: as measured by the Petrosian 20 percent radius, RP20, Lya radii are larger than those of Halpha by factors ranging from 1 to 3.6, with an average of 2.4. The average ratio of Lya-to-FUV radii is 2.9. This suggests that much of the Lya light is pushed to large radii by resonance scattering. Defining the "Relative Petrosian Extension" of Lya compared to Halpha, \xi_ext = RP20_Lya / RP20_Ha, we find \xi_ext to be uncorrelated with total Lya luminosity. However \xi_ext is strongly correlated with quantities that scale with dust content, in the sense that a low dust abundance is a necessary requirement (although not the only one) in order to spread Lya photons throughout the interstellar medium and drive a large extended Lya halo.

We present the Snapshot Hubble U-band Cluster Survey (SHUCS), a project aimed at characterizing the star cluster populations of ten nearby galaxies (d<23 Mpc, half within 12 Mpc) through new F336W (U band equivalent) imaging from WFC3, and archival BVI-equivalent data with HST. Completing the UBVI baseline reduces the age-extinction degeneracy of optical colours, thus enabling the measurement of reliable ages and masses for the thousands of clusters covered by our survey. The sample consists chiefly of face-on spiral galaxies at low inclination, in various degrees of isolation (isolated, in group, merging), and includes two AGN hosts. This first paper outlines the survey itself, the observational datasets, the analysis methods, and presents a proof-of-concept study of the large-scale properties and star cluster population of NGC 4041, a massive SAbc galaxy at a distance of 23 Mpc, and part of a small grouping of six giant members. We resolve two structural components with distinct stellar populations, a morphology more akin to merging and interacting systems. We also find strong evidence of a truncated, Schechter-type mass function, and a similarly segmented luminosity function. These results indicate that binning must erase much of the substructure present in the mass and luminosity functions, and might account for the conflicting reports on the intrinsic shape of these functions in the literature. We also note a tidal feature in the outskirts of the galaxy in GALEX UV imaging, and follow it up with a comprehensive multi-wavelength study of NGC 4041 and its parent group. We deduce a minor merger as a likely cause of its segmented structure and the observed pattern of a radially decreasing star formation rate. We propose that combining the study of star cluster populations with broad-band metrics is not only advantageous, but often easily achievable through archival datasets.

We analyse the clump population of the spiral galaxy Sp 1149 at redshift 1.5. Located behind the galaxy cluster MACS J1149.5+2223, Sp 1149 has been significantly magnified allowing us to study the galaxy on physical scales down to ~100 pc. We have used the publicly available multi-band imaging dataset (CLASH) to reconstruct the spectral energy distributions (SEDs) of the clumps in Sp 1149, and derive, by means of stellar evolutionary models, their physical properties. We found that 40% of the clumps observed in Sp 1149 are older than 30 Myr and can be as old as 300 Myr. These are also the more massive (luminous) clumps in the galaxy. Among the complexes in the local reference sample, the star-forming knots in luminous blue compact galaxies could be considered progenitor analogs of these long-lived clumps. The remaining 60% of clumps have colors comparable to local cluster complexes, suggesting a similar young age. We observe that the Sp 1149 clumps follow the M ~ R^2 relation similar to local cluster complexes, suggesting similar formation mechanisms although they may have different initial conditions (e.g. higher gas surface densities). We suggest that the galaxy is experiencing a slow decline in star formation rate and a likely transitional phase toward a more quiescent star-formation mode. The older clumps have survived between 6 and 20 dynamical times and are all located at projected distances smaller than 4 kpc from the centre. Their current location suggests migration toward the centre and the possibility to be the building blocks of the bulge. On the other hand, the dynamical timescale of the younger clumps are significantly shorter, meaning that they are quite close to their birthplace. We show that the clumps of Sp 1149 may account for the expected metal-rich globular cluster population usually associated with the bulge and thick disk components of local spirals.

The determination of age is a critical component in the study of a population of stellar clusters. In this letter we present a new absolute age indicator for young massive star clusters based on J-H colour. This novel method identifies clusters as older or younger than 5.7 +/- 0.8 Myr based on the appearance of the first population of red supergiant stars. We test the technique on the stellar cluster population of the nearby spiral galaxy, M83, finding good agreement with the theoretical predictions. The localisation of this technique to the near-IR promises that it may be used well into the future with space-- and ground--based missions optimised for near-IR observations.

We present the Snapshot Hubble U-band Cluster Survey (SHUCS), an ongoing deep U-band imaging survey of nearby star-forming galaxies. Thanks to the information provided by the U band, together with archival Hubble Space Telescope (HST) optical data, we are able to constrain reliable ages, masses, and extinctions of the cluster populations of these galaxies. We show some preliminary results from the study of one of the SHUCS galaxies, NGC 2146. Using the recovered cluster ages we try to understand the propagation of the star formation in one of the tidal streams where a ring-like cluster complex has been found. The Ruby Ring, so named due to its appearance, shows a clear ring-like distribution of star clusters around a central object. We find evidence of a spatial and temporal correlation between the central cluster and the clusters in the ring. The Ruby Ring is the product of an intense and localised burst of star formation, similar to the extended cluster complexes observed in M 51 and the Antennae, but more impressive because is quite isolated. We discuss the formation of the Ruby Ring in a "collect & collapse" framework. The predictions made by this model agree quite well with the estimated bubble radius and expansion velocity produced by the feedback from the central cluster, making the Ruby Ring an interesting case of triggered star formation.

We report the discovery of a ring-like cluster complex in the starburst galaxy NGC 2146. The Ruby Ring, so named due to its appearance, shows a clear ring-like distribution of star clusters around a central object. It is located in one of the tidal streams which surround the galaxy. NGC 2146 is part of the Snapshot Hubble U-band Cluster Survey (SHUCS). The WFC3/F336W data has added critical information to the available archival Hubble Space Telescope imaging set of NGC 2146, allowing us to determine ages, masses, and extinctions of the clusters in the Ruby Ring. These properties have then been used to investigate the formation of this extraordinary system. We find evidence of a spatial and temporal correlation between the central cluster and the clusters in the ring. The latter are about 4 Myr younger than the central cluster, which has an age of 7 Myr. This result is supported by the H alpha emission which is strongly coincident with the ring, and weaker at the position of the central cluster. From the derived total H alpha luminosity of the system we constrain the star formation rate density to be quite high, e.g. ~ 0.47 Msun/yr/kpc^2. The Ruby Ring is the product of an intense and localised burst of star formation, similar to the extended cluster complexes observed in M51 and the Antennae, but more impressive because is quite isolated. The central cluster contains only 5 % of the total stellar mass in the clusters that are determined within the complex. The ring-like morphology, the age spread, and the mass ratio support a triggering formation scenario for this complex. We discuss the formation of the Ruby Ring in a "collect & collapse" framework. The predictions made by this model agree quite well with the estimated bubble radius and expansion velocity produced by the feedback from the central cluster, making the Ruby Ring an interesting case of triggered star formation.

We study the stellar cluster population in two adjacent fields in the nearby, face-on spiral galaxy, M83, using WFC3/HST imaging. The clusters are selected through visual inspection to be centrally concentrated, symmetric, and resolved on the images, which allows us to differentiate between clusters and likely unbound associations. We compare our sample with previous studies and show that the differences between the catalogues are largely due to the inclusion of large numbers of diffuse associations within previous catalogues. The luminosity function of the clusters is well approximated by a power-law with index, -2, over most of the observed range, however a steepening is seen at M_V = -9.3 and -8.8 in the inner and outer fields, respectively. Additionally, we show that the cluster population is inconsistent with a pure power-law mass distribution, but instead exhibits a truncation at the high mass end. If described as a Schechter function, the characteristic mass is 1.6 and 0.5 * 10^5 Msun, for the inner and outer fields, respectively, in agreement with previous estimates of other cluster populations in spiral galaxies. Comparing the predictions of the mass independent disruption (MID) and mass dependent disruption (MDD) scenarios with the observed distributions, we find that both models can accurately fit the data. However, for the MID case, the fraction of clusters destroyed (or mass lost) per decade in age is dependent on the environment, hence, the age/mass distributions of clusters are not universal. In the MDD case, the disruption timescale scales with galactocentric distance (being longer in the outer regions of the galaxy) in agreement with analytic and numerical predictions. Finally, we discuss the implications of our results on other extragalactic surveys, focussing on the fraction of stars that form in clusters and the need (or lack thereof) for infant mortality.

The numerous and massive young star clusters in blue compact galaxies (BCGs) are used to investigate the properties of their hosts. We test whether BCGs follow claimed relations between cluster populations and their hosts, such as the the fraction of the total luminosity contributed by the clusters as function of the mean star formation rate density; the $V$ band luminosity of the brightest youngest cluster as related to the mean host star formation rate; and the cluster formation efficiency (i.e., the fraction of star formation happening in star clusters) versus the density of the SFR. We find that BCGs follow the trends, supporting a scenario where cluster formation and environmental properties of the host are correlated. They occupy, in all the diagrams, the regions of higher SFRs, as expected by the extreme nature of the starbursts operating in these systems. We find that the star clusters contribute almost to the 20 % of the UV luminosity of the hosts. We suggest that the BCG starburst environment has most likely favoured the compression and collapse of the giant molecular clouds, enhancing the local star formation efficiency, so that massive clusters have been formed. The estimated cluster formation efficiency supports this scenario. BCGs have a cluster formation efficiency comparable to luminous IR galaxies and spiral starburst nuclei (the averaged value is about 35 %) which is much higher than the 8 - 10 % reported for quiescent spirals and dwarf star-forming galaxies.

We report the results of a recent study of young star clusters (YSCs) in luminous blue compact galaxies (BCGs). The age distributions of the YSCs suggest that the starburst episode in Haro 11, ESO 185-IG13, and Mrk 930 started not more than 30-40 Myr ago. A peak of cluster formation only 3 - 4 Myr old is observed, unveiling a unique sample of clusters still partially embedded. A considerable fraction of clusters (30 - 50%), mainly younger than 10 Myr, shows an observed flux excess between 0.8 and 2.2 micron. This so-called near-infrared (NIR) excess is impossible to reproduce even with the most recent spectral synthesis models (that include a self-consistent treatment of the photoionized gas). We have used these YSCs to probe the very early evolution phase of star clusters. In all the three host galaxies, the analysis is limited to the optically brightest objects, i.e., systems that are only partially embedded by their natal cocoons (since deeply embedded clusters are probably too faint to be detected). We discuss possible explanations for this NIR excess, in the context of IR studies of both extragalactic young star clusters and resolved massive star forming regions in the Milky Way and in the nearby Magellanic Clouds.

Using multi-wavelength imaging from the Wide Field Camera 3 on the Hubble Space Telescope we study the stellar cluster populations of two adjacent fields in the nearby face-on spiral galaxy, M83. The observations cover the galactic centre and reach out to ~6 kpc, thereby spanning a large range of environmental conditions, ideal for testing empirical laws of cluster disruption. The clusters are selected by visual inspection to be centrally concentrated, symmetric, and resolved on the images. We find that a large fraction of objects detected by automated algorithms (e.g. SExtractor or Daofind) are not clusters, but rather are associations. These are likely to disperse into the field on timescales of tens of Myr due to their lower stellar densities and not due to gas expulsion (i.e. they were never gravitationally bound). We split the sample into two discrete fields (inner and outer regions of the galaxy) and search for evidence of environmentally dependent cluster disruption. Colour-colour diagrams of the clusters, when compared to simple stellar population models, already indicate that a much larger fraction of the clusters in the outer field are older by tens of Myr than in the inner field. This impression is quantified by estimating each cluster's properties (age, mass, and extinction) and comparing the age/mass distributions between the two fields. Our results are inconsistent with "universal" age and mass distributions of clusters, and instead show that the ambient environment strongly affects the observed populations.

We present a HST multiband analysis of the large population of star clusters in the blue compact galaxy (BCG) Mrk 930. We have reconstructed the spectral energy distributions of the star clusters and estimated age, mass, and extinction for a representative sample. We observe a very young cluster population with 70% of the systems formed less than 10 Myr ago. The peak in the star cluster age distribution at 4 Myr is corroborated by the presence of Wolf-Rayet spectral features, and by the observed optical and IR lines ratios [OIII]/H_beta and [NeIII]/[NeII]. The recovered extinction in these very young clusters shows large variations, with a decrease at older ages. It is likely that our analysis is limited to the optically brightest objects (i.e. systems only partially embedded in their natal cocoons). We map the extinction across the galaxy using low-resolution spectra and the H_alpha/H_beta ratio, as obtained from ground-based narrow band imaging. We find that the mean optical extinction derived in the starburst regions is close to the averaged value observed in the clusters, but locally, do not trace the more extinguished clusters. Previous HST studies of BCGs have revealed a population of young and extremely red super star clusters. We detect a considerable fraction of clusters affected by a red excess also in Mrk 930. The nature of the red excess, which turns up at near-IR wavelengths remains unknown. We compare the cluster and the star formation history, the latter derived from the fit of spectral population synthesis models to the spectra. We find a general agreement between the two independently estimated quantities. Using the cluster properties we perform a study of the host environmental properties. We find that the cluster formation efficiency is significantly higher, suggesting a key role of the environment for the formation of these massive objects.[Abridged]

We have investigated the starburst properties of the luminous blue compact galaxy ESO 185-IG13. The galaxy has been imaged with the high resolution cameras onboard to the Hubble Space Telescope. From the UV to the IR, the data reveal a system shaped by hundreds of young star clusters, and fine structures, like a tidal stream and a shell. The presence of numerous clusters and the perturbed morphology indicate that the galaxy has been involved in a recent merger event. Using previous simulations of shell formation in galaxy mergers we constrain potential progenitors of ESO 185-IG13. The analysis of the star cluster population is used to investigate the properties of the present starburst and to date the final merger event, which has produced hundreds of clusters younger than 100 Myr. We have found a peak of cluster formation only 3.5 Myr old. A large fraction of these clusters will not survive after 10-20 Myr, due to the "infant mortality" caused by gas expulsion. However, this sample of clusters represents an unique chance to investigate the youngest phases of cluster evolution. As already observed in the analog blue compact galaxy Haro 11, a fraction of young clusters are affected by a flux excess at wavelengths longer than 8000 \AA. Ages, masses, and extinctions of clusters with this NIR excess are estimated from UV and optical data. We discuss similarities and differences of the observed NIR excess in ESO 185-IG13 clusters with other cases in the literature. The cluster ages and masses are used to distinguish among the potential causes of the excess. We observe, as in Haro 11, that the use of the IR and the (commonly used) I band data results in overestimates of age and mass in clusters affected by the NIR excess. This has important implications for a number of related studies of star clusters.

The spectral energy distribution analysis of very young unresolved star clusters challenges our understanding of the cluster formation process. Studies of resolved massive clusters in the Milky Way and in the nearby Magellanic Clouds show us that the contribution from photoionized gas is very important during the first Myr of cluster evolution. We present our models which include both a self-consistent treatment of the photoionized gas and the stellar continuum and quantify the impact of such nebular component on the total flux of young unresolved star clusters. A comparison with other available models is considered. The very young star clusters in the SBS 0335-052E dwarf starburst galaxy are used as a test for our models. Due to the low metallicity of the galactic medium our models predict a longer lasted nebular phase which contributes between 10-40% of the total near infrared (NIR) fluxes at around 10 Myr. We propose thus a possible solution for the observed flux excess in the 6 bright super star clusters of SBS 0335-052E. Reines et al. showed that the observed cluster fluxes, in the red-optical and NIR range, sit irreconcilably above the provided stellar continuum models. We find that in the age range estimated from the H_alpha emission we can explain the red excess in all the 6 super star clusters as due to nebular emission, which at cluster ages around 10 Myr still affects the NIR wavebands substantially.

We have used multi-band imaging to investigate the nature of the extreme starburst environment in Haro 11 galaxy. The central starburst region has been observed in 8 HST wavebands and at 2.16 micron at the ESO-VLT. We constructed integrated spectral energy distributions (SEDs) for about 200 star clusters and compared them with single stellar population models in order to derive ages, masses and extinctions of thestar clusters. The present starburst has lasted for 40 Myr, and shows a peak of cluster formation only 3.5 Myr old. With such an extremely young cluster population, Haro 11 represents a unique opportunity to investigate the youngest phase of the cluster formation process and evolution in starburst systems. Extinction tends to diminish as function of the cluster age, but the spread is large and for clusters in partial embedded phases (< 5 Myr). A fraction of low-mass (> 10^4 Msun), very young (1-3 Myr) clusters is missing, either because they are embedded inthe parental molecular cloud and heavily extinguished, or because of blending. Almost half of the cluster sample is affected by flux excesses at wavelengths 8000 \AA which cannot be explained by simple stellar evolutionary models. Fitting SED models over all wavebands leads to systematic overestimates of cluster ages and incorrect masses for the stellar population supplying the light in these clusters. We show that the red excess affects also the HST F814W filter, which is typically used to constrain cluster physical properties. The clusters which show the red excess are younger than 40 Myr; we propose possible physical explanations for the phenomenon. Finally, we estimate that Haro 11 hasproduced bound clusters at a rate almost a factor of 10 higher than the massive and regular spirals, like the Milky Way. (Abriged)

We have performed a deep multi-band photometric analysis of the star cluster population of Haro 11. This starburst galaxy (log L_FUV = 10.3 L_sun) is considered a nearby analogue of Lyman break galaxies (LBGs) at high redshift. The study of the numerous star clusters in the systems is an effective way to investigate the formation and evolution of the starburst phase. In fact, the SED fitting models have revealed a surprisingly young star cluster population, with ages between 0.5 and 40 Myr, and estimated masses between 10^3 and 10^7 solar masses. An independent age estimation has been done with the EW(Halpha) analysis of each cluster. This last analysis has confirmed the young ages of the clusters. We noticed that the clusters with ages between 1 and 10 Myr show a flux excess in H (NIC3/F160W) and/or I (WFPC2/F814W) bands with respect to the evolutionary models. Once more Haro 11 represents a challenge to our understanding.

We present preliminary results from a study of the compact star cluster populations in three local luminous blue compact galaxies: ESO 185-IG 013, ESO 350-IG 038 (a.k.a. Haro 11), and MRK 930. These systems show peculiar morphologies and the presence of hundreds of SCs that have been produced by the past, recent, and/or current starburst phases. We use a complete set of HST images ranging from the UV to IR for each galaxy. Deep images in V (WFPC2/f606w) and I (WFPC2/f814w) are used to capture most of the star cluster candidates up to the old ones (fainter) which have had, in the past, less possibility to be detected. The other bands are used in the SED fitting technique for constraining ages and masses. Our goals are to investigate the evolution of these three blue compact galaxies and the star cluster formation impact on their star formation history.

In a programme of observations of local luminous blue compact galaxies (BCGs), we are investigating kinematics by using tracers of both stars and ionized gas. Here we summarise our program and present new data on the local Lyman break galaxy analogue Haro 11. From spatially-resolved spectroscopy around the near-infrared Ca II triplet, we find that its stars and ionized gas have similar velocity fields. Our programme so far indicates however that emission line velocities can differ locally by a few tens of km/s from the Ca II values. Comparing our data to simple stellar population models, we assess which stellar population the Ca II triplet traces and its potential beyond the local universe.

(abridged) We present near-infrared (NIR) adaptive optics imaging obtained with VLT/NACO and optical spectroscopy from the Southern African Large Telescope (SALT) of a luminous infrared galaxy (LIRG) IRAS 19115-2124. These data are combined with archival HST imaging and Spitzer imaging and spectroscopy, allowing us to study this disturbed interacting/merging galaxy, dubbed the Bird, in extraordinary detail. In particular, the data reveal a triple system where the LIRG phenomenon is dominated by the smallest of the components. One nucleus is a regular barred spiral with significant rotation, while another is highly disturbed with a surface brightness distribution intermediate to that of disk and bulge systems, and hints of remaining arm/bar structure. We derive dynamical masses in the range 3-7x10^10 M_solar for both. The third component appears to be a 1-2x10^10 M_solar mass irregular galaxy. The total system exhibits HII galaxy-like optical line ratios and strengths, and no evidence for AGN activity is found from optical or mid-infrared data. The star formation rate is estimated to be 190 M_solar/yr. We search for SNe, super star clusters, and detect 100-300 km/s outflowing gas from the Bird. Overall, the Bird shows kinematic, dynamical, and emission line properties typical for cool ultra luminous IR galaxies. However, the interesting features setting it apart for future studies are its triple merger nature, and the location of its star formation peak - the strongest star formation does not come from the two major K-band nuclei, but from the third irregular component. Aided by simulations, we discuss scenarios where the irregular component is on its first high-speed encounter with the more massive components.

(Abridged) Lyman-alpha (Lya) is a dominant probe of the galaxy population at high-z. However, interpretation of data drawn from Lya alone hinges on the Lya escape fraction which, due to the complex radiative transport, may vary greatly. Here we map the Lya emission from local starburst Haro 11, a Lya emitter and the only known candidate for low-z Lyman continuum emission (LyC). To aid in the interpretation we perform a detailed multi-wavelength analysis and model the stellar population, dust distribution, ionising photon budget, and star-cluster population. We use archival X-ray observations to further constrain properties of the starburst and estimate the HI column density. The Lya morphology is found to be strongly decoupled from stellar and nebular (H-alpha) morphologies. General surface photometry finds only very slight correlation between Lya and H-halpha, E(B-V), and stellar age. Only around the central Lya-bright cluster do we find the Lya/Ha ratio at values predicted by recombination theory. The total Lya escape fraction is found to be just 3%. We compute that ~90% of the Lya photons that escape do so after undergoing multiple resonance scattering events, masking their point of origin. This leads to a largely symmetric distribution and, by increasing the distance that photons must travel to escape, decreases the escape probability significantly. While dust must ultimately be responsible for the destruction of Lya, it plays little role in governing the observed morphology, which is regulated more by ISM kinematics and geometry. We find tentative evidence for local Lya equivalent width in the immediate vicinity of star-clusters being a function of cluster age, consistent with hydrodynamic studies. We estimate the ionising photon production and further constrain the escape fraction at 900 AA to <~9% .