Joakim Rosdahl, Jérémy Blaizot, Harley Katz, Taysun Kimm, Thibault Garel, Martin Haehnelt, Laura C. Keating, Sergio Martin-Alvarez, Léo Michel-Dansac, Pierre Ocvirk
5 min
This study investigates the escape fraction of ionizing Lyman-continuum (LyC) radiation from galaxies during the Epoch of Reionization (EoR). Using the SPHINX suite of cosmological radiation-hydrodynamical simulations, the authors aim to understand how radiation escapes the interstellar medium (ISM) of high-redshift galaxies and enters the intergalactic medium (IGM). By resolving halos down to with a 10 pc resolution, the researchers track the interplay between star formation, stellar feedback, and radiative transfer to determine which galaxies were responsible for reionizing the Universe.
The simulations reveal that the escape fraction () is not a static property but fluctuates significantly in individual galaxies over timescales of a few million years. This variability is driven by supernova and radiation feedback, which periodically clears gas from the ISM, creating temporary channels for ionizing radiation to escape. Statistically, is highest in intermediate-mass, low-metallicity galaxies (). The study finds that a large portion of the ionizing budget—approximately 55%—is provided by galaxies that are currently too dim to be detected by existing observational surveys. Furthermore, the global average decreases over cosmic time as specific star formation rates decline.
Understanding the escape of LyC photons is critical for modeling the reionization of the Universe. Because direct observations of LyC escape are nearly impossible at high redshifts due to IGM absorption, researchers must rely on simulations to constrain the ionizing budget. This study provides a robust statistical framework for identifying the primary drivers of reionization, suggesting that low-mass, metal-poor galaxies are the dominant contributors, rather than the brighter, more massive galaxies often targeted in observational campaigns.
We measure escape fractions, $f_{\rm esc}$, of ionizing radiation from galaxies in the SPHINX suite of cosmological radiation-hydrodynamical simulations of reionization, resolving halos with $M_{\rm vir} \gtrapprox 7.5 \times 10^7 \ M_{\odot}$ with a minimum cell width of $\approx 10$ pc. Our new and largest $20$ co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times $10^{11} \ M_{\odot}$. The simulated galaxies agree well with observational constraints of the UV luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over timescales of a few Myrs, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the inter-galactic medium. Statistically, $f_{\rm esc}$ peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies ($M_{*} \approx 10^7 \ M_{\odot}$, $M_{1500} \approx -17$, $Z\lesssim 5 \times 10^{-3} \ Z_{\odot}$), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to SPHINX, galaxies too dim to be yet observed, with $M_{1500} \gtrapprox -17$, provide about $55$ percent of the photons contributing to reionization. The global averaged $f_{\rm esc}$ naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time.
Sam: But simulating every single atom in the universe would be impossible, right? So how do you manage the scale?
Alex: That's where a technique called adaptive mesh refinement comes in. Imagine drawing a grid over your simulation. In the empty regions of space where not much is happening, the grid squares are large and coarse. But wherever there's a dense cloud of gas or a forming star — wherever the interesting physics is — the grid automatically zooms in and becomes much finer. It lets you focus your computing power exactly where it matters, without wasting it on empty space.
Sam: So you get the fine detail inside the galaxy and the broader cosmic picture at the same time.
Alex: Exactly. You can see the small "windows" opening and closing in the gas, while still tracking how the escaped light spreads through the wider universe.
Sam: So what does the simulation actually show about why reionization wound down?
Alex: The most significant finding is that reionization isn't a smooth, steady process. At any given moment, only a tiny fraction of galaxies are actively punching light through their gas. The rest have their windows closed. So the universe's clearing process looks less like a lamp being slowly turned up, and more like a strobe light — flickering, intermittent, driven by whichever galaxies happen to be in the middle of a burst at that moment.
Sam: A strobe light. So it's not that galaxies fundamentally change — it's that the *timing* of their bursts changes?
Alex: Precisely. What the simulations track is something called the escape fraction — the proportion of ionising light that actually makes it out of a galaxy rather than being absorbed by the surrounding gas. That fraction is closely tied to how bursty a galaxy's star formation is. As galaxies age and their star formation becomes more stable, the escape fraction drops. Not because they've run out of fuel, but because they've stopped being chaotic enough to blast open the channels.
Sam: So the universe didn't dim because galaxies ran out of energy. It dimmed because they grew up.
Alex: That's a good way to put it. The early universe was dominated by small, volatile galaxies — constantly flaring, constantly punching light through their gas. Over time, those galaxies settled. The violence that drove reionization gradually faded, not with a dramatic cutoff, but as the natural consequence of galaxies becoming more organised. It's a meaningful shift in how we understand the process — from thinking about which galaxies are biggest, to asking which ones are still chaotic enough to keep the windows open.
Sam: It makes you think differently about what "growing up" means for a galaxy.
Alex: It does. And it's a reminder that some of the largest transformations in cosmic history were driven not by size or power alone, but by timing — by whether the right kind of chaos was happening in the right place at the right moment. Thanks for listening to ResearchPod.