pith. machine review for the scientific record. sign in

arxiv: 2605.10354 · v1 · submitted 2026-05-11 · 🌌 astro-ph.CO

Recognition: 2 theorem links

· Lean Theorem

STEPSIC: Initial condition generator for stereographic cosmological simulations

Bal\'azs P\'al , G\'abor R\'acz , Istv\'an Csabai , Istv\'an Szapudi
Authors on Pith no claims yet

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

classification 🌌 astro-ph.CO
keywords initial conditionsLagrangian perturbation theorycosmological simulationsstereographic projectionN-bodypower spectrumStePS
0
0 comments X

The pith

STEPSIC generates initial conditions for stereographic and cylindrical cosmological volumes using extended Lagrangian perturbation theory.

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

The paper introduces STEPSIC to produce initial conditions for N-body simulations in spherical and cylindrical geometries that stereographic codes like StePS require, rather than restricting to periodic cubic boxes. It builds Gaussian random fields on cubic Fourier grids, applies first- and second-order LPT, interpolates displacements with B-spline kernels, and maps fields across radially varying resolutions via a multiresolution scheme. This matters because conventional generators cannot support observer-centric or irregular domains without introducing boundary artifacts. Tests confirm power spectrum recovery to better than 0.5% of linear theory up to half the Nyquist scale, independent of aspect ratio, and full evolution runs show second-order LPT removes transient excess power.

Core claim

STEPSIC constructs Gaussian random density fields on anisotropy-free Fourier grids with cubic voxels, applies first- and second-order LPT to obtain displacement and velocity fields, interpolates these onto particles via B-spline mass-assignment kernels with Fourier-space deconvolution, and for stereographic geometries employs a multiresolution scheme to map displacement fields across the radially varying particle mass resolution. In periodic cubic boxes the recovered matter power spectrum agrees with the input linear theory prediction to better than 0.5% up to half the Nyquist wavenumber independent of box aspect ratio, and full N-body evolution of matched cylindrical StePS runs confirms the

What carries the argument

Multiresolution mapping scheme that transfers displacement fields across radially varying particle mass resolution in stereographic projections, combined with B-spline interpolation after LPT.

If this is right

  • Power spectrum recovery matches linear theory to better than 0.5% up to half Nyquist wavenumber and holds independent of box aspect ratio up to 10:1.
  • Second-order LPT suppresses the 2-3% transient power excess that appears when first-order initial conditions are evolved.
  • Paired-and-fixed variance-reduced realizations are generated with the same accuracy as standard ones.
  • Identical white-noise fields produce sub-percent power-spectrum agreement with independent generators such as monofonic.

Where Pith is reading between the lines

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

  • This method opens the door to observer-centric simulations that avoid periodic boundary artifacts when studying large-scale structure around a chosen location.
  • The same multiresolution approach could be tested in full spherical StePS volumes to verify performance beyond the periodic validation cases.
  • Higher-order perturbation schemes or different mass-assignment kernels could be added to further reduce small-scale transients in non-cubic domains.

Load-bearing premise

The multiresolution scheme for mapping displacement fields across radially varying particle mass resolution in stereographic geometries introduces no significant artifacts or biases.

What would settle it

Compare the early-time power spectrum measured in an N-body run started from STEPSIC cylindrical initial conditions against the linear theory input and check whether any deviation exceeds 0.5% below half the Nyquist scale.

Figures

Figures reproduced from arXiv: 2605.10354 by Bal\'azs P\'al, G\'abor R\'acz, Istv\'an Csabai, Istv\'an Szapudi.

Figure 1
Figure 1. Figure 1: Two-dimensional slices through the four particle load configurations supported by stepsic, shown before any Lagrangian displacement is applied or glass relaxation is performed. Each panel displays particles within a thin slab of thickness 25 h −1Mpc centred on the midplane and oriented perpendicular to the x-axis, projected onto the y–z plane. (a) Uniform random (Poisson) sampling within a cubical T 3 doma… view at source ↗
Figure 2
Figure 2. Figure 2: Particle mass distribution as a function of radial distance for the constant stereographic angle ω (solid) and constant compact-space volume (dashed) binning modes in the spherical R 3 geometry (top) and the cylindrical S 1×R 2 geometry along the non-periodic R 2 directions (bottom). The constant ∆ω profile includes a constant-resolution inner zone with rc = 50 h −1Mpc. The compactification diameter is D4D… view at source ↗
Figure 3
Figure 3. Figure 3: Distribution of per-component absolute displacements in a slab geometry with dimensions Lx = Ly = 1000 h −1Mpc and Lz = 200 h −1Mpc, averaged over 200 realisations using second-order LPT at z = 31. The transverse curve |Ψ⊥| is the mean of the x and y compo￾nents, exploiting the exact x ↔ y symmetry of the slab. The longitudinal component |Ψ∥ | ≡ |Ψz | is strongly suppressed, reflecting the anisotropic infr… view at source ↗
Figure 5
Figure 5. Figure 5: Starting from a cubic box of side L = 1000 h −1Mpc, the z￾dimension is progressively reduced to Lz = 100 h −1Mpc while the trans￾verse dimensions remain fixed. To ensure a fair comparison, the physi￾cal cell size ∆ = Lz,min/Ninput is held constant across all boxes; larger boxes receive proportionally more mesh cells so that spatial resolu￾tion is identical in every case. All realizations use second-order L… view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of different sources of systematic effects: (a) Res￾olution dependence: comparison of multiple mesh resolutions (N 3 ) at a fixed redshift of z = 31, using second-order LPT. (b) Redshift de￾pendence: comparison across different target redshifts at the fixed res￾olution of 5123 , using second-order LPT. (c) LPT order: comparison between first- and second-order LPT at the fixed resolution of 5123 … view at source ↗
Figure 6
Figure 6. Figure 6: Qualitative visualization of the projected density field at z = 0 from the second-order LPT cylindrical S 1×R 2 StePS simulation used for the validation presented in [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Ratio of the matter power spectra at z = 0 from two matched cylindrical S 1×R 2 StePS simulations evolved from first-order and second-order LPT initial conditions, respectively. Both runs share the same white noise realization (∼ 2 × 106 particles, cylindrical glass, R3D = 750 Mpc, D4D = 35 Mpc, Lz = 200 Mpc, zinit = 31). Power spectra were measured with the non-periodic Feldman–Kaiser–Peacock estimator of… view at source ↗
Figure 8
Figure 8. Figure 8: Cross-validation of stepsic against monofonic, both initialised from the same white noise realization in a fully periodic cubic box of side L = 1000 h −1Mpc with 5123 particles at z = 31, using second￾order LPT. Shaded bands in all panels indicate ± 1% deviation from unity. (a) Ratio of the matter power spectra P(k), measured with inter￾laced CIC deposition and deconvolution. The median ratio is ≈ 0.994; t… view at source ↗
read the original abstract

Conventional cosmological initial condition generators are designed exclusively for fully periodic cubic domains and cannot produce the non-periodic, observer-centric configurations required by stereographically projected N-body codes such as StePS. We present STEPSIC, an open-source initial condition generator that extends Lagrangian perturbation theory-based initial conditions to the spherical and cylindrical geometries used by StePS, while also supporting cuboid domains with arbitrary aspect ratios. The code constructs Gaussian random density fields on anisotropy-free Fourier grids with cubic voxels, applies first- and second-order LPT to obtain displacement and velocity fields, and interpolates these onto particles via B-spline mass-assignment kernels with Fourier-space deconvolution. For stereographic geometries, a multiresolution scheme maps displacement fields across the radially varying particle mass resolution intrinsic to the projection. Both standard and paired-and-fixed variance-reduced realizations are supported. In periodic cubic boxes, the recovered matter power spectrum agrees with the input linear theory prediction to better than 0.5% up to half the Nyquist wavenumber, independent of box aspect ratio (tested up to 10:1). Cross-validation against monofonic using identical white noise fields yields sub-percent power spectrum agreement, with a small residual offset consistent with differences between two independent implementations. Full N-body evolution of matched cylindrical StePS runs confirms that second-order LPT correctly suppresses the 2-3% transient power excess present in first-order initial conditions.

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

1 major / 0 minor

Summary. The manuscript introduces STEPSIC, an open-source initial condition generator extending Lagrangian perturbation theory (first- and second-order) to non-periodic stereographic and cylindrical geometries required by StePS N-body codes, while also supporting cuboid domains of arbitrary aspect ratio. It constructs Gaussian random fields on cubic-voxel Fourier grids, applies LPT, interpolates displacements/velocities via B-spline kernels with Fourier deconvolution, and uses a multiresolution scheme to handle radially varying particle masses in stereographic projections. Reported results include matter power spectrum recovery to better than 0.5% versus linear theory up to k_Nyq/2 in periodic cubic boxes (independent of aspect ratio up to 10:1), sub-percent agreement with monofonic on identical white-noise fields, and N-body confirmation that 2LPT suppresses the 2-3% transient excess seen in 1LPT for cylindrical StePS runs. Both standard and paired-and-fixed realizations are supported.

Significance. If the multiresolution mapping performs as intended, STEPSIC would provide a needed capability for accurate initial conditions in observer-centric stereographic simulations, which are otherwise unsupported by standard cubic periodic generators. The cubic-box power-spectrum recovery and monofonic cross-check supply quantitative evidence of correctness for periodic cases, while the cylindrical N-body test directly demonstrates the benefit of 2LPT for transient control. These elements strengthen the paper's contribution for the supported geometries.

major comments (1)
  1. [Abstract and validation results] The multiresolution displacement-mapping scheme for stereographic geometries (described in the abstract as mapping LPT fields across radially varying particle mass resolution) is central to the code's purpose yet receives no direct validation. All quantitative tests (power-spectrum recovery to <0.5%, monofonic comparison, and cylindrical N-body evolution) are performed on periodic cubic or uniform-resolution cylindrical domains and therefore do not exercise the radial shell transitions or B-spline interpolation across resolution boundaries. Any mismatch in deconvolution or mass assignment at these interfaces could introduce localized density biases or spurious power that would not appear in the reported tests.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comment on validation. We address the major comment below and agree that strengthening the validation for the multiresolution scheme is warranted. The revised manuscript will incorporate additional tests.

read point-by-point responses

  1. Referee: [Abstract and validation results] The multiresolution displacement-mapping scheme for stereographic geometries (described in the abstract as mapping LPT fields across radially varying particle mass resolution) is central to the code's purpose yet receives no direct validation. All quantitative tests (power-spectrum recovery to <0.5%, monofonic comparison, and cylindrical N-body evolution) are performed on periodic cubic or uniform-resolution cylindrical domains and therefore do not exercise the radial shell transitions or B-spline interpolation across resolution boundaries. Any mismatch in deconvolution or mass assignment at these interfaces could introduce localized density biases or spurious power that would not appear in the reported tests.

    Authors: We agree with the referee that the quantitative tests reported in the manuscript are performed on periodic cubic boxes and uniform-resolution cylindrical domains, which do not directly exercise the radial resolution transitions in the stereographic multiresolution scheme. The core B-spline interpolation with Fourier deconvolution is validated through the power-spectrum recovery and monofonic cross-check in uniform cases, and the multiresolution mapping applies the same kernels across shells. However, we acknowledge that explicit checks for interface consistency (e.g., density continuity or absence of spurious power at shell boundaries) are not presented. We will revise the manuscript to include a dedicated validation subsection for stereographic geometries. This will consist of (i) generating matched initial conditions in a stereographic domain and a high-resolution equivalent uniform grid, (ii) measuring the matter power spectrum in thin radial shells to verify sub-percent agreement with linear theory across transitions, and (iii) confirming that mass assignment and deconvolution remain consistent at resolution boundaries. These additions will be placed in Section 4 or a new appendix, with corresponding updates to the abstract and conclusions if needed. revision: yes

Circularity Check

0 steps flagged

No circularity: validations use external linear theory and independent code cross-checks

full rationale

The paper extends standard LPT IC generation to stereographic/cylindrical geometries via explicit B-spline interpolation and a multiresolution displacement mapping. All reported results compare the generated fields to external linear theory power spectra or to monofonic outputs on shared white noise; the N-body test confirms the known 2LPT transient suppression. No equation or claim reduces by construction to a fitted parameter, self-citation, or renamed input. The multiresolution scheme is described procedurally without any self-referential validation loop.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work relies on standard cosmological assumptions and prior LPT methods without new free parameters or invented entities; all core techniques are drawn from established literature.

axioms (2)
  • domain assumption Initial density fluctuations form a Gaussian random field drawn from linear theory power spectrum
    Invoked in the construction of the Fourier-space density field on the grid.
  • domain assumption Lagrangian perturbation theory at first and second order provides accurate displacements and velocities at early times
    Used to obtain displacement and velocity fields before interpolation.

pith-pipeline@v0.9.0 · 5570 in / 1388 out tokens · 40650 ms · 2026-05-12T05:16:02.848350+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

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

  1. [1]

    Angulo, R. E. & Pontzen, A. 2016, MNRAS, 462, L1

    work page 2016

  2. [2]

    Bagla, J. S. & Padmanabhan, T. 1997, Pramana, 49, 161 Benítez-Llambay, A. 2025, Research Notes of the AAS, 9, 122

    work page 1997

  3. [3]

    2002, Phys

    Bernardeau, F., Colombi, S., Gaztañaga, E., & Scoccimarro, R. 2002, Phys. Rep., 367, 1

    work page 2002

  4. [4]

    R., Colombi, S., Hivon, E., & Juszkiewicz, R

    Bouchet, F. R., Colombi, S., Hivon, E., & Juszkiewicz, R. 1995, A&A, 296, 575

    work page 1995

  5. [5]

    Bouchet, F. R. & Hernquist, L. 1992, ApJ, 400, 25

    work page 1992

  6. [6]

    2006, MNRAS, 373, 369

    Crocce, M., Pueblas, S., & Scoccimarro, R. 2006, MNRAS, 373, 369

    work page 2006

  7. [7]

    S., & White, S

    Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, Astrophysical Journal, 292, 371

    work page 1985

  8. [8]

    2018, ApJS, 239, 35

    Diemer, B. 2018, ApJS, 239, 35

    work page 2018

  9. [9]

    Eisenstein, D. J. & Hu, W. 1998, ApJ, 496, 605 Euclid Collaboration, Castander, F. J., Fosalba, P., et al. 2025, A&A, 697, A5

    work page 1998

  10. [10]

    & Abel, T

    Hahn, O. & Abel, T. 2011, MNRAS, 415, 2101

    work page 2011

  11. [11]

    Hahn, O., List, F., & Porqueres, N. 2024, J. Cosmology Astropart. Phys., 2024, 063

    work page 2024

  12. [12]

    2021, MNRAS, 503, 426

    Hahn, O., Rampf, C., & Uhlemann, C. 2021, MNRAS, 503, 426

    work page 2021

  13. [13]

    Hockney, R. W. & Eastwood, J. W. 1988, Computer simulation using particles (Bristol and Philadelphia: Institute of Physics Publishing)

    work page 1988

  14. [14]

    Jing, Y . P. 2005, ApJ, 620, 559

    work page 2005

  15. [15]

    & Holtzman, J

    Klypin, A. & Holtzman, J. 1997, arXiv e-prints, astro

    work page 1997

  16. [16]

    2020, MNRAS, 496, 3862

    Klypin, A., Prada, F., & Byun, J. 2020, MNRAS, 496, 3862

    work page 2020

  17. [17]

    2000, ApJ, 538, 473

    Lewis, A., Challinor, A., & Lasenby, A. 2000, ApJ, 538, 473

    work page 2000

  18. [18]

    2017, arXiv e-prints, arXiv:1709.04703

    Llinares, C. 2017, arXiv e-prints, arXiv:1709.04703

    work page arXiv 2017

  19. [19]

    Michaux, M., Hahn, O., Rampf, C., & Angulo, R. E. 2021, MNRAS, 500, 663

    work page 2021

  20. [20]

    R., Pellat, R., & Ramani, A

    Moutarde, F., Alimi, J.-M., Bouchet, F. R., Pellat, R., & Ramani, A. 1991, ApJ, 382, 377

    work page 1991

  21. [21]

    Peebles, P. J. E. 1980, The large-scale structure of the universe (Princeton, NJ: Princeton University Press)

    work page 1980

  22. [22]

    Pontzen, A., Slosar, A., Roth, N., & Peiris, H. V . 2016, Phys. Rev. D, 93, 103519

    work page 2016

  23. [23]

    H., Teukolsky, S

    Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 2007, Nu- merical Recipes: The Art of Scientific Computing, 3rd edn. (Cambridge, UK: Cambridge University Press) Rácz, G., Szapudi, I., Csabai, I., & Dobos, L. 2018, Monthly Notices of the Royal Astronomical Society, 477, 1949 Rácz, G., Szapudi, I., Csabai, I., & Dobos, L. 2021, MNRAS, 5...

    work page internal anchor Pith review arXiv 2007

  24. [24]

    Springel, V ., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629

    work page 2005

  25. [25]

    R., & White, S

    Tormen, G., Bouchet, F. R., & White, S. D. M. 1997, MNRAS, 286, 865

    work page 1997

  26. [26]

    White, S. D. M. 1994, arXiv:astro-ph/9410043 Zel’dovich, Y . B. 1970, A&A, 5, 84 Article number, page 10 of 10

    work page arXiv 1994