Search

We present the first direct spectroscopic measurement of the stellar velocity dispersion function (VDF) for massive quiescent and star-forming galaxies at $0.6 < z \leq 1.0$. For this analysis we use individual measurements of stellar velocity dispersion from high-S/N spectra from the public Large Early Galaxy Astrophysics Census (LEGA-C) survey. We report a remarkable stability of the VDF for both quiescent and star-forming galaxies within this redshift range, though we note the presence of weak evolution in the number densities of star-forming galaxies. We compare both VDFs with previous direct and inferred measurements at local and intermediate redshifts, with the caveat that previous measurements of the VDF for star-forming galaxies are poorly constrained at all epochs. We emphasize that this work is the first to directly push to low-stellar velocity dispersion ($σ_\star > 100$ km s$^{-1}$) and extend to star-forming galaxies. We are largely consistent with the high-sigma tail measured from BOSS, and we find that the VDF remains constant from the median redshift of LEGA-C, $z\sim0.8$, to the present day., 11 pages, 3 figures, accepted for publication in ApJ

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) is expected to map thousands of square degrees of the northern sky with 56 narrowband filters in the upcoming years. This will make J-PAS a very competitive and unbiased emission line survey compared to spectroscopic or narrowband surveys with fewer filters. The miniJPAS survey covered 1 deg$^2$, and it used the same photometric system as J-PAS, but the observations were carried out with the pathfinder J-PAS camera. In this work, we identify and characterize the sample of emission line galaxies (ELGs) from miniJPAS with a redshift lower than $0.35$. Using a method based on artificial neural networks, we detect the ELG population and measure the equivalent width and flux of the $Hα$, $Hβ$, [OIII], and [NII] emission lines. We explore the ionization mechanism using the diagrams [OIII]/H$β$ versus [NII]/H$α$ (BPT) and EW(H$α$) versus [NII]/H$α$ (WHAN). We identify 1787 ELGs ($83$%) from the parent sample (2154 galaxies) in the AEGIS field. For the galaxies with reliable EW values that can be placed in the WHAN diagram (2000 galaxies in total), we obtained that $72.8 \pm 0.4$%, $17.7 \pm 0.4$% , and $9.4 \pm 0.2$% are star-forming (SF), active galactic nucleus (Seyfert), and quiescent galaxies, respectively. Based on the flux of $Hα$ we find that the star formation main sequence is described as $\log$ SFR $[M_\mathrm{\odot} \mathrm{yr}^{-1}] = 0.90^{+ 0.02}_{-0.02} \log M_{\star} [M_\mathrm{\odot}] -8.85^{+ 0.19}_{-0.20}$ and has an intrinsic scatter of $0.20^{+ 0.01}_{-0.01}$. The cosmic evolution of the SFR density ($ρ_{\text{SFR}}$) is derived at three redshift bins: $0 < z \leq 0.15$, $0.15 < z \leq 0.25$, and $0.25 < z \leq 0.35$, which agrees with previous results that were based on measurements of the $Hα$ emission line., 22 pages, 19 figures

We present the third and final data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO/VLT public spectroscopic survey targeting $0.6 < z < 1.0$, Ks-selected galaxies. The data release contains 3528 spectra with measured stellar velocity dispersions and stellar population properties, a 25-fold increase in sample size compared to previous work. This $K_s$-selected sample probes the galaxy population down to $\sim0.3 L^*$, for all colors and morphological types. Along with the spectra we publish a value-added catalog with stellar and ionized gas velocity dispersions, stellar absorption line indices, emission line fluxes and equivalent widths, complemented with structural parameters measured from HST/ACS imaging. With its combination of high precision and large sample size, LEGA-C provides a new benchmark for galaxy evolution studies., Accepted for publication in ApJS. Data available at http://archive.eso.org/cms/eso-archive-news/Third-and-final-release-of-the-Large-Early-Galaxy-Census-LEGA-C-Spectroscopic-Public-Survey-published.html or https://users.ugent.be/\string~avdrwel/research.html\#legac

We present the first results from the X-SHOOTER Lyman-$��$ survey at $z=2$ (XLS-$z2$). XLS-$z2$ is a deep spectroscopic survey of 35 Lyman-$��$ emitters (LAEs) utilising $\approx90$ hours of exposure time with VLT/X-SHOOTER and covers rest-frame Ly$��$ to H$��$ emission with R$\approx4000$. We present the sample selection, the observations and the data reduction. Systemic redshifts are measured from rest-frame optical lines for 33/35 sources. In the stacked spectrum, our LAEs are characterised by an interstellar medium with little dust, a low metallicity and a high ionisation state. The ionising sources are young hot stars that power strong emission-lines in the optical and high ionisation lines in the UV. The LAEs exhibit clumpy UV morphologies and have outflowing kinematics with blue-shifted SiII absorption, a broad [OIII] component and a red-skewed Ly$��$ line. Typically 30 % of the Ly$��$ photons escape, of which one quarter on the blue side of the systemic velocity. A fraction of Ly$��$ photons escapes directly at the systemic suggesting clear channels enabling a $\approx10$ % escape of ionising photons, consistent with an inference based on MgII. A combination of a low effective HI column density, a low dust content and young star-burst determine whether a star forming galaxy is observed as a LAE. The first is possibly related to outflows and/or a fortunate viewing angle, while we find that the latter two in LAEs are typical for their stellar mass of 10$^9$ M$_{\odot}$., Main text 26 pages. Accepted for publication in MNRAS

Galaxy evolution is driven by a variety of physical processes which are predicted to proceed at different rates for different dark matter haloes and environments across cosmic times. A record of this evolution is preserved in galaxy stellar populations, which we can access using absorption-line spectroscopy. Here we explore the large LEGA-C survey (DR3) to investigate the role of the environment and stellar mass on stellar populations at z~0.6-1.0 in the COSMOS field. Leveraging the statistical power and depth of LEGA-C, we reveal significant gradients in D4000 and H-delta equivalent widths (EWs) distributions over the stellar mass vs environment 2D spaces for the massive galaxy population (M>10^10 M$_{\odot}$) at z~0.6-1.0. D4000 and H-delta EWs primarily depend on stellar mass, but they also depend on environment at fixed stellar mass. By splitting the sample into centrals and satellites, and in terms of star-forming galaxies and quiescent galaxies, we reveal that the significant environmental trends of D4000 and H-delta EW when controlling for stellar mass are driven by quiescent galaxies. Regardless of being centrals or satellites, star-forming galaxies reveal D4000 and H-delta EWs which depend strongly on their stellar mass and are completely independent of the environment at 0.6, Comment: Accepted for publication in ApJ. 18 pages, 9 figures

Using deep rest-frame optical spectroscopy from the Large Early Galaxy Astrophysical Census (LEGA-C) survey, conducted using VIMOS on the ESO Very Large Telescope, we systematically search for low-ionization [OII] 3726,3729 emission in the spectra of a mass-complete sample of z~0.85 galaxies. Intriguingly, we find that 59 percent of UVJ-quiescent (i.e. non star-forming) galaxies in the sample have ionized gas, as traced by [OII] emission, detected above our completeness limit of 1.5 Angstroms. The median stacked spectrum of the lowest equivalent width quiescent galaxies also shows [OII] emission. The overall fraction of sources with [OII] above our equivalent width limit is comparable to what we find in the low-redshift Universe from GAMA and MASSIVE, except perhaps at the highest stellar masses (log Mstar/Msol > 11.5). However, stacked spectra for the individual low-equivalent width systems uniquely indicates ubiquitous [OII] emission in the higher-z LEGA-C sample, with typical [OII] luminosities per unit stellar mass that are a factor of 3 larger than the lower-z GAMA sample. Star formation in these otherwise quiescent galaxies could play a role in producing the [OII] emission at higher-z, although it is unlikely to provide the bulk of the ionizing photons. More work is required to fully quantify the contributions of evolved stellar populations or active galactic nuclei to the observed spectra., ApJ accepted: 16 pages, 11 figures

The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here we introduce an emissivity framework based on Ly$α$ emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between the production and escape of Ly$α$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Lyα}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$α$ profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at $z\approx2-6$, which is non-trivial given the strong evolution in both the star-formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between $z\approx6-9$ under plausible assumptions. Within this framework, the $>10\times$ rise in the UV population-averaged $f_{\rm{esc}}$ between $z\approx3-7$ naturally arises due to the same phenomena that drive the growing Ly$α$ emitter fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm{esc}}=50\%$, $ξ_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), a highly ionizing minority of galaxies with $M_{\rm UV}, Accepted for publication in MNRAS. Figure 2 shows the main result -- the comoving emissivity due to bright LAEs. Our fiducial model is based on results presented in our companion paper -- Naidu & Matthee et al. 2022, arXiv: 2110.11961

The ionizing photon escape fraction (LyC $f_{\rm{esc}}$) of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshift. Here we circumvent this uncertainty by inferring $f_{\rm{esc}}$ with resolved (R>4000) LyA profiles from the X-SHOOTER LyA survey at z=2 (XLS-z2). We select leakers ($f_{\rm{esc}}>20$%) and non-leakers ($f_{\rm{esc}}0.2 L^{*}$ LyA emitters (LAEs). With median stacked spectra of these subsets covering 1000-8000 �� (rest-frame) we investigate the conditions for LyC $f_{\rm{esc}}$. We find the following differences between leakers vs. non-leakers: (i) strong nebular CIV and HeII emission vs. non-detections, (ii) O32~8.5 vs. ~3, (iii) Ha/Hb indicating no dust vs. E(B-V)~0.3, (iv) MgII emission close to the systemic velocity vs. redshifted, optically thick MgII, (v) LyA $f_{\rm{esc}}$ of ~50% vs. ~10%. The extreme EWs in leakers (O3+Hb~1100 ��) constrain the characteristic timescale of LyC escape to ~3-10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers -- extremely ionizing stellar populations, low column densities, a dust-free, high ionization state ISM -- occur simultaneously in the $f_{\rm{esc}}>20\%$ stack, suggesting they are causally connected, and motivating why indicators like O32 may suffice to constrain $f_{\rm{esc}}$ at z>6 with JWST. The leakers comprise half our sample, have a median LyC $f_{\rm{esc}}$~50%, and an ionising production efficiency $\log({��_{\rm{ion}}/\rm{Hz\ erg^{-1}}})$~25.9. These results show LAEs -- the type of galaxies rare at z=2, but that become the norm at higher redshift -- are highly efficient ionizers, with extreme $��_{\rm{ion}}$ and prolific $f_{\rm{esc}}$ occurring in sync. (ABRIDGED), Submitted to MNRAS. Figure 4 and Table 2 show key results from the stacks. Figure 8 presents a unifying scheme for the LyC duty cycle. Our companion paper (Matthee & Naidu et al. 2021) works out implications of these results for the cosmic UV background at $z\approx2-8$. Comments warmly welcomed and greatly appreciated!

We present the photometric determination of the bright-end (L_Lya>10^43.5 erg/s) of the Lya luminosity function (LF) within four redshifts windows in the interval 2.210^43.3 erg/s) Lya-emitting candidates. We test our selection with two spectroscopic follow-up programs at the GTC telescope, confirming ~89% of the targets as line-emitting sources, with ~64% being genuine z~2.2 QSOs. We extend the 2.22 AGN/QSOs (likely dominant in our samples) are described by a structurally different LF than z>2 star-forming LAEs, namely with L^*_QSOs ~ 100 L^*_LAEs and Phi^*_QSOs ~ 10^-3 Phi^*_LAEs. Finally, our method identifies very efficiently as high-z line-emitters sources without previous spectroscopic confirmation, currently classified as stars (~2000 objects in each redshift bin, on average). Assuming a large predominance of Lya-emitting AGN/QSOs in our samples, this supports the scenario by which these are the most abundant class of z>2 Lya emitters at L_Lya>10^43.3 erg/s., 27 pages (6 in appendix), 21 figures (5 in appendix), submitted to A&A

CR7 is among the most luminous Lyman-$��$ emitters (LAEs) known at $z = 6.6$ and consists of at least three UV components that are surrounded by Lyman-$��$ (Ly$��$) emission. Previous studies have suggested that it may host an extreme ionising source. Here, we present deep integral field spectroscopy of CR7 with VLT/MUSE. We measure extended emission with a similar halo scale length as typical LAEs at $z\approx5$. CR7's Ly$��$ halo is clearly elongated along the direction connecting the multiple components, likely tracing the underlying gas distribution. The Ly$��$ emission originates almost exclusively from the brightest UV component, but we also identify a faint kinematically distinct Ly$��$ emitting region nearby a fainter component. Combined with new near-infrared data, the MUSE data show that the rest-frame Ly$��$ equivalent width (EW) is $\approx100$ ��. This is a factor four higher than the EW measured in low-redshift analogues with carefully matched Ly$��$ profiles (and thus arguably HI column density), but this EW can plausibly be explained by star formation. Alternative scenarios requiring AGN powering are also disfavoured by the narrower and steeper Ly$��$ spectrum and much smaller IR to UV ratio compared to obscured AGN in other Ly$��$ blobs. CR7's Ly$��$ emission, while extremely luminous, resembles the emission in more common LAEs at lower redshifts very well and is likely powered by a young metal poor starburst., Main text 16 pages, 9 figures. Accepted for publication by MNRAS

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., 5 pages, 3 figures. Proceeding of IAU Symposium #352, 2019, "Uncovering early galaxy evolution in the ALMA and JWST era"

The observed properties of the Lyman-$��$ (Ly$��$) emission line are a powerful probe of neutral gas in and around galaxies. We present spatially resolved Ly$��$ spectroscopy with VLT/MUSE targeting VR7, a UV-luminous galaxy at $z=6.532$ with moderate Ly$��$ equivalent width (EW$_0\approx38$ ��). These data are combined with deep resolved [CII]$_{\rm 158 ��m}$ spectroscopy obtained with ALMA and UV imaging from {\it HST} and we also detect UV continuum with MUSE. Ly$��$ emission is clearly detected with S/N $\approx40$ and FWHM of 374 km s$^{-1}$. Ly$��$ 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��, halo} = 3.45^{+1.08}_{-0.87}$ kpc when measured jointly with the UV continuum. The Ly$��$ 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$��$ 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$��$ properties in VR7 resemble the conditions in post-re-ionisation galaxies., 13 pages, 10 figures. Accepted for publication in MNRAS

[ABRIDGED] We present spectroscopic observations of 466 galaxies in and around a superstructure at $z\sim0.84$ targeted by the VIMOS Spectroscopic Survey of a Supercluster in the COSMOS field (VIS$^{3}$COS). We use [OII]$��$3727, H$��$, and $D_n4000$ to trace the recent, mid-, and long-term star formation histories and investigate how stellar mass and the local environment impacts those. By studying trends both in individual and composite galaxy spectra, we find that both stellar mass and environment play a role in the observed galactic properties. We find that the median [OII] equivalent width (|EW$_\mathrm{[OII]}|$) decreases from $27\pm2$ �� to $2.0_{-0.4}^{+0.5}$ �� and $D_n4000$ increases from $1.09\pm0.01$ to $1.56\pm0.03$ with increasing stellar mass (from $\sim10^{9.25}$ to $\sim10^{11.35}\ \mathrm{M_\odot}$). Concerning the dependence on the environment, we find that at fixed stellar mass |EW$_\mathrm{[OII]}|$ is tentatively lower in higher density environments. Regarding $D_n4000$, we find that the increase with stellar mass is sharper in denser environments, hinting that such environments may accelerate galaxy evolution. Moreover, we find larger $D_n4000$ values in denser environments at fixed stellar mass, suggesting that galaxies are on average older and/or more metal-rich in such dense environments. This set of tracers depicts a scenario where the most massive galaxies have, on average, the lowest sSFRs and the oldest stellar populations (age $\gtrsim1$ Gyr, showing a mass-downsizing effect). We also hypothesize that the observed increase in star formation (higher EW$_\mathrm{[OII]|}$, higher sSFR) at intermediate densities may lead to quenching since we find the quenched fraction to increase sharply from the filament to cluster-like regions at similar stellar masses., 20 pages, 13 figures, accepted for publication in A&A

We investigate the clustering and halo properties of $\sim 5000$ Ly$��$-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}��} \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}��} \sim 10^{43 - 43.6}$ erg s$^{-1}$. We find a strong, redshift-independent correlation between halo mass and Ly$��$ luminosity normalized by the characteristic Ly$��$ 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 ��~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., 19 pages, 13 figures, 3 tables, submitted to MNRAS. Comments are welcomed

Using high resolution spectra from the VLT LEGA-C program, we reconstruct the star formation histories (SFHs) of 607 galaxies at redshifts $z = 0.6-1.0$ and stellar masses $\gtrsim10^{10}$M$_{\odot}$ using a custom full spectrum fitting algorithm that incorporates the emcee and FSPS packages. We show that the mass-weighted age of a galaxy correlates strongly with stellar velocity dispersion ($��_*$) and ongoing star-formation (SF) activity, with the stellar content in higher-$��_*$ galaxies having formed earlier and faster. The SFHs of quiescent galaxies are generally consistent with passive evolution since their main SF epoch, but a minority show clear evidence of a rejuvenation event in their recent past. The mean age of stars in galaxies that are star-forming is generally significantly younger, with SF peaking after $z100$Myrs. This indicates that $z>2$ progenitors of $z\sim1$ star-forming galaxies are generally far less massive. Finally, despite considerable variance in the individual SFHs, we show that the current SF activity of massive galaxies ($>$L$_*$) at $z\sim1$ correlates with SF levels at least $3$Gyrs prior: SFHs retain `memory' on a large fraction of the Hubble time. Our results illustrate a novel approach to resolve the formation phase of galaxies, and in identifying their individual evolutionary paths, connects progenitors and descendants across cosmic time. This is uniquely enabled by the high-quality continuum spectroscopy provided by the LEGA-C survey., 18 pages, 9 figures, accepted by ApJ

We present the second data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO 130-night public spectroscopic survey conducted with VIMOS on the Very Large Telescope. We release 1988 spectra with typical continuum S / N ~= 20 /Angstrom of galaxies at 0.6 ~< z ~< 1.0, each observed for ~20 hours and fully reduced with a custom-built pipeline. We also release a catalog with spectroscopic redshifts, emission line fluxes, Lick/IDS indices, and observed stellar and gas velocity dispersions that are spatially integrated quantities including both rotational motions and genuine dispersion. To illustrate the new parameter space in the intermediate redshift regime probed by LEGA-C we explore relationships between dynamical and stellar population properties. The star-forming galaxies typically have observed stellar velocity dispersions of ~150 km/s and strong Hdelta absorption (Hd_A ~ 5 Angstrom), while passive galaxies have higher observed stellar velocity dispersions ~200 km/s and weak Hdelta absortion (Hd_A ~ 0 Angstrom). Strong [O III]5007 / Hbeta ratios tend to occur mostly for galaxies with weak Hd_A or galaxies with higher observed velocity dispersion. Beyond these broad trends, we find a large diversity of possible combinations of rest-frame colors, absorption line strengths and emission line detections, illustrating the utility of spectroscopic measurements to more accurately understand galaxy evolution. By making the spectra and value-added catalogs publicly available we encourage the community to take advantage of this very substantial investment in telescope time provided by ESO., accepted for publication by ApJ; LEGA-C data-release http://www.mpia.de/home/legac/

Drawing from the LEGA-C dataset, we present the spectroscopic view of the stellar population across a large volume- and mass-selected sample of galaxies at large lookback time. We measure the 4000�� break (D$_n$4000) and Balmer absorption line strengths (probed by H$��$) from 1019 high-quality spectra of $z=0.6 - 1.0$ galaxies with $M_\ast = 2 \times 10^{10} M_\odot - 3 \times 10^{11} M_\odot$. Our analysis serves as a first illustration of the power of high-resolution, high-S/N continuum spectroscopy at intermediate redshifts as a qualitatively new tool to constrain galaxy formation models. The observed D$_n$4000-EW(H$��$) distribution of our sample overlaps with the distribution traced by present-day galaxies, but $z\sim 0.8$ galaxies populate that locus in a fundamentally different manner. While old galaxies dominate the present-day population at all stellar masses $> 2\times10^{10} M_\odot$, we see a bimodal D$_n$4000-EW(H$��$) distribution at $z\sim0.8$, implying a bimodal light-weighted age distribution. The light-weighted age depends strongly on stellar mass, with the most massive galaxies $>1\times10^{11}M_\odot$ being almost all older than 2 Gyr. At the same time we estimate that galaxies in this high mass range are only $\sim3$ Gyr younger than their $z\sim0.1$ counterparts, at odd with pure passive evolution given a difference in lookback time of $>5$ Gyr; younger galaxies must grow to $>10^{11}M_\odot$ in the meantime, and/or small amounts of young stars must keep the light-weighted ages young. Star-forming galaxies at $z\sim0.8$ have stronger H$��$ absorption than present-day galaxies with the same D$_n$4000, implying larger short-term variations in star-formation activity., 17 pages, 9 figures, accepted by ApJ

We present spectroscopic follow-up of candidate luminous Ly$��$ 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��} \approx 3\times10^{43}$ erg s$^{-1}$, very high rest-frame equivalent widths of EW$_0\gtrsim 200$ �� and narrow Ly$��$ 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$��$, 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$��$ equivalent widths. We confirm that Ly$��$ line-widths increase with Ly$��$ luminosity at $z=5.7$, while there are indications that Ly$��$ 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$��$ 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$��$ escape fractions $\gtrsim100$ \%, indicating bursty star-formation histories, alternative Ly$��$ production mechanisms, or dust attenuating Ly$��$ emission differently than UV emission. Finally, we present a method to compute $��_{ion}$, the production efficiency of ionising photons, and find that LAEs at $z\approx6-7$ have high values of log$_{10}(��_{ion}$/Hz erg$^{-1}) \approx 25.51\pm0.09$ that may alleviate the need for high Lyman-Continuum escape fractions required for reionisation., 18 pages, 8 figures. Accepted for publication by MNRAS. Added new HST observations of VR7

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 ($, Comment: Accepted for publication to MNRAS. Comments welcome! 17 pages, 12 figures, 8 tables. Paper II can be found here: https://arxiv.org/abs/1704.01124

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}, Comment: Accepted for publication to MNRAS. Comments welcome! 11 pages, 7 figures, 5 tables. Paper I can be found here: https://arxiv.org/abs/1703.10169

We investigate the clustering properties of $\sim 7000$ H$��$+[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$��$+[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$��$+[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$��$+[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., 19 pages, 11 figures, 4 tables, submitted to MNRAS

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��tes field identified with a suite of six narrow-band filters at $\approx 0.4-2.1$ $��$m. These galaxies have been selected on their Ly$��$ (73), {\sc [Oii]} (285), H$��$/{\sc [Oiii]} (387) or H$��$ (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$��$) detections and/or spectroscopy. In this survey paper, we present the observations, selection and catalogs of emitters. We measure number densities of Ly$��$, [O{\sc ii}], H$��$/{\sc [Oiii]} and H$��$ 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$��$ emitters to measure the observed [O{\sc ii}]/H$��$ ratio at $z=1.47$. The observed [O{\sc ii}]/H$��$ 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$��$ and Ly$��$ 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}$., Accepted for publication in MNRAS. Main text 16 pages, 7 Figures. Catalogues of line-emitters available

We study the effects of local environment and stellar mass on galaxy properties using a mass complete sample of quiescent and star-forming systems in the COSMOS field at $z\lesssim$ 3. We show that at $z\lesssim$ 1, the median star-formation rate (SFR) and specific SFR (sSFR) of all galaxies depend on environment, but they become independent of environment at $z\gtrsim$ 1. However, we find that only for \textit{star-forming} galaxies, the median SFR and sSFR are similar in different environments, regardless of redshift and stellar mass. We find that the quiescent fraction depends on environment at $z\lesssim$ 1, and on stellar mass out to $z\sim$ 3. We show that at $z\lesssim$ 1, galaxies become quiescent faster in denser environments and that the overall environmental quenching efficiency increases with cosmic time. Environmental and mass quenching processes depend on each other. At $z\lesssim$ 1, denser environments more efficiently quench galaxies with higher masses (log($M/M_{\odot}$)$\gtrsim$ 10.7), possibly due to a higher merger rate of massive galaxies in denser environments, and that mass quenching is more efficient in denser regions. We show that the overall mass quenching efficiency ($\epsilon_{mass}$) for more massive galaxies (log($M/M_{\odot}$)$\gtrsim$ 10.2) rises with cosmic time until $z\sim$ 1 and flattens out since then. However, for less massive galaxies, the rise in $\epsilon_{mass}$ continues to the present time. Our results suggest that environmental quenching is only relevant at $z\lesssim$ 1, likely a fast process, whereas mass quenching is the dominant mechanism at $z\gtrsim$ 1, with a possible stellar feedback physics., Comment: 17 pages, 10 figures, Accepted for publication in the ApJ

`Direct collapse black holes' (DCBHs) provide possible seeds for supermassive black holes that exist at redshifts as high as z~7. We study Lyman Alpha (Lya) radiative transfer through simplified representations of the DCBH-scenario. We find that gravitational heating of the collapsing cloud gives rise to a Lya cooling luminosity of up to ~ 1e38(M_gas/1e6 Msun)^2 erg/s. The Lya production rate can be significantly larger during the final stages of collapse, but collisional deexcitation efficiently suppresses the emerging Lya flux. Photoionization by a central source boosts the Lya luminosity to L~1e43(M_BH/1e6 M_sun) erg/s during specific evolutionary stages of the cloud, where M_BH denotes the mass of the black hole powering this source. We predict that the width and velocity off-set of the Lya spectral line range from a few tens to few thousands km/s, depending sensitively on the evolutionary state of the cloud. We also compare our predictions to observations of CR7 (Sobral et al. 2015), a luminous Lya emitter at z~7, which is potentially associated with a DCBH. If CR7 is powered by a black hole, then its Lya flux alone requires that M_BH> 1e7 M_sun, which exceeds the mass of DCBHs when they first form. The observed width of the Lya spectrum favors the presence of only a low column density of hydrogen, log [N_HI/cm^-2]~19-20. The shape of the Lya spectrum indicates that this gas is outflowing. These requirements imply that if CR7 harbors a DCBH, then the physical conditions that enabled its formation have been mostly erased, which is in agreement with theoretical expectations., 16 pages, 7 figures, submitted to ApJ following revision based on referee report

We use new near-infrared spectroscopic observations to investigate the nature and evolution of the most luminous H��(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��}>L_{H��}^*$, 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��}/L_{H��}^*(z)$. While $L_{H��}80%), the most luminous Ha emitters ($L_{H��}>10L_{H��}^*(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��}>L^*_{H��}(z)$) contribute the same fractional amount of ~15% to the total star-formation rate density, at least up to z=2.23., 16 pages, 10 figures. Accepted for publication in MNRAS

We investigate the properties of $\sim7000$ narrow-band selected galaxies with strong H$��+$[OIII] and [OII] nebular emission lines from the High-$z$ Emission Line Survey (HiZELS) between $z \sim 0.8 - 5.0$. Our sample covers a wide range in stellar mass ($M\sim10^{7.5 - 12.0}$ M$_\odot$), rest-frame equivalent widths (EW$_\mathrm{rest}\sim 10 - 10^5$ ��), and line luminosities ($L\sim10^{40.5 - 43.2}$ erg s$^{-1}$). We measure the H$��+$[OIII] and [OII]-selected stellar mass functions out to $z \sim 3.5$ and find that both $M_\star$ and $��_\star$ increases with cosmic time, which may be due to the [OIII] selection including an increasing fraction of AGN at lower redshifts. The [OII]-selected stellar mass functions show a constant $M_\star\sim10^{11.6}$ M$_\odot$ and a strong, increasing evolution with cosmic time in $��_\star$ in line with H$��$ studies. We also investigate the EW$_\mathrm{rest}$ evolution as a function of redshift with a fixed mass range (10$^{9.5 - 10.0}$ M$_\odot$) and find an increasing trend best represented by $(1+z)^{3.81\pm0.14}$ and $(1+z)^{2.72\pm0.19}$ up to $z\sim 2$ and $\sim 3$ for H$��$+[OIII] and [OII] emitters, respectively. This is the first time that the EW evolution has been directly measured for H$��+$[OIII] and [OII] emitters up to these redshifts. There is evidence for a slower evolution for $z > 2$ in the H$��+$[OIII] EW and a decreasing trend for $z > 3$ in the [OII] EW evolution, which would imply low [OII] EW at the highest redshifts and higher [OIII]/[OII] line ratios. This suggests that the ionization parameter at higher redshift may be significantly higher than the local Universe. Our results set the stage for future near-IR space-based spectroscopic surveys to test our extrapolated predictions and also produce $z > 5$ measurements to constrain the high-$z$ end of the EW$_\mathrm{rest}$ and [OIII]/[OII] evolution., 20 pages, 10 figures, 6 tables, submitted to MNRAS

We present the first results from the largest H$��$ survey of star formation and AGN activity in galaxy clusters. Using 9 different narrow band filters, we select $>3000$ H$��$ emitters within $19$ clusters and their larger scale environment over a total volume of $1.3\times10^5$ Mpc$^3$. The sample includes both relaxed and merging clusters, covering the $0.15-0.31$ redshift range and spanning from $5\times10^{14}$ $M_{\odot}$ to $30\times10^{14}$ $M_{\odot}$. We find that the H$��$ luminosity function (LF) for merging clusters has a higher characteristic density $��^*$ compared to relaxed clusters. $��^*$ drops from cluster core to cluster outskirts for both merging and relaxed clusters, with the merging cluster values $\sim0.3$ dex higher at each projected radius. The characteristic luminosity $L^*$ drops over the $0.5-2.0$ Mpc distance from the cluster centre for merging clusters and increases for relaxed objects. Among disturbed objects, clusters hosting large-scale shock waves (traced by radio relics) are overdense in H$��$ emitters compared to those with turbulence in their intra-cluster medium (traced by radio haloes). We speculate that the increase in star formation activity in disturbed, young, massive galaxy clusters can be triggered by interactions between gas-rich galaxies, shocks and/or the intra-cluster medium, as well as accretion of filaments and galaxy groups. Our results indicate that disturbed clusters represent vastly different environments for galaxy evolution compared to relaxed clusters or average field environments., Accepted for publication in MNRAS. Comments welcome. 20 pages, 10 figures

We present the first results from our CAlibrating LYMan-$��$ with H$��$ (CALYMHA) pilot survey at the Isaac Newton Telescope. We measure Ly$��$ emission for 488 H$��$ selected galaxies at $z=2.23$ from HiZELS in the COSMOS and UDS fields with a specially designed narrow-band filter ($��_c$ = 3918 ��, $����$= 52 ��). We find 17 dual H$��$-Ly$��$ emitters ($f_{\rm Ly��} >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$��$ 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$��$ 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$��$ is significantly more extended than H$��$ and the UV. f$_{\rm esc}$ continues to increase up to at least 20 kpc (3$��$, 40 kpc [2$��$]) for typical SFGs and thus the aperture is the most important predictor of f$_{\rm esc}$., Accepted for publication by MNRAS

We carried out the largest ($>3.5\times10^5$ Mpc$^3$, 26 deg$^2$) H$��$ narrow band survey to date at $z\sim0.2$ in the SA22, W2 and XMMLSS extragalactic fields. Our survey covers a large enough volume to overcome cosmic variance and to sample bright and rare H$��$ emitters up to an observed luminosity of $\sim10^{42.4}$ erg s$^{-1}$, equivalent to $\sim11 M_\odot$ yr$^{-1}$. Using our sample of $220$ sources brighter than $>10^{41.4}$ erg s$^{-1}$ ($>1 M_\odot$ yr$^{-1}$), we derive H$��$ luminosity functions, which are well described by a Schechter function with $��^* = 10^{-2.85\pm0.03}$ Mpc$^{-3}$ and $L^*_{H��} = 10^{41.71\pm0.02}$ erg s$^{-1}$ (with a fixed faint end slope $��=-1.35$). We find that surveys probing smaller volumes ($\sim3\times10^4$ Mpc$^3$) are heavily affected by cosmic variance, which can lead to errors of over $100$ per cent in the characteristic density and luminosity of the H$��$ luminosity function. We derive a star formation rate density of $��_\mathrm{SFRD} = 0.0094\pm0.0008$ $M_\odot$ yr$^{-1}$, in agreement with the redshift-dependent H$��$ parametrisation from Sobral et al. (2013). The two-point correlation function is described by a single power law $��(��) = (0.159\pm0.012) ��^{(-0.75\pm0.05)}$, corresponding to a clustering length of $r_0 = 3.3\pm0.8$ Mpc/h. We find that the most luminous H$��$ emitters at $z\sim0.2$ are more strongly clustered than the relatively fainter ones. The $L^*_{H��}$ H$��$ emitters at $z\sim0.2$ in our sample reside in $\sim10^{12.5-13.5}$ $M_\odot$ dark matter haloes. This implies that the most star forming galaxies always reside in relatively massive haloes or group-like environments and that the typical host halo mass of star-forming galaxies is independent of redshift if scaled by $L_\mathrm{H��}/L^*_{H��}(z)$, as proposed by Sobral et al. (2010)., Accepted for publication in MNRAS; 18 pages, 19 figures, 6 tables

We investigate the evolution of the H$��$+[OIII] and [OII] luminosity functions from $z \sim 0.8$ to $\sim5$ in four redshift slices per emission line using data from the High-{\it z} Emission Line Survey (HiZELS). This is the first time that the H$��$+[OIII] and [OII] luminosity functions have been studied at these redshifts in a self-consistent analysis. This is also the largest sample of [OII] and H$��$+[OIII] emitters (3475 and 3298 emitters, respectively) in this redshift range, with large co-moving volumes $\sim 1 \times 10^6$ Mpc$^{-3}$ in two independent volumes (COSMOS and UDS), greatly reducing the effects of cosmic variance. The emitters were selected by a combination of photometric redshift and color-color selections, as well as spectroscopic follow-up, including recent spectroscopic observations using DEIMOS and MOSFIRE on the Keck Telescopes and FMOS on Subaru. We find a strong increase in $L_\star$ and a decrease in $��_\star$ for both H$��$+[OIII] and [OII] emitters. We derive the [OII] star-formation history of the Universe since $z\sim5$ and find that the cosmic SFRD rises from $z \sim 5$ to $\sim 3$ and then drops towards $z \sim 0$. We also find that our star-formation history is able to reproduce the evolution of the stellar mass density up to $z\sim 5$ based only on a single tracer of star-formation. When comparing the H$��$+[OIII] SFRDs to the [OII] and H$��$ SFRD measurements in the literature, we find that there is a remarkable agreement, suggesting that the H$��$+[OIII] sample is dominated by star-forming galaxies at high-$z$ rather than AGNs., 24 pages, 14 figures, 7 Tables, MNRAS accepted

Faint Lyman-$��$ (Ly$��$) emitters become increasingly rarer towards the re-ionisation epoch (z~6-7). However, observations from a very large (~5deg$^2$) Ly$��$ 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$��$ 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$��$ with a line width of $386\pm30$ km/s (FWHM) and with high EW$_0$ (>200 ��), but it is undetected in the continuum. CR7, with an observed Ly$��$ luminosity of $10^{43.93\pm0.05}$erg/s is the most luminous Ly$��$ emitter ever found at z>6. CR7 reveals a narrow Ly$��$ line with $266\pm15$ km/s FWHM, being detected in the NIR (rest-frame UV, with $��=-2.3\pm0.1$) with an excess in $J$, and also strongly detected in IRAC/Spitzer. We detect a narrow HeII1640$��$ emission line ($6��$) which explains the excess seen in the $J$ band photometry (EW$_0$~80 ��). 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$��$ 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., 15 pages, 8 figures. Accepted for publication in ApJ

Using wide field narrow-band surveys, we provide a new measurement of the $z=6.6$ Lyman-$��$ 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$��$ 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$��$ 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$��$ 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., 20 pages, main results shown in Fig. 6 and Fig. 7, accepted by MNRAS

We investigate resolved kpc-scale stellar and nebular dust distribution in eight star-forming galaxies at z~0.4 in the GOODS fields. Constructing the observed Spectral Energy Distributions (SEDs) per pixel, based on seven bands photometric data from HST/ACS and WFC3, we performed pixel-by-pixel SED fits to population synthesis models and estimated small-scale distribution of stellar dust extinction. We use Halpha / Hbeta nebular emission line ratios from Keck/DEIMOS high resolution spectra at each spatial resolution element to measure the amount of attenuation faced by ionized gas at different radii from the center of galaxies. We find a good agreement between the integrated and median of resolved color excess measurements in our galaxies. The ratio of integrated nebular to stellar dust extinction is always greater than unity, but does not show any trend with stellar mass or star formation rate. We find that inclination plays an important role in the variation of the nebular to stellar excess ratio. The stellar color excess profiles are found to have higher values at the center compared to outer parts of the disk. However, for lower mass galaxies, a similar trend is not found for the nebular color excess. We find that the nebular color excess increases with stellar mass surface density. This explains the absence of radial trend in the nebular color excess in lower mass galaxies which lack a large radial variation of stellar mass surface density. Using standard conversions of star formation rate surface density to gas mass surface density, and the relation between dust mass surface density and color excess, we find no significant variation in the dust to gas ratio in regions with high gas mass surface densities, over the scales probed in this study., Comment: 10 pages, 8 figures, accepted for publication in the ApJ

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 $��$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$��$ emitters at z=0.8, 1343 [OIII]+H$��$ 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$��$ (spectroscopically confirmed) over-density of H$��$ emitters at z=0.81. We derive H$��$, [OIII]+H$��$ 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 5x10$^{5}$ Mpc$^{3}$; errors on $L^*$ and $��^*$ 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., 22 pages, 15 figures, accepted for publication in MNRAS. Full catalogue of line emitters available in FITS format with final MNRAS published paper

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-$��$ (Ly$��$) 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$��$ luminosity limit of $6.3\times10^{43}$ erg s$^{-1}$. We look for Ly$��$ 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$��$ 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$��$ emitters. We put a strong constraint on the Ly$��$ 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., 15 pages, 9 figures, accepted for publication in MNRAS

The most extreme cluster mergers can lead to massive cluster-wide travelling shock waves. The CIZA J2242.8+5301 ('sausage') and 1RXS J0603.3+4213 (`toothbrush') clusters ($z\sim0.2$) host enormous radio-emitting shocks with simple geometry. We investigate the role of mergers and shocks in shaping the H$��$ luminosity function, using custom-made narrow-band filters matching the cluster redshifts mounted on the INT. We surveyed $\sim0.28$ deg$^2$ for each cluster and found $181$ line emitters in the `sausage' (volume of $3.371\times10^3$ Mpc$^3$ for H$��$ at $z=0.1945$) and $141$ in the `toothbrush' ($4.546\times10^3$ Mpc$^3$ for H$��$ at $z=0.225$), out of which $49$ (`sausage') and $30$ (`toothbrush') are expected to be H$��$. We build luminosity functions for the field-of-view down to an average limiting star formation rate of $0.14$ M$_{\odot}$ yr$^{-1}$, find good agreement with field luminosity functions at $z=0.2$, but significant differences between the shapes of the luminosity functions for the two clusters. We discover extended, tens-of-kpc-wide H$��$ haloes in galaxies neighbouring relics, which were possibly disrupted by the passage of the shock wave. By comparing the `sausage' cluster with blank fields and other clusters, we also uncover an order of magnitude boost (at $9��$ level) in the normalisation $��^*$ of the luminosity function in the relic areas. Our results suggest that cluster mergers may play an important role in the evolution of cluster galaxies through shock-induced star formation., Accepted for publication in MNRAS, 14 pages, 9 figures

We present adaptive optics assisted, spatially resolved spectroscopy of a sample of nine H-alpha-selected galaxies at z=0.84--2.23 drawn from the HiZELS narrow-band survey. These galaxies have star-formation rates of 1-27Mo/yr and are therefore representative of the typical high-redshift star-forming population. Our ~kpc-scale resolution observations show that approximately half of the sample have dynamics suggesting that the ionised gas is in large, rotating disks. We model their velocity fields to infer the inclination-corrected, asymptotic rotational velocities. We use the absolute B-band magnitudes and stellar masses to investigate the evolution of the B-band and stellar mass Tully-Fisher relationships. By combining our sample with a number of similar measurements from the literature, we show that, at fixed circular velocity, the stellar mass of star-forming galaxies has increased by a factor 2.5 between z=2 and z=0, whilst the rest-frame B-band luminosity has decreased by a factor ~6 over the same period. Together, these demonstrate a change in mass-to-light ratio in the B-band of Delta(M/L_B)/(M/L_B)_(z=0) \sim 3.5 between z=1.5 and z=0, with most of the evolution occurring below z=1. We also use the spatial variation of [NII]/Halpha to show that the metallicity of the ionised gas in these galaxies declines monotonically with galacto-centric radius, with an average Delta(log O/H)/DeltaR=-0.027+/-0.005dex/kpc. This gradient is consistent with predictions for high-redshift disk galaxies from cosmologically based hydrodynamic simulations., Comment: Accepted for publication in MNRAS. 18 pages, 10 figures, 3 tables

We investigate the relationship between H�� and [OII](��3727) emission in faint star-forming galaxies at z=1.47 with dust uncorrected star formation rates (SFRs) down to 1.4 Msun/yr, using data in two narrow-bands from WFCAM/UKIRT and Suprime-Cam/Subaru. A stacking analysis allows us to investigate H�� emission flux from bright [OII] emitters as well as faint ones for which H�� is not individually detected, and to compare them with a large sample of local galaxies. We find that there is a clear, positive correlation between the average H�� and [OII] luminosities for [OII] emitters at z=1.47, with its slope being consistent with the local relation. [OII] emitters at z=1.47 have lower mean observed ratios of H��/[OII] suggesting a small but systematic offset (at 2.8�� significance) towards lower values of dust attenuation, A(H��)~0.35, than local galaxies. This confirms that [OII] selection tends to pick up galaxies which are significantly less dusty on average than H�� selected ones, with the difference being higher at z=1.47 than at z=0. The discrepancy of the observed line ratios between [OII] emitters at z=1.47 and the local galaxies may in part be due to the samples having different metallicities. However, we demonstrate that metallicity is unlikely to be the main cause. Therefore, it is important to take into account that the relations for the dust correction which are derived using H�� emitter samples, and frequently used in many studies of high-z galaxies, may overestimate the intrinsic SFRs of [OII]-selected galaxies, and that surveys of [OII] emission galaxies are likely to miss dusty populations., 10 pages, 7 figures, 1 table, accepted for publication in MNRAS

This paper presents the results from the first wide and deep dual narrow-band survey to select H-alpha (Ha) and [OII] line emitters at z=1.47+-0.02 (using matched narrow-band filters in the H and z' bands), exploiting synergies between the UKIRT and Subaru telescopes. The Ha survey at z=1.47 reaches a flux limit of ~7x10^-17 erg/s/cm^2 and detects ~200 Ha emitters over 0.7deg^2, while the much deeper [OII] survey reaches an effective flux of ~7x10^-18 erg/s/cm^2, detecting ~1400 z=1.47 [OII] emitters in a matched co-moving volume of ~2.5x10^5 Mpc^3. The combined survey results in the identification of 190 simultaneous Ha and [OII] emitters at z=1.47. Ha and [OII] luminosity functions are derived and both are shown to evolve significantly from z~0 in a consistent way. The star formation rate density of the Universe at z=1.47 is evaluated, with the Ha analysis yielding 0.16+-0.05 M_sun/yr/Mpc^3 and the [OII] analysis 0.17+-0.04 M_sun/yr/Mpc^3. The measurements are combined with other studies, providing a self-consistent measurement of the star formation history of the Universe over the last ~11Gyrs. By using a large comparison sample at z~0.1 (from the SDSS), [OII]/Ha line ratios are calibrated as probes of dust-extinction. Ha emitters at z~1.47 show on average 1 mag of extinction at Ha, similar to the SDSS sources at z~0. Although we find that dust extinction correlates with SFR, the relation evolves by about ~0.5 mag from z~1.5 to z~0, with z~0 relations over-predicting the dust extinction corrections at high-z by that amount. Stellar mass is found to be a much more fundamental extinction predictor, with the relation between mass and extinction being valid at both z~0 and z~1.5. Dust extinction corrections as a function of optical colours are also derived, offering simpler mechanisms for estimating extinction in moderately star-forming systems over the last ~9Gyrs [Abridged]., Accepted for publication in MNRAS. Revised version following the reviewer's report. 21 pages, 23 figures