---
title: 'Symfind: Particle-Tracking Subhalo Finder'
url: https://www.emergentmind.com/topics/symfind
type: topic
---

# Symfind: Particle-Tracking Subhalo Finder

Searching arXiv for the cited Symfind papers and closely related context.
Symfind is a particle-tracking-based subhalo finder developed for cosmological $N$-body simulations to address known failure modes of conventional single-epoch halo-finding pipelines, especially artificial subhalo disruption and false convergence in inferred subhalo statistics. In its original $\Lambda$CDM formulation, Symfind begins from a conventional halo catalog and merger tree, tags the most-bound particles of a subhalo at first infall, and subsequently identifies descendants by locating the density peak containing the majority of those tagged particles, followed by iterative unbinding [2308.10926]. Subsequent work has shown that the same fixed-core strategy that is advantageous in cold dark matter can become problematic in self-interacting dark matter (SIDM), where self-interactions and tides can diffuse originally tagged core particles outward and strip them from the remnant, leading to premature loss of otherwise surviving subhalos [2507.09799].

## 1. Conceptual basis and motivation

Symfind was introduced to overcome two specific failings attributed to traditional “single-epoch” pipelines such as Rockstar + consistent-trees (RCT): artificial disruption of resolved subhalos at masses orders of magnitude above those predicted by idealized, high-resolution simulations, and false convergence of subhalo statistics with increasing simulation resolution [2308.10926]. The central premise is that a subhalo should be identified by temporal continuity of its most tightly bound material rather than solely by instantaneous phase-space overdensity.

In this framework, the defining object is not merely a density peak at a given snapshot, but a tracked remnant anchored to a pre-selected “core” of bound particles. This design is intended to reduce confusion between a genuine subhalo remnant and surrounding tidal debris. In cold dark matter simulations, the method is reported to be insensitive to temporary confusion between the true subhalo core and tidal material and to recover more highly stripped subhalos at small host-centric radii than conventional phase-space methods such as Rockstar [2507.09799].

A plausible implication is that Symfind shifts the ontology of subhalo identification from snapshot-local classification to branch-continuous remnant tracking. That distinction is central to why its convergence properties differ from those of RCT in the published comparisons [2308.10926].

## 2. Algorithmic pipeline

Symfind begins with a conventional halo catalog and merger tree, specifically Rockstar + consistent-trees in the published implementation, and corrects merger-tree errors such as spurious short-lived branches and aphysical central-subhalo label swapping during mergers [2308.10926]. For each branch that eventually becomes a subhalo, it walks the branch backward in time and collects every $N$-body particle that was ever within the branch’s virial radius, $R_{\rm vir}$, before first infall, separating smoothly accreted from non-smoothly accreted particles [2308.10926].

At the snapshot of first infall, Symfind computes binding energies and flags the $N_{\rm core}$ most-bound particles as the subhalo core. The fiducial choice is $N_{\rm core}=32$ [2308.10926]. The SIDM analysis describes the same step as ranking all particles belonging to a newly detected subhalo by their binding energy and tagging the top $N_{\rm core}=32$ most-bound particles at the “infall snapshot” [2507.09799].

At each later snapshot, Symfind re-identifies candidate density peaks and uses the pre-selected core particles to decide which peak is the descendant. In the 2023 description, candidate peaks are identified among smoothly accreted particles using Subfind with kernel-smoothed density over the $k$ nearest neighbors, and the peak containing the largest number of core particles is taken as the subhalo center of mass and velocity [2308.10926]. The fiducial neighbor count is $k=16$, and the pair $(k,N_{\rm core})=(16,32)$ yields an estimated instantaneous error rate below $1\%$ for subhalos with $n_{\rm peak}>300$ [2308.10926]. The SIDM study describes an equivalent descendant-identification stage as a local-density-peak search using a small Friends-of-Friends run applied to the tagged particles and their neighbors [2507.09799]. In both descriptions, the descendant is defined by the density peak containing the largest fraction of the original tagged core.

After the center is identified, Symfind iteratively unbinds tracked particles and computes halo properties from the resulting bound set. The SIDM description gives the unbinding loop explicitly: gather all particles within $3\times$ the tidal radius or within $R_{\rm vir}$ if available, compute each particle’s total energy relative to the peak, remove particles with $E_{\rm tot}>0$, and repeat until convergence [2507.09799]. The tracked bound set is then used to derive quantities such as mass $m$, half-mass radius $r_{1/2}$, $v_{\rm max}$, and $r_{\rm max}$ [2308.10926].

Symfind declares a subhalo disrupted if any of the following hold: all core particles lie outside $r_{1/2}$, the subhalo’s distance to its host drops below its own $r_{1/2}$, or it is completely unbound [2308.10926]. It also continues searching in later snapshots to recover transient misassignments and interpolates properties across short gaps [2308.10926].

By construction, the tagged core set never grows. New tightly bound particles are incorporated only if they lie near the density peak containing the tagged core particles [2507.09799]. This design is essential to both the method’s robustness in CDM and its vulnerability in some SIDM regimes.

## 3. Quantitative behavior in cold dark matter

The principal quantitative claim for Symfind in $\Lambda$CDM is that it tracks subhalos to substantially lower surviving masses than commonly used tools such as Rockstar and consistent-trees [2308.10926]. In the Symphony dark-matter-only simulations, the reported gains are approximately $15\%-40\%$ more subhalos within the virial radius and approximately $35\%-120\%$ more subhalos within $R_{\rm vir}/4$ at fixed peak subhalo mass [2308.10926]. At eight-times higher resolution in SymphonyMilkyWayHR, the gains increase to approximately $25\%-130\%$ inside $R_{\rm vir}$ and approximately $60\%-250\%$ inside $R_{\rm vir}/4$ [2308.10926].

The published comparison also emphasizes changes in inferred radial structure. Symfind returns subhalo cumulative radial distributions that steepen with resolution and approach the host dark-matter particle cumulative distribution, whereas RCT shows no change with resolution [2308.10926]. In the same study, the low-mass slope of the $m_{\rm peak}$ mass function is shallower for Symfind than for Rockstar by approximately $-0.04$ inside $R_{\rm vir}$ and $-0.11$ inside $R_{\rm vir}/4$ [2308.10926].

The interpretation advanced in the source is that single-epoch finders frequently miss low-mass remnants or misidentify them during disruption because subhalos can lose large fractions of their mass before disappearing [2308.10926]. Symfind’s particle-tracking scheme is intended to maintain continuity through precisely this regime. The same study argues that Symfind can trace resolved subhalos until the point of typical galaxy disruption without invoking orphan modeling, provided the numerical resolution is sufficient [2308.10926].

## 4. Resolution thresholds, convergence, and validation

A major contribution of the original Symfind paper is the distinction between resolving mass loss and resolving internal structural quantities. Idealized simulations are used to argue that resolving subhalo mass loss rates requires $n_{\rm peak}\gtrsim4\times10^3$ particles at peak, whereas structural properties such as $v_{\rm max}$ require $n_{\rm peak}\gtrsim3\times10^4$ [2308.10926]. These thresholds are reiterated as recommended best practice for sample selection [2308.10926].

The same work formulates numerical disruption thresholds through fitted functions
\[
n_\mathrm{lim,mass}(n_{\rm peak})=n_\mathrm{lim,ideal}(8\,n_{\rm peak}), \qquad
n_\mathrm{lim,vmax}(n_{\rm peak})=n_\mathrm{lim,ideal}(n_{\rm peak}),
\]
with
\[
n_\mathrm{lim,ideal}(n_\star;q)=10^{\,b_2(\log_{10}n_\star)^2+b_1\log_{10}n_\star+b_0},
\]
where the best-fit coefficients are given in Table 4 of the source paper [2308.10926]. The disruption analysis uses survival-analysis methodology. Letting $\mu=m/m_{\rm peak}$ and $\mu_{\rm disrupt}$ denote the minimum $\mu$ before a branch disappears, the Kaplan-Meier estimator is written as
\[
\widehat{\Pr}(\mu_{\rm disrupt}<\mu_i)=\prod_{\mu_j\ge \mu_i}\left(1-\frac{d_j}{N(\le \mu_j)}\right),
\]
with Greenwood variance
\[
\widehat{\Var}[\widehat{\Pr}(<\mu)]=\widehat{\Pr}(<\mu)^2\sum_j\frac{d_j}{N(\le\mu_j)\,[N(\le\mu_j)-d_j]}.
\]
These formulae are used to account for censored subhalos surviving to $z=0$ [2308.10926].

In the reported comparison for subhalos with $10^{4.5}<n_{\rm peak}<10^5$, Symfind’s median $\mu_{\rm disrupt}$ is approximately $30$–$100$ times smaller than Rockstar’s [2308.10926]. Rockstar’s disruption fraction is described as essentially independent of $n_{\rm peak}$, producing a flat $\mu_{\rm disrupt}$; this is identified as a textbook example of false convergence [2308.10926]. Symfind, by contrast, is reported not to falsely converge because $\mu_{\rm disrupt}(n_{\rm peak})$ decreases with $n_{\rm peak}$ until numerical limits are reached [2308.10926].

The Symfind framework also proposes four validation tests for any subhalo finder: Kaplan-Meier survival curves for disruption thresholds, structural convergence using the $v_{\rm max}/v_{\rm max,inf}$ versus $m/m_{\rm infall}$ relation, mass-loss convergence at fixed $(t-t_{\rm infall})/t_{\rm cross}$, and core-particle consistency during disruption [2308.10926]. This suggests that Symfind is not only an algorithm but also a methodological program for certifying subhalo catalogs against numerical and physical criteria.

## 5. Reformulation in self-interacting dark matter

The 2025 SIDM study re-examines Symfind in a regime where the physics alters the behavior of the very particles used to define the tracked core [2507.09799]. The simulations adopt a velocity-dependent scattering cross section with Rutherford-like angular dependence,
\[
\frac{d\sigma}{d\cos\theta}=
\frac{\sigma_0\cdot w^4}{2\,[w^2+v^2\sin^2(\theta/2)]^2},
\]
where $\sigma_0/m$ is the normalization, $w$ is the turnover velocity, and $v$ is the relative speed of colliding particles [2507.09799]. An approximate phenomenological form is also given:
\[
\sigma(v)m^{-1}\simeq
\sigma_0m^{-1}\cdot\frac{w^4}{(w^2+v^2)^2}.
\]
The study compares Symfind with RCT on Milky-Way- and Group-mass zoom hosts in four models: SIDM70, SIDM147-Group, SIDM147-MW, and a CDM baseline [2507.09799].

Several auxiliary quantities are defined to characterize why Symfind can fail in SIDM. The instantaneous tidal radius is
\[
r_t=r\left[\frac{M_{\rm sub}}{3\,M_{\rm host}(<r)}\right]^{1/3},
\]
for subhalo mass $M_{\rm sub}$ at host-centric radius $r$ [2507.09799]. A diffusion radius is introduced to describe heat-conduction-driven migration of core particles,
\[
r_{\rm diff}\simeq \sqrt{6D\Delta t},
\]
with diffusion coefficient
\[
D\sim \frac{1}{3}(\sigma/m)\rho v^3 l_{\rm scatt}\sim \frac{1}{3}v^2 l_{\rm scatt}/\tau_{\rm heat},
\]
and $l_{\rm scatt}\sim(\sigma\rho)^{-1}$ [2507.09799]. In practice, $r_{\rm diff}$ is measured as the radius enclosing half of the originally tagged cores after time $\Delta t$ [2507.09799]. The study also defines a core-particle-loss fraction after pericenter,
\[
f_{\rm loss}\equiv 1-\frac{N_{\rm core}(t_{\rm after})}{N_{\rm core}(t_{\rm before})}.
\]

The key physical result is that the fixed-core technique does not always yield accurate results in SIDM [2507.09799]. Self-interactions transfer heat from hotter outer regions to colder inner regions, causing originally tagged core particles to move onto larger orbits; if those particles reach radii $r\gtrsim r_t$, the host tidal field can strip them away [2507.09799]. The same mechanism can cause Symfind to lose the tracked core even though a bound remnant survives.

The paper distinguishes two regimes. In the core-expansion phase, subhalos develop large constant-density inner cores, with $r_{\rm core}\sim10$–$20\,{\rm kpc}$ for massive $M_{\rm peak}\sim10^{11}\,M_\odot$ systems; after tidal pruning of the outer NFW wings, the remaining isothermal core particles diffuse outward and are rapidly stripped, leading Symfind to lose the object prematurely [2507.09799]. In the core-collapse phase, central densities rise and diffusion is subdominant, so core particles remain bound through repeated pericenters; here Symfind can track farther into the host center than RCT, which often confuses the subhalo with host outskirts at small radii [2507.09799].

## 6. Comparative performance and hybrid catalog construction

The SIDM study condenses the comparison between Symfind and RCT into four metrics: mass completeness, tracking longevity, false-loss rate, and model dependence on $\sigma/m$ and $w$ [2507.09799]. The reported behavior is summarized below.

| Metric | Reported result | Regime |
|---|---|---|
| Mass completeness | Symfind recovers $\sim20\%$ more subhalos within $0.5\,R_{\rm vir}$ than RCT | CDM |
| Mass completeness | Symfind’s gain falls to $\sim5\%$ | SIDM70 |
| Mass completeness | Symfind under-performs RCT by $\sim10\%$ for massive subhalos in the inner $0.25\,R_{\rm vir}$ | SIDM147-Group |
| Mass completeness | Symfind outperforms by $\sim10\%$ | SIDM147-MW |
| Tracking longevity | Symfind is $\sim3\,{\rm Gyr}$ longer than RCT for $M_{\rm peak}>10^{10}\,M_\odot$ | CDM |
| Tracking longevity | Symfind can be lost $\sim1\,{\rm Gyr}$ earlier than RCT for core-expansion subhalos | SIDM147-Group |
| False-loss rate | $<5\%$ | CDM |
| False-loss rate | rises to $\sim10\%$ | SIDM70 |
| False-loss rate | can exceed $20\%$ for the most massive core-expansion subhalos | SIDM147-Group |

The model dependence is stated explicitly: higher $\sigma/m$ and higher $w$ produce larger isothermal cores, more core-particle migration, larger $r_{\rm diff}$, larger $f_{\rm loss}$, and degraded Symfind performance, whereas smaller $w$ places many low-mass subhalos in core collapse, where tightly bound cores suffer little diffusion and Symfind excels [2507.09799].

Because neither finder is uniformly superior across all SIDM regimes, the study recommends a hybrid workflow that merges Symfind and RCT outputs [2507.09799]. The procedure is: run both finders; pre-match subhalos one-to-one by particle IDs or most-bound center; if Symfind’s core-loss fraction after infall exceeds approximately $30\%$ and $M_{\rm peak}>10^{10}\,M_\odot$, adopt the RCT descendant; else if $r_{\rm pericenter}<0.2\,R_{\rm vir}$ and the subhalo is in core collapse with rising $\rho_{\rm inner}$, adopt the Symfind result; otherwise use whichever finder reports the later disruption time; finally, include subhalos found in only one catalog but flag the finder of origin [2507.09799]. The resulting combined catalog recovers approximately $2$–$15\%$ more total subhalos and up to approximately $30\%$ more subhalos at $r<0.25\,R_{\rm vir}$ than either method alone across a wide range of SIDM parameters [2507.09799].

This suggests that Symfind’s methodological value persists outside CDM, but only when paired with diagnostics sensitive to the underlying microphysics. In SIDM, the core particles are not merely tracers; their transport properties become part of the failure analysis.

## 7. Relation to orphan modeling and broader significance

A recurring implication of the Symfind literature is that the longevity of tracked subhalos bears directly on the need for orphan prescriptions in satellite-galaxy modeling. The 2023 study argues that, for $\Lambda$CDM applications, one can avoid orphan modeling whenever $n_{\rm peak}\gtrsim4\times10^3$, and that studies requiring reliable $v_{\rm max}$ should ensure $n_{\rm peak}\gtrsim3\times10^4$ [2308.10926]. If orphans are still required, the recommendation is to generate them at the simulation’s numerical disruption threshold rather than at the last surviving snapshot reported by the halo finder, and to use multiple core particles to track positions so as to mitigate diffusion errors [2308.10926].

The SIDM results qualify that conclusion. In CDM, fixed core-particle tracking is a source of robustness because it resists confusion with tidal debris and host background [2507.09799]. In SIDM, the same fixed core set can become physically non-representative of the surviving remnant after heat conduction and tidal stripping [2507.09799]. Thus, the general lesson is not that particle tracking is universally superior, but that its validity depends on whether the tagged particles remain a stable proxy for remnant identity.

Within the subhalo-finding literature, Symfind is therefore significant in two senses. First, it provides a concrete alternative to single-epoch finders, with explicit numerical tests, survival-analysis diagnostics, and published resolution criteria [2308.10926]. Second, its later SIDM reassessment clarifies the boundary conditions of that alternative: the durability of the tracked core is model-dependent, and catalog construction may require hybridization with phase-space methods when new physics alters core transport [2507.09799].

Source: https://www.emergentmind.com/topics/symfind