Pith. Sign in

REVIEW

Non-equilibrium Ionization State of Warm-Hot Intergalactic Medium

Astrophysics · Submitted 3 Mar 2006 · 40 references

Pith's verdict: A solid, honest early study showing WHIM oxygen is out of ionization equilibrium and line ratios are biased, but the headline quantitative claim is conditional on the acknowledged Te=Ti assumption.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/0603072 (version 2)
Pith ID2DEGH43H, Pith's permanent record of this paper and its review
CitationKohji Yoshikawa, Shin Sasaki. "Non-equilibrium Ionization State of Warm-Hot Intergalactic Medium." arXiv:astro-ph/0603072 (2006). https://pith.science/paper/2DEGH43H
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

At z ~ 0, a significant fraction of O VII and O VIII in the warm-hot intergalactic medium is not in ionization equilibrium, making O VIII/O VII line ratios systematically different from equilibrium predictions.

The reading

This paper drops the assumption of ionization equilibrium for the warm-hot intergalactic medium (WHIM)—the diffuse $10^5$–$10^7$ K gas thought to hold a large share of the universe's missing baryons—and follows the time evolution of oxygen ionization along realistic thermal histories. It finds that at $z \simeq 0$ a significant fraction of O VII and O VIII has not reached equilibrium: recombination timescales are comparable to the Hubble time in low-density WHIM plasma, and the photoionizing background has faded rapidly since $z<2$, leaving these ions over-abundant. Consequently, O VIII/O VII emission and absorption line ratios are systematically lower than equilibrium predicts, even though overall detectability through oxygen lines is barely changed. A reader should care because those line ratios are used as a WHIM temperature diagnostic; if the paper is right, that diagnostic is biased.

What carries the argument

The load-bearing object is the coupled rate equation $$\frac{d f_j}{dt} = \sum_{k=1}^{j-1} S_{j-k,k} f_k - \sum_{i=j+1}^{Z+1} S_{i-j,j} f_j - \alpha_j f_j + \alpha_{j+1} f_{j+1},$$ integrated with a backward difference formula—a stable implicit update that advances ionization stages sequentially—along each gas element's thermal history, using the paper's adopted ionization and recombination rates and a redshift-dependent UV/X-ray photoionizing background. The diagnostic carrying the argument is the ratio $R \equiv f_{\rm OVIII}/f_{\rm OVII}$ and its equilibrium counterpart $R_{\rm eq}$; the paper shows $R/R_{\rm eq} < 1$ ubiquitously in filamentary WHIM, which is why line-ratio diagnostics shift while overall line fluxes do not.

Extended reading notes

Core claim

Tracing the ionization balance equation along both idealized shock-heating histories and the actual density–temperature tracks of gas particles in a cosmological smoothed-particle hydrodynamics (SPH) simulation, the paper establishes that at $z=0$ the oxygen ionization state in typical WHIM conditions is not in equilibrium. O VI is close to its equilibrium fraction where its abundance matters, but O VII and O VIII are over-abundant relative to equilibrium, and the ratio $R \equiv f_{\rm OVIII}/f_{\rm OVII}$ is systematically below its equilibrium value in low-overdensity, $10^6$–$10^7$ K gas at the outskirts of filaments and clusters. The deviation is largest for gas shocked recently, around $z \simeq 0.1$, and is amplified by the rapid decline of the UV/X-ray background at $z<2$ because recombination cannot keep pace. Convolving the results into mock observations, the non-equilibrium state leaves O VII and O VIII surface brightnesses and column densities within a factor of order unity of equilibrium predictions, so detectability is not significantly changed; but the ratios $S_{\rm OVIII}/S_{\rm OVII}$ and $N_{\rm OVIII}/N_{\rm OVII}$ are systematically lower, so equilibrium-based temperature estimates are biased.

Load-bearing premise

The calculation assumes electrons and ions are instantly at the same temperature, although cosmological shocks heat ions first and Coulomb relaxation is slow in low-density gas, and the ionization and recombination rates depend on the electron temperature.

Editorial extensions

If this is right

  • Temperature estimates of WHIM gas that assume ionization equilibrium and use O VIII/O VII line ratios will be biased in filament outskirts, because $R/R_{\rm eq}<1$ there shifts the inferred temperature.
  • Detectability forecasts for oxygen emission and absorption in future X-ray surveys remain approximately valid; the main change is in line ratios, not line strengths.
  • O VI remains a usable equilibrium tracer of WHIM, while O VII and O VIII do not, so multi-ion analyses should treat the latter with time-dependent ionization.
  • Line-ratio maps trace recent shock heating in the cosmic web, since the deviation is concentrated at the edges of filaments where shocks occurred around $z \simeq 0.1$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If electron temperature lags ion temperature by a factor of a few—the paper's own caveat—the bias in O VIII/O VII ratios would likely be larger, because the collisional ionization that produces O VIII would be slower; a two-temperature simulation is the natural test.
  • The same timescale argument should apply to other hydrogen- and helium-like metal ions, such as Ne IX, C V, and Fe XVII, so non-equilibrium effects may contaminate other WHIM diagnostics, not just oxygen.
  • Mapping $R/R_{\rm eq}$ across filament edges could act as a chronometer for shock heating: regions with the largest suppression were shocked most recently, complementing temperature measurements.
  • Because the paper neglects galaxy feedback and metal inhomogeneity, and argues its deviations are a lower limit, a future simulation with feedback could predict larger non-equilibrium signatures, giving observers a way to constrain feedback strength from line-ratio data.
Original abstract

Time evolution of the ionization state of metals in the cosmic baryons is investigated in a cosmological context without the assumption of ionization equilibrium. We find that a significant fraction of ionized oxygen ions (OVII and OVIII) in the warm-hot intergalactic medium (WHIM) is not in the ionization equilibrium state at a redshift of z~0. We also investigate the effect on the detectability and observables of WHIM as a consequence of such deviation from ionization equilibrium. It is found that the detectability of WHIM is not altered very much both through its emission and absorption signatures, but line ratios between OVII and OVIII are significantly different from those in the ionization equilibrium state.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 3,712 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.15 · human participants 0.01 · controlled comparison 0.45 · randomized assignment 0.01 · states quantitative result 0.10 · compares with prior results 0.56 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

Time evolution of the ionization state of metals in the cosmic baryons is investigated in a cosmological context without the assumption of ionization equilibrium. We find that a significant fraction of ionized oxygen ions (OVII and OVIII) in the warm-hot intergalactic medium (WHIM) is not in the ionization equilibrium state at a redshift of z~0. We also investigate the effect on the detectability and observables of WHIM as a consequence of such deviation from ionization equilibrium. It is found that the detectability of WHIM is not altered very much both through its emission and absorption signatures, but line ratios between OVII and OVIII are significantly different from those in the ionization equilibrium state.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Time evolution of the ionization state of metals in the cosmic baryons is investigated in a cosmological context without the assumption of ionization equilibrium. contribution 0.96 0.76 0.05 asserted 0.70
s1 We find that a significant fraction of ionized oxygen ions (OVII and OVIII) in the warm-hot intergalactic medium (WHIM) is not in the ionization equilibrium state at a redshift of z~0. result 1.00 0.96 0.84 asserted 1.00
s2 We also investigate the effect on the detectability and observables of WHIM as a consequence of such deviation from ionization equilibrium. contribution 0.99 0.89 0.04 unclear 0.68
s3 It is found that the detectability of WHIM is not altered very much both through its emission and absorption signatures, but line ratios between OVII and OVIII are significantly different from those in the ionization equilibrium state. result 0.99 0.95 0.25 asserted 1.00

What kind of paper the classifier reports this is: simulation study 0.91 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Redshift dependence of the cosmic UV and X-ray background radiation at photon energy of 0.5 keV (left panel) and its spectrum at redshift z = 0, 0.5, 1, and 3 (right panel). In the right panel, we also show the observed background spectra in UV and X-ray bands obtained by Shull et al. (1999) and Miyaji et al. (1998), respectively. 3
Figure 2. Ionization fractions of O vi, O vii and O viii in ionization equilibrium as a function of temperature for the low density plasma with nH = 10−5 and 10−6 cm−3 under the presence of photoionizing background radiation at a redshift of z = 0. Ionization fractions in collisional ionization equilibrium are also shown for the comparison. τion =   Z X +1 i=j+1 Si−j,j   −1 (3) and τrec = α −1 j . (4)
Figure 3. Timescales of ionization and recombination for O vi, O vii and O viii ions in hydrogen number density of nH = 10−5 cm−3 (left) and 10−6 cm−3 (right) under the presence of UV and X-ray background radiation at a redshift of z = 0. Horizontal dotted lines indicate the hubble time H −1 0 (≃ 13 Gyr). Since the timescales for the ionization and recombination differ by many order of mag￾nitude depending on ionic species and ionization stages, equation (1) is a stiff set of differential equations. In numerically solving a stiff set of equations, explicit integration schemes are unstable u
Figure 4. Time evolution of density and temperature of baryons for the model calculations. We have three free parameters δs, Ts and zs. (see text for detail.) At the initial redshift zi = 10, we assume that all baryons are in the ionization equilibrium state. This assumption can be justified by the fact that physical gas density is so high at such a high redshift that ionization equilibrium can be quickly realized. Equation (1) is integrated along the thermal histories described above till z = 0 for a given set of parameters δs , Ts , and zs
Figure 5. Time evolution of temperature (left panel) and hydrogen number density (right panel) of gas particles in a cosmological hydrodynamic simulation presented in section 4. See section 4 for details of the simulation. and O ix on a (1 + δs)–Ts plane. White contour lines are separated by a factor of 1.1 in ratios of ionization fractions, fj/feq j , where the bold lines are contours for fj/feq j = 1. Black contours with numeric labels indicate the ionization fractions in the ionization equilibrium state on (1 +δs)−T plane. One can see that, for O vi ion, the ratios of ionization fractions…
Figure 6. Maps of ratios of ionization fractions between non-equilibrium and equilibrium ionization states for O vi, O vii, O viii and O ix at redshift z = 0 on a (1 + δs)–Ts plane, where shock heating epoch is set to zs = 1.0. White contours are separated by a factor of 1.1 in the ratios of the ionization fractions. Black lines with numeric labels are the contours of the ionization fractions in the ionization equilibrium states. in the initial epoch. Thermal history of each gas particle is interpolated if an time interval between two continuous outputs is longer than the timestep ∆t in…
Figure 7. Time evolution of ionization fractions of O vii, O viii and O ix for four sets of the densities and temperatures at redshift z = 0 marked in figure 6. Bold and thin lines indicate the non-equilibrium and equilibrium states, respectively. non-equilibrium state at z = 0 relative to those in the equilibrium state on a (1 + δb)–T plane, where δb is the baryon overdensity. Black contours in each panel indicate the mass distribution of the corresponding ion, and 25%, 50% and 75% of its mass are enclosed by three contour lines from inside to outside. One can see that most of O vi ions are…
Figure 8. Same as figure 6 except that shock heating epoch is set to zs = 0.5. equilibrium state in a 7.5h −1 thick slice taken from the simulation box at z = 0. In this slice, we have a galaxy cluster at (X, Y ) = (32, 25)h −1 Mpc, and a prominent filamentary structure associated with this galaxy cluster. A close look at figure 11 shows that we have significant de￾viation from the ionization equilibrium state at the outskirts of galaxy clusters and filamentary structures, where cosmological shocks occur (Kang et al. 2005). 5. EFFECT ON DETECTABILITY AND OBSERVABLES OF WHIM The deviations of i
Figure 9. Same as figure 6 except that shock heating epoch is set to zs = 0.1. Note that the ranges of the ratios are different from those in figure 6 and 8. 5.1. Detection through emission Surface brightness of emission lines of O vii at 574eV and O viii at 653eV (SOVII and SOVIII, respectively) are computed for both of the equilibrium and non-equilibrium states. Here, we assume the latest specification of the detector proposed for DIOS mission (Ohashi et al. 2004), namely a spectroscopic energy resolution of 2 eV and a spatial resolution of 5 arcmin
Figure 10. Mean ratios of ionization fractions at z = 0 between non-equilibrium and equilibrium ionization states for O vi, O vii, O viii and O ix on the (1 + δb)–T plane. White contours are drawn in the regions with 0.5 < fj/f eq j < 2 and separated by a factor of 1.2 in the ratios of ionization fractions. The white solid lines indicate the contour for fj/f eq j = 1. Three black contours from inside to outside enclose the 25%, 50% and 75% of mass of each ionic species. in the ionization equilibrium state above the detection limit of DIOS mission, irrespective of redshift of the line…
Figure 11. Spatial distribution of hydrogen number density, and ratios of ionization fractions of O vi, O vii and O viii between non-equilibrium and equilibrium states in the simulation volume at z = 0. The thickness of the slice is 7.5h −1Mpc. Contours of baryon mass density for 1 + δb = 1 are also shown in the top-right, bottom-left and bottom right panels. 5.2. Detection through absorption Column densities of O vi, O vii and O viii (NOVI, NOVII, and NOVIII, respectively) and those in ionization equilibrium (N eq OVI, N eq OVII and N eq OVIII) are calculated for 104 randomly se￾lected…
Figure 12. Ratios of surface brightness of O vii (574eV) and O viii (653eV) emission lines for non-equilibrium states relative to equilibrium states as a function of surface brightness in equilibrium states. Upper and lower panels are for line emissions at redshift ranges of 0.15 < z < 0.3 and 0.0 < z < 0.15, respectively. Symbols and error bars indicate the averages and the standard deviations for each horizontal bin. Vertical dotted lines in both panels indicate the nominal detection limit of DIOS mission. absorbers. On the other hand, we have systematically larger O vii and O viii column…
Figure 13. Same as figure 12, except that ratios relative to surface brightness in the collisional ionization equilibrium state are shown
Figure 14. Ratios of O vi, O vii and O viii column densities of absorbers at redshift 0.0 < z < 0.15 (lower panels) and 0.15<z <0.3 (upper panels) in non-equilibrium states relative to those in ionization equilibrium states identified in the lightcone data as a function of column densities in ionization equilibrium states. Symbols and error bars indicate the averages and standard deviations for corresponding horizontal bins, respectively
Figure 15. Relations of line emission ratio SOVIII/SOVII between non-equilibrium and equilibrium states. Symbols and error bars indicate the averages and standard deviations for corresponding horizontal bins, respectively. 19
Figure 16. Relations of column density ratios NOVII/NOVI (filled squares) and NOVIII/NOVII (open trian￾gles) between non-equilibrium and equilibrium states. Symbols and error bars indicate the averages and standard deviations for corresponding horizontal bins, respectively. exceed 1014 cm−2 . Although NOVII/NOVI in the non-equilibrium state statistically agrees well with that in the ionization equilibrium state in both redshift ranges, the non-equilibrium state exhibits systematically lower NOVIII/NOVII than the ionization equilibrium state. It is found that these differences in NOVIII/NOVII…
Figure 17. Spatial distribution of R/Req at z = 0 in the same slice of the simulation volume as in figure 11. Contours for baryon mass density of 1 + δb = 1 are also shown
Figure 18. Maps of R/Req at z = 0 on a (1 + δb)–T plane in a numerical simulation (left panel) and in a simple model with zs = 0.1 (right panel). White contours are separated by a difference of 0.25 in R/Req and bold lines indicate R/Req = 1. Three black contours in each panel are for Req = 0.1, 1 and 10 from bottom to top. 21

Figures reproduced from arXiv:astro-ph/0603072 by the authors.

References
How the literature uses it

1 arXiv paper cite this in 1 sentence

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.