In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

High-resolution multi-wavelength photometry is crucial to explore the spatial distribution of star formation in galaxies and understand how these evolve. To this aim, in this paper we exploit the deep, multi-wavelength Hubble Space Telescope (HST) data available in the central parts of the GOODS fields and study the distribution of star formation activity and mass in galaxies located at different positions with respect to the Main Sequence (MS) of star-forming galaxies. Our sample consists of galaxies with stellar mass $\geq 10^{9.5} M_{\odot}$ in the redshift range 0.2 $ \leq z \leq 1.2$. Exploiting 10-band photometry from the UV to the near-infrared at HST resolution, we derive spatially resolved maps of galaxies properties, such as stellar mass and star formation rate and specific star formation rate, with a resolution of $\sim 0.16$ arcsec. We find that the star formation activity is centrally enhanced in galaxies above the MS and centrally suppressed below the MS, with quiescent galaxies (1 dex below the MS) characterised by the highest suppression. The sSFR in the outer region does not show systematic trends of enhancement or suppression above or below the MS. The distribution of mass in MS galaxies indicates that bulges are growing when galaxies are still on the MS relation. Galaxies below the MS are more bulge-dominated with respect to MS counterparts at fixed stellar mass, while galaxies in the upper envelope are more extended and have S\'ersic indexes that are always smaller than or comparable to their MS counterparts. The suppression of star formation activity in the central region of galaxies below the MS hints at \textit{inside-out} quenching, as star formation is still ongoing in the outer regions.

We present the Hubble Deep UV Legacy Survey (HDUV), a 132 orbit imaging program with the WFC3/UVIS camera onboard the Hubble Space Telescope (HST). The HDUV extends and builds on the few previous UV imaging surveys in the two GOODS/CANDELS-Deep fields to provide deep images over a total area of ~100 arcmin2 in the two filters F275W and F336W. Our release also includes all the F275W imaging data taken by the CANDELS survey, which were aligned using a novel approach and combined with the HDUV survey data. By reaching depths of 27.5-28.0 mag (5sigma, in 0.4" apertures), these are the deepest high-resolution UV data over such a large area taken to date. Such unique UV imaging enables a wide range of science by the community. Among the main goals of the HDUV survey are: (1) provide a complete sample of faint star-forming galaxies at z~1-3, (2) constrain the ionizing photon escape fraction from galaxies at z~2-3, and (3) track the build-up of bulges and the disappearance of clumpy disk galaxies through reliable internal stellar population properties at sub-kpc resolution out to z~3. The addition of the HDUV data further enhances the legacy value of the two GOODS/CANDELS-Deep fields, which now include deep 11-band HST imaging as well as very deep ancillary data from X-ray to radio, enabling unique multi-wavelength studies. Here, we provide an overview of the survey design, describe the data reduction, and highlight a few basic analyses on the images which are released to the community as high level science products via the Mikulski Archive for Space Telescopes (MAST).

With the Hubble Frontier Fields program, gravitational lensing has provided a powerful way to extend the study of the ultraviolet luminosity function (LF) of galaxies at $z \sim 6$ down to unprecedented magnitude limits. At the same time, significant discrepancies between different studies were found at the very faint end of the LF. In an attempt to understand such disagreements, we present a comprehensive assessment of the uncertainties associated with the lensing models and the size distribution of galaxies. We use end-to-end simulations from the source plane to the final LF that account for all lensing effects and systematic uncertainties by comparing several mass models. In addition to the size distribution, the choice of lens model leads to large differences at magnitudes fainter than $M_{UV} = -15~$ AB mag, where the magnification factor becomes highly uncertain. We perform MCMC simulations that include all these uncertainties at the individual galaxy level to compute the final LF, allowing, in particular, a crossover between magnitude bins. The best LF fit, using a modified Schechter function that allows for a turnover at faint magnitudes, gives a faint-end slope of $\alpha = -2.01_{-0.14}^{+0.12}$, a curvature parameter of $\beta = 0.48_{-0.25}^{+0.49}$, and a turnover magnitude of $M_{T} = -14.93_{-0.52}^{+0.61}$. Most importantly our procedure shows that robust constraints on the LF at magnitudes fainter than $M_{UV} = -15~$ AB remain unrealistic. More accurate lens modeling and future observations of lensing clusters with the James Webb Space Telescope can reliably extend the UV LF to fainter magnitudes.

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.

The first billion years of the Universe is a pivotal time: stars, black holes (BHs) and galaxies form and assemble, sowing the seeds of galaxies as we know them today. Detecting, identifying and understand- ing the first galaxies and BHs is one of the current observational and theoretical challenges in galaxy formation. In this paper we present a population synthesis model aimed at galaxies, BHs and Active Galactic Nuclei (AGNs) at high redshift. The model builds a population based on empirical relations. Galaxies are characterized by a spectral energy distribution determined by age and metallicity, and AGNs by a spectral energy distribution determined by BH mass and accretion rate. We validate the model against observational constraints, and then predict properties of galaxies and AGN in other wavelength and/or luminosity ranges, estimating the contamination of stellar populations (normal stars and high-mass X-ray binaries) for AGN searches from the infrared to X-rays, and vice-versa for galaxy searches. For high-redshift galaxies, with stellar ages < 1 Gyr, we find that disentangling stellar and AGN emission is challenging at restframe UV/optical wavelengths, while high-mass X-ray binaries become more important sources of confusion in X-rays. We propose a color-color selection in JWST bands to separate AGN vs star-dominated galaxies in photometric observations. We also esti- mate the AGN contribution, with respect to massive, hot, metal-poor stars, at driving high ionization lines, such as C IV and He II. Finally, we test the influence of the minimum BH mass and occupa- tion fraction of BHs in low mass galaxies on the restframe UV/near-IR and X-ray AGN luminosity function.

We provide the first observational constraints on the sizes of the faintest galaxies lensed by the Hubble Frontier Fields (HFF) clusters. Ionizing radiation from faint galaxies likely drives cosmic reionization, and the HFF initiative provides a key opportunity to find such galaxies. Yet, we cannot really assess their ionizing emissivity without a robust measurement of their sizes, since this is key to quantifying both their prevalence and the faint-end slope to the UV luminosity function. Here we provide the first such size constraints with 2 new techniques. The first utilizes the fact that the detectability of highly-magnified galaxies as a function of shear is very dependent on a galaxy's size. Only the most compact galaxies will remain detectable in regions of high shear (vs. a larger detectable size range for low shear), a phenomenon we carefully quantify using simulations. Remarkably, however, no correlation is found between the surface density of faint galaxies and the predicted shear, using 87 faint high-magnification mu>10 z~2-8 galaxies seen behind the first 4 HFF clusters. This can only be the case if such faint (~-15 mag) galaxies have significantly smaller sizes than luminous galaxies. We constrain their half-light radii to be <~30 mas (<160-240 pc). As a 2nd size probe, we rotate and stack 26 faint high-magnification sources along the major shear axis. Less elongation is found than even for objects with an intrinsic half-light radius of 10 mas. Together these results indicate that extremely faint z~2-8 galaxies have near point-source profiles in the HFF dataset (half-light radii conservatively <30 mas and likely 5-10 mas). These results suggest smaller completeness corrections and hence much lower volume densities for faint z~2-8 galaxies and shallower faint-end slopes than have been derived in many recent studies (by factors of ~2-3 and by dalpha>~0.1-0.3).

We use a newly assembled large sample of 3,545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z=1.5-2.5 to constrain the relationship between UV slope (beta) and dust attenuation (L(IR)/L(UV)=IRX). Our sample benefits from the combination of deep Hubble WFC3/UVIS photometry from the Hubble Deep UV (HDUV) Legacy survey and existing photometric data compiled in the 3D-HST survey, and extends the range of UV luminosity and beta probed in previous UV-selected samples. IRX is measured using stacks of deep Herschel/PACS 100 and 160 micron data, and the results are compared with predictions of the IRX-beta relation for different assumptions of the stellar population model and obscuration curve. We find that z=1.5-2.5 galaxies have an IRX-beta relation that is consistent with the predictions for an SMC extinction curve if we invoke sub-solar metallicity models that are currently favored for high-redshift galaxies, while the commonly assumed starburst attenuation curve over-predicts the IRX at a given beta by a factor of ~3. The IRX of high-mass (M*>10^9.75 Msun) galaxies is a factor of >4 larger than that of low-mass galaxies, lending support for the use of stellar mass as a proxy for attenuation. The commonly observed trend of fainter galaxies having bluer beta may simply reflect bluer intrinsic UV slopes for such galaxies, rather than lower obscurations. The IRX-beta for young/low-mass galaxies implies a dust curve that is steeper than the SMC, suggesting a lower attenuation at a given beta relative to older/more massive galaxies. The lower attenuations and higher ionizing photon output implied by low metallicity stellar population models point to Lyman continuum production efficiencies, xi_ion, that may be elevated by a factor of ~2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z~2. [Abridged]

We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations.

We present the results of a systematic search for Lyman-alpha emitters (LAEs) at $6 \lesssim z \lesssim 7.6$ using the HST WFC3 Infrared Spectroscopic Parallel (WISP) Survey. Our total volume over this redshift range is $\sim 8 \times10^5$ Mpc$^3$, comparable to many of the narrowband surveys despite their larger area coverage. We find two LAEs at $z=6.38$ and $6.44$ with line luminosities of L$_{\mathrm{Ly}\alpha} \sim 4.7 \times 10^{43}$ erg s$^{-1}$, putting them among the brightest LAEs discovered at these redshifts. Taking advantage of the broad spectral coverage of WISP, we are able to rule out almost all lower-redshift contaminants. The WISP LAEs have a high number density of $7.7\times10^{-6}$ Mpc$^{-3}$. We argue that the LAEs reside in Mpc-scale ionized bubbles that allow the Lyman-alpha photons to redshift out of resonance before encountering the neutral IGM. We discuss possible ionizing sources and conclude that the observed LAEs alone are not sufficient to ionize the bubbles.

We present six galaxies at z~2 that show evidence of Lyman continuum (LyC) emission based on the newly acquired UV imaging of the Hubble Deep UV legacy survey (HDUV) conducted with the WFC3/UVIS camera on the Hubble Space Telescope (HST). At the redshift of these sources, the HDUV F275W images partially probe the ionizing continuum. By exploiting the HST multi-wavelength data available in the HDUV/GOODS fields, models of the UV spectral energy distributions, and detailed Monte Carlo simulations of the intergalactic medium absorption, we estimate the absolute ionizing photon escape fractions of these galaxies to be very high -- typically >60% (>13% for all sources at 90% likelihood). Our findings are in broad agreement with previous studies that found only a small fraction of galaxies to show high escape fraction. These six galaxies comprise the largest sample yet of LyC leaking candidates at z~2 whose inferred LyC flux has been cleanly observed at HST resolution. While three of our six candidates show evidence of hosting an active galactic nucleus (AGN), two of these are heavily obscured and their LyC emission appears to originate from star-forming regions rather than the central nucleus. This suggests an AGN-aided pathway for LyC escape from these sources. Extensive multi-wavelength data in the GOODS fields, especially the near-IR grism spectra from the 3D-HST survey, enable us to study the candidates in detail and tentatively test some recently proposed indirect methods to probe LyC leakage -- namely, the [OIII]/[OII] line ratio and the H$\beta-$UV slope diagram. High-resolution spectroscopic followup of our candidates will help constrain such indirect methods which are our only hope of studying $f_{esc}$ at z~5-9 in the fast-approaching era of the James Webb Space Telescope.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and could thus be a potentially powerful tool to trace cosmic star formation. However, especially at low redshifts (z < 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the impact of the environment on GRB formation. We study host galaxy spectra of the Swift/BAT6 complete sample of 14 z < 1 bright LGRBs. We use the detected nebular emission lines to measure the dust extinction, star formation rate (SFR) and nebular metallicity (Z) of the hosts and supplement the data set with previously measured stellar masses M$_{\star}$. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M$_{\star}$ relations) are compared to samples of field star-forming galaxies.We find that LGRB hosts at z < 1 have on average lower SFRs than if they were direct star-formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to Z $\sim 8.4-8.5$, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are not inconsistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of increased incidence of starburst galaxies in the Swift/BAT6 z < 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevent us from reaching more quantitative conclusions.

We analyze how passive galaxies at z $\sim$ 1.5 populate the mass-size plane as a function of their stellar age, to understand if the observed size growth with time can be explained with the appearance of larger quenched galaxies at lower redshift. We use a sample of 32 passive galaxies extracted from the Wide Field Camera 3 Infrared Spectroscopic Parallel (WISP) survey with spectroscopic redshift 1.3 $\lesssim$ z $\lesssim$ 2.05, specific star-formation rates lower than 0.01 Gyr$^{-1}$, and stellar masses above 4.5 $\times$ 10$^{10}$ M$_\odot$. All galaxies have spectrally determined stellar ages from fitting of their rest-frame optical spectra and photometry with stellar population models. When dividing our sample into young (age $\leq$ 2.1 Gyr) and old (age $>$ 2.1 Gyr) galaxies we do not find a significant trend in the distributions of the difference between the observed radius and the one predicted by the mass-size relation. This result indicates that the relation between the galaxy age and its distance from the mass-size relation, if it exists, is rather shallow, with a slope alpha $\gtrsim$ -0.6. At face value, this finding suggests that multiple dry and/or wet minor mergers, rather than the appearance of newly quenched galaxies, are mainly responsible for the observed time evolution of the mass-size relation in passive galaxies.

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

[abridged] We present a strong-lensing analysis of MACSJ0717.5+3745, based on the full depth of the Hubble Frontier Field (HFF) observations, which brings the number of multiply imaged systems to 61, ten of which are spectroscopically confirmed. The total number of images comprised in these systems rises to 165. Our analysis uses a parametric mass reconstruction technique, as implemented in the Lenstool software, to constrain a mass distribution composed of four large-scale mass components + galaxy-scale perturbers. We find a superposition of cored isothermal mass components to provide a good fit to the observational constraints, resulting in a very shallow mass distribution for the smooth (large-scale) component. Given the implications of such a flat mass profile, we investigate whether a model composed of "peaky" non-cored mass components can also reproduce the observational constraints. We find that such a non-cored mass model reproduces the observational constraints equally well. Although the total mass distributions of both models are consistent, as well as the integrated two dimensional mass profiles, we find that the smooth and the galaxy-scale components are very different. We conclude that, even in the HFF era, the generic degeneracy between smooth and galaxy-scale components is not broken, in particular in such a complex galaxy cluster. Consequently, insights into the mass distribution of MACS J0717 remain limited, underlining the need for additional probes beyond strong lensing. Our findings also have implications for estimates of the lensing magnification: we show that the amplification difference between the two models is larger than the error associated with either model. This uncertainty decreases the area of the image plane where we can reliably study the high-redshift Universe by 50 to 70%.

We present a high-precision mass model of the galaxy cluster MACSJ1149.6+2223, based on a strong-gravitational-lensing analysis of Hubble Space Telescope Frontier Fields (HFF) imaging data and spectroscopic follow-up with Gemini/GMOS and VLT/MUSE. Our model includes 12 new multiply imaged galaxies, bringing the total to 22, comprised of 65 individual lensed images. Unlike the first two HFF clusters, Abell 2744 and MACSJ0416.1-2403, MACSJ1149 does not reveal as many multiple images in the HFF data. Using the Lenstool software package and the new sets of multiple images, we model the cluster with several cluster-scale dark-matter halos and additional galaxy-scale halos for the cluster members. Consistent with previous analyses, we find the system to be complex, composed of five cluster-scale halos. Their spatial distribution and lower mass, however, makes MACSJ1149 a less powerful lens. Our best-fit model predicts image positions with an RMS of 0.91". We measure the total projected mass inside a 200~kpc aperture as ($1.840\pm 0.006$)$\times 10^{14}$M$_{\odot}$, thus reaching again 1% precision, following our previous HFF analyses of MACSJ0416.1-2403 and Abell 2744. In light of the discovery of the first resolved quadruply lensed supernova, SN Refsdal, in one of the multiply imaged galaxies identified in MACSJ1149, we use our revised mass model to investigate the time delays and predict the rise of the next image between November 2015 and January 2016.

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 use deep Hubble Space Telescope imaging of the Frontier Fields to accurately measure the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range $z \sim 6-8$. We combine observations in three lensing clusters A2744, MACS0416, MACS0717 and their associated parallels fields to select high-redshift dropout candidates. We use the latest lensing models to estimate the flux magnification and the effective survey volume in combination with completeness simulations performed in the source plane. We report the detection of 227 galaxy candidates at $z=6-7$ and 25 candidates at $z \sim 8$. While the total survey area is about 4 arcmin$^{2}$ in each parallel field, it drops to about 0.6 to 1 arcmin$^{2}$ in the cluster core fields because of the strong lensing. We compute the UV luminosity function at $z \sim 7$ using the combined galaxy sample and perform Monte Carlo simulations to determine the best fit Schechter parameters. We are able to reliably constrain the LF down to an absolute magnitude of $M_{UV}=-15.25$, which corresponds to 0.005$L^{\star}$. More importantly, we find that the faint-end slope remains steep down to this magnitude limit with $\alpha=-2.04_{-0.17}^{+0.13}$. Our results confirm the most recent results in deep blank fields but extend the LF measurements more than two magnitudes deeper. The UV LF at $z \sim 8$ is not very well constrained below $M_{UV}=-18$ due to the small number statistics and incompleteness uncertainties. To assess the contribution of galaxies to cosmic reionization we derive the UV luminosity density at $z\sim7$ by integrating the UV LF down to an observationally constrained limit of $M_{UV} = -15$. We show that our determination of Log($\rho_{UV}$)=$26.2\pm0.13$ (erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$) can be sufficient to maintain the IGM ionized.

We present a high-precision mass model of galaxy cluster Abell 2744, based on a strong-gravitational-lensing analysis of the \emph{Hubble Space Telescope Frontier Fields} (HFF) imaging data, which now include both \emph{Advanced Camera for Surveys} and \emph{Wide-Field Camera 3} observations to the final depth. Taking advantage of the unprecedented depth of the visible and near-infrared data, we identify 34 new multiply imaged galaxies, bringing the total to 61, comprising 181 individual lensed images. In the process, we correct previous erroneous identifications and positions of multiple systems in the northern part of the cluster core. With the \textsc{Lenstool} software and the new sets of multiple images, we model the cluster using two cluster-scale dark matter halos plus galaxy-scale halos for the cluster members. Our best-fit model predicts image positions with an \emph{RMS} error of 0.69$\arcsec$, which constitutes an improvement by almost a factor of two over previous parametric models of this cluster. We measure the total projected mass inside a 200~kpc aperture as ($2.162\pm 0.005$)$\times 10^{14}M_{\odot}$, thus reaching 1\% level precision for the second time, following the recent HFF measurement of MACSJ0416.1-2403. Importantly, the higher quality of the mass model translates into an overall improvement by a factor of 4 of the derived magnification factor. % for the high-redshift lensed background galaxies. Together with our previous HFF gravitational lensing analysis, this work demonstrates that the HFF data enables high-precision mass measurements for massive galaxy clusters and the derivation of robust magnification maps to probe the early Universe.

We present photometry and derived redshifts from up to eleven bandpasses for 9927 galaxies in the Hubble Ultra Deep field (UDF), covering an observed wavelength range from the near-ultraviolet (NUV) to the near-infrared (NIR) with Hubble Space Telescope observations. Our Wide Field Camera 3 (WFC3)/UV F225W, F275W, and F336W image mosaics from the ultra-violet UDF (UVUDF) imaging campaign are newly calibrated to correct for charge transfer inefficiency, and use new dark calibrations to minimize background gradients and pattern noise. Our NIR WFC3/IR image mosaics combine the imaging from the UDF09 and UDF12 campaigns with CANDELS data to provide NIR coverage for the entire UDF field of view. We use aperture-matched point-spread function corrected photometry to measure photometric redshifts in the UDF, sampling both the Lyman break and Balmer break of galaxies at z~0.8-3.4, and one of the breaks over the rest of the redshift range. Our comparison of these results with a compilation of robust spectroscopic redshifts shows an improvement in the galaxy photometric redshifts by a factor of two in scatter and a factor three in outlier fraction over previous UDF catalogs. The inclusion of the new NUV data is responsible for a factor of two decrease in the outlier fraction compared to redshifts determined from only the optical and NIR data, and improves the scatter at z<0.5 and at z>2. The panchromatic coverage of the UDF from the NUV through the NIR yields robust photometric redshifts of the UDF, with the lowest outlier fraction available.

We present high-resolution far-UV spectroscopy of the 14 galaxies of the Lyman Alpha Reference Sample; a sample of strongly star-forming galaxies at low redshifts ($0.028 < z < 0.18$). We compare the derived properties to global properties derived from multi band imaging and 21 cm HI interferometry and single dish observations, as well as archival optical SDSS spectra. Besides the Lyman $\alpha$ line, the spectra contain a number of metal absorption features allowing us to probe the kinematics of the neutral ISM and evaluate the optical depth and and covering fraction of the neutral medium as a function of line-of-sight velocity. Furthermore, we show how this, in combination with precise determination of systemic velocity and good Ly$\alpha$ spectra, can be used to distinguish a model in which separate clumps together fully cover the background source, from the "picket fence" model named by Heckman et al. (2011). We find that no one single effect dominates in governing Ly$\alpha$ radiative transfer and escape. Ly$\alpha$ escape in our sample coincides with a maximum velocity-binned covering fraction of $\lesssim 0.9$ and bulk outflow velocities of $\gtrsim 50$ km s$^{-1}$, although a number of galaxies show these characteristics and yet little or no Ly$\alpha$ escape. We find that Ly$\alpha$ peak velocities, where available, are not consistent with a strong backscattered component, but rather with a simpler model of an intrinsic emission line overlaid by a blueshifted absorption profile from the outflowing wind. Finally, we find a strong anticorrelation between H$\alpha$ equivalent width and maximum velocity-binned covering factor, and propose a heuristic explanatory model.

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

Exploiting the power of gravitational lensing, the Hubble Frontier Fields (HFF) program aims at observing six massive galaxy clusters to explore the distant Universe far beyond the depth limits of blank field surveys. Using the complete Hubble Space Telescope observations of the first HFF cluster Abell 2744, we report the detection of 50 galaxy candidates at $z \sim 7$ and eight candidates at $z \sim 8$ in a total survey area of 0.96 arcmin$^{2}$ in the source plane. Three of these galaxies are multiply-imaged by the lensing cluster. Using an updated model of the mass distribution in the cluster we were able to calculate the magnification factor and the effective survey volume for each galaxy in order to compute the ultraviolet galaxy luminosity function at both redshifts 7 and 8. Our new measurements extend the $z \sim 7$ UV LF down to an absolute magnitude of $M_{UV} \sim -15.5$. We find a characteristic magnitude of $M^{\star}_{UV}=-20.63^{+0.69}_{-0.56}$ mag and a faint-end slope $\alpha = -1.88^{+0.17}_{-0.20}$ close to previous determinations in blank fields. We show here for the first time that this slope remains steep down to very faint luminosities of 0.01$L^{\star}$. Although prone to large uncertainties, our results at $z \sim 8$ also seem to confirm a steep faint-end slope below 0.1$L^{\star}$. The HFF program is therefore providing an extremely efficient way to study the faintest galaxy populations at $z > 7$ that would otherwise be inaccessible with current instrumentation. The full sample of six galaxy clusters will provide yet better constraints on the build-up of galaxies at early epochs and their contribution to cosmic reionization.

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)

Extending over three Hubble Space Telescope (HST) cycles, the Hubble Frontier Fields (HFF) initiative constitutes the largest commitment ever of HST time to the exploration of the distant Universe via gravitational lensing by massive galaxy clusters. We here present models of the mass distribution in the six HFF cluster lenses, derived from a joint strong- and weak-lensing analysis anchored by a total of 88 multiple-image systems identified in existing HST data. The resulting maps of the projected mass distribution and of the gravitational magnification effectively calibrate the HFF clusters as gravitational telescopes. Allowing the computation of search areas in the source plane, these maps are provided to the community to facilitate the exploitation of forthcoming HFF data for quantitative studies of the gravitationally lensed population of background galaxies. Our models of the gravitational magnification afforded by the HFF clusters allow us to quantify the lensing-induced boost in sensitivity over blank-field observations and predict that galaxies at $z>10$ and as faint as m(AB)=32 will be detectable, up to 2 magnitudes fainter than the limit of the Hubble Ultra Deep Field.

Near infrared slitless spectroscopy with the Wide Field Camera 3, onboard the Hubble Space Telescope, offers a unique opportunity to study low-mass galaxy populations at high-redshift ($z\sim$1-2). While most high$-z$ surveys are biased towards massive galaxies, we are able to select sources via their emission lines that have very-faint continua. We investigate the star formation rate (SFR)-stellar mass ($M_{\star}$) relation for about 1000 emission-line galaxies identified over a wide redshift range of $0.3 \lesssim z \lesssim 2.3$. We use the H$_{\alpha}$ emission as an accurate SFR indicator and correct the broadband photometry for the strong nebular contribution to derive accurate stellar masses down to $M_{\star} \sim 10^{7} M_{\odot}$. We focus here on a subsample of galaxies that show extremely strong emission lines (EELGs) with rest-frame equivalent widths ranging from 200 to 1500 \AA. This population consists of outliers to the normal SFR-$M_{\star}$ sequence with much higher specific SFRs ($> 10$ Gyr$^{-1}$). While on-sequence galaxies follow a continuous star formation process, EELGs are thought to be caught during an extreme burst of star formation that can double their stellar mass in less than $100$ Myr. The contribution of starbursts to the total star formation density appears to be larger than what has been reported for more massive galaxies in previous studies. In the complete mass range $8.2 <$ log($M_{\star}/M_{\odot}$) $< 10$ and a SFR lower completeness limit of about 2 $M_{\odot}$ yr$^{-1}$ (10 $M_{\odot}$ yr$^{-1}$) at $z\sim1$ ($z \sim 2$), we find that starbursts having EW$_{rest}$(H$_{\alpha}$)$>$ 300, 200, and 100 A contribute up to $\sim13$, 18, and 34 %, respectively, to the total SFR of emission-line selected sample at $z\sim1-2$. The comparison with samples of massive galaxies shows an increase in the contribution of starbursts towards lower masses.

We present a high-precision mass model of the galaxy cluster MACSJ0416.1-2403, based on a strong-gravitational-lensing analysis of the recently acquired Hubble Space Telescope Frontier Fields (HFF) imaging data. Taking advantage of the unprecedented depth provided by HST/ACS observations in three passbands, we identify 51 new multiply imaged galaxies, quadrupling the previous census and bringing the grand total to 68, comprising 194 individual lensed images. Having selected a subset of the 57 most securely identified multiply imaged galaxies, we use the Lenstool software package to constrain a lens model comprised of two cluster-scale dark-matter halos and 98 galaxy-scale halos. Our best-fit model predicts image positions with an $RMS$ error of 0.68'', which constitutes an improvement of almost a factor of two over previous, pre-HFF models of this cluster. We find the total projected mass inside a 200~kpc aperture to be $(1.60\pm0.01)\times 10^{14}\ M_\odot$, a measurement that offers a three-fold improvement in precision, reaching the percent level for the first time in any cluster. Finally, we quantify the increase in precision of the derived gravitational magnification of high-redshift galaxies and find an improvement by a factor of $\sim$2.5 in the statistical uncertainty. Our findings impressively confirm that HFF imaging has indeed opened the domain of high-precision mass measurements for massive clusters of galaxies.

Using three magnified Type Ia supernovae (SNe Ia) detected behind CLASH clusters, we perform a first pilot study to see whether standardizable candles can be used to calibrate cluster mass maps created from strong lensing observations. Such calibrations will be crucial when next generation HST cluster surveys (e.g. FRONTIER) provide magnification maps that will, in turn, form the basis for the exploration of the high redshift Universe. We classify SNe using combined photometric and spectroscopic observations, finding two of the three to be clearly of type SN Ia and the third probable. The SNe exhibit significant amplification, up to a factor of 1.7 at $\sim5\sigma$ significance (SN-L2). We conducted this as a blind study to avoid fine tuning of parameters, finding a mean amplification difference between SNe and the cluster lensing models of $0.09 \pm 0.09^{stat} \pm 0.05^{sys}$ mag. This impressive agreement suggests no tension between cluster mass models and high redshift standardized SNe Ia. However, the measured statistical dispersion of $\sigma_{\mu}=0.21$ mag appeared large compared to the dispersion expected based on statistical uncertainties ($0.14$). Further work with the supernova and cluster lensing models, post unblinding, reduced the measured dispersion to $\sigma_{\mu}=0.12$. An explicit choice should thus be made as to whether SNe are used unblinded to improve the model, or blinded to test the model. As the lensed SN samples grow larger, this technique will allow improved constraints on assumptions regarding e.g. the structure of the dark matter halo.

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

The Hubble Frontier Fields (HFF) program combines the capabilities of the Hubble Space Telescope (HST) with the gravitational lensing of massive galaxy clusters to probe the distant Universe to an unprecedented depth. Here, we present the results of the first combined HST and Spitzer observations of the cluster Abell 2744. We combine the full near-infrared data with ancillary optical images to search for gravitationally lensed high-redshift (z > 6) galaxies. We report the detection of 15 I814-dropout candidates at z ~ 6-7 and one Y105-dropout at z ~ 8 in a total survey area of 1.43 arcmin^2 in the source plane. The predictions of our lens model allow us to also identify five multiply-imaged systems lying at redshifts between z ~ 6 and z ~ 8. Thanks to constraints from the mass distribution in the cluster, we were able to estimate the effective survey volume corrected for completeness and magnification effects. This was in turn used to estimate the rest-frame ultraviolet luminosity function (LF) at z ~ 6-8. Our LF results are generally in agreement with the most recent blank field estimates, confirming the feasibility of surveys through lensing clusters. Although based on a shallower observations than what will be achieved in the final dataset including the full ACS observations, the LF presented here extends down to Muv ~ -18.5 at z ~ 7 with one identified object at Muv ~ -15 thanks to the highly-magnified survey areas. This early study forecasts the power of using massive galaxy clusters as cosmic telescopes and its complementarity to blank fields.

We present results from near-infrared spectroscopy of 26 emission-line galaxies at z ~ 2 obtained with the FIRE spectrometer on the Magellan Baade telescope. The sample was selected from the WISP survey, which uses the near-infrared grism of the Hubble Space Telescope Wide Field Camera 3 to detect emission-line galaxies over 0.3 < z < 2.3. Our FIRE follow-up spectroscopy (R~5000) over 1.0-2.5 micron permits detailed measurements of physical properties of the z~2 emission-line galaxies. Dust-corrected star formation rates for the sample range from ~5-100 M_sun yr-1. We derive a median metallicity for the sample of ~0.45 Z_sun, and the estimated stellar masses range from ~10^8.5 - 10^9.5 M_sun. The average ionization parameters measured for the sample are typically much higher than what is found for local star-forming galaxies. We derive composite spectra from the FIRE sample, from which we infer typical nebular electron densities of ~100-400 cm^-3. Based on the location of the galaxies and composite spectra on BPT diagrams, we do not find evidence for significant AGN activity in the sample. Most of the galaxies as well as the composites are offset in the BPT diagram toward higher [O III]/H-beta at a given [N II]/H-alpha, in agreement with other observations of z > 1 star-forming galaxies, but composite spectra derived from the sample do not show an appreciable offset from the local star-forming sequence on the [O III]/H-beta versus [S II]/H-alpha diagram. We infer a high nitrogen-to-oxygen abundance ratio from the composite spectrum, which may contribute to the offset of the high-redshift galaxies from the local star-forming sequence in the [O III]/H-beta versus [N II]/H-alpha diagram. We speculate that the elevated nitrogen abundance could result from substantial numbers of Wolf-Rayet stars in starbursting galaxies at z~2. (Abridged)

We present in this letter the first analysis of a z~8 galaxy candidate found in the Hubble and Spitzer imaging data of Abell 2744, as part of the Hubble Frontier Fields legacy program. We applied the most commonly-used methods to select exceptionally high-z galaxies by combining non-detection and color-criteria using seven HST bands. We used GALFIT on IRAC images for fitting and subtracting contamination of bright nearby sources. The physical properties have been inferred from SED-fitting using templates with and without nebular emission. This letter is focussed on the brightest candidate we found (m$_{F160W}$=26.2) over the 4.9 arcmin$^2$ field of view covered by the WFC3. It shows a non-detection in the ACS bands and at 3.6{\mu}m whereas it is clearly detected at 4.5{\mu}m with rather similar depths. This break in the IRAC data could be explained by strong [OIII]+H{\beta} lines at z~8 which contribute to the 4.5{\mu}m photometry. The best photo-z is found at z~8.0$^{+0.2}_{-0.5}$, although solutions at low-redshift (z~1.9) cannot be completely excluded, but they are strongly disfavoured by the SED-fitting work. The amplification factor is relatively small at {\mu}=1.49$\pm$0.02. The Star Formation Rate in this object is ranging from 8 to 60 Mo/yr, the stellar mass is in the order of M$_{\star}$=(2.5-10) x 10$^{9}$Mo and the size is r~0.35$\pm$0.15 kpc. This object is one of the first z~8 LBG candidates showing a clear break between 3.6{\mu}m and 4.5{\mu}m which is consistent with the IRAC properties of the first spectroscopically confirmed galaxy at a similar redshift. Due to its brightness, the redshift of this object could potentially be confirmed by near infrared spectroscopy with current 8-10m telescopes. The nature of this candidate will be revealed in the coming months with the arrival of new ACS and Spitzer data, increasing the depth at optical and near-IR wavelengths.

We present near-infrared emission line counts and luminosity functions from the HST WFC3 Infrared Spectroscopic Parallels (WISP) program for 29 fields (0.037 deg^2) observed using both the G102 and G141 grisms. Altogether we identify 1048 emission line galaxies with observed equivalent widths greater than 40 Angstroms, 467 of which have multiple detected emission lines. The WISP survey is sensitive to fainter flux levels (3-5x10^-17 ergs/s/cm^2) than the future space near-infrared grism missions aimed at baryonic acoustic oscillation cosmology (1-4x10^-16 ergs/s/cm^2), allowing us to probe the fainter emission line galaxies that the shallower future surveys may miss. Cumulative number counts of 0.7<z<1.5 galaxies reach 10,000 deg^-2 above an H-alpha flux of 2x10^-16 ergs/s/cm^2. H-alpha-emitting galaxies with comparable [OIII] flux are roughly 5 times less common than galaxies with just H-alpha emission at those flux levels. Galaxies with low H-alpha/[OIII] ratios are very rare at the brighter fluxes that future near-infrared grism surveys will probe; our survey finds no galaxies with H-alpha/[OIII] < 0.95 that have H-alpha flux greater than 3x10^-16 ergs/s/cm^2. Our H-alpha luminosity function contains a comparable number density of faint line emitters to that found by the NICMOS near-infrared grism surveys, but significantly fewer (factors of 3-4 less) high luminosity emitters. We also find that our high redshift (z=0.9-1.5) counts are in agreement with the high redshift (z=1.47) narrow band H-alpha survey of HiZELS (Sobral et al. 2013), while our lower redshift luminosity function (z=0.3-0.9) falls slightly below their z=0.84 result. The evolution in both the H-alpha luminosity function from z=0.3--1.5 and the [OIII] luminosity function from z=0.7-2.3 is almost entirely in the L* parameter, which steadily increases with redshift over those ranges.

We present an overview of a 90-orbit Hubble Space Telescope treasury program to obtain near ultraviolet imaging of the Hubble Ultra Deep Field using the Wide Field Camera 3 UVIS detector with the F225W, F275W, and F336W filters. This survey is designed to: (i) Investigate the episode of peak star formation activity in galaxies at 1<z<2.5; (ii) Probe the evolution of massive galaxies by resolving sub-galactic units (clumps); (iii) Examine the escape fraction of ionizing radiation from galaxies at z~2-3; (iv) Greatly improve the reliability of photometric redshift estimates; and (v) Measure the star formation rate efficiency of neutral atomic-dominated hydrogen gas at z~1-3. In this overview paper, we describe the survey details and data reduction challenges, including both the necessity of specialized calibrations and the effects of charge transfer inefficiency. We provide a stark demonstration of the effects of charge transfer inefficiency on resultant data products, which when uncorrected, result in uncertain photometry, elongation of morphology in the readout direction, and loss of faint sources far from the readout. We agree with the STScI recommendation that future UVIS observations that require very sensitive measurements use the instrument's capability to add background light through a "post-flash". Preliminary results on number counts of UV-selected galaxies and morphology of galaxies at z~1 are presented. We find that the number density of UV dropouts at redshifts 1.7, 2.1, and 2.7 is largely consistent with the number predicted by published luminosity functions. We also confirm that the image mosaics have sufficient sensitivity and resolution to support the analysis of the evolution of star-forming clumps, reaching 28-29th magnitude depth at 5 sigma in a 0.2 arcsecond radius aperture depending on filter and observing epoch.

We present the first robust measurement of the high redshift mass-metallicity (MZ) relation at 10^{8}< M/M_{\sun} < 10^{10}, obtained by stacking spectra of 83 emission-line galaxies with secure redshifts between 1.3 < z < 2.3. For these redshifts, infrared grism spectroscopy with the Hubble Space Telescope Wide Field Camera 3 is sensitive to the R23 metallicity diagnostic: ([OII]3726,3729 + [OIII] 4959,5007)/H\beta. Using spectra stacked in four mass quartiles, we find a MZ relation that declines significantly with decreasing mass, extending from 12+log(O/H) = 8.8 at M=10^{9.8} M_{\sun} to 12+log(O/H)= 8.2 at M=10^{8.2} M_{\sun}. After correcting for systematic offsets between metallicity indicators, we compare our MZ relation to measurements from the stacked spectra of galaxies with M>10^{9.5} M_{\sun} and z~2.3. Within the statistical uncertainties, our MZ relation agrees with the z~2.3 result, particularly since our somewhat higher metallicities (by around 0.1 dex) are qualitatively consistent with the lower mean redshift z=1.76 of our sample. For the masses probed by our data, the MZ relation shows a steep slope which is suggestive of feedback from energy-driven winds, and a cosmological downsizing evolution where high mass galaxies reach the local MZ relation at earlier times. In addition, we show that our sample falls on an extrapolation of the star-forming main sequence (the SFR-M_{*} relation) at this redshift. This result indicates that grism emission-line selected samples do not have preferentially high SFRs. Finally, we report no evidence for evolution of the mass-metallicity-SFR plane; our stack-averaged measurements show excellent agreement with the local relation.

We combine Hubble Space Telescope (HST) G102 & G141 NIR grism spectroscopy with HST/WFC3-UVIS, HST/WFC3-IR and Spitzer/IRAC[3.6\mu m] photometry to assemble a sample of massive (log(M_star/M_sun) ~ 11) and quenched galaxies at z~1.5. Our sample of 41 galaxies is the largest with G102+G141 NIR spectroscopy for quenched sources at these redshifts. In contrast to the local Universe, z~1.5 quenched galaxies in the high-mass range have a wide range of stellar population properties. We find their SEDs are well fitted with exponentially decreasing SFHs, and short star-formation time-scales (\tau<100Myr). Quenched galaxies also show a wide distribution in ages, between 1-4Gyr. In the (u-r)_0-versus-mass space quenched galaxies have a large spread in rest-frame color at a given mass. Most quenched galaxies populate the z~1.5 red-sequence (RS), but an important fraction of them (32%) have substantially bluer colors. Although with a large spread, we find that the quenched galaxies ON the RS have older median ages (3.1Gyr) than the quenched galaxies OFF the RS (1.5Gyr). We also show that a rejuvenated SED cannot reproduce the observed stacked spectra of (the bluer) quenched galaxies OFF the RS. We derive the upper limit on the fraction of massive galaxies ON the RS at z~1.5 to be <43%. We speculate that the young quenched galaxies OFF the RS are in a transition phase between vigorous star formation at z>2 and the z~1.5 RS. According to their estimated ages, the time required for quenched galaxies OFF the RS to join their counterparts ON the z~1.5 RS is of the order of ~1Gyr.

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

We report upon new results regarding the Lya output of galaxies, derived from the Lyman alpha Reference Sample (LARS), focusing on Hubble Space Telescope imaging. For 14 galaxies we present intensity images in Lya, Halpha, and UV, and maps of Halpha/Hbeta, Lya equivalent width (EW), and Lya/Halpha. We present Lya and UV light profiles and show they are well-fitted by S\'ersic profiles, but Lya profiles show indices systematically lower than those of the UV (n approx 1-2 instead of >~4). This reveals a general lack of the central concentration in Lya that is ubiquitous in the UV. Photometric growth curves increase more slowly for Lya than the FUV, showing that small apertures may underestimate the EW. For most galaxies, however, flux and EW curves flatten by radii ~10 kpc, suggesting that if placed at high-z, only a few of our galaxies would suffer from large flux losses. We compute global properties of the sample in large apertures, and show total luminosities to be independent of all other quantities. Normalized Lya throughput, however, shows significant correlations: escape is found to be higher in galaxies of lower star formation rate, dust content, mass, and several quantities that suggest harder ionizing continuum and lower metallicity. Eight galaxies could be selected as high-z Lya emitters, based upon their luminosity and EW. We discuss the results in the context of high-z Lya and UV samples. A few galaxies have EWs above 50 AA, and one shows f_escLya of 80%; such objects have not previously been reported at low-z.

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.

Spectroscopic observations of Halpha and Hbeta emission lines of 128 star-forming galaxies in the redshift range 0.75<z<1.5 are presented. These data were taken with slitless spectroscopy using the G102 and G141 grisms of the Wide-Field-Camera 3 (WFC3) on board the Hubble Space Telescope as part of the WFC3 Infrared Spectroscopic Parallel (WISP) survey. Interstellar dust extinction is measured from stacked spectra that cover the Balmer decrement (Halpha/Hbeta). We present dust extinction as a function of Halpha luminosity (down to 3 x 10^{41} erg/s), galaxy stellar mass (reaching 4 x 10^{8} Msun), and rest-frame Halpha equivalent width. The faintest galaxies are two times fainter in Halpha luminosity than galaxies previously studied at z~1.5. An evolution is observed where galaxies of the same Halpha luminosity have lower extinction at higher redshifts, whereas no evolution is found within our error bars with stellar mass. The lower Halpha luminosity galaxies in our sample are found to be consistent with no dust extinction. We find an anti-correlation of the [OIII]5007/Halpha flux ratio as a function of luminosity where galaxies with L_{Halpha}<5 x 10^{41} erg/s are brighter in [OIII]5007 than Halpha. This trend is evident even after extinction correction, suggesting that the increased [OIII]5007/Halpha ratio in low luminosity galaxies is likely due to lower metallicity and/or higher ionization parameters.

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

We present the discovery of three late type (>T4) brown dwarfs, including a probable Y dwarf, in the WFC3 Infrared Spectroscopic Parallels (WISP) Survey. We use the G141 grism spectra to determine the spectral types of the dwarfs and derive distance estimates based on a comparison with nearby T dwarfs with known parallaxes. These are the most distant spectroscopically confirmed T/Y dwarfs, with the farthest at an estimated distance of ~400 pc. We compare the number of cold dwarfs found in the WISP survey with simulations of the brown dwarf mass function. The number found is generally consistent with an initial stellar mass function dN/dM \propto M^{-\alpha} with \alpha = 0.0--0.5, although the identification of a Y dwarf is somewhat surprising and may be indicative of either a flatter absolute magnitude/spectral type relation than previously reported or an upturn in the number of very late type brown dwarfs in the observed volume.

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission-line galaxies at z~2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with 1-sigma detections of emission lines to f > 2.5x10^{-18} erg/s/cm^2, means that the galaxies in the sample are typically ~7 times less massive (median M_* = 10^{9.5} M_sun) than previously studied z~2 emission-line galaxies. Despite their lower mass, the galaxies have OIII/Hb ratios which are very similar to previously studied z~2 galaxies and much higher than the typical emission-line ratios of local galaxies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the two-dimensional spectra of the galaxies for this purpose. In the stacked data the OIII emission line is more spatially concentrated than the Hb emission line with 98.1 confidence. We additionally stack the X-ray data (all sources are individually undetected), and find that the average L(OIII)/L(0.5-10 keV) ratio is intermediate between typical z~0 obscured active galaxies and star-forming galaxies. Together the compactness of the stacked OIII spatial profile and the stacked X-ray data suggest that at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei.

The WFC3 Infrared Spectroscopic Parallel Survey (WISP) uses the Hubble Space Telescope (HST) infrared grism capabilities to obtain slitless spectra of thousands of galaxies over a wide redshift range including the peak of star formation history of the Universe. We select a population of very strong emission-line galaxies with rest-frame equivalent widths higher than 200 A. A total of 176 objects are found over the redshift range 0.35 < z < 2.3 in the 180 arcmin^2 area we analyzed so far. After estimating the AGN fraction in the sample, we show that this population consists of young and low-mass starbursts with higher specific star formation rates than normal star-forming galaxies at any redshift. After spectroscopic follow-up of one of these galaxies with Keck/LRIS, we report the detection at z = 0.7 of an extremely metal-poor galaxy with 12+Log(O/H)= 7.47 +- 0.11. The nebular emission-lines can substantially affect the broadband flux density with a median brightening of 0.3 mag, with examples producing brightening of up to 1 mag. The presence of strong emission lines in low-z galaxies can mimic the color-selection criteria used in the z ~ 8 dropout surveys. In order to effectively remove low redshift interlopers, deep optical imaging is needed, at least 1 mag deeper than the bands in which the objects are detected. Finally, we empirically demonstrate that strong nebular lines can lead to an overestimation of the mass and the age of galaxies derived from fitting of their SED. Without removing emission lines, the age and the stellar mass estimates are overestimated by a factor of 2 on average and up to a factor of 10 for the high-EW galaxies. Therefore the contribution of emission lines should be systematically taken into account in SED fitting of star-forming galaxies at all redshifts.

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

We present the WFC3 Infrared Spectroscopic Parallel (WISP) Survey. WISP is obtaining slitless, near-infrared grism spectroscopy of ~ 90 independent, high-latitude fields by observing in the pure parallel mode with Wide Field Camera-3 on the Hubble Space Telescope for a total of ~ 250 orbits. Spectra are obtained with the G102 (lambda=0.8-1.17 microns, R ~ 210) and G141 grisms (lambda=1.11-1.67 microns, R ~ 130), together with direct imaging in the J- and H-bands (F110W and F140W, respectively). In the present paper, we present the first results from 19 WISP fields, covering approximately 63 square arc minutes. For typical exposure times (~ 6400 sec in G102 and ~ 2700 sec in G141), we reach 5-sigma detection limits for emission lines of 5 x 10^(-17) ergs s^(-1) cm^(-2) for compact objects. Typical direct imaging 5sigma-limits are 26.8 and 25.0 magnitudes (AB) in F110W and F140W, respectively. Restricting ourselves to the lines measured with highest confidence, we present a list of 328 emission lines, in 229 objects, in a redshift range 0.3 < z < 3. The single-line emitters are likely to be a mix of Halpha and [OIII]5007,4959 A, with Halpha predominating. The overall surface density of high-confidence emission-line objects in our sample is approximately 4 per arcmin^(2).These first fields show high equivalent width sources, AGN, and post starburst galaxies. The median observed star formation rate of our Halpha selected sample is 4 Msol/year. At intermediate redshifts, we detect emission lines in galaxies as faint as H_140 ~ 25, or M_R < -19, and are sensitive to star formation rates down to less than 1 Msol/year. The slitless grisms on WFC3 provide a unique opportunity to study the spectral properties of galaxies much fainter than L* at the peak of the galaxy assembly epoch.

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

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

(Abridged) Despite the privileged position that Lyman-alpha (Lya) emission line holds in the exploration of the distant universe and modern observational cosmology, the origin of the observed diversity of lya profiles remains to be thoroughly explained. Observations of nearby star forming galaxies bring their batch of apparent contradictions between Lya emission and their physical parameters, and call for a detailed understanding of the physical processes at work. IZw 18, one of the most metal-poor galaxies known is of particular interest in this context. We use a 3D Lya radiation transfer code to model Hubble Space Telescope (HST) observations of IZw 18 and to fit its Lya spectrum. Different geometrical configurations of the source and the neutral gas are explored. The integrated Lya profile of NW region of IZw 18 is reproduced using the observed small amount of dust (E(B-V) ~ 0.05) and a spherical HI shell with N(HI) = 6.5 x 10^(21) cm^(-2). Such a high column density makes it possible to transform a strong Lya emission (EW(Lya) = 60 A) into a damped absorption even with a small extinction. When a slab geometry is applied and a given line of sight is chosen, the Lya profile can be successfully reproduced with no dust at all and N(HI) = 3 x 10^(21) cm^(-2). The spatial variations of the profile shape are naturally explained by radiation transfer effects. In the case of outflowing Inter Stellar Medium (ISM), as commonly observed in Lyman Break Galaxies (LBGs), a high N(H) and dust content are required to observe Lya in absorption. For nearly static neutral gas as observed in IZw 18 and other low luminosity galaxies only a small amount of dust is required provided a sufficiently high N(H) covers the galaxy.

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]

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

(Abridged) Lyman-alpha (Lya) is a dominant probe of the galaxy population at high-z. However, interpretation of data drawn from Lya alone hinges on the Lya escape fraction which, due to the complex radiative transport, may vary greatly. Here we map the Lya emission from local starburst Haro 11, a Lya emitter and the only known candidate for low-z Lyman continuum emission (LyC). To aid in the interpretation we perform a detailed multi-wavelength analysis and model the stellar population, dust distribution, ionising photon budget, and star-cluster population. We use archival X-ray observations to further constrain properties of the starburst and estimate the HI column density. The Lya morphology is found to be strongly decoupled from stellar and nebular (H-alpha) morphologies. General surface photometry finds only very slight correlation between Lya and H-halpha, E(B-V), and stellar age. Only around the central Lya-bright cluster do we find the Lya/Ha ratio at values predicted by recombination theory. The total Lya escape fraction is found to be just 3%. We compute that ~90% of the Lya photons that escape do so after undergoing multiple resonance scattering events, masking their point of origin. This leads to a largely symmetric distribution and, by increasing the distance that photons must travel to escape, decreases the escape probability significantly. While dust must ultimately be responsible for the destruction of Lya, it plays little role in governing the observed morphology, which is regulated more by ISM kinematics and geometry. We find tentative evidence for local Lya equivalent width in the immediate vicinity of star-clusters being a function of cluster age, consistent with hydrodynamic studies. We estimate the ionising photon production and further constrain the escape fraction at 900 AA to <~9% .