## Charting the Lymanα escape fraction in the range 2 9 . &lt; z &lt; 6 7 . and consequences for the LAE reionisation contribution

I. Goovaerts 1 , T. T. Thai 1 2 3 , , , R. Pello 1 , P. Tuan-Anh 2 3 , , N. Laporte 1 , J. Matthee , T. Nanayakkara 5 4 , and J. Pharo 6

- 1 Aix Marseille Université, CNRS, CNES, LAM (Laboratoire d'Astrophysique de Marseille), UMR 7326, 13388 Marseille, France e-mail: ilias.goovaerts@lam.fr; tranthai.physics@gmail.com
- 2 Department of Astrophysics, Vietnam National Space Center, Vietnam Academy of Science and Technology, 18 Hoang Quoc Viet, Hanoi, Vietnam
- 3 Graduate University of Science and Technology, VAST, 18 Hoang Quoc Viet, Cau Giay, Vietnam
- 4 Institute of Science and Technology Austria (ISTA), Am Campus 1, 3400 Klosterneuburg, Austria
- 5 Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Hawthorn, VIC 3122, Australia
- 6 Leibniz Institut für Astrophysik (AIP), An der Sternwarte 16, 14482 Potsdam, Germany

Received 9th July 2024; accepted 13th July 2024

## ABSTRACT

Context. The escape of Lymanα photons at redshifts greater than two is an ongoing subject of study and an important quantity to further understanding of Lymanα emitters (LAEs), the transmission of Lymanα photons through the interstellar medium and intergalactic medium, and the impact these LAEs have on cosmic reionisation.

Aims. This study aims to assess the Lymanα escape fraction, fesc Ly , α , over the redshift range 2 9 . &lt; z &lt; 6 7, focusing on Very Large . Telescope Multi Unit Spectroscopic Explorer (VLT MUSE) selected, gravitationally lensed, intrinsically faint LAEs. These galaxies / / are of particular interest as the potential drivers of cosmic reionisation.

Methods. We assessed fesc Ly , α in two ways: through an individual study of 96 LAEs behind the A2744 lensing cluster, with James Webb Space Telescope Near-Infrared Camera (JWST NIRCam) and HST data, and through a study of the global evolution of fesc Ly / / , α using the state-of-the-art luminosity functions for LAEs and the UV-selected 'parent' population (dust-corrected). We compared these studies to those in the literature based on brighter samples.

Results. We find a negligible redshift evolution of fesc Ly , α for our individual galaxies; it is likely that it was washed out by significant intrinsic scatter. We observed a more significant evolution towards higher escape fractions with decreasing UV magnitude and fit this relation. When comparing the two luminosity functions to derive fesc Ly , α in a global sense, we saw agreement with previous literature when integrating the luminosity functions to a bright limit. However, when integrating using a faint limit equivalent to the observational limits of our samples, we observed enhanced values of fesc Ly , α , particularly around z ∼ 6, where fesc Ly , α becomes consistent with 100% escape. This indicates for the faint regimes we sampled that galaxies towards reionisation tend to allow very large fractions of Lymanα photons to escape. We interpret this as evidence of a lack of any significant dust in these populations; our sample is likely dominated by young, highly star-forming chemically unevolved galaxies. Finally, we assessed the contribution of the LAE population to reionisation using our latest values for fesc Ly , α and the LAE luminosity density. The dependence on the escape fraction of Lyman continuum photons is strong, but for values similar to those observed recently in z ∼ 3 LAEs and high-redshift analogues, LAEs could provide all the ionising emissivity necessary for reionisation.

Key words. galaxies: high redshift - galaxies: evolution - gravitational lensing: strong

## 1. Introduction

For more than 20 years, the study of high-redshift galaxies has been facilitated by the Lymanα line (Cowie &amp; Hu 1998; Hu et al. 1998). Galaxies for which this line is observable, Lymanα emitters (LAEs; equivalent width, EWLy α &gt; 25Å), have formed an essential part of the drive to characterise star-forming galaxies (SFGs) at increasingly higher redshifts (the most recent examples include, Witstok et al. 2024b; Iani et al. 2023; Chen et al. 2024, see also Ferrara 2024 and references therein). The Lymanα line, while intrinsically very strong and therefore an attractive target for observations, undergoes a very complicated process escaping a galaxy due to the resonant nature of the line (Verhamme et al. 2008; Hayes et al. 2010; Zheng et al. 2010; Dijkstra &amp; Kramer 2012; Matthee et al. 2016b) as well as dust within the galaxy (Atek et al. 2008; Finkelstein et al. 2009; Hayes et al. 2013a), inflows and outflows (Hansen &amp; Oh 2006;

Gronke &amp; Dijkstra 2016), and the intergalactic medium (IGM) (Stark et al. 2010, 2011). Understanding the escape of these Lymanα photons is of paramount importance to the characterisation of these high-redshift SFGs as well as to understanding their impact on extragalactic processes such as cosmic reionisation.

The escape of Lymanα photons is also expected to be connected to the escape of ionising Lyman continuum photons ( λ LyC &lt; 912 Å) (Dijkstra et al. 2016; Verhamme et al. 2017; Steidel et al. 2018; Izotov et al. 2022; Yuan et al. 2024; Choustikov et al. 2024a; Pahl et al. 2024; Gazagnes et al. 2024), although any relation between these two photon escapes tends to be reported with significant scatter between individual galaxies. Since the escape of Lyman continuum photons is extremely challenging to observe directly at redshifts higher than z ∼ 3 due to the intervening IGM, the escape of Lymanα photons can be an important proxy.

Many e ff orts have been made towards a better understanding, both on the theoretical and observational front, and a better quantification of the escape fraction (henceforth fesc Ly , α ). The wellknown relations in Kennicutt (1998a) quantify the basic relationship between Lymanα luminosity and the star-formation rate (SFR) expected to produce it, subject to assumptions on the initial mass function (IMF), stellar metallicity, and star-formation history (SFH). However, the resonance of the Lymanα line means that Lymanα photons will scatter in a galaxy's neutral hydrogen. During this process, dust may absorb these photons (Schaerer &amp; Verhamme 2008; Ciardullo et al. 2014; Dijkstra 2016 and references therein), so the Lymanα escape fraction is expected to decrease with increasing dust content (Atek et al. 2008; Runnholm et al. 2020). The size and stellar mass of galaxies may also play a role, with fesc Ly , α generally increasing with lower stellar mass (Oyarzún et al. 2017; Yang et al. 2017; Goovaerts et al. 2024).

The increasing neutrality of the IGM at redshifts above z ∼ 5 -6 (McGreer et al. 2018; Planck Collaboration et al. 2020) is expected to play a significant role in Lymanα visibility. Neutral hydrogen around the galaxy and along the line of sight absorbs Lymanα emission, meaning that the number density of LAEs decreases around these redshifts (Schaerer et al. 2011; Pentericci et al. 2011; De Barros et al. 2017; Arrabal Haro et al. 2018; Pentericci et al. 2018). However, this e ff ect is subject to significant uncertainties, and for individual cases, it is highly dependant on the physical conditions and environment of the galaxy in question (Goovaerts et al. 2023; Witten et al. 2024).

Despite these challenges, strides have been made in recent years towards quantifying fesc Ly , α and its evolution with redshift. Blanc et al. (2011), henceforth B11, studied the evolution of dust properties and fesc Ly , α across a sample of ∼ 100 blank-field LAEs in the range 1 9 . &lt; z &lt; 3 8 and compared star-formation rate den-. sities (SFRDs) derived from LAE and UV luminosity functions (henceforth LFs) in order to quantify the evolution of fesc Ly , α with redshift. They found a sample median of ∼ 29% when considering their LAEs individually and a negligible redshift evolution.

Hayes et al. (2011) found a significantly lower f esc Ly , α , around 5% at z = 2 2 using a blind narrow-band survey search-. ing for H α and Lymanα . The authors made use of the H α line to estimate un-obscured star formation and compared this to the star formation derived from Lymanα . The H α line originates from the same process as the Lymanα line and is less sensitive to dust attenuation; therefore, it is often a more reliable indicator of star formation.

One general agreement in the literature is that f esc Ly , α increases with redshift. Low values such as 0 -5%are found in the aforementioned study as well as Hayes et al. (2013b), Verhamme et al. (2017), Sobral et al. (2017), Sobral &amp; Matthee (2019), and Runnholm et al. (2020). Larger values such as 30 -40% and even higher have been found at redshifts towards the epoch of reionisation (B11; Hayes et al. 2011; Chen et al. 2024; Witstok et al. 2024a; Lin et al. 2024; Napolitano et al. 2024).

Another general agreement found in the literature relates to sample selection. Samples that are selected by their Lymanα emission typically display higher values of fesc Ly , α when compared to samples selected by H α (Song et al. 2014; Trainor et al. 2016; Sobral &amp; Matthee 2019; Matthee et al. 2021).

By contrast, significant scatter in fesc Ly , α is expected in all samples, likely due to the large number of factors that a ff ect Lymanα photon escape and the di ff ering conditions present in any sample of galaxies. This can be seen in all studies on fesc Ly , α , from low redshift (Hayes et al. 2010, 2013b; Ciardullo et al. 2014; Matthee et al. 2016b; Runnholm et al. 2020) to the higher redshift, such as B11, Chen et al. (2024), Lin et al. (2024), and Napolitano et al. (2024).

Currently, the most powerful estimators for determining f esc Ly , α in an indirect manner include EWLy α (Sobral &amp; Matthee 2019) and Lymanα line separation (Yang et al. 2017; Matthee et al. 2021). The correlation with Lymanα line separation is likely due to the presence of outflows, to which the Lymanα line shape is sensitive (Gronke &amp; Dijkstra 2016; Blaizot et al. 2023). These outflows produce channels of low HI column density in the interstellar medium (ISM) and circumgalactic medium (CGM), which allows for more Lymanα escape (Hashimoto et al. 2015).

We can also constrain fesc Ly , α in a global way by comparing the SFRDs derived from the LAE and general UV-selected populations (with SFR derived from the UV emission). The ratio of these gives the escape fraction of Lymanα photons over a specific survey volume. Employing this method, Konno et al. (2016), B11 and Hayes et al. (2011) found a gentle evolution over their redshift ranges that could be fitted with a power law of the form: fesc Ly , α = C × (1 + z ) ξ , where C = (4 79 . + 5 68 . -0 69 . ) × 10 -4 , ξ = 3 38 . + 0 10 . -0 37 . at z ∼ 2 2 . -6 7 (Hayes et al. 2011), . C = 3 98 . + 2 32 . -1 46 . × 10 -3 , ξ = 2 2 . ± 0 3 at . z = 1 9 . -3 8 (B11) and . C = 5 0 . × 10 -4 , ξ = 2 8 . in z = 0 -6 (Konno et al. 2016).

Hayes et al. (2011) mentioned that fesc Ly , α may reach unity at z ∼ 11 + 0 8 . -0 6 . , while Cassata et al. (2011) argued that fesc Ly , α would reach unity around a redshift of six. In general, such research has been done with a few LAEs fainter than log LLy α [erg s] / &lt; 40, hence without the ability to fully constrain the faint part of the LAE LF and its impact on the SFRD.

Any di ff erence in the evolution of fesc Ly , α may come from the integration limits during the estimation of the Lymanα and UV luminosity densities. Hence, data on the fainter luminosity regime is necessary in order to integrate to fainter luminosity limits and better understand these faint populations' e ff ect on the evolution of the LF as well as fesc Ly , α with redshift.

Fortunately, the recent VLT Multi Unit Spectroscopic Ex-/ plorer (MUSE) observations of lensing clusters (Claeyssens et al. 2022; Thai et al. 2023) and blank fields (Vitte et al. 2024) have allowed us to observe Lymanα down to log LLy α [erg s] / ∼ 39, 39 5, respectively. .

In this paper, we seek to combine two state-of-the-art LFs of intrinsically faint lensed galaxies, the latest UV LF from Bouwens et al. (2022) and the latest LAE LF from Thai et al. (2023), allowing us to compare the SFRD derived from Lymanα emission, SFRDLy α , to the total SFR density, SFRDtotal, and hence track the redshift evolution of fesc Ly , α over a wide redshift range: 2 9 . &lt; z &lt; 6 7. In Hayes et al. (2011), the independent LF . comparison is advocated as a method to calculate fesc Ly , α , though it is subject to cosmic-variance related errors. We can significantly reduce the uncertainties related to cosmic variance by using six lensing clusters for the UV LF and 17 lensing clusters for the LAE LF. The LFs we used provide this study with a unique opportunity to explore the faint galaxy regime (MUV &lt; -13 and log LLy α [erg / s] &gt; 39) without the need for extrapolation in the LFs.

We also calculated fesc Ly , α on an object-per-object basis in order to compare galaxy properties as well as the global fesc Ly , α evolution, based on spectral energy distribution (SED) fitting of HSTand JWST NIRCam data of / ∼ 100 LAEs in the A2744 lensing cluster. The NIRCam data, extending to ∼ 5 µ m , allowed us to access, for the first time, the rest-frame optical emission of these LAEs and greatly improves the reliability of the properties extracted with SED fitting. The NIRCam imaging notably

also covers the entire MUSE field of view, meaning the full sample size of MUSE LAEs can be explored. This would not be the case with calculations of SFR from H α and dust attenuation from the Balmer decrement obtained from spectroscopy, for example from JWST NIRSpec's MSA mode. This would, in principle, be / preferable, due to the fewer assumptions needed and our lack of photometric coverage of the infrared for our LAEs. However, samples of LAEs with Balmer line detections are rare at the redshifts concerned by this study, and sample sizes are often small. Using photometry, we could gain an appreciable sample size using just one cluster, although more would be preferable (and soon possible with the current and upcoming JWST observations).

Studies of emission lines with slit spectroscopy for these galaxies (especially lensed LAEs and therefore potentially extended sources, both by lensing and due to the extended nature of Lymanα haloes) are not without di ffi culties themselves. Parts of flux can be missed in slits, especially slits as small as JWST NIRSpec's MSA (0 2 to 0 45; Napolitano et al. (2024) / ′′ . ′′ . provide a good discussion of these di ffi culties, see also Jung et al. (2024)), and sometimes, flux can be missed altogether, such as in Jiang et al. (2023).

Section 2 outlines the determination of fesc Ly , α for individual galaxies, describing the data and process used. In Section 3, we describe the global evolution of fesc Ly , α , the two LFs used for this calculation, and the results thereof. A comparison of the two methods to each other and to studies in the literature is presented in Section 4 along with consequences for reionisation and a discussion of the validity of comparing the UV and LAE LFs. Finally, Section 5 o ff ers conclusions and an outlook for fesc Ly , α determination in the near future.

Throughout this paper, we adopt a value for the Hubble constant of H 0 = 70 km s -1 Mpc -1 , and the cosmology used is ΩΛ = 0 7 and . Ω m = 0 3. The adopted IMF is that of Salpeter . (1955). All values of Lymanα luminosity and absolute UV magnitude are given corrected for magnification.

## 2. Escape fraction of individual LAEs

## 2.1. Data: Combining spectroscopy and photometry

In order to perform analysis of this genre, photometry and spectroscopy of su ffi cient quality and depth is essential. For our purposes we combined public MUSE integral field unit (IFU) spectroscopic data (Richard et al. 2021) (094.A-0115, 095.A0181, 096.A-0496) with the latest, deepest combination of Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) photometry from the UNCOVER survey (PIs Labbé and Bezanson, JWST-GO-2561, Bezanson et al. (2022)) of the lensing cluster A2744. This cluster is extremely well studied and benefits from a great number of multiple images to constrain its lens models. Below, we detail the two sets of observations and how they are used for this study.

## 2.1.1. MUSE IFU spectroscopy

The use of an IFU allowed us to blindly select LAEs, rather than performing spectroscopic follow-up on UV-selected targets. This ensures our sample has a simple reproducible selection function, namely Lymanα luminosity limited. MUSE has a field of view of 1 × 1 arcmin 2 and a spectral resolving power of R ∼ 3000 at λ = 808 nm (Bacon et al. 2015). Wavelength coverage ranges between 4750 Å and 9350 Å, meaning the Lymanα line can be detected between redshifts of 2.9 and 6.7.

The A2744 lensing cluster used in this work has integration times varying between 3.5 and 7 hours. It forms part of a data release by Richard et al. (2021) 1 and the LAEs used in this work are comprehensively described as part of the Lensed Lymanα MUSE Arcs (LLAMA) Sample, in Claeyssens et al. (2022).The full process of LAE detection in MUSE data is detailed in Richard et al. (2021) and Weilbacher et al. (2020). We give an abridged version here.

Emission line sources are identified in MUSE narrow-band datacubes using the MUSELET software (Piqueras et al. 2019) . 2 Thereafter the Source Inspection package (Bacon et al. 2023) is used to identify line emitters as LAEs and assign redshifts. During this process, the redshift confidence of each LAE is determined, based on the emission lines seen, the shape of the Lymanα line, ancillary HST data and lensing considerations. The confidence scale ranges from 1 for a tentative redshift attribution to 3 for a very secure redshift attribution. In this work, only sources that have been identified with redshift confidence levels of two and three are used, as well as a S N greater than / three.

The flux of these LAEs is derived using SExtractor (Bertin &amp; Arnouts 1996) on a continuum-subtracted Narrow-Band subcube with a size of 15' × 15' × 20Å. The full process is described in de La Vieuville et al. (2019).

## 2.1.2. HST + JWST photometry and galaxy properties

The photometric catalogues of the UNCOVER project (Weaver et al. 2023) comprise photometry in 15 bands ranging from the F435W HST filter to the F444W JWST NIRCam filter. This / gives excellent, deep coverage up to 4.4 µ m, which corresponds to ∼ 580 nm in the rest frame at z = 6 7, the highest redshift of . the MUSE LAEs. We used the catalogue optimised for faint and compact sources (0 32 . ′′ apertures). 5 σ depths in the 15 filters range between 27.16 and 29.56 in AB magnitudes. Of the 154 total LAE images, 121 were detected, which corresponds to 99 individual LAEs. We rejected two of these LAEs as being clearly contaminated by bright cluster galaxies and a further LAE, as it lies on a critical line and has a magnification factor of 137 ± 1500. As such, its properties would likely not be well constrained (see, for example, Limousin et al. 2016). Having thus blindly matched these Lymanα -selected galaxies to their UV counterparts, their SEDs were fitted with the SED-fitting code CIGALE (Boquien et al. 2019). The details of the fitting procedure are outlined in Goovaerts et al. (2024); here we give a brief outline.

Taking advantage of the secure spectroscopic redshifts of these LAEs, we fit their SEDs using two di ff erent SFHs: a single exponentially decaying burst with a range of timescales and a double burst model with a second, delayed burst at a time τ after the initial burst, to be fit (Małek et al. 2018). Ages for the bursts range from 10 to 700 Myr and mass fractions in the second burst range from 0.001 to 0.65. The stellar template library of Bruzual &amp; Charlot (2003) is employed with metallicities ranging from 0.001 to 0.02, together with the nebular emission models from Theulé et al. (2024). Gas metallicities range between 0.0004 and 0.02. Dust attenuation is taken into account by means of the Calzetti attenuation law (Calzetti et al. 2000).

The best-fit SFH and model was carried forward in this analysis (see Goovaerts et al. 2024 for comparisons between the two SFHs). One of the parameters in the CIGALE output is the intrin-

1 https://cral-perso.univ-lyon1.fr/labo/perso/johan.

richard/MUSE\_data\_release/

2 https://mpdaf.readthedocs.io/en/latest/muselet.html

Fig. 1. Redshift evolution of the escape fraction of Lymanα photons from LAEs in the A2744 lensing cluster colour-coded by the absolute UV magnitude of each LAE. The global medians of our sample and that in B11, are shown with red and cyan dashed lines respectively. The median uncertainty of the sample ( ∼ 18%) is shown by the black error bar, as if it was at the level of the global median, o ff set for clarity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00517v1-figure-416e77a5-4e78-4d7e-bbcc-786ae2f60614.jpg)

sic dust-corrected star formation, which is used in the computation of fesc Ly , α .

## 2.1.3. Lensing models

For studies involving lensing clusters, the lens models used are pivotal to the reliability of the results, as these govern the lensing magnifications that are assigned to each galaxy and therefore the correction made to their observed luminosity. The lens models used for this work for the A2744 lensing cluster are described in Richard et al. (2021). As the Hubble Frontier Fields (HFF) and A2744 in particular have been so well studied over the years, these clusters have very well constrained mass models, drawing on many multiple images identified behind each cluster. While systematic uncertainties are still associated to the particular lens model used (Acebron et al. 2017; Furtak et al. 2021), the quantity of multiple images for these particular clusters render the statistical error on the mass distribution as low as 1% (Richard et al. 2021). Typical, accepted magnifications in the sample range from 1.5 to 25, with the maximum magnification being 137. We note that escape fraction measurements for individual galaxies (as in the following Section), are not a ff ected by magnification.

## 2.2. Escape fraction

To calculate the escape fraction on an object-per-object basis for the 96 individual LAEs behind the A2744 lensing cluster, from the data described in Section 2.1, we compare the SFR inferred from the Lymanα flux, to the dust-corrected SFR from the CIGALE SED fitting process. The SFR inferred from the Lymanα flux is calculated using the prescription in Kennicutt (1998a) and the conversion factor of 8.7 between Lymanα and H α luminosities, using a Case B recombination scenario (Osterbrock 1989) ( T = 10 4 K).

The ratio of these two SFRs gives the escape fraction of Lymanα photons, fesc Ly , α . The redshift evolution of this quantity is shown in Fig. 1. The sample median is shown as the red dashed line and the sample median in B11 is shown as the cyan dashed line. They are consistent to within each other's errors, even though B11 use a di ff erent prescription for finding fesc Ly , α . The median escape fraction uncertainty of our sample (18%) is shown by the black error bar. Each LAE is colour-coded by its absolute UV magnitude, calculated from the filter that sees the rest frame 1500 Å emission. One can clearly see that Lymanα photons escape more easily from UV-fainter objects, a result already remarked upon by Lin et al. (2024), Napolitano et al. (2024). To clarify this further, we can fit a relation between the absolute UV magnitude and fesc Ly , α for each galaxy. This is shown in Fig. 2, where the galaxies are now colour-coded by redshift-bin and median error bars for each bin are shown in the corresponding colour. The best-fit to this data is

$$\Delta ^ { - \varepsilon \varepsilon } \quad \log ( f _ { \text{esc}, L y \alpha } ) = ( 0. 2 7 \pm 0. 0 2 ) \, M _ { \text{UV} } + ( 4. 2 \pm 0. 5 ). \quad \quad ( 1 )$$

This is clearly a tighter relation than the one shown by the redshift versus fesc Ly , α , with an uncertainty on the slope and normalisation of about 10 %. These quantities, MUV and fesc Ly , α are not independent as MUV is related to the intrinsic SFR of a galaxy, so this relation is to be expected, but it is nevertheless noticeable that below MUV = -18 there are few galaxies with f esc Ly , α values less than 10%. However, care must be taken when interpreting this result, as incompleteness can play a role as we start to probe fainter galaxies (towards the right of Fig. 2). To better assess this, in Fig. 2 we plot three regions where we expect to be incomplete as a function of UV magnitude and lensing magnification. This is based on the limiting flux for a MUSE detection (of an LAE). The three magnification regions plotted are µ = 1 5, the minimum magnification included in the sample . (the area covered by MUSE observations do not reach the areas with less magnification than this, most of the area has µ &gt; 2), µ = 9, the maximum magnification for an image included in the sample (we note that there are LAEs not included and images not included with µ &gt; 9) and µ = 137, the maximum magnification of any LAE image in the A2744 lensing cluster. These three regions delineate the area of this plot where we expect to become incomplete in terms of detecting a faint galaxy with a small fesc Ly , α value with MUSE. We can see that this effect likely becomes important towards the faint end of this plot, at MUV ≳ -17. However, at the brighter end of this graph, we do not expect the incompleteness to shape the relation we observe.

We did see LAEs undetected in the continuum that might be expected to fill this region of the graph. However, they are uniquely objects of high escape fraction (based on upper limits of the UV continuum, and no dust correction). These objects are denoted by the black triangles in Fig. 2: lower limits on fesc Ly , α and upper limits on MUV.

We note also that around this MUV value of -18 that we start to see objects with fesc Ly , α &gt; 100%. fesc Ly , α &gt; 100% means that we are seeing more Lymanα emission than we should, given the standard conversions between Lymanα and UV luminosities and SFR detailed above. Reasons for this may include our assumptions in making those conversions, such as constant SFH, ages greater than ∼ 50 Myr and the case B recombination scenario (Osterbrock 1989), becoming less valid (see, for example, McClymont et al. 2024). This may point concurrently to our models failing to correctly estimate the dust content of these galaxies and thereby underestimating the intrinsic SFR. However, this would mean that some property of these galaxies is heavily attenuating the UV emission, but not attenuating the Lymanα emission accordingly. Atek et al. (2008) postulate that geometry or the presence of large outflows may play a role in disconnecting Lymanα emission and dust such that one may see a high fesc Ly , α from a heavily dust-attenuated system.

We also observed LAEs for which the continuum is completely undetected in the deepest photometry to date (Bacon et al. 2015; Maseda et al. 2020; de La Vieuville et al. 2020; Goovaerts et al. 2023; Maseda et al. 2023), and fesc Ly , α for these galaxies is, naturally, very high, often greater than 100% (as calculated from 2 σ upper limits on the UV continuum). This suggests that there is more than simply model inadequacy at play here. We discuss these objects further in Section 4.

Unlike MUV, we see no significant evolution of fesc Ly , α with redshift, although we have relatively few objects at redshifts higher than 6. It is clear from Fig. 1 that we are a ff ected by missed objects above z = 6, both in the high-fesc Ly , α regime (as these objects have low MUV) and the low-escape fraction regime, as objects with low fesc Ly , α are harder to detect at higher redshifts. The greater presence of skylines and decreased MUSE sensitivity towards longer wavelengths are additional challenges to detect LAEs towards the high end of our redshift range.

The comparison to the results in B11 suggest that fesc Ly , α does not evolve significantly over the redshift range 1 8 . &lt; z &lt; 6. This is in tension with certain literature results such as Hayes et al. (2011), who find that fesc Ly , α increases with redshift (0 6 . &lt; z &lt; 6). However, this study was performed in a similar way to the global analysis of fesc Ly , α that we present in the following sections. Hence, we return to this point after having taken both methods into consideration.

The scatter in the individual results shown in Fig. 1 is to be expected, given the wide variety of conditions within galaxies, even LAEs, expected to be mostly young and highly star forming. Dust content of these galaxies likely varies, which plays a significant role in the escape of Lymanα photons (Verhamme et al. 2008). Additionally and as mentioned, geometry and lineof-sight e ff ects are also expected to have a hand in the amount of Lymanα photons that reach the observer (Verhamme et al. 2012; Dijkstra &amp; Kramer 2012; Gronke &amp; Dijkstra 2016; Blaizot et al. 2023; Giovinazzo et al. 2024).

While we see consistency between the sample median of B11 and ours, the o ff set, as well as the smaller scatter of our results, is to be expected based on the di ff ering methods for finding f esc Ly , α . B11 assume an intrinsic UV slope of -2 23 in order to . calculate the dust reddening E B ( -V ). As this UV slope likely depends on factors other than dust, such as metallicity, age and mass of the stellar population (Buat et al. 2012; Wilkins et al. 2013; Bouwens et al. 2016; Reddy et al. 2018; Topping et al. 2022), and as there is evidence for slopes bluer than this value (De Barros et al. 2017; Naidu et al. 2022; Goovaerts et al. 2023; Iani et al. 2023; Nanayakkara et al. 2023; Napolitano et al. 2024), we prefer not to use this assumption. If there are bluer intrinsic slopes than this, we would expect the corresponding E B ( -V ) values from B11 to be higher, and therefore, the intrinsic SFR would be higher, and the calculated escape fractions would be lower. This is what we observed in our results. We also saw less scatter in fesc Ly , α , as we do not depend on estimations of the UV slope, which are known to have significant scatter (as well as uncertainty and possibly bias, e.g. Austin et al. 2024).

Our sample median, 22 5 %, agrees with determinations at . slightly lower redshift from Sobral et al. (2017); Matthee et al. (2021), who both quote escape fractions of ∼ 30 % for LAEs. Most recently, also using JWST, Chen et al. (2024) have found a sample median of 28 % for 10 LAEs at 5 &lt; z &lt; 8.

However, we note that the scatter on individual determinations of fesc Ly , α remains a very significant consideration when assessing samples of the kind in Fig. 1. It is therefore unlikely that we will be able to specify a precise fraction for Lymanα escape at any given redshift. It is clear that the trend of fesc Ly , α

Fig. 2. Evolution of fesc Ly , α with absolute UV magnitude MUV. Photometrically detected galaxies are split into three redshift bins to stress the lack of evolution in this sense, while the evolution with MUV is fit by the relation in cyan (details in text). Median error bars for the three redshift bins are shown in the corresponding colour and fesc Ly , α error bars are also listed here for clarity (relevant redshift bin in brackets): 0 18 (2 9 . . &lt; z &lt; 4), 0 12 (4 . &lt; z &lt; 5) and 0 25 (5 . &lt; z &lt; 6 7). Abso-. lute magnitude values and their errors are adjusted for magnification. Black triangles denote continuum-undetected LAEs. MUV values are upper limits and escape fraction values are lower limits. These objects are not included in the fit. Red shaded areas denote regions of incompleteness based on MUSE line-detections: a darker colour representing increasing lensing magnification (see text).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00517v1-figure-39fd36b9-7289-4dee-8751-1360c6b46644.jpg)

with MUV is more useful.

We note also that while individual fesc Ly , α values are better determined using H α fluxes and the Balmer decrement, the intrinsic scatter seen remains similar, with values ranging from &lt; 1% escape to 100% (see, for example, Chen et al. 2024; Napolitano et al. 2024). For our sample, we can take additional encouragement from the fact that the dust content of these lowmass LAEs, at these redshifts, is both expected to be, and observed to be, very low. We therefore do not expect many cases where our SED fitting significantly underestimates the dust attenuation of a galaxy. Both methods of determining fesc Ly , α are valuable then, and a more detailed comparison between them for the same sample of galaxies is desirable.

## 3. Evolution of the global Lymanα escape fraction

One can derive estimates of fesc Ly , α in a global manner by comparing the LF of LAEs and that of the 'total' galaxy population. Integrating both of these LFs gives the Lymanα luminosity density ρ Ly α and the total luminosity density ρ total . From this one can derive the respective SFRDs. Then taking the ratio of these SFRDs gives fesc Ly , α (Hayes et al. 2011; Blanc et al. 2011; Zheng et al. 2013; Konno et al. 2016). In order to perform this calculation, we start from the two state-of-the-art LFs in lensing fields, the LAE LF from Thai et al. (2023) and the UV LF from Bouwens et al. (2022), assumed here to be the 'total' galaxy population. This assumption and its consequences for the results of this study are discussed in detail in Sect. 4.3.

Before presenting the derivation of fesc Ly , α , we briefly detail both LFs used, starting with the LAE LF from Thai et al. (2023).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00517v1-figure-39699c53-dae1-417b-8f17-a8f77a95ef17.jpg)

Fig. 3. Ultraviolet and LAE LFs included in our analysis on a common x-axis of SFR. The LFs of Bouwens et al. (2022) and Thai et al. (2023) are shown in bold and the literature LFs are shown in coloured dash-dotted lines. The LF integration limits in terms of L ∗ that we use throughout this work are denoted by the pink shaded regions. Left: Results from the redshift range 3 to 4. Right: Results from the redshift range 5 to 6. The literature LFs in the graph are shown in shorthand in the caption and here in full, in order: Ouchi et al. (2008); Cassata et al. (2011); Sobral et al. (2018); Thai et al. (2023); Drake et al. (2017); Bouwens et al. (2022, 2009).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00517v1-figure-df85c9b3-4754-4091-a94c-fbf5d60b507e.jpg)

## 3.1. Lymanα luminosity function

The LAE LF is computed by using the non-parametric V max method to estimate the volume of the survey in which we can detect sources individually. The inverse of this parameter is the contribution of the source to the numerical density of the galaxies. The procedure to compute this parameter using data collected by lensing observations by MUSE VLT was developed by / de La Vieuville et al. (2019) and presented in Thai et al. (2023) with some improvements related to the magnification and completeness values of the sources. These completeness values reflect the expected number of galaxies actually in the field based on the observed galaxies and the noise statistics of the field. The sample in Thai et al. (2023) includes 600 LAEs behind 17 lensing clusters observed by MUSE VLT, identified with highly / secure redshifts, zconf = 2 3 (high Lyman-, α line S N, clearly / asymmetric line profile, etc) as defined in Richard et al. (2021). These data are within a redshift interval 2 9 . &lt; z &lt; 6 7 and cover-. ing four orders of magnitude in luminosity 39 &lt; log L [erg s -1 ] &lt; 43. These luminosity values were corrected for the magnification of the lensing e ff ect using models from the work of Pello et al. (1991); Covone et al. (2006); Richard et al. (2010); Mahler et al. (2018); Lagattuta et al. (2019); Beauchesne et al. (2023); Richard et al. (2021).

exp -( LT / L m ) to the traditional Schechter function. The bestfit results suggest that m is about unity and log LT [erg s -1 ] is ∼ 40 and ∼ 40 7 for the redshift intervals . z 45 and z 56 respectively. Accounting for the e ff ects of source selection, completeness cut, and flux measurement, the final slope values in the four redshift intervals are -2 00 . ± 0 50, . -1 97 . ± 0 50, . -2 28 . ± 0 50 and . -2 06 . ± 0 60. These values are in good agreement with results . from other literature studies within a 1 σ deviation. The best-fit value of the characteristic luminosity values are log L ∗ [erg s -1 ] = 42 85 . + 0 10 . -0 10 . , 42 87 . + 0 11 . -0 1 . , 42 97 . + 0 13 . -0 11 . , 43 09 . + 0 10 . -0 08 . and the normalisation factors are ϕ ∗ [10 -4 Mpc -3 ] = 7 41 . + 2 70 . 2 70 . , 6 56 . + 3 20 . -2 40 . , 4 06 . + 2 70 . -1 70 . and 3 49 . + 2 11 . + 1 50 in redshift bins zall , z 35 , z 45 and z 60 , respectively.

The evolution of the LF with redshift was investigated in four redshift intervals 2 9 . &lt; z &lt; 6 7, 2 9 . . &lt; z &lt; 4 0, 4 0 . . &lt; z &lt; 5 0 and . 5 0 . &lt; z &lt; 6 7 (labelled . zall , z 35 , z 45 and z 56 ) and the Schechter function used to fit the LF points. Using the data sample collected with the help of gravitational lensing is e ffi cient to study the LF at faint Lymanα luminosities (log L [erg s -1 ] &lt; 41). However, it becomes less sensitive when examining the function surrounding the characteristic luminosity value L ∗ of the Schechter function. To tackle this problem, Thai et al. (2023) used the average LF values obtained from the literature in the same luminosity and redshift bins as a constraint at the bright end of the LF. A turnover in the shape of the LF was observed at luminosities log L [erg s -1 ] &lt; 41 when studying the evolution in the z 45 and z 56 redshift intervals. This was explained by inefficient cooling of gas in small dark matter haloes (Jaacks et al. 2013; Gnedin 2016; Yue et al. 2016). To accommodate the potential turnover in these two redshift bins, a modified Schechter function was introduced by multiplying an exponential term

. In our analysis we also include results for the LAE LF at di ff erent redshift intervals from literature studies to build a general picture, and provide brief descriptions of these works here. de La Vieuville et al. (2019) studied the LF evolution of ∼ 120 LAEs behind four lensing clusters observed by MUSE, with the same redshift range as Thai et al. (2023) and a luminosity range of 39 ≲ log LLy α [erg / s] ≲ 43. In the blank field observations, Ouchi et al. (2008); Blanc et al. (2011); Cassata et al. (2011); Drake et al. (2017); Konno et al. (2016) presented the LF evolution with redshift of LAEs with luminosities greater than log L [erg s] / = 41. Due to lacking data at faint luminosities, Ouchi et al. (2008) fixed the slope value of the Schechter function at α = -1 5 ( . z = 3 1 . , 3 7 . , 5 7). B11 fixed this parameter at . α = -1 7 . ( z = 1 9 . -3 8), while Cassata et al. (2011) found . α = -1 6 . + 0 12 . -0 12 . , -1 78 . + 0 10 . -0 12 . , -1 69 (fixed) at . z = 1 95 . -3 0 . , 3 0 . -4 55 4 55 . , . -6 6, . respectively. However, at a higher redshift and a brighter luminosity range, this value has been found to be much steeper, varying from -2 26 to . -2 56 at redshift . z = 5 7, and from . -1 86 . to -2 5 at redshift . z = 6 6 (Konno et al. 2016); or . ∼ -1 98 at . redshift z = 5 4 . ± 0 4 (Sobral et al. 2018), depending on the . luminosity limitations used for fitting.

## 3.2. Ultraviolet luminosity function

Bouwens et al. (2022) presented a new derivation of the UV LF in the redshift range 2 0 . &lt; z &lt; 9 0 including over 2500 lensed . galaxies behind the HFF clusters using HST observations from the HFF (Lotz et al. 2017), Grism Lens-Amplified Survey from Space (GLASS) (Treu et al. 2015) and Cluster Lensing And Su-

pernova survey with Hubble (CLASH) (Postman et al. 2012) programmes. This LF probes down to MUV &lt; -16 over the full redshift range we are considering in this study. A smooth flattening of α was found, from α = -2 28 . ± 0 10 at redshift . z = 9 to α = -1 53 . ± 0 03 at redshift . z ≈ 2. Over our redshift range, 3 ≲ z ≲ 7, the faint-end slope is found to evolve from α ∼ -1 60 . to α ∼ -2 05. The shape of the LF and its evolution agrees well . with other results in the literature, including determinations from blank fields (see Parsa et al. 2016, as well as other references in Bouwens et al. 2022).

We treat this LF as the LF of the total 'parent' population in order to calculate the total SFRD. In this way, inclusive of completeness corrections, we assume that both LFs represent the 'total' of their respective populations within the quoted Lymanα luminosity and absolute UV magnitude limits. This is an assumption that we return to and discuss in Section 4.3.

All LFs, UV and LAE, are shown on a common scale (SFR) in Fig. 3. The LFs of Thai et al. (2023) and Bouwens et al. (2022) are highlighted in thicker, solid lines.

## 3.3. Global evolution results

The results from the Lymanα LF in each redshift interval are integrated with luminosity to map out the evolution of Lymanα luminosity density with redshift:

$$\rho _ { L y \alpha } & = \, \int L _ { L y \alpha } \, \phi _ { L y \alpha } \, d L _ { L y \alpha }. & ( 2 ) \, \begin{array} { c } \rho _ { 1 1 } \, \text{and} \\ C \text{ is con} \\ \text{two red} \end{array}$$

This is then converted to the Lymanα SFRD based on the calibration of Kennicutt (1998a), assuming the factor of 8.7 between intrinsic Lymanα and H α flux and a case B recombination (Osterbrock 1989):

$$S F R D _ { L y \alpha } [ M _ { \odot } y r ^ { - 1 } M p c ^ { - 3 } ] = 7. 9 \times 1 0 ^ { - 4 2 } \times \rho _ { L y \alpha } / 8. 7. \quad \ ( 3 ) \quad \text{from } \mathcal { O } u \ \infty.$$

The UV luminosity density is converted into the equivalent SFRD using the form SFRD = κ × ρ UV (Madau &amp; Dickinson 2014), in which κ is the conversion factor, sensitive to the recent SFH, metal-enrichment history and the IMF. The term ρ UV is the UV luminosity density, which is expressed in units of erg s -1 Hz -1 , estimated by integration the UV LF with MUV. The widely used value of the conversion factor from Kennicutt (1998b) is κ = 1 4 . × 10 -28 , where a Salpeter IMF (Salpeter 1955) is assumed. Bouwens et al. (2022) used a Chabrier IMF (Chabrier 2003) for this conversion, so we correct for this, using a factor of 1.96. Going forward, the SFRD obtained from UV LF is treated as the total SFRD of the entire galaxy population. We estimated the total UV SFRD down to MUV = -13 at redshifts z = 2 3 4 5 6 7 (see Table 1 for the values used, after , , , , , correction to the Salpeter IMF). Subsequently we interpolated the SFRD to redshifts z = 2 5 . , 3 1 . , 3 5 . , 3 7 4 5 . , . , 4 7 5 4 . , . , 5 7 to . adapt to redshift bins used in the literature studies.

The ratio between the LAE SFRD and the total, dustcorrected SFRD in a given redshift interval gives the global Lymanα escape fraction. This definition was applied in the works of Hayes et al. (2011), Cassata et al. (2011), B11, Zheng et al. (2013) and Konno et al. (2016). The results for the global escape fraction derived in this work are shown in Figs 4 and 5. We compared our results to the results derived using several LAE LFs in the literature, for which we adopted an identical approach to that described above (Ouchi et al. 2008; Blanc et al. 2011; Cassata et al. 2011; de La Vieuville et al. 2019; Thai et al. 2023; Sobral et al. 2018). We used the same UV LF for each comparison.

The lower limits of the respective LF integrations are worth noting. We used two di ff erent limits, the first was 0 04 . L ∗ and the second was 0 0006 . L ∗ (where L ∗ is the characteristic galaxy luminosity of the Schechter function). The first limit, as suggested in Hayes et al. (2011) is designed to compare our study to previous results using brighter galaxies. It corresponds to limits of log LLy α [erg s -1 ] ∼ 41 and MUV ∼ -17 respectively. The results for f esc Ly , α calculated with this limit are shown in Fig. 4. The second limit, 0 0006 . L ∗ , corresponds to the faint regimes that we reach with the LAE and UV LFs (log LLy α [erg s -1 ] ∼ 39 5 . and MUV ∼ -13) and is shown in Fig. 5. The vertical error bars are estimated from published uncertainty values in ϕ ∗ , α , L ∗ (LAE LFs) and ϕ ∗ , α , M ∗ (UV LF). The horizontal error bars display the redshift intervals of each survey. As de La Vieuville et al. (2019) and Thai et al. (2023) probed the LAE LF at the same redshift intervals, the data points derived from the LF of de La Vieuville et al. (2019) have been shifted by 0.1 in redshift. Table A1 presents the results for each study using this integration limit, as well as the SFRD value derived. By comparing the results obtained with the di ff erent limits we can assess the impact of the faint galaxies that we reach with the latest LAE and UV LFs.

To investigate the global evolution of escape fraction with redshift, we recall the power-law form that was proposed by B11 and Hayes et al. (2011): f Ly α escp = C × (1 + z ) ξ . The coe ffi cient C is constrained strongly by the evolution of escape fraction at low redshift while the exponent ξ is defined by the evolution at high redshift. In Fig. 4, the best fits calculated for this work are shown as red and blue dash-dotted lines, fitting all the points in the graph and only those of Thai et al. (2023), respectively. In Fig. 5, the red fit is to all the data. The green dashed line shows the fit from Hayes et al. (2011), which used the UV LF obtained from Ouchi et al. (2004); Reddy et al. (2008); Bouwens et al. (2009) and an LF integration down to 0 04 . L ∗ (fitted to data over the range 2 &lt; z &lt; 6). We leave this fit in both plots to guide the eye as to the di ff erence when integrating to a fainter limit. B11 use the results of Bouwens et al. (2010) and an integration limit of 0. The fit of B11 is shown in Fig. 5 to provide the comparison to earlier results which have been integrated to 0.

We find C = 14 3 . ± 14 0 . × 10 -4 and ξ = 2 3 . ± 0 6 when fitting . to all the literature and C = 12 ± 43 × 10 -4 and ξ = 2 8 . ± 2 1 . when only fitting the results from Thai et al. (2023) with the integration limit at 0 04 . L ∗ . With the lower integration limit of 0 0006 . L ∗ , we find C = 4 1 . ± 7 7 . × 10 -4 and ξ = 3 6 . ± 1 2 . when fitting our results together with all the ancillary data in Fig. 5. Because the other studies included in this graph do not probe the faint regime that our LFs do (equivalent to the integration limit of 0 0006 . L ∗ ), and they therefore require an extrapolation outside the limits of the datasets used. During this fitting process, we set the weight of each point equal to the weight of the equivalent point from our study (i.e. the cyan point in the same redshift bin). This relation is qualitatively consistent with the equivalent relation derived in Hayes et al. (2011) (green dashed line, noting that this relation is derived using a di ff erent UV LF and a brighter integration limit) and suggests that fesc Ly , α will reach 100% around z = 8.

## 4. Discussion

With the individual and global results now in hand, we can compare them and consider what these can tell us about fesc Ly , α and its evolution.

Table 1. Ultraviolet SFRD used to calculate fesc Ly , α .

| Redshift   | SFRD unobscured              | SFRD obscured                | SFRD total                   |
|------------|------------------------------|------------------------------|------------------------------|
|            | log 10 [M ⊙ Mpc - 3 yr - 1 ] | log 10 [M ⊙ Mpc - 3 yr - 1 ] | log 10 [M ⊙ Mpc - 3 yr - 1 ] |
| z = 2 . 0  | - 1 . 28 ± 0 . 03            | - 0 . 8 ± 0 . 1              | 0 . 220 ± 0 . 043            |
| z = 3 . 0  | - 1 . 11 ± 0 . 07            | - 1 . 1 ± 0 . 1              | 0 . 163 ± 0 . 026            |
| z = 4 . 0  | - 1 . 12 ± 0 . 12            | - 1 . 2 ± 0 . 1              | 0 . 189 ± 0 . 031            |
| z = 5 . 0  | - 1 . 33 ± 0 . 09            | - 1 . 5 ± 0 . 1              | 0 . 078 ± 0 . 013            |
| z = 6 . 0  | - 1 . 46 ± 0 . 05            | - 2 . 0 ± 0 . 1              | 0 . 045 ± 0 . 005            |
| z = 7 . 0  | - 1 . 56 ± 0 . 06            | - 2 . 6 ± 0 . 1              | 0 . 030 ± 0 . 004            |

Notes. Values adapted from (Bouwens et al. 2022), adjusted for the conversion from the Chabrier IMF to the Salpeter IMF. Each value was calculated with an integration limit of MUV = -13.

Fig. 4. Global redshift evolution of fesc Ly , α with an integration limit on both LFs of 0 04 . L ∗ . Our results, using the LFs from Thai et al. (2023) are shown in cyan for the modified Schechter function and pink for the general Schechter function. The dotted horizontal line denotes f esc = 100%. The dash-dotted green curve is from Hayes et al. (2011) in the redshift range of 2 2 . -6 6 and calculated with the integration limit 0 04 . . L ∗ . The red dash-dotted curve is the best fit to all the data shown and the blue dash-dotted curve fits just the data of Thai et al. (2023) (cyan points: modified Schechter function).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00517v1-figure-ffb92c0d-4b2c-4a8f-82b4-841a77481076.jpg)

## 4.1. Global f esc Ly , α evolution

Fig. 4 shows the determination of fesc Ly , α using the LF integration limit 0 04 . L ∗ . We plotted this graph in order to compare our results to previous literature determinations of the escape fraction, which used samples of galaxies brighter than this work, and to assess the impact of the bright (MUV &lt; -17) and faint populations (MUV &lt; -13) on the evolution of fesc Ly , α . We plotted results derived using both the modified (cyan) and the general Schechter (pink) functions in order to show the di ff erence in these results created by the presence or lack of a turnover in the LAE LF. However, at this integration limit, the behaviour of the modified Schechter function is close to that of the traditional Schechter function, so we did not see a significant di ff erence in SFRD Ly α . As the current best derivations of the LF come from the modified Schechter functions, going forth we discuss the results related to these (cyan points).

results derived in this work and the results from the literature. We saw little evolution between z ∼ 3 5 and . z ∼ 4 5, followed . by a jump in fesc Ly , α at z ∼ 6, but at each redshift, our results agree well with the established literature. When fitting our data (blue line) compared to all the results (red line), we observed an enhancement of fesc Ly , α , likely due to the constraints placed by the literature in the redshift range 2 &lt; z &lt; 4. Our points do lie systematically above the points from the literature, indicating a greater SFRD based on the latest LAE LF of Thai et al. (2023). We note, however, that the LFs of Konno et al. (2016) and Drake et al. (2017) have not been included in Figs 4 and 5 due to large uncertainties, especially at the highest redshifts and faint luminosities. These uncertainties contribute to unrealistically high (and poorly constrained) escape fractions based on these LFs.

We observed a good agreement at all redshifts between the

Additionally, the choice of completeness cut has an impact. In Thai et al. (2023), a completeness cut of 1% was chosen (i.e. LAEs with completeness values &lt; 1% are not taken into account

in the LF determination), in order to maximise the sample size, while still negating the e ff ect of the most incomplete and therefore uncertain sources. Studies have typically chosen di ff erent completeness cuts, or none at all. For example, de La Vieuville et al. (2019) chose 10%, rejecting a larger number of sources. In order to check the e ff ect this has, we recomputed fesc Ly , α using the LF cut at 10% completeness. This yielded values similar to those seen in the literature: 10 -15% escape. Therefore, we consider our results completely consistent with what has been found previously when considering a brighter integration limit (0 04 . L ∗ ).

The green dashed fit of Hayes et al. (2011) lies slightly higher than the literature data. We attribute this to the significant update in the UV LF since the publishing of that study (the UV LFs of Bouwens et al. (2009) and Reddy et al. (2008) were used for this fit). This emphasises the importance of the UV LF in these calculations and the value of the latest determination from Bouwens et al. (2022). The values used in this work, inclusive of dust and IMF correction are shown in Table 1.

The values we derive for f esc Ly , α are considerably higher when we extend the LF integrations down to 0 0006 . L ∗ , particularly for the highest redshift bin, z ∼ 6, where we see an increase of 0.4 dex with respect to the 0 04 . L ∗ integration. In fact, this point is consistent with the fit to the data in B11, which is integrated to zero, that is, taking into account the full population by assuming the validity of the two LFs down to zero. This large escape fraction value is unsurprising, as we see a high number density of faint galaxies in the LAE LF, and these populations typically have higher fesc Ly , α (see Fig. 2, as well as Goovaerts et al. 2023; Chen et al. 2024; Goovaerts et al. 2024; Napolitano et al. 2024; Nakane et al. 2024). In the z ∼ 6 bin, f esc Ly , α is consistent with 100%, i.e. a similar SFRD derived from LAEs as from the UV-selected population. While still consistent at 1 σ with the fit derived by B11 (using an LF integration to 0), our point lies above points from the literature. We note that no other studies have had access to the faint regime of LAEs used in this study and as such, the literature LFs we compare to are integrated beyond the limits of their samples to create this graph. Reasons for this enhanced fesc Ly , α value for faint LAEs at z ∼ 6 include the e ff ect of the LAE incompleteness treatment, which acts most strongly for faint objects at the highest redshifts. E ff ectively, this would mean that we are counting in the LAE LF a significant number of faint LAEs not seen in the continuum.

This result may also point to a general scenario where faint galaxies (i.e. MUV ∼ -13) generally see an extremely high production and escape fraction of Lymanα photons (see also Atek et al. 2024). This appears to be a clear e ff ect when integrating down to the L ∗ limit our sample reaches: faint galaxies experience a higher escape of Lymanα photons, especially at higher redshifts. This result is also in agreement with the results presented in Fig. 2 as in this plot, between MUV ∼ -17 and MUV ∼ -13, we see LAEs with exclusively high escape fractions (almost all above 10% and many around 100%). This result implies that at z ∼ 6 and beyond, faint galaxies tend to be highly star forming and chemically unevolved, which should allow high escape fractions. On the other hand, seen in the uncertainties in the faint end of the LAE LF, it is clear that there is still a large intrinsic scatter in the populations (once more in agreement with the results in Section 2).

It is worth noting, however, that if the IMF of these faint, reionisation-era LAEs is top-heavy with respect to the assumed Salpeter IMF, this may also explain an elevated fesc Ly , α . Such an IMF would increase the number of ionising photons, from mas- sive stars, which would in turn increase the Lymanα emission of the galaxy. A higher ionising photon production e ffi ciency, ξ ion has been observed for some LAEs (Matthee et al. 2021; Saxena et al. 2024, see also Simmonds et al. 2023) although there is no firm consensus on how much higher. If ξ ion were indeed to be 25 5 . -25 6, this would mean a factor of 2 . -3 more ionising photons produced. This would in turn mean more Lymanα photons produced and an overestimation of fesc Ly , α in our results by a similar factor (we note that this would apply to all the results in Figs. 4 and 5). This is an important dependance on the high-mass slope of the IMF and the production of ionising photons, and could be particularly relevant to the faintest galaxies towards the epoch of reionisation, which are more likely to host exotic stellar populations.

In summary, we find that the LAE SFRD matches the UV SFRD during the final stages of reionisation, and is ∼ 20% of it at redshifts between 3 and 5. For both bright and faint integration limits, we see this jump in the highest redshift bin ( z ∼ 6) compared to little evolution between redshifts of 3 and 5. This reflects a significant evolution of the LAE LF in this redshift bin, therefore suggesting that the LAE population evolves significantly between the end of reionisation ( z ∼ 6) and later times ( z ∼ 4 5). . Although it is worth noting the uncertainty on our z ∼ 6 result that we attribute to a combination of several factors: a lower dust content in galaxies at this redshift and above combined with higher specific SFRs in our sample of faint galaxies. As mentioned, the completeness correction for the LAE LF acts most strongly on the faint, high-redshift LAEs, making this a possible additional contributor.

## 4.2. Comparison of the global and individual f esc Ly , α evolutions

It is immediately clear that the global evolution of fesc Ly , α shows a more significant trend than the individual evolution, despite significant uncertainties, especially in the highest redshift bin. This becomes additionally clear when incorporating the agreement in the literature for the global evolution.

For our datasets, this is likely due at least in part to the far greater sample size used in the global evolution: 2500 UVselected galaxies and 600 LAEs, compared to ∼ 100 LAEs for which fesc Ly , α is calculated (Fig. 1). For a smaller sample such as the one used for the individual tracking of fesc Ly , α , particularly considering the lack of su ffi cient galaxies at redshifts greater than 6, any redshift evolution is washed out by the scatter in individual results. The results of B11 display a similar lack of evolution over 1 9 . &lt; z &lt; 3 8 when taking into account the signif-. icant scatter of individual galaxies. The same is even true for the more recent studies of Chen et al. (2024) and Napolitano et al. (2024), using Balmer decrements to calculate fesc Ly , α , albeit with smaller sample sizes. The intrinsic variety of escape scenarios for Lymanα photons from di ff erent galaxies (even bearing in mind that all these systems are selected as LAEs) is evident.

The trend seen in the global evolution of fesc Ly , α (integration limit 0 0006 . L ∗ ), towards fesc Ly , α ∼ 80% around z = 6, is not seen in the individual results. Likely this is at least in part due to the incompleteness of the sample used to calculate fesc Ly , α for individual galaxies at this redshift, as both LAEs with very small escape fractions and LAEs with very large escape fractions are harder to observe. For the LAEs with large escape fractions this di ffi culty comes from the faintness of the associated continuum: the galaxies with the highest fesc Ly , α values tend to be faint in terms of UV magnitude (Fig. 2). This incompleteness should be accounted for in the two LFs used, so we do not expect this to

Fig. 5. Global redshift evolution of fesc Ly , α with an integration limit on both LFs of 0 0006 . L ∗ . Results from this study are shown by cyan (modified Schechter function) and pink (general Schechter function) points. Power-law fits from the literature are shown in grey and green (B11, Hayes et al. 2011) and the equivalent fit to all the data is shown in red. The fit from Hayes et al. (2011) is derived from results integrated to 0 04 . L ∗ , and the fit from B11 from results integrated to 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00517v1-figure-56e67f72-52f9-4aa1-ae1b-0df782a1cf2b.jpg)

play a significant role for the global fesc Ly , α evolution.

Contrasting these two results, global and individual, it is clear to see that a large sample size, as well as considered completeness corrections, are necessary to observe any significant redshift evolution of fesc Ly , α . Indeed, further exploration of the faintest regimes of LAEs (and UV-selected galaxies) would be desirable in order to better constrain the faint end of the LFs and reduce the uncertainties on the highest-redshift points. For example, it is still very challenging to observe more than a few galaxies with log LLy α &lt; 40 at z &gt; 5 5 (see, for example, Fig. 3 . of Thai et al. 2023).

ditionally, and as mentioned, the completeness correction added to the LAE LF acts strongly in the faint, high-redshift regime so we also expect this to play a role. The results in Fig. 5 suggest that this may be a stronger role than the equivalent completeness correction added to the UV LF.

We note also that the global comparison is e ff ectively the Lymanα escape fraction of a given volume of the Universe, bounded by the luminosity and redshift limits of the surveys used. This calculation takes into account all galaxies, not just those actually selected as LAEs, as is the case for the individual results. We might therefore reasonably expect to see a smaller escape fraction due to the inclusion of galaxies not selected as LAEs. In the lower and intermediate redshift bins, the results are consistent with the sample median shown by the red dashed line in Fig. 1; however in the highest redshift bin, the global calculation of fesc Ly , α appears to yield a larger value. Taking the median of the galaxies from Fig. 1 above z = 5 would yield a higher result (36% escape), and it is likely that we simply do not detect enough LAEs of high escape fraction at higher redshifts in our limited sample behind A2744. As this is just one cluster, extending these measurements to other clusters e ffi cient in lensing galaxies above z = 5 is very desirable.

The galaxies detected by their Lymanα emission but not in the continuum (the black triangles in Fig. 2) naturally all have high escape fractions and are counted in the LAE LF, so this may go some way towards explaining the discrepancy we observe between our di ff erent methods of calculating fesc Ly , α . Ad-

Contrary to the lack of similarity in redshift evolution in the global and individual fesc Ly , α results, there is a similar trend in terms of an evolution with the UV magnitude of the host galaxy. In Sect.2 this is established as a much more significant trend and the results when comparing LFs support this. When integrating to fainter limits we see significantly higher values of fesc Ly , α across the redshift range probed.

## 4.3. The overlap of the LAE and ultraviolet luminosity functions

Here we revisit the assumption made in order to compare the two LFs using single integration limits in terms of L ∗ and treat them as representative of their respective 'total' populations. This, in principle, is troublesome as there is no direct mapping between Lymanα luminosity and UV magnitude; for every given Lymanα luminosity bin, there are a range of UV magnitudes that galaxies display. In fact, this encodes exactly what we want to measure, the escape of Lymanα photons. It is clearly invalid to compare a LF based on just the brightest LAEs to the complete UV LF, including the faintest sources. Therefore, in order to make the comparison discussed above and shown in Figs. 4 and 5, we have to assume that we adequately represent the 'total' populations with both LFs, within the UV and Lymanα luminosity integration limits used. This approach is discussed at length in Hayes et al. (2011) and B11 and we revisit this in the context of our dataset.

That we adequately represent both Lymanα and UV popu-

Fig. 6. Redshift and MUV distribution of LAEs and Lyman-break galaxies (UV-selected galaxies) from the sample in Goovaerts et al. (2023). Grey points denote the Lyman-break galaxies and are thus analogous to the UV-selected sample. Orange points denote objects that were selected as LBGs and LAEs and upwards-facing triangles indicate LAEs that were not detected at all in the HST imaging (Shipley et al. 2018). The 2 σ upper limits for the continuum of these objects have been estimated using the local noise in the region where they would have been detected in the HST image that would see the 1500 Å emission. The horizontal dashed line denotes the integration limit in the ρ total calculation in Bouwens et al. (2022).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00517v1-figure-3dbfad84-70c6-4d2c-9c4d-a052638bb18d.jpg)

lations is reasonably evident for the brightest objects, they are certain to be included in their respective surveys. This is less evident for the fainter regimes explored. For example, it is useful to ask whether the faintest LAEs included in the LAE LF are adequately represented by the faint end of the UV luminosity function. In order to assess this, we can look at the faintest LAEs and compare them to the UV magnitude regimes included in the UV LF of Bouwens et al. (2022).

It is by now well established that certain LAEs remain undetected in even the deepest photometric searches to date (see, for example, Bacon et al. 2015; Maseda et al. 2018; Sobral &amp; Matthee 2019; Maseda et al. 2020; de La Vieuville et al. 2020; Goovaerts et al. 2023). Additionally, some LAEs display faint continuum that would not pass the criteria for many UV-based selections, be they based on colour-colour cuts (Stark et al. 2010; Pentericci et al. 2011, 2018) or photometric redshifts (Caruana et al. 2018; Kusakabe et al. 2020; Goovaerts et al. 2023). We note that preliminary searches conducted by JWST NIRCam in / lensing clusters yields almost exactly the same number of undetected LAEs (Goovaerts et al. 2024). On the other hand, nebular lines are detected in these systems, so they are likely to be real detections (Maseda et al. 2023). With future, deep JWST (and HST) surveys, we will be able to shed more light on this issue.

This raises concerns as to whether the UV LF truly represents the entirety of the galaxy population, or whether there are a significant number of 'UV-dark' LAEs that are counted in the LAE LF but not in the UV LF. This is important to consider given the high values of fesc Ly , α we derive when integrating to 0 0006 . L ∗ (Fig. 5).

A way that we can assess which of these two cases is closer to reality is to compare the UV magnitude regime where these HST-dark LAEs are expected to lie to the regime covered by the UV LF in Bouwens et al. (2022) and the integration limits we use. Upper limits on the UV magnitude of these LAEs can be calculated by manually extracting the continuum level from the

HST images, as in Goovaerts et al. (2023). These objects, denoted by upward-pointing triangles to indicate upper estimates on MUV, are displayed in Fig. 6. This shows that for the clusters considered by the UV LF in question, and for the photometry in question, the UV-dark LAEs would be expected to appear around MUV = -16, with many residing between this and MUV = -14, and some down to MUV ∼ -12.

The UV LF in Bouwens et al. (2022), inclusive of the correction for selection incompleteness, extends past MUV = -16 at all redshifts relevant for this work, and extends past MUV = -14 at z ∼ 3 and z ∼ 6, and close to it at z ∼ 4 and z ∼ 5. At z ∼ 3, the coverage extends all the way to MUV = -12. The integration limit (0 0006 . L ∗ ) used for the calculation of ρ total from the UV LF corresponds to MUV = -13, as denoted by the horizontal dashed line in Fig. 6.

There are four LAEs estimated to have upper UV magnitude limits fainter than this in the sample of four lensing clusters in Goovaerts et al. (2023). As the non-detection of these LAEs likely depends at least in part on the relative depths of the spectroscopy and photometry, it is di ffi cult to estimate what e ff ect this has across our sample, but we can reasonably expect this to have a minor e ff ect, considering the small number of LAEs in this position in comparison to the sample size used for the two LFs.

We can additionally reassure ourselves by comparing what the UV luminosities should be from the faintest LAEs assuming di ff erent levels of Lymanα escape. The faintest LAEs included in the LAE LF are no fainter than LLy α = 10 39 erg s -1 , which equates to an absolute UV magnitude of MUV ∼ -10 4 . assuming 100% Lymanα escape (meaning this value is a lower (faint) limit on the absolute UV magnitude). If we assume a more likely fesc Ly , α value of our sample median from Section 2 of 22.5%, this UV magnitude becomes MUV = -12 2. This is . very faint compared to the reach of the UV LF we use, below the integration limit of MUV = -13. However, using this integration limit and an escape fraction of 22.5%, we recover a value LLy α = 2 4 . × 10 39 erg s -1 . There is only one LAE in our sample close to this luminosity (see Figure 3 of Thai et al. 2023). On the other hand, following this analysis shows that faint LAEs with LLy α &lt; 10 40 erg s -1 are not accounted for in the UV LF with its integration limit set at 0 0006 . L ∗ as in Fig. 5 if they have fesc Ly , α &gt; 95%. Across the board, this e ff ect will bias our results for f esc Ly , α low.

We can perform the same exercise in an inverse manner, to assess whether we are likely to see the Lymanα emission from the faintest UV objects included in the UV LF. The LAE luminosity expected from a galaxy of MUV = -12 is 10 39 6 . erg s -1 , assuming 100% Lymanα escape. Assuming 22.5% escape, this value will be 10 39 0 . erg s -1 . So for galaxies of this absolute magnitude, we can expect to account for all the LAEs in the population (when taking into account the completeness correction). However the escape fraction can be lower (down to 0%), meaning that we will miss the objects of very low UV magnitude and very low Lymanα escape. Unlike missing objects in the UV selection, this biases our results towards high Lymanα escape values. This counter-bias can alleviate the bias introduced by missing UV-faint objects as described in the previous paragraph, and indeed this is the motivation behind choosing one value of L ⋆ down to which the two LF integrations are calculated.

In summary, and in light of the high values for fesc Ly , α found in Fig. 5, for these results to be overestimated, a large population of continuum-undetected LAEs must be playing a role (together with the e ff ects introduced by their incompleteness) in the LAE LF that are not accounted for in the UV LF, even with the com-

pleteness treatment described in Bouwens et al. (2022). This potential e ff ect is important for studies looking to evaluate the total SFRD of the Universe.

Alternatively, if we assume that this e ff ect is negligible and the high fesc Ly , α value derived at z ∼ 6 when including our whole faint sample is reflecting a real similarity between SFRDLy α and SFRDUV, this indicates that galaxies at z ≳ 6 are emitting significantly higher amounts of Lymanα photons than later in the Universe. This supports a picture of galaxies towards reionisation as chemically unevolved and highly star forming, creating a lot of Lymanα emission and allowing it to escape.

## 4.4. Revisiting the contribution of LAEs to cosmic reionisation

Seeing as the escape of Lymanα photons is a crucial quantity for the assessment of the LAE contribution to cosmic reionisation, we can now revisit this in light of our results. There is significant debate in the literature about the relative contribution of the LAE population to reionisation, (see, for example Cassata et al. 2011; Drake et al. 2017; de La Vieuville et al. 2019; Matthee et al. 2022; Thai et al. 2023). All of these results are, naturally, dependent on the fesc Ly , α values assumed in the various studies. The f esc Ly , α value used to calculate this contribution is often pivotal in the classification of this contribution. We can now recalculate the LAE contribution taking into account our results.

We started by considering the ionising emissivity from SFGs, using the following formalism (e.g. Madau et al. 1999; Robertson et al. 2013; Duncan &amp; Conselice 2015):

$$\dot { n } _ { i o n } ( z ) & = \rho _ { U V } ( z ) \xi _ { i o n } \, \mathfrak { f } _ { \text{esc}, \text{LyC} }, & ( 4 ) & \quad A \text{ v}$$

where ˙ nion is expressed as the product of the UV luminosity density ρ UV ( z ) (Bouwens et al. 2015b; Finkelstein 2016; Oesch et al. 2018), global UV luminosity to ionising photon conversion factor ξ ion (Matthee et al. 2016a; Lam et al. 2019), and the ionising photon escape fraction fesc LyC. We reformulated Eq. 4, as we are , interested in the quantity ˙ nion ( z ) for the LAE population only:

$$i _ { i o n, L y \alpha } ( z ) & = \rho _ { L y \alpha } ( z ) \, \xi _ { i o n } ^ { L y \alpha } \, f _ { e s c, L y C }, & ( 5 ) \quad \text{value to} \\ & \quad \text{rived by}$$

where ρ Ly α is the Lymanα luminosity density, ξ Ly α ion is the conversion factor of the observed Lymanα luminosity to the number of ionising photons, and fesc LyC , is the Lyman continuum escape fraction of the LAEs. Ignoring collisionally excited emission and using Case B recombination (Brocklehurst 1971; Osterbrock &amp; Ferland 2006; see also Matthee et al. 2022) as well as the intrinsic flux ratio between Lymanα and H α of 8.7, the ionising emissivity of LAEs, corrected for fesc Ly , α , is written as:

$$\dot { n } _ { i o n, L y \alpha } ( z ) & = \frac { \rho _ { L y \alpha } ( z ) \times f _ { \text{esc,LyC} } } { 8. 7 \, c _ { H \alpha } \times f _ { \text{esc,Ly\alpha} } \, ( 1 - f _ { \text{esc,LyC} } ) }, & ( 6 ) & \quad \text{cases we} \\ & \quad \text{on} \, f _ { \text{esc,Ly} } \\ & \quad \text{not repre}$$

where cH α is the coe ffi cient factor of the H α emission line, ∼ 1 25 . -1 35 . × 10 -12 for Case B recombination. It remains to estimate the individual parameters appearing in the above function.

- - The LAE luminosity density ( ρ Ly α ): this parameter was presented in Sect. 3.3.

- The escape fraction of Lymanα photons, fesc Ly , α . We cannot correct ρ Ly α by the escape fraction derived in Sect. 3.3 as this, by construction, corrects ρ Ly α to the total value of the parent population ( ρ total ) that we take from Bouwens et al. (2022). We can nonetheless use the value we have derived in Sect. 2. We discuss both our sample median, fesc Ly , α = 22 5% and the sample . median for the LAEs at z &gt; 5, f esc Ly , α = 36%.

- LyC escape fraction, fesc LyC: This parameter is poorly un-, derstood, in particular at high redshift, due to the increasing IGM neutrality, which absorbs Lyman continuum photons along the light of sight. However, recent progress has been made in local galaxies (high-redshift analogues) (Flury et al. 2022; Izotov et al. 2024) and in LAEs up to z ∼ 3 (Liu et al. 2023; Kerutt et al. 2024; Pahl et al. 2024; Jung et al. 2024). However, results from these studies, as well as simulations (Choustikov et al. 2024a), suggest that there is a significant scatter in the fesc LyC , across the LAE population. Therefore, placing one value on fesc LyC, such , as fesc Ly , α , is a simplification. Many observational studies point to values between 10% and 15%: Kerutt et al. (2024) derive an underlying population value of 12%, Jung et al. (2024) obtain values between 3% and 15%, Liu et al. (2023) give an upper limit of 16% and Izotov et al. (2024) find values between 1% and 34%, with most of their galaxies having fesc LyC , &lt; 15%. Based on the SPHINX simulations (Garel et al. 2021) and the model for Lyman continuum escape developed by Rosdahl et al. (2022), fesc LyC , is 1 17% 0 4% at MUV . , . = -17 , -13 (we note that this is for all galaxies, not just those selected as LAEs).

The ability to turn Eq. 6 into a summation, with each LAE contributing its own ρ Ly α , f esc Ly , α and with fesc LyC , consequently following a relation such as those presented in Rosdahl et al. (2022), Choustikov et al. (2024a) or Izotov et al. (2024) would be an improvement. Increased statistics, especially for faint LAEs and at redshifts greater than 5 would be necessary to make this approach feasible. This is within reach for JWST in the coming years, as observed samples increase in size and more faint galaxies are uncovered.

A very recent approach by Choustikov et al. (2024b) attempts this individual approach for a large sample of JWST NIRCam-observed galaxies (not just LAEs) (observed by / the JADES programme, Eisenstein et al. 2023). The authors conclude that faint galaxies (MUV ≥ -18 5) are potentially respon-. sible for ∼ 80% of the reionisation budget. The majority of our sample of LAEs fall within this bracket (see, e.g. Fig. 6).

With ˙ nion Ly , α , at z = 6, from Eq. 6, we can compare this value to the ionising emissivity required for reionisation as derived by Bouwens et al. (2015a), ˙ nion crit , = 50 92 (we note that . this value uses a fiducial clumping factor C HII = 3). In Table 2, we show the range of results, applying di ff erent values for ρ Ly α , f esc Ly , α and fesc LyC. For , ρ Ly α derived from our results, (the LF of Thai et al. 2023) we find LAE contributions ranging from minor (5-10%) for a low fesc LyC , (0 4%) to dominant . &gt; 100% for f esc LyC , = 10%. Changing fesc Ly , α between the sample median and z &gt; 5 median does not qualitatively change the outcome. As fesc LyC , changes more significantly depending on the chosen value, this represents the main e ff ect on the outcome. In several cases we see unrealistic contributions ( &gt; 100%) for LAEs based on fesc LyC , = 10%, suggesting that this value may be too high and not representative of our whole sample. Crucial work remains to correlate fesc LyC , with observable parameters (such as Lymanα parameters) at higher redshifts and constrain the escape of ionising photons for the LAE population at z ∼ 6 and above.

We investigate the e ff ect on the LAE reionisation contribution in the case of a completeness cut in the calculation of the LAE LF at 10%. We take the results of de La Vieuville et al. (2019) as emblematic of this scenario (as mentioned, the LF of Thai et al. (2023) is consistent with that of de La Vieuville et al. (2019) if this cut is applied). In this case, the contribution to reionisation derived for the LAE population is significantly less (1 -2% for fesc LyC , = 0 4% and 30 . -50% for fesc LyC , = 10%). This emphasises the sensitivity of these results to the completeness treatment applied for the LAE population.

Table 2. LAE contribution to reionisation.

| Key             | ρ Ly α (10 39 )[erg s - 1 Mpc - 3 ]   |   f esc , Ly α [%] | Contribution [%] (f esc , LyC = 0 . 4%)   | Contribution [%] (f esc , LyC = 10%)   |
|-----------------|---------------------------------------|--------------------|-------------------------------------------|----------------------------------------|
| T23, median     | 47 . 31 ± 37 . 63                     |               22.5 | 9 + 16 - 7                                | 260 + 210 - 50                         |
| T23, z > 5      | 47 . 31 ± 37 . 63                     |               36   | 6 + 10 - 4                                | 160 + 130 - 30                         |
| dLV19, median   | 9 . 52 ± 5 . 65                       |               22.5 | 1 . 9 + 1 . 1 - 0 . 8                     | 52 + 31 - 21                           |
| dLV19, z > 5    | 9 . 52 ± 5 . 65                       |               36   | 1 . 2 + 0 . 7 - 0 . 4                     | 52 + 31 - 21                           |
| T23 GSF, median | 103 . 2 ± 76 . 5                      |               22.5 | 20 + 15 - 5                               | 560 + 420 - 150                        |

Notes. The key indicates the provenance of the values used in columns two and three, in the order they appear in the table. GSF indicates where the general Schechter function was used rather than the modified Schechter function. Values of ρ Ly α come from the LAE LF integrated to 0 0006 . L ∗ (log LLy α ∼ 39). T23, dLV19 refer to Thai et al. (2023) and de La Vieuville et al. (2019) respectively.

We additionally show the di ff erence between the modified Schechter function, which we have chosen to adopt throughout this study, and the general Schechter function (noted GSF in Table 2. Using the general Schechter function roughly doubles the LAE contribution to reionisation.

detections to play a role here, but only for the faintest objects, MUV ≳ -14.

Highlighting the e ff ects of both the function chosen to fit the LAE LF and the completeness cut applied serves to reinforce the need for a larger sample at z &gt; 5 and for faint galaxies. Better statistics in this parameter space will constrain the faint end of the LF more strongly, discriminating between the modified and general Schechter functions. With deeper observations of faint galaxies, the completeness for the faint LAEs will also improve (i.e. the sample will be more complete at a given Lymanα luminosity) and the chosen cut will become a less important factor.

The chosen value of fesc LyC , is clearly determinant in this calculation with LAEs contributing very little for small values of f esc LyC , and providing the dominant contribution for reasonable values of fesc LyC , observed in many lower-redshift LAEs and lowmass galaxies acting as high-redshift analogs. Using the ρ Ly α from Thai et al. (2023) and fesc Ly , α = 36%, a sample average f esc LyC , of just 6 5% would be enough for LAEs to provide all . the necessary ionising photons at z = 6. If the values of fesc LyC , observed in these low-redshift LAEs are maintained across the LAE population into reionisation, it is certainly conceivable that these galaxies could have been responsible for reionisation in its entirety. We note that the results in this study do not constrain the reionisation scenario beyond z = 6 7 and that LAE LFs at . z ≥ 7 are needed to track the contribution of LAEs further into the EoR.

## 5. Conclusions

We have derived the escape fraction of Lymanα photons across the redshift range 2 9 . &lt; z &lt; 6 7 using two methods, namely, . comparing dust-corrected star formation to star formation derived from Lymanα emission and comparing the SFRD from the UV 'parent' population LF and the Lymanα SFRD derived from the LAE LF. Here we summarise the main conclusions of the work.

- - We determined the values of f esc Ly , α for 96 individual galaxies behind the A2744 lensing cluster, finding a significant scatter and a sample median of 22 5 %. .
- -We observed a significant evolution of f esc Ly , α with absolute UV magnitude, following the relation log(fesc Ly , α ) = (0 28 . ± 0 03) MUV . + (4 3 . ± 0 6). . Notably, we saw very few LAEs with fesc Ly , α &lt; 10 % among objects fainter than MUV = -18.
- - We expect the line-incompleteness of MUSE Lymanα
- - For the global redshift evolution of f esc Ly , α , we derived results in agreement with previous determinations in the literature when using a bright LF integration limit (0 04 . L ∗ ). f esc Ly , α appears to reach 100% around z = 11 -12.
- - When we assessed fesc Ly , α with an integration limit that reflects the faint limits of our sample (0 0006 . L ∗ ), we found that f esc Ly , α becomes consistent with 100% escape by z ∼ 6. This result supports galaxies being chemically unevolved and highly star forming at these redshifts, producing a high quantity of Lymanα photons and allowing them to escape.
- - With both bright and faint integration limits, we observed little evolution between z = 3 -5 and a jump at z ∼ 6, suggesting a rapid evolution of the LAE population towards reionisation, although this e ff ect is at &lt; 1 σ for our dataset.
- - We saw more significant redshift evolution in our global assessment than for individual galaxies due to the large intrinsic scatter in the LAE population combined with the far greater statistics in the respective LFs used for the global comparison.
- -We revisited the LAE contribution to the reionisation process, taking into account our new values for fesc Ly , α and a range of motivated values for fesc LyC. The results depend heavily , on the assumed values for f esc LyC , . Taking the latest values observed in LAEs at z ∼ 3, we find that LAEs can provide all the ionising emissivity needed for reionisation.
- - We demonstrated the additional strong dependence of the LAE contribution on the completeness cuts and fitting form of the LAE LF. This stresses the urgent need for better statistics at z &gt; 5 and log LLy α [erg / s] &lt; 40.

At the moment, it is challenging to improve the two luminosity functions used in this work, as they already use the full extent of deep observations in lensing clusters from HFF and MUSE. Further clusters to be included in this sample are desirable, particularly to better constrain the high-redshift ( z ∼ 6) faint end of the LAE LF. This will be made possible through further MUSE-observed lensing clusters, such as those presented in Claeyssens et al. (2024). This will also help bin the results into smaller redshift bins to better constrain the precise evolution in the range 3 &lt; z &lt; 6. Observations of LAEs at higher redshift from JWST, such as Chen et al. (2024); Tang et al. (2024); Lin et al. (2024); Witstok et al. (2024a); Nakane et al. (2024) are starting to build larger sample sizes, but they are not yet at the quantity of LAEs and UV-selected galaxies used for the results in this work. In particular, intrinsically faint galaxies have poor statistics, although this can be improved with the continuing observations of lensing clusters by JWST NIRCam / and JWST NIRSpec. /

On the other hand, it is possible, with further JWST observations, photometry, and spectroscopy in order to improve the fesc Ly , α determinations for individual galaxies. The Balmer decrement provides accurate estimations of dust attenuation along with SED fitting utilising a greater number of photometric bands and a more extended coverage. Studying these systems in more detail will reveal the range of di ff erent escape scenarios leading to the varied fesc Ly , α values that we have observed. This can, in turn, lead to better understanding of galaxy evolution around the epoch of reionisation and how LAEs a ff ect this process.

Acknowledgements. This work is done based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme IDs 060.A-9345, 092.A-0472, 094.A-0115, 095.A-0181, 096.A-0710, 097.A0269, 100.A-0249, and 294.A-5032. Also based on observations obtained with the NASA ESA / Hubble Space Telescope, retrieved from the Mikulski Archive for Space Telescopes (MAST) at the Space Telescope Science Institute (STScI). STScI is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555. All plots in this paper were created using Matplotlib (Hunter 2007). Part of this work was supported by the French CNRS, the Aix-Marseille University, the French Programme National de Cosmologie et Galaxies (PNCG) of CNRS INSU with INP and IN2P3, co-/ funded by CEA and CNES. This work also received support from the French government under the France 2030 investment plan, as part of the Excellence Initiative of Aix-Marseille University - A*MIDEX (AMX-19-IET-008 - IPhU). Financial support from the World Laboratory, the Odon Vallet Foundation and VNSC is gratefully acknowledged. Tran Thi Thai was funded by Vingroup JSC and supported by the Master, PhD Scholarship Programme of Vingroup Innovation Foundation (VINIF), Institute of Big Data, code VINIF.2023.TS.108. This research was funded by Vingroup Innovation Foundation under project code VINIF.2023.DA.057.

## References

Acebron, A., Jullo, E., Limousin, M., et al. 2017, MNRAS, 470, 1809

- Arrabal Haro, P., Rodríguez Espinosa, J., Muñoz-Tuñón, C., et al. 2018, Monthly Notices of the Royal Astronomical Society, 478, 3740
- Atek, H., Kunth, D., Hayes, M., Östlin, G., &amp; Mas-Hesse, J. M. 2008, Astronomy &amp;Astrophysics, 488, 491
- Atek, H., Labbé, I., Furtak, L. J., et al. 2024, Nature, 626, 975
- Austin, D., Conselice, C. J., Adams, N. J., et al. 2024, arXiv e-prints, arXiv:2404.10751
- Bacon, R., Brinchmann, J., Conseil, S., et al. 2023, Astronomy and Astrophysics, 670, A4
- Bacon, R., Brinchmann, J., Richard, J., et al. 2015, Astronomy &amp; Astrophysics, 575, A75
- Beauchesne, B., Clément, B., Hibon, P., et al. 2023, Monthly Notices of the Royal Astronomical Society, 527, 3246
- Bertin, E. &amp; Arnouts, S. 1996, Astronomy and astrophysics supplement series, 117, 393
- Bezanson, R., Labbe, I., Whitaker, K. E., et al. 2022, arXiv e-prints, arXiv:2212.04026
- Blaizot, J., Garel, T., Verhamme, A., et al. 2023, MNRAS, 523, 3749
- Blanc, G. A., Adams, J. J., Gebhardt, K., et al. 2011, The Astrophysical Journal, 736, 31
- Boquien, M., Burgarella, D., Roehlly, Y., et al. 2019, Astronomy &amp; Astrophysics, 622, A103
- Bouwens, R. J., Aravena, M., Decarli, R., et al. 2016, The Astrophysical Journal, 833, 72
- Bouwens, R. J., Illingworth, G., Ellis, R. S., Oesch, P., &amp; Stefanon, M. 2022, ApJ, 940, 55
- Bouwens, R. J., Illingworth, G., Oesch, P., et al. 2015a, The Astrophysical Journal, 803, 34
- Bouwens, R. J., Illingworth, G. D., Franx, M., et al. 2009, ApJ, 705, 936
- Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2015b, The Astrophysical Journal, 811, 140
- Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2010, ApJ, 709, L133 Brocklehurst, M. 1971, MNRAS, 153, 471
- Bruzual, G. &amp; Charlot, S. 2003, Monthly Notices of the Royal Astronomical Society, 344, 1000
- Buat, V., Noll, S., Burgarella, D., et al. 2012, Astronomy &amp; Astrophysics, 545, A141
- Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, The Astrophysical Journal, 533, 682
- Caruana, J., Wisotzki, L., Herenz, E. C., et al. 2018, Monthly Notices of the Royal Astronomical Society, 473, 30
- Cassata, P., Le Fèvre, O., Garilli, B., et al. 2011, Astronomy &amp; Astrophysics, 525, A143
- Chabrier, G. 2003, Publications of the Astronomical Society of the Pacific, 115, 763
- Chen, Z., Stark, D. P., Mason, C., et al. 2024, MNRAS, 528, 7052
- Choustikov, N., Katz, H., Saxena, A., et al. 2024a, MNRAS[ arXiv:2401.09557 ]
- Choustikov, N., Stiskalek, R., Saxena, A., et al. 2024b, arXiv e-prints, arXiv:2405.09720
- Ciardullo, R., Zeimann, G. R., Gronwall, C., et al. 2014, The Astrophysical Journal, 796, 64
- Claeyssens, A., Richard, J., Blaizot, J., et al. 2022, Astronomy and Astrophysics, 666, A78
- Claeyssens, A., Richard, J., &amp; Runnholm, A. 2024, Astronomy &amp; Astrophysics [ submitted ]
- Covone, G., Kneib, J. P., Soucail, G., et al. 2006, A&amp;A, 456, 409
- Cowie, L. L. &amp; Hu, E. M. 1998, The Astronomical Journal, 115, 1319
- De Barros, S., Pentericci, L., Vanzella, E., et al. 2017, Astronomy &amp; Astrophysics, 608, A123
- de La Vieuville, G., Bina, D., Pello, R., et al. 2019, Astronomy &amp; Astrophysics, 628, A3
- de La Vieuville, G., Pelló, R., Richard, J., et al. 2020, Astronomy &amp; Astrophysics, 644, A39
- Dijkstra, M. 2016, Understanding the Epoch of Cosmic Reionization: Challenges and Progress, 145
- Dijkstra, M., Gronke, M., &amp; Venkatesan, A. 2016, ApJ, 828, 71
- Dijkstra, M. &amp; Kramer, R. 2012, Monthly Notices of the Royal Astronomical Society, 424, 1672
- Drake, A.-B., Garel, T., Wisotzki, L., et al. 2017, Astronomy &amp; Astrophysics, 608, A6
- Duncan, K. &amp; Conselice, C. J. 2015, Monthly Notices of the Royal Astronomical Society, 451, 2030
- Eisenstein, D. J., Willott, C., Alberts, S., et al. 2023, arXiv e-prints, arXiv:2306.02465
- Ferrara, A. 2024, A&amp;A, 684, A207
- Finkelstein, S. L. 2016, Publications of the Astronomical Society of Australia, 33, e037
- Finkelstein, S. L., Malhotra, S., Rhoads, J. E., Hathi, N. P., &amp; Pirzkal, N. 2009, Monthly Notices of the Royal Astronomical Society, 393, 1174
- Flury, S. R., Jaskot, A. E., Ferguson, H. C., et al. 2022, ApJS, 260, 1
- Furtak, L. J., Atek, H., Lehnert, M. D., Chevallard, J., &amp; Charlot, S. 2021, MNRAS, 501, 1568
- Garel, T., Blaizot, J., Rosdahl, J., et al. 2021, Monthly Notices of the Royal Astronomical Society, 504, 1902-1926
- Gazagnes, S., Cullen, F., Mauerhofer, V., et al. 2024, ApJ, 969, 50
- Giovinazzo, E., Trebitsch, M., Mauerhofer, V., Dayal, P., &amp; Oesch, P. A. 2024 [ Arxiv:2402.17635v1 ]
- Gnedin, N. Y. 2016, ApJ, 825, L17
- Goovaerts, I., Pello, R., Burgarella, D., et al. 2024, A&amp;A, 683, A184
- Goovaerts, I., Pello, R., Thai, T. T., et al. 2023, A&amp;A, 678, A174
- Gronke, M. &amp; Dijkstra, M. 2016, The Astrophysical Journal, 826, 14
- Hansen, M. &amp; Oh, S. P. 2006, Monthly Notices of the Royal Astronomical Society, 367, 979
- Hashimoto, T., Verhamme, A., Ouchi, M., et al. 2015, The Astrophysical Journal, 812, 157
- Hayes, M., Östlin, G., Schaerer, D., et al. 2010, Nature, 464, 562
- Hayes, M., Östlin, G., Schaerer, D., et al. 2013a, The Astrophysical Journal Letters, 765, L27
- Hayes, M., Östlin, G., Schaerer, D., et al. 2013b, The Astrophysical Journal Letters, 765, L27
- Hayes, M., Schaerer, D., Östlin, G., et al. 2011, The Astrophysical Journal, 730, 8
- Hu, E. M., Cowie, L. L., &amp; McMahon, R. G. 1998, The Astrophysical Journal, 502, L99
- Iani, E., Caputi, K. I., Rinaldi, P., et al. 2023, arXiv e-prints, arXiv:2309.08515 Izotov, Y. I., Chisholm, J., Worseck, G., et al. 2022, MNRAS, 515, 2864
- Izotov, Y. I., Thuan, T. X., Guseva, N. G., et al. 2024, MNRAS, 527, 281
- Jaacks, J., Thompson, R., &amp; Nagamine, K. 2013, The Astrophysical Journal, 766, 94
- Jiang, H., Wang, X., Cheng, C., et al. 2023, arXiv e-prints, arXiv:2312.04151 Jung, I., Ferguson, H. C., Hayes, M. J., et al. 2024 [ Arxiv:2403.02388v1 ] Jung, I., Finkelstein, S. L., Arrabal Haro, P., et al. 2024, ApJ, 967, 73 Kennicutt, R. C. 1998a, The astrophysical journal, 498, 541 Kennicutt, R. C. 1998b, Annual Review of Astronomy and Astrophysics, 36, 189 Kerutt, J., Oesch, P. A., Wisotzki, L., et al. 2024, A&amp;A, 684, A42 Konno, A., Ouchi, M., Nakajima, K., et al. 2016, The Astrophysical Journal, 823, 20

Article number, page 14 of 16

Kusakabe, H., Blaizot, J., Garel, T., et al. 2020, Astronomy &amp; Astrophysics, 638, A12

Lagattuta, D. J., Richard, J., Bauer, F. E., et al. 2019, Monthly Notices of the

Royal Astronomical Society Lam, Bouwens, Rychard J., Labbé, Ivo, et al. 2019, A&amp;A, 627, A164 Limousin, M., Richard, J., Jullo, E., et al. 2016, A&amp;A, 588, A99 Lin, X., Cai, Z., Wu, Y., et al. 2024, ApJS, 272, 33

- Liu, Y., Jiang, L., Windhorst, R. A., Guo, Y., &amp; Zheng, Z.-Y. 2023, ApJ, 958, 22 Lotz, J. e., Koekemoer, A., Coe, D., et al. 2017, The Astrophysical Journal, 837, 97
- Madau, P. &amp; Dickinson, M. 2014, Annual Review of Astronomy and Astrophysics, 52, 415-486

Madau, P., Haardt, F., &amp; Rees, M. J. 1999, The Astrophysical Journal, 514, 648 Mahler, G., Richard, J., Clément, B., et al. 2018, Monthly Notices of the Royal Astronomical Society, 473, 663

- Małek, K., Buat, V., Roehlly, Y., et al. 2018, Astronomy &amp; Astrophysics, 620, A50
- Maseda, M. V., Bacon, R., Franx, M., et al. 2018, ApJ, 865, L1
- Maseda, M. V., Bacon, R., Lam, D., et al. 2020, MNRAS, 493, 5120
- Maseda, M. V., Lewis, Z., Matthee, J., et al. 2023, ApJ, 956, 11
- Matthee, J., Naidu, R. P., Pezzulli, G., et al. 2022, Monthly Notices of the Royal Astronomical Society, 512, 5960
- Matthee, J., Sobral, D., Best, P., et al. 2016a, Monthly Notices of the Royal Astronomical Society, 465, 3637
- Matthee, J., Sobral, D., Hayes, M., et al. 2021, Monthly Notices of the Royal Astronomical Society, 505, 1382
- Matthee, J., Sobral, D., Oteo, I., et al. 2016b, Monthly Notices of the Royal Astronomical Society, 458, 449
- McClymont, W., Tacchella, S., D'Eugenio, F., et al. 2024, arXiv e-prints, arXiv:2405.15859
- McGreer, I. D., Fan, X., Jiang, L., &amp; Cai, Z. 2018, The Astronomical Journal, 155, 131

Naidu, R. P., Oesch, P. A., van Dokkum, P., et al. 2022, ApJ, 940, L14 Nakane, M., Ouchi, M., Nakajima, K., et al. 2024, ApJ, 967, 28 Nanayakkara, T., Glazebrook, K., Jacobs, C., et al. 2023, ApJ, 947, L26 Napolitano, L., Pentericci, L., Santini, P., et al. 2024 [ Arxiv:2402.11220v1 ] Oesch, P. A., Bouwens, R. J., Illingworth, G. D., Labbé, I., &amp; Stefanon, M. 2018,

The Astrophysical Journal, 855, 105

- Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and active galactic nuclei
- Osterbrock, D. E. &amp; Ferland, G. J. 2006, Astrophysics Of Gas Nebulae and Active Galactic Nuclei (University science books)
- Ouchi, M., Shimasaku, K., Akiyama, M., et al. 2008, ApJS, 176, 301
- Ouchi, M., Shimasaku, K., Okamura, S., et al. 2004, The Astrophysical Journal, 611, 660
- Oyarzún, G. A., Blanc, G. A., González, V., Mateo, M., &amp; Bailey, J. I. 2017, The Astrophysical Journal, 843, 133
- Pahl, A. J., Shapley, A. E., Steidel, C. C., et al. 2024 [ Arxiv:2401.09526v1 ]
- Parsa, S., Dunlop, J. S., McLure, R. J., &amp; Mortlock, A. 2016, MNRAS, 456, 3194
- Pello, R., Le Borgne, J.-F., Soucail, G., Mellier, Y., &amp; Sanahuja, B. 1991, ApJ, 366, 405

Pentericci, L., Fontana, A., Vanzella, E., et al. 2011, The Astrophysical Journal, 743, 132

Pentericci, L., Vanzella, E., Castellano, M., et al. 2018, Astronomy &amp; Astrophysics, 619, A147

Piqueras, L., Conseil, S., Shepherd, M., et al. 2019, ASP Conf. Ser. Vol. 521, Astronomical Data Analysis Software and Systems XXVI

Planck Collaboration, Aghanim, N., Akrami, Y., et al. 2020, A&amp;A, 641, A6

Postman, M., Coe, D., Benítez, N., et al. 2012, ApJS, 199, 25

Reddy, N. A., Oesch, P. A., Bouwens, R. J., et al. 2018, The Astrophysical Journal, 853, 56

- Reddy, N. A., Steidel, C. C., Pettini, M., et al. 2008, ApJS, 175, 48
- Richard, J., Claeyssens, A., Lagattuta, D., et al. 2021, Astronomy &amp; Astrophysics, 646, A83

Richard, J., Smith, G. P., Kneib, J.-P., et al. 2010, Monthly Notices of the Royal Astronomical Society, 404, 325

- Robertson, B. E., Furlanetto, S. R., Schneider, E., et al. 2013, The Astrophysical Journal, 768, 71

Rosdahl, J., Blaizot, J., Katz, H., et al. 2022, Monthly Notices of the Royal Astronomical Society, 515, 2386-2414

Runnholm, A., Hayes, M., Melinder, J., et al. 2020, ApJ, 892, 48

Salpeter, E. E. 1955, The Astrophysical Journal, 121, 161

- Saxena, A., Bunker, A. J., Jones, G. C., et al. 2024, A&amp;A, 684, A84
- Schaerer, D., de Barros, S., &amp; Stark, D. P. 2011, Astronomy &amp; Astrophysics, 536, A72
- Schaerer, D. &amp; Verhamme, A. 2008, Astronomy &amp; Astrophysics, 480, 369 Shipley, H. V., Lange-Vagle, D., Marchesini, D., et al. 2018, The Astrophysical Journal Supplement Series, 235, 14
- Simmonds, C., Tacchella, S., Maseda, M., et al. 2023, Monthly Notices of the Royal Astronomical Society, 523, 5468-5486
- Sobral, D. &amp; Matthee, J. 2019, Astronomy &amp; Astrophysics, 623, A157
- Sobral, D., Matthee, J., Best, P., et al. 2017, Monthly Notices of the Royal Astronomical Society, 466, 1242
- Sobral, D., Santos, S., Matthee, J., et al. 2018, MNRAS, 476, 4725
- Song, M., Finkelstein, S. L., Gebhardt, K., et al. 2014, The Astrophysical Journal, 791, 3
- Stark, D. P., Ellis, R. S., Chiu, K., Ouchi, M., &amp; Bunker, A. 2010, Monthly Notices of the Royal Astronomical Society, 408, 1628
- Stark, D. P., Ellis, R. S., &amp; Ouchi, M. 2011, The Astrophysical Journal Letters, 728, L2

Steidel, C. C., Bogosavljevi´ c, M., Shapley, A. E., et al. 2018, ApJ, 869, 123 Tang, M., Stark, D. P., Ellis, R. S., et al. 2024, MNRAS, 531, 2701 Thai, T., Tuan-Anh, P., Pello, R., et al. 2023, A&amp;A, 678, A139 Theulé, P., Burgarella, D., Buat, V., et al. 2024, A&amp;A, 682, A119 Topping, M. W., Stark, D. P., Endsley, R., et al. 2022, ApJ, 941, 153 Trainor, R. F., Strom, A. L., Steidel, C. C., &amp; Rudie, G. C. 2016, The Astrophysical Journal, 832, 171

- Treu, T., Schmidt, K., Brammer, G., et al. 2015, The Astrophysical Journal, 812, 114
- Verhamme, Dubois, Y., Blaizot, J., et al. 2012, A&amp;A, 546, A111
- Verhamme, A., Orlitová, I., Schaerer, D., et al. 2017, Astronomy &amp; Astrophysics, 597, A13
- Verhamme, A., Schaerer, D., Atek, H., &amp; Tapken, C. 2008, Astronomy &amp; Astrophysics, 491, 89
- Vitte, E., Verhamme, A., &amp; Hibon, P. 2024, Astronomy &amp; Astrophysics [ submitted ]
- Weaver, J. R., Cutler, S. E., Pan, R., et al. 2023, arXiv e-prints, arXiv:2301.02671 Weilbacher, P. M., Palsa, R., Streicher, O., et al. 2020, Astronomy &amp; Astrophysics, 641, A28
- Wilkins, S. M., Bunker, A., Coulton, W., et al. 2013, Monthly Notices of the Royal Astronomical Society, 430, 2885
- Witstok, J., Maiolino, R., Smit, R., et al. 2024a, arXiv e-prints, arXiv:2404.05724 Witstok, J., Smit, R., Saxena, A., et al. 2024b, A&amp;A, 682, A40
- Witten, C., Laporte, N., Martin-Alvarez, S., et al. 2024, Nature Astronomy, 8, 384
- Yang, H., Malhotra, S., Gronke, M., et al. 2017, The Astrophysical Journal, 844, 171
- Yuan, Y., Martin-Alvarez, S., Haehnelt, M. G., Garel, T., &amp; Sijacki, D. 2024, MNRAS
- Yue, B., Ferrara, A., &amp; Xu, Y. 2016, MNRAS, 463, 1968
- Zheng, Z., Cen, R., Trac, H., &amp; Miralda-Escudé, J. 2010, The Astrophysical Journal, 716, 574
- Zheng, Z.-Y., Finkelstein, S. L., Finkelstein, K., et al. 2013, Monthly Notices of the Royal Astronomical Society, 431, 3589

## Appendix A: LAE luminosity density and its star-forming rate density

Table A1. LAE luminosity density and its SFRD.

| Redshift           | Reference                                        | α              | log L ∗   |   log 0 . 0006 L ∗ | ρ Ly α (10 39 )   | LAE SFRD (10 - 3 )   | f escp [%]          |
|--------------------|--------------------------------------------------|----------------|-----------|--------------------|-------------------|----------------------|---------------------|
| z = 3 . 1          | Ouchi et al. (2008)                              | - 1 . 5 (fix)  | 42.76     |              39.54 | 9 . 21 ± 2 . 89   | 8 . 37 ± 2 . 63      | 5 . 21 ± 1 . 85     |
| z = 3 . 7          | Ouchi et al. (2008)                              | - 1 . 5(fix)   | 43.0      |              39.78 | 5 . 99 ± 2 . 06   | 5 . 44 ± 1 . 87      | 3 . 66 ± 1 . 44     |
| z = 5 . 7          | Ouchi et al. (2008)                              | - 1 . 5 (fix)  | 42.83     |              39.61 | 9 . 04 ± 6 . 11   | 8 . 21 ± 5 . 55      | 15 . 06 ± 10 . 39   |
| z = 1 . 95 - 3     | Cassata et al. (2011)                            | - 1 . 6        | 42.74     |              39.52 | 7 . 45 ± 3 . 15   | 6 . 76 ± 2 . 86      | 3 . 57 ± 1 . 63     |
| z = 3 - 4 . 55     | Cassata et al. (2011)                            | - 1 . 78       | 42.83     |              39.61 | 7 . 86 ± 2 . 60   | 7 . 14 ± 2 . 36      | 4 . 86 ± 1 . 87     |
| z = 4 . 55 - 6 . 6 | Cassata et al. (2011)                            | - 1 . 69 (fix) | 43 . 0    |              39.78 | 12 . 42 ± 3 . 10  | 11 . 28 ± 2 . 82     | 18 . 5 ± 5 . 4      |
| z = 1 . 9 - 3 . 8  | Blanc et al. (2011)                              | - 1 . 7 (fix)  | 43.08     |              39.86 | 6 . 96 ± 5 . 17   | 6 . 32 ± 4 . 70      | 6 . 42 ± 2 . 95     |
| z = 3 . 1 ± 0 . 3  | Sobral et al. (2018)                             | - 1 . 63       | 42.77     |              39.55 | 11 . 44 ± 6 . 73  | 10 . 39 ± 6 . 11     | 6 . 47 ± 3 . 95     |
| z = 3 . 9 ± 0 . 2  | Sobral et al. (2018)                             | - 2 . 26       | 42.93     |              39.71 | 40 . 56 ± 42 . 59 | 36 . 83 ± 38 . 68    | 25 . 42 ± 27 . 2    |
| z = 4 . 7 ± 0 . 1  | Sobral et al. (2018)                             | - 2 . 35       | 43.28     |              40.06 | 36 . 78 ± 44 . 43 | 33 . 40 ± 40 . 35    | 34 . 46 ± 42 . 14   |
| z = 5 . 4 ± 0 . 4  | Sobral et al. (2018)                             | - 1 . 98       | 43.28     |              40.06 | 17 . 78 ± 11 . 75 | 16 . 14 ± 10 . 67    | 25 . 12 ± 17 . 07   |
| z = 3 . 5          | de La Vieuville et al. (2019)                    | - 1 . 58       | 42.77     |              39.55 | 9 . 80 ± 4 . 85   | 8 . 90 ± 4 . 41      | 5 . 83 ± 3 . 08     |
| z = 4 . 5          | de La Vieuville et al. (2019)                    | - 1 . 72       | 42.96     |              39.74 | 9 . 29 ± 4 . 86   | 8 . 43 ± 4 . 42      | 7 . 66 ± 4 . 29     |
| z = 6 . 0          | de La Vieuville et al. (2019)                    | - 1 . 87       | 43 . 16   |              39.94 | 9 . 52 ± 5 . 65   | 8 . 64 ± 5 . 13      | 15 . 85 ± 9 . 66    |
| z = 3 . 5          | Thai et al. (2023)                               | - 2 . 0        | 42.87     |              39.65 | 33 . 34 ± 15 . 19 | 30 . 27 ± 13 . 79    | 19 . 82 ± 9 . 73    |
| z = 3 . 5          | Thai et al. (2023) (modified Schechter function) | - 2 . 03       | 42.90     |              39.68 | 27 . 69 ± 11 . 24 | 25 . 15 ± 10 . 21    | 16 . 47 ± 7 . 33    |
| z = 4 . 5          | Thai et al. (2023)                               | - 1 . 97       | 42.97     |              39.75 | 23 . 17 ± 11 . 44 | 21 . 04 ± 10 . 39    | 19 . 12 ± 10 . 16   |
| z = 4 . 5          | Thai et al. (2023) (modified Schechter function) | - 2 . 06       | 43.04     |              39.82 | 17 . 47 ± 11 . 01 | 15 . 9 ± 10 . 0      | 14 . 41 ± 9 . 52    |
| z = 6 . 0          | Thai et al. (2023)                               | - 2 . 28       | 43.03     |              39.81 | 103 . 23 ± 76 . 5 | 93 . 74 ± 69 . 47    | 171 . 94 ± 129 . 61 |
| z = 6 . 0          | Thai et al. (2023) (modified Schechter function) | - 2 . 67       | 43.26     |              40.04 | 47 . 31 ± 37 . 63 | 42 . 96 ± 34 . 17    | 78 . 8 ± 63 . 61    |

Notes. The SFRD was calculated by integrating each LF to 0 0006 . L ∗ . The annotation 'fix' denotes slopes that were fixed by the authors of the study in question. The units of ρ Ly α and the LAE SFRD are erg s -1 Mpc 3 -and M yr ⊙ -1 Mpc 3 -respectively and the units of L ∗ are erg s -1 .
Ilias Goovaerts Tran Thi Thai Roser Pello Pham Tuan-Anh Nicolas Laporte Jorryt Matthee Themiya Nanayakkara John Pharo
Charting the Lyman-alpha escape fraction in the range 2.9<z<6.7 and consequences for the LAE reionisation contribution
The escape of Lyman-alpha photons at redshifts greater than two is an ongoing
subject of study and an important quantity to further understanding of
Lyman-alpha emitters (LAEs), the transmission of Lyman-alpha photons through
the interstellar medium and intergalactic medium, and the impact these LAEs
have on cosmic reionisation. This study aims to assess the Lyman-alpha escape
fraction over the redshift range 2.9<z<6.7, focusing on VLT/MUSE-selected,
gravitationally lensed, intrinsically faint LAEs. These galaxies are of
particular interest as the potential drivers of cosmic reionisation. We
assessed the Lyman-alpha escape fraction in two ways: through an individual
study of 96 LAEs behind the A2744 lensing cluster, with JWST/NIRCam and HST
data, and through a study of the global evolution of the escape fraction using
the state-of-the-art luminosity functions for LAEs and the UV-selected `parent'
population (dust-corrected). We compared these studies to those in the
literature based on brighter samples. We find a negligible redshift evolution
of the escape fraction for our individual galaxies; it is likely that it was
washed out by significant intrinsic scatter. We observed a more significant
evolution towards higher escape fractions with decreasing UV magnitude and fit
this relation. When comparing the two luminosity functions to derive the escape
fraction in a global sense, we saw agreement with previous literature when
integrating the luminosity functions to a bright limit. However, when
integrating using a faint limit equivalent to the observational limits of our
samples, we observed enhanced escape fraction values, particularly around z~6,
where the escape fraction becomes consistent with 100%. This indicates for the
faint regimes we sampled that galaxies towards reionisation tend to allow very
large fractions of Lyman-alpha photons to escape. (shortened)