A Random Walk Model for Halo Triaxiality

This paper presents a semi-analytic random walk model that utilizes merger trees and relaxation dynamics to successfully predict the mass-dependent triaxial shapes of dark matter halos, showing strong agreement with N-body simulations while offering insights into the physical mechanisms driving their evolution.

The Gist

Imagine the universe as a giant, cosmic construction site. For decades, astronomers have been trying to understand the "buildings" on this site: dark matter halos. These are invisible, massive clouds of dark matter that act as the scaffolding for galaxies.

For a long time, scientists assumed these halos were perfect, round balls, like billiard balls. But when they looked closer using supercomputer simulations, they realized the truth: these halos are actually lumpy, stretched-out blobs. Some look like rugby balls (long and thin), others like pancakes (flat and wide), and most are somewhere in between, looking like slightly squashed eggs. This shape is called triaxiality.

The paper you provided, "A Random Walk Model for Halo Triaxiality," by Paul Menker and Andrew Benson, tries to answer a simple question: Why do these cosmic blobs have the shapes they do, and can we predict them without running massive, expensive supercomputer simulations every time?

Here is the breakdown of their discovery, explained with everyday analogies.

1. The Problem: The "Black Box" of Simulations

Currently, to figure out the shape of a dark matter halo, scientists run massive N-body simulations. Think of this like trying to predict the weather by building a giant, perfect scale model of the atmosphere in a room and watching how the wind blows. It works, but it takes weeks of supercomputer time and is incredibly expensive.

The authors wanted a faster way. They wanted a "recipe" or a "rule of thumb" that could predict the shape of a halo just by looking at its history.

2. The Solution: The "Energy Ledger" (The Random Walk)

The authors built a model based on Merger Trees.

• The Analogy: Imagine a family tree. Every person (halo) is born from two parents (smaller halos merging).
• The Innovation: In previous models, scientists tracked things that are "conserved" (like money in a bank account that doesn't disappear). But the shape of a halo isn't like money; it changes constantly.

Instead, the authors decided to track the Energy Tensor.

• The Metaphor: Think of the Energy Tensor as a 3D "stress map" or a ledger of forces inside the halo. It records how much energy is pushing the halo in the X, Y, and Z directions.
• The "Random Walk": As a halo grows, it eats smaller halos. Every time it eats one, it gets a "kick" or a "push." The authors treat these kicks like steps in a random walk. If you take a step forward, then a step left, then a step forward again, you end up in a specific spot. Similarly, every merger adds a specific "push" to the halo's energy map, changing its shape.

3. The Two Forces at Play

The model relies on two competing forces that determine the final shape:

A. The "Kicks" (Mergers)

When two halos crash into each other, they smash together.

• Head-on collision: If they hit straight on, the result tends to stretch out, making a long, thin (prolate) shape, like a rugby ball.
• Glancing blow: If they hit at an angle, they tend to flatten out, making a pancake-like (oblate) shape.
  The model calculates these "kicks" based on the history of the merger tree.

B. The "Relaxation" (Sphericalization)

Here is the tricky part. If you just kept adding kicks, the halo would get crazier and crazier. But in reality, halos tend to settle down and become more round over time.

• The Analogy: Think of a spinning pizza dough. When you first toss it, it's wobbly and irregular. But as it spins, centrifugal force and internal friction smooth it out into a circle.
• The Science: The authors call this sphericalization. Over billions of years, the chaotic orbits of particles inside the halo slowly smooth out, trying to turn the lumpy blob back into a sphere. The authors added a "damping factor" to their model to simulate this smoothing process.

4. The Results: Does the Recipe Work?

The authors tested their "recipe" against the massive supercomputer simulations (the "real" universe).

• The Calibration: They had to tweak one knob in their recipe (called the "sphericalization parameter") to make sure the average shape matched reality. Once they did that, they let the model run on its own.
• The Success:
  • Mass Matters: The model correctly predicted that heavy halos are more stretched out (rugby balls) and light halos are rounder (like billiard balls).
  • History Matters: It correctly predicted that halos in empty space (voids) look different than those in crowded clusters.
  • Speed: Instead of taking weeks of supercomputer time, their model calculates the shape in a split second.

5. Why Does This Matter?

You might ask, "Why do we care if a dark matter halo is a rugby ball or a pancake?"

• Satellite Galaxies: The shape of the halo dictates how smaller galaxies orbit inside it. If the halo is a rugby ball, the satellites line up in a specific way.
• Gravitational Lensing: Dark matter bends light. If the halo is lumpy, the way it bends light from distant galaxies changes. This affects how we measure the mass of the universe.
• Understanding the Cosmos: By linking the shape directly to the merger history, this model connects the "fingerprint" of a halo to the fundamental laws of the universe (the Big Bang and the power spectrum of density).

Summary

In simple terms, Menker and Benson created a fast, semi-analytic calculator for the shape of dark matter.

They realized that a halo's shape is the result of a tug-of-war:

1. Mergers constantly kick it into weird, lumpy shapes.
2. Time slowly smooths it back out into a sphere.

By tracking the "energy kicks" from every merger in a halo's family tree and applying a smoothing factor, they can predict the shape of a galaxy's invisible home with surprising accuracy, without needing to run a billion-year simulation. It's a "random walk" through the history of the universe that tells us exactly what the future (or present) shape of the cosmos should look like.

Technical Summary

1. Problem Statement

Dark matter halos, the gravitational scaffolding of the universe, are not perfectly spherical but are significantly triaxial (ellipsoidal). While $N$-body simulations have extensively characterized these shapes (finding them mostly prolate with a mass-dependent distribution), there is a lack of semi-analytic models that can predict halo shapes directly from their formation histories.

• Current Limitations: Existing analytic approaches often rely on Gaussian random field statistics, while $N$-body simulations are computationally expensive and do not easily isolate the causal link between specific merger events and the resulting shape.
• Goal: The authors aim to construct a semi-analytic model that predicts the triaxial shape of dark matter halos by tracking the evolution of the energy tensor through a halo's merger tree, similar to how previous works tracked angular momentum (spin) and scalar energy (concentration).

2. Methodology

The model extends the semi-analytic framework of Benson et al. (2020) and Johnson et al. (2021). Instead of tracking scalar quantities, it tracks the Energy Tensor ($\mathbf{E}$), defined as the sum of the kinetic energy tensor ($\mathbf{K}$) and the potential energy tensor ($\mathbf{W}$).

A. Core Physics & Assumptions

1. Energy Tensor Conservation: The model assumes that during a merger, the energy tensor is approximately conserved. For two merging halos (1 and 2), the total energy tensor is calculated as:
   $$\mathbf{E}^{(tot)} = \frac{\mathbf{W}^{(1)}}{2} + \left[ \frac{\mathbf{W}^{(2)}}{2} + \mathbf{V}^{(12)} + \mathbf{K}^{(orb)} \right] \left(1 + \frac{M_2}{M_1}\right)^{-\gamma}(1+b)$$
   • $\mathbf{V}^{(12)}$ and $\mathbf{K}^{(orb)}$ represent the interaction potential and orbital kinetic energy, respectively.
   • The terms $(1 + M_2/M_1)^{-\gamma}$ and $(1+b)$ are correction factors accounting for mass/energy loss to the environment and non-conservation during the merger (calibrated from Johnson et al. 2021).
2. Sphericalization: Unlike scalar energy, the energy tensor is not strictly conserved over cosmological timescales; halos tend to relax toward sphericality. The model introduces a "sphericalization" process between merger events where the eigenvalues of the energy tensor evolve toward equality:
   $$\dot{\lambda}_i = -\alpha \omega (\lambda_i - \bar{\lambda})$$
   • $\lambda_i$ are the eigenvalues.
   • $\omega$ is the precession frequency derived from the halo's flattening.
   • $\alpha$ is a dimensionless free parameter controlling the timescale of sphericalization.

B. Implementation Steps

1. Merger Trees: The model utilizes merger trees (generated via the Parkinson et al. 2008 algorithm) to track the hierarchical assembly of halos.
2. Initialization: Unresolved progenitor halos are initialized as spherical with NFW profiles.
3. Recursive Evolution:
   • At each merger node, the energy tensors of progenitors are rotated to a common frame and summed (with corrections).
   • The resulting tensor's eigenvalues are used to calculate the new axis ratios ($a_2, a_3$) of the merged halo.
   • Between mergers, the sphericalization equation is applied to evolve the tensor.
4. Calibration: The model is calibrated against $N$-body simulations (Bonamigo et al. 2015; Schneider et al. 2012) by tuning the sphericalization parameter $\alpha$ to match the median axis ratio $a_3$ at a specific halo mass ($10^{14} M_\odot$).

C. Validation

The model was first tested against idealized $N$-body simulations of isolated merging halo pairs (using Gadget-2). These simulations varied mass ratios, approach angles, and initial triaxialities to verify the model's ability to reproduce the time evolution of energy tensor components and axis ratios.

3. Key Contributions

• Novel Framework: First semi-analytic model to predict halo triaxiality by tracking the energy tensor rather than scalar energy or angular momentum.
• Physical Insight: Demonstrates that halo triaxiality is a "random walk" driven by the stochastic accumulation of energy tensors during mergers, counteracted by a relaxation process (sphericalization).
• Analytic Efficiency: Provides a rapid method to predict halo shapes directly from merger histories, linking shape statistics to the underlying cosmological power spectrum without running full $N$-body simulations.

4. Results

• Idealized Simulations: The model reproduces the evolution of axis ratios in idealized merger simulations with a median absolute deviation (MAD) of $\Delta a \approx 0.04$. It successfully captures the initial violent relaxation phase and the subsequent slow sphericalization.
• Cosmological Distributions:
  • After calibrating $\alpha = 1.75$ to match the median $a_3$ of $10^{14} M_\odot$halos, the model reproduces the **median** axis ratios ($a_2 \approx 0.69$,$a_3 \approx 0.53$) found in$N$-body simulations.
  • The model correctly predicts the mass dependence of halo shapes (higher mass halos are more prolate) and the conditioned distributions (the correlation between $a_2$ and $a_3$).
• Discrepancies:
  • The model predicts a broader distribution of axis ratios than $N$-body simulations, particularly overpredicting the number of highly triaxial halos (very low $a_3 < 0.3$).
  • The model struggles to capture complex late-time behaviors in some idealized simulations where halos temporarily become less spherical (potentially due to fallback of ejected particles), a phenomenon not included in the simple sphericalization term.

5. Significance

• Predictive Power: The model successfully connects the statistical properties of halo shapes to the fundamental cosmological model and the specific formation history of individual halos.
• Astrophysical Applications: Accurate shape predictions are crucial for understanding:
  • The anisotropic distribution of satellite galaxies and subhalos.
  • Gravitational lensing strength (which depends on halo ellipticity).
  • Cluster mass estimates and correlation functions.
• Future Directions: The authors suggest that while the current model is a robust first-order approximation, future improvements could involve:
  • Refining the sphericalization physics (e.g., accounting for tidal heating).
  • Extending the model to subhalos (affected by host tidal forces).
  • Applying the framework to self-interacting dark matter (SIDM) scenarios.
  • Generating joint distributions of spin, concentration, and triaxiality.

In conclusion, Menker & Benson provide a compelling "random walk" model that successfully explains the origins of halo triaxiality through merger history, offering a computationally efficient alternative to $N$-body simulations for studying the shape statistics of dark matter halos.