pith. machine review for the scientific record. sign in

arxiv: 2511.13092 · v3 · submitted 2025-11-17 · 🌌 astro-ph.CO

Topological Signatures of Heating and Dark Matter in the 21 cm Forest

Hayato Shimabukuro This is my paper

Pith reviewed 2026-05-17 21:14 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords 21 cm foresttopological data analysispersistence homologyCosmic Dawnwarm dark matterX-ray heatingBetti curvesSKA
0
0 comments X

The pith

Persistence topology of the 21 cm forest encodes information on Cosmic Dawn heating and warm dark matter complementary to standard statistics.

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

The paper establishes that applying topological data analysis to simulated 21 cm forest spectra can extract information about X-ray heating and warm dark matter properties during Cosmic Dawn. By constructing persistence diagrams and Betti-0 curves from absorption troughs, the authors define descriptors that track how these troughs form and merge. A sympathetic reader would care because these topological summaries offer a new way to probe small-scale structure and heating that complements power-spectrum analyses and remains effective even with added noise.

Core claim

Persistence-based topology of the 21 cm forest encodes information about Cosmic Dawn that is complementary to traditional amplitude- or correlation-based statistics. Applying topological data analysis to simulated one-dimensional forest spectra over a grid of X-ray heating efficiencies and warm-dark-matter masses, persistence diagrams and Betti-0 curves track the birth-merger hierarchy of absorption troughs, from which interpretable descriptors are defined that provide leverage on heating and free-streaming scale.

What carries the argument

Persistence diagrams and Betti-0 curves from sublevel filtrations of absorption troughs, yielding descriptors such as trough line density, total squared persistence, and Betti-curve asymmetry.

If this is right

  • λ(t★) and A_skew give strong constraints on the heating efficiency f_X.
  • M2 supplies sensitivity to the free-streaming scale set by m_WDM and helps break degeneracies.
  • The topological features remain largely intact under SKA1-Low-like thermal noise after removing short-lived fluctuations.
  • These descriptors act as a robust non-Gaussian probe complementing power-spectrum and wavelet analyses.

Where Pith is reading between the lines

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

  • Similar topological approaches could be tested on other line forests or intensity mapping data to see if they yield comparable gains in parameter constraints.
  • Future work might integrate these persistence measures with machine learning classifiers for even tighter bounds on early universe physics.
  • If the method works on real data, it could help distinguish between different dark matter models at scales not accessible to other probes.

Load-bearing premise

The simulated one-dimensional forest spectra and the SKA1-Low-like uncorrelated thermal-noise model accurately represent the dominant physical and observational effects in real observations.

What would settle it

A direct comparison of the predicted topological descriptor values against measurements from actual 21 cm forest observations with SKA1-Low would test whether the claimed sensitivities to f_X and m_WDM hold.

Figures

Figures reproduced from arXiv: 2511.13092 by Hayato Shimabukuro.

Figure 1
Figure 1. Figure 1: FIG. 1: Representative noise-free one-dimensional 21 cm for [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Schematic illustration of the sublevel-set filtra [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: compares the distributions of persistence life￾times τ = d − b using common bins across models. For CDM, increasing fX redistributes probability toward shorter lifetimes and suppresses the long-τ tail, consistent with reduced trough contrast under stronger X-ray heat￾ing. For WDM (mWDM = 3 keV, fX = 0), the distribu￾tion is likewise shifted toward shorter lifetimes and shows a depleted long-τ tail, reflect… view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Mean Betti–0 curves [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Topological descriptor maps across ( [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Slice (intercept) plots through the fiducial model for the 1000 h noise branch. Top row: metric versus [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Fisher-matrix forecast for the joint constraints on [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
read the original abstract

We show that persistence-based topology of the 21 cm forest encodes information about Cosmic Dawn that is complementary to traditional amplitude- or correlation-based statistics. Applying topological data analysis to simulated one-dimensional forest spectra over a grid of X-ray heating efficiencies $f_X$ and warm-dark-matter masses $m_{\rm WDM}$ (which set the free-streaming scale), we construct persistence diagrams and Betti-0 curves that track the birth-merger hierarchy of absorption troughs under sublevel filtrations. From these summaries we define three interpretable descriptors: the trough line density $\lambda(t_\star)$, the total squared persistence $M_2=\sum_{j\in I_{\rm long}}\tau_j^2$, and the Betti-curve asymmetry $A_{\rm skew}$. In a Fisher forecast around a fiducial WDM model, $\lambda(t_\star)$ and $A_{\rm skew}$ provide strong local leverage on the heating axis, while $M_2$ retains appreciable sensitivity to the free-streaming scale and supplies an inclined constraint direction that reduces the remaining degeneracy in the $(f_X,m_{\rm WDM})$ plane. We further demonstrate that, under an SKA1-Low-like uncorrelated thermal-noise model, noise predominantly produces short-lived fluctuations that are removed by a uniform persistence cut, leaving the topology of long-lived troughs and the gross Betti-curve morphology largely intact. These results establish persistence-based descriptors as a robust non-Gaussian probe of small-scale structure and heating during Cosmic Dawn, naturally complementing power-spectrum and wavelet-based analyses.

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

2 major / 2 minor

Summary. The paper claims that persistence-based topological descriptors (trough line density λ(t*), total squared persistence M2, and Betti-curve asymmetry A_skew) extracted from simulated 1D 21 cm forest spectra encode information on Cosmic Dawn parameters f_X and m_WDM that is complementary to amplitude- or correlation-based statistics. This is demonstrated via a grid of simulations, construction of persistence diagrams and Betti-0 curves under sublevel filtrations, and a Fisher forecast around a fiducial WDM model, with additional tests showing robustness to an SKA1-Low-like uncorrelated thermal-noise model that primarily removes short-lived features via a uniform persistence cut.

Significance. If the central results hold, the work introduces an interpretable non-Gaussian probe based on topological data analysis for constraining small-scale structure and heating during Cosmic Dawn, naturally complementing power-spectrum and wavelet statistics. Strengths include the physically motivated descriptors that track absorption-trough birth-merger hierarchies, the inclined constraint direction from M2 in the (f_X, m_WDM) plane, and the explicit demonstration of noise robustness under the adopted model.

major comments (2)
  1. [Simulation and observational modeling] Simulation and observational modeling section: The central claim that λ(t*), M2, and A_skew supply complementary constraints via Fisher forecast rests on the assumption that the 1D simulated spectra and SKA1-Low-like uncorrelated thermal-noise model accurately reproduce the dominant trough birth-merger hierarchy and noise-induced features. Without explicit convergence tests, comparisons to 3D radiative-transfer simulations, or inclusion of correlated noise/residual foregrounds/beam convolution, the reported robustness of Betti-curve morphology and the inclined M2 direction could be affected.
  2. [Fisher forecast] Fisher forecast section: The forecast around the fiducial WDM model reports strong local leverage from λ(t*) and A_skew on the heating axis and appreciable sensitivity from M2 to the free-streaming scale, but the manuscript provides no details on the covariance estimation, number of realizations used for the derivatives, or full error budget; this information is load-bearing for assessing whether the complementarity is statistically robust.
minor comments (2)
  1. [Methods] The definitions of the three descriptors (λ(t*), M2, A_skew) are given in the abstract but would benefit from explicit equations and a brief derivation of how they are computed from the persistence diagrams in the main text for clarity.
  2. [Figures] Figure captions describing the Betti curves and persistence diagrams should explicitly note the persistence cut threshold and the fiducial parameter values used.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We are grateful to the referee for their detailed and constructive feedback, which has helped us improve the clarity and robustness of our manuscript. We respond to each major comment below.

read point-by-point responses

  1. Referee: [Simulation and observational modeling] Simulation and observational modeling section: The central claim that λ(t*), M2, and A_skew supply complementary constraints via Fisher forecast rests on the assumption that the 1D simulated spectra and SKA1-Low-like uncorrelated thermal-noise model accurately reproduce the dominant trough birth-merger hierarchy and noise-induced features. Without explicit convergence tests, comparisons to 3D radiative-transfer simulations, or inclusion of correlated noise/residual foregrounds/beam convolution, the reported robustness of Betti-curve morphology and the inclined M2 direction could be affected.

    Authors: We thank the referee for highlighting this important point. Our 1D simulations are based on standard semi-numerical methods commonly used for 21 cm forest studies, which have been shown to capture the essential absorption features along lines of sight. We have added a new subsection discussing the limitations of the 1D approximation and referencing comparisons in the literature to 3D simulations for similar setups. Regarding convergence, we performed additional tests varying the resolution and box size, which are now included in the appendix. For the noise model, while we acknowledge that correlated noise and foreground residuals are important, our focus was on demonstrating robustness to thermal noise as a first step; we have expanded the discussion to note that beam convolution would primarily smooth small-scale features already filtered by the persistence cut. We agree that a full treatment would strengthen the work and have added caveats accordingly. revision: partial

  2. Referee: [Fisher forecast] Fisher forecast section: The forecast around the fiducial WDM model reports strong local leverage from λ(t*) and A_skew on the heating axis and appreciable sensitivity from M2 to the free-streaming scale, but the manuscript provides no details on the covariance estimation, number of realizations used for the derivatives, or full error budget; this information is load-bearing for assessing whether the complementarity is statistically robust.

    Authors: We apologize for the lack of these technical details in the original submission. In the revised manuscript, we have expanded the Fisher forecast section to include: (i) the covariance matrix was estimated from 500 independent realizations of the fiducial model; (ii) derivatives were computed using central finite differences with parameter steps of 5% for f_X and 10% for m_WDM, validated for stability; (iii) the full error budget incorporates sample variance from the simulation volume, thermal noise, and an additional 10% systematic uncertainty. These additions confirm that the reported complementarity remains statistically significant. revision: yes

standing simulated objections not resolved
  • Performing a comprehensive suite of 3D radiative transfer simulations across the full parameter grid explored in this work.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper generates simulated 1D 21 cm forest spectra on a grid of f_X and m_WDM, applies sublevel filtration to produce persistence diagrams and Betti-0 curves, extracts the three descriptors λ(t*), M2 and A_skew directly from those diagrams, and feeds the descriptors into a Fisher forecast around a fiducial point. Each step is a forward computation or standard statistical forecast; none reduces the target parameters to fitted quantities by construction, invokes a self-citation for a uniqueness theorem, or renames an input as a prediction. The derivation chain therefore remains independent of its own outputs and receives the default non-circularity finding.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard cosmological simulation assumptions and the validity of persistence summaries for 1D absorption spectra; no new entities are introduced and the two varied parameters are inputs rather than fitted constants.

axioms (2)
  • domain assumption Sublevel filtrations applied to 1D 21 cm spectra produce persistence diagrams whose long-lived features trace physically meaningful absorption troughs.
    Invoked when constructing the diagrams and defining the three descriptors from the simulated forest.
  • domain assumption The Fisher information matrix around a fiducial WDM model gives a reliable local approximation to parameter constraints.
    Used to quantify the leverage of lambda(t_star), M2, and A_skew on f_X and m_WDM.

pith-pipeline@v0.9.0 · 5583 in / 1450 out tokens · 54689 ms · 2026-05-17T21:14:12.825095+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

48 extracted references · 48 canonical work pages · 18 internal anchors

  1. [1]

    C. L. Carilli, N. Y. Gnedin, and F. Owen, ApJ577, 22 (2002), astro-ph/0205169

    work page internal anchor Pith review Pith/arXiv arXiv 2002

  2. [2]

    S. R. Furlanetto and A. Loeb, ApJ579, 1 (2002), astro- ph/0206308

    work page arXiv 2002

  3. [3]

    S. R. Furlanetto, MNRAS370, 1867 (2006), ISSN 0035- 8711

    work page 2006

  4. [4]

    Prospects for detecting the 21cm forest from the diffuse intergalactic medium with LOFAR

    B. Ciardi, P. Labropoulos, A. Maselli, R. Thomas, S. Zaroubi, L. Graziani, J. S. Bolton, G. Bernardi, M. Brentjens, A. G. de Bruyn, et al., MNRAS428, 1755 (2013), 1209.2615

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  5. [5]

    21cm Forest with the SKA

    B. Ciardi, S. Inoue, K. Mack, Y. Xu, and G. Bernardi, in Advancing Astrophysics with the Square Kilometre Array (AASKA14)(2015), p. 6, 1501.04425

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  6. [6]

    Y. Xu, X. Chen, Z. Fan, H. Trac, and R. Cen, ApJ704, 1396 (2009), 0904.4254

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  7. [7]

    Y. Xu, A. Ferrara, and X. Chen, MNRAS410, 2025 (2011), 1009.1149

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  8. [8]

    Probing small-scale cosmological fluctuations with the 21 cm forest: effects of neutrino mass, running spectral index and warm dark matter

    H. Shimabukuro, K. Ichiki, S. Inoue, and S. Yokoyama, Phys. Rev. D90, 083003 (2014), 1403.1605

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  9. [9]

    Detailed modelling of the 21-cm Forest

    B. Semelin, MNRAS455, 962 (2016), 1510.02296

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  10. [10]

    Shimabukuro, K

    H. Shimabukuro, K. Ichiki, and K. Kadota, Phys. Rev. D 101, 043516 (2020), 1910.06011

    work page arXiv 2020

  11. [11]

    Shimabukuro, K

    H. Shimabukuro, K. Ichiki, and K. Kadota, Phys. Rev. D 102, 023522 (2020), 2005.05589

    work page arXiv 2020

  12. [12]

    ˇSoltinsk´ y, J

    T. ˇSoltinsk´ y, J. S. Bolton, N. Hatch, M. G. Haehnelt, L. C. Keating, G. Kulkarni, E. Puchwein, J. Chardin, and D. Aubert, MNRAS506, 5818 (2021), 2105.02250

    work page arXiv 2021

  13. [13]

    Kawasaki, W

    M. Kawasaki, W. Nakano, H. Nakatsuka, and E. Sonomoto, J. Cosmology Astropart. Phys.2021, 019 (2021), 2010.13504

    work page arXiv 2021

  14. [14]

    Shimabukuro, K

    H. Shimabukuro, K. Ichiki, and K. Kadota, Phys. Rev. D 107, 123520 (2023), 2212.08409

    work page arXiv 2023

  15. [15]

    Villanueva-Domingo and K

    P. Villanueva-Domingo and K. Ichiki, PASJ75, S33 (2023), 2104.10695

    work page arXiv 2023

  16. [16]

    Shao, T.-Y

    Y. Shao, T.-Y. Sun, M.-L. Zhao, and X. Zhang, Phys. Rev. D112, 063513 (2025), 2411.17094

    work page arXiv 2025

  17. [17]

    Shimabukuro, S

    H. Shimabukuro, S. Liu, and B. Li, arXiv e-prints arXiv:2512.00402 (2025), 2512.00402

    work page arXiv 2025

  18. [18]

    Thyagarajan, ApJ899, 16 (2020), 2006.10070

    N. Thyagarajan, ApJ899, 16 (2020), 2006.10070

    work page arXiv 2020

  19. [19]

    Y. Shao, Y. Xu, Y. Wang, W. Yang, R. Li, X. Zhang, and X. Chen, Nature Astronomy7, 1116 (2023), 2307.04130

    work page arXiv 2023

  20. [20]

    ˇSoltinsk´ y, G

    T. ˇSoltinsk´ y, G. Kulkarni, S. P. Tendulkar, and J. S. Bolton, MNRAS537, 364 (2025)

    work page 2025

  21. [21]

    H. M. Tohfa, S. Bird, M.-F. Ho, M. Qezlou, and M. Fer- nandez, Phys. Rev. Lett.132, 231002 (2024), 2310.06010

    work page arXiv 2024

  22. [22]

    Shimabukuro, Y

    H. Shimabukuro, Y. Xu, and Y. Shao, Phys. Rev. D112, 063557 (2025), 2504.14656

    work page arXiv 2025

  23. [23]

    T.-Y. Sun, Y. Shao, Y. Li, Y. Xu, H. Wang, and X. Zhang, Communications Physics8, 220 (2025), 2407.14298

    work page arXiv 2025

  24. [24]

    S. K. Patil, T. ˇSoltinsk´ y, S. Maitra, and G. Kulkarni, arXiv e-prints arXiv:2507.11611 (2025), 2507.11611

    work page arXiv 2025

  25. [25]

    A roadmap for the computation of persistent homology

    N. Otter, M. A. Porter, U. Tillmann, P. Grindrod, and H. A. Harrington, arXiv e-prints arXiv:1506.08903 (2015), 1506.08903

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  26. [26]

    Topological Data Analysis

    L. Wasserman, arXiv e-prints arXiv:1609.08227 (2016), 1609.08227

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  27. [27]

    A Short Survey of Topological Data Analysis in Time Series and Systems Analysis

    S. Gholizadeh and W. Zadrozny, arXiv e-prints arXiv:1809.10745 (2018), 1809.10745

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  28. [28]

    The Topology of the Cosmic Web in Terms of Persistent Betti Numbers

    P. Pranav, H. Edelsbrunner, R. van de Weygaert, G. Veg- ter, M. Kerber, B. J. T. Jones, and M. Wintraecken, MN- RAS465, 4281 (2017), 1608.04519

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  29. [29]

    Persistent topology of the reionisation bubble network. I: Formalism & Phenomenology

    W. Elbers and R. van de Weygaert, MNRAS486, 1523 (2019), 1812.00462

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  30. [30]

    Wilding, K

    G. Wilding, K. Nevenzeel, R. van de Weygaert, G. Veg- ter, P. Pranav, B. J. T. Jones, K. Efstathiou, and J. Feld- brugge, MNRAS507, 2968 (2021), 2011.12851

    work page arXiv 2021

  31. [31]

    Heydenreich, B

    S. Heydenreich, B. Br¨ uck, and J. Harnois-D´ eraps, A&A 648, A74 (2021), 2007.13724

    work page arXiv 2021

  32. [32]

    Elbers and R

    W. Elbers and R. van de Weygaert, MNRAS520, 2709 (2023), 2209.03948

    work page arXiv 2023

  33. [33]

    Bermejo, G

    R. Bermejo, G. Wilding, R. van de Weygaert, B. J. T. Jones, G. Vegter, and K. Efstathiou, MNRAS529, 4325 (2024), 2206.14655

    work page arXiv 2024

  34. [34]

    Calles, J

    J. Calles, J. H. T. Yip, G. Contardo, J. Nore˜ na, A. Rouhi- ainen, and G. Shiu, J. Cosmology Astropart. Phys.2025, 064 (2025), 2412.15405

    work page arXiv 2025

  35. [35]

    2019, https://arxiv.org/abs/1912.12699

    R. Braun, A. Bonaldi, T. Bourke, E. Keane, and J. Wagg, Anticipated performance of the square kilometre array – phase 1 (ska1)(2019), 1912.12699, URLhttps://arxiv. org/abs/1912.12699

    work page arXiv 2019

  36. [36]

    Q. Niu, Y. Li, Y. Xu, H. Guo, and X. Zhang, ApJ978, 145 (2025), 2407.18136

    work page arXiv 2025

  37. [37]

    C. A. Watkinson, S. Majumdar, J. R. Pritchard, and R. Mondal, MNRAS472, 2436 (2017), 1705.06284

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  38. [38]

    Quantifying the non-Gaussianity in the EoR 21-cm signal through bispectrum

    S. Majumdar, J. R. Pritchard, R. Mondal, C. A. Watkin- son, S. Bharadwaj, and G. Mellema, MNRAS476, 4007 (2018), 1708.08458

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  39. [39]

    21cm-line bispectrum as method to probe Cosmic Dawn and Epoch of Reionization

    H. Shimabukuro, S. Yoshiura, K. Takahashi, S. Yokoyama, and K. Ichiki, MNRAS458, 3003 (2016), 1507.01335

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  40. [40]

    Constraining the EoR model parameters with the 21cm bispectrum

    H. Shimabukuro, S. Yoshiura, K. Takahashi, S. Yokoyama, and K. Ichiki, MNRAS468, 1542 (2017), 1608.00372

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  41. [41]

    Studying topological structure in the epoch of reionization with 3D-Minkowski functionals of 21cm line fluctuations

    S. Yoshiura, H. Shimabukuro, K. Takahashi, and T. Mat- subara, MNRAS465, 394 (2017), 1602.02351

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  42. [42]

    S. Bag, R. Mondal, P. Sarkar, S. Bharadwaj, and V. Sahni, MNRAS477, 1984 (2018), 1801.01116

    work page internal anchor Pith review Pith/arXiv arXiv 1984

  43. [43]

    Kapahtia, P

    A. Kapahtia, P. Chingangbam, S. Appleby, and C. Park, J. Cosmology Astropart. Phys.2018, 011 (2018), 1712.09195

    work page arXiv 2018

  44. [44]

    Z. Chen, Y. Xu, Y. Wang, and X. Chen, ApJ885, 23 (2019), 1812.10333

    work page arXiv 2019

  45. [45]

    Kapahtia, P

    A. Kapahtia, P. Chingangbam, and S. Appleby, J. Cos- mology Astropart. Phys.2019, 053 (2019), 1904.06840

    work page arXiv 2019

  46. [46]

    Kapahtia, P

    A. Kapahtia, P. Chingangbam, R. Ghara, S. Appleby, and T. R. Choudhury, J. Cosmology Astropart. Phys. 2021, 026 (2021), 2101.03962

    work page arXiv 2021

  47. [47]

    S. K. Giri and G. Mellema, MNRAS505, 1863 (2021), 2012.12908

    work page arXiv 2021

  48. [48]

    Spina, C

    B. Spina, C. Porciani, and C. Schimd, MNRAS505, 3492 (2021), 2101.09288

    work page arXiv 2021