pith. sign in

arxiv: 2604.09539 · v1 · submitted 2026-04-10 · 🌌 astro-ph.GA

The missing ultra-faint satellites of the Milky Way

Isabel M.E. Santos-Santos , Carlos S. Frenk , Julio F. Navarro This is my paper

Pith reviewed 2026-05-10 16:39 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords Milky Way satellitesultra-faint dwarfsdark matter halostidal strippingcuspy profilesgalaxy formationsatellite mass function
0
0 comments X

The pith

A simulation of a Milky Way-like halo predicts many ultra-faint dwarf satellites remain undiscovered and share properties with recently found compact systems.

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

The paper combines a high-resolution N-body simulation of a Milky Way-mass halo with a semianalytic galaxy formation model to forecast the entire satellite population. It shows that once tidal stripping is accounted for, the predicted masses of the brighter satellites match observations. This leads to the conclusion that a large population of ultra-faint dwarfs must still be missing from current surveys. The missing satellites are expected to lie outside the innermost halo and to have densities and other traits that overlap with ultra-faint compact satellites, which the model treats as dark-matter-dominated objects that survived in dense subhalo cores.

Core claim

The central claim is that the total number of luminous satellites ever accreted by a Milky Way-like halo is limited to a few hundred down to arbitrarily low luminosity because galaxies form only in subhalos above the hydrogen-cooling threshold. After applying a treatment of dark matter and stellar tidal stripping that preserves all cuspy subhalos to the present day, the model reproduces the massive end of the observed satellite mass function while predicting a substantial population of missing ultra-faint dwarfs. These missing systems avoid the central regions and possess mean densities of order 10^10 to 10^11 solar masses per cubic kiloparsec, properties that align with those of ultra-faint

What carries the argument

The novel treatment of dark matter and stellar tidal stripping that assumes all dark matter subhalos survive until the present because of their cuspy inner density profiles, together with the tracking of orphan galaxies after their subhalos become unresolved.

If this is right

  • The current Milky Way satellite census is incomplete, with a large number of ultra-faint dwarfs still undetected.
  • Missing ultra-faint dwarfs avoid the innermost regions of the host halo.
  • Ultra-faint compact satellites are dark matter dominated with velocity dispersions of 1 to 3 km/s and mean densities of 10^10 to 10^11 solar masses per cubic kiloparsec.
  • Confirmation would support both the cuspy inner structure of dark matter halos and the hydrogen cooling threshold as the limit for galaxy formation.

Where Pith is reading between the lines

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

  • Kinematic data on ultra-faint compact satellites could distinguish whether they are truly dark-matter dominated or star clusters.
  • Future surveys sensitive to faint objects at moderate distances from the galactic center could directly test the spatial distribution prediction.
  • The high predicted densities imply that some objects currently classified near the globular-cluster boundary may instead be the densest surviving dwarfs.

Load-bearing premise

Galaxies form only in subhalos whose peak circular velocity exceeds the hydrogen cooling threshold and all dark matter subhalos survive to the present day because of their cuspy central density profiles.

What would settle it

A complete census of Milky Way satellites that finds either far fewer or far more ultra-faint systems with mean densities of 10^10 to 10^11 solar masses per cubic kiloparsec, or that finds such systems concentrated in the innermost halo regions, would contradict the predictions.

Figures

Figures reproduced from arXiv: 2604.09539 by Carlos S. Frenk, Isabel M.E. Santos-Santos, Julio F. Navarro.

Figure 1
Figure 1. Figure 1: Enclosed mass profile and circular velocity of the Aq-A-L1 dark matter halo, for the DMO and Disc potential models. Results at 𝑧 = 0 are shown in black and at 𝑧 = 2 in grey. precise value of the critical mass and its redshift dependence are still debated, this overall scheme is widely accepted at present. On the observational side, the advent of digital sky surveys has led to the discovery of new MW satell… view at source ↗
Figure 2
Figure 2. Figure 2: 𝑀∗,peak-𝑉peak relation predicted by GALFORM for Aq-A-L1 satel￾lites at infall time. A grey dashed line shows an extrapolation of the Moster et al. (2013) abundance-matching model. A green line indicates a power-law fit to the data (see SS25). kpc, a value roughly consistent with the region of the Milky Way where the circular velocity profile is observed to be reasonably flat (Lian et al. 2024). 2.2 The GAL… view at source ↗
Figure 3
Figure 3. Figure 3: Evolution of an example type-1 (“surviving”) Aq-A-L1 satellite. Left: Radial distance from the host halo (top panel), and 𝑉max, and 𝑟max as a function of time (bottom panel). The grey dashed line in the top panel shows the evolution of the host’s 𝑟200. The period of time before infall, when the galaxy is considered a central, is marked with circles. The dashed cyan line shows the orbit integrated backwards… view at source ↗
Figure 4
Figure 4. Figure 4: Subhalo characteristic radii and velocities (𝑟max and 𝑉max) of Aq-A-L1 satellites at 𝑧 = 0, according to the SUBFIND halo finder catalogue, versus their values as predicted by the tidal tracks of Errani & Navarro (2021). Left: 𝑉max vs 𝑟max relation in physical units. Grey symbols (circles for type-1s and stars for type-2s) show the 𝑧 = 0 values as measured by SUBFIND. In the case of type-2s, values corresp… view at source ↗
Figure 5
Figure 5. Figure 5: Left and Middle panels: Pericentric vs apocentric distances computed for Aq-A-L1 satellites using the DMO (left) or Disc (middle) potentials. Points are colored by the fraction of dark matter mass bound to the subhalo by 𝑧 = 0, after stripping, computed using the tidal track model described in Sec.2.4. Right: Points show the mean dark matter density within 𝑟peak of Aq-A-L1 satellites at infall time, versus… view at source ↗
Figure 6
Figure 6. Figure 6: Stellar tidal tracks from Errani et al. (2022) used in this work. Lines of constant color indicate the track of a model with fixed initial radial segregation values, where segregation is defined as 𝑆 = 𝑅half/𝑟peak. The plot shows the change in stellar mass (top) and change in size (bottom), relative to their peak values at infall time, as a function of remnant subhalo mass fraction, log 𝑀max/𝑀peak (see Sec… view at source ↗
Figure 7
Figure 7. Figure 7: 𝑀∗,peak-𝑅half,peak relations assumed in this work for Aq-A-L1 satellites. 𝑀∗,peak values are predicted by GALFORM ( [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: illustrates the tidal evolution of the stellar component of the same satellite chosen to illustrate the dark matter evolution in [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Left: Satellite stellar mass function in Aq-A-L1. The black line shows the total (type-1s and type-2s/orphans) “unabridged” stellar mass function of Aq-A-L1 satellites, at infall time (𝑀∗,peak, see fig.4 in SS25). Coloured lines show the predicted 𝑧 = 0 functions after accounting for stellar tidal stripping (see details in Secs. 2.4 and 2.5). Cyan lines correspond to results assuming the "Scaled" luminosit… view at source ↗
Figure 10
Figure 10. Figure 10: The predicted 𝑧 = 0 cumulative 𝑉-band luminosity function of MW satellites. The red (cyan) line shows the prediction assuming the "Fiducial" ("Scaled") luminosity-size prescription for our stellar stripping model. We assume the "Disc" host potential. Absolute magnitudes 𝑀𝑉 have been estimated from stellar masses assuming a fixed mass-to-light ratio of 𝑀∗/𝐿𝑉 = 2. Grey lines with symbols show currently know… view at source ↗
Figure 11
Figure 11. Figure 11: Predicted 𝑧 = 0 size-luminosity relation for MW satellites. Size is expressed as the projected half-light stellar radii of satellite galaxies, 𝑅half. The left (right) panel shows predictions from our stellar stripping model assuming the "Scaled" ("Fiducial") initial luminosity-size prescription. We adopt the "Disc" host potential. Points are colored with different shades according to their 𝑀𝑉 . Observatio… view at source ↗
Figure 12
Figure 12. Figure 12: shows the predicted radial distribution of satellites, shown as a scatter plot of absolute magnitude, 𝑀𝑉 , as a function of Galac￾tocentric distance, 𝐷MW, at 𝑧 = 0. As in [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Left: Predicted 𝑧 = 0 mean density within the half-mass radius, 𝜌¯(< 𝑟half ) vs 𝑟half for MW satellites. As in previous figures, predictions from Aq-A-L1 are shown in shades of red while observational data are shown in in grey. Densities for Aq-A-L1 satellites are computed as described in Sec. 3.5. The "Fiducial" luminosity-size model with Disc potential is assumed to estimate stellar stripping. For obser… view at source ↗
Figure 14
Figure 14. Figure 14: Predicted 𝑧 = 0 velocity dispersion (𝜎) vs projected half-light radius (𝑅half) for MW satellites. The Fiducial Disc luminosity-size prescription is assumed. As in previous figures, predictions from Aq-A-L1 are shown in red and observational data in grey (see App. A). Small black dots mark the model galaxies that overlap with observed UFCS systems in luminosity and size (see box in [PITH_FULL_IMAGE:figure… view at source ↗
read the original abstract

We combine the highest resolution N-body simulation of a $\sim 10^{12}\, M_\odot$ $\Lambda$CDM halo (Aquarius-A) with the {\sc GALFORM} galaxy formation semianalytic model to study the full satellite population expected in a MW-like galaxy. The model assumes that galaxies only form in subhalos whose peak circular velocity exceeds the H-cooling threshold, all of which are well resolved in the simulation. The number of luminous subhalos ever accreted into the main halo is thus well defined, and implies that the total number of MW satellites, down to arbitrarily low luminosity, should not exceed a few hundred. The model tracks satellites even after their halos cease to be resolved ("orphan" galaxies), and includes a novel treatment of dark matter and stellar tidal stripping which takes into account that all $\Lambda$CDM subhalos survive until the present because of their cuspy inner density profiles. After accounting for tides, our results match well the massive end of the observed MW satellite mass function and predict that a large number of ultra-faint dwarfs are missing from the current MW satellite census. The missing UFDs are predicted to avoid the innermost regions of the host, and to have properties that overlap with those of the many ultra-faint compact MW satellites (UFCSs) discovered recently, with properties intermediate between globular clusters and dwarf galaxies. Our results suggest that many UFCS systems are dark matter-dominated dwarfs with velocity dispersions between $1-3$km/s, which have survived disruption because they reside in the dense cusp of $\Lambda$CDM subhalos. UFCSs should have mean densities of order $10^{10}$-$10^{11}\,M_\odot/$kpc$^3$, higher than those of more extended ultra-faint systems. If confirmed, our results would provide support for the cuspy nature of $\Lambda$CDM dark matter halos and for the hydrogen-cooling threshold for galaxy formation.

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 / 3 minor

Summary. The manuscript combines the high-resolution Aquarius-A N-body simulation of a ~10^12 M_⊙ ΛCDM halo with the GALFORM semianalytic galaxy formation model. Galaxies are assumed to form only in subhalos whose peak circular velocity exceeds the hydrogen-cooling threshold; the model tracks 'orphan' galaxies after subhalos drop below resolution and applies a novel prescription for dark-matter and stellar tidal stripping that assumes all cuspy ΛCDM subhalos survive to z=0. After accounting for tides, the results match the massive end of the observed Milky Way satellite mass function and predict a large population of missing ultra-faint dwarfs whose properties overlap with recently discovered ultra-faint compact satellites (UFCSs), with mean densities 10^{10}–10^{11} M_⊙ kpc^{-3}. The work concludes that many UFCSs are dark-matter-dominated systems with velocity dispersions 1–3 km s^{-1} and that the findings support cuspy CDM profiles and the H-cooling threshold.

Significance. If the survival assumption and stripping prescription hold, the paper supplies a concrete, observationally testable prediction for the number, radial distribution, and internal densities of undetected ultra-faint satellites, potentially identifying many UFCSs as the missing population. The combination of a well-resolved simulation with a calibrated semianalytic model that explicitly follows orphans to arbitrarily low luminosity is a clear strength, as is the derivation of falsifiable density and velocity-dispersion ranges for the predicted objects. The result, if confirmed, would strengthen the case for cuspy inner profiles in low-mass halos within ΛCDM.

major comments (3)
  1. [§3] §3 (model description) and the novel tidal-stripping prescription: the central claim that a large number of UFDs survive with densities 10^{10}–10^{11} M_⊙ kpc^{-3} and overlap with UFCS properties rests on the assumption that all cuspy ΛCDM subhalos persist to z=0; no quantitative comparison is presented to hydrodynamical simulations that include baryonic feedback capable of softening cusps or enhancing disruption via disk shocking. This directly determines the predicted survival fraction at the lowest masses.
  2. [Model assumptions] Model assumptions (H-cooling threshold): the total number of luminous satellites is fixed by the chosen peak-circular-velocity threshold for hydrogen cooling, a free parameter calibrated inside GALFORM to reproduce observed galaxy populations; the statement that the total number 'should not exceed a few hundred' therefore follows from the input calibration rather than constituting an independent prediction from the simulation alone.
  3. [Results section] Comparison to observations (results section): while a match to the massive end of the MW satellite mass function is asserted, no quantitative metrics (e.g., number counts above a stated mass threshold, Kolmogorov–Smirnov statistic, or tabulated residuals) are supplied, making it difficult to judge how tightly the model is constrained before the extrapolation to the ultra-faint regime is performed.
minor comments (3)
  1. The mean-density range quoted for UFCSs (10^{10}–10^{11} M_⊙ kpc^{-3}) would benefit from an explicit formula or reference to the exact definition used in the stripping calculation.
  2. A supplementary table listing predicted numbers of missing UFDs in radial and luminosity bins would improve clarity for future observational tests.
  3. A few references to the most recent UFCS discoveries appear slightly dated; updating the citation list would strengthen the discussion of observational overlap.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their insightful comments on our manuscript. We address each of the major comments below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: §3 (model description) and the novel tidal-stripping prescription: the central claim that a large number of UFDs survive with densities 10^{10}–10^{11} M_⊙ kpc^{-3} and overlap with UFCS properties rests on the assumption that all cuspy ΛCDM subhalos persist to z=0; no quantitative comparison is presented to hydrodynamical simulations that include baryonic feedback capable of softening cusps or enhancing disruption via disk shocking. This directly determines the predicted survival fraction at the lowest masses.

    Authors: We agree that our results depend on the survival of cuspy subhalos to z=0, as assumed in the Aquarius simulation and our stripping prescription. A direct quantitative comparison to hydrodynamical simulations is not included in this work, as our approach combines N-body with semianalytic modeling. However, we will add a paragraph in the revised manuscript discussing this assumption, citing hydrodynamical studies on subhalo disruption and baryonic effects, to provide better context for readers. revision: partial

  2. Referee: Model assumptions (H-cooling threshold): the total number of luminous satellites is fixed by the chosen peak-circular-velocity threshold for hydrogen cooling, a free parameter calibrated inside GALFORM to reproduce observed galaxy populations; the statement that the total number 'should not exceed a few hundred' therefore follows from the input calibration rather than constituting an independent prediction from the simulation alone.

    Authors: The H-cooling threshold is calibrated in GALFORM, but the Aquarius simulation supplies the detailed accretion history and the exact number of subhalos that exceed this threshold. This makes the prediction of no more than a few hundred satellites a specific outcome of applying the calibrated model to this high-resolution simulation. We will revise the text to clarify this point and highlight the role of the simulation in determining the count. revision: yes

  3. Referee: Comparison to observations (results section): while a match to the massive end of the MW satellite mass function is asserted, no quantitative metrics (e.g., number counts above a stated mass threshold, Kolmogorov–Smirnov statistic, or tabulated residuals) are supplied, making it difficult to judge how tightly the model is constrained before the extrapolation to the ultra-faint regime is performed.

    Authors: We accept this criticism and will strengthen the results section by adding quantitative metrics. Specifically, we will report the number of model satellites above certain mass thresholds and include a Kolmogorov-Smirnov test to assess the goodness of fit to the observed massive satellite mass function. revision: yes

Circularity Check

0 steps flagged

No significant circularity; model predictions are independent computations from stated inputs

full rationale

The paper's central results follow from applying the Aquarius N-body simulation (which resolves all subhalos above the H-cooling threshold) plus the GALFORM model with its explicit assumptions on galaxy formation and a novel tidal-stripping prescription. The total satellite count is obtained by direct counting of resolved subhalos in the simulation that satisfy the threshold; this count is not set by fitting to the MW satellite census. After applying the stripping model, the massive-end mass function is shown to match observations (a validation step), and the number, locations, and densities of fainter systems are computed outputs. No equation or result is shown to equal its input by construction, and the survival assumption is stated openly rather than smuggled in via self-citation. The derivation therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 1 invented entities

The central claim rests on the hydrogen cooling threshold as the sole limit on galaxy formation, the assumption that all subhalos survive due to cusps, and the GALFORM model parameters; no new entities are postulated beyond standard orphan galaxies.

free parameters (1)
  • H-cooling threshold (peak circular velocity)
    The minimum halo size required for galaxy formation; its specific value determines the total number of luminous satellites.
axioms (2)
  • domain assumption LambdaCDM cosmology with cuspy inner density profiles in all subhalos
    Invoked to argue that subhalos survive tidal disruption until the present day.
  • domain assumption GALFORM semianalytic model accurately captures galaxy formation physics above the cooling threshold
    The model is used without re-derivation to assign luminosities and track stripping.
invented entities (1)
  • orphan galaxies no independent evidence
    purpose: Track satellites whose dark matter halos are no longer resolved in the simulation
    Standard bookkeeping device in semianalytic models; no independent evidence provided here.

pith-pipeline@v0.9.0 · 5677 in / 1836 out tokens · 54585 ms · 2026-05-10T16:39:09.445986+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

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

  • IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    The model assumes that galaxies only form in subhalos whose peak circular velocity exceeds the H-cooling threshold... novel treatment of dark matter and stellar tidal stripping which takes into account that all ΛCDM subhalos survive until the present because of their cuspy inner density profiles.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    After accounting for tides, our results match well the massive end of the observed MW satellite mass function and predict... UFCSs should have mean densities of order 10^10-10^11 M_⊙/kpc^3.

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.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. The Hubble Missing Globular Clusters Survey IV. Ultra-faint compact satellites of the Milky Way. The case of Koposov 2

    astro-ph.GA 2026-06 unverdicted novelty 6.0

    Koposov 2 is shown to be an old (13.7 Gyr) star cluster with half-light radius 2.7 pc, absolute magnitude -0.95, and stellar mass 372 solar masses, supporting a star cluster classification over a dwarf galaxy.

  2. Discovery and Spectroscopic Characterization of a Distant, Compact Milky Way Satellite in Gemini

    astro-ph.GA 2026-06 unverdicted novelty 6.0

    Discovery and spectroscopic confirmation of DELVE 8/Gemini I as an ultra-faint compact Milky Way satellite at 120 kpc with low metallicity upper limit.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages · cited by 2 Pith papers

  1. [1]

    B., Garling, C

    Abazajian K. N., et al., 2009, ApJS, 182, 543 Abbott T. M. C., et al., 2018, ApJS, 239, 18 Ahvazi N., et al., 2025, arXiv e-prints, p. arXiv:2511.15808 Aihara H., et al., 2018, PASJ, 70, S4 Allen C. W., 1973, Astrophysical quantities Applebaum E., Brooks A. M., Christensen C. R., Munshi F., Quinn T. R., Shen S., Tremmel M., 2021, ApJ, 906, 96 Baumgardt H....

    work page arXiv 2009

  2. [2]

    Local Group Field

    and assuming a solar 𝑉-band absolute magnitude of𝑀𝑉 ,⊙ =4.83(Allen 1973). 3D half- mass radii𝑟half are estimated from observed 2D projected half-light radii𝑅 half by multiplying the latter by a factor of4/3to correct for spherical projection effects. Central mean densities are derived from the observed velocity dispersions using the dynamical mass estimat...

    work page 1973