We study the projected spatial offset between the ultraviolet continuum and Ly$\alpha$ emission for a sample of 65 lensed and unlensed galaxies in the Epoch of Reionization (5<z<7), the first such study at these redshifts, in order to understand the potential for these offsets to confuse estimates of the Ly$\alpha$ properties of galaxies observed in slit spectroscopy. While we find that ~40% of galaxies in our sample show significant projected spatial offsets ($|\Delta_{Ly\alpha-UV}|$), we find a modest average offset of 0.61$\pm$0.08 kpc. A small fraction of our sample, ~10%, exhibits offsets of 2-4 kpc, sizes that are larger than the effective radii of typical galaxies at these redshifts. An internal comparison and a comparison to studies at lower redshift yielded no significant evidence of evolution of $|\Delta_{Ly\alpha-UV}|$ with redshift. In our own sample, UV-bright galaxies showed offsets a factor of three greater than their fainter counterparts, 0.89$\pm$0.18 vs. 0.27$\pm$0.05 kpc, respectively. We argue that offsets are likely not the result of merging processes, but are rather due to anisotropic internal processes internal to the galaxies as a result of stellar feedback, which facilitates Ly$\alpha$ fluorescence and/or backscattering from nearby or outflowing gas. The reduction in the Ly$\alpha$ flux due to offset effects for various observational setups was quantified through mock observations of simple simulations. It was found that the loss of Ly$\alpha$ photons for galaxies with average offsets is not, if corrected for, a limiting factor for all but the narrowest slit widths (<0.4''). However, for the largest offsets, if such offsets are mostly perpendicular to the slit major axis, slit losses were found to be extremely severe in cases where slit widths of $\le$1'' were employed, such as those planned for James Webb Space Telescope/NIRSpec observations. (abridged)

The Lyman-$\alpha$ (Ly$\alpha$) emission line has been ubiquitously used to confirm and study high redshift galaxies. We report on the line morphology as seen in the 2D spectra from the VIMOS Ultra Deep Survey in a sample of 914 Ly$\alpha$ emitters from a parent sample of 4192 star-forming galaxies at $2<z_\mathrm{spec}\lesssim6$. The study of the spatial extent of Ly$\alpha$ emission provides insight into the escape of Ly$\alpha$ photons from galaxies. We classify the line emission as either non-existent, coincident, projected spatial offset, or extended with respect to the observed 2D UV continuum emission. The line emitters in our sample are classified as ~45% coincident, ~24% extended and ~11% offset emitters. For galaxies with detected UV continuum, we show that extended Ly$\alpha$ emitters (LAEs) correspond to the highest equivalent width galaxies (with an average $W_\mathrm{Ly\alpha}\sim-22${\AA}). This means that this class of objects is the most common in narrow-band selected samples, which usually select high equivalent width LAEs, $<-20${\AA}. Extended Ly$\alpha$ emitters are found to be less massive, less star-forming, with lower dust content, and smaller UV continuum sizes ($r_{50}\sim0.9$kpc) of all the classes considered here. We also find that galaxies with larger UV-sizes have lower fractions of Ly$\alpha$ emitters. By stacking the spectra per emitter class we find that the weaker Ly$\alpha$ emitters have stronger low ionization inter-stellar medium (ISM) absorption lines. Interestingly, we find that galaxies with Ly$\alpha$ offset emission (median separation of $1.1_{-0.8}^{+1.3}$kpc from UV continuum) show similar velocity offsets in the ISM as those with no visible emission (and different from other Ly$\alpha$ emitting classes). This class of objects may hint at episodes of gas accretion, bright offset clumps, or on-going merging activity into the larger galaxies.

We present the results of a new study investigating the relationship between observed Ly$\alpha$ equivalent width ($W_{\lambda}$(Ly$\alpha$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_{\lambda}$(Ly$\alpha$) across the range $-58 \unicode{xC5} \lesssim$ $W_{\lambda}$(Ly$\alpha$) $\lesssim 110 \unicode{xC5}$ we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_{\lambda}$(Ly$\alpha$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=-18\unicode{xC5}$) and the highest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_{\lambda}$(Ly$\alpha$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_{\lambda}$(Ly$\alpha$) $\leq20\unicode{xC5}$) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_{\lambda}$(Ly$\alpha$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies.

Strong He II emission is produced by low-metallicity stellar populations. Here, we aim to identify and study a sample of He II $\lambda 1640$-emitting galaxies at redshifts of $z \sim 2.5-5$ in the deep VANDELS spectroscopic survey.. We identified a total of 33 Bright He II emitters (S/N > 2.5) and 17 Faint emitters (S/N < 2.5) in the VANDELS survey and used the available deep multi-wavelength data to study their physical properties. After identifying seven potential AGNs in our sample and discarding them from further analysis, we divided the sample of \emph{Bright} emitters into 20 \emph{Narrow} (FWHM < 1000 km s$^{-1}$) and 6 \emph{Broad} (FWHM > 1000 km s$^{-1}$) He II emitters. We created stacks of Faint, Narrow, and Broad emitters and measured other rest-frame UV lines such as O III] and C III] in both individual galaxies and stacks. We then compared the UV line ratios with the output of stellar population-synthesis models to study the ionising properties of He II emitters. We do not see a significant difference between the stellar masses, star-formation rates, and rest-frame UV magnitudes of galaxies with He II and no He II emission. The stellar population models reproduce the observed UV line ratios from metals in a consistent manner, however they under-predict the total number of \heii ionising photons, confirming earlier studies and suggesting that additional mechanisms capable of producing He II are needed, such as X-ray binaries or stripped stars. The models favour subsolar metallicities ($\sim0.1Z_\odot$) and young stellar ages ($10^6 - 10^7$ years) for the He II emitters. However, the metallicity measured for He II emitters is comparable to that of non-He II emitters at similar redshifts. We argue that galaxies with He II emission may have undergone a recent star-formation event, or may be powered by additional sources of He II ionisation.

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works.

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models.

We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range $8.5 < \mathrm{log}(\langle M_{\ast} \rangle/\rm{M}_{\odot}) < 10.2$ we find a strong correlation between stellar metallicity and stellar mass, with stellar metallicity monotonically increasing from $Z_{\ast}/\mathrm{Z}_{\odot} < 0.09$ at $\langle M_{\ast} \rangle = 3.2 \times 10^{8} \rm{M}_{\odot}$ to $Z_{\ast}/Z_{\odot} = 0.27$ at $\langle M_{\ast} \rangle = 1.7 \times 10^{10} \rm{M}_{\odot}$. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the $z=0$ stellar mass-metallicity relation, we find that the $\langle z \rangle = 3.5$ relation is consistent with being shifted to lower metallicities by $\simeq 0.6$ dex. Contrasting our derived stellar metallicities with estimates of gas-phase metallicities at similar redshifts, we find evidence for enhanced $\rm{O}/\rm{Fe}$ ratios of the order (O/Fe) $\gtrsim 1.8$ $\times$ (O/Fe)$_{\odot}$. Finally, by comparing our results to simulation predictions, we find that the $\langle z \rangle = 3.5$ stellar mass-metallicity relation is consistent with current predictions for how outflow strength scales with galaxy mass. This conclusion is supported by an analysis of analytic models, and suggests that the mass loading parameter ($\eta=\dot{M}_{\mathrm{outflow}}/M_{\ast}$) scales as $\eta \propto M_{\ast}^{\beta}$ with $\beta \simeq -0.4$.

We investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.

Detections and non-detections of Lyman alpha (Ly$\alpha$) emission from $z>6$ galaxies ($<1$ Gyr after the Big Bang) can be used to measure the timeline of cosmic reionization. Of key interest to measuring reionization's mid-stages, but also increasing observational challenge, are observations at z > 7, where Ly$\alpha$ redshifts to near infra-red wavelengths. Here we present a search for z > 7.2 Ly$\alpha$ emission in 53 intrinsically faint Lyman Break Galaxy candidates, gravitationally lensed by massive galaxy clusters, in the KMOS Lens-Amplified Spectroscopic Survey (KLASS). With integration times of ~7-10 hours, we detect no Ly$\alpha$ emission with S/N>5 in our sample. We determine our observations to be 80% complete for 5$\sigma$ spatially and spectrally unresolved emission lines with integrated line flux $>5.7\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. We define a photometrically selected sub-sample of 29 targets at $z=7.9\pm0.6$, with a median 5$\sigma$ Ly$\alpha$ EW limit of 58A. We perform a Bayesian inference of the average intergalactic medium (IGM) neutral hydrogen fraction using their spectra. Our inference accounts for the wavelength sensitivity and incomplete redshift coverage of our observations, and the photometric redshift probability distribution of each target. These observations, combined with samples from the literature, enable us to place a lower limit on the average IGM neutral hydrogen fraction of $> 0.76 \; (68\%), \; > 0.46 \; (95\%)$ at z ~ 8, providing further evidence of rapid reionization at z~6-8. We show that this is consistent with reionization history models extending the galaxy luminosity function to $M_\textrm{UV} \lesssim -12$, with low ionizing photon escape fractions, $f_\textrm{esc} \lesssim 15\%$.

We present the results of CANDELSz7, an ESO large program aimed at confirming spectroscopically a homogeneous sample of z~6 and z~7 star forming galaxies. The candidates were selected in the GOODS-South, UDS and COSMOS fields using the official CANDELS catalogs based on H160-band detections. Standard color criteria, which were tailored depending on the ancillary multi-wavelength data available for each field, were applied to select more than 160 candidate galaxies at z~6 and z~7. Deep medium resolution FORS2 spectroscopic observations were then conducted with integration times ranging from 12 to 20 hours, to reach a Lyalpha flux limit of approximately 1-3x 10-18 erg/s/cm^2 at 3sigma. For about 40% of the galaxies we could determine a spectroscopic redshift, mainly through the detection of a single emission line that we interpret as Lyalpha emission, or for some of the brightest objects (H160< 25.5) from the presence of faint continuum and sharp drop that we interpret as a Lyman break. In this paper we present the redshifts and main properties of 65 newly confirmed high redshift galaxies. Adding previous proprietary and archival data we assemble a sample of ~260 galaxies that we use to explore the evolution of the Lyalpha fraction in Lyman break galaxies and the change in the shape of the emission line between z~6 and z~7. We also discuss the accuracy of the CANDELS photometric redshifts in this redshift range.

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection.

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged)

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]

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

Context: Ultraviolet (UV) emission-line spectra are used to spectroscopically confirm high-z galaxies and increasingly also to determine their physical properties. Aims: To interpret the observed UV spectra of distant galaxies in terms of the dominant radiation field and the physical condition of the interstellar medium. Methods: We construct a large grid of photoionization models and derive new spectral UV line diagnostics using equivalent widths (EWs) of CIII]1909, CIV1549 and the line ratios of CIII], CIV, and HeII1640 recombination lines. We apply these diagnostics to a sample of 450 CIII]-emitting galaxies at z=2-4 previously identified in the VIMOS Ultra Deep Survey. Results: We show that the average star-forming galaxy (EW(CIII])~2A) is well described by stellar photoionization from single and binary stars. The inferred metallicity and ionization parameter is typically Z=0.3-0.5Zsun and logU=-2.7 to -3, in agreement with earlier works at similar redshifts. The models also indicate an average age of 50-200Myr since the beginning of the current star-formation, and an ionizing photon production rate, xi_ion, of log(xi_ion/[Hz/erg])~25.3-25.4. Among the sources with EW(CIII])=10-20A, ~30% are likely dominated by AGNs. Their derived metallicity is low, Z=0.02-0.2Zsun, and the ionization parameter higher (logU=-1.7). To explain the average UV observations of the strongest but rarest CIII] emitters (EW(CIII])>20A), we find that stellar photoionization is clearly insufficient. A radiation field consisting of a mix of a young stellar population (log(xi_ion/[Hz/erg])~25.7) plus an AGN component is required. Furthermore an enhanced C/O abundance ratio is needed for metallicities Z=0.1-0.2Zsun and logU=-1.7 to -1.5. Conclusions: The UV diagnostics we propose should serve as an important basis for the interpretation of observations of high-redshift galaxies. [abridged]

We investigate the morphology of the [CII] emission in a sample of "normal" star-forming galaxies at $5<z<7.2$ in relation to their UV (rest-frame) counterpart. We use new ALMA observations of galaxies at $z\sim6-7$ as well as a careful re-analysis of archival ALMA data. In total 29 galaxies were analysed, 21 of which are detected in [CII]. For several of the latter the [CII] emission breaks into multiple components. Only a fraction of these [CII] components, if any, is associated with the primary UV systems, while the bulk of the [CII] emission is associated either with fainter UV components, or not associated with any UV counterpart at the current limits. By taking into account the presence of all these components, we find that the L$_{\rm [CII]}$-SFR relation at early epochs is fully consistent with the local relation, but it has a dispersion of 0.48$\pm$0.07 dex, which is about two times larger than observed locally. We also find that the deviation from the local L$_{\rm [CII]}$-SFR relation has a weak anti-correlation with the EW(Ly$\alpha$). The morphological analysis also reveals that [CII] emission is generally much more extended than the UV emission. As a consequence, these primordial galaxies are characterised by a [CII] surface brightness generally much lower than expected from the local ${\rm \Sigma_{[CII}-\Sigma_{SFR}}$ relation. These properties are likely a consequence of a combination of different effects, namely: gas metallicity, [CII] emission from obscured star-forming regions, strong variations of the ionisation parameter, and circumgalactic gas in accretion or ejected by these primeval galaxies.

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 analyze the CIII]-1908\AA emission properties in a sample of 3899 star-forming galaxies (SFGs) at 2<z<3.8 from the VIMOS Ultra-Deep Survey (VUDS). We find a median EW(CIII])=2.0$\pm$0.2 to 2.2$\pm$0.2\AA ~for the SFG population at 2<z<3 and 3<z<4, resp. About 24% of SFGs are showing EW(CIII])>3\AA, including 20% with 3<EW(CIII])<10\AA ~and 4% with strong emission EW(CIII])>10\AA. A significant fraction of 1.2% of SFGs presents strong CIII] emission 20<EW(CIII])<40\AA. This makes CIII] the second most-frequent emission line in the UV-rest spectra of SFGs after Lyman-$\alpha$. We find a large dispersion in the weak EW(CIII]) - EW(Ly\alpha) correlation, with galaxies showing strong CIII] and no Ly$\alpha$, and vice-versa. SFGs with 10<EW(CIII])<0\AA ~present strong emission lines including CIV-1549, HeII-1640, OIII-1664, but also weaker emission of highly ionized elements SiIV-1403, NIV-1485, NIII-1750, or SiIII-1888, indicating a hard radiation field. We present a broad range of observational evidence supporting the presence of AGN in the strong CIII] emitting population. As EW(CIII]) is rising, we identify powerful outflows with velocities up to 1014 km/s. The strongest CIII] emitters are preferentially located below the SFGs main sequence, with the SFR reduced by x2. In addition, the median stellar age of the strongest emitters is ~0.8 Gyr, three times that of galaxies with EW(CIII])<10\AA. Spectral line analysis presented in a joint paper by Nakajima et al. (2017) confirms that the strongest emitters require an AGN. We conclude that these properties are indicative of AGN feedback in SFGs at 2<z<3.8, contributing to star-formation quenching. We find that quenching timescales of 0.25-0.5x10^9 years are necessary for such AGN feedback to turn part of the star-forming galaxy population with Mstar>10^10 MSun at z~3 into the population of quiescent galaxies observed at redshift z~1-2.

Fast and energetic winds are invoked by galaxy formation models as essential processes in the evolution of galaxies. These outflows can be powered either by star-formation and/or AGN activity, but the relative dominance of the two mechanisms is still under debate. We use spectroscopic stacking analysis to study the properties of the low-ionization phase of the outflow in a sample of 1330 star-forming galaxies (SFGs) and 79 X-ray detected (42<log(L_X)<45 erg/s) Type 2 AGN at 1.7<z<4.6 selected from a compilation of deep optical spectroscopic surveys, mostly zCOSMOS-Deep and VUDS. We measure mean velocity offsets of -150 km/s in the SFGs while in the AGN sample the velocity is much higher (-950 km/s), suggesting that the AGN is boosting the outflow up to velocities that could not be reached only with the star- formation contribution. The sample of X-ray AGN has on average a lower SFR than non-AGN SFGs of similar mass: this, combined with the enhanced outflow velocity in AGN hosts, is consistent with AGN feedback in action. We further divide our sample of AGN into two X-ray luminosity bins: we measure the same velocity offsets in both stacked spectra, at odds with results reported for the highly ionized phase in local AGN, suggesting that the two phases of the outflow may be mixed only up to relatively low velocities, while the highest velocities can be reached only by the highly ionized phase.

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

The reionization of the Universe is one of the most important topics of present day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the HI ionizing background at z~3. We aim at measuring the Lyman continuum escape fraction, which is one of the key parameters to compute the contribution of star-forming galaxies to the UV background. We have used ultra-deep U-band imaging (U=30.2mag at 1sigma) by LBC/LBT in the CANDELS/GOODS-North field, as well as deep imaging in COSMOS and EGS fields, in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27<z<3.40 to faint magnitude limits (L=0.2L*, or equivalently M1500~-19). We have measured through stacks a stringent upper limit (<1.7% at 1sigma) for the relative escape fraction of HI ionizing photons from bright galaxies (L>L*), while for the faint population (L=0.2L*) the limit to the escape fraction is ~10%. We have computed the contribution of star-forming galaxies to the observed UV background at z~3 and we have found that it is not enough to keep the Universe ionized at these redshifts, unless their escape fraction increases significantly (>10%) at low luminosities (M1500>-19). We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental to measure the Lyman continuum escape fraction down to faint magnitudes (M1500~-16) which are inaccessible with the present instrumentation on blank fields.

Determining the average fraction of Lyman continuum (LyC) photons escaping high redshift galaxies is essential for understanding how reionization proceeded in the z>6 Universe. We want to measure the LyC signal from a sample of sources in the Chandra Deep Field South (CDFS) and COSMOS fields for which ultra-deep VIMOS spectroscopy as well as multi-wavelength Hubble Space Telescope (HST) imaging are available. We select a sample of 46 galaxies at $z\sim 4$ from the VIMOS Ultra Deep Survey (VUDS) database, such that the VUDS spectra contain the LyC part of the spectra, that is, the rest-frame range $880-910\AA$. Taking advantage of the HST imaging, we apply a careful cleaning procedure and reject all the sources showing nearby clumps with different colours, that could potentially be lower-redshift interlopers. After this procedure, the sample is reduced to 33 galaxies. We measure the ratio between ionizing flux (LyC at $895\AA$) and non-ionizing emission (at $\sim 1500 \AA$) for all individual sources. We also produce a normalized stacked spectrum of all sources. Assuming an intrinsic average $L_{\nu}(1470)/L_{\nu}(895)$ of 3, we estimate the individual and average relative escape fraction. We do not detect ionizing radiation from any individual source, although we identify a possible LyC emitter with very high Ly$\alpha$ equivalent width (EW). From the stacked spectrum and assuming a mean transmissivity for the sample, we measure a relative escape fraction $f_{esc}^{rel}=0.09\pm0.04$. We also look for correlations between the limits in the LyC flux and source properties and find a tentative correlation between LyC flux and the EW of the Ly$\alpha$ emission line. Our results imply that the LyC flux emitted by $V=25-26$ star-forming galaxies at z$\sim$4 is at most very modest, in agreement with previous upper limits from studies based on broad and narrow band imaging.

Deep observations are revealing a growing number of young galaxies in the first billion year of cosmic time. Compared to typical galaxies at later times, they show more extreme emission-line properties, higher star formation rates, lower masses, and smaller sizes. However, their faintness precludes studies of their chemical abundances and ionization conditions, strongly limiting our understanding of the physics driving early galaxy build-up and metal enrichment. Here we study a rare population of UV-selected, sub$-L^{*}$(z=3) galaxies at redshift 2.4$<z<$3.5 that exhibit all the rest-frame properties expected from primeval galaxies. These low-mass, highly-compact systems are rapidly-forming galaxies able to double their stellar mass in only few tens million years. They are characterized by very blue UV spectra with weak absorption features and bright nebular emission lines, which imply hard radiation fields from young hot massive stars. Their highly-ionized gas phase has strongly sub-solar carbon and oxygen abundances, with metallicities more than a factor of two lower than that found in typical galaxies of similar mass and star formation rate at $z\lesssim$2.5. These young galaxies reveal an early and short stage in the assembly of their galactic structures and their chemical evolution, a vigorous phase which is likely to be dominated by the effects of gas-rich mergers, accretion of metal-poor gas and strong outflows.

We present the discovery and spectrophotometric characterization of a large sample of 164 faint ($i_{AB}$ $\sim$ $23$-$25$ mag) star-forming dwarf galaxies (SFDGs) at redshift $0.13$ $\leq z \leq$ $0.88$ selected by the presence of bright optical emission lines in the VIMOS Ultra Deep Survey (VUDS). We investigate their integrated physical properties and ionization conditions, which are used to discuss the low-mass end of the mass-metallicity relation (MZR) and other key scaling relations. We use optical VUDS spectra in the COSMOS, VVDS-02h, and ECDF-S fields, as well as deep multiwavelength photometry, to derive stellar masses, star formation rates (SFR) and gas-phase metallicities. The VUDS SFDGs are compact (median $r_{e}$ $\sim$ $1.2$ kpc), low-mass ($M_{*}$ $\sim$ $10^7-10^9$ $M_{\odot}$) galaxies with a wide range of star formation rates (SFR($H\alpha$) $\sim 10^{-3}-10^{1}$ $M_{\odot}/yr$) and morphologies. Overall, they show a broad range of subsolar metallicities (12+log(O/H)=$7.26$-$8.7$; $0.04$ $\lesssim Z/Z_{\odot} \lesssim$ $1$). The MZR of SFDGs shows a flatter slope compared to previous studies of galaxies in the same mass range and redshift. We find the scatter of the MZR partly explained in the low mass range by varying specific SFRs and gas fractions amongst the galaxies in our sample. Compared with simple chemical evolution models we find that most SFDGs do not follow the predictions of a "closed-box" model, but those from a gas regulating model in which gas flows are considered. While strong stellar feedback may produce large-scale outflows favoring the cessation of vigorous star formation and promoting the removal of metals, younger and more metal-poor dwarfs may have recently accreted large amounts of fresh, very metal-poor gas, that is used to fuel current star formation.

We present new results on [CII]158$\mu$ m emission from four galaxies in the reionization epoch. These galaxies were previously confirmed to be at redshifts between 6.6 and 7.15 from the presence of the Ly$\alpha$ emission line in their spectra. The Ly$\alpha$ emission line is redshifted by 100-200 km/s compared to the systemic redshift given by the [CII] line. These velocity offsets are smaller than what is observed in z~3 Lyman break galaxies with similar UV luminosities and emission line properties. Smaller velocity shifts reduce the visibility of Ly$\alpha$ and hence somewhat alleviate the need for a very neutral IGM at z~7 to explain the drop in the fraction of Ly$\alpha$ emitters observed at this epoch. The galaxies show [CII] emission with L[CII]=0.6-1.6 x10$^8 L_\odot$: these luminosities place them consistently below the SFR-L[CII] relation observed for low redshift star forming and metal poor galaxies and also below z =5.5 Lyman break galaxies with similar star formation rates. We argue that previous undetections of [CII] in z~7 galaxies with similar or smaller star formation rates are due to selection effects: previous targets were mostly strong Ly$\alpha$ emitters and therefore probably metal poor systems, while our galaxies are more representative of the general high redshift star forming population .

In this paper we investigate the impact of different star formation histories (SFHs) on the relation between stellar mass M$_{*}$ and star formation rate (SFR) using a sample of galaxies with reliable spectroscopic redshift zspec>2 drawn from the VIMOS Ultra-Deep Survey (VUDS). We produce an extensive database of dusty model galaxies, calculated starting from the new library of single stellar population (SSPs) models presented in Cassara' et al. 2013 and weighted by a set of 28 different SFHs based on the Schmidt function, and characterized by different ratios of the gas infall time scale $\tau_{infall}$ to the star formation efficiency $\nu$. The treatment of dust extinction and re-emission has been carried out by means of the radiative transfer calculation. The spectral energy distribution (SED) fitting technique is performed by using GOSSIP+, a tool able to combine both photometric and spectroscopic information to extract the best value of the physical quantities of interest, and to consider the Intergalactic Medium (IGM) attenuation as a free parameter. We find that the main contribution to the scatter observed in the $SFR-M_{*}$ plane is the possibility of choosing between different families of SFHs in the SED fitting procedure, while the redshift range plays a minor role. The majority of the galaxies, at all cosmic times, are best-fit by models with SFHs characterized by a high $\tau_{\rm infall}/\nu$ ratio. We discuss the reliability of the presence of a small percentage of dusty and highly star forming galaxies, in the light of their detection in the FIR.

This paper describes the first data release (DR1) of the VIMOS Ultra Deep Survey (VUDS). The DR1 includes all low-resolution spectroscopic data obtained in 276.9 arcmin2 of the CANDELS-COSMOS and CANDELS-ECFDS survey areas, including accurate spectroscopic redshifts z_spec and individual spectra obtained with VIMOS on the ESO-VLT. A total of 698 objects have a measured redshift, with 677 galaxies, two type-I AGN and a small number of 19 contaminating stars. The targets of the spectroscopic survey are selected primarily on the basis of their photometric redshifts to ensure a broad population coverage. About 500 galaxies have z_spec>2, 48 with z_spec>4, and the highest reliable redshifts reach beyond z_spec=6. This dataset approximately doubles the number of galaxies with spectroscopic redshifts at z>3 in these fields. We discuss the general properties of the sample in terms of the spectroscopic redshift distribution, the distribution of Lyman-alpha equivalent widths, and physical properties including stellar masses M_star and star formation rates (SFR) derived from spectral energy distribution fitting with the knowledge of z_spec. We highlight the properties of the most massive star-forming galaxies, noting the large range in spectral properties, with Lyman-alpha in emission or in absorption, and in imaging properties with compact, multi-component or pair morphologies. We present the catalogue database and data products. All data are publicly available and can be retrieved from a dedicated query-based database available at http://cesam.lam.fr/vuds.

We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization.

Aim: We aim to measure the LyC signal from a sample of sources in the Chandra deep field south. We collect star-forming galaxies (SFGs) and active galactic nuclei (AGN) with accurate spectroscopic redshifts, for which Hubble Space Telescope (HST) coverage and multi-wavelength photometry are available. Method: We selected a sample of about 200 sources at z~3. Taking advantage of HST resolution, we applied a careful cleaning procedure and rejected sources showing nearby clumps with different colours, which could be lower-z interlopers. Our clean sample consisted of 86 SFGs (including 19 narrow-band selected Lya emitters) and 8 AGN (including 6 detected in X-rays). We measured the LyC flux from aperture photometry in four narrow-band filters covering wavelengths below a 912 A rest frame (3.11<z<3.53). We estimated the ratio between ionizing (LyC flux) and 1400 A non-ionizing emissions for AGN and galaxies. Results: By running population synthesis models, we assume an average intrinsic L(1400 A)/L(900 A) ratio of 5 as the representative value for our sample. With this value and an average treatment of the lines of sight of the inter-galactic medium, we estimate the LyC escape fraction relative to the intrinsic value (fesc_rel(LyC)). We do not directly detect ionizing radiation from any individual SFG, but we are able to set a 1(2)sigma upper limit of fesc_rel(LyC)<12(24)%. This result is consistent with other non-detections published in the literature. No meaningful limits can be calculated for the sub-sample of Lya emitters. We obtain one significant direct detection for an AGN at z=3.46, with fesc_rel(LyC) = (72+/-18)%. Conclusions: Our upper limit on fescrel(LyC) implies that the SFGs studied here do not present either the physical properties or the geometric conditions suitable for efficient LyC-photon escape.

The aim of this paper is to investigate spectral and photometric properties of 854 faint ($i_{AB}$<~25 mag) star-forming galaxies (SFGs) at 2<z<2.5 using the VIMOS Ultra-Deep Survey (VUDS) spectroscopic data and deep multi-wavelength photometric data in three extensively studied extragalactic fields (ECDFS, VVDS, COSMOS). These SFGs were targeted for spectroscopy based on their photometric redshifts. The VUDS spectra are used to measure the UV spectral slopes ($\beta$) as well as Ly$\alpha$ equivalent widths (EW). On average, the spectroscopically measured $\beta$ (-1.36$\pm$0.02), is comparable to the photometrically measured $\beta$ (-1.32$\pm$0.02), and has smaller measurement uncertainties. The positive correlation of $\beta$ with the Spectral Energy Distribution (SED)-based measurement of dust extinction, E$_{\rm s}$(B-V), emphasizes the importance of $\beta$ as an alternative dust indicator at high redshifts. To make a proper comparison, we divide these SFGs into three subgroups based on their rest-frame Ly$\alpha$ EW: SFGs with no Ly$\alpha$ emission (SFG$_{\rm N}$; EW$\le$0\AA), SFGs with Ly$\alpha$ emission (SFG$_{\rm L}$; EW$>$0\AA), and Ly$\alpha$ emitters (LAEs; EW$\ge$20\AA). The fraction of LAEs at these redshifts is $\sim$10%, which is consistent with previous observations. We compared best-fit SED-estimated stellar parameters of the SFG$_{\rm N}$, SFG$_{\rm L}$ and LAE samples. For the luminosities probed here ($\sim$L$^*$), we find that galaxies with and without Ly$\alpha$ in emission have small but significant differences in their SED-based properties. We find that LAEs have less dust, and lower star-formation rates (SFR) compared to non-LAEs. We also find that LAEs are less massive compared to non-LAEs, though the difference is smaller and less significant compared to the SFR and E$_{\rm s}$(B-V). [abridged]

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.

[abbreviated] We present a census of Ly\alpha\ emission at $z\gtrsim7$ utilizing deep near infrared HST grism spectroscopy from the first six completed clusters of the Grism Lens-Amplified Survey from Space (GLASS). In 24/159 photometrically selected galaxies we detect emission lines consistent with Ly\alpha\ in the GLASS spectra. Based on the distribution of signal-to-noise ratios and on simulations we expect the completeness and the purity of the sample to be 40-100% and 60-90%, respectively. For the objects without detected emission lines we show that the observed (not corrected for lensing magnification) 1$\sigma$ flux limits reaches $5\times10^{-18}$erg/s/cm$^{2}$ per position angle over the full wavelength range of GLASS (0.8-1.7$\mu$m). Based on the conditional probability of Ly\alpha\ emission measured from the ground at $z\sim7$ we would have expected 12-18 Ly\alpha\ emitters. This is consistent with the number of detections, within the uncertainties, confirming the drop in Ly\alpha\ emission with respect to $z\sim6$. These candidates include a promising source at $z=8.1$. The spatial extent of Ly\alpha\ in a deep stack of the most convincing Ly\alpha\ emitters with $\langle z\rangle=7.2$ is consistent with that of the rest-frame UV continuum. Extended Ly$\alpha$ emission, if present, has a surface brightness below our detection limit, consistent with the properties of lower redshift comparison samples. From the stack we estimate upper limits on rest-frame UV emission line ratios and find $f_\textrm{CIV} / f_\textrm{Ly${\alpha}$} \lesssim 0.32$ and $f_\textrm{CIII]} / f_\textrm{Ly$\alpha$} \lesssim 0.23$ in good agreement with other values published in the literature.

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 Lyman continuum (LyC) flux escaping from high-z galaxies into the IGM is a fundamental quantity to understand the physical processes involved in the reionization epoch. We have investigated a sample of star-forming galaxies at z~3.3 in order to search for possible detections of LyC photons escaping from galaxy halos. UV deep imaging in the COSMOS field obtained with the prime focus camera LBC at the LBT telescope was used together with a catalog of spectroscopic redshifts obtained by the VIMOS Ultra Deep Survey (VUDS) to build a sample of 45 galaxies at z~3.3 with L>0.5L*. We obtained deep LBC images of galaxies with spectroscopic redshifts in the interval 3.27<z<3.40 both in the R and deep U bands. A sub-sample of 10 galaxies apparently shows escape fractions>28% but a detailed analysis of their properties reveals that, with the exception of two marginal detections (S/N~2) in the U band, all the other 8 galaxies are most likely contaminated by the UV flux of low-z interlopers located close to the high-z targets. The average escape fraction derived from the stacking of the cleaned sample was constrained to fesc_rel<2%. The implied HI photo-ionization rate is a factor two lower than that needed to keep the IGM ionized at z~3, as observed in the Lyman forest of high-z QSO spectra or by the proximity effect. These results support a scenario where high redshift, relatively bright (L>0.5L*) star-forming galaxies alone are unable to sustain the level of ionization observed in the cosmic IGM at z~3. Star-forming galaxies at higher redshift and at fainter luminosities (L<<L*) can be the major contributors to the reionization of the Universe only if their physical properties are subject to rapid changes from z~3 to z~6-10. Alternatively, ionizing sources could be discovered looking for fainter sources among the AGN population at high-z.

We present ultraviolet (UV) observations of six nearby Type~Ia supernovae (SNe~Ia) obtained with the Hubble Space Telescope, three of which were also observed in the near-IR (NIR) with Wide-Field Camera~3. UV observations with the Swift satellite, as well as ground-based optical and near-infrared data provide complementary information. The combined data-set covers the wavelength range $0.2$--$2~\mu$m. By also including archival data of SN 2014J, we analyse a sample spanning observed colour excesses up to $E(B-V)=1.4~$mag. We study the wavelength dependent extinction of each individual SN and find a diversity of reddening laws when characterised by the total-to-selective extinction $R_V$. In particular, we note that for the two SNe with $E(B-V)\gtrsim1~$mag, for which the colour excess is dominated by dust extinction, we find $R_V=1.4\pm0.1$ and $R_V=2.8\pm0.1$. Adding UV photometry reduces the uncertainty of fitted $R_V$ by $\sim50\,$% allowing us to also measure $R_V$ of individual low-extinction objects which point to a similar diversity, currently not accounted for in the analyses when SNe~Ia are used for studying the expansion history of the universe.

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

We 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 aim to investigate the effect of the escaping ionizing radiation on the color selection of high redshift galaxies and identify candidate Lyman continuum (LyC) emitters. The intergalactic medium prescription of Inoue et al.(2014) and galaxy synthesis models of Bruzual&Charlot (2003) have been used to properly treat the ultraviolet stellar emission, the stochasticity of the intergalactic transmission and mean free path in the ionizing regime. Color tracks are computed by turning on/off the escape fraction of ionizing radiation. At variance with recent studies, a careful treatment of IGM transmission leads to no significant effects on the high-redshift broad-band color selection. The decreasing mean free path of ionizing photons with increasing redshift further diminishes the contribution of the LyC to broad-band colors. We also demonstrate that prominent LyC sources can be selected under suitable conditions by calculating the probability of a null escaping ionizing radiation. The method is applied to a sample of galaxies extracted from the GOODS-S field. A known LyC source at z=3.795 is successfully recovered as a LyC emitter candidate and another convincing candidate at z=3.212 is reported. A detailed analysis of the two sources (including their variability and morphology) suggests a possible mixture of stellar and non-stellar (AGN) contribution in the ultraviolet. Conclusions: Classical broad-band color selection of 2.5<z<4.5 galaxies does not prevent the inclusion of LyC emitters in the selected samples. Large fesc in relatively bright galaxies (L>0.1L*) could be favored by the presence of a faint AGN not easily detected at any wavelength. A hybrid stellar and non-stellar (AGN) ionizing emission could coexist in these systems and explain the tensions found among the UV excess and the stellar population synthesis models reported in literature.

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 use the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) GOODS-S multi-wavelength catalog to identify counterparts for 20 Lya Emitting (LAE) galaxies at z=2.1. We build several types of stacked Spectral Energy Distributions (SEDs) of these objects. We combine photometry to form average and median flux-stacked SEDs, and postage stamp images to form average and median image-stacked SEDs. We also introduce scaled flux stacks that eliminate the influence of variation in overall brightness. We use the SED fitting code SpeedyMC to constrain the physical properties of individual objects and stacks. Our LAEs at z = 2.1 have stellar masses ranging from 2x10^7 Msun - 8x10^9 Msun (median = 3x10^8 Msun), ages ranging from 4 Myr to 500 Myr (median =100 Myr), and E(B-V) between 0.02 and 0.24 (median = 0.12). We do not observe strong correlations between Lya equivalent width (EW) and stellar mass, age, or E(B-V). The Lya radiative transfer (q) factors of our sample are predominantly close to one and do not correlate strongly with EW or E(B-V), implying that Lya radiative transfer prevents Lya photons from resonantly scattering in dusty regions. The SED parameters of the flux stacks match the average and median values of the individual objects, with the flux-scaled median SED performing best with reduced uncertainties. Median image-stacked SEDs provide a poor representation of the median individual object, and none of the stacking methods captures the large dispersion of LAE properties.

We search for evidence of diffuse Ly-alpha emission from extended neutral hydrogen surrounding Ly-alpha emitting galaxies (LAEs) using deep narrow-band images of the Extended Chandra Deep Field South. By stacking the profiles of 187 LAEs at z = 2.06, 241 LAEs at z = 3.10, and 179 LAEs at z = 3.12, and carefully performing low-surface brightness photometry, we obtain mean surface brightness maps that reach 9.9, 8.7, and 6.2 * 10^{-19} ergs cm^{-2} s^{-1} arcsec^{-2} in the emission line. We undertake a thorough investigation of systematic uncertainties in our surface brightness measurements, and find that our limits are 5--10 times larger than would be expected from Poisson background fluctuations; these uncertainties are often underestimated in the literature. At z ~ 3.1, we find evidence for extended halos with small scale lengths of 5--8 kpc in some, but not all of our sub-samples. We demonstrate that sub-samples of LAEs with low equivalent widths and brighter continuum magnitudes are more likely to possess such halos. At z ~ 2.1, we find no evidence of extended Ly-alpha emission down to our detection limits. Through Monte-Carlo simulations, we also show that we would have detected large diffuse LAE halos if they were present in our data sets. We compare these findings to other measurements in the literature, and discuss possible instrumental and astrophysical reasons for the discrepancies.

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 analyse the clump population of the spiral galaxy Sp 1149 at redshift 1.5. Located behind the galaxy cluster MACS J1149.5+2223, Sp 1149 has been significantly magnified allowing us to study the galaxy on physical scales down to ~100 pc. We have used the publicly available multi-band imaging dataset (CLASH) to reconstruct the spectral energy distributions (SEDs) of the clumps in Sp 1149, and derive, by means of stellar evolutionary models, their physical properties. We found that 40% of the clumps observed in Sp 1149 are older than 30 Myr and can be as old as 300 Myr. These are also the more massive (luminous) clumps in the galaxy. Among the complexes in the local reference sample, the star-forming knots in luminous blue compact galaxies could be considered progenitor analogs of these long-lived clumps. The remaining 60% of clumps have colors comparable to local cluster complexes, suggesting a similar young age. We observe that the Sp 1149 clumps follow the M ~ R^2 relation similar to local cluster complexes, suggesting similar formation mechanisms although they may have different initial conditions (e.g. higher gas surface densities). We suggest that the galaxy is experiencing a slow decline in star formation rate and a likely transitional phase toward a more quiescent star-formation mode. The older clumps have survived between 6 and 20 dynamical times and are all located at projected distances smaller than 4 kpc from the centre. Their current location suggests migration toward the centre and the possibility to be the building blocks of the bulge. On the other hand, the dynamical timescale of the younger clumps are significantly shorter, meaning that they are quite close to their birthplace. We show that the clumps of Sp 1149 may account for the expected metal-rich globular cluster population usually associated with the bulge and thick disk components of local spirals.

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

Small galaxies consisting entirely of population III (pop III) stars may form at high redshifts, and could constitute one of the best probes of such stars. Here, we explore the prospects of detecting gravitationally lensed pop III galaxies behind the galaxy cluster J0717.5+3745 (J0717) with both the Hubble Space Telescope (HST) and the upcoming James Webb Space Telescope (JWST). By projecting simulated catalogs of pop III galaxies at z~7-15 through the J0717 magnification maps, we estimate the lensed number counts as a function of flux detection threshold. We find that the ongoing HST survey CLASH, targeting a total of 25 galaxy clusters including J0717, potentially could detect a small number of pop III galaxies if ~1% of the baryons in these systems have been converted into pop III stars. Using JWST exposures of J0717, this limit can be pushed to ~0.1% of the baryons. Ultra-deep JWST observations of unlensed fields are predicted to do somewhat worse, but will be able to probe pop III galaxies with luminosities intermediate between those detectable in HST/CLASH and in JWST observations of J0717. We also explain how current measurements of the galaxy luminosity function at z=7-10 can be used to constrain pop III galaxy models with very high star formation efficiencies (~10% of the baryons converted into pop III stars).

Lyman Alpha Emitting (LAE) galaxies are thought to be progenitors of present-day L* galaxies. Clustering analyses have suggested that LAEs at z ~ 3 might evolve into LAEs at z ~ 2, but it is unclear whether the physical nature of these galaxies is compatible with this hypothesis. Several groups have investigated the properties of LAEs using spectral energy distribution (SED) fitting, but direct comparison of their results is complicated by inconsistencies in the treatment of the data and in the assumptions made in modeling the stellar populations, which are degenerate with the effects of galaxy evolution. By using the same data analysis pipeline and SED fitting software on two stacked samples of LAEs at z = 3.1 and z = 2.1, and by eliminating several systematic uncertainties that might cause a discrepancy, we determine that the physical properties of these two samples of galaxies are dramatically different. LAEs at z = 3.1 are found to be old (age ~ 1 Gyr) and metal-poor (Z < 0.2 Z_Sun), while LAEs at z = 2.1 appear to be young (age ~ 50 Myr) and metal-rich (Z > Z_Sun). The difference in the observed stellar ages makes it very unlikely that z = 3.1 LAEs evolve directly into z = 2.1 LAEs. Larger samples of galaxies, studies of individual objects and spectroscopic measurements of metallicity at these redshifts are needed to confirm this picture, which is difficult to reconcile with the effects of 1 Gyr of cosmological evolution.

We present properties of individual and composite rest-UV spectra of continuum- and narrowband-selected star-forming galaxies (SFGs) at a redshift of 2<z<3.5 discovered by the MUSYC collaboration in the ECDF-S. Among our sample of 81 UV-bright SFGs, 59 have R<25.5, of which 32 have rest-frame equivalent widths W_{Ly{\alpha}}>20 {\AA}, the canonical limit to be classified as a LAE. We divide our dataset into subsamples based on properties we are able to measure for each individual galaxy: Ly{\alpha} equivalent width, rest-frame UV colors, and redshift. Among our subsample of galaxies with R<25.5, those with rest-frame W_{Ly{\alpha}}>20 {\AA} have bluer UV continua, weaker low-ionization interstellar absorption lines, weaker C IV absorption, and stronger Si II* nebular emission than those with W_{Ly{\alpha}}<20 {\AA}. We measure a typical velocity offset of {\Delta}v~600 km s$^{-1}$ between Ly{\alpha} emission and low-ionization absorption among our subsamples. We find that the interstellar component, as opposed to the stellar component, dominates the high-ionization absorption line profiles. We find the low- and high-ionization Si ionization states have similar kinematic properties, yet the low-ionization absorption is correlated with Ly$\alpha$ emission and the high-ionization absorption is not. These trends are consistent with outflowing neutral gas being in the form of neutral clouds embedded in ionized gas as previously suggested by \cite{Steidel2010}. Moreover, our galaxies with bluer UV colors have stronger Ly{\alpha} emission, weaker low-ionization absorption and more prominent nebular emission line profiles. Among our dataset, UV-bright galaxies with W_{Ly{\alpha}}>20 {\AA} exhibit weaker Ly{\alpha} emission at lower redshifts, although we caution that this could be caused by spectroscopic confirmation of low Ly{\alpha} equivalent width galaxies being harder at z~3 than z~2.

We present GalMC (Acquaviva et al 2011), our publicly available Markov Chain Monte Carlo algorithm for SED fitting, show the results obtained for a stacked sample of Lyman Alpha Emitting galaxies at z ~ 3, and discuss the dependence of the inferred SED parameters on the assumptions made in modeling the stellar populations. We also introduce SpeedyMC, a version of GalMC based on interpolation of pre-computed template libraries. While the flexibility and number of SED fitting parameters is reduced with respect to GalMC, the average running time decreases by a factor of 20,000, enabling SED fitting of each galaxy in about one second on a 2.2GHz MacBook Pro laptop, and making SpeedyMC the ideal instrument to analyze data from large photometric galaxy surveys.

We study the physical properties of 216 z ~ 2.1 LAEs discovered in an ultra-deep narrow-band MUSYC image of the ECDF-S. We fit their stacked Spectral Energy Distribution (SED) using Charlot & Bruzual templates. We consider star formation histories parametrized by the e-folding time parameter tau, allowing for exponentially decreasing (tau>0), exponentially increasing (tau<0), and constant star formation rates. These LAEs are characterized by best fit parameters and 68% confidence intervals of log(M_*/M_sun)=8.6[8.4-9.1], E(B-V)=0.22[0.00-0.31], tau=-0.02[(-4)-18] Gyr, and age_ SF=0.018[0.009-3] Gyr. Thus, we obtain robust measurements of low stellar mass and dust content, but we cannot place meaningful constraints on the age or star formation history of the LAEs. We also calculate the instantaneous SFR to be 35[0.003-170] M_sun/yr, with its average over the last 100 Myr before observation giving <SFR>_100=4[2-30] M_sun/yr. When we compare the results for the same star formation history, LAEs at z~2.1 are dustier and show higher instantaneous SFRs than z~3.1 LAEs, while the observed stellar masses of the two samples seem consistent. LAEs appear to occupy the low-mass end of the distribution of star forming galaxies at z~2. We perform SED fitting on several sub-samples selected based on photometric properties and find that LAE sub-samples at z~2.1 exhibit heterogeneous properties. The IRAC-bright, UV-bright and red LAEs have the largest stellar mass and dust reddening. The UV-faint, IRAC-faint, and high equivalent width LAE sub-samples appear less massive (<10^9 M_sun) and less dusty, with E(B-V) consistent with zero.

We have coadded X-ray flux of all known Lyman Alpha Emitters (LAEs) in the 4 Msec Chandra Deep Field South (CDF-S) region, to place sensitive upper limits on the average unobscured star-formation rate (SFR_X) in these galaxies. A very small fraction of Lyman-Alpha galaxies in the field are individually detected in the X-rays, implying a low fraction of AGN activity. After excluding the few X-ray detected LAEs, we stack the undetected LAEs located in the 4 Ms CDF-S data and 250 ks Extended CDF-S (ECDFS) data, and compute a 1-\sigma upper limit on SFR_X < 14, 28, 28, 140, 440, 880 M$_{\sun}$ yr$^{-1}$ for LAEs located at z = 2.1, 3.1, 3.2, 4.5, 5.7 and 6.5, respectively. The upper limit of SFR_X in LAEs can be then be compared to SFR$_{Ly\alpha}$ derived from Lyman-Alpha line and thus can constrain on the Lyman-Alpha escape fraction ($f^{Esc}_{Ly\alpha}$). We derive a lower limit on f(Lyman-Alpha Escape) > 14% (84 % confidence level, 1-\sigma lower limit) for LAEs at redshift z ~ 2.1 and z ~ 3.1-3.2. At z > 4, the current LAE samples are not of sufficient size to constrain SFR_X well. By averaging all the LAEs at z> 2, the X-ray non-detection constrains f(Lyman-Alpha Escape) > 17% (84 % confidence level, 1-\sigma lower limit), and rejects f(Lyman-Alpha Escape) < 5.7% at the 99.87% confidence level from 2.1 < z < 6.5.

We describe the results of a new, wide-field survey for z=3.1 Ly-alpha emission-line galaxies (LAEs) in the Extended Chandra Deep Field South (ECDF-S). By using a nearly top-hat 5010 Angstrom filter and complementary broadband photometry from the MUSYC survey, we identify a complete sample of 141 objects with monochromatic fluxes brighter than 2.4E-17 ergs/cm^2/s and observers-frame equivalent widths greater than ~ 80 Angstroms (i.e., 20 Angstroms in the rest-frame of Ly-alpha). The bright-end of this dataset is dominated by x-ray sources and foreground objects with GALEX detections, but when these interlopers are removed, we are still left with a sample of 130 LAE candidates, 39 of which have spectroscopic confirmations. This sample overlaps the set of objects found in an earlier ECDF-S survey, but due to our filter's redder bandpass, it also includes 68 previously uncataloged sources. We confirm earlier measurements of the z=3.1 LAE emission-line luminosity function, and show that an apparent anti-correlation between equivalent width and continuum brightness is likely due to the effect of correlated errors in our heteroskedastic dataset. Finally, we compare the properties of z=3.1 LAEs to LAEs found at z=2.1. We show that in the ~1 Gyr after z~3, the LAE luminosity function evolved significantly, with L* fading by ~0.4 mag, the number density of sources with L > 1.5E42 ergs/s declining by ~50%, and the equivalent width scale-length contracting from 70^{+7}_{-5} Angstroms to 50^{+9}_{-6} Angstroms. When combined with literature results, our observations demonstrate that over the redshift range z~0 to z~4, LAEs contain less than ~10% of the star-formation rate density of the universe.

We present GalMC, a MCMC algorithm designed to fit the spectral energy distributions (SED) of galaxies to infer physical properties such as age, stellar mass, dust reddening, metallicity, redshift, and star formation rate. We describe the features of the code and the extensive tests conducted to ensure that our procedure leads to unbiased parameter estimation and accurate evaluation of uncertainties. We compare its performance to grid-based algorithms, showing that the efficiency in CPU time is ~ 100 times better for MCMC for a three dimensional parameter space and increasing with the number of dimensions. We use GalMC to fit the stacked SEDs of two samples of Lyman Alpha Emitters (LAEs) at redshift z=3.1. Our fit reveals that the typical LAE detected in the IRAC 3.6 micron band has age = 0.67 [0.37 - 1.81] Gyr and stellar mass = 3.2 [2.5 - 4.2] x 10^9 M_Sun, while the typical LAE not detected at 3.6 micron has age = 0.06 [0.01-0.2] Gyr and stellar mass = 2 [1.1 - 3.4] x 10^8 M_Sun. The SEDs of both stacks are consistent with the absence of dust. The data do not significantly prefer exponential with respect to constant star formation history. The stellar populations of these two samples are consistent with the previous study by Lai et al, with some differences due to the improved modeling of the stellar populations. A constraint on the metallicity of z=3.1 LAEs from broad-band photometry, requiring Z < Z_Sun at 95% confidence, is found here for the first time.

We present a rest-frame ultraviolet morphological analysis of 108 z=2.1 Lyman Alpha Emitters (LAEs) in the Extended Chandra Deep Field South (ECDF-S) and compare it to a similar sample of 171 LAEs at z=3.1. Using Hubble Space Telescope (HST) images from the Galaxy Evolution from Morphology and SEDs survey, Great Observatories Origins Deep Survey, and Hubble Ultradeep Field, we measure size and photometric component distributions, where photometric components are defined as distinct clumps of UV-continuum emission. At both redshifts, the majority of LAEs have observed half-light radii <~ 2 kpc, but the median half-light radius rises from 1.0 kpc at z=3.1 to 1.4 kpc at z=2.1. A similar evolution is seen in the sizes of individual rest-UV components, but there is no evidence for evolution in the number of multi-component systems. In the z=2.1 sample, we see clear correlations between the size of an LAE and other physical properties derived from its SED. LAEs are found to be larger for galaxies with higher stellar mass, star formation rate, and dust obscuration, but there is no evidence for a trend between equivalent width and half-light radius at either redshift. The presence of these correlations suggests that a wide range of objects are being selected by LAE surveys at z~2, including a significant fraction of objects for which a massive and moderately extended population of old stars underlies the young starburst giving rise to the Lyman alpha emission.

We present MIPS 24 micron observations of the Hubble Deep Field South taken with the Spitzer Space Telescope. The resulting image is 254 arcmin^2 in size and has a sensitivity ranging between ~12 to ~30 microJy rms, with a median sensitivity of ~20 microJy rms. A total of 495 sources have been cataloged with a signal-to-noise ratio greater than 5 sigma. The source catalog is presented as well as source counts which have been corrected for completeness and flux boosting. The IR sources are then combined with MUSYC optical/NIR and ATHDFS radio observations to obtain redshifts and radio flux densities of the sample. We use the IR/radio flux density ratio (q_24) to explore the IR-radio correlation for this IR sample and find q_24 = 0.71 +- 0.31 for sources detected in both IR and radio. The results are extended by stacking IR sources not detected in the radio observations and we derive an average q_24 for redshift bins between 0 < z < 2.5. We find the high redshift (z > 1) sources have an average q_{24} ratio which is better fit by local LIRG SEDs rather than local ULIRG SEDs, indicating that high redshift ULIRGs differ in their IR/radio properties. So ULIRGs at high redshift have SEDs different from those found locally. Infrared faint radio sources are examined, and while nine radio sources do not have a MIPS detection and are therefore radio-loud AGN, only one radio source has an extreme IRAC 3.6 micron to radio flux density ratio indicating it is a radio-loud AGN at z > 1.

We measure the evolution of galaxy clustering out to a redshift of z~1.5 using data from two MUSYC fields, the Extended Hubble Deep Field South (EHDF-S) and the Extended Chandra Deep Field South (ECDF-S). We use photometric redshift information to calculate the projected-angular correlation function, omega(sigma), from which we infer the projected correlation function Xi(sigma). We demonstrate that this technique delivers accurate measurements of clustering even when large redshift measurement errors affect the data. To this aim we use two mock MUSYC fields extracted from a LambdaCDM simulation populated with GALFORM semi-analytic galaxies which allow us to assess the degree of accuracy of our estimates of Xi(sigma) and to identify and correct for systematic effects in our measurements. We study the evolution of clustering for volume limited subsamples of galaxies selected using their photometric redshifts and rest-frame r-band absolute magnitudes. We find that the real-space correlation length r_0 of bright galaxies, M_r<-21 (rest-frame) can be accurately recovered out to z~1.5, particularly for ECDF-S given its near-infrared photometric coverage. There is mild evidence for a luminosity dependent clustering in both fields at the low redshift samples (up to <z>=0.57), where the correlation length is higher for brighter galaxies by up to 1Mpc/h between median rest-frame r-band absolute magnitudes of -18 to -21.5. As a result of the photometric redshift measurement, each galaxy is assigned a best-fit template; we restrict to E and E+20%Sbc types to construct subsamples of early type galaxies (ETGs). Our ETG samples show a strong increase in r_0 as the redshift increases, making it unlikely (95% level) that ETGs at median redshift z_med=1.15 are the direct progenitors of ETGs at z_med=0.37 with equivalent passively evolved luminosities. (ABRIDGED)

We discovered a sample of 250 Ly-Alpha emitting (LAE) galaxies at z=2.1 in an ultra-deep 3727 A narrow-band MUSYC image of the Extended Chandra Deep Field-South. LAEs were selected to have rest-frame equivalent widths (EW) > 20 A and emission line fluxes > 2.0 x 10^(-17)erg /cm^2/s, after carefully subtracting the continuum contributions from narrow band photometry. The median flux of our sample is 4.2 x 10^(-17)erg/cm^2/s, corresponding to a median Lya luminosity = 1.3 x 10^(42) erg/s at z=2.1. At this flux our sample is > 90% complete. Approximately 4% of the original NB-selected candidates were detected in X-rays by Chandra, and 7% were detected in the rest-frame far-UV by GALEX. At luminosity>1.3 x 10^42 erg/s, the equivalent width distribution is unbiased and is represented by an exponential with scale-length of 83+/-10 A. Above this same luminosity threshold, we find a number density of 1.5+/-0.5 x 10^-3 Mpc^-3. Neither the number density of LAEs nor the scale-length of their EW distribution show significant evolution from z=3 to z=2. We used the rest frame UV luminosity to estimate a median star formation rate of 4 M_(sun) /yr. The median rest frame UV slope, parametrized by B-R, is that typical of dust-free, 0.5-1 Gyr old or moderately dusty, 300-500 Myr old populations. Approximately 40% of the sample occupies the z~2 star-forming galaxy locus in the UVR two color diagram. Clustering analysis reveals that LAEs at z=2.1 have r_0=4.8+/-0.9 Mpc and a bias factor b=1.8+/-0.3. This implies that z=2.1 LAEs reside in dark matter halos with median masses Log(M/M_(sun))=11.5^(+0.4)_(-0.5), which are among of the lowest-mass halos yet probed at this redshift. We used the Sheth-Tormen conditional mass function to study the descendants of these LAEs and found that their typical present-day descendants are local galaxies with L* properties, like the Milky Way.