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.

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 Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

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.

Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decreasing metal content mimicking the environment at the earliest epochs. Ultimately, a physical model that could be extrapolated to zero metallicity would enable tackling long-standing questions such as "What did the First, very massive stars of the Universe look like?" or "What was their role in the re-ionization of the Universe?". Yet, most our knowledge of metal-poor massive stars is drawn from one single point in metallicity. Massive stars in the Small Magellanic Cloud (SMC, $\sim 1/5 Z_{\odot}$) currently serve as templates for low-metallicity objects in the early Universe, even though significant differences with respect to massive stars with poorer metal content have been reported. This White Paper summarizes the current knowledge on extremely (sub-SMC) metal poor massive stars, highlighting the most outstanding open questions and the need to supersede the SMC as standard. A new paradigm can be built from nearby extremely metal-poor galaxies that make a new metallicity ladder, but massive stars in these galaxies are out of reach to current observational facilities. Such task would require an L-size mission, consisting of a 10m-class space telescope operating in the optical and the ultraviolet ranges. Alternatively, we propose that ESA unites efforts with NASA to make the LUVOIR mission concept a reality, thus continuing the successful partnership that made Hubble Space Telescope one of the greatest observatories of all time.

The prioritization and improvement of ethics, planetary protection, and safety standards in the astro-sciences is the most critical priority as our scientific and exploratory capabilities progress, both within government agencies and the private sector. These priorities lie in the belief that every single science mission - crewed or non-crewed, ground-based or not - should heed strict ethical and safety standards starting at the very beginning of a mission. Given the inevitability of the private sector in influencing future crewed missions both in and beyond low-Earth orbit, it is essential to the science community to agree on universal standards of safety, mission assurance, planetary protection, and especially anti-colonization. These issues will impact all areas of space science. Examples that are particularly relevant to the Astro2020 Decadal Survey include but are not limited to: light pollution from satellites, the voices and rights of Native people when constructing telescopes on their lands, and the need to be cognizant of contamination when searching for and exploring habitable environments beyond Earth. Ultimately, moving international space law and domestic space policy from a reactive nature to a proactive one will ensure the future of space exploration is one that is safe, transparent, and anti-imperialist.

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.

How do stars manage to form within low-density, HI-dominated gas? Such environments provide a laboratory for studying star formation with physical conditions distinct from starbursts and the metal-rich disks of spiral galaxies where most effort has been invested. Here we outline fundamental open questions about the nature of star formation at low-density. We describe the wide-field, high-resolution UV-optical-IR-radio observations of stars, star clusters and gas clouds in nearby galaxies needed in the 2020's to provide definitive answers, essential for development of a complete theory of star formation.

Multiple generations of massive stars have lived and died during Cosmic History, invigorating host galaxies with ionizing photons, kinetic energy, fresh material and stellar-size black holes. At present, massive stars in the Small Magellanic Cloud (SMC) serve as templates for low-metallicity objects in the early Universe. However, recent results have highlighted important differences in the evolution, death and feedback of massive stars with poorer metal content that better matches the extremely low metallicity of previous Cosmic epochs. This paper proposes to supersede the SMC standard with a new metallicity ladder built from very metal-poor galaxies, and provides a brief overview of the technological facilities needed to this aim.

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.

Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale.

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.

Measurements of the galaxy UV luminosity function at z>6 suggest that young stars hosted in low-mass star-forming galaxies produced the bulk of hydrogen-ionizing photons necessary to reionize the intergalactic medium (IGM) by redshift z~6. Whether star-forming galaxies dominated cosmic reionization, however, also depends on their stellar populations and interstellar medium properties, which set, among other things, the production rate of H-ionizing photons, $\xi_\text{ion}^\star$, and the fraction of these escaping into the IGM. Given the difficulty of constraining with existing observatories the physical properties of z>6 galaxies, in this work we focus on a sample of ten nearby objects showing UV spectral features comparable to those observed at z>6. We use the new-generation Beagle tool to model the UV-to-optical photometry and UV/optical emission lines of these Local 'analogues' of high-redshift galaxies, finding that our relatively simple, yet fully self-consistent, physical model can successfully reproduce the different observables considered. Our galaxies span a broad range of metallicities and are characterised by high ionization parameters, low dust attenuation, and very young stellar populations. Through our analysis, we derive a novel diagnostic of the production rate of H-ionizing photons per unit UV luminosity, $\xi_\text{ion}^\star$, based on the equivalent width of the bright $[\text{OIII}] \lambda 4959,5007$ doublet, which does not require measurements of H-recombination lines. This new diagnostic can be used to estimate $\xi_\text{ion}^\star$ from future direct measurements of the $[\text{OIII}] \lambda 4959,5007$ line using JWST/NIRSpec (out to z~9.5), and by exploiting the contamination by $\text{H}\beta + [\text{OIII}] \lambda 4959,5007$ of photometric observations of distant galaxies, for instance from existing Spitzer/IRAC data and from future ones with JWST/NIRCam.

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.

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.

The cosmological origin of carbon, the fourth most abundant element in the Universe, is not well known and matter of heavy debate. We investigate the behavior of C/O to O/H in order to constrain the production mechanism of carbon. We measured emission-line intensities in a spectral range from 1600 to 10000 \AA\ on Space Telescope Imaging Spectrograph (STIS) long-slit spectra of 18 starburst galaxies in the local Universe. We determined chemical abundances through traditional nebular analysis and we used a Markov Chain Monte Carlo (MCMC) method to determine where our carbon and oxygen abundances lie in the parameter space. We conclude that our C and O abundance measurements are sensible. We analyzed the behavior of our sample in the [C/O] vs. [O/H] diagram with respect to other objects such as DLAs, neutral ISM measurements, and disk and halo stars, finding that each type of object seems to be located in a specific region of the diagram. Our sample shows a steeper C/O vs. O/H slope with respect to other samples, suggesting that massive stars contribute more to the production of C than N at higher metallicities, only for objects where massive stars are numerous; otherwise intermediate-mass stars dominate the C and N production.

Nearby dwarf galaxies provide a unique laboratory in which to test stellar population models below $Z_\odot/2$. Such tests are particularly important for interpreting the surprising high-ionization UV line emission detected at $z>6$ in recent years. We present HST/COS ultraviolet spectra of ten nearby metal-poor star-forming galaxies selected to show He II emission in SDSS optical spectra. The targets span nearly a dex in gas-phase oxygen abundance ($7.8<12+\log\mathrm{O/H}<8.5$) and present uniformly large specific star formation rates (sSFR $\sim 10^2$ $\mathrm{Gyr}^{-1}$). The UV spectra confirm that metal-poor stellar populations can power extreme nebular emission in high-ionization UV lines, reaching C III] equivalent widths comparable to those seen in systems at $z\sim 6-7$. Our data reveal a marked transition in UV spectral properties with decreasing metallicity, with systems below $12+\log\mathrm{O/H}\lesssim 8.0$ ($Z/Z_\odot \lesssim 1/5$) presenting minimal stellar wind features and prominent nebular emission in He II and C IV. This is consistent with nearly an order of magnitude increase in ionizing photon production beyond the $\mathrm{He^+}$-ionizing edge relative to H-ionizing flux as metallicity decreases below a fifth solar, well in excess of standard stellar population synthesis predictions. Our results suggest that often neglected sources of energetic radiation such as stripped binary products and very massive O-stars produce a sharper change in the ionizing spectrum with decreasing metallicity than expected. Consequently, nebular emission in C IV and He II powered by these stars may provide useful metallicity constraints in the reionization era.

We studied Lyman-$\alpha$ (Ly$\alpha$) escape in a statistical sample of 43 Green Peas with HST/COS Ly$\alpha$ spectra. Green Peas are nearby star-forming galaxies with strong [OIII]$\lambda$5007 emission lines. Our sample is four times larger than the previous sample and covers a much more complete range of Green Pea properties. We found that about 2/3 of Green Peas are strong Ly$\alpha$ line emitters with rest-frame Ly$\alpha$ equivalent width $>20$ \AA. The Ly$\alpha$ profiles of Green Peas are diverse. The Ly$\alpha$ escape fraction, defined as the ratio of observed Ly$\alpha$ flux to intrinsic Ly$\alpha$ flux, shows anti-correlations with a few Ly$\alpha$ kinematic features -- both the blue peak and red peak velocities, the peak separations, and FWHM of the red portion of the Ly$\alpha$ profile. Using properties measured from SDSS optical spectra, we found many correlations -- Ly$\alpha$ escape fraction generally increases at lower dust reddening, lower metallicity, lower stellar mass, and higher [OIII]/[OII] ratio. We fit their Ly$\alpha$ profiles with the HI shell radiative transfer model and found Ly$\alpha$ escape fraction anti-correlates with the best-fit $N_{HI}$. Finally, we fit an empirical linear relation to predict Ly$\alpha$ escape fraction from the dust extinction and Ly$\alpha$ red peak velocity. The standard deviation of this relation is about 0.3 dex. This relation can be used to isolate the effect of IGM scatterings from Ly$\alpha$ escape and to probe the IGM optical depth along the line of sight of each $z>7$ Ly$\alpha$ emission line galaxy in the JWST era.

Green Peas are nearby analogs of high-redshift Ly$\alpha$-emitting galaxies. To probe their Ly$\alpha$ escape, we study the spatial profiles of Ly$\alpha$ and UV continuum emission of 24 Green Pea galaxies using the Cosmic Origins Spectrograph (COS) on Hubble Space Telescope (HST). We extract the spatial profiles of Ly$\alpha$ emission from their 2D COS spectra, and of UV continuum from both the 2D spectra and NUV images. The Ly$\alpha$ emission shows more extended spatial profiles than the UV continuum in most Green Peas. The deconvolved Full Width Half Maximum (FWHM) of the Ly$\alpha$ spatial profile is about 2 to 4 times that of the UV continuum in most cases. Since Green Peas are analogs of high-z LAEs, it suggests that most high-z LAEs likely have larger Ly$\alpha$ sizes than UV sizes. We also compare the spatial profiles of Ly$\alpha$ photons at blueshifted and redshifted velocities in eight Green Peas with sufficient data quality, and find the blue wing of the Ly$\alpha$ line has a larger spatial extent than the red wing in four Green Peas with comparatively weak blue Ly$\alpha$ line wings. We show that Green Peas and MUSE $z=3-6$ LAEs have similar Ly$\alpha$ and UV continuum sizes, which probably suggests starbursts in both low-z and high-z LAEs drive similar gas outflows illuminated by Ly$\alpha$ light. Five Lyman continuum (LyC) leakers in this sample have similar Ly$\alpha$ to UV continuum size ratios (~1.4-4.3) to the other Green Peas, indicating their LyC emission escape through ionized holes in the interstellar medium.

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 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 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 present the science cases and technological discussions that came from the workshop entitled "Finding the UV-Visible Path Forward" held at NASA GSFC June 25-26, 2015. The material presented outlines the compelling science that can be enabled by a next generation space-based observatory dedicated for UV-visible science, the technologies that are available to include in that observatory design, and the range of possible alternative launch approaches that could also enable some of the science. The recommendations to the Cosmic Origins Program Analysis Group from the workshop attendees on possible future development directions are outlined.

We study the ionization structure of galactic outflows in 37 nearby, star forming galaxies with the Cosmic Origins Spectrograph on the Hubble Space Telescope. We use the O I, Si II, Si III, and Si IV ultraviolet absorption lines to characterize the different ionization states of outflowing gas. We measure the equivalent widths, line widths, and outflow velocities of the four transitions, and find shallow scaling relations between them and galactic stellar mass and star formation rate. Regardless of the ionization potential, lines of similar strength have similar velocities and line widths, indicating that the four transitions can be modeled as a co-moving phase. The Si equivalent width ratios (e.g. Si IV/Si II) have low dispersion, and little variation with stellar mass; while ratios with O I and Si vary by a factor of 2 for a given stellar mass. Photo-ionization models reproduce these equivalent width ratios, while shock models under predict the relative amount of high ionization gas. The photo-ionization models constrain the ionization parameter (U) between -2.25 < log(U) < -1.5, and require that the outflow metallicities are greater than 0.5 Z$_\odot$. We derive ionization fractions for the transitions, and show that the range of ionization parameters and stellar metallicities leads to a factor of 1.15-10 variation in the ionization fractions. Historically, mass outflow rates are calculated by converting a column density measurement from a single metal ion into a total Hydrogen column density using an ionization fraction, thus mass outflow rates are sensitive to the assumed ionization structure of the outflow.

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.

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.

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.

We report on a sample of 51 nearby, star-forming galaxies observed with the Cosmic Origin Spectrograph on the Hubble Space Telescope. We calculate Si II kinematics and densities arising from warm gas entrained in galactic outflows. We use multi-wavelength ancillary data to estimate stellar masses (M$_\ast$), star-formation rates (SFR), and morphologies. We derive significant correlations between outflow velocity and SFR$^{\sim 0.1}$, M$_\ast^{\sim 0.1}$ and v$_\text{circ}^{\sim 1/2}$. Some mergers drive outflows faster than these relations prescribe, launching the outflow faster than the escape velocity. Calculations of the mass outflow rate reveal strong scaling with SFR$^{\sim 1/2}$ and M$_\ast^{\sim 1/2}$. Additionally, mass-loading efficiency factors (mass outflow rate divided by SFR) scale approximately as M$_\ast^{-1/2}$. Both the outflow velocity and mass-loading scaling suggest that these outflows are powered by supernovae, with only 0.7% of the total supernovae energy converted into the kinetic energy of the warm outflow. Galaxies lose some gas if log(M$_\ast$/M$_\odot$) < $9.5$, while more massive galaxies retain all of their gas, unless they undergo a merger. This threshold for gas loss can explain the observed shape of the mass-metallicity relation.

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.

Deep spectroscopic observations of z~6.5 galaxies have revealed a marked decline with increasing redshift in the detectability of Lyman-alpha emission. While this may offer valuable insight into the end of the reionisation process, it presents a fundamental challenge to the detailed spectroscopic study of the many hundreds of photometrically-selected distant sources now being found via deep HST imaging, and particularly those bright sources viewed through foreground lensing clusters. In this paper we demonstrate the validity of a new way forward via the convincing detection of an alternative diagnostic line, CIII]1909, seen in spectroscopic exposures of two star forming galaxies at z=6.029 and 7.213. The former detection is based on a 3.5 hour X-shooter spectrum of a bright (J=25.2) gravitationally-lensed galaxy behind the cluster Abell 383. The latter detection is based on a 4.2 hour MOSFIRE spectra of one of the most distant spectroscopically confirmed galaxies, GN-108036, with J=25.2. Both targets were chosen for their continuum brightness and previously-known redshift (based on Lyman-alpha), ensuring that any CIII] emission would be located in a favorable portion of the near-infrared sky spectrum. We compare our CIII] and Lyman-alpha equivalent widths in the context of those found at z~2 from earlier work and discuss the motivation for using lines other than Lyman-alpha to study galaxies in the reionisation era.

We present deep spectroscopy of 17 very low mass (M* ~ 2.0x10^6 Msun to 1.4x10^9 Msun) and low luminosity (M_UV ~ -13.7 to -19.9) gravitationally lensed galaxies in the redshift range z~1.5-3.0. Deep rest-frame ultraviolet spectra reveal large equivalent width emission from numerous lines (NIV], OIII], CIV, Si III], CIII]) which are rarely seen in individual spectra of more massive star forming galaxies. CIII] is detected in 16 of 17 low mass star forming systems with rest-frame equivalent widths as large as 13.5 Angstroms. Nebular CIV emission is present in the most extreme CIII] emitters, requiring an ionizing source capable of producing a substantial component of photons with energies in excess of 47.9 eV. Photoionization models support a picture whereby the large equivalent widths are driven by the increased electron temperature and enhanced ionizing output arising from metal poor gas and stars, young stellar populations, and large ionization parameters. The young ages implied by the emission lines and continuum SEDs indicate that the extreme line emitters in our sample are in the midst of a significant upturn in their star formation activity. The low stellar masses, blue UV colors, and large sSFRs of our sample are similar to those of typical z>6 galaxies. Given the strong attenuation of Ly-alpha in z>6 galaxies we suggest that CIII] is likely to provide our best probe of early star forming galaxies with ground-based spectrographs and one of the most efficient means of confirming z>10 galaxies with the James Webb Space Telescope.

Super star cluster A1 in the nearby starburst galaxy NGC 3125 is characterized by broad He\ii \lam1640 emission (full width at half maximum, $FWHM\sim1200$ km s$^{-1}$) of unprecedented strength (equivalent width, $EW=7.1\pm0.4$ \AA). Previous attempts to characterize the massive star content in NGC 3125-A1 were hampered by the low resolution of the UV spectrum and the lack of co-spatial panchromatic data. We obtained far-UV to near-IR spectroscopy of the two principal emitting regions in the galaxy with the Space Telescope Imaging Spectrograph (STIS) and the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope (\hst). We use these data to study three clusters in the galaxy, A1, B1, and B2. We derive cluster ages of 3-4 Myr, intrinsic reddenings of $E(B-V)=0.13$, 0.15, and 0.13, and cluster masses of $1.7\times10^5$, $1.4\times10^5$, and $1.1\times10^5$ M$_\odot$, respectively. A1 and B2 show O\vb \lam1371 absorption from massive stars, which is rarely seen in star-forming galaxies, and have Wolf-Rayet (WR) to O star ratios of $N(WN5-6)/N(O)=0.23$ and 0.10, respectively. The high $N(WN5-6)/N(O)$ ratio of A1 cannot be reproduced by models that use a normal IMF and generic WR star line luminosities. We rule out that the extraordinary He\ii \lam1640 emission and O\vb \lam1371 absorption of A1 are due to an extremely flat upper IMF exponent, and suggest that they originate in the winds of very massive ($>120\,M_\odot$) stars. In order to reproduce the properties of peculiar clusters such as A1, the present grid of stellar evolution tracks implemented in Starburst99 needs to be extended to masses $>120\,M_\odot$.

We obtained medium-resolution ultraviolet (UV) spectra between 1150 and 1450 Angstroms of the four UV-bright, infrared (IR)-luminous starburst galaxies IRAS F08339+6517, NGC 3256, NGC 6090, and NGC 7552 using the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. The selected sightlines towards the starburst nuclei probe the properties of the recently formed massive stars and the physical conditions in the starburst-driven galactic superwinds. Despite being metal-rich and dusty, all four galaxies are strong Lyman-alpha emitters with equivalent widths ranging between 2 and 13 Angstroms. The UV spectra show strong P Cygni-type high-ionization features indicative of stellar winds and blueshifted low-ionization lines formed in the interstellar and circumgalactic medium. We detect outflowing gas with bulk velocities of about 400 km/s and maximum velocities of almost 900 km/s. These are among the highest values found in the local universe and comparable to outflow velocities found in luminous Lyman-break galaxies at intermediate and high redshift. The outflow velocities are unlikely to be high enough to cause escape of material from the galactic gravitational potential. However, the winds are significant for the evolution of the galaxies by transporting heavy elements from the starburst nuclei and enriching the galaxy halos. The derived mass outflow rates of ~100 Msol/yr are comparable to, or even higher than the star-formation rates. The outflows can quench star formation and ultimately regulate the starburst as has been suggested for high-redshift galaxies.

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

The usefulness of H I Lyman-alpha photons for characterizing star formation in the distant universe is limited by our understanding of the astrophysical processes that regulate their escape from galaxies. These processes can only be observed in detail out to a few x100 Mpc. Past nearby (z<0.3) spectroscopic studies are based on small samples and/or kinematically unresolved data. Taking advantage of the high sensitivity of HST's COS, we observed the Lyman-alpha lines of 20 H-alpha-selected galaxies located at <z>=0.03. The galaxies cover a broad range of luminosity, oxygen abundance, and reddening. In this paper, we characterize the observed Lyman-alpha lines and establish correlations with fundamental galaxy properties. We find seven emitters. These host young (\le 10 Myr) stellar populations, have rest-frame equivalent widths in the range 1-12 \AA, and have Lyman-alpha escape fractions within the COS aperture in the range 1-12 %. One emitter has a double-peaked Lyman-alpha with peaks 370 km/s apart and a stronger blue peak. Excluding this object, the emitters have Lyman-alpha and O I \lambda 1302 offsets from H-alpha in agreement with expanding shell models and LBG observations. The absorbers have offsets that are almost consistent with a static medium. We find no one-to-one correspondence between Lyman-alpha emission and age, metallicity, or reddening. Thus, we confirm that Lyman-alpha is enhanced by outflows and is regulated by the dust and H I column density surrounding the hot stars.

We describe observations in the nearby universe (<100 Mpc) with a 10-m or larger space-based telescope having imaging and spectral capabilities in the range 912-9000 \AA that would enable advances in the fields of massive stars, young populations, and star-forming galaxies, that are essential for achieving the Cosmic Origins Program objectives i) how are the chemical elements distributed in galaxies and dispersed in the circumgalactic and intergalactic medium; and ii) when did the first stars in the universe form, and how did they influence their environments. We stress the importance of observing hundreds of massive stars and their descendants individually, which will make it possible to separate the many competing factors that influence the observed properties of these systems (mass, composition, convection, mass-loss, rotation rate, binarity, magnetic fields, and cluster mass).

We analyze archival HST/STIS/FUV-MAMA imaging and spectroscopy of 13 compact star clusters within the circumnuclear starburst region of M83, the closest such example. We compare the observed spectra with semi-empirical models, which are based on an empirical library of Galactic O and B stars observed with IUE, and with theoretical models, which are based on a new theoretical UV library of hot massive stars computed with WM-Basic. The models were generated with Starburst99 for metallicities of Z=0.020 and Z=0.040, and for stellar IMFs with upper mass limits of 10, 30, 50, and 100 M_sol. We estimate the ages and masses of the clusters from the best fit model spectra, and find that the ages derived from the semi-empirical and theoretical models agree within a factor of 1.2 on average. A comparison of the spectroscopic age estimates with values derived from HST/WFC3/UVIS multi-band photometry shows a similar level of agreement for all but one cluster. The clusters have a range of ages from about 3 to 20 Myr, and do not appear to have an age gradient along M83's starburst. Clusters with strong P-Cygni profiles have masses of a few times 10^4 M_sol, seem to have formed stars more massive than 30 M_sol, and are consistent with a Kroupa IMF from 0.1-100 M_sol. Field regions in the starburst lack P-Cygni profiles and are dominated by B stars.

The circumnuclear starburst of M83 (NGC 5236), the nearest such example (4.6 Mpc), constitutes an ideal site for studying the massive star IMF at high metallicity (12+log[O/H]=9.1$\pm$0.2, Bresolin & Kennicutt 2002). We analyzed archival HST/STIS FUV imaging and spectroscopy of 13 circumnuclear star clusters in M83. We compared the observed spectra with two types of single stellar population (SSP) models, semi-empirical models, which are based on an empirical library of Galactic O and B stars observed with IUE (Robert et al. 1993), and theoretical models, which are based on a new theoretical UV library of hot massive stars described in Leitherer et al. (2010) and computed with WM-Basic (Pauldrach et al. 2001). The models were generated with Starburst99 (Leitherer & Chen 2009). We derived the reddenings, the ages, and the masses of the clusters from model fits to the FUV spectroscopy, as well as from optical HST/WFC3 photometry.

In general, HII regions do not show clear signs of self-enrichment in products from massive stars (M > 8 M_sun). In order to explore why, I modeled the contamination with Wolf-Rayet star ejecta of metal-poor (Z=0.001) HII regions, ionised either by a 10^6 M_sun cluster of coeval stars (cluster 1), or a cluster resulting from continuous star formation at a rate of 1 M_sun yr^-1 (cluster 2). The clusters have Z=0.001 and a Salpeter initial mass function (IMF) from 0.1 to 120 M_sun. Independent one dimensional constant density simulations of the emission-line spectra of unenriched HII regions were computed at the discrete ages 1, 2, 3, 4, and 5 Myr, with the photoionisation code CLOUDY, using as input, radiative and mechanical stellar feedbacks predicted by the evolutionary synthesis code STARBURST99. Each HII region was placed at the outer radius of the adiabatically expanding superbubble of Mac Low and McCray (1988). For models with thermal and ionisation balance time-scales of less than 1 Myr, and with oxygen emission-line ratios in agreement with observations, the interior of the superbubble and the HII region were uniformly and instantaneously polluted with stellar ejecta predicted by STARBURST99. I obtained a maximum oxygen abundance enhancement of 0.025 dex, with cluster 1, at 4 Myr. It would be unobservable.