@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{11490,
  abstract     = {Directly characterising the first generations of stars in distant galaxies is a key quest of observational cosmology. We present a case study of ID53 at z = 4.77, the UV-brightest (but L⋆) star-forming galaxy at z > 3 in the MUSE eXtremely Deep Field with a mass of ≈109 M⊙. In addition to very strong Lyman-α (Lyα) emission, we clearly detect the (stellar) continuum and an N V P Cygni feature, interstellar absorption, fine-structure emission and nebular C IV emission lines in the 140 h spectrum. Continuum emission from two spatially resolved components in Hubble Space Telescope data are blended in the MUSE data, but we show that the nebular C IV emission originates from a subcomponent of the galaxy. The UV spectrum can be fit with recent BPASS stellar population models combined with single-burst or continuous star formation histories (SFHs), a standard initial mass function, and an attenuation law. Models with a young age and low metallicity (log10(age/yr) = 6.5–7.6 and [Z/H] = −2.15 to −1.15) are preferred, but the details depend on the assumed SFH. The intrinsic Hα luminosity of the best-fit models is an order of magnitude higher than the Hα luminosity inferred from Spitzer/IRAC data, which either suggests a high escape fraction of ionising photons, a high relative attenuation of nebular to stellar dust, or a complex SFH. The metallicity appears lower than the metallicity in more massive galaxies at z = 3 − 5, consistent with the scenario according to which younger galaxies have lower metallicities. This chemical immaturity likely facilitates Lyα escape, explaining why the Lyα equivalent width is anti-correlated with stellar metallicity. Finally, we stress that uncertainties in SFHs impose a challenge for future inferences of the stellar metallicity of young galaxies. This highlights the need for joint (spatially resolved) analyses of stellar spectra and photo-ionisation models.},
  author       = {Matthee, Jorryt J and Feltre, Anna and Maseda, Michael and Nanayakkara, Themiya and Boogaard, Leindert and Bacon, Roland and Verhamme, Anne and Leclercq, Floriane and Kusakabe, Haruka and Urrutia, Tanya and Wisotzki, Lutz},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / techniques: spectroscopic / galaxies: stellar content / galaxies: formation},
  publisher    = {EDP Sciences},
  title        = {{Deciphering stellar metallicities in the early universe: Case study of a young galaxy at z = 4.77 in the MUSE eXtremely Deep Field}},
  doi          = {10.1051/0004-6361/202142187},
  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{11522,
  abstract     = {The decline in abundance of Lyman-α (Lyα) emitting galaxies at z ≳ 6 is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CODAII simulation, which is a radiation hydrodynamic simulation featuring a box of ∼94 comoving Mpc side length, to compute the Lyα transmission properties of the intergalactic medium (IGM) at z ∼ 5.8 to 7. Our results mainly confirm previous studies, i.e. we find a declining Lyα transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Lyα peaks at high redshift, we also analyse the IGM transmission on the blue side, which shows a rapid decline at z ≳ 6 of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as nHI∼10−9cm−3 (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that ∼1 per cent of sightlines towards M1600AB ∼ −21 galaxies at z ∼ 7 are transparent enough to allow a transmission of a blue Lyα peak. We discuss our results in the context of the interpretation of observations.},
  author       = {Gronke, Max and Ocvirk, Pierre and Mason, Charlotte and Matthee, Jorryt J and Bosman, Sarah E I and Sorce, Jenny G and Lewis, Joseph and Ahn, Kyungjin and Aubert, Dominique and Dawoodbhoy, Taha and Iliev, Ilian T and Shapiro, Paul R and Yepes, Gustavo},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {dark ages, reionization, first stars, intergalactic medium, galaxies: formation},
  number       = {3},
  pages        = {3697--3709},
  publisher    = {Oxford University Press},
  title        = {{Lyman-α transmission properties of the intergalactic medium in the CoDaII simulation}},
  doi          = {10.1093/mnras/stab2762},
  volume       = {508},
  year         = {2021},
}

@article{11523,
  abstract     = {We present the first results from the X-SHOOTER Lyman α survey at z = 2 (XLS-z2). XLS-z2 is a deep spectroscopic survey of 35 Lyman α emitters (LAEs) utilizing ≈90 h of exposure time with Very Large Telescope/X-SHOOTER and covers rest-frame Ly α to H α emission with R ≈ 4000. We present the sample selection, the observations, and the data reduction. Systemic redshifts are measured from rest-frame optical lines for 33/35 sources. In the stacked spectrum, our LAEs are characterized by an interstellar medium with little dust, a low metallicity, and a high ionization state. The ionizing sources are young hot stars that power strong emission lines in the optical and high-ionization lines in the ultraviolet (UV). The LAEs exhibit clumpy UV morphologies and have outflowing kinematics with blueshifted Si II absorption, a broad [O III] component, and a red-skewed Ly α line. Typically, 30 per cent of the Ly α photons escape, of which one quarter on the blue side of the systemic velocity. A fraction of Ly α photons escape directly at the systemic suggesting clear channels enabling an ≈10 per cent escape of ionizing photons, consistent with an inference based on Mg II. A combination of a low effective H I column density, a low dust content, and young starburst determines whether a star-forming galaxy is observed as an LAE. The first is possibly related to outflows and/or a fortunate viewing angle, while we find that the latter two in LAEs are typical for their stellar mass of 109 M⊙.},
  author       = {Matthee, Jorryt J and Sobral, David and Hayes, Matthew and Pezzulli, Gabriele and Gronke, Max and Schaerer, Daniel and Naidu, Rohan P and Röttgering, Huub and Calhau, João and Paulino-Afonso, Ana and Santos, Sérgio and Amorín, Ricardo},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: formation, galaxies: ISM, galaxies: starburst, dark ages, reionization, first stars},
  number       = {1},
  pages        = {1382--1412},
  publisher    = {Oxford University Press},
  title        = {{The X-SHOOTER Lyman α survey at z = 2 (XLS-z2) I: What makes a galaxy a Lyman α emitter?}},
  doi          = {10.1093/mnras/stab1304},
  volume       = {505},
  year         = {2021},
}

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

@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{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{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{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{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{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{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{11575,
  abstract     = {We investigate correlations between different physical properties of star-forming galaxies in the ‘Evolution and Assembly of GaLaxies and their Environments’ (EAGLE) cosmological hydrodynamical simulation suite over the redshift range 0 ≤ z ≤ 4.5. A principal component analysis reveals that neutral gas fraction (fgas,neutral), stellar mass (Mstellar) and star formation rate (SFR) account for most of the variance seen in the population, with galaxies tracing a two-dimensional, nearly flat, surface in the three-dimensional space of fgas, neutral–Mstellar–SFR with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their fgas, neutral and SFR drop with redshift. The positions of galaxies along the plane are highly correlated with gas metallicity. The metallicity can therefore be robustly predicted from fgas, neutral, or from the Mstellar and SFR. We argue that the appearance of this ‘Fundamental Plane of star formation’ is a consequence of self-regulation, with the plane's curvature set by the dependence of the SFR on gas density and metallicity. We analyse a large compilation of observations spanning the redshift range 0 ≲ z ≲ 3, and find that such a plane is also present in the data. The properties of the observed Fundamental Plane of star formation are in good agreement with EAGLE's predictions.},
  author       = {Lagos, Claudia del P. and Theuns, Tom and Schaye, Joop and Furlong, Michelle and Bower, Richard G. and Schaller, Matthieu and Crain, Robert A. and Trayford, James W. and Matthee, Jorryt J},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics  stars: formation, ISM: evolution, galaxies: evolution, galaxies: formation, galaxies: ISM},
  number       = {3},
  pages        = {2632--2650},
  publisher    = {Oxford University Press},
  title        = {{The Fundamental Plane of star formation in galaxies revealed by the EAGLE hydrodynamical simulations}},
  doi          = {10.1093/mnras/stw717},
  volume       = {459},
  year         = {2016},
}

@article{11580,
  abstract     = {We present results from the largest contiguous narrow-band survey in the near-infrared. We have used the wide-field infrared camera/Canada–France–Hawaii Telescope and the lowOH2 filter (1.187 ± 0.005 μm) to survey ≈10 deg2 of contiguous extragalactic sky in the SA22 field. A total of ∼6000 candidate emission-line galaxies are found. We use deep ugrizJK data to obtain robust photometric redshifts. We combine our data with the High-redshift(Z) Emission Line Survey (HiZELS), explore spectroscopic surveys (VVDS, VIPERS) and obtain our own spectroscopic follow-up with KMOS, FMOS and MOSFIRE to derive large samples of high-redshift emission-line selected galaxies: 3471 Hα emitters at z = 0.8, 1343 [O III] + Hβ emitters at z = 1.4 and 572 [O II] emitters at z = 2.2. We probe comoving volumes of >106 Mpc3 and find significant overdensities, including an 8.5σ (spectroscopically confirmed) overdensity of Hα emitters at z = 0.81. We derive Hα, [O III] + Hβ and [O II] luminosity functions at z = 0.8, 1.4, 2.2, respectively, and present implications for future surveys such as Euclid. Our uniquely large volumes/areas allow us to subdivide the samples in thousands of randomized combinations of areas and provide a robust empirical measurement of sample/cosmic variance. We show that surveys for star-forming/emission-line galaxies at a depth similar to ours can only overcome cosmic-variance (errors <10 per cent) if they are based on volumes >5 × 105 Mpc3; errors on L* and ϕ* due to sample (cosmic) variance on surveys probing ∼104 and ∼105 Mpc3 are typically very high: ∼300 and ∼40–60 per cent, respectively.},
  author       = {Sobral, D. and Matthee, Jorryt J and Best, P. N. and Smail, I. and Khostovan, A. A. and Milvang-Jensen, B. and Kim, J.-W. and Stott, J. and Calhau, J. and Nayyeri, H. and Mobasher, B.},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: formation, galaxies: luminosity function, mass function, cosmology: observations, early Universe, large-scale structure of Universe},
  number       = {3},
  pages        = {2303--2323},
  publisher    = {Oxford University Press},
  title        = {{CF-HiZELS, an ∼10 deg2 emission-line survey with spectroscopic follow-up: Hα, [O III] + Hβ and [O II] luminosity functions at z = 0.8, 1.4 and 2.2 }},
  doi          = {10.1093/mnras/stv1076},
  volume       = {451},
  year         = {2015},
}