Dark Matter’s Internal Life: Shaping the Universe We See

Author: Denis Avetisyan

New research explores how interactions within dark matter itself influence the growth of cosmic structures and leaves detectable signatures in the universe.

The evolution of density contrasts-illustrated for a mode <span class="katex-eq" data-katex-display="false"> k=60\,h\,\mathrm{Mpc}^{-1} </span>-demonstrates how differing velocity-dependent scattering cross-sections-defined by parameters ranging from <span class="katex-eq" data-katex-display="false"> n=-1 </span> to <span class="katex-eq" data-katex-display="false"> n=-2 </span> or incorporating a velocity scale <span class="katex-eq" data-katex-display="false"> v_{0} </span> between 0.01 and 0.1 km/s-shape the transition from superhorizon scales through phases of early subhorizon growth, conversion-driven behavior, acoustic oscillation, damping, and ultimately, late-time gravitational dominance. — The evolution of density contrasts-illustrated for a mode $k=60\,h\,\mathrm{Mpc}^{-1}$ -demonstrates how differing velocity-dependent scattering cross-sections-defined by parameters ranging from $n=-1$ to $n=-2$ or incorporating a velocity scale $v_{0}$ between 0.01 and 0.1 km/s-shape the transition from superhorizon scales through phases of early subhorizon growth, conversion-driven behavior, acoustic oscillation, damping, and ultimately, late-time gravitational dominance.

This review details the linear evolution and observational constraints on cosmological models featuring inelastic self-interacting dark matter with a mass splitting.

The enduring mystery of dark matter composition necessitates exploration beyond purely gravitational interactions. This is the central motivation behind ‘Cosmology of Inelastic Self-Interacting Dark Matter: Linear Evolution and Observational Constraints’, which investigates the cosmological effects of a two-component dark sector undergoing inelastic scattering. The authors demonstrate that exothermic conversions between dark matter species generate pressure support suppressing small-scale structure and induce observable acoustic oscillations, leaving detectable imprints on the matter power spectrum. Can these subtle signatures, constrained by Lyman-α forest data and high-redshift luminosity functions, offer a novel pathway to probe the internal thermodynamics of secluded dark sectors and unveil the true nature of dark matter?

The Illusion of Simplicity: Unveiling Dark Matter’s Hidden Layers

Contemporary cosmological simulations, built upon the premise of cold, collisionless dark matter, increasingly diverge from observed galactic structures. Specifically, simulations predict an overabundance of dwarf galaxies and overly dense galactic cores-discrepancies that suggest the standard model of dark matter is incomplete. These tensions arise from the assumption that dark matter particles interact only gravitationally, a simplification that may not reflect reality. Evidence hints at the necessity of incorporating self-interactions, or interactions within the dark sector itself, to accurately model the formation and evolution of galaxies. The inability of current models to fully align with observational data strongly motivates exploration into more complex dark matter candidates and interaction mechanisms, pushing the boundaries of our understanding of the universe’s hidden mass.

The limitations of current dark matter models, which often fail to accurately reproduce observed galactic structures, motivate exploration beyond the Standard Model and into the realm of self-interacting dark matter. Inelastic Self-Interacting Dark Matter, or Inelastic SIDM, proposes a dark sector comprised of particles that not only interact with each other, but also possess internal energy states. This means collisions between dark matter particles aren’t simply about changing direction; they can involve transitions between these energy levels, absorbing or releasing energy in the process. Such interactions, unlike those predicted by collisionless cold dark matter, can significantly alter the density profiles of dark matter halos, potentially resolving discrepancies between simulations and observations.

The Inelastic Self-Interacting Dark Matter (Inelastic SIDM) model hinges on a subtle, yet crucial, premise: a mass difference between dark matter particles. This seemingly small variation allows for inelastic scattering events – collisions where kinetic energy is exchanged, effectively ‘up-scattering’ or ‘down-scattering’ dark matter particles between different mass states. Unlike purely elastic interactions, inelastic scattering is velocity-dependent, significantly altering the velocity distribution of dark matter in galactic halos. This modification addresses persistent discrepancies between simulations employing standard collisionless dark matter and observed galactic structures, particularly in dwarf galaxies where the core-cusp problem and ‘too-big-to-fail’ issues arise. By providing a mechanism for efficient dark matter self-interaction at relatively low velocities, Inelastic SIDM can produce cored density profiles and suppress the formation of overly massive subhalos, offering a compelling pathway towards a more accurate cosmological model.

The squared transfer functions reveal that benchmark cosmological models with varying spectral indices <span class="katex-eq" data-katex-display="false">n</span> and warm dark matter velocities <span class="katex-eq" data-katex-display="false">v_0</span> exhibit distinct matter power spectrum ratios at redshifts <span class="katex-eq" data-katex-display="false">z=0</span> and <span class="katex-eq" data-katex-display="false">z=50</span>, and can be effectively approximated by equivalent thermal warm dark matter models with corresponding masses indicated at the half-mode scale <span class="katex-eq" data-katex-display="false">k_{1/2}</span>. — The squared transfer functions reveal that benchmark cosmological models with varying spectral indices $n$ and warm dark matter velocities $v_0$ exhibit distinct matter power spectrum ratios at redshifts $z=0$ and $z=50$ , and can be effectively approximated by equivalent thermal warm dark matter models with corresponding masses indicated at the half-mode scale $k_{1/2}$ .

Deconstructing the Darkness: The Architecture of a Hidden Sector

The dark sector is hypothesized to comprise two distinct dark matter particle species, denoted as $\chi_l$ and $\chi_h$ , representing the light and heavy components, respectively. These components are characterized by differing masses; $\chi_l$ possesses a smaller mass, while $\chi_h$ is significantly heavier. This mass differentiation is crucial, as it enables potential interactions within the dark sector, specifically through processes involving energy transfer between the two particle types. The existence of these two components provides a framework for modeling a self-interacting dark matter scenario, differing from single-component models and opening possibilities for explaining observed astrophysical phenomena.

The difference in mass, denoted as $\Delta m$ , between the light ( $\chi_l$ ) and heavy ( $\chi_h$ ) dark matter components directly influences the energy available during inelastic scattering events. Specifically, the magnitude of $\Delta m$ determines the energy transfer possible when a dark matter particle scatters with another, impacting the cross-section for self-interaction. A larger $\Delta m$ implies a greater energy exchange, potentially increasing the self-interaction rate, while a smaller difference limits the available energy and consequently reduces the interaction probability. This relationship is crucial as the self-interaction rate affects the dark matter’s distribution within galactic halos and its potential detectability through indirect means.

A thermally isolated dark sector is characterized by weak couplings to Standard Model particles, effectively limiting energy exchange between the dark and visible sectors. This isolation is not absolute; rather, it’s maintained while allowing for internal energy redistribution through interactions within the dark sector, specifically between the light component $χ_l$ and the heavy component $χ_h$ . This internal energy transfer occurs without significant leakage to Standard Model particles, preserving the dark sector’s thermal properties and enabling processes like dark matter self-interaction. The degree of thermal isolation is determined by the strength of the couplings to Standard Model particles, while the efficiency of internal energy transfer is dictated by the mass splitting $\Delta m$ and interaction rates between $χ_l$ and $χ_h$ .

The evolution of the heavy-component fraction <span class="katex-eq" data-katex-display="false">r_h = n_h / (n_h + n_l)</span> and light-component temperature <span class="katex-eq" data-katex-display="false">T_l</span> with redshift <span class="katex-eq" data-katex-display="false">1+z</span> for different parameterizations (BP1, BP2, BP3) reveals sensitivity to both the mass splitting <span class="katex-eq" data-katex-display="false">\Delta m / m_l</span> at a fixed <span class="katex-eq" data-katex-display="false">m_l = 100\\,\\mathrm{MeV}</span> and the light mass <span class="katex-eq" data-katex-display="false">m_l</span> at a fixed <span class="katex-eq" data-katex-display="false">\Delta m / m_l = 10^{-2}</span>, with shaded bands representing a factor-of-ten cross-section variation. — The evolution of the heavy-component fraction $r_h = n_h / (n_h + n_l)$ and light-component temperature $T_l$ with redshift $1+z$ for different parameterizations (BP1, BP2, BP3) reveals sensitivity to both the mass splitting $\Delta m / m_l$ at a fixed $m_l = 100\\,\\mathrm{MeV}$ and the light mass $m_l$ at a fixed $\Delta m / m_l = 10^{-2}$ , with shaded bands representing a factor-of-ten cross-section variation.

Echoes of Expansion: Modeling Cosmological Evolution

Background Evolution Equations are utilized to model the expansion of the dark sector under the assumption of spatial homogeneity. These equations govern the time evolution of scale factor $a(t)$ and are derived from the Friedmann equations, incorporating contributions from dark energy and dark matter. Crucially, the model allows for energy transfer between “light” components, such as dark energy with an equation of state $w \neq -1$ , and “heavy” components, typically dark matter, through a transfer term $Q$ . This transfer term represents the rate of energy exchange per unit volume and is a key parameter in determining the dark sector’s overall evolution and its impact on cosmological observables. The equations account for both energy density and pressure contributions from each component, providing a self-consistent framework for investigating dark sector interactions.

Perturbation equations model the evolution of small deviations from the homogeneous background cosmology, specifically tracking the growth of density and velocity fluctuations in the early universe. These equations are crucial for understanding structure formation because they account for the influence of inelastic scattering processes – collisions between particles that do not conserve kinetic energy – on the dynamics of these fluctuations. Inelastic scattering introduces a damping effect, altering the amplitude and timescale of structure growth; this is particularly relevant in scenarios involving dark matter self-interactions or interactions with other relativistic species. The solutions to these perturbation equations predict how initial fluctuations evolve into the large-scale structures observed today, and are fundamental for interpreting observational probes of cosmic structure.

Predictions from our cosmological models regarding the growth of large-scale structure are quantitatively compared with observations of the Lyman-alpha forest and UV luminosity functions. The Lyman-alpha forest, resulting from absorption by intervening neutral hydrogen clouds, provides a mapping of the distribution of matter along the line of sight. Similarly, the UV luminosity function traces the abundance of galaxies at high redshifts, reflecting the underlying matter distribution. By analyzing the statistical properties of these datasets – specifically the power spectrum of density fluctuations – we can constrain the parameters governing the dark sector’s evolution. This work demonstrates that these observational constraints are sensitive to the internal thermodynamics of the dark sector, allowing for probes of energy transfer rates and equation of state parameters beyond those accessible from cosmic microwave background observations alone.

The ratio of inelastic reaction rate to the Hubble rate, <span class="katex-eq" data-katex-display="false">\mathcal{R}/(aH)</span>, reveals that inelastic conversions exceed the expansion rate (marked by the dotted line at <span class="katex-eq" data-katex-display="false">\mathcal{R}/(aH)=1</span>) for power-law (left panel) and saturated (right panel) cross-section parametrizations, with varying parameters as detailed in the text. — The ratio of inelastic reaction rate to the Hubble rate, $\mathcal{R}/(aH)$ , reveals that inelastic conversions exceed the expansion rate (marked by the dotted line at $\mathcal{R}/(aH)=1$ ) for power-law (left panel) and saturated (right panel) cross-section parametrizations, with varying parameters as detailed in the text.

Beyond the Horizon: Refining the Dark Matter Narrative

To truly understand how dark matter organizes itself in the universe, researchers are turning to N-body simulations that model the gravitational interactions of millions of particles. These aren’t simply linear approximations; they delve into the non-linear regime of structure formation – the chaotic, complex stage where gravity’s effects become dominant and simple calculations fail. By tracking the evolution of these simulated universes, scientists can generate realistic depictions of dark matter halo formation, the gravitational ‘wells’ where galaxies ultimately reside. These simulations allow for a detailed examination of halo shapes, internal structure, and the distribution of matter within them, offering crucial insights that are unattainable through analytical methods and providing a powerful means to test the predictions of Inelastic Self-Interacting Dark Matter against observations of the cosmos.

Rigorous statistical analysis, specifically employing Markov Chain Monte Carlo (MCMC) methods, is paramount to precisely determine the values of the free parameters within the Inelastic Self-Interacting Dark Matter (SIDM) model. These analyses don’t simply yield best-fit values; they meticulously map the probability distributions of each parameter, revealing the uncertainties and correlations that govern the model’s behavior. Crucially, MCMC techniques allow researchers to quantitatively assess how well the model reproduces observed astronomical data – from the distribution of galaxies to the properties of dark matter halos – by calculating statistical measures of goodness-of-fit. This process moves beyond qualitative agreement, providing a robust framework for determining whether Inelastic SIDM is a viable explanation for the observed universe or requires further refinement, and ultimately, distinguishes it from competing dark matter theories with statistical confidence.

The convergence of sophisticated N-body simulations and rigorous Markov Chain Monte Carlo (MCMC) analyses represents a pivotal step towards validating or refuting the Inelastic Self-Interacting Dark Matter (SIDM) model. Current dark matter paradigms struggle to reconcile simulations with observed galactic structures, prompting exploration of alternative interactions beyond simple gravitational effects. This combined approach allows for a comprehensive examination of SIDM’s parameter space, testing its ability to accurately reproduce observed galactic dynamics and large-scale structure formation. Success in this endeavor could not only resolve existing discrepancies but also fundamentally reshape cosmological understanding, offering a more complete picture of the universe’s composition and evolution, and potentially revealing the nature of dark matter itself.

Varying parameters around fiducial values modulates the linear power spectra of light (solid lines) and heavy (dashed lines) components, expressed as ratios to the ΛCDM matter power spectrum at <span class="katex-eq" data-katex-display="false">z=0</span>, as described by <span class="katex-eq" data-katex-display="false">P\\_{\\chi\\_{i}}(k)\\sim eq P\\_{\\Lambda\\mathrm{CDM}}(k)\\,[\\delta\\_{i}(k)/\\delta\\_{\\Lambda\\mathrm{CDM}}(k)]^{2}</span>. — Varying parameters around fiducial values modulates the linear power spectra of light (solid lines) and heavy (dashed lines) components, expressed as ratios to the ΛCDM matter power spectrum at $z=0$ , as described by $P\\_{\\chi\\_{i}}(k)\\sim eq P\\_{\\Lambda\\mathrm{CDM}}(k)\\,[\\delta\\_{i}(k)/\\delta\\_{\\Lambda\\mathrm{CDM}}(k)]^{2}$ .

The pursuit of dark matter’s true nature feels increasingly like charting a course toward an event horizon. This study, detailing the cosmological evolution within an inelastic self-interacting dark matter (SIDM) framework, meticulously maps the interplay between internal thermodynamic processes and large-scale structure formation. It’s a precise calculation, yet one must remember, as Erwin Schrödinger observed, “The total number of states of a system is finite.” The elegance of the perturbation equations, detailing how these interactions imprint themselves on the Lyman-alpha forest, is undeniable. However, the model, however sophisticated, exists until it collides with data – and the universe has a habit of obscuring its secrets.

What Lies Beyond the Horizon?

The exploration of inelastic self-interacting dark matter, as detailed in this work, reveals not so much a path to understanding, but a refinement of the questions. The cosmos generously shows its secrets to those willing to accept that not everything is explainable. To model dark matter’s internal dynamics is to invite increasingly complex thermodynamic landscapes, each capable of obscuring the initial conditions with its own entropic haze. The precision with which structure formation can constrain such models will undoubtedly improve, but the fundamental challenge – discerning signal from noise, theory from self-deception – will only sharpen.

Future investigations will likely focus on bridging the gap between simplified collision kernels and the full complexity of multi-component dark sectors. However, it is crucial to remember that increasingly sophisticated models do not necessarily equate to increasing proximity to truth. The Lyman-alpha forest, while a powerful probe, offers only a limited window onto the universe’s dark underbelly. Other observational avenues, such as gravitational lensing and the cosmic microwave background, must be leveraged to create a more holistic picture.

Ultimately, this line of inquiry serves as a potent reminder: black holes are nature’s commentary on our hubris. Each successful constraint, each refined parameter space, merely illuminates the vastness of what remains unknown. The pursuit of dark matter’s secrets, therefore, is not about solving a mystery, but about acknowledging the exquisite limitations of the tools-and the minds-employed in its unraveling.

Original article: https://arxiv.org/pdf/2604.15006.pdf

Contact the author: https://www.linkedin.com/in/avetisyan/

The Illusion of Simplicity: Unveiling Dark Matter’s Hidden Layers

Deconstructing the Darkness: The Architecture of a Hidden Sector

Echoes of Expansion: Modeling Cosmological Evolution

Beyond the Horizon: Refining the Dark Matter Narrative

What Lies Beyond the Horizon?

See also: