We study the dynamical evolution of globular clusters using our H\'enon-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 interactions 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 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 stellar 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.

The planet formation process and subsequent planet migration may lead to configurations resulting in strong dynamical interactions among the various planets. Well-studied possible outcomes include collisions between planets, scattering events that eject one or more of the planets, and a collision of one or more of the planets with the parent star. In this work we consider one other possibility that has seemingly been overlooked in the various scattering calculations presented in the literature: the tidal capture of two planets which leads to the formation of a binary planet (or binary brown dwarf) in orbit about the parent star. We carry out extensive numerical simulations of such dynamical and tidal interactions to explore the parameter space for the formation of such binary planets. We show that tidal formation of binary planets is possible for typical planet masses and distances from the host star. The detection (or lack thereof) of planet-planet binaries can thus be used to constrain the properties of planetary systems, including their mutual spacing during formation, and the fraction of close planets in very eccentric orbits which are believed to form by a closely related process.

Inspired by the recent identification of the first candidate BH-WD X-ray binaries, where the compact accretors may be stellar-mass black hole candidates in extragalactic globular clusters, we explore how such binaries could be formed in a dynamical environment. We provide analyses of the formation rates via well known formation channels like binary exchange and physical collisions and propose that the only possibility to form BH-WD binaries is via coupling these usual formation channels with subsequent hardening and/or triple formation. Indeed, we find that the most important mechanism to make a BH-WD X-ray binary from an initially dynamically formed BH-WD binary is triple induced mass transfer via the Kozai mechanism. Even using the most optimistic estimates for the formation rates, we cannot match the observationally inferred production rates if black holes undergo significant evaporation from the cluster or form a completely detached subcluster of black holes. We estimate that at least 1% of all formed black holes, or presumably 10% of the black holes present in the core now, must be involved in interactions with the rest of the core stellar population.

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 evolution 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.

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 with 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.

In this paper we address the consumption of energy by astronomers while performing their professional duties. Although we find that astronomy uses a negligible fraction of the US energy budget, the rate at which energy is consumed by an average astronomer is similar to that of a typical high-flying businessperson. We review some of the ways in which astronomers are already acting to reduce their energy consumption. In the coming decades, all citizens will have to reduce their energy consumption to conserve fossil fuel reserves and to help avert a potentially catastrophic change in the Earth's climate. The challenges are the same for astronomers as they are for everyone: decreasing the distances we travel and investing in energy-efficient infrastructure. The high profile of astronomy in the media, and the great public interest in our field, can play a role in promoting energy-awareness to the wider population. Our specific recommendations are therefore to 1) reduce travel when possible, through efficient meeting organization, and by investing in high-bandwidth video conference facilities and virtual-world software, 2) create energy-efficient observatories, computing centers and workplaces, powered by sustainable energy resources, and 3) actively publicize these pursuits.

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 treatments 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 present the first results of a series of Monte-Carlo simulations investigating the imprint of a central black hole on the core structure of a globular cluster. We investigate the three-dimensional and the projected density profile of the inner regions of idealized as well as more realistic globular cluster models, taking into account a stellar mass spectrum, 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 simulations published previously in the literature.

In this paper, the second of a series, we study the stellar dynamical and evolutionary processes leading to the formation of compact binaries containing neutron stars (NSs) in dense globular clusters (GCs). For this study, 70 dense clusters were simulated independently, with a total stellar mass ~2x10^7Msun, exceeding the total mass of all dense GCs in our Galaxy. We find that, in order to reproduce the empirically derived formation rate of low-mass X-ray binaries (LMXBs), we must assume that NSs can be formed via electron-capture supernovae (ECS) with typical natal kicks smaller than in core-collapse supernovae. Our results explain the observed dependence of the number of LMXBs on ``collision number'' as well as the large scatter observed between different GCs. We predict that the number of quiescent LMXBs in different GCs should not have a strong metallicity dependence. In our cluster model the following mass-gaining events create populations of MSPs that do not match the observations: (i) accretion during a common envelope event with a NS formed through ECS, and (ii) mass transfer (MT) from a WD donor. Some processes lead only to a mild recycling. In addition, for MSPs, we distinguish low-magnetic-field (long-lived) and high-magnetic-field (short-lived) populations. With this distinction and by considering only those mass-gaining events that appear to lead to NS recycling, we obtain good agreement of our models with the numbers and characteristics of observed MSPs in 47 Tuc and Terzan 5, as well as with the cumulative statistics for MSPs detected in GCs of different dynamical properties. We find that significant production of merging double NSs potentially detectable as short gamma-ray bursts occurs only in very dense, most likely core-collapsed GCs. (abridged)

We consider the formation of double white dwarfs (DWDs) through dynamical interactions in globular clusters. Such interactions can give rise to eccentric DWDs, in contrast to the exclusively circular population expected to form in the Galactic disk. We show that for a 5-year Laser Interferometer Space Antenna (LISA) mission and distances as far as the Large Magellanic Cloud, multiple harmonics from eccentric DWDs can be detected at a signal-to-noise ratio higher than 8 for at least a handful of eccentric DWDs, given their formation rate and typical lifetimes estimated from current cluster simulations. Consequently the association of eccentricity with stellar-mass LISA sources does not uniquely involve neutron stars, as is usually assumed. Due to the difficulty of detecting (eccentric) DWDs with present and planned electromagnetic observatories, LISA could provide unique dynamical identifications of these systems in globular clusters.

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 driven 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 developed, with a clear correlation between the number of X-ray sources in a cluster, $N_X$, and the cluster's 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.

We systematically study the effects of collisions on the overall dynamical evolution of dense star clusters using Monte Carlo simulations over many relaxation times. We derive many observable properties of these clusters, including their core radii and the radial distribution of collision products. We also study different aspects of collisions in a cluster taking into account the shorter lifetimes of more massive stars, which has not been studied in detail before. Depending on the lifetimes of the significantly more massive collision products, observable properties of the cluster can be modified qualitatively; for example, even without binaries, core collapse can sometimes be avoided simply because of stellar collisions.

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 material 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.

We study the dynamical evolution of globular clusters containing populations of primordial binaries, using our newly updated Monte Carlo cluster evolution code with the inclusion of direct integration of binary scattering interactions. We describe the modifications we have made to the code, as well as improvements we have made to the core Monte Carlo method. We present several test calculations to verify the validity of the new code, and perform many comparisons with previous analytical and numerical work in the literature. We simulate the evolution of a large grid of models, with a wide range of initial cluster profiles, and with binary fractions ranging from 0 to 1, and compare with observations of Galactic globular clusters. We find that our code yields very good agreement with direct N-body simulations of clusters with primordial binaries, but yields some results that differ significantly from other approximate methods. Notably, the direct integration of binary interactions reduces their energy generation rate relative to the simple recipes used in Paper III, and yields smaller core radii. Our results for the structural parameters of clusters during the binary-burning phase are now in the tail of the range of parameters for observed clusters, implying that either clusters are born significantly more or less centrally concentrated than has been previously considered, or that there are additional physical processes beyond two-body relaxation and binary interactions that affect the structural characteristics of clusters.

This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discussion at the meeting. Our meeting, which had close to 300 registered participants, was broadly aimed at the large community of astronomers around the world working on the formation and evolution of compact objects and interacting binary systems in dense star clusters, such as globular clusters and galactic nuclei. Great advances have occurred in this field during the past few years, including the introduction of fundamentally new theoretical paradigms for the formation and evolution of compact objects in binaries as well as countless new discoveries by astronomers that have challenged many accepted models. Some of the highlights include: a nearly complete census of all the millisecond pulsars in 47 Tuc; first detections of many new radio pulsars in other clusters, particularly Terzan 5; detailed studies of X-ray binary populations and their luminosity functions in many galaxies and extragalactic globular clusters; increasing evidence for intermediate-mass black holes in clusters and greatly improved theoretical understanding of their possible formation processes.

Recent numerical simulations have suggested the possibility of forming double intermediate mass black holes (IMBHs) via the collisional runaway scenario in young dense star clusters. The two IMBHs formed would exchange into a common binary shortly after their birth, and quickly inspiral and merge. Since space-borne gravitational wave (GW) observatories such as LISA will be able to see the late phases of their inspiral out to several Gpc, and LIGO will be able to see the merger and ringdown out to similar distances, they represent potentially significant GW sources. In this Letter we estimate the rate at which LISA and LIGO will see their inspiral and merger in young star clusters, and discuss the information that can be extracted from the observations. We find that LISA will likely see tens of IMBH--IMBH inspirals per year, while advanced LIGO could see ~10 merger and ringdown events per year, with both rates strongly dependent on the distribution of cluster masses and densities.

In this paper, the first of a series, we study the stellar dynamical and evolutionary processes leading to the formation of compact binaries containing white dwarfs in dense globular clusters. We examine the processes leading to the creation of X-ray binaries such as cataclysmic variables and AM CVn systems. Using numerical simulations, we identify the dominant formation channels and we predict the expected numbers and characteristics of detectable systems, emphasizing how the cluster sources differ from the field population. We explore the dependence of formation rates on cluster properties and we explain in particular why the distribution of cataclysmic variables has only a weak dependence on cluster density. We also discuss the frequency of dwarf nova outbursts in globular clusters and their connection with moderately strong white dwarf magnetic fields. We examine the rate of Type Ia supernovae via both single and double degenerate channels in clusters and we argue that those rates may contribute to the total SN Ia rate in elliptical galaxies. Considering coalescing white dwarf binaries we discuss possible constraints on the common envelope evolution of their progenitors and we derive theoretical expectations for gravitational wave detection by LISA.

Recent theoretical work has solidified the viability of the collisional runaway scenario in young dense star clusters for the formation of very massive stars (VMSs), which may be precursors to intermediate-mass black holes (IMBHs). We present first results from a numerical study of the collisional runaway process in dense star clusters containing primordial binaries. Stellar collisions during binary scattering encounters offer an alternate channel for runaway growth, somewhat independent of direct collisions between single stars. We find that clusters with binary fractions >~10% yield two VMSs via collisional runaways, presenting the exotic possibility of forming IMBH--IMBH binaries in star clusters. We discuss the implications for gravitational wave observations, and the impact on cluster structure.

Observations and theoretical work suggest that globular clusters may be born with initially very large binary fractions. We present first results from our newly modified Monte-Carlo cluster evolution code, which treats binary interactions exactly via direct N-body integration. It is shown that binary scattering interactions generate significantly less energy than predicted by the recipes that have been used in the past to model them in approximate cluster evolution methods. The new result that the cores of globular clusters in the long-lived binary-burning phase are smaller than previously predicted weakens the agreement with observations, thus implying that more than simply stellar dynamics is at work in shaping the globular clusters we observe today.

We study dynamical interactions of star--planet binaries with other single stars. We derive analytical cross sections for all possible outcomes, and confirm them with numerical scattering experiments. We find that a wide mass ratio in the binary introduces a region in parameter space that is inaccessible to comparable-mass systems, in which the nature of the dynamical interaction is fundamentally different from what has traditionally been considered in the literature on binary scattering. We study the properties of the planetary systems that result from the scattering interactions for all regions of parameter space, paying particular attention to the location of the "hard--soft" boundary. The structure of the parameter space turns out to be significantly richer than a simple statement of the location of the "hard--soft" boundary would imply. We consider the implications of our findings, calculating characteristic lifetimes for planetary systems in dense stellar environments, and applying the results to previous analytical studies, as well as past and future observations. Recognizing that the system PSR B1620-26 in the globular cluster M4 lies in the "new" region of parameter space, we perform a detailed analysis quantifying the likelihood of different scenarios in forming the system we see today.

We model the dynamical evolution of primordial black holes (BHs) in dense star clusters using a simplified treatment of stellar dynamics in which the BHs are assumed to remain concentrated in an inner core, completely decoupled from the background stars. Dynamical interactions involving BH binaries are computed exactly and are generated according to a Monte Carlo prescription. Recoil and ejections lead to complete evaporation of the BH core on a timescale ~10^9 yr for typical globular cluster parameters. Orbital decay driven by gravitational radiation can make binaries merge and, in some cases, successive mergers can lead to significant BH growth. Our highly simplified treatment of the cluster dynamics allows us to study a large number of models and to compute statistical distributions of outcomes, such as the probability of massive BH growth and retention in a cluster. We find that, in most models, there is a significant probability (~20-80%) of BH growth with final masses > 100 M_{\sun}. In one case, a final BH formed with mass ~ 620 M_{\sun}. However, if the typical merger recoil speed (due to asymmetric emission of gravitational radiation) significantly exceeds the cluster escape speed, no growth ever occurs. Independent of the recoil speed, we find that BH-BH mergers enhanced by dynamical interactions in cluster cores present an important source of gravitational waves for ground--based laser interferometers. Under optimistic conditions, the total rate of detections by Advanced LIGO, could be as high as a few tens of events per year from inspiraling BHs from clusters.

While many observed ultra-luminous X-ray sources (ULXs, Lx > 10^39 erg s^-1) could be extragalactic X-ray binaries (XRBs) emitting close to the Eddington limit, the highest-luminosity ULXs (Lx > 3x10^39 erg s^-1) exceed the isotropic Eddington luminosity for even high-stellar-mass accreting black hole XRBs. It has been suggested that these highest-luminosity ULXs may contain accreting intermediate-mass black hole (IMBH) binaries. We consider this hypothesis for dense, young (about 100 Myr) stellar clusters where we assume that a 50-500 solar mass central IMBH has formed through runaway growth of a massive star. Using numerical simulations of the dynamics and evolution of the central black hole's captured companions, we obtain estimates of the incidence of mass transfer phases and possible ULX activity throughout the IMBH's evolutionary history. We find that, although it is common for the central black hole to acquire binary companions, there is a very low probability that these interacting binaries will become observable ULX sources.

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code (StarTrack), and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing 3-body and 4-body interactions (Fewbody). We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster like 47 Tuc, assuming an initial binary fraction of 100%, is only about 5-10%. We show that this is in good agreement with recent HST observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modeling of globular clusters, suggesting that ``realistic models'' should incorporate much larger initial binary fractions than has usually been done in the past.

Physical collisions between stars occur frequently in dense star clusters, either via close encounters between two single stars, or during strong dynamical interactions involving binary stars. Here we study stellar collisions that occur during binary-single and binary-binary interactions, by performing numerical scattering experiments. Our results include cross sections, branching ratios, and sample distributions of parameters for various outcomes. For interactions of hard binaries containing main-sequence stars, we find that the normalized cross section for at least one collision to occur (between any two of the four stars involved) is essentially unity, and that the probability of collisions involving more than two stars is significant. Hydrodynamic calculations have shown that the effective radius of a collision product can be 2-30 times larger than the normal main-sequence radius for a star of the same total mass. We study the effect of this expansion, and find that it increases the probability of further collisions considerably. We discuss these results in the context of recent observations of blue stragglers in globular clusters with masses exceeding twice the main-sequence turnoff mass. We also present Fewbody, a new, freely available numerical toolkit for simulating small-N gravitational dynamics that is particularly suited to performing scattering experiments.

We investigate the dynamical formation and evolution of binaries containing neutron stars in dense globular clusters. Our numerical simulations combine a simple Monte Carlo prescription for stellar dynamics, a sophisticated binary population synthesis code, and a small-N-body integrator for computing 3-body and 4-body interactions. Our results suggest that there is no ``retention problem,'' i.e., that, under standard assumptions, globular clusters can retain enough neutron stars to produce the observed numbers of millisecond pulsars. We also identify the dominant evolutionary and dynamical channels through which globular clusters produce their two main types of binary millisecond pulsars

We study the evolution of binary stars in globular clusters using a novel approach combining a state-of-the-art population synthesis code with a simple treatment of dynamical interactions in the dense cluster core. We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster like 47 Tuc, assuming an initial binary fraction of 100%, is only about 5%. We show that this is in good agreement with recent HST observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modeling of globular clusters, suggesting that ``realistic models'' should incorporate much larger initial binary fractions than has usually been done in the past.

The recent emergence of a new class of accretion-powered, transient, millisecond X-ray pulsars presents some difficulties for the conventional picture of accretion onto rapidly rotating magnetized neutron stars and their spin behavior during outbursts. In particular, it is unclear from the standard paradigm how these systems manage to accrete over such a wide range in dM/dt (i.e., >~ a factor of 50), and why the neutron stars exhibit a high rate of spindown in at least a number of cases. Following up on prior suggestions, we propose that `fast' X-ray pulsars can continue to accrete, and that their accretion disks terminate at approximately the corotation radius. We demonstrate the existence of such disk solutions by modifying the Shakura-Sunyaev equations with a simple magnetic torque prescription. The solutions are completely analytic, and have the same dependence on dM/dt and alpha (the viscosity parameter) as the original Shakura-Sunyaev solutions; but, the radial profiles can be considerably modified, depending on the degree of fastness. We apply these results to compute the torques expected during the outbursts of the transient millisecond pulsars, and find that we can explain the large spindown rates that are observed for quite plausible surface magnetic fields of ~10^9 G.

We study the dynamical evolution of globular clusters using our 2D Monte Carlo code with the inclusion of primordial binary interactions for equal-mass stars. We use approximate analytical cross sections for energy generation from binary-binary and binary-single interactions. After a brief period of slight contraction or expansion of the core over the first few relaxation times, all clusters enter a much longer phase of stable "binary burning" lasting many tens of relaxation times. The structural parameters of our models during this phase match well those of most observed globular clusters. At the end of this phase, clusters that have survived tidal disruption undergo deep core collapse, followed by gravothermal oscillations. Our results clearly show that the presence of even a small fraction of binaries in a cluster is sufficient to support the core against collapse significantly beyond the normal core collapse time predicted without the presence of binaries. For tidally truncated systems, collapse is easily delayed sufficiently that the cluster will undergo complete tidal disruption before core collapse. As a first step toward the eventual goal of computing all interactions exactly using dynamical three- and four-body integration, we have incorporated an exact treatment of binary-single interactions in our code. We show that results using analytical cross sections are in good agreement with those using exact three-body integration, even for small binary fractions where binary-single interactions are energetically most important.

We present the results of a new study of mass segregation in two-component star clusters, based on a large number of numerical N-body simulations using our recently developed dynamical Monte Carlo code. Specifically, we follow the dynamical evolution of clusters containing stars with individual masses m_1 as well as a tracer population of objects with individual masses m_2=\mu m_1, using N=10^5 total stars. For heavy tracers, which could represent stellar remnants such as neutron stars or black holes in a globular cluster, we characterize in a variety of ways the tendency for these objects to concentrate in or near the cluster core. In agreement with simple theoretical arguments, we find that the characteristic time for this mass segregation process varies as 1/\mu. For models with very light tracers (\mu <~ 10^-2), which could represent free-floating planets or brown dwarfs, we find the expected depletion of light objects in the cluster core, but also sometimes a significant enhancement in the halo. Using these results we estimate the optical depth to gravitational microlensing by planetary mass objects or brown dwarfs in typical globular clusters. For some initial conditions, the optical depth in the halo due to very low-mass objects could be much greater than that of luminous stars. If we apply our results to M22, using the recent null detection of Sahu, Anderson, & King (2001), we find an upper limit of ~25% at the 63% confidence level for the current mass fraction of M22 in the form of very low-mass objects.

We estimate the rate of collisions between stars and free-floating planets (FFPs) in globular clusters, in particular the collision of FFPs with red giant branch (RGB) stars. Recent dynamical simulations imply that the density of such objects could exceed million per cubic parsec near the cores of rich globular clusters. We show that in these clusters 5-10 per cents of all RGB stars near the core would suffer a collision with a FFP, and that such a collision can spin up the RGB star's envelope by an order of magnitude. In turn, the higher rotation rates may lead to enhanced mass-loss rates on the RGB, which could result in bluer horizontal branch (HB) stars. Hence, it is plausible that the presence of a large population of FFPs in a globular cluster can influence the distribution of stars on the HB of that cluster to a detectable degree.

We summarize the results of recent theoretical work on the dynamical evolution of globular clusters containing primordial binaries. Even a very small initial binary fraction (e.g., 10%) can play a key role in supporting a cluster against gravothermal collapse for many relaxation times. Inelastic encounters between binaries and single stars or other binaries provide a very significant energy source for the cluster. These dynamical interactions also lead to the production of large numbers of exotic systems such as ultracompact X-ray binaries, recycled radio pulsars, double degenerate systems, and blue stragglers. Our work is based on a new parallel supercomputer code implementing Henon's Monte Carlo method for simulating the dynamical evolution of dense stellar systems in the Fokker-Planck approximation. This new code allows us to calculate very accurately the evolution of a cluster containing a realistic number of stars (N ~ 10^5 - 10^6) in typically a few hours to a few days of computing time. The discrete, star-by-star representation of the cluster in the simulation makes it possible to treat naturally a number of important processes, including single and binary star evolution, all dynamical interactions of single stars and binaries, and tidal interactions with the Galaxy.

We improve and extend Shapiro's model of a relativistic, compact object which is stable in isolation but is driven dynamically unstable by the tidal field of a binary companion. Our compact object consists of a dense swarm of test particles moving in randomly-oriented, initially circular, relativistic orbits about a nonrotating black hole. The binary companion is a distant, slowly inspiraling point mass. The tidal field of the companion is treated as a small perturbation on the background Schwarzschild geometry near the hole; the resulting metric is determined by solving the perturbation equations of Regge and Wheeler and Zerilli in the quasi-static limit. The perturbed spacetime supports Bekenstein's conjecture that the horizon area of a near-equilibrium black hole is an adiabatic invariant. We follow the evolution of the system and confirm that gravitational collapse can be induced in a compact collisionless cluster by the tidal field of a binary companion.