We study the dynamical evolution of globular clusters using our Henon-type Monte Carlo code for stellar dynamics including all relevant physics such as two-body relaxation, single and binary stellar evolution, Galactic tidal stripping, and strong int
eractions such as physical collisions and binary mediated scattering. We compute a large database of several hundred models starting from broad ranges of initial conditions guided by observations of young and massive star clusters. We show that these initial conditions very naturally lead to present day clusters with properties including the central density, core radius, half-light radius, half-mass relaxation time, and cluster mass, that match well with those of the old Galactic globular clusters. In particular, we can naturally reproduce the bimodal distribution in observed core radii separating the core-collapsed vs the non core-collapsed clusters. We see that the core-collapsed clusters are those that have reached or are about to reach the equilibrium binary burning phase. The non core-collapsed clusters are still undergoing gravo-thermal contraction.
We study the dynamical evolution of globular clusters containing primordial binaries, including full single and binary stellar evolution using our Monte Carlo cluster evolution code updated with an adaptation of the single and binary stellar evolutio
n codes SSE/BSE from Hurley et. al (2000, 2002). We describe the modifications we have made to the code. We present several test calculations and comparisons with existing studies to illustrate the validity of the code. We show that our code finds very good agreement with direct N-body simulations including primordial binaries and stellar evolution. We find significant differences in the evolution of the global properties of the simulated clusters using stellar evolution compared to simulations without any stellar evolution. In particular, we find that the mass loss from stellar evolution acts as a significant energy production channel simply by reducing the total gravitational binding energy and can significantly prolong the initial core contraction phase before reaching the binary-burning quasi steady state of the cluster evolution as noticed in Paper IV. We simulate a large grid of clusters varying the initial cluster mass, binary fraction, and concentration and compare properties of the simulated clusters with those of the observed Galactic globular clusters (GGCs). We find that our simulated cluster properties agree well with the observed GGC properties. We explore in some detail qualitatively different clusters in different phases of their evolution, and construct synthetic Hertzprung-Russell diagrams for these clusters.
We present results of a series of Monte Carlo simulations investigating the imprint of a central intermediate-mass black hole (IMBH) on the structure of a globular cluster. We investigate the three-dimensional and projected density profiles, and stel
lar disruption rates for idealized as well as realistic cluster models, taking into account a stellar mass spectrum and stellar evolution, and allowing for a larger, more realistic, number of stars than was previously possible with direct N-body methods. We compare our results to other N-body and Fokker-Planck simulations published previously. We find, in general, very good agreement for the overall cluster structure and dynamical evolution between direct N-body simulations and our Monte Carlo simulations. Significant differences exist in the number of stars that are tidally disrupted by the IMBH, which is most likely an effect of the wandering motion of the IMBH, not included in the Monte Carlo scheme. These differences, however, are negligible for the final IMBH masses in realistic cluster models as the disruption rates are generally much lower than for single-mass clusters. As a direct comparison to observations we construct a detailed model for the cluster NGC 5694, which is known to possess a central surface brightness cusp consistent with the presence of an IMBH. We find that not only the inner slope but also the outer part of the surface brightness profile agree well with observations. However, there is only a slight preference for models harboring an IMBH compared to models without.
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.
