Globular cluster distributions as a dynamical probe of dark matter

This paper utilizes N-body and semianalytic simulations to demonstrate that the orbital contraction of globular clusters due to dynamical friction serves as a robust, independent probe for confirming the presence of dark matter halos in ultradiffuse and dwarf galaxies like NGC5846-UDG1 and Fornax.

The Gist

Imagine a galaxy as a giant, invisible ocean of dark matter, and inside this ocean float heavy, glowing islands called Globular Clusters (GCs). These clusters are like massive ships sailing through the water.

This paper is about figuring out how thick and heavy that invisible ocean is by watching how the ships move.

The Problem: The Invisible Ocean

We know galaxies are held together by gravity, but most of that gravity comes from Dark Matter, which we can't see. Usually, astronomers try to measure this invisible stuff by watching how fast stars move (like watching cars on a highway to guess the road's width). But this paper uses a different trick.

It looks at Dynamical Friction. Think of this like a swimmer moving through a pool.

• If the pool is filled with thick honey (a lot of dark matter), the swimmer moves fast but gets slowed down quickly by the sticky fluid.
• If the pool is just thin air (no dark matter), the swimmer doesn't get slowed down much by the air, but they might crash into other things or sink faster due to their own weight.

In a galaxy, the "swimmer" is a Globular Cluster. As it moves through the dark matter "honey," it drags the dark matter behind it, creating a wake that pulls the cluster backward. This causes the cluster to lose energy and spiral inward toward the center of the galaxy.

The Detective Work

The authors, Nativ Ben-Yeda, Kfir Blum, and Inbar Havilio, acted like detectives trying to solve a mystery: How much dark matter is in these galaxies?

They picked three specific galaxies to investigate:

1. UDG1: A very faint, spread-out galaxy.
2. Fornax: A small dwarf galaxy near our own Milky Way.
3. DF44: Another very faint, spread-out galaxy.

They used supercomputers to run thousands of simulations. They asked: "If we start these clusters in different places and with different masses, where will they end up after 10 billion years?"

The Findings

1. The "Honey" vs. "Air" Test

• The "No Dark Matter" Scenario: If there were no dark matter (just the visible stars), the "honey" would be very thin. The heavy clusters would spiral inward very quickly, crashing into the center and forming a giant, dense ball of stars.
• The "Dark Matter" Scenario: If there is a lot of dark matter, the "honey" is thick. The clusters get slowed down gently, and they stay spread out over a wider area.

2. The Results for UDG1 and Fornax
When the authors compared their computer simulations to real telescope photos, they found a clear pattern:

• UDG1 and Fornax: The clusters in these galaxies are spread out exactly as if they were swimming through thick honey. If there were no dark matter, the clusters would have already crashed into the center. The fact that they are still spread out is strong proof that a massive, invisible dark matter halo is holding them back.
• The Surprise: This is a new way to prove dark matter exists. It doesn't rely on measuring the speed of stars (kinematics); it relies on the position of the clusters. It's like knowing a room is full of people not by hearing them talk, but by seeing how hard it is to walk through the crowd.

3. The Result for DF44

• DF44: This galaxy is so diffuse (spread out) that the "honey" is very thin, or the clusters are so far out that the friction is too weak to give a clear answer. The data here is a bit too fuzzy to say for sure if there is dark matter or not, though it doesn't rule it out.

The "What If" Scenarios

The authors were careful. They knew that maybe the clusters didn't start where the stars are now.

• The "Stretched" Start: What if the clusters started way further out and just drifted in? They tested this. Even if they started further out, the "No Dark Matter" models still predicted that the clusters would crash into the center too fast. The "Dark Matter" models were the only ones that matched the real data.
• The "Heavy" vs. "Light" Clusters: They also tested if the clusters were losing mass (getting lighter) over time. Even with different assumptions about how heavy the clusters are, the conclusion for UDG1 and Fornax remained the same: They need dark matter to explain why the clusters haven't crashed into the center yet.

The Bottom Line

This paper argues that Globular Clusters are excellent "probes" for dark matter.

• In UDG1 and Fornax, the clusters are acting like buoys in a thick ocean. They haven't sunk to the bottom because the ocean (dark matter) is heavy and thick.
• This confirms that these galaxies are dominated by dark matter, using a method that is completely different from the usual way astronomers measure it.
• It suggests that the standard theory of cold dark matter works perfectly well to explain these galaxies, without needing any "exotic" or strange new physics.

In short: The clusters are still floating where they should be, proving that the invisible ocean of dark matter is real and heavy.

Technical Summary

Technical Summary: Globular Cluster Distributions as a Dynamical Probe of Dark Matter

Problem Statement
Globular clusters (GCs) act as massive probe particles traversing the dark matter (DM) halos of their host galaxies. Over time, gravitational dynamical friction (DF) arising from scattering between GCs and halo particles causes GC orbits to contract. This effect provides a "beyond-mean-field" test of the cold dark matter paradigm, distinct from traditional stellar kinematics which rely on the mean gravitational field. While previous studies have suggested that GC mass segregation in dwarf and ultradiffuse galaxies (UDGs) could constrain DM properties, the interpretation of such data is complicated by significant observational and theoretical uncertainties. Key challenges include robust GC identification, mass estimation, and, most critically, the unknown initial conditions of the GC radial distribution and the initial mass function (GCIMF). Without a systematic survey of plausible initial conditions, it is difficult to distinguish between a DM-induced signal and alternative dynamical histories.

Methodology
The authors employ a mixed computational approach to survey the space of initial conditions and DM halo models:

1. Semianalytic Simulations: The primary workhorse involves orbit integration using the Chandrasekhar formula for DF, treating the halo potential as a mean field but modeling the GC system as an N-body problem to include GC-GC interactions. This method allows for the simulation of thousands of systems.
2. N-body Simulations: To validate the semianalytic approach, the authors use a modified Barnes-Hut tree code ("GONBY") with a "live" halo. In these simulations, DM and stellar particles are initialized using the Eddington formalism to ensure statistical ergodicity. GCs are treated as point particles that can merge.
3. Physical Modeling:
   • Halo Profiles: The study compares Burkert (core) and Navarro-Frenk-White (NFW, cusp) DM profiles against a DM-free stellar-only model.
   • GC Evolution: The models incorporate mass loss via stellar evolution, two-body evaporation, and tidal stripping. GC mergers are modeled based on energy and proximity criteria.
   • Initial Conditions: The authors systematically vary the initial radial distribution of GCs (ranging from matching the current stellar distribution to being "stretched" by factors of 2 or 3) and the maximum mass in the GCIMF ($M_{max}$).
4. Observational Targets: The analysis focuses on three specific systems: NGC5846-UDG1 (UDG1), the Fornax dwarf spheroidal (dSph), and the Coma cluster UDG-DF44.

Key Results

• NGC5846-UDG1 (UDG1): The GC luminosity cumulative distribution function (CDF) strongly favors models containing a substantial DM halo. DM-free models predict rapid orbital contraction and excessive mass segregation that contradicts the observed gentle trend. While the specific DM halo mass required depends on the initial radial stretch of the GCs, a DM halo is robustly required to explain the data within reasonable limits. The results are consistent with, though slightly weaker than, constraints derived from stellar kinematics.
• Fornax dSph: Despite having only six GCs, the system provides a strong indication of a DM halo. Naive estimates suggest that without a massive DM halo, Fornax's GCs would have spiraled into the center long ago. The authors demonstrate that standard DM-induced DF can explain the current GC morphology without invoking exotic DM physics, provided the initial GC distribution was somewhat stretched relative to the current stellar body. Both core and cusp halo models are consistent with the data under these conditions.
• UDG-DF44: This galaxy is too diffuse for DF to provide strong constraints. The stellar radius is significantly larger than that of UDG1, and reasonable initial conditions (including significant radial stretch) can match the observed GC distribution even in DM-free models. Consequently, the GC-based constraints for DF44 are much weaker than those from kinematics.
• Model Robustness: The authors find that while the specific DM mass constraints depend on initial conditions (particularly the radial stretch and mass loss rates), the qualitative conclusion—that UDG1 and Fornax require DM halos to prevent excessive mass segregation—remains robust across a wide parameter space.

Significance and Claims
The paper claims to provide a "beyond-mean-field" positive signature of dark matter. Unlike kinematic analyses that infer DM from the mean gravitational potential, this work utilizes the non-linear, scattering-driven process of dynamical friction. The authors argue that the absence of extreme mass segregation in UDG1 and Fornax is a "null signal" indicating the presence of a high-velocity dispersion DM halo, which suppresses DF efficiency.

The study emphasizes that while GC morphology cannot yet replace stellar kinematics as a precise quantitative tracer of DM due to systematic uncertainties (initial conditions, mass loss, mergers), it offers a complementary and independent test. The results suggest that exotic DM properties are not required to explain the GC data in Fornax or UDG1; standard cold dark matter models with induced dynamical friction are sufficient, provided the initial conditions of the GC systems are accounted for. The authors explicitly state they do not claim to definitively discriminate between core and cusp profiles for these galaxies given current uncertainties, but they establish the viability of using GC distributions as a dynamical probe for DM halos in dwarf and ultradiffuse galaxies.