We report on the results of a direct N-body simulation of a star cluster that started with N = 200 000, comprising 195 000 single stars and 5000 primordial binaries. The code used for the simulation includes stellar evolution, binary evolution, an external tidal field and the effects of two-body relaxation. The model cluster is evolved to 12 Gyr, losing more than 80 per cent of its stars in the process. It reaches the end of the main core-collapse phase at 10.5 Gyr and experiences core oscillations from that point onwards -direct numerical confirmation of this phenomenon. However, we find that after a further 1 Gyr the core oscillations are halted by the ejection of a massive binary comprised of two black holes from the core, producing a core that shows no signature of the prior core-collapse. We also show that the results of previous studies with N ranging from 500 to 100 000 scale well to this new model with larger N. In particular, the time-scale to core collapse (in units of the relaxation time-scale), mass segregation, velocity dispersion and the energies of the binary population all show similar behaviour at different N.Key words: methods: numerical -binaries: close -stars: evolution -stars: kinematics and dynamics -globular clusters: general -galaxies: star clusters: general.
I N T RO D U C T I O NThe increasing ability of the direct N-body method to provide reliable models of the dynamical evolution of star clusters has closely mirrored increases in computing power (Heggie 2011). The community has progressed from small-N models performed on workstations (e.g. In this paper we present an N-body simulation of star cluster evolution that begins with N = 200 000 stars and binaries. This extends the N parameter space covered by direct N-body models and performs two important functions. First it provides a new calibration point for the MC method -this statistical method is increasingly valid for increasing N, so calibrations at higher N are more reliable. It also allows us to further develop our theoretical understanding of star cluster evolution and investigate how well inferences drawn from models of smaller N scale to larger values. The latter is the focus of this paper. A good example of the small-N models that we wish to compare with is the comprehensive study of star cluster evolution presented by Giersz & Heggie (1997) using models that included a mass function, stellar evolution and the tidal field of a point mass galaxy, albeit starting with 500 stars instead of 200 000. More recent examples for comparison include Baumgardt & Makino (2003) and Küpper, Kroupa & Baumgardt (2008). We were also motivated to produce a model that exhibited core collapse close to a Hubble time without dissolving by that time. What we find when interpreting this model is that much of the behaviour reported previously for smaller N-body models stands up well in comparison but that the actions of a binary comprised of two black holes (BHs) provide a late twist to the evolution of the cluster core.In Section 2 we describ...