🔭 astrophysics

Growth of Structure in Multi-species Wave Dark Matter

This paper presents a comprehensive theoretical framework for analyzing the growth of structure and deriving power spectra in multi-species wave and particle dark matter, accounting for arbitrary component mixtures, initial conditions, and scale-dependent effects like free-streaming and Jeans scales.

Original authors: Mustafa A. Amin, M. Sten Delos

Published 2026-05-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Mustafa A. Amin, M. Sten Delos

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 Big Picture: A Cosmic Orchestra

Imagine the universe's dark matter isn't just one single type of invisible stuff. Instead, think of it as a cosmic orchestra made of many different instruments (species) playing together. Some instruments are heavy and slow (like a tuba), some are light and fast (like a flute), and some are so light they act more like waves rippling across a pond than solid particles.

This paper provides a new "score" (a mathematical framework) to predict how this orchestra creates the structure of the universe—how galaxies and clusters form—when these different instruments are playing at the same time.

The Main Characters: Particles vs. Waves

The authors study two main types of dark matter "instruments":

The Particle Players (The "Warm" Ones): Imagine a crowd of people running around a field. They have some speed (warmth) and bump into each other. If they run too fast, they can't stick together to form a tight group. In physics, this is called "free-streaming."
The Wave Players (The "Fuzzy" Ones): Imagine a giant, invisible ocean wave. Because they are waves, they have a natural "quantum pressure" that keeps them from collapsing into a tiny point. They act like a fuzzy cloud rather than a sharp dot.

The paper looks at scenarios where you have a mix of these. Maybe 90% of the dark matter is a slow, heavy particle (like standard Cold Dark Matter), and 10% is a light, wavy, "warm" field.

The Problem: How Do They Grow Together?

In the early universe, tiny ripples in density (like small bumps on a trampoline) started to grow. Gravity pulled matter into these bumps to eventually form galaxies.

The authors asked: What happens when a "wave" species and a "particle" species are mixed?

Does the "warmth" of the running particles wash out the small bumps?
Does the "fuzziness" of the waves prevent them from clumping?
Can a small, noisy "wave" component actually help the big "particle" component grow, even if the wave part is weak?

The Solution: A New Mathematical Recipe

The team derived a set of equations (a recipe) that calculates exactly how the "power spectrum" (a measure of how much clumping happens at different sizes) evolves over time.

Think of the Power Spectrum as a graph showing how much "clumpiness" exists at different scales:

Big scales: Galaxy clusters.
Small scales: Individual galaxies or dwarf galaxies.

Their recipe accounts for three main forces that fight against clumping:

Free-streaming: Particles running away too fast to stick together.
Thermal Pressure: The "warmth" of the particles pushing them apart.
Quantum Pressure: The "fuzziness" of the waves keeping them spread out.

Key Findings (The "Plot Twist")

1. The "Noise" Can Be Helpful
Usually, we think of "noise" (random fluctuations) as bad. But in this paper, the authors found that even a tiny, sub-dominant "wave" component can create its own random "noise" (called isocurvature fluctuations).

Analogy: Imagine a quiet room (the main dark matter) where a few people are clapping randomly (the wave component). Even though the clappers are few, their clapping can start a rhythm that the whole room eventually follows.
Result: This "noise" from the small wave component can actually seed growth in the dominant particle component, making small structures form more easily than expected.

2. Different Scales, Different Rules
The paper identifies specific "cutoff" sizes where things change:

The "Free-Stream" Cutoff: If particles are too fast, they erase large-scale bumps.
The "Jeans" Cutoff: If the pressure is too high (either from heat or quantum waves), small bumps get smoothed out.
The "Fuzzy" Cutoff: For wave dark matter, there is a specific size below which waves simply refuse to clump, no matter how much gravity pulls.

3. The "Mixed" Universe Looks Different
When you mix a cold particle with a warm wave, the resulting pattern of galaxies isn't just a simple average of the two.

Scenario A: If the main dark matter is a "fuzzy" wave and the small part is warm, the small part gets erased quickly, but the main wave part still forms a specific pattern of clumps.
Scenario B: If the main dark matter is standard particles and the small part is a warm wave, the small part's "noise" can boost the formation of tiny structures, but the "warmth" eventually stops them from growing too big.

The Tool

The authors didn't just write equations; they built a computer code (available online) that lets other scientists plug in different masses and mixtures of dark matter to see what the universe would look like. This allows researchers to compare their theories with real telescope data (like the Lyman-alpha forest or galaxy surveys) to figure out what dark matter actually is.

Summary

This paper is a "user manual" for a universe with a complex, multi-species dark matter sector. It explains how different types of dark matter (particles vs. waves, cold vs. warm) interact, how they suppress or boost the formation of galaxies, and provides the tools to calculate exactly what we should see if our universe is made of this complex mixture.

Technical Summary: Growth of Structure in Multi-species Wave Dark Matter

Problem Statement
The nature of the dark sector remains unknown, with possibilities ranging from a single bosonic field (e.g., the QCD axion) to a multi-species ensemble containing bosonic and fermionic components. While previous works have addressed the growth of structure in single-species wave dark matter or multi-species particle dark matter, a comprehensive framework for calculating the gravitational growth of structure in a multi-species system comprising both wave-like (bosonic fields) and particle-like components was lacking. Such systems may arise in high-energy theories like the "string axiverse," where different fields possess distinct masses, density fractions, and initial conditions (e.g., globally misaligned scalars vs. locally misaligned fields or vector dark photon components). The challenge lies in modeling how these diverse components, potentially exhibiting "warmth" (random motion) and intrinsic Poisson-like fluctuations, interact gravitationally to shape the matter power spectrum.

Methodology
The authors develop a theoretical framework based on the Schrödinger–Poisson system for $N$ non-relativistic scalar fields $\phi_s(t, \mathbf{x})$ . The evolution is analyzed using Fourier transforms of the fields and their complex conjugates, treating them as random variables with specified initial statistics.

Key methodological steps include:

Bilinear Definition: The authors define a bilinear quantity $\hat{f}^s_k(p)$ related to the phase-space distribution, allowing the fractional density perturbation $\delta$ to be expressed as an integral over momentum space.
Evolution Equations: By applying the Schrödinger–Poisson system, they derive evolution equations for the correlation functions of these bilinears. These equations account for the coupling between different species via the gravitational potential.
Decomposition of Perturbations: The total density power spectrum is decomposed into two distinct contributions:
- Adiabatic Perturbations: Arising from primordial curvature perturbations, evolving similarly to standard Cold Dark Matter (CDM) but modified by the multi-species dynamics.
- Isocurvature Perturbations: Arising from intrinsic Poisson-like fluctuations (shot noise) inherent to the wave nature of the fields, particularly relevant for post-inflationary production mechanisms or locally misaligned fields.
Volterra Equations: The time evolution of the power spectrum is governed by a set of coupled Volterra integral equations. These equations involve "free-streaming kernels" ( $T^{fs}$ ) that depend on the initial field power spectra ( $f^s_0(p)$ ) and the mass of each species. The kernels capture effects such as free-streaming, effective pressure from random motion ("warmth"), and wave-dynamical "quantum pressure."
Numerical Implementation: The authors provide a numerical algorithm to solve these coupled Volterra equations iteratively. A public code is released to calculate density power spectra for arbitrary numbers of species, density fractions, and initial field power spectra.

Key Contributions and Results
The paper derives general expressions for the total density contrast power spectrum $P_\delta(k)$ and the cross-spectra between different species ( $P^{ss'}_\delta(k)$ ).

General Formalism: The total power spectrum is expressed as the sum of an adiabatic term and an isocurvature term, each evolved via specific transfer functions ( $T^{(ad)}_k$ and $T^{(iso)}_k$ ). These transfer functions are solutions to Volterra equations that incorporate the weighted sum of free-streaming kernels from all species.
Relevant Scales: The authors identify and analyze several critical scales that suppress structure growth:
- Free-streaming scale ( $k_{fs}$ ): Erases initial adiabatic perturbations.
- "Warm" Jeans scale ( $k_J$ ): Suppresses growth due to effective pressure from random motion.
- "Fuzzy" Jeans scale ( $k_j$ ): Suppresses growth due to wave-dynamical quantum pressure.
- Initial Isocurvature Scale ( $k_*$ ): The scale at which the initial Poisson-like power spectrum peaks.
Interplay of Components: The framework demonstrates that a subdominant warm wave component can significantly influence the total power spectrum.
- Isocurvature Boost: Even if a component is subdominant in mass, its intrinsic isocurvature fluctuations can seed perturbations in the dominant component, potentially boosting the total power spectrum at small scales.
- Suppression Mechanisms: The combination of free-streaming, warm Jeans suppression, and fuzzy Jeans suppression leads to complex, scale-dependent features in the power spectrum. For instance, in scenarios with a dominant globally misaligned scalar and a subdominant warm wave component, the total spectrum shows suppression at high $k$ due to the dominant component's fuzzy pressure, while the subdominant component's isocurvature seeds growth that is subsequently suppressed by its own warmth.
Limits: The formalism naturally recovers known limits:
- Globally Misaligned Wave DM: Corresponds to $k_* \to 0$ , yielding negligible isocurvature and warmth, behaving like standard fuzzy dark matter.
- Particle DM: Corresponds to $k_* \to \infty$ with fixed velocity dispersion, recovering standard warm/cold particle dark matter results.
Numerical Examples: The paper presents numerical results for two-component scenarios (e.g., CDM + warm wave DM, or globally misaligned wave DM + warm wave DM). These examples illustrate how varying the mass fraction of the warm component alters the power spectrum, showing that subdominant components can leave distinct imprints on the matter power spectrum at scales accessible to Lyman- $\alpha$ forest and galaxy surveys.

Significance
The paper claims to provide a necessary and general framework for studying structure formation in complex dark matter scenarios that cannot be reduced to a single effective fluid. By unifying wave-dynamical calculations with multi-species interactions, it allows for the investigation of:

Scenarios where dark matter consists of multiple fields with different masses and production mechanisms (e.g., the string axiverse).
The impact of intrinsic Poisson fluctuations (isocurvature) on structure growth, which can be significant even for subdominant components.
The interplay between different suppression mechanisms (free-streaming, thermal pressure, quantum pressure) in multi-component systems.

The authors emphasize that their results are general enough to include various dark matter candidates, from ultralight scalars to vector fields and classical particles. They provide a publicly available code to facilitate the calculation of power spectra for these complex models, enabling comparisons with observational data to constrain the nature of the dark sector. The work serves as a foundation for future studies on nonlinear structure formation and the observational signatures of multi-species dark matter.