🔭 astrophysics

Low-ionization Metal Absorption at $0.7 \lesssim z \lesssim 2$ Confronting Cosmological Simulations with Observations

This paper benchmarks the IllustrisTNG cosmological simulation against observational data of low-ionization metal absorption in quasar spectra, finding that while the simulation broadly reproduces column density distributions, it underestimates the cosmic incidence of MgII systems at higher redshifts unless photo-ionization by the ultraviolet background is included.

Original authors: Ivan Rapoport, Ehud Behar, Vincent Desjacques

Published 2026-04-14

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ivan Rapoport, Ehud Behar, Vincent Desjacques

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

Imagine the universe as a vast, cosmic ocean. Most of the water in this ocean is invisible to our eyes—it's the "dark" gas floating between galaxies. But sometimes, this gas isn't just empty space; it's rich with heavy elements like magnesium and iron, the "minerals" of the cosmos. These elements are the leftovers from dead stars, scattered across the universe.

This paper is like a detective story where scientists are trying to figure out if their computer simulations of the universe are accurate enough to explain what we actually see in the sky.

Here is the breakdown of their investigation, using simple analogies:

1. The Mystery: The "Cosmic Fog"

When light from a distant, bright lighthouse (a quasar) travels through the universe to reach us, it passes through this cosmic fog. As the light passes through clouds of gas, the gas absorbs specific colors of the light, leaving dark "fingerprints" (absorption lines) in the spectrum.

The scientists are looking specifically for low-ionization metals (like Magnesium II). Think of these as the "cool" gas clouds—gas that isn't super hot and ionized, but is relatively cool and dense, sitting right around galaxies or in the space between them.

2. The Simulation: The "Virtual Universe"

The researchers used a massive, super-computer model called IllustrisTNG. Imagine this as a giant, 3D video game of the universe. It simulates gravity, gas physics, and how galaxies form.

However, there's a catch: The computer doesn't naturally "know" how to turn gas into specific chemical fingerprints (like Magnesium II). It just knows the gas exists. So, the scientists had to build a translator. They created two different "rulesets" to guess how the gas behaves:

Rule Set A (The Brawler): Gas particles just crash into each other (collisions) to change their state.
Rule Set B (The Sunlight): Gas particles crash into each other AND get hit by a constant, invisible rain of ultraviolet light from the early universe (the UV Background).

3. The Experiment: Comparing the Game to Reality

The team took their "Virtual Universe" and projected it onto a grid, creating millions of fake "lines of sight" (like looking through a window). They calculated how much Magnesium and Iron should be visible in these fake views.

Then, they compared these fake views to real data from powerful telescopes (like HIRES, UVES, and the new DESI instrument). They looked at two main things:

The "Density Histogram" (PDF): They asked, "How many gas clouds are there with a little bit of metal, a medium amount, or a huge amount?"
- Result: The computer simulation was surprisingly good! It matched the real world almost perfectly for most of the gas clouds. It showed that the basic physics of how gas clumps together is working well in the model.
- The Glitch: The model struggled with the "neutral" (non-ionized) gas, likely because the computer's model for the coldest gas (inside star-forming regions) is a bit too simple.
The "Frequency Count" (Cosmic Incidence): They asked, "How often do we see these gas clouds as we look further back in time (higher redshift)?"
- Result: This is where the models started to stumble.
- For weak, faint clouds (the "drizzle" of gas), the model that included the Ultraviolet Background (Rule Set B) was spot on. It proved that the "sunlight" from the early universe is crucial for keeping these low-density clouds visible.
- For strong, dense clouds (the "storms"), the computer failed. It didn't produce enough of them, especially as we look back to when the universe was younger (around redshift 2). The simulation seems to be missing something that creates these heavy, dense metal clouds.

4. The New Data: The "DESI" Problem

The paper also looked at data from DESI, a new, massive survey that is mapping the universe.

The Analogy: Imagine looking at a forest through a high-powered telescope (HIRES/UVES). You can see every single leaf and branch. Now, imagine looking at the same forest through a foggy window (DESI). You can see the forest is there, but the leaves blur together into big blobs.
The Issue: Because DESI has lower resolution, it tends to "smear" distinct gas clouds together. It sees one big, strong cloud when there might actually be several smaller ones. This makes it hard to compare DESI data directly with the sharp, detailed computer simulations.

The Bottom Line

The scientists concluded that:

The Big Picture is Good: Our current computer models of how galaxies and gas clouds form are doing a great job at the general level.
The "Sunlight" Matters: To get the faint, diffuse gas right, you must include the background ultraviolet light in your model.
The "Heavy" Stuff is Missing: The models are still missing the physics that creates the strongest, densest metal clouds. We don't fully understand how galaxies pump out these heavy elements yet.
Resolution is Key: To fix these models, we need to match the sharpness of our computer simulations with the sharpness of our telescopes. New, sharper data will help us solve the mystery of the missing heavy clouds.

In short: The map is mostly accurate, but we need to zoom in closer to find the missing treasure.

1. Problem Statement

Low-ionization metal absorption lines (specifically Mg I, Mg II, and Fe II) serve as primary probes for the cool, enriched gas in the circumgalactic medium (CGM) and intergalactic medium (IGM). While cosmological simulations like IllustrisTNG successfully reproduce high-ionization species (e.g., O VI, C IV) and hot gas properties, their ability to accurately model low-ionization absorption remains a critical test for galaxy formation physics.

The specific challenges addressed in this work include:

Ionization Physics: Determining whether collisional ionization alone or the inclusion of a cosmic ultraviolet background (UVB) is necessary to reproduce observed column densities.
Statistical Discrepancies: Benchmarking simulation predictions against high-resolution observational data (HIRES/UVES) and large-scale surveys (SDSS, DESI) to see if current models capture the correct incidence rates ($dN/dz$) and column density distributions.
Data Limitations: Assessing the feasibility of using lower-resolution data (DESI) for direct comparison with high-resolution simulations.

2. Methodology

Observational Data

The authors utilize two primary datasets:

HIRES/UVES Catalog (Churchill et al. 2020): A high-resolution ( $R \approx 40,000-50,000$ ) sample of 2,989 kinematic components from 422 Mg II systems ( $0.19 \le z \le 2.55$ ). This provides precise column densities ( $N$ ) and Doppler parameters ( $b$ ).
DESI Catalog (Napolitano et al. 2025): A large, homogeneous sample from the Dark Energy Spectroscopic Instrument (DR1), though with lower spectral resolution ( $R \sim 2000$ ).

Simulation Framework

Simulation: IllustrisTNG50-1, a high-resolution magnetohydrodynamic simulation ( $L=35 h^{-1}$ Mpc) with a baryonic mass resolution of $\sim 8 \times 10^4 M_\odot$ .
Ionization Modeling: Since TNG does not track individual ionization states, the authors implemented a post-processing grid-based framework. They solved the steady-state ionization equilibrium equations for Mg and Fe, considering two models:
- Coll: Purely collisional ionization/recombination.
- Coll+UVB: Collisional processes plus photoionization from a uniform, redshift-dependent cosmic UVB (Faucher-Giguère et al. 2009).
Temperature Handling: Star-forming gas cells were manually set to $T = 10^3$ K to match the TNG two-phase ISM model, as the simulation does not resolve gas below $T \sim 10^4$ K accurately.
Column Density Calculation: Ionic masses were projected onto 2D grids using a Smoothed Particle Hydrodynamics (SPH) kernel to generate synthetic column density maps ( $N_{\text{Mg I}}, N_{\text{Mg II}}, N_{\text{Fe II}}$ ).

Statistical Analysis

Column Density PDFs: Probability Distribution Functions ( $\phi(N)$ ) were constructed for the simulation and compared to observations, normalized over specific ranges to mitigate low-column-density incompleteness biases.
Cosmic Incidence ($dN/dz$): The number of absorbers per unit redshift was calculated for specific Equivalent Width (EW) bins ( $W_{2796}^0$ ) and compared to observational measurements from SDSS and VLT/UVES.

3. Key Contributions

Grid-Based Ionization Post-Processing: The authors developed a computationally efficient method to map TNG gas properties to low-ionization metal column densities using pre-computed ionization grids, avoiding the need for full spectral synthesis in every cell.
Benchmarking Ionization Prescriptions: A rigorous comparison between collisional-only and UVB-included models against high-resolution data.
DESI Feasibility Study: A critical assessment of the DESI Mg II catalog, highlighting the impact of spectral resolution on absorption line blending and component identification.
Redshift Evolution Analysis: An investigation into the evolution of absorption systems between $z \approx 0.7$ and $z \approx 2$ .

4. Key Results

Column Density PDFs

Agreement: The TNG simulation broadly reproduces the observed column density PDFs for Mg II and Fe II across the range $10^{11.4} \lesssim N \lesssim 10^{16}$ cm $^{-2}$ .
Model Divergence:
- At low-to-moderate column densities ( $N \lesssim 10^{14}$ cm $^{-2}$ ), both Coll and Coll+UVB models yield similar results.
- At high column densities ( $N \gtrsim 10^{15.5}$ cm $^{-2}$ ), the Coll+UVB model maintains a continuous slope, while the Coll model shows a sharp break. The Coll model underestimates high-density absorbers because it lacks UVB to sustain Mg II in cooler, denser regions (or conversely, over-ionizes in the absence of UVB shielding logic in the specific temperature regime).
Mg I Discrepancy: TNG fails to reproduce the Mg I PDF slope, predicting a shallower distribution. This is attributed to the TNG multiphase ISM model's inability to accurately resolve gas temperatures below $10^4$ K, where neutral species dominate.
Redshift Independence: The PDFs show little evolution between $z=0.7$ and $z=2$ , suggesting these systems form early and remain stable.

Cosmic Incidence ($dN/dz$)

Low EW Systems ( $W_{2796}^0 < 0.6$ Å): The Coll+UVB model accurately reproduces the observed incidence rate. This confirms that weak absorbers are dominated by low-density, photoionized gas in the outer CGM. The Coll-only model overestimates the incidence of these systems.
High EW Systems ( $W_{2796}^0 > 1$ Å): Both models underestimate the incidence of strong absorbers. TNG fails to capture the observed rise in incidence toward $z \sim 2$ $z \sim 2$ .
- Potential Causes: The simulation may underestimate metal content in cold star-forming regions (where $T=10^3$ K is an oversimplification) or the conversion from column density to EW may be flawed if high-EW systems possess unusually high turbulence ( $b$ -parameters).

DESI vs. HIRES/UVES

Blending Effects: The DESI catalog shows a systematic excess of high-EW systems compared to HIRES/UVES. This is attributed to DESI's lower resolution, which blends multiple narrow kinematic components into single, broader, high-EW features.
Identification Issues: While the total integrated EW along a sightline is generally conserved between the two catalogs, DESI identifies fewer distinct components and exhibits mis-identifications at $z \gtrsim 1.5$ .
Conclusion: Direct comparison of DESI PDFs with TNG is limited by resolution effects; detailed spectral synthesis matching DESI's resolution would be required for robust comparison.

5. Significance and Conclusions

Validation of Physics: The success of the Coll+UVB model in reproducing the incidence of weak absorbers validates the importance of including the cosmic UVB in modeling the outer CGM.
Limitations of Current Simulations: The failure of TNG to reproduce strong absorbers ( $W_{2796}^0 > 1$ Å) and their redshift evolution suggests that current feedback mechanisms or the treatment of cold, dense gas in galaxy centers are insufficient.
Diagnostic Power: The authors conclude that cosmic incidence ($dN/dz$) is a more sensitive diagnostic for testing ionization prescriptions and galaxy formation physics than column density PDFs, which appear largely redshift-agnostic in the current models.
Future Outlook: While grid-based modeling is efficient for broad statistics, future work must combine full spectral synthesis with higher-resolution simulations to resolve the physics of the strongest, most metal-rich absorbers. The study also cautions that lower-resolution surveys like DESI and Euclid NISP require careful handling of line blending before they can be used to constrain cosmological simulations.