The role of energy shear in the collapse of protohaloes

This paper demonstrates that strong correlations between the trace and traceless parts of the energy shear tensor in positive-definite matrices allow for a refined parameterization of the initial energy overdensity, thereby explaining scatter in collapse thresholds and enabling accurate predictions of protohalo properties.

The Gist

Imagine the universe as a giant, invisible ocean of dark matter. In the beginning, this ocean was mostly calm, but with tiny ripples and waves. Over billions of years, gravity acted like a giant magnet, pulling the water in these ripples together to form islands. These islands are dark matter haloes, the invisible scaffolding that holds galaxies together.

This paper is a detective story about how these islands form. The authors, Marcello Musso and Ravi Sheth, are trying to figure out exactly what makes a specific spot in that early ocean decide, "Yes, I will collapse into a galaxy cluster today," while its neighbors just float along.

Here is the breakdown of their discovery, using some everyday analogies.

1. The Three-Directional Squeeze

Usually, we think of gravity as a simple squeeze from all sides, like a ball being squashed. But the universe is messy. The authors focus on a concept called Energy Shear.

Think of a piece of dough on a table.

• If you just press down on it, it spreads out.
• But if you have a specific force pushing it from the top, bottom, left, and right all at once, it collapses into a tight ball.

For a protohalo (a baby galaxy cluster) to form, the forces acting on it must be strong enough to crush it in all three dimensions (up-down, left-right, front-back). The authors call this "positive definiteness." If the forces aren't strong enough in even one direction, the dough just stretches out into a long noodle (a filament) instead of a ball.

The Discovery: They found that for a spot to become a halo, the "squeeze" must be positive in all three directions. This is a strict rule that filters out most random spots in the universe.

2. The Mystery of the "Special" Spots

In a completely random, calm ocean (what physicists call a "Gaussian random field"), the strength of the squeeze in one direction has nothing to do with the squeeze in another. They are independent.

However, the authors looked at the actual "baby haloes" in computer simulations and found something weird: The squeezes were strongly correlated. If the squeeze was strong in one direction, it was almost always strong in the others.

The Analogy: Imagine a lottery where you usually pick numbers completely at random. But suddenly, you notice that every time the winning ticket has a "7" in the first spot, it also has a "7" in the second spot. That's suspicious! It means the winners aren't random; they are special.

The authors realized this correlation happens simply because of the rule of collapse. You can't have a collapse in all three directions unless the forces are balanced and strong enough. This mathematical "rule" forces the numbers to line up, creating the correlation we see. It's not magic; it's just the physics of squeezing dough.

3. The "Threshold" Problem

Knowing that a spot can collapse isn't enough. We need to know when it collapses.

• Some spots collapse quickly.
• Some take a long time.
• Some never collapse at all.

The authors wanted to find a "Threshold"—a specific number that tells us, "If the squeeze is stronger than this, you become a galaxy today."

Previously, scientists thought this threshold was just a single number (like a speed limit). But the data was messy. Some spots had the same "squeeze" but collapsed at different times.

The New Insight: The authors realized the threshold isn't just one number; it's a combination of numbers.

• Imagine you are trying to fill a bucket (the halo).
• The "main number" is how much water you have (the total energy).
• But the "shape" of the bucket matters too. If the bucket is lopsided (uneven squeeze), you need more water to fill it to the brim than if it were a perfect cylinder.

They found a clever way to combine the "shape" factors into a new variable (which they call $v_+$). When they plotted the data using this new variable, the messy scatter disappeared. The "speed limit" became a clear, straight line.

4. The "Magic" Rescaling

When they compared their math to the computer simulations, they found the math was slightly off. The simulations showed haloes forming with slightly more "squeeze" than the math predicted.

The Fix: They realized their math assumed the baby haloes were perfect spheres (like a marble). But in reality, they are lumpy and irregular (like a potato).

• The Analogy: If you try to measure the volume of a potato using the formula for a marble, you'll get the wrong answer.
• They applied a simple "correction factor" (multiplying by 1.35) to account for the fact that haloes are lumpy potatoes, not perfect marbles. Suddenly, their math matched the simulations perfectly.

Why Does This Matter?

This paper is important because it simplifies a very complex problem.

1. It explains the "Assembly Bias": It helps us understand why some galaxies form in crowded neighborhoods and others in lonely places. The "shape" of the initial squeeze determines the final location.
2. It's a New Tool: Instead of needing super-complex computer simulations to guess where galaxies form, we can now use these simple mathematical rules (the "squeeze" and the "shape") to make very accurate predictions.
3. Future Tech: The authors mention that this could help train Artificial Intelligence to predict where galaxies will form, which is crucial for mapping the universe in the future.

In a Nutshell:
The universe is full of random fluctuations, but the ones that turn into galaxies are special because they are squeezed hard enough in all three directions. By understanding the exact shape of that squeeze, the authors found a simple "recipe" to predict exactly where and when these cosmic islands will form. They turned a chaotic mess of numbers into a clear, predictable pattern.

Technical Summary

Here is a detailed technical summary of the paper "The role of energy shear in the collapse of protohaloes" by Musso and Sheth.

1. Problem Statement

Dark matter haloes form from the gravitational collapse of matter around specific locations in the initial density field where flows converge. While the standard model of halo formation often relies on the initial matter overdensity ($\delta$) or the deformation tensor (derived from the gravitational potential), recent work suggests the energy overdensity tensor (or energy shear tensor, $u_{ij}$) is more physically relevant for collapse dynamics.

The central problem addressed in this paper is twofold:

1. Correlation of Invariants: In a generic Gaussian random field, the rotational invariants of a tensor (trace, traceless shear amplitude, and traceless determinant) are statistically independent. However, observations of protohaloes show strong correlations between these invariants. The authors seek to determine if these correlations arise simply from the physical requirement that protohaloes must collapse along all three axes (positive definiteness of the energy shear).
2. Collapse Threshold: While positive definiteness ensures an object will collapse, it does not determine when. Previous work established a critical threshold for the trace (energy overdensity, $\epsilon$) required for collapse by the present day ($z=0$), but this threshold exhibits significant scatter. The authors aim to explain this scatter and formulate a more precise, parameterized threshold condition.

2. Methodology

The study combines analytical statistical mechanics with numerical validation using cosmological simulations.

• Theoretical Framework:

  • Energy Shear Tensor ($u_{ij}$): Defined as the source of the evolution of the inertia tensor. For a protohalo to collapse triaxially, $u_{ij}$ must be positive definite (all three eigenvalues $\lambda_1, \lambda_2, \lambda_3 > 0$).
  • Rotational Invariants: The authors analyze three invariants:
    • $\epsilon = \text{tr}(u)$ (Energy overdensity).
    • $q = \sqrt{3/2 \sum (\lambda_i - \epsilon/3)^2}$ (Traceless shear amplitude).
    • $\mu = (9/2)(U_3/q^3)$, where $U_3$ is the traceless determinant.
  • Constraint Analysis: They derive the conditional probability distributions of these invariants given the constraint $\lambda_3 > 0$ (positive definiteness). They introduce new variables $v_+ = \lambda_1 + \lambda_2 - 2\lambda_3$ and $v_- = \lambda_1 - \lambda_2$ to simplify the positivity constraint ($\epsilon \ge v_+$).
  • Threshold Ansatz: They propose a new collapse threshold condition: $\epsilon = \sqrt{\epsilon_c^2 + v_+^2}$, where $\epsilon_c \approx 2.2$ is a critical value and $v_+$ accounts for the scatter.

• Numerical Validation:

  • Simulations: Data from the SBARBINE simulation suite (specifically the Bice and Flora simulations) containing $1024^3$ particles.
  • Samples: Protohaloes were identified at $z=0$ using a Spherical Overdensity finder. The study analyzed four mass bins ranging from $10^{11} h^{-1}M_\odot$to$>10^{15} h^{-1}M_\odot$.
  • Measurement: For each protohalo, the energy shear tensor was reconstructed from initial conditions, diagonalized, and the invariants ($\epsilon, q, \mu$) and derived variables ($v_+, v_-$) were calculated.

3. Key Contributions

1. Origin of Invariant Correlations: The paper demonstrates that the strong correlations observed between the trace ($\epsilon$) and the shear ($q$ or $\mu$) in protohaloes are a direct mathematical consequence of the positive definiteness constraint ($\lambda_3 > 0$). In a generic Gaussian field, these are independent; in positive definite matrices, they are strongly coupled.
2. Efficient Parameterization of Collapse: The authors introduce $v_+$ (a linear combination of eigenvalues) as a superior parameter for describing the collapse threshold compared to the standard invariants $q$ and $\mu$. The scatter in the relationship between $\epsilon$ and $v_+$ is significantly narrower than in the $\epsilon$-$q$ plane.
3. Analytical Threshold Model: They derive an analytical expression for the collapse threshold:
   $$\epsilon \approx \sqrt{\epsilon_c^2 + v_+^2}$$
   This formula explains the stochasticity of the collapse threshold, linking the energy overdensity required for collapse to the anisotropy of the initial patch.
4. Mass Dependence and Rescaling: The study identifies that while the theoretical distributions depend on mass via the smoothing scale $\sigma_0^2$, the observed distributions of $q$ require a rescaling factor ($\sigma_\epsilon \approx 1.35 \sigma_0$). This is attributed to the non-spherical nature of protohalo patches compared to the spherical volumes assumed in the analytical derivation.

4. Key Results

• Correlation Structure: The scatter plot of $\epsilon$ vs. $q$ shows a lower bound determined by $\epsilon > q X(\mu)$, where $X(\mu)$ is a function derived from the positivity constraint. This bound explains why protohaloes with lower $\mu$ (more prolate) require higher $\epsilon/q$ ratios.
• Distribution of Invariants:
  • $\epsilon$: The distribution of the energy overdensity for protohaloes is shifted to higher values compared to a Gaussian distribution. The mean is $\approx 1.66$ (in units of $\sigma_0$).
  • $q$: The distribution of the shear amplitude is less mass-dependent than $\epsilon$ but is systematically higher than predicted by the unconstrained positive-definite theory. A rescaling of the variance by 1.35 brings the theory into excellent agreement with simulation data.
  • $\mu$: The distribution of the third invariant $\mu$ shifts towards 1 (spherical symmetry) for lower mass haloes, consistent with the constraint $\epsilon/q$ being tighter for small objects.
• Threshold Validation: The proposed ansatz $\epsilon = \sqrt{\epsilon_c^2 + v_+^2}$ (with $\epsilon_c \approx 2.2$) successfully reproduces the measured distribution of protohaloes across different mass bins. The residual scatter in $\epsilon$ at fixed $v_+$ is small, suggesting $v_+$ captures the dominant physics of the collapse threshold.
• Mass Regimes:
  • Low Mass: Dominated by the positivity constraint ($\epsilon > q$).
  • High Mass: Dominated by the collapse-time constraint ($\epsilon > \epsilon_c \approx 2.2$), making them less sensitive to the shear anisotropy.

5. Significance

• First-Principles Derivation: The paper moves away from complex peak theory and excursion set calculations for predicting secondary properties of protohaloes. Instead, it shows that conditioning on a physically motivated threshold (positivity + collapse time) is sufficient to predict the statistical properties of protohaloes analytically.
• Assembly Bias: By clarifying the role of the third invariant ($\mu$) and the variable $v_+$, the work provides a better framework for understanding "assembly bias" (how halo properties depend on environment beyond mass).
• Cosmological Modeling: The results offer a more accurate way to model the initial conditions of dark matter haloes, which is crucial for interpreting large-scale structure surveys and constraining cosmological parameters.
• Machine Learning Applications: The analytical mapping from initial field properties to collapse outcomes provides a robust physical prior for training machine learning algorithms to predict halo formation from initial conditions.

In summary, Musso and Sheth establish that the statistical properties of protohaloes are largely dictated by the mathematical requirement of positive definiteness and a simple geometric threshold for collapse time, providing a unified and analytically tractable description of halo formation physics.