Search

The Lyman-$\alpha$ (Lya) line of hydrogen is a well-known tracer of galaxies at high-z. However, the connection between Lya observables and galaxy properties has not fully been established, limiting the use of the line to probe the physics of galaxies. Here, we derive global neutral hydrogen gas (HI) properties of nearby Lya-emitting galaxies to assess the impact of HI on the Lya output of galaxies. We observed 21cm line emission using the VLA in D-array configuration (~55" resolution, ~38 kpc) for 37 star-forming galaxies with available Lya imaging from the Lyman Alpha Reference Samples (LARS and eLARS). We detect 21cm emission for 33/37 galaxies observed. We find no significant correlation of global HI properties with Lya luminosity, escape fraction or equivalent width derived with HST photometry. Additionally, both Lya-emitters and weak or non-emitters are distributed evenly along the HI parameter space of optically-selected z=0 galaxies. Around 74% of the sample is undergoing galaxy interaction, this fraction is higher for Lya-emitters (83% for galaxies with EW$\geq$20\r{A}) than for non or weak emitters (70%). Nevertheless, galaxies identified as interacting have Lya and HI properties statistically consistent with those of non-interacting galaxies. Our results show that global HI properties (on scales > 30kpc) have little direct impact on the Lya output from galaxies. Instead, HI likely regulates Lya emission on small scales: statistical comparisons of Lya and high angular resolution 21cm observations are required to fully assess the role of HI in Lya radiative transfer. While our study indicates that galaxy mergers could play a role in the emission of Lya photons in the local universe, especially for galaxies with high HI fractions, the line-of-sight through which a system is observed ultimately determines Lya observables., Comment: 15+35 pages, 5+33 figures, Accepted for publication in A&A

Gas kinematics affect the radiative transfer and escape of hydrogen Lyman-$\alpha$ (Ly$\alpha$) emission from galaxies. We investigate this interplay empirically by relating the ionised gas kinematics of 42 galaxies in the extended Ly$\alpha$ Reference Sample with their Ly$\alpha$ escape fractions, $f_\rm{esc}$, Ly$\alpha$ equivalent widths, $\rm{EW}_\rm{Ly\alpha}$, and Ly$\alpha$ luminosities, $L_\rm{Ly\alpha}$. To this aim we use PMAS integral-field spectroscopic observations of the Balmer-$\alpha$ line. We calculate shearing velocities, $v_\rm{shear}$, and intrinsic velocity dispersions, $\sigma_0^\rm{obs}$ (empirically corrected for beam-smearing effects), as global kinematical measures for each galaxy. The sample is characterised by highly turbulent motions and more than half of the sample shows dispersion dominated kinematics. We uncover clear trends between Ly$\alpha$ observables and global kinematical statistics. We discuss statistically the importance of $v_\rm{shear}$, $\sigma_0^\rm{obs}$, and $v_\rm{shear}/\sigma_0^\rm{obs}$ for regulating the Ly$\alpha$ observables in comparison to other galaxy parameters. It emerges that $\sigma_0^\rm{obs}$ is the dominating parameter for $\rm{EW}_\rm{Ly\alpha}$ and that is as important as nebular extinction, gas covering fraction, and ionising photon production efficiency in regulating $f_\rm{esc}$. A simple scenario where the starburst age is simultaneously regulating turbulence, $\rm{EW}_\rm{Ly\alpha}$, and $f_\rm{esc}$ does not find support by our observations. However, we show that the small scale distribution of dust appears to be influenced by turbulence in some galaxies. In support of our observational result we discuss how turbulence is theoretically expected to play a significant role in modulating $f_\rm{esc}$. (abridged), Comment: Accepted version for publication in A&A

Feedback driven massive outflows play a crucial role in galaxy evolution by regulating star formation and influencing the dynamics of surrounding media. Extracting outflow properties from spectral lines is a notoriously difficult process for a number of reasons, including the possibility that a substantial fraction of the outflow is carried by dense gas in a very narrow range in velocity. This gas can hide in spectra with insufficient resolution. Empirically motivated analysis based on the Apparent Optical Depth method, commonly used in the literature, neglects the contribution of this gas, and may therefore underestimate the true gas column density. More complex semi-analytical line transfer (e.g., SALT) models, on the other hand, allow for the presence of this gas by modeling the radial density and velocity of the outflows as power laws. Here we compare the two approaches to quantify the uncertainties in the inferences of outflow properties based on 1-D "down-the-barrel" using the UV spectra of the CLASSY galaxy sample. We find that empirical modeling may significantly underestimate the column densities relative to SALT analysis, particularly in the optically thick regime. We use simulations to show that the main reason for this discrepancy is the presence of large amount of dense material at low velocities, which can be hidden by the finite spectral resolution of the data. The SALT models in turn could over-estimate the column densities if the assumed power laws of the density profiles strong are not a property of actual outflows., Comment: Accepted to the Astrophysical Journal

"Godzilla" is a peculiar object within the gravitationally lensed Sunburst Arc at $z=2.37$. Despite being very bright, it appears in only one of the twelve lensed images of the source galaxy, and shows exotic spectroscopic properties not found in any other clumps. We use JWST's unique combination of spatial resolution and spectroscopic sensitivity to provide a unified, coherent explanation of the physical nature of Godzilla. We measure fluxes and kinematic properties of rest-optical emission lines in Godzilla and surrounding regions. Using standard line ratio-based diagnostic methods in combination with NIRCam imaging and ground based rest-UV spectra, we characterize Godzilla and its surroundings. Among around 60 detected lines, we find a cascade of strong O I lines pumped by intense Ly$\beta$ emission, as well as Ly$\alpha$-pumped rest-optical Fe II lines, reminiscent of the Weigelt blobs in the local LBV star Eta Carinae. Spectra and images of Godzilla and two faint adjacent images, and the detection of a low-surface brightness foreground galaxy in the NIRCam data, support the interpretation that Godzilla is an extremely magnified object due to the alignment with lensing caustics. We find that Godzilla is part of a previously identified clump, comprising $\sim10 - 25$ % of it, with magnifications in the range of $\approx600-25,000$ depending on the models and images in comparison. The unique O I source in Godzilla is well explained by a non-erupting LBV accompanied by a hotter companion and/or gas condensations exposed to more intense radiation compared to the Weigelt blobs. If Godzilla is confirmed to contain an LBV star, it expands the distance to the furthest known LBV from a dozen Mpc to several Gpc., Comment: 22 pages, 19 figures. Submitted to A&A

In an ongoing search for low-mass extreme emission line galaxies, we identified a galaxy with a Ha/Hb Balmer line ratio of 2.620 +- 0.078. Ha/Hb Balmer ratios lower than the dust-free Case~B value appear relatively frequently in extreme emission line galaxies. These low values suggest that the Case~B assumption may not be valid in these objects. After ruling out the possibility that the low Ha/Hb ratio is due to systematic errors introduced by observational effects, we use constraints from the total Hb luminosity, the [OIII]/[OII] line ratio and the Balmer line equivalent widths, to suggest that the gas is optically thick to both Ha and Lya photons, and the geometry and orientation of the scattering gas causes Ha photons to be preferentially removed from the line of sight with respect to higher order Balmer series photons. Finally, we use data from the SDSS survey to show that Balmer self-absorption may be more important than previously assumed in high excitation emission line galaxies, where Lya pumping of the hydrogen excited state can be effective. If not recognized, Balmer self-absorption could lead to inaccurate estimates of galaxy physical properties. As an example, the effect of dust extinction could be over-estimated, for spherically symmetric scattering medium, or under-estimated, for a not spherically-symmetric distribution.

We report the detection of a population of Wolf-Rayet (WR) stars in the Sunburst Arc, a strongly gravitationally lensed galaxy at redshift $z=2.37$. As the brightest known lensed galaxy, the Sunburst Arc has become an important cosmic laboratory for studying star and cluster formation, Lyman $\alpha$ radiative transfer, and Lyman Continuum (LyC) escape. Here, we present the first results of JWST/NIRSpec IFU observations of the Sunburst Arc, focusing on a stacked spectrum of the 12-fold imaged LyC-emitting (Sunburst LCE) cluster. In agreement with previous studies, we find that the cluster is massive and compact, with $M_{\text{dyn}} = (9\pm1) \times 10^{6} M_{\odot}$, Our age estimate of 4.2--4.5 Myr is much larger than the crossing time of $t_{\text{cross}} = 183 \pm 9 $ kyr, indicating that the cluster is dynamically evolved and consistent with being gravitationally bound. We find a significant nitrogen enhancement of the low ionization state ISM, with $\log(N/O) = -0.74 \pm 0.09$, which is $\approx 0.8$ dex above typical values for H II regions of similar metallicity in the local Universe. We find broad stellar emission complexes around He II$\lambda 4686$ and C IV$\lambda 5808$ with associated nitrogen emission -- this is the first time WR signatures have been directly observed at redshifts above $\sim 0.5$. The strength of the WR signatures cannot be reproduced by stellar population models that only include single-star evolution. While models with binary evolution better match the WR features, they still struggle to reproduce the nitrogen-enhanced WR features. JWST reveals the Sunburst LCE to be a highly ionized, proto-globular cluster with low oxygen abundance and extreme nitrogen enhancement that hosts a population of Wolf-Rayet stars, and possibly Very Massive stars (VMSs), which are rapidly enriching the surrounding medium., Comment: 8 pages, 4 figures, 3 tables. Accepted for publication in A&A

One of the surprising early findings with JWST has been the discovery of a strong "roll-over" or a softening of the absorption edge of Ly$\alpha$ in a large number of galaxies at ($z\gtrsim 6$), in addition to systematic offsets from photometric redshift estimates and fundamental galaxy scaling relations. This has been interpreted as damped Ly$\alpha$ absorption (DLA) wings from high column densities of neutral atomic hydrogen (HI), signifying major gas accretion events in the formation of these galaxies. To explore this new phenomenon systematically, we assemble the JWST/NIRSpec PRImordial gas Mass AssembLy (PRIMAL) legacy survey of 494 galaxies at $z=5.5-13.4$. We characterize this benchmark sample in full and spectroscopically derive the galaxy redshifts, metallicities, star-formation rates, and ultraviolet slopes. We define a new diagnostic, the Ly$\alpha$ damping parameter $D_{\rm Ly\alpha}$ to measure and quantify the Ly$\alpha$ emission strength, HI fraction in the IGM, or local HI column density for each source. The JWST-PRIMAL survey is based on the spectroscopic DAWN JWST Archive (DJA-Spec). All the software, reduced spectra, and spectroscopically derived quantities and catalogs are made publicly available in dedicated repositories. The fraction of strong galaxy DLAs are found to be in the range $65-95\%$ at $z>5.5$. The fraction of strong Ly$\alpha$ emitters (LAEs) is found to increase with decreasing redshift, in qualitative agreement with previous observational results, and are predominantly associated with low-metallicity and UV faint galaxies. By contrast, strong DLAs are observed in galaxies with a variety of intrinsic physical properties. Our results indicate that strong DLAs likely reflect a particular early assembly phase of reionization-era galaxies, at which point they are largely dominated by pristine HI gas accretion. [abridged], Comment: 20 pages, 18 figures, 3 tables. Submitted to A&A. Comments welcome! All data and catalogs are available through the DAWN JWST Archive (DJA): https://dawn-cph.github.io/dja/ and https://github.com/keheintz/jwst-primal

The radio-quiet quasar SDSS J1240+1455 lies at a redshift of z=3.11, is surrounded by a Ly-alpha blob (LAB), and is absorbed by a proximate damped Ly-alpha system. In order to better define the morphology of the blob and determine its emission mechanism, we gathered deep narrow-band images isolating the Ly-alpha line of this object in linearly polarized light. We provide a deep intensity image of the blob, showing a filamentary structure extending up to 16'' (or ~122 physical kpc) in diameter. No significant polarization signal could be extracted from the data, but 95% probability upper limits were defined through simulations. They vary between ~3% in the central 0.75'' disk (after subtraction of the unpolarized quasar continuum) and ~10% in the 3.8-5.5'' annulus. The low polarization suggests that the Ly-alpha photons are emitted mostly in situ, by recombination and de-excitation in a gas largely ionized by the quasar ultraviolet light, rather than by a central source and scattered subsequently by neutral hydrogen gas. This blob shows no detectable polarization signal, contrary to LAB1, a brighter and more extended blob that is not related to the nearby active galactic nucleus (AGN) in any obvious way, and where a significant polarization signal of about 18% was detected., Comment: 15 pages, 8 figures; accepted for publication in A&A

The physical processes by which gas is accreted onto galaxies, transformed into stars and then expelled from galaxies are of paramount importance to galaxy evolution studies. Observationally constraining each of these baryonic components in the same systems however, is challenging. Furthermore, simulations indicate that the stellar mass of galaxies is a key factor influencing CGM properties. Indeed, absorption lines detected against background quasars offer the most compelling way to study the cold gas in the circumgalactic medium (CGM). The MUSE-ALMA Haloes survey is composed of quasar fields covered with VLT/MUSE observations, comprising 32 \ion{H}{i} absorbers at 0.2 $<$ $z$ $<$ 1.4 and 79 associated galaxies, with available or upcoming molecular gas measurements from ALMA. We use a dedicated 40-orbit HST UVIS and IR WFC3 broad-band imaging campaign to characterise the stellar content of these galaxies. By fitting their spectral energy distribution, we establish they probe a wide range of stellar masses: 8.1 $<$ log($M_*$/M$_{\odot}$) $<$ 12.4. Given their star-formation rates, most of these objects lie on the main sequence of galaxies. We also confirm a previously reported anti-correlation between the stellar masses and CGM hydrogen column density N(\ion{H}{i}), indicating an evolutionary trend where higher mass galaxies are less likely to host large amounts of \ion{H}{i} gas in their immediate vicinity up to 120 kpc. Together with other studies from the MUSE-ALMA Haloes survey, these data provide stellar masses of absorber hosts, a key component of galaxy formation and evolution, and observational constraints on the relation between galaxies and their surrounding medium., Comment: accepted in MNRAS

We report observational evidence of highly turbulent ionized gas kinematics in a sample of 20 Lyman continuum (LyC) emitters (LCEs) at low redshift ($z\sim 0.3$). Detailed Gaussian modeling of optical emission line profiles in high-dispersion spectra consistently shows that both bright recombination and collisionally excited lines can be fitted as one or two narrow components with intrinsic velocity dispersion of $\sigma$ $\sim$ 40-100 km s$^{-1}$, in addition to a broader component with $\sigma \sim$ 100-300 km s$^{-1}$, which contributes up to $\sim$40% of the total flux and is preferentially blueshifted from the systemic velocity. We interpret the narrow emission as highly ionized gas close to the young massive star clusters and the broader emission as a signpost of unresolved ionized outflows, resulting from massive stars and supernova feedback. We find a significant correlation between the width of the broad emission and the LyC escape fraction, with strong LCEs exhibiting more complex and broader line profiles than galaxies with weaker or undetected LyC emission. We provide new observational evidence supporting predictions from models and simulations; our findings suggest that gas turbulence and outflows resulting from strong radiative and mechanical feedback play a key role in clearing channels through which LyC photons escape from galaxies. We propose that the detection of blueshifted broad emission in the nebular lines of compact extreme emission-line galaxies can provide a new indirect diagnostic of Lyman photon escape, which could be useful to identify potential LyC leakers in the epoch of reionization with the JWST., Comment: 11 Pages, 7 Figures. A&A in Press

Green Pea galaxies are remarkable for their intense star formation and serve as a window into the early universe. In our study, we used integral field spectroscopy to examine 24 of these galaxies in the optical spectrum. We focused on the interaction between their ionized interstellar medium and the star formation processes within them. Our research generated spatial maps of emission lines and other properties like ionization structures and chemical conditions. These maps showed that areas with higher levels of excitation are usually located where starbursts are occurring. Continuum maps displayed more intricate structures than emission line maps and hinted at low brightness ionized gas in the galaxies' outer regions. We also analyzed integrated spectra from selected areas within these galaxies to derive physical properties like electron densities and temperatures. In some galaxies, we were able to determine metallicity levels. Our observations revealed the presence of high-ionizing lines in three galaxies, two of which had extremely high rates of star formation. Our findings provide valuable insights into the properties and star-forming processes in Green Pea galaxies, contributing to our broader understanding of galactic evolution in the early universe.

Context: Outflows from low-mass star-forming galaxies are a fundamental ingredient for models of galaxy evolution and cosmology. Aims: The onset of kpc-scale ionised filaments in the halo of the metal-poor compact dwarf SBS 0335-052E was previously not linked to an outflow. We here we investigate whether these filaments provide evidence for an outflow. Methods: We obtained new VLT/MUSE WFM and deep NRAO/VLA B-configuration 21cm data of the galaxy. The MUSE data provide morphology, kinematics, and emission line ratios H$\beta$/H$\alpha$ and [\ion{O}{iii}]$\lambda5007$/H$\alpha$ of the low surface-brightness filaments, while the VLA data deliver morphology and kinematics of the neutral gas in and around the system. Both datasets are used in concert for comparisons between the ionised and the neutral phase. Results: We report the prolongation of a lacy filamentary ionised structure up to a projected distance of 16 kpc at $\mathrm{SB}_\mathrm{H\alpha} = 1.5\times10^{-18}$erg s$^{-1}$ cm$^{-2}$arcsec$^{-2}$. The filaments exhibit unusual low H$\alpha$/H$\beta \approx 2.4$ and low [\ion{O}{iii}]/H$\alpha \sim 0.4 - 0.6$ typical of diffuse ionised gas. They are spectrally narrow ($\sim 20$ km s$^{-1}$) and exhibit no velocity sub-structure. The filaments extend outwards of the elongated \ion{H}{I} halo. On small scales the $N_\mathrm{HI}$ peak is offset from the main star-forming sites. Morphology and kinematics of \ion{H}{I} and \ion{H}{II} reveal how star-formation driven feedback interacts differently with the ionised and the neutral phase. Conclusions: We reason that the filaments are a large scale manifestation of star-formation driven feedback, namely limb-brightened edges of a giant outflow cone that protrudes through the halo of this gas-rich system. A simple toy model of such a conical-structure is found to be commensurable with the observations., Comment: Accepted version in A&A after language editing. 22 pages, 24 figures

The CLusters in the Uv as EngineS (CLUES) survey is a Cosmic Origins Spectrograph (COS) campaign aimed at acquiring the 1130 to 1770 {\AA}, restframe spectroscopy of very young (<20 Myr) and massive (>10^4 solar masses) star clusters in galaxies that are part of the Hubble treasury program Legacy ExtraGalactic Uv Survey (LEGUS). In this first paper of a series, we describe the CLUES sample consisting of 20 young star clusters and report their physical properties as derived by both multi-wavelength photometry and far-UV (FUV) spectroscopy with Hubble Space Telescope (HST). Thanks to the synergy of the two different datasets we build a coherent picture of the diverse stellar populations found in each region (with sizes of 40 to 160 pc). We associate the FUV-brightest stellar population to the central targeted star cluster and the other modeled population to the diffuse stars that are included in the COS aperture. We observe better agreement between photometric and spectroscopic ages for star clusters younger than 5 Myr. For clusters older than 5 Myr, photometry and spectroscopy measurements deviate, with the latter producing older ages, due to the degeneracy of photometric models. FUV spectroscopy enables us to better constrain the stellar metallicities, a parameter that optical colors are insensitive to. Finally, the derived E(B-V) are quite similar, with a tendency for FUV spectroscopy to favor solutions with higher extinctions. The recovered masses are in agreement within a factor of 2 for all the clusters., Comment: Accepted for publication on the Astronomical Journal on 16th September 2022. 26 pages, 12 figures

Most of the hydrogen in the intergalactic medium (IGM) was rapidly ionized at high-redshifts. While observations have established that reionization occurred, observational constraints on the emissivity of ionizing photons at high-redshift remains elusive. Here, we present a new analysis of the Low-redshift Lyman Continuum Survey (LzLCS) and archival observations, a combined sample of 89 star-forming galaxies at z~0.3 with Hubble Space Telescope observations of their ionizing continua (or Lyman Continuum, LyC). We find a strong (6$\sigma$ significant) inverse correlation between the continuum slope at 1550\r{A} (defined as F$_\lambda\propto\lambda^{\beta}$) and both the LyC escape fraction (f$_{esc}$) and f$_{esc}$ times the ionizing photon production efficiency ($\xi_{ ion}$). On average, galaxies with redder continuum slopes have smaller f$_{esc}$ than galaxies with bluer slopes due to higher dust attenuation. More than 5% (20%) of the LyC emission escapes galaxies with $\beta$<-2.1 (-2.6). We find strong correlations between $\beta$ and the gas-phase ionization ([OIII]/[OII] flux ratio; at 7.5$\sigma$ significance), galaxy stellar mass (at 5.9$\sigma$), the gas-phase metallicity (at 4.6$\sigma$), and the observed FUV absolute magnitude (M$_{UV}$ at 3.4$\sigma$). Using previous observations of $\beta$ at high-redshift, we estimate the evolution of f$_{esc}$ with both $z$ and M$_{UV}$. The LzLCS suggest that fainter and lower mass galaxies dominate the ionizing photon budget at higher redshift, possibly due to their rapidly evolving metal and dust content. Finally, we use our correlation between $\beta$ and f$_{ esc}\times\xi_{ion}$ to predict the ionizing emissivity of galaxies during the epoch of reionization. Our estimated emissivities match IGM observations, and suggest that star-forming galaxies emit sufficient LyC photons into the IGM to exceed recombinations near redshifts of 7-8., Comment: 14 pages plus appendix. Submitted to MNRAS. Comments encouraged

The Ly$\alpha$ emission line is one of the main observables of galaxies at high redshift, but its output depends strongly on the neutral gas distribution and kinematics around the star-forming regions where UV photons are produced. We present observations of Ly$\alpha$ and 21-cm HI emission at comparable scales with the goal to qualitatively investigate how the neutral interstellar medium (ISM) properties impact Ly$\alpha$ transfer in galaxies. We have observed 21-cm HI at the highest angular resolution possible (~ 3" beam) with the VLA in two local galaxies from the Lyman Alpha Reference Sample. We contrast this data with HST Ly$\alpha$ imaging and spectroscopy, and MUSE and PMAS ionized gas observations. In LARS08, high intensity Ly$\alpha$ emission is co-spatial with high column density HI where dust content is the lowest. The Ly$\alpha$ line is strongly redshifted, consistent with velocity redistribution which allows Ly$\alpha$ escape from high column density neutral medium with low dust content. In eLARS01, high intensity Ly$\alpha$ emission is located in regions of low column density HI, below the HI data sensitivity limit ($<2\times10^{20}\,$cm$^{-2}$). The perturbed ISM distribution with low column density gas in front of the Ly$\alpha$ emission region plays an important role in the escape. In both galaxies, the faint Ly$\alpha$ emission ($\sim 1\times10^{-16}$erg.s$^{-1}$cm$^{-2}$arcsec$^{-2}$) traces intermediate H$\alpha$ emission regions where HI is found, regardless of the dust content. Dust seems to modulate, but not prevent, the formation of a faint Ly$\alpha$ halo. This study suggests the existence of scaling relations between dust, H$\alpha$, HI, and Ly$\alpha$ emission in galaxies., Comment: 28 pages, 14 figures, accepted for publication in ApJ

The properties that govern the production and escape of hydrogen ionizing photons (Lyman continuum, LyC; with energies >13.6 eV) in star-forming galaxies are still poorly understood, but they are key to identifying and characterizing the sources that reionized the Universe. Here we empirically explore the relationship between the hardness of ionizing radiation and the LyC leakage in a large sample of low-$z$ star-forming galaxies from the recent Hubble Space Telescope Low-$z$ Lyman Continuum Survey. Using Sloan Digital Sky Survey stacks and deep XShooter observations, we investigate the hardness of the ionizing spectra ($Q_{\rm He^+}/Q_{\rm H}$) between 54.4 eV (He$^{+}$) and 13.6 eV (H) from the optical recombination lines HeII 4686A and H$\beta$ 4861A for galaxies with LyC escape fractions spanning a wide range, $f_{\rm esc} \rm (LyC) \simeq 0 - 90\%$. We find that the observed intensity of HeII/H$\beta$ is primarily driven by variations in the metallicity, but is not correlated with LyC leakage. Both very strong ($ \simeq 0.5$) and nonleakers ($ < f_{\rm esc} \rm (LyC) > \simeq 0$) present similar observed intensities of HeII/H$\beta$ at comparable metallicity, between $\simeq 0.01$ and $\simeq 0.02$ for $12 + \log({\rm O/H}) > 8.0$ and $<8.0$, respectively. Our results demonstrate that $Q_{\rm He^+}/Q_{\rm H}$ does not correlate with $f_{\rm esc} \rm (LyC)$, which implies that strong LyC emitters do not show harder ionizing spectra than nonleakers at similar metallicity., Comment: 6 pages, 2 figures, 1 tables. Accepted for publication in A&A Letters. v2: updated to match published version, including new title

Star formation is a clustered process that regulates the structure and evolution of galaxies. We investigate this process in the dwarf galaxy Haro 11, forming stars in three knots (A, B, C). The exquisite resolution of HST imaging allows us to resolve the starburst into tens of bright star clusters. We derive masses between $10^5$ and $10^7\,\rm M_{\odot}$ and ages younger than 20 Myr, using photometric modeling. We observe that the clustered star formation has propagated from knot C (the oldest) through knot A (in between) towards knot B (the youngest). We use aperture-matched ultraviolet and optical spectroscopy (HST + MUSE) to independently study the stellar populations of Haro 11 and determine the physical properties of the stellar populations and their feedback in 1 kpc diameter regions. We discuss these results in light of the properties of the ionised gas within the knots. We interpret the broad blue-shifted components of the optical emission lines as outflowing gas ($v_{max} \sim 400$ km/s). The strongest outflow is detected in knot A with a mass-rate of $\dot{M}_{out}\sim 10\,\rm M_{\odot}/yr$, ten times higher than the star-formation in the same region. Knot B hosts a young and not fully developed outflow, whereas knot C has likely been already evacuated. Because Haro 11 has properties similar to high-redshift unresolved galaxies, our work can additionally aid the understanding of star formation at high redshift, a window that will be opened by upcoming facilities., Comment: Accepted for publication in MNRAS on 22nd Dec 2022: 18 pages, 9 figures