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

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

We have used the capability of the MUSE instrument to explore the impact of stellar feedback at large scales in Haro 11, a galaxy under extreme starburst condition and one of the first galaxies where Lyman continuum (LyC) has been detected. Using Ha, [OIII] and [OI] emission lines from deep MUSE observations, we have constructed a sequence of velocity-dependent maps of the Ha emission, the state of the ionised gas and a tracer of fast shocks. These allowed us to investigate the ionisation structure of the galaxy in 50 kms^2 bins over a velocity range of -400 to 350 kms. The ionised gas in Haro 11 is assembled by a rich arrangement of structures, such as superbubbles, filaments, arcs and galactic ionised channels, whose appearances change drastically with velocity. The central star forming knots and the star forming dusty arm are the main engines that power the strong mechanical feedback in this galaxy, although with different impact on the ionisation structure. Haro 11 appears to leak LyC radiation in many directions. We found evidence of a kpc-scale fragmented superbubble, that may have cleared galactic-scale channels in the ISM. Additionally, the southwestern hemisphere is highly ionised in all velocities, hinting at a density bound scenario. A compact kpc-scale structure of lowly ionised gas coincides with the diffuse Lya emission and the presence of fast shocks. Finally, we find evidence that a significant fraction of the ionised gas mass may escape the gravitational potential of the galaxy.

We study young star-forming clumps on physical scales of 10-500 pc in the Lyman-Alpha Reference Sample (LARS), a collection of low-redshift (z = 0.03-0.2) UV-selected star-forming galaxies. In each of the 14 galaxies of the sample, we detect clumps for which we derive sizes and magnitudes in 5 UV-optical filters. The final sample includes $\sim$1400 clumps, of which $\sim$600 have magnitude uncertainties below 0.3 in all filters. The UV luminosity function for the total sample of clumps is described by a power-law with slope $\alpha = -2.03^{+0.11}_{-0.13}$. Clumps in the LARS galaxies have on average $\Sigma_{SFR}$ values higher than what observed in HII regions of local galaxies and comparable to typical SFR densities of clumps in z = 1-3 galaxies. We derive the clumpiness as the relative contribution from clumps to the UV emission of each galaxy, and study it as a function of galactic-scale properties, i.e. $\Sigma_{SFR}$ and the ratio between rotational and dispersion velocities of the gas ($v_s/\sigma_0$). We find that in galaxies with higher $\Sigma_{SFR}$ or lower $v_s/\sigma_0$, clumps dominate the UV emission of their host systems. All LARS galaxies with Ly$\alpha$ escape fractions larger than 10% have more than 50% of the UV luminosity from clumps. We tested the robustness of these results against the effect of different physical resolutions. At low resolution, the measured clumpiness appears more elevated than if we could resolve clumps down to single clusters. This effect is small in the redshift range covered by LARS, thus our results are not driven by the physical resolution.

Stellar feedback strongly affects the interstellar medium (ISM) of galaxies. Stellar feedback in the first galaxies likely plays a major role in enabling the escape of LyC photons, which contribute to the re-ionization of the Universe. Nearby starburst galaxies serve as local analogues allowing for a spatially resolved assessment of the feedback processes in these galaxies. We characterize the feedback effects from the star clusters in the local high-redshift analogue ESO 338-IG04 on the ISM and compare the results with the properties of the most massive clusters. We use high quality VLT/MUSE optical integral field data to derive the physical properties of the ISM such as ionization, density, shocks, and perform new fitting of the spectral energy distributions of the brightest clusters in ESO 338-IG04 from HST imaging. ESO 338-IG04 has a large ionized halo which we detect to a distance of 9 kpc. We identify 4 Wolf-Rayet (WR) clusters based on the blue and red WR bump. We follow previously identified ionization cones and find that the ionization of the halo increases with distance. Analysis of the galaxy kinematics shows two complex outflows driven by the numerous young clusters in the galaxy. We find a ring of shocked emission traced by an enhanced [OI]/H$\alpha$ ratio surrounding the starburst and at the end of the outflow. Finally we detect nitrogen enriched gas associated with the outflow, likely caused by the WR stars in the massive star clusters. Photo-ionization dominates the central starburst and sets the ionization structure of the entire halo, resulting in a density bounded halo, facilitating the escape of LyC photons. Outside the central starburst, shocks triggered by an expanding super bubble become important. The shocks at the end of the outflow suggest interaction between the hot outflowing material and the more quiescent halo gas.

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

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

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

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

While the Lyman-$\alpha$ ($\mathrm{Ly}\alpha$) emission line serves as an important tool in the study of galaxies at $z\lesssim 6$, finding Ly$\alpha$ emitters (LAE) at significantly higher redshifts has been more challenging, probably because of the increasing neutrality of the intergalactic medium above $z\sim6$. Galaxies with extremely high rest-frame Ly$\alpha$ equivalent widths, EW(Ly$\alpha$) $\gtrsim 150$ \AA{}, at $z>6$ are good candidates for Ly$\alpha$ follow-up observations, and can stand out in multiband imaging surveys because of their unusual colors. We have conducted a photometric search for such objects in the Cluster Lensing And Supernova survey with Hubble (CLASH), and report here the identification of three likely gravitationally-lensed images of a single LAE candidate at $z\sim6.3$, behind the galaxy cluster Abell 2261($z = 0.225$). In the process, we also measured with Keck/MOSFIRE the first spectroscopic redshift of a multiply-imaged galaxy behind Abell 2261, at $z = 3.337$. This allows us to calibrate the lensing model, which in turn is used to study the properties of the candidate LAE. Population III galaxy spectral energy distribution (SED) model fits to the CLASH broadband photometry of the possible LAE provide a slightly better fit than Population I/II models. The best fitted model suggests intrinsic EW(Ly$\alpha$) $\approx 160$ \AA{} after absorption in the interstellar and intergalactic medium. Future spectroscopic observations will examine this prediction as well as shed more light on the morphology of this object, which indicates it may be a merger of two smaller galaxies.

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

While numerical simulations suggest that the strength of the Lyman alpha (Lya) line of star-forming disk galaxies strongly depends on the inclination at which they are observed (i.e. from edge-on to face-on, we expect to see a change from an attenuated Lya line to a strong Lya emission line), recent observations with the Hubble space telescope (HST) have highlighted few low-redshift highly inclined (edge-on) disk galaxies that breaks this trend. We aim to understand how a strong Lya emission line is able to escape from one of those inclined disk galaxies, named Mrk1486 (z=0.0338). For that purpose we used a large set of HST imaging and spectroscopic data to investigate both the ISM structure and the dominant source of Lya radiation inside Mrk1486. Moreover, we used a 3D Monte Carlo Lya radiation transfer code to study the radiative transfer of Lya and UV continuum photons inside a 3D geometry of neutral hydrogen (HI) and dust that models the ISM structure at the galaxy center. The analysis of IFU Halpha spectroscopic data of Mrk1486 indicates the presence of two bipolar galactic winds of HI gas above and bellow the disk plane of Mrk1486. Furthermore, comparing different diagnostic diagrams (such as [OIII]5007/Hbeta versus [OI]6300/Halpha) to photo- and shock-ionization models, we find that the Lya production of Mrk1486 is dominated by photoionization inside the galaxy disk. From this perspective, our numerical simulations succeed in reproducing the strength and spectral shape of the observed Lya line of Mrk1486 by assuming a scenario in which the Lya photons are produced inside the disk, travel along the galactic winds and scatter on cool HI materials toward the observer. As bipolar galactic winds are ubiquitous in star-forming disk galaxies, this mechanism may explain the origin of strong Lya emission lines recently observed from highly inclined galaxies at high-redshift.

We present integral field spectroscopic observations with the Potsdam Multi Aperture Spectrophotometer of all 14 galaxies in the $z\sim 0.1$ Lyman Alpha Reference Sample (LARS). We produce 2D line of sight velocity maps and velocity dispersion maps from the Balmer $\alpha$ (H$\alpha$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$\alpha$ radiation field. We show our kinematic maps spatially registered onto the Hubble Space Telescope H$\alpha$ and Lyman $\alpha$ (Ly$\alpha$) images. Only for individual galaxies a causal connection between spatially resolved H$\alpha$ kinematics and Ly$\alpha$ photometry can be conjectured. However, no general trend can be established for the whole sample. Furthermore, we compute non-parametric global kinematical statistics -- intrinsic velocity dispersion $\sigma_0$, shearing velocity $v_\mathrm{shear}$, and the $v_\mathrm{shear}/\sigma_0$ ratio -- from our kinematic maps. In general LARS galaxies are characterised by high intrinsic velocity dispersions (54\,km\,s$^{-1}$ median) and low shearing velocities (65\,km\,s$^{-1}$ median). $v_\mathrm{shear}/\sigma_0$ values range from 0.5 to 3.2 with an average of 1.5. Noteworthy, five galaxies of the sample are dispersion dominated systems with $v_\mathrm{shear}/\sigma_0 <1$ and are thus kinematically similar to turbulent star forming galaxies seen at high redshift. When linking our kinematical statistics to the global LARS Ly$\alpha$ properties, we find that dispersion dominated systems show higher Ly$\alpha$ equivalent widths and higher Ly$\alpha$ escape fractions than systems with $v_\mathrm{shear}/\sigma_0 > 1$. Our result indicates that turbulence in actively star-forming systems is causally connected to interstellar medium conditions that favour an escape of Ly$\alpha$ radiation.

The leakage of Lyman continuum photons from star forming galaxies is an elusive parameter. When observed, it provides a wealth of information on star formation in galaxies and the geometry of the interstellar medium, and puts constraints on the role of star forming galaxies in the reionization of the universe. H-alpha-selected galaxies at z~2 trace the highest star formation population at the peak of cosmic star formation history, providing a base for directly measuring Lyman continuum escape. Here we present this method, and highlight its benefits as well as caveats. We also use the method on 10 H-alpha emitters in the Chandra Deep Field South at z=2.2, also imaged with the Hubble Space Telescope in the ultraviolet. We find no individual Lyman continuum detections, and our stack puts a 5 sigma upper limit on the average absolute escape fraction of <24%, consistent with similar studies. With future planned observations, the sample sizes would rapidly increase and the method presented here would provide very robust constraints on the escape fraction.

The stellar content of young massive star clusters emit large amounts of Lyman continuum photons and inject momentum into the inter stellar medium (ISM) by the strong stellar winds of the most massive stars in the cluster. When the most massive stars explode as supernovae, large amounts of mechanical energy are injected in the ISM. A detailed study of the ISM around these massive cluster provides insights on the effect of cluster feedback. We present high quality integral field spectroscopy taken with VLT/MUSE of two starburst galaxies: ESO 338-IG04 and Haro 11. Both galaxies contain a significant number of super star clusters. The MUSE data provide us with an unprecedented view of the state and kinematics of the ionized gas in the galaxy allowing us to study the effect of stellar feedback on small and large spatial scales. We present our recent results on studying the ISM state of these two galaxies. The data of both galaxies show that the mechanical and ionization feedback of the super star clusters in the galaxy modify the state and kinematics of the ISM substancially by creating highly ionized bubbles around the cluster, making the central part of the galaxy highly ionized. This shows that the HII regions around the individual clusters are density bounded, allowing the ionizing photons to escape and ionize the ISM further out.

The Ly$\alpha$ line is an important diagnostic for star formation at high redshift, but interpreting its flux and line profile is difficult because of the resonance nature of Ly$\alpha$. Trends between the escape of Ly$\alpha$ photons and dust and properties of the interstellar medium (ISM) have been found, but detailed comparisons between Ly$\alpha$ emission and the properties of the gas in local high-redshift analogs are vital for understanding the relation between Ly$\alpha$ emission and galaxy properties. For the first time, we can directly infer the properties of the ionized gas at the same location and similar spatial scales of the extended Ly$\alpha$ halo around ESO 338-IG04. We obtained VLT/MUSE integral field spectra. We used ionization parameter mapping of the [SII]/[OIII] line ratio and the kinematics of H$\alpha$ to study the ionization state and kinematics of the ISM of ESO338-IG04. The velocity map reveals two outflows. The entire central area of the galaxy is highly ionized by photons leaking from the HII regions around the youngest star clusters. Three highly ionized cones have been identified, of which one is associated with an outflow. We propose a scenario where the outflows are created by mechanical feedback of the older clusters, while the highly ionized gas is caused by the hard ionizing photons emitted by the youngest clusters. A comparison with the Ly$\alpha$ map shows that the (approximately bipolar) asymmetries observed in the Ly$\alpha$ emission are consistent with the base of the outflows detected in H$\alpha$. No clear correlation with the ionization cones is found. The mechanical and ionization feedback of star clusters significantly changes the state of the ISM by creating ionized cones and outflows. The comparison with Ly$\alpha$ suggests that especially the outflows could facilitate the escape of Ly$\alpha$ photons [Abridged].

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

Population III galaxies are predicted to exist at high redshifts and may be rendered sufficiently bright for detection with current telescopes when gravitationally lensed by a foreground galaxy cluster. Population III galaxies that exhibit strong Lya emission should furthermore be identifiable from broadband photometry because of their unusual colors. Here, we report on a search for such objects at z > 6 in the imaging data from the Cluster Lensing And Supernova survey with Hubble (CLASH), covering 25 galaxy clusters in 16 filters. Our selection algorithm returns five singly-imaged candidates with Lya-like color signatures, for which ground-based spectroscopy with current 8-10 m class telescopes should be able to test the predicted strength of the Lya line. None of these five objects have been included in previous CLASH compilations of high-redshift galaxy candidates. However, when large grids of spectral synthesis models are applied to the study of these objects, we find that only two of these candidates are significantly better fitted by Population III models than by more mundane, low-metallicity stellar populations.

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

The Lyman-alpha reference sample (LARS) is a program with the Hubble Space Telescope (HST) that provides a sample of local universe laboratory galaxies in which to study the astrophysics of the visibility and strength of the Lyman-alpha (Lya) line of hydrogen. This article presents an overview of the survey, its selection function and HST imaging observations. The sample was selected from the GALEX+SDSS catalogue at z=0.028-0.19, in order to allow Lya to be captured with combinations of long pass filters in the Solar Blind Channel (SBC) of HST/ACS. In addition, LARS utilises Halpha and Hbeta narrow, and U, B, i broad-band imaging with ACS and WFC3. In order to study galaxies in which large numbers of Lya photons are produced we demanded an Halpha equivalent width > 100{\AA}. The sample of 14 galaxies covers far UV (FUV) luminosities that overlaps with those of high-z Lya emitters and Lyman Break Galaxies (LBGs). We present the reduction steps used to obtain the Lya images, including our LARS eXtraction software (LaXs) which utilises pixel-by-pixel spectral synthesis fitting to subtract the continuum at Lya. We demonstrate that the use of SBC long pass filters dramatically increase the signal to noise compared to the nominal Lya filter in SBC. To exemplify, we present results for Mrk 259. This irregular galaxy shows extended but strongly asymmetric Lya emission. Spectroscopy from HST/COS show a moderate outflow in the neutral ISM and an asymmetric Lya profile. Radiative transfer modeling is able to reproduce the Lya line profile, confirming the presence of an outflow. The integrated photometry provides a Lya luminosity of 1.3E42 erg/s an equivalent width W(Lya) =45{\AA} and a FUV absolute magnitude -19.2. Mrk 259 would hence be detectable in high-z Lya and LBG surveys. The total Lya escape fraction is 12%, similar to that observed at z>4. (shortened)

We present new H I imaging and spectroscopy of the 14 UV-selected star-forming galaxies in the Lyman Alpha Reference Sample (LARS), aimed for a detailed study of the processes governing the production, propagation, and escape of Ly$\alpha$ photons. New H I spectroscopy, obtained with the 100m Green Bank Telescope (GBT), robustly detects the H I spectral line in 11 of the 14 observed LARS galaxies (although the profiles of two of the galaxies are likely confused by other sources within the GBT beam); the three highest redshift galaxies are not detected at our current sensitivity limits. The GBT profiles are used to derive fundamental H I line properties of the LARS galaxies. We also present new pilot H I spectral line imaging of 5 of the LARS galaxies obtained with the Karl G. Jansky Very Large Array (VLA). This imaging localizes the H I gas and provides a measurement of the total H I mass in each galaxy. In one system, LARS 03 (UGC 8335 or Arp 238), VLA observations reveal an enormous tidal structure that extends over 160 kpc from the main interacting systems and that contains $>$10$^9$ M$_{\odot}$ of H I. We compare various H I properties with global Ly$\alpha$ quantities derived from HST measurements. The measurements of the Ly$\alpha$ escape fraction are coupled with the new direct measurements of H I mass and significantly disturbed H I velocities. Our robustly detected sample reveals that both total H I mass and linewidth are tentatively correlated with key Ly$\alpha$ tracers. Further, on global scales, these data support a complex coupling between Ly$\alpha$ propagation and the H I properties of the surrounding medium.

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

We estimate the fraction of core-collapse supernovae (CCSNe) that remain undetected by optical SN searches due to obscuration by large amounts of dust in their host galaxies. This effect is especially important in luminous and ultraluminous infrared galaxies, which are locally rare but dominate the star formation at redshifts of z~1-2. We perform a detailed investigation of the SN activity in the nearby luminous infrared galaxy Arp 299 and estimate that up to 83% of the SNe in Arp 299 and in similar galaxies in the local Universe are missed by observations at optical wavelengths. For rest-frame optical surveys we find the fraction of SNe missed due to high dust extinction to increase from the average local value of ~19% to ~38% at z~1.2 and then stay roughly constant up to z~2. It is therefore crucial to take into account the effects of obscuration by dust when determining SN rates at high redshift and when predicting the number of CCSNe detectable by future high-z surveys such as LSST, JWST, and Euclid. For a sample of nearby CCSNe (distances 6-15 Mpc) detected during the last 12 yr, we find a lower limit for the local CCSN rate of 1.5 +0.4/-0.3 x 10^-4 yr^-1 Mpc^-3, consistent with that expected from the star formation rate. Even closer, at distances less than ~6 Mpc, we find a significant increase in the CCSN rate, indicating a local overdensity of star formation caused by a small number of galaxies that have each hosted multiple SNe.

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

We report the discovery of two consecutive supernovae (SNe), 2010cu and 2011hi, located at 0.37" (180 pc) and 0.79" (380 pc) projected distance respectively from the centre of the K-band nucleus of the luminous infrared galaxy IC 883. The SNe were discovered in an ongoing near-infrared K-band search for core-collapse SNe in such galaxies using the ALTAIR/NIRI adaptive optics system with laser guide star at the Gemini-North Telescope. These are thus the closest SNe yet discovered to a LIRG nucleus in optical or near-infrared wavelengths. The near-infrared light curves and colours of both SNe are consistent with core-collapse events. Both SNe seem to suffer from relatively low host galaxy extinction suggesting that regardless of their low projected galactocentric distances, they are not deeply buried in the nuclear regions of the host galaxy.

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

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 the work presented in this paper is to test and optimise supernova detection methods based on the optimal image subtraction technique. The main focus is on applying the detection methods to wide field supernova imaging surveys and in particular to the Stockholm VIMOS Supernova Survey (SVISS). We have constructed a supernova detection pipeline for imaging surveys. The core of the pipeline is image subtraction using the ISIS 2.2 package. Using real data from the SVISS we simulate supernovae in the images, both inside and outside galaxies. The detection pipeline is then run on the simulated frames and the effects of image quality and subtraction parameters on the detection efficiency and photometric accuracy are studied. The pipeline allows efficient detection of faint supernovae in the deep imaging data. It also allows controlling and correcting for possible systematic effects in the SN detection and photometry. We find such a systematic effect in the form of a small systematic flux offset remaining at the positions of galaxies in the subtracted frames. This offset will not only affect the photometric accuracy of the survey, but also the detection efficiencies. Our study has shown that ISIS 2.2 works well for the SVISS data. We have found that the detection efficiency and photometric accuracy of the survey are affected by the stamp selection for the image subtraction and by host galaxy brightness. With our tools the subtraction results can be further optimised, any systematic effects can be controlled and photometric errors estimated, which is very important for the SVISS, as well as for future SN searches based on large imaging surveys such as Pan-STARRS and LSST.

An extensive set of optical and NIR photometry and low-resolution spectra the Type Ia supernova (SN Ia) 2003du was obtained using a number of facilities. The observations started 13 days before B-band maximum light and continued for 480 days with exceptionally good time sampling. The optical photometry was calibrated through the S-correction technique. The UBVRIJHK light curves and the color indices of SN 2003du closely resemble those of normal SNe Ia. SN 2003du reached a B-band maximum of 13.49 (+/-0.02) mag on JD2452766.38 (+/-0.5). We derive a B-band stretch parameter of 0.988 (+/-0.003), which corresponds to dM15=1.02 (+/-0.05), indicative of a SN Ia of standard luminosity. The reddening in the host galaxy was estimated by three methods, and was consistently found to be negligible. We estimate a distance modulus of 32.79 (+/-0.15) mag to the host galaxy, UGC 9391. The peak UVOIR bolometric luminosity of 1.35(+/-0.20) 10^43 erg/s and Arnett's rule implies that M(Ni56)=0.68 (+/-0.14)M_sun of Ni56 was synthesized during the explosion. Modeling of the UVOIR bolometric light curve also indicates M(Ni56) in the range 0.6-0.8 M_sun. The spectral evolution of SN 2003du at both optical and NIR wavelengths also closely resembles normal SNe Ia. In particular, the Si II ratio at maximum R(Si II)=0.22 (+/-0.02) and the time evolution of the blueshift velocities of the absorption line minima are typical. The pre-maximum spectra of SN 2003du showed conspicuous high-velocity features in the Ca II H&K doublet and infrared triplet, and possibly in Si II 6355, lines. We compare the time evolution of the profiles of these lines with other well-observed SNe Ia and we suggest that the peculiar pre-maximum evolution of Si II 6355 line in many SNe Ia is due to the presence of two blended absorption components.

We have obtained deep optical images with the Very Large Telescope at ESO of the first well-localized short-duration gamma-ray burst, GRB 050509B. From V and R imaging, initiated ~2 days after the GRB trigger and lasting up to three weeks, we detect no variable object inside the small Swift/XRT X-ray error circle down to 2sigma limits of V = 26.5 and R = 25.1. The X-ray error circle includes a giant elliptical galaxy at z = 0.225, which has been proposed as the likely host of this GRB. Our limits indicate that if the GRB originated at z = 0.225, any supernova-like event accompanying the GRB would have to be over 100 times fainter than normal Type Ia SNe or Type Ic hypernovae, 5 times fainter than the faintest known Ia or Ic SNe, and fainter than the faintest known Type II SNe. Moreover, we use the optical limits to constrain the energetics of the GRB outflow. Simple models indicate that, unless the intrinsic energy in the outflow from GRB 050509B was << 10^51 erg, there was very little radioactive material with efficient decay timescales for generating a large luminosity. These limits strongly constrain progenitor models for this short GRB.