Pith. sign in

REVIEW 4 major objections 4 minor 74 references

The Sahyadri simulation suite resolves dark matter halos down to 3.2e9 solar masses, a factor of about 25 better than the AbacusSummit suite, and opens the faint-galaxy regime of the DESI Bright Galaxy Survey to cosmological parameter studi

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-03 10:59 UTC pith:43ZYX42Z

load-bearing objection The core resolution claim checks out, but the abstract promises six parameters and derivative-ready data that the body does not yet deliver. the 4 major comments →

arxiv 2601.07924 v2 pith:43ZYX42Z submitted 2026-01-12 astro-ph.CO

Sahyadri: A simulation suite for the cosmology dependence of the Cosmic Web

classification astro-ph.CO
keywords cosmological N-body simulationscosmic webDESI Bright Galaxy Surveyhalo mass functionVoronoi volume functionk-nearest neighbour statisticsseed-matched initial conditionsOmega_m sensitivity
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper presents Sahyadri, a suite of dark-matter N-body simulations that evolve 2048^3 particles in a 200 h^-1 Mpc box, reaching a particle mass of 8.1e7 h^-1 M_sun. This resolves halos down to 3.2e9 h^-1 M_sun, about 25 times lower mass than the previous best parameter-varying suite, AbacusSummit. The suite varies cosmological parameters with matched initial-condition seeds, so differences between runs isolate the effect of each parameter, enabling derivative and Fisher-matrix calculations. The authors estimate this resolution brings roughly 40% of DESI BGS galaxies at z<0.15 into reach, and they demonstrate clear Omega_m sensitivity in both standard clustering and cosmic-web statistics such as the Voronoi volume function and k-nearest-neighbour distributions.

Core claim

The central claim is that Sahyadri supplies the mass resolution and parameter coverage needed to model the faint galaxy population targeted by low-redshift spectroscopic surveys, especially DESI BGS. By resolving halos at 3.2e9 h^-1 M_sun, the suite captures roughly 1.6 million BGS galaxies at z<0.15 that AbacusSummit cannot reach. The paper shows that the simulations reproduce theoretical matter power spectra and halo mass functions, and that the Voronoi volume function and k-nearest-neighbour statistics respond strongly to Omega_m variations, making these high-density cosmic-web statistics usable for cosmological inference.

What carries the argument

The carrying element is the combination of a compact box (200 h^-1 Mpc) with a very large particle count (2048^3), yielding a particle mass of 8.1e7 h^-1 M_sun and resolved halos down to 3.2e9 h^-1 M_sun (40 particles). Seed-matched initial conditions across parameter variations isolate cosmological responses, while a custom compression scheme reduces the data footprint by about a factor of three, keeping the ~117 TB total storage manageable.

Load-bearing premise

The estimate that about 40% of DESI BGS galaxies at z<0.15 become accessible relies on the assumed relation between galaxy stellar mass and halo mass, and between stellar mass and r-band luminosity, calibrated in an unpublished companion methodology; if that galaxy-halo connection is wrong, the headline accessibility gain over AbacusSummit is unsupported even though the mass-resolution claim stands.

What would settle it

Compare the same halo finder on the Sahyadri snapshots with a version degraded to AbacusSummit particle mass (2e9 h^-1 M_sun) and verify that the 40-particle minimum halo mass shifts by the expected factor of ~25. Separately, generate mock BGS galaxy catalogs from Sahyadri using the paper's assumed stellar-mass-to-halo-mass relation and check whether the predicted BGS stellar mass function at z=0.15 matches observations; a significant mismatch would invalidate the 40% figure.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • DESI BGS galaxies at z<0.15 can be modeled at the halo level for roughly 40% more galaxies than with AbacusSummit, enabling mocks and systematics studies for the faint end of the survey.
  • Seed-matched runs allow cosmological parameter derivatives for beyond-two-point statistics, opening the door to Fisher forecasts that combine VVF, kNN, and the power spectrum.
  • The demonstrated Omega_m sensitivity of high-density VVF and kNN statistics suggests these observables can contribute constraints from the non-linear regime where two-point statistics are cosmic-variance limited.
  • The same simulation infrastructure can be extended to 4MOST and other low-redshift surveys, and to higher-order statistics like Minkowski functionals.

Where Pith is reading between the lines

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

  • The 40% BGS accessibility figure depends on assumed stellar-to-halo and mass-to-light relations applied to BGS targets; if those calibrations shift, the percentage changes even though the raw mass-resolution improvement stands.
  • The small 200 h^-1 Mpc box limits access to large-scale modes; the paper itself shows finite-volume truncation of halo-by-halo bias tails, so users should correct for this when interpreting environment statistics.
  • Only Omega_m, h, and n_s variations are complete so far; the planned full six-parameter derivative set is needed before the suite can be used as a general emulator across the full parameter space.
  • The compression scheme's sub-percent clustering accuracy would be worth independent verification with a different power-spectrum estimator, since the claim rests on a single implementation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper presents Sahyadri, a suite of dark-matter-only N-body simulations with 2048^3 particles in 200 h^-1 Mpc boxes, reaching a particle mass of 8.1e7 h^-1 M_sun and halo masses down to 3.2e9 h^-1 M_sun with 40 particles. The suite is designed around Planck 2018 cosmology with parameter variations, seed-matched initial conditions, and a claimed capability for cosmological parameter derivatives. The paper showcases results for the matter power spectrum, halo mass function, halo power spectrum, Voronoi volume function, kNN statistics, and environment-dependent halo correlations, and it argues that the resolution opens up access to low-mass halos relevant to DESI BGS at z<0.15. Data products are advertised as publicly available with a custom compression scheme.

Significance. If the numerical specifications hold, Sahyadri is a genuinely useful resource: its particle mass is roughly a factor 25 lower than AbacusSummit and more than two orders of magnitude below Quijote/Aemulus, while retaining a non-trivial simulation volume. The paper ships a public pipeline, includes consistency checks against halofit and Tinker, and demonstrates that high-density tracer samples enable measurements of beyond-2-point statistics such as the VVF and kNN-CDFs. These are concrete strengths. However, the two headline claims that would elevate the paper's significance—(i) reliable cosmological parameter derivatives for small-scale halo statistics and (ii) access to ~40% of DESI BGS galaxies—are not demonstrated in the manuscript: the first relies on an unpublished companion paper, and the second relies on an unpublished methodology. The central simulation product itself appears sound, but the advertised capabilities are currently only partly substantiated.

major comments (4)
  1. [§2.1.1, §3.2] The abstract and §2.1.1 state that seed-matched initial conditions enable reliable estimation of derivatives of small-scale observables involving halos, but the same paragraph concedes that seed matching is not sufficient to control gravitationally induced stochasticity for halo statistics and defers to the unpublished post-processing method of [24]. In the results sections, no derivative is computed for any statistic; Figures 8, 10, and 11 show raw Ωm responses with jackknife errors described as 'for illustrative purpose only.' The central advertised capability is therefore not demonstrated in this paper and depends entirely on an external, unpublished method. Please either include derivative estimates validated at Sahyadri's resolution/volume, or explicitly recast the derivative claim as a design goal to be established in a companion paper.
  2. [Abstract, §2, §4] The abstract and storage estimates describe the suite as containing 'systematic variations of six cosmological parameters,' and the 117 TB footprint is computed for the full six-parameter set. However, §2 and §4 state that only Ωm, h, and ns variations are complete, with As, w∥, and Ωk 'in progress.' This inconsistency overstates the presently available data products and should be corrected in the abstract, the storage calculations, and the summary section.
  3. [§1, Figure 1] The headline estimate that ~40% of DESI BGS galaxies at z<0.15 are accessible to Sahyadri but not AbacusSummit is based on the unpublished 'Alam et al. (in prep)' methodology, applied to the UniverseMachine SHMR and the Bell et al. mass-to-light mapping. The paper gives no details of the BGS selection model, the abundance matching implementation, or validation against observed BGS clustering/luminosity functions. Because this estimate is used to justify the entire survey-motivation of the suite, it needs to be either fully specified in an appendix or released as a public companion document before it can be assessed.
  4. [§4 vs Appendix C] The summary in §4 states that finite-volume effects on linear halo bias become negligible (≲1%) for halos with M≳10^12 h^-1 M_sun, but Appendix C shows that while the mean halo-by-halo bias is unbiased, the median bias is underestimated by Δb∼0.2 and the scatter σ_b by ∼0.7 in the fiducial 200 h^-1 Mpc box. These effects truncate the tails of the b1 distribution that are used in the environment-correlation analyses (§3.2, Figures 12–13). The text acknowledges that Spearman correlations 'may slightly weaken' the trends, but this is in tension with the 'negligible' summary statement. Please reconcile these statements or quantify the impact on the environment results.
minor comments (4)
  1. [§2.1.2] The compression scheme is claimed to maintain 'sub-percent clustering accuracy,' but no direct test of compressed versus uncompressed power spectra or correlation functions is presented. A small validation figure or quantitative statement would substantiate this claim.
  2. [§3.2, Appendix A] The text in §3.2 says the matter power spectrum shows 'good agreement' with halofit and deviations begin for k≳4 h/Mpc, but Appendix A reports deviations up to ∼20% at z=0 over a wide range of scales. Please make the characterization consistent and state explicitly whether this level of agreement is expected for the simulation resolution and halofit's calibration.
  3. [§2.2, Appendix C] Typos and minor wording issues: 'reproduceability' (§2.2), 'while while' (Appendix C), and 'spectrua' (§1).
  4. [References] Reference [24] and 'Alam et al. (in prep)' are central to the paper's claims but appear only as unpublished works. If they are not yet on arXiv, please provide versions or at least a detailed public methodological appendix for the referee and readers.

Circularity Check

0 steps flagged

No significant circularity: simulation outputs are benchmarked against external fitting functions and prior suites, and the Omega_m responses are direct measurements rather than reductions of their own inputs.

full rationale

Sahyadri's central results are numerical experiments: 2048^3 N-body runs whose matter/halo power spectra and mass functions are compared with halofit and Tinker, and whose Omega_m sensitivities are obtained by re-running with varied input cosmologies. No parameter is fitted to the same statistic it is used to predict; the Omega_m responses in Figs. 7-11 are raw simulation ratios, not outputs of a model calibrated on those ratios. The resolution claims (m_p = 8.1e7 Msun/h, M_min = 3.2e9 Msun/h at 40 particles) follow by arithmetic from box size, particle count and a standard 40-particle threshold. The BGS fraction (~40%) is an application estimate using external relations (UniverseMachine, Bell et al.) plus the unpublished Alam et al. methodology; it is not in the simulation's derivation chain, though it is a verifiability concern. The paper does lean on prior self-authored work: [24] is invoked for stabilizing small-scale derivative estimates, and the text explicitly concedes in Sec. 2.1.1 that seed matching alone 'is not sufficient to control the gravitationally induced stochasticity due to phase-mixing at small scales, which affects statistical probes involving halos', deferring to [24]. [48,49,63] provide the VVF/tidal-anisotropy methods. These are method and companion-paper citations, not self-defining reductions of central measurements to their own inputs. The deferred derivative capability is a limitation—no derivative is actually computed here—but absence of demonstration is not circularity. Core validation is self-contained against external benchmarks, so circularity score is 0.

Axiom & Free-Parameter Ledger

6 free parameters · 8 axioms · 0 invented entities

The central claims rest on standard N-body modeling choices, hand-chosen parameter steps and sample cuts, external fitting functions used as benchmarks, and unpublished galaxy-halo connection models for the BGS-impact estimate. No new physical entities are introduced; the only new object is the simulation suite itself, which is a data product rather than an invented physical entity.

free parameters (6)
  • Omega_m variation step = dOmega_m = 0.05 Omega_m,f ~ 0.0157
    Hand-chosen finite-difference step in Table 2; all Omega_m sensitivity plots and derivative claims depend on this step size.
  • 40-particle halo mass threshold = M_min = 40 mp = 3.2e9 h^-1 M_sun
    Hand-chosen standard resolution threshold; defines the headline mass-resolution claim and all halo samples.
  • Halo relaxation (QE) cut = eta = 2T/|U| in [0.5, 1.5]
    Adopted from Bett et al. (2007); every analysis uses halos satisfying this cut, so it shapes all quantitative results.
  • Tracer number densities for Vpeak-selected samples = 2e-4, 2e-3, 2e-2 Mpc^-3
    Chosen by hand in Section 2.3; the VVF and kNN analyses are performed for these fixed densities.
  • Vpeak thresholds at z=0 (fiducial) = 422.5, 196.6, 74.6 km/s
    Obtained by inverting the Vpeak distribution to match the chosen number densities; in effect fitted to target sample densities.
  • Tidal-anisotropy smoothing scale = R = 4 R200b
    Hand-chosen scale defining alpha; the environment-property correlation results in Figures 12-13 depend on it.
axioms (8)
  • domain assumption FLRW background plus Newtonian N-body dynamics as implemented in gadget-4 is a valid model of dark-matter structure formation from z=49 to z=0.
    Used throughout Section 2; all measurements inherit this modeling framework.
  • domain assumption 2LPT initial conditions at z=49 with the CLASS linear power spectrum are sufficient for the resolved scales.
    Section 2.1.1; no convergence test for initial-condition transients is presented.
  • domain assumption halofit and the Tinker mass function are appropriate external benchmarks for these cosmologies and mass ranges.
    Appendix A compares against them; deviations are attributed to the fits rather than to the simulation.
  • domain assumption Rockstar and consistent-trees correctly identify halos and merger trees at 40+ particles, and the QE/parent-halo cuts do not bias the samples.
    Section 2.2-2.3; halo samples are the basis of the mass-function and assembly-bias results.
  • domain assumption UniverseMachine abundance matching and the Bell et al. mass-to-light relation describe DESI BGS galaxies at z~0.15.
    Section 1 and Figure 1; underlies the 40%-of-BGS claim and depends on the unpublished Alam et al. companion paper.
  • domain assumption Seed-matched initial conditions isolate cosmological parameter dependence despite small-scale gravitational stochasticity.
    Section 2.1.1 acknowledges the limitation; the Omega_m comparisons rely on this assumption.
  • ad hoc to paper The w_parallel parametrization along the CMB degeneracy with A=-3.68 is a suitable dark-energy model for derivative calculations.
    Section 2 defines it via Lodha et al. (in prep); it is not independently justified in the paper.
  • domain assumption The custom compression scheme preserves clustering statistics at the sub-percent level.
    Section 2.1.2 asserts sub-percent accuracy but provides no detailed numerical validation of this claim.

pith-pipeline@v1.3.0-alltime-deepseek · 21929 in / 20162 out tokens · 184648 ms · 2026-08-03T10:59:19.531943+00:00 · methodology

0 comments
read the original abstract

We present Sahyadri, a suite of cosmological $N$-body simulations designed to enable precision studies of the low-redshift Universe with next-generation spectroscopic surveys. Sahyadri includes systematic variations of four cosmological parameters around Planck 2018 constraints, with seed-matched initial conditions enabling cosmological parameter derivatives. It is planned to ultimately extend to six parameters. Each simulation evolves $2048^3$ particles in a periodic box of side length $200$ $h^{-1}$ Mpc, yielding a particle mass of $m_{\rm{p}} = 8.1 \times 10^{7}\,h^{-1}\,M_{\odot}$ in the fiducial Planck 2018 cosmology. This resolution enables robust identification of dark matter halos down to $M_{\rm min} = 3.2 \times 10^{9}$ $h^{-1}$ $M_\odot$, which represents a factor of $\sim$25 improvement over the AbacusSummit suite, and is over two orders of magnitude better than the Quijote and Aemulus suites. We estimate that approximately 40% of DESI BGS galaxies at redshift $z < 0.15$ - roughly 1.6 million objects - reside in halos accessible to Sahyadri but beyond the reach of existing parameter-varying simulation suites. We demonstrate Sahyadri's capabilities through measurements of the matter power spectrum, halo mass function and power spectrum, and beyond 2-point statistics such as the Voronoi volume function and $k^{\rm th}$ nearest neighbour statistics, showing excellent agreement with theoretical predictions and significant sensitivity to $\Omega_{\rm m}$ variations. We implement a custom compression scheme reducing storage requirements by a factor of $\sim$3 while maintaining sub-percent clustering accuracy. Key data products are publicly available.

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

74 extracted references · 31 linked inside Pith

  1. [1]

    Aghamousa, J

    DESI Collaboration, A. Aghamousa, J. Aguilar, S. Ahlen et al.,The desi experiment part i: Science, targeting, and survey design,arXiv e-prints(2016) [1611.00036]

  2. [2]

    de Jong, O

    R.S. de Jong, O. Agertz, A.A. Berbel et al.,4most: Project overview and information for the first call for proposals,The Messenger175(2019) 3

  3. [3]

    Takada, R.S

    M. Takada, R.S. Ellis, M. Chiba et al.,Extragalactic science, cosmology, and galactic archaeology with the subaru prime focus spectrograph,Publications of the Astronomical Society of Japan66(2014) R1

  4. [4]

    Mellier, F.B

    Euclid Collaboration, Y. Mellier, F.B. Abdalla et al.,Euclid: Overview of the mission and survey,arXiv e-prints(2022) [2209.02534]

  5. [5]

    Ivezi´ c, S.M

    ˇZ. Ivezi´ c, S.M. Kahn, J.A. Tyson et al.,Lsst: From science drivers to reference design and anticipated data products,The Astrophysical Journal873(2019) 111

  6. [6]

    Ishiyama, F

    T. Ishiyama, F. Prada, A.A. Klypin et al.,The uchuu simulations: Data release 1 and dark matter halo concentrations,Monthly Notices of the Royal Astronomical Society506(2021) 4210

  7. [7]

    Dav´ e, D

    R. Dav´ e, D. Angl´ es-Alc´ azar, D. Narayanan et al.,Simba: Cosmological simulations with black hole growth and feedback,Monthly Notices of the Royal Astronomical Society486(2019) 2827

  8. [8]

    Schaye, R.A

    J. Schaye, R.A. Crain, R.G. Bower et al.,The eagle project: simulating the evolution and assembly of galaxies and their environments,Monthly Notices of the Royal Astronomical Society446(2015) 521

  9. [9]

    Vogelsberger, S

    M. Vogelsberger, S. Genel, V. Springel et al.,Introducing the illustris project: simulating the coevolution of dark and visible matter in the universe,Monthly Notices of the Royal Astronomical Society444(2014) 1518

  10. [10]

    Nelson, A

    D. Nelson, A. Pillepich, V. Springel et al.,The illustristng simulations: public data release, Computational Astrophysics and Cosmology6(2019) 2

  11. [11]

    Pillepich, D

    A. Pillepich, D. Nelson, L. Hernquist et al.,First results from the illustristng simulations: the galaxy colour bimodality,Monthly Notices of the Royal Astronomical Society475(2018) 648

  12. [12]

    Springel, S.D.M

    V. Springel, S.D.M. White, A. Jenkins et al.,Simulations of the formation, evolution and clustering of galaxies and quasars,Nature435(2005) 629

  13. [13]

    Boylan-Kolchin, V

    M. Boylan-Kolchin, V. Springel, S.D.M. White, A. Jenkins and G. Lemson,Resolving cosmic structure formation with the millennium-ii simulation,Monthly Notices of the Royal Astronomical Society398(2009) 1150

  14. [14]

    Klypin, S

    A.A. Klypin, S. Trujillo-Gomez and J. Primack,Dark matter halos in the standard cosmological model: Results from the bolshoi simulation,The Astrophysical Journal740(2011) 102

  15. [15]

    Klypin, G

    A. Klypin, G. Yepes, S. Gottl¨ ober, F. Prada and S. Heß,Multidark simulations: the story of dark matter halo concentrations and density profiles,Monthly Notices of the Royal Astronomical Society457(2016) 4340

  16. [16]

    Villaescusa-Navarro, D

    F. Villaescusa-Navarro, D. Angl´ es-Alc´ azar, S. Genel, D.N. Spergel, R.S. Somerville, R. Dave et al.,The CAMELS Project: Cosmology and Astrophysics with Machine-learning Simulations, ApJ915(2021) 71 [2010.00619]. 10https://matplotlib.org/ 11https://jupyter.org – 22 –

  17. [17]

    Villaescusa-Navarro, C

    F. Villaescusa-Navarro, C. Hahn, E. Massara et al.,The quijote simulations,The Astrophysical Journal Supplement Series250(2020) 2

  18. [18]

    DeRose, R.H

    J. DeRose, R.H. Wechsler, J.L. Tinker et al.,The aemulus project. ii. emulating the halo mass function,The Astrophysical Journal875(2019) 69

  19. [19]

    McClintock, T.N

    T. McClintock, T.N. Varga, D. Gruen et al.,Dark energy survey year 1 results: Weak lensing mass calibration of redmapper galaxy clusters,Monthly Notices of the Royal Astronomical Society482(2019) 1352

  20. [20]

    Nishimichi, A

    T. Nishimichi, A. Shirata, A. Taruya et al.,Dark quest. i. fast and accurate emulation of halo clustering statistics and its application to galaxy clustering,The Astrophysical Journal884 (2019) 29

  21. [21]

    Garrison, D.J

    L.H. Garrison, D.J. Eisenstein, D. Ferrer et al.,The abacus cosmos: A suite of cosmological n-body simulations,The Astrophysical Journal Supplement Series236(2018) 43

  22. [22]

    Maksimova, L.H

    N.A. Maksimova, L.H. Garrison, B. Hadzhiyska et al.,Abacussummit: a massive set of high-accuracy, high-resolution n-body simulations,Monthly Notices of the Royal Astronomical Society508(2021) 4017

  23. [23]

    Alam, J.A

    S. Alam, J.A. Peacock, D.J. Farrow, J. Loveday and A.M. Hopkins,Using gama to probe the impact of small-scale galaxy physics on nonlinear redshift-space distortions,Monthly Notices of the Royal Astronomical Society503(2021) 59–76

  24. [24]

    Dhawalikar, A

    S. Dhawalikar, A. Paranjape and S. Alam,Stabilizing simulation-based cosmological fisher forecasts: a case study using the voronoi volume function, 2025

  25. [25]

    Hahn, M.J

    C. Hahn, M.J. Wilson, O. Ruiz-Macias et al.,The desi bright galaxy survey: Final target selection, design, and validation,The Astronomical Journal165(2023) 253

  26. [26]

    Yildiz, M

    E. Yildiz, M. Biagetti, U. Seljak et al.,Towards testing the theory of gravity with desi: summary statistics, model predictions and future simulation requirements,Monthly Notices of the Royal Astronomical Society493(2020) 4783

  27. [27]

    Behroozi, R.H

    P. Behroozi, R.H. Wechsler, A.P. Hearin and C. Conroy,Universemachine: The correlation between galaxy growth and dark matter halo assembly from z = 0-10,Monthly Notices of the Royal Astronomical Society488(2019) 3143

  28. [28]

    Bell, D.H

    E.F. Bell, D.H. McIntosh, N. Katz and M.D. Weinberg,The optical and near-infrared properties of galaxies. i. luminosity and stellar mass functions,The Astrophysical Journal Supplement Series149(2003) 289

  29. [29]

    Driver, D.T

    S.P. Driver, D.T. Hill, L.S. Kelvin et al.,Galaxy and mass assembly (gama): survey diagnostics and core data release,Monthly Notices of the Royal Astronomical Society413(2011) 971

  30. [30]

    Baldry, S.P

    I.K. Baldry, S.P. Driver, J. Loveday et al.,Galaxy and mass assembly (gama): the galaxy stellar mass function at z ¡ 0.06,Monthly Notices of the Royal Astronomical Society421(2012) 621

  31. [31]

    Wetzel, J.L

    A.R. Wetzel, J.L. Tinker, C. Conroy and F.C. van den Bosch,Galaxy evolution in groups and clusters: Satellite star formation histories and quenching time-scales in a hierarchical universe, Monthly Notices of the Royal Astronomical Society432(2013) 336

  32. [32]

    Reines, J.E

    A.E. Reines, J.E. Greene and M. Geha,Dwarf galaxies with optical signatures of active massive black holes,The Astrophysical Journal775(2013) 116

  33. [33]

    Riess, W

    A.G. Riess, W. Yuan, L.M. Macri et al.,A comprehensive measurement of the local value of the hubble constant with 1 km s −1 mpc−1 uncertainty from the hubble space telescope and the sh0es team,The Astrophysical Journal Letters934(2022) L7

  34. [34]

    Aghanim, Y

    Planck Collaboration, N. Aghanim, Y. Akrami et al.,Planck 2018 results. vi. cosmological parameters,Astronomy & Astrophysics641(2020) A6. – 23 –

  35. [35]

    Adame, J

    DESI Collaboration, A.G. Adame, J. Aguilar, S. Ahlen, S. Alam, D.M. Alexander et al.,DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations, arXiv e-prints(2024) arXiv:2404.03002 [2404.03002]

  36. [36]

    Abdul Karim, J

    M. Abdul Karim, J. Aguilar, S. Ahlen, S. Alam, L. Allen, C.A. Prieto et al.,Desi dr2 results. ii. measurements of baryon acoustic oscillations and cosmological constraints,Physical Review D112(2025)

  37. [37]

    Wadekar and R

    D. Wadekar and R. Scoccimarro,Galaxy power spectrum multipoles covariance in perturbation theory,Physical Review D102(2020) 123521

  38. [38]

    Q. Hang, S. Alam, J.A. Peacock and Y.-C. Cai,Galaxy clustering in the desi legacy survey and its imprint on the cmb,Monthly Notices of the Royal Astronomical Society501(2020) 1481–1498

  39. [39]

    Ishiyama, F

    T. Ishiyama, F. Prada and A.A. Klypin,Evolution of clustering in cosmological models with time-varying dark energy,Phys. Rev. D112(2025) 043504

  40. [40]

    Springel, R

    V. Springel, R. Pakmor, O. Zier and M. Reinecke,Simulating cosmic structure formation with the GADGET-4 code, MNRAS506(2021) 2871 [2010.03567]

  41. [41]

    Scoccimarro,Transients from initial conditions: a perturbative analysis, MNRAS299 (1998) 1097 [astro-ph/9711187]

    R. Scoccimarro,Transients from initial conditions: a perturbative analysis, MNRAS299 (1998) 1097 [astro-ph/9711187]

  42. [42]

    D. Blas, J. Lesgourgues and T. Tram,The cosmic linear anisotropy solving system (class). part ii: Approximation schemes,Journal of Cosmology and Astroparticle Physics2011(2011) 034

  43. [43]

    Lesgourgues,The cosmic linear anisotropy solving system (class) i: Overview,arXiv e-prints (2011) [1104.2932]

    J. Lesgourgues,The cosmic linear anisotropy solving system (class) i: Overview,arXiv e-prints (2011) [1104.2932]

  44. [44]

    Coulton, F

    W.R. Coulton, F. Villaescusa-Navarro, D. Jamieson, M. Baldi, G. Jung, D. Karagiannis et al., Quijote-PNG: The Information Content of the Halo Power Spectrum and Bispectrum, ApJ943 (2023) 178 [2206.15450]

  45. [45]

    Maksimova, L.H

    N.A. Maksimova, L.H. Garrison, B. Hadzhiyska, S. Bose and D.J. Eisenstein,Abacussummit: a massive set of high-accuracy, high-resolution n-body simulations,Monthly Notices of the Royal Astronomical Society508(2021) 4017

  46. [46]

    Behroozi, R.H

    P.S. Behroozi, R.H. Wechsler and H.-Y. Wu,The ROCKSTAR Phase-space Temporal Halo Finder and the Velocity Offsets of Cluster Cores, ApJ762(2013) 109 [1110.4372]

  47. [47]

    Behroozi, R.H

    P.S. Behroozi, R.H. Wechsler, H.-Y. Wu, M.T. Busha, A.A. Klypin and J.R. Primack, Gravitationally Consistent Halo Catalogs and Merger Trees for Precision Cosmology, ApJ763 (2013) 18 [1110.4370]

  48. [48]

    Paranjape, O

    A. Paranjape, O. Hahn and R.K. Sheth,Halo assembly bias and the tidal anisotropy of the local halo environment, MNRAS476(2018) 3631 [1706.09906]

  49. [49]

    Paranjape and S

    A. Paranjape and S. Alam,Voronoi volume function: a new probe of cosmology and galaxy evolution, MNRAS495(2020) 3233 [2001.08760]

  50. [50]

    P. Bett, V. Eke, C.S. Frenk, A. Jenkins, J. Helly and J. Navarro,The spin and shape of dark matter haloes in the Millennium simulation of aΛcold dark matter universe, MNRAS376 (2007) 215 [astro-ph/0608607]

  51. [51]

    Reddick, R.H

    R.M. Reddick, R.H. Wechsler, J.L. Tinker and P.S. Behroozi,The Connection between Galaxies and Dark Matter Structures in the Local Universe, ApJ771(2013) 30 [1207.2160]

  52. [52]

    Carretero, F.J

    J. Carretero, F.J. Castander, E. Gazta˜ naga, M. Crocce and P. Fosalba,An algorithm to build mock galaxy catalogues using MICE simulations, MNRAS447(2015) 646 [1411.3286]

  53. [53]

    Chaves-Montero, R.E

    J. Chaves-Montero, R.E. Angulo, J. Schaye, M. Schaller, R.A. Crain, M. Furlong et al.,Subhalo abundance matching and assembly bias in the EAGLE simulation, MNRAS460(2016) 3100 [1507.01948]. – 24 –

  54. [54]

    DeRose, M.R

    J. DeRose, M.R. Becker and R.H. Wechsler,Modeling Redshift-space Clustering with Abundance Matching, ApJ940(2022) 13 [2105.12104]

  55. [55]

    Takahashi, M

    R. Takahashi, M. Sato, T. Nishimichi, A. Taruya and M. Oguri,Revising the Halofit Model for the Nonlinear Matter Power Spectrum, ApJ761(2012) 152 [1208.2701]

  56. [56]

    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, ApJ688(2008) 709 [0803.2706]

  57. [57]

    Voronoi,Nouvelles applications des param` etres continus ` a la th´ eorie des formes quadratiques

    G. Voronoi,Nouvelles applications des param` etres continus ` a la th´ eorie des formes quadratiques. deuxi` eme m´ emoire. recherches sur les parall´ ello` edres primitifs.,Journal f¨ ur die reine und angewandte Mathematik134(1908) 198

  58. [58]

    Banerjee and T

    A. Banerjee and T. Abel,Nearest neighbour distributions: New statistical measures for cosmological clustering,Mon. Not. Roy. Astron. Soc.500(2020) 5479 [2007.13342]

  59. [59]

    Banerjee and T

    A. Banerjee and T. Abel,Cosmological cross-correlations and nearest neighbour distributions, Mon. Not. Roy. Astron. Soc.504(2021) 2911 [2102.01184]

  60. [60]

    Banerjee and T

    A. Banerjee and T. Abel,Tracer-field cross-correlations with k-nearest neighbour distributions, Mon. Not. Roy. Astron. Soc.519(2023) 4856 [2210.05140]

  61. [61]

    Gangopadhyay, A

    K. Gangopadhyay, A. Banerjee and T. Abel,Geometric Interpretations of thek-Nearest Neighbour Distributions,Mon. Not. Roy. Astron. Soc.3409(2025) 3427 [2502.07713]

  62. [62]

    Banerjee, N

    A. Banerjee, N. Kokron and T. Abel,Modelling nearest neighbour distributions of biased tracers using hybrid effective field theory,Mon. Not. Roy. Astron. Soc.511(2022) 2765 [2107.10287]

  63. [63]

    Ramakrishnan, A

    S. Ramakrishnan, A. Paranjape, O. Hahn and R.K. Sheth,Cosmic web anisotropy is the primary indicator of halo assembly bias, MNRAS489(2019) 2977 [1903.02007]

  64. [64]

    Ramakrishnan, V

    S. Ramakrishnan, V. Gonzalez-Perez, G. Parimbelli and G. Yepes,The multi-dimensional halo assembly bias can be preserved when enhancing halo properties with HALOSCOPE, A&A697 (2025) A70 [2410.07361]

  65. [65]

    Ramakrishnan and P

    S. Ramakrishnan and P. Velmani,Properties beyond mass for unresolved haloes across redshift and cosmology using correlations with local halo environment, MNRAS516(2022) 5849 [2112.15305]

  66. [66]

    Ramakrishnan, A

    S. Ramakrishnan, A. Paranjape and R.K. Sheth,Mock halo catalogues: assigning unresolved halo properties using correlations with local halo environment, MNRAS503(2021) 2053 [2012.10170]

  67. [67]

    Van Der Walt, S.C

    S. Van Der Walt, S.C. Colbert and G. Varoquaux,The NumPy array: a structure for efficient numerical computation,ArXiv e-prints(2011) [1102.1523]

  68. [68]

    Virtanen, R

    P. Virtanen, R. Gommers, T.E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau et al., SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python,Nature Methods17 (2020) 261

  69. [69]

    56 – 61, 2010, DOI

    Wes McKinney,Data Structures for Statistical Computing in Python, inProceedings of the 9th Python in Science Conference, St´ efan van der Walt and Jarrod Millman, eds., pp. 56 – 61, 2010, DOI

  70. [70]

    pandas development team,pandas-dev/pandas: Pandas, Feb., 2020

    T. pandas development team,pandas-dev/pandas: Pandas, Feb., 2020. 10.5281/zenodo.3509134

  71. [71]

    Hunter,Matplotlib: A 2d graphics environment,Computing In Science & Engineering9 (2007) 90

    J.D. Hunter,Matplotlib: A 2d graphics environment,Computing In Science & Engineering9 (2007) 90

  72. [72]

    van Daalen, J

    M.P. van Daalen, J. Schaye, C.M. Booth and C. Dalla Vecchia,The effects of galaxy formation on the matter power spectrum: a challenge for precision cosmology, MNRAS415(2011) 3649 [1104.1174]. – 25 –

  73. [73]

    Tinker, B.E

    J.L. Tinker, B.E. Robertson, A.V. Kravtsov, A. Klypin, M.S. Warren, G. Yepes et al.,The Large-scale Bias of Dark Matter Halos: Numerical Calibration and Model Tests, ApJ724 (2010) 878 [1001.3162]

  74. [74]

    N. Paul, I. Pahwa and A. Paranjape,Global analysis of luminosity- and colour-dependent galaxy clustering in the Sloan Digital Sky Survey, MNRAS488(2019) 1220 [1902.01566]. A Comparison with literature fitting functions As a consistency check, we compare the matter power spectrum and halo mass function measured from our simulations with widely used fitting...