@inproceedings{11586,
  abstract     = {Distant luminous Lyman-α 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-α 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.},
  author       = {Matthee, Jorryt J and Sobral, David},
  booktitle    = {Proceedings of the International Astronomical Union},
  issn         = {1743-9221},
  keywords     = {Astronomy and Astrophysics, Space and Planetary Science, galaxies: formation, galaxies: evolution, galaxies: high-redshift},
  number       = {S352},
  pages        = {21--25},
  publisher    = {Cambridge University Press},
  title        = {{Unveiling the most luminous Lyman-α emitters in the epoch of reionisation}},
  doi          = {10.1017/s1743921319009451},
  volume       = {15},
  year         = {2020},
}

@article{11499,
  abstract     = {Deep optical spectroscopic surveys of galaxies provide a unique opportunity to investigate rest-frame ultra-violet (UV) emission line properties of galaxies at z ∼ 2 − 4.5. Here we combine VLT/MUSE Guaranteed Time Observations of the Hubble Deep Field South, Ultra Deep Field, COSMOS, and several quasar fields with other publicly available data from VLT/VIMOS and VLT/FORS2 to construct a catalogue of He II λ1640 emitters at z ≳ 2. The deepest areas of our MUSE pointings reach a 3σ line flux limit of 3.1 × 10−19 erg s−1 cm−2. After discarding broad-line active galactic nuclei, we find 13 He II λ1640 detections from MUSE with a median MUV = −20.1 and 21 tentative He II λ1640 detections from other public surveys. Excluding Lyα, all except two galaxies in our sample show at least one other rest-UV emission line, with C III] λ1907, λ1909 being the most prominent. We use multi-wavelength data available in the Hubble legacy fields to derive basic galaxy properties of our sample through spectral energy distribution fitting techniques. Taking advantage of the high-quality spectra obtained by MUSE (∼10 − 30 h of exposure time per pointing), we use photo-ionisation models to study the rest-UV emission line diagnostics of the He II λ1640 emitters. Line ratios of our sample can be reproduced by moderately sub-solar photo-ionisation models, however, we find that including effects of binary stars lead to degeneracies in most free parameters. Even after considering extra ionising photons produced by extreme sub-solar metallicity binary stellar models, photo-ionisation models are unable to reproduce rest-frame He II λ1640 equivalent widths (∼0.2 − 10 Å), thus additional mechanisms are necessary in models to match the observed He II λ1640 properties.},
  author       = {Nanayakkara, Themiya and Brinchmann, Jarle and Boogaard, Leindert and Bouwens, Rychard and Cantalupo, Sebastiano and Feltre, Anna and Kollatschny, Wolfram and Marino, Raffaella Anna and Maseda, Michael and Matthee, Jorryt J and Paalvast, Mieke and Richard, Johan and Verhamme, Anne},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: ISM / galaxies: star formation / galaxies: evolution / galaxies: high-redshift},
  publisher    = {EDP Sciences},
  title        = {{Exploring He II λ1640 emission line properties at z ∼2−4}},
  doi          = {10.1051/0004-6361/201834565},
  volume       = {648},
  year         = {2019},
}

@article{11505,
  abstract     = {Contact. This paper presents the results obtained with the Multi-Unit Spectroscopic Explorer (MUSE) at the ESO Very Large Telescope on the faint end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The goal of our project is to set strong constraints on the relative contribution of the Lyman-alpha emitter (LAE) population to cosmic reionization.

Aims. The precise aim of the present study is to further constrain the abundance of LAEs by taking advantage of the magnification provided by lensing clusters to build a blindly selected sample of galaxies which is less biased than current blank field samples in redshift and luminosity. By construction, this sample of LAEs is complementary to those built from deep blank fields, whether observed by MUSE or by other facilities, and makes it possible to determine the shape of the LF at fainter levels, as well as its evolution with redshift.

Methods. We selected a sample of 156 LAEs with redshifts between 2.9 ≤ z ≤ 6.7 and magnification-corrected luminosities in the range 39 ≲ log LLyα [erg s−1] ≲43. To properly take into account the individual differences in detection conditions between the LAEs when computing the LF, including lensing configurations, and spatial and spectral morphologies, the non-parametric 1/Vmax method was adopted. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure to properly determine the effective volume Vmax for each galaxy. In this paper we present a complete procedure for the determination of the LF based on IFU detections in lensing clusters. This procedure, including some new methods for masking, effective volume integration and (individual) completeness determinations, has been fully automated when possible, and it can be easily generalized to the analysis of IFU observations in blank fields.

Results. As a result of this analysis, the Lyman-alpha LF has been obtained in four different redshift bins: 2.9 <  z <  6, 7, 2.9 <  z <  4.0, 4.0 <  z <  5.0, and 5.0 <  z <  6.7 with constraints down to log LLyα = 40.5. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis also including data from the literature to complete the present sample towards the brightest luminosities, a steep faint end slope was measured varying from α = −1.69−0.08+0.08 to α = −1.87−0.12+0.12 between the lowest and the highest redshift bins.

Conclusions. The contribution of the LAE population to the star formation rate density at z ∼ 6 is ≲50% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. Depending on the intersection between the LAE/LBG populations, the contribution of the observed galaxies to the ionizing flux may suffice to keep the universe ionized at z ∼ 6.},
  author       = {de La Vieuville, G. and Bina, D. and Pello, R. and Mahler, G. and Richard, J. and Drake, A. B. and Herenz, E. C. and Bauer, F. E. and Clément, B. and Lagattuta, D. and Laporte, N. and Martinez, J. and Patrício, V. and Wisotzki, L. and Zabl, J. and Bouwens, R. J. and Contini, T. and Garel, T. and Guiderdoni, B. and Marino, R. A. and Maseda, M. V. and Matthee, Jorryt J and Schaye, J. and Soucail, G.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, gravitational lensing: strong / galaxies: high-redshift / dark ages, reionization, first stars / galaxies: clusters: general / galaxies: luminosity function, mass function},
  publisher    = {EDP Sciences},
  title        = {{Faint end of the z ∼ 3–7 luminosity function of Lyman-alpha emitters behind lensing clusters observed with MUSE}},
  doi          = {10.1051/0004-6361/201834471},
  volume       = {628},
  year         = {2019},
}

@article{11507,
  abstract     = {Lyman-α (Lyα) is intrinsically the brightest line emitted from active galaxies. While it originates from many physical processes, for star-forming galaxies the intrinsic Lyα luminosity is a direct tracer of the Lyman-continuum (LyC) radiation produced by the most massive O- and early-type B-stars (M⋆ ≳ 10 M⊙) with lifetimes of a few Myrs. As such, Lyα luminosity should be an excellent instantaneous star formation rate (SFR) indicator. However, its resonant nature and susceptibility to dust as a rest-frame UV photon makes Lyα very hard to interpret due to the uncertain Lyα escape fraction, fesc, Lyα. Here we explore results from the CAlibrating LYMan-α with Hα (CALYMHA) survey at z = 2.2, follow-up of Lyα emitters (LAEs) at z = 2.2 − 2.6 and a z ∼ 0−0.3 compilation of LAEs to directly measure fesc, Lyα with Hα. We derive a simple empirical relation that robustly retrieves fesc, Lyα as a function of Lyα rest-frame EW (EW0): fesc,Lyα = 0.0048 EW0[Å] ± 0.05 and we show that it constrains a well-defined anti-correlation between ionisation efficiency (ξion) and dust extinction in LAEs. Observed Lyα luminosities and EW0 are easy measurable quantities at high redshift, thus making our relation a practical tool to estimate intrinsic Lyα and LyC luminosities under well controlled and simple assumptions. Our results allow observed Lyα luminosities to be used to compute SFRs for LAEs at z ∼ 0−2.6 within ±0.2 dex of the Hα dust corrected SFRs. We apply our empirical SFR(Lyα,EW0) calibration to several sources at z ≥ 2.6 to find that star-forming LAEs have SFRs typically ranging from 0.1 to 20 M⊙ yr−1 and that our calibration might be even applicable for the most luminous LAEs within the epoch of re-ionisation. Our results imply high ionisation efficiencies (log10[ξion/Hz erg−1] = 25.4−25.6) and low dust content in LAEs across cosmic time, and will be easily tested with future observations with JWST which can obtain Hα and Hβ measurements for high-redshift LAEs.},
  author       = {Sobral, David and Matthee, Jorryt J},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / galaxies: star formation / galaxies: statistics / galaxies: evolution / galaxies: formation / galaxies: ISM},
  publisher    = {EDP Sciences},
  title        = {{Predicting Lyα escape fractions with a simple observable: Lyα in emission as an empirically calibrated star formation rate indicator}},
  doi          = {10.1051/0004-6361/201833075},
  volume       = {623},
  year         = {2019},
}

@article{11535,
  abstract     = {We investigate the clustering and halo properties of ∼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 ∼ 2.5 and 6. We measure clustering lengths of r0 ∼ 3–6 h−1 Mpc and typical halo masses of ∼1011 M⊙ for our narrowband-selected LAEs with typical LLy α ∼ 1042–43 erg s−1. The intermediate-band-selected LAEs are observed to have r0 ∼ 3.5–15 h−1 Mpc with typical halo masses of ∼1011–12 M⊙ and typical LLy α ∼ 1043–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⋆(z). The faintest LAEs (L ∼ 0.1 L⋆(z)) typically identified by deep narrowband surveys are found in 1010 M⊙ haloes and the brightest LAEs (L ∼ 7 L⋆(z)) are found in ∼5 × 1012 M⊙ haloes. A dependency on the rest-frame 1500 Å UV luminosity, MUV, is also observed where the halo masses increase from 1011 to 1013 M⊙ for MUV ∼ −19 to −23.5 mag. Halo mass is also observed to increase from 109.8 to 1012 M⊙ for dust-corrected UV star formation rates from ∼0.6 to 10 M⊙ yr−1 and continues to increase up to 1013 M⊙ in halo mass, where the majority of those sources are active galactic nuclei. 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.},
  author       = {Khostovan, A A and Sobral, D and Mobasher, B and Matthee, Jorryt J and Cochrane, R K and Chartab, N and Jafariyazani, M and Paulino-Afonso, A and Santos, S and Calhau, J},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: haloes, galaxies: high-redshift, galaxies: star formation, cosmology: observations, large-scale structure of Universe},
  number       = {1},
  pages        = {555--573},
  publisher    = {Oxford University Press},
  title        = {{The clustering of typical Ly α emitters from z ∼ 2.5–6: Host halo masses depend on Ly α and UV luminosities}},
  doi          = {10.1093/mnras/stz2149},
  volume       = {489},
  year         = {2019},
}

@article{11540,
  abstract     = {Observations have revealed that the star formation rate (SFR) and stellar mass (Mstar) of star-forming galaxies follow a tight relation known as the galaxy main sequence. However, what physical information is encoded in this relation is under debate. Here, we use the EAGLE cosmological hydrodynamical simulation to study the mass dependence, evolution, and origin of scatter in the SFR–Mstar relation. At z = 0, we find that the scatter decreases slightly with stellar mass from 0.35 dex at Mstar ≈ 109 M⊙ to 0.30 dex at Mstar ≳ 1010.5 M⊙. The scatter decreases from z = 0 to z = 5 by 0.05 dex at Mstar ≳ 1010 M⊙ and by 0.15 dex for lower masses. We show that the scatter at z = 0.1 originates from a combination of fluctuations on short time-scales (ranging from 0.2–2 Gyr) that are presumably associated with self-regulation from cooling, star formation, and outflows, but is dominated by long time-scale (∼10 Gyr) variations related to differences in halo formation times. Shorter time-scale fluctuations are relatively more important for lower mass galaxies. At high masses, differences in black hole formation efficiency cause additional scatter, but also diminish the scatter caused by different halo formation times. While individual galaxies cross the main sequence multiple times during their evolution, they fluctuate around tracks associated with their halo properties, i.e. galaxies above/below the main sequence at z = 0.1 tend to have been above/below the main sequence for ≫1 Gyr.},
  author       = {Matthee, Jorryt J and Schaye, Joop},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics : galaxies: evolution, galaxies: formation, galaxies: star formation, cosmology: theory},
  number       = {1},
  pages        = {915--932},
  publisher    = {Oxford University Press},
  title        = {{The origin of scatter in the star formation rate–stellar mass relation}},
  doi          = {10.1093/mnras/stz030},
  volume       = {484},
  year         = {2019},
}

@article{11541,
  abstract     = {We present new Hubble Space Telescope (HST)/WFC3 observations and re-analyse VLT data to unveil the continuum, variability, and rest-frame ultraviolet (UV) lines of the multiple UV clumps of the most luminous Lyα emitter at z = 6.6, CR7 (COSMOS Redshift 7). Our re-reduced, flux-calibrated X-SHOOTER spectra of CR7 reveal an He II emission line in observations obtained along the major axis of Lyα emission with the best seeing conditions. He II is spatially offset by ≈+0.8 arcsec from the peak of Lyα emission, and it is found towards clump B. Our WFC3 grism spectra detects the UV continuum of CR7’s clump A, yielding a power law with β=−2.5+0.6−0.7 and MUV=−21.87+0.25−0.20⁠. No significant variability is found for any of the UV clumps on their own, but there is tentative (≈2.2 σ) brightening of CR7 in F110W as a whole from 2012 to 2017. HST grism data fail to robustly detect rest-frame UV lines in any of the clumps, implying fluxes ≲2×10−17 erg s−1 cm−2 (3σ). We perform CLOUDY modelling to constrain the metallicity and the ionizing nature of CR7. CR7 seems to be actively forming stars without any clear active galactic nucleus activity in clump A, consistent with a metallicity of ∼0.05–0.2 Z⊙. Component C or an interclump component between B and C may host a high ionization source. Our results highlight the need for spatially resolved information to study the formation and assembly of early galaxies.},
  author       = {Sobral, David and Matthee, Jorryt J and Brammer, Gabriel and Ferrara, Andrea and Alegre, Lara and Röttgering, Huub and Schaerer, Daniel and Mobasher, Bahram and Darvish, Behnam},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: high-redshift, galaxies: ISM, cosmology: observations, dark ages, reionization, first stars, early Universe},
  number       = {2},
  pages        = {2422--2441},
  publisher    = {Oxford University Press},
  title        = {{On the nature and physical conditions of the luminous Ly α emitter CR7 and its rest-frame UV components}},
  doi          = {10.1093/mnras/sty2779},
  volume       = {482},
  year         = {2019},
}

@article{11508,
  abstract     = {Distant luminous Lyman-α 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. We rule out the possibility 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 ± 0.03) and ii) an unphysical [OII]/Hα ratio ([OII]/Hα >  22). We show that COLA1’s observed B-band flux is explained by a faint extended foreground LAE, for which we detect Lyα 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. COLA1 is UV luminous (M1500 = −21.6 ± 0.3), has a high equivalent width (EW0,Lyα = 120−40+50 Å) and very compact Lyα emission (r50,Lyα = 0.33−0.04+0.07 kpc). Relatively weak inferred Hβ+[OIII] line-emission from Spitzer/IRAC indicates an extremely low metallicity of Z <  1/20 Z⊙ or reduced strength of nebular lines due to high escape of ionising photons. The small Lyα peak separation of 220 ± 20 km s−1 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 (2.3 cMpc) bubble, allowing the blue Lyα 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.},
  author       = {Matthee, Jorryt J and Sobral, David and Gronke, Max and Paulino-Afonso, Ana and Stefanon, Mauro and Röttgering, Huub},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / galaxies: formation / dark ages / reionization / first stars / techniques: spectroscopic / intergalactic medium},
  publisher    = {EDP Sciences},
  title        = {{Confirmation of double peaked Lyα emission at z = 6.593: Witnessing a galaxy directly contributing to the reionisation of the universe}},
  doi          = {10.1051/0004-6361/201833528},
  volume       = {619},
  year         = {2018},
}

@article{11549,
  abstract     = {We investigate the clustering properties of ∼7000 H β + [O III] and [O II] narrowband-selected emitters at z ∼ 0.8–4.7 from the High-z Emission Line Survey. We find clustering lengths, r0, of 1.5–4.0 h−1 Mpc and minimum dark matter halo masses of 1010.7–12.1 M⊙ for our z = 0.8–3.2 H β + [O III] emitters and r0 ∼ 2.0–8.3 h−1 Mpc and halo masses of 1011.5–12.6 M⊙ for our z = 1.5–4.7 [O II] emitters. We find r0 to strongly increase both with increasing line luminosity and redshift. By taking into account the evolution of the characteristic line luminosity, L⋆(z), and using our model predictions of halo mass given r0, we find a strong, redshift-independent increasing trend between L/L⋆(z) and minimum halo mass. The faintest H β + [O III] emitters are found to reside in 109.5 M⊙ haloes and the brightest emitters in 1013.0 M⊙ haloes. For [O II] emitters, the faintest emitters are found in 1010.5 M⊙ haloes and the brightest emitters in 1012.6 M⊙ haloes. A redshift-independent stellar mass dependency is also observed where the halo mass increases from 1011 to 1012.5 M⊙ for stellar masses of 108.5 to 1011.5 M⊙, respectively. We investigate the interdependencies of these trends by repeating our analysis in a Lline−Mstar grid space for our most populated samples (H β + [O III] z = 0.84 and [O II] z = 1.47) and find that the line luminosity dependency is stronger than the stellar mass dependency on halo mass. For L > L⋆ emitters at all epochs, we find a relatively flat trend with halo masses of 1012.5–13 M⊙, which may be due to quenching mechanisms in massive haloes that is consistent with a transitional halo mass predicted by models.},
  author       = {Khostovan, A A and Sobral, D and Mobasher, B and Best, P N and Smail, I and Matthee, Jorryt J and Darvish, B and Nayyeri, H and Hemmati, S and Stott, J P},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: haloes, galaxies: high-redshift, galaxies: star formation, cosmology: observations, large-scale structure of Universe},
  number       = {3},
  pages        = {2999--3015},
  publisher    = {Oxford University Press},
  title        = {{The clustering of H β + [O III] and [O II] emitters since z ∼ 5: Dependencies with line luminosity and stellar mass}},
  doi          = {10.1093/mnras/sty925},
  volume       = {478},
  year         = {2018},
}

@article{11555,
  abstract     = {We investigate the morphology of the [C II] emission in a sample of ‘normal’ star-forming galaxies at 5 < z < 7.2 in relation to their UV (rest-frame) counterpart. We use new Atacama Large Millimetre/submillimetre Array (ALMA) observations of galaxies at z ∼ 6–7, as well as a careful re-analysis of archival ALMA data. In total 29 galaxies were analysed, 21 of which are detected in [C II]. For several of the latter the [C II] emission breaks into multiple components. Only a fraction of these [C II] components, if any, is associated with the primary UV systems, while the bulk of the [C II] emission is associated either with fainter UV components, or not associated with any UV counterpart at the current limits. By taking into account the presence of all these components, we find that the L[CII]–SFR (star formation rate) relation at early epochs is fully consistent with the local relation, but it has a dispersion of 0.48 ± 0.07 dex, which is about two times larger than observed locally. We also find that the deviation from the local L[CII]–SFR relation has a weak anticorrelation with the EW(Ly α). The morphological analysis also reveals that [C II] emission is generally much more extended than the UV emission. As a consequence, these primordial galaxies are characterized by a [C II] surface brightness generally much lower than expected from the local Σ[CII]−ΣSFR relation. These properties are likely a consequence of a combination of different effects, namely gas metallicity, [C II] emission from obscured star-forming regions, strong variations of the ionization parameter, and circumgalactic gas in accretion or ejected by these primeval galaxies.},
  author       = {Carniani, S and Maiolino, R and Amorin, R and Pentericci, L and Pallottini, A and Ferrara, A and Willott, C J and Smit, R and Matthee, Jorryt J and Sobral, D and Santini, P and Castellano, M and De Barros, S and Fontana, A and Grazian, A and Guaita, L},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: high-redshift, galaxies: ISM, galaxies: formation},
  number       = {1},
  pages        = {1170--1184},
  publisher    = {Oxford University Press},
  title        = {{Kiloparsec-scale gaseous clumps and star formation at z = 5–7}},
  doi          = {10.1093/mnras/sty1088},
  volume       = {478},
  year         = {2018},
}

@article{11557,
  abstract     = {Deep narrow-band surveys have revealed a large population of faint Ly α emitters (LAEs) in the distant Universe, but relatively little is known about the most luminous sources (⁠LLyα≳1042.7 erg s−1; LLyα≳L∗Lyα⁠). 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 (≈4 × 106 Mpc3 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 ionization UV lines. Luminous LAEs at z ∼ 2–3 have blue UV slopes (⁠β=−2.0+0.3−0.1⁠) and high Ly α escape fractions (⁠50+20−15 per cent) and span five orders of magnitude in UV luminosity (MUV ≈ −19 to −24). Many (70 per cent) show at least one high ionization rest-frame UV line such as C IV, N V, C III], He II or O III], typically blue-shifted by ≈100–200 km s−1 relative to Ly α. Their Ly α profiles reveal a wide variety of shapes, including significant blue-shifted components and widths from 200 to 4000 km s−1. Overall, 60 ± 11  per cent appear to be active galactic nucleus (AGN) dominated, and at LLyα > 1043.3 erg s−1 and/or MUV < −21.5 virtually all LAEs are AGNs with high ionization parameters (log U = 0.6 ± 0.5) and with metallicities of ≈0.5 − 1 Z⊙. Those lacking signatures of AGNs (40 ± 11  per cent) have lower ionization parameters (⁠logU=−3.0+1.6−0.9 and log ξion = 25.4 ± 0.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 2×L∗Lyα and 2×M∗UV mark a sharp transition in the nature of LAEs, from star formation dominated to AGN dominated.},
  author       = {Sobral, David and Matthee, Jorryt J and Darvish, Behnam and Smail, Ian and Best, Philip N and Alegre, Lara and Röttgering, Huub and Mobasher, Bahram and Paulino-Afonso, Ana and Stroe, Andra and Oteo, Iván},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: active, galaxies: evolution, galaxies: high-redshift, galaxies: ISM, galaxies: starburst, cosmology: observations},
  number       = {2},
  pages        = {2817--2840},
  publisher    = {Oxford University Press},
  title        = {{The nature of luminous Ly α emitters at z ∼ 2–3: Maximal dust-poor starbursts and highly ionizing AGN}},
  doi          = {10.1093/mnras/sty782},
  volume       = {477},
  year         = {2018},
}

@article{11558,
  abstract     = {We present and explore deep narrow- and medium-band data obtained with the Subaru and the Isaac Newton Telescopes in the ∼2 deg2 COSMOS field. We use these data as an extremely wide, low-resolution (R ∼ 20–80) Integral Field Unit survey to slice through the COSMOS field and obtain a large sample of ∼4000 Ly α emitters (LAEs) from z ∼ 2 to 6 in 16 redshift slices (SC4K). We present new Ly α luminosity functions (LFs) covering a comoving volume of ∼108 Mpc3. SC4K extensively complements ultradeep surveys, jointly covering over 4 dex in Ly α luminosity and revealing a global (2.5 < z < 6) synergy LF with α=−1.93+0.12−0.12⁠, log10Φ∗Lyα=−3.45+0.22−0.29 Mpc−3, and log10L∗Lyα=42.93+0.15−0.11 erg s−1. The Schechter component of the Ly α LF reveals a factor ∼5 rise in L∗Lyα and a ∼7 × decline in Φ∗Lyα from z ∼ 2 to 6. The data reveal an extra power-law (or Schechter) component above LLy α ≈ 1043.3 erg s−1 at z ∼ 2.2–3.5 and we show that it is partially driven by X-ray and radio active galactic nucleus (AGN), as their Ly α LF resembles the excess. The power-law component vanishes and/or is below our detection limits above z > 3.5, likely linked with the evolution of the AGN population. The Ly α luminosity density rises by a factor ∼2 from z ∼ 2 to 3 but is then found to be roughly constant (⁠1.1+0.2−0.2×1040 erg s−1 Mpc−3) to z ∼ 6, despite the ∼0.7 dex drop in ultraviolet (UV) luminosity density. The Ly α/UV luminosity density ratio rises from 4 ± 1 per cent to 30 ± 6 per cent from z ∼ 2.2 to 6. Our results imply a rise of a factor of ≈2 in the global ionization efficiency (ξion) and a factor ≈4 ± 1 in the Ly α escape fraction from z ∼ 2 to 6, hinting for evolution in both the typical burstiness/stellar populations and even more so in the typical interstellar medium conditions allowing Ly α photons to escape.},
  author       = {Sobral, David and Santos, Sérgio and Matthee, Jorryt J and Paulino-Afonso, Ana and Ribeiro, Bruno and Calhau, João and Khostovan, Ali A},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: formation, galaxies: high-redshift, galaxies: luminosity function, mass function, galaxies: statistics},
  number       = {4},
  pages        = {4725--4752},
  publisher    = {Oxford University Press},
  title        = {{Slicing COSMOS with SC4K: The evolution of typical Ly α emitters and the Ly α escape fraction from z ∼ 2 to 6}},
  doi          = {10.1093/mnras/sty378},
  volume       = {476},
  year         = {2018},
}

@article{11584,
  abstract     = {Observations show that star-forming galaxies reside on a tight 3D plane between mass, gas-phase metallicity, and star formation rate (SFR), which can be explained by the interplay between metal-poor gas inflows, SFR and outflows. However, different metals are released on different time-scales, which may affect the slope of this relation. Here, we use central, star-forming galaxies with Mstar = 109.0–10.5 M⊙ from the EAGLE hydrodynamical simulation to examine 3D relations between mass, SFR, and chemical enrichment using absolute and relative C, N, O, and Fe abundances. We show that the scatter is smaller when gas-phase α-enhancement is used rather than metallicity. A similar plane also exists for stellar α-enhancement, implying that present-day specific SFRs are correlated with long time-scale star formation histories. Between z = 0 and 1, the α-enhancement plane is even more insensitive to redshift than the plane using metallicity. However, it evolves at z > 1 due to lagging iron yields. At fixed mass, galaxies with higher SFRs have star formation histories shifted towards late times, are more α-enhanced, and this α-enhancement increases with redshift as observed. These findings suggest that relations between physical properties inferred from observations may be affected by systematic variations in α-enhancements.},
  author       = {Matthee, Jorryt J and Schaye, Joop},
  issn         = {1745-3933},
  journal      = {Monthly Notices of the Royal Astronomical Society: Letters},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: abundances, galaxies: evolution, galaxies: formation, galaxies: star formation},
  number       = {1},
  pages        = {L34 -- L39},
  publisher    = {Oxford University Press},
  title        = {{Star-forming galaxies are predicted to lie on a fundamental plane of mass, star formation rate, and α-enhancement}},
  doi          = {10.1093/mnrasl/sly093},
  volume       = {479},
  year         = {2018},
}

@article{11518,
  abstract     = {We present spectroscopic follow-up observations of CR7 with ALMA, targeted at constraining the infrared (IR) continuum and [C II]158 mm line-emission at high spatial resolution matched to the HST/WFC3 imaging. CR7 is a luminous Lyα emitting galaxy at z = 6.6 that consists of three separated UV-continuum components. Our observations reveal several well-separated components of [C II] emission. The two most luminous components in [C II] coincide with the brightest UV components (A and B), blueshifted by »150 km s−1 with respect to the
peak of Lyα emission. Other [C II] components are observed close to UV clumps B and C and are blueshifted by »300 and ≈80 km s−1 with respect to the systemic redshift. We do not detect FIR continuum emission due to dust with a 3σ limiting luminosity LIR T L d 35 K 3.1 10 = <´ 10 ( ) . This allows us to mitigate uncertainties in the dust-corrected SFR and derive SFRs for the three UV clumps A, B, and C of 28, 5, and 7 M yr−1. All clumps have [C II] luminosities consistent within the scatter observed in the local relation between SFR and L[ ] C II , implying that strong Lyα emission does not necessarily anti-correlate with [C II] luminosity. Combining
our measurements with the literature, we show that galaxies with blue UV slopes have weaker [C II] emission at fixed SFR, potentially due to their lower metallicities and/or higher photoionization. Comparison with hydrodynamical simulations suggests that CR7ʼs clumps have metallicities of 0.1 Z Z 0.2 < < . The observed ISM structure of CR7 indicates that we are likely witnessing the build up of a central galaxy in the early universe through complex accretion of satellites.},
  author       = {Matthee, Jorryt J and Sobral, D. and Boone, F. and Röttgering, H. and Schaerer, D. and Girard, M. and Pallottini, A. and Vallini, L. and Ferrara, A. and Darvish, B. and Mobasher, B.},
  issn         = {1538-4357},
  journal      = {The Astrophysical Journal},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, dark ages, reionization, first stars – galaxies: formation – galaxies: high-redshift – galaxies: ISM – galaxies: kinematics and dynamics},
  number       = {2},
  publisher    = {IOP Publishing},
  title        = {{ALMA reveals metals yet no dust within multiple components in CR7}},
  doi          = {10.3847/1538-4357/aa9931},
  volume       = {851},
  year         = {2017},
}

@article{11561,
  abstract     = {We present a sample of ∼1000 emission-line galaxies at z = 0.4–4.7 from the ∼0.7deg2 High-z Emission-Line Survey in the Boötes field identified with a suite of six narrow-band filters at ≈0.4–2.1 μm. These galaxies have been selected on their Ly α (73), [O II] (285), H β/[O III] (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 II]–H α) detections and/or spectroscopy. In this survey paper, we present the observations, selection and catalogues of emitters. We measure number densities of Ly α, [O II], H β/[O III] and H α and confirm strong luminosity evolution in star-forming galaxies from z ∼ 0.4 to ∼5, in agreement with previous results. To demonstrate the usefulness of dual-line emitters, we use the sample of dual [O II]–H α emitters to measure the observed [O II]/H α ratio at z = 1.47. The observed [O II]/H α ratio increases significantly from 0.40 ± 0.01 at z = 0.1 to 0.52 ± 0.05 at z = 1.47, which we attribute to either decreasing dust attenuation with redshift, or due to a bias in the (typically) fibre measurements in the local Universe that only measure the central kpc regions. At the bright end, we find that both the H α and Ly α number densities at z ≈ 2.2 deviate significantly from a Schechter form, following a power law. We show that this is driven entirely by an increasing X-ray/active galactic nucleus fraction with line luminosity, which reaches ≈100 per cent at line luminosities L ≳ 3 × 1044 erg s−1.},
  author       = {Matthee, Jorryt J and Sobral, David and Best, Philip and Smail, Ian and Bian, Fuyan and Darvish, Behnam and Röttgering, Huub and Fan, Xiaohui},
  issn         = {0035-8711},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics galaxies, active, galaxies, evolution, galaxies, high-redshift, galaxies, luminosity function, mass function, galaxies: star formation},
  number       = {1},
  pages        = {629--649},
  publisher    = {Oxford University Press},
  title        = {{Boötes-HiZELS: An optical to near-infrared survey of emission-line galaxies at z = 0.4–4.7}},
  doi          = {10.1093/mnras/stx1569},
  volume       = {471},
  year         = {2017},
}

@article{11562,
  abstract     = {We present the CAlibrating LYMan-α with Hα (CALYMHA) pilot survey and new results on Lyman α (Lyα) selected galaxies at z ∼ 2. We use a custom-built Lyα narrow-band filter at the Isaac Newton Telescope, designed to provide a matched volume coverage to the z = 2.23 Hα HiZELS survey. Here, we present the first results for the COSMOS and UDS fields. Our survey currently reaches a 3σ line flux limit of ∼4 × 10−17 erg s−1 cm−2, and a Lyα luminosity limit of ∼1042.3 erg s−1. We find 188 Lyα emitters over 7.3 × 105 Mpc3, but also find significant numbers of other line-emitting sources corresponding to He II, C III] and C IV emission lines. These sources are important contaminants, and we carefully remove them, unlike most previous studies. We find that the Lyα luminosity function at z = 2.23 is very well described by a Schechter function up to LLy α ≈ 1043 erg s−1 with L∗=1042.59+0.16−0.08 erg s−1, ϕ∗=10−3.09+0.14−0.34 Mpc−3 and α = −1.75 ± 0.25. Above LLy α ≈ 1043 erg s−1, the Lyα luminosity function becomes power-law like, driven by X-ray AGN. We find that Lyα-selected emitters have a high escape fraction of 37 ± 7 per cent, anticorrelated with Lyα luminosity and correlated with Lyα equivalent width. Lyα emitters have ubiquitous large (≈40 kpc) Lyα haloes, ∼2 times larger than their Hα extents. By directly comparing our Lyα and Hα luminosity functions, we find that the global/overall escape fraction of Lyα photons (within a 13 kpc radius) from the full population of star-forming galaxies is 5.1 ± 0.2 per cent at the peak of the star formation history. An extra 3.3 ± 0.3 per cent of Lyα photons likely still escape, but at larger radii.},
  author       = {Sobral, David and Matthee, Jorryt J and Best, Philip and Stroe, Andra and Röttgering, Huub and Oteo, Iván and Smail, Ian and Morabito, Leah and Paulino-Afonso, Ana},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: haloes, galaxies: high-redshift, galaxies: luminosity function, mass function, galaxies: statistics, cosmology: observations},
  number       = {1},
  pages        = {1242--1258},
  publisher    = {Oxford University Press},
  title        = {{The CALYMHA survey: Lyα luminosity function and global escape fraction of Lyα photons at z = 2.23}},
  doi          = {10.1093/mnras/stw3090},
  volume       = {466},
  year         = {2017},
}

@article{11564,
  abstract     = {We study the production rate of ionizing photons of a sample of 588 Hα emitters (HAEs) and 160 Lyman-α emitters (LAEs) at z = 2.2 in the COSMOS field in order to assess the implied emissivity from galaxies, based on their ultraviolet (UV) luminosity. By exploring the rest-frame Lyman Continuum (LyC) with GALEX/NUV data, we find fesc < 2.8 (6.4) per cent through median (mean) stacking. By combining the Hα luminosity density with intergalactic medium emissivity measurements from absorption studies, we find a globally averaged 〈fesc〉 of 5.9+14.5−4.2 per cent at z = 2.2 if we assume HAEs are the only source of ionizing photons. We find similarly low values of the global 〈fesc〉 at z ≈ 3–5, also ruling out a high 〈fesc〉 at z < 5. These low escape fractions allow us to measure ξion, the number of produced ionizing photons per unit UV luminosity, and investigate how this depends on galaxy properties. We find a typical ξion ≈ 1024.77 ± 0.04 Hz erg−1 for HAEs and ξion ≈ 1025.14 ± 0.09 Hz erg−1 for LAEs. LAEs and low-mass HAEs at z = 2.2 show similar values of ξion as typically assumed in the reionization era, while the typical HAE is three times less ionizing. Due to an increasing ξion with increasing EW(Hα), ξion likely increases with redshift. This evolution alone is fully in line with the observed evolution of ξion between z ≈ 2 and 5, indicating a typical value of ξion ≈ 1025.4 Hz erg−1 in the reionization era.},
  author       = {Matthee, Jorryt J and Sobral, David and Best, Philip and Khostovan, Ali Ahmad and Oteo, Iván and Bouwens, Rychard and Röttgering, Huub},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: high-redshift, cosmology: observations, dark ages, reionization, first stars},
  number       = {3},
  pages        = {3637--3655},
  publisher    = {Oxford University Press},
  title        = {{The production and escape of Lyman-Continuum radiation from star-forming galaxies at z ∼ 2 and their redshift evolution}},
  doi          = {10.1093/mnras/stw2973},
  volume       = {465},
  year         = {2017},
}

@article{11565,
  abstract     = {We use the hydrodynamical EAGLE simulation to study the magnitude and origin of the scatter in the stellar mass–halo mass relation for central galaxies. We separate cause and effect by correlating stellar masses in the baryonic simulation with halo properties in a matched dark matter only (DMO) simulation. The scatter in stellar mass increases with redshift and decreases with halo mass. At z = 0.1, it declines from 0.25 dex at M200, DMO ≈ 1011 M⊙ to 0.12 dex at M200, DMO ≈ 1013 M⊙, but the trend is weak above 1012 M⊙. For M200, DMO < 1012.5 M⊙ up to 0.04 dex of the scatter is due to scatter in the halo concentration. At fixed halo mass, a larger stellar mass corresponds to a more concentrated halo. This is likely because higher concentrations imply earlier formation times and hence more time for accretion and star formation, and/or because feedback is less efficient in haloes with higher binding energies. The maximum circular velocity, Vmax, DMO, and binding energy are therefore more fundamental properties than halo mass, meaning that they are more accurate predictors of stellar mass, and we provide fitting formulae for their relations with stellar mass. However, concentration alone cannot explain the total scatter in the Mstar−M200,DMO relation, and it does not explain the scatter in Mstar–Vmax, DMO. Halo spin, sphericity, triaxiality, substructure and environment are also not responsible for the remaining scatter, which thus could be due to more complex halo properties or non-linear/stochastic baryonic effects.},
  author       = {Matthee, Jorryt J and Schaye, Joop and Crain, Robert A. and Schaller, Matthieu and Bower, Richard and Theuns, Tom},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: formation, galaxies: haloes, cosmology: theory},
  number       = {2},
  pages        = {2381--2396},
  publisher    = {Oxford University Press},
  title        = {{The origin of scatter in the stellar mass–halo mass relation of central galaxies in the EAGLE simulation}},
  doi          = {10.1093/mnras/stw2884},
  volume       = {465},
  year         = {2017},
}

@article{11566,
  abstract     = {While traditionally associated with active galactic nuclei (AGN), the properties of the C II] (λ = 2326 Å), C III] (λ, λ = 1907, 1909 Å) and C IV (λ, λ = 1549, 1551 Å) emission lines are still uncertain as large, unbiased samples of sources are scarce. We present the first blind, statistical study of C II], C III] and C IV emitters at z ∼ 0.68, 1.05, 1.53, respectively, uniformly selected down to a flux limit of ∼4 × 10−17 erg s−1 cm−1 through a narrow-band survey covering an area of ∼1.4 deg2 over COSMOS and UDS. We detect 16 C II], 35 C III] and 17 C IV emitters, whose nature we investigate using optical colours as well as Hubble Space Telescope (HST), X-ray, radio and far-infrared data. We find that z ∼ 0.7 C II] emitters are consistent with a mixture of blue (UV slope β = −2.0 ± 0.4) star-forming (SF) galaxies with discy HST structure and AGN with Seyfert-like morphologies. Bright C II] emitters have individual X-ray detections as well as high average black hole accretion rates (BHARs) of ∼0.1 M⊙ yr−1. C III] emitters at z ∼ 1.05 trace a general population of SF galaxies, with β = −0.8 ± 1.1, a variety of optical morphologies, including isolated and interacting galaxies and low BHAR (<0.02 M⊙ yr−1). Our C IV emitters at z ∼ 1.5 are consistent with young, blue quasars (β ∼ −1.9) with point-like optical morphologies, bright X-ray counterparts and large BHAR (0.8  M⊙ yr−1). We also find some surprising C II], C III] and C IV emitters with rest-frame equivalent widths (EWs) that could be as large as 50–100 Å. AGN or spatial offsets between the UV continuum stellar disc and the line-emitting regions may explain the large EW. These bright C II], C III] and C IV emitters are ideal candidates for spectroscopic follow-up to fully unveil their nature.},
  author       = {Stroe, Andra and Sobral, David and Matthee, Jorryt J and Calhau, João and Oteo, Ivan},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: active, galaxies: high-redshift, quasars: emission lines, galaxies: star formation, cosmology: observations},
  number       = {3},
  pages        = {2558--2574},
  publisher    = {Oxford University Press},
  title        = {{A 1.4 deg2 blind survey for C II], C III] and C IV at z ∼ 0.7–1.5 – I. Nature, morphologies and equivalent widths }},
  doi          = {10.1093/mnras/stx1712},
  volume       = {471},
  year         = {2017},
}

@article{11567,
  abstract     = {Recently, the C III] and C IV 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 active galactic nuclei (AGN). We explore the first blind C II], C III] and C IV survey (z ∼ 0.68, 1.05, 1.53, respectively) presented in Stroe et al. (2017). We derive luminosity functions (LF) and study properties of C II], C III] and C IV line emitters through comparisons to the LFs of H α and Ly α emitters, UV selected star-forming (SF) galaxies and quasars at similar redshifts. The C II] LF at z ∼ 0.68 is equally well described by a Schechter or a power-law LF, characteristic of a mixture of SF and AGN activity. The C III] LF (z ∼ 1.05) is consistent to a scaled down version of the Schechter H α and Ly α LF at their redshift, indicating a SF origin. In stark contrast, the C IV LF at z ∼ 1.53 is well fit by a power-law, quasar-like LF. We find that the brightest UV sources (MUV < −22) will universally have C III] and C IV emission. However, on average, C III] and C IV are not as abundant as H α or Ly α emitters at the same redshift, with cosmic average ratios of ∼0.02–0.06 to H α and ∼0.01–0.1 to intrinsic Ly α. We predict that the C III] and C IV lines can only be truly competitive in confirming high-redshift candidates when the hosts are intrinsically bright and the effective Ly α escape fraction is below 1 per cent. While C III] and C IV were proposed as good tracers of young, relatively low-metallicity galaxies typical of the early Universe, we find that, at least at z ∼ 1.5, C IV is exclusively hosted by AGN/quasars, especially at large line equivalent widths.},
  author       = {Stroe, Andra and Sobral, David and Matthee, Jorryt J and Calhau, João and Oteo, Ivan},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: active, galaxies: high redshift, galaxies: luminosity function, mass function, quasars: emission lines, star formation, cosmology: observations},
  number       = {3},
  pages        = {2575--2586},
  publisher    = {Oxford University Press},
  title        = {{A 1.4 deg2 blind survey for C II], C III] and C IV at z ∼ 0.7–1.5 – II. Luminosity functions and cosmic average line ratios}},
  doi          = {10.1093/mnras/stx1713},
  volume       = {471},
  year         = {2017},
}

