REVIEW 4 major objections 4 minor 4 references
The splashback radius of galaxy clusters, measured two independent ways, agrees with cold-dark-matter predictions at 10–14% precision.
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 →
The splashback radius of AMICO KiDS-1000 clusters, measured with weak lensing and cluster-galaxy clustering, agrees with ΛCDM predictions, with clustering achieving 10% precision.
T0 review reviewed 2026-08-03 challenge →
load-bearing objection Solid two-probe splashback measurement on AMICO/KiDS-1000, but the quoted 10–14% precision is partly prior-driven and the model is not yet mock-tested. the 4 major comments →
AMICO galaxy clusters in KiDS-1000: Splashback radius from weak lensing and cluster-galaxy correlation function
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
For 8,730 rich clusters in the redshift range 0.1–0.8, the paper models the stacked reduced shear and the projected cluster-galaxy correlation function with the same truncated Einasto-plus-outer-power-law density profile, marginalising over selection effects, photometric-redshift scatter, miscentring, and the richness–mass relation. The fitted profiles place the average splashback radius at rsp/r200m at values consistent with the ΛCDM theoretical relation as a function of peak height, and give a mass accretion rate consistent with simulation predictions; the two probes agree within 1σ, with the clustering probe yielding a slightly smaller rsp, interpreted as dynamical friction on satellite g
What carries the argument
A single-cusp density profile with a transition term and an outer power law (Eq. 8), evaluated as an ensemble average over richness and redshift bins. The splashback radius is not measured directly but defined as the minimum of the logarithmic slope of the total density profile, and the ensemble-averaging step (Eq. 34) converts individual halo profiles into predicted observables free of radial binning effects.
Load-bearing premise
The analytic density profile (Eq. 8) describes the true stacked mass distribution from 0.4 to 5 h−1 Mpc and stays valid when extrapolated outward, so that the radius of minimum logarithmic slope of the fitted profile equals the true splashback radius of the halo population.
What would settle it
Measure the stacked density profile of the same clusters non-parametrically (e.g., by Abel-inverting the lensing signal with no assumed shape) and compare the radius of steepest slope with the profile-fitted rsp; if they disagree by more than the quoted 10–14% uncertainty, the result is model-dependent. Alternatively, apply the same pipeline to mock observations with known true splashback radii and check recovery within 1σ.
If this is right
- If true, optically selected clusters have a universal splashback boundary that tracks the ΛCDM Rsp–ν200m relation, with no residual trend in redshift or richness beyond mass.
- Galaxy clustering is a sharper splashback probe than lensing (10% vs 14% precision) and also constrains the amplitude and slope of the infalling-matter profile, which lensing alone cannot.
- The small but systematic offset between lensing and galaxy-traced rsp implies that galaxies trace a splashback boundary biased slightly inward by dynamical friction — a bias that will matter for any cluster-based cosmology using galaxy positions.
- The inferred mass accretion rates agree with simulation expectations, supporting the use of rsp as a mass-accretion probe with cosmological sensitivity to Ωm and σ8.
- The results confirm earlier X-ray, SZ, and optical measurements within 1–2σ, so the splashback radius is now measured consistently across cluster selection methods.
Where Pith is reading between the lines
- If the fitted-profile rsp depends on cosmology as the Planck18 test suggests (lower rsp for higher Ωm), then with the statistical power of upcoming surveys, rsp–ν200m constraints could become a competitive standalone cosmological probe.
- The lensing-vs-clustering offset can be tested directly in simulations with galaxy formation: if dynamical friction is the cause, the offset should grow with satellite galaxy mass and with the magnitude gap between the brightest galaxy and its satellites.
- Replacing the analytic outer power law with a matter-power-spectrum two-halo term, or measuring the profile non-parametrically, would test whether the extrapolation beyond 5 h−1 Mpc biases rsp; this is the most natural next step to check the model dependence.
- A joint fit of lensing and clustering with an explicit dynamical-friction parameter, rather than separate fits, could break the degeneracy between infall profile shape and galaxy bias and yield a single, robust rsp measurement.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes stacked weak-lensing reduced shear g_t and cluster-galaxy correlation function w_cg for 8730 AMICO clusters in KiDS-1000, binned in richness and redshift. Using a Diemer & Kravtsov (2014) profile with a transition factor and an outer infall term, the authors model the stacked observables and define the splashback radius as the minimum-log-slope radius of the 3D density profile. They report constraints on r_sp, the mass accretion rate Gamma, and the relation between R_sp = r_sp/r_200m and peak height nu_200m, claiming that the two probes are mutually consistent and agree with LambdaCDM simulation predictions, with per-stack precision of 14% (g_t) and 10% (w_cg).
Significance. If the model-derived r_sp indeed recovers the true splashback radius of the underlying halo population, this is a competitive and observationally valuable measurement, leveraging a large optical cluster sample and combining two complementary probes. The analysis is thorough in its treatment of known systematics: bootstrap and jackknife covariance matrices, propagation of shear-calibration, SOM redshift-distribution, miscentering, and mass-richness relation uncertainties, plus robustness tests against cosmology and model choices in Appendix C. The explicit comparison with the L25 mass calibration is a strength. However, the central claim depends on the DK14 profile shape and on priors calibrated with the same simulations used for the theoretical comparison; the paper's own appendices and discussion concede that several key parameters are prior-dominated and that end-to-end simulation validation is deferred. The significance of the quoted precision therefore remains conditional.
major comments (4)
- [§4.1–4.5, Eq. (36), Table 2, Fig. C.5] r_sp is not directly observed; it is the minimum-log-slope radius of the DK14 profile. The steepening that sets this radius is controlled by F_t, beta, and gamma_0, which are assigned Gaussian priors from DK14 (Sect. 4.6). Table 2 and Fig. C.5 show that the posteriors of these parameters remain close to the priors, with gamma_0 effectively unconstrained by either probe. Since the data only cover R in [0.4,5] h^-1 Mpc and the model is integrated to 40 h^-1 Mpc in Eq. (13), the claimed 10–14% precision per stack is substantially inherited from the simulation-calibrated profile shape rather than demonstrated by the data. The manuscript itself states in Sect. 6 that the impact of the DK14 extrapolation to very large scales 'shall be assessed through simulations.' An end-to-end mock validation—injecting clusters with known r_sp and verifying posterior recovery—is required to support the preci
- [§5, Eq. (45)] The mass accretion rate Gamma is obtained by inserting the model-derived R_sp into the More et al. (2015) fitting formula, and the comparison model of Diemer (2020) is calibrated on the same simulation suite. The text explicitly says the agreement 'is expected' for this reason. Therefore the Gamma constraints reported in Table 1 and Fig. 3 are not an independent test of LambdaCDM; they are a consistency check that is partly circular. The abstract and results should state this limitation, or the analysis should derive Gamma through an independent route before claiming a constraint.
- [§4.6] The w_cg analysis uses the g_t posteriors on the log lambda*–log M_200m relation (A, B, C, sigma_intr) as priors, and the text assumes M_200m = M_gt = M_wcg. The two probes are therefore not independent: a systematic error in the mass-richness calibration, or in the lensing masses, would shift both r_sp estimates coherently. The 'consistent results' claim in the abstract is weakened by this shared calibration. The authors should either run w_cg with uninformative mass-richness priors as a robustness check, or quantify the correlation between the two probe results.
- [§6] The paper acknowledges that anisotropic projection and selection effects may bias the w_cg measurements and that the impact on r_sp 'will be tested' with future dedicated mocks. Since w_cg provides the tighter constraint (10%), the central w_cg-based result—including the R_sp–nu relation and the possible dynamical-friction offset—rests on an unquantified systematic. A first-order assessment using the existing L25 anisotropic-boost model, or a simple test with mock galaxy catalogues, should be included before asserting that the w_cg constraints are unbiased.
minor comments (4)
- [§4.6, Eq. (40)] Typo: 'statical part of the covariance' should read 'statistical part.'
- [Appendix C] In the Planck18 robustness paragraph, there is a duplicated phrase: 'than the one assumed in our baseline analysis, analysis, namely Omega_m = 0.22.'
- [§4.2, Eq. (13)] The choice R_max = 40 h^-1 Mpc for the surface-density integration is not justified in the text. Given that the fitting range is [0.4,5] h^-1 Mpc and that Section 6 flags the large-scale extrapolation as a concern, a sentence explaining why 40 h^-1 Mpc is sufficient (or a convergence test) would help.
- [Fig. 2 caption] The caption states that error bars include 'residual uncertainties coming from systematic errors,' but the text (Sect. 4.6) models these as an additive covariance term rather than as error bars. Consider aligning the caption phrasing with the covariance treatment.
Circularity Check
Γ agreement is inherited from the same More+15/Diemer+20 simulation relation used to convert Rsp; r_sp model-prior dependence further weakens the ΛCDM comparison.
specific steps
-
self definitional
[Section 5, Eq. (45) and following paragraph]
"More et al. (2015) derived the following expression for the dimensionless mass accretion rate, Γ≡Δlog Mvir/Δlog a=0.935−3.04 ln(...)−1 ... We computed Γ posteriors by injecting ⟨Rsp(Δλ∗ob,Δzob)⟩ values, derived from Eq. (34) at each MCMC step, into Eq. (45). ... This agreement is expected, as our Rsp results are consistent with the ΛCDM predictions from More et al. (2015) and Diemer (2020a), and Eq. (45) is itself derived from those same simulations."
The paper reports a ΛCDM-agreeing Γ as one of its constraints, but Γ is not measured independently: it is obtained by inserting the fitted Rsp into Eq. (45), a mapping calibrated on the same More+15/Diemer+20 simulations that define the theoretical Γ prediction. The Γ agreement is therefore a restatement of the Rsp agreement through a simulation-derived conversion, not a new test of ΛCDM. The text itself concedes that the agreement is 'expected' because Eq. (45) comes from those simulations.
full rationale
The splashback radii and the Rsp–ν relation are obtained by fitting a DK14 profile (Eq. 8) to the stacked g_t and w_cg data and then locating the minimum-log-slope radius (Eq. 36); this is a model-dependent measurement, but the comparison with the More+15 and Diemer+20 predictions is not formally circular because the data can in principle shift the profile parameters away from the DK14 priors. The one step that does reduce by construction is the Γ constraint: Γ is computed by inserting the fitted Rsp into Eq. (45), a relation calibrated on the same ΛCDM simulations (More+15) that define the theoretical Γ prediction (Diemer+20), and the paper explicitly says the agreement is 'expected'. This makes the Γ agreement a restatement of the Rsp agreement rather than an independent confirmation. Two further factors lower the evidential weight without being formal circularity: the DK14 transition parameters (Ft, β, γ0) are prior-dominated (posteriors align with priors; Fig. C.5), and the w_cg analysis uses the g_t posteriors for the mass-richness relation, so the two probes are not fully independent. The central 14%/10% r_sp precision claims and the Rsp–ν comparison still contain independent data content.
Axiom & Free-Parameter Ledger
free parameters (18)
- F_t (truncation factor) =
1.40 (gt), 1.51 (wcg)
- β (transition sharpness) =
3.35 (wcg)
- γ0 (steepness normalization)
- b_e,0 (outer profile amplitude) =
1.97 (wcg)
- b_e,z (outer amplitude redshift scaling)
- s_e (outer profile slope) =
1.461 (wcg)
- f_off (miscentering fraction)
- σ_off (miscentering scale) =
0.21 h^-1 Mpc
- A (richness-mass amplitude) =
-0.27 (gt)
- B (richness-mass slope) =
0.59
- C (richness-mass redshift evolution) =
0.28
- σ_intr (richness intrinsic scatter) =
0.07
- log c0 (concentration amplitude) =
0.68 (gt), 0.82 (wcg)
- s, q (mass function correction)
- b_g,0 (galaxy bias amplitude) =
1.17 (wcg)
- b_g,λ* (galaxy bias richness trend) =
-0.10
- b_g,z (galaxy bias redshift trend) =
0.31
- A_sp, B_sp (Rsp-vs-ν relation) =
A=0.91, B=0.49 (gt); A=0.87, B=0.30 (wcg)
axioms (10)
- domain assumption DK14 profile (Eq. 8) describes the total mass profile of the stacked halo population, including the infalling region, over 0.4–40 h−1 Mpc.
- domain assumption The splashback radius is the radius of steepest slope of the total density profile (Eq. 36).
- domain assumption Tinker et al. (2008) mass function and Costanzi et al. (2019) bias correction are valid for this sample.
- domain assumption Richness–mass relation is lognormal with mean Eq. (26) and the observational scattering P(λ*_ob|λ*_tr) from L25 is correct.
- domain assumption Cluster photo-z uncertainty is Gaussian with σ=0.014(1+z).
- domain assumption SOM-reconstructed galaxy redshift distributions n(z_g) are unbiased.
- domain assumption Galaxy bias is constant along the line of sight and independent of projected radius.
- domain assumption Miscentering follows a Rayleigh distribution with scale σ_off.
- domain assumption The Γ–Rsp relation (Eq. 45) from More et al. (2015) is valid.
- domain assumption The c200m–M200m slope and redshift evolution are fixed to Duffy et al. (2008).
Cite this review
Pith. "Pith review of AMICO galaxy clusters in KiDS-1000: Splashback radius from weak lensing and cluster-galaxy correlation function." pith.science (2026). https://pith.science/paper/H23RX2HY
@misc{pith2026260207114,
author = {Pith},
title = {Pith review of: AMICO galaxy clusters in KiDS-1000: Splashback radius from weak lensing and cluster-galaxy correlation function},
year = {2026},
howpublished = {\url{https://pith.science/paper/H23RX2HY}},
note = {Machine review of arXiv:2602.07114}
}
abstract
We present the splashback radius analysis of the Adaptive Matched Identifier of Clustered Objects (AMICO) galaxy cluster sample in the fourth data release of the Kilo Degree Survey (KiDS). The sample contains 9049 rich galaxy clusters within $z\in[0.1,0.8]$, with shear measurements available for 8730 of them. We measure and model the stacked reduced shear, $g_{\rm t}$, and the cluster-galaxy correlation function, $w_{\rm cg}$, in bins of observed intrinsic richness, $\lambda^*$, and redshift, $z$. Building on the methods employed in recent cosmological analyses, we model the average splashback radius, $r_{\rm sp}$, of the underlying dark matter halo distribution, accounting for the known systematic uncertainties affecting measurements and theoretical models. By modelling $g_{\rm t}$ and $w_{\rm cg}$ separately, in the cluster-centric radial range $R\in[0.4,5]$ $h^{-1}$Mpc, we constrain $r_{\rm sp}$, the mass accretion rate, $\Gamma$, and the relation between $\mathcal{R}_{\rm sp}\equiv r_{\rm sp}/r_{200\rm m}$ and the peak height, $\nu_{200\rm m}$, over the mass range $M_{200\rm m}\in[0.4,20]$ $10^{14}h^{-1}$M$_\odot$. The two probes provide consistent results that also agree with $\Lambda$-cold dark matter model predictions. Our $\mathcal{R}_{\rm sp}$ constraints are consistent with those from previous observations. For $g_{\rm t}$ and $w_{\rm cg}$, we achieve a precision of 14% and 10% per cluster stack, respectively. The higher precision of $w_{\rm cg}$, enabled by its combination with weak-lensing constraints on the mass-richness relation, highlights the complementarity of lensing and clustering in measuring $r_{\rm sp}$ and constraining the properties of the infalling material region.
Figures
Reference graph
Works this paper leans on
-
[1]
Adelberger, K. L., Steidel, C. C., Pettini, M., et al. 2005, ApJ, 619, 697 Adhikari, S., Dalal, N., & Chamberlain, R. T. 2014, J. Cosm. Astro-Particle Phys., 2014, 019 Aihara, H., Allende Prieto, C., An, D., et al. 2011, ApJ Suppl., 193, 29 Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJ Suppl., 219, 12 Baltz, E. A., Marshall, P., & Oguri,...
Pith/arXiv arXiv 2005
-
[3]
54 in L25) and those from L25
Figure C.1 displays the comparison between the individual cluster mass estimates derived from theg t modelling presented in this work (based on Eq. 54 in L25) and those from L25. De- spite L25 adopting the truncated NFW profile by Baltz, Mar- shall, & Oguri (2009, BMO) rather than the DK14 profile in Eq. (8), and modelling theg t profiles only up to 3.5h−...
2009
-
[4]
Figure C.3 shows that the two probes yield similar constraints on these parameters
by modellingg t, subsequently using the resulting posteriors as priors for thew cg modelling. Figure C.3 shows that the two probes yield similar constraints on these parameters. The most significant differences appear in the pos- teriors for the amplitude,A, and the intrinsic scatter,σ intr. For wcg, theAposterior is shifted to lower values and is more he...
2008
-
[25]
Shaded areas represent 68% confidence regions
from the modelling ofg t (blue solid lines) andw cg (orange dashed lines). Shaded areas represent 68% confidence regions. 1014 1015 M200m [h 1M ] 10 2 3 4 6 c200m Duffy+08 gt wcg Fig. C.4.logc 200m−logM 200 relation constrained byg t (blue band) and wcg (orange band) measurements. The width of the bands represents the 68% confidence of the models, while t...
2008
This paper was first reviewed by deepseek-v4-flash on August 3, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.