Intracluster light as a dark matter tracer: how their spatial and kinematic relationship is shaped by satellite demographics

2); (2) Laboratory of Astrophysics; (3) Departamento de Fisica Teorica; 4; (4) CIAFF; 5); (5) Institute for Astronomy; (6) Institut d Astrophysique de Paris; 7) ((1) School of Physics; (7) ICRAR

arxiv: 2601.18693 · v2 · submitted 2026-01-26 · 🌌 astro-ph.GA · astro-ph.CO

Intracluster light as a dark matter tracer: how their spatial and kinematic relationship is shaped by satellite demographics

G Martin (1) , F R Pearce (1) , N A Hatch (1) , H J Brown (1) , J Butler (1) , Y M Bahe (1 , 2) , W Cui (3

show 24 more authors

4 5) Y Dubois (6) A Knebe (3 7) ((1) School of Physics Astronomy University of Nottingham UK (2) Laboratory of Astrophysics EPFL Switzerland (3) Departamento de Fisica Teorica Universidad Autonoma de Madrid Spain (4) CIAFF (5) Institute for Astronomy Royal Observatory Edinburgh (6) Institut d Astrophysique de Paris CNRS Sorbonne Universite France (7) ICRAR University of Western Australia Australia)

This is my paper

Pith reviewed 2026-05-16 11:07 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO

keywords intracluster lightdark mattergalaxy clusterstidal strippingsatellite galaxiesN-body simulationsphase space

0 comments

The pith

The radial relationship between intracluster light and dark matter is governed by the masses of infalling satellites.

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

The paper uses N-body simulations to track how stars stripped from satellite galaxies form intracluster light and sit relative to the cluster dark matter halo. Stripped stars consistently reach lower-energy and lower-angular-momentum orbits than the dark matter from the same satellites, with the size of the difference growing for higher satellite-to-host mass ratios. A model built from an entire population of infalling satellites shows that the resulting phase-space offset and radial concentration difference are set mainly by the characteristic mass of the satellite stellar mass function. When that mass function is matched to full hydrodynamical simulations, the predicted ICL-to-DM radial offsets agree within the scatter among the simulations. This indicates that ICL density profiles can trace the underlying dark matter distribution once the infalling satellite population is known.

Core claim

In controlled N-body simulations that vary satellite-to-host mass ratio and orbital circularity, the stripped stellar material from infalling satellites occupies lower specific orbital energy and angular momentum regions than the stripped dark matter. The authors construct a predictive model for the phase-space properties of both components across a whole satellite population and find that the offsets are driven primarily by the characteristic mass of the infalling satellite stellar mass function. The resulting ICL is always more centrally concentrated than the DM, with the magnitude of the offset increasing toward higher characteristic masses. Comparisons to four independent cosmological,

What carries the argument

Predictive model for the phase-space properties of stripped stars and dark matter from the infalling satellite population, with the characteristic mass of the satellite stellar mass function as the controlling parameter.

If this is right

The ICL is always more centrally concentrated than the DM halo.
The magnitude of the ICL-DM radial offset increases with the characteristic mass of the infalling satellites.
Orbital circularity has only weak influence on the phase-space offset compared with mass ratio.
ICL density profiles can serve as tracers of the radial DM distribution once the infalling satellite population is adequately constrained.
Matching the satellite stellar mass function alone reproduces the radial offsets seen in full hydrodynamical simulations within inter-simulation scatter.

Where Pith is reading between the lines

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

ICL observations in clusters could be used to infer properties of the satellite population that assembled the system.
The same relation might allow ICL measurements to help map DM profiles in clusters where direct mass estimates are limited.
Extending the model with baryonic physics would test whether gas and star formation change the predicted offsets.
Analogous stripping processes could produce usable intragroup light tracers in lower-mass systems.

Load-bearing premise

That dark-matter-only N-body simulations without gas or star formation accurately reproduce the tidal stripping and phase-space evolution of satellites once only the stellar mass function is matched.

What would settle it

Independently measure the radial ICL profile and the radial DM profile in a real galaxy cluster, determine the infalling satellite stellar mass function from observations or other data, and test whether the observed ICL-DM offset matches the model's prediction for that mass function.

read the original abstract

We investigate how the orbital evolution and mass distribution of infalling satellite galaxies shape the phase-space and radial distributions of intracluster light (ICL) relative to the underlying cluster dark matter (DM) halo. Using N-body simulations, we follow the tidal stripping and orbital evolution of satellite galaxies as they are accreted into a live cluster halo, systematically varying satellite-to-host mass ratio and orbital circularity. We measure the specific orbital energy and angular momentum of stripped stellar and DM material, finding that the stripped stars consistently occupy lower-energy and lower-angular momentum regions of phase-space than the stripped DM. The magnitude of this difference increases strongly towards more equal satellite--to--host mass ratios, while the dependence on orbital circularity is weak. We construct a predictive model for the phase-space properties of stripped stars and DM from a whole infalling satellite population and find that the resulting phase-space difference between the components are driven primarily by the characteristic mass of the infalling satellite stellar mass function. We find that the ICL is always more centrally concentrated than the DM. The magnitude of this offset depends on the characteristic mass and increases towards higher characteristic masses. Comparisons with four independent cosmological hydrodynamical simulations show that, once the infalling satellite stellar mass function is matched, the model reproduces the radial stellar-to-DM density profile offsets to better than the inter-simulation scatter. This demonstrates that the radial relationship between the ICL and the DM distribution is largely governed by satellite demographics. With adequate constraints on the infalling satellite population, ICL density profiles can therefore be used as informative tracers of the underlying radial DM distribution in clusters.

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.

Desk Editor's Note private letter to a colleague

Satellite demographics control the ICL-DM offset, with the N-body model holding up against hydro simulations.

read the letter

The main point is that the offset between intracluster light and dark matter in clusters is set mostly by the mass distribution of infalling satellites. Their N-body simulations show stars from more massive satellites end up more centrally concentrated and lower in phase space than the dark matter, and a model based on the characteristic satellite mass reproduces the radial differences seen in four separate hydrodynamical simulations once the mass function is matched. What the paper does well is the systematic scan over satellite-to-host mass ratio and orbital circularity. The mass ratio drives the effect strongly while circularity does not. Building the predictive model from the population and validating it directly against independent hydro runs is a solid check that the N-body setup captures the essential dynamics for this question. The soft spots are minor. The simulations omit gas and star formation, but the hydro comparison indicates this does not change the main trends. I would want to see the exact details of how they match the satellite mass function and any error bars on the offsets, since the abstract is light on those. The result also depends on the definition of stripped material, which could be checked for robustness. This work is aimed at cluster simulators and observers interested in ICL as a dark matter tracer. It deserves a serious referee because the multi-code validation gives it a concrete foundation.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that the radial and phase-space relationship between intracluster light (ICL) and dark matter (DM) in galaxy clusters is largely governed by the demographics of infalling satellites. Using N-body simulations that systematically vary satellite-to-host mass ratio and orbital circularity, the authors show that stripped stars occupy lower-energy, lower-angular-momentum regions than stripped DM, with the offset increasing for more equal-mass satellites. They construct a predictive model driven primarily by the characteristic mass of the infalling satellite stellar mass function, demonstrate that ICL is always more centrally concentrated than DM, and show that this model reproduces the stellar-to-DM density offsets seen in four independent hydrodynamical simulations once the satellite mass function is matched, to better than the inter-simulation scatter.

Significance. If the central result holds, the work provides a physically motivated framework for interpreting ICL density profiles as tracers of the underlying DM distribution once the infalling satellite population is adequately constrained, with direct relevance to cluster cosmology and galaxy evolution studies. The manuscript is strengthened by its systematic N-body parameter exploration and by the explicit cross-validation against multiple independent hydrodynamical simulations, which reproduces observed offsets without requiring gas or star-formation physics once the satellite mass function is matched.

major comments (2)

[Abstract and comparison section] Abstract and comparison section: the claim that the predictive model reproduces the radial stellar-to-DM density offsets 'to better than the inter-simulation scatter' is central to the validation but is presented without quantitative values for the offsets, the inter-simulation scatter, error bars on the profiles, or the precise matching procedure and data-exclusion criteria applied to the hydrodynamical runs.
[Predictive-model construction] Predictive-model construction (inferred §4): the statement that phase-space differences are 'driven primarily by the characteristic mass of the infalling satellite stellar mass function' requires an explicit derivation or equation showing how the population-level averages are computed from the individual satellite runs, including any weighting by the mass function and assumptions about orbital-parameter independence.

minor comments (2)

[Figures] All figures showing density profiles or phase-space distributions should include uncertainty estimates or bootstrap errors to allow direct assessment of the reported offsets.
[Methods] Clarify the exact number of N-body realizations per parameter combination and any convergence tests performed on the stripped-material measurements.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive assessment of its significance. We address each major comment below and will make the suggested revisions to improve the clarity and rigor of the presentation.

read point-by-point responses

Referee: [Abstract and comparison section] Abstract and comparison section: the claim that the predictive model reproduces the radial stellar-to-DM density offsets 'to better than the inter-simulation scatter' is central to the validation but is presented without quantitative values for the offsets, the inter-simulation scatter, error bars on the profiles, or the precise matching procedure and data-exclusion criteria applied to the hydrodynamical runs.

Authors: We acknowledge that the abstract and the comparison section would be strengthened by the inclusion of quantitative metrics. In the revised manuscript, we will add specific values for the stellar-to-DM density offsets, the inter-simulation scatter, error bars on the profiles, and a precise description of the matching procedure and data-exclusion criteria used for the hydrodynamical simulations. These additions will provide a more transparent validation of the model. revision: yes
Referee: [Predictive-model construction] Predictive-model construction (inferred §4): the statement that phase-space differences are 'driven primarily by the characteristic mass of the infalling satellite stellar mass function' requires an explicit derivation or equation showing how the population-level averages are computed from the individual satellite runs, including any weighting by the mass function and assumptions about orbital-parameter independence.

Authors: We agree that an explicit derivation is needed. We will include in the revised §4 a detailed equation and explanation of how the population-level phase-space averages are obtained by weighting the results from individual satellite runs according to the satellite stellar mass function. We will also state the assumption of independence from orbital parameters, justified by the weak dependence on circularity found in our simulations. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper derives phase-space offsets between stripped stars and DM directly from N-body simulations that systematically vary satellite-to-host mass ratio and orbital circularity. It then constructs a predictive model whose output is governed by the characteristic mass of the infalling satellite stellar mass function and validates that model against four independent hydrodynamical simulations by matching only the satellite SMF. The central claim that ICL radial profiles trace DM via satellite demographics is therefore supported by this external cross-validation rather than by any reduction of the result to fitted inputs or self-citations by construction. No load-bearing step matches the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

3 free parameters · 1 axioms · 0 invented entities

The central claim rests on N-body simulation results for tidal stripping and the assumption that satellite demographics alone explain the ICL-DM offset when matched across simulation types.

free parameters (3)

satellite-to-host mass ratio
Systematically varied in simulations to measure dependence of energy and angular momentum differences
orbital circularity
Systematically varied in simulations to measure dependence of energy and angular momentum differences
characteristic mass of infalling satellite stellar mass function
Identified as the primary driver of phase-space and radial offsets in the predictive model

axioms (1)

domain assumption N-body simulations accurately capture tidal stripping and orbital evolution of stellar and DM material from satellites
Core method used to follow material as satellites are accreted into the live cluster halo

pith-pipeline@v0.9.0 · 5753 in / 1390 out tokens · 43181 ms · 2026-05-16T11:07:30.332729+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/DimensionForcing.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We construct a predictive model for the phase-space properties of stripped stars and DM from a whole infalling satellite population and find that the resulting phase-space difference between the components are driven primarily by the characteristic mass of the infalling satellite stellar mass function.

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.