🔭 astrophysics

Non-linear Evolution of Dark Plasma Subhalos

Using particle-in-cell simulations, this study demonstrates that non-linear electrostatic plasma instabilities in self-interacting dark matter subhalos orbiting the Milky Way can cause substantial turbulent heating and mass loss, potentially suppressing the low-mass subhalo population by up to 97% depending on orbital eccentricity.

Original authors: Andrew Liu, Anirudh Prabhu, Akaxia Cruz, Mariangela Lisanti

Published 2026-07-31

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Andrew Liu, Anirudh Prabhu, Akaxia Cruz, Mariangela Lisanti

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Invisible Storm in the Dark

Imagine the universe is filled with a mysterious, invisible substance called dark matter. We know it's there because its gravity acts like a cosmic glue, holding galaxies together and shaping the vast web of the cosmos. For decades, scientists thought of this dark matter as a cold, silent crowd of particles that only ever bumped into each other by accident, like ghosts passing through walls. But what if this crowd isn't so silent? What if, deep down, these particles can talk to each other, pushing and pulling with invisible forces?

This is the world of "dark plasma." Just like the air in a room can be filled with sound waves or the water in a pond can ripple when you throw a stone, a dark plasma can have its own internal storms. If dark matter particles carry a tiny electric charge (but a "dark" one), they can create waves and currents that heat them up, much like how a microwave heats food. This paper explores a wild idea: what happens when a small clump of this dark matter, called a subhalo, dives into the giant dark matter cloud surrounding our Milky Way galaxy? Does it stay intact, or does the invisible storm rip it apart?

The Cosmic Surfing Disaster

In this study, a team of physicists decided to play a high-stakes game of cosmic surfing. They imagined a small, dense clump of dark matter (a subhalo) zooming through the vast, diffuse ocean of dark matter that makes up the Milky Way. In the old, boring version of the story, this clump would just slowly lose a few particles as the galaxy's gravity tugged on it, like a snowball slowly melting as it rolls down a hill. This is called "tidal stripping," and it's a well-known process.

But the authors of this paper asked: What if the dark matter isn't just a passive snowball? What if it's a charged plasma that can get excited? They used powerful computer simulations to watch what happens when this subhalo "surfs" through the Milky Way's dark ocean. They found that when the subhalo moves fast enough, it creates a massive traffic jam in the invisible dark matter sea. This jam triggers a chaotic, invisible storm called a "streaming instability."

Think of it like this: Imagine a quiet crowd of people (the Milky Way's dark matter) standing still. Suddenly, a fast-moving group of skaters (the subhalo) zooms through them. In a normal world, the skaters just push past. But in this "dark plasma" world, the skaters' speed creates a shockwave that makes the crowd start vibrating and swirling wildly. This vibration acts like a giant, invisible blender. It heats up the dark matter particles inside the subhalo, making them move so fast that they fly off into space, leaving the subhalo shrunken and broken.

The Results: A Mass Extinction for Small Clumps

The simulations revealed some dramatic results. The team found that this "plasma heating" is a much more effective destroyer than the slow, gentle pull of gravity.

The Big Difference: For a small subhalo with a mass of 10 million suns ( $10^7 M_\odot$ ) on a very stretched-out, oval-shaped orbit, the plasma storm stripped away about 97% of its mass just after it made its closest pass to the center of the galaxy. In contrast, if only gravity were at work (the old way), it would have lost only about 30%.
The Size Matters: The smaller the subhalo, the more vulnerable it is. A slightly larger clump of 1 billion suns ( $10^9 M_\odot$ ) lost about 84% of its mass under the same conditions, compared to the 30% it would lose from gravity alone.
The "Outside-In" Effect: The destruction doesn't happen all at once. It starts at the edges. As the subhalo speeds up on its way to the center of the galaxy, the invisible storm gets stronger, peeling away the outer layers first. Only when the subhalo is moving very fast does the storm get strong enough to eat into the core.

The authors note that this effect depends heavily on how "eccentric" (oval-shaped) the orbit is. If the subhalo is on a nice, circular path, it moves slower and stays in the calmer, outer regions, so the storm is weaker. But on a wild, oval orbit that dives deep into the galaxy, the storm is fierce.

Why This Changes the Map of the Galaxy

The most exciting part of this discovery is what it implies for the population of dark matter clumps in our galaxy. Because this plasma storm is so good at destroying small, cold clumps, it suggests that the Milky Way might have far fewer tiny subhalos than we thought.

If this theory is correct, the "subhalo mass function" (a fancy way of counting how many clumps of different sizes exist) would look different. Instead of having a huge number of tiny, invisible clumps, we might see a sharp drop-off, with very few small clumps surviving at all. This could explain why we don't see as many tiny satellite galaxies as some older theories predicted.

The authors are careful to say that these results come from computer simulations, not direct observations yet. They used a specific model where dark matter particles interact through a "dark photon" (a hidden version of light). They found that for certain strengths of this interaction, the plasma effect is the dominant force, completely reshaping how these cosmic structures evolve.

In short, this paper suggests that the dark universe might be a much more turbulent and violent place than we imagined. Instead of a quiet, cold collection of ghosts, the dark matter halo of our galaxy might be a churning sea of invisible storms that can wipe out entire worlds before they even have a chance to shine.

Technical Summary: Non-linear Evolution of Dark Plasma Subhalos

Problem Statement
This paper investigates the dynamical evolution of dark matter (DM) subhalos orbiting within the Milky Way (MW) potential under the assumption that DM constitutes a "dark plasma." In models where DM self-interacts via long-range forces mediated by a hidden photon (a dark $U(1)$ sector), the dynamics resemble those of Standard Model plasmas. Specifically, when a DM subhalo moves through the ambient MW dark plasma, the relative drift velocity between the two populations can exceed their internal velocity dispersions. This configuration is susceptible to plasma instabilities (e.g., two-stream and beam-plasma instabilities), which generate dark electrostatic and electromagnetic fields. These fields scatter particles, modify the phase-space distribution, and drive the system toward thermal equilibrium. The central problem addressed is the non-linear evolution of these instabilities and their impact on subhalo mass loss, specifically comparing plasma-induced heating against standard gravitational tidal stripping.

Methodology
The authors employ a hybrid semi-analytic and numerical approach to model subhalo evolution:

Semi-Analytic Framework:
- Initial Conditions: Subhalos are modeled as spherical, net-neutral plasmas with NFW density profiles and Maxwell-Boltzmann velocity distributions. Two benchmark masses are considered: $10^7 M_\odot$ and $10^9 M_\odot$ .
- Orbital Dynamics: Subhalo orbits are integrated using the galpy package within a MW potential (MWPotential2014). Orbits are initialized at the virial radius with varying eccentricities ( $\epsilon = 0.1$ to $0.9$).
- Mass Loss Coupling: The subhalo density profile is discretized into concentric radial shells. The evolution of each shell is tracked over time steps ( $\delta t = 7$ Myr) up to $1.2 \times t_{\text{peri}}$ (time of first pericenter passage).
- Mechanisms: Mass loss is calculated via two distinct channels:
  - Tidal Stripping: Modeled using a standard efficiency parameter ( $\alpha = 0.55$ ) applied to shells outside the instantaneous tidal radius.
  - Plasma Heating: The fraction of bound particles ( $f_{\text{bound}}$ ) in each shell is determined by interpolating results from a suite of Particle-in-Cell (PIC) simulations. Particles heated above the local escape velocity ( $v_{\text{esc}} \approx 4\sigma_{\text{SH}}$ ) are removed, and the remaining bound mass re-virializes.
Particle-in-Cell (PIC) Simulations:
- Code: The study utilizes the TRISTAN-MP code, adapted for dark plasmas via coordinate rescaling invariance of the Vlasov-Poisson equations.
- Setup: One-dimensional (1D) electrostatic simulations are performed with periodic boundary conditions. The setup models counter-streaming beams representing the subhalo and MW background.
- Parameters: Simulations vary the density ratio ( $\rho_{\text{SH}}/\rho_{\text{MW}} \in [10^{-3}, 10^2]$ ) and relative streaming speed ( $v_{\text{SH}}/\sigma_{\text{MW}}$ ). The benchmark charge-to-mass ratio is $q_\chi/m_\chi = 10^{-14} e/m_p$ .
- Dimensionality: While the primary results are 1D, the authors perform a limited 2D-2V (two spatial, two velocity dimensions) simulation to verify that heating efficiency increases in higher dimensions, suggesting the 1D results are conservative.

Key Contributions

First Non-Linear Study: This work presents the first study of the non-linear evolution of electrostatic instabilities in DM subhalos using fully kinetic PIC simulations.
Quantification of Plasma Heating: The authors quantify the efficiency of plasma-induced mass loss, demonstrating that it can significantly exceed tidal stripping, particularly for low-mass subhalos on eccentric orbits.
Semi-Analytic Model: A novel framework is developed that couples local plasma instability growth (derived from PIC data) with global orbital evolution to predict subhalo mass and density profile evolution.

Results

Mass Loss Efficiency: For subhalos on highly eccentric orbits ( $\epsilon = 0.9$ $ϵ = 0.9$ ), plasma heating is the dominant mass loss mechanism. By the time of first pericenter ( $t = 1.2 t_{\text{peri}}$ $t = 1.2 t_{peri}$ ):
- A $10^7 M_\odot$ subhalo loses $\sim 97\%$ of its initial mass.
- A $10^9 M_\odot$ subhalo loses $\sim 84\%$ of its initial mass.
- In contrast, tidal stripping alone removes only $\sim 30\%$ of the mass for both cases.
Radial Dependence ("Outside-In" Stripping): Mass loss proceeds from the outside in. Instabilities grow most efficiently where the density contrast between the subhalo and the MW is modest ( $\rho_{\text{SH}}/\rho_{\text{MW}} \sim 1$ ) and the relative velocity is high. Consequently, outer shells are heated and stripped before the inner core. For highly eccentric orbits, the core is eventually affected near pericenter; for circular orbits, mass loss is confined to the outer regions.
Mass Function Implications: The preferential stripping of low-mass subhalos (due to their lower internal velocity dispersions and escape velocities) suggests a suppression of the low-mass end of the Subhalo Mass Function (SHMF). The authors estimate that for the benchmark parameters, the SHMF could exhibit a shallower slope or a sharp cutoff below $M_{\text{SH}} \sim 10^7 M_\odot$ .
Orbital Eccentricity: The effect is highly sensitive to eccentricity. For circular orbits ( $\epsilon = 0.1$ ), plasma heating removes only $\sim 50\%$ ( $10^7 M_\odot$ ) and $\sim 30\%$ ( $10^9 M_\odot$ ) of the mass, comparable to or slightly exceeding tidal effects.

Significance and Claims
The paper claims that dark plasma physics offers a distinct mechanism for subhalo evolution that can leave observable signatures in the MW satellite population. Specifically, the authors argue that plasma-induced heating and stripping can:

Reshape the SHMF: By preferentially destroying low-mass, cold subhalos, dark plasma models may resolve or alter the "missing satellites" problem or the "too big to fail" problem, potentially leading to a suppressed low-mass slope in the observed subhalo population.
Provide Observational Constraints: The suppression of low-mass subhalos and the resulting modification of the SHMF provide a new avenue for constraining dark sector parameters (specifically the charge-to-mass ratio $q_\chi/m_\chi$ ) using data from surveys like ELVES, SPARC, SAGA, and Merian.
Highlight Non-Linear Dynamics: The study emphasizes that collective plasma effects, driven by streaming instabilities, can operate on timescales comparable to galactic dynamical times, leading to rapid evaporation that standard collisionless CDM simulations do not capture.

The authors maintain a modest tone regarding the scope of their findings, noting that their results are based on 1D simulations (though 2D tests suggest 1D is conservative) and that a full population-level prediction requires future work incorporating self-consistent gravity and plasma effects in higher dimensions.

The Invisible Storm in the Dark

The Cosmic Surfing Disaster

The Results: A Mass Extinction for Small Clumps

Why This Changes the Map of the Galaxy

More like this