A stochastic Monte Carlo approach to model real star cluster evolution; 2, Self-consistent models and primordial binaries

The new approach outlined in Paper I (Spurzem \& Giersz 1996) to follow the individual formation and evolution of binaries in an evolving, equal point-mass star cluster is extended for the self-consistent treatment of relaxation and close three- and four-body encounters for many binaries (typically a few percent of the initial number of stars in the cluster). The distribution of single stars is treated as a conducting gas sphere with a standard anisotropic gaseous model. A Monte Carlo technique is used to model the motion of binaries, their formation and subsequent hardening by close encounters, and their relaxation (dynamical friction) with single stars and other binaries. The results are a further approach towards a realistic model of globular clusters with primordial binaries without using special hardware. We present, as our main result, the self-consistent evolution of a cluster consisting of 300.000 equal point-mass stars, plus 30.000 equal mass binaries over several hundred half-mass relaxation times, well into the phase where most of the binaries have been dissolved and evacuated from the core. In a self-consistent model it is the first time that such a realistically large number of binaries is evolving in a cluster with an even ten times larger number of single stars. Due to the Monte Carlo treatment of the binaries we can at every moment analyze their external and internal parameters in the cluster as in an N-body simulation