pith. sign in

arxiv: 2509.13690 · v2 · submitted 2025-09-17 · 🌌 astro-ph.CO

Conditioning halos on the tidal environment for fast and accurate HI power spectra during reionization

Gaurav Pundir , Tirthankar Roy Choudhury , Aseem Paranjape This is my paper

Pith reviewed 2026-05-18 16:46 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords reionizationneutral hydrogenHI power spectrumtidal tensorcollapse fractionsemi-numerical models21 cm cosmology
0
0 comments X

The pith

Conditioning halo collapse on tidal eigenvalues recovers large-scale HI power spectra to within 2-5 percent.

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

The paper develops a method to predict the neutral hydrogen distribution during cosmic reionization more accurately and efficiently. It conditions the mean collapse fraction of halos on both the local density and the eigenvalues of the tidal tensor, drawing this conditional mean from small high-resolution simulations. This conditioned field is then used in larger, lower-resolution simulations to generate the HI density and its power spectrum. The resulting large-scale power spectrum matches reference high-resolution runs far better than previous approaches that used density alone. This matters for making reliable forecasts for 21 cm observations that will probe the epoch of reionization.

Core claim

By computing the conditional mean of f_coll given the dark matter overdensity and the eigenvalues of the tidal tensor from high-resolution small-volume simulations, and transferring this to sample f_coll in low-resolution large-volume simulations, the subsequent semi-numerical reionization modeling yields HI power spectra accurate to ≲2%–5% for k ≲ 0.5 h Mpc^{-1} across a range of redshifts, ionized fractions, grid resolutions, and minimum halo masses, substantially improving on the ∼10% accuracy of density-only conditioning.

What carries the argument

Conditional mean of the collapse fraction f_coll on density and tidal tensor eigenvalues, transferred from small high-resolution volumes to large low-resolution ones.

If this is right

  • Large-volume HI power spectra can be computed accurately using only low-resolution density fields augmented with tidal info.
  • The improvement holds across variations in redshift, ionization level, resolution, and halo mass threshold.
  • Large-scale HI power spectra become reliable enough for direct comparison with upcoming 21 cm data.
  • Semi-numerical codes can now incorporate environmental effects beyond density with minimal added cost.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This conditioning technique could be applied to other tracers like galaxies or Lyman-alpha to improve clustering predictions.
  • Generating large ensembles of such HI fields would allow better estimation of covariances for observational analyses.
  • Extending the conditioning to include velocity fields or higher-order statistics might further reduce errors on smaller scales.

Load-bearing premise

The conditional mean of f_coll given density and tidal eigenvalues from small high-resolution volumes transfers directly to large low-resolution volumes without needing volume-dependent corrections.

What would settle it

A high-resolution large-volume simulation in a different cosmological model or at higher redshift where the error exceeds 5% at k=0.5 h/Mpc would falsify the transferability of the conditional means.

read the original abstract

Predicting the statistical properties of the neutral hydrogen (HI) density field during reionization is an important step in using upcoming 21 cm observations to constrain models of reionization. Semi-numerical models of reionization are often coupled with the collapse fraction field $f_{\text{coll}}(\mathbf{x})$, which determines the fraction of dark matter within halos. In this work, we improve upon earlier prescriptions that compute $f_{\text{coll}}$ based on the dark matter overdensity $\delta(\mathbf{x})$ alone, to include more information about the environment in the form of eigenvalues of the tidal tensor. We compute the mean of the $f_{\text{coll}}$ conditioned on these eigenvalues from a set of high-resolution, small-volume simulations and use them to sample the $f_{\text{coll}}$ field of a low-resolution, large-volume simulation. We subsequently use a semi-numerical code for reionization to compute the HI density field and its power spectrum, and benchmark our results against a reference high-resolution, large-volume simulation. Across variations in redshift, ionized fraction, grid resolution, and minimum halo mass, our method recovers the large-scale HI power spectrum with errors at the $\lesssim 2\%-5\%$ level for $k \lesssim 0.5~ h~ \text{Mpc}^{-1}$, providing a substantial improvement over the $\sim 10\%$ results previously obtained using density-only conditioning. Overall, this makes our method a simple yet efficient tool for forward modeling HI maps during reionization.

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 manuscript presents a method to compute the collapse fraction f_coll by conditioning on both the dark matter overdensity δ and the eigenvalues of the tidal tensor. Conditional means are measured from high-resolution small-volume simulations and applied to populate f_coll in low-resolution large-volume simulations. These f_coll fields are then used in a semi-numerical reionization model to generate the HI density field and power spectrum, which is benchmarked against a reference high-resolution large-volume simulation. The central claim is that this tidal conditioning recovers the large-scale HI power spectrum to ≲2–5% accuracy for k ≲ 0.5 h Mpc^{-1} across redshifts, ionized fractions, resolutions, and halo masses, improving on prior density-only results at the ~10% level.

Significance. If the transfer of conditional means is robust, the approach offers an efficient route to accurate large-volume HI maps without full high-resolution runs, which is valuable for 21 cm cosmology forecasts. The direct benchmarking against reference simulations and the quantified error levels across multiple parameters provide concrete support for the improvement over density-only conditioning. The physical motivation from tidal environment for halo formation is a clear strength.

major comments (2)
  1. [§2] §2 (conditioning procedure): The central method transfers the conditional mean f_coll(δ, λ) measured in small high-resolution volumes directly to large low-resolution volumes. No explicit test or correction is presented for differences in long-wavelength tidal mode sampling or grid smoothing between the volumes; if the conditional mean depends on these, the f_coll field could be biased even when the downstream HI power spectrum appears accurate. The validation is performed only on the final power spectrum after semi-numerical reionization, not on the fidelity of the conditional means themselves.
  2. [§4] §4 (results): While errors are reported at the 2–5% level on large scales, the manuscript does not show a direct comparison of the measured conditional means f_coll(δ, λ) between the small-box conditioning set and what would be obtained in the target large volume, leaving open whether the accuracy is coincidental or robust to the volume/resolution mismatch.
minor comments (2)
  1. [Abstract] The abstract and introduction should include a specific citation to the prior density-only work that achieved ~10% accuracy for context.
  2. [Figures] Figure legends and axis labels in the power-spectrum comparison panels should be made fully consistent across redshifts and resolutions to aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript, positive assessment of the work, and recommendation for minor revision. The comments on the conditioning procedure and validation are well taken, and we have revised the manuscript to address them directly.

read point-by-point responses

  1. Referee: §2 (conditioning procedure): The central method transfers the conditional mean f_coll(δ, λ) measured in small high-resolution volumes directly to large low-resolution volumes. No explicit test or correction is presented for differences in long-wavelength tidal mode sampling or grid smoothing between the volumes; if the conditional mean depends on these, the f_coll field could be biased even when the downstream HI power spectrum appears accurate. The validation is performed only on the final power spectrum after semi-numerical reionization, not on the fidelity of the conditional means themselves.

    Authors: We agree that a more explicit validation of the conditional means themselves would strengthen the presentation. The method assumes that the local dependence of f_coll on δ and the tidal eigenvalues λ is sufficiently universal that it can be transferred between volumes, with the large-volume simulation supplying the correct long-wavelength modes. To address the referee's concern, we have added to the revised §2 both a comparison of the tidal eigenvalue distributions sampled in the small high-resolution conditioning boxes versus the target large volume and a direct plot of the conditional means f_coll(δ, λ) measured in the conditioning set and as applied to the large-volume grid. These additions confirm that the relevant ranges are well sampled and that the means agree within the expected cosmic variance of the small boxes. While the primary scientific validation remains the recovered HI power spectrum, we now also show the intermediate fidelity of the f_coll field. revision: yes

  2. Referee: §4 (results): While errors are reported at the 2–5% level on large scales, the manuscript does not show a direct comparison of the measured conditional means f_coll(δ, λ) between the small-box conditioning set and what would be obtained in the target large volume, leaving open whether the accuracy is coincidental or robust to the volume/resolution mismatch.

    Authors: We thank the referee for this suggestion. In the revised §4 we now include a direct side-by-side comparison of the conditional means f_coll(δ, λ) extracted from the small-box conditioning simulations and the values realized on the large-volume grid. The two agree to within ~5% across the populated range of δ and λ, with residual differences consistent with the finite volume of the conditioning runs. This comparison demonstrates that the reported 2–5% accuracy on the large-scale HI power spectrum is robust to the volume and resolution differences rather than coincidental. revision: yes

Circularity Check

0 steps flagged

No significant circularity; conditional statistics measured independently and validated externally

full rationale

The derivation measures conditional means of f_coll on density and tidal eigenvalues from independent high-resolution small-volume simulations, then applies those means to sample f_coll in a separate low-resolution large-volume run. The HI power spectrum is subsequently computed via semi-numerical reionization and benchmarked directly against a reference high-resolution large-volume simulation. This external validation step ensures no claimed accuracy reduces by construction to the measured conditional means themselves. No self-citations, fitted parameters renamed as predictions, self-definitional equations, or ansatzes imported via prior work appear in the methodology.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method rests on standard cosmological assumptions about halo collapse and tidal fields plus the transferability of conditional statistics across simulation volumes.

axioms (1)
  • domain assumption Halo collapse fraction depends on local density and tidal environment in a manner that can be captured by conditional means from smaller simulations.
    Invoked when the authors compute mean f_coll conditioned on eigenvalues and apply it to large-volume runs.

pith-pipeline@v0.9.0 · 5821 in / 1266 out tokens · 32787 ms · 2026-05-18T16:46:05.768508+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

40 extracted references · 40 canonical work pages · 1 internal anchor

  1. [1]

    Gnedin and P

    N.Y. Gnedin and P. Madau,Modeling cosmic reionization,Living Reviews in Computational Astrophysics8(2022) 3

    work page 2022

  2. [2]

    Choudhury,A short introduction to reionization physics,General Relativity and Gravitation54(2022)

    T.R. Choudhury,A short introduction to reionization physics,General Relativity and Gravitation54(2022)

    work page 2022

  3. [3]

    Furlanetto, S

    S.R. Furlanetto, S. Peng Oh and F.H. Briggs,Cosmology at low frequencies: The 21cm transition and the high-redshift universe,Physics Reports433(2006) 181

    work page 2006

  4. [4]

    Pritchard and A

    J.R. Pritchard and A. Loeb,21 cm cosmology in the 21st century,Reports on Progress in Physics75(2012) 086901

    work page 2012

  5. [5]

    2019, The Cosmic 21-Cm Revolution, 2514-3433 (IOP Publishing), doi: 10.1088/2514-3433/ab4a73

    A. Mesinger, ed.,The Cosmic 21-cm Revolution, 2514-3433, IOP Publishing (2019), 10.1088/2514-3433/ab4a73

    work page doi:10.1088/2514-3433/ab4a73 2019

  6. [6]

    Gnedin,Cosmological reionization by stellar sources,The Astrophysical Journal535 (2000) 530

    N.Y. Gnedin,Cosmological reionization by stellar sources,The Astrophysical Journal535 (2000) 530

    work page 2000

  7. [7]

    Ciardi, A

    B. Ciardi, A. Ferrara and S.D.M. White,Early reionization by the first galaxies,Monthly Notices of the Royal Astronomical Society344(2003) L7 [https://academic.oup.com/mnras/article-pdf/344/1/L7/18652206/344-1-L7.pdf]

    work page 2003

  8. [8]

    Mellema, I.T

    G. Mellema, I.T. Iliev, U.-L. Pen and P.R. Shapiro,Simulating cosmic reionization at large scales – ii. the 21-cm emission features and statistical signals,Monthly Notices of the Royal Astronomical Society372(2006)679 [https://academic.oup.com/mnras/article-pdf/372/2/679/2986158/mnras0372-0679.pdf]

    work page 2006

  9. [9]

    Trac and R

    H. Trac and R. Cen,Radiative transfer simulations of cosmic reionization. i. methodology and initial results,The Astrophysical Journal671(2007) 1

    work page 2007

  10. [10]

    Gnedin,Cosmic reionization on computers

    N.Y. Gnedin,Cosmic reionization on computers. i. design and calibration of simulations,The Astrophysical Journal793(2014) 29

    work page 2014

  11. [11]

    Rosdahl, H

    J. Rosdahl, H. Katz, J. Blaizot, T. Kimm, L. Michel-Dansac, T. Garel et al.,The sphinx cosmological simulations of the first billion years: the impact of binary stars on reionization, Monthly Notices of the Royal Astronomical Society479(2018) 994 [https://academic.oup.com/mnras/article-pdf/479/1/994/25129300/sty1655.pdf]

    work page 2018

  12. [12]

    Lewis, P

    J.S.W. Lewis, P. Ocvirk, J.G. Sorce, Y. Dubois, D. Aubert, L. Conaboy et al.,The short ionizing photon mean free path at z = 6 in cosmic dawn iii, a new fully coupled radiation-hydrodynamical simulation of the epoch of reionization,Monthly Notices of the Royal Astronomical Society516(2022) 3389 [https://academic.oup.com/mnras/article-pdf/516/3/3389/458827...

    work page 2022

  13. [13]

    Iliev, G

    I.T. Iliev, G. Mellema, K. Ahn, P.R. Shapiro, Y. Mao and U.-L. Pen,Simulating cosmic reionization: how large a volume is large enough?,Monthly Notices of the Royal Astronomical Society439(2014) 725 [https://academic.oup.com/mnras/article-pdf/439/1/725/5599101/stt2497.pdf]

    work page 2014

  14. [14]

    H.D. Kaur, N. Gillet and A. Mesinger,Minimum size of 21-cm simulations,Monthly Notices of the Royal Astronomical Society495(2020) 2354 [https://academic.oup.com/mnras/article-pdf/495/2/2354/33323159/staa1323.pdf]

    work page 2020

  15. [15]

    J.R. Bond, S. Cole, G. Efstathiou and N. Kaiser,Excursion Set Mass Functions for Hierarchical Gaussian Fluctuations, ApJ379(1991) 440

    work page 1991

  16. [16]

    Furlanetto, M

    S.R. Furlanetto, M. Zaldarriaga and L. Hernquist,The growth of HII regions during reionization,The Astrophysical Journal613(2004) 1

    work page 2004

  17. [17]

    Mesinger and S

    A. Mesinger and S. Furlanetto,Efficient simulations of early structure formation and reionization,The Astrophysical Journal669(2007) 663

    work page 2007

  18. [18]

    O. Zahn, A. Lidz, M. McQuinn, S. Dutta, L. Hernquist, M. Zaldarriaga et al.,Simulations and analytic calculations of bubble growth during hydrogen reionization,The Astrophysical Journal 654(2007) 12

    work page 2007

  19. [19]

    Choudhury, M.G

    T.R. Choudhury, M.G. Haehnelt and J. Regan,Inside-out or outside-in: the topology of reionization in the photon-starved regime suggested by lyαforest data,Monthly Notices of the Royal Astronomical Society394(2009) 960 [https://academic.oup.com/mnras/article-pdf/394/2/960/3710519/mnras0394-0960.pdf]

    work page 2009

  20. [20]

    Mesinger, S

    A. Mesinger, S. Furlanetto and R. Cen,21cmfast: a fast, seminumerical simulation of the high-redshift 21-cm signal,Monthly Notices of the Royal Astronomical Society411(2011) 955 [https://academic.oup.com/mnras/article-pdf/411/2/955/4099991/mnras0411-0955.pdf]

    work page 2011

  21. [21]

    Lin, S.P

    Y. Lin, S.P. Oh, S.R. Furlanetto and P.M. Sutter,The distribution of bubble sizes during reionization,Monthly Notices of the Royal Astronomical Society461(2016) 3361 [https://academic.oup.com/mnras/article-pdf/461/3/3361/8112292/stw1542.pdf]

    work page 2016

  22. [22]

    Choudhury and A

    T.R. Choudhury and A. Paranjape,Photon number conservation and the large-scale 21 cm power spectrum in seminumerical models of reionization,Monthly Notices of the Royal Astronomical Society481(2018) 3821 [https://academic.oup.com/mnras/article-pdf/481/3/3821/25844366/sty2551.pdf]

    work page 2018

  23. [23]

    Press and P

    W.H. Press and P. Schechter,Formation of Galaxies and Clusters of Galaxies by Self-Similar Gravitational Condensation, ApJ187(1974) 425

    work page 1974

  24. [24]

    Large scale bias and the peak background split

    R.K. Sheth and G. Tormen,Large-scale bias and the peak background split, MNRAS308(1999) 119 [astro-ph/9901122]

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  25. [25]

    Sheth and G

    R.K. Sheth and G. Tormen,An excursion set model of hierarchical clustering: ellipsoidal collapse and the moving barrier,Monthly Notices of the Royal Astronomical Society329(2002) 61 [https://academic.oup.com/mnras/article-pdf/329/1/61/3882215/329-1-61.pdf]

    work page 2002

  26. [26]

    D.S. Reed, R. Bower, C.S. Frenk, A. Jenkins and T. Theuns,The halo mass function from the dark ages through the present day,Monthly Notices of the Royal Astronomical Society374 (2006) 2 [https://academic.oup.com/mnras/article-pdf/374/1/2/2835466/mnras0374-0002.pdf]

    work page 2006

  27. [27]

    Tinker, A.V

    J. Tinker, A.V. Kravtsov, A. Klypin, K. Abazajian, M. Warren, G. Yepes et al.,Toward a halo mass function for precision cosmology: The limits of universality,The Astrophysical Journal 688(2008) 709

    work page 2008

  28. [28]

    Courtin, Y

    J. Courtin, Y. Rasera, J.-M. Alimi, P.-S. Corasaniti, V. Boucher and A. Füzfa,Imprints of dark energy on cosmic structure formation – ii. non-universality of the halo mass function, – 12 – Monthly Notices of the Royal Astronomical Society410(2011) 1911 [https://academic.oup.com/mnras/article-pdf/410/3/1911/2864964/mnras0410-1911.pdf]

    work page 2011

  29. [29]

    Crocce, P

    M. Crocce, P. Fosalba, F.J. Castander and E. Gaztañaga,Simulating the universe with mice: the abundance of massive clusters,Monthly Notices of the Royal Astronomical Society403 (2010) 1353 [https://academic.oup.com/mnras/article-pdf/403/3/1353/6170753/mnras0403-1353.pdf]

    work page 2010

  30. [30]

    Ahn, I.T

    K. Ahn, I.T. Iliev, P.R. Shapiro, G. Mellema, J. Koda and Y. Mao,Detecting the rise and fall of the first stars by their impact on cosmic reionization,The Astrophysical Journal Letters756 (2012) L16

    work page 2012

  31. [31]

    Pundir, A

    G. Pundir, A. Paranjape and T.R. Choudhury,Accelerating hi density predictions during the epoch of reionization using a gpr-based emulator on n-body simulations,Journal of Cosmology and Astroparticle Physics2025(2025) 045

    work page 2025

  32. [32]

    O. Hahn, C. Porciani, C.M. Carollo and A. Dekel,Properties of dark matter haloes in clusters, filaments, sheets and voids,Monthly Notices of the Royal Astronomical Society375(2007) 489 [https://academic.oup.com/mnras/article-pdf/375/2/489/4244383/mnras0375-0489.pdf]

    work page 2007

  33. [33]

    Forero–Romero, Y

    J.E. Forero–Romero, Y. Hoffman, S. Gottlöber, A. Klypin and G. Yepes,A dynamical classification of the cosmic web,Monthly Notices of the Royal Astronomical Society396(2009) 1815 [https://academic.oup.com/mnras/article-pdf/396/3/1815/5804803/mnras0396-1815.pdf]

    work page 2009

  34. [34]

    real-space power spectra,A&A661(2022) A146

    Bonnaire, Tony, Aghanim, Nabila, Kuruvilla, Joseph and Decelle, Aurélien,Cosmology with cosmic web environments - i. real-space power spectra,A&A661(2022) A146

    work page 2022

  35. [35]

    Barsode and T.R

    A. Barsode and T.R. Choudhury,Efficient hybrid technique for generating sub-grid haloes in reionization simulations,Journal of Cosmology and Astroparticle Physics2024(2024) 036

    work page 2024

  36. [36]

    V. Springel,The cosmological simulation code gadget-2,Monthly Notices of the Royal Astronomical Society364(2005) 1105 [https://academic.oup.com/mnras/article-pdf/364/4/1105/18657201/364-4-1105.pdf]

    work page 2005

  37. [37]

    et al.,Planck 2018 results - vi

    Planck Collaboration, Aghanim, N., Akrami, Y., Ashdown, M., Aumont, J., Baccigalupi, C. et al.,Planck 2018 results - vi. cosmological parameters,A&A641(2020) A6

    work page 2018

  38. [38]

    Davis, G

    M. Davis, G. Efstathiou, C.S. Frenk and S.D.M. White,The evolution of large-scale structure in a universe dominated by cold dark matter, ApJ292(1985) 371

    work page 1985

  39. [39]

    McQuinn, A

    M. McQuinn, A. Lidz, O. Zahn, S. Dutta, L. Hernquist and M. Zaldarriaga,The morphology of h ii regions during reionization,Monthly Notices of the Royal Astronomical Society377(2007) 1043 [https://academic.oup.com/mnras/article-pdf/377/3/1043/5678910/mnras0377-1043.pdf]

    work page 2007

  40. [40]

    Doussot, Aristide and Semelin, Benoît,A bubble size distribution model for the epoch of reionization, A&A667(2022) A118. – 13 –

    work page 2022