🔭 astrophysics

Mass and energy cascade in self-gravitating collisionless dark matter flow

This paper utilizes Illustris and Virgo simulations to demonstrate that self-gravitating collisionless dark matter flow exhibits a fundamental mass and energy cascade from large to small scales, which drives a statistically steady state characterized by scale-independent rates, self-similar halo structures, and universal density profiles derived from Fokker-Planck equations describing random walks in mass and particle space.

Original authors: Zhijie Xu

Published 2026-07-30

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Zhijie Xu

Original paper licensed under CC BY 4.0 (https://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

Imagine the universe not as a static stage, but as a churning, invisible ocean. This isn't water, but dark matter, a mysterious substance that doesn't shine, reflect light, or bump into normal matter. It only interacts through gravity, acting like a ghostly scaffolding that holds galaxies together. For decades, scientists have tried to understand how this invisible stuff moves and organizes itself. They often compare it to turbulence, like the swirling eddies in a rushing river or the steam rising from a hot cup of coffee. In a normal river, big swirls break apart into smaller ones, passing energy down until friction (viscosity) turns that motion into heat. But dark matter is "collisionless," meaning its particles never actually crash into each other like billiard balls; they just glide past one another, pulled only by gravity. The big question has been: if there's no friction to stop them, how does this cosmic ocean settle into the shapes we see today? Why do dark matter clumps (called haloes) form in specific sizes, and why do they have the specific density patterns they do?

This paper, titled "Mass and energy cascade in self-gravitating collisionless dark matter flow," proposes a fascinating new way to look at this cosmic dance. The author, Zhijie Xu, suggests that even without collisions, dark matter behaves like a turbulent fluid, but with a twist. Instead of energy just getting lost to heat, the universe is constantly "dissipating" energy because space itself is expanding. Think of the universe as a giant, stretching rubber sheet. As it stretches, the dark matter particles lose energy, much like a spinning top slowing down as the floor beneath it moves away. The paper argues that this energy loss drives a cascade: a continuous flow of mass and energy moving between different scales.

On the largest scale, small clumps of dark matter merge to form bigger ones, passing their mass "up" the ladder in an inverse mass cascade. It's like a game of musical chairs where the music never stops, and the chairs keep getting bigger. This process creates a "random walk" for the haloes, hopping from one mass size to the next, which naturally explains the number of haloes we see at different sizes (the halo mass function).

On the smaller scale, inside a single dark matter halo, particles are constantly migrating. They drift outward due to the expansion of the universe but are pulled back inward by gravity. This tug-of-war creates an energy cascade. The paper finds that this flow establishes a "statistically steady state," a kind of cosmic balance where energy is constantly flowing from the center of the halo to the edges and then "lost" to the expanding background. This flow dictates the density of the halo, leading to a specific rule: the density drops off as you move away from the center following a -4/3 power law (meaning if you double the distance, the density drops by a specific, predictable amount).

The most surprising twist comes from the "rules" of this game. By analyzing how fast particles must move and how long they wait before jumping to a new spot, the paper suggests that dark matter particles must be incredibly heavy—far heavier than the standard "WIMP" (Weakly Interacting Massive Particle) candidates that many physicists have been hunting for. The calculations suggest these particles might be superheavy, perhaps related to right-handed neutrinos, with masses around 10¹² GeV (about 10⁻¹⁵ kg). While this is a bold suggestion based on simulations and theoretical models, it offers a fresh perspective on the invisible architecture of our universe, turning the mystery of dark matter into a story of cascading flows and cosmic random walks.

Technical Summary: Mass and Energy Cascade in Self-Gravitating Collisionless Dark Matter Flow

Problem Statement
The paper addresses the fundamental physical principles governing the motion and evolution of self-gravitating collisionless dark matter (SG-CFD), a non-equilibrium system responsible for large-scale structure formation. While the abundance of haloes is typically described by the halo mass function (e.g., Press-Schechter formalism) and internal structures by density profiles (e.g., NFW, Einasto), the exact origins of these universal distributions and the variations in inner density slopes remain not fully understood. The author posits that these phenomena can be better understood by drawing an analogy to turbulence, specifically the cascade of mass and energy across scales, despite the absence of viscous forces in collisionless dark matter.

Methodology
The study utilizes large-scale cosmological N-body simulations, primarily the Illustris-1-Dark (dark matter only) and the Virgo consortium (SCDM) simulations. The methodology involves:

Cosmic Energy Evolution: Deriving the energy evolution equations for N-body systems in an expanding background, demonstrating that total specific energy decreases linearly with time due to Hubble expansion, acting as an effective "dissipation" mechanism.
Cascade Formulation: Defining mass and energy fluxes in two distinct domains:
- Halo Mass Space: Analyzing the transfer of mass and kinetic energy between haloes of different sizes via hierarchical merging.
- Radial Space within Haloes: Analyzing the flow of mass and energy along the radial direction of individual haloes, driven by particle migration in fluctuating gravitational potentials.
Stochastic Modeling: Modeling the migration of haloes in mass space and particles in radial space as random walks. The waiting times for these transitions are derived based on gravitational potentials and density contrasts.
Analytical Derivation: Solving the corresponding Fokker-Planck equations for these random walks to derive analytical forms for the halo mass function and density profiles.
Scaling Analysis: Identifying universal scaling laws (e.g., $-4/3$ law for density) and validating them against simulation data across various redshifts and mass scales.

Key Contributions and Results

Inverse Mass Cascade: The paper identifies an inverse mass cascade where mass flows from small to large scales via halo merging. This leads to a "random walk" of haloes in mass space.
- Halo Mass Function: The distribution of haloes is analytically derived from the Fokker-Planck equation for this random walk. The result is a "double- $\lambda$ " mass function that naturally fits simulation data across the entire mass range, outperforming standard Press-Schechter and Jenkins-Keel (JK) models, particularly for large haloes.
- Universal Scaling: A "small-scale permanence" is identified where the group mass of haloes follows a power law ( $m_g \propto m_h^{-\lambda}$ ) that is independent of redshift in the propagation range.
Energy Cascade and "Dissipation": The total energy of the dark matter system decreases continuously over time due to the expanding background. This decrease is facilitated by a dual cascade:
- Kinetic Energy: Cascades inversely (small to large scales) at a rate $\varepsilon_u$ .
- Potential Energy: Cascades directly (large to small scales) at a rate of $-1.4\varepsilon_u$ .
- Total Energy: Cascades directly to small scales at a rate of $-0.4\varepsilon_u$ , where it is effectively "dissipated" (decreased) by the Hubble expansion.
Radial Energy Flow and Density Profiles: Within individual haloes, a statistically steady state is established where the inward mass flux (accretion) balances the outward mass flux (Hubble expansion-induced migration), resulting in zero net mass flux. However, a net energy flux exists.
- Inner Density Slope: The scale-independent rate of energy flow leads to a universal inner density slope of $\rho \propto r^{-4/3}$ . This is derived both from dimensional analysis of the energy cascade and from the solution of the Fokker-Planck equation for particle random walk (yielding a "double- $\gamma$ " profile).
- Explanation of Slope Variations: The paper suggests that the wide range of inner density slopes observed in simulations ( $s \approx -1.0$ to $-1.3$) arises from non-zero net mass and energy fluxes in individual, non-equilibrium haloes, whereas the $-4/3$ slope represents the limit for composite haloes with vanishing net mass flux.
Constraints on Dark Matter Particle Mass: By applying physical constraints (waiting times and jumping lengths must exceed Planck time and Planck length, respectively) to the derived scaling laws for particle migration, the paper proposes a lower bound for the dark matter particle mass.
- Result: The analysis suggests a particle mass $m_p \gtrsim 10^{-15}$ kg ( $\approx 10^{12}$ GeV).
- Implication: This constraint excludes standard Weakly Interacting Massive Particles (WIMPs) and points toward a heavy dark matter scenario, such as superheavy right-handed neutrinos.

Significance
The paper claims to provide a unified theoretical framework for the abundance and internal structure of dark matter haloes based on the concepts of mass and energy cascades. By establishing a "statistically steady state" characterized by scale-independent cascade rates, the theory explains the self-similar and scale-free nature of halo structures without relying on specific collapse models (spherical or ellipsoidal). Furthermore, it offers a novel perspective on the "cusp-core" problem and the variations in inner density slopes, attributing them to the dynamic balance of mass and energy fluxes. Finally, the work proposes a method to constrain the fundamental properties of dark matter particles (specifically their mass) by extending the identified scaling laws to the smallest physical scales, suggesting that dark matter may be significantly heavier than currently assumed in standard cosmological models.

More like this