SCYON Abstract

Received on February 21 2006

Star Clusters with Primordial Binaries: II. Dynamical Evolution of Models in a Tidal Field

AuthorsM. Trenti(1,2), D.C. Heggie(3), P. Hut(4)
Affiliation
(1) Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA
(2) Yukawa Institute for Theoretical Physics, Kyoto University, 606-8502 Kyoto, Japan
(3) School of Mathematics, University of Edinburgh, King's Buildings, Edinburgh EH9 3JZ, Scotland, U.K.
(4) Institute for Advanced Study, Princeton, NJ 08540, USA
Submitted toMonthly Notices of the Royal Astronomical Society
Contacttrenti@stsci.edu
URLhttp://arxiv.org/abs/astro-ph/0602409
Links

Abstract

[abridged] We extend our analysis of the dynamical evolution of simple star cluster models, in order to provide comparison standards that will aid in interpreting the results of more complex realistic simulations. We augment our previous primordial-binary simulations by introducing a tidal field, and starting with King models of different central concentrations. We present the results of N-body calculations of the evolution of equal-mass models, starting with primordial binary fractions of 0 - 100 %, and N values from 512 to 16384. We also attempt to extrapolate some of our results to the larger number of particles that are necessary to model globular clusters. We characterize the steady-state "deuterium main sequence" phase in which primordial binaries are depleted in the core in the process of "gravitationally burning". In this phase we find that the ratio of the core to half-mass radius, rc/rh, is similar to that measured for isolated systems. In addition to the generation of energy due to hardening and depletion of the primordial binary population, the overall evolution of the star clusters is driven by a competing process: the tidal disruption of the system. We find that the depletion of primordial binaries before tidal dissolution of the system is possible only if the initial number is below 0.05 N, in the case of a King model with W0=7 and N=4096 (which is one of our longest living models). We compare our findings, obtained by means of direct N-body simulations but scaled, where possible, to larger N, with similar studies carried out by means of Monte Carlo methods.