We present the Lya luminosity functions (LFs) derived by our deep Subaru narrowband survey that identifies a total of 3,137 Lya emitters (LAEs) at $z = 2.2$ in five independent blank fields. The sample of these LAEs is the largest, to date, and covers a very wide Lya luminosity range of $\log L_{Ly\alpha} = 41.7-44.4$ erg s$^{-1}$. We determine the Lya LF at $z = 2.2$ with unprecedented accuracies, and obtain the best-fit Schechter parameters of $L^{*}_{Ly\alpha} = 5.29^{+1.67}_{-1.13} \times 10^{42}$ erg s$^{-1}$, $\phi^{*}_{Ly\alpha} = 6.32^{+3.08}_{-2.31} \times 10^{-4}$ Mpc$^{-3}$, and $\alpha = -1.75^{+0.10}_{-0.09}$ showing a steep faint-end slope. We identify a significant hump at the LF bright end ($\log L_{Ly\alpha} > 43.4$ erg s$^{-1}$). Because all of the LAEs in the bright-end hump have (a) bright counterpart(s) either in the X-ray, UV, or radio data, this bright-end hump is not made by gravitational lensing magnification bias but AGNs. These AGNs allow us to derive the AGN UV LF at $z \sim 2$ down to the faint magnitude limit of $M_{UV} \simeq -22.5$, and to constrain the faint-end slope of AGN UV LF, $\alpha_{AGN}=-1.2 \pm 0.1$, that is flatter than those at $z > 4$. Based on the Lya and UV LFs from our and previous studies, we find the increase of Lya escape fraction $f^{Ly\alpha}_{esc}$ from $z \sim 0$ to $6$ by two orders of magnitude. This large $f^{Ly\alpha}_{esc}$ increase can be explained neither by the evolution of stellar population nor outflow alone, but the evolution of neutral hydrogen HI density in inter-stellar medium that enhances dust attenuation for Lya by resonance scattering. Our uniform expanding shell models suggest that the typical HI column density decreases from $N_{HI} \sim 7 \times 10^{19}$ ($z \sim 0$) to $\sim 1 \times 10^{18}$ cm$^{-2}$ ($z \sim 6$) to explain the large $f^{Ly\alpha}_{esc}$ increase.

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.

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

Aims: We aim to study the effects of an inhomogeneous interstellar medium (ISM) on the strength and the shape of the Lyman alpha (Lya) line in starburst galaxies. Methods: Using our 3D Monte Carlo Lya radiation transfer code, we study the radiative transfer of Lya, UV and optical continuum photons in homogeneous and clumpy shells of neutral hydrogen and dust surrounding a central source. Our simulations predict the Lya and continuum escape fraction, the Lya equivalent width EW(Lya), the Lya line profile and their dependence on the gas geometry and the main input physical parameters. Results: The ISM clumpiness is found to have a strong impact on the Lya line radiative transfer, entailing a strong dependence of the emergent features of the Lya line (escape fraction, EW(Lya)) on the ISM morphology. Although a clumpy and dusty ISM appears more transparent to radiation (both line and continuum) compared to an equivalent homogeneous ISM of equal dust optical depth, we find that the Lya photons are, in general, still more attenuated than UV continuum radiation. As a consequence, the observed Lya equivalent width (EWobs(Lya)) is lower than the intrinsic one (EWint(Lya)) for nearly all clumpy ISM configurations considered. There are, however, special conditions under which Lya photons escape more easily than the continuum, resulting in an enhanced EWobs(Lya). The requirement for this to happen is that the ISM is almost static (galactic outflows < 200 km/s), extremely clumpy (with density contrasts >10^7 in HI between clumps and the interclump medium), and very dusty (E(B-V) > 0.30). When these conditions are fulfilled the emergent Lya line profile shows no velocity shift and little asymmetry. Given the asymmetry and velocity shifts generally observed in starburst galaxies with Lya emission, we conclude that clumping is unlikely to significantly enhance their relative Lya/UV transmission.

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.

It has been suggested that radiative transfer effects may explain the unusually high equivalent widths (EWs) of the Lya line, observed occasionally from starburst galaxies, especially at high redshifts. If the dust is locked up inside high-density clouds dispersed in an empty intercloud medium, the Lya photons could scatter off of the surfaces of the clouds, effectively having their journey confined to the dustless medium. The continuum radiation, on the other hand, does not scatter, and would thus be subject to absorption inside the clouds. This scenario is routinely invoked when Lya EWs higher than what is expected theoretically are observed, although the ideal conditions under which the results are derived usually are not considered. Here we systematically examine the relevant physical parameters in this idealized framework, testing whether any astrophysically realistic scenarios may lead to such an effect. It is found that although clumpiness indeed facilitates the escape of Lya, it is highly unlikely that any real interstellar media should result in a preferential escape of Lya over continuum radiation. Other possible causes are discussed, and it is concluded that the observed high EWs are more likely to be caused by cooling radiation from cold accretion and/or anisotropic escape of the Lya radiation.

Population III galaxies, made partly or exclusively of metal-free stars, are predicted to exist at high redshifts and may produce very strong Lya emission. A substantial fraction of these Lya photons are likely absorbed in the intergalactic medium at z>6, but recent simulations suggest that significant Lya emission may be detectable up to z~8.5, i.e. well into the reionization epoch. Here, we argue that high-redshift population III galaxies with strong Lya emission can be identified in Hubble Space Telescope imaging data because of their unusual colours. We quantify this effect in some of the filters used in Y-band dropout searches for galaxies at z~8 and find that population III galaxies with high Lya fluxes may exhibit much bluer J-H colours at z=8-10 than any normal type of galaxy at these redshifts. This colour signature can arise even if pop III stars account for as little as ~1e-3 to ~1e-2 of the stellar mass in these galaxies. Some of the anomalously blue objects reported in current Y-band dropout samples do in fact meet the colour criteria for Lya-emitting population III galaxies.