Influence of tides and self-gravity on Ultra-Light Dark Matter Bounds from Dwarf Galaxies

This paper investigates how tidal interactions with the Milky Way and stellar self-gravity affect ultra-light dark matter constraints derived from dwarf galaxies, finding that despite these systematics, dark matter masses between $5\times 10^{-22}$ and $5\times 10^{-21}$ eV remain in tension with current observational data.

The Gist

Imagine the universe is filled with a mysterious, invisible fog called Dark Matter. For decades, scientists have been trying to figure out what this fog is made of. One popular idea is that it's made of "Ultra-Light Dark Matter" (ULDM)—particles so incredibly light and wavy that they act more like a giant, shimmering ocean than a collection of solid rocks.

This paper is like a team of detectives (Andrea Caputo and Luca Teodori) re-investigating a cold case: Can this wavy fog exist in the tiny, lonely galaxies orbiting our own Milky Way?

Here is the story of their investigation, broken down into simple concepts.

1. The Crime Scene: The "Dwarf" Galaxies

Our Milky Way is a massive city of stars. Orbiting it are tiny, dim "dwarf galaxies" (like Fornax, Carina, and Leo II). These are the perfect crime scenes because they are small and quiet.

• The Theory: If ULDM exists, it creates ripples and waves in the fog. These waves should act like a gentle but constant shaking, causing the stars in these dwarf galaxies to heat up and spread out over time (like popcorn kernels popping and scattering).
• The Problem: In a previous study, the detectives found that for certain weights of these ULDM particles, the stars should have spread out much more than they actually have. The galaxies looked too "tight" and calm. This suggested that ULDM of that specific weight might not exist.

2. The Suspects: Two New Alibis

The paper's authors decided to check if there were any "loopholes" that could explain why the galaxies looked so calm, even if the wavy fog was there. They looked at two main suspects:

Suspect A: The "Tidal Stripper" (The Milky Way's Gravity)

Imagine a small boat (the dwarf galaxy) sailing near a giant, turbulent ocean liner (the Milky Way). The giant ship's wake creates huge waves that can rip parts of the small boat off.

• The Analogy: As these dwarf galaxies orbit the Milky Way, the Milky Way's gravity pulls on the outer edges of the dark matter fog. This "tidal stripping" tears away the outer, wavy parts of the fog.
• The Effect: If the outer waves are ripped away, the remaining fog is calmer. It stops shaking the stars as violently. The authors asked: Did the Milky Way strip the fog so much that it stopped heating the stars, making the ULDM theory look innocent again?
• The Verdict: They simulated the orbits of these galaxies for billions of years. Even with the "ripping" effect, the fog was still too wavy. The stars should still be shaking more than they are.

Suspect B: The "Compact Crowd" (Stellar Self-Gravity)

Imagine a crowd of people in a room. If the room is empty, a gentle breeze (the ULDM waves) can easily push people around. But if the people are huddled together in a tight, dense group, they hold onto each other, and the breeze can't move them as easily.

• The Analogy: The authors considered that in the past, the stars in these dwarf galaxies might have been packed much tighter than they are today. If the stars were a "dense crowd," their own gravity would hold them together, making them resistant to the shaking caused by the ULDM waves.
• The Verdict: They ran simulations where the stars were packed very tightly. While this did slow down the shaking a bit, it wasn't enough to save the ULDM theory. The stars still shouldn't be as calm as they are.

3. The Final Ruling: The Fog is Still Too Wavy

After running thousands of computer simulations—accounting for the Milky Way ripping the fog apart and the stars huddling together tightly—the detectives reached a conclusion:

The "loopholes" don't work.

Even with the best possible excuses (strong tides and dense stars), the Ultra-Light Dark Matter particles with masses between $5 \times 10^{-22}$ and $5 \times 10^{-21}$ eV simply don't fit the evidence. If these particles existed, the dwarf galaxies would look much more chaotic and spread out than they do.

The Takeaway

Think of it like trying to fit a square peg into a round hole. The "square peg" is the theory of Ultra-Light Dark Matter with those specific masses. The "round hole" is the actual data we see from the dwarf galaxies.

The authors tried to sand down the square peg (by adding tidal stripping and stellar gravity) to make it fit, but it still didn't work. The hole is just too round.

In short: The universe might still be full of dark matter, but if it is this specific type of "ultra-light, wavy" stuff, it's likely not the kind that weighs between those two tiny numbers. The search for the true nature of dark matter must continue, perhaps looking for heavier particles or different types of fog entirely.

Technical Summary

1. Problem Statement

Ultra-Light Dark Matter (ULDM), also known as Fuzzy Dark Matter (FDM), is a candidate for dark matter with particle masses $m \sim 10^{-22} - 10^{-21}$ eV. Its defining feature is a large de Broglie wavelength that induces quantum interference effects, leading to dynamical heating of stellar systems. This heating causes the expansion of stellar distributions (increasing the half-light radius) and alters velocity dispersion profiles.

Previous work (specifically Ref. [38] by the same authors) found that ULDM in the mass range $5 \times 10^{-22} \lesssim m/\text{eV} \lesssim 5 \times 10^{-21}$ is in tension with observational data from dwarf spheroidal galaxies (dSphs) like Fornax, Carina, and Leo II. However, this tension relied on simulations that neglected two critical physical effects which could potentially suppress dynamical heating and reconcile the model with data:

1. Stellar Self-Gravity: The gravitational potential generated by the stars themselves, which could stabilize the system against ULDM fluctuations, especially if the stellar component was more compact in the past.
2. Tidal Stripping: The removal of the outer halo of the dwarf galaxy due to the Milky Way's gravitational potential. Since ULDM heating arises from interference between high-energy eigenmodes (which reside in the outer halo), stripping these modes could significantly reduce heating.

This paper aims to rigorously test whether including these two effects invalidates the previous exclusion of the ULDM mass range.

2. Methodology

The authors performed numerical simulations using a 3D Schrödinger-Poisson (SP) solver coupled with $N$-body stellar dynamics.

A. Simulation Setup

• ULDM Evolution: The SP solver evolves the ULDM field $\psi$ using a pseudo-spectral method. The field is treated as a superposition of energy eigenfunctions.
• Stellar Dynamics: Stars are modeled as massive particles evolving in the total gravitational potential $\Phi_{\text{tot}} = \Phi_{\text{ULDM}} + \Phi_{\text{stars}}$. The code explicitly includes stellar self-gravity (the potential generated by the stars acting back on themselves and the ULDM field), correcting a limitation of previous test-particle approximations.
• Target Systems: Simulations were run for five dwarf galaxies: Fornax, Carina, Leo II, Leo I, and Draco.
• Initial Conditions:
  • No Tides: Standard soliton + NFW profiles.
  • With Tides: To model tidal stripping, the authors used an eigenfunction truncation method (following Refs. [44, 46]). They decomposed the ULDM halo into eigenmodes and imposed a hard cutoff on the maximum radial ($n_{\text{max}}$) and angular ($l_{\text{max}}$) quantum numbers. This effectively removes high-energy, spatially extended modes, simulating a halo truncated at a specific tidal radius $r_t$.

B. Tidal History Reconstruction

To determine the appropriate tidal radius $r_t$ for each galaxy, the authors:

1. Reconstructed orbital histories for the dwarfs over the last 10 Gyr using the McMillan17 Milky Way potential (and cross-checked with MWPotential2014).
2. Propagated observational uncertainties (distance, proper motion, line-of-sight velocity) via Monte Carlo sampling ($N_{\text{MC}}$ realizations).
3. Calculated the instantaneous tidal radius $r_t(t)$ along each orbit using the full tidal tensor of the host potential.
4. Adopted a conservative estimate for the initial conditions: the 2.5th percentile of the minimum tidal radius-to-half-light ratio ($r_t/r_{1/2}$) encountered over the orbit. This ensures the simulations assume the maximum plausible stripping consistent with data.

C. Stellar Self-Gravity Scenarios

The authors tested scenarios where the stellar component was significantly more compact in the past (smaller initial half-light radius) or more massive than currently inferred, to see if this could suppress heating enough to save the ULDM model.

3. Key Contributions

• Explicit Self-Gravity Implementation: Unlike previous studies that treated stars as massless test particles, this work fully couples stellar self-gravity to the ULDM field, allowing for a realistic assessment of how a dense stellar core resists heating.
• Conservative Tidal Modeling: The authors developed a framework to reconstruct tidal histories and implement stripping via eigenmode truncation. They argue this "hard cutoff" approach is conservative (i.e., it likely overestimates the suppression of heating compared to realistic, time-dependent stripping) because it removes all power above the cutoff immediately, whereas real stripping might retain some residual power.
• Multi-Galaxy Analysis: The study extends the analysis beyond the original three galaxies (Fornax, Carina, Leo II) to include Leo I and Draco, providing a broader statistical basis for constraints.

4. Key Results

A. Stellar Self-Gravity

• Effect: Including stellar self-gravity does slow down the growth of the half-light radius ($r_{\text{half}}$) compared to simulations without it. A more compact or massive stellar core provides a deeper potential well, making the system more "adiabatic" against ULDM fluctuations.
• Magnitude: However, the suppression is insufficient to resolve the tension. Even in extreme scenarios where the stellar mass was 5 times larger than current estimates (or the initial radius was 5 times smaller), the dynamical heating for $m = 5 \times 10^{-21}$ eV remained too high to match observations for Carina and Leo II.
• Conclusion: Stellar self-gravity cannot save the ULDM model in the $10^{-21}$ eV range for these systems.

B. Tidal Stripping

• Effect: Tidal stripping significantly reduces dynamical heating by removing the outer halo modes responsible for the interference fluctuations.
• Fornax: For $m = 10^{-21}$ eV and $5 \times 10^{-22}$ eV, tidal stripping slows the heating of the half-light radius. However, the constraint for Fornax is driven by the inner velocity dispersion peak caused by the growing soliton core, not just heating. This peak persists even in tidally stripped halos, keeping the model in tension with data.
• Carina/Leo II: For $m = 5 \times 10^{-21}$ eV, the tension is driven by heating. While extreme tidal stripping (e.g., $r_t \sim 0.5$ kpc for Leo II) could theoretically make the data compatible with the model, such small tidal radii are statistically unlikely (lying in the extreme tail of the Monte Carlo distribution) and are inconsistent with the lack of observed tidal features in these galaxies today.
• Leo I/Draco: These systems also show strong tension due to both heating and soliton-induced kinematic peaks, even with tidal stripping included.

C. Final Constraints

The study confirms that the mass range $5 \times 10^{-22} \lesssim m/\text{eV} \lesssim 5 \times 10^{-21}$ remains in tension with current dwarf galaxy data, even when accounting for:

1. Conservative estimates of tidal stripping.
2. Explicit stellar self-gravity.
3. Variations in initial stellar density and mass.

5. Significance and Conclusion

• Robustness of Exclusion: The exclusion of the ULDM mass range $5 \times 10^{-22} - 5 \times 10^{-21}$ eV is robust against the two most significant systematic uncertainties identified in the field: tidal stripping and stellar self-gravity.
• Methodological Benchmark: The paper establishes a framework for including stellar self-gravity and tidal effects in ULDM simulations. The authors argue that while their tidal implementation is a simplified "hard cutoff," it provides a conservative lower bound on heating, making the resulting constraints reliable.
• Future Directions: The authors note that future work should focus on ultra-faint dwarf galaxies (like Segue I) to probe higher masses ($10^{-20} - 10^{-19}$ eV) and on more sophisticated cosmological zoom-in simulations that evolve the halo and stars self-consistently within the Milky Way potential over cosmic time.

In summary, this work strengthens the case against Ultra-Light Dark Matter in the $10^{-21}$ eV regime by demonstrating that realistic astrophysical systematics cannot reconcile the model's predictions with the observed stability and kinematics of local dwarf galaxies.