The origin and formation of cuspy density profiles through violent relaxation of stellar sy.pdf

a r X i v : a s t r o - p h / 9 9 0 6 3 5 4 v 1 2 2 J u n 1 9 9 9 Mon. Not. R. Astron. Soc. 000, 1–6 (1999) Printed 5 February 2008 (MN LaTEX style file v1.4) The origin and formation of cuspy density profiles through violent relaxation of stellar systems S. Hozumi,1? A. Burkert2? and T. Fujiwara3? 1Faculty of Education, Shiga University, 2-5-1 Hiratsu, Otsu, Shiga 520-0862, Japan 2Max-Planck-Institut fu?r Astronomie, Ko?nigstuhl 17, D-69117 Heidelberg, Germany 3Kyoto City University of Arts, Nishikyo-ku, Kyoto 610-1197, Japan Accepted 1999 xxx xx. Received 1999 January 4 ABSTRACT It is shown that the cuspy density distributions observed in the cores of elliptical galaxies can be realized by dissipationless gravitational collapse. The initial models consist of power-law density spheres such as ρ ∝ r?1 with anisotropic velocity disper- sions. Collapse simulations are carried out by integrating the collisionless Boltzmann equation directly, on the assumption of spherical symmetry. From the results obtained, the extent of constant density cores, formed through violent relaxation, decreases as the velocity anisotropy increases radially, and practically disappears for extremely ra- dially anisotropic models. As a result, the relaxed density distributions become more cuspy with increasing radial velocity anisotropy. It is thus concluded that the velocity anisotropy could be a key ingredient for the formation of density cusps in a dissipa- tionless collapse picture. The velocity dispersions increase with radius in the cores according to the nearly power-law density distributions. The power-law index, n, of the density profiles, defined as ρ ∝ r?n, changes from n ≈ 2.1 at intermediate radii, to a shallower power than n ≈ 2.1 toward the centre. This density bend can be explained from our postulated local phase-space constraint that the phase-space density acces- sible to the relaxed state is determined at each radius by the maximum phase-space density of the initial sta