Pith. sign in

REVIEW 2 major objections 5 minor 2 cited by

The minimum spanning tree of galaxy positions carries a neutrino-mass signal that two-point statistics miss.

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 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Massive neutrinos shrink clusters and voids, thicken filaments, and leave distinct signatures in minimum-spanning-tree edge lengths that differ from baryonic feedback effects.

T0 review reviewed 2026-08-03 challenge →

load-bearing objection Solid proof-of-concept: MST edge-length sensitivity to neutrinos is localized to filaments, but single-realization errors make the quantitative claims provisional. the 2 major comments →

arxiv 2512.16517 v2 pith:W7YKOOMA submitted 2025-12-18 astro-ph.CO

How massive neutrinos reshape the cosmic web

Leonor N. L. Sim\~oes , Krishna Naidoo , Benjamin Joachimi , Willem Elbers , Carlos S. Frenk This is my paper
classification astro-ph.CO
keywords massive neutrinoscosmic webminimum spanning treefilamentslarge-scale structureNEXUS+neutrino mass constraints
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 reading

The paper argues that the minimum spanning tree (MST) of subhalo positions — the shortest network connecting all points — is sensitive to the sum of neutrino masses because roughly 70% of its edges lie in cosmic-web filaments. It shows that neutrinos and baryons imprint opposite, distinguishable signatures on the MST edge-length distribution: neutrinos suppress short edges near 1 Mpc, boost intermediate edges near 3–4 Mpc, and cut long edges near 9 Mpc, while baryons mainly boost the shortest edges. Using a NEXUS+ classification of the density field, the paper traces each scale of the signal to a specific environment — clusters at small scales, voids at large scales, filaments in between. If correct, the MST offers a complementary, computationally cheap probe for breaking the neutrino–baryon degeneracy that plagues power-spectrum analyses, with direct application to galaxy-survey-like densities.

Core claim

The central discovery is that the MST edge-length statistic encodes neutrino mass through its cosmic-web environment content. Increasing neutrino mass reduces MST edges at ~1 Mpc (linked to clusters and sheets), increases them at ~3–4 Mpc (linked to filaments), and decreases them at ~9 Mpc (linked to voids), while baryonic physics increases edges at ~1 Mpc — the opposite of the neutrino effect at small scales. The paper attributes this to neutrinos delaying structure formation, making clusters less dense and voids shallower, and shows that these distinct signatures are the likely origin of the MST's previously reported neutrino-mass sensitivity.

What carries the argument

The argument runs on two coupled objects: (1) the NEXUS+ cosmic-web classification, which assigns each volume element to void, sheet, filament, or cluster by taking the maximum Hessian-based shape signature over five Gaussian smoothing scales; and (2) the minimum spanning tree (MST) built on 500,000 randomly selected subhaloes above a mass threshold, whose edge-length distribution is the observable. To connect them, each MST edge is split into segments and assigned to environments by straight-line sampling through the NEXUS+ field. The edge-length histogram, split by environment, is what carries the neutrino-baryon discrimination.

Load-bearing premise

The NEXUS+ cosmic-web classifier, with its chosen smoothing scales and threshold criteria, is assumed to capture the true physical environments in a way that makes straight-line MST edge segments reliably trace those same environments; if a different classification or a non-straight path assignment changed the ~70% filament fraction or the scale-environment mapping, the link between MST edge lengths and neutrino mass could weaken.

What would settle it

A direct test: repeat the MST edge-length analysis on smaller neutrino masses (e.g., 0.06 eV vs 0.0 eV) in simulations or on survey-like mock catalogs, and check whether the 1 Mpc suppression and 3–4 Mpc excess survive above jackknife errors. Alternatively, rerun the NEXUS+ classification with a different smoothing-scale hierarchy (e.g., R0=2 Mpc/h or n=0..5) and see whether the ~70% filament fraction and the environment-split edge-length signals remain unchanged; if they shift significantly, the claimed environment-MST link is not robust.

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

If this is right

  • MST edge lengths can be added to two-point clustering in galaxy surveys as an independent summary statistic, potentially tightening neutrino-mass constraints by a factor of several.
  • Because neutrinos and baryons move the edge-length distribution in opposite directions at small scales, the MST can help break a well-known degeneracy in the matter power spectrum.
  • The ~70% filament-edge fraction means the MST is effectively a filament tracer; its sensitivity to neutrino mass is a sensitivity to how neutrinos diffuse filamentary structure.
  • The environment-split edge-length signal gives a physical explanation for why earlier MST-based likelihoods improved neutrino constraints, and identifies which scales to weight in a future analysis.
  • The method is ready for survey-like number densities in real space; accounting for redshift-space distortions and survey geometry is the natural next step.

Where Pith is reading between the lines

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

  • The straight-line assignment of MST edges to environments treats an edge as a Euclidean bridge between subhaloes; real filament geometry can curve on scales comparable to the longest edges, so the reported scale–environment mapping may blur in more tortuous regions. A test would repeat the split with curved or density-weighted paths.
  • The simulations use neutrino masses up to 0.48 eV, far above the current upper bound. The signal is monotonic in mass, but the small- and large-scale features may drop below noise for realistic masses near 0.06 eV; a forecast of the required survey volume would tell whether the MST can actually compete with power-spectrum constraints at allowed masses.
  • The paper works in 3D real space with true subhalo positions. Survey data include redshift-space distortions and selection effects; a mock-observation pipeline that applies RSDs and a survey mask would test whether the environmental signatures survive, and could reveal that the MST's filament sensitivity is partially an artifact of geometric distortions.
  • The NEXUS+ classification itself has tunable smoothing scales and thresholds; if these are varied, the '70% in filaments' and the scale-environment linkage could shift. A robustness scan over smoothing hierarchy would establish whether the neutrino signature is tied to the classifier or to the underlying matter distribution.
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.

Referee Report

2 major / 5 minor

Summary. The paper uses FLAMINGO L1000 simulations with three neutrino masses (0.06, 0.24, and 0.48 eV), the NEXUS+ cosmic-web classifier, and minimum spanning trees (MSTs) constructed from 500,000 subhaloes selected to match a DESI-like number density. It reports that increasing neutrino mass reduces the volume occupied by clusters and voids, shifts and narrows the volume-weighted density distribution, and that roughly 70% of MST edges lie in filaments. The central quantitative claim is that MST edge lengths show a characteristic neutrino-mass signature—suppression of edges at l≈1 Mpc, enhancement at l≈3–4 Mpc, and suppression at l≈9 Mpc—while baryonic physics produces a distinct small-scale enhancement. The paper interprets these features as evidence that MST statistics are sensitive to neutrino mass and can complement two-point statistics for future surveys, and it attributes the scale-dependent signal to clusters, filaments, and voids respectively.

Significance. If the claimed edge-length signatures are robust, this is a useful step toward understanding why MST statistics improve neutrino-mass constraints relative to two-point statistics (Naidoo et al. 2022). The explicit decomposition of MST edges by cosmic-web environment, the use of a DESI-like subhalo number density, and the comparison of neutrino versus baryonic effects are all valuable and forward-model, non-circular: the neutrino mass is an input, and the MST/environment statistics are outputs. The authors are also careful to acknowledge that the simulated neutrino masses exceed current upper limits. However, the quantitative conclusions currently rest on a single simulation realization per cosmology, which is a serious limitation that needs to be addressed before the central claim can be accepted as established.

major comments (2)
  1. [Sections 2.1, 2.4, Fig. 5, Table 2] All quantitative claims—the ≈70% filament fraction, the scale-dependent edge-length differences (suppression at ~1 Mpc, enhancement at ~3–4 Mpc, deficit at ~9 Mpc), and the baryon-versus-neutrino separation—are drawn from exactly one FLAMINGO L1000 realization per neutrino mass. The jackknife errors from Eq. (8) partition this single box into 125 subboxes; these subboxes share long-wavelength density modes, so jackknife scatter severely underestimates cosmic variance. A single seed's particular initial conditions can produce fluctuations in the MST edge-length PDF that mimic or mask the ~1–3% differences shown in Fig. 5. The paper nowhere tests seed-to-seed reproducibility. Because the headline claim is that MST statistics are a reliable neutrino-mass probe, this is load-bearing: please provide multiple independent realizations for at least the 0.06 and 0.48 eV cases, or otherwise quanti
  2. [Section 2.2, Table 2, Fig. 4] The environment decomposition depends on the NEXUS+ smoothing scales (R0=1 Mpc/h, n=0..4) and the thresholds adopted from Cautun et al. (2013). The paper mentions that cosmic-web classifications are scheme-dependent (Libeskind et al. 2018) but does not test robustness to these choices. The interpretation that MST edges preferentially trace filaments, and the scale-dependent attributions (clusters at small scales, filaments at intermediate scales, voids at large scales), could shift if different smoothing scales or thresholds were used. Please add a sensitivity test—for example varying R0 or the thresholds, or comparing with an alternative classifier such as the T-web—to confirm that the ≈70% filament fraction and the environment–scale decomposition are not artefacts of the specific NEXUS+ settings adopted.
minor comments (5)
  1. [Section 3.2.1] Please define 'minimum resolution of the NEXUS+ classification' explicitly (presumably the grid cell size) and justify the oversampling choice for the edge-splitting procedure. The method assumes a straight line between subhaloes is representative of the environments the edge traverses; a brief note on this approximation would aid reproducibility.
  2. [Section 3.1.2] Typo: 'density biy' should read 'density bin'.
  3. [Figure 5 caption] 'different neutrinos masses' should be 'different neutrino masses'.
  4. [Section 3.2.2] The claim that baryons 'increase the number of edges at the smallest scales of l≈1 Mpc' is made from a single comparison (HYDRO_0.06 versus DMO_0.06). Given the jackknife errors, please state the significance of this difference explicitly rather than relying on visual inspection of the PDFs.
  5. [Discussion, final paragraph] The last paragraph says the results 'demonstrate' that MST statistics provide complementary information. Given the single-realization limitation and the unaccounted-for redshift-space distortions (noted in the abstract), 'suggest' or 'indicate' would be more appropriate.

Circularity Check

0 steps flagged

No significant circularity: the MST and cosmic-web results are measured outputs of a forward-model simulation, and the only self-citations are motivational rather than load-bearing.

full rationale

The paper is a simulation-forward analysis: neutrino mass is an input cosmology, and the NEXUS+ cosmic-web classifications, volume/mass fractions, and MST edge-length distributions are all measured outputs. No parameter is fitted to the target MST signal and then renamed a prediction. The NEXUS+ smoothing scales and thresholds are adopted from Cautun et al. (2013), an external public method, and the subhalo density is chosen to match DESI-like surveys (Section 2.3), not to optimize the neutrino signal. The central claim that MST edge lengths are sensitive to neutrino mass is supported directly by the measured differences among the 0.06, 0.24 and 0.48 eV FLAMINGO simulations (Fig. 5a), not only by the earlier Naidoo et al. (2022) result. That earlier work is used as motivation ('The origin of the MST's sensitivity to neutrino mass is unclear, but is believed to come from its ability to trace filaments'), and the paper independently quantifies environment fractions (Table 2) and scale-dependent attributions (Section 3.2.3). The explicit statement in Section 3.1.2 that the density hierarchy of environments is 'by construction' is an honest admission about that particular hierarchy, and it is not used as the basis for a predictive claim. The single-realization cosmic-variance issue raised by the skeptic is a statistical robustness concern, not circularity: jackknife resampling does not enter the definition of the measured MST signal. No equation or fitted input reduces to another by construction, and no load-bearing argument depends on a self-citation chain. Hence the circularity score is minimal.

Axiom & Free-Parameter Ledger

4 free parameters · 4 axioms · 0 invented entities

The central analysis depends on a small number of hand-chosen modeling parameters (subhalo mass cut, number density, NEXUS+ settings, jackknife partition) and on domain assumptions about the fidelity of FLAMINGO and the validity of NEXUS+ thresholds. No new physical entities are introduced.

free parameters (4)
  • Subhalo mass cut = 10^12 Msun
    Chosen to reduce MST sensitivity to density tracer and particle mass resolution; affects which subhaloes are used as MST nodes.
  • Number of MST subhaloes / number density = 500,000 (≈0.0016 h^3 Mpc^-3)
    Chosen to match DESI-like galaxy density; MST statistics can depend on tracer density.
  • NEXUS+ smoothing scales and thresholds = R0=1 Mpc/h, n=0..4; cluster volume <0.5 h^-3 Mpc^3, filament/wall <10 h^-3 Mpc^3
    Adopted from Cautun et al. (2013); classification results depend on these choices, which are not re-calibrated or tested for robustness.
  • Jackknife partition = 5^3 = 125 subboxes
    Choice of partition size for error estimation; affects error estimates.
axioms (4)
  • domain assumption FLAMINGO simulations with δf neutrino method and subgrid baryonic feedback accurately model the effect of neutrino mass on large-scale structure.
    The paper uses these simulations as ground truth to measure neutrino effects; no independent validation of the neutrino implementation is given here (see Section 2.1).
  • domain assumption NEXUS+ thresholds and smoothing scales from Cautun et al. (2013) are valid for the FLAMINGO density fields.
    Section 2.2 adopts the Cautun et al. thresholds without re-derivation; classification of environments depends on this.
  • domain assumption A random subset of 500,000 subhaloes (mass cut 10^12 Msun) is representative of the full subhalo population for MST topology.
    Section 3.2 selects subhaloes at random; the MST is a function of the point set, so undersampling could alter edge-length statistics.
  • standard math MST algorithms (Kruskal) and Euclidean edge lengths are appropriate statistics for cosmic web filaments.
    Standard graph theory; cited in Section 2.3.

reviewed 2026-08-03 · how reviews work

0 comments
Cite this review

Pith. "Pith review of How massive neutrinos reshape the cosmic web." pith.science (2026). https://pith.science/paper/W7YKOOMA

@misc{pith2026251216517,
  author       = {Pith},
  title        = {Pith review of: How massive neutrinos reshape the cosmic web},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/W7YKOOMA}},
  note         = {Machine review of arXiv:2512.16517}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

We explore the effects of massive neutrinos on the cosmic web using the FLAMINGO simulations. We classify the cosmic web into voids, sheets, filaments, and clusters, and find that massive neutrinos affect the environment by decreasing the volume occupied by clusters and voids. We find that increasing the neutrino mass shifts the volume-weighted density distribution towards higher densities and leads to a more narrow density distribution, which we interpret as neutrinos delaying structure formation. We construct the minimum spanning tree (MST) graph from the subhaloes, adopting a number density chosen to match that expected for DESI-like observations. We show that most MST edges lie in filaments, approximately 70% throughout different simulations, which we link to its sensitivity to neutrino mass. We also link the MST's edge length signal at different scales to different cosmic web environments, with clusters dominating the signal at small scales, voids at longer scales, and filaments at intermediate scales. The strong correlation between MST edges and cosmic web environments reinforces the MST's potential to be used as a classifier for large-scale structure in galaxy surveys. We compare the effects of baryonic physics and massive neutrinos and find that each produces distinct signatures in MST edge lengths. This analysis is performed in 3D space, using the true positions of subhaloes and not accounting for redshift space distortions. Nevertheless, these results emphasise the MST's capability to go beyond two-point statistics, motivating future applications to real observational data.

Figures

Figures reproduced from arXiv: 2512.16517 by Benjamin Joachimi, Carlos S. Frenk, Krishna Naidoo, Leonor N. L. Sim\~oes, Willem Elbers.

Figure 1
Figure 1. Figure 1: The same slice (projection of 155 × 155 × 5Mpc3 subvolume) of the total density contrast (top row, where 1 + 𝛿𝑡𝑜𝑡𝑎𝑙 = 𝜌𝑡𝑜𝑡𝑎𝑙 𝜌 ¯ 𝑡𝑜𝑡𝑎𝑙 ), the neutrino density contrast (centre row, where 1 + 𝛿𝜈 = 𝜌𝜈 𝜌¯𝜈 ), and the NEXUS+ classification (bottom row) for each simulation with different neutrino masses (each column is one mass). This highlights the subtlety of neutrino effects: although in the centre row we ca… view at source ↗
Figure 2
Figure 2. Figure 2: Volume (top) and mass (bottom) fractions, relative to the values for the DMO+𝜈 simulation with 𝑀𝜈 = 0.06 eV. sities when we increase neutrino mass, which we link again to the neutrinos washing out power on small scales. In the bottom panel of [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The normalised density (top) and mass-weighted density (middle) distributions for each neutrino mass. The top panels show the total density whereas in the bottom panels we subtract the 𝑀𝜈 = 0.06 eV distribution from the other two. In the bottom panel we split by cosmic web environment to show which components dominate at different densities. which we construct the MST. This corresponds to a density of ap￾p… view at source ↗
Figure 4
Figure 4. Figure 4: Histograms of the distribution of the edge lengths (𝑙) in the MST separated into each classification of the cosmic web that they are tracing. In Fig. 4a we plot the total probability distribution function (PDF), normalised to total counts. In Fig. 4b we plot the individual PDFs, normalised for each classification. These plots come from the analysis of the DMO+𝜈 simulation with 𝑀𝜈 = 0.06 eV , but the relati… view at source ↗
Figure 5
Figure 5. Figure 5: A comparison between the histogram distributions of the MST statistics in different models, using 500 000 haloes for the MST. In Fig. 5a we compare different neutrinos masses, and in Fig. 5b we compare the effect of baryons and the effect of massive neutrinos. In the top panels we show the distributions of the degrees (d), edge lengths (l), branch lengths (b), and branch shapes (s). The bottom panels show … view at source ↗
Figure 6
Figure 6. Figure 6: The edge length statistic (𝑙) split by each cosmic web environ￾ment, compared to the total statistic from [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗

discussion (0)

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

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Dark Energy in the DESI Era: A Brief Review of Evidence, Beyond-$\Lambda$CDM Interpretations, and Tensions

    astro-ph.CO 2026-06 conditional novelty 3.0

    DESI data hint that dark energy evolves (w0>−1, wa<0, phantom crossing) at 2.5–4.2σ, but the same expansion history can also be produced by interacting dark sectors, modified gravity, or non-cold dark matter.

  2. Dark Energy in the DESI Era: A Brief Review of Evidence, Beyond-$\Lambda$CDM Interpretations, and Tensions

    astro-ph.CO 2026-06 unverdicted novelty 1.0

    Review of DESI evidence for dynamical dark energy, its dependence on parametrization and datasets, and alternative beyond-LambdaCDM interpretations that may address cosmological tensions.

Reference graph

Works this paper leans on

1 extracted references · 1 linked inside Pith · cited by 1 Pith paper

  1. [1]

    Abbott T. M. C., et al., 2018, Phys. Rev. D, 98, 043526 Abdul Karim M., et al., 2025, Phys. Rev. D, 112, 083515 Agarwal S., Feldman H. A., 2011, MNRAS, 410, 1647 Alam S., et al., 2017, MNRAS, 470, 2617 Alpaslan M., et al., 2014, MNRAS, 438, 177 Aragón-CalvoM.A.,vandeWeygaertR.,JonesB.J.T.,2010,MNRAS,408, 2163 Bahé Y. M., et al., 2022, MNRAS, 516, 167 Barr...

    Pith/arXiv arXiv 2018

This paper was first reviewed by deepseek-v4-flash on August 3, 2026.