@article{11488,
  abstract     = {Hydrogen Lyα haloes (LAHs) are commonly used as a tracer of the circumgalactic medium (CGM) at high redshifts. In this work, we aim to explore the existence of Lyα haloes around individual UV-selected galaxies, rather than around Lyα emitters (LAEs), at high redshifts. Our sample was continuum-selected with F775W ≤ 27.5, and spectroscopic redshifts were assigned or constrained for all the sources thanks to the deepest (100- to 140-h) existing Very Large Telescope (VLT)/Multi-Unit Spectroscopic Explorer (MUSE) data with adaptive optics. The final sample includes 21 galaxies that are purely F775W-magnitude selected within the redshift range z ≈ 2.9 − 4.4 and within a UV magnitude range −20 ≤ M1500 ≤ −18, thus avoiding any bias toward LAEs. We tested whether galaxy’s Lyα emission is significantly more extended than the MUSE PSF-convolved continuum component. We find 17 LAHs and four non-LAHs. We report the first individual detections of extended Lyα emission around non-LAEs. The Lyα halo fraction is thus as high as 81.0−11.2+10.3%, which is close to that for LAEs at z = 3 − 6 in the literature. This implies that UV-selected galaxies generally have a large amount of hydrogen in their CGM. We derived the mean surface brightness (SB) profile for our LAHs with cosmic dimming corrections and find that Lyα emission extends to 5.4 arcsec (≃40 physical kpc at the midpoint redshift z = 3.6) above the typical 1σ SB limit. The incidence rate of surrounding gas detected in Lyα per one-dimensional line of sight per unit redshift, dn/dz, is estimated to be 0.76−0.09+0.09 for galaxies with M1500 ≤ −18 mag at z ≃ 3.7. Assuming that Lyα emission and absorption arise in the same gas, this suggests, based on abundance matching, that LAHs trace the same gas as damped Lyα systems (DLAs) and sub-DLAs.},
  author       = {Kusakabe, Haruka and Verhamme, Anne and Blaizot, Jérémy and Garel, Thibault and Wisotzki, Lutz and Leclercq, Floriane and Bacon, Roland and Schaye, Joop and Gallego, Sofia G. and Kerutt, Josephine and Matthee, Jorryt J and Maseda, Michael and Nanayakkara, Themiya and Pelló, Roser and Richard, Johan and Tresse, Laurence and Urrutia, Tanya and Vitte, Eloïse},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / galaxies: formation / galaxies: evolution / cosmology: observations},
  publisher    = {EDP Sciences},
  title        = {{The MUSE eXtremely Deep Field: Individual detections of Ly<i>α</i> haloes around rest-frame UV-selected galaxies at <i>z</i> ≃ 2.9–4.4}},
  doi          = {10.1051/0004-6361/202142302},
  volume       = {660},
  year         = {2022},
}

@article{11497,
  abstract     = {Context. The hydrogen Lyman α line is often the only measurable feature in optical spectra of high-redshift galaxies. Its shape and strength are influenced by radiative transfer processes and the properties of the underlying stellar population. High equivalent widths of several hundred Å are especially hard to explain by models and could point towards unusual stellar populations, for example with low metallicities, young stellar ages, and a top-heavy initial mass function. Other aspects influencing equivalent widths are the morphology of the galaxy and its gas properties.
Aims. The aim of this study is to better understand the connection between the Lyman α rest-frame equivalent width (EW0) and spectral properties as well as ultraviolet (UV) continuum morphology by obtaining reliable EW0 histograms for a statistical sample of galaxies and by assessing the fraction of objects with large equivalent widths.
Methods. We used integral field spectroscopy from the Multi Unit Spectroscopic Explorer (MUSE) combined with broad-band data from the Hubble Space Telescope (HST) to measure EW0. We analysed the emission lines of 1920 Lyman α emitters (LAEs) detected in the full MUSE-Wide (one hour exposure time) and MUSE-Deep (ten hour exposure time) surveys and found UV continuum counterparts in archival HST data. We fitted the UV continuum photometric images using the Galfit software to gain morphological information on the rest-UV emission and fitted the spectra obtained from MUSE to determine the double peak fraction, asymmetry, full-width at half maximum, and flux of the Lyman α line.
Results. The two surveys show different histograms of Lyman α EW0. In MUSE-Wide, 20% of objects have EW0 > 240 Å, while this fraction is only 11% in MUSE-Deep and ≈16% for the full sample. This includes objects without HST continuum counterparts (one-third of our sample), for which we give lower limits for EW0. The object with the highest securely measured EW0 has EW0 = 589 ± 193 Å (the highest lower limit being EW0 = 4464 Å). We investigate the connection between EW0 and Lyman α spectral or UV continuum morphological properties.
Conclusions. The survey depth has to be taken into account when studying EW0 distributions. We find that in general, high EW0 objects can have a wide range of spectral and UV morphological properties, which might reflect that the underlying causes for high EW0 values are equally varied.},
  author       = {Kerutt, J. and Wisotzki, L. and Verhamme, A. and Schmidt, K. B. and Leclercq, F. and Herenz, E. C. and Urrutia, T. and Garel, T. and Hashimoto, T. and Maseda, M. and Matthee, Jorryt J and Kusakabe, H. and Schaye, J. and Richard, J. and Guiderdoni, B. and Mauerhofer, V. and Nanayakkara, T. and Vitte, E.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / galaxies: formation / galaxies: evolution / cosmology: observations},
  publisher    = {EDP Sciences},
  title        = {{Equivalent widths of Lyman α emitters in MUSE-Wide and MUSE-Deep}},
  doi          = {10.1051/0004-6361/202141900},
  volume       = {659},
  year         = {2022},
}

@article{11524,
  abstract     = {We measure the evolution of the rest-frame UV luminosity function (LF) and the stellar mass function (SMF) of Lyman-α (Ly α) emitters (LAEs) from z ∼ 2 to z ∼ 6 by exploring ∼4000 LAEs from the SC4K sample. We find a correlation between Ly α luminosity (LLy α) and rest-frame UV (MUV), with best fit MUV=−1.6+0.2−0.3log10(LLyα/ergs−1)+47+12−11 and a shallower relation between LLy α and stellar mass (M⋆), with best fit log10(M⋆/M⊙)=0.9+0.1−0.1log10(LLyα/ergs−1)−28+4.0−3.8⁠. An increasing LLy α cut predominantly lowers the number density of faint MUV and low M⋆ 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 (Φ*). For the SMF, we measure a characteristic stellar mass (⁠M∗⋆/M⊙⁠) increase with increasing redshift, and a Φ* decline. However, if we apply a uniform luminosity cut of log10(LLyα/ergs−1)≥43.0⁠, we find much milder to no evolution in the UV and SMF of LAEs. The UV luminosity density (ρUV) of the full sample of LAEs shows moderate evolution and the stellar mass density (ρM) decreases, with both being always lower than the total ρUV and ρM of more typical galaxies but slowly approaching them with increasing redshift. Overall, our results indicate that both ρUV and ρ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 reionization.},
  author       = {Santos, S and Sobral, D and Butterworth, J and Paulino-Afonso, A and Ribeiro, B and da Cunha, E and Calhau, J and Khostovan, A A and Matthee, Jorryt J and Arrabal Haro, P},
  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: luminosity function, mass function},
  number       = {1},
  pages        = {1117--1134},
  publisher    = {Oxford University Press},
  title        = {{The evolution of the UV luminosity and stellar mass functions of Lyman-α emitters from z ∼ 2 to z ∼ 6}},
  doi          = {10.1093/mnras/stab1218},
  volume       = {505},
  year         = {2021},
}

@article{11501,
  abstract     = {We investigated the ultraviolet (UV) spectral properties of faint Lyman-α emitters (LAEs) in the redshift range 2.9 ≤ z ≤ 4.6, and we provide material to prepare future observations of the faint Universe. We used data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median MUV of −18 and down to MUV of −16), low stellar mass (median value of 108.4 M⊙ and down to 107 M⊙) LAEs at redshift z ≳ 3. We computed various averaged spectra of LAEs, subsampled on the basis of their observational (e.g., Lyα strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We searched for UV spectral features other than Lyα, such as higher ionization nebular emission lines and absorption features. We successfully observed the O III]λ1666 and [C III]λ1907+C III]λ1909 collisionally excited emission lines and the He IIλ1640 recombination feature, as well as the resonant C IVλλ1548,1551 doublet either in emission or P-Cygni. We compared the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Lyα equivalent widths, fainter UV magnitudes, bluer UV spectral slopes, and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z >  2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass, and faint LAEs. We publicly released the stacked spectra, as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.},
  author       = {Feltre, Anna and Maseda, Michael V. and Bacon, Roland and Pradeep, Jayadev and Leclercq, Floriane and Kusakabe, Haruka and Wisotzki, Lutz and Hashimoto, Takuya and Schmidt, Kasper B. and Blaizot, Jeremy and Brinchmann, Jarle and Boogaard, Leindert and Cantalupo, Sebastiano and Carton, David and Inami, Hanae and Kollatschny, Wolfram and Marino, Raffaella A. and Matthee, Jorryt J and Nanayakkara, Themiya and Richard, Johan and Schaye, Joop and Tresse, Laurence and Urrutia, Tanya and Verhamme, Anne and Weilbacher, Peter M.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution / galaxies: high-redshift / ISM: lines and bands / ultraviolet: ISM / ultraviolet: galaxies},
  publisher    = {EDP Sciences},
  title        = {{The MUSE Hubble Ultra Deep Field Survey: XV. The mean rest-UV spectra of Lyα emitters at z > 3}},
  doi          = {10.1051/0004-6361/202038133},
  volume       = {641},
  year         = {2020},
}

@article{11503,
  abstract     = {Context. The Lyα emitter (LAE) fraction, XLAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of XLAE are still under debate.
Aims. Thanks to deep data obtained with the integral field spectrograph Multi Unit Spectroscopic Explorer (MUSE), we can measure the evolution of the LAE fraction homogeneously over a wide redshift range of z ≈ 3–6 for UV-faint galaxies (down to UV magnitudes of M1500 ≈ −17.75). This is a significantly fainter range than in former studies (M1500 ≤ −18.75) and it allows us to probe the bulk of the population of high-redshift star-forming galaxies.
Methods. We constructed a UV-complete photometric-redshift sample following UV luminosity functions and measured the Lyα emission with MUSE using the latest (second) data release from the MUSE Hubble Ultra Deep Field Survey.
Results. We derived the redshift evolution of XLAE for M1500 ∈ [ − 21.75; −17.75] for the first time with a equivalent width range EW(Lyα) ≥ 65 Å and found low values of XLAE ≲ 30% at z ≲ 6. The best-fit linear relation is XLAE = 0.07+0.06−0.03z − 0.22+0.12−0.24. For M1500 ∈ [ − 20.25; −18.75] and EW(Lyα) ≥ 25 Å, our XLAE values are consistent with those in the literature within 1σ at z ≲ 5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of XLAE on M1500 for EW(Lyα) ≥ 50 Å at z ≈ 3–4, in contrast with previous work. The differences in XLAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lyα emission, hence XLAE is biased towards higher values when those samples are used.
Conclusions. Our results suggest either a lower increase of XLAE towards z ≈ 6 than previously suggested, or even a turnover of XLAE at z ≈ 5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed LAE fractions much better than models where SF is a smooth function of time.},
  author       = {Kusakabe, Haruka and Blaizot, Jérémy and Garel, Thibault and Verhamme, Anne and Bacon, Roland and Richard, Johan and Hashimoto, Takuya and Inami, Hanae and Conseil, Simon and Guiderdoni, Bruno and Drake, Alyssa B. and Christian Herenz, Edmund and Schaye, Joop and Oesch, Pascal and Matthee, Jorryt J and Anna Marino, Raffaella and Borello Schmidt, Kasper and Pelló, Roser and Maseda, Michael and Leclercq, Floriane and Kerutt, Josephine and Mahler, Guillaume},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, dark ages / reionization / first stars / early Universe / cosmology: observations / galaxies: evolution / galaxies: high-redshift / intergalactic medium},
  publisher    = {EDP Sciences},
  title        = {{The MUSE Hubble Ultra Deep Field Survey: XIV. Evolution of the Lyα emitter fraction from z = 3 to z = 6}},
  doi          = {10.1051/0004-6361/201937340},
  volume       = {638},
  year         = {2020},
}

@article{11504,
  abstract     = {We present spatially resolved maps of six individually-detected Lyman α haloes (LAHs) as well as a first statistical analysis of the Lyman α (Lyα) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies (−17.5 >  MUV >  −21.5) using the Multi-Unit Spectroscopic Explorer. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lyα line profile. Using a three-dimensional two-component model for the Lyα emission, we measured the full width at half maximum (FWHM), the peak velocity shift, and the asymmetry of the Lyα line in the core and in the halo of 19 galaxies. We find that the Lyα line shape is statistically different in the halo compared to the core (in terms of width, peak wavelength, and asymmetry) for ≈40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lyα line both at the interstellar and circum-galactic scales. This trend has been predicted by radiative transfer simulations of galactic winds as a result of resonant scattering in outflows. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies (MUV <  −20) show broader, more redshifted, and less asymmetric Lyα lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lyα lines. The generally broad and red line shapes found in the halo component suggest that the Lyα haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lyα line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.},
  author       = {Leclercq, Floriane and Bacon, Roland and Verhamme, Anne and Garel, Thibault and Blaizot, Jérémy and Brinchmann, Jarle and Cantalupo, Sebastiano and Claeyssens, Adélaïde and Conseil, Simon and Contini, Thierry and Hashimoto, Takuya and Herenz, Edmund Christian and Kusakabe, Haruka and Marino, Raffaella Anna and Maseda, Michael and Matthee, Jorryt J and Mitchell, Peter and Pezzulli, Gabriele and Richard, Johan and Schmidt, Kasper Borello and Wisotzki, Lutz},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics galaxies: high-redshift / galaxies: formation / galaxies: evolution / cosmology: observations},
  publisher    = {EDP Sciences},
  title        = {{The MUSE Hubble Ultra Deep field survey: XIII. Spatially resolved spectral properties of Lyman α haloes around star-forming galaxies at z > 3}},
  doi          = {10.1051/0004-6361/201937339},
  volume       = {635},
  year         = {2020},
}

@article{11529,
  abstract     = {CR7 is among the most luminous Ly α emitters (LAEs) known at z = 6.6 and consists of at least three UV components that are surrounded by Ly α emission. Previous studies have suggested that it may host an extreme ionizing source. Here, we present deep integral field spectroscopy of CR7 with VLT/Multi Unit Spectroscopic Explorer (MUSE). We measure extended emission with a similar halo scale length as typical LAEs at z ≈ 5. CR7’s Ly α halo is clearly elongated along the direction connecting the multiple components, likely tracing the underlying gas distribution. The Ly α emission originates almost exclusively from the brightest UV component, but we also identify a faint kinematically distinct Ly α emitting region nearby a fainter component. Combined with new near-infrared data, the MUSE data show that the rest-frame Ly α equivalent width (EW) is ≈100 Å. This is a factor 4 higher than the EW measured in low-redshift analogues with carefully matched Ly α profiles (and thus arguably H I column density), but this EW can plausibly be explained by star formation. Alternative scenarios requiring active galactic nucleus (AGN) powering are also disfavoured by the narrower and steeper Ly α spectrum and much smaller IR to UV ratio compared to obscured AGN in other Ly α blobs. CR7’s Ly α emission, while extremely luminous, resembles the emission in more common LAEs at lower redshifts very well and is likely powered by a young metal-poor starburst.},
  author       = {Matthee, Jorryt J and Pezzulli, Gabriele and Mackenzie, Ruari and Cantalupo, Sebastiano and Kusakabe, Haruka and Leclercq, Floriane and Sobral, David and Richard, Johan and Wisotzki, Lutz and Lilly, Simon and Boogaard, Leindert and Marino, Raffaella and Maseda, Michael and Nanayakkara, Themiya},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: high-redshift, dark ages, reionization, first stars, cosmology: observations},
  number       = {2},
  pages        = {3043--3059},
  publisher    = {Oxford University Press},
  title        = {{The nature of CR7 revealed with MUSE: A young starburst powering extended Ly α emission at z = 6.6}},
  doi          = {10.1093/mnras/staa2550},
  volume       = {498},
  year         = {2020},
}

@article{11531,
  abstract     = {While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z ≈ 4–5 continuum-faint Lyman-α emitters from deep VLT/MUSE spectroscopy and directly measure their H α emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than −16 (≈ 0.01 L⋆), implying a median Lyman-α equivalent width of 259 Å. By combining the H α measurement with the UV magnitude, we determine the ionizing photon production efficiency, ξion, a first for such faint galaxies. The measurement of log10 (ξion [Hz erg−1]) = 26.28 (⁠+0.28−0.40⁠) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower luminosity, Lyman-α-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4 × 10−4 and 0.008, with light-weighted ages less than 3 Myr.},
  author       = {Maseda, Michael V and Bacon, Roland and Lam, Daniel and Matthee, Jorryt J and Brinchmann, Jarle and Schaye, Joop and Labbe, Ivo and Schmidt, Kasper B and Boogaard, Leindert and Bouwens, Rychard and Cantalupo, Sebastiano and Franx, Marijn and Hashimoto, Takuya and Inami, Hanae and Kusakabe, Haruka and Mahler, Guillaume and Nanayakkara, Themiya and Richard, Johan and Wisotzki, Lutz},
  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},
  number       = {4},
  pages        = {5120--5130},
  publisher    = {Oxford University Press},
  title        = {{Elevated ionizing photon production efficiency in faint high-equivalent-width Lyman-α emitters}},
  doi          = {10.1093/mnras/staa622},
  volume       = {493},
  year         = {2020},
}

@article{11533,
  abstract     = {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.) Lyman α (Ly α) emitters (LAEs) at z ∼ 2–6. We find typical stellar masses of 109.3 ± 0.6 M⊙ and star formation rates (SFR) of SFRSED=4.4+10.5−2.4 M⊙ yr−1 and SFRLyα=5.9+6.3−2.6 M⊙ yr−1, combined with very blue UV slopes of β=−2.1+0.5−0.4⁠, but with significant variations within the population. MUV and β 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 Ly α rest-frame equivalent width (EW0) and find 45 ‘extreme’ LAEs with EW0 > 240 Å (3σ), implying a low number density of (7 ± 1) × 10−7 Mpc−3. Overall, we measure little to no evolution of the Ly α EW0 and scale length parameter (w0), which are consistently high (EW0=140+280−70 Å, w0=129+11−11 Å) from z ∼ 6 to z ∼ 2 and below. However, w0 is anticorrelated with MUV and stellar mass. Our results imply that sources selected as LAEs have a high Ly α escape fraction (fesc,Ly α) irrespective of cosmic time, but fesc,Ly α is still higher for UV-fainter and lower mass LAEs. The least massive LAEs (<109.5 M⊙) are typically located above the star formation ‘main sequence’ (MS), but the offset from the MS decreases towards z ∼ 6 and towards 1010 M⊙. 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.},
  author       = {Santos, S and Sobral, D and Matthee, Jorryt J and Calhau, J and da Cunha, E and Ribeiro, B and Paulino-Afonso, A and Arrabal Haro, P and Butterworth, J},
  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: star formation},
  number       = {1},
  pages        = {141--160},
  publisher    = {Oxford University Press},
  title        = {{The evolution of rest-frame UV properties, Ly α EWs, and the SFR–stellar mass relation at z ∼ 2–6 for SC4K LAEs}},
  doi          = {10.1093/mnras/staa093},
  volume       = {493},
  year         = {2020},
}

@article{11534,
  abstract     = {The observed properties of the Lyman-α (Ly α) emission line are a powerful probe of neutral gas in and around galaxies. We present spatially resolved Ly α spectroscopy with VLT/MUSE targeting VR7, a UV-luminous galaxy at z = 6.532 with moderate Ly α equivalent width (EW0 ≈ 38 Å). These data are combined with deep resolved [CII]158μm spectroscopy obtained with ALMA and UV imaging from HST and we also detect UV continuum with MUSE. Ly α emission is clearly detected with S/N ≈ 40 and FWHM of 374 km s−1. Ly α and [C II] are similarly extended beyond the UV, with effective radius reff = 2.1 ± 0.2 kpc for a single exponential model or reff,Lyα,halo=3.45+1.08−0.87 kpc when measured jointly with the UV continuum. The Ly α profile is broader and redshifted with respect to the [C II] line (by 213 km s−1), but there are spatial variations that are qualitatively similar in both lines and coincide with resolved UV components. This suggests that the emission originates from two components with plausibly different H I column densities. We place VR7 in the context of other galaxies at similar and lower redshift. The Ly α halo scale length is similar at different redshifts and velocity shifts with respect to the systemic are typically smaller. Overall, we find little indications of a more neutral vicinity at higher redshift. This means that the local (∼10 kpc) neutral gas conditions that determine the observed Ly α properties in VR7 resemble the conditions in post-reionization galaxies.},
  author       = {Matthee, Jorryt J and Sobral, David and Gronke, Max and Pezzulli, Gabriele and Cantalupo, Sebastiano and Röttgering, Huub and Darvish, Behnam and Santos, Sérgio},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: high-redshift, dark ages, reionization, first stars, cosmology: observations},
  number       = {2},
  pages        = {1778--1790},
  publisher    = {Oxford University Press},
  title        = {{Resolved Lyman-α properties of a luminous Lyman-break galaxy in a large ionized bubble at z = 6.53 }},
  doi          = {10.1093/mnras/stz3554},
  volume       = {492},
  year         = {2020},
}

@article{11539,
  abstract     = {Despite recent progress in understanding Ly α emitters (LAEs), relatively little is known regarding their typical black hole activity across cosmic time. Here, we study the X-ray and radio properties of ∼4000 LAEs at 2.2 < z < 6 from the SC4K survey in the COSMOS field. We detect 254 (⁠6.8per cent±0.4per cent⁠) LAEs individually in the X-rays (S/N > 3) with an average luminosity of 1044.31±0.01ergs−1 and average black hole accretion rate (BHAR) of 0.72±0.01 M⊙ yr−1, consistent with moderate to high accreting active galactic neuclei (AGNs). We detect 120 sources in deep radio data (radio AGN fraction of 3.2per cent±0.3per cent⁠). The global AGN fraction (⁠8.6per cent±0.4per cent⁠) rises with Ly α luminosity and declines with increasing redshift. For X-ray-detected LAEs, Ly α luminosities correlate with the BHARs, suggesting that Ly α luminosity becomes a BHAR indicator. Most LAEs (⁠93.1per cent±0.6per cent⁠) at 2 < z < 6 have no detectable X-ray emission (BHARs < 0.017 M⊙ yr−1). The median star formation rate (SFR) of star-forming LAEs from Ly α and radio luminosities is 7.6+6.6−2.8 M⊙ yr−1. The black hole to galaxy growth ratio (BHAR/SFR) for LAEs is <0.0022, consistent with typical star-forming galaxies and the local BHAR/SFR relation. We conclude that LAEs at 2 < z < 6 include two different populations: an AGN population, where Ly α luminosity traces BHAR, and another with low SFRs which remain undetected in even the deepest X-ray stacks but is detected in the radio stacks.},
  author       = {Calhau, João and Sobral, David and Santos, Sérgio and Matthee, Jorryt J and Paulino-Afonso, Ana and Stroe, Andra and Simmons, Brooke and Barlow-Hall, Cassandra and Adams, Benjamin},
  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, quasars: supermassive black holes, galaxies: star formation, cosmology: observations, X-rays: galaxies},
  number       = {3},
  pages        = {3341--3362},
  publisher    = {Oxford University Press},
  title        = {{The X-ray and radio activity of typical and luminous Ly α emitters from z ∼ 2 to z ∼ 6: Evidence for a diverse, evolving population}},
  doi          = {10.1093/mnras/staa476},
  volume       = {493},
  year         = {2020},
}

@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{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{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},
}