STATISTICS OF N-BODY SIMULATIONS .2. EQUAL MASSES AFTER CORE COLLAPSE

Citation
M. Giersz et Dc. Heggie, STATISTICS OF N-BODY SIMULATIONS .2. EQUAL MASSES AFTER CORE COLLAPSE, Monthly Notices of the Royal Astronomical Society, 270(2), 1994, pp. 298-324
Citations number
29
Categorie Soggetti
Astronomy & Astrophysics
ISSN journal
00358711
Volume
270
Issue
2
Year of publication
1994
Pages
298 - 324
Database
ISI
SICI code
0035-8711(1994)270:2<298:SONS.E>2.0.ZU;2-0
Abstract
This paper presents and analyses statistical results from a large numb er of N-body simulations of isolated systems with equal masses, in whi ch 250 less than or equal to N less than or equal to 2000. It concentr ates on the phase starting around the end of core collapse. Binaries p lay a crucial role, and we find that the total energy of bound pairs i s in line with theoretical expectations. Interpretation of the total n umber is complicated by the presence of a number of binaries on the ha rd/soft threshold. Interactions of hard binaries are consistent with t he Spitzer cross-section. The spatial evolution of the half-mass radiu s after core collapse nearly follows classical theory, and, by compari son with Fokker-Planck and gas models, allows a redetermination of the effective thermal conductivity and the argument of the Coulomb logari thm in the expression for the relaxation time. The evolution of the in ner parts of the system around the time of core bounce is consistent w ith these simplified models, provided that the continuous production o f energy, as is usually assumed, is replaced by a model of stochastic energy production. Similarly, postcollapse evolution of the core requi res a modest recalibration of the coefficient of energy generation, es pecially for small N. These remarks refer to the behaviour averaged ov er many models; individual cases show alternate and irregular phases o f expansion and recollapse. The distributions of velocity dispersion a nd anisotropy become remarkably homologous soon after core bounce. The bound mass of the systems very nearly follows a power law with time. A small number of escapers, presumed to be those associated with binar y activity, dominate the energy that is carried off: the distribution of energies of escapers changes abruptly at the end of core collapse. The 'internal' energy of escaping binaries is consistent with theoreti cal expectations, and again supports Spitzer's reaction cross-section for hard binaries.