pith. machine review for the scientific record. sign in

arxiv: 2604.05033 · v1 · submitted 2026-04-06 · 🌌 astro-ph.CO · astro-ph.GA

Recognition: 2 theorem links

· Lean Theorem

Universal Dark-matter Density Profiles of Cosmic Filaments

Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:21 UTC · model grok-4.3

classification 🌌 astro-ph.CO astro-ph.GA
keywords cosmic filamentsdark matter density profilesuniversal profilescosmic webself-similarityfilament spinesterminal nodes
0
0 comments X

The pith

Scaling cosmic filament radii by the virial radii of their terminal nodes produces a nearly universal dark-matter density profile.

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

The analysis extracts filaments from the TNG50 simulation using DisPerSE and applies a shrinking-cylinder re-centering step to trace true density ridges more accurately. Once the radial coordinate is scaled by the virial radii of the nodes at each filament's ends, the density profiles collapse to a common shape with only weak dependence on redshift, node mass, or filament length. This indicates that filaments, like dark-matter halos, follow structural self-similarity when the right characteristic scale is chosen. Separating bound and unbound particles shows that the central cusp arises mainly from embedded low-mass halos while the smooth component develops a flat core inside one-tenth of the scaled radius, with the smooth part becoming less dominant at lower redshifts.

Core claim

After applying the shrinking-cylinder re-centering algorithm to correct the filament spines returned by DisPerSE, the radial dark-matter density profiles of cosmic filaments, when the radial coordinate is scaled by the virial radii of the terminal nodes, exhibit a nearly universal form that depends only weakly on redshift, node mass, and filament length. This result suggests that cosmic filaments obey a form of structural self-similarity once an appropriate characteristic scale is introduced. The apparent central cusp of the full profile is primarily produced by low-mass halos embedded along the filament spines, while the smooth component develops a flat core within R/R_vir ≲ 0.1. The redoxh

What carries the argument

The shrinking-cylinder re-centering algorithm that corrects DisPerSE spine locations to better follow density ridges, combined with scaling the radial coordinate by the virial radii of the filament's terminal nodes to reveal the universal profile shape.

If this is right

  • Filaments exhibit structural self-similarity comparable to that of dark-matter halos once the terminal-node virial radius is adopted as the scale.
  • The smooth unbound component develops a flat core at small scaled radii while embedded halos produce the central cusp.
  • The redshift evolution of the smooth component indicates a shift from predominantly smooth accretion at high redshift to increasingly clumpy accretion at late times.
  • The universal profile is accurately described by a generalized triple-power-law functional form.

Where Pith is reading between the lines

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

  • If the universality holds in observations, stacked weak-lensing or galaxy-count profiles around filaments could be used to infer average filament properties without detailed modeling of each structure.
  • The transition from smooth to clumpy accretion implied by the smooth-component evolution offers a testable prediction for how gas and galaxies are distributed along filaments at different epochs.
  • Adopting the same scaling in other large-scale structure analyses might reduce scatter in filament property measurements and improve comparisons between simulations and data.

Load-bearing premise

The shrinking-cylinder re-centering algorithm correctly recovers the true density ridges of the filaments and scaling by the virial radii of the terminal nodes supplies the appropriate characteristic length for revealing self-similarity.

What would settle it

Extracting filament density profiles from an independent high-resolution simulation or from galaxy survey data, applying the same re-centering and scaling procedure, and finding that the profiles fail to collapse onto a single curve across redshifts and masses would falsify the universality claim.

Figures

Figures reproduced from arXiv: 2604.05033 by David Koo, Douglas Hellinger, Fangzhou Jiang, Farhanul Hasan, Joanna Woo, Joel R. Primack, Peng Xu, Sandra M. Faber.

Figure 1
Figure 1. Figure 1: Illustration of the cosmic filament network identi￾fied by DisPerSE at z = 1 in TNG50-1. Blue lines show the filaments in a 20 cMpc-thick slice, while purple and green points mark maxima and saddle points, respectively. The background image shows the DM surface density projected over the same slice. Our analysis is restricted to the area enclosed by the white box, excluding the outer ∼ 5% of the simulation… view at source ↗
Figure 2
Figure 2. Figure 2: Left: Illustration of the method used to refine the spine location of a filament segment. The sequence of “shrinking cylinder” operations is indicated by the green, blue, and black cylinders. At the (k + 1)th iteration, the spine is repositioned to pass through the center of mass of particles enclosed within the cylindrical region of radius Rk from the kth iteration, after which the radius is reduced by 10… view at source ↗
Figure 3
Figure 3. Figure 3: Left: DM density profiles of filaments in physical radius at different redshifts. The unresolved region is shown in gray shades with the same definition as [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Dark-matter density profiles around filaments as functions of node mass (top) and filament length (bottom). In each panel, the length-weighted mean profiles across redshifts are shown for the mass bin or length bin indicated. The best-fit gNFW (dotted) and PL3 (dashed) models for the full sample, as in [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Impact of halos near the filament spine on the DM density profiles. Left: Comparison of length-weighted mean filament density profiles calculated using all DM particles (dashed) and unbound DM particles (solid). Here the “unbound DM” profile is constructed by removing all DM particles bound to central halos with mass above 108M⊙ together with the volumes occupied by these halos. The resulting smooth filame… view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of length-weighted mean filament profiles from this work with those reported by W. Wang et al. (2024), which used galaxies in the MTNG simulation to identify and characterize filaments. The upper and lower rows present the density profiles and the logarithmic density slopes, respectively. The left column shows the W. Wang et al. results. The middle column shows our results measured around the ra… view at source ↗
read the original abstract

We present a comprehensive analysis of the radial dark-matter (DM) density profiles of cosmic filaments in the hydrodynamical simulation TNG50. The cosmic web is extracted from high-resolution density grids at redshifts $z =$ 0, 0.5, 1, 2 and 3 using the DisPerSE algorithm. We show that the filament spine locations returned directly by DisPerSE do not accurately reflect the true density ridges. To address this issue, we introduce a "shrinking-cylinder" re-centering algorithm, which significantly increases the inferred central densities and restores the inner power-law behavior of the profiles. When the radial coordinate is scaled by the virial radii of the terminal nodes, the filament density profiles exhibit a nearly universal form that depends only weakly on redshift, node mass, and filament length. This result suggests that cosmic filaments, much like dark-matter halos, obey a form of structural self-similarity once an appropriate characteristic scale is introduced. By repeating the measurement using only smoothly distributed, unbound DM particles, we find that the apparent central cusp of the full profile is primarily produced by low-mass halos embedded along the filament spines, while the smooth component develops a flat core within $R/R_{\rm vir}\lesssim0.1$. The redshift evolution of this smooth component further suggests a transition from predominantly smooth filamentary accretion at high redshift to increasingly clumpy accretion at late times. Finally, we show that the universal filament profile is accurately described by a generalized triple-power-law model.

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

3 major / 3 minor

Summary. The paper analyzes radial dark-matter density profiles of cosmic filaments extracted via DisPerSE from TNG50 hydrodynamical simulations at z=0,0.5,1,2,3. It introduces a shrinking-cylinder re-centering algorithm to correct spine locations, reports that scaling radii by terminal-node virial radii yields a nearly universal profile with weak dependence on redshift, node mass and filament length, fits the profiles with a generalized triple-power-law, and separates the smooth unbound component (showing a core) from the clumpy halo contribution (producing the cusp), with implications for accretion history.

Significance. If validated, the result would establish structural self-similarity for filaments analogous to NFW halos, offering a new characteristic scale and a quantitative description of the transition from smooth to clumpy accretion. Strengths include the use of high-resolution TNG50 data across multiple redshifts with explicit methodology and the decomposition into smooth versus embedded-halo components, which provides falsifiable predictions for future observations or simulations.

major comments (3)
  1. [§3.2] §3.2 (shrinking-cylinder algorithm description): the re-centering procedure is shown to raise central densities and restore an inner power-law, but no quantitative validation on synthetic filaments with a priori known ridge locations and radial profiles is presented. Recovery accuracy metrics (e.g., positional offset distributions or profile reconstruction error) on controlled mocks are required to demonstrate that the restored cusp is not an algorithmic artifact; this step is load-bearing for the universality claim.
  2. [§4.3] §4.3 and Figure 7 (universality tests): scaling by terminal-node R_vir is asserted to reveal self-similarity with only weak residual dependence on redshift, mass and length, yet no systematic comparison to alternative characteristic scales (filament length, mean inter-node distance, or local density) is shown. Without such tests it remains possible that the apparent universality is specific to the chosen scale rather than intrinsic.
  3. [§5.2] §5.2 (smooth-component analysis): the claim that the central cusp arises primarily from low-mass embedded halos while the smooth component develops a flat core at R/R_vir ≲ 0.1 rests on the re-centering step; any residual bias in ridge recovery would propagate directly into the reported redshift evolution of the smooth profile and the inferred transition from smooth to clumpy accretion.
minor comments (3)
  1. The functional form and free parameters of the generalized triple-power-law model should be written explicitly (with equation number) in the text where it is first introduced, rather than only in the abstract.
  2. Figure captions for the profile stacks (e.g., Figures 4–6) should state the number of filaments in each stack and the bootstrap or jackknife method used for the shaded uncertainty regions.
  3. A short paragraph comparing the derived filament profile to existing literature parametrizations (e.g., those based on N-body or other hydro runs) would help place the triple-power-law result in context.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment point by point below, indicating where the manuscript will be revised to incorporate the suggestions.

read point-by-point responses
  1. Referee: [§3.2] §3.2 (shrinking-cylinder algorithm description): the re-centering procedure is shown to raise central densities and restore an inner power-law, but no quantitative validation on synthetic filaments with a priori known ridge locations and radial profiles is presented. Recovery accuracy metrics (e.g., positional offset distributions or profile reconstruction error) on controlled mocks are required to demonstrate that the restored cusp is not an algorithmic artifact; this step is load-bearing for the universality claim.

    Authors: We agree that quantitative validation on synthetic filaments would strengthen the case that the restored inner power-law is not an artifact of the shrinking-cylinder procedure. The current manuscript shows the effect by direct before-and-after comparison on the TNG50 data, but does not include controlled mocks. In the revised manuscript we will add an appendix with synthetic tests: we will generate mock filaments with prescribed ridge lines and radial density profiles, apply the algorithm, and report recovery metrics including positional offset distributions and profile reconstruction errors. This will directly address the concern. revision: yes

  2. Referee: [§4.3] §4.3 and Figure 7 (universality tests): scaling by terminal-node R_vir is asserted to reveal self-similarity with only weak residual dependence on redshift, mass and length, yet no systematic comparison to alternative characteristic scales (filament length, mean inter-node distance, or local density) is shown. Without such tests it remains possible that the apparent universality is specific to the chosen scale rather than intrinsic.

    Authors: The scaling by terminal-node virial radii is chosen because of the physical role of the nodes in setting the filament potential; the manuscript reports only weak residual trends with redshift, mass and length under this scaling. We acknowledge that alternative normalizations were not systematically tested. In the revised version we will add a comparison (new figure or panels in Figure 7) showing the scatter in the stacked profiles when radii are instead scaled by filament length, mean inter-node distance, and local density. This will allow readers to evaluate whether the terminal R_vir scaling indeed minimizes dispersion. revision: yes

  3. Referee: [§5.2] §5.2 (smooth-component analysis): the claim that the central cusp arises primarily from low-mass embedded halos while the smooth component develops a flat core at R/R_vir ≲ 0.1 rests on the re-centering step; any residual bias in ridge recovery would propagate directly into the reported redshift evolution of the smooth profile and the inferred transition from smooth to clumpy accretion.

    Authors: We recognize that the smooth-component results depend on the accuracy of the spine locations after re-centering. The decomposition itself removes bound halo particles, but the radial coordinate uses the corrected ridge. In the revised §5.2 we will add explicit discussion of possible residual biases, include a direct comparison of smooth profiles computed before and after re-centering, and note that any uniform bias would affect the absolute normalization but not necessarily the reported redshift trend in the core development. The transition from smooth to clumpy accretion is inferred from the increasing relative contribution of the clumpy component at lower redshift, which is robust to small radial shifts. revision: partial

Circularity Check

0 steps flagged

No significant circularity; universality is an empirical observation from simulation data after re-centering and scaling

full rationale

The paper extracts cosmic filaments from TNG50 density grids using DisPerSE, introduces a shrinking-cylinder re-centering procedure to adjust spine locations, scales the radial coordinate by the virial radii of terminal nodes, and reports that the resulting DM density profiles collapse to a nearly universal form with only weak dependence on redshift, node mass, and filament length. This is presented as a direct measurement from simulation snapshots rather than a derivation from equations. The generalized triple-power-law model is introduced as a descriptive fit to the observed profiles. No steps reduce to self-definition, fitted inputs renamed as predictions, or load-bearing self-citations/uniqueness theorems. The chain is self-contained empirical analysis, consistent with the default expectation of no significant circularity.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the fidelity of the TNG50 hydrodynamical simulation and on the validity of the newly introduced shrinking-cylinder re-centering procedure; the descriptive triple-power-law model introduces fitted parameters but the universality itself is an observed scaling relation.

free parameters (1)
  • triple-power-law slopes and break radii
    The generalized triple-power-law model used to describe the universal profile contains free parameters that are fitted to the measured densities.
axioms (2)
  • domain assumption TNG50 accurately reproduces the dark-matter distribution and halo properties in a Lambda-CDM universe.
    All density measurements and the separation into smooth versus bound particles rely on the simulation's output.
  • domain assumption DisPerSE plus shrinking-cylinder correctly identifies filament spines.
    The paper introduces the re-centering step to correct DisPerSE, making this assumption load-bearing for the central densities.

pith-pipeline@v0.9.0 · 5595 in / 1654 out tokens · 60620 ms · 2026-05-10T19:21:56.497467+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

33 extracted references · 29 canonical work pages · 2 internal anchors

  1. [1]

    2006, MNRAS, 373, 1619, doi: 10.1111/j.1365-2966.2006.11113.x 3I/ATLAS Observed by Tianwen-1 Orbiter17 Salazar Manzano, L

    Allgood, B., Flores, R. A., Primack, J. R., et al. 2006, Monthly Notices of the Royal Astronomical Society, 367, 1781, doi: 10.1111/j.1365-2966.2006.10094.x

  2. [2]

    M., & Jablonka, P

    Bahe, Y. M., & Jablonka, P. 2025, Galaxies in the Simulated Cosmic Web: I. Filament Identification and Their Properties, arXiv, doi: 10.48550/arXiv.2502.06484 11 Figure A1.Left: Mass function of the FoF groups associated with nodes at different redshifts.Right: Filament length distribution at different redshifts. By definition, filaments are sets of conti...

  3. [3]

    Richard and Kofman, Lev and Pogosyan, Dmitri

    Bond, J. R., Kofman, L., & Pogosyan, D. 1996, Nature, 380, 603, doi: 10.1038/380603a0

  4. [4]

    Cautun, M., Van De Weygaert, R., & Jones, B. J. T. 2013, Monthly Notices of the Royal Astronomical Society, 429, 1286, doi: 10.1093/mnras/sts416

  5. [5]

    G., Baugh, C

    Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, Monthly Notices of the Royal Astronomical Society

  6. [6]

    2009, MNRAS, 398, 515, doi: 10.1111/j.1365-2966.2009.14895.x

    Dolag, K., Borgani, S., Murante, G., & Springel, V. 2009, Monthly Notices of the Royal Astronomical Society, 399, 497, doi: 10.1111/j.1365-2966.2009.15034.x

  7. [7]

    A., & Macci` o, A

    Dutton, A. A., & Macci` o, A. V. 2014, Monthly Notices of the Royal Astronomical Society, 441, 3359, doi: 10.1093/mnras/stu742 Gal´ arraga-Espinosa, D., Aghanim, N., Langer, M., Gouin, C., & Malavasi, N. 2020, Astronomy & Astrophysics, 641, A173, doi: 10.1051/0004-6361/202037986 Gal´ arraga-Espinosa, D., Cadiou, C., Gouin, C., et al. 2024, Astronomy & Ast...

  8. [8]

    2016, Monthly Notices of the Royal Astronomical Society, 462, 448, doi: 10.1093/mnras/stw1595

    Gheller, C., Vazza, F., Br¨ uggen, M., et al. 2016, Monthly Notices of the Royal Astronomical Society, 462, 448, doi: 10.1093/mnras/stw1595

  9. [9]

    N., Hellinger, D., et al

    Hasan, F., Burchett, J. N., Hellinger, D., et al. 2024, The Astrophysical Journal, 970, 177, doi: 10.3847/1538-4357/ad4ee2 Hern´ andez-Aguayo, C., Springel, V., Pakmor, R., et al. 2023, Monthly Notices of the Royal Astronomical Society, 524, 2556, doi: 10.1093/mnras/stad1657

  10. [10]

    Jiang, F., & Van Den Bosch, F. C. 2014, Monthly Notices of the Royal Astronomical Society, 440, 193, doi: 10.1093/mnras/stu280

  11. [11]

    2008, Riemannian Geometry and Geometric

    Jost, J., ed. 2008, Riemannian Geometry and Geometric

  12. [12]

    Analysis, Universitext (Berlin, Heidelberg: Springer-Verlag Berlin Heidelberg), doi: 10.1007/978-3-540-77341-2

  13. [13]

    , keywords =

    Kannan, R., Springel, V., Hernquist, L., et al. 2023, Monthly Notices of the Royal Astronomical Society, 524, 2594, doi: 10.1093/mnras/stac3743

  14. [14]

    Knebe, A., Gill, S. P. D., Gibson, B. K., et al. 2004, The Astrophysical Journal, 603, 7, doi: 10.1086/381306

  15. [15]

    I., et al., 2018, MNRAS, 473, 1195, doi:10.1093/mnras/stx1976 https://doi.org/10.1093/mnras/stx1976, arXiv:1705.03021 https://arxiv.org/abs/1705.03021

    Libeskind, N. I., Van De Weygaert, R., Cautun, M., et al. 2018, Monthly Notices of the Royal Astronomical Society, 473, 1195, doi: 10.1093/mnras/stx1976

  16. [16]

    Milnor, J. W. 2016, Morse Theory: AM-51, Annals of Mathematics Studies No. 51 (Princeton, NJ: Princeton University Press), doi: 10.1515/9781400881802

  17. [17]

    F., Frenk, C

    Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, The Astrophysical Journal

  18. [18]

    First results from the IllustrisTNG simulations: the galaxy color bimodality

    Nelson, D., Pillepich, A., Springel, V., et al. 2018, Monthly Notices of the Royal Astronomical Society, 475, 624, doi: 10.1093/mnras/stx3040

  19. [19]

    2019, Computational Astrophysics and Cosmology, 6, 2, doi: 10.1186/s40668-019-0028-x

    Nelson, D., Springel, V., Pillepich, A., et al. 2019, Computational Astrophysics and Cosmology, 6, 2, doi: 10.1186/s40668-019-0028-x

  20. [20]

    doi:10.1093/mnras/stac3620 , archiveprefix =

    Pakmor, R., Springel, V., Coles, J. P., et al. 2023, Monthly Notices of the Royal Astronomical Society, 524, 2539, doi: 10.1093/mnras/stac3620

  21. [21]

    2007, MNRAS, 382, 1833, doi: 10.1111/j.1365-2966.2007.12492.x

    Parkinson, H., Cole, S., & Helly, J. 2008, Monthly Notices of the Royal Astronomical Society, 383, 557, doi: 10.1111/j.1365-2966.2007.12517.x

  22. [22]

    2018b, MNRAS, 473, 4077, doi: 10.1093/mnras/stx2656

    Pillepich, A., Springel, V., Nelson, D., et al. 2018, Monthly Notices of the Royal Astronomical Society, 473, 4077, doi: 10.1093/mnras/stx2656 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, Astronomy & Astrophysics, 594, A13, doi: 10.1051/0004-6361/201525830

  23. [23]

    S., Shaviv, N

    Power, C., Navarro, J. F., Jenkins, A., et al. 2003, Monthly Notices of the Royal Astronomical Society, 338, 14, doi: 10.1046/j.1365-8711.2003.05925.x

  24. [24]

    2025, The Astrophysical Journal Supplement Series, 280, 53, doi: 10.3847/1538-4365/adfa26

    Shi, F., Wang, Z., Yang, X., et al. 2025, The Astrophysical Journal Supplement Series, 280, 53, doi: 10.3847/1538-4365/adfa26

  25. [25]

    W., & Ebeling, H

    Sousbie, T. 2011, Monthly Notices of the Royal Astronomical Society, 414, 350, doi: 10.1111/j.1365-2966.2011.18394.x

  26. [26]

    W., & Ebeling, H

    Sousbie, T., Pichon, C., & Kawahara, H. 2011, Monthly Notices of the Royal Astronomical Society, Volume 414, Issue 1, pp. 384-403., 414, 384, doi: 10.1111/j.1365-2966.2011.18395.x

  27. [27]

    2009, MNRAS, 398, 515, doi: 10.1111/j.1365-2966.2009.14895.x

    Springel, V. 2010, Monthly Notices of the Royal Astronomical Society, 401, 791, doi: 10.1111/j.1365-2966.2009.15715.x

  28. [28]

    Springel, V., White, S. D. M., Tormen, G., & Kauffmann, G. 2001, Monthly Notices of the Royal Astronomical Society, 328, 726, doi: 10.1046/j.1365-8711.2001.04912.x van Haarlem, M., & van de Weygaert, R. 1993, The Astrophysical Journal, 418, 544, doi: 10.1086/173416

  29. [29]

    2024, Monthly Notices of the Royal Astronomical Society, 532, 4604, doi: 10.1093/mnras/stae1801

    Wang, W., Wang, P., Guo, H., et al. 2024, Monthly Notices of the Royal Astronomical Society, 532, 4604, doi: 10.1093/mnras/stae1801

  30. [30]

    2023, Science China

    Wang, Z., Shi, F., Yang, X., et al. 2023, Science China

  31. [31]

    Physics, Mechanics & Astronomy, 67, 219513, doi: 10.1007/s11433-023-2192-9

  32. [32]

    C., et al

    Yang, S., Du, X., Zeng, Z. C., et al. 2023, The Astrophysical Journal, 946, 47, doi: 10.3847/1538-4357/acbd49

  33. [33]

    1996, MNRAS, 278, 488, doi: 10.1093/mnras/278.2.488

    Zhao, H. 1996, Monthly Notices of the Royal Astronomical Society, 278, 488, doi: 10.1093/mnras/278.2.488