The Accuracy of Subhalo Detection

With the ever increasing resolution of N-body simulations, accurate subhalo detection is becoming essential in the study of the formation of structure, the production of merger trees and the seeding of semi-analytic models. To investigate the state of halo finders, we compare two different approaches to detecting subhaloes; the first based on overdensities in a halo and the second being adaptive mesh refinement. A set of stable mock NFW dark matter haloes were produced and a subhalo was placed at different radii within a larger halo. SUBFIND (a Friends-of-Friends based finder) and AHF (an adaptive mesh based finder) were employed to recover the subhalo. As expected, we found that the mass of the subhalo recovered by SUBFIND has a strong dependence on the radial position and that neither halo finder can accurately recover the subhalo when it is very near the centre of the halo. This radial dependence is shown to be related to the subhalo being truncated by the background density of the halo and originates due to the subhalo being defined as an overdensity. If the subhalo size is instead determined using the peak of the circular velocity profile, a much more stable value is recovered. The downside to this is that the maximum circular velocity is a poor measure of stripping and is affected by resolution. For future halo finders to recover all the particles in a subhalo, a search of phase space will need to be introduced.Comment: 9 pages, 7 figures, accepted for publication in MNRA

The Impact of Box Size on the Properties of Dark Matter Haloes in Cosmological Simulations

We investigate the impact finite simulation box size has on the structural and kinematic properties of Cold Dark Matter haloes forming in cosmological simulations. Our approach involves generating a single realisation of the initial power spectrum of density perturbations and studying how truncation of this power spectrum on scales larger than L_cut affects the structure of dark matter haloes at z=0. In particular, we have examined the cases of L_cut = f_cut L_box with f_cut=1 (i.e. no truncation), 1/2, 1/3 and 1/4. In common with previous studies, we find that the suppression of long wavelength perturbations reduces the strength of clustering, as measured by a suppression of the 2-point correlation function xi(r), and reduces the numbers of the most massive haloes, as reflected in the depletion of the high mass end of the mass function n(M). Interestingly, we find that truncation has little impact on the internal properties of haloes. The masses of high mass haloes decrease in a systematic manner as L_cut is reduced, but the distribution of concentrations is unaffected. On the other hand, the median spin parameter is ~50% lower in runs with f_cut<1. We argue that this is an imprint of the linear growth phase of the halo's angular momentum by tidal torquing, and that the absence of any measurable trend in concentration and the weak trend observed in halo shape reflect the importance of virialisation and complex mass accretion histories for these quantities. These results are of interest for studies that require high mass resolution and statistical samples of simulated haloes, such as simulations of the population of first stars. Our analysis shows that large-scale tidal fields have relatively little effect on the internal properties of Cold Dark Matter haloes and hence may be ignored in such studies.Comment: Submitted to MNRAS; 10 pages, 12 figure