Using our recently improved Monte Carlo evolution code, we study the evolution of the binary fraction in globular clusters. In agreement with previous N-body simulations, we find generally that the hard binary fraction in the core tends to increase w
ith time over a range of initial cluster central densities for initial binary fractions <~ 90%. The dominant processes driving the evolution of the core binary fraction are mass segregation of binaries into the cluster core and preferential destruction of binaries there. On a global scale, these effects and the preferential tidal stripping of single stars tend to roughly balance, leading to overall cluster binary fractions that are roughly constant with time. Our findings suggest that the current hard binary fraction near the half-mass radius is a good indicator of the hard primordial binary fraction. However, the relationship between the true binary fraction and the fraction of main-sequence stars in binaries (which is typically what observers measure) is non-linear and rather complicated. We also consider the importance of soft binaries, which not only modify the evolution of the binary fraction, but can drastically change the evolution of the cluster as a whole. Finally, we describe in some detail the recent addition of single and binary stellar evolution to our cluster evolution code.
Recent observations of the white dwarf (WD) populations in the Galactic globular cluster NGC 6397 suggest that WDs receive a kick of a few km/s shortly before they are born. Using our Monte Carlo cluster evolution code, which includes accurate treatm
ents of all relevant physical processes operating in globular clusters, we study the effects of the kicks on their host cluster and on the WD population itself. We find that in clusters whose velocity dispersion is comparable to the kick speed, WD kicks are a significant energy source for the cluster, prolonging the initial cluster core contraction phase significantly so that at late times the cluster core to half-mass radius ratio is a factor of up to ~ 10 larger than in the no-kick case. WD kicks thus represent a possible resolution of the large discrepancy between observed and theoretically predicted values of this key structural parameter. Our modeling also reproduces the observed trend for younger WDs to be more extended in their radial distribution in the cluster than older WDs.
We summarize and discuss recent work (Fregeau 2007) that presents the confluence of three results suggesting that most Galactic globular clusters are still in the process of core contraction, and have not yet reached the thermal equilibrium phase dri
ven by binary scattering interactions: that 1) the three clusters that appear to be overabundant in X-ray binaries per unit encounter frequency are observationally classified as core-collapsed, 2) recent numerical simulations of cluster evolution with primordial binaries show that structural parameters of clusters in the binary-burning phase agree only with core-collapsed clusters, and 3) a cluster in the binary-burning phase for the last few Gyr should have about 5 times more dynamically formed X-ray sources than if it were in the core contraction phase for the same time.
It has been known for over 30 years that Galactic globular clusters (GCs) are overabundant by orders of magnitude in bright X-ray sources per unit mass relative to the disk population. Recently a quantitative understanding of this phenomenon has deve
loped, with a clear correlation between the number of X-ray sources in a cluster, $N_X$, and the clusters encounter frequency, $Gamma$, becoming apparent. We derive a refined version of $Gamma$ that incorporates the finite lifetime of X-ray sources and the dynamical evolution of clusters. With it we find we are able to explain the few clusters that lie off the $N_X$--$Gamma$ correlation, and resolve the discrepancy between observed GC core radii and the values predicted by theory. Our results suggest that most GCs are still in the process of core contraction and have not yet reached the thermal equilibrium phase driven by binary scattering interactions.
In this brief proceedings article I summarize the review talk I gave at the IAU 246 meeting in Capri, Italy, glossing over the well-known results from the literature, but paying particular attention to new, previously unpublished material. This new m
aterial includes a careful comparison of the apparently contradictory results of two independent methods used to simulate the evolution of binary populations in dense stellar systems (the direct N-body method of Hurley, et al. 2007 and the approximate Monte Carlo method of Ivanova, et al. 2005), that shows that the two methods may not actually yield contradictory results, and suggests future work to more directly compare the two methods.
