pith. sign in

arxiv: 2604.22920 · v1 · submitted 2026-04-24 · 🌌 astro-ph.GA

The abundance and radial distribution of faint and ultra-faint dwarfs in galaxy clusters

Pradyumna Sadhu , Laura V. Sales , Julio F. Navarro , Rapha\"el Errani , Jose Benavides , Eric W. Peng This is my paper

Pith reviewed 2026-05-08 10:43 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords galaxy clustersdwarf galaxiesultra-faint dwarfstidal evolutioncosmological simulationsradial distributionsatellite galaxies
0
0 comments X

The pith

Galaxy clusters with virial mass around 10^14 solar masses contain 2000 to 7000 faint dwarf galaxies with stellar masses above 100 solar masses.

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

The paper overcomes the resolution limits of cosmological simulations by combining the TNG50 hydrodynamical run with an empirical model of tidal evolution calibrated on high-resolution N-body simulations. This combination reveals a large population of dwarf galaxies well below the simulation's formal limit. The dwarfs as a group follow a radial distribution that matches the dark matter profile of the host cluster. However, realistic luminosity or mass thresholds for detection exclude the most heavily stripped galaxies, which are concentrated toward the center. This prediction can be tested with deeper observations of cluster environments.

Core claim

At z=0, clusters with virial mass M200 around 10^14 solar masses host 2000 to 7000 dwarf galaxies with stellar mass above 100 solar masses within the virial radius. These satellites together follow a radial distribution that matches the underlying dark matter profile of the host. Applying a minimum mass or luminosity threshold, as expected in observational studies, tends to exclude the most heavily stripped objects that populate the inner regions.

What carries the argument

An empirical model of tidal evolution calibrated on high-resolution idealized N-body simulations, applied to sub-resolution dwarf galaxies in the TNG50 hydrodynamical simulation to predict their stellar masses and radial positions.

If this is right

  • Clusters with M200 around 10^14 solar masses host 2000-7000 systems with stellar mass above 100 solar masses within the virial radius.
  • The combined population of these satellites follows a radial distribution that matches the host dark matter profile.
  • Luminosity or mass thresholds for detection exclude the most stripped inner objects and bias the observed distribution.

Where Pith is reading between the lines

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

  • Deep surveys could detect a higher central concentration of ultra-faint galaxies than current observations show.
  • The total number of satellites in clusters is likely underestimated when relying only on resolved simulation objects.
  • The same modeling approach could be used to predict dwarf populations in groups of lower mass.

Load-bearing premise

The empirical model of tidal evolution calibrated on idealized N-body simulations applies without major systematic bias to the sub-resolution dwarf galaxies in the TNG50 simulation.

What would settle it

A deep survey that counts ultra-faint dwarfs down to stellar masses of 100 solar masses across a galaxy cluster and checks whether their radial distribution matches the dark matter profile or shows a central deficit from selection effects.

Figures

Figures reproduced from arXiv: 2604.22920 by Eric W. Peng, Jose Benavides, Julio F. Navarro, Laura V. Sales, Pradyumna Sadhu, Rapha\"el Errani.

Figure 2
Figure 2. Figure 2: Tidal tracks for the evolution of baryonic properties – projected (2D) half-light radius 𝑅h (upper panel) and stellar mass 𝑀∗ (lower panel) – of satellite galaxies as a function of amount of stripping measured in terms of the fraction of mass remaining inside 𝑟mx, 𝑀mx/𝑀mx0. 𝑀mx here denotes the mass enclosed within the radius of peak circular velocity 𝑟mx, 𝑀mx = 𝑟mx𝑉 2 mx/𝐺 (while the subscript 0 refers to… view at source ↗
Figure 3
Figure 3. Figure 3: Top Left: The abundance matching relation for TNG50 simulation along with two chosen extrapolations adopted for the initial stellar mass - halo mass relation. Grey dots show maximum stellar mass (𝑀∗,max ) vs peak maximum circular velocity (𝑉peak ) for Type-1 satellites from FoFs-0, 1 and 2 in TNG50. Subhalos without stellar particles in the simulation are artificially shown at 𝑀∗ = 103M⊙. Horizontal dotted… view at source ↗
Figure 4
Figure 4. Figure 4: The initial stellar segregation, defined as the ratio between the projected half light radius and 𝑟mx. for the stars in their dark matter subhalos as a function of the infall stellar mass (𝑀∗0 ). Here, unresolved subhalos are populated using the power-law extrapolation. 𝑅 pl h0 is the projected infall stellar half-light radius, whereas 𝑟mx0 is the radius at which the circular velocity profile of the subhal… view at source ↗
Figure 5
Figure 5. Figure 5: Projected thin slice of FoF-0 from TNG50-1 in the xy-plane, excluding the central subhalo for the simulation versus power-law and cutoff models. The region shown corresponds to √︁ 𝑥 2 + 𝑦 2 < 300 kpc and |𝑧 | < 50 kpc in cluster-centric coordinates. This region is encircled by a green dashed circle. Left panel: Stellar particles at 𝑧 = 0 are binned to represent the surface brightness of the satellite galax… view at source ↗
Figure 6
Figure 6. Figure 6: 2D histogram for the combined size-mass relation for subhalos from all three FoF halos in TNG. The left and right panels correspond to the power-law and cutoff models, respectively. The dashed grey lines show lines of constant surface brightness. Yellow solid line shows the best-fit double power-law to observations of the galaxies belonging to Local Group and Virgo clusters. Most of the subhalos evolve alo… view at source ↗
Figure 7
Figure 7. Figure 7: Left:The projected satellite mass function for all three FoF halos. Power-law models are shown in red, and cutoff models are shown in blue. Dotted lines represent the satellite mass functions for the resolved Type-1 satellites in the TNG50 simulation. Different transparencies represent each FoF halo. Right: projected median satellite stellar mass function for the three FoF groups considered within 309 kpc … view at source ↗
Figure 8
Figure 8. Figure 8: Upper panel: predictions for the satellite mass function in FoF￾0 with surface brightness cuts at 24, 28 and 35 mag/arcsec2 (dash-dotted, dashed and solid lines, respectively). Lower panel: Completeness, defined as the fraction of satellites at a given 𝑀∗ that satisfies the surface brightness cut. Current surveys such as the SDSS can probe up to 24 mag/arcsec2 , and are thus complete down to a stellar mass… view at source ↗
Figure 9
Figure 9. Figure 9 view at source ↗
read the original abstract

Cosmological simulations of galaxy clusters are unable to resolve dwarf galaxies due to limited numerical resolution which drives the artificial disruption of dark matter substructures. We address these limitations by combining the results of the cosmological hydrodynamical simulation TNG50 in $\Lambda$CDM with an empirical model of tidal evolution of cluster galaxies calibrated using high-resolution idealized N-body simulations. Applied to the three most massive clusters in TNG50, our model allows us to study the stellar mass and radial distribution of dwarfs well below the formal resolution limit of the parent simulation. We find that, at $z=0$, clusters with virial mass $M_{200} \sim 10^{14}~\mathrm{M_\odot}$ host a vast population of dwarf galaxies within the virial radius, amounting to $2000$-$7000$ systems with $M_* > 100~\mathrm{M_\odot}$. Taken together, these satellites follow a radial distribution that matches the underlying dark matter profile of the host. However, applying a minimum mass or luminosity threshold for detection, as expected in observational studies, tends to exclude the most heavily-stripped objects, which tend to populate the inner regions. Future surveys targeting ultra-faint galaxies in group and cluster environments, such as those made possible by the Euclid, Rubin, or Roman telescopes, will be fundamental to refute or confirm this prediction.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript presents a hybrid approach combining the TNG50 cosmological hydrodynamical simulation with an empirical model of tidal evolution calibrated on high-resolution idealized N-body simulations to predict the abundance and radial distribution of faint and ultra-faint dwarf galaxies in galaxy clusters below the resolution limit of TNG50. For clusters with M_200 ~ 10^14 M_sun, it finds 2000-7000 dwarfs with M_* > 100 M_sun within the virial radius, whose collective radial distribution matches the host dark matter profile, though luminosity thresholds would preferentially exclude the most stripped inner satellites.

Significance. If the central results hold, this work provides a valuable prediction for the population of ultra-faint dwarfs in clusters, which has implications for the interpretation of current and future observations with telescopes like Euclid, Rubin, and Roman. The method addresses a key limitation in cosmological simulations by leveraging external high-resolution calibrations, and the external grounding of the tidal model on separate N-body runs is a strength that reduces circularity concerns.

major comments (3)
  1. [Methods] Methods section (description of empirical model application): The tidal evolution prescription is calibrated exclusively on idealized dark-matter-only N-body simulations and applied directly to sub-resolution galaxies in the hydrodynamical TNG50 run. No quantitative validation or comparison is shown demonstrating that the model reproduces the stellar mass function, survival rates, or radial positions of the resolved (above-resolution) dwarfs in TNG50 prior to extrapolation. This assumption is load-bearing for both the 2000-7000 abundance estimate and the claim that the satellites match the dark matter radial profile.
  2. [Results] Results section (radial distribution analysis): The conclusion that the satellites follow a radial distribution matching the host dark matter profile depends on the tidal model not introducing systematic biases in orbital decay or inner survival. The manuscript does not test or discuss whether baryonic processes present in TNG50 (gas stripping, feedback) alter mass-loss efficiency or positions relative to the dark-matter-only calibration, particularly for the most stripped objects that dominate the inner regions.
  3. [Results] Results or Discussion (abundance and uncertainty): The quoted range of 2000-7000 systems with M_* > 100 M_sun is derived from the three most massive TNG50 clusters, but the manuscript provides no details on error propagation from the empirical model parameters, variation across the three clusters, or sensitivity to the minimum mass threshold. This weakens the robustness of the quantitative prediction.
minor comments (2)
  1. [Abstract/Introduction] The abstract and introduction could more explicitly reference prior works on subhalo disruption in clusters to better contextualize the novelty of the hybrid method.
  2. [Figures] Figure captions for the radial distribution plots should clarify whether the dark matter profile is normalized to the same total mass or number of objects as the satellite sample.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us identify areas for improvement. We respond to each major comment below.

read point-by-point responses

  1. Referee: [Methods] Methods section (description of empirical model application): The tidal evolution prescription is calibrated exclusively on idealized dark-matter-only N-body simulations and applied directly to sub-resolution galaxies in the hydrodynamical TNG50 run. No quantitative validation or comparison is shown demonstrating that the model reproduces the stellar mass function, survival rates, or radial positions of the resolved (above-resolution) dwarfs in TNG50 prior to extrapolation. This assumption is load-bearing for both the 2000-7000 abundance estimate and the claim that the satellites match the dark matter radial profile.

    Authors: We agree that a direct comparison would strengthen the case for extrapolation. The empirical model is applied exclusively to sub-resolution objects; resolved satellites in TNG50 are taken directly from the simulation outputs. To address the concern, we will add a new validation subsection in the Methods that applies the tidal model to the initial conditions of resolved TNG50 satellites and compares the resulting stellar mass function, survival rates, and radial distributions against the native TNG50 results at z=0, quantifying the level of agreement. revision: yes

  2. Referee: [Results] Results section (radial distribution analysis): The conclusion that the satellites follow a radial distribution matching the host dark matter profile depends on the tidal model not introducing systematic biases in orbital decay or inner survival. The manuscript does not test or discuss whether baryonic processes present in TNG50 (gas stripping, feedback) alter mass-loss efficiency or positions relative to the dark-matter-only calibration, particularly for the most stripped objects that dominate the inner regions.

    Authors: We acknowledge that the DM-only calibration omits baryonic effects present in TNG50. For the ultra-faint regime, however, the dominant process is tidal stripping by the cluster potential, and baryonic mass loss is sub-dominant once the galaxies are below ~10^7 M_sun. We will expand the Discussion to include a dedicated paragraph on this approximation, citing literature on baryon-tide interactions in dwarfs, and note that any residual bias would primarily affect the innermost, most-stripped objects without altering the overall conclusion that the population traces the dark-matter profile. revision: partial

  3. Referee: [Results] Results or Discussion (abundance and uncertainty): The quoted range of 2000-7000 systems with M_* > 100 M_sun is derived from the three most massive TNG50 clusters, but the manuscript provides no details on error propagation from the empirical model parameters, variation across the three clusters, or sensitivity to the minimum mass threshold. This weakens the robustness of the quantitative prediction.

    Authors: The quoted range already encodes cluster-to-cluster variation. In the revised manuscript we will tabulate the individual counts for each of the three clusters, add a sensitivity plot showing how the total number changes with the minimum stellar-mass threshold, and include a Monte-Carlo error band that propagates the principal uncertainties in the empirical model's mass-loss and disruption parameters. revision: yes

Circularity Check

0 steps flagged

No significant circularity; external calibration grounds the model

full rationale

The derivation applies an empirical tidal-evolution model (calibrated on separate idealized high-resolution N-body simulations) to sub-resolution galaxies identified in the independent TNG50 hydrodynamical runs. The reported abundances (2000-7000 systems with M*>100 Msun) and the radial-distribution claim both emerge from this forward application rather than from any fit to the same TNG50 data or from self-referential definitions. No equations or self-citations reduce the central results to tautologies; the calibration set is distinct from the target clusters, satisfying the criteria for non-circular external grounding.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The prediction rests on the transferability of a tidal-stripping model calibrated in idealized simulations to the cosmological TNG50 clusters and on the assumption that TNG50's large-scale tidal fields remain accurate below its resolution limit.

free parameters (1)
  • parameters of the empirical tidal evolution model
    The model is calibrated on high-resolution N-body simulations, implying fitted parameters that control mass loss rates under tides.
axioms (2)
  • domain assumption The tidal evolution model calibrated on idealized N-body simulations applies accurately to sub-resolution dwarfs inside TNG50 clusters.
    Invoked when the empirical model is combined with TNG50 output to reach below the simulation resolution.
  • domain assumption TNG50 correctly captures the large-scale gravitational tidal fields and cluster assembly history despite limited resolution.
    Required for the post-processing step to be valid.

pith-pipeline@v0.9.0 · 5571 in / 1590 out tokens · 77002 ms · 2026-05-08T10:43:30.536082+00:00 · methodology

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

1 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    J., Du X., 2022, MNRAS, 517, 1398 Boylan-KolchinM.,SpringelV.,WhiteS.D.M.,JenkinsA.,2010,MNRAS, 406, 896 Brough S., et al., 2020, arXiv e-prints, p

    Benitez-Llambay A., Frenk C., 2020, MNRAS, 498, 4887 Benson A. J., Du X., 2022, MNRAS, 517, 1398 Boylan-KolchinM.,SpringelV.,WhiteS.D.M.,JenkinsA.,2010,MNRAS, 406, 896 Brough S., et al., 2020, arXiv e-prints, p. arXiv:2001.11067 Davis M., Efstathiou G., Frenk C. S., White S. D. M., 1985, ApJ, 292, 371 Diemer B., Behroozi P., Mansfield P., 2024, MNRAS, 533...

    work page internal anchor Pith review doi:10.48550/arxiv.astro-ph/0307362 2020