Search

The code FADO is the first publicly available population spectral synthesis tool that treats the contribution from ionised gas to the observed emission self-consistently. We study the impact of the nebular contribution on the determination of the star formation rate (SFR), stellar mass, and consequent effect on the star-forming main sequence (SFMS) at low redshift. We applied FADO to the spectral database of the SDSS to derive the physical properties of galaxies. As a comparison, we used the data in the MPA-JHU catalogue, which contains the properties of SDSS galaxies derived without the nebular contribution. We selected a sample of SF galaxies with H$\alpha$ and H$\beta$ flux measurements, and we corrected the fluxes for the nebular extinction through the Balmer decrement. We then calculated the H$\alpha$ luminosity to estimate the SFR. Then, by combining the stellar mass and SFR estimates from FADO and MPA-JHU, the SFMS was obtained. The H$\alpha$ flux estimates are similar between FADO and MPA-JHU. Because the H$\alpha$ flux was used as tracer of the SFR, FADO and MPA-JHU agree in their SFR. The stellar mass estimates are slightly higher for FADO than for MPA-JHU on average. However, considering the uncertainties, the differences are negligible. With similar SFR and stellar mass estimates, the derived SFMS is also similar between FADO and MPA-JHU. Our results show that for SDSS normal SF galaxies, the additional modelling of the nebular contribution does not affect the retrieved fluxes and consequentially also does not influence SFR estimators based on the extinction-corrected H$\alpha$ luminosity. For the stellar masses, the results point to the same conclusion. These results are a consequence of the fact that the vast majority of normal SF galaxies in the SDSS have a low nebular contribution., Comment: 13 pages, 10 figures and 3 tables. Forthcoming article in Astronomy and Astrophysics

Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions into the main luminosity contributors (e.g. stellar populations or ionised gas). However, the impact nebular emission has on the inferred properties of star-forming (SF) galaxies has been largely overlooked over the years. The objective of this work is to estimate the relations between stellar properties of SF galaxies from SDSS DR7 by simultaneously fitting the stellar and nebular continua with FADO and comparing them to the results derived using STARLIGHT, a representative of purely stellar population synthesis codes. Differences between codes regarding average mass, mean age and mean metallicity values can go as high as $\sim$0.06 dex for the overall population of galaxies and $\sim$0.12 dex for SF galaxies (galaxies with EW(H$\alpha$)>3 \AA), with the most prominent difference between both codes in the light-weighted mean stellar age. A closer look into the average light- and mass-weighted star formation histories of intensively SF galaxies (EW(H$\alpha$)>75 \AA) suggests that STARLIGHT is underestimating the average light-weighted age of intensively SF galaxies by up to $\sim$0.17 dex and overestimating the light-weighted metallicity by up to $\sim$0.13 dex compared to FADO (or vice versa). The comparison between the average stellar properties of passive, SF and intensively SF galaxy samples also reveals that differences between codes increase with increasing EW(H$\alpha$) and decreasing total stellar mass. This work finds indirect evidence that a purely stellar population synthesis approach negatively impacts the inferred stellar properties of galaxies with relatively high star formation rates. In turn, this can bias interpretations of fundamental relations such as the mass-age or mass-metallicity., Comment: 27 pages, 25 figures. Accepted for publication in Astronomy & Astrophysics

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

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$\alpha$ 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$\alpha$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Ly\alpha}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$\alpha$ 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$\alpha$ 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\%$, $\xi_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), a highly ionizing minority of galaxies with $M_{\rm UV}<-17$ accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a "disco" ionizing background. We conclude proposing tests to further develop our suggested Ly$\alpha$-anchored formalism., Comment: 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}}<5$%) from a representative sample of $>0.2 L^{*}$ LyA emitters (LAEs). With median stacked spectra of these subsets covering 1000-8000 {\AA} (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 {\AA}) 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({\xi_{\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 $\xi_{\rm{ion}}$ and prolific $f_{\rm{esc}}$ occurring in sync. (ABRIDGED), Comment: 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!

The field of galaxy evolution will make a great leap forward in the next decade as a consequence of the huge effort by the scientific community in multi-object spectroscopic facilities. To maximise the impact of such incoming data, the analysis methods must also step up, extracting reliable information from the available spectra. In this paper, we aim to investigate the limits and the reliability of different spectral synthesis methods in the estimation of the mean stellar age and metallicity. The main question this work aims to address is which signal-to-noise ratios (S/N) are needed to reliably determine the mean stellar age and metallicity from a galaxy spectrum and how this depends on the tool used to model the spectra. To address this question we built a set of realistic simulated spectra containing stellar and nebular emission, reproducing the evolution of a galaxy in two limiting cases: a constant star formation rate and an exponentially declining star formation. We degraded the synthetic spectra built from these two star formation histories (SFHs) to different S/N and analysed them with three widely used spectral synthesis codes, namely FADO, STECKMAP, and STARLIGHT. For S/N < 5 all three tools show a large diversity in the results. The FADO and STARLIGHT tools find median differences in the light-weighted mean stellar age of ~0.1 dex, while STECKMAP shows a higher value of ~0.2 dex. Detailed investigations of the best-fit spectrum for galaxies with overestimated mass-weighted quantities point towards the inability of purely stellar models to fit the observed spectra around the Balmer jump. Our results imply that when a galaxy enters a phase of high specific star formation rate the neglect of the nebular continuum emission in the fitting process has a strong impact on the estimation of its SFH when purely stellar fitting codes are used, even in presence of high S/N spectra., Comment: New version matching the published one. 19 pages, 16 figures, accepted for publication in A&A

We measure the evolution of the rest-frame UV luminosity function (LF) and the stellar mass function (SMF) of Lyman-alpha (Lya) emitters (LAEs) from z~2 to z~6 by exploring ~4000 LAEs from the SC4K sample. We find a correlation between Lya luminosity (LLya) and rest-frame UV (M_UV), with best-fit M_UV=-1.6+-0.2 log10(LLya/erg/s)+47+-12 and a shallower relation between LLya and stellar mass (Mstar), with best-fit log10( Mstar/Msun)=0.9+-0.1 log10(LLya/erg/s)-28+-4.0. An increasing LLya cut predominantly lowers the number density of faint M_UV and low Mstar LAEs. We estimate a proxy for the full UV LFs and SMFs of LAEs with simple assumptions of the faint end slope. For the UV LF, we find a brightening of the characteristic UV luminosity (M_UV*) with increasing redshift and a decrease of the characteristic number density (Phi*). For the SMF, we measure a characteristic stellar mass (Mstar*/Msun) increase with increasing redshift, and a Phi* decline. However, if we apply a uniform luminosity cut of log10 (LLya/erg/s) >= 43.0, we find much milder to no evolution in the UV and SMF of LAEs. The UV luminosity density (rho_UV) of the full sample of LAEs shows moderate evolution and the stellar mass density (rho_M) decreases, with both being always lower than the total rho_UV and rho_M of more typical galaxies but slowly approaching them with increasing redshift. Overall, our results indicate that both rho_UV and rho_M of LAEs slowly approach the measurements of continuum-selected galaxies at z>6, which suggests a key role of LAEs in the epoch of reionisation., Comment: Accepted for publication in MNRAS. 20 pages, 13 figures, 4 tables

We present the first results from the X-SHOOTER Lyman-$\alpha$ survey at $z=2$ (XLS-$z2$). XLS-$z2$ is a deep spectroscopic survey of 35 Lyman-$\alpha$ emitters (LAEs) utilising $\approx90$ hours of exposure time with VLT/X-SHOOTER and covers rest-frame Ly$\alpha$ to H$\alpha$ 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$\alpha$ line. Typically 30 % of the Ly$\alpha$ photons escape, of which one quarter on the blue side of the systemic velocity. A fraction of Ly$\alpha$ 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}$., Comment: Main text 26 pages. Accepted for publication in MNRAS

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.

Taking advantage of HST CANDELS data, we analyze the lowest redshift (z<0.5) massive galaxies in order to disentangle their structural constituents and study possible faint non-axis-symmetric features. Due to the excellent HST spatial resolution for intermediate-z objects, they are hard to model by purely automatic parametric fitting algorithms. We performed careful single and double S\'ersic fits to their galaxy surface brightness profiles. We also compare the model color profiles with the observed ones and also derive multi-component global effective radii attempting to obtain a better interpretation of the mass-size relation. Additionally, we test the robustness of our measured structural parameters via simulations. We find that the S\'ersic index does not offer a good proxy for the visual morphological type for our sample of massive galaxies. Our derived multi-component effective radii give a better description of the size of our sample galaxies than those inferred from single S\'ersic models with GALFIT. Our galaxy population lays on the scatter of the local mass-size relation, indicating that these massive galaxies do not experience a significant growth in size since z~0.5. Interestingly the few outliers are late-type galaxies, indicating that spheroids must reach the local mass-size relation earlier. For most of our sample galaxies, both single and multi-component S\'ersic models with GALFIT show substantial systematic deviations from the observed SBPs in the outskirts. These residuals may be partly due to several factors, namely a non-optimal data reduction for low surface brightness features, the existence of prominent stellar haloes for massive galaxies and could also arise from conceptual shortcomings of parametric 2D image decomposition tools. They consequently propagate into galaxy color profiles.

[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]$\lambda$3727, H$\delta$, 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$ {\AA} to $2.0_{-0.4}^{+0.5}$ {\AA} 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., Comment: 20 pages, 13 figures, accepted for publication in A&A

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., Comment: Accepted for publication in MNRAS. 23 pages, 11 figures, 3 tables + appendices. The full SC4K catalogue of LAEs with PSF photometry, stellar masses, SFRs and other properties for individual LAEs, is available at https://goo.gl/q9yfKo

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.23) 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

We study the impact of local density and stellar mass on the morphology of ~500 quiescent (Q) and SFGs from the VIS3COS survey. We perform B/D decomposition of the SBPs and find ~41% of 10^10 Msun galaxies to be best fitted with 2 components. We complement our analysis with non-parametric measurements and qualitative visual classifications. We find that galaxy morphology depends on stellar mass and environment for our sample as a whole. We only find an impact of the environment on galaxy sizes for galaxies more massive than 10^11 Msun. We find higher n and B/T in high-density regions when compared to low-density counterparts at similar stellar masses. We also find that higher stellar mass galaxies have steeper light profiles compared to the lower ones. Using visual classifications, we find a morphology-density relation at z~0.84 for galaxies more massive than 10^10 Msun, with elliptical galaxies being dominant at high-density regions and disks more common in low-density regions. However, when splitting the sample into colour-colour selected SF and Q sub-populations, there are no statistically significant differences between low- and high-density regions. We find that Q galaxies are smaller, have higher n, and higher B/T when compared to SF counterparts. We confirm these trends with non-parametric quantities, finding Q galaxies to be smoother and having most of their light over smaller areas than SFGs. Overall, we find that the dependence of morphology on stellar mass is stronger than with local density and these relations are strongly correlated with the quenching fraction. The change in average morphology corresponds to a change in the relative fractions of blue disk-like galaxies and red elliptical galaxies with stellar mass and environment. We hypothesize that the processes responsible for the quenching of SF must also affect the galaxy morphology on similar timescales., Comment: 25 pages, 14 figures, change in the series number, to publish in Astronomy & Astrophysics

The galaxy catalogs generated from low-resolution emission line surveys often contain both foreground and background interlopers due to line misidentification, which can bias the cosmological parameter estimation. In this paper, we present a method for correcting the interloper bias by using the joint-analysis of auto- and cross-power spectra of the main and the interloper samples. In particular, we can measure the interloper fractions from the cross-correlation between the interlopers and survey galaxies, because the true cross-correlation must be negligibly small. The estimated interloper fractions, in turn, remove the interloper bias in the cosmological parameter estimation. For example, in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift ($z<0.5$) [O II] $\lambda3727${\AA} emitters contaminate high-redshift ($1.9

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., Comment: 19 pages, 13 figures, 3 tables, submitted to MNRAS. Comments are welcomed

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., Comment: 14 pages, 10 figures. Accepted for publication by A&A. Data are fully public

We present the VIMOS Spectroscopic Survey of a Supercluster in the COSMOS field (VIS3COS) at z~0.84. We use VIMOS high-resolution spectra (GG475 filter) to spectroscopically select 490 galaxies in and around the super-structure and an additional 481 galaxies in the line of sight. We present the redshift distribution, the catalog to be made public and the first results on the properties of individual galaxies and stacked spectra (3500{\AA}<\lambda<4200{\AA} rest-frame). We probe a wide range of densities and environments (from low-density field to clusters and rich groups). We find a decrease in the median star formation rate from low to high density environments in all bins of stellar mass and a sharp rise of the quenched fraction (from ~10% to ~40-60%) of intermediate stellar mass galaxies $(10<\log_{10}\left(M_\star/\mathrm{M_\odot}\right)<10.75)$ from filaments to clusters. The quenched fraction for massive galaxies shows little dependence on the environment being constant at ~30-40%. We find a break in the strength of the [OII] emission, with nearly constant line equivalent widths at lower densities (~-11{\AA}) and then a drop to ~-2.5{\AA} towards higher densities. The break in the [OII] line strength happens at similar densities ($\log_{10}(1+\delta)\sim0.0-0.5$) as the observed rise in the quenched fraction. Our results may provide further clues regarding the different environmental processes affecting galaxies with different stellar masses and highlight the advantages of a single dataset in the COSMOS field probing a wide range of stellar masses and environments. We hypothesize that quenching mechanisms are enhanced in high-density regions., Comment: 16 pages, 14 figures, accepted for publication on A&A

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., Comment: Accepted for publication in MNRAS. 1D spectra available with the final published version on MNRAS (Supplementary data)

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.53.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., Comment: Accepted for publication in MNRAS. 24 pages + Appendices with SC4K catalogue and LFs. SC4K catalogue of LAEs, SC4K LFs and the S-SC4K compilation are available (https://goo.gl/q9yfKo) prior to publication

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 255 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., Comment: Accepted for publication on MNRAS

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., Comment: 18 pages, 8 figures. Accepted for publication by MNRAS. Added new HST observations of VR7

We present new results on the evolution of rest-frame blue/UV sizes and Sersic indices of H$\alpha$-selected star-forming galaxies over the last 11 Gyrs. We investigate how the perceived evolution can be affected by a range of biases and systematics such as cosmological dimming and resolution effects. We use GALFIT and an artificial redshifting technique, which includes the luminosity evolution of H$\alpha$-selected galaxies, to quantify the change on the measured structural parameters with redshift. We find typical sizes of 2 to 3 kpc and Sersic indices of n~1.2, close to pure exponential disks all the way from z=2.23 to z=0.4. At z=0 we find typical sizes of 4-5 kpc. Our results show that, when using GALFIT, cosmological dimming has a negligible impact on the derived effective radius for galaxies with <10 kpc, but we find a ~20% bias on the estimate of the median Sersic indices, rendering galaxies more disk-like. Star-forming galaxies have grown on average by a factor of 2-3 in the last 11 Gyrs with $r_e\propto(1+z)^{-0.75}$. By exploring the evolution of the stellar mass-size relation we find evidence for a stronger size evolution of the most massive star-forming galaxies since z~2, as they grow faster towards z~0 when compared to the lower stellar mass counterparts. As we are tracing the rest-frame blue/UV, we are likely witnessing the growth of disks where star formation is ongoing in galaxies while their profiles remain close to exponential disks, n<1.5, across the same period., Comment: 20 pages, 18 figures, Accepted for publication in MNRAS

We present the first results from the largest H$\alpha$ survey of star formation and AGN activity in galaxy clusters. Using 9 different narrow band filters, we select $>3000$ H$\alpha$ 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$\alpha$ luminosity function (LF) for merging clusters has a higher characteristic density $\phi^*$ compared to relaxed clusters. $\phi^*$ 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$\alpha$ 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., Comment: Accepted for publication in MNRAS. Comments welcome. 20 pages, 10 figures

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., Comment: 19 pages, 11 figures. Accepted for publication in MNRAS. Catalogue of line emitters available with MNRAS refereed paper

Cluster star-forming galaxies are found to have an excess of Far-Infrared emission relative to H-alpha (Ha), when compared to those in the field, which could be caused by intense AGN activity, dust and/or declining star formation histories. Here we present spectroscopic observations of Ha emitters in the Cl 0939+4713 (Abell 851) super-cluster at z=0.41, using AF2+WYFFOS on the WHT. We measure [OII], Hbeta (Hb), [OIII], Ha and [NII] for a sample of 119 Ha emitters in and around the cluster. We find that 17+-5% of the Ha emitters are AGN, irrespective of environment. For star-forming galaxies, we obtain Balmer decrements, metallicities and ionisation parameters with different methods, individually and by stacking. We find a strong mass-metallicity relation at all environments, with no significant dependence on environment. The ionisation parameter declines with increasing stellar mass for low-mass galaxies. Ha emitters residing in intermediate environments show the highest ionisation parameters (along with high [OIII]/Ha and high [OIII]/[OII] line ratios, typically twice as large as in the highest and lowest densities), which decline with increasing environmental density. Dust extinction (A$_{H\alpha}$) correlates strongly with stellar mass, but also with environmental density. Star-forming galaxies in the densest environments are found to be significantly dustier (A$_{H\alpha}$~1.5-1.6) than those residing in the lowest density environments (A$_{H\alpha}$~0.6), deviating significantly from what would be predicted given their stellar masses., Comment: Accepted for publication in MNRAS. 13 pages, 11 figures

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}$., Comment: Accepted for publication by MNRAS