Galaxies in N-Body Simulations: Overcoming the Overmerging Problem

We present analysis of the evolution of dark matter halos in dense environments of groups and clusters in dissipationless cosmological simulations. The premature destruction of halos in such environments, known as "the overmerging," reduces the predictive power of N-body simulations and makes difficult any comparison between models and observations. We analyze the possible processes that cause the overmerging and assess the extent to which this problem can be cured with current computer resources and codes. Using both analytic estimates and high-resolution numerical simulations, we argue that the overmerging is mainly due to the lack of numerical resolution. We find that the force and mass resolution required for a simulated halo to survive in galaxy groups and clusters is extremely high and was almost never reached before: ~1-3 kpc and 108-109 M☉, respectively. We use the high-resolution Adaptive Refinement Tree (ART) N-body code to run cosmological simulations with particle mass ≈2 × 108 h-1 M☉ and spatial resolution ≈1-2 h-1 kpc and show that in these simulations the halos do survive in regions that would appear overmerged with lower force resolution. Nevertheless, the halo identification in very dense environments remains a challenge even with resolution this high. We present two new halo-finding algorithms developed to identify both isolated and satellite halos that are stable (existed at previous moments) and gravitationally bound. To illustrate the use of the satellite halos that survive the overmerging, we present a series of halo statistics, which can be compared with those of observed galaxies. Particularly, we find that, on average, halos in groups have the same velocity dispersion as the dark matter particles; i.e., they do not exhibit significant velocity bias. The small-scale (100 kpc to 1 Mpc) halo correlation function in both models is well described by the power law ξ ∝ r-1.7 and is in good agreement with observations. It is slightly antibiased (b≈0.7-0.9) relative to the dark matter. To test other galaxy statistics, we use the maximum of the halo rotation velocity and the Tully-Fisher relation to assign luminosity to the halos. For two cosmological models, a flat model with the cosmological constant and Ω0=1-ΩΛ=0.3,h=0.7 and a model with a mixture of cold and hot dark matter and Ω0=1.0,Ων=0.2,h=0.5, we construct luminosity functions and evaluate mass-to-light ratios in groups. Both models produce luminosity functions and mass-to-light ratios ( ~200-400) that are in reasonable agreement with observations. The latter implies that the mass-to-light ratio in galaxy groups (at least for Mvir ≲ 3 × 1013 h-1 M☉ analyzed here) is not a good indicator of Ω0.