pith. sign in

arxiv: 2604.14676 · v1 · submitted 2026-04-16 · 🌌 astro-ph.CO

Tracking ionization balance in intergalactic medium and its implications towards metallicity

Bhaskar Arya , Kartick C. Sarkar , Shiv K. Sethi This is my paper

Pith reviewed 2026-05-10 10:21 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords intergalactic mediumionization balancemetal ionsC IVreionizationultraviolet backgroundnon-equilibrium chemistry
0
0 comments X

The pith

A zero-dimensional model evolves the ionization and temperature of intergalactic gas including 107 metal ions while matching full simulations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a simplified framework that follows a single parcel of gas as the universe expands, tracking how hydrogen, helium, and many metal ions change their ionization states over time under a changing ultraviolet background. This matters because quasar light absorbed by the IGM reveals ionic columns that depend on the gas density, temperature, metallicity, and radiation history, and those timescales overlap with cosmic evolution. The framework solves the coupled rate equations and temperature evolution self-consistently, then checks the results against three-dimensional hydrodynamical runs. It reproduces the overall thermal and ionization paths, including the extra heating when helium becomes doubly ionized. As a demonstration, the model calculates the cosmic density of C IV ions and converts observed values into estimates of IGM metallicity that align with existing constraints.

Core claim

The authors present a fast, metals-inclusive, zero-dimensional framework for modeling the redshift evolution of the IGM that solves stiff time-dependent rate equations for H, He, and 107 metal ions while self-consistently evolving the temperature through photoheating and cooling; the framework reproduces the thermal and ionization histories from full three-dimensional hydrodynamical non-equilibrium calculations over a wide redshift range, including He II reionization heating, and predicts the cosmic C IV density parameter to infer metallicities consistent with observations.

What carries the argument

The zero-dimensional Lagrangian gas parcel that solves time-dependent ionization rate equations for hydrogen, helium, and 107 metal ions while evolving temperature from photoheating and cooling.

Load-bearing premise

A single average gas parcel evolving in a uniform radiation field can represent the overall thermal and ionization behavior of the entire intergalactic medium despite real variations in density and radiation across different regions.

What would settle it

A mismatch larger than reported between the model's predicted ionization fractions or C IV density and either new high-precision quasar absorption measurements or independent three-dimensional simulations that include spatial inhomogeneities.

Figures

Figures reproduced from arXiv: 2604.14676 by Bhaskar Arya, Kartick C. Sarkar, Shiv K. Sethi.

Figure 1
Figure 1. Figure 1: Evolution of𝑇0 w.r.t.redshift for our zero dimension NEI simulation (shown in red curve) compared to full 3D NEI hydro simulation (H/He NEI only; shown in black dashed line). Overlaid are observational constraints on 𝑇0 as shown in colored circles (Schaye et al. 2000; Becker et al. 2011; Bolton et al. 2012; Boera et al. 2014). We find that the zero-dimensional model reproduces the 3D NEI thermal history to… view at source ↗
Figure 3
Figure 3. Figure 3: Evolution of ion fractions of H i and He ii w.r.t. redshift at mean density, for our zero dimension NEI simulation (shown in red curves) com￾pared to full 3D NEI hydro simulation (shown in black curves). We also show He i fraction for our model (dotted red). We again see that the zero￾dimensional model recovers the ion fractions in very good agreement with 3D NEI simulations. that the simplified one-zone t… view at source ↗
Figure 4
Figure 4. Figure 4: Quadratic fit to the observed C iv density parameter, log ΩCIV, as a function of redshift. Coloured symbols show individual observational constraints: Yu et al. (2025) (teal circle), Davies et al. (2023) (red square), Codoreanu et al. (2018) (blue down triangle), Boksenberg & Sargent (2015) (magenta right triangle), Ryan-Weber et al. (2009) (yellow cross), Simcoe (2011) (cyan up triangle), D’Odorico et al.… view at source ↗
Figure 5
Figure 5. Figure 5: Top: C iv density parameter, ΩCIV, obtained using our 0D framework in the IGM using full NEI thermochemistry with Illustris-TNG weights (red solid) and Sherwood weights (red dashed). The black dashed curve represents the empirical fit for various observations of C iv density parameter as mentioned in fig. 4. The gray shaded region denotes the 1-𝜎 uncertainty. Bottom panel: Metallicities inferred for the ob… view at source ↗
Figure 7
Figure 7. Figure 7: C iv density parameter for three choices of maximum overdensity: (i) the full case including gas up to 𝛿 ∼ 1000 (solid), (ii) restricted to 𝛿 ≲ 600 (dashed), and (iii) restricted to 𝛿 ≲ 300 (dotted). a factor of ∼ 15. Thus, Civ-bearing gas would only become self￾shielded at 𝛿 ∼ 3000. Gas at such high densities, however, is unlikely to host a significant fraction of C iv, since rapid recombination drives ca… view at source ↗
Figure 8
Figure 8. Figure 8: Normalized C iv abundance (see text for definition) in our 0D model w.r.t. overdensity at given redshifts. Black vertical line denotes virial overdensity for a Lagrangian density point (refer to AppendixBfor calculation and assumptions). spectra, typically outside obvious DLAs and sub-DLAs. Yet they re￾cover essentially the full ΩC IV quoted in the literature from 𝑧 ≳ 1.5 up to the reionization epoch. This… view at source ↗
Figure 9
Figure 9. Figure 9: Left: Colormap of the ratio between our ionization correction and the standard PIE-based correction used in observational metallicity estimates, (𝑥CIV/𝑥HI)0D/(𝑥CIV/𝑥HI)PIE, shown in the 𝑧–Δ phase space. The dotted black curves trace the trajectories of a representative set of density points used in our analysis. The gray dashed curves show overdensities containing 25% and 75% of global C iv content, showin… view at source ↗
read the original abstract

Ionization balance in the intergalactic medium (IGM) is central to the interpretation of quasar absorption spectra, linking observed ionic columns to the underlying gas density, temperature, metallicity, and ionizing radiation field. Because ionization, recombination, and cooling timescales can be comparable to the timescales over which the ultraviolet background (UVB) and gas thermodynamic state evolve, ion populations may retain a strong memory of their past history. To this end, we present a fast, metals-inclusive, zero-dimensional framework for modeling the redshift evolution of the IGM. The model follows the coupled thermal and ionization evolution of a Lagrangian gas parcel in a redshift-dependent UVB, solving stiff, time-dependent rate equations for H, He, and 107 metal ions while self-consistently evolving the temperature through photoheating and standard cooling processes. We validate the framework against full three-dimensional hydrodynamical non-equilibrium calculations and find that it reproduces the thermal and ionization histories of the IGM with good accuracy over a wide redshift range, including the heating associated with $\rm He_{\,\rm II}$ reionization. As an application, we predict the cosmic $\rm C_{\,\rm IV}$ density parameter, $\Omega_{\rm CIV}$, and use it to infer the origin of metal ions in the IGM and the corresponding metallicities from observational measurements, obtaining values broadly consistent with literature constraints. The framework is well suited for rapid parameter studies of how reionization timing, UVB spectral hardness, self-shielding, and UVB inhomogeneity shape the thermal and ionization history of the IGM and the resulting metal-line observables.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper presents a zero-dimensional Lagrangian framework that evolves the coupled thermal and ionization state of a gas parcel in a redshift-dependent UVB, solving stiff rate equations for H, He, and 107 metal ions while self-consistently computing temperature via photoheating and cooling. It validates the model against full 3D non-equilibrium hydrodynamical simulations, claiming good accuracy in reproducing IGM thermal and ionization histories (including He II reionization heating) over a wide redshift range. As an application, the framework is used to predict the cosmic C IV density parameter Ω_CIV and to infer IGM metallicities from observational data, obtaining values broadly consistent with literature constraints.

Significance. If the validation holds quantitatively, the framework provides a fast, metals-inclusive tool for parameter studies of reionization timing, UVB spectral shape, self-shielding, and inhomogeneity effects on IGM observables. This is useful for interpreting quasar absorption spectra and constraining metal enrichment histories. The self-consistent treatment of 107 metal ions and temperature evolution, together with explicit comparison to independent 3D calculations, are notable strengths.

major comments (2)
  1. [Validation section] Validation section: The manuscript asserts that the 0D model 'reproduces the thermal and ionization histories of the IGM with good accuracy' including He II reionization heating, but does not report quantitative metrics (e.g., fractional RMS error in T(z), x_HII(z), or C IV fraction versus redshift) comparing the single-parcel trajectory to the mass- or volume-weighted averages from the 3D simulations. Without these, it is difficult to assess whether the accuracy is sufficient to support the subsequent Ω_CIV prediction.
  2. [Application to Ω_CIV] § on application to Ω_CIV: The central claim that the 0D model enables reliable inference of metallicities from observed Ω_CIV rests on the assumption that a single mean-density parcel with uniform UVB captures the volume-averaged C IV fraction. The paper should demonstrate (or quantify the error from) how well the 0D C IV fraction matches the 3D distribution when local density variations and inhomogeneous UVB/self-shielding are present, as these are known to affect metal-ion fractions.
minor comments (2)
  1. [Methods] Notation for metal ions (e.g., C IV vs. C^{3+}) should be defined consistently in the methods section and used uniformly in figures and tables.
  2. [Introduction] The abstract and introduction would benefit from a brief statement of the specific UVB model (e.g., which spectral shape and normalization) adopted for the fiducial run.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight important aspects of the validation and application sections that we address below. We have revised the manuscript to incorporate quantitative metrics and additional discussion as outlined in our point-by-point responses.

read point-by-point responses

  1. Referee: [Validation section] Validation section: The manuscript asserts that the 0D model 'reproduces the thermal and ionization histories of the IGM with good accuracy' including He II reionization heating, but does not report quantitative metrics (e.g., fractional RMS error in T(z), x_HII(z), or C IV fraction versus redshift) comparing the single-parcel trajectory to the mass- or volume-weighted averages from the 3D simulations. Without these, it is difficult to assess whether the accuracy is sufficient to support the subsequent Ω_CIV prediction.

    Authors: We agree that quantitative metrics would allow a more rigorous evaluation of the claimed accuracy. In the revised manuscript we have added a new subsection (and accompanying table) that reports the fractional RMS errors between the 0D trajectories and the mass- and volume-weighted averages extracted from the 3D simulations for temperature, x_HII, x_HeII, and the C IV fraction over z = 2–6. These metrics show typical deviations of order 10 % or less for the key quantities, thereby substantiating the statement of good accuracy and the subsequent use of the model for Ω_CIV predictions. revision: yes

  2. Referee: [Application to Ω_CIV] § on application to Ω_CIV: The central claim that the 0D model enables reliable inference of metallicities from observed Ω_CIV rests on the assumption that a single mean-density parcel with uniform UVB captures the volume-averaged C IV fraction. The paper should demonstrate (or quantify the error from) how well the 0D C IV fraction matches the 3D distribution when local density variations and inhomogeneous UVB/self-shielding are present, as these are known to affect metal-ion fractions.

    Authors: We acknowledge that the 0D framework adopts a mean-density parcel with a uniform UVB and that local density variations and UVB inhomogeneities can in principle affect ion fractions. Our existing validation already compares the 0D results directly to volume-averaged quantities from the 3D runs, which incorporate these effects. In the revised manuscript we have added a dedicated paragraph that quantifies the residual uncertainty by examining the scatter in C IV fractions across density bins in the 3D simulations; we find that the mean-density approximation recovers the volume-averaged C IV fraction to within ~15–20 % over the relevant redshift range. This additional analysis supports the robustness of the metallicity inferences while making the approximation explicit. revision: yes

Circularity Check

0 steps flagged

No significant circularity: 0D IGM model validated externally and applied to independent observables

full rationale

The paper defines a zero-dimensional Lagrangian parcel that solves time-dependent rate equations for H, He, and 107 metal ions coupled to temperature evolution under a redshift-dependent UVB. Validation consists of direct numerical comparison to separate three-dimensional hydrodynamical non-equilibrium calculations, reproducing thermal and ionization histories including He II reionization heating. The derived Ω_CIV is obtained from the model's ion fractions and then compared to separate observational data for metallicity inference, without any parameter fitting to the target quantities or reduction of the central result to self-citations, self-defined quantities, or ansatzes imported from prior author work. All load-bearing steps remain independent of the final predictions.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The framework rests on standard atomic physics rates and cosmological assumptions rather than new postulates. No new particles or forces are introduced.

free parameters (1)
  • UVB spectral shape and normalization parameters
    The redshift-dependent UVB is an external input whose detailed functional form and amplitude are chosen or fitted from literature models; these choices affect the ionization and heating rates.
axioms (2)
  • domain assumption Atomic ionization, recombination, and cooling rate coefficients from standard databases are accurate and complete for the included ions
    The stiff rate equations for 107 metal ions rely on these pre-tabulated coefficients without re-derivation in the paper.
  • domain assumption A single Lagrangian gas parcel evolving in a uniform UVB adequately represents the volume-averaged IGM thermal and ionization state
    This zero-dimensional approximation is the core modeling choice and is validated only by comparison to 3D simulations.

pith-pipeline@v0.9.0 · 5591 in / 1681 out tokens · 36716 ms · 2026-05-10T10:21:40.058255+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    AguirreA.,SchayeJ.,TheunsT.,2002,TheAstrophysicalJournal,576,1–20 AguirreA.,Dow-HygelundC.,SchayeJ.,TheunsT.,2008,TheAstrophysical Journal, 689, 851–864 Banerjee E., Muzahid S., Schaye J., Johnson S. D., Cantalupo S., 2023, Monthly Notices of the Royal Astronomical Society, 524, 5148 BeckerG.D.,BoltonJ.S.,HaehneltM.G.,SargentW.L.W.,2011,MNRAS, 410, 1096 B...

    work page doi:10.1007/978-3- 2002

  2. [2]

    B1 NFW profile and the virial density For an NFW halo (Navarro et al

    We define the virial radius of a halo at redshift𝑧 by requiring that the mean density inside𝑅vir is a fixed overdensity Δvir times some reference density𝜌ref (𝑧), usually either the critical density𝜌 crit (𝑧)or the mean background matter density¯𝜌m (𝑧): ¯𝜌(< 𝑅vir) ≡ 3𝑀vir 4𝜋𝑅 3 vir = Δvir 𝜌ref (𝑧).(B1) Foragivenhalomass𝑀 vir andchoiceofΔ vir,thisimmediate...

    work page 1997