2004/03/05 by Itai Arad, I. Arad, Avishai Dekel +3
Physics and Astronomy · #Astronomy and Astrophysical Research #Galaxies: Formation, Evolution, Phenomena #Scientific Research and Discoveries #astro-ph
paper · pdf · doi:10.1111/j.1365-2966.2004.08045.x
published as Mon.Not.Roy.Astron.Soc. 353 (2004) 15 · 16 pages, 11 figures, submitted to mnras
arxiv created 2004/03/05 · openalex publication_date 2004/08/18 · arxiv updated 2009/12/01 · openalex created_date 2016/06/24 · openalex updated_date 2026/07/28
We present a method for computing the six-dimensional coarse-grained phase-space density f(x, v) in an N-body system, and derive its distribution function v(f). The method is based on Delaunay tessellation, where v(f) is obtained with an effective fixed smoothing window over a wide f range. The errors are estimated, and v(f) is found to be insensitive to the sampling resolution or the simulation technique. We find that in gravitationally relaxed haloes built by hierarchical clustering, v(f) is well approximated by a robust power law, v(f) ∝f−2.5±0.05, over more than four decades in f, from its virial level to the numerical resolution limit. This is tested to be valid in the Λ cold dark matter cosmology for haloes with masses 109–1015 M⊙, indicating insensitivity to the slope of the initial fluctuation power spectrum. By mapping the phase-space density in position space, we find that the high-f end of v(f) is dominated by the ‘cold’ subhaloes rather than the parent-halo central region and its global spherical profile. The value of f in subhaloes near the virial radius is typically >100 times higher than its value at the halo centre, and it decreases gradually from the outside in toward its value at the halo centre. This seems to reflect phase mixing due to mergers and tidal effects involving puffing up and heating. The phase-space density can thus provide a sensitive tool for studying the evolution of subhaloes during the hierarchical build-up of haloes. It remains to be understood why the evolved substructure adds up to the actual universal power law of v(f) ∝f−5/2. It seems that this behaviour results from the hierarchical clustering process and is not a general result of violent relaxation.