We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$\alpha$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$\alpha$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$\alpha$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.

Identifying the mechanisms driving the escape of Lyman Continuum (LyC) photons is crucial to find Lyman Continuum Emitter (LCE) candidates. To understand the physical properties involved in the leakage of LyC photons, we investigate the connection between the HI covering fraction, HI velocity width, the Lyman alpha (LyA) properties and escape of LyC photons in a sample of 22 star-forming galaxies including 13 LCEs. We fit the stellar continua, dust attenuation, and absorption lines between 920 and 1300 A to extract the HI covering fractions and dust attenuation. Additionally, we measure the HI velocity widths of the optically thick Lyman series and derive the LyA equivalent widths (EW), escape fractions (fesc), peak velocities and fluxes at the minimum of the LyA profiles. Overall, we highlight strong correlations between the presence of low HI covering fractions and (1) low LyA peak velocities; (2) more flux at the profile minimum; and (3) larger EW(LyA), fesc(LyA), and fesc(LyC). Hence, low column density channels are crucial ISM ingredients for the leakage of LyC and LyA photons. Additionally, galaxies with narrower HI absorption velocity widths have higher LyA equivalent widths, larger LyA escape fractions, and lower LyA peak velocity separations. This suggests that these galaxies have low HI column density. Finally, we find that dust regulates the amount of LyA and LyC radiation that actually escapes the ISM. Overall, the ISM porosity is one origin of strong LyA emission and enables the escape of ionizing photons in low-z leakers. However, this is not enough to explain the largest fesc(LyC) observed, which indicates that the most extreme LCEs are likely density-bounded along all lines of sight to the observer. Overall, the neutral gas porosity constrains a lower limit to the escape fraction of LyC and LyA photons, providing a key estimator of the leakage of ionizing photons.

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

We present spatially resolved maps of six individually-detected Lyman alpha haloes (LAHs) as well as a first statistical analysis of the Lyman alpha (Lya) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies using MUSE. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lya line profile. Using a three-dimensional two-component model for the Lya emission, we measure the full width at half maximum (FWHM), the peak velocity shift and the asymmetry of the Lya line in the core and in the halo of 19 galaxies. We find that the Lya line shape is statistically different in the halo compared to the core for ~40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lya line both at the interstellar and circum-galactic scales. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies show broader, more redshifted and less asymmetric Lya lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lya lines. The generally broad and red line shapes found in the halo component suggests that the Lya haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lya line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.

Resonant lines are powerful probes of the interstellar and circumgalactic medium of galaxies. Their transfer in gas being a complex process, the interpretation of their observational signatures, either in absorption or in emission, is often not straightforward. Numerical radiative transfer simulations are needed to accurately describe the travel of resonant line photons in real and in frequency space, and to produce realistic mock observations. This paper introduces RASCAS, a new public 3D radiative transfer code developed to perform the propagation of any resonant line in numerical simulations of astrophysical objects. RASCAS was designed to be easily customisable and to process simulations of arbitrarily large sizes on large supercomputers. RASCAS performs radiative transfer on an adaptive mesh with an octree structure using the Monte Carlo technique. RASCAS features full MPI parallelisation, domain decomposition, adaptive load-balancing, and a standard peeling algorithm to construct mock observations. The radiative transport of resonant line photons through different mixes of species (e.g. \ion{H}{i}, \ion{Si}{ii}, \ion{Mg}{ii}, \ion{Fe}{ii}), including their interaction with dust, is implemented in a modular fashion to allow new transitions to be easily added to the code. RASCAS is very accurate and efficient. It shows perfect scaling up to a minimum of a thousand cores. It has been fully tested against radiative transfer problems with analytic solutions and against various test cases proposed in the literature. Although it was designed to describe accurately the many scatterings of line photons, RASCAS may also be used to propagate photons at any wavelength (e.g. stellar continuum or fluorescent lines), or to cast millions of rays to integrate the optical depths of ionising photons, making it highly versatile.

Galactic outflows are thought to eject baryons back out to the circum-galactic medium (CGM). Studies based on metal absorption lines (MgII in particular) in the spectra of background quasars indicate that the gas is ejected anisotropically, with galactic winds likely leaving the host in a bi-conical flow perpendicular to the galaxy disk. In this paper, we present a detailed analysis of an outflow from a z = 0.7 "green-valley" galaxy (log($M_*$/$\mathrm{M}_\odot$) = 9.9; SFR = 0.5 $\mathrm{M}_\odot\,\mathrm{yr}^{-1}$) probed by two background sources part of the MUSE Gas Flow and Wind (MEGAFLOW) survey. Thanks to a fortuitous configuration with a background quasar (SDSSJ1358+1145) and a bright background galaxy at $z = 1.4$, both at impact parameters of $\approx 15\,\mathrm{kpc}$, we can - for the first time - probe both the receding and approaching components of a putative galactic outflow around a distant galaxy. We measure a significant velocity shift between the MgII absorption from the two sightlines ($84\pm17\,\mathrm{km}\,\mathrm{s}^{-1}$), which is consistent with the expectation from our simple fiducial wind model, possibly combined with an extended disk contribution.

We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of eight compact star-forming galaxies at redshifts z=0.02811-0.06540, with low oxygen abundances 12+log(O/H)=7.43-7.82 and extremely high emission-line flux ratios O32=[OIII]5007/[OII]3727~22-39, aiming to study the properties of Ly-alpha emission in such conditions. We find a diversity in Ly-alpha properties. In five galaxies Ly-alpha emission line is strong, with equivalent width (EW) in the range 45-190A. In the remaining galaxies, weak Ly-alpha emission with EW(Ly-alpha)~2-7A is superposed on a broad Ly-alpha absorption line, indicating a high neutral hydrogen column density N(HI)~(1-3)x10^21 cm^-2. We examine the relation between the Ly-alpha escape fraction fesc(Ly-alpha) and the Lyman continuum escape fraction fesc(LyC), using direct measures of the latter in eleven low-redshift LyC leakers, to verify whether fesc(Ly-alpha) can be an indirect measure of escaping LyC radiation. The usefulness of O32, of the Ly-alpha equivalent width EW(Ly-alpha) and of the Ly-alpha peak separation Vsep as indirect indicators of Ly-alpha leakage is also discussed. It is shown that there is no correlation between O32 and fesc(Ly-alpha). We find an increase of fesc(Ly-alpha) with increasing EW(Ly-alpha) for EW(Ly-alpha)<100A, but for higher EW(Ly-alpha)>150A the fesc(Ly-alpha) is nearly constant attaining the value of ~0.25. We find an anticorrelation between fesc(Ly-alpha) and Vsep, though not as tight as the one found earlier between fesc(LyC) and Vsep. This finding makes Vsep a promising indirect indicator of both the Ly-alpha and ionizing radiation leakage.

Lyman-\alpha (Ly\alpha) emission lines are typically found to be redshifted with respect to the systemic redshifts of galaxies, likely due to resonant scattering of Ly\alpha photons. Here we measure the average velocity offset for a sample of 96 z ~ 3.3 Ly\alpha emitters (LAEs) with a median Ly\alpha flux (luminosity) of ~10^{-17} erg cm^{-2} s^{-1} (~10^{42} erg s^{-1}) and a median unobscured star formation rate (SFR) of ~1.3 Msun/yr, detected by the Multi-Unit Spectroscopic Explorer as part of our MUSEQuBES circumgalactic medium (CGM) survey. By postulating that the stacked CGM absorption profiles of these LAEs, probed by 8 background quasars, must be centered on the systemic redshift, we measure an average velocity offset, Voffset = 169 \pm 10 km/s, between the Ly\alpha emission peak and the systemic redshift. The observed Voffset is lower by factors of ~1.4 and ~2.6 compared to the velocity offsets measured for narrow-band selected LAEs and Lyman break galaxies, respectively, which probe galaxies with higher masses and SFRs. Consistent with earlier studies based on direct measurements for individual objects, we find that the Voffset is correlated with the full width at half-maximum of the red peak of the Ly\alpha line, and anti-correlated with the rest-frame equivalent width. Moreover, we find that Voffset is correlated with SFR with a sub-linear scaling relation, Voffset ~ SFR^(0.16 \pm 0.03). Adopting the mass scaling for main sequence galaxies, such a relation suggests that Voffset scales with the circular velocity of the dark matter halos hosting the LAEs.

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 present an analysis of HI Lyman-alpha emission in deep VLT/MUSE observations of two highly magnified and extended galaxies at z=3.5 and 4.03, including a newly discovered, almost complete Einstein ring. While these Lyman-alpha haloes are intrinsically similar to the ones typically seen in other MUSE deep fields, the benefits of gravitational lensing allows us to construct exceptionally detailed maps of Lyman-alpha line properties at sub-kpc scales. By combining all multiple images, we are able to observe complex structures in the Lyman-alpha emission and uncover small (~ 120 km/s in Lyman-alpha peak shift), but significant at > 4 sigma, systematic variations in the shape of the Lyman-alpha line profile within each halo. Indeed, we observe a global trend for the line peak shift to become redder at large radii, together with a strong correlation between the peak wavelength and line width. This systematic intrahalo variation is markedly similar to the object-to-object variations obtained from the integrated properties of recent large samples. Regions of high surface brightness correspond to relatively small line shifts, which could indicate that Lyman-alpha emission escapes preferentially from regions where the line profile has been less severely affected by scattering of Lyman-alpha photons.

Understanding the escape of Lyman continuum (LyC) and Lyman alpha (Lya) photons from molecular clouds is one of the keys to constraining the reionization history of the Universe and the evolution of galaxies at high redshift. Using a set of radiation-hydrodynamic simulations with adaptive mesh refinement, we investigate how photons propagate and escape from turbulent clouds with different masses, star formation efficiencies (SFEs), and metallicities, as well as with different models of stellar spectra and supernova feedback. We find that the escape fractions in both LyC and Lya are generally increasing with time if the cloud is efficiently dispersed by radiation and supernova feedback. When the total SFE is low (1% of the cloud mass), 0.1-5% of LyC photons leave the metal-poor cloud, whereas the fractions increase to 20-70% in clouds with a 10% SFE. LyC photons escape more efficiently if gas metallicity is lower, if the upper mass limit in the stellar initial mass function is higher, if binary interactions are allowed in the evolution of stars, or if additional strong radiation pressure, such as Lya pressure, is present. As a result, the number of escaping LyC photons can easily vary by a factor of $\sim4$ on cloud scales. The escape fractions of Lya photons are systemically higher (60-80%) than those of LyC photons despite large optical depths at line centre ($\tau_0\sim10^6-10^9$). Scattering of Lya photons is already significant on cloud scales, leading to double-peaked profiles with peak separations of $v_{\rm sep}\sim400\,{\rm km\,s^{-1}}$ during the initial stage of the cloud evolution, while it becomes narrower than $v_{\rm sep} \le 150 \, {\rm km\,s^{-1}}$ in the LyC bright phase. Comparisons with observations of low-redshift galaxies suggest that Lya photons require further interactions with neutral hydrogen to reproduce their velocity offset for a given LyC escape fraction.

Deep optical spectroscopic surveys of galaxies provide us a unique opportunity to investigate rest-frame ultra-violet (UV) emission line properties of galaxies at ${z \sim 2-4.5}$. Here we combine VLT/MUSE Guaranteed Time Observations of the Hubble Deep Field South, Ultra Deep Field, COSMOS, and several quasar fields with other publicly available data from VLT/VIMOS and VLT/FORS2 to construct a catalogue of HeII${\lambda}$1640 emitters at ${z\sim2}$. The deepest areas of our MUSE pointings reach a ${3\sigma}$ line flux limit of 3.1${\times}$10-19 erg/ s/ cm$^2$. After discarding broad line active galactic nuclei we find 13 HeII${\lambda}$1640 detections from MUSE with a median MUV = $-20.1$ and 21 tentative HeII${\lambda}$1640 detections from other public surveys. Excluding Ly${\alpha}$, all except two galaxies in our sample show at least one other rest-UV emission line, with CIII]${\lambda}$1907,${\lambda}$1909 being the most prominent. We use multi-wavelength data available in the Hubble legacy fields to derive basic galaxy properties of our sample via spectral energy distribution fitting techniques. Taking advantage of the high quality spectra obtained by MUSE (${\sim10 - 30}$h of exposure time per pointing), we use photo-ionisation models to study the rest-UV emission line diagnostics of the HeII${\lambda}$1640 emitters. Line ratios of our sample can be reproduced by moderately sub-solar photo-ionisation models, however, we find that including effects of binary stars lead to degeneracies in most free parameters. Even after considering extra ionising photons produced by extreme sub-solar metallicity binary stellar models, photo-ionisation models are unable to reproduce rest-frame HeII${\lambda}$1640 equivalent widths (${\sim}$ 0.2 - 10 A), thus additional mechanisms are necessary in models to match the observed HeII${\lambda}$1640 properties.

Using an ultra-deep blind survey with the MUSE integral field spectrograph on the ESO Very Large Telescope, we obtain spectroscopic redshifts to a depth never explored before: galaxies with observed magnitudes m > 30 - 32. Specifically, we detect objects via Lyman-alpha emission at 2.9 < z < 6.7 without individual continuum counterparts in areas covered by the deepest optical/near-infrared imaging taken by the Hubble Space Telescope, the Hubble Ultra Deep Field. In total, we find 102 such objects in 9 square arcminutes at these redshifts. Detailed stacking analyses confirm the Lyman-alpha emission as well as the 1216 Angstrom-breaks and faint UV continua (M_UV ~ -15). This makes them the faintest spectroscopically-confirmed objects at these redshifts, similar to the sources believed to reionize the universe. A simple model for the expected fraction of detected/undetected Lyman-alpha emitters as a function of luminosity is consistent with these objects being the high-equivalent width tail of the normal Lyman-alpha-emitter population at these redshifts.

We present a study of the [OIII]5007/[OII]3727 (O32) ratios of star-forming galaxies drawn from MUSE data spanning a redshift range 0.28<z<0.85. Recently discovered Lyman continuum (LyC) emitters have extremely high oxygen line ratios: O32>4. Here we aim to understand the properties and the occurrences of galaxies with such high line ratios. Combining data from several MUSE GTO programmes, we select a population of star-forming galaxies with bright emission lines, from which we draw 406 galaxies for our analysis based on their position in the z-dependent star formation rate (SFR) - stellar mass (M*) plane. Out of this sample 15 are identified as extreme oxygen emitters based on their O32 ratios (3.7%) and 104 galaxies have O32>1 (26%). Our analysis shows no significant correlation between M*, SFR, and the distance from the SFR-M* relation with O32. We find a decrease in the fraction of galaxies with O32>1 with increasing M*, however, this is most likely a result of the relationship between O32 and metallicity, rather than between O32 and M*. We draw a comparison sample of local analogues with <z>~0.03 from SDSS, and find similar incidence rates for this sample. In order to investigate the evolution in the fraction of high O32 emitters with redshift, we bin the sample into three redshift subsamples of equal number, but find no evidence for a dependence on redshift. Furthermore, we compare the observed line ratios with those predicted by nebular models with no LyC escape and find that most of the extreme oxygen emitters can be reproduced by low metallicity models. The remaining galaxies are likely LyC emitter candidates. Finally, based on a comparison between electron temperature estimates from the [OIII4363]/[OIII]5007 ratio of the extreme oxygen emitters and nebular models, we argue that the galaxies with the most extreme O32 ratios have young light-weighted ages.

We have obtained the first complete ultraviolet (UV) spectrum of a strong Lyman continuum(LyC) emitter at low redshift -- the compact, low-metallicity, star-forming galaxy J1154+2443 -- with a Lyman continuum escape fraction of 46% discovered recently. The Space Telescope Imaging Spectrograph spectrum shows strong Lya and CIII] 1909 emission, as well as OIII] 1666. Our observations show that strong LyC emitters can have UV emission lines with a high equivalent width (e.g. EW(CIII])$=11.7 \pm 2.9 \AA$ rest-frame), although their equivalent widths should be reduced due to the loss of ionizing photons. The intrinsic ionizing photon production efficiency of J1154+2443 is high, $\log(\xi_{\rm ion}^0)=25.56$ erg$^{-1}$ Hz, comparable to that of other recently discovered $z \sim 0.3-0.4$ LyC emitters. Combining our measurements and earlier determinations from the literature, we find a trend of increasing $\xi_{\rm ion}^0$ with increasing CIII] 1909 equivalent width, which can be understood by a combination of decreasing stellar population age and metallicity. Simple ionization and density-bounded photoionization models can explain the main observational features including the UV spectrum of J1154+2443.

The physical origin of the near-ultraviolet MgII emission remains an under-explored domain, contrary to more typical emission lines detected in the spectra of star-forming galaxies. We explore the nebular and physical properties for a sample of 381 galaxies between 0.70 < z < 2.34 drawn from the MUSE Hubble Ultra Deep Survey. The spectra of these galaxies show a wide variety of profiles of the MgII 2796,2803 resonant doublet, from absorption to emission. We present a study on the main drivers for the detection of MgII emission in galaxy spectra. By exploiting photoionization models we verified that the emission-line ratios observed in galaxies with MgII in emission are consistent with nebular emission from HII regions. From a simultaneous analysis of MUSE spectra and ancillary HST information via spectral energy distribution (SED) fitting, we find that galaxies with MgII in emission have lower stellar masses, smaller sizes, bluer spectral slopes and lower optical depth than those with absorption. This leads us to suggest that MgII emission is a potential tracer of physical conditions not merely related to those of the ionized gas. We show that these differences in MgII emission/absorption can be explained in terms of a higher dust and neutral gas content in the interstellar medium (ISM) of galaxies showing MgII in absorption, confirming the extreme sensitivity of MgII to the presence of the neutral ISM. We conclude with an analogy between the MgII doublet and the Ly-alpha line, due to their resonant nature. Further investigations with current and future facilities, including JWST, are promising as the detection of MgII emission and its potential connection with Ly-alpha could provide new insights on the ISM content in the early Universe.

Recent observations have discovered an escape of ionizing Lyman-continuum radiation from a population of compact, actively star-forming galaxies at redshift z~0.3, also known as "green peas". We here study the peculiar green pea Lyman-alpha (Lya) spectra, which are mostly double-peaked, unlike in any other galaxy sample. We select a sample of twelve archival green peas and we apply numerical radiative transfer models to reproduce the observed Lya spectral profiles, using the geometry of expanding, homogeneous spherical shells. We use ancillary optical and ultraviolet data to constrain the model parameters, and we evaluate the match between the models and the observed Lya spectra. As a second step, we allow all the fitting parameters to be free, and examine the agreement between the interstellar medium parameters derived from the models and those from ancillary data. The green pea double-peaked Lya line profiles are not correctly reproduced by the constrained shell models. Conversely, unconstrained models fit the spectra, but parameters derived from the best-fitting models are not in agreement with the ancillary data. In particular: 1) the best-fit systemic redshifts are larger by 10 - 250 km/s than those derived from optical emission lines, 2) the double-peaked Lya profiles are best reproduced with low-velocity (<150 km/s) outflows that contradict the observed ultraviolet absorption lines of low-ionization-state elements with characteristic velocities as large as 300 km/s, and 3) the models need to consider intrinsic Lya profiles that are on average three times broader than the observed Balmer lines. Differences between the modelled and observed velocities are larger for targets with prominent Lya blue peaks. The blue peak position and flux appear to be connected to low column densities of neutral hydrogen, leading to Lya and Lyman-continuum escape.

We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of five star-forming galaxies at redshifts z in the range 0.2993-0.4317 and with high emission-line flux ratios O32=[OIII]5007/[OII]3727 ~ 8-27 aiming to detect the Lyman continuum (LyC) emission. We detect LyC emission in all galaxies with the escape fractions fesc(LyC) in a range of 2-72 per cent. A narrow Ly-alpha emission line with two peaks in four galaxies and with three peaks in one object is seen in medium-resolution COS spectra with a velocity separation between the peaks Vsep varying from ~153 km/s to ~345 km/s. We find a general increase of the LyC escape fraction with increasing O32 and decreasing stellar mass M*, but with a large scatter of fesc(LyC). A tight anti-correlation is found between fesc(LyC) and Vsep making Vsep a good parameter for the indirect determination of the LyC escape fraction. We argue that one possible source driving the escape of ionizing radiation is stellar winds and radiation from hot massive stars.

The processes allowing the escape of ionizing photons from galaxies into the intergalactic medium are poorly known. To understand how Lyman continuum (LyC) photons escape galaxies, we constrain the HI covering fractions and column densities using ultraviolet HI and metal absorption lines of 18 star-forming galaxies which have Lyman series observations. Nine of these galaxies are confirmed LyC emitters. We fit the stellar continuum, dust attenuation, metal, and HI properties to consistently determine the UV attenuation, as well as the column densities and covering factors of neutral hydrogen and metals. We use synthetic interstellar absorption lines to explore the systematics of our measurements. Then we apply our method to the observed UV spectra of low-redshift and z-2 galaxies. The observed HI lines are found to be saturated in all galaxies. An indirect approach using OI column densities and the observed O/H abundances yields HI column densities of 18.6 to 20 cm-2. These columns are too high to allow the escape of ionizing photons. We find that the known LyC leakers have HI covering fractions less than unity. Ionizing photons escape through optically thin holes/channels in a clumpy interstellar medium. Our simulations confirm that the HI covering fractions are accurately recovered. The SiII and HI covering fractions scale linearly, in agreement with observations from stacked Lyman break galaxy spectra at z-3. Thus, with an empirical correction, the SiII absorption lines can also be used to determine the HI coverage. Finally, we show that a consistent fitting of dust attenuation, continuum and absorption lines is required to properly infer the covering fraction of neutral gas and subsequently to infer the escape fraction of ionizing radiation. These measurements can estimate the LyC escape fraction, as we demonstrate in a companion paper.

The Lyman alpha (lya) line of Hydrogen is a prominent feature in the spectra of star-forming galaxies, usually redshifted by a few hundreds of km/s compared to the systemic redshift. This large offset hampers follow-up surveys, galaxy pair statistics and correlations with quasar absorption lines when only lya is available. We propose diagnostics that can be used to recover the systemic redshift directly from the properties of the lya line profile. We use spectroscopic observations of Lyman-Alpha Emitters (LAEs) for which a precise measurement of the systemic redshift is available. Our sample contains 13 sources detected between z~3 and z~6 as part of various Multi Unit Spectroscopic Explorer (MUSE) Guaranteed Time Observations (GTO). We also include a compilation of spectroscopic lya data from the literature spanning a wide redshift range (z~0-8). First, restricting our analysis to double-peaked lya spectra, we find a tight correlation between the velocity offset of the red peak with respect to the systemic redshift, Vpeak, and the separation of the peaks. Secondly, we find a correlation between Vpeak and the full width at half maximum of the lya line. Fitting formulas, to estimate systemic redshifts of galaxies with an accuracy of +-100 km/s when only the lya emission line is available, are given for the two methods.

The fraction of ionizing photons that escape high-redshift galaxies sensitively determines whether galaxies reionized the early universe. However, this escape fraction cannot be measured from high-redshift galaxies because the opacity of the intergalactic medium is large at high redshifts. Without methods to indirectly measure the escape fraction of high-redshift galaxies, it is unlikely that we will know what reionized the universe. Here, we analyze the far-ultraviolet (UV) H I (Lyman series) and low-ionization metal absorption lines of nine low-redshift, confirmed Lyman continuum emitting galaxies. We use the H I covering fractions, column densities, and dust attenuations measured in a companion paper to predict the escape fraction of ionizing photons. We find good agreement between the predicted and observed Lyman continuum escape fractions (within $1.4\sigma$) using both the H I and ISM absorption lines. The ionizing photons escape through holes in the H I, but we show that dust attenuation reduces the fraction of photons that escape galaxies. This means that the average high-redshift galaxy likely emits more ionizing photons than low-redshift galaxies. Two other indirect methods accurately predict the escape fractions: the Ly$\alpha$ escape fraction and the optical [O III]/[O II] flux ratio. We use these indirect methods to predict the escape fraction of a sample of 21 galaxies with rest-frame UV spectra but without Lyman continuum observations. Many of these galaxies have low escape fractions ($f_{\rm esc} \le 1$\%), but 11 have escape fractions $>1$\%. The methods presented here will measure the escape fractions of high-redshift galaxies, enabling future telescopes to determine whether star-forming galaxies reionized the early universe.

The Mg II 2796,2803 doublet is often used to measure interstellar medium absorption in galaxies, thereby serving as a diagnostic for feedback and outflows. However, the interpretation of Mg II remains confusing, due to resonant trapping and re-emission of the photons, analogous to Lyman Alpha. Therefore, in this paper, we present new MMT Blue Channel Spectrograph observations of Mg II for a sample of 10 Green Pea galaxies at z~0.2-0.3, where Lyman Alpha was previously observed with the Cosmic Origins Spectrograph on Hubble Space Telescope. With strong, (mostly) double-peaked Lyman Alpha profiles, these galaxies allow us to observe Mg II in the limit of low H I column density. We find strong Mg II emission and little-to-no absorption. We use photoionization models to show that nebular Mg II from H II regions is non-negligible, and the ratios of Mg II/[OIII] 5007 vs. [OIII] 5007/[OII] 3727 form a tight sequence. Using this relation, we predict intrinsic Mg II flux, and show that Mg II escape fractions range from 0 to 0.9. We find that the Mg II escape fraction correlates tightly with the Lyman Alpha escape fraction, and the Mg II line profiles show evidence for broader and more redshifted emission when the escape fractions are low. These trends are expected if the escape fractions and velocity profiles of Lyman Alpha and Mg II are shaped by resonant scattering in the same low column density gas. As a consequence of a close relation with Lyman Alpha, Mg II may serve as a useful diagnostic in the epoch of reionization, where Lyman Alpha and Lyman continuum photons are not easily observed.

To identify the galaxies responsible for the reionization of the Universe, we must rely on the investigation of the Lyman Continuum (LyC) properties of z<5 star-forming galaxies, where we can still directly observe their ionizing radiation. We selected a sample of 201 star-forming galaxies from the Vimos Ultra Deep Survey at 3.5<z<4.3 to explore the validity of some of the proposed indirect indicators of LyC radiation. We created subsamples of galaxies with EWLya>70{\AA}, Lya(ext)<5.7kpc, rUV<0.3kpc and |Dv Lya|<200km/s, stacked all the galaxies in each subsample and measured the flux density ratio fnu(895)/fnu(1470), that we consider to be a proxy for LyC emission. We compared these ratios to those obtained for the complementary samples. We find that the stacks of galaxies which are UV compact (rUV<0.3kpc) and have bright Lya emission (EWLya>70{\AA}), have much higher LyC fluxes compared to the rest of the galaxy population in agreement with theoretical studies and previous observational works. We also find that galaxies with low Lya(ext) have the highest LyC flux: this new correlation seems even stronger than the correlations with high EWLya and small rUV. These results assume that the stacks from all the subsamples present the same statistical contamination from lower redshift interlopers. If we subtract a statistical contamination from low redshift interlopers obtained with dedicated Monte Carlo simulations, from the flux density ratios (fnu(895)/fnu(1470)) of the significant subsamples we find that these samples contain real LyC leaking flux with a very high probability, but the true average escape fractions remain uncertain. Our work indicates that galaxies with very high EWLya, small Lya(ext) and small rUV are very likely the best candidates to show LyC radiation at z=4 and could therefore be the galaxies that contributed more to reionization.

The "Antennae Galaxy" (NGC 4038/39) is the closest major interacting galaxy system and therefore often taken as merger prototype. We present the first comprehensive integral field spectroscopic dataset of this system, observed with the MUSE instrument at the ESO VLT. We cover the two regions in this system which exhibit recent star-formation: the central galaxy interaction and a region near the tip of the southern tidal tail. In these fields, we detect HII regions and diffuse ionized gas to unprecedented depth. About 15% of the ionized gas was undetected by previous observing campaigns. This newly detected faint ionized gas is visible everywhere around the central merger, and shows filamentary structure. We estimate diffuse gas fractions of about 60% in the central field and 10% in the southern region. We are able to show that the southern region contains a significantly different population of HII regions, showing fainter luminosities. By comparing HII region luminosities with the HST catalog of young star clusters in the central field, we estimate that there is enough Lyman-continuum leakage in the merger to explain the amount of diffuse ionized gas that we detect. We compare the Lyman-continuum escape fraction of each HII region against ionization-parameter sensitive emission line ratios. While we find no systematic trend between these properties, the most extreme line ratios seem to be strong indicators of density bounded ionization. Extrapolating the Lyman-continuum escape fractions to the southern region, we conclude that just from the comparison of the young stellar populations to the ionized gas there is no need to invoke other ionization mechanisms than Lyman-continuum leaking HII regions for the diffuse ionized gas in the Antennae.

We report the detection of the Lyman continuum (LyC) radiation of the compact star-forming galaxy (SFG) J1154+2443 observed with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope. This galaxy, at a redshift of z=0.3690, is characterized by a high emission-line flux ratio O32=[OIII]5007/[OII]3727=11.5. The escape fraction of the LyC radiation fesc(LyC) in this galaxy is 46 per cent, the highest value found so far in low-redshift SFGs and one of the highest values found in galaxies at any redshift. The narrow double-peaked Lya emission line is detected in the spectrum of J1154+2443 with a separation between the peaks Vsep of 199 km/s, one of the lowest known for Lya-emitting galaxies, implying a high fesc(Lya). Comparing the extinction-corrected Lya/Hb flux ratio with the case B value we find fesc(Lya) = 98 per cent. Our observations, combined with previous detections in the literature, reveal an increase of O32 with increasing fesc(LyC). We also find a tight anticorrelation between fesc(LyC) and Vsep. The surface brightness profile derived from the COS acquisition image reveals a bright star-forming region in the centre and an exponential disc in the outskirts with a disc scale length alpha=1.09 kpc. J1154+2443, compared to other known low-redshift LyC leakers, is characterized by the lowest metallicity, 12+logO/H=7.65+/-0.01, the lowest stellar mass M*=10^8.20 Msun, a similar star formation rate SFR=18.9 Msun/yr and a high specific SFR of 1.2x10^-7 yr^-1.

We present the deepest study to date of the Lya luminosity function (LF) in a blank field using blind integral field spectroscopy from MUSE. We constructed a sample of 604 Lya emitters (LAEs) across the redshift range 2.91 < z < 6.64 using automatic detection software in the Hubble Ultra Deep Field. We calculate accurate total Lya fluxes capturing low surface brightness extended Lya emission now known to be a generic property of high-redshift star-forming galaxies. We simulated realistic extended LAEs to characterise the selection function of our samples, and performed flux-recovery experiments to test and correct for bias in our determination of total Lya fluxes. We find an accurate completeness correction accounting for extended emission reveals a very steep faint-end slope of the LF, alpha, down to luminosities of log10 L erg s^-1< 41.5, applying both the 1/Vmax and maximum likelihood estimators. Splitting the sample into three broad redshift bins, we see the faint-end slope increasing from -2.03+1.42-inf at z ~ 3.44 to -2.86+0.76-inf at z ~ 5.48, however no strong evolution is seen between the 68% confidence regions in L*-alpha parameter space. Using the Lya line flux as a proxy for star formation activity, and integrating the observed LFs, we find that LAEs' contribution to the cosmic SFRD rises with redshift until it is comparable to that from continuum-selected samples by z ~ 6. This implies that LAEs may contribute more to the star-formation activity of the early Universe than previously thought - any additional interglactic medium correction would act to further boost the Lya luminosities. Finally, assuming fiducial values for the escape of Lya and LyC radiation, and the clumpiness of the IGM, we integrated the maximum likelihood LF at 5.00 < z < 6.64 and find we require only a small extrapolation beyond the data (< 1 dex in L) for LAEs alone to maintain an ionised IGM at z ~ 6.

We present rest-frame Ly$\alpha$ equivalent widths (EW) of 417 Ly$\alpha$ emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at $2.9 < z < 6.6$ in the Hubble Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured EW values taking into account extended Ly$\alpha$ emission and UV continuum slopes ($\beta$). Our LAEs reach unprecedented depths, both in Ly$\alpha$ luminosities and UV absolute magnitudes, from log($L_{\rm Ly\alpha}$/erg s$^{-1}$) $\sim$41.0 to 43.0 and from Muv $\sim$ -16 to -21 (0.01-1.0 $L^{*}_{\rm z=3}$). The EW values span the range of $\sim$ 5 to 240 \AA\ or larger, and their distribution can be well fitted by an exponential law $N = N_{\rm 0}$ exp($-$EW/$w_{\rm 0}$). Owing to the high dynamic range in Muv, we find that the scale factor, $w_{\rm 0}$, depends on Muv in the sense that including fainter Muv objects increases $w_{\rm 0}$, i.e., the Ando effect. The results indicate that selection functions affect the EW scale factor. Taking these effects into account, we find that our $w_{\rm 0}$ values are consistent with those in the literature within $1\sigma$ uncertainties at $2.9 < z < 6.6$ at a given threshold of Muv and $L_{\rm Ly\alpha}$. Interestingly, we find 12 objects with EW $>200$ \AA\ above $1\sigma$ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak CIV $\lambda 1549$ emission line. For the remaining 10 very large EW LAEs, we find that the EW values can be reproduced by young stellar ages ($< 100$ Myr) and low metallicities ($\lesssim 0.02$ $Z_{\rm \odot}$). Otherwise, at least part of the Ly$\alpha$ emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

We present a measurement of the fraction of Lyman $\alpha$ (Ly$\alpha$) emitters ($X_{\rm{Ly} \alpha}$) amongst HST continuum-selected galaxies at $3<z<6$ with the Multi-Unit Spectroscopic Explorer (MUSE) on the VLT. Making use of the first 24 MUSE-Wide pointings in GOODS-South, each having an integration time of 1 hour, we detect 100 Ly$\alpha$ emitters and find $X_{\rm{Ly} \alpha}\gtrsim0.5$ for most of the redshift range covered, with 29 per cent of the Ly$\alpha$ sample exhibiting rest equivalent widths (rest-EWs) $\leq$ 15\AA. Adopting a range of rest-EW cuts (0 - 75\AA), we find no evidence of a dependence of $X_{\rm{Ly} \alpha}$ on either redshift or UV luminosity.

We report the detection of extended Lyman-alpha (Lya) haloes around 145 individual star-forming galaxies at redshifts 3<z<6 in the Hubble Ultra Deep Field observed with the Multi-Unit Spectroscopic Explorer at ESO-VLT. Our sample consists of continuum-faint (-15> M_{UV}> -22) Lya emitters (LAEs). Using a 2D, two-component decomposition of Lya emission assuming circular exponential distributions, we measure scale lengths and luminosities of Lya haloes. We find that 80% of our objects having reliable Lya halo measurements show Lya emission that is significantly more extended than the UV continuum detected by HST (by a factor ~4 to >20). The median exponential scale length of the Lya haloes in our sample is ~4.5 kpc. By comparing the maximal detected extent of the Lya emission with the predicted dark matter halo virial radii of simulated galaxies, we show that the detected Lya emission of our selected sample of LAEs probes a significant portion of the cold circum-galactic medium (CGM) of these galaxies (>50% in average). This result shows that there must be significant HI reservoirs in the CGM and reinforces the idea that Lya haloes are ubiquitous around high-redshift Lya emitting galaxies. Our characterization of the Lya haloes indicates that the majority of the Lya flux comes from the halo (~65%) and that their scale lengths seem to be linked to the UV properties of the galaxies. We do not observe a significant Lya halo size evolution with redshift. We also find that the Lya lines cover a large range of full width at half maximum (FWHM) from 118 to 512 km/s. While the FWHM does not seem to be correlated to the Lya scale length, most compact Lya haloes and those that are not detected with high significance tend to have narrower Lya profiles. Finally, we investigate the origin of the extended Lya emission but we conclude that our data do not allow us to disentangle the possible processes.

Non-resonant FeII* 2365, 2396, 2612, 2626 emission can potentially trace galactic winds in emission and provide useful constraints to wind models. From the 3'x3' mosaic of the Hubble Ultra Deep Field (UDF) obtained with the VLT/MUSE integral field spectrograph, we identify a statistical sample of 40 FeII* emitters and 50 MgII 2796, 2803 emitters from a sample of 271 [OII] 3726, 3729 emitters with reliable redshifts from z = 0.85 - 1.5 down to 2E-18 (3 sigma) ergs/s/cm^2 (for [OII]), covering the stellar mass range 10^8 - 10^11 Msun. The FeII* and MgII emitters follow the galaxy main sequence, but with a clear dichotomy. Galaxies with masses below 10^9 Msun and star formation rates (SFRs) of <1 Msun/year have MgII emission without accompanying FeII* emission, whereas galaxies with masses above 10^10 Msun and SFRs >10 Msun/year have FeII* emission without accompanying MgII emission. Between these two regimes, galaxies have both MgII and FeII* emission, typically with MgII P-Cygni profiles. Indeed, the MgII profile shows a progression along the main sequence from pure emission to P-Cygni profiles to strong absorption, due to resonant trapping. Combining the deep MUSE data with HST ancillary information, we find that galaxies with pure MgII emission profiles have lower star formation rate surface densities than those with either MgII P-Cygni profiles or FeII* emission. These spectral signatures produced through continuum scattering and fluorescence, MgII P-Cygni profiles and FeII* emission, are better candidates for tracing galactic outflows than pure MgII emission, which may originate from HII regions. We compare the absorption and emission rest-frame equivalent widths for pairs of FeII transitions to predictions from outflow models and find that the observations consistently have less total re-emission than absorption, suggesting either dust extinction or non-isotropic outflow geometries.

The C III] 1907,1909 emission doublet has been proposed as an alternative to Lyman-alpha in redshift confirmations of galaxies at z > 6 since it is not attenuated by the largely neutral intergalactic medium at these redshifts and is believed to be strong in the young, vigorously star-forming galaxies present at these early cosmic times. We present a statistical sample of 17 C III]-emitting galaxies beyond z~1.5 using 30 hour deep VLT/MUSE integral field spectroscopy covering 2 square arcminutes in the Hubble Deep Field South (HDFS) and Ultra Deep Field (UDF), achieving C III] sensitivities of ~2e-17 erg/s/cm^2 in the HDFS and ~7e-18 erg/s/cm^2 in the UDF. The rest-frame equivalent widths range from 2 to 19 Angstroms. These 17 galaxies represent ~3% of the total sample of galaxies found between 1.5 < z < 4. They also show elevated star formation rates, lower dust attenuation, and younger mass-weighted ages than the general population of galaxies at the same redshifts. Combined with deep slitless grism spectroscopy from the HST/WFC3 in the UDF, we can tie the rest-frame ultraviolet C III] emission to rest-frame optical emission lines, namely [O III] 5007, finding a strong correlation between the two. Down to the flux limits that we observe (~1e-18 erg/s/cm^2 with the grism data in the UDF), all objects with a rest-frame [O III] 4959,5007 equivalent width in excess of 250 Angstroms, the so-called Extreme Emission Line Galaxies, have detections of C III] in our MUSE data. More detailed studies of the C III]-emitting population at these intermediate redshifts will be crucial to understand the physical conditions in galaxies at early cosmic times and to determine the utility of C III] as a redshift tracer.

We tested the performance of photometric redshifts for galaxies in the Hubble Ultra Deep field down to 30th magnitude. We compared photometric redshift estimates from three spectral fitting codes from the literature (EAZY, BPZ and BEAGLE) to high quality redshifts for 1227 galaxies from the MUSE integral field spectrograph. All these codes can return photometric redshifts with bias |Dzn|=|z-z_phot|/(1+z)<0.05 down to F775W=30 and spectroscopic incompleteness is unlikely to strongly modify this statement. We have, however, identified clear systematic biases in the determination of photometric redshifts: in the 0.4<z<1.5 range, photometric redshifts are systematically biased low by as much as Dzn=-0.04 in the median, and at z>3 they are systematically biased high by up to Dzn = 0.05, an offset that can in part be explained by adjusting the amount of intergalactic absorption applied. In agreement with previous studies we find little difference in the performance of the different codes, but in contrast to those we find that adding extensive ground-based and IRAC photometry actually can worsen photo-z performance for faint galaxies. We find an outlier fraction, defined through |Dzn|>0.15, of 8% for BPZ and 10% for EAZY and BEAGLE, and show explicitly that this is a strong function of magnitude. While this outlier fraction is high relative to numbers presented in the literature for brighter galaxies, they are very comparable to literature results when the depth of the data is taken into account. Finally, we demonstrate that while a redshift might be of high confidence, the association of a spectrum to the photometric object can be very uncertain and lead to a contamination of a few percent in spectroscopic training samples that do not show up as catastrophic outliers, a problem that must be tackled in order to have sufficiently accurate photometric redshifts for future cosmological surveys.

To reionize the early universe, high-energy photons must escape the galaxies that produce them. It has been suggested that stellar feedback drives galactic outflows out of star-forming regions, creating low density channels through which ionizing photons escape into the inter-galactic medium. We compare the galactic outflow properties of confirmed Lyman continuum (LyC) leaking galaxies to a control sample of nearby star-forming galaxies to explore whether the outflows from leakers are extreme as compared to the control sample. We use data from the Cosmic Origins Spectrograph on the Hubble Space Telescope to measure the equivalent widths and velocities of Si II and Si III absorption lines, tracing neutral and ionized galactic outflows. We find that the Si II and Si III equivalent widths of the LyC leakers reside on the low-end of the trend established by the control sample. The leakers' velocities are not statistically different than the control sample, but their absorption line profiles have a different asymmetry: their central velocities are closer to their maximum velocities. The outflow kinematics and equivalent widths are consistent with the scaling relations between outflow properties and host galaxy properties -- most notably metallicity -- defined by the control sample. Additionally, we use the Ly\alpha\ profiles to show that the Si II equivalent width scales with the Ly\alpha\ peak velocity separation. We determine that the low equivalent widths of the leakers are likely driven by low metallicities and low H I column densities, consistent with a density-bounded ionization region, although we cannot rule out significant variations in covering fraction. While we do not find that the LyC leakers have extreme outflow velocities, the low maximum-to-central velocity ratios demonstrate the importance of the acceleration and density profiles for LyC and Ly\alpha\ escape. [abridged]

We selected a sample of 76 Lya emitting galaxies from the VIMOS Ultra Deep Survey (VUDS) at 2<z<4. We estimated the velocity of the neutral gas flowing out of the interstellar medium as the velocity offset, Deltav, between the systemic redshift (zsys) and the center of low-ionization absorption line systems (LIS). To increase the SN of VUDS spectra, we stacked subsamples. We measured the systemic redshift from the rest-frame UV spectroscopic data using the CIII]1908 nebular emission line, and we considered SiII1526 as the highest signal-to-noise LIS line. We calculated the Lya peak shift with respect to the zsys, the EW(Lya), and the Lya spatial extension, Ext(Lya-C), from the profiles in the 2D stacked spectra. The galaxies that are faint in the rest-frame UV continuum, strong in Lya and CIII], with compact UV morphology, and localized in an underdense environment are characterized by outflow velocities of the order of a few hundreds of km/sec. The subsamples with smaller Deltav are characterized by larger Lya peak shifts, larger Ext(Lya-C), and smaller EW(Lya). In general we find that EW(Lya) anti-correlates with Ext(Lya-C) and Lya peak shift. We interpret these trends using a radiative-transfer shell model. The model predicts that an HI gas with a column density larger than 10^20/cm^2 is able to produce Lya peak shifts larger than >300km/sec. An ISM with this value of NHI would favour a large amount of scattering events, especially when the medium is static, so it can explain large values of Ext(Lya-C) and small EW(Lya). On the contrary, an ISM with a lower NHI, but large velocity outflows would lead to a Lya spatial profile peaked at the galaxy center (i.e. low values of Ext(Lya-C)) and to a large EW(Lya), as we see in our data. Our results and their interpretation via radiative-transfer models tell us that it is possible to use Lya to study the properties of the HI gas.

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.

Emission signatures from galactic winds provide an opportunity to directly map the outflowing gas, but this is traditionally challenging because of the low surface brightness. Using deep observations (27 hours) of the Hubble Deep Field South from the Multi Unit Spectroscopic Explorer (MUSE) instrument, we identify signatures of an outflow in both emission and absorption from a spatially resolved galaxy at z = 1.29 with a stellar mass M* = 8 x 10^9 Msun, star formation rate SFR = 77 Msun/yr, and star formation rate surface brightness 1.6 Msun/kpc^2 within the [OII] half-light radius R_1/2,[OII] = 2.76 +- 0.17 kpc. From a component of the strong resonant MgII and FeII absorptions at -350 km/s, we infer a mass outflow rate that is comparable to the star formation rate. We detect non-resonant FeII* emission, at lambda 2626, 2612, 2396, and 2365, at 1.2-2.4-1.5-2.7 x 10^-18 egs s-1 cm-2 respectively. These flux ratios are consistent with the expectations for optically thick gas. By combining the four non-resonant FeII* emission lines, we spatially map the FeII* emission from an individual galaxy for the first time. The FeII* emission has an elliptical morphology that is roughly aligned with the galaxy minor kinematic axis, and its integrated half-light radius R_1/2,FeII* = 4.1 +- 0.4 kpc is 50% larger than the stellar continuum (R_1/2,* = 2.34 +- 0.17 kpc) or the [OII] nebular line. Moreover, the FeII* emission shows a blue wing extending up to -400 km/s, which is more pronounced along the galaxy minor kinematic axis and reveals a C-shaped pattern in a p-v diagram along that axis. These features are consistent with a bi-conical outflow.

We have recently reported the discovery of five low redshift Lyman continuum (LyC) emitters (LCEs, hereafter) with absolute escape fractions fesc(LyC) ranging from 6 to 13%, higher than previously found, and which more than doubles the number of low redshift LCEs.We use these observations to test theoretical predictions about a link between the characteristics of the Lyman-alpha (Lya) line from galaxies and the escape of ionising photons. We analyse the Lya spectra of eight LCEs of the local Universe observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope (our five leakers and three galaxies from the litterature), and compare their strengths and shapes to the theoretical criteria and comparison samples of local galaxies: the Lyman Alpha Reference Survey, Lyman Break Analogs, Green Peas, and the high-redshift strong LyC leaker Ion2. Our LCEs are found to be strong Lya emitters, with high equivalent widths, EW(Lya)> 70 {\AA}, and large Lya escape fractions, fesc(Lya) > 20%. The Lya profiles are all double-peaked with a small peak separation, in agreement with our theoretical expectations. They also have no underlying absorption at the Lya position. All these characteristics are very different from the Lya properties of typical star-forming galaxies of the local Universe. A subset of the comparison samples (2-3 Green Pea galaxies) share these extreme values, indicating that they could also be leaking. We also find a strong correlation between the star formation rate surface density and the escape fraction of ionising photons, indicating that the compactness of star-forming regions plays a role in shaping low column density paths in the interstellar medium of LCEs. The Lya properties of LCEs are peculiar: Lya can be used as a reliable tracer of LyC escape from galaxies, in complement to other indirect diagnostics proposed in the literature.

Direct Ly $\alpha$ imaging of intergalactic gas at $z\sim2$ has recently revealed giant cosmological structures around quasars, e.g. the Slug Nebula (Cantalupo et al. 2014). Despite their high luminosity, the detection rate of such systems in narrow-band and spectroscopic surveys is less than 10%, possibly encoding crucial information on the distribution of gas around quasars and the quasar emission properties. In this study, we use the MUSE integral-field instrument to perform a blind survey for giant Ly $\alpha$ nebulae around 17 bright radio-quiet quasars at $3<z<4$ that does not suffer from most of the limitations of previous surveys. After data reduction and analysis performed with specifically developed tools, we found that each quasar is surrounded by giant Ly $\alpha$ nebulae with projected sizes larger than 100 physical kpc and, in some cases, extending up to 320 kpc. The circularly averaged surface brightness profiles of the nebulae appear very similar to each other despite their different morphologies and are consistent with power laws with slopes $\approx-1.8$. The similarity between the properties of all these nebulae and the Slug Nebula suggests a similar origin for all systems and that a large fraction of gas around bright quasars could be in a relatively "cold" (T$\sim$10$^4$K) and dense phase. In addition, our results imply that such gas is ubiquitous within at least 50 kpc from bright quasars at $3<z<4$ independently of the quasar emission opening angle, or extending up to 200 kpc for quasar isotropic emission.

We have recently discovered five Lyman continuum leaking galaxies at z~0.3, selected for their compactness, intense star-formation, and high [OIII]/[OII] ratio (Izotov et al. 2016ab). Here we derive their ionizing photon production efficiency, a fundamental quantity for inferring the number of photons available to reionize the Universe, for the first time for galaxies with confirmed strong Lyman continuum escape (fesc~6-13%). We find an ionizing photon production per unit UV luminosity, which is a factor 2-6 times higher than the canonical value when reported to their observed UV luminosity. After correction for extinction this value is close to the canonical value. The properties of our five Lyman continuum leakers are found to be very similar to those of the confirmed z=3.218 leaker Ion2 from de Barros et al. (2016) and very similar to those of typical star-forming galaxies at z>~6. Our results suggest that UV bright galaxies at high-z such as Lyman break galaxies can be Lyman continuum leakers and that their contribution to cosmic reionization may be underestimated.

Following our first detection reported in Izotov et al. (2016), we present the detection of Lyman continuum (LyC) radiation of four other compact star-forming galaxies observed with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope (HST). These galaxies, at redshifts of z~0.3, are characterized by high emission-line flux ratios [OIII]5007/[OII]3727 > 5. The escape fractions of the LyC radiation fesc(LyC) in these galaxies are in the range of ~6%-13%, the highest values found so far in low-redshift star-forming galaxies. Narrow double-peaked Lyalpha emission lines are detected in the spectra of all four galaxies, compatible with predictions for Lyman continuum leakers. We find escape fractions of Lyalpha, fesc(Lyalpha) ~20%-40%, among the highest known for Lyalpha emitters (LAEs). Surface brightness profiles produced from the COS acquisition images reveal bright star-forming regions in the center and exponential discs in the outskirts with disc scale lengths alpha in the range ~0.6-1.4 kpc. Our galaxies are characterized by low metallicity, ~1/8-1/5 solar, low stellar mass ~(0.2 - 4)e9 Msun, high star formation rates SFR~14-36 Msun/yr, and high SFR densities Sigma~2-35 Msun/yr/kpc^2. These properties are comparable to those of high-redshift star-forming galaxies. Finally, our observations, combined with our first detection reported in Izotov et al. (2016), reveal that a selection for compact star-forming galaxies showing high [OIII]5007/[OII]3727 ratios appears to pick up very efficiently sources with escaping Lyman continuum radiation: all five of our selected galaxies are LyC leakers.

We present VLT/X-Shooter and MUSE spectroscopy of an faint F814W=28.60+/-0.33 (Muv=-17.0), low mass (~<10^7 Msun) and compact (Reff=62pc) freshly star-forming galaxy at z=3.1169 magnified (16x) by the Hubble Frontier Fields galaxy cluster Abell S1063. Gravitational lensing allows for a significant jump toward low-luminosity regimes, in moderately high resolution spectroscopy (R=lambda/dlambda ~ 3000-7400). We measured CIV1548,1550, HeII1640, OIII]1661,1666, CIII]1907,1909, Hbeta, [OIII]4959,5007, emission lines with FWHM< 50 km/s and (de-lensed) fluxes spanning the interval 1.0x10^-19 - 2.0x10^-18 erg/s/cm2 at S/N=4-30. The double peaked Lya emission with Delta_v(red-blue) = 280(+/-7)km/s and de-lensed fluxes 2.4_(blue)|8.5_(red)x10^-18 erg/s/cm2 (S/N=38_(blue)|110_(red)) indicate a low column density of neutral hydrogen gas consistent with a highly ionized interstellar medium as also inferred from the large [OIII]5007/[OII]3727>10 ratio. We detect CIV1548,1550 resonant doublet in emission, each component with FWHM ~< 45 km/s, and redshifted by +51(+/-10)km/s relative to the systemic redshift. We interpret this as nebular emission tracing an expanding optically-thin interstellar medium. Both CIV1548,1550 and HeII1640 suggest the presence of hot and massive stars (with a possible faint AGN). The ultraviolet slope is remarkably blue, beta =-2.95 +/- 0.20 (F_lambda=lambda^beta), consistent with a dust-free and young ~<20 Myr galaxy. Line ratios suggest an oxygen abundance 12+log(O/H)<7.8. We are witnessing an early episode of star-formation in which a relatively low NHI and negligible dust attenuation might favor a leakage of ionizing radiation. This galaxy currently represents a unique low-luminosity reference object for future studies of the reionization epoch with JWST.

Lyman-$\alpha$ nebulae are usually found in massive environments at high redshift ($z > 2$). The origin of their Lyman-$\alpha$ (Lya) emission remains debated. Recent polarimetric observations showed that at least some Lya sources are polarized. This is often interpreted as a proof that the photons are centrally produced, and opposed to the scenario in which the Lya emission is the cooling radiation emitted by gas heated during the accretion onto the halo. We suggest that this scenario is not incompatible with the polarimetric observations. In order to test this idea, we post-process a radiative hydrodynamics simulation of a blob with the MCLya Monte Carlo transfer code. We compute radial profiles for the surface brightness and the degree of polarization and compare them to existing observations. We find that both are consistent with a significant contribution of the extragalactic gas to the Lya emission. Most of the photons are centrally emitted and scattered inside the filament afterwards, producing the observed high level of polarization. We argue that the contribution of the extragalactic gas to the Lya emission does not prevent polarization to arise. On the contrary, we find that pure galactic emission causes the polarization profile to be too steep to be consistent with observations.

We present the results obtained with MUSE on the core of the lensing cluster A1689. Integral-field observations with MUSE provide a unique view of the central region, allowing us to conduct a complete census on both cluster galaxies and lensed background sources, identified based on their spectral features without preselection. We investigate the multiple-image configuration for all known sources in the field. Previous to our survey, 28 different lensed galaxies displaying 46 multiple images were known in the MUSE field of view, most of them based on photometric redshifts and lensing considerations. Among them, we spectroscopically confirm 12 images based on their emission-lines, corresponding to 7 different lensed galaxies between z = 0.95 and 5.0. In addition, 14 new galaxies have been spectroscopically identified in this area, with redshifts ranging between 0.8 and 6.2. All background sources within the MUSE field of view correspond to multiple-imaged systems lensed by A1689. 17 sources in total are found at z > 3 based on their Lyman-alpha emission, with Lyman-alpha luminosities ranging between 40.5 < log(Ly{\alpha}) < 42.5 after correction for magnification. This sample is particularly sensitive to the slope of the LF toward the faintest-end. The density of sources obtained in this survey is consistent with a steep value of {\alpha} < -1.5, although this result still needs further investigation. These results illustrate the efficiency of MUSE in the characterization of lensing clusters on one hand, and the study of faint and distant populations of galaxies on the other hand. In particular, our current survey of lensing clusters should provide a unique census of sources responsible for the reionization in a representative volume at z ~ 4-7.

We present the results of a Lya profile analysis of 12 Lya emitters (LAEs) at z = 2.2 with high-resolution Lya spectra. We find that all 12 objects have a Lya profile with the main peak redward of the systemic redshift defined by nebular lines, and five have a weak, secondary peak blueward of the systemic redshift (blue bump). The average velocity offset of the red main peak (the blue bump, if any) with respect to the systemic redshift is Delta_v_Lya,r = 174+/- 19 km s-1 (Delta_v_Lya,b = -316+/-45 km s-1), which is smaller than (comparable to) that of Lyman-break galaxies (LBGs). The outflow velocities inferred from metal absorption lines in three individual and one stacked spectra are comparable to those of LBGs. The expanding shell model constructed by Verhamme et al. (2006) reproduces not only the Lya profiles but also other observed quantities including the outflow velocity and the FWHM of nebular lines for the non-blue bump objects. On the other hand, the model predicts too high FWHMs of nebular lines for the blue bump objects, although this discrepancy may disappear if we introduce additional Lya photons produced by gravitational cooling. We show that the small Delta_v_Lya,r values of our sample can be explained by low neutral-hydrogen column densities of log(NHI) = 18.9 cm-2 on average. This value is more than one order of magnitude lower than those of LBGs but is consistent with recent findings that LAEs have high ionization parameters and low Hi gas masses. This result suggests that low NHI values, giving reduced numbers of resonant scattering of Lya photons, are the key to the strong Lya emission of LAEs.

One of the key questions in observational cosmology is the identification of the sources responsible for ionisation of the Universe after the cosmic Dark Ages, when the baryonic matter was neutral. The currently identified distant galaxies are insufficient to fully reionise the Universe by redshift z~6, but low-mass star-forming galaxies are thought to be responsible for the bulk of the ionising radiation. Since direct observations at high redshift are difficult for a variety of reasons, one solution is to identify local proxies of this galaxy population. However, starburst galaxies at low redshifts are generally opaque to their ionising radiation. This radiation with small escape fractions of 1-3% is directly detected only in three low-redshift galaxies. Here we present far-ultraviolet observations of a nearby low-mass star-forming galaxy, J0925+1403, selected for its compactness and high excitation. The galaxy is leaking ionising radiation, with an escape fraction of ~8%. The total number of photons emitted during the starburst phase is sufficient to ionize intergalactic medium material, which is about 40 times more massive than the stellar mass of the galaxy.

We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.

Spatially resolved studies of high redshift galaxies, an essential insight into galaxy formation processes, have been mostly limited to stacking or unusually bright objects. We present here the study of a typical (L$^{*}$, M$_\star$ = 6 $\times 10^9$ $M_\odot$) young lensed galaxy at $z=3.5$, observed with MUSE, for which we obtain 2D resolved spatial information of Ly$\alpha$ and, for the first time, of CIII] emission. The exceptional signal-to-noise of the data reveals UV emission and absorption lines rarely seen at these redshifts, allowing us to derive important physical properties (T$_e\sim$15600 K, n$_e\sim$300 cm$^{-3}$, covering fraction f$_c\sim0.4$) using multiple diagnostics. Inferred stellar and gas-phase metallicities point towards a low metallicity object (Z$_{\mathrm{stellar}}$ = $\sim$ 0.07 Z$_\odot$ and Z$_{\mathrm{ISM}}$ $<$ 0.16 Z$_\odot$). The Ly$\alpha$ emission extends over $\sim$10 kpc across the galaxy and presents a very uniform spectral profile, showing only a small velocity shift which is unrelated to the intrinsic kinematics of the nebular emission. The Ly$\alpha$ extension is $\sim$4 times larger than the continuum emission, and makes this object comparable to low-mass LAEs at low redshift, and more compact than the Lyman-break galaxies and Ly$\alpha$ emitters usually studied at high redshift. We model the Ly$\alpha$ line and surface brightness profile using a radiative transfer code in an expanding gas shell, finding that this model provides a good description of both observables.

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.

[Abridged] We investigate the physical properties of a Lyman continuum emitter candidate at $z=3.212$ with photometric coverage from $U$ to MIPS 24$\mu$m band and VIMOS/VLT and MOSFIRE/Keck spectroscopy. Investigation of the UV spectrum confirms a direct spectroscopic detection of the Lyman continuum emission with $S/N>5$. Non-zero Ly$\alpha$ flux at the systemic redshift and high Lyman-$\alpha$ escape fraction suggest a low HI column density. The weak C and Si low-ionization absorption lines are also consistent with a low covering fraction along the line of sight. The [OIII]$\lambda\lambda4959,5007+\mathrm{H}\beta$ equivalent width is one of the largest reported for a galaxy at $z>3$ ($\mathrm{EW}([\mathrm{OIII}]\lambda\lambda4959,5007+\mathrm{H}\beta) \simeq 1600\AA$, rest-frame) and the NIR spectrum shows that this is mainly due to an extremely strong [OIII] emission. The large observed [OIII]/[OII] ratio ($>10$) and high ionization parameter are consistent with prediction from photoionization models in case of a density-bounded nebula scenario. Furthermore, the $\mathrm{EW}([\mathrm{OIII}]\lambda\lambda4959,5007+\mathrm{H}\beta)$ is comparable to recent measurements reported at $z\sim7-9$, in the reionization epoch. We also investigate the possibility of an AGN contribution to explain the ionizing emission but most of the AGN identification diagnostics suggest that stellar emission dominates instead. This source is currently the first high-$z$ example of a Lyman continuum emitter exhibiting indirect and direct evidences of a Lyman continuum leakage and having physical properties consistent with theoretical expectation from Lyman continuum emission from a density-bounded nebula.

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.

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.

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

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.

[abridged] Aims. The aim of this paper is to assess the impact of the interstellar medium (ISM) physics on Lyman-alpha (Lya) radiation transfer and to quantify how galaxy orientation with respect to the line of sight alters observational signatures. Methods. We compare the results of Lya radiation transfer calculations through the ISM of a couple of idealized galaxy simulations with different ISM models. Results. First, the small-scale structuration of the ISM plays a determinant role in shaping a galaxys Lya properties.The artificially warm, and hence smooth, ISM of G1 yields an escape fraction of 50 percent at the Lya line center, and produces symmetrical double-peak profiles. On the contrary, in G2, most young stars are embedded in thick star-forming clouds, and the result is a 10 times lower escape fraction. G2 also displays a stronger outflowing velocity field, which favors the escape of red-shifted photons, resulting in an asymmetric Lya line. Second, the Lya properties of G2 strongly depend on the inclination at which it is observed: From edge-on to face-on, the line goes from a double-peak profile with an equivalent width of -5 Angstrom to a 15 times more luminous red-shifted asymmetric line with EW 90 Angstrom. Conclusions. Lya radiation transfer calculations can only lead to realistic properties in simulations where galaxies are resolved into giant molecular clouds, putting these calculations out of reach of current large scale cosmological simulations. Finally, we find inclination effects to be much stronger for Lya photons than for continuum radiation. This could potentially introduce severe biases in the selection function of narrow-band Lya emitter surveys, which could indeed miss a significant fraction of the high-z galaxy population.

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.

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.

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

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.

The Aim of our study is to understand the variety of observed Lyman-alpha (Lya) line profiles and strengths in Lyman Break Galaxies (LBGs) and Lya emitters (LAEs), the physical parameters governing them, and hence deriving constraints on the gas and dust content and stellar populations of these objects. Using our 3D Lya radiation transfer code including gas and dust (Verhamme et al. 2006), we fit 11 LBGs from the FORS Deep Field with redshifts between 2.8 and 5 observed by Tapken et al. (2007). A simple geometry of a spherically expanding shell of HI is adopted. RESULTS : The variety of observed Lya profiles is successfully reproduced. Most objects show outflow velocities of 150-200 km/s; two objects are most likely quasi-static. The radial HI column density ranges from NH=2.10^{19} to 7.10^{20} cm^{-2}. Our Lya profile fits yield values of E(B-V)~0.05-0.2 for the gas extinction. We find indications for a dust-to-gas ratio higher than the Galactic value, and for a substantial scatter. The escape fraction of Lya photons is found to be determined primarily by the extinction, and a simple fit formula is proposed. Intrinsic EW(Lya)~50-100 Angstroms are found for 8/11 objects, as expected for stellar populations forming constantly over long periods (> 10-100 Myr). In three cases we found indications for younger populations. Correlations between the observed EW(Lya) and other observables such as FWHM(Lya), E(B-V),SFR(UV) etc, are reproduced. We also show that there is a clear overlap between LBGs and LAEs. Radiation transfer and dust effects explain the increase of the LAE/LBG ratio, and a higher percentage of LBGs with strong Lya emission with increasing redshift. [shortened]

Using our 3D Lya radiation transfer code, we compute the radiation transfer of Lya and UV continuum photons including dust. Observational constraints on the neutral gas (column density, kinematics, etc.) are taken from other analysis of this object. RESULTS: The observed Lya profile of MS 1512--cB58 is reproduced for the first time taking radiation transfer and all observational constraints into account. The observed absorption profile is found to result naturally from the observed amount of dust and the relatively high HI column density. Radiation transfer effects and suppresion by dust transform a strong intrinsic Lya emission with EW(Lya)>~ 60 Ang into the observed faint superposed Lya emission peak. We propose that the vast majority of LBGs have intrinsically EW(Lya)~60-80 Ang or larger, and that the main physical parameter responsible for the observed variety of Lya strengths and profiles in LBGs is N_H and the accompanying variation of the dust content. Observed EW(Lya) distributions, Lya luminosity functions, and related quantities must therefore be corrected for radiation transfer and dust effects. The implications from our scenario on the duty-cycle of Lya emitters are also discussed.

Using a Monte Carlo technique, we have developed a 3D lyman-alpha radiation transfer code allowing for prescribed arbitrary hydrogen density, ionisation, temperature structures, and dust distribution, and arbitrary velocity fields and UV photon sources. We have examined the lyman-alpha line profiles predicted for several simple geometrical configurations and their dependence on the main input parameters. Overall, we find line profiles reaching from doubly peaked symmetric emission to symmetric Voigt (absorption) in static configurations with increasing dust content, and asymmetric red-(blue-) shifted emission lines with a blue (red) counterpart ranging from absorption to emission (with increasing line/continuum strength) in expanding (infalling) media. The following results are of interest for the interpretation of lya profiles from galaxies. 1/ Standard lya absorption line fitting of global spectra of galaxies may lead to an underestimate of the true hydrogen column density in certain geometrical conditions. 2/ Normal (inverted) P-Cygni like lya profiles can be obtained in expanding (infalling) media from objects without any intrinsic lya emission, as a natural consequence of radiation transfer effects. 3/ The formation and the detailed shape of lya profiles resulting from expanding shells has been thoroughly revised: for sufficiently large column densities, the position of the main lya emission peak is redshifted by twice the expansion velocity.This is in excellent agreement with the observations of z~3 LBGs. This finding indicates also that large scale, fairly symmetric shell structures must be a good description for the outflows in LBGs.(shortened abstract)

This series of papers is intended to present astrocladistics in some detail and evaluate this methodology in reconstructing phylogenies of galaxies. Being based on the evolution of all the characters describing galaxies, it is an objective way of understanding galaxy diversity through evolutionary relationships. In this first paper, we present the basic steps of a cladistic analysis and show both theoretically and practically that it can be applied to galaxies. For illustration, we use a sample of 50 simulated galaxies taken from the GALICS database, which are described by 91 observables (dynamics, masses and luminosities). These 50 simulated galaxies are indeed 10 different galaxies taken at 5 cosmological epochs, and they are free of merger events. The astrocladistic analysis easily reconstructs the true chronology of evolution relationships within this sample. It also demonstrates that burst characters are not relevant for galaxy evolution as a whole. A companion paper is devoted to the formalization of the concepts of formation and diversification in galaxy evolution.

This series of papers is intended to evaluate astrocladistics in reconstructing phylogenies of galaxies. The objective of this second paper is to formalize the concept of galaxy formation and to identify the processes of diversification. We show that galaxy diversity can be expected to organize itself in a hierarchy. In order to better understand the role of mergers, we have selected a sample of 43 galaxies from the GALICS database built from simulations with a hybrid model for galaxy formation studies. These simulated galaxies, described by 119 characters and considered as representing still undefined classes, have experienced different numbers of merger events during evolution. Our cladistic analysis yields a robust tree that proves the existence of a hierarchy. Mergers, like interactions (not taken into account in the GALICS simulations), are probably a strong driver for galaxy diversification. Our result shows that mergers participate in a branching type of evolution, but do not seem to play the role of an evolutionary clock.