Identifying the mechanisms driving the escape of Lyman Continuum (LyC) photons is crucial to find Lyman Continuum Emitter (LCE) candidates. To understand the physical properties involved in the leakage of LyC photons, we investigate the connection between the HI covering fraction, HI velocity width, the Lyman alpha (LyA) properties and escape of LyC photons in a sample of 22 star-forming galaxies including 13 LCEs. We fit the stellar continua, dust attenuation, and absorption lines between 920 and 1300 A to extract the HI covering fractions and dust attenuation. Additionally, we measure the HI velocity widths of the optically thick Lyman series and derive the LyA equivalent widths (EW), escape fractions (fesc), peak velocities and fluxes at the minimum of the LyA profiles. Overall, we highlight strong correlations between the presence of low HI covering fractions and (1) low LyA peak velocities; (2) more flux at the profile minimum; and (3) larger EW(LyA), fesc(LyA), and fesc(LyC). Hence, low column density channels are crucial ISM ingredients for the leakage of LyC and LyA photons. Additionally, galaxies with narrower HI absorption velocity widths have higher LyA equivalent widths, larger LyA escape fractions, and lower LyA peak velocity separations. This suggests that these galaxies have low HI column density. Finally, we find that dust regulates the amount of LyA and LyC radiation that actually escapes the ISM. Overall, the ISM porosity is one origin of strong LyA emission and enables the escape of ionizing photons in low-z leakers. However, this is not enough to explain the largest fesc(LyC) observed, which indicates that the most extreme LCEs are likely density-bounded along all lines of sight to the observer. Overall, the neutral gas porosity constrains a lower limit to the escape fraction of LyC and LyA photons, providing a key estimator of the leakage of ionizing photons.

A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$\sigma$ upper limit of $\Delta^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power.

We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space.

The processes allowing the escape of ionizing photons from galaxies into the intergalactic medium are poorly known. To understand how Lyman continuum (LyC) photons escape galaxies, we constrain the HI covering fractions and column densities using ultraviolet HI and metal absorption lines of 18 star-forming galaxies which have Lyman series observations. Nine of these galaxies are confirmed LyC emitters. We fit the stellar continuum, dust attenuation, metal, and HI properties to consistently determine the UV attenuation, as well as the column densities and covering factors of neutral hydrogen and metals. We use synthetic interstellar absorption lines to explore the systematics of our measurements. Then we apply our method to the observed UV spectra of low-redshift and z-2 galaxies. The observed HI lines are found to be saturated in all galaxies. An indirect approach using OI column densities and the observed O/H abundances yields HI column densities of 18.6 to 20 cm-2. These columns are too high to allow the escape of ionizing photons. We find that the known LyC leakers have HI covering fractions less than unity. Ionizing photons escape through optically thin holes/channels in a clumpy interstellar medium. Our simulations confirm that the HI covering fractions are accurately recovered. The SiII and HI covering fractions scale linearly, in agreement with observations from stacked Lyman break galaxy spectra at z-3. Thus, with an empirical correction, the SiII absorption lines can also be used to determine the HI coverage. Finally, we show that a consistent fitting of dust attenuation, continuum and absorption lines is required to properly infer the covering fraction of neutral gas and subsequently to infer the escape fraction of ionizing radiation. These measurements can estimate the LyC escape fraction, as we demonstrate in a companion paper.

The fraction of ionizing photons that escape high-redshift galaxies sensitively determines whether galaxies reionized the early universe. However, this escape fraction cannot be measured from high-redshift galaxies because the opacity of the intergalactic medium is large at high redshifts. Without methods to indirectly measure the escape fraction of high-redshift galaxies, it is unlikely that we will know what reionized the universe. Here, we analyze the far-ultraviolet (UV) H I (Lyman series) and low-ionization metal absorption lines of nine low-redshift, confirmed Lyman continuum emitting galaxies. We use the H I covering fractions, column densities, and dust attenuations measured in a companion paper to predict the escape fraction of ionizing photons. We find good agreement between the predicted and observed Lyman continuum escape fractions (within $1.4\sigma$) using both the H I and ISM absorption lines. The ionizing photons escape through holes in the H I, but we show that dust attenuation reduces the fraction of photons that escape galaxies. This means that the average high-redshift galaxy likely emits more ionizing photons than low-redshift galaxies. Two other indirect methods accurately predict the escape fractions: the Ly$\alpha$ escape fraction and the optical [O III]/[O II] flux ratio. We use these indirect methods to predict the escape fraction of a sample of 21 galaxies with rest-frame UV spectra but without Lyman continuum observations. Many of these galaxies have low escape fractions ($f_{\rm esc} \le 1$\%), but 11 have escape fractions $>1$\%. The methods presented here will measure the escape fractions of high-redshift galaxies, enabling future telescopes to determine whether star-forming galaxies reionized the early universe.