Where the Forest Goes Dark: Halo-centered Characterization of High Column Density Systems with IllustrisTNG
Abstract
High column density systems of neutral hydrogen (HCDs) --- Lyman limit systems (LLSs), sub-damped and damped Lyα absorbers (DLAs) --- contaminate both the 3D and the 1D statistics of the Lyα forest. Most DLAs are masked out from the analyses, but cleaning algorithms are not perfect, and weaker LLSs are individually undetectable, so the residual contamination is absorbed into nuisance parameters that dilute the cosmological constraining power. Simulations provide the ideal setting to study this effect, allowing absorption features to be directly traced back to the gas producing them. We identify high-column-density structures on the gas cells of the TNG50 hydrodynamical simulation, deblend them in velocity space, and split every sightline into HCD-only and forest-only spectra. Applied to halos at z ≃ 3 out to 50 virial radii, it yields each class's covering fraction against impact parameter b and halo mass. The impact parameter defines the type of absorption appearing in the spectrum: DLAs give way to sub-DLAs at b ≈ 0.11 R200c, sub-DLAs to LLSs at $0.3$, and LLSs to the Lyα-forest at $0.5$, which covers ≃ 78% of sightlines at R200c. In units of R200c the sequence is nearly mass-independent over three decades in mass, so the virial radius sets the scale of the neutral gas distribution. A single Voigt component recovers the column density of essentially every damped system, but not of LLSs, whose absorption features often arise from several separate contributions from gas structures along the sightline. Our detailed characterization of HCDs gives the first step to the construction of advanced techniques to directly forward-model their contribution to Lyα spectra.
Related
- The connection between surface brightness and satellite systems for central galaxies through IllustrisTNG (Corrigendum)
- Uncertainty-Aware Deep Learning for the Lyα Forest: CNN-Based Absorber Detection and Characterization
- Something green beneath the surface: The dynamical nature of fossil structures in IllustrisTNG-100
- The stability of radially anisotropic rotating stellar systems with a central density cusp
- Turbulent Transport of Galactic Magnetic Fields into Cosmic Voids: Insights from IllustrisTNG
BibTeX
@misc{silva2026where,
title = {{Where the Forest Goes Dark: Halo-centered Characterization of High Column Density Systems with IllustrisTNG}},
author = {Victor R. Soares da Silva and Daniel López-Cano and L. Raul Abramo and J. Chaves-Montero and Francisco Maion},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.19235},
eprint = {2609.19235},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.19235},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Where the Forest Goes Dark: Halo-centered Characterization of High Column Density Systems with IllustrisTNG
AU - Victor R. Soares da Silva
AU - Daniel López-Cano
AU - L. Raul Abramo
AU - J. Chaves-Montero
AU - Francisco Maion
PY - 2026
DA - 2026/09/16
PB - arXiv
DO - 10.48550/arxiv.2609.19235
UR - https://arxiv.org/abs/2609.19235
AB - High column density systems of neutral hydrogen (HCDs) --- Lyman limit systems (LLSs), sub-damped and damped Lyα absorbers (DLAs) --- contaminate both the 3D and the 1D statistics of the Lyα forest. Most DLAs are masked out from the analyses, but cleaning algorithms are not perfect, and weaker LLSs are individually undetectable, so the residual contamination is absorbed into nuisance parameters that dilute the cosmological constraining power. Simulations provide the ideal setting to study this effect, allowing absorption features to be directly traced back to the gas producing them. We identify high-column-density structures on the gas cells of the TNG50 hydrodynamical simulation, deblend them in velocity space, and split every sightline into HCD-only and forest-only spectra. Applied to halos at z ≃ 3 out to 50 virial radii, it yields each class's covering fraction against impact parameter b and halo mass. The impact parameter defines the type of absorption appearing in the spectrum: DLAs give way to sub-DLAs at b ≈ 0.11 R_200c, sub-DLAs to LLSs at $0.3$, and LLSs to the Lyα-forest at $0.5$, which covers ≃ 78% of sightlines at R_200c. In units of R_200c the sequence is nearly mass-independent over three decades in mass, so the virial radius sets the scale of the neutral gas distribution. A single Voigt component recovers the column density of essentially every damped system, but not of LLSs, whose absorption features often arise from several separate contributions from gas structures along the sightline. Our detailed characterization of HCDs gives the first step to the construction of advanced techniques to directly forward-model their contribution to Lyα spectra.
ER -