@article{14852, abstract = {The physical conditions giving rise to high escape fractions of ionizing radiation (LyC fesc) in star-forming galaxies – most likely protagonists of cosmic reionization – are not yet fully understood. Using the VLT/MUSE observations of ∼1400 Ly α emitters at 2.9 &lt; z &lt; 6.7, we compare stacked rest-frame UV spectra of candidates for LyC leakers and non-leakers selected based on their Ly α profiles. We find that the stacks of potential LyC leakers, i.e. galaxies with narrow, symmetric Ly α profiles with small peak separation, generally show (i) strong nebular O iii]λ1666, [Si iii]λ1883, and [C iii]λ1907 +C iii]λ1909 emission, indicating a high-ionization state of the interstellar medium (ISM); (ii) high equivalent widths of He iiλ1640 (∼1 − 3 Å), suggesting the presence of hard ionizing radiation fields; (iii) Si ii*λ1533 emission, revealing substantial amounts of neutral hydrogen off the line of sight; (iv) high C ivλλ1548,1550 to [C iii]λ1907 +C iii]λ1909 ratios (C iv/C iii] ≳0.75) , signalling the presence of low column density channels in the ISM. In contrast, the stacks with broad, asymmetric Ly α profiles with large peak separation show weak nebular emission lines, low He iiλ1640 equivalent widths (≲1 Å), and low C iv/C iii] (≲0.25), implying low-ionization states and high-neutral hydrogen column densities. Our results suggest that C iv/C iii] might be sensitive to the physical conditions that govern LyC photon escape, providing a promising tool for identification of ionizing sources among star-forming galaxies in the epoch of reionization.}, author = {Kramarenko, Ivan and Kerutt, J and Verhamme, A and Oesch, P A and Barrufet, L and Matthee, Jorryt J and Kusakabe, H and Goovaerts, I and Thai, T T}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {4}, pages = {9853--9871}, publisher = {Oxford University Press}, title = {{Linking UV spectral properties of MUSE Ly α emitters at z ≳ 3 to Lyman continuum escape}}, doi = {10.1093/mnras/stad3853}, volume = {527}, year = {2024}, } @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α haloes around rest-frame UV-selected galaxies at z ≃ 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{11509, abstract = {We present Keck/LRIS follow-up spectroscopy for 13 photometric candidates of extremely metal-poor galaxies (EMPGs) selected by a machine-learning technique applied to the deep (∼26 AB mag) optical and wide-area (∼500 deg2) Subaru imaging data in the EMPRESS survey. Nine out of the 13 candidates are EMPGs with an oxygen abundance (O/H) less than ∼10% solar value (O/H)⊙, and four sources are contaminants of moderately metal-rich galaxies or no emission-line objects. Notably, two out of the nine EMPGs have extremely low stellar masses and oxygen abundances of 5 × 10⁴x–7 × -10⁵ M⊙ and 2%–3% (O/H)⊙, respectively. With a sample of five EMPGs with (Fe/O) measurements, two (three) of which are taken from this study (the literature), we confirm that two EMPGs with the lowest (O/H) ratios of ∼2% (O/H)⊙ show high (Fe/O) ratios of ∼0.1, close to the solar abundance ratio. Comparing galaxy chemical enrichment models, we find that the two EMPGs cannot be explained by a scenario of metal-poor gas accretion/episodic star formation history due to their low (N/O) ratios. We conclude that the two EMPGs can be reproduced by the inclusion of bright hypernovae and/or hypothetical pair-instability supernovae (SNe) preferentially produced in a metal-poor environment. This conclusion implies that primordial galaxies at z ∼ 10 could have a high abundance of Fe that did not originate from Type Ia SNe with delays and that Fe may not serve as a cosmic clock for primordial galaxies.}, author = {Isobe, Yuki and Ouchi, Masami and Suzuki, Akihiro and Moriya, Takashi J. and Nakajima, Kimihiko and Nomoto, Ken’ichi and Rauch, Michael and Harikane, Yuichi and Kojima, Takashi and Ono, Yoshiaki and Fujimoto, Seiji and Inoue, Akio K. and Kim, Ji Hoon and Komiyama, Yutaka and Kusakabe, Haruka and Lee, Chien-Hsiu and Maseda, Michael and Matthee, Jorryt J and Michel-Dansac, Leo and Nagao, Tohru and Nanayakkara, Themiya and Nishigaki, Moka and Onodera, Masato and Sugahara, Yuma and Xu, Yi}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {2}, publisher = {IOP Publishing}, title = {{EMPRESS. IV. Extremely metal-poor galaxies including very low-mass primordial systems with M∗= 10⁴-10⁵⊙ and 2%–3% (O/H): High (Fe/O) suggestive of metal enrichment by hypernovae/pair-instability supernovae}}, doi = {10.3847/1538-4357/ac3509}, volume = {925}, year = {2022}, } @article{11510, abstract = {Galaxy evolution is driven by a variety of physical processes that are predicted to proceed at different rates for different dark matter haloes and environments across cosmic times. A record of this evolution is preserved in galaxy stellar populations, which we can access using absorption-line spectroscopy. Here we explore the large LEGA-C survey (DR3) to investigate the role of the environment and stellar mass on stellar populations at z ∼ 0.6–1 in the COSMOS field. Leveraging the statistical power and depth of LEGA-C, we reveal significant gradients in Dn4000 and Hδ equivalent widths (EWs) distributions over the stellar mass versus environment 2D spaces for the massive galaxy population (M > 1010 M⊙) at z ∼ 0.6–1.0. Dn4000 and Hδ EWs primarily depend on stellar mass, but they also depend on environment at fixed stellar mass. By splitting the sample into centrals and satellites, and in terms of star-forming galaxies and quiescent galaxies, we reveal that the significant environmental trends of Dn4000 and Hδ EW, when controlling for stellar mass, are driven by quiescent galaxies. Regardless of being centrals or satellites, star-forming galaxies reveal Dn4000 and Hδ EWs, which depend strongly on their stellar mass and are completely independent of the environment at 0.6 < z < 1.0. The environmental trends seen for satellite galaxies are fully driven by the trends that hold only for quiescent galaxies, combined with the strong environmental dependency of the quiescent fraction at fixed stellar mass. Our results are consistent with recent predictions from simulations that point toward massive galaxies forming first in overdensities or the most compact dark matter haloes.}, author = {Sobral, David and van der Wel, Arjen and Bezanson, Rachel and Bell, Eric and Muzzin, Adam and D’Eugenio, Francesco and Darvish, Behnam and Gallazzi, Anna and Wu, Po-Feng and Maseda, Michael and Matthee, Jorryt J and Paulino-Afonso, Ana and Straatman, Caroline and van Dokkum, Pieter G.}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {2}, publisher = {IOP Publishing}, title = {{The LEGA-C of nature and nurture in stellar populations at z ∼ 0.6–1.0: Dn4000 and Hδ reveal different assembly histories for quiescent galaxies in different environments}}, doi = {10.3847/1538-4357/ac4419}, volume = {926}, year = {2022}, } @article{11511, abstract = {The ratio of α-elements to iron in galaxies holds valuable information about the star formation history (SFH) since their enrichment occurs on different timescales. The fossil record of stars in galaxies has mostly been excavated for passive galaxies, since the light of star-forming galaxies is dominated by young stars, which have much weaker atmospheric absorption features. Here we use the largest reference cosmological simulation of the EAGLE project to investigate the origin of variations in stellar α-enhancement among star-forming galaxies at z = 0, and their impact on integrated spectra. The definition of α-enhancement in a composite stellar population is ambiguous. We elucidate two definitions—termed “mean” and “galactic” α-enhancement—in more detail. While a star-forming galaxy has a high “mean” α-enhancement when its stars formed rapidly, a galaxy with a large “galactic” α-enhancement generally had a delayed SFH. We find that absorption-line strengths of Mg and Fe correlate with variations in α-enhancement. These correlations are strongest for the “galactic” α-enhancement. However, we show that these are mostly caused by other effects that are cross-correlated with α-enhancement, such as variations in the light-weighted age. This severely complicates the retrieval of α-enhancements in star-forming galaxies. The ambiguity is not severe for passive galaxies, and we confirm that spectral variations in these galaxies are caused by measurable variations in α-enhancements. We suggest that this more complex coupling between α-enhancement and SFHs can guide the interpretation of new observations of star-forming galaxies.}, author = {Gebek, Andrea and Matthee, Jorryt J}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {2}, publisher = {IOP Publishing}, title = {{On the variation in stellar α-enhancements of star-forming galaxies in the EAGLE simulation}}, doi = {10.3847/1538-4357/ac350b}, volume = {924}, year = {2022}, } @article{11521, abstract = {The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here, we introduce an emissivity framework based on Lyα emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between production and escape of Lyα and LyC photons. Using constraints on the escape fractions of bright LAEs (LLyα > 0.2L*) at z ≈ 2 obtained from resolved Lyα profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at z ≈ 2–6, which is non-trivial given the strong evolution in both the star formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between z ≈ 6−9 under plausible assumptions. Within this framework, the >10 × rise in the UV population-averaged fesc between z ≈ 3–7 naturally arises due to the same phenomena that drive the growing LAE fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters (⁠fesc=50 per cent⁠, ξion = 1025.9 Hz erg−1 for half the bright LAEs), a highly ionizing minority of galaxies with MUV < −17 accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a ‘disco’ ionizing background. We conclude proposing tests to further develop our suggested Lyα-anchored formalism.}, author = {Matthee, Jorryt J and Naidu, Rohan P. and Pezzulli, Gabriele and Gronke, Max and Sobral, David and Oesch, Pascal A. and Hayes, Matthew and Erb, Dawn and Schaerer, Daniel and Amorín, Ricardo and Tacchella, Sandro and Ana Paulino-Afonso, Ana Paulino-Afonso and Llerena, Mario and Calhau, João and Röttgering, Huub}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {galaxies: high-redshift, intergalactic medium, cosmology: observations, dark ages, reionization, first stars, ultraviolet: galaxies}, number = {4}, pages = {5960--5977}, publisher = {Oxford University Press}, title = {{(Re)Solving reionization with Lyα: How bright Lyα emitters account for the z ≈ 2 − 8 cosmic ionizing background}}, doi = {10.1093/mnras/stac801}, volume = {512}, year = {2022}, } @article{11498, abstract = {Rest-frame ultraviolet (UV) emission lines probe electron densities, gas-phase abundances, metallicities, and ionization parameters of the emitting star-forming galaxies and their environments. The strongest main UV emission line, Lyα, has been instrumental in advancing the general knowledge of galaxy formation in the early universe. However, observing Lyα emission becomes increasingly challenging at z ≳ 6 when the neutral hydrogen fraction of the circumgalactic and intergalactic media increases. Secondary weaker UV emission lines provide important alternative methods for studying galaxy properties at high redshift. We present a large sample of rest-frame UV emission line sources at intermediate redshift for calibrating and exploring the connection between secondary UV lines and the emitting galaxies’ physical properties and their Lyα emission. The sample of 2052 emission line sources with 1.5 < z < 6.4 was collected from integral field data from the MUSE-Wide and MUSE-Deep surveys taken as part of Guaranteed Time Observations. The objects were selected through untargeted source detection (i.e., no preselection of sources as in dedicated spectroscopic campaigns) in the three-dimensional MUSE data cubes. We searched optimally extracted one-dimensional spectra of the full sample for UV emission features via emission line template matching, resulting in a sample of more than 100 rest-frame UV emission line detections. We show that the detection efficiency of (non-Lyα) UV emission lines increases with survey depth, and that the emission line strength of He IIλ1640 Å, [O III] λ1661 + O III] λ1666, and [Si III] λ1883 + Si III] λ1892 correlate with the strength of [C III] λ1907 + C III] λ1909. The rest-frame equivalent width (EW0) of [C III] λ1907 + C III] λ1909 is found to be roughly 0.22 ± 0.18 of EW0(Lyα). We measured the velocity offsets of resonant emission lines with respect to systemic tracers. For C IVλ1548 + C IVλ1551 we find that ΔvC IV ≲ 250 km s−1, whereas ΔvLyα falls in the range of 250−500 km s−1 which is in agreement with previous results from the literature. The electron density ne measured from [Si III] λ1883 + Si III] λ1892 and [C III] λ1907 + C III] λ1909 line flux ratios is generally < 105 cm−3 and the gas-phase abundance is below solar at 12 + log10(O/H)≈8. Lastly, we used “PhotoIonization Model Probability Density Functions” to infer physical parameters of the full sample and individual systems based on photoionization model parameter grids and observational constraints from our UV emission line searches. This reveals that the UV line emitters generally have ionization parameter log10(U) ≈ −2.5 and metal mass fractions that scatter around Z ≈ 10−2, that is Z ≈ 0.66 Z⊙. Value-added catalogs of the full sample of MUSE objects studied in this work and a collection of UV line emitters from the literature are provided with this paper.}, author = {Schmidt, K. B. and Kerutt, J. and Wisotzki, L. and Urrutia, T. and Feltre, A. and Maseda, M. V. and Nanayakkara, T. and Bacon, R. and Boogaard, L. A. and Conseil, S. and Contini, T. and Herenz, E. C. and Kollatschny, W. and Krumpe, M. and Leclercq, F. and Mahler, G. and Matthee, Jorryt J and Mauerhofer, V. and Richard, J. and Schaye, J.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, ultraviolet: galaxies / galaxies: high-redshift / galaxies: ISM / ISM: lines and bands / methods: observational / techniques: imaging spectroscopy}, publisher = {EDP Sciences}, title = {{Recovery and analysis of rest-frame UV emission lines in 2052 galaxies observed with MUSE at 1.5 < z < 6.4}}, doi = {10.1051/0004-6361/202140876}, volume = {654}, year = {2021}, } @article{11500, abstract = {We report the discovery of diffuse extended Lyα emission from redshift 3.1 to 4.5, tracing cosmic web filaments on scales of 2.5−4 cMpc. These structures have been observed in overdensities of Lyα emitters in the MUSE Extremely Deep Field, a 140 h deep MUSE observation located in the Hubble Ultra-Deep Field. Among the 22 overdense regions identified, five are likely to harbor very extended Lyα emission at high significance with an average surface brightness of 5 × 10−20 erg s−1 cm−2 arcsec−2. Remarkably, 70% of the total Lyα luminosity from these filaments comes from beyond the circumgalactic medium of any identified Lyα emitter. Fluorescent Lyα emission powered by the cosmic UV background can only account for less than 34% of this emission at z ≈ 3 and for not more than 10% at higher redshift. We find that the bulk of this diffuse emission can be reproduced by the unresolved Lyα emission of a large population of ultra low-luminosity Lyα emitters (< 1040 erg s−1), provided that the faint end of the Lyα luminosity function is steep (α ⪅ −1.8), it extends down to luminosities lower than 1038 − 1037 erg s−1, and the clustering of these Lyα emitters is significant (filling factor < 1/6). If these Lyα emitters are powered by star formation, then this implies their luminosity function needs to extend down to star formation rates < 10−4 M⊙ yr−1. These observations provide the first detection of the cosmic web in Lyα emission in typical filamentary environments and the first observational clue indicating the existence of a large population of ultra low-luminosity Lyα emitters at high redshift.}, author = {Bacon, R. and Mary, D. and Garel, T. and Blaizot, J. and Maseda, M. and Schaye, J. and Wisotzki, L. and Conseil, S. and Brinchmann, J. and Leclercq, F. and Abril-Melgarejo, V. and Boogaard, L. and Bouché, N. F. and Contini, T. and Feltre, A. and Guiderdoni, B. and Herenz, C. and Kollatschny, W. and Kusakabe, H. and Matthee, Jorryt J and Michel-Dansac, L. and Nanayakkara, T. and Richard, J. and Roth, M. and Schmidt, K. B. and Steinmetz, M. and Tresse, L. and Urrutia, T. and Verhamme, A. and Weilbacher, P. M. and Zabl, J. and Zoutendijk, S. L.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / galaxies: groups: general / cosmology: observations}, publisher = {EDP Sciences}, title = {{The MUSE Extremely Deep Field: The cosmic web in emission at high redshift}}, doi = {10.1051/0004-6361/202039887}, volume = {647}, year = {2021}, } @article{11512, abstract = {We study the molecular gas content of 24 star-forming galaxies at z = 3–4, with a median stellar mass of 109.1 M⊙, from the MUSE Hubble Ultra Deep Field (HUDF) Survey. Selected by their Lyα λ1216 emission and HF160W-band magnitude, the galaxies show an average $\langle {\mathrm{EW}}_{\mathrm{Ly}\alpha }^{0}\rangle \approx 20$ Å, below the typical selection threshold for Lyα emitters (${\mathrm{EW}}_{\mathrm{Ly}\alpha }^{0}\gt 25$ Å), and a rest-frame UV spectrum similar to Lyman-break galaxies. We use rest-frame optical spectroscopy from KMOS and MOSFIRE, and the UV features observed with MUSE, to determine the systemic redshifts, which are offset from Lyα by 〈Δv(Lyα)〉 = 346 km s−1, with a 100 to 600 km s−1 range. Stacking 12CO J = 4 → 3 and [C i]3P1 → 3P0 (and higher-J CO lines) from the ALMA Spectroscopic Survey of the HUDF, we determine 3σ upper limits on the line luminosities of 4.0 × 108 K km s−1pc2 and 5.6 × 108 K km s−1pc2, respectively (for a 300 km s−1 line width). Stacking the 1.2 mm and 3 mm dust-continuum flux densities, we find a 3σ upper limits of 9 μJy and 1.2 μJy, respectively. The inferred gas fractions, under the assumption of a "Galactic" CO-to-H2 conversion factor and gas-to-dust ratio, are in tension with previously determined scaling relations. This implies a substantially higher αCO ≥ 10 and δGDR ≥ 1200, consistent with the subsolar metallicity estimated for these galaxies ($12+\mathrm{log}({\rm{O}}/{\rm{H}})\approx 7.8\pm 0.2$). The low metallicity of z ≥ 3 star-forming galaxies may thus make it very challenging to unveil their cold gas through CO or dust emission, warranting further exploration of alternative tracers, such as [C ii].}, author = {Boogaard, Leindert A. and Bouwens, Rychard J. and Riechers, Dominik and van der Werf, Paul and Bacon, Roland and Matthee, Jorryt J and Stefanon, Mauro and Feltre, Anna and Maseda, Michael and Inami, Hanae and Aravena, Manuel and Brinchmann, Jarle and Carilli, Chris and Contini, Thierry and Decarli, Roberto and González-López, Jorge and Nanayakkara, Themiya and Walter, Fabian}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {1}, publisher = {IOP Publishing}, title = {{Measuring the average molecular gas content of star-forming galaxies at z = 3–4}}, doi = {10.3847/1538-4357/ac01d7}, volume = {916}, year = {2021}, } @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{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{11525, abstract = {The intensity of the Cosmic UV background (UVB), coming from all sources of ionizing photons such as star-forming galaxies and quasars, determines the thermal evolution and ionization state of the intergalactic medium (IGM) and is, therefore, a critical ingredient for models of cosmic structure formation. Most of the previous estimates are based on the comparison between observed and simulated Lyman-α forest. We present the results of an independent method to constrain the product of the UVB photoionization rate and the covering fraction of Lyman limit systems (LLSs) by searching for the fluorescent Lyman-α emission produced by self-shielded clouds. Because the expected surface brightness is well below current sensitivity limits for direct imaging, we developed a new method based on 3D stacking of the IGM around Lyman-α emitting galaxies (LAEs) between 2.9 < z < 6.6 using deep MUSE observations. Combining our results with covering fractions of LLSs obtained from mock cubes extracted from the EAGLE simulation, we obtain new and independent constraints on the UVB at z > 3 that are consistent with previous measurements, with a preference for relatively low UVB intensities at z = 3, and which suggest a non-monotonic decrease of ΓH I with increasing redshift between 3 < z < 5. This could suggest a possible tension between some UVB models and current observations which however require deeper and wider observations in Lyman-α emission and absorption to be confirmed. Assuming instead a value of UVB from current models, our results constrain the covering fraction of LLSs at 3 < z < 4.5 to be less than 25 per cent within 150 kpc from LAEs.}, author = {Gallego, Sofia G and Cantalupo, Sebastiano and Sarpas, Saeed and Duboeuf, Bastien and Lilly, Simon and Pezzulli, Gabriele and Marino, Raffaella Anna and Matthee, Jorryt J and Wisotzki, Lutz and Schaye, Joop and Richard, Johan and Kusakabe, Haruka and Mauerhofer, Valentin}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {1}, pages = {16--32}, publisher = {Oxford University Press}, title = {{Constraining the cosmic UV background at z > 3 with MUSE Lyman-α emission observations}}, doi = {10.1093/mnras/stab796}, volume = {504}, year = {2021}, } @article{11526, abstract = {We present the results from a MUSE survey of twelve z ≃ 3.15 quasars, which were selected to be much fainter (20 < iSDSS < 23) than in previous studies of giant Ly α nebulae around the brightest quasars (16.6 < iAB < 18.7). We detect H I Ly α nebulae around 100 per cent of our target quasars, with emission extending to scales of at least 60 physical kpc, and up to 190 pkpc. We explore correlations between properties of the nebulae and their host quasars, with the goal of connecting variations in the properties of the illuminating QSO to the response in nebular emission. We show that the surface brightness profiles of the nebulae are similar to those of nebulae around bright quasars, but with a lower normalization. Our targeted quasars are on average 3.7 mag (≃30 times) fainter in UV continuum than our bright reference sample, and yet the nebulae around them are only 4.3 times fainter in mean Ly α surface brightness, measured between 20 and 50 pkpc. We find significant correlations between the surface brightness of the nebula and the luminosity of the quasar in both UV continuum and Ly α. The latter can be interpreted as evidence for a substantial contribution from unresolved inner parts of the nebulae to the narrow components seen in the Ly α lines of some of our faint quasars, possibly from the inner circumgalactic medium or from the host galaxy’s interstellar medium.}, author = {Mackenzie, Ruari and Pezzulli, Gabriele and Cantalupo, Sebastiano and Marino, Raffaella A and Lilly, Simon and Muzahid, Sowgat and Matthee, Jorryt J and Schaye, Joop and Wisotzki, Lutz}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, techniques: imaging spectroscopy, intergalactic medium, quasars: emission lines, quasars: general}, number = {1}, pages = {494--509}, publisher = {Oxford University Press}, title = {{Revealing the impact of quasar luminosity on giant Lyα nebulae}}, doi = {10.1093/mnras/staa3277}, volume = {502}, year = {2021}, } @article{11585, abstract = {Observations show that star-forming galaxies reside on a tight three-dimensional 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 three-dimensional 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 toward 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}, issn = {2397-3366}, journal = {Nature Astronomy}, keywords = {Astronomy and Astrophysics, galaxies, formation - galaxies, evolution - galaxies, star formation - galaxies, abundances}, pages = {984--985}, publisher = {Springer Nature}, title = {{Differences in galaxy colours are not just about the mass}}, doi = {10.1038/s41550-021-01415-y}, volume = {5}, 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}, }