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

We report the spectroscopic confirmation of a new protocluster in the COSMOS field at $z$ $\sim$ 2.2, COSMOS Cluster 2.2 (CC2.2), originally identified as an overdensity of narrowband selected H$\alpha$ emitting candidates. With only two masks of Keck/MOSFIRE near-IR spectroscopy in both $H$ ($\sim$ 1.47-1.81 $\mu$m) and $K$ ($\sim$ 1.92-2.40 $\mu$m) bands ($\sim$ 1.5 hour each), we confirm 35 unique protocluster members with at least two emission lines detected with S/N $>$ 3. Combined with 12 extra members from the zCOSMOS-deep spectroscopic survey (47 in total), we estimate a mean redshift and a line-of-sight velocity dispersion of $z_{mean}$=2.23224 $\pm$ 0.00101 and $\sigma_{los}$=645 $\pm$ 69 km s$^{-1}$ for this protocluster, respectively. Assuming virialization and spherical symmetry for the system, we estimate a total mass of $M_{vir}$ $\sim$ $(1-2) \times$10$^{14}$ $M_{\odot}$ for the structure. We evaluate a number density enhancement of $\delta_{g}$ $\sim$ 7 for this system and we argue that the structure is likely not fully virialized at $z$ $\sim$ 2.2. However, in a spherical collapse model, $\delta_{g}$ is expected to grow to a linear matter enhancement of $\sim$ 1.9 by $z$=0, exceeding the collapse threshold of 1.69, and leading to a fully collapsed and virialized Coma-type structure with a total mass of $M_{dyn}$($z$=0) $\sim$ 9.2$\times$10$^{14}$ $M_{\odot}$ by now. This observationally efficient confirmation suggests that large narrowband emission-line galaxy surveys, when combined with ancillary photometric data, can be used to effectively trace the large-scale structure and protoclusters at a time when they are mostly dominated by star-forming galaxies.

A prediction of the classic active galactic nuclei (AGN) unification model is the presence of ionisation cones with different orientations depending on the AGN type. Confirmations of this model exist for present times, but it is less clear in the early Universe. Here, we use the morphology of giant Ly$\alpha$ nebulae around AGNs at redshift z$\sim$3 to probe AGN emission and therefore the validity of the AGN unification model at this redshift. We compare the spatial morphology of 19 nebulae previously found around type I AGNs with a new sample of 4 Ly$\alpha$ nebulae detected around type II AGNs. Using two independent techniques, we find that nebulae around type II AGNs are more asymmetric than around type I, at least at radial distances $r>30$~physical kpc (pkpc) from the ionizing source. We conclude that the type I and type II AGNs in our sample show evidence of different surrounding ionising geometries. This suggests that the classical AGN unification model is also valid for high-redshift sources. Finally, we discuss how the lack of asymmetry in the inner parts (r$\lesssim$30 pkpc) and the associated high values of the HeII to Ly$\alpha$ ratios in these regions could indicate additional sources of (hard) ionizing radiation originating within or in proximity of the AGN host galaxies. This work demonstrates that the morphologies of giant Ly$\alpha$ nebulae can be used to understand and study the geometry of high redshift AGNs on circum-nuclear scales and it lays the foundation for future studies using much larger statistical samples.

The decline in abundance of Lyman-$\alpha$ (Ly$\alpha$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$\alpha$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$\alpha$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$\alpha$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 5\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$\alpha$ peak. We discuss our results in the context of the interpretation of observations.

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

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

Despite recent progress in understanding Lyman-alpha (Lya) emitters (LAEs), relatively little is known regarding their typical black hole activity across cosmic time. Here, we study the X-ray and radio properties of ~4000 LAEs at 2.2<z<6 from the SC4K survey in the COSMOS field. We detect 254 (6.8+-0.4%) of the LAEs individually in the X-rays (S/N>3) and find an average luminosity of 10^{44.31+-0.01} erg/s and an average black hole accretion rate (BHAR) of 0.72+-0.01 Msun/yr, consistent with moderate to high accreting AGN. We detect 120 sources in deep radio data (radio AGN fraction of 3.2%+-0.3%). The global AGN fraction (8.6%+-0.4%) rises with Lya luminosity and declines with increasing redshift. For X-ray detected LAEs, Lya luminosities correlate with the BHARs, suggesting that Lya luminosity becomes an accretion rate indicator. Most LAEs (93.1+-0.6%) at 2<z<6 have no detectable X-ray emission (BHARs<0.017 Msun/yr). The median star formation rate (SFR) of star-forming LAEs from Lya and radio luminosities is 7.6^{+6.6}_{-2.8} Msun/yr. The black hole to galaxy growth ratio (BHAR/SFR) for LAEs is <0.0022, consistent with typical star forming galaxies and the local BHAR/SFR relation. We conclude that LAEs at 2<z<6 include two different populations: an AGN population, where Lya luminosity traces BHAR, making them bright in Lya, and another with low SFRs which remain undetected in even the deepest X-ray stacks but is detected in the radio stacks.

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

We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.4^{+10.5}_{-2.4}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr, combined with very blue UV slopes of beta=-2.1$^{+0.5}_{-0.4}$, but with significant variations within the population. M$_{UV}$ and beta are correlated in a similar way to UV-selected sources, but LAEs are consistently bluer. This suggests that LAEs are the youngest and/or most dust-poor subset of the UV-selected population. We also study the Lya rest-frame equivalent width (EW$_0$) and find 45 "extreme" LAEs with EW$_0>240$ A (3 $\sigma$), implying a low number density of $(7\pm1)\times10^{-7}$ Mpc$^{-3}$. Overall, we measure little to no evolution of the Lya EW$_0$ and scale length parameter ($w_0$) which are consistently high (EW$_0=140^{+280}_{-70}$ A, $w_0=129^{+11}_{-11}$ A) from z~6 to z~2 and below. However, $w_0$ is anti-correlated with M$_{UV}$ and stellar mass. Our results imply that sources selected as LAEs have a high Lya escape fraction (f$_{esc, Lya}$) irrespective of cosmic time, but f$_{esc, Lya}$ is still higher for UV-fainter and lower mass LAEs. The least massive LAEs ($<10^{9.5}$ M$_{\odot}$) are typically located above the star formation "Main Sequence" (MS), but the offset from the MS decreases towards z~6 and towards $10^{10}$ M$_{\odot}$. Our results imply a lack of evolution in the properties of LAEs across time and reveals the increasing overlap in properties of LAEs and UV-continuum selected galaxies as typical star-forming galaxies at high redshift effectively become LAEs.

We present new deep ALMA and HST/WFC3 observations of MASOSA and VR7, two luminous Ly$\alpha$ emitters (LAEs) at $z=6.5$, for which the UV continuum level differ by a factor four. No IR dust continuum emission is detected in either, indicating little amounts of obscured star formation and/or high dust temperatures. MASOSA, with a UV luminosity M$_{1500}=-20.9$, compact size and very high Ly$\alpha$ EW$_{0}\approx145$ A, is undetected in [CII] to a limit of L$_{\rm [CII]}<2.2\times10^7$ L$_{\odot}$ implying a metallicity $Z\lesssim0.07 Z_{\odot}$. Intriguingly, our HST data indicates a red UV slope $\beta=-1.1\pm0.7$, at odds with the low dust content. VR7, which is a bright (M$_{1500}=-22.4$) galaxy with moderate color ($\beta=-1.4\pm0.3$) and Ly$\alpha$ EW$_0 = 34$ A, is clearly detected in [CII] emission (S/N=15). VR7's rest-frame UV morphology can be described by two components separated by $\approx1.5$ kpc and is globally more compact than the [CII] emission. The global [CII]-UV ratio indicates $Z\approx0.2 Z_{\odot}$, but there are large variations in the UV-[CII] ratio on kpc scales. We also identify diffuse, possibly outflowing, [CII]-emitting gas at $\approx 100$ km s$^{-1}$ with respect to the peak. VR7 appears assembling its components at a slightly more evolved stage than other luminous LAEs, with outflows already shaping its direct environment at $z\sim7$. Our results further indicate that the global [CII]-UV relation steepens at SFR $<30$ M$_{\odot}$ yr$^{-1}$, naturally explaining why the [CII]-UV ratio is anti-correlated with Ly$\alpha$ EW in many, but not all, observed LAEs.

The observed properties of the Lyman-$\alpha$ (Ly$\alpha$) emission line are a powerful probe of neutral gas in and around galaxies. We present spatially resolved Ly$\alpha$ spectroscopy with VLT/MUSE targeting VR7, a UV-luminous galaxy at $z=6.532$ with moderate Ly$\alpha$ equivalent width (EW$_0\approx38$ {\AA}). These data are combined with deep resolved [CII]$_{\rm 158 \mu m}$ spectroscopy obtained with ALMA and UV imaging from {\it HST} and we also detect UV continuum with MUSE. Ly$\alpha$ emission is clearly detected with S/N $\approx40$ and FWHM of 374 km s$^{-1}$. Ly$\alpha$ and [CII] are similarly extended beyond the UV, with effective radius r$_{\rm eff} = 2.1\pm0.2$ kpc for a single exponential model or r$_{\rm eff, Ly\alpha, halo} = 3.45^{+1.08}_{-0.87}$ kpc when measured jointly with the UV continuum. The Ly$\alpha$ profile is broader and redshifted with respect to the [CII] line (by 213 km s$^{-1}$), but there are spatial variations that are qualitatively similar in both lines and coincide with resolved UV components. This suggests that the emission originates from two components with plausibly different HI column densities. We place VR7 in the context of other galaxies at similar and lower redshift. The Ly$\alpha$ halo scale length is similar at different redshifts and velocity shifts with respect to the systemic are typically smaller. Overall, we find little indications of a more neutral vicinity at higher redshift. This means that the local ($\sim 10$ kpc) neutral gas conditions that determine the observed Ly$\alpha$ properties in VR7 resemble the conditions in post-re-ionisation galaxies.

Distant luminous Lyman-alpha emitters are excellent targets for detailed observations of galaxies in the epoch of reionisation. Spatially resolved observations of these galaxies allow us to simultaneously probe the emission from young stars, partially ionised gas in the interstellar medium and to constrain the properties of the surrounding hydrogen in the circumgalactic medium. We review recent results from (spectroscopic) follow-up studies of the rest-frame UV, Lyman-alpha and [CII] emission in luminous galaxies observed ~500 Myr after the Big Bang with ALMA, HST/WFC3 and VLT/X-SHOOTER. These galaxies likely reside in early ionised bubbles and are complex systems, consisting of multiple well separated and resolved components where traces of metals are already present.

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

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged)

The well-known quasar SDSS J095253.83+011421.9 (J0952+0114) at z=3.02 has one of the most peculiar spectra discovered so far, showing the presence of narrow Ly$\alpha$ and broad metal emission lines. Although recent studies have suggested that a Proximate Damped Ly$\alpha$ system (PDLA) causes this peculiar spectrum, the origin of the gas associated with the PDLA is unknown. Here we report the results of MUSE observations that reveal a new giant ($\approx$ 100 physical kpc) Lyman $\alpha$ nebula. The detailed analysis of the Ly$\alpha$ velocity, velocity dispersion, and surface brightness profiles suggests that the J0952+0114 Ly$\alpha$ nebula shares similar properties of other QSO nebulae previously detected with MUSE, implying that the PDLA in J0952+0144 is covering only a small fraction of the QSO emission solid angle. We also detected bright and spectrally narrow CIV$\lambda$1550 and HeII$\lambda$1640 extended emission around J0952+0114 with velocity centroids similar to the peak of the extended and central narrow Ly$\alpha$ emission. The presence of a peculiarly bright, unresolved, and relatively broad HeII$\lambda$1640 emission in the central region at exactly the same PDLA redshift hints at the possibility that the PDLA originates in a clumpy outflow with a bulk velocity of about 500 km/s. The smaller velocity dispersion of the large scale Ly$\alpha$ emission suggests that the high-speed outflow is confined to the central region. Lastly, the derived spatially resolved HeII/Ly$\alpha$ and CIV/Ly$\alpha$ maps show a positive gradient with the distance to the QSO hinting at a non-homogeneous ionization parameter distribution.

To understand star formation in galaxies, we investigate the star formation rate (SFR) surface density ($\Sigma_{\rm SFR}$) profiles for galaxies, based on a well-defined sample of 976 star-forming MaNGA galaxies. We find that the typical $\Sigma_{\rm SFR}$ profiles within 1.5Re of normal SF galaxies can be well described by an exponential function for different stellar mass intervals, while the sSFR profile shows positive gradients, especially for more massive SF galaxies. This is due to the more pronounced central cores or bulges rather than the onset of a `quenching' process. While galaxies that lie significantly above (or below) the star formation main sequence (SFMS) show overall an elevation (or suppression) of $\Sigma_{\rm SFR}$ at all radii, this central elevation (or suppression) is more pronounced in more massive galaxies. The degree of central enhancement and suppression is quite symmetric, suggesting that both the elevation and suppression of star formation are following the same physical processes. Furthermore, we find that the dispersion in $\Sigma_{\rm SFR}$ within and across the population is found to be tightly correlated with the inferred gas depletion time, whether based on the stellar surface mass density or the orbital dynamical time. This suggests that we are seeing the response of a simple gas-regulator system to variations in the accretion rate. This is explored using a heuristic model that can quantitatively explain the dependence of $\sigma(\Sigma_{\rm SFR})$ on gas depletion timescale. Variations in accretion rate are progressively more damped out in regions of low star-formation efficiency leading to a reduced amplitude of variations in star-formation.

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation.

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

Lyman-alpha (Lya) is intrinsically the brightest line emitted from active galaxies. While it originates from many physical processes, for star-forming galaxies the intrinsic Lya luminosity is a direct tracer of the Lyman-continuum (LyC) radiation produced by the most massive O- and early-type B-stars with lifetimes of a few Myrs. As such, Lya luminosity should be an excellent instantaneous star formation rate (SFR) indicator. However, its resonant nature and susceptibility to dust as a rest-frame UV photon makes Lya very hard to interpret due to the uncertain Lya escape fraction, f$_{esc,Ly\alpha}$. Here we explore results from the CALYMHA survey at z=2.2, follow-up of Lya emitters (LAEs) at z=2.2-2.6 and a z~0-0.3 compilation of LAEs to directly measure f$_{esc,Ly\alpha}$ with H-alpha (Ha). We derive a simple empirical relation that robustly retrieves f$_{esc,Ly\alpha}$ as a function of Lya rest-frame EW (EW$_0$): f$_{esc,Ly\alpha}$= 0.0048 EW$_0$[A]$\pm0.05$ and we show that it constrains a well-defined anti-correlation between ionisation efficiency and dust extinction in LAEs. Observed Lya luminosities and EW$_0$ are easy measurable quantities at high redshift, thus making our relation a practical tool to estimate intrinsic Lya and LyC luminosities under well controlled and simple assumptions. Our results allow observed Lya luminosities to be used to compute SFRs for LAEs at z~0-2.6 within 0.2 dex of the Ha dust corrected SFRs. We apply our empirical SFR(Lya,EW$_0$) calibration to several sources at z>2.6 to find that star-forming LAEs have SFRs typically ranging from 0.1 to 20 M$_{\odot}$/yr and that our calibration might be even applicable for the most luminous LAEs within the epoch of re-ionisation. Our results imply high ionisation efficiencies and low dust content in LAEs across cosmic time, and will be easily tested with future observations with JWST.

Observations have revealed that the star formation rate (SFR) and stellar mass (M$_{\rm star}$) of star-forming galaxies follow a tight relation known as the galaxy main sequence. However, what physical information is encoded in this relation is under debate. Here, we use the EAGLE cosmological hydrodynamical simulation to study the mass dependence, evolution and origin of scatter in the SFR-M$_{\rm star}$ relation. At $z=0$, we find that the scatter decreases slightly with stellar mass from 0.35 dex at M$_{\rm star} \approx 10^9$ M$_{\odot}$ to 0.30 dex at M$_{\rm star} \gtrsim 10^{10.5}$ M$_{\odot}$. The scatter decreases from $z=0$ to $z=5$ by 0.05 dex at M$_{\rm star} \gtrsim 10^{10}$ M$_{\odot}$ and by 0.15 dex for lower masses. We show that the scatter at $z=0.1$ originates from a combination of fluctuations on short time-scales (ranging from 0.2-2 Gyr) that are presumably associated with self-regulation from cooling, star formation and outflows, but is dominated by long time-scale ($\sim 10$ Gyr) variations related to differences in halo formation times. Shorter time-scale fluctuations are relatively more important for lower-mass galaxies. At high masses, differences in black hole formation efficiency cause additional scatter, but also diminish the scatter caused by different halo formation times. While individual galaxies cross the main sequence multiple times during their evolution, they fluctuate around tracks associated with their halo properties, i.e. galaxies above/below the main sequence at $z = 0.1$ tend to have been above/below the main sequence for $\gg1$ Gyr.

We investigate the clustering and halo properties of $\sim 5000$ Ly$\alpha$-selected emission line galaxies (LAEs) from the Slicing COSMOS 4K (SC4K) and from archival NB497 imaging of SA22 split in 15 discrete redshift slices between $z \sim 2.5 - 6$. We measure clustering lengths of $r_0 \sim 3 - 6\ h^{-1}$ Mpc and typical halo masses of $\sim 10^{11}$ M$_\odot$ for our narrowband-selected LAEs with typical $L_{\rm{Ly}\alpha} \sim 10^{42 - 43}$ erg s$^{-1}$. The intermediate band-selected LAEs are observed to have $r_0 \sim 3.5 - 15\ h^{-1}$ Mpc with typical halo masses of $\sim 10^{11 - 12}$ M$_\odot$ and typical $L_{\rm{Ly}\alpha} \sim 10^{43 - 43.6}$ erg s$^{-1}$. We find a strong, redshift-independent correlation between halo mass and Ly$\alpha$ luminosity normalized by the characteristic Ly$\alpha$ luminosity, $L^\star(z)$. The faintest LAEs ($L \sim 0.1\ L^\star(z)$) typically identified by deep narrowband surveys are found in $10^{10}$ M$_\odot$ halos and the brightest LAEs ($L \sim 7\ L^\star(z)$) are found in $\sim 5 \times 10^{12}$ M$_\odot$ halos. A dependency on the rest-frame 1500 \AA~UV luminosity, M$_\rm{UV}$, is also observed where the halo masses increase from $10^{11}$ to $10^{13}$ M$_\odot$ for M$_\rm{UV} \sim -19$ to $-23.5$ mag. Halo mass is also observed to increase from $10^{9.8}$ to $10^{12.3}$ M$_\odot$ for dust-corrected UV star formation rates from $\sim 0.6$ to $10$ M$_\odot$ yr$^{-1}$ and continues to increase up to $10^{13.5}$ M$_\odot$ in halo mass, where the majority of those sources are AGN. All the trends we observe are found to be redshift-independent. Our results reveal that LAEs are the likely progenitors of a wide range of galaxies depending on their luminosity, from dwarf-like, to Milky Way-type, to bright cluster galaxies. LAEs therefore provide unique insight into the early formation and evolution of the galaxies we observe in the local Universe.

We present new HST/WFC3 observations and re-analyse VLT data to unveil the continuum, variability and rest-frame UV lines of the multiple UV clumps of the most luminous Ly$\alpha$ emitter at z=6.6, CR7. Our re-reduced, flux calibrated X-SHOOTER spectra of CR7 reveal a HeII emission line in observations obtained along the major axis of Lyman-alpha (Lya) emission with the best seeing conditions. HeII is spatially offset by +0.8'' from the peak of Lya emission, and it is found towards clump B. Our WFC3 grism spectra detects the UV continuum of CR7's clump A, yielding a power law with $\beta=-2.5^{+0.6}_{-0.7}$ and $M_{UV}=-21.87^{+0.25}_{-0.20}$. No significant variability is found for any of the UV clumps on their own, but there is tentative (~2.2$\sigma$) brightening of CR7 in F110W as a whole from 2012 to 2017. HST grism data fail to robustly detect rest-frame UV lines in any of the clumps, implying fluxes <2x10$^{-17}$ erg s$^{-1}$ cm$^{-2}$ (3 $\sigma$). We perform CLOUDY modelling to constrain the metallicity and the ionising nature of CR7. CR7 seems to be actively forming stars without any clear AGN activity in clump A, consistent with a metallicity of ~0.05-0.2 Z$_{\odot}$. Component C or an inter-clump component between B and C may host a high ionisation source. Our results highlight the need for spatially resolved information to study the formation and assembly of early galaxies.

Distant luminous Lyman-alpha emitters (LAEs) are excellent targets for spectroscopic observations of galaxies in the epoch of reionisation (EoR). We present deep high-resolution (R=5000) VLT/X-SHOOTER observations, along with an extensive collection of photometric data of `COLA1', a proposed double peaked LAE at z=6.6 (Hu et al. 2016). We rule out that COLA1's emission line is an [OII] doublet at z=1.475 on the basis of i) the asymmetric red line-profile and flux ratio of the peaks (blue/red=$0.31\pm0.03$) and ii) an unphysical [OII]/Ha ratio ([OII]/Ha > 22). We show that COLA1's observed B-band flux is explained by a faint extended foreground LAE, for which we detect Lya and [OIII] at z=2.142. We thus conclude that COLA1 is a real double-peaked LAE at z=6.593, the first discovered at z>6, confirming the result from Hu et al. (2016). COLA1 is UV luminous (M$_{1500}=-21.6\pm0.3$), has a high equivalent width (EW$_{0}$~120 \AA) and very compact Lya emission (r$_{50} = 0.3$ kpc). Relatively weak inferred Hb+[OIII] line-emission from Spitzer/IRAC indicates an extremely low metallicity of Z<1/20 Z$_{\odot}$ or reduced strength of nebular lines due to high escape of ionising photons. The small Lya peak separation of $220\pm20$ km/s implies a low HI column density and an ionising photon escape fraction of ~15-30 %, providing the first direct evidence that such galaxies contribute actively to the reionisation of the Universe at z>6. Based on simple estimates, we find that COLA1 could have provided just enough photons to reionise its own ~0.3 pMpc bubble, allowing the blue Lya line to be observed. However, we also discuss alternative scenarios explaining the detected double peaked nature of COLA1. Our results show that future high-resolution observations of statistical samples of double peaked LAEs at z>5 are a promising probe of the occurrence of ionised regions around galaxies in the EoR.

Observations show that star-forming galaxies reside on a tight three-dimensional plane between mass, gas-phase metallicity and star formation rate (SFR), which can be explained by the interplay between metal-poor gas inflows, SFR and outflows. However, different metals are released on different time-scales, which may affect the slope of this relation. Here, we use central, star-forming galaxies with M$_{\rm star}=10^{9.0-10.5}$ M$_{\odot}$ from the EAGLE hydrodynamical simulation to examine three-dimensional relations between mass, SFR and chemical enrichment using absolute and relative C, N, O and Fe abundances. We show that the scatter is smaller when gas-phase $\alpha$-enhancement is used rather than metallicity. A similar plane also exists for stellar $\alpha$-enhancement, implying that present-day specific SFRs are correlated with long time-scale star formation histories. Between $z=0$ and 1, the $\alpha$-enhancement plane is even more insensitive to redshift than the plane using metallicity. However, it evolves at $z>1$ due to lagging iron yields. At fixed mass, galaxies with higher SFRs have star formation histories shifted toward late times, are more $\alpha$-enhanced and this $\alpha$-enhancement increases with redshift as observed. These findings suggest that relations between physical properties inferred from observations may be affected by systematic variations in $\alpha$-enhancements.

Deep narrow-band surveys have revealed a large population of faint Lyman-alpha (Lya) emitters (LAEs) in the distant Universe, but relatively little is known about the most luminous sources ($L_{Ly\alpha}>10^{42.7}$ erg/s; $L_{Ly\alpha}>L^*_{Ly\alpha}$). Here we present the spectroscopic follow-up of 21 luminous LAEs at z~2-3 found with panoramic narrow-band surveys over five independent extragalactic fields (~4x10$^6$ Mpc$^{3}$ surveyed at z~2.2 and z~3.1). We use WHT/ISIS, Keck/DEIMOS and VLT/X-SHOOTER to study these sources using high ionisation UV lines. Luminous LAEs at z~2-3 have blue UV slopes ($\beta=-2.0^{+0.3}_{-0.1}$), high Lya escape fractions ($50^{+20}_{-15}$%) and span five orders of magnitude in UV luminosity ($M_{UV}\approx-19$ to -24). Many (70%) show at least one high ionisation rest-frame UV line such as CIV, NV, CIII], HeII or OIII], typically blue-shifted by ~100-200 km/s relative to Lya. Their Lya profiles reveal a wide variety of shapes, including significant blue-shifted components and widths from 200 to 4000 km/s. Overall, 60+-11% appear to be AGN dominated, and at $L_{Ly\alpha}>10^{43.3}$ erg/s and/or $M_{UV}<-21.5$ virtually all LAEs are AGN with high ionisation parameters (log U=0.6+-0.5) and with metallicities of ~0.5-1 Zsun. Those lacking signatures of AGN (40+-11%) have lower ionisation parameters ($\log U=-3.0^{+1.6}_{-0.9}$ and $\log\xi_{\rm ion}=25.4\pm0.2$) and are apparently metal-poor sources likely powered by young, dust-poor "maximal" starbursts. Our results show that luminous LAEs at z~2-3 are a diverse population and that 2xL$^*_{Ly\alpha}$ and 2xM$_{UV}^*$ mark a sharp transition in the nature of LAEs, from star formation dominated to AGN dominated.

We investigate the rest-frame UV morphologies of a large sample of Lyman-a emitters (LAEs) from z~2 to z~6, selected in a uniform way with 16 different narrow- and medium-bands over the full COSMOS field. We use 3045 LAEs with HST coverage in a stacking analysis and find that they have M_UV~-20, below M*_UV at these redshifts. We also focus our analysis on a subsample of 780 individual galaxies with i_AB<25 for which GALFIT converges for 429 of them. The individual median size (re~1 kpc), ellipticities (slightly elongated with (b/a)~0.45), S\'ersic index (disk-like with n<2) and light concentration (comparable to that of disk or irregular galaxies, with C~2.7) of LAEs show mild evolution from z~2 to z~6. LAEs with the highest rest-frame equivalent widths (EW) are the smallest/most compact (re~0.8 kpc, compared to re~1.5 kpc for the lower EW LAEs). When stacking our samples in bins of fixed Lya luminosity and Lya EW we find evidence for redshift evolution in n and C, but not in galaxy sizes. The evolution seems to be stronger for LAEs with 25<EW<100 {\AA}. When compared to other SFGs, LAEs are found to be smaller at all redshifts. The difference between the two populations changes with redshift, from a factor of ~1 at z>5 to SFGs being a factor of ~2-4 larger than LAEs for z<2. This means that at the highest redshifts, where typical sizes approach those of LAEs, the fraction of galaxies showing Lya in emission (and with a high Lya escape fraction) should be much higher, consistent with observations.

We present and explore deep narrow- and medium-band data obtained with the Subaru and the Isaac Newton telescopes in the ~2 deg$^2$ COSMOS field. We use these data as an extremely wide, low-resolution (R~20-80) IFU survey to slice through the COSMOS field and obtain a large sample of ~4000 Lyman-$\alpha$ (Lya) emitters from z~2 to z~6 in 16 redshift slices (SC4K). We present new Lya luminosity functions (LFs) covering a co-moving volume of ~10$^8$Mpc$^3$. SC4K extensively complements ultra-deep surveys, jointly covering over 4 dex in Lya luminosity and revealing a global (2.5<z<6) synergy LF with $\alpha=-1.93\pm0.12$, $\log_{10}\Phi^*=-3.45^{+0.22}_{-0.29}$ Mpc$^{-3}$ and $\log_{10}L^*=42.93^{+0.15}_{-0.11}$ erg/s. The Schechter component of the Lya LF reveals a factor ~5 rise in $L^*$ and a ~7x decline in $\Phi^*$ from z~2 to z~6. The data reveal an extra power-law (or Schechter) component above L~10$^{43.3}$ erg/s at z~2.2-3.5 and we show that it is partially driven by X-ray and radio AGN, as their Lya LF resembles the excess. The power-law component vanishes and/or is below our detection limits above z>3.5, likely linked with the evolution of the AGN population. The Lya luminosity density rises by a factor ~2 from z~2 to z~3 but is then found to be roughly constant (~$1.1\times10^{40}$ erg s$^{-1}$ Mpc$^{-3}$) to z~6, despite the ~0.7 dex drop in UV luminosity density. The Lya/UV luminosity density ratio rises from $4\pm1$% to $30\pm6$% from z~2.2 to z~6. Our results imply a rise of a factor of ~2 in the global ionisation efficiency ($\xi_{\rm ion}$) and a factor ~$4\pm1$ in the Lya escape fraction from z~2 to z~6, hinting for evolution in both the typical burstiness/stellar populations and even more so in the typical ISM conditions allowing Ly$\alpha$ photons to escape.

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 present spectroscopic follow-up observations of CR7 with ALMA, targeted at constraining the infrared (IR) continuum and [CII]$_{158 \mu \rm m}$ line-emission at high spatial resolution matched to the HST/WFC3 imaging. CR7 is a luminous Ly$\alpha$ emitting galaxy at $z=6.6$ that consists of three separated UV-continuum components. Our observations reveal several well-separated components of [CII] emission. The two most luminous components in [CII] coincide with the brightest UV components (A and B), blue-shifted by $\approx 150$ km s$^{-1}$ with respect to the peak of Ly$\alpha$ emission. Other [CII] components are observed close to UV clumps B and C and are blue-shifted by $\approx300$ and $\approx80$ km s$^{-1}$ with respect to the systemic redshift. We do not detect FIR continuum emission due to dust with a 3$\sigma$ limiting luminosity L$_{\rm IR} (T_d = 35 \rm \, K) < 3.1\times10^{10}$ L$_{\odot}$. This allows us to mitigate uncertainties in the dust-corrected SFR and derive SFRs for the three UV clumps A, B and C of 28, 5 and 7 M$_{\odot}$ yr$^{-1}$. All clumps have [CII] luminosities consistent within the scatter observed in the local relation between SFR and L$_{\rm [CII]}$, implying that strong Ly$\alpha$ emission does not necessarily anti-correlate with [CII] luminosity. Combining our measurements with the literature, we show that galaxies with blue UV slopes have weaker [CII] emission at fixed SFR, potentially due to their lower metallicities and/or higher photoionisation. Comparison with hydrodynamical simulations suggests that CR7's clumps have metallicities of $0.1<\rm Z/Z_{\odot}<0.2$. The observed ISM structure of CR7 indicates that we are likely witnessing the build up of a central galaxy in the early Universe through complex accretion of satellites.

We present spectroscopic follow-up of candidate luminous Ly$\alpha$ emitters (LAEs) at $z=5.7-6.6$ in the SA22 field with VLT/X-SHOOTER. We confirm two new luminous LAEs at $z=5.676$ (SR6) and $z=6.532$ (VR7), and also present {\it HST} follow-up of both sources. These sources have luminosities L$_{\rm Ly\alpha} \approx 3\times10^{43}$ erg s$^{-1}$, very high rest-frame equivalent widths of EW$_0\gtrsim 200$ {\AA} and narrow Ly$\alpha$ lines (200-340 km s$^{-1}$). VR7 is the most UV-luminous LAE at $z>6.5$, with M$_{1500} = -22.5$, even brighter in the UV than CR7. Besides Ly$\alpha$, we do not detect any other rest-frame UV lines in the spectra of SR6 and VR7, and argue that rest-frame UV lines are easier to observe in bright galaxies with low Ly$\alpha$ equivalent widths. We confirm that Ly$\alpha$ line-widths increase with Ly$\alpha$ luminosity at $z=5.7$, while there are indications that Ly$\alpha$ lines of faint LAEs become broader at $z=6.6$, potentially due to reionisation. We find a large spread of up to 3 dex in UV luminosity for $>L^{\star}$ LAEs, but find that the Ly$\alpha$ luminosity of the brightest LAEs is strongly related to UV luminosity at $z=6.6$. Under basic assumptions, we find that several LAEs at $z\approx6-7$ have Ly$\alpha$ escape fractions $\gtrsim100$ \%, indicating bursty star-formation histories, alternative Ly$\alpha$ production mechanisms, or dust attenuating Ly$\alpha$ emission differently than UV emission. Finally, we present a method to compute $\xi_{ion}$, the production efficiency of ionising photons, and find that LAEs at $z\approx6-7$ have high values of log$_{10}(\xi_{ion}$/Hz erg$^{-1}) \approx 25.51\pm0.09$ that may alleviate the need for high Lyman-Continuum escape fractions required for reionisation.

While traditionally associated with active galactic nuclei (AGN), the properties of the CII] ($\lambda=2326$\,{\AA}), CIII] ($\lambda,\lambda=1907,1909$\,{\AA}) and CIV ($\lambda,\lambda=1549, 1551$\,{\AA}) emission lines are still uncertain as large, unbiased samples of sources are scarce. We present the first blind, statistical study of CII], CIII] and CIV emitters at $z\sim0.68,1.05,1.53$, respectively, uniformly selected down to a flux limit of $\sim4\times10^{-17}$ erg s$^{-1}$ cm$^{-1}$ through a narrow band survey covering an area of $\sim1.4$ deg$^2$ over COSMOS and UDS. We detect 16 CII], 35 CIII] and 17 CIV emitters, whose nature we investigate using optical colours as well as HST, X-ray, radio and far infra-red data. We find that $z\sim0.7$ CII] emitters are consistent with a mixture of blue (UV slope $\beta=-2.0\pm0.4$) star forming galaxies with disky HST structure and AGN with Seyfert-like morphologies. Bright CII] emitters have individual X-ray detections as well as high average black hole accretion rates (BHAR) of $\sim0.1$ $M_{\odot}$ yr$^{-1}$. CIII] emitters at $z\sim1.05$ trace a general population of SF galaxies, with $\beta=-0.8\pm1.1$, a variety of optical morphologies, including isolated and interacting galaxies and low BHAR ($<0.02$ $M_{\odot}$ yr$^{-1}$). Our CIV emitters at $z\sim1.5$ are consistent with young, blue quasars ($\beta\sim-1.9$) with point-like optical morphologies, bright X-ray counterparts and large BHAR ($0.8$ $M_{\odot}$ yr$^{-1}$). We also find some surprising CII], CIII] and CIV emitters with rest-frame equivalent widths which could be as large as $50-100$ {\AA}. AGN or spatial offsets between the UV continuum stellar disk and the line emitting regions may explain the large EW. These bright CII], CIII] and CIV emitters are ideal candidates for spectroscopic follow up to fully unveil their nature.

Recently, the CIII] and CIV emission lines have been observed in galaxies in the early Universe ($z>5$), providing new ways to measure their redshift and study their stellar populations and AGN. We explore the first blind CII], CIII] and CIV survey ($z\sim0.68, 1.05, 1.53$, respectively) presented in Stroe et al. (2017). We derive luminosity functions (LF) and study properties of CII], CIII] and CIV line emitters through comparisons to the LFs of H$\alpha$ and Ly$\alpha$ emitters, UV selected star forming (SF) galaxies and quasars at similar redshifts. The CII] LF at $z\sim0.68$ is equally well described by a Schechter or a power-law LF, characteristic of a mixture of SF and AGN activity. The CIII] LF ($z\sim1.05$) is consistent to a scaled down version of the Schechter H$\alpha$ and Ly$\alpha$ LF at their redshift, indicating a SF origin. In stark contrast, the CIV LF at $z\sim1.53$ is well fit by a power-law, quasar-like LF. We find that the brightest UV sources ($M_{UV}<-22$) will universally have CIII] and CIV emission. However, on average, CIII] and CIV are not as abundant as H$\alpha$ or Ly$\alpha$ emitters at the same redshift, with cosmic average ratios of $\sim0.02-0.06$ to H$\alpha$ and $\sim0.01-0.1$ to intrinsic Ly$\alpha$. We predict that the CIII] and CIV lines can only be truly competitive in confirming high redshift candidates when the hosts are intrinsically bright and the effective Ly$\alpha$ escape fraction is below 1 per cent. While CIII] and CIV were proposed as good tracers of young, relatively low-metallicity galaxies typical of the early Universe, we find that, at least at $z\sim1.5$, CIV is exclusively hosted by AGN/quasars, especially at large line equivalent widths.

We present a sample of $\sim 1000$ emission line galaxies at $z=0.4-4.7$ from the $\sim0.7$deg$^2$ High-$z$ Emission Line Survey (HiZELS) in the Bo\"otes field identified with a suite of six narrow-band filters at $\approx 0.4-2.1$ $\mu$m. These galaxies have been selected on their Ly$\alpha$ (73), {\sc [Oii]} (285), H$\beta$/{\sc [Oiii]} (387) or H$\alpha$ (362) emission-line, and have been classified with optical to near-infrared colours. A subsample of 98 sources have reliable redshifts from multiple narrow-band (e.g. [O{\sc ii}]-H$\alpha$) detections and/or spectroscopy. In this survey paper, we present the observations, selection and catalogs of emitters. We measure number densities of Ly$\alpha$, [O{\sc ii}], H$\beta$/{\sc [Oiii]} and H$\alpha$ and confirm strong luminosity evolution in star-forming galaxies from $z\sim0.4$ to $\sim 5$, in agreement with previous results. To demonstrate the usefulness of dual-line emitters, we use the sample of dual [O{\sc ii}]-H$\alpha$ emitters to measure the observed [O{\sc ii}]/H$\alpha$ ratio at $z=1.47$. The observed [O{\sc ii}]/H$\alpha$ ratio increases significantly from 0.40$\pm0.01$ at $z=0.1$ to 0.52$\pm0.05$ at $z=1.47$, which we attribute to either decreasing dust attenuation with redshift, or due to a bias in the (typically) fiber-measurements in the local Universe which only measure the central kpc regions. At the bright end, we find that both the H$\alpha$ and Ly$\alpha$ number densities at $z\approx2.2$ deviate significantly from a Schechter form, following a power-law. We show that this is driven entirely by an increasing X-ray/AGN fraction with line-luminosity, which reaches $\approx 100$ \% at line-luminosities $L\gtrsim3\times10^{44}$ erg s$^{-1}$.

We investigate the clustering properties of $\sim 7000$ H$\beta$+[OIII] and [OII] narrowband-selected emitters at $z \sim 0.8 - 4.7$ from the High-$z$ Emission Line Survey. We find clustering lengths, $r_0$, of $1.5 - 4.0h^{-1}$ Mpc and minimum dark matter halo masses of $10^{10.7 - 12.1}\rm{M}_\odot$ for our $z = 0.8 - 3.2$ H$\beta$+[OIII] emitters and $r_0 \sim 2.0 - 8.3h^{-1}$ Mpc and halo masses of $10^{11.5 - 12.6}\rm{M}_\odot$ for our $z = 1.5 - 4.7$ [OII] emitters. We find $r_0$ to strongly increase both with increasing line luminosity and redshift. By taking into account the evolution of the characteristic line luminosity, $L^\star(z)$, and using our model predictions of halo mass given $r_0$, we find a strong, redshift-independent increasing trend between $L/L^\star(z)$ and minimum halo mass. The faintest H$\beta$+[OIII] emitters are found to reside in $10^{9.5}\rm{M}_\odot$ halos and the brightest emitters in $10^{13.0}\rm{M}_\odot$ halos. For [OII] emitters, the faintest emitters are found in $10^{10.5} \rm{M}_\odot$ halos and the brightest emitters in $10^{12.6}\rm{M}_\odot$ halos. A redshift-independent stellar mass dependency is also observed where the halo mass increases from $10^{11}\rm{M}_\odot$ to $10^{12.5} \rm{M}_\odot$ for stellar masses of $10^{8.5}\rm{M}_\odot$ to $10^{11.5}\rm{M}_\odot$, respectively. We investigate the interdependencies of these trends by repeating our analysis in a $L_\textrm{line} - \rm{M}_\textrm{star}$ grid space for our most populated samples (H$\beta$+[OIII] $z = 0.84$ and [OII] $z = 1.47$) and find that the line luminosity dependency is stronger than the stellar mass dependency on halo mass. For $L > L^\star$ emitters at all epochs, we find a relatively flat trend with halo masses of $10^{12.5 - 13}\rm{M}_\odot$ which may be due to quenching mechanisms in massive halos which is consistent with a transitional halo mass predicted by models.

We present dynamical measurements from the KMOS (K-band Multi-Object Spectrograph) Deep Survey (KDS), which is comprised of 77 typical star-forming galaxies at z $\simeq$ 3.5 in the mass range 9.0 < log(M$_{\star}$/M$_{\odot}$) < 10.5. These measurements constrain the internal dynamics, the intrinsic velocity dispersions ({\sigma}$_{int}$) and rotation velocities (V$_{C}$) of galaxies in the high redshift Universe. The mean velocity dispersion of the galaxies in our sample is {\sigma}$_{int}$ = $70.8^{+3.3}_{-3.1}$ km s$^{-1}$, revealing that the increasing average {\sigma}$_{int}$ with increasing redshift, reported for z $\lesssim2$, continues out to z $\simeq$ 3.5. Only 34 $\pm$ 8% of our galaxies are rotation-dominated (V$_{C}$/{\sigma}$_{int}$ > 1), with the sample average V$_{C}$/{\sigma}$_{int}$ value much smaller than at lower redshift. After carefully selecting comparable star-forming samples at multiple epochs, we find that the rotation-dominated fraction evolves with redshift with a z$^{-0.2}$ dependence. The rotation-dominated KDS galaxies show no clear offset from the local rotation velocity-stellar mass (i.e. V$_{C}$-M$_{\star}$) relation, although a smaller fraction of the galaxies are on the relation due to the increase in the dispersion-dominated fraction. These observations are consistent with a simple equilibrium model picture, in which random motions are boosted in high redshift galaxies by a combination of the increasing gas fractions, accretion efficiency, specific star-formation rate and stellar feedback and which may provide significant pressure support against gravity on the galactic disk scale.

We present the CAlibrating LYMan-$\alpha$ with H$\alpha$ (CALYMHA) pilot survey and new results on Lyman-$\alpha$ (Lya) selected galaxies at z~2. We use a custom-built Lya narrow-band filter at the Isaac Newton Telescope, designed to provide a matched volume coverage to the z=2.23 Ha HiZELS survey. Here we present the first results for the COSMOS and UDS fields. Our survey currently reaches a 3$\sigma$ line flux limit of ~4x10$^{-17}$ erg/s/cm$^{2}$, and a Lya luminosity limit of ~10$^{42.3}$ erg/s. We find 188 Lya emitters over 7.3x10$^5$ Mpc$^{3}$, but also find significant numbers of other line emitting sources corresponding to HeII, CIII] and CIV emission lines. These sources are important contaminants, and we carefully remove them, unlike most previous studies. We find that the Lya luminosity function at z=2.23 is very well described by a Schechter function up to L~10$^{43}$ erg/s with L$^*=10^{42.59+-0.05}$ erg/s, $\phi^*=10^{-3.09+-0.08}$ Mpc$^{-3}$ and $\alpha$=-1.75+-0.15. Above L~10$^{43}$ erg/s the Lya luminosity function becomes power-law like, driven by X-ray AGN. We find that Lya-selected emitters have a high escape fraction of 37+-7%, anti-correlated with Lya luminosity and correlated with Lya equivalent width. Lya emitters have ubiquitous large (~40 kpc) Lya haloes, 2x larger than their Ha extents. By directly comparing our Lya and Ha luminosity functions we find that the global/overall escape fraction of Lya photons (within a 13 kpc radius) from the full population of star-forming galaxies is 5.1+-0.2% at the peak of the star formation history. An extra 3.3+-0.3% of Lya photons likely still escape, but at larger radii.

We study the production rate of ionizing photons of a sample of 588 H$\alpha$ emitters (HAEs) and 160 Lyman-$\alpha$ emitters (LAEs) at $z=2.2$ in the COSMOS field in order to assess the implied emissivity from galaxies, based on their UV luminosity. By exploring the rest-frame Lyman Continuum (LyC) with GALEX/$NUV$ data, we find f$_{\rm esc} < 2.8\, (6.4)$% through median (mean) stacking. By combining the H$\alpha$ luminosity density with IGM emissivity measurements from absorption studies, we find a globally averaged $\langle$f$_{\rm esc}\rangle$ of $5.9^{+14.5}_{-4.2}$ % at $z=2.2$ if we assume HAEs are the only source of ionizing photons. We find similarly low values of the global $\langle$f$_{\rm esc}\rangle$ at $z\approx3-5$, also ruling out a high $\langle$f$_{\rm esc}\rangle$ at $z<5$. These low escape fractions allow us to measure $\xi_{ion}$, the number of produced ionizing photons per unit UV luminosity, and investigate how this depends on galaxy properties. We find a typical $\xi_{ion} \approx 10^{24.77\pm0.04}$ Hz erg$^{-1}$ for HAEs and $\xi_{ion} \approx 10^{25.14\pm0.09}$ Hz erg$^{-1}$ for LAEs. LAEs and low mass HAEs at $z=2.2$ show similar values of $\xi_{ion}$ as typically assumed in the reionization era, while the typical HAE is three times less ionizing. Due to an increasing $\xi_{ion}$ with increasing EW(H$\alpha$), $\xi_{ion}$ likely increases with redshift. This evolution alone is fully in line with the observed evolution of $\xi_{ion}$ between $z\approx2-5$, indicating a typical value of $\xi_{ion} \approx 10^{25.4}$ Hz erg$^{-1}$ in the reionization era.

We use the hydrodynamical EAGLE simulation to study the magnitude and origin of the scatter in the stellar mass - halo mass relation for central galaxies. We separate cause and effect by correlating stellar masses in the baryonic simulation with halo properties in a matched dark matter only (DMO) simulation. The scatter in stellar mass increases with redshift and decreases with halo mass. At $z = 0.1$ it declines from 0.25 dex at $M_{200, \rm DMO} \approx 10^{11}$ M$_{\odot}$ to 0.12 dex at $M_{200, \rm DMO} \approx 10^{13}$ M$_{\odot}$, but the trend is weak above $10^{12}$ M$_{\odot}$. For $M_{200, \rm DMO} < 10^{12.5}$ M$_{\odot}$ up to 0.04 dex of the scatter is due to scatter in the halo concentration. At fixed halo mass, a larger stellar mass corresponds to a more concentrated halo. This is likely because higher concentrations imply earlier formation times and hence more time for accretion and star formation, and/or because feedback is less efficient in haloes with higher binding energies. The maximum circular velocity, $V_{\rm max, DMO}$, and binding energy are therefore more fundamental properties than halo mass, meaning that they are more accurate predictors of stellar mass, and we provide fitting formulae for their relations with stellar mass. However, concentration alone cannot explain the total scatter in the $M_{\rm star} - M_{200, \rm DMO}$ relation, and it does not explain the scatter in $M_{\rm star} -V_{\rm max, DMO}$. Halo spin, sphericity, triaxiality, substructure and environment are also not responsible for the remaining scatter, which thus could be due to more complex halo properties or non-linear/stochastic baryonic effects.

We present new results from the widest narrow band survey search for Lyman-alpha (Lya) emitters at z=5.7, just after reionization. We survey a total of 7 deg$^2$ spread over the COSMOS, UDS and SA22 fields. We find over 11,000 line emitters, out of which 514 are robust Lya candidates at z=5.7 within a volume of 6.3x10$^6$ Mpc$^3$. Our Lya emitters span a wide range in Lya luminosities, from faint to bright (L$_{\rm Ly\alpha}\sim10^{42.5-44}$ erg s$^{-1}$) and rest-frame equivalent widths (EW$_0$~25-1000 \AA) in a single, homogeneous data-set. By combining all our fields we find that the faint end slope of the z=5.7 Lya luminosity function is very steep, with $\alpha=-2.3^{+0.4}_{-0.3}$. We also present an updated z=6.6 Lya luminosity function, based on comparable volumes and obtained with the same methods, which we directly compare with that at z=5.7. We find a significant decline of the number density of faint Lya emitters from z=5.7 to z=6.6 (by $0.5\pm0.1$ dex), but no evolution at the bright end/no evolution in L*. Faint Lya emitters at z=6.6 show much more extended haloes than those at z=5.7, suggesting that neutral Hydrogen plays an important role, increasing the scattering and leading to observations missing faint Lya emission within the epoch of reionization. All together, our results suggest that we are observing patchy reionization which happens first around the brightest Lya emitters, allowing the number densities of those sources to remain unaffected by the increase of neutral Hydrogen fraction from z~5 to z~7.

We present the first results from our CAlibrating LYMan-$\alpha$ with H$\alpha$ (CALYMHA) pilot survey at the Isaac Newton Telescope. We measure Ly$\alpha$ emission for 488 H$\alpha$ selected galaxies at $z=2.23$ from HiZELS in the COSMOS and UDS fields with a specially designed narrow-band filter ($\lambda_c$ = 3918 {\AA}, $\Delta\lambda$= 52 {\AA}). We find 17 dual H$\alpha$-Ly$\alpha$ emitters ($f_{\rm Ly\alpha} >5\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, of which 5 are X-ray AGN). For star-forming galaxies, we find a range of Ly$\alpha$ escape fractions (f$_{\rm esc}$, measured with 3$"$ apertures) from $2$\%$-30$\%. These galaxies have masses from $3\times10^8$ M$_{\odot}$ to 10$^{11}$ M$_{\odot}$ and dust attenuations E$(B-V)=0-0.5$. Using stacking, we measure a median escape fraction of $1.6\pm0.5$\% ($4.0\pm1.0$\% without correcting H$\alpha$ for dust), but show that this depends on galaxy properties. The stacked f$_{\rm esc}$ tends to decrease with increasing SFR and dust attenuation. However, at the highest masses and dust attenuations, we detect individual galaxies with f$_{\rm esc}$ much higher than the typical values from stacking, indicating significant scatter in the values of f$_{\rm esc}$. Relations between f$_{\rm esc}$ and UV slope are bimodal, with high f$_{\rm esc}$ for either the bluest or reddest galaxies. We speculate that this bimodality and large scatter in the values of f$_{\rm esc}$ is due to additional physical mechanisms such as outflows facilitating f$_{\rm esc}$ for dusty/massive systems. Ly$\alpha$ is significantly more extended than H$\alpha$ and the UV. f$_{\rm esc}$ continues to increase up to at least 20 kpc (3$\sigma$, 40 kpc [2$\sigma$]) for typical SFGs and thus the aperture is the most important predictor of f$_{\rm esc}$.

We investigate correlations between different physical properties of star-forming galaxies in the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) cosmological hydrodynamical simulation suite over the redshift range $0\le z\le 4.5$. A principal component analysis reveals that neutral gas fraction ($f_{\rm gas, neutral}$), stellar mass ($M_{\rm stellar}$) and star formation rate (SFR) account for most of the variance seen in the population, with galaxies tracing a two-dimensional, nearly flat, surface in the three-dimensional space of $f_{\rm gas, neutral}-M_{\rm stellar}-\rm SFR$ with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their $f_{\rm gas, neutral}$ and SFR drop with redshift. The positions of galaxies along the plane are highly correlated with gas metallicity. The metallicity can therefore be robustly predicted from $f_{\rm gas, neutral}$, or from the $M_{\rm stellar}$ and SFR. We argue that the appearance of this "fundamental plane of star formation" is a consequence of self-regulation, with the plane's curvature set by the dependence of the SFR on gas density and metallicity. We analyse a large compilation of observations spanning the redshift range $0\lesssim \rm z\lesssim 2.5$, and find that such a plane is also present in the data. The properties of the observed fundamental plane of star formation are in good agreement with EAGLE's predictions.

We use new near-infrared spectroscopic observations to investigate the nature and evolution of the most luminous H\alpha (Ha) emitters at z~0.8-2.23, which evolve strongly in number density over this period, and compare them to more typical Ha emitters. We study 59 luminous Ha emitters with $L_{H\alpha}>L_{H\alpha}^*$, roughly equally split per redshift slice at z~0.8, 1.47 and 2.23 from the HiZELS and CF-HiZELS surveys. We find that, overall, 30$\pm$8% are AGN (80$\pm$30% of these AGN are broad-line AGN, BL-AGN), and we find little to no evolution in the AGN fraction with redshift, within the errors. However, the AGN fraction increases strongly with Ha luminosity and correlates best with $L_{H\alpha}/L_{H\alpha}^*(z)$. While $L_{H\alpha}<L_{\rm H\alpha}^*(z)$ Ha emitters are largely dominated by star-forming galaxies (>80%), the most luminous Ha emitters ($L_{H\alpha}>10L_{H\alpha}^*(z)$) at any cosmic time are essentially all BL-AGN. Using our AGN-decontaminated sample of luminous star-forming galaxies, and integrating down to a fixed Ha luminosity, we find a factor of ~1300x evolution in the star formation rate density from z=0 to z=2.23. This is much stronger than the evolution from typical Ha star-forming galaxies and in line with the evolution seen for constant luminosity cuts used to select "Ultra-Luminous" Infrared Galaxies and/or sub-millimetre galaxies. By taking into account the evolution in the typical Ha luminosity, we show that the most strongly star-forming Ha-selected galaxies at any epoch ($L_{H\alpha}>L^*_{H\alpha}(z)$) contribute the same fractional amount of ~15% to the total star-formation rate density, at least up to z=2.23.

CR7 is the brightest $z=6.6 \, {\rm Ly}\alpha$ emitter (LAE) known to date, and spectroscopic follow-up by Sobral et al. (2015) suggests that CR7 might host Population (Pop) III stars. We examine this interpretation using cosmological hydrodynamical simulations. Several simulated galaxies show the same "Pop III wave" pattern observed in CR7. However, to reproduce the extreme CR7 ${\rm Ly}\alpha$/HeII1640 line luminosities ($L_{\rm \alpha/He II}$) a top-heavy IMF and a massive ($>10^{7}{\rm M}_{\odot}$) PopIII burst with age $<2$ Myr are required. Assuming that the observed properties of ${\rm Ly}\alpha$ and HeII emission are typical for Pop III, we predict that in the COSMOS/UDS/SA22 fields, 14 out of the 30 LAEs at $z=6.6$ with $L_{\alpha} >10^{43.3}{\rm erg}\,{\rm s}^{-1}$ should also host Pop III stars producing an observable $L_{\rm He II}>10^{42.7}{\rm erg}\,{\rm s}^{-1}$. As an alternate explanation, we explore the possibility that CR7 is instead powered by accretion onto a Direct Collapse Black Hole (DCBH). Our model predicts $L_{\alpha}$, $L_{\rm He II}$, and X-ray luminosities that are in agreement with the observations. In any case, the observed properties of CR7 indicate that this galaxy is most likely powered by sources formed from pristine gas. We propose that further X-ray observations can distinguish between the two above scenarios.

Faint Lyman-$\alpha$ (Ly$\alpha$) emitters become increasingly rarer towards the re-ionisation epoch (z~6-7). However, observations from a very large (~5deg$^2$) Ly$\alpha$ survey at z=6.6 (Matthee et al. 2015) show that this is not the case for the most luminous emitters. Here we present follow-up observations of the two most luminous z~6.6 Ly$\alpha$ candidates in the COSMOS field: `MASOSA' and `CR7'. We used X-SHOOTER, SINFONI and FORS2 (VLT), and DEIMOS (Keck), to confirm both candidates beyond any doubt. We find redshifts of z=6.541 and z=6.604 for MASOSA and CR7, respectively. MASOSA has a strong detection in Ly$\alpha$ with a line width of $386\pm30$ km/s (FWHM) and with high EW$_0$ (>200 \AA), but it is undetected in the continuum. CR7, with an observed Ly$\alpha$ luminosity of $10^{43.93\pm0.05}$erg/s is the most luminous Ly$\alpha$ emitter ever found at z>6. CR7 reveals a narrow Ly$\alpha$ line with $266\pm15$ km/s FWHM, being detected in the NIR (rest-frame UV, with $\beta=-2.3\pm0.1$) with an excess in $J$, and also strongly detected in IRAC/Spitzer. We detect a narrow HeII1640$\AA$ emission line ($6\sigma$) which explains the excess seen in the $J$ band photometry (EW$_0$~80 \AA). We find no other emission lines from the UV to the NIR in our X-SHOOTER spectra (HeII/OIII]1663A>3 and HeII/CIII]1908A>2.5). We find that CR7 is best explained by a combination of a PopIII-like population which dominates the rest-frame UV and the nebular emission, and a more normal stellar population which dominates the mass. HST/WFC3 observations show that the light is indeed spatially separated between a very blue component, coincident with Ly$\alpha$ and HeII emission, and two red components (~5 kpc away), which dominate the mass. Our findings are consistent with theoretical predictions of a PopIII wave, with PopIII star formation migrating away from the original sites of star formation.

Using wide field narrow-band surveys, we provide a new measurement of the $z=6.6$ Lyman-$\alpha$ Emitter (LAE) luminosity function (LF), which constraints the bright end for the first time. We use a combination of archival narrow-band NB921 data in UDS and new NB921 measurements in SA22 and COSMOS/UltraVISTA, all observed with the Subaru telescope, with a total area of $\sim 5$ deg$^2$. We exclude lower redshift interlopers by using broad-band optical and near-infrared photometry and also exclude three supernovae with data split over multiple epochs. Combining the UDS and COSMOS samples we find no evolution of the bright end of the Ly$\alpha$ LF between $z=5.7$ and $6.6$, which is supported by spectroscopic follow-up, and conclude that sources with \emph{Himiko}-like luminosity are not as rare as previously thought, with number densities of $\sim 1.5\times10^{-5}$ Mpc$^{-3}$. Combined with our wide-field SA22 measurements, our results indicate a non-Schechter-like bright end of the LF at $z=6.6$ and a different evolution of \emph{observed} faint and bright LAEs. This differential evolution is also seen in the spectroscopic follow-up of UV selected galaxies and is now also confirmed for Ly$\alpha$ emitters, and we argue that it may be an effect of re-ionisation. Using a toy-model, we show that such differential evolution of the LF is expected, since brighter sources are able to ionise their surroundings earlier, such that Ly$\alpha$ photons are able to escape. Our targets are excellent candidates for detailed follow-up studies and provide the possibility to give a unique view on the earliest stages in the formation of galaxies and re-ionisation process.

We present results from the largest contiguous narrow-band survey in the near-infrared. We have used WIRCam/CFHT and the lowOH2 filter (1.187$\pm$0.005 $\mu$m) to survey ~10 deg$^2$ of contiguous extragalactic sky in the SA22 field. A total of ~6000 candidate emission-line galaxies are found. We use deep ugrizJK data to obtain robust photometric redshifts. We combine our data with the High-redshift Emission Line Survey (HiZELS), explore spectroscopic surveys (VVDS, VIPERS) and obtain our own spectroscopic follow-up with KMOS, FMOS and MOSFIRE to derive large samples of high-redshift emission-line selected galaxies: 3471 H$\alpha$ emitters at z=0.8, 1343 [OIII]+H$\beta$ emitters at z=1.4 and 572 [OII] emitters at z=2.2. We probe co-moving volumes of >10$^6$ Mpc$^3$ and find significant over-densities, including an 8.5$\sigma$ (spectroscopically confirmed) over-density of H$\alpha$ emitters at z=0.81. We derive H$\alpha$, [OIII]+H$\beta$ and [OII] luminosity functions at z=0.8,1.4,2.2, respectively, and present implications for future surveys such as Euclid. Our uniquely large volumes/areas allow us to sub-divide the samples in thousands of randomised combinations of areas and provide a robust empirical measurement of sample/cosmic variance. We show that surveys for star-forming/emission-line galaxies at a depth similar to ours can only overcome cosmic-variance (errors <10%) if they are based on volumes >5x10$^{5}$ Mpc$^{3}$; errors on $L^*$ and $\phi^*$ due to sample (cosmic) variance on surveys probing ~10$^4$ Mpc$^{3}$ and ~10$^5$ Mpc$^{3}$ are typically very high: ~300% and ~40-60%, respectively.

We have observed a sample of typical z=1 star forming galaxies, selected from the HiZELS survey, with the new KMOS near-infrared, multi-IFU instrument on the VLT, in order to obtain their dynamics and metallicity gradients. The majority of our galaxies have a metallicity gradient consistent with being flat or negative (i.e. higher metallicity cores than outskirts). Intriguingly, we find a trend between metallicity gradient and specific star formation rate (sSFR), such that galaxies with a high sSFR tend to have relatively metal-poor centres, a result which is strengthened when combined with datasets from the literature. This result appears to explain the discrepancies reported between different high redshift studies and varying claims for evolution. From a galaxy evolution perspective, the trend we see would mean that a galaxy's sSFR is governed by the amount of metal poor gas that can be funnelled into its core, triggered either by merging or through efficient accretion. In fact merging may play a significant role as it is the starburst galaxies at all epochs, which have the more positive metallicity gradients. Our results may help to explain the origin of the fundamental metallicity relation, in which galaxies at a fixed mass are observed to have lower metallicities at higher star formation rates, especially if the metallicity is measured in an aperture encompassing only the central regions of the galaxy. Finally, we note that this study demonstrates the power of KMOS as an efficient instrument for large scale resolved galaxy surveys.

Candidate galaxies at redshifts of $z \sim 10$ are now being found in extremely deep surveys, probing very small areas. As a consequence, candidates are very faint, making spectroscopic confirmation practically impossible. In order to overcome such limitations, we have undertaken the CF-HiZELS survey, which is a large area, medium depth near infrared narrow-band survey targeted at $z=8.8$ Lyman-$\alpha$ (Ly$\alpha$) emitters (LAEs) and covering 10 deg$^2$ in part of the SSA22 field with the Canada-France-Hawaii Telescope. We surveyed a comoving volume of $4.7\times 10^6$ Mpc$^3$ to a Ly$\alpha$ luminosity limit of $6.3\times10^{43}$ erg s$^{-1}$. We look for Ly$\alpha$ candidates by applying the following criteria: i) clear emission line source, ii) no optical detections ($ugriz$ from CFHTLS), iii) no visible detection in the optical stack ($ugriz > 27$), iv) visually checked reliable NB$_J$ and $J$ detections and v) $J-K \leq 0$. We compute photometric redshifts and remove a significant amount of dusty lower redshift line-emitters at $z \sim 1.4 $ or $2.2$. A total of 13 Ly$\alpha$ candidates were found, of which two are marked as strong candidates, but the majority have very weak constraints on their SEDs. Using follow-up observations with SINFONI/VLT we are able to exclude the most robust candidates as Ly$\alpha$ emitters. We put a strong constraint on the Ly$\alpha$ luminosity function at $z \sim 9$ and make realistic predictions for ongoing and future surveys. Our results show that surveys for the highest redshift LAEs are susceptible of multiple contaminations and that spectroscopic follow-up is absolutely necessary.

We present the spatially resolved H-alpha (Ha) dynamics of sixteen star-forming galaxies at z~0.81 using the new KMOS multi-object integral field spectrograph on the ESO VLT. These galaxies were selected using 1.18 um narrow-band imaging from the 10 deg^2 CFHT-HiZELS survey of the SA22hr field, are found in a ~4Mpc over-density of Ha emitters and likely reside in a group/intermediate environment, but not a cluster. We confirm and identify a rich group of star-forming galaxies at z=0.813+-0.003, with thirteen galaxies within 1000 km/s of each other, and 7 within a diameter of 3Mpc. All our galaxies are "typical" star-forming galaxies at their redshift, 0.8+-0.4 SFR*(z=0.8), spanning a range of specific star formation rate of sSFR=0.2-1.1 Gyr^-1 and have a median metallicity very close to solar of 12+log(O/H)=8.62+-0.06. We measure the spatially resolved Ha dynamics of the galaxies in our sample and show that thirteen out of sixteen galaxies can be described by rotating disks and use the data to derive inclination corrected rotation speeds of 50-275 km/s. The fraction of disks within our sample is 75+-8, consistent with previous results based on HST morphologies of Ha selected galaxies at z~1 and confirming that disks dominate the star formation rate density at z~1. Our Ha galaxies are well fitted by the z~1-2 Tully-Fisher relation, confirming the evolution seen in the zero-point. Apart from having, on average, higher stellar masses and lower sSFRs, our group galaxies at z=0.813 present the same mass-metallicity and TF relation as z~1 field galaxies, and are all disk galaxies.