The Dependence of Halo Clustering on Subhalo Anisotropy and Planarity

This paper demonstrates that the clustering strength of cold dark matter host haloes depends on the anisotropy and planarity of their subhaloes, revealing a distinct environmental effect where systems with less anisotropic and planar subhaloes located further from their centers exhibit stronger clustering, independent of other known halo properties like concentration.

The Gist

Imagine the universe as a giant, invisible web of dark matter. Within this web, massive clumps called haloes form, acting like the "cities" where galaxies live. Inside these dark matter cities, there are smaller, satellite cities called subhaloes (which often host dwarf galaxies).

For a long time, astronomers have been asking: How are these satellite cities arranged inside their host cities? Do they scatter randomly like marbles in a jar, or do they line up in specific patterns, like planets in a solar system or stars in a galaxy?

This paper, written by Johnson and Zentner, investigates exactly that. They used a supercomputer simulation to look at thousands of these dark matter "cities" and discovered two major things.

1. The "Dance Floor" is Tilted

First, they confirmed what many suspected: The satellites don't scatter randomly.

Imagine a host galaxy as a giant, slightly squashed balloon (an ellipsoid). If you threw marbles (satellites) into it randomly, they would be everywhere. But the researchers found that the marbles actually prefer to sit on the "equator" of the balloon or line up along its longest axis. They are anisotropic (not the same in all directions) and planar (they like to form flat sheets or discs).

However, there's a twist. The satellites aren't as perfectly aligned as the shape of the dark matter balloon itself. They are a bit more "messy" than the balloon's shape suggests. It's like if the balloon was a perfect oval, but the marbles inside were arranged in a slightly wobbly, less perfect oval.

2. The "Social Clue": Who Hangs Out With Whom?

This is the paper's biggest new discovery. The authors asked: Does the way these satellites are arranged tell us anything about where the host galaxy lives in the universe?

Think of it like a high school party.

• The "Popular" Group: Some galaxies are very "social." They hang out in big, dense clusters with lots of other galaxies nearby.
• The "Loner" Group: Other galaxies are more isolated, hanging out in the quiet suburbs of the universe.

The researchers found a surprising link between a galaxy's "social life" (how many neighbors it has) and how its satellites are arranged:

• The "Messy" Hosts: Galaxies whose satellites are less aligned and less flat (more random, more "messy") tend to be the most popular. They are found in the densest, most crowded clusters of galaxies.
• The "Organized" Hosts: Galaxies whose satellites are very aligned and form perfect, thin planes tend to be more isolated. They hang out in quieter, less crowded areas.

The Analogy:
Imagine a group of people at a party.

• If a person's friends are standing in a tight, perfect circle around them, that person is likely standing alone in a quiet corner of the room.
• If a person's friends are scattered all over the place, bumping into others, and not forming any clear shape, that person is likely in the middle of the crowded dance floor, surrounded by other groups.

The paper found that galaxies with "scattered" satellites are the ones living in the crowded dance floors of the universe.

Why Does This Matter?

For decades, astronomers have been puzzled by the "Local Group" (our neighborhood of galaxies, including the Milky Way and Andromeda). We see our satellite galaxies arranged in thin, flat planes, which some people thought was impossible according to the standard model of the universe (Cold Dark Matter).

This paper suggests that both scenarios are possible.

• If you live in a quiet neighborhood (low density), your satellites might form a perfect, thin plane.
• If you live in a busy city (high density), your satellites might get jumbled up and become less flat.

So, the fact that we see flat planes in our neighborhood doesn't necessarily break the laws of physics; it might just mean we live in a relatively quiet part of the cosmic city.

The "Secret Sauce" Check

The authors were careful to make sure this wasn't just a side effect of something else. They checked if this pattern was just caused by:

• How heavy the galaxy is.
• How "spun up" the galaxy is (angular momentum).
• How many satellites it has.

They found that no, this is a unique effect. The way satellites are arranged (their "personality") is a direct clue to where the galaxy lives in the cosmic web, independent of its size or spin.

In a Nutshell

The universe is a giant social network. The way a galaxy's "friends" (satellite galaxies) are arranged around it acts like a social badge. If the friends are messy and scattered, the galaxy is likely in a crowded, busy cluster. If the friends are perfectly organized in a flat line, the galaxy is likely living a quieter, more isolated life. This helps astronomers understand how galaxies form and evolve in different environments.

Technical Summary

1. Problem Statement

The paper addresses two interconnected issues in the standard Cold Dark Matter (CDM) cosmological model:

1. Subhalo Anisotropy and Planarity: While it is established that satellite galaxies and subhaloes are not distributed isotropically around their host haloes (often forming planar structures), the degree to which this anisotropy depends on the environment of the host halo remains unclear.
2. Secondary Assembly Bias: Host haloes are known to cluster based on properties other than mass (e.g., concentration, spin, formation time). The authors investigate whether the clustering of host haloes depends on the spatial configuration of their subhaloes (anisotropy, alignment, planarity) and whether this dependence is a distinct effect or merely a byproduct of correlations between subhalo distribution and other known secondary halo properties (like concentration).

The study aims to determine if the spatial distribution of subhaloes is an independent predictor of host halo clustering, which would have implications for testing the CDM model against observations of satellite planes (e.g., the "Vast Polar Structures" in the Local Group).

2. Methodology

Simulation and Data

• Simulation: The authors utilized the Small MultiDark-Planck (SMDPL) simulation, a gravity-only $N$-body simulation with a box size of $400 \, h^{-1}$ Mpc and $3840^3$ particles.
• Cosmology: Planck 2013 $\Lambda$CDM parameters ($h=0.6777$, $\Omega_m=0.307$, $\sigma_8=0.8228$).
• Halo Finding: Haloes were identified using the Rockstar halo finder at redshift $z=0$.
• Sample Selection:
  • Host Haloes: Selected with a minimum mass of $1.52 \times 10^{13} \, h^{-1} M_\odot$.
  • Subhaloes: Minimum mass of $6.07 \times 10^9 \, h^{-1} M_\odot$ (63 particles).
  • Mass Ratio Cut: Subhaloes must be $>0.05\%$ of the host mass to ensure a self-similar subhalo count across host masses.
  • Final Sample: 18,441 host haloes, with 95% containing $\ge 10$ subhaloes.

Quantifying Subhalo Distributions (Marks)

The authors defined six "marks" to quantify the spatial distribution of subhaloes relative to their hosts. To isolate environmental effects from mass bias, these marks were mass-normalized (converted to percentile ranks within mass bins).

1. Alignment ($C_{50}, C_{90}$): The 50th and 90th percentiles of the absolute cosine of the angle between subhalo position vectors and the host halo's major axis.
2. Planarity ($D_{rms}$): The root-mean-square distance of subhaloes from a best-fitting plane (mass-weighted and number-weighted), normalized by the median subhalo radius.
3. Plane Orientation ($|\cos \theta_{plane}|$): The cosine of the angle between the normal of the best-fitting subhalo plane and the host's major axis.
4. Shape ($(c/a)_{sub}$): The ratio of the minor to major axes of the subhalo inertia tensor.
5. Inertia Alignment ($|\cos \theta_I|$): The cosine of the angle between the subhalo distribution's major axis and the host's major axis.
6. Radial Concentration ($Med(r_{sub})$): The median radial distance of subhaloes (normalized by virial radius).

Statistical Analysis

• Two-Point Correlation Function (TPCF): Used to visualize clustering differences between subsamples (e.g., high vs. low alignment).
• Marked Correlation Functions (MCF): The primary statistical tool. MCFs measure the correlation between the spatial separation of haloes and the values of their marks.
  • Formula: $M_m(r) = \frac{\langle m_i m_j \rangle(r) - \langle m \rangle^2}{\text{Var}(m)}$.
  • A non-zero $M_m(r)$ indicates that haloes with specific mark values cluster differently than the average.
• Control for Secondary Biases: To ensure the results were not driven by known correlations (e.g., between subhalo alignment and host concentration), the authors performed double normalization. They removed the dependence on mass and a specific secondary property (concentration, spin, shape, or subhalo count) simultaneously before calculating MCFs.

3. Key Contributions

1. Novel Quantification of Anisotropy: The paper provides a comprehensive statistical characterization of subhalo anisotropy and planarity, comparing real simulation data against two types of mocks:
   • Isotropic Mocks: Subhaloes placed randomly on spheres.
   • Ellipsoidal Mocks: Subhaloes distributed to trace the triaxial shape of the host halo's dark matter distribution.
2. Discovery of Environmental Dependence: The authors demonstrate for the first time that host halo clustering strength is a function of the subhalo spatial distribution.
3. Independence from Known Biases: Through rigorous double-normalization, the study proves that this clustering dependence is not primarily induced by correlations with standard secondary biases (concentration, spin, shape, or subhalo count). It represents a distinct, novel effect.

4. Key Results

Subhalo Distribution Properties

• Anisotropy: Subhaloes are significantly anisotropic and aligned with the major axes of their host haloes.
• Planarity: Subhaloes exhibit significant planarity (thinner planes) compared to an isotropic distribution.
• Deviation from Host Shape: Crucially, subhaloes are more isotropic and less planar than the underlying triaxial ellipsoidal shape of the host halo's dark matter distribution. They do not simply trace the host's mass ellipsoid.

Clustering Dependence (MCF Results)

The study found strong, statistically significant correlations between subhalo properties and host clustering:

• Alignment: Host haloes with less aligned subhaloes (subhaloes poorly aligned with the host major axis) cluster more strongly than those with well-aligned subhaloes.
• Planarity: Host haloes with less planar (thicker) subhalo distributions cluster more strongly.
• Radial Distribution: Host haloes where subhaloes reside at larger radii (less concentrated) cluster more strongly.

Robustness Checks

• Secondary Bias Removal: When the authors removed the effects of host concentration, spin, shape, and subhalo count from the analysis, the clustering signals for subhalo alignment ($C_{50}$) and radial position ($Med(r_{sub})$) remained statistically significant.
• Conclusion: The dependence of clustering on subhalo spatial distribution is a distinct physical effect, not an artifact of known assembly biases.

Local Group Analogs

The authors tested "Milky Way-Andromeda" (MWA) analog systems in the simulation. They found that the subhalo distributions in these specific analogs are statistically consistent with the global population, suggesting that the Local Group's satellite planes are not outliers in the context of the SMDPL simulation, though resolution limits prevent a definitive test of the most extreme thin planes.

5. Significance and Implications

• Testing CDM: The discovery that subhalo spatial configuration drives host clustering offers a new, powerful test for the standard cosmological model. If observational data (e.g., from galaxy surveys) shows different clustering trends based on satellite alignment than predicted, it could challenge CDM or the galaxy-halo connection.
• Observational Feasibility: The magnitude of the effect (e.g., $\sim 0.5\sigma$ to $0.7\sigma$ deviations in clustering at $10 \, h^{-1}$ Mpc) suggests it is measurable with current and future large-scale structure surveys.
• Galaxy Formation Physics: The result implies that the environmental history of a halo (which dictates its clustering) is intimately linked to the dynamical state and accretion history of its substructure.
• Systematic Errors: Ignoring subhalo anisotropy in clustering analyses could introduce systematic errors in interpreting large-scale structure data, particularly for weak lensing and galaxy-galaxy lensing studies where subhalo distributions bias the inferred mass profiles.

In summary, Johnson and Zentner establish that subhalo anisotropy and planarity are environmental markers that predict how strongly a host halo clusters, independent of its mass or other standard formation properties. This adds a new layer of complexity to the understanding of structure formation in the $\Lambda$CDM paradigm.