We demonstrate the redshift-evolution of the spectral profile of HI Lyman-alpha (Ly{\alpha}) emission from star-forming galaxies. In this first study we pay special attention to the contribution of blueshifted emission. At redshift z=2.9-6.6, we compile spectra of a sample of 229 Ly{\alpha}-selected galaxies identified with the Multi-Unit Spectroscopic Explorer at the Very Large Telescope, while at low-z (< 0.44) we use a sample of 74 ultraviolet-selected galaxies observed with the Cosmic Origin Spectrograph onboard the Hubble Space Telescope. At low-z, where absorption from the intergalactic medium (IGM) is negligible, we show that the relative contribution of blueshifted emission (blue/red) increases rapidly with increasing equivalent width ($W_{Ly{\alpha}}$). This result does not, however, emerge at z=3-4, and we use bootstrap simulations to demonstrate the blue emission to be suppressed by stochastic IGM absorption. Our main result is that the blue-peak contribution evolves rapidly downwards with increasing redshift: the blue/red flux ratio is $\approx 30$% at $z\approx 0$, but drops to 15% at $z\approx 3$, and to below 3% by $z\approx 6$. Applying further simulations of the IGM absorption to the unabsorbed COS spectrum, we demonstrate that this decrease in the blue-wing contribution can be entirely attributed to the thickening of intervening Ly{\alpha} absorbing systems, with no need for additional HI opacity from local structure, companion galaxies, or cosmic infall. We discuss our results in light of the numerical radiative transfer simulations, the evolving total Ly{\alpha} output of galaxies, and the utility of resolved Ly{\alpha} spectra in the reionization epoch.

We study the global dust and (molecular) gas content in the Lyman Alpha Reference Sample (LARS), i.e. 14 local star-forming galaxies. We characterize their interstellar medium and relate newly derived properties to quantities relevant for Ly$\alpha$ escape. We observed LARS galaxies with Herschel/PACS, SOFIA/FIFI-LS, the IRAM 30m telescope and APEX, targeting far-infrared (FIR) continuum and emission lines of [C II]158$\mu$m, [O I]63$\mu$m, [O III]88$\mu$m and low-J CO lines. Using Bayesian methods we derive dust model parameters and estimate total gas masses for all LARS galaxies, taking into account a metallicity-dependent gas-to-dust ratio. Star formation rates were estimated from FIR, [C II]158$\mu$m and [O I]63$\mu$m luminosities. LARS covers a wide dynamic range in the derived properties, with FIR-based star formation rates from $\sim$0.5-100 $M_{\odot}\ yr^{-1}$, gas fractions between $\sim$15-80% and gas depletion times ranging from a few hundred Myr up to more than 10 Gyr. The distribution of LARS galaxies in the $\Sigma_{gas}$ vs. $\Sigma_{SFR}$ (Kennicutt-Schmidt plane) is thus quite heterogeneous. However, we find that LARS galaxies with the longest gas depletion times, i.e. relatively high gas surface densities ($\Sigma_{gas}$) and low star formation rate densities ($\Sigma_{SFR}$), have by far the highest Ly$\alpha$ escape fraction. A strong $\sim$linear relation is found between Ly$\alpha$ escape fraction and the total gas (HI+H$_2$) depletion time. We argue that the Ly$\alpha$ escape in those galaxies is driven by turbulence in the star-forming gas that shifts the Ly$\alpha$ photons out of resonance close to the places where they originate. We further report on an extreme [C II]158$\mu$m excess in LARS 5, corresponding to $\sim$14$\pm$3% of the FIR luminosity, i.e. the most extreme [C II]-to-FIR ratio observed in a non-AGN galaxy to date.

Understanding the production and escape of Lyman $\alpha$ (Ly$\alpha$) radiation from star-forming galaxies is a long standing problem in astrophysics. The ability to predict the Ly$\alpha$ luminosity of galaxies would open up new ways of exploring the Epoch of Reionization (EoR), and to estimate Ly$\alpha$ emission from galaxies in cosmological simulations where radiative transfer calculations cannot be done. We apply multivariate regression methods to the Lyman Alpha Reference Sample dataset to obtain a relation between the galaxy properties and the emitted Ly$\alpha$. The derived relation predicts the Ly$\alpha$ luminosity of our galaxy sample to good accuracy, regardless of whether we consider only direct observables (root-mean-square (RMS) dispersion around the relation of $\sim 0.19$ dex) or derived physical quantities (RMS $\sim 0.27$ dex). We confirm the predictive ability on a separate sample of compact star-forming galaxies and find that the prediction works well, but that aperture effects on measured Ly$\alpha$ luminosity may be important, depending on the redshift of the galaxy. We apply statistical feature selection techniques to determine an order of importance of the variables in our dataset, enabling future observations to be optimized for predictive ability. When using physical variables, we are able to determine that the most important predictive parameters are, in order, star formation rate, dust extinction, compactness and the gas covering fraction. We discuss the application of our results in terms of studying the EoR and intensity mapping experiments.

Ly$\alpha$ emission is a standard tracer of starburst galaxies at high redshift. However, a number of local Ly$\alpha$ emitters (LAEs) are X-ray sources, suggesting a possible origin of Ly$\alpha$ photons other than young, hot stars, and which may be active at much later ages relative to the parent starburst. Resolved, nearby LAEs offer the opportunity to discriminate between diffuse X-ray emission arising from supernova-heated gas, high-mass X-ray binaries (HMXBs), or low-luminosity active galactic nuclei (LLAGN). We examine archival X-ray imaging from Chandra and XMM-Newton for 11 galaxies with spatially resolved Ly$\alpha$ imaging to determine the luminosity, morphology, and spectral hardness of the X-ray sources. The data are consistent with 9 of the 12, bright Ly$\alpha$ sources being driven by luminous, $10^{40}$ erg s$^{-1}$ X-ray sources. Half of the 8 Chandra sources are unresolved. The data suggest that nuclear activity, whether from LLAGN or nuclear starbursts, may play an important role in Ly$\alpha$ emission. Our results also suggest a significant link between Ly$\alpha$ emission and HMXBs, ULXs, and/or LLAGN, which would imply that Ly$\alpha$ may be generated over timescales 1 - 2 orders of magnitude longer than produced by photoionization from OB stars. This highlights a critical need to quantify the relative contributions of different sources across cosmic time, to interpret Ly$\alpha$ observations and the resulting properties of distant galaxies.

Context: Lyman $\alpha$ blobs (LABs) are large-scale radio-quiet Lyman $\alpha$ (Ly$\alpha$) nebula at high-$z$ that occur predominantly in overdense proto-cluster regions. Especially the prototypical SSA22a-LAB1 at $z=3.1$ has become an observational reference for LABs across the electromagnetic spectrum. Aims: We want to understand the powering mechanisms that drive the LAB to gain empirical insights into galaxy formation processes within a rare dense environment at high-$z$. Methods: LAB 1 was observed for 17.5h with the VLT/MUSE integral-field spectrograph. We produced optimally extracted narrow band images in Ly$\alpha$ $\lambda1216$, HeII $\lambda1640$, and we tried to detect CIV $\lambda1549$ emission. By using a moment based analysis we mapped the kinematics of the blob. Results: We detect Ly$\alpha$ emission to surface-brightness limits of $10^{-19}$erg s$^{-1}$cm$^{-2}$arcsec$^{-2}$. At this depth we reveal a bridge between LAB 1 and its northern neighbour LAB 8, as well as a shell-like filament towards the south of LAB 1. We find a coherent large scale east-west $\sim$1000 km s$^{-1}$ velocity gradient that is aligned perpendicular to the major axis of the blob. We detect HeII emission in three distinct regions, but we can only provide upper limits for CIV. Conclusions: Various gas excitation mechanisms are at play in LAB 1: Ionising radiation and feedback effects dominate near the embedded galaxies, while Ly$\alpha$ scattering is contributing at larger distances. The HeII/Ly$\alpha$ ratios combined with the upper limits on CIV/Ly$\alpha$ ratios favour shock models over photo-ionisation. The alignment of the angular momentum vector parallel to the morphological principal axis appears odds with the predicted norm for high-mass halos, but likely reflects that LAB\,1 resides at a node of multiple intersecting filaments of the cosmic web. (Abridged)

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

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.

Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.

The ionizing photons produced by massive stars are key actors in galaxy evolution. Ionizing photon production and escape is poorly understood. Improved space-based, spatially-resolved, multiplexed spectroscopic capabilities covering observed wavelengths of 1000 to 3000 Angstrom, in concert with spectroscopy from the ELTs and JWST, would lead to definitive answers as to how ionizing photons are produced and leaked, what populations of galaxies are responsible for ionizing photon leakage, what determines whether escape is possible, and how ionizing galaxy populations evolve over cosmic time.

Diagnosing the physical and chemical conditions within star-forming galaxies (SFGs) is of paramount importance to understanding key components of galaxy formation and evolution: star-formation, gas enrichment, outflows, and accretion. Well established optical emission-line diagnostics used to discern such properties (i.e., metal content, density, strength/shape of ionizing radiation) will be observationally inaccessible for the earliest galaxies, emphasizing the need for robust, reliable interstellar medium (ISM) diagnostics at ultraviolet (UV) wavelengths. Calibrating these UV diagnostics requires a comprehensive comparison of the UV and optical emission lines in nearby SFGs. Optical integral field unit (IFU) surveys have revealed the inhomogeneous nature of the ISM in SFGs, which leads to non-systematic biases in the interpretation of unresolved sources. Spatial variations are especially important to consider at UV wavelengths, where the strongest emission features originate from only the highest excitation regions of the nebula and are challenging to distinguish from competing high-ionization sources (e.g., shocks, AGN, etc.). Since surveys collecting large-scale optical integral field unit (IFU) spectroscopy are already underway, this white paper calls for an IFU or multi-object far-UV (FUV) spectroscopic instrument with high sensitivity, high spatial resolution, and large field of view (FoV). Given the impact of large-scale optical IFU surveys over the past decade, this white paper emphasizes the scientific need for a comparable foundation of spatially-resolved far-UV spectroscopy survey of nearby galaxies that will lay the foundation of diagnostics critical to the interpretation of the distant universe.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The aim of this paper is to introduce the WIRCam Ultra Deep Survey (WUDS), a near-IR photometric public survey carried out at the CFH Telescope in the field of the CFHTLS-D3 field (Groth Strip). WUDS includes four near-IR bands (Y, J, H and K_s) over a field of view of ~400 arcmin^2. The typical depth of WUDS data reaches between ~26.8 in Y and J, and ~26 in H and K_s (AB, 3 sigma in 1.3 arcsec aperture). The area and depth of this survey were specifically tailored to set strong constraints on the cosmic star formation rate and the luminosity function brighter or around L* in the z~6-10 redshift domain, although these data are also useful for a variety of extragalactic projects.This first paper is intended to present the properties of WUDS: catalog building, completeness and depth, number counts, photometric redshifts, and global properties of the galaxy population. We have also concentrated on the study of galaxy samples at z~[4.5-7] in this field. UV luminosity functions were derived at z~5 and z~6 taking advantage from the fact that WUDS covers a particularly interesting regime at intermediate luminosities, which allows a combined determination of M* and Phi* with increased accuracy. Our results on the luminosity function are consistent with a small evolution of both M* and Phi* between z=5 and z=6, irrespective of the method used to derive them, either photometric redshifts applied to blindly-selected dropout samples or the classical Lyman Break Galaxy color-preselected samples. Our results lend support to higher Phi* determinations at z=6 than usually reported. The selection and combined analysis of different galaxy samples at z>7 will be presented in a forthcoming paper. WUDS is intended to provide a robust database in the near-IR for the selection of targets for detailed spectroscopic studies, in particular for the EMIR/GTC GOYA Survey (Abridged)

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 CO J=4-3 line and 3 mm dust continuum observations of a 100 kpc-scale filamentary Ly{\alpha} nebula (SSA22 LAB18) at z=3.1 using the Atacama Large Millimeter/submillimeter Array (ALMA). We detected the CO J=4-3 line at a systemic z(CO)=3.093 {\pm} 0.001 at 11 {\sigma} from one of the ALMA continuum sources associated with the Ly{\alpha} filament. We estimated the CO J=4-3 luminosity of L'CO(4-3)=(2.3\pm0.2)x10^9 K km s^{-1} pc^2 for this CO source, which is one order of magnitude smaller than those of typical z>1 dusty star-forming galaxies (DSFGs) of similar far-infrared luminosity L(IR)~10^{12} Lsun. We derived a molecular gas mass of Mgas=(4.4^{+0.9}_{-0.6})x10^9 Msun and a star-formation rate of SFR=270\pm160 Msun yr^{-1}. We also estimated a gas depletion time of {\tau}(dep)=17\pm10 Myr, being shorter than those of typical DSFGs. It is suggested that this source is in a transition phase from DSFG to a gas-poor, early-type galaxy. From ALMA to Herschel multi-band dust continuum observations, we measured a dust emissivity index {\beta}=2.3\pm0.2, which is similar to those of local gas-poor, early-type galaxies. Such a high {\beta} can be reproduced by specific chemical compositions for interstellar dust at the submillimeter wavelengths from recent laboratory experiments. ALMA CO and multi-band dust continuum observations can constrain the evolutionary stage of high-redshift galaxies through {\tau}(dep) and {\beta}, and thus we can investigate dust chemical compositions even in the early Universe.

Lyman-$\alpha$ (Ly$\alpha$) is a powerful astrophysical probe. Not only is it ubiquitous at high redshifts, it is also a resonant line, making Ly$\alpha$ photons scatter. This scattering process depends on the physical conditions of the gas through which Ly$\alpha$ propagates, and these conditions are imprinted on observables such as the Ly$\alpha$ spectrum and its surface brightness profile. In this work, we focus on a less-used observable capable of probing any scattering process: polarization. We implement the density matrix formalism of polarization into the Monte Carlo radiative transfer code tlac. This allows us to treat it as a quantum mechanical process where single photons develop and lose polarization from scatterings in arbitrary gas geometries. We explore static and expanding ellipsoids, biconical outflows, and clumpy multiphase media. We find that photons become increasingly polarized as they scatter and diffuse into the wings of the line profiles, making scattered Ly$\alpha$ polarized in general. The degree and orientation of Ly$\alpha$ polarization depends on the kinematics and distribution of the scattering HI gas. We find that it generally probes spatial or velocity space asymmetries and aligns itself tangentially to the emission source. We show that the mentioned observables, when studied separately, can leave similar signatures for different source models. We conclude by revealing how a joint analysis of the Ly$\alpha$ spectra, surface brightness profiles, and polarization can break these degeneracies and help us extract unique physical information on galaxies and their environments from their strongest, most prominent emission line.

The Lyman Alpha Reference Sample (LARS) of 14 star-forming galaxies offers a wealth of insight into the workings of these local analogs to high-redshift star-forming galaxies. The sample has been well-studied in terms of LyA and other emission line properties, such as HI mass, gas kinematics, and morphology. We analyze deep surface photometry of the LARS sample in UBIK broadband imaging obtained at the Nordic Optical Telescope and the Canada-France-Hawaii Telescope, and juxtaposition their derived properties with a sample of local high-redshift galaxy analogs, namely, with blue compact galaxies (BCGs). We construct radial surface brightness and color profiles with both elliptical and isophotal integration, as well as RGB images, deep contours, color maps, a burst fraction estimate, and a radial mass-to-light ratio profile for each LARS galaxy. Standard morphological parameters like asymmetry, clumpiness, the Gini and M20 coefficients are [...] analyzed, as well as isophotal asymmetry profiles for each galaxy. [...] We compare the LARS to the properties of the BCG sample and highlight the differences. Several diagnostics indicate that the LARS galaxies have highly disturbed morphologies even at the level of the faintest isophotes [...]. The ground-based photometry [...] reveals previously unexplored isophotes [...]. The burst fraction estimate suggests a spatially more extended burst region in LARS than in the BCGs. [...] The galaxies in the LARS sample appear to be in earlier stages of a merger event compared to the BCGs. Standard morphological diagnostics like asymmetry, clumpiness, Gini and M20 coefficients cannot separate the two samples, although an isophotal asymmetry profile successfully captures the average difference in morphology. These morphological diagnostics do not show any correlation with the equivalent width or the escape fraction of Lyman Alpha. [abridged]

We examine the dust geometry and Ly{\alpha} scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Ly{\alpha} emitting and starbursting systems with Hubble Space Telescope Ly{\alpha}, H{\alpha}, and H{\beta} imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media while others of them show significantly lower Ly{\alpha} emission compared to their Balmer decrement. With the LARS imaging, we present Ly{\alpha}/H{\alpha} and H{\alpha}/H{\beta} maps with spatial resolutions as low as $\sim$ 40 pc, and use these data to show that in most galaxies, the dust geometry is best modeled by three distinct regions: a central core where dust acts as a screen, an annulus where dust is distributed in clumps, and an outer envelope where Ly{\alpha} photons only scatter. We show that the dust that affects the escape of Ly{\alpha} is more restricted to the galaxies' central regions, while the larger Ly{\alpha} halos are generated by scattering at large radii. We present an empirical modeling technique to quantify how much Ly{\alpha} scatters in the halo, and find that this "characteristic" scattering distance correlates with the measured size of the Ly{\alpha} halo. We note that there exists a slight anti-correlation between the scattering distance of Ly{\alpha} and global dust properties.

We study the heating mechanisms and Ly{\alpha} escape fractions of 35 Ly{\alpha} blobs (LABs) at z = 3.1 in the SSA22 field. Dust continuum sources have been identified in 11 of the 35 LABs, all with star formation rates (SFRs) above 100 Msun/yr. Likely radio counterparts are detected in 9 out of 29 investigated LABs. The detection of submm dust emission is more linked to the physical size of the Ly{\alpha} emission than to the Ly{\alpha} luminosities of the LABs. A radio excess in the submm/radio detected LABs is common, hinting at the presence of active galactic nuclei. Most radio sources without X-ray counterparts are located at the centers of the LABs. However, all X-ray counterparts avoid the central regions. This may be explained by absorption due to exceptionally large column densities along the line-of-sight or by LAB morphologies, which are highly orientation dependent. The median Ly{\alpha} escape fraction is about 3\% among the submm-detected LABs, which is lower than a lower limit of 11\% for the submm-undetected LABs. We suspect that the large difference is due to the high dust attenuation supported by the large SFRs, the dense large-scale environment as well as large uncertainties in the extinction corrections required to apply when interpreting optical data.

We report on the discovery of ionised gas filaments in the circum-galactic halo of the extremely metal-poor compact starburst SBS 0335-052E in a 1.5h integration with the MUSE integral-field spectrograph. We detect these features in H${\alpha}$ and [OIII] emission down to surface-brightness levels of $5 \times 10^{-19}$erg s$^{-1}$cm$^{-2}$arcsec$^{-2}$. The filaments have projected diameters of 2.1 kpc and extend more than 9 kpc to the north and north-west from the main stellar body. We also detect extended nebular HeII $\lambda$4686 emission that brightens towards the north-west at the rim of a star-burst driven super-shell, suggestive of a locally enhanced UV radiation field due to shocks. We also present a velocity field of the ionised gas. The filaments appear to connect seamlessly in velocity space to the kinematical disturbances caused by the shell. Similar to high-$z$ star-forming galaxies, the ionised gas in this galaxy is dispersion dominated. We argue that the filaments were created via feedback from the starburst and that these ionised structures in the halo may act as escape channels for Lyman continuum radiation in this gas-rich system.

Context. Type Ia supernovae (SNe Ia) can be used to address numerous questions in astrophysics and cosmology. Due to their well known spectral and photometric properties, SNe Ia are well suited to study gas and dust along the lines-of-sight to the explosions. For example, narrow Na I D and Ca II H&K absorption lines can be studied easily, because of the well-defined spectral continuum of SNe Ia around these features. Aims. We study the gas and dust along the line-of-sight to iPTF16abc, which occurred in an unusual location, in a tidal arm, 80 kpc from centre of the galaxy NGC 5221. Methods. Using a time-series of high-resolution spectra, we examine narrow Na I D and Ca II H&K absorption features for variations in time, which would be indicative for circumstellar (CS) matter. Furthermore, we take advantage of the well known photometric properties of SNe Ia to determine reddening due to dust along the line-of-sight. Results. From the lack of variations in Na I D and Ca II H&K, we determine that none of the detected absorption features originate from the CS medium of iPTF16abc. While the Na I D and Ca II H&K absorption is found to be optically thick, a negligible amount of reddening points to a small column of interstellar dust. Conclusions. We find that the gas along the line-of-sight to iPTF16abc is typical of what might be found in the interstellar medium (ISM) within a galaxy. It suggests that we are observing gas that has been tidally stripped during an interaction of NGC 5221 with one of its neighbouring galaxies in the past $\sim10^9$ years. In the future, the gas clouds could become the locations of star formation. On a longer time scale, the clouds might diffuse, enriching the circum-galactic medium (CGM) with metals. The gas profile along the line-of-sight should be useful for future studies of the dynamics of the galaxy group containing NGC 5221.

Recent work suggests that strong emission line, star-forming galaxies may be significant Lyman Continuum leakers. We combine archival HST broadband ultraviolet and optical imaging (F275W and F606W, respectively) with emission line catalogs derived from WFC3 IR G141 grism spectroscopy to search for escaping Lyman Continuum (LyC) emission from homogeneously selected $z\simeq$2.5 SFGs. We detect no escaping Lyman Continuum from SFGs selected on [OII] nebular emission (N=208) and, within a narrow redshift range, on [OIII]/[OII]. We measure 1$\sigma$ upper limits to the LyC escape fraction relative to the non-ionizing UV continuum from [OII] emitters, $f_{esc}<$5.6%, and strong [OIII]/[OII]$>$5 ELGs, $f_{esc}<$14.0%. Our observations are not deep enough to detect $f_{esc}\lesssim$10% typical of the low redshift Lyman continuum emitters. However, we find that this population represents a small fraction of the star-forming galaxy population at $z\simeq$2. Thus, unless the number of extreme emission line galaxies grows substantially to z$>$6, such galaxies may be insufficient for reionization. Deeper survey data in the rest-frame ionizing UV will be necessary to determine whether strong line ratios could be useful for pre-selecting LyC leakers at high redshift.

The radio-quiet quasar Q2059-360 at redshift $z=3.08$ is known to be close to a small Lyman $\alpha$ blob (LAB) and to be absorbed by a proximate damped Ly$\alpha$ (PDLA) system. Here, we present the Multi Unit Spectroscopic Explorer (MUSE) integral field spectroscopy follow-up of this quasi-stellar object (QSO). Our primary goal is to characterize this LAB in detail by mapping it both spatially and spectrally using the Ly$\alpha$ line, and by looking for high-ionization lines to constrain the emission mechanism. Combining the high sensitivity of the MUSE integral field spectrograph mounted on the Yepun telescope at ESO-VLT with the natural coronagraph provided by the PDLA, we map the LAB down to the QSO position, after robust subtraction of QSO light in the spectral domain. In addition to confirming earlier results for the small bright component of the LAB, we unveil a faint filamentary emission protruding to the south over about 80 pkpc (physical kpc); this results in a total size of about 120 pkpc. We derive the velocity field of the LAB (assuming no transfer effects) and map the Ly$\alpha$ line width. Upper limits are set to the flux of the N V $\lambda 1238-1242$, C IV $\lambda 1548-1551$, He II $\lambda 1640$, and C III] $\lambda 1548-1551$ lines. We have discovered two probable Ly$\alpha$ emitters at the same redshift as the LAB and at projected distances of 265 kpc and 207 kpc from the QSO; their Ly$\alpha$ luminosities might well be enhanced by the QSO radiation. We also find an emission line galaxy at $z=0.33$ near the line of sight to the QSO. This LAB shares the same general characteristics as the 17 others surrounding radio-quiet QSOs presented previously. However, there are indications that it may be centered on the PDLA galaxy rather than on the QSO.

Low- and intermediate mass galaxies are widely discussed as cause of reionization at redshift $z\sim10-6$. However, observational proof of galaxies that are leaking ionizing radiation (Lyman continuum; LyC) is a currently ongoing challenge and the list of LyC emitting candidates is still short. Tololo 1247-232 is among those very few galaxies with recently reported leakage. We performed intermediate resolution ultraviolet (UV) spectroscopy with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope and confirm ionizing radiation emerging from Tololo 1247-232. Adopting an improved data reduction procedure, we find that LyC escapes from the central stellar clusters, with an escape fraction of 1.5$\pm$0.5% only, i.e. the lowest value reported for the galaxy so far. We further make use of FUV absorption lines of Si II and Si IV as a probe of the neutral and ionized interstellar medium. We find that most of the ISM gas is ionized, likely facilitating LyC escape from density bounded regions. Neutral gas covering as a function of line-of-sight velocity is derived using the apparent optical depth method. The ISM is found to be sufficiently clumpy, supporting the direct escape of LyC photons. We further report on broadband UV and optical continuum imaging as well as narrowband imaging of Ly$\alpha$, H$\alpha$ and H$\beta$. Using stellar population synthesis, a Ly$\alpha$ escape fraction of 8% was derived. We also performed VLA 21cm imaging. The hydrogen hyperfine transition was not detected, but a deep upper limit atomic gas mass of $\lesssim10^9 M_{\odot}$ could be derived. The upper limit gas fraction defined as $\frac{M_{H I}}{M_*}$ is only 20%. Evidence is found that the H I gas halo is relatively small compared to other Lyman Alpha emitters.

Star forming galaxies are believed to be a major source of Lyman Continuum (LyC) radiation responsible for reionizing the early Universe. Direct observations of escaping ionizing radiation have however been few and with low escape fractions. In the local Universe, only ~10 emitters have been observed, with typical escape fractions of a few percent. The mechanisms regulating this escape need to be strongly evolving with redshift in order to account for the Epoch of Reionization. Gas content and star formation feedback are among the main suspects, known to both regulate neutral gas coverage and evolve with cosmic time. In this paper, we reanalyze HST-COS spectrocopy of the first detected local LyC leaker, Haro 11. We examine the connection between LyC leakage and Lyman-$\alpha$ line shape, and feedback-influenced neutral ISM properties like kinematics and gas distribution. We discuss the two extremes of an optically thin, density bounded ISM and a riddled, optically thick, ionization bounded ISM, and how Haro 11 fits into their theoretical predictions. We find that the most likely ISM model is a clumpy neutral medium embedded in a highly ionized medium with a combined covering fraction of unity and a residual neutral gas column density in the ionized medium high enough to be optically thick to Lyman-$\alpha$, but low enough to be at least partly transparent to Lyman continuum and undetected in Si II. This suggests that SF feedback and galaxy-scale interaction events play a major role in opening passageways for ionizing radiation through the neutral medium.

We report the result from observations conducted with the Atacama Large Millimeter/submillimeter Array (ALMA) to detect [CII] 158 um fine structure line emission from galaxies embedded in one of the most spectacular Lyman-alpha blobs (LABs) at z=3.1, SSA22-LAB1. Of three dusty star-forming galaxies previously discovered by ALMA 860 um dust continuum survey toward SSA22-LAB1, we detected the [CII] line from one, LAB1-ALMA3 at z=3.0993+/-0.0004. No line emission was detected, associated with the other ALMA continuum sources or from three rest-frame UV/optical selected z_spec~3.1 galaxies within the field of view. For LAB1-ALMA3, we find relatively bright [CII] emission compared to the infrared luminosity (L_[CII]/L_[CII]) and an extremely high [CII] 158 um and [NII] 205 um emission line ratio (L_[CII]/L_[NII]>55). The relatively strong [CII] emission may be caused by abundant photodissociation regions and sub-solar metallicity, or by shock heating. The origin of the unusually strong [CII] emission could be causally related to the location within the giant LAB, although the relationship between extended Lyman-alpha emission and ISM conditions of associated galaxies is yet to be understand.

We present the first robust detection of HI 21 cm emission in the blue compact galaxy Haro 11 using the 100m Robert C. Byrd Green Bank Telescope (GBT). Haro 11 is a luminous blue compact galaxy with emission in both Lyman Alpha and the Lyman continuum. We detect (5.1 $\pm$ 0.7 $\times$10$^8$) M$_{\odot}$ of HI gas at an assumed distance of 88 Mpc, making this galaxy HI deficient compared to other local galaxies with similar optical properties. Given this small HI mass, Haro 11 has an elevated M$_{H2}$/M$_{HI}$ ratio and a very low gas fraction compared to most local galaxies, and contains twice as much mass in ionized hydrogen as in neutral hydrogen. The HI emission has a linewidth of 71 kms$^{-1}$ and is offset 60 kms$^{-1}$ redward of the optical line center. It is undergoing a starburst after a recent merger which has elevated the star formation rate, and will deplete the gas supply in $<$ 0.2 Gyr. Although this starburst has elevated the SFR compared to galaxies with similar HI masses and linewidths, Haro 11 matches a trend of lower gas fractions toward higher star formation rates and is below the general trend of increasing HI mass with increasing luminosity. Taken together, our results paint Haro 11 as a standard low-mass galaxy that is undergoing an unusually efficient star formation episode.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

We report results from a new HST study of the OVI 1032,1038\AA\ doublet in emission around intensely star-forming galaxies. The programme aims to characterize the energy balance in starburst galaxies and gas cooling in the difficult-to-map coronal temperature regime of 2-5 x $10^5$K. We present the first resolved image of gas emission in the OVI line. Our target, SDSS J1156+5008, is very compact in the continuum but displays OVI emission to radii of 23 kpc. The surface brightness profile is well fit by an exponential with a scale of 7.5kpc. This is ten times the size of the photoionized gas, and we estimate that 1/6 the total OVI luminosity comes from resonantly scattered continuum radiation. Spectroscopy - which closely resembles a stacked sample of archival spectra - confirms the OVI emission, and determines the column density and outflow velocity from blueshifted absorption. The combination of measurements enables several new calculations with few assumptions. The OVI regions fill only ~$10^{-3}$ of the volume. By comparing the cooling time with the cloud sound-crossing time, the cooling distance with the size, and the pressure in the OVI and nebular gas, we conclude that the OVI-bearing gas cannot have been lifted to the scale height at this temperature, and must be cooling in situ through this coronal temperature regime. The coronal phase contains ~1% of the ionized mass, and its kinetic energy is currently ~1% of the budget set by supernova feedback. However a much larger amount of the gas must have cooled through this phase during the star formation episode. The outflow exceeds the escape velocity and the gas may become unbound, but it will recombine before it escapes and become visible to Lyman (and OI) spectroscopy. The mapping of this gas represents a crucial step in further constraining galaxy formation scenarios and guiding the development of future satellites.

We present the Spitzer-IRAC/MIPS Extragalactic survey (SIMES) in the South Ecliptic Pole (SEP) field. The large area covered (7.7 deg$^2$), together with one of the lowest Galactic cirrus emissions in the entire sky and a very extensive coverage by Spitzer, Herschel, Akari, and GALEX, make the SIMES field ideal for extragalactic studies. The elongated geometry of the SIMES area ($\approx$4:1), allowing for a significant cosmic variance reduction, further improves the quality of statistical studies in this field. Here we present the reduction and photometric measurements of the Spitzer/IRAC data. The survey reaches a depth of 1.93 and 1.75 $\mu$Jy (1$\sigma$) at 3.6 and 4.5 $\mu$m, respectively. We discuss the multiwavelength IRAC--based catalog, completed with optical, mid-- and far--IR observations. We detect 341,000 sources with F$_{3.6\mu m} \geq 3\sigma$. Of these, 10% have an associated 24 $\mu$m counterpart, while 2.7% have an associated SPIRE source. We release the catalog through the NASA/IPAC Infrared Science Archive (IRSA). Two scientific applications of these IRAC data are presented in this paper: first we compute integral number counts at 3.6 $\mu$m. Second, we use the [3.6]--[4.5] color index to identify galaxy clusters at z$>$1.3. We select 27 clusters in the full area, a result consistent with previous studies at similar depth.

We present a follow-up study to the imaging polarimetry performed by Hayes et. al. 2011 on LAB1 in the SSA22 protocluster region. Arguably the most well-known Lyman-$\alpha$ "blob", this radio-quiet emission-line nebula likely hosts a galaxy which is either undergoing significant star formation or hosts an AGN, or both. We obtain deep, spatially resolved spectro-polarimetry of the Ly$\alpha$ emission and detect integrated linear polarization of $9$-$13\%\pm2$-$3\%$ at a distance of approximately 15 kpc north and south of the peak of the Lyman-$\alpha$ surface brightness with polarization vectors lying tangential to the galactic central source. In these same regions, we also detect a wavelength dependence in the polarization which is low at the center of the Ly$\alpha$ line profile and rises substantially in the wings of the profile. These polarization signatures are easily explained by a weak out-flowing shell model. The spectral dependence of the polarization presented here provide a framework for future observations and interpretations of the southern portion of LAB1 in that any model for this system must be able to reproduce this particular spectral dependence. However, questions still remain for the northern-most spur of LAB1. In this region we detect total linear polarization of between $3$ and $20\%$ at the $5\%$ significance level. Simulations predict that polarization should increase with radius for a symmetric geometry. That the northern spur does not suggests either that this region is not symmetric (which is likely) and exhibits variations in columns density, or that it is kinematically distinct from the rest of LAB1 and powered by another mechanism altogether.

To date no direct detection of Lyman continuum emission has been measured for intermediate--redshift z~1 star-forming galaxies . We combine HST grism spectroscopy with GALEX UV and ground--based optical imaging to extend the search for escaping Lyman continuum to a large (~600) sample of z~1 low-mass, moderately star-forming galaxies selected initially on H$\alpha$ emission. The characteristic escape fraction of LyC from SFGs that populate this parameter space remains weakly constrained by previous surveys, but these faint SFGs are assumed to play a significant role in the reionization of neutral hydrogen in the intergalactic medium (IGM) at high redshift (z>6). We do not make an unambiguous detection of escaping LyC radiation from this $z\sim1$ sample, individual non--detections to constrain the absolute Lyman continuum escape fraction, $f_{esc}$<2.1% (3$\sigma$). We measure upper limits of $f_{esc}$<9.6% from a sample of SFGs selected on high H$\alpha$ equivalent width (EW>200\AA), which are thought to be close analogs of high redshift sources of reionization. For reference, we also present an emissivity--weighted escape fraction which is useful as a measurement of the general contribution of the SFGs to the z~1 ionizing UV background. In the discussion, we consider the implications of these intermediate redshift constraints for the reionization of hydrogen in the intergalactic medium at high (z > 6) redshift. If we assume our $z\sim1$ SFGs, for which we measure this emissivity-weighted $f_{esc}$, are analogs to the high redshift sources of reionization, we find is difficult reconcile reionization by faint (M$_{UV}<-13$) SFGs with a low escape fraction ($f_{esc}<$3%), with constraints from independent high redshift observations. If $f_{esc}$ evolves with redshift, reionization by SFGs may be consistent with observations from Planck.

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

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.

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

The Lya emission line of HI is intrinsically the brightest feature in the spectrum of astrophysical nebulae, making it a very attractive observational tool with which to survey galaxies. Moreover as a UV resonance line, Lya possesses several unique characteristics that make it useful to study the ISM and ionizing stellar population at all cosmic epochs. In this review I present a summary of Lya observations of galaxies in the nearby universe. At UV magnitudes reachable with current facilities, only ~5% of the local galaxy population shows a Lya equivalent width (EW_Lya) that exceeds 20\AA. This fraction increases dramatically at higher z, but only in the local universe can we study galaxies in detail and assemble unprecedented multi-wavelength datasets. I discuss many local Lya observations, showing that when galaxies show net Lya emission, they ubiquitously produce large halos of scattered Lya, that dominate the integrated luminosity. We discuss how global EW_Lya and the Lya escape fraction (fescLya) are higher (EW_Lya >~ 20\AA\ and fescLya> 10%) in galaxies that represent the less massive and younger end of the distributions for local objects. This is connected with various properties, such that Lya-emitters have lower metallicities (median value of 12+log(O/H) ~ 8.1) and dust reddening. However, the presence of galactic outflows is also vital to Doppler shift the Lya line out of resonance with the HI, as high EW_Lya is found only among galaxies with winds faster than ~50km/s. The evidence is then assembled into a coherent picture, and the requirement for star formation driven feedback is discussed in the context of an evolutionary sequence where the ISM is accelerated and/or subject to fluid instabilities, which reduce the scattering of Lya. Concluding remarks take the form of perspectives upon the most pressing questions that can be answered by observation.

Context. Lyman alpha emitting galaxies (LAEs) are used to probe the distant universe and are therefore important for galaxy evolution studies and for providing clues to the nature of the epoch of reionization, but the exact circumstances under which Lyman alpha escapes a galaxy are still not fully understood. Aims. The Trident project is designed to simultaneously examine Lyman alpha, H-alpha and Lyman Continuum emission from galaxies at redshift z~2, thus linking together these three aspects of ionising radiation in galaxies. In this paper, we outline the strategy of this project and examine the properties of LAEs in the GOODS North field. Methods. We performed a narrowband LAE survey in GOODS North using existing and two custom made filters at the Nordic Optical Telescope with MOSCA. We use complementary broad band archival data in the field to make a careful candidate selection and perform optical to near-IR SED fitting. We also estimate far-infrared luminosities by matching our candidates to detections in Spitzer/MIPS 24{\mu}m and Herschel/PACS catalogs. Results. We find a total of 25 LAE candidates, probing mainly the bright end of the LAE luminosity function with L_Ly {\alpha} ~ 1-15e42 erg/s. They display a range of masses of ~0.5-50e9 M_solar, and average ages from a few tens of Myr to 1 Gyr when assuming a constant star formation history. The majority of our candidates also show signs of recent elevated star formation. Three candidates have counterparts in the GOODS-Herschel far-IR catalogue, with luminosities consistent with ultra-luminous infrared galaxies (ULIRGs). Conclusions. The wide range of parameters derived from our SED fitting, as well as part of our sample being detected as ULIRGs, seems to indicate that at these Lyman alpha luminosities, LAEs do not necessarily have to be young dwarfs, and that a lack of dust is not required for Lyman alpha to escape.

We propose to infer ionising continuum leaking properties of galaxies by looking at their Lyman-alpha line profiles. We carry out Lyman-alpha radiation transfer calculations in two models of HII regions which are porous to ionising continuum escape: 1) the so-called "density bounded" media, in which massive stars produce enough ionising photons to keep the surrounding interstellar medium transparent to the ionising continuum, i.e almost totally ionised, and 2) "riddled ionisation-bounded" media, surrounded by neutral interstellar medium, but with holes, i.e. with a covering factor lower than unity. The Lyman-alpha spectra emergent from these configurations have distinctive features: 1) a "classical" asymmetric redshifted profile in the first case, but with a small shift of the maximum of the profile compare to the systemic redshift (Vpeak < 150 km/s); 2) a main peak at the systemic redshift in the second case (Vpeak = 0 km/s), with, as a consequence, a non-zero Lyman-alpha flux bluewards the systemic redshift. Assuming that in a galaxy leaking ionising photons, the Lyman-alpha component emerging from the leaking star cluster(s) dominates the total Lyman-alpha spectrum, the Lyman-alpha shape may be used as a pre-selection tool to detect Lyman continuum (LyC) leaking galaxies, in objects with well determined systemic redshift, and high spectral resolution Lyman-alpha spectra (R >= 4000). The examination of a sample of 10 local starbursts with high resolution HST-COS Lyman-alpha spectra and known in the literature as LyC leakers or leaking candidates, corroborates our predictions. Observations of Lyman-alpha profiles at high resolution should show definite signatures revealing the escape of Lyman continuum photons from star-forming galaxies.

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 explore the properties of high-redshift Lyman-alpha emitters (LAE), and their link with the Lyman-Break galaxy population (LBG), using a semi-analytic model of galaxy formation that takes into account resonant scattering of Lya photons in gas outflows. We can reasonably reproduce the abundances of LAEs and LBGs from redshift 3 to 7, as well as most UV LFs of LAEs. The stronger dust attenuation for (resonant) Lya photons compared to UV continuum photons in bright LBGs provides a natural interpretation to the increase of the LAE fraction in LBG samples, X_LAE, towards fainter magnitudes. The redshift evolution of X_LAE seems however very sensitive to UV magnitudes limits and EW cuts. In spite of the apparent good match between the statistical properties predicted by the model and the observations, we find that the tail of the Lya equivalent width distribution (EW > 100 A) cannot be explained by our model, and we need to invoke additional mechanisms. We find that LAEs and LBGs span a very similar dynamical range, but bright LAEs are about 4 times rarer than LBGs in massive halos. Moreover, massive halos mainly contain weak LAEs in our model, which might introduce a bias towards low-mass halos in surveys which select sources with high EW cuts. Overall, our results are consistent with the idea that LAEs and LBGs make a very similar galaxy population. Their apparent differences seem mainly due to EW selections, UV detection limits, and a decreasing Lya-to-UV escape fraction ratio in high SFR galaxies.

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.

We measured the sizes and morphological parameters of LARS galaxies in the continuum, Lya, and Ha images. We studied morphology by using the Gini coefficient vs M20 and asymmetry vs concentration diagrams. We then simulated LARS galaxies at z~2 and 5.7, performing the same morphological measurements. We also investigated the detectability of LARS galaxies in current deep field observations. The subsample of LAEs within LARS (LARS-LAEs) was stacked to provide a comparison to stacking studies performed at high redshift. LARS galaxies have continuum size, stellar mass, and rest-frame absolute magnitude typical of Lyman break analogues in the local Universe and also similar to 2<z<3 star-forming galaxies and massive LAEs. LARS optical morphology is consistent with the one of merging systems, and irregular or starburst galaxies. For the first time we quantify the morphology in Lya images: even if a variety of intrinsic conditions of the interstellar medium can favour the escape of Lya photons, LARS-LAEs appear small in the continuum, and their Lya is compact. LARS galaxies tend to be more extended in Lya than in the rest-frame UV. It means that Lya photons escape by forming haloes around HII regions of LARS galaxies. The stack of LARS-LAE Lya images is peaked in the centre, indicating that the conditions, which make a galaxy an LAE, tend to produce a concentrated surface brightness profile. On the other hand, the stack of all LARS galaxies is shallower and more extended. This can be caused by the variety of dust and HI amount and distribution, which produces a more complex, patchy, and extended profile, like the one observed for Lyman break galaxies that can contribute to the stack. We cannot identify a single morphological property that controls whether a galaxy emits a net positive Lya flux. However, the LARS-LAEs have continuum properties consistent with merging systems.

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.

We study the prospects for constraining the ionized fraction of the intergalactic medium (IGM) at $z>6$ with the next generation of large Ly$\alpha$ emitter surveys. We make predictions for the upcoming Subaru Hyper Suprime-Cam (HSC) Ly$\alpha$ survey and a hypothetical spectroscopic survey performed with the James Webb Space Telescope (JWST). Considering various scenarios where the observed evolution of the Ly$\alpha$ luminosity function of Ly$\alpha$ emitters at $z>6$ is explained partly by an increasingly neutral IGM and partly by intrinsic galaxy evolution, we show how clustering measurements will be able to distinguish between these scenarios. We find that the HSC survey should be able to detect the additional clustering induced by a neutral IGM if the global IGM neutral fraction is greater than $\sim$20 per cent at $z=6.5$. If measurements of the Ly$\alpha$ equivalent widths (EWs) are also available, neutral fractions as small as 10 per cent may be detectable by looking for correlation between the EW and the local number density of objects. In this case, if it should turn out that the IGM is significantly neutral at $z=6.5$ and the intrinsic EW distribution is relatively narrow, the observed EWs can also be used to construct a map of the locations and approximate sizes of the largest ionized regions. For the JWST survey, the results appear a bit less optimistic. Since such surveys probe a large range of redshifts, the effects of the IGM will be mixed up with any intrinsic galaxy evolution that is present, making it difficult to disentangle the effects. However, we show that a survey with the JWST will have a possibility of observing a large group of galaxies at $z\sim7$, which would be a strong indication of a partially neutral IGM.

Though Lyman alpha emission (Lya) is one of the most used tracers of massive star formation at high redshift, a correct understanding of radiation transfer effects by neutral gas is required to properly quantify the star formation rate along the history of the Universe. We are embarked in a program to study the properties of the Lya emission (spectral profile, spatial distribution, relation to Balmer lines intensity,...) in several local starburst galaxies. We present here the results obtained for IRAS 08339+6517. Using evolutionary population synthesis models, we have characterized the properties of the starburst (UV continuum, Halpha, total infrared and X-ray emissions, etc.), which transformed 1.4e+8 Mo of gas into stars around 5-6 Myr ago. In addition to the central compact emission blob, we have identified a diffuse Lya emission component smoothly distributed over the whole central area of IRAS 08339+6517. This diffuse emission is spatially decoupled from the UV continuum, the Halpha emission or the Halpha/Hbeta ratio. Both locally and globally, the Lya/Halpha ratio is lower than the Case B predictions, even after reddening correction, with an overall Lya escape fraction of only 4%. We conclude that in IRAS 08339+6517 the resonant scattering of Lya photons by an outflowing shell of neutral gas causes their highly-efficient destruction by dust, which explains the low Lya escape fraction measured. These results stress again the importance of a proper correction of scattering and transfer effects when using Lya to derive the star formation rate in high-redshift galaxies.

[abridged] Among the different observational techniques used to select high-redshift galaxies, the hydrogen recombination line Lyman-alpha (Lya) is of particular interest as it gives access to the measurement of cosmological quantities such as the star formation rate of distant galaxy populations. However, the interpretation of this line and the calibration of such observables is still subject to serious uncertainties. Therefore, it important to understand under what conditions the Lya line can be used as a reliable star formation diagnostic tool. We use a sample of 24 Lya emitters at z ~ 0.3 with an optical spectroscopic follow-up to calculate the Lya escape fraction and its dependency upon different physical properties. We also examine the reliability of Lya as a star formation rate indicator. We combine these observations with a compilation of Lya emitters selected at z = 0 - 0.3 to assemble a larger sample. The Lya escape fraction depends clearly on the dust extinction following the relation fesc(Lya) = C(Lya) x 10^(-0.4 E(B-V) k(Lya)), but with a shallower slope than previously reported, with k(Lya) ~ 6.67 and C(Lya) = 0.22. However, the correlation does not follow the expected curve for a simple dust attenuation. We explore the various mechanisms than lead to fesc(Lya) values above the continuum extinction curve, i.e. to an enhancement of the Lya output. We also observe that the strength of Lya and the escape fraction appear unrelated to the galaxy metallicity. Regarding the reliability of Lya as a star formation rate (SFR) indicator, we show that the deviation of SFR(Lya) from the true SFR (as traced by the UV continuum) is a function of the observed SFR(UV), which can be seen as the decrease of fesc(Lya) with increasing UV luminosity. Moreover, we observe a redshift-dependence of this relationship revealing the underlying evolution of fesc(Lya) with redshift.

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.

Context. The Lyman alpha emission line of galaxies is an important tool for finding galaxies at high redshift, and thus probe the structure of the early universe. However, the resonance nature of the line and its sensitivity to dust and neutral gas is still not fully understood. Aims. We present measurements of the velocity, covering fraction and optical depth of neutral gas in front of two well known local blue compact galaxies that show Lyman alpha in emission: ESO 338-IG 04 and Haro 11. We thus test observationally the hypothesis that Lyman alpha can escape through neutral gas by being Doppler shifted out of resonance. Methods. We present integral field spectroscopy from the GIRAFFE/Argus spectrograph at VLT/FLAMES in Paranal, Chile. The excellent wavelength resolution allows us to accurately measure the velocity of the ionized and neutral gas through the H-alpha emission and Na D absorption, which traces the ionized medium and cold interstellar gas, respectively. We also present independent measurements with the VLT/X-shooter spectrograph which confirm our results. Results. For ESO 338-IG04, we measure no significant shift of neutral gas. The best fit velocity is -15 (16) km/s. For Haro 11, we see an outflow from knot B at 44 (13) km/s and infalling gas towards knot C with 32 (12) km/s. Based on the relative strength of the Na D absorption lines, we estimate low covering fractions of neutral gas (down to 10%) in all three cases. Conclusions. The Na D absorption likely occurs in dense clumps with higher column densities than where the bulk of the Ly-alpha scattering takes place. Still, we find no strong correlation between outflowing neutral gas and a high Lyman alpha escape fraction. The Lyman alpha photons from these two galaxies are therefore likely escaping due to a low column density and/or covering fraction.

The escape fraction of hydrogen ionizing photons (f_esc) from galaxies has been suggested to be evolving with time, but the picture is far from clear. While evidence for significant escape fractions has been found at high redshifts in several studies, the picture looks different in the more nearby universe. Here, we apply a new background subtraction routine on archival data from the Far Ultraviolet Spectroscopic Explorer (FUSE), in order to study local galaxies in search for possible Lyman Continuum (LyC) leakage. In the process, for the first time a stacked spectrum in the LyC is produced for local galaxies. With this small sample, we also make a more tentative approach to look for possible correlations between f_esc and physical parameters such as internal absorption E(B-V)_i, mass, H I mass, specific star formation rate (SSFR), metallicity, and Ly-alpha emission. Eight star forming galaxies with redshifts z > 0.015 from the FUSE archive were re-examined. Also, a sub-sample of an additional four galaxies with lower redshifts were included, for which the escape fraction was estimated from residual flux in the low ionization interstellar C II(1036{\AA}) line. Out of the eight galaxies, only one was found to have significant LyC leakage, Tol 1247-232 (S/N=5.2). This is the second detection of a leaking galaxy in the local universe. We find an absolute escape fraction for Tol 1247-232 of f_esc = 2.4(+0.9/-0.8) %. The stacked sample show an excess in the LyC with f_esc = 1.4(+0.6/-0.5) %, but we note that there might be important selection biases involved. With the small sample, we suggest a possible trend for higher f_esc with lower mass and with enhanced SSFR. None of the galaxies with high values of E(B-V)_i were found to show any sign of leakage.

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.

To investigate the ingredients, which allow star-forming galaxies to present Lyalpha line in emission, we studied the kinematics and gas phase metallicity (Z) of the interstellar medium. We used multi-object NIR spectroscopy with Magellan/MMIRS to study nebular emission from z=2-3 star-forming galaxies discovered in 3 MUSYC fields. We detected emission lines from four active galactic nuclei and 13 high-z star-forming galaxies, including Halpha lines down to a flux of 4.E-17 erg/sec/cm^2. This yielded 7 new redshifts. The most common emission line detected is [OIII]5007, which is sensitive to Z. We were able to measure Z for 2 galaxies and to set upper(lower) limits for another 2(2). The Z values are consistent with 0.3<Z/Zsun<1.2. Comparing the Lyalpha central wavelength with the systemic redshift, we find Delta_v(Lyalpha-[OIII])=70-270 km/sec. High-redshift star-forming galaxies, Lyalpha emitting (LAE) galaxies, and Halpha emitters appear to be located in the low mass, high star-formation rate (SFR) region of the SFR versus stellar mass diagram, confirming that they are experiencing burst episodes of star formation, which are building up their stellar mass. Their Zs are consistent with the relation found for z<2.2 galaxies in the Z versus stellar mass plane. The measured Delta_v(Lyalpha-[OIII]) values imply that outflows of material, driven by star formation, could be present in the z=2-3 LAEs of our sample. Comparing with the literature, we note that galaxies with lower Z than ours are also characterized by similar Delta_v(Lyalpha-[OIII]) velocity offsets. Strong [OIII] is detected in many Lyalpha emitters. Therefore, we propose the Lyalpha/[OIII] flux ratio as a tool for the study of high-z galaxies; while influenced by Z, ionization, and Lyalpha radiative transfer in the ISM, it may be possible to calibrate this ratio to primarily trace one of these effects.

The timing and duration of the reionization epoch is crucial to the emergence and evolution of structure in the universe. The relative roles that star-forming galaxies, active galactic nuclei and quasars play in contributing to the metagalactic ionizing background across cosmic time remains uncertain. Deep quasar counts provide insights into their role, but the potentially crucial contribution from star-formation is highly uncertain due to our poor understanding of the processes that allow ionizing radiation to escape into the intergalactic medium (IGM). The fraction of ionizing photons that escape from star-forming galaxies is a fundamental free parameter used in models to "fine-tune" the timing and duration of the reionization epoch that occurred somewhere between 13.4 and 12.7 Gyrs ago (redshifts between 12 > z > 6). However, direct observation of Lyman continuum (LyC) photons emitted below the rest frame \ion{H}{1} ionization edge at 912 \AA\ is increasingly improbable at redshifts z > 3, due to the steady increase of intervening Lyman limit systems towards high z. Thus UV and U-band optical bandpasses provide the only hope for direct, up close and in depth, observations of the types of environment that favor LyC escape. By quantifying the evolution over the past 11 billion years (z < 3) of the relationships between LyC escape and local and global parameters ..., we can provide definitive information on the LyC escape fraction that is so crucial to answering the question of, how did the universe come to be ionized? Here we provide estimates of the ionizing continuum flux emitted by "characteristic" (L_{uv}^*) star-forming galaxies as a function of look back time and escape fraction, finding that at z = 1 (7.6 Gyrs ago) L_{uv}^* galaxies with an escape fraction of 1% have a flux of 10^{-19} ergs cm^{-2} s^{-1} \AA^{-1}.

In order to use Lyman alpha (Lya) emission as star formation tracer in cosmological studies, we must understand how the resonant scattering affects the escape fraction of the Lya photons. Thus, high spatial resolution multiwavelength studies of nearby Lya emitters, like Haro 2, are highly needed. For that purpose, we have used Chandra X-ray and HST (UV, optical and NIR) images of Haro 2, and STIS and ground-based spectral images along its major and minor axes, to characterize the Lya emission and the properties of the stellar population. The UV, Ha (Halpha) and FIR luminosities of the Haro 2 nuclear starburst are reproduced using evolutionary synthesis models assuming a young stellar population with ages ~3.5-5.0 Myr, affected by differential interstellar extinctions. The observed X-ray emission is attributed to gas heated by the mechanical energy released by the starburst (soft component) and a Ultra-Luminous X-ray source candidate (hard). Both compact and diffuse Lya components are observed. Whereas Lya is spatially decoupled from Balmer lines emission, Balmer decrement and UV continuum, the diffuse Lya component is spatially correlated with the diffuse soft X-ray emission. Moreover, unlike the compact one, diffuse Lya shows luminosities larger than predicted from Ha, assuming case B recombination and dust extinction as derived from Ha/Hbeta. We propose that, whereas the compact Lya emission is associated to the massive stellar clusters and is affected by outflows and dust extinction, the diffuse Lya originates in gas ionized by the hot plasma responsible for the soft X-ray radiation, as suggested by their spatial correlation and by the measured L(Ha)/LsoftX ratios. Calibration of Lya as star formation rate tracer should therefore include both effects (destruction vs. enhancement) to avoid biases in the study of galaxies at cosmological distances.

We present supernova rate measurements at redshift 0.1-1.0 from the Stockholm VIMOS Supernova Survey (SVISS). The sample contains 16 supernovae in total. The discovered supernovae have been classified as core collapse or type Ia supernovae (9 and 7, respectively) based on their light curves, colour evolution and host galaxy photometric redshift. The rates we find for the core collapse supernovae are 3.29 (-1.78,-1.45)(+3.08,+1.98) x 10^-4 yr^-1 Mpc^-3 h70^3 (with statistical and systematic errors respectively) at average redshift 0.39 and 6.40 (-3.12,-2.11)(+5.30,+3.65) x 10^-4 yr^-1 Mpc^-3 h70^3 at average redshift 0.73. For the type Ia supernovae we find a rate of 1.29 (-0.57,-0.28)(+0.88,+0.27) x 10^-4 yr^-1 Mpc^-3 h70^3 at average redshift 0.62. All of these rate estimates have been corrected for host galaxy extinction, using a method that includes supernovae missed in infrared bright galaxies at high redshift. We use Monte Carlo simulations to make a thorough study of the systematic effects from assumptions made when calculating the rates and find that the most important errors come from misclassification, the assumed mix of faint and bright supernova types and uncertainties in the extinction correction. We compare our rates to other observations and to the predicted rates for core collapse and type Ia supernovae based on the star formation history and different models of the delay time distribution. Overall, our measurements, when taking the effects of extinction into account, agree quite well with the predictions and earlier results. Our results highlight the importance of understanding the role of systematic effects, and dust extinction in particular, when trying to estimate the rates of supernovae at moderate to high redshift.

Laporte et al. (2011) reported a very high redshift galaxy candidate: a lensed J-band dropout (A2667-J1). J1 has a photometric redshift of z=9.6-12, the probability density function for which permits no low or intermediate z solution. We here report new spectroscopic observations of this galaxy with VLT/XShooter, which show clear [OIII]5007AA, Ly-alpha, H-alpha, and H-beta emission and place the galaxy firmly at z=2.082. The oxygen lines contribute only ~25% to the H-band flux, and do not significantly affect the dropout selection of J1. After correcting the broadband fluxes for line emission, we identify two roughly equally plausible natures for A2667-J1: either it is young heavily reddened starburst, or a maximally old system with a very pronounced 4000AA break, upon which a minor secondary burst of star formation is superimposed. Fits show that to make a 3 sigma detection of this object in the B-band (V-band), imaging of depth AB=30.2 (29.5) would be required - despite the relatively bright NIR magnitude, we would need optical data of equivalent depth to the Hubble Ultra Deep Field to rule out the mid-z solution on purely photometric grounds. Assuming that this stellar population can be scaled to the NIR magnitudes of recent HST/WFC3 IR-selected galaxies, we conclude that infeasibly deep optical data AB~32 would be required for the same level of security. There is a population of galaxies at z~2 with continuum colours alone that mimic those of our z=7-12 candidates.

We report new constraints on the galaxy luminosity function at z~9 based on observations carried out with ESO/VLT FORS2, HAWK-I and X-Shooter around the lensing cluster A2667, as part of our project aimed at selecting z~7-10 candidates accessible to spectroscopy. Only one J-dropout source was selected in this field fulfilling the color and magnitude criteria. This source was recently confirmed as a mid-z interloper based on X-Shooter spectroscopy. The depth and the area covered by our survey are well suited to set strong constraints on the bright-end of the galaxy luminosity function and hence on the star formation history at very high redshift. The non-detection of reliable J-dropout sources over the ~36arcmin2 field of view towards A2667 was used to carefully determine the lens-corrected effective volume and the corresponding upper-limit on the density of sources. The strongest limit is obtained for Phi(M_{1500}=-21.4+/-0.50)<6.70x10^{-6}Mpc^{-3}mag^{-1} at z~9. A maximum-likelihood fit of the luminosity function using all available data points including the present new result yields M*>-19.7 with fixed alpha=-1.74 and Phi*=1.10x10^{-3}Mpc^{-3}. The corresponding star formation rate density should be rho_{SFR}<5.97x10^{-3}M_{solar}/yr/Mpc^{3} at z~9. These results are in good agreement with the most recent estimates already published in this range of redshift and for this luminosity domain. This new result confirms the decrease in the density of luminous galaxies at very high-redshift, hence providing strong constraints for the design of future surveys aiming to explore the very high-redshift Universe.

We present a new model for high redshift Lyman-Alpha Emitters (LAEs) in the cosmological context which takes into account the resonant scattering of Ly-a photons through expanding gas. The GALICS semi-analytic model provides us with the physical properties of a large sample of high redshift galaxies. We implement a gas outflow model for each galaxy based on simple scaling arguments. The coupling with a library of numerical experiments of Ly-a transfer through expanding or static dusty shells of gas allows us to derive the Ly-a escape fractions and profiles. The predicted distribution of Ly-a photons escape fraction shows that galaxies with a low star formation rate have a f_esc of the order of unity, suggesting that, for those objects, Ly-a may be used to trace the star formation rate assuming a given conversion law. In galaxies forming stars intensely, the escape fraction spans the whole range from 0 to 1. The model is able to get a good match to the UV and Ly-a luminosity function (LF) data at 3 < z < 5. We find that we are in good agreement with both the bright Ly-a data and the faint population observed by Rauch et al. (2008) at z=3. Most of the Ly-a profiles of our LAEs are redshifted by the diffusion in the outflow which suppresses IGM absorption. The bulk of the observed Ly-a equivalent width (EW) distribution is recovered by our model, but we fail to obtain the very large values sometimes detected. Predictions for stellar masses and UV LFs of LAEs show a satisfactory agreement with observational estimates. The UV-brightest galaxies are found to show only low Ly-a EWs in our model, as it is reported by many observations of high redshift LAEs. We interpret this effect as the joint consequence of old stellar populations hosted by UV-bright galaxies, and high HI column densities that we predict for these objects, which quench preferentially resonant Ly-a photons via dust extinction.

(Abridged) We use deep VLT optical and near-IR spectroscopy and deep Spitzer/IRAC imaging to examine the properties of two of the most luminous Lyman alpha emitters at z=5.7. The continuum red-ward of the Lyman alpha line is clearly detected in both objects, thus facilitating a relatively accurate measurement (10--20% uncertainties) of the observed rest-frame equivalent widths, which are around 160 Angstroms for both objects. Through detailed modelling of the profile of the Lyman alpha line with a 3-D Monte-Carlo radiative transfer code, we estimate the intrinsic rest-frame equivalent width of Lyman alpha and find values that are around 300 Angstroms, which is at the upper end of the range allowed for very young, moderately metal-poor star-forming galaxies. However, the uncertainties are large and values as high as 700 Angstroms are permitted by the data. Both Lyman alpha emitters are detected at 3.6 micron in deep images taken with the Spitzer Space Telescope. We use these measurements, the measurement of the continuum red-ward of Lyman alpha and other photometry to constrain the spectral energy distributions of these very luminous Lyman alpha emitters and to compare them to three similar Lyman alpha emitters from the literature. The contribution from nebular emission is included in our models: excluding it results in significantly higher masses. Evidence for dust is found in all objects, and emission from nebular lines often make a dominant contribution to the rest frame 3.6 micron flux. Rich in emission lines, these objects are prime targets for the next generation of extremely large telescopes, JWST and ALMA.

High-redshift Lyman-alpha blobs are extended, luminous, but rare structures that appear to be associated with the highest peaks in the matter density of the Universe. Their energy output and morphology are similar to powerful radio galaxies, but the source of the luminosity is unclear. Some blobs are associated with ultraviolet or infrared bright galaxies, suggesting an extreme starburst event or accretion onto a central black hole. Another possibility is gas that is shock excited by supernovae. However some blobs are not associated with galaxies, and may instead be heated by gas falling into a dark matter halo. The polarization of the Ly-alpha emission can in principle distinguish between these options, but a previous attempt to detect this signature returned a null detection. Here we report on the detection of polarized Ly-alpha from the blob LAB1. Although the central region shows no measurable polarization, the polarized fraction (P) increases to ~20 per cent at a radius of 45 kpc, forming an almost complete polarized ring. The detection of polarized radiation is inconsistent with the in situ production of Ly-alpha photons, and we conclude that they must have been produced in the galaxies hosted within the nebula, and re-scattered by neutral hydrogen.

Supernova surveys can be used to study a variety of subjects such as: (i) cosmology through type Ia supernovae (SNe), (ii) star-formation rates through core-collapse SNe, and (iii) supernova properties and their connection to host galaxy characteristics. The Stockholm VIMOS Supernova Survey (SVISS) is a multi-band imaging survey aiming to detect supernovae at redshift ~0.5 and derive thermonuclear and core-collapse supernova rates at high redshift. In this paper we present the supernovae discovered in the survey along with light curves and a photometric classification into thermonuclear and core-collapse types. To detect the supernovae in the VLT/VIMOS multi-epoch images, we used difference imaging and a combination of automatic and manual source detection to minimise the number of spurious detections. Photometry for the found variable sources was obtained and careful simulations were made to estimate correct errors. The light curves were typed using a Bayesian probability method and Monte Carlo simulations were used to study misclassification. We detected 16 supernovae, nine of which had a core-collapse origin and seven had a thermonuclear origin. The estimated misclassification errors are quite small, in the order of 5%, but vary with both redshift and type. The mean redshift of the supernovae is 0.58. Additionally, we found a variable source with a very extended light curve that could possibly be a pair instability supernova.

Lya is a key diagnostic for numerous observations of distant star-forming galaxies. It's interpretation requires, however, detailed radiation transfer models. We provide an extensive grid of 3D radiation transfer models simulating the Lya and UV continuum radiation transfer in the interstellar medium of star-forming galaxies. We have improved our Monte Carlo MCLya code, and have used it to compute a grid of 6240 radiation transfer models for homogeneous spherical shells containing HI and dust surrounding a central source. The simulations cover a wide range of parameter space. We present the detailed predictions from our models including in particular the Lya escape fraction fesc, the continuum attenuation, and detailed Lya line profiles. The Lya escape fraction is shown to depend strongly on dust content, but also on other parameters (HI column density and radial velocity). The predicted line profiles show a great diversity of morphologies ranging from broad absorption lines to emission lines with complex features. The results from our simulations are distributed in electronic format. Our models should be of use for the interpretation of observations from distant galaxies, for other simulations, and should also serve as an important base for comparison for future, more refined, radiation transfer models.

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 Ly-alpha emission line has been proven a powerful tool by which to study evolving galaxies at the highest redshifts. However, in order to use Lya as a physical probe of galaxies, it becomes vital to know the Lya escape fraction (fescLya). Unfortunately, due to the resonant nature of Lya, fescLya may vary unpredictably and requires empirical measurement. Here we compile Lya luminosity functions between redshift z=0 and 8 and, combined with H-alpha and ultraviolet data, assess how fescLya evolves with redshift. We find a strong upwards evolution in fescLya over the range z=0.3-6, which is well-fit by the power-law fescLya \propto (1+z)^\xi, with \xi =(2.57_-0.12^+0.19). This predicts that fescLya should reach unity at z=11.1. By comparing fescLya and E(B-V) in individual galaxies we derive an empirical relationship between fescLya and E(B-V), which includes resonance scattering and can explain the redshift evolution of fescLya between z=0 and 6 purely as a function of the evolution in the dust content of galaxies. Beyond z~6.5, fescLya drops more substantially; an effect attributed to either ionizing photon leakage, or an increase in the neutral gas fraction of the intergalactic medium. While distinguishing between those two scenarios may be extremely challenging, by framing the problem this way we remove the uncertainty of the halo mass from Lya-based tests of reionization. We finally derive a new method by which to estimate the dust content of galaxies based purely upon the observed Lya and UV LFs. These data are characterized by an exponential with an e-folding redshift of ~3.5.

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.

Aims: We investigate the nature of the source GDS J033218.92-275302.7at redshift ~ 5.56. Methods: The SED of the source is well sampled by 16 bands photometry, from UV-optical, near infrared and mid-infrared (MID-IR).The detection of signal in the MID-IR Spitzer/IRAC bands 5.8, 8.0 um -- where the nebular emission contribution is less effective -- suggests the presence of a Balmer break, signature of an underlying stellar population formed at earlier epochs. The optical spectrum shows a clear Lya emission line together with semi-forbidden NIV] 1483.3-1486.5 also in emission. Results: From the SED fitting and the Lya modelling it turns out that the source seems to have an evolved component with stellar mass of ~5 x10^(10) Msolar and age ~ 0.4 Gyrs, and a young component with an age of ~ 0.01 Gyrs and SFR in the range of 30-200 Msolar yr^(-1). The limits on the effective radius derived from the ACS/z850 and VLT/Ks bands indicate that this galaxy is denser than the local ones with similar mass. A relatively high nebular gas column density is favored from the Lya line modelling (NHI>=10^(21) cm^(-2)). A vigorous outflow (~ 450 km/s) has been measured from the optical spectrum,consistent with the Lya modelling. From ACS observations it turns out that the region emitting Lya photons is spatially compact and of the same order of the effective radius estimated at the ~1400A rest-frame wavelength, whose emission is dominated by the stellar continuum and/or AGN. The gas is blown out from the central region,but given the mass of the galaxy it is uncertain whether it will pollute the IGM to large distances. We argue that a burst of star formation in a dense gas environment is active (possibly containing hot and massive stars and/or a low luminosity AGN), superimposed to an already formed fraction of stellar mass (abridged).

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)

The Lyman-alpha (Lya) emission line is the primary observational signature of star-forming galaxies at the highest redshifts, and has enabled the compilation of large samples of galaxies with which to study cosmic evolution. The resonant nature of the line, however, means that Lya photons scatter in the neutral interstellar medium of their host galaxies, and their sensitivity to absorption by interstellar dust may therefore be enhanced greatly. This implies that the Lya luminosity may be significantly reduced, or even completely suppressed. Hitherto, no unbiased empirical test of the escaping fraction (f_esc) of Lya photons has been performed at high redshifts. Here we report that the average fesc from star-forming galaxies at redshift z = 2.2 is just 5 per cent by performing a blind narrowband survey in Lya and Ha. This implies that numerous conclusions based on Lya-selected samples will require upwards revision by an order of magnitude and we provide a benchmark for this revision. We demonstrate that almost 90 per cent of star-forming galaxies emit insufficient Lya to be detected by standard selection criteria. Both samples show an anti-correlation of fesc with dust content, and we show that Lya- and Ha-selection recovers populations that differ substantially in dust content and fesc.

We aim to place new, strengthened constraints on the luminosity function (LF) of H-alpha emitting galaxies at redshift z=2.2, and to further constrain the instantaneous star-formation rate density of the universe (rho*). We have used the new HAWK-I instrument at ESO-VLT to obtain extremely deep narrow-band (line; NB2090) and broad-band (continuum; Ks) imaging observations. The target field is in the GOODS-South, providing us with a rich multi-wavelength auxiliary data set, which we utilise for redshift confirmation and to estimate dust content. We use this new data to measure the faint-end slope (alpha) of LF(H-alpha) with unprecedented precision. The data are well fit by a Schechter function and also a single power-law, yielding alpha=(-1.72 +/- 0.20) and (-1.77 +/- 0.21), respectively. Thus we are able to confirm the steepening of alpha from low- to high-z predicted by a number of authors and observed at other wavelengths. We combine our LF data-points with those from a much shallower but wider survey at z=2.2 (Geach et al. 2008), constructing a LF spanning a factor of 50 in luminosity. Re-fitting the Schechter parameters, we obtain log L*=(43.07+/-0.22)erg s^-1 ; log phi*=(-3.45+/-0.52)Mpc^-3 ; alpha=(-1.60+/-0.15). We integrate over LF(Halpha) and apply a correction for dust attenuation to determine the instantaneous cosmic star-formation rate density at z=2 without assuming alpha or extrapolating it from lower-z. Our measurement of rho* is (0.215+/-0.090) Msun yr^-1 Mpc^-3, integrated over a range of 37 <log(LHhalpha / erg s^-1) < 47.

We present the optical spectroscopic follow-up of 31 z=0.3 Lyman-alpha (Lya) emitters, previously identified by Deharveng et al. (2008). We find that 17% of the Lya emitters have line ratios that require the hard ionizing continuum produced by an AGN. The uniform dust screen geometry traditionally used in studies similar to ours is not able to simultaneously reproduce the observed high Lya/Halpha and Halpha/Hbeta line ratios. We consider different possibilities for the geometry of the dust around the emitting sources. We find that also a uniform mixture of sources and dust does not reproduce the observed line ratios. Instead, these are well reproduced by a clumpy dust screen. This more realistic treatment of the geometry results in extinction corrected (Lya/Halpha)_C values consistent with Case B recombination theory, whereas a uniform dust screen model would imply values (Lya/Halpha)_C higher than 8.7. Our analysis shows that there is no need to invoke "ad-hoc" multi phase media in which the Lya photons only scatter between the dusty clouds and eventually escape.

The Lyman-alpha (Lya) recombination line is a fundamental tool for galaxy evolution studies and modern observational cosmology. However, subsequent interpretations are still prone to a number of uncertainties. Besides numerical efforts, empirical data are urgently needed for a better understanding of Lya escape process. We empirically estimate the Lyman-alpha escape fraction fesc(Lya) in a statistically significant sample of z ~ 0 - 0.3 galaxies in order to calibrate high-redshift Lyman-alpha observations. An optical spectroscopic follow-up of a sub-sample of 24 Lyman-alpha emitters (LAEs) detected by GALEX at z ~ 0.2-0.3, combined with a UV-optical sample of local starbursts, both with matched apertures, allow us to quantify the dust extinction through Balmer lines, and to estimate the Lyman-alpha escape fraction from the Halpha flux corrected for extinction in the framework of the recombination theory. The global escape fraction of Lyman-alpha radiation spans nearly the entire range of values, from 0.5 to 100 %, and fesc(Lya) clearly decreases with increasing nebular dust extinction E(B-V). Several objects show fesc(Lya) greater than fesc(continuum) which may be an observational evidence for clumpy ISM geometry or for an aspherical ISM. Selection biases and aperture size effects may still prevail between z ~ 0.2-0.3 LAEs and local starbursts, which may explain the difference observed for fesc(Lya).

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.

[ABRIDGED] We are undertaking an imaging study of local star-forming galaxies in the Lyman-alpha (Lya) emission line using the Solar Blind Channel (SBC) of the ACS onboard the HST. Observations have been obtained in Lya and H-alpha (Ha) and 6 continuum filters between ~1500AA and the I-band. Previously (Hayes et al 2005) we demonstrated that the production of Lya line-only images in the SBC-only data-set is non-trivial and that supporting data is a requirement. We here develop various methods of continuum subtraction and assess their relative performance using a variety of spectral energy distributions (SED) as input. We conclude that simple assumptions about the behavior of the ultraviolet continuum consistently lead to results that are wildly erroneous, and determine that an SED fitting approach is essential. Moreover, a single component stellar or stellar+nebular spectrum is not always sufficient for realistic template SEDs and, in order to successfully recover the input observables, care must be taken to control contributions of nebular gas and underlying stellar populations. Independent measurements of the metallicity must first be obtained, but details of the initial mass function play only a small role. We identify the need to bin together pixels in our data to obtain signal-to-noise (S/N) of ~10 in each band before processing. At S/N=10 we can recover Lya fluxes accurate to ~30% for Lya lines with intrinsic equivalent width [W(Lya)] of 10AA. This accuracy improves to ~10% for W(Lya)=100AA. We describe the image processing applied to the observations presented in \"Ostlin et al. (2009) and the associated data-release. We also present simulations for an observing strategy for an alternative low-redshift SBC Lya imaging campaign using adjacent combinations of long-pass filters at slightly higher redshift.

We present reduced and calibrated high resolution Lyman-alpha (Lya) images for a sample of 6 local star forming galaxies. Targets were selected to represent a range in luminosity and metallicity and to include both known Lya emitters and non-emitters. Far ultraviolet imaging was carried out with the Solar Blind Channel of the ACS on HST in the F122M (Lya on-line) and F140LP (continuum) filters. The resulting Lya images are the product of careful modeling of both the stellar and nebular continua, facilitated by supporting HST imaging in Ha and 5 continuum bands, combined with Starburst99 evolutionary synthesis models, and prescriptions for dust extinction on the continuum. In all, the resulting morphologies in Lya, Ha, and UV-continuum are qualitatively very different and we show that the bulk of Lya emerges in a diffuse component resulting from resonant scattering events. Lya escape fractions, computed from integrated Ha luminosities and recombination theory, are found never to exceed 14%. Internal dust extinction is estimated in each pixel and used to correct Lya fluxes. However, the extinction corrections are far too small (factors from 2.6 to infinity) to reconcile the global Lya luminosities with recombination theory. Surprisingly, when comparing the global equivalent widths of Lya and Ha, the two quantities appear anti-correlated, which may be due to the evolution of mechanical feedback. This calls for caution in the interpretation of Lya observations. The images presented have a physical resolution 3 orders of magnitude better than attainable at high-z from the ground with current instrumentation and our images may therefore serve as useful templates for comparing with observations and modeling of primeval galaxy formation. We therefore provide the reduced Lya, Ha, and continuum images to the community.

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.

Lyman-alpha is now widely used to investigate the galaxy formation and evolution in the high redshift universe. However, without a rigorous understanding of the processes which regulate the Lya escape fraction, physical interpretations of high-z observations remain questionable. We examine six nearby star-forming galaxies to disentangle the role of the dust from other parameters such as gas kinematics, geometry and ISM morphology in the obscuration of Ly-alpha. Thereby we aim to understand the Ly-a escape physics and infer the implications for high-redshift studies. We use HST/ACS to produce continuum-subtracted Lya maps, and ground-based observations (ESO/NTT and NOT) to map the Halpha emission and the extinction E(B-V) in the gas phase derived from the Balmer decrement Halpha/Hbeta. When large outflows are present, the Lya emission appears not to correlate with the dust content, confirming the role of the HI kinematics in the escape of Lya photons. In the case of a dense, static HI covering, we observe a damped absorption with a declining relationship between Lya and E(B-V). We found that the Lya escape fraction does not exceed 10% in all our galaxies and is mostly about 3% or below. Finally, because of the radiative transfer complexity of the Lya line, star formation rate based on Lya luminosity is underestimated with respect to that derived from UV luminosity. The failure of simple dust correction to recover the intrinsic Lya/Ha ratio or the total star formation rate should prompt us to be more cautious when interpreting high-z observations and related properties, such as SFRs based on Lya alone. To this end we propose a more realistic calibration for SFR(Lya) which accounts for dust attenuation and resonant scattering effects via the Lya escape fraction.

The long-duration gamma-ray burst (GRB) 030329, associated with supernova (SN) 2003dh, occurred inside a star-forming dwarf galaxy at redshift $z=0.1685$. The low redshift, and a rich set of archival Hubble Space Telescope (HST) images, makes this GRB well-suited for a detailed study of the stellar population in the immediate vicinity of the explosion. Since the lifetime of a star is directly tied to its mass, the age of the stellar population can be used to put constraints on the GRB and SN progenitor mass. From the HST images we extract the colours of the precise site from which the GRB originated, and find that the colours are significantly different from those of the overall host galaxy and the surrounding starburst environment. We have used spectral evolutionary models, including nebular emission, to carefully constrain the age of the stellar population, and hence the progenitor, at the very explosion site. For instantaneous burst models we find that a population age of 5 Myr best matches the data, suggesting a very massive (M > 50 M_sun) star as the progenitor, with an upper limit of 8 Myr (M > 25 M_sun). For more extended star formation scenarios, the inferred progenitor age is in most cases still very young (age < 7 Myr, implying M > 25 M_sun), with an upper limit of 20 Myr (M > 12 M_sun). These age estimates are an order of magnitude lower than the ages inferred from the overall host galaxy colours, indicating that progenitor mass estimates based on data for spatially unresolved GRB host galaxies will in general be of limited use. Our results are consistent with the collapsar scenario.

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

Numerous surveys are underway or planned that aim to exploit the Lyman-alpha emission line as a probe of the distant universe. We investigate the reliability of such high-z Lya studies and the validity of the conclusions that are based upon them. We examine whether reliable Lya fluxes (F_Lya) and equivalent widths (W_Lya) can be estimated from narrowband imaging surveys and whether any observational biases may be present. We have developed software to simulate the observed line and continuum properties of synthetic LAEs in the distant universe by adopting various typical observational techniques. This was used to investigate how detected F_Lya and W_Lya vary with properties of the host or IGM: internal dust reddening; intervening Lya absorption systems; the presence of underlying stellar populations. None of the techniques studied are greatly susceptible to underlying stellar populations or the contribution of nebular gas. We find that techniques that use one off-line filter on the red side of Lya result in highly inaccurate measurements of W_Lya. Adopting two off-line filters to estimate continuum at Lya is an improvement but is still unreliable when dust extinction is considered. Techniques that employ single narrow- and broad-band filters with the same central wavelength are not susceptible to internal dust, but IGM Lya absorption can cause W_Lya to be overestimated by factors of up to 2: at z=6, the median W_Lya is overestimated by ~25%. The most robust approach is a SED fitting technique that fits E(B-V) and burst-age from synthetic models - broadband observations are needed that sample the UV continuum slope, 2175AA dust feature, and the 4000AA discontinuity. We also notice a redshift-dependent incompleteness that results from DLA systems close to the target LAEs, amounting to ~10% at z=6.

We present a satellite mission concept to measure the dark energy equation of state parameter w with percent-level precision. The Very Ambitious Dark Energy Research satellite (VADER) is a multi-wavelength survey mission joining X-ray, optical, and IR instruments for a simultaneous spectral coverage from 4microns (0.3eV) to 10keV over a field of view (FoV) of 1 square degree. VADER combines several clean methods for dark energy studies, the baryonic acoustic oscillations in the galaxy and galaxy cluster power spectrum and weak lensing, for a joint analysis over an unrivalled survey volume. The payload consists of two XMM-like X-ray telescopes with an effective area of 2,800cm^2 at 1.5keV and state-of-the-art wide field DEPFET pixel detectors (0.1-10keV) in a curved focal plane configuration to extend the FoV. The X-ray telescopes are complemented by a 1.5m optical/IR telescope with 8 instruments for simultaneous coverage of the same FoV from 0.3 to 4 microns. The 8 dichroic-separated bands (u,g,r,z,J,H,K,L) provide accurate photometric galaxy redshifts, whereas the diffraction-limited resolution of the central z-band allows precise shape measurements for cosmic shear analysis. The 5 year VADER survey will cover a contiguous sky area of 3,500 square degrees to a depth of z~2 and will yield accurate photometric redshifts and multi-wavelength object parameters for about 175,000 galaxy clusters, one billion galaxies, and 5 million AGN. VADER will not only provide unprecedented constraints on the nature of dark energy, but will additionally extend and trigger a multitude of cosmic evolution studies to very large (>10 Gyrs) look-back times.

We review the history and the current status of the understanding of the processes that regulate Lyman alpha emission from star-forming galaxies. We present some of the most recent results of our study to image local starburst galaxies in the Lyman alpha emission line using the Advanced Camera for Surveys on the Hubble Space Telescope. Particular attention is dedicated to our study of the low-metallicity, dust-poor Blue Compact Galaxy ESO338-IG04. We discuss some of our local observational results with reference to the interpretation of results of high-redshift Lyman alpha surveys and recent simulations of the detection properties of high-z Lyman alpha emitting objects performed by our group.

ESO338-IG04 is a well known, local, luminous Blue Compact Galaxy. It has a complex morphology and contains a central young starburst with compact young star clusters. The galaxy was imaged using HST/ACS in the Lyman alpha line and continuum. Using the Starburst99 synthetic spectra and other imaging data from HST/WFPC2 we developed a technique that allows us to make the first photometrically valid continuum subtraction for the LYa line. Our results are in qualitative agreement with the models of LYa escape being regulated by kinematical properties of the interstellar medium. The line-only image shows LYa in both emission and absorption. Notably, emission is seen from central bright young clusters and is in spatial agreement with the features present in a longslit spectrum taken with HST/STIS. Diffuse LYa emission regions are seen surrounding the central starburst where photons escape after one or more resonant scatterings in HI. These regions dominate the total LYa output. Photometry reveals a total LYa flux of 194 x 10^-14 erg/s/cm^2. Inside an artificial aperture matching that of the IUE satelite, we measure a LYa flux and equivalent width consistent with previous studies using the IUE. Analysis of parameter dependencies show our technique to be largely parameter independent, producing LYa maps indistinguishable from one another by eye and with fluxes consistent with one another to better than 50%. By comparison of LYa fluxes with Ha fluxes of a previous study, we estimate that each observed LYa photon has undergone ~2 additional scatterings. We see that LYa line closely correlate with other tracers of kinematics.

We review the imaging and spectroscopic properties of the Lyman alpha emission in starburst galaxies. The Lyman alpha photons escape is largely driven by kinematical and orientation effects related to the presence of large scale expanding shells of neutral hydrogen. The various Lyman alpha profiles can be linked to the evolutionary state of the starburst. Our recent Lyman alpha imaging study using the HST-ACS is presented here: it highlights the presence of diffuse extended emission of low intensity leaking out the diffuse HI components. The similarities with the line profiles observed on high redshift galaxies is an important fact that cautions the use of Lyman alpha to derive the cosmic star formation rate and constrains the search for high redshift galaxies.