STRAWBERRY: Finding haloes in the gravitational potential

The paper introduces STRAWBERRY, a novel algorithm that identifies dark matter haloes by distinguishing bound and unbound particles through an accelerated reference frame's gravitational potential, revealing that haloes function as two-component systems comprising a stable, virialized core and a rapidly evolving unbound component.

The Gist

Here is an explanation of the paper "STRAWBERRY: Finding haloes in the gravitational potential," translated into simple language with everyday analogies.

The Big Picture: The Cosmic Neighborhood Problem

Imagine the universe as a giant, expanding city. In this city, Dark Matter Haloes are like neighborhoods or cities themselves. They are massive clumps of invisible matter that hold galaxies together.

For decades, astronomers have struggled with a simple question: "Where does a neighborhood end, and the open countryside begin?"

• The Old Way: Previously, scientists drew a circle around a neighborhood. If the density of people (matter) inside the circle was, say, 200 times higher than the average density of the whole country, they called it a "city."

  • The Problem: This is like saying a city ends exactly where the streetlights stop being bright. It's arbitrary. Real neighborhoods aren't perfect circles; they are messy, lumpy, and have suburbs that fade out slowly. Also, as the country expands, the definition of "bright streetlights" changes, making it hard to track how the city grows over time.

• The New Way (STRAWBERRY): The authors of this paper created a new tool called STRAWBERRY. Instead of looking at how crowded a place is, they look at energy and gravity. They ask: "If you were a particle (a tiny speck of matter) here, could you run away and never come back, or are you stuck in this neighborhood?"

The Core Idea: The "Boosted" Potential

To understand STRAWBERRY, imagine you are standing on a moving train.

• The Problem: If you try to measure the shape of the train car while the train is speeding up, the floor feels like it's tilting. It's hard to tell where the car actually ends because the whole world around you is accelerating.
• The Solution: STRAWBERRY puts you in a "free-falling" reference frame. It's like imagining you are floating right next to the train, moving at the exact same speed. Suddenly, the train feels stationary. You can see its true shape without the "tilt" of the universe's expansion messing up your view.

In this "boosted" view, the scientists look for a Saddle Point.

• The Analogy: Imagine a landscape of hills and valleys. A "halo" is a deep valley.
  • Bound Particles: These are like balls sitting at the bottom of the valley. They don't have enough energy to roll up the hill and escape.
  • Unbound Particles: These are balls rolling down the hill from a higher mountain. They might dip into the valley for a moment, but they have too much speed to stay. They will roll right back out.
  • The Saddle Point: This is the lowest point on the ridge separating your valley from a deeper valley next door. If a particle has enough energy to get over this ridge, it belongs to the other valley (or is escaping).

STRAWBERRY finds this "ridge" automatically and says, "Everything inside this ridge is part of the city. Everything outside is just passing through."

How the Algorithm Works (The "Strawberry" Method)

The name "STRAWBERRY" stands for Structure Tracking Refining And Binding Refining Refining Y... well, actually, it's just a catchy name, but here is how it works step-by-step:

1. Pick a Seed: Start at the very bottom of the deepest valley (the center of the halo).
2. Fill the Bowl: Imagine pouring water into the valley. The water spreads out, filling the bottom and moving up the sides.
3. Find the Spill: Keep filling until the water reaches the "saddle point" (the lowest point of the rim).
4. Stop: Anything the water touches is "bound" (part of the halo). Anything the water didn't reach is "unbound" (just passing through).

The algorithm does this mathematically for millions of particles, checking if they have enough "speed" to climb out of the valley.

What They Discovered

By using this new method, the team found some fascinating things about how these cosmic neighborhoods behave:

1. The "First Dip" Rule
They found that particles usually become "official" members of the neighborhood after they swing around the center for the first time.

• Analogy: Imagine a person moving to a new city. They might visit for a day (unbound), but they don't really "belong" until they've had their first full day of work and sleep (the first orbit). Once they've done that, they are part of the community.

2. Two Distinct Groups
The halo isn't a uniform blob. It has two distinct populations:

• The Residents (Bound): These are the particles that have settled down. They are "virialized" (a fancy word for "calm and stable"). They orbit the center in a predictable way, like planets around a sun. They form a distinct, finite shape with a clear edge.
• The Commuters (Unbound): These are particles currently falling in or just bouncing off the edge (splashback). They are chaotic, fast, and haven't settled yet. They form a fuzzy cloud around the residents.

3. The Edge is Real
Old methods suggested haloes fade out forever. STRAWBERRY shows that the "Residents" have a sharp edge. Beyond a certain point, the particles are no longer part of the halo; they are just falling in from the outside.

Why This Matters

This new way of defining haloes is like switching from measuring a city by a fixed map grid to measuring it by where the actual people live and work.

• Better Science: It helps astronomers understand how galaxies form and evolve without the confusion of "fake" boundaries.
• Universal Truths: It suggests that if we look at the "Residents" only, the rules of the universe might look the same everywhere, regardless of when in history we are looking.
• Future Tools: This method can help simulate the universe more accurately, helping us understand dark matter, dark energy, and the ultimate fate of our cosmic neighborhood.

In summary: STRAWBERRY is a smart tool that stops guessing where a galaxy cluster ends and starts. Instead, it asks, "Who is stuck here, and who is just passing through?" By answering that, it reveals the true, stable heart of the universe's structure.

Technical Summary

Here is a detailed technical summary of the paper "STRAWBERRY: Finding haloes in the gravitational potential" by Richardson, Stücker, and Angulo.

1. Problem Statement

Defining the boundaries and mass of dark matter haloes in cosmological N-body simulations is a fundamental challenge. Current standard methods suffer from significant theoretical and practical limitations:

• Spherical Overdensity (SO): Defines haloes based on a fixed overdensity threshold (e.g., $M_{200b}$). This enforces spherical symmetry on inherently ellipsoidal structures, ignores substructure, and suffers from "pseudo-evolution" (mass increases due to background density changes rather than accretion).
• Friends-of-Friends (FoF): Links particles based on a linking length. While efficient, it lacks a sharp physical boundary, can artificially bridge distinct haloes, and often misidentifies chance alignments or numerical fragments as haloes.
• Dynamical Definitions: Methods based on orbital history (e.g., particles that have passed apocenter) are physically motivated but computationally expensive as they require tracking particle trajectories over time.

There is a need for a definition that is parameter-free, relies on instantaneous information, accounts for the global gravitational potential (not just self-potential), and distinguishes between bound (virialized) and unbound (accreting) matter.

2. Methodology: The STRAWBERRY Algorithm

The authors propose a novel algorithm, STRAWBERRY (Structure assignment within boosted reference frames in cython), which identifies haloes by analyzing the boosted gravitational potential.

Theoretical Framework

• Boosted Potential: The algorithm operates in an accelerated reference frame that removes large-scale gradients. The potential is defined as:
  $$\phi_b(\mathbf{x}) = \phi(\mathbf{x}) - (\mathbf{x} - \mathbf{x}_0) \cdot \nabla \phi|_{\mathbf{x}_0}$$
  This is equivalent to studying the halo in a frame free-falling with the local acceleration $\mathbf{a}_0$, isolating the local potential well from the cosmic expansion and large-scale tidal fields.
• Turn-Around Potential: To ensure the existence of a saddle point (necessary for defining a boundary), the authors introduce a "turn-around" correction based on spherical collapse theory. This forces a saddle point at the turn-around radius ($\tilde{\delta} \approx 4.55$), ensuring a well-defined escape energy even in pathological cases.
• Binding Criterion: A particle is considered bound if its total energy (kinetic + boosted potential) is less than the energy at the saddle point ($\phi_{sad}$) that connects the local potential minimum to a deeper minimum.
  $$E_{particle} < E_{\infty} = \phi_{sad} + \text{kinetic terms at saddle}$$

Algorithmic Steps

1. Connectivity: Establishes a graph of nearest neighbors (default $N=10$) for all particles to define local topology.
2. Initialization: Uses an initial FoF catalog to locate potential wells. It identifies the particle with the lowest potential as the seed.
3. Well Filling: Iteratively grows a group of "internal particles" from the minimum. It checks surface particles; if a surface particle has neighbors with lower potential that are not yet in the group, it indicates a path to a deeper valley.
4. Saddle Detection: If a path to a deeper valley is found, the algorithm traces it to find the saddle point. The difference in potential between the minimum and the saddle is the persistence ($\Delta \phi$).
5. Binding Check: Calculates the energy of all particles within the saddle contour. Those with $E < E_{sad}$ are bound.
6. Substructure Handling: Subgroups (subhaloes) are treated as single objects with their own bulk velocity to prevent particles from "cheating" the binding check due to the motion of the substructure relative to the host.

3. Key Contributions

• Parameter-Free Definition: Unlike SO or FoF, STRAWBERRY requires no arbitrary thresholds (like linking lengths or overdensity factors $\Delta$). The boundary is determined dynamically by the potential landscape.
• Global Potential Inclusion: It naturally includes the gravitational influence of surrounding structures, improving upon "self-binding" checks that ignore external tides.
• Instantaneous Analysis: Unlike orbital history methods, it determines bound status using only the snapshot's instantaneous positions, velocities, and potentials, making it computationally efficient.
• Two-Component Halo Model: The algorithm explicitly separates haloes into a bound, virialized core and an unbound, rapidly evolving envelope.

4. Key Results

Using a $\Lambda$CDM simulation (Illustris TNG300 cosmology, $N=625^3$ particles), the authors analyzed the properties of the identified haloes:

• Binding Timescales: Particles typically become bound between their first pericenter and apocenter passage.
  • ~50% of particles are bound at the first pericenter.
  • ~75% are bound at the first apocenter.
  • ~90% are bound after two dynamical times.
• Phase Space Structure:
  • Bound Population: Occupies a triangular region in phase space ($r$ vs. $v_r$), is symmetric in radial velocity, and is virialized.
  • Unbound Population: Consists of infalling streams (negative $v_r$) and splashback streams (positive $v_r$). These are not virialized and dominate the outer regions.
• Density Profiles:
  • The bound mass profile follows an NFW/Einasto shape but exhibits a sharp exponential cutoff at the boundary (approx. $r_{70\%}$), leading to a finite, convergent mass.
  • The total profile (bound + unbound) does not show this sharp cutoff, explaining why traditional definitions struggle to define a finite halo edge.
  • The bound profile is remarkably stable across redshifts, while the unbound envelope evolves rapidly with the accretion rate and background density.
• Virialization:
  • The virial ratio ($G/T$) for STRAWBERRY-selected bound particles is consistent with theoretical expectations ($\approx 1$).
  • In contrast, standard methods (like SUBFIND) include unbound particles, resulting in a virial ratio biased low by 10–25%.
• Mass Comparison:
  • At $z=0$, the bound mass ($M_{bound}$) tracks the standard spherical overdensity mass ($M_{200b}$).
  • At higher redshifts, $M_{bound}$ becomes more extended relative to $M_{200b}$. The physical radius $r_{200b}$ contains a smaller fraction of bound matter at high $z$ because the bound region expands while the background density drops.

5. Significance and Future Prospects

• Universal Halo Mass Function: The stability of the bound profile suggests that $M_{bound}$ may provide a more universal definition of halo mass, potentially resolving the non-universality observed in standard $M_{200b}$ mass functions.
• Physical Insight: The method provides a rigorous way to distinguish between virialized equilibrium and non-equilibrium accretion, crucial for studying gas physics, tidal stripping, and the assembly history of haloes.
• Efficiency: By avoiding time-consuming orbit integration, STRAWBERRY enables the analysis of large-scale simulations with a physically robust definition of halo boundaries.
• Code Availability: The algorithm is open-source and publicly available, facilitating its adoption in the community for future cosmological studies.

In conclusion, STRAWBERRY offers a paradigm shift from density-based to energy-based halo definitions, providing a physically motivated, parameter-free, and computationally efficient tool to study the true dynamical boundaries of dark matter structures.