The Rockstar Phase-Space Temporal Halo Finder and the Velocity Offsets of Cluster Cores
read the original abstract
We present a new algorithm for identifying dark matter halos, substructure, and tidal features. The approach is based on adaptive hierarchical refinement of friends-of-friends groups in six phase-space dimensions and one time dimension, which allows for robust (grid-independent, shape-independent, and noise-resilient) tracking of substructure; as such, it is named Rockstar (Robust Overdensity Calculation using K-Space Topologically Adaptive Refinement). Our method is massively parallel (up to 10^5 CPUs) and runs on the largest current simulations (>10^10 particles) with high efficiency (10 CPU hours and 60 gigabytes of memory required per billion particles analyzed). A previous paper (Knebe et al 2011) has shown Rockstar to have class-leading recovery of halo properties; we expand on these comparisons with more tests and higher-resolution simulations. We show a significant improvement in substructure recovery as compared to several other halo finders and discuss the theoretical and practical limits of simulations in this regard. Finally, we present results which demonstrate conclusively that dark matter halo cores are not at rest relative to the halo bulk or satellite average velocities and have coherent velocity offsets across a wide range of halo masses and redshifts. For massive clusters, these offsets can be up to 350 km/s at z=0 and even higher at high redshifts. Our implementation is publicly available at http://code.google.com/p/rockstar .
This paper has not been read by Pith yet.
Forward citations
Cited by 43 Pith papers
-
Dark Matter in Draco and Bo\"otes I: Hints of a Core in an Ultra-Faint Dwarf from Simulation-Based Inference
GraphNPE recovers a significantly lower central density for Boötes I consistent with a core while Draco remains marginally cuspy, and demonstrates that higher-order velocity moments reduce bias in dynamical modeling.
-
Cluster Mass Inference from Galaxy Kinematics
Simulation-based Deep Sets model with neural posterior estimation halves scatter in cluster mass estimates from galaxy kinematics compared to the M-sigma relation.
-
Mergers Matter: Gravothermal Collapse in Dwarf Halos with Self-Interacting Dark Matter
SIDM simulations of dwarf halos show that quiescent merger histories produce gravothermal core collapse while sustained mergers prevent collapse and can yield central densities below gravothermal fluid model predictions.
-
The AGORA High-resolution Galaxy Simulations Comparison Project. XI: Solving the Non-Spherical Morphology and Evolution of Dark Matter Halos with Haskap Pie
New non-spherical halo solving procedure applied to AGORA simulations shows merger-responsive morphology measures and redshift trends in spin, axis ratio, and overdensity.
-
The pre-infall bias of subhalos
Future subhalos show a pre-infall bias to higher progenitor masses, modeled in extended Press-Schechter theory by multiplying the collapse barrier by β(x,a)=(1-x)^{1.20+0.14a} for M200c, leading to 10-15% higher centr...
-
Informative Priors on Primordial Non-Gaussianity Bias $b_{\phi}$ From Galaxy Formation
Conditioning CAMELS-SAM simulations on the stellar mass function or stellar-to-halo mass relation reduces uncertainty in b_phi by 88-97% for DESI emission line galaxy samples while remaining consistent across galaxy f...
-
Measuring neutrino mass and asymmetry through galaxy pairwise peculiar velocity
Galaxy pairwise peculiar velocities from Cosmicflows-4 yield M_ν = 0.24^{+0.34}_{-0.18} eV and η² = 2.14^{+0.30}_{-0.32} (7σ non-zero asymmetry) in the CMB framework, consistent with prior Planck results.
-
Tomographic Alcock-Paczynski Test with Marked Correlation Functions
First integration of tomographic AP tests with MCFs and PCA compression yields 48% and 45% tighter errors on Ω_m and w versus standard two-point functions.
-
Hot or Cold? Radial Redistribution of Stars in FIRE Simulations of Milky Way-Mass Galaxies and the Asymmetry of Inward versus Outward Migrators
Simulations of Milky Way-mass galaxies reveal direction-dependent dynamical heating during stellar radial migration, with inward migrators heating, some cooling, and cold preservation uncommon especially for older stars.
-
The impact of evolving cosmic filaments on mass and spin evolution of dark matter halos
A new algorithm traces filament evolution in simulations and finds suppressed halo mass accretion plus non-random spin alignment changes near dense filaments.
-
A Unified Halo Mass Function Across Dark Matter Models from High-Resolution Multi-Scale Simulations
A calibrated fitting function for the halo mass function that unifies predictions across CDM and non-CDM models over 10 orders of magnitude in mass with typical 12% precision after modeling systematics.
-
Learning the Universe with cosmological rescaling of merger trees and semi-analytic galaxy formation models
Rescaling merger trees with a halo-profile correction enables cheap generation of galaxy summary statistics across cosmologies using semi-analytic models, matching dedicated simulation accuracy with far fewer base runs.
-
A universal model for the accretion rates and formation times of dark matter halos
A six-parameter function of peak height ν, power spectrum slope n_eff, and growth rate α_eff accurately describes median halo mass accretion rates from simulations in ΛCDM and Einstein-de Sitter cosmologies at z=0-14.
-
Euclid preparation. CosmoPostProcess: A simulation calibrated framework for weak lensing selection bias in richness-selected galaxy clusters
CosmoPostProcess delivers simulation-calibrated radial corrections for projection-induced selection bias (20-40% amplitude near 1 h^{-1} Mpc) and baryonic effects in Euclid richness-selected cluster weak lensing profiles.
-
NEFERTITI: Linking early galaxy formation to the assembly of the Milky Way
NEFERTITI simulations show that the Milky Way's most metal-poor stars largely come from a handful of accreted massive dwarf galaxies, while reproducing the JWST Hebe galaxy at z~11 as a pure Population III system.
-
Cosmological constraints from the small scale clustering of Emission Line Galaxies
SHAMe-SF modeling of small-scale DESI ELG clustering delivers 6% precision on σ8 and Ωm h², matching full DR1 results with 1% volume.
-
Constraining Neutrino Mass with the Void Weak Lensing Effect
Simulations of void-shear cross-correlation demonstrate that void lensing can constrain total neutrino mass to σ(M_ν)=0.096 eV without shape noise and 0.340 eV with Stage-III-like noise.
-
Impact of projection-induced optical selection bias on the weak lensing mass calibration of galaxy clusters
Projection-induced selection bias causes 20-50% overestimation of weak lensing masses for optically selected galaxy clusters, larger on scales >3 Mpc.
-
Universal Fitting Formulae for the Peak Concentration of Dark Matter Halos
Empirical universal fitting formula for the peak (most probable) concentration of dark matter halos derived from lognormal fits to simulation distributions and shown to be invariant across cosmologies.
-
TNG SAM: Bridging Hydrodynamical Complexity and Semi-Analytic Efficiency to Model Galaxy Formation
The TNG SAM reproduces TNG hydro simulation gas and metal flows plus galaxy and halo properties within 30% accuracy out to z=6 via five targeted updates to the Santa Cruz SAM calibrated on stellar feedback-dominated galaxies.
-
Halo mass functions in mixed cold and fuzzy dark matter models
Simulations show fuzzy dark matter fraction up to 0.3 suppresses low-mass halos in mixed DM models, and a redshift- and fraction-dependent suppression function maps CDM HMFs to MDM HMFs within 0.1-0.2 dex accuracy for z=1-4.
-
The impact of source and survey modelling on the connection between [O III] emitters and Ly $\alpha$ forest transmission at z ~ 6
Empirical halo-to-[O III] emitter modeling with realistic JWST survey mocks produces cross-correlations consistent with z~6 data within large scatter, but with a ~10 cMpc offset in the 1D peak.
-
Resolving galaxy formation in the early Universe with BonFIRE and CampFIRE
BonFIRE and CampFIRE simulations show bursty clustered star formation in early galaxies and predict UV luminosity functions matching observations at faint magnitudes with a turnover at M_UV approx -14 but overpredicti...
-
Cosmological constraints from neighbor-density-weighted marked correlation functions
Neighbor-density-weighted marked correlation functions improve FoM for Ωm–σ8 by 1.7–2.5× over standard 2PCF using Gaussian-process emulators on 129 w0waCDM+∑mν simulations.
-
Cosmological constraints from neighbor-density-weighted marked correlation functions
Neighbor-density-weighted marked correlation functions improve FoM for Ω_m–σ_8 by 1.7–2.5× over standard 2PCF using emulators from 129 w0waCDM+∑m_ν simulations.
-
Euclid preparation: Testing multi-field inflation with galaxy power spectrum and bispectrum
Validates redshift-space power spectrum and bispectrum analysis on Abacus-PNG mocks to recover unbiased f_NL constraints for Euclid spectroscopic sample.
-
Correlation between baryonic process and galaxy assembly bias
Simulations show gas cooling and stellar feedback dominate assembly bias for stellar-mass selected galaxies while star formation gives way to gas cooling for SFR-selected galaxies as number density rises.
-
How High-Specific-Energy Winds Regulate the Circumgalactic Medium of Dwarf Galaxies
Simulations demonstrate that high-specific-energy supernova outflows sustain hot CGM at virial temperature, raise t_cool/t_ff above 10, and transition dwarf galaxy feedback from ejective to preventive mode around 5 Gyr.
-
The Simulated Oxygen Shortage (SOS): Mapping the Missing OVI in Simulated Dwarf Galaxies to Subgrid Physics
The paper identifies underproduction of oxygen in low-mass simulated dwarf galaxies as the likely cause of missing OVI in the CGM, based on comparisons across two simulation suites.
-
The LISA Astrophysics MBHcatalogues Project: A comparison of predictions of simulated massive black hole binaries
A large collaboration compiles and compares merger rate predictions for massive black holes across multiple galaxy formation models to forecast LISA detections and quantify uncertainties.
-
A graph-based Neural Network surrogate model for accelerating semi-analytical model of galaxy formation and evolution
A conditional graph neural network serves as an accurate and fast surrogate for semi-analytic galaxy formation models, predicting key properties across cosmic time.
-
An Extended Parametric Model for Self-interacting Dark Matter Halos
An extended parametric model for SIDM halos incorporates mass accretion to delay core collapse and reduces errors in predicted V_max at z=0 relative to the prior Yang et al. (2024b) model.
-
Inferring Halo Mass and Scale Radius of Galaxy Clusters Using Convolutional Neural Networks and Uchuu-UniverseMachine Catalogs
Convolutional neural networks can infer galaxy cluster virial masses and scale radii from 2D projected position and line-of-sight velocity distributions with nearly unbiased results and reduced scatter when richness i...
-
ELG$\times$LRG distribution through dark matter halo dynamics
A new halo occupation model called HOMe reproduces the anisotropic clustering of ELGs and LRGs down to 200 h^{-1} kpc scales by sampling satellites from dark matter particle positions and fitting parameters to two-poi...
-
FOGGIE: Figuring Out Gas & Galaxies In Enzo XII. The Formation and Evolution of Extended HI Galactic Disks and Warps with a Dynamic Circumgalactic medium
FOGGIE cosmological simulations find that less HI-populated CGMs produce thin, coherently rotating extended disks while more populated CGMs yield irregular morphologies, with all systems settling kinematically by z=0 ...
-
The evolution of the galaxy stellar mass function and star formation rates in the COLIBRE simulations from redshift 17 to 0
COLIBRE simulations match observed galaxy stellar mass functions, star formation rates, and quenched fractions from z=17 to z=0, including JWST massive quiescent galaxies at high redshift.
-
DESI DR2 Reference Mocks: Clustering results from UCHUU ELGs and QSOs
Mock catalogs for DESI DR2 ELGs and QSOs are constructed via modified subhalo abundance matching on the Uchuu simulation to reproduce observed number density and clustering statistics.
-
Galaxy formation in modified gravity -- II. galaxy halo connection and assembly bias
Simulations in f(R) gravity show that environment density in the HOD model reduces assembly bias effects to 2-3% at z≲0.5 in both ΛCDM and modified gravity.
-
Learning the Universe with the 2nd Generation of CAMELS: Varying 35 parameters of the IllustrisTNG model in (50Mpc/h)^3 boxes
New CAMELS simulations in larger (50 Mpc/h)^3 boxes with 35 varied parameters produce tighter neural-network constraints on model parameters than prior smaller-volume runs, with public data release.
-
Little Red Dot progenitors from Compact Starbursts: A Natural Path to Early AGN Formation
High-resolution simulations produce compact galaxies where gas inflows and dynamical processes accumulate enough mass in 10 Myr to form ~10^6 solar mass central black holes under 10% feedback efficiency.
-
JWST Predictions for $z > 10$ Galaxies from the Renaissance Simulations -- I: Photometry and Sizes
Renaissance Simulations produce galaxies at z>10 whose stellar masses, star formation rates, sizes, and colors overlap with JWST observations while extending to lower masses.
-
An Improved Fit for Linear Halo Bias at High Redshift
An updated linear halo bias fit calibrated on high-redshift simulations reduces systematic offsets in early-universe clustering predictions to under 1%.
-
Galaxy and halo angular clustering in LCDM and Modified Gravity cosmologies
N-body light-cone mocks show 2-4 sigma deviations in third-order angular clustering between LCDM and f(R)/nDGP models at z=0.15-0.3 for halos and galaxies, with stronger signals in the dark-matter field.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.