Abundance and phase-space distribution of subhalos in cosmological N-body simulations: testing numerical convergence and correction methods

Abstract

Subhalos play a crucial role in accurately modeling galaxy formation and galaxy-based cosmological probes within the highly nonlinear, virialized regime. However, numerical convergence of subhalo evolution is difficult to achieve, especially in the inner regions of host halos where tidal forces are strongest. I investigate the numerical convergence and correction methods for the abundance, spatial, and velocity distributions of subhalos using two $6144^3$-particle cosmological N-body simulations with different mass resolutions -- Jiutian-300 ($1.0\times 10^7{h}^{-1}{M}_{\odot }$) and Jiutian-1G ($3.7\times 10^8{h}^{-1}{M}_{\odot }$) -- with subhalos identified by HBT+. My study shows that the Surviving subhalo Peak Mass Function (SPMF) converges only for subhalos with $m_{\mathrm{peak}}$ above $5000$ particles but can be accurately recovered by including orphan subhalos that survive according to the merger timescale model of Jiang et al., which outperforms other models. Including orphan subhalos also enables recovery of the real-space spatial and velocity distributions to $5$--$10\%$ accuracy down to scales of $0.1$--$0.2h^{-1}\mathrm{Mpc}$. The remaining differences are likely due to cosmic variance and finite-box effects in the smaller Jiutian-300 simulation. Convergence below $0.1h^{-1}\mathrm{Mpc}$ remains challenging and requires more sophisticated modeling of orphan subhalos. I further highlight that redshift-space multipoles are more difficult to recover even at larger scales because unreliable small-scale pairs at $r_{\mathrm{p}}<0.1h^{-1}\mathrm{Mpc}$ in real space affect scales of tens of $\mathrm{Mpc}$ in redshift space due to elongated Fingers-of-God effects. Therefore, for redshift-space statistics, I recommend using modified or alternative measures that reduce sensitivity to small projected separations in subhalo-based studies.

arXiv comment24 pages, 15 figures. Published in JCAP

References in corpus (47)

• A Universal Density Profile from Hierarchical Clustering
• Dark Matter Substructure in Galactic Halos
• Where are the missing galactic satellites?
• Populating a cluster of galaxies - I. Results at z=0
• Statistical Properties of X-ray Clusters: Analytic and Numerical Comparisons
• Simulating Galaxy Formation with the IllustrisTNG Model
• Hierarchical Galaxy Formation
• The Rockstar Phase-Space Temporal Halo Finder and the Velocity Offsets of Cluster Cores
• UniverseMachine: The Correlation between Galaxy Growth and Dark Matter Halo Assembly from z=0-10
• The Diversity and Similarity of Simulated Cold Dark Matter Halos
• The Connection between Galaxies and their Dark Matter Halos
• How do galaxies populate Dark Matter halos?
• A measurement of the cosmological mass density from clustering in the 2dF Galaxy Redshift Survey
• The subhalo populations of LCDM dark haloes
• Destruction of the Galactic Globular Cluster System
• Scaling relations for galaxy clusters in the Millennium-XXL simulation
• Rotation-Dependent Catastrophic Disruption of Gravitational Aggregates
• Dependence of the outer density profiles of halos on their mass accretion rate
• Dynamical Friction and Galaxy Merging Timescales
• Simulating cosmic structure formation with the GADGET-4 code
• Dark matter and cosmic structure
• Cosmology from the Chinese Space Station Optical Survey (CSS-OS)
• Dark Matter Substructure in Numerical Simulations: A Tale of Discreteness Noise, Runaway Instabilities, and Artificial Disruption
• A fitting formula for the merger timescale of galaxies in hierarchical clustering
• ECOSMOG: An Efficient Code for Simulating Modified Gravity
• The Dynamics of Sinking Satellites Around Disk Galaxies: A Poor Man's Alternative to High-Resolution Numerical Simulations
• Precise Strong Lensing Mass Modeling of Four Hubble Frontier Fields Clusters and a Sample of Magnified High-Redshift Galaxies
• An excess of small-scale gravitational lenses observed in galaxy clusters
• ELUCID - Exploring the Local Universe with reConstructed Initial Density field I: Hamiltonian Markov Chain Monte Carlo Method with Particle Mesh Dynamics
• HBT+: an improved code for finding subhalos and building merger trees in cosmological simulations
• A unified model for the spatial and mass distribution of subhaloes
• Resolving Subhaloes' Lives with the Hierarchical Bound-Tracing Algorithm
• The Preferential Tidal Stripping of Dark Matter versus Stars in Galaxies
• Orbital parameters of infalling satellite haloes in the hierarchical $\Lambda$CDM model
• Constraining Extended Gamma-ray Emission from Galaxy Clusters
• Evidence for extended gamma-ray emission from galaxy clusters
• The Aemulus Project V: Cosmological constraint from small-scale clustering of BOSS galaxies
• An analytical model for the accretion of dark matter subhalos
• The APM cluster-galaxy cross-correlation function : Constraints on Omega and galaxy bias
• Modelling Galaxy Merger Timescales and Tidal Destruction
• Internal kinematics of groups of galaxies in the Sloan Digital Sky Survey data release 7
• A Fully Self-Consistent Cosmological Simulation Pipeline for Interacting Dark Energy Models
• Influence of baryonic physics on the merger time-scale of galaxies in N-body/hydrodynamical simulations
• CSST Cosmological Emulator I: Matter Power Spectrum Emulation with one percent accuracy to k = 10 h/Mpc
• The Universal Specific Merger Rate of Dark Matter Halos
• Tracking the orbit of unresolved subhalos for semi-analytic models
• Describing the Nonuniversal Galaxy Merger Timescales in IllustrisTNG: Effects of Host Halo Mass, Baryons, and Sample Selection