Pith. sign in

REVIEW 3 major objections 5 minor 100 references

Weak-lensing peak statistics can outperform shear two-point functions for key cosmological parameters.

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 13:12 UTC pith:L3MCKFDN

load-bearing objection WL peak redshift distribution is a genuinely strong probe, and the joint peak forecast beats a weakened shear-2PCF baseline; the paper is honest but the abstract overstates and the covariance/emulator-error treatment needs work. the 3 major comments →

arxiv 2607.29128 v1 pith:L3MCKFDN submitted 2026-07-31 astro-ph.CO

Cosmological constraining power of the redshifts, heights, and angular clustering of weak gravitational lensing peaks

classification astro-ph.CO
keywords gravitational lensing: weakcosmological parameterslarge-scale structurecosmological simulationsnon-Gaussian statisticspeaksdark energycosmic shear
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.

This paper forecasts what a Euclid-like weak-lensing survey would learn from the peaks of the lensing convergence field—the tallest bumps in the maps—compared with the standard shear two-point correlation functions. It argues that when peak heights, peak redshifts, and peak clustering are combined, the constraints on the matter density Ωm, the baryon density Ωb, and the dark-energy equation-of-state parameters w0 and wa are all tighter than those from shear two-point functions alone. The redshift distribution of high-signal peaks is the single most powerful statistic, driving the Ωm, w0, and wa constraints and breaking degeneracies that let the combination also probe Ωb and h. If correct, future lensing surveys can extract substantially more cosmological information from the same maps by adding peak statistics to their standard analysis.

Core claim

For the parameters Ωm, Ωb, w0, and wa, the combination of WL peak statistics provides a tighter constraint than the shear two-point correlation functions ξ±, with the redshift distribution of SNR>3 peaks being the most powerful individual statistic—driving the Ωm, w0, and wa constraints and breaking degeneracies that let the combination also probe Ωb and h. The height distribution and angular clustering are most sensitive to the amplitude of the primordial power spectrum, ln(10^10 A_s). Combining the three statistics also breaks degeneracies that no single statistic can resolve on its own, and the paper finds that smaller smoothing scales typically yield better constraints, with the figure o

What carries the argument

The central objects are weak lensing peaks—local maxima in the smoothed convergence field that mark massive collapsed structures. The analysis uses three summary statistics of these peaks: the height distribution dN/dκ, the redshift distribution dN/dz (for peaks above signal-to-noise 3), and the angular two-point correlation function ω(θ). The argument is carried by a forward-modelling pipeline: 100 dark-matter-only simulations spanning a 10-dimensional cosmological parameter space are used to train Gaussian-process emulators for each statistic, and a Bayesian likelihood with a covariance estimated from 80 fixed-cosmology simulations produces the forecast. The key mechanism in the result is

Load-bearing premise

The load-bearing assumption is that the covariance matrix measured from 80 simulations at one fixed cosmology (the D3A flat ΛCDM cosmology) accurately describes the scatter of the peak statistics across the whole 10-dimensional parameter space; the paper states this has not been tested for non-Gaussian statistics.

What would settle it

Estimate the covariance matrix from hypercube simulations at several distinct cosmologies (e.g., high and low Ωm, evolving dark energy, and non-zero DDM) and re-run the likelihood forecast; if the reported credible-interval sizes or the ranking of the three peak statistics change by more than the statistical error bars, the fixed-covariance assumption is invalidated.

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

If this is right

  • A Euclid-like survey that includes peak statistics in its likelihood can obtain tighter dark-energy constraints than cosmic-shear two-point analysis alone.
  • The redshift distribution of high-signal peaks should be treated as a first-class observable in lensing analyses, not just an auxiliary check.
  • Combining peak statistics with two-point functions improves constraints further, especially for Ωm, where the two probes have different degeneracies.
  • The smallest smoothing scale that can be reliably modelled gives the best constraints, so surveys should push to smaller angular scales where feasible.
  • Baryon density and h, which are usually inaccessible to lensing alone, can be constrained by the combined peak statistics.

Where Pith is reading between the lines

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

  • If the covariance's dependence on cosmology turns out to be significant, the ranking of the three peak statistics and the absolute widths of the constraints could shift; the paper itself flags this as an untested assumption.
  • The method's real-world power depends on accurate photometric redshifts for massive clusters, so its performance in practice hinges on cluster photo-z precision staying near the assumed 0.02(1+z) level.
  • The same forward-modelling framework could be extended to other non-Gaussian statistics, such as peak steepness or void statistics, to see whether further gains come from combining even more summary statistics.
  • The complementarity between the redshift and height distributions may help calibrate baryonic feedback in future hydrodynamical analyses, since the paper notes these statistics respond differently to baryonic physics.

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

3 major / 5 minor

Summary. The paper presents an idealized cosmological forecast for a Euclid-like weak lensing survey using three WL peak statistics: the height distribution, redshift distribution, and angular clustering of SNR>3 peaks. The analysis is based on a new 10D dark-matter-only hypercube of simulations (varying Ωm, fb, h, ns, As, w0, wa, Σmν, αs, Γdcdm), from which Gaussian-process emulators are trained. A multivariate t-distribution likelihood uses a covariance matrix estimated from 80 flat-ΛCDM simulations at the D3A cosmology. The reference data vector is an additional D3A simulation. The authors find that the three statistics are complementary, with the redshift distribution driving constraints on Ωm, w0, and wa; the height distribution and angular clustering constraining ln(10^10 As); and their combination probing Ωb and h. They report that the combined WL peak statistics outperform the (non-tomographic) shear 2PCFs for several parameters, and that combining multiple smoothing scales improves 2D FoMs by a factor ≈2. The paper is explicit about its idealized nature, omitting baryonic physics, intrinsic alignments, source clustering, and nuisance-parameter marginalization.

Significance. If the central claim holds, the paper provides a concrete motivation for Stage IV surveys to extract peak heights, peak redshifts, and peak clustering jointly from lensing maps, rather than relying solely on two-point shear statistics. The 10D hypercube and the emulator-based forecasting pipeline are valuable contributions, and the internal validation (mean FoB = 0.59σ for the well-constrained parameters in §5.2) demonstrates that the pipeline is self-consistent. The paper is honest about many of its limitations. However, the headline 'outperform' claim rests on two assumptions that are acknowledged but untested for non-Gaussian statistics: the cosmology-independence of the covariance matrix (§4.4.3) and the omission of emulator-error covariance in the fiducial 10D forecast (§5.2 vs §5.5). These assumptions affect the relative widths and FoMs that support the central conclusion, so the significance is conditional on addressing them.

major comments (3)
  1. [§5.2–§5.5, Appendix A, Appendix C] The fiducial 10D forecasts (§5.2–§5.4) use only the sample covariance from 80 D3A simulations and do not include emulator-error covariance, whereas the 2D smoothing-scale forecasts in §5.5 explicitly add an emulator-error covariance because the results were 'less stable to small variations in the emulation set-up.' Appendix A reports leave-one-out fractional errors up to ≈3% for angular clustering, and Appendix C states that individual nodes 'can be misestimated on average by several percent' and that multimodalities 'may indicate an insufficiently regular emulator.' Since the central claim is a relative comparison of peak statistics versus 2PCFs, omitting emulator error in the peak likelihood preferentially shrinks the peak contours. The 2PCF data vector is also emulated but is smoother and likely easier to emulate accurately, so the asymmetry could change the ranking. The authors shoul
  2. [§4.4.3] The covariance matrix is estimated from 80 flat-ΛCDM simulations at the fixed D3A cosmology and applied across the full 10D hypercube parameter space. The paper acknowledges that this approximation 'has not yet been studied for the non-Gaussian statistics we consider.' This is a load-bearing assumption: every reported constraint width, FoM, and the peak-versus-2PCF ranking derive from this covariance. The non-monotonic Ωm dependence of the angular clustering (§5.2.1, Fig. 3) already suggests that the scatter of the clustering statistic may vary across the hypercube. A concrete test — e.g., estimating the covariance from hypercube nodes in different parts of parameter space, or using an analytic cosmology-dependent covariance — is needed to justify the assumption or to quantify its impact. Without this, the absolute and relative constraint widths are uncertain.
  3. [§5.4, Abstract] The comparison to ξ± uses non-tomographic shear two-point correlation functions measured from the same maps (9 bins between 0.1 and 3 deg). The authors note in §5.4 that 'in reality, the cosmic shear signal is measured across different tomographic bins, meaning our forecast on the constraining power of the 2PCFs is a conservative estimate.' Yet the abstract states without qualification that 'WL peaks alone outperform the commonly used statistics.' In practice, 'commonly used statistics' include tomographic 2PCFs and 3×2pt analyses, which are substantially more powerful than the non-tomographic case considered here. The abstract and conclusions should be tempered to reflect the scope of the comparison, or the comparison should be extended to a tomographic 2PCF forecast. The body is transparent about this limitation, but the headline overstates the result.
minor comments (5)
  1. [§5.5, Table 2] The factor-of-≍2 FoM improvement from combining smoothing scales is obtained from 2D forecasts with eight parameters fixed at their true values. The conclusion (§6) states 'When combining different smoothing scales... the FoMs to typically improve by a factor ≈2' without repeating this caveat. Please qualify the conclusion or verify the improvement in the full 10D inference.
  2. [Figure 5 and Figure 6 captions] The caption of Fig. 5(c) says 'Angular clustering (top)' but the panel is the bottom one; the same appears in Fig. 6(c). Please correct.
  3. [§5.1] Typo: 'suppresses the growth of of matter perturbations' should read 'suppresses the growth of matter perturbations.'
  4. [§4.4.1 vs Appendix A] §4.4.1 states 'sub-per-cent-level relative emulator accuracy' but Appendix A reports individual angular-clustering nodes mispredicted by up to ≈3%. Please qualify the main-text statement to avoid contradiction.
  5. [Appendix A2] The height distribution is truncated at the 99th percentile of the κ distribution across training simulations. This is an ad hoc choice that removes the highest peaks; the paper justifies it by emulator stability but does not test the sensitivity of the cosmological constraints to this cut. A short robustness test or a caveat would be helpful.

Circularity Check

0 steps flagged

No circular derivation: the forecast is a self-consistency test of a forward model, not a definitional tautology.

full rationale

The paper's inference chain is a forward-modeled forecast. Emulators are trained on 100 hypercube simulations, a target data vector is generated from a separate D3A simulation using the same simulation framework, and the likelihood (Eq. 14) is evaluated against that target. There is no step in which a fitted parameter is renamed as a prediction. The covariance matrix is estimated from 80 simulations at the D3A cosmology and assumed approximately cosmology-independent; the paper explicitly states this is an approximation in Sec. 4.4.3 and flags that it has not been studied for non-Gaussian statistics. This is a modeling assumption, not a circular reduction. The self-citations to Broxterman et al. (2024, 2025) motivate the peak-selection and redshift-assignment methods; they are methodological inputs, and the current paper's cosmological constraints do not reduce to those citations. The leave-one-out emulator validation in Appendix A is an out-of-sample test, and the paper itself notes in Appendix C that emulator inaccuracies may affect some contours. The qualitative claim that peak statistics outperform the 2PCFs is derived from the relative widths of posteriors computed from the same simulated maps and covariance; it is a property of the statistics in this idealized model, not an identity. No equation or fitted parameter is equivalent to the claimed output by construction.

Axiom & Free-Parameter Ledger

7 free parameters · 7 axioms · 0 invented entities

No new physics entities are introduced; the ledger is dominated by statistical and modeling choices. The load-bearing assumptions are the cosmology-independent covariance (§4.4.3), the fidelity of the 100-node 10D emulators, the DMO reduction, and the idealized peak-redshift assignment. The paper is transparent about most of these, which is why they appear as stated axioms rather than hidden parameters.

free parameters (7)
  • SNR > 3 peak-selection threshold = 3
    Hand-chosen cut on the smoothed noise-realisations (Appendix A1) that selects which peaks enter all three statistics; justified by Broxterman et al. (2025), but not re-optimized for this 10D analysis.
  • Height-distribution upper cut at 99th percentile of κ = κ(p99) across training simulations
    The highest 1% of κ peaks are excluded from emulator training and inference because the emulator 'performed better without them' (Appendix A2). Post-hoc and affects the height-distribution constraints.
  • Photo-z scatter for peak redshifts = σ_z = 0.02(1+z), no catastrophic outliers
    Assumed precision for massive-cluster photo-zs, motivated by the KiDS-1000 bright sample (§4.2); this smoothing directly shapes the redshift distribution, the paper's most powerful statistic.
  • Hypercube w0–wa sampling strip = w0 ∈ [-1.2,-0.3]; wa = -3.7(w0+1.0) + [0,0.5]
    The dark-energy plane is sampled along a degeneracy direction aligned with CMB+DESI BAO+SNe (§3, Table 1); the posterior's w0–wa degeneracy direction is partly inherited from this imposed sampling, bracketing but also restricting the emulator's training support.
  • Priors on beyond-w0waCDM parameters = αs: Gaussian σ = 0.02; Γdcdm: half-Gaussian width 1
    Adopted 'to penalise strong deviations from these less common extensions' (§4.4.2); both posteriors largely recover the prior, so the 10D claim is effectively weaker than the framing suggests.
  • Emulator architecture = 15 PCA components, Matérn-3/2 kernel, white-noise kernel
    Fixed by hand for the fiducial setup (Appendix A); the response surface, and hence all reported constraint widths, depends on these choices, and the authors note sensitivity to small changes in §5.5.
  • Binning choices = 25 z-bins over [0,1.25]; 10 SNR bins; 12 angular bins [0.1,3] deg
    Hand-chosen bins for the three statistics (Appendix A); the reported FoMs are binning-dependent and the paper does not explore binning sensitivity.
axioms (7)
  • domain assumption Spatially flat universe (Ωk = 0)
    Assumed in §2 so that the comoving angular diameter distance reduces to the comoving distance; the lensing kernel equations (Eqs. 5–10) depend on it.
  • domain assumption Born and Limber approximations hold
    Lensing quantities are evaluated on unperturbed photon paths and line-of-sight correlations are factorized (§2, Eqs. 2 and 6). Standard in cosmic shear, but an approximation.
  • domain assumption Dark-matter-only simulations are adequate for comparing statistics
    Baryonic feedback is ignored although the paper calls it 'one of the dominant systematics in WL inferences' (§3); the claim that it does not affect relative differences between statistics is stated, not tested.
  • domain assumption Covariance matrix is cosmology-independent
    The covariance is estimated from 80 flat-ΛCDM realisations at D3A and applied across the 10D space (§4.4.3); the paper notes this 'has not yet been studied for the non-Gaussian statistics we consider.'
  • domain assumption Emulators trained on 100 hypercube nodes are faithful in 10D
    Leave-one-out accuracy averages are sub-percent, but individual angular-clustering nodes are mispredicted by up to ~3% (Appendix A4), and posterior multimodalities 'may indicate an insufficiently regular emulator' (Appendix C); emulator error is only propagated in the §5.5 2D forecasts.
  • domain assumption A unique 'primary redshift' can be assigned to each high-SNR peak
    The redshift distribution statistic assigns each peak the redshift of the overdensity dominating Eq. 9 (§4.2), a procedure that is only possible inside the simulation; its fidelity to real Euclid observations is untested.
  • domain assumption Gaussian shape noise and Euclid-like source distribution
    Noise is drawn from a Gaussian with σ_ε = 0.26√2 and n_gal = 30 arcmin⁻² (§4, Eq. 11), and the source n(z) is a Euclidean parameterization with mean redshift 0.9 (Eq. 8); intrinsic alignments, source clustering, and reduced shear are ignored.

pith-pipeline@v1.3.0-daily-deepseek · 28481 in / 22983 out tokens · 220584 ms · 2026-08-03T13:12:05.516621+00:00 · methodology

0 comments
read the original abstract

Weak gravitational lensing (WL) peaks probe non-Gaussian information of the large-scale distribution of matter that is not captured by two-point lensing statistics. We study the cosmological potential of the height distribution, redshift distribution, and angular clustering of high-valued WL peaks using a Bayesian inference approach that mimics a $\textit{Euclid}$ analysis. We use a forthcoming dark-matter-only hypercube, which varies cosmology in a ten-dimensional space, including evolving dark energy, neutrino mass, decaying dark matter, and the running of the scalar spectral index. We find the individual WL peak statistics to be complementary, as the redshift distribution best constrains the matter density, $\Omega_\mathrm{m}$, and the dark energy equation-of-state parameters, $w_0$, and $w_a$; the height distribution and angular clustering are most sensitive to the amplitude of the primordial power spectrum, $\ln(10^{10}A_\mathrm{s})$; while combining the three statistics allows us to also probe the baryon density, $\Omega_\mathrm{b}$, and the Hubble parameter, $h$. A comparison to the shear two-point correlation function demonstrates that WL peaks alone outperform the commonly used statistics, while even tighter constraints are obtained when combining all. Considering the 2-dimensional $\Omega_\mathrm{m}-\ln(10^{10}A_\mathrm{s})$ and $w_0-w_a$ parameter planes, we find that the redshift distribution outperforms the angular clustering and peak-height distributions. We study the impact of the smoothing scale and find that, typically, the smallest scales yield the best results and the figure of merit improves by a factor $\approx2$ when combining multiple scales.

Figures

Figures reproduced from arXiv: 2607.29128 by Elena Sellentin, Henk Hoekstra, Ian G. McCarthy, Jaime Salcido, Jeger C. Broxterman, John Helly, Joop Schaye, Konrad Kuijken, Matthieu Schaller, Naomi Schutte, Willem Elbers.

Figure 1
Figure 1. Figure 1: WL peak statistics used in the fiducial forecast with 1 arcmin smoothing. Left: redshift distribution of SNR > 3 peaks. Middle: height distribution. Right: angular correlation function of SNR > 3 peaks. The measurements from the hypercube simulations are shown as dashed light-grey curves, the covariance matrix realisations as red lines and the measurements from the additional simulation at the FLAMINGO D3A… view at source ↗
Figure 2
Figure 2. Figure 2: One-dimensional posteriors of the fiducial forecast using WL peak angular clustering (blue), height distribution (red), redshift distribution (orange), or the three statistics combined (black) for SNR > 3 WL peaks with 1 arcmin smoothing. The prior is indicated by the dotted grey curve and the shaded grey area, which uniformly spans the parameter range for all panels except the final two. The true values a… view at source ↗
Figure 3
Figure 3. Figure 3: Emulator predictions for the redshift distribution (left), height distribution (middle), and angular clustering (right) of SNR > 3 WL peaks at 1 arcmin smoothing for varying Ωm, as indicated by the colour bar, while keeping all other parameters fixed to their mean values within the hypercube. All three statistics are impacted. Larger values of Ωm increase the overall number of peaks at fixed WL convergence… view at source ↗
Figure 4
Figure 4. Figure 4: Corner plot of the cosmological parameter forecast using the WL two-point correlation functions (blue) or angular clustering, redshifts and heights of WL peaks (black) of SNR > 3 peaks or a combined inference using all 4 statistics (orange) at 1 arcmin smoothing. The contours indicate the 1- and 2𝜎 credible intervals. The combination of WL peak statistics outperforms the 𝜉+/− 2PCFs for Ωm, Ωb, 𝑤0, and 𝑤𝑎, … view at source ↗
Figure 5
Figure 5. Figure 5: Posteriors for the 2 dimensional forecast for Ωm − ln(1010𝐴s) using (a) the redshift distribution (top), (b) the height distribution (middle) or (c) the angular clustering (top) for 0.5′ (grey), 1.0′ (blue), 2.0′ (green), or the combination of all three (black) smoothing scales. The redshift distribution shows the clearest improvement when combining multiple smoothing scales. The height distribution shows … view at source ↗
Figure 6
Figure 6. Figure 6: Same as [PITH_FULL_IMAGE:figures/full_fig_p013_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.

Reference graph

Works this paper leans on

100 extracted references · 3 canonical work pages · 2 internal anchors

  1. [1]

    , keywords =

    The FLAMINGO project: baryonic impact on weak gravitational lensing convergence peak counts. , keywords =. doi:10.1093/mnras/stae698 , archivePrefix =. 2312.08450 , primaryClass =

  2. [2]

    Euclid preparation. XCVI. Cosmology Likelihood for Observables in Euclid (CLOE). 3. Inference and Forecasts. A&A, in press,. doi:10.48550/arXiv.2510.09153 , archivePrefix =. 2510.09153 , primaryClass =

  3. [3]

    arXiv e-prints , keywords =

    Comparing explicit likelihood and likelihood-free simulation-based inference for weak lensing cosmic shear. arXiv e-prints , keywords =

  4. [4]

    , keywords =

    KiDS-Legacy: Constraining dark energy, neutrino mass, and curvature. , keywords =. doi:10.1051/0004-6361/202558581 , archivePrefix =. 2512.11041 , primaryClass =

  5. [5]

    , keywords =

    The FLAMINGO project: cosmology with the redshift dependence of weak gravitational lensing peaks. , keywords =. doi:10.1093/mnras/staf357 , archivePrefix =. 2412.02736 , primaryClass =

  6. [6]

    , keywords =

    The FLAMINGO project: cosmological hydrodynamical simulations for large-scale structure and galaxy cluster surveys. , keywords =. doi:10.1093/mnras/stad2419 , archivePrefix =. 2306.04024 , primaryClass =

  7. [7]

    , keywords =

    FLAMINGO: calibrating large cosmological hydrodynamical simulations with machine learning. , keywords =. doi:10.1093/mnras/stad2540 , archivePrefix =. 2306.05492 , primaryClass =

  8. [8]

    , keywords =

    CORRFUNC - a suite of blazing fast correlation functions on the CPU. , keywords =. 2020. doi:10.1093/mnras/stz3157 , adsurl =

  9. [9]

    CORRFUNC: Blazing Fast Correlation Functions with AVX512F SIMD Intrinsics

    Sinha, Manodeep and Garrison, Lehman. CORRFUNC: Blazing Fast Correlation Functions with AVX512F SIMD Intrinsics. Software Challenges to Exascale Computing. 2019

  10. [10]

    , keywords =

    Bias and Variance of Angular Correlation Functions. , keywords =. doi:10.1086/172900 , adsurl =

  11. [11]

    10.1051/0004-6361/202450810

    Euclid - I. Overview of the Euclid mission , DOI= "10.1051/0004-6361/202450810", url= "https://doi.org/10.1051/0004-6361/202450810", journal =

  12. [12]

    preparation. XXVIII. Forecasts for ten different higher-order weak lensing statistics. , keywords =. doi:10.1051/0004-6361/202346017 , archivePrefix =. 2301.12890 , primaryClass =

  13. [13]

    , keywords =

    Big-Bang Nucleosynthesis after Planck. , keywords =. doi:10.1088/1475-7516/2020/03/010 , archivePrefix =. 1912.01132 , primaryClass =

  14. [14]

    , keywords =

    B-modes in cosmic shear from source redshift clustering. , keywords =. doi:10.1051/0004-6361:20020626 , archivePrefix =. astro-ph/0112441 , primaryClass =

  15. [15]

    , keywords =

    Constraints on neutrino physics from DESI DR2 BAO and DR1 full shape. , keywords =. doi:10.1103/w9pk-xsk7 , archivePrefix =. 2503.14744 , primaryClass =

  16. [16]

    , keywords =

    Three-point intrinsic alignments of galaxies and haloes in the FLAMINGO simulations. , keywords =. doi:10.1093/mnras/stag1008 , archivePrefix =. 2601.17914 , primaryClass =

  17. [17]

    , keywords =

    GetDist: a Python package for analysing Monte Carlo samples. , keywords =. doi:10.1088/1475-7516/2025/08/025 , archivePrefix =. 1910.13970 , primaryClass =

  18. [18]

    , keywords =

    SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications. , keywords =. doi:10.1093/mnras/stae922 , archivePrefix =. 2305.13380 , primaryClass =

  19. [19]

    MUSIC2-monofonIC: 3LPT initial condition generator

  20. [20]

    , keywords =

    Cosmology constraints from the weak lensing peak counts and the power spectrum in CFHTLenS data. , keywords =. doi:10.1103/PhysRevD.91.063507 , archivePrefix =. 1412.0757 , primaryClass =

  21. [21]

    , keywords =

    Revealing the neutrino mass through persistent homology of the cosmic web. , keywords =. doi:10.1093/mnras/stag937 , archivePrefix =. 2604.02300 , primaryClass =

  22. [22]

    arXiv e-prints , keywords =

    On combining estimated and analytic covariance matrices. arXiv e-prints , keywords =. doi:10.48550/arXiv.2604.19463 , archivePrefix =. 2604.19463 , primaryClass =

  23. [23]

    , keywords =

    Accurate initial conditions for cosmological N-body simulations: minimizing truncation and discreteness errors. , keywords =. doi:10.1093/mnras/staa3149 , archivePrefix =. 2008.09588 , primaryClass =

  24. [24]

    Euclid preparation. VII. Forecast validation for Euclid cosmological probes. , keywords =. doi:10.1051/0004-6361/202038071 , archivePrefix =. 1910.09273 , primaryClass =

  25. [25]

    , keywords =

    Lifting weak lensing degeneracies with a field-based likelihood. , keywords =. doi:10.1093/mnras/stab3234 , archivePrefix =. 2108.04825 , primaryClass =

  26. [26]

    , keywords =

    Effect of baryonic feedback on two- and three-point shear statistics: prospects for detection and improved modelling. , keywords =. doi:10.1093/mnras/stt1013 , archivePrefix =. 1210.7303 , primaryClass =

  27. [27]

    , keywords =

    Quantifying the effect of baryon physics on weak lensing tomography. , keywords =. doi:10.1111/j.1365-2966.2011.19385.x , archivePrefix =. 1105.1075 , primaryClass =

  28. [28]

    Strong and weak lensing united. I. The combined strong and weak lensing cluster mass reconstruction method. , keywords =. doi:10.1051/0004-6361:20042233 , adsurl =

  29. [29]

    , keywords =

    Detection of the significant impact of source clustering on higher order statistics with DES Year 3 weak gravitational lensing data. , keywords =. doi:10.1093/mnrasl/slad143 , archivePrefix =. 2307.13860 , primaryClass =

  30. [30]

    , keywords =

    The Effect of Intrinsic Alignments on Weak-lensing Statistics in Hydrodynamical Simulations. , keywords =. doi:10.3847/1538-4357/ae1ca7 , archivePrefix =. 2504.12460 , primaryClass =

  31. [31]

    , keywords =

    The impact of baryons on cosmological inference from weak lensing statistics. , keywords =. doi:10.1093/mnras/stab1978 , archivePrefix =. 2104.04165 , primaryClass =

  32. [32]

    International Journal of Modern Physics D , keywords =

    Accelerating Universes with Scaling Dark Matter. International Journal of Modern Physics D , keywords =. doi:10.1142/S0218271801000822 , archivePrefix =. gr-qc/0009008 , primaryClass =

  33. [33]

    , keywords =

    Exploring the Expansion History of the Universe. , keywords =. doi:10.1103/PhysRevLett.90.091301 , archivePrefix =. astro-ph/0208512 , primaryClass =

  34. [34]

    , keywords =

    The FLAMINGO project: the coupling between baryonic feedback and cosmology in light of the S _ 8 tension. , keywords =. doi:10.1093/mnras/staf093 , archivePrefix =. 2403.12967 , primaryClass =

  35. [35]

    , keywords =

    HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere. , keywords =. doi:10.1086/427976 , archivePrefix =. astro-ph/0409513 , primaryClass =

  36. [36]

    Computing in Science and Engineering , keywords =

    Matplotlib: A 2D Graphics Environment. Computing in Science and Engineering , keywords =. doi:10.1109/MCSE.2007.55 , adsurl =

  37. [37]

    Physics of the Dark Universe , keywords =

    Status of the S _ 8 tension: A 2026 review of probe discrepancies. Physics of the Dark Universe , keywords =. doi:10.1016/j.dark.2026.102286 , archivePrefix =. 2602.12238 , primaryClass =

  38. [38]

    , keywords =

    Ten-parameter simulation suite for cosmological emulation beyond CDM. , keywords =. doi:10.1103/PhysRevD.111.083529 , archivePrefix =. 2501.06296 , primaryClass =

  39. [39]

    The Coyote Universe. II. Cosmological Models and Precision Emulation of the Nonlinear Matter Power Spectrum. , keywords =. doi:10.1088/0004-637X/705/1/156 , archivePrefix =. 0902.0429 , primaryClass =

  40. [40]

    , keywords =

    A new approach to observational cosmology using the scattering transform. , keywords =. doi:10.1093/mnras/staa3165 , archivePrefix =. 2006.08561 , primaryClass =

  41. [41]

    arXiv e-prints , keywords =

    Constraints on Dynamical Dark Energy from Multiple Probes in the Full Dark Energy Survey. arXiv e-prints , keywords =. doi:10.48550/arXiv.2605.27221 , archivePrefix =. 2605.27221 , primaryClass =

  42. [42]

    The Journal of Open Source Software , keywords =

    healpy: equal area pixelization and spherical harmonics transforms for data on the sphere in Python. The Journal of Open Source Software , keywords =. doi:10.21105/joss.01298 , adsurl =

  43. [43]

    , keywords =

    Array programming with NumPy. , keywords =. doi:10.1038/s41586-020-2649-2 , archivePrefix =. 2006.10256 , primaryClass =

  44. [44]

    The Journal of Open Source Software , keywords =

    swiftsimio: A Python library for reading SWIFT data. The Journal of Open Source Software , keywords =. doi:10.21105/joss.02430 , adsurl =

  45. [45]

    , keywords =

    The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package. , keywords =. doi:10.3847/1538-4357/ac7c74 , archivePrefix =. 2206.14220 , primaryClass =

  46. [46]

    , keywords =

    emcee: The MCMC Hammer. , keywords =. doi:10.1086/670067 , archivePrefix =. 1202.3665 , primaryClass =

  47. [47]

    , keywords =

    Impact of baryonic feedback on HSC-Y1 weak lensing non-Gaussian statistics. , keywords =. doi:10.1103/PhysRevD.110.103539 , archivePrefix =. 2403.03807 , primaryClass =

  48. [48]

    , keywords =

    Effects of baryons on weak lensing peak statistics. , keywords =. doi:10.1088/1475-7516/2019/10/011 , archivePrefix =. 1905.11636 , primaryClass =

  49. [49]

    Experimental Astronomy , keywords =

    The Kilo-Degree Survey. Experimental Astronomy , keywords =. doi:10.1007/s10686-012-9306-1 , archivePrefix =. 1206.1254 , primaryClass =

  50. [50]

    Ground-based and Airborne Instrumentation for Astronomy IV , year = 2012, editor =

    Hyper Suprime-Cam. Ground-based and Airborne Instrumentation for Astronomy IV , year = 2012, editor =. doi:10.1117/12.926844 , adsurl =

  51. [51]

    arXiv e-prints , keywords =

    The Dark Energy Survey. arXiv e-prints , keywords =. doi:10.48550/arXiv.astro-ph/0510346 , archivePrefix =. astro-ph/0510346 , primaryClass =

  52. [52]

    , keywords =

    Cosmological forecast for non-Gaussian statistics in large-scale weak lensing surveys. , keywords =. doi:10.1088/1475-7516/2021/01/028 , archivePrefix =. 2006.12506 , primaryClass =

  53. [53]

    , keywords =

    KiDS+VIKING-450: Cosmic shear tomography with optical and infrared data. , keywords =. doi:10.1051/0004-6361/201834878 , archivePrefix =. 1812.06076 , primaryClass =

  54. [54]

    , keywords =

    Impact of magnification and size bias on the weak lensing power spectrum and peak statistics. , keywords =. doi:10.1103/PhysRevD.89.023515 , archivePrefix =. 1310.7517 , primaryClass =

  55. [55]

    , keywords =

    Parameter inference with estimated covariance matrices. , keywords =. doi:10.1093/mnrasl/slv190 , archivePrefix =. 1511.05969 , primaryClass =

  56. [56]

    arXiv e-prints , keywords =

    Dark Energy Survey Year 6 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing. arXiv e-prints , keywords =. doi:10.48550/arXiv.2601.14559 , archivePrefix =. 2601.14559 , primaryClass =

  57. [57]

    arXiv e-prints , keywords =

    KiDS-Legacy: Constraining dark energy, neutrino mass, and curvature. arXiv e-prints , keywords =. doi:10.48550/arXiv.2512.11041 , archivePrefix =. 2512.11041 , primaryClass =

  58. [58]

    , keywords =

    Weak gravitational lensing. , keywords =. doi:10.1016/S0370-1573(00)00082-X , archivePrefix =. astro-ph/9912508 , primaryClass =

  59. [59]

    Reports on Progress in Physics , keywords =

    Cosmology with cosmic shear observations: a review. Reports on Progress in Physics , keywords =. doi:10.1088/0034-4885/78/8/086901 , archivePrefix =. 1411.0115 , primaryClass =

  60. [60]

    , keywords =

    The fifth data release of the Kilo Degree Survey: Multi-epoch optical/NIR imaging covering wide and legacy-calibration fields. , keywords =. doi:10.1051/0004-6361/202346730 , archivePrefix =. 2503.19439 , primaryClass =

  61. [61]

    DESI DR2 results. II. Measurements of baryon acoustic oscillations and cosmological constraints. , keywords =. doi:10.1103/tr6y-kpc6 , archivePrefix =. 2503.14738 , primaryClass =

  62. [62]

    , keywords =

    DESI 2024 VII: cosmological constraints from the full-shape modeling of clustering measurements. , keywords =. doi:10.1088/1475-7516/2025/07/028 , archivePrefix =. 2411.12022 , primaryClass =

  63. [63]

    , keywords =

    KiDS-Legacy: Cosmological constraints from cosmic shear with the complete Kilo-Degree Survey. , keywords =. doi:10.1051/0004-6361/202554908 , archivePrefix =. 2503.19441 , primaryClass =

  64. [64]

    Planck 2018 results. VI. Cosmological parameters. , keywords =. doi:10.1051/0004-6361/201833910 , archivePrefix =. 1807.06209 , primaryClass =

  65. [65]

    , keywords =

    Persistent homology in cosmic shear: Constraining parameters with topological data analysis. , keywords =. doi:10.1051/0004-6361/202039048 , archivePrefix =. 2007.13724 , primaryClass =

  66. [66]

    Annual Review of Nuclear and Particle Science , keywords =

    Weak Gravitational Lensing and Its Cosmological Applications. Annual Review of Nuclear and Particle Science , keywords =. doi:10.1146/annurev.nucl.58.110707.171151 , archivePrefix =. 0805.0139 , primaryClass =

  67. [67]

    , keywords =

    New interpretable statistics for large-scale structure analysis and generation. , keywords =. doi:10.1103/PhysRevD.102.103506 , archivePrefix =. 2006.06298 , primaryClass =

  68. [68]

    Euclid preparation. LXXXV. Towards a DR1 application of higher-order weak lensing statistics. , keywords =. doi:10.1051/0004-6361/202557573 , archivePrefix =. 2510.04953 , primaryClass =

  69. [69]

    arXiv e-prints , keywords =

    FLAMINGO: Baryonic effects on the weak lensing scattering transform. arXiv e-prints , keywords =. doi:10.48550/arXiv.2510.09761 , archivePrefix =. 2510.09761 , primaryClass =

  70. [70]

    , keywords =

    Constraining neutrino mass with weak lensing Minkowski Functionals. , keywords =. doi:10.1088/1475-7516/2019/06/019 , archivePrefix =. 1812.08206 , primaryClass =

  71. [71]

    non-Gaussian random fields: graphs towards Betti numbers and persistence diagrams

    Stochastic homology of Gaussian vs. non-Gaussian random fields: graphs towards Betti numbers and persistence diagrams. , keywords =. doi:10.1088/1475-7516/2019/09/052 , archivePrefix =. 1908.01619 , primaryClass =

  72. [72]

    , keywords =

    Statistics of Dark Matter Halos from Gravitational Lensing. , keywords =. doi:10.1086/312480 , archivePrefix =. astro-ph/9910459 , primaryClass =

  73. [73]

    Gaussian Processes for Machine Learning

  74. [74]

    , keywords =

    Constraining Cosmology with High-Convergence Regions in Weak Lensing Surveys. , keywords =. doi:10.1088/0004-637X/691/1/547 , archivePrefix =. 0809.4052 , primaryClass =

  75. [75]

    , keywords =

    Searching for massive clusters in weak lensing surveys. , keywords =. doi:10.1111/j.1365-2966.2004.07691.x , archivePrefix =. astro-ph/0310607 , primaryClass =

  76. [76]

    , keywords =

    Probing cosmology with weak lensing peak counts. , keywords =. doi:10.1103/PhysRevD.81.043519 , archivePrefix =. 0907.0486 , primaryClass =

  77. [77]

    , keywords =

    Cosmology with the shear-peak statistics. , keywords =. doi:10.1111/j.1365-2966.2009.15948.x , archivePrefix =. 0906.3512 , primaryClass =

  78. [78]

    , keywords =

    Weak-lensing peak statistics - steepness versus height. , keywords =. doi:10.1093/mnras/stad534 , archivePrefix =. 2302.08255 , primaryClass =

  79. [79]

    , keywords =

    The cosmological information of shear peaks: beyond the abundance. , keywords =. doi:10.1093/mnras/stt552 , archivePrefix =. 1301.5001 , primaryClass =

  80. [80]

    , keywords =

    Cosmological forecasts with the clustering of weak lensing peaks. , keywords =. doi:10.1093/mnras/stac1204 , archivePrefix =. 2110.10164 , primaryClass =

Showing first 80 references.