Star clusters with primordial binaries – I. Dynamical evolution of isolated models

In order to interpret the results of complex realistic star cluster simulations, which rely on many simplifying approximations and assumptions, it is essential to study the behaviour of even more idealized models, which can highlight the essential physical effects and are amenable to more exact methods. With this aim, we present the results of N -body calculations of the evolution of equal-mass models, starting with primordial binary fractions of 0–100 per cent, with values of N ranging from 256 to 16 384. This allows us to extrapolate the main features of the evolution to systems comparable in particle number with globular clusters.
In this range, we find that the steady-state 'deuterium main sequence' is characterized by a ratio of the core radius to half-mass radius that follows qualitatively the analytical estimate by Vesperini & Chernoff, although the N dependence is steeper than expected. Interestingly, for an initial binary fraction f greater than 10 per cent, the binary heating in the core during the post-collapse phase almost saturates (becoming nearly independent of f ), and so little variation in the structural properties is observed. Thus, although we observe a significantly lower binary abundance in the core with respect to the Fokker–Planck simulations by Gao et al., this is of little dynamical consequence.
At variance with the study of Gao et al., we see no sign of gravothermal oscillations before 150 half-mass relaxation times. At later times, however, oscillations become prominent. We demonstrate the gravothermal nature of these oscillations.     

On the onset of runaway stellar collisions in dense star clusters – I. Dynamics of the first collision

We study the circumstances under which first collisions occur in young and dense star clusters. The initial conditions for our direct N -body simulations are chosen such that the clusters experience core collapse within a few million years, before the most massive stars have left the main sequence. It turns out that the first collision is typically driven by the most massive stars in the cluster. Upon arrival in the cluster core, by dynamical friction, massive stars tend to form binaries. The enhanced cross-section of the binary compared to a single star causes other stars to engage the binary. A collision between one of the binary components and the incoming third star is then mediated by the encounters between the binary and other cluster members. Due to the geometry of the binary–single star engagement the relative velocity at the moment of impact is substantially different than in a two-body encounter. This may have profound consequences for the further evolution of the collision product.     

Ratios of star cluster core and half-mass radii: a cautionary note on intermediate-mass black holes in star clusters

There is currently much interest in the possible presence of intermediate-mass black holes (IMBHs) in the cores of globular clusters (GCs). Based on theoretical arguments and simulation results it has previously been suggested that a large core radius – or particularly a large ratio of the core radius to half-mass radius – is a promising indicator for finding such a black hole (BH) in a star cluster. In this study N -body models of 100 000 stars with and without primordial binaries are used to investigate the long-term structural evolution of star clusters. Importantly, the simulation data are analysed using the same processes by which structural parameters are extracted from observed star clusters. This gives a ratio of the core and half-mass (or half-light) radii that are directly comparable to the Galactic GC sample. As a result, it is shown that the ratios observed for the bulk of this sample can be explained without the need for an IMBH. Furthermore, it is possible that clusters with large core to half-light radius ratios harbour a BH binary (comprising stellar mass BHs) rather than a single massive BH. This work does not rule out the existence of IMBHs in the cores of at least some star clusters.     

A stochastic Monte Carlo approach to model real star cluster evolution – II. Self-consistent models and primordial binaries

The new approach outlined in Paper I 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 per cent of the initial number of stars in the cluster mass). 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. The cluster evolution is about three times slower than found by Gao et al. Other features are rather comparable. At every moment we are able to show the individual distribution of binaries in the cluster.     

Limits on initial mass segregation in young clusters

Mass segregation is observed in many star clusters, including several that are less than a few Myr old. Time-scale arguments are frequently used to argue that these clusters must be displaying primordial segregation, because they are too young to be dynamically relaxed. Looking at this argument from the other side, the youth of these clusters and the limited time available to mix spatially distinct populations of stars can provide constraints on the amount of initial segregation that is consistent with current observations. We present n -body experiments testing this idea, and discuss the implications of our results for theories of star formation. For system ages less than a few crossing times, we show that star formation scenarios predicting general primordial mass segregation are inconsistent with observed segregation levels.     

Mass transfer from a giant star to a main-sequence companion and its contribution to long-orbital-period blue stragglers

Binary population synthesis shows that mass transfer from a giant star to a main-sequence (MS) companion may account for some observed long-orbital-period blue stragglers. However, little attention is paid to this blue straggler formation scenario as dynamical instability often happens when the mass donor is a giant star. In this paper, we have studied the critical mass ratio, q _c, for dynamically stable mass transfer from a giant star to a MS companion using detailed evolution calculations. The results show that a more evolved star is generally less stable for Roche lobe overflow. Meanwhile,   q _c  almost linearly increases with the amount of the mass and angular momentum lost during mass transfer, but has little dependance on stellar wind. To conveniently use the result, we give a fit of q _c as a function of the stellar radius at the onset of Roche lobe overflow and of the mass-transfer efficiency during the Roche lobe overflow.
To examine the formation of blue stragglers from mass transfer between giants and MS stars, we have performed Monte Carlo simulations with various q _c. The simulations show that some binaries with the mass donor on the first giant branch may contribute to blue stragglers with q _c obtained in this paper but will not from previous q _c. Meanwhile, from our q _c, blue stragglers from the mass transfer between an asymptotic giant branch star and a MS companion may be more numerous and have a wider range of orbital periods than those from the other q _c.     

Effects of external tidal field on the evolution of the outer regions of multi-mass star clusters

We present N -body simulations (including an initial mass function) of globular clusters in the Galaxy in order to study effects of the tidal field systematically on the properties of the outer parts of globular clusters. Using nbody6 , which correctly takes into account the two-body relaxation, we investigate the development of tidal tails of globular clusters in the Galactic tidal field. For simplicity, we have employed only the spherical components (bulge and halo) of the Galaxy, and ignored the effects of stellar evolution which could have been important in the very early phase of the cluster evolution. The total number of stars in our simulations is about 20 000, which is much smaller than the realistic number of stars. All simulations had been done for several orbital periods in order to understand the development of the tidal tails. In our scaled-down models, the relaxation time is sufficiently short to show the mass segregation effect, but we did not go far enough to see the core collapse, and the fraction of stars lost from the cluster at the end of the simulations is only ∼10 per cent. The radial distribution of extra-tidal stars can be described by a power law with a slope around −3 in surface density. The directions of tidal tails are determined by the orbits and locations of the clusters. We find that the length of tidal tails increases towards the apogalacticon and decreases towards the perigalacticon. This is an anti-correlation with the strength of the tidal field, caused by the fact that the time-scale for the stars to respond to the potential is similar to the orbital time-scale of the cluster. The escape of stars in the tidal tails towards the pericentre could be another reason for the decrease of the length of tidal tails. We find that the rotational angular velocity of tidally induced clusters shows quite different behaviour from that of initially rotating clusters.