🔭 astrophysics

Galaxy Proximate Damped Lyman-Alpha Systems and HI Reionization Topology in TECHNICOLOR DAWN

Using the TECHNICOLOR DAWN cosmological hydrodynamical simulation, this study demonstrates that while circumgalactic media reionize later than the intergalactic medium, the neutral hydrogen column density of proximate damped Lyman- $\alpha$ systems is primarily determined by halo mass, suggesting they can effectively trace reionization progress at high redshifts if stellar or halo masses are precisely estimated.

Original authors: Maya Steen, Kristian Finlator, Samir Kušmić, Ezra Huscher

Published 2026-05-20

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Maya Steen, Kristian Finlator, Samir Kušmić, Ezra Huscher

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: A Cosmic "Fog" Check

Imagine the early universe as a giant, dark room filled with a thick, invisible fog made of neutral hydrogen gas. For a long time, this fog was so dense that light couldn't travel through it easily. Then, the first stars and galaxies turned on like lightbulbs. Their light started to "burn away" the fog, turning the gas from neutral to ionized (a process called reionization).

Astronomers want to know: When exactly did the fog clear up? Was it a slow, uneven process? Did some parts of the room clear before others?

Usually, scientists look at distant quasars (super-bright black holes) to see how much fog is left in the space between galaxies. But quasars are rare. Recently, the James Webb Space Telescope (JWST) found a new way to look: by studying the "fog" right around the galaxies themselves. These are called Proximate Damped Lyman-Alpha Systems (PDLAs). Think of them as looking at the mist swirling directly around a streetlamp, rather than looking at the fog in the distance.

The Experiment: A Cosmic Simulation

The authors didn't just look at the sky; they built a giant, 3D video game simulation of the universe called TECHNICOLOR DAWN. This simulation includes gas, dark matter, and stars, and it runs the laws of physics forward in time to see how the universe evolved from redshift 10 (very early) to 5.5 (a bit later).

They used this simulation to answer two main questions:

Does the fog around a galaxy (the CGM) clear up at the same time as the fog between galaxies (the IGM)?
Can we measure how "foggy" the universe is just by looking at how thick the fog is around a galaxy?

Key Finding 1: The "Inside-Out-Middle" Clearing

The simulation revealed that the universe didn't clear up in a simple way. It followed a pattern the authors call "Inside-Out-Middle."

The Analogy: Imagine a party in a crowded house.
1. Inside: The people right next to the music (the bright galaxies) clear the air first because the light is strongest there.
2. Outside: The light then shoots out through the empty hallways (low-density space) very quickly because there's nothing to block it.
3. Middle: Finally, the light has to work its way back into the crowded rooms (medium-density filaments) to clear the last bits of fog.

The Result: The fog between galaxies (IGM) cleared up first. The fog around the galaxies (CGM) stayed thick and neutral for a long time, even after the rest of the universe was mostly clear. By the time the universe was "halfway" through clearing, the gas right next to galaxies was still very foggy.

Key Finding 2: It's About Size, Not Just "Fogginess"

The researchers wanted to know: If we see a lot of fog (a high column density of hydrogen) around a galaxy, does that mean the whole universe is still foggy?

They found a surprising twist. The amount of fog you see around a galaxy depends mostly on how heavy the galaxy is, not just on how foggy the universe is.

The Analogy: Think of two different-sized buckets.
- Bucket A (Small Galaxy): Even if the air is very foggy, a small bucket can only hold a little bit of water (fog).
- Bucket B (Massive Galaxy): This bucket is huge. Even if the air is only slightly foggy, the bucket is so big that it still catches a massive amount of water.

The Result: A massive galaxy will always have a thick layer of neutral gas around it, even if the universe is mostly clear. A small galaxy might have very little gas around it, even if the universe is still very foggy. Therefore, you can't just look at the fog and guess the age of the universe; you have to know the size (mass) of the galaxy first.

The Conclusion: A New Tool for Astronomers

The paper concludes that galaxy PDLAs can be used to trace the progress of reionization, but with a catch.

If you know the mass of the galaxy (how heavy it is), you can use the amount of fog around it to figure out how much fog is left in the whole universe.

High Redshift (Early Universe): The universe is very foggy. Massive galaxies have huge amounts of fog.
Low Redshift (Later Universe): The universe is clearing up. Massive galaxies still have fog, but less of it.

The Catch: To use this method, you need to know the galaxy's mass very precisely. If you do, this new method is a powerful way to map out the history of the universe's "fog clearing," especially for the very early times when other methods struggle.

Summary in One Sentence

By simulating the early universe, the authors found that the gas around galaxies clears up later than the gas between them, and that the thickness of this gas is mostly a sign of the galaxy's size, meaning we can use it to track the universe's history if we know how big the galaxies are.

Technical Summary: Galaxy Proximate Damped Lyman-Alpha Systems and HI Reionization Topology in TECHNICOLOR DAWN

Problem Statement
Recent observations from the James Webb Space Telescope (JWST) have identified proximate damped Lyman- $\alpha$ systems (PDLAs) in the foreground of high-redshift galaxies ( $z > 5$ ). These systems are interpreted as neutral circumgalactic media (CGM). While PDLAs have historically been used to study quasar environments, their presence around high-redshift galaxies offers a potential new probe for tracing the progress of cosmic reionization. However, the utility of this method is uncertain. It remains unclear whether the ionization history of the CGM mirrors that of the intergalactic medium (IGM) given the complex topology of reionization. Furthermore, it is unknown if the observed column density of neutral hydrogen ( $N_{\text{HI}}$ ) in PDLAs correlates primarily with the global neutral fraction (a proxy for reionization progress) or if it is driven by other factors such as halo mass, gas mass fraction, or local neutral fraction.

Methodology
The authors utilize the TECHNICOLOR DAWN cosmological hydrodynamical simulation to investigate these questions. The simulation employs a $\Lambda$ CDM model with Planck cosmological parameters, utilizing smoothed particle hydrodynamics (SPH) and on-the-fly radiative transfer via the moment method. The simulation box has side lengths of 16.5 cMpc/ $h$ and includes $2 \times 704^3$ particles.

The analysis proceeds in two stages:

Topology Analysis: The authors compute the overdensity ( $\Delta$ ) and neutral hydrogen fraction ( $x_{\text{HI}}$ ) for individual gas particles across redshift snapshots ( $z = 10, 8, 6, 5.5$ ). They analyze the ratio of mass-averaged to volume-averaged neutral fractions to characterize the reionization topology.
Halo-by-Halo Analysis: Galaxies are identified using the SKID algorithm, and their host dark matter halos are defined via spherical overdensity. The study focuses on resolved halos with masses $M > 10^{8.5} M_\odot$ . The authors compute the column density of neutral hydrogen ( $N_{\text{HI}}$ ) in the foreground of these halos, weighting the calculation by the local star formation rate (SFR) to mimic observational sightlines through star-forming regions. They then correlate $N_{\text{HI}}$ with halo mass, stellar mass, gas mass fraction ( $f_g$ ), and mass-weighted local neutral fraction ( $\langle f_{\text{HI}} \rangle$ ).

Finally, the authors fit a functional relationship between $\log_{10}(N_{\text{HI}})$ , halo/stellar mass ( $M$ ), and the volume-averaged neutral fraction ( $\langle x_{\text{HI}} \rangle$ ) to determine the feasibility of using PDLAs to infer reionization progress.

Key Results

Reionization Topology: The simulation reveals an "inside-out-middle" (IOM) reionization topology. Initially, reionization proceeds in an inside-out manner, ionizing dense regions immediately surrounding sources. This is followed by a dominant outside-in phase where ionizing flux escapes to low-density regions, ionizing them rapidly due to low recombination rates. Finally, the ionization fronts propagate back toward intermediate-density regions (filaments). Consequently, the CGM reionizes after the IGM and remains partially neutral even at $z = 5.5$ .
Drivers of Column Density: The study finds that the column density of neutral hydrogen ( $N_{\text{HI}}$ ) is not primarily driven by the gas mass fraction or the local neutral fraction. Instead, $N_{\text{HI}}$ shows a clear, dominant correlation with halo mass (and stellar mass). Massive halos exhibit higher $N_{\text{HI}}$ even when their local neutral fractions are low, simply because they contain more total gas mass.
Correlation with Reionization Progress: Despite the dominance of mass, a statistically significant relationship exists between $N_{\text{HI}}$ and the global neutral fraction. The authors establish the following functional form:
$\log_{10}(N_{\text{HI}}) = a + b \log_{10}(M) + c \langle x_{\text{HI}} \rangle$
where $b$ (mass dependence) is the dominant term, but $c$ (neutral fraction dependence) is non-zero.
Observational Feasibility: The authors perform an error propagation analysis to estimate the number of galaxies required to constrain $\langle x_{\text{HI}} \rangle$ to a specific fractional precision. They find that the method is most sensitive at high redshifts where $\langle x_{\text{HI}} \rangle$ is large. As reionization nears completion and $\langle x_{\text{HI}} \rangle$ drops, the required number of observations scales as $\langle x_{\text{HI}} \rangle^{-2}$ , making the method less sensitive at lower redshifts.

Significance and Claims
The paper claims that galaxy PDLAs represent a viable, new observational method for tracing the progress of cosmic reionization, provided that precise estimates of halo or stellar mass are available to disentangle the mass dependence from the ionization state.

The authors emphasize that their findings support the interpretation that high column densities in PDLAs are not merely indicators of high gas fractions or locally neutral environments, but are strongly tied to the host galaxy's mass. By accounting for this mass dependence, PDLAs can serve as a tracer for the global neutral fraction. The study suggests this method is particularly well-suited for constraining early reionization (high $z$ ) using current and projected JWST data, though it notes that sightline-to-sightline variations and limited simulation box sizes may introduce scatter that requires further investigation. The results align with previous simulations (e.g., Finlator et al. 2009) regarding IOM topology and offer consistency with observational trends reported in recent literature (e.g., Heintz et al. 2024; Umeda et al. 2024).