Aspects of gravitational clustering and structure formation in the Universe

Swati Gavas

Aspects of gravitational clustering and structure formation in the Universe

T0 review · 3 major / 2 minor · reviewed 2026-05-08 · grok-4.3

Challenge this review Download PDF Read on arXiv ↗

Pith's one-line read The halo mass function depends explicitly on the slope of the input power spectrum, and local density variations bias Hubble constant estimates by up to 5 percent in Milky Way-sized halos.

desk verdict The paper shows explicit slope dependence for the HMF in scale-free EdS runs and up to 5% local H0 shifts in MW halos, but the effective-index bridge to full ΛCDM is untested. read the letter →

arxiv 2604.21634 v1 submitted 2026-04-23 astro-ph.CO

Swati Gavas This is my paper · ORCID

classification astro-ph.CO

keywords halomassfunctionpowerspectrumslopeN-bodysimulationsHubbletensionstructureformationmodecouplingcosmologicalinhomogeneities

topics The Hubble Tension

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

N-body simulations of structure formation show that the halo mass function deviates systematically from universal theoretical predictions because it carries an explicit dependence on the slope of the initial power spectrum. In Einstein-de Sitter cosmologies with scale-free spectra, this dependence allows an effective index to serve as a first approximation to the full Lambda CDM case. Non-linear clustering also produces deviations from homogeneity on scales of at least 100 Mpc/h that correlate with errors in local Hubble constant measurements. Mock observations placed in Milky Way-sized halos yield shifts of up to 5 percent, suggesting the effect could account for part of the Hubble tension. Insufficient resolution of small-scale halos in the simulations further introduces errors in root-mean-square fluctuations and the mass function through incomplete mode coupling between scales.

What carries the argument

The explicit dependence of the halo mass function on the slope of the input power spectrum, combined with the correlation between local density and errors in Hubble constant estimates.

What would settle it

Compare halo mass functions measured in N-body simulations that use power spectra with different slopes but the same effective index to test whether the claimed explicit dependence appears; or measure the correlation between local galaxy density and Hubble constant estimates in real surveys.

Watch

Extended reading notes

Core claim

For scale-free power spectra in an Einstein-de Sitter cosmology, the halo mass function is not universal but depends explicitly on the slope of the input power spectrum. An effective index drawn from the Lambda CDM power spectrum can therefore reproduce the scale-free halo mass function as a first approximation. Structure formation drives deviations from homogeneity and isotropy out to at least 100 Mpc/h, which produce errors in Hubble constant estimates that correlate with the local density around the observer. Mock observations show deviations reaching 5 percent in Milky Way-sized halos, indicating that this local effect may explain part of the Hubble tension. Missing small-scale power in

Load-bearing premise

Scale-free power spectra in an Einstein-de Sitter cosmology provide a valid first approximation to the halo mass function of the full Lambda CDM model via an effective index, and the N-body simulations resolve enough small-scale halos for accurate mode coupling.

Editorial extensions

If this is right

Theoretical models of halo abundance must incorporate the slope of the input power spectrum rather than assuming universality.
An effective spectral index allows scale-free simulations to approximate halo statistics in the full Lambda CDM cosmology.
Local density variations around observers can introduce systematic errors of several percent in Hubble constant measurements.
Accurate mode coupling between small and large scales requires simulations to resolve collapsed halos at small scales.

Reading between the lines

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

Observer location within the cosmic web may contribute to apparent discrepancies between local and global Hubble constant measurements.
Higher-resolution simulations could test whether the effective-index approximation continues to hold when the power spectrum includes baryonic features and a cosmological constant.
Accounting for local density in analyses of large-scale structure surveys could reduce reported tensions without invoking new physics.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

The main takeaways are straightforward. In scale-free Einstein-de Sitter simulations the halo mass function deviates from universality in a way that tracks the input power-spectrum slope, and mock local observations find density variations around Milky Way-sized halos that shift inferred H0 by as much as 5%. An analytical mode-coupling calculation is used to argue that small-scale resolution matters for accurate large-scale power transfer. These are concrete outputs from the N-body work rather than abstract claims. The slope dependence and the local H0 correlation are the clearest new pieces; both extend earlier simulation studies but give specific numbers for these setups. The paper also flags the practical limit that missing small-scale power in the runs affects rms fluctuations and the mass function itself. That acknowledgment is useful. The weakest link is the step from scale-free results to ΛCDM. The authors propose that an effective index can serve as a first approximation, yet no side-by-side comparison to a full ΛCDM run at matched mass and redshift is reported. Without that check, differences from the actual transfer function, baryon effects, or post-z=1 growth stay unquantified. The abstract gives 5-20% deviation ranges without error bars or resolution details, though the full text may supply them. This is a moderate rather than fatal gap for a methods-focused paper. The work is aimed at people already running or analyzing N-body simulations who care about HMF non-universality and possible local contributions to the Hubble tension. It is not a new framework, so most readers will treat it as an incremental check rather than a must-read. I would bring the H0 section to a reading group for discussion but would not cite the paper in the next year unless the effective-index validation is added. A serious editor should still send it to referees; the topic is relevant and the simulation results are worth checking, provided the authors address the direct ΛCDM comparison and convergence tests.

Referee Report

3 major / 2 minor

Summary. The manuscript uses N-body simulations of scale-free power spectra in an Einstein-de Sitter cosmology to study non-linear gravitational clustering. It claims that the halo mass function (HMF) shows an explicit dependence on the slope of the input power spectrum, that an effective spectral index supplies a first approximation to the ΛCDM HMF, and that mock observations of local structure yield up to 5% deviations in H0 estimates for Milky Way-sized halos due to density variations around the observer, potentially contributing to the Hubble tension. The work also analyzes limitations of N-body simulations arising from unresolved small-scale power and its effect on mode coupling.

Significance. If the effective-index mapping and the quantified H0 bias hold under direct validation, the results would strengthen understanding of HMF non-universality and local cosmological biases. The controlled use of scale-free EdS runs combined with an analytical mode-coupling calculation is a methodological strength that isolates spectral-index effects cleanly. The 5% H0 deviation estimate, if robust, would be a concrete, falsifiable prediction relevant to ongoing Hubble-tension analyses.

major comments (3)

[Abstract and simulation section] Abstract and simulation section: The reported 5-20% systematic deviations of the simulated HMF from theoretical predictions are presented without error bars, convergence tests, or resolution details (particle number, box size, force softening), preventing assessment of whether the claimed explicit dependence on power-spectrum slope is statistically significant or resolution-limited.
[Effective-index mapping] Effective-index mapping: The statement that an effective index of the ΛCDM model can correspond to the HMF from scale-free cosmologies is offered as a first approximation, yet the manuscript contains no side-by-side comparison of this mapped HMF against a full ΛCDM N-body run at matched halo masses and redshifts; differences due to the k-dependent transfer function, baryonic suppression, or post-z≈1 growth remain unquantified and directly affect both the universality-breaking claim and the subsequent 5% H0 bias estimate.
[H0 mock-observation analysis] H0 mock-observation analysis: The correlation between H0 errors and local density, together with the 5% deviation claim for Milky Way-sized halos, is presented without an explicit description of the mock-observation pipeline, the fitting procedure for H0, or the sample variance across multiple observer positions, making it impossible to judge whether the reported bias is load-bearing or an artifact of the chosen observer selection.

minor comments (2)

[Abstract] The abstract refers to 'the concordance model' without a brief definition or reference in the introduction, which could confuse readers unfamiliar with the specific context of H0 measurements.
[Figures and tables] Figure captions and table legends should explicitly list the spectral indices (n) used in each scale-free run and the precise definition of the effective index adopted for the ΛCDM comparison.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We will revise the manuscript to improve clarity on methodological details and strengthen the presentation of our results. We respond to each major comment below.

read point-by-point responses

Referee: [Abstract and simulation section] Abstract and simulation section: The reported 5-20% systematic deviations of the simulated HMF from theoretical predictions are presented without error bars, convergence tests, or resolution details (particle number, box size, force softening), preventing assessment of whether the claimed explicit dependence on power-spectrum slope is statistically significant or resolution-limited.

Authors: We agree that error bars, convergence tests, and full resolution specifications are required to establish the statistical significance of the reported HMF deviations and their dependence on the input power-spectrum slope. In the revised manuscript we will add error bars computed from multiple independent realizations, present convergence tests varying particle number and box size, and explicitly state the simulation parameters (particle number, box size, force softening). These additions will confirm that the slope dependence is not an artifact of limited resolution. revision: yes
Referee: [Effective-index mapping] Effective-index mapping: The statement that an effective index of the ΛCDM model can correspond to the HMF from scale-free cosmologies is offered as a first approximation, yet the manuscript contains no side-by-side comparison of this mapped HMF against a full ΛCDM N-body run at matched halo masses and redshifts; differences due to the k-dependent transfer function, baryonic suppression, or post-z≈1 growth remain unquantified and directly affect both the universality-breaking claim and the subsequent 5% H0 bias estimate.

Authors: The effective-index mapping is presented strictly as a first approximation that isolates the leading effect of the local spectral slope within our controlled scale-free EdS simulations. A direct side-by-side comparison against a full ΛCDM run that includes the k-dependent transfer function, baryons, and post-z≈1 growth lies outside the scope of the present work. In the revision we will expand the discussion to explicitly quantify the limitations of the approximation and to state that the primary evidence for non-universality is the explicit slope dependence demonstrated in the scale-free runs themselves. revision: partial
Referee: [H0 mock-observation analysis] H0 mock-observation analysis: The correlation between H0 errors and local density, together with the 5% deviation claim for Milky Way-sized halos, is presented without an explicit description of the mock-observation pipeline, the fitting procedure for H0, or the sample variance across multiple observer positions, making it impossible to judge whether the reported bias is load-bearing or an artifact of the chosen observer selection.

Authors: We acknowledge that the H0 analysis section requires a more complete description of the methodology. In the revised manuscript we will provide a detailed account of the mock-observation pipeline, the precise fitting procedure used to extract H0, the criteria for selecting observer positions inside Milky Way-sized halos, and quantitative estimates of sample variance obtained from multiple observer locations and independent realizations. These additions will allow readers to assess the robustness of the reported correlation and the 5% deviation amplitude. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: HMF slope dependence and effective-index mapping derived from scale-free EdS simulations plus analytical mode-coupling, not by construction from inputs

full rationale

The paper's core results on explicit power-spectrum slope dependence in the HMF and the effective-index approximation to ΛCDM are obtained from direct N-body runs in scale-free EdS cosmologies plus an analytical calculation of mode coupling. These outputs do not reduce to fitted parameters renamed as predictions, self-definitions, or self-citation chains; the mapping to ΛCDM is presented as a first approximation without claiming it is forced by the scale-free data alone. The derivation chain remains self-contained against the reported simulation outputs and mode-coupling analysis.

Assumptions & free parameters 1 free parameters · 2 assumptions · 0 invented entities

The central claims rest on standard cosmological N-body assumptions and the validity of scale-free approximations; no new entities are postulated.

free parameters (1)

effective index
Chosen to match ΛCDM halo mass function with scale-free results; appears fitted rather than derived from first principles.

assumptions (2)

domain assumption Einstein-de Sitter background for scale-free power spectra
Invoked to isolate power-spectrum-slope effects without dark energy.
domain assumption N-body simulations faithfully capture non-linear gravitational clustering and mode coupling when halos are resolved
Underlies all quantitative HMF and H0 claims.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Aspects of gravitational clustering and structure formation in the Universe." pith.science (2026). https://pith.science/paper/2604.21634

@misc{pith2026260421634,
  author       = {Pith},
  title        = {Pith review of: Aspects of gravitational clustering and structure formation in the Universe},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2604.21634}},
  note         = {Machine review of arXiv:2604.21634}
}

abstract

The distribution of galaxies, halo abundance, and peculiar velocities are influenced by non-linear gravitational interactions, making the study of non-linear evolution crucial for accurate cosmological predictions. We explore these aspects using N-body simulations. Theoretical models of the halo mass function (HMF) can be formulated without referencing a cosmological model or input power spectrum. HMF obtained from N-body simulations show systematic deviations of 5-20\% from theoretical predictions. The physical origin of deviations may result from cosmology, the power spectrum, or both. We examine HMF deviations from universality for scale-free power spectra with an Einstein-de Sitter cosmology. We demonstrate that the mass function exhibits an explicit dependence on the slope of the input power spectrum. We find that an effective index of the $\Lambda$CDM model can correspond to the HMF from scale-free cosmologies as a first approximation. Furthermore, structure formation has led to deviations from homogeneity and isotropy on scales up to at least $100$ Mpc/h, expected to affect measurements of $H_0$. We revisit this issue of the concordance model. We find a correlation between errors in $H_0$ estimates and the density around the observer. Further, our mock observations reveal that deviations of up to 5\% can occur in Milky Way-sized halos. While this finding alone does not fully resolve the Hubble tension, it may account for part of it. It is essential to understand the limitations of N-body simulations to avoid misinterpreting data. We show that the missing power at small scales introduces errors in the root-mean-square fluctuations and in the simulated mass function. Our analytical calculation indicates that mode coupling between small and large scales depends on resolving collapsed halos. Therefore, accurate mode coupling estimates require sufficient halos in the simulation.

Discussion (0). Continue with ORCID to comment.

Pith tools

Reviewed May 8, 2026 · model on record in the stance chip above.