Recent numerical simulations of globular clusters (GCs) have shown that stellar-mass black holes (BHs) play a fundamental role in driving cluster evolution and shaping their present-day structure. Rapidly mass-segregating to the center of GCs, BHs act as a dynamical energy source via repeated super-elastic scattering, delaying onset of core collapse and limiting mass segregation for visible stars. While recent discoveries of BH candidates in Galactic and extragalactic GCs have further piqued interest in BH-mediated cluster dynamics, numerical models show that even if significant BH populations remain in today's GCs, they are typically in configurations that are not directly detectable. We demonstrated in Weatherford et al. (2018) that an anti-correlation between a suitable measure of mass segregation ($\Delta$) in observable stellar populations and the number of retained BHs in GC models can be applied to indirectly probe BH populations in real GCs. Here, we estimate the number and total mass of BHs presently retained in 50 Milky Way GCs from the ACS Globular Cluster Survey by measuring $\Delta$ between populations of main sequence stars, using correlations found between $\Delta$ and BH retention in the CMC Cluster Catalog models. We demonstrate that the range in $\Delta$'s distribution from our models matches that for observed GCs to a remarkable degree. Our results further provide the narrowest constraints to-date on the retained BH populations in the GCs analyzed. Of these 50 GCs, we identify NGCs 2808, 5927, 5986, 6101, and 6205 to presently contain especially large BH populations, each with total BH mass exceeding $10^3\,\rm{M_{\odot}}$.

Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae. However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called "upper-mass gap" would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through pre-collapse stellar collisions in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as $\sim20\%$ of all BH progenitors undergo one or more collisions prior to stellar collapse and up to $\sim1\%$ of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the "multiple-generation" merger channel described in other analyses. Finally, we describe how stellar collisions in clusters may provide a unique pathway to pair-instability supernovae and briefly discuss the expected rate of these events and other electromagnetic transients.

We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics.

Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt, despite their potential role as massive seeds for quasars and sources of tidal disruption events, ultra-luminous X-ray sources, dwarf galaxy feedback, and hypervelocity stars. Gravitational wave (GW) observations can help to find and confirm the existence of IMBHs. Current and upcoming detectors, such as LIGO, Virgo, KAGRA, LISA, ET, and DECIGO promise to identify the full range from stellar-mass to supermassive black holes (SMBHs). In this paper, we address the question of whether IMBHs can produce GWs in galactic nuclei. We consider the possibility that stellar black holes (SBHs) form bound systems and later coalesce with an IMBH through gravitational captures in the dense nucleus. We show that this mechanism is efficient for IMBH masses in the range $\sim 3\times 10^3$ M$_\odot$--$2\times 10^4$ M$_\odot$. We find that the typical distributions of peak frequencies and merger timescales depend mainly on the IMBH mass. In particular, the typical peak frequency is about $0.2\,$Hz, $0.1\,$Hz, $0.09\,$Hz, and $0.05\,$Hz for $M_{\rm IMBH}=5\times 10^3$ M$_\odot$, $8\times 10^3$ M$_\odot$, $1\times 10^4$ M$_\odot$, and $2\times 10^4$ M$_\odot$, respectively. Our results show that, at design sensitivity, both DECIGO and ET should be able to detect these IMBH--SBH mergers. Furthermore, most of the mergers will appear eccentric ($e \gtrsim 0.1$), providing an indication of their dynamical origin.

The origin of the black hole (BH) binary mergers observed by LIGO-Virgo is still uncertain, as are the boundaries of the stellar BH mass function. Stellar evolution models predict a dearth of BHs both at masses $\gtrsim 50$ M$_\odot$ and $\lesssim 5$ M$_\odot$, thus leaving low- and high-mass gaps in the BH mass function. A natural way to form BHs of these masses is through mergers of neutron stars (NSs; for the low-mass gap) or lower-mass BHs (for the high-mass gap); the low- or high-mass-gap BH produced as a merger product can then be detected by LIGO-Virgo if it merges again with a new companion. We show that the evolution of a 2+2 quadruple system can naturally lead to BH mergers with component masses in the low- or high-mass gaps. In our scenario, the BH in the mass gap originates from the merger of two NSs, or two BHs, in one of the two binaries and the merger product is imparted a recoil kick (from anisotropic gravitational wave emission), which triggers its interaction with the other binary component of the quadruple system. The outcome of this three-body interaction is usually a new eccentric compact binary containing the BH in the mass gap, which can then merge again. The merger rate is $\sim 10^{-7} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ and $\sim 10^{-3} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ for BHs in the low-mass and high-mass gap, respectively. As the sensitivity of gravitational wave detectors improves, tighter constraints will soon be placed on the stellar BH mass function.

We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters.

The first detection of gravitational waves from a neutron star - neutron star (NS-NS) merger, GW170817, and the increasing number of observations of short gamma-ray bursts (SGRBs) have greatly motivated studies of the origins of NS-NS and neutron star - black hole (NS-BH) binaries. We calculate the merger rates of NS-NS and NS-BH binaries from globular clusters (GCs) using realistic GC simulations with the \texttt{CMC} cluster catalog. We use a large sample of models with a range of initial numbers of stars, metallicities, virial radii and galactocentric distances, representative of the present-day Milky Way GCs, to quantify the inspiral times and volumetric merger rates as a function of redshift, both inside and ejected from clusters. We find that over the complete lifetime of most GCs, stellar BHs dominate the cluster cores and prevent the mass segregation of NSs, thereby reducing the dynamical interaction rates of NSs so that at most a few NS binary mergers are ever produced. We estimate the merger rate in the local universe to be $\sim\rm{0.02\,Gpc^{-3}\,yr^{-1}}$ for both NS-NS and NS-BH binaries, or a total of $\sim 0.04$~Gpc$^{-3}$~yr$^{-1}$ for both populations. These rates are about 5 orders of magnitude below the current empirical merger rate from LIGO/Virgo. We conclude that dynamical interactions in GCs do not play a significant role in enhancing the NS-NS and NS-BH merger rates.

The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA.

Among the growing list of confirmed exoplanets, the number of planets identified in dense star clusters remains sparse. Previous analyses have suggested this may be due in part to dynamical interactions that unbind large fractions of planets from their host stars, limiting the survival of planetary systems in clusters. Thus, alternative detection strategies may be necessary to study planets in clusters that may no longer be bound to a host star. Here, we use the cluster Monte Carlo code CMC to explore the evolution of planetary systems in dense star clusters. Depending on a number of initial conditions, we show that $10-50\%$ of primordial planetary systems are broken through dynamical encounters over a cluster's full lifetime, populating clusters with "free-floating" planets. Furthermore, a large number ($30-80\%$) of planets are ejected from their host cluster through strong dynamical encounters and/or tidal loss. Additionally, we show that planets naturally mix with stellar-mass black holes (BHs) in the central regions of their host cluster. As a consequence, up to a few hundreds of planets will be tidally disrupted through close passages of BHs. We show these BH-planet tidal disruption events (TDEs) occur in clusters at a rate of up to $10^{-5}\,\rm{yr}^{-1}$ in a Milky Way-type galaxy. We predict BH-planet TDEs should be detected by upcoming transient surveys such as LSST at a rate of roughly a few events per year. The observed rate of BH-planet TDEs would place new constraints upon the formation and survival of both planetary systems and BHs in dense star clusters. Additionally, depending on various assumptions including the initial number of planets and their orbital properties, we predict that typical globular clusters may contain a few dynamically-formed NS-planet systems at present as well as up to roughly 100 dynamically-formed WD-planet systems.

When two black holes merge in a dense star cluster, they form a new black hole with a well-defined mass and spin. If that "second-generation" black hole remains in the cluster, it will continue to participate in dynamical encounters, form binaries, and potentially merge again. Using a grid of 96 dynamical models of dense star clusters and a cosmological model of cluster formation, we explore the production of binary black hole mergers where at least one component of the binary was forged in a previous merger. We create four hypothetical universes where every black hole born in the collapse of a massive star has a dimensionless Kerr spin parameter, $\chi_{\rm birth}$, of 0.0, 0.1, 0.2, or 0.5. We show that if all stellar-born black holes are non-spinning ($\chi_{\rm birth} = 0.0$), then more than 10% of merging binary black holes from clusters have components formed from previous mergers, accounting for more than 20% of the mergers from globular clusters detectable by LIGO/Virgo. Furthermore, nearly 7% of detectable mergers would have a component with a mass $\gtrsim 55M_{\odot}$, placing them clearly in the mass "gap" region where black holes cannot form from isolated collapsing stars due to the pulsational-pair instability mechanism. On the other hand, if black holes are born spinning, then the contribution from these second-generation mergers decreases, making up as little as 1% of all detections from globular clusters when $\chi_{\rm birth} = 0.5$. We make quantitative predictions for the detected masses, mass ratios, and spin properties of first- and second-generation mergers from dense star clusters, and show how these distributions are highly sensitive to the birth spins of black holes.

As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, which leads to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward so-called "core collapse," the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the core-collapsed and non-core-collapsed clusters. This suggests an internal energy source is at work, delaying the onset of core collapse in many clusters. Primordial binary stars have been thought for a long time to provide this energy source, but recent analyses have cast doubt upon the corresponding "binary-burning" mechanism as a viable explanation. Over the past decade, a large amount of both observational and theoretical work has suggested that many stellar-mass black holes (BHs) are retained in typical clusters today and that they play a dynamically-significant role in these clusters throughout their entire lifetimes. Here we review our latest understanding of the formation and evolution of BH populations in GCs and demonstrate that, through their dynamical interaction with their host cluster, BHs can naturally explain the distinction between core-collapsed and non-core-collapsed clusters through a process we call "black hole burning."

In a dense stellar environment, such as the core of a globular cluster (GC), dynamical interactions with black holes (BHs) are expected to lead to a variety of astrophysical transients. Here we explore tidal disruption events (TDEs) of stars by stellar-mass BHs through collisions and close encounters. Using state-of-the-art $N$-body simulations, we show that these TDEs occur at significant rates throughout the evolution of typical GCs and we study how their relative rates relate to cluster parameters such mass and size. By incorporating a realistic cosmological model of GC formation, we predict a BH - main-sequence-star TDE rate of approximately $3\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ in the local universe ($z<0.1$) and a cosmological rate that peaks at roughly $25\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ for redshift 3. Furthermore, we show that the ejected mass associated with these TDEs could produce optical transients of luminosity $\sim 10^{41} - 10^{44}\rm\,erg\,s^{-1}$ with timescales of about a day to a month. These should be readily detectable by optical transient surveys such as the Zwicky Transient Facility. Finally, we comment briefly on BH - giant encounters and discuss how these events may contribute to the formation of BH - white-dwarf binaries.

Over a hundred millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. Here we study MSPs in GCs using realistic $N$-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae (including both accretion-induced and merger-induced collapse of white dwarfs) can be spun up through mass transfer to form MSPs. Both NS formation and spin-up through accretion are greatly enhanced through dynamical interaction processes. We find that our models for average GCs at the present day with masses $\approx 2 \times 10^5\,M_\odot$ can produce up to $10-20$ MSPs, while a very massive GC model with mass $\approx 10^6\,M_\odot$ can produce close to $100$. We show that the number of MSPs is anti-correlated with the total number of stellar-mass black holes (BHs) retained in the host cluster. The radial distributions are also affected: MSPs are more concentrated towards the center in a host cluster with a smaller number of retained BHs. As a result, the number of MSPs in a GC could be used to place constraints on its BH population. Some intrinsic properties of MSP systems in our models (such as the magnetic fields and spin periods) are in good overall agreement with observations, while others (such as the distribution of binary companion types) less so, and we discuss the possible reasons for such discrepancies. Interestingly, our models also demonstrate the possibility of dynamically forming NS--NS and NS--BH binaries in GCs, although the predicted numbers are very small.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

The dynamical processing of black holes in the dense cores of globular clusters (GCs), makes them efficient factories for producing binary black holes (BBHs). Here we explore the population of BBHs that form dynamically in GCs and may be observable at mHz frequencies or higher with LISA. We use our Monte Carlo stellar dynamics code, which includes gravitational radiation reaction effects for all BH encounters. By creating a representative local universe of GCs, we show that up to dozens of these systems may be resolvable by LISA with signal-to-noise ratios of at least 5. Approximately one third of these binaries will have measurable eccentricities ($e > 10^{-3}$) in the LISA band and a small number ($\lesssim 5$) may evolve from the LISA band to the LIGO band during the LISA mission.

About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary stars. However, observations suggest that the eccentricity distributions of most observed binaries, and particularly those with masses < 5 Msun, are flatter than thermal and follow more closely to a uniform distribution. Nonetheless, it is often argued that dynamical interactions in a star cluster would quickly thermalize the binaries, which could justify imposing a thermal eccentricity distribution at birth for all binaries. In this paper we investigate the validity of this assumption. We develop our own rapid semi-analytic model for binary evolution in star clusters, and also compare with detailed N-body and Monte Carlo star cluster models. We show that, for nearly all binaries, dynamical encounters fail to convert an initially uniform eccentricity distribution to thermal within a star cluster's lifetime. Thus, if a thermal eccentricity distribution is observed, it is likely imprinted upon formation rather than through subsequent long-term dynamical processing. Theoretical investigations that initialize all binaries with a thermal distribution will make incorrect predictions for the evolution of the binary population. Such models may overpredict the merger rate for binaries with modest orbital separations by a factor of about two.

Globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the so-called "core-collapsed" and "non-core-collapsed" clusters. Here, we use our H\'enon-type Monte Carlo code, CMC, to explore initial cluster parameters that map into this bimodality. Remarkably, we find that by varying the initial size of clusters (specified in our initial conditions in terms of the initial virial radius, $r_v$) within a relatively narrow range consistent with the measured radii of young star clusters in the local universe ($r_v \approx 0.5-5$ pc), our models reproduce the variety of present-day cluster properties. Furthermore, we show that stellar-mass black holes (BHs) play an intimate role in this mapping from initial conditions to the present-day structural features of GCs. We identify "best-fit" models for three GCs with known observed BH candidates, NGC 3201, M22, and M10, and show that these clusters harbor populations of $\sim 50-100$ stellar-mass BHs at present. As an alternative case, we also compare our models to the core-collapsed cluster NGC 6752 and show that this cluster likely contains few BHs at present. Additionally, we explore the formation of BH binaries in GCs and demonstrate that these systems form naturally in our models in both detached and mass-transferring configurations with a variety of companion stellar types, including low-mass main sequence stars, white dwarfs, and sub-subgiants.

Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic population of dynamically-formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission between unbound black holes during chaotic resonant encounters, with roughly half of those events having eccentricities detectable by current ground-based gravitational-wave detectors. The mergers that occur inside clusters typically have lower masses than binaries that were ejected from the cluster many Gyrs ago. Gravitational-wave captures from globular clusters contribute 1-2 Gpc^-3 yr^-1 to the binary merger rate in the local universe, increasing to ~10 Gpc^-3 yr^-1 at z~3. Finally, we discuss some of the technical difficulties associated with post-Newtonian scattering encounters, and how care must be taken when measuring the binary parameters during a dynamical capture.

The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics. Unfortunately, current techniques to model such systems suffer from severe limitations. A direct N-body simulation with more than 10^5 particles can require months or even years to complete, while an orbit-sampling Monte Carlo approach cannot adequately model the dynamics in a dense cluster core, particularly in the presence of many black holes. We have developed a new technique combining the precision of a direct N-body integration with the speed of a Monte Carlo approach. Our Rapid And Precisely Integrated Dynamics code, the RAPID code, statistically models interactions between neighboring stars and stellar binaries while integrating directly the orbits of stars or black holes in the cluster core. This allows us to accurately simulate the dynamics of the black holes in a realistic globular cluster environment without the burdensome N^2 scaling of a full N-body integration. We compare RAPID models of idealized globular clusters to identical models from the direct N-body and Monte Carlo methods. Our tests show that RAPID can reproduce the half-mass radii, core radii, black hole ejection rates, and binary properties of the direct N-body models far more accurately than a standard Monte Carlo integration while remaining significantly faster than a full N-body integration. With this technique, it will be possible to create more realistic models of Milky Way globular clusters with sufficient rapidity to explore the full parameter space of dense stellar clusters.

Planetary systems with more than two bodies will experience orbital crossings at a time related to the initial orbital separations of the planets. After a crossing, the system enters a period of chaotic evolution ending in the reshaping of the system's architecture via planetary collisions or ejections. We carry out N-body integrations on a large number of systems with equally-spaced planets (in units of the Hill radius) to determine the distribution of instability times for a given planet separation. We investigate both the time to the initiation of instability through a close encounter and the time to a planet-planet collision. We find that a significant portion of systems with non-zero mutual inclinations survive after a close encounter and do not promptly experience a planet-planet collision. Systems with significant inclinations can continue to evolve for over 1,000 times longer than the encounter time. The fraction of long lived systems is dependent on the absolute system scale and the initial inclination of the planets. These results have implications to the assumed stability of observed planetary systems.

Recent discoveries of black hole (BH) candidates in Galactic and extragalactic globular clusters (GCs) have ignited interest in understanding how BHs dynamically evolve in a GC and the number of BHs ($N_{\rm{BH}}$) that may still be retained by today's GCs. Numerical models show that even if stellar-mass BHs are retained in today's GCs, they are typically in configurations that are not directly detectable. We show that a suitably defined measure of mass segregation ($\Delta$) between, e.g., giants and low-mass main-sequence stars, can be an effective probe to indirectly estimate $N_{\rm{BH}}$ in a GC aided by calibrations from numerical models. Using numerical models including all relevant physics we first show that $N_{\rm{BH}}$ is strongly anticorrelated with $\Delta$ between giant stars and low-mass main-sequence stars. We apply the distributions of $\Delta$ vs $N_{\rm{BH}}$ obtained from models to three Milky Way GCs to predict the $N_{\rm{BH}}$ retained by them at present. We calculate $\Delta$ using the publicly available ACS survey data for 47 Tuc, M 10, and M 22, all with identified stellar-mass BH candidates. Using these measured $\Delta$ and distributions of $\Delta$ vs $N_{\rm{BH}}$ from models as calibration we predict distributions for $N_{\rm{BH}}$ expected to be retained in these GCs. For 47 Tuc, M 10, and M 22 our predicted distributions peak at $N_{\rm{BH}}\approx20$, $24$, and $50$, whereas, within the $2\sigma$ confidence level, $N_{\rm{BH}}$ can be up to $\sim150$, $50$, and $200$, respectively.

We explore the formation of double-compact-object binaries in Milky Way (MW) globular clusters (GCs) that may be detectable by the Laser Interferometer Space Antenna (LISA). We use a set of 137 fully evolved GC models that, overall, effectively match the properties of the observed GCs in the MW. We estimate that, in total, the MW GCs contain $\sim21$ sources that will be detectable by LISA. These detectable sources contain all combinations of black hole (BH), neutron star, and white dwarf components. We predict $\sim7$ of these sources will be BH-BH binaries. Furthermore, we show that some of these BH-BH binaries can have signal-to-noise ratios large enough to be detectable at the distance of the Andromeda galaxy or even the Virgo cluster.

The MeerKAT telescope represents an outstanding opportunity for radio pulsar timing science with its unique combination of a large collecting area and aperture efficiency (effective area $\sim$7500 m$^2$), system temperature ($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and ability to form up to four sub-arrays. The MeerTime project is a five-year program on the MeerKAT array by an international consortium that will regularly time over 1000 radio pulsars to perform tests of relativistic gravity, search for the gravitational wave signature induced by supermassive black hole binaries in the timing residuals of millisecond pulsars, explore the interiors of neutron stars through a pulsar glitch monitoring programme, explore the origin and evolution of binary pulsars, monitor the swarms of pulsars that inhabit globular clusters and monitor radio magnetars. In addition to these primary programmes, over 1000 pulsars will have their arrival times monitored and the data made immediately public. The MeerTime pulsar backend comprises two server-class machines each of which possess four Graphics Processing Units. Up to four pulsars can be coherently dedispersed simultaneously up to dispersion measures of over 1000 pc cm$^{-3}$. All data will be provided in psrfits format. The MeerTime backend will be capable of producing coherently dedispersed filterbank data for timing multiple pulsars in the cores of globular clusters that is useful for pulsar searches of tied array beams. All MeerTime data will ultimately be made available for public use, and any published results will include the arrival times and profiles used in the results.

We present models of realistic globular clusters with post-Newtonian dynamics for black holes. By modeling the relativistic accelerations and gravitational-wave emission in isolated binaries and during three- and four-body encounters, we find that nearly half of all binary black hole mergers occur inside the cluster, with about 10% of those mergers entering the LIGO/Virgo band with eccentricities greater than 0.1. In-cluster mergers lead to the birth of a second generation of black holes with larger masses and high spins, which, depending on the black hole natal spins, can sometimes be retained in the cluster and merge again. As a result, globular clusters can produce merging binaries with detectable spins regardless of the birth spins of black holes formed from massive stars. These second-generation black holes would also populate any upper mass gap created by pair-instability supernovae.

Numerical simulations have shown that black holes (BHs) can strongly influence the evolution and present-day observational properties of globular clusters (GCs). Using a Monte Carlo code, we construct GC models that match the Milky Way (MW) cluster NGC 3201, the first cluster in which a stellar-mass BH was identified through radial-velocity measurements. We predict that NGC 3201 contains $\gtrsim 200$ stellar-mass BHs. Furthermore, we explore the dynamical formation of main sequence-BH binaries and demonstrate that systems similar to the observed BH binary in NGC 3201 are produced naturally. Additionally, our models predict the existence of bright blue-straggler-BH binaries unique to core-collapsed clusters, which otherwise retain few BHs.

Nuclear star clusters around massive black holes are expected to be abundant in stellar mass black holes and black hole binaries. These binaries form a hierarchical triple system with the massive black hole at the center. Gravitational perturbations from the massive black hole can cause high eccentricity excitation. During this process, the eccentricity may approach unity, and the pericenter distance may become sufficiently small that gravitational wave emission drives the binary to merge. In this paper, we consider a simple proof of concept and explore the effect of the eccentric Kozai-Lidov mechanism for unequal mass binaries. We perform a set of Monte Carlo simulations on BH-BH binaries in galactic nuclei with quadrupole and octupole-level secular perturbations, general relativistic precession, and gravitational wave emission. For a nominal number of steady-state BH-BH binaries, our model gives a total merger rate $\sim 1 - 3$$Gpc^{-3} yr^{-1}$, depending on the assumed density profile. Thus, our model potentially competes with other dynamical mechanisms, such as the dynamical formations and mergers of BH binaries in globular clusters or dense nuclear clusters without a massive black hole. We provide predictions for the distributions of these LIGO sources in galactic nuclei.

We explore the population of mass-transferring binaries ejected from globular clusters (GCs) with both black hole (BH) and neutron star (NS) accretors. We use a set of 137 fully evolved globular cluster models which span a large range in cluster properties and, overall, match very well the properties of old GCs observed in the Milky Way. We identify all binaries ejected from our set of models that eventually undergo mass-transfer. These binaries are ejected from their host clusters over a wide range of ejection times and include white dwarf, giant, and main sequence donors. We calculate the orbits of these ejected systems in the Galactic potential to determine their present-day positions in the Galaxy and compare to the distribution of observed low-mass X-ray binaries (XRBs) in the Milky Way. We estimate $\sim 300$ mass-transferring NS binaries and $\sim 180$ mass-transferring BH binaries may currently be present in the Milky Way that originated from within GCs. Of these, we estimate, based on mass-transfer rates and duty cycles at the present time, at most a few would be observable as BH--XRBs and NS--XRBs at the present day. Based on our results, XRBs that originated from GCs are unlikely to contribute significantly to the total population of low-mass XRBs in the Galactic field.

Studies of high-multiplicity, tightly-packed planetary systems suggest that dynamical instabilities are common and affect both the orbits and planet structures, where the compact orbits and typically low densities make physical collisions likely outcomes. Since the structure of many of these planets is such that the mass is dominated by a rocky core, but the volume is dominated by a tenuous gas envelope, the sticky-sphere approximation, used in dynamical integrators, may be a poor model for these collisions. We perform five sets of collision calculations, including detailed hydrodynamics, sampling mass ratios and core mass fractions typical in Kepler Multis. In our primary set of calculations, we use Kepler-36 as a nominal remnant system, as the two planets have a small dynamical separation and an extreme density ratio. We use an N-body code, Mercury 6.2, to integrate initially unstable systems and study the resultant collisions in detail. We use these collisions, focusing on grazing collisions, in combination with realistic planet models created using gas profiles from Modules for Experiments in Stellar Astrophysics and core profiles using equations of state from Seager et al. (2007), to perform hydrodynamic calculations, finding scatterings, mergers, and even a potential planet-planet binary. We dynamically integrate the remnant systems, examine the stability, and estimate the final densities, finding the remnant densities are sensitive to the core masses, and collisions result in generally more stable systems. We provide prescriptions for predicting the outcomes and modeling the changes in mass and orbits following collisions for general use in dynamical integrators.

We explore the formation of mass-transferring binary systems containing black holes within globular clusters. We show that it is possible to form mass-transferring black hole binaries with main sequence, giant, and white dwarf companions with a variety of orbital parameters in globular clusters spanning a large range in present-day properties. All mass-transferring black hole binaries found in our models at late times are dynamically created. The black holes in these systems experienced a median of $\sim 30$ dynamical encounters within the cluster before and after acquiring the donor. Furthermore, we show that the presence of mass-transferring black hole systems has little correlation with the total number of black holes within the cluster at any time. This is because the net rate of formation of black hole-non-black hole binaries in a cluster is largely independent of the total number of retained black holes. Our results suggest that the detection of a mass-transferring black hole binary in a globular cluster does not necessarily indicate that the host cluster contains a large black hole population.

Gravitational waves (GWs) from binary black hole (BBH) mergers provide a new probe of massive-star evolution and the formation channels of binary compact objects. By coupling the growing sample of BBH systems with population synthesis models, we can begin to constrain the parameters of such models and glean unprecedented knowledge about the inherent physical processes that underpin binary stellar evolution. In this study, we apply a hierarchical Bayesian model to mass measurements from a synthetic GW sample to constrain the physical prescriptions in population models and the relative fraction of systems generated from various channels. We employ population models of two canonical formation scenarios in our analysis --- isolated binary evolution involving a common-envelope phase and dynamical formation within globular clusters --- with model variations for different black hole natal kick prescriptions. We show that solely with chirp mass measurements, it is possible to constrain natal kick prescriptions and the relative fraction of systems originating from each formation channel with $\mathcal{O}(100)$ of confident detections. This framework can be extended to include additional formation scenarios, model parameters, and measured properties of the compact binary.

High-multiplicity Kepler systems (referred to as Kepler multis) are often tightly packed and may be on the verge of instability. Many systems of this type could have experienced past instabilities, where the compact orbits and often low densities make physical collisions likely outcomes. We use numerical simulations to study the dynamical instabilities and planet-planet interactions in a synthetically generated sample of closely-packed, high-multiplicity systems. We focus specifically on systems resembling Kepler-11, a Kepler multi with six planets, and run a suite of dynamical integrations, sampling the initial orbital parameters around the nominal values reported in Lissauer et al. (2011a), finding that most of the realizations are unstable, resulting in orbit crossings and, eventually, collisions and mergers. We study in detail the dependence of stability on the orbital parameters of the planets and planet-pair characteristics to identify possible precursors to instability, compare the systems that emerge from dynamical instabilities to the observed Kepler sample (after applying observational corrections), and propose possible observable signatures of these instabilities. We examine the characteristics of each planet-planet collision, categorizing collisions by the degree of contact and collision energy, and find that grazing collisions are more common than direct impacts. Since the structure of many planets found in Kepler multis is such that the mass is dominated by a rocky core, but the volume is dominated by a low-density gaseous envelope, the sticky-sphere approximation may not be valid, and we present hydrodynamic calculations of planet-planet collisions clearly deviating from this approximation. Finally, we rerun a subset of our dynamical calculations using instead a modified prescription to handle collisions, finding, in general, higher multiplicity remnant systems.

Disk migration and high-eccentricity migration are two well-studied theories to explain the formation of hot Jupiters. The former predicts that these planets can migrate up until the planet-star Roche separation ($a_{Roche}$) and the latter predicts they will tidally circularize at a minimum distance of 2$a_{Roche}$. Considering long-running radial velocity and transit surveys have identified a couple hundred hot Jupiters to date, we can revisit the classic question of hot Jupiter formation in a data-driven manner. We approach this problem using data from several exoplanet surveys (radial velocity, Kepler, HAT, and WASP) allowing for either a single population or a mixture of populations associated with these formation channels, and applying a hierarchical Bayesian mixture model of truncated power laws of the form $x^{\gamma-1}$ to constrain the population-level parameters of interest (e.g., location of inner edges, $\gamma$, mixture fractions). Within the limitations of our chosen models, we find the current radial velocity and Kepler sample of hot Jupiters can be well explained with a single truncated power law distribution with a lower cutoff near 2$a_{Roche}$, a result that still holds after a decade, and $\gamma=-0.51\pm^{0.19}_{0.20}$. However, the HAT and WASP data show evidence for multiple populations (Bayes factor $\approx 10^{21}$). We find that $15\pm^{9}_{6}\%$ reside in a component consistent with disk migration ($\gamma=-0.04\pm^{0.53}_{1.27}$) and $85\pm^{6}_{9}\%$ in one consistent with high-eccentricity migration ($\gamma=-1.38\pm^{0.32}_{0.47}$). We find no immediately strong connections with some observed host star properties and speculate on how future exoplanet surveys could improve upon hot Jupiter population inference.

Using numerical models for star clusters spanning a wide range in ages and metallicities (Z) we study the masses of binary black holes (BBHs) produced dynamically and merging in the local universe ($z\lesssim0.2$). After taking into account cosmological constraints on star-formation rate and metallicity evolution, which realistically relate merger delay times obtained from models with merger redshifts, we show here for the first time that while old, metal-poor globular clusters can naturally produce merging BBHs with heavier components, as observed in GW150914, lower-mass BBHs like GW151226 are easily formed dynamically in younger, higher-metallicity clusters.More specifically, we show that the mass of GW151226 is well within $1\sigma$ of the mass distribution obtained from our models for clusters with Z/Zsun$\gtrsim0.5$. Indeed dynamical formation of a system like GW151226 likely requires a cluster that is younger and has a higher metallicity than typical Galactic globular clusters. The LVT151012 system, if real, could have been created in any cluster with Z/Zsun$\lesssim0.25$. On the other hand, GW150914 is more massive (beyond $1\sigma$) than typical BBHs from even the lowest-metallicity (Z/Zsun=0.005) clusters we consider, but is within $2\sigma$ of the intrinsic mass distribution from our cluster models with Z/Zsun$\lesssim0.05$; of course detection biases also push the observed distributions towards higher masses.

Recent N-body simulations predict that large numbers of stellar black holes (BHs) could remain bound to globular clusters (GCs) at present, and merging BH--BH binaries are produced dynamically in significant numbers. We systematically vary "standard" assumptions made by numerical simulations related to, e.g., BH formation, stellar winds, binary properties of high-mass stars, and IMF within existing uncertainties, and study the effects on the evolution of the structural properties of GCs, and the BHs in GCs. We find that variations in initial assumptions can set otherwise identical initial clusters on completely different evolutionary paths, significantly affecting their present observable properties, or even affecting the cluster's very survival to the present. However, these changes usually do not affect the numbers or properties of local BH--BH mergers. The only exception is that variations in the assumed winds and IMF can change the masses and numbers of local BH--BH mergers, respectively. All other variations (e.g., in initial binary properties and binary fraction) leave the masses and numbers of locally merging BH--BH binaries largely unchanged. This is in contrast to binary population synthesis models for the field, where results are very sensitive to many uncertain parameters in the initial binary properties and binary stellar-evolution physics. We find that weak winds are required for producing GW150914-like mergers from GCs at low redshifts. LVT151012 can be produced in GCs modeled both with strong and weak winds. GW151226 is lower-mass than typical mergers from GCs modeled with weak winds, but is similar to mergers from GCs modeled with strong winds.

The recent detections of the binary black hole mergers GW150914 and GW151226 have inaugurated the field of gravitational-wave astronomy. For the two main formation channels that have been proposed for these sources, isolated binary evolution in galactic fields and dynamical formation in dense star clusters, the predicted masses and merger rates overlap significantly, complicating any astrophysical claims that rely on measured masses alone. Here, we examine the distribution of spin- orbit misalignments expected for binaries from the field and from dense star clusters. Under standard assumptions for black-hole natal kicks, we find that black-hole binaries similar to GW150914 could be formed with significant spin-orbit misalignment only through dynamical processes. In particular, these heavy-black-hole binaries can only form with a significant spin-orbit anti-alignment in the dynamical channel. Our results suggest that future detections of merging black hole binaries with measurable spins will allow us to identify the main formation channel for these systems.

The recent detections of GW150914 and GW151226 imply an abundance of stellar-mass binary-black-hole mergers in the local universe. While ground-based gravitational-wave detectors are limited to observing the final moments before a binary merges, space-based detectors, such as the Laser Interferometer Space Antenna (LISA), can observe binaries at lower orbital frequencies where such systems may still encode information about their formation histories. In particular, the orbital eccentricity and mass of binary black holes in the LISA frequency band can be used together to discriminate between binaries formed in isolation in galactic fields and those formed in dense stellar environments such as globular clusters. In this letter, we explore the orbital eccentricity and mass of binary-black-hole populations as they evolve through the LISA frequency band. Overall we find that there are two distinct populations discernible by LISA. We show that up to ~90% of binaries formed either dynamically or in isolation have eccentricities measurable by LISA. Finally, we note how measured eccentricities of low-mass binary black holes evolved in isolation could provide detailed constraints on the physics of black-hole natal kicks and common-envelope evolution.

Motivated by the recent detection of gravitational waves from the black hole binary merger GW150914, we study the dynamical evolution of black holes in galactic nuclei where massive star clusters reside. With masses of ~10^7M_Sun and sizes of only a few parsecs, nuclear star clusters are the densest stellar systems observed in the local universe and represent a robust environment where (stellar mass) black hole binaries can dynamically form, harden and merge. We show that due to their large escape speeds, nuclear star clusters can keep a large fraction of their merger remnants while also evolving rapidly enough that the holes can sink back to the central regions where they can swap in new binaries that can subsequently harden and merge. This process can repeat several times and produce black hole mergers of several tens of solar masses similar to GW150914 and up to a few hundreds of solar masses, without the need of invoking extremely low metallicity environments or implausible initial conditions. We use a semi-analytical approach to describe the formation and dynamics of black holes in massive star clusters. We find a black hole binary merger rate per volume from nuclear star clusters of ~1.5 Gpc^-3 yr^-1, implying up to a few tens of possible detections per year with Advanced LIGO. Our models suggest a local merger rate of 0.3- 1 Gpc^-3 yr^-1 for high mass black hole binaries similar to GW150914 (total mass >~ 50 M_Sun, redshift z< 0.3); a merger rate comparable to that of high mass black hole binaries that are dynamically assembled in globular clusters. Finally, we show that if all black holes receive high natal kicks, >~50km s^-1, then nuclear star clusters could dominate the local merger rate of binary black holes compared to the merger rate of similar binaries produced in either globular clusters or through isolated binary evolution.

We explore the possibility that GW150914, the binary black hole merger recently detected by Advanced LIGO, was formed by gravitational interactions in the core of a dense star cluster. Using models of globular clusters with detailed $N$-body dynamics and stellar evolution, we show that a typical cluster with a mass of $3\times10^5M_{\odot}$ to $6\times10^5M_{\odot}$ is optimal for forming GW150914-like binary black holes that will merge in the local universe. We identify the most likely dynamical processes for forming GW150914 in such a cluster, and we show that the detection of GW150914 is consistent with the masses and merger rates expected for binary black hole mergers from globular clusters. Our results show that dynamical processes provide a significant and well-understood pathway for forming binary black hole mergers in the local universe. Understanding the contribution of dynamics to the binary black hole merger problem is a critical step in unlocking the full potential of gravitational-wave astronomy.

The recent discovery of GW150914, the binary black hole merger detected by Advanced LIGO, has the potential to revolutionize observational astrophysics. But to fully utilize this new window into the universe, we must compare these new observations to detailed models of binary black hole formation throughout cosmic time. Expanding upon our previous work (Rodriguez et al., 2015), we study merging binary black holes formed in globular clusters using our Monte Carlo approach to stellar dynamics. We have created a new set of 52 cluster models with different masses, metallicities, and radii to fully characterize the binary black hole merger rate. These models include all the relevant dynamical processes (such as two-body relaxation, strong encounters, and three-body binary formation) and agree well with detailed direct N-body simulations. In addition, we have enhanced our stellar evolution algorithms with updated metallicity-dependent stellar wind and supernova prescriptions, allowing us to compare our results directly to the most recent population synthesis predictions for merger rates from isolated binary evolution. We explore the relationship between a cluster's global properties and the population of binary black holes that it produces. In particular, we derive a numerically calibrated relationship between the merger times of ejected black hole binaries and a cluster's mass and radius. We explore the masses and mass ratios of these binaries as a function of redshift, and find a merger rate of ~5 Gpc$^{-3}$ yr$^{-1}$ in the local universe, with 80% of sources having total masses from $32M_{\odot}$ to $64M_{\odot}$. Under standard assumptions, approximately 1 out of every 7 binary black hole mergers in the local universe will have originated in a globular cluster, but we also explore the sensitivity of this result to different assumptions for binary stellar evolution. [Abridged]

We present the first detailed comparison between million-body globular cluster simulations computed with a H\'enon-type Monte Carlo code, CMC, and a direct $N$-body code, NBODY6++GPU. Both simulations start from an identical cluster model with $10^6$ particles, and include all of the relevant physics needed to treat the system in a highly realistic way. With the two codes "frozen" (no fine-tuning of any free parameters or internal algorithms of the codes) we find excellent agreement in the overall evolution of the two models. Furthermore, we find that in both models, large numbers of stellar-mass black holes (> 1000) are retained for 12 Gyr. Thus, the very accurate direct $N$-body approach confirms recent predictions that black holes can be retained in present-day, old globular clusters. We find only minor disagreements between the two models and attribute these to the small-$N$ dynamics driving the evolution of the cluster core for which the Monte Carlo assumptions are less ideal. Based on the overwhelming general agreement between the two models computed using these vastly different techniques, we conclude that our Monte Carlo approach, which is more approximate, but dramatically faster compared to the direct $N$-body, is capable of producing a very accurate description of the long-term evolution of massive globular clusters even when the clusters contain large populations of stellar-mass black holes.

The predicted rate of binary black hole mergers from galactic fields can vary over several orders of magnitude and is extremely sensitive to the assumptions of stellar evolution. But in dense stellar environments such as globular clusters, binary black holes form by well-understood gravitational interactions. In this letter, we study the formation of black hole binaries in an extensive collection of realistic globular cluster models. By comparing these models to observed Milky Way and extragalactic globular clusters, we find that the mergers of dynamically-formed binaries could be detected at a rate of ~100 per year, potentially dominating the binary black hole merger rate. We also find that a majority of cluster-formed binaries are more massive than their field-formed counterparts, suggesting that Advanced LIGO could identify certain binaries as originating from dense stellar environments.

Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions between the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the inner orbit ("Lidov-Kozai cycles"), which can lead to a collision and merger of the two inner components. In this paper we use Monte Carlo models of dense star clusters to identify all triple systems formed dynamically and we compute their evolution using a highly accurate three-body integrator which incorporates relativistic and tidal effects. We find that a large fraction of these triples evolve through a non-secular dynamical phase which can drive the inner binary to higher eccentricities than predicted by the standard secular perturbation theory (even including octupole-order terms). We place constraints on the importance of Lidov-Kozai-induced mergers for producing: (i) gravitational wave sources detectable by Advanced LIGO (aLIGO), for triples with an inner pair of stellar black holes; and (ii) blue straggler stars, for triples with main-sequence-star components. We find a realistic aLIGO detection rate of black hole mergers due to the Lidov-Kozai mechanism of 1yr^-1, with about 20% of these having a finite eccentricity when they first chirp into the aLIGO frequency band. While rare, these events are likely to dominate among eccentric compact object inspirals that are potentially detectable by aLIGO. For blue stragglers, we find that the Lidov-Kozai mechanism can contribute only up to ~10% of their total numbers in globular clusters.

The formation of Black Hole (BH) Low-Mass X-ray Binaries (LMXB) poses a theoretical challenge, as low-mass companions are not expected to survive the common-envelope scenario with the BH progenitor. Here we propose a formation mechanism that skips the common-envelope scenario and relies on triple-body dynamics. We study the evolution of hierarchical triples, following the secular dynamical evolution up to the octupole-level of approximation, including general relativity, tidal effects and post-main-sequence evolution, such as mass loss, changes to stellar radii and supernovae. During the dynamical evolution of the triple system, the "eccentric Kozai-Lidov" mechanism can cause large eccentricity excitations in the LMXB progenitor, resulting in three main BH-LMXB formation channels. Here we define BH-LMXB candidates as systems where the inner BH companion star crosses its Roche limit. In the "eccentric" channel (~ 81% of the LMXBs in our simulations), the donor star crosses its Roche limit during an extreme eccentricity excitation, while still on a wide orbit. Second, we find a "giant" LMXB channel (~ 11%), where a system undergoes only moderate eccentricity excitations, but the donor star fills its Roche lobe after evolving toward the giant branch. Third, we identify a "classical" channel (~8%), where tidal forces and magnetic braking shrink and circularize the orbit to short periods, triggering mass transfer. Finally, for the giant channel, we predict an eccentric ($\sim 0.3-0.6$), preferably inclined (~40, ~140 degreed) tertiary, typically on a wide enough orbit (~10^4AU), to potentially become unbound later in the triple evolution.

Many exoplanets have now been detected in orbits with ultra-short periods, very close to the Roche limit. Building upon our previous work, we study the possibility that mass loss through Roche lobe overflow (RLO) may affect the evolution of these planets, and could possibly transform a hot Jupiter into a lower-mass planet (hot Neptune or super-Earth). We focus here on systems in which the mass loss occurs slowly ("stable mass transfer" in the language of binary star evolution) and we compute their evolution in detail with the binary evolution code MESA. We include the effects of tides, RLO, irradiation and photo-evaporation of the planet, as well as the stellar wind and magnetic braking. Our calculations all start with a hot Jupiter close to its Roche limit, in orbit around a sun-like star. The initial orbital decay and onset of RLO are driven by tidal dissipation in the star. We confirm that such a system can indeed evolve to produce lower-mass planets in orbits of a few days. The RLO phase eventually ends and, depending on the details of the mass transfer and on the planetary core mass, the orbital period can remain around a few days for several Gyr. The remnant planets have a rocky core and some amount of envelope material, which is slowly removed via photo-evaporation at nearly constant orbital period; these have properties resembling many of the observed super-Earths and sub-Neptunes. For these remnant planets we also predict an anti-correlation between mass and orbital period; very low-mass planets ($M_{\rm pl}\,\lesssim\,5\,M_{\oplus}$) in ultra-short periods ($P_{\rm orb}$<1d) cannot be produced through this type of evolution.

Massive stars are usually found in binaries, and binaries with periods less than 10 days may have a preference for near equal component masses. In this paper we investigate the evolution of these binaries all the way to contact and the possibility that these systems can be progenitors of double neutron star binaries. The small orbital separations of observed double neutron star binaries suggest that the progenitor systems underwent a common envelope phase at least once during their evolution. Bethe & Brown (1998) proposed that massive binary twins will undergo a common envelope evolution while both components are ascending the red giant branch or asymptotic giant branch simultaneously, also known as double-core evolution. Using models generated from the stellar evolution code Evolve Zero Age Main Sequence, we determine the range of mass ratios resulting in both components simultaneously ascending the RGB or AGB as a function of the difference in birth times, t. We find that, even for a generous t=5 Myr, the minimum mass ratio qmin=0.933 for an 8 Solar Mass primary and increases for larger primaries. We use a hydrodynamics code, StarSmasher, to study specifically the evolution of q=1 common envelope systems as a function of initial component mass, age, and orbital separation. We find the dynamical stability limit, the largest orbital separation where the binary becomes dynamically unstable, as a function of the component mass and age. Finally, we calculate the efficiency of ejecting matter during the inspiral phase to extrapolate the properties of the remnant binary from our numerical results, assuming the common envelope is completely ejected. We find that for the nominal core masses, there is a minimum orbital separation for a given component mass such that the helium cores survive common envelope evolution in a tightly bound binary and are viable progenitors for double neutron stars.

Our current understanding of the stellar initial mass function and massive star evolution suggests that young globular clusters may have formed hundreds to thousands of stellar-mass black holes, the remnants of stars with initial masses from $\sim 20 - 100\, M_\odot$. Birth kicks from supernova explosions may eject some black holes from their birth clusters, but most should be retained. Using a Monte Carlo method we investigate the long-term dynamical evolution of globular clusters containing large numbers of stellar black holes. We describe numerical results for 42 models, covering a range of realistic initial conditions, including up to $1.6\times10^6$ stars. In almost all models we find that significant numbers of black holes (up to $\sim10^3$) are retained all the way to the present. This is in contrast to previous theoretical expectations that most black holes should be ejected dynamically within a few Gyr. The main reason for this difference is that core collapse driven by black holes (through the Spitzer "mass segregation instability") is easily reverted through three-body processes, and involves only a small number of the most massive black holes, while lower-mass black holes remain well-mixed with ordinary stars far from the central cusp. Thus the rapid segregation of stellar black holes does not lead to a long-term physical separation of most black holes into a dynamically decoupled inner core, as often assumed previously. Combined with the recent detections of several black hole X-ray binary candidates in Galactic globular clusters, our results suggest that stellar black holes could still be present in large numbers in many globular clusters today, and that they may play a significant role in shaping the long-term dynamical evolution and the present-day dynamical structure of many clusters.

Kepler-56 is a multi-planet system containing two coplanar inner planets that are in orbits misaligned with respect to the spin axis of the host star, and an outer planet. Various mechanisms have been proposed to explain the broad distribution of spin-orbit angles among exoplanets, and these theories fall under two broad categories. The first is based on dynamical interactions in a multi-body system, while the other assumes that disk migration is the driving mechanism in planetary configuration and that the star (or disk) is titled with respect to the planetary plane. Here we show that the large observed obliquity of Kepler-56 system is consistent with a dynamical origin. In addition, we use observations by Huber et al. (2013) to derive the obliquity's probability distribution function, thus improving the constrained lower limit. The outer planet may be the cause of the inner planets' large obliquities, and we give the probability distribution function of its inclination, which depends on the initial orbital configuration of the planetary system. We show that even in the presence of precise measurement of the true obliquity, one cannot distinguish the initial configurations. Finally we consider the fate of the system as the star continues to evolve beyond the main sequence, and we find that the obliquity of the system will not undergo major variations as the star climbs the red giant branch. We follow the evolution of the system and find that the innermost planet will be engulfed in ~129 Myr. Furthermore we put an upper limit of ~155 Myr for the engulfment of the second planet. This corresponds to ~ 3% of the current age of the star.

Through tidal dissipation in a slowly spinning host star the orbits of many hot Jupiters may decay down to the Roche limit. We expect that in most cases the ensuing mass transfer will be stable. Using detailed numerical calculations we find that this evolution is quite rapid, potentially leading to complete removal of the gaseous envelope in a few Gyr, and leaving behind an exposed rocky core ("hot super-Earth"). Final orbital periods are quite sensitive to the details of the planet's mass-radius relation, and to the effects of irradiation and photo-evaporation, but could be as short as a few hours, or as long as several days. Our scenario predicts the existence of planets with intermediate masses ("hot Neptunes") that should be found precisely at their Roche limit and in the process of losing mass through Roche lobe overflow. The observed excess of small single-planet candidate systems observed by Kepler may also be the result of this process. If so, the properties of their host stars should track those of the hot Jupiters. Moreover, the number of systems that produced hot Jupiters could be 2-3 times larger than one would infer from contemporary observations.

Hot Jupiters formed through circularization of high-eccentricity orbits should be found at orbital separations $a$ exceeding $twice$ that of their Roche limit $a_{\rm R}$. Nevertheless, about a dozen giant planets have now been found well within this limit ($a_{\rm R}<a<2a_{\rm R}$), with one coming as close as 1.2$a_{\rm R}$. In this Letter, we show that orbital decay (starting beyond 2$a_{\rm R}$) driven by tidal dissipation in the star can naturally explain these objects. For a few systems (WASP-4 and 19), this explanation requires the linear reduction in convective tidal dissipation proposed originally by Zahn (1966, 1989) and verified by recent numerical simulations Penev et al. (2007), but rules out the quadratic prescription proposed by Goldreich & Nicholson (1977). Additionally, we find that WASP-19-like systems could potentially provide direct empirical constraints on tidal dissipation, as we could soon be able to measure their orbital decay through high precision transit timing measurements.

Two formation scenarios have been proposed to explain the tight orbits of hot Jupiters. They could be formed in orbits with a small inclination (with respect to the stellar spin) via disk migration, or in more highly inclined orbits via high-eccentricity migration, where gravitational interactions with a companion and tidal dissipation are at play. Here we target hot Jupiter systems where the misalignment $\lambda$ has been inferred observationally and we investigate whether their properties are consistent with high-eccentricity migration. Specifically, we study whether stellar tides can be responsible for the observed distribution of $\lambda$ and orbital separations. Improving on previous studies, we use detailed models for each star, thus accounting for how convection (and tidal dissipation) depends on stellar properties. In line with observations suggesting that hotter stars have higher $\lambda$, we find that $\lambda$ increases as the amount of stellar surface convection decreases. This trend supports the hypothesis that tides are the mechanism shaping the observed distribution of $\lambda$. Furthermore, we study the past orbital evolution of five representative systems, chosen to cover a variety of temperatures and misalignments. We consider various initial orbital configurations and integrate the equations describing the coupled evolution of the orbital separation, stellar spin, and misalignment. We account for stellar tides and wind mass loss, stellar evolution, and magnetic braking. We show that the current properties of these five representative systems can be explained naturally, given our current understanding of tidal dissipation and with physically motivated assumptions for the effects driving the orbital evolution.

Interacting galaxies often have complexes of hundreds of young stellar clusters of individual masses $\sim 10^{4-6}~M_\odot$ in regions that are a few hundred parsecs across. These cluster complexes interact dynamically, and their coalescence is a candidate for the origin of some ultracompact dwarf galaxies (UCDs). Individual clusters with short relaxation times are candidates for the production of intermediate-mass black holes of a few hundred solar masses, via runaway stellar collisions prior to the first supernovae in a cluster. It is therefore possible that a cluster complex hosts multiple intermediate-mass black holes that may be ejected from their individual clusters due to mergers or binary processes, but bound to the complex as a whole. Here we explore the dynamical interaction between initially free-flying massive black holes and clusters in an evolving cluster complex. We find that, after hitting some clusters, it is plausible that the massive black hole will be captured in an ultracompact dwarf forming near the center of the complex. In the process, the hole typically triggers electromagnetic flares via stellar disruptions, and is also likely to be a prominent source of gravitational radiation for the advanced ground-based detectors LIGO and VIRGO. We also discuss other implications of this scenario, notably that the central black hole could be considerably larger than expected in other formation scenarios for ultracompact dwarfs.

Observations of exoplanets over the last two decades have revealed a new class of Jupiter-size planets with orbital periods of a few days, the so-called "hot Jupiters". Recent measurements using the Rossiter-McLaughlin effect have shown that many (~ 50%) of these planets are misaligned; furthermore, some (~ 15%) are even retrograde with respect to the stellar spin axis. Motivated by these observations, we explore the possibility of forming retrograde orbits in hierarchical triple configurations consisting of a star-planet inner pair with another giant planet, or brown dwarf, in a much wider orbit. Recently Naoz et al. (2011) showed that in such a system, the inner planet's orbit can flip back and forth from prograde to retrograde, and can also reach extremely high eccentricities. Here we map a significant part of the parameter space of dynamical outcomes for these systems. We derive strong constraints on the orbital configurations for the outer perturber that could lead to the formation of hot Jupiters with misaligned or retrograde orbits. We focus only on the secular evolution, neglecting other dynamical effects such as mean-motion resonances, as well as all dissipative forces. For example, with an inner Jupiter-like planet initially on a nearly circular orbit at 5 AU, we show that a misaligned hot Jupiter is likely to be formed in the presence of a more massive planetary companion (> 2 MJ) within 140 AU of the inner system, with mutual inclination 50 degrees and eccentricity above 0.25. This is in striking contrast to the test-particle approximation, where an almost perpendicular configuration can still cause large eccentricity excitations, but flips of an inner Jupiter-like planet are much less likely to occur. The constraints we derive can be used to guide future observations, and, in particular, searches for more distant companions in systems containing a hot Jupiter.

We present theoretical models for stellar black hole (BH) properties in young, massive star clusters. Using a Monte Carlo code for stellar dynamics, we model realistic star clusters with $N\simeq 5\times10^5$ stars and significant binary fractions (up to 50%) with self-consistent treatments of stellar dynamics and stellar evolution. We compute the formation rates and characteristic properties of single and binary BHs for various representative ages, cluster parameters, and metallicities. Because of dynamical interactions and supernova (SN) kicks, more single BHs end up retained in clusters compared to BHs in binaries. We also find that the ejection of BHs from a cluster is a strong function of initial density. In low-density clusters (where dynamical effects are negligible), it is mainly SN kicks that eject BHs from the cluster, whereas in high-density clusters (initial central density $\rho_c(0) \sim 10^5 \, M_\odot\, {\rm pc}^{-3} $ in our models) the BH ejection rate is enhanced significantly by dynamics. Dynamical interactions of binary systems in dense clusters also modify the orbital period and eccentricity distributions while also increasing the probability of a BH having a more massive companion.

Blue straggler stars (BSS) are abundantly observed in all Galactic globular clusters (GGC) where data exist. However, observations alone cannot reveal the relative importance of various formation channels or the typical formation times for this well studied population of anomalous stars. Using a state-of-the-art H\'enon-type Monte Carlo code that includes all relevant physical processes, we create 128 models with properties typical of the observed GGCs. These models include realistic numbers of single and binary stars, use observationally motivated initial conditions, and span large ranges in central density, concentration, binary fraction, and mass. Their properties can be directly compared with those of observed GGCs. We can easily identify the BSSs in our models and determine their formation channels and birth times. We find that for central densities above ~10^3 Msun/pc^3 the dominant formation channel is stellar collisions while for lower density clusters, mass transfer in binaries provides a significant contribution (up to ~ 60% in our models). The majority of these collisions are binary-mediated, occurring during 3-body and 4-body interactions. As a result a strong correlation between the specific frequency of BSSs and the binary fraction in a cluster can be seen in our models. We find that the number of BSSs in the core shows only a weak correlation with the collision rate estimator \Gamma traditionally used by observers, in agreement with the latest Hubble Space Telescope (ACS) data. Using an idealized "full mixing" prescription for collision products, our models indicate that the BSSs observed today may have formed several Gyrs ago. However, denser clusters tend to have younger (~1 Gyr) BSSs.

We created artificial color-magnitude diagrams of Monte Carlo dynamical models of globular clusters, and then used observational methods to determine the number of blue stragglers in those clusters. We compared these blue stragglers to various cluster properties, mimicking work that has been done for blue stragglers in Milky Way globular clusters to determine the dominant formation mechanism(s) of this unusual stellar population. We find that a mass-based prescription for selecting blue stragglers will choose approximately twice as many blue stragglers than a selection criterion that was developed for observations of real clusters. However, the two numbers of blue stragglers are well-correlated, so either selection criterion can be used to characterize the blue straggler population of a cluster. We confirm previous results that the simplified prescription for the evolution of a collision or merger product in the BSE code overestimates their lifetimes. We show that our model blue stragglers follow similar trends with cluster properties (core mass, binary fraction, total mass, collision rate) as the true Milky Way blue stragglers, as long as we restrict ourselves to model clusters with an initial binary fraction higher than 5%. We also show that, in contrast to earlier work, the number of blue stragglers in the cluster core does have a weak dependence on the collisional parameter Gamma in both our models and in Milky Way globular clusters.

Here we present CAFein, a new computational tool for investigating radiative dissipation of dynamic tides in close binaries and of non-adiabatic, non-radial stellar oscillations in isolated stars in the linear regime. For the latter, CAFein computes the non-adiabatic eigenfrequencies and eigenfunctions of detailed stellar models. The code is based on the so-called Riccati method, a numerical algorithm that has been successfully applied to a variety of stellar pulsators, and which doesn't suffer of the major drawbacks of commonly used shooting and relaxation schemes. Here we present an extension of the Riccati method to investigate dynamic tides in close binaries. We demonstrate CAFein's capabilities as a stellar pulsation code both in the adiabatic and non-adiabatic regime, by reproducing previously published eigenfrequencies of a polytrope, and by successfully identifying the unstable modes of a stellar model in the $\beta$ Cephei/SPB region of the Hertzsprung-Russell diagram. Finally, we verify CAFein's behavior in the dynamic tides regime by investigating the effects of dynamic tides on the eigenfunctions and orbital and spin evolution of massive Main Sequence stars in eccentric binaries, and of hot Jupiter host stars. The plethora of asteroseismic data provided by the NASA's Kepler satellite, some of which include the direct detection of tidally excited stellar oscillations, make CAFein quite timely. Furthermore, the increasing number of observed short-period detached double white dwarfs (WD) and the observed orbital decay in the tightest of such binaries open up a new possibility of investigating WD interiors through the effects of tides on their orbital evolution

The secular approximation for the evolution of hierarchical triple configurations has proven to be very useful in many astrophysical contexts, from planetary to triple-star systems. In this approximation the orbits may change shape and orientation, on time scales longer than the orbital time scales, but the semi major axes are constant. For example, for highly inclined triple systems, the Kozai-Lidov mechanism can produce large-amplitude oscillations of the eccentricities and inclinations. Here we revisit the secular dynamics of hierarchical triple systems. We derive the secular evolution equations to octupole order in Hamiltonian perturbation theory. Our derivation corrects an error in some previous treatments of the problem that implicitly assumed a conservation of the z-component of the angular momentum of the inner orbit (i.e., parallel to the total angular momentum of the system). Already to quadrupole order, our results show new behaviors including the possibility for a system to oscillate from prograde to retrograde orbits. At the octupole order, for an eccentric outer orbit, the inner orbit can reach extremely high eccentricities and undergo chaotic flips in its orientation. We discuss applications to a variety of astrophysical systems, from stellar triples to merging compact binaries and planetary systems. Our results agree with those of previous studies done to quadrupole order only in the limit in which one of the inner two bodies is a massless test particle and the outer orbit is circular;our results agree with previous studies at octupole order for the eccentricity evolution, but not for the inclination evolution.

Globular clusters should be born with significant numbers of stellar-mass black holes (BHs). It has been thought for two decades that very few of these BHs could be retained through the cluster lifetime. With masses ~10 MSun, BHs are ~20 times more massive than an average cluster star. They segregate into the cluster core, where they may eventually decouple from the remainder of the cluster. The small-N core then evaporates on a short timescale. This is the so-called Spitzer instability. Here we present the results of a full dynamical simulation of a globular cluster containing many stellar-mass BHs with a realistic mass spectrum. Our Monte Carlo simulation code includes detailed treatments of all relevant stellar evolution and dynamical processes. Our main finding is that old globular clusters could still contain many BHs at present. In our simulation, we find no evidence for the Spitzer instability. Instead, most of the BHs remain well-mixed with the rest of the cluster, with only the innermost few tens of BHs segregating significantly. Over the 12 Gyr evolution, fewer than half of the BHs are dynamically ejected through strong binary interactions in the cluster core. The presence of BHs leads to long-term heating of the cluster, ultimately producing a core radius on the high end of the distribution for Milky Way globular clusters (and those of other galaxies). A crude extrapolation from our model suggests that the BH--BH merger rate from globular clusters could be comparable to the rate in the field.

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 a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parallel random number generation scheme, as well as a parallel sorting algorithm, required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. The implementation uses the Message Passing Interface (MPI) library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude, from 10^5 to 10^7. We find that our results are in good agreement with self-similar core-collapse solutions, and the core collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within less than 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N=10^5, 128 for N=10^6 and 256 for N=10^7. The runtime reaches a saturation with the addition of more processors beyond these limits which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60x, 100x, and 220x, respectively.

While hundreds of planets have been discovered around field stars, only a few are known in star clusters. To explain the lack of short-period giant planets in globular clusters (GC), such as 47 Tucane and \omega\ Centauri, it has been suggested that their low metallicities may have prevented planet formation. Alternatively, the high rates of close stellar encounters in these clusters may have influenced the formation and subsequent evolution of planetary systems. How common are planets in clusters around normal main-sequence stars? Here we consider whether this question can be addressed using data from the Kepler mission. The Kepler field of view contains 4 low-density (relative to GCs) open clusters where the metallicities are about solar (or even higher) and stellar encounters are much less frequent than in typical GCs. We provide detailed $N$-body models and show that most planets in Kepler-detectable orbits are not significantly perturbed by stellar encounters in these open clusters. We focus on the most massive cluster, NGC 6791, which has super-solar metallicity, and find that if planets formed in this cluster at the same frequency as observed in the field, Kepler could detect 1 -- 20 transiting planets depending on the planet-size distribution and the duration of data collection. Due to the large distance to NGC 6791 Kepler will have to search relatively faint ($K_p<20$) stars for the full extended mission to achieve such a yield.

Decades after the first predictions of intermediate-mass black holes (IMBHs) in globular clusters (GCs) there is still no unambiguous observational evidence for their existence. The most promising signatures for IMBHs are found in the cores of GCs, where the evidence now comes from the stellar velocity distribution, the surface density profile, and, for very deep observations, the mass-segregation profile near the cluster center. However, interpretation of the data, and, in particular, constraints on central IMBH masses, require the use of detailed cluster dynamical models. Here we present results from Monte Carlo cluster simulations of GCs that harbor IMBHs. As an example of application, we compare velocity dispersion, surface brightness and mass-segregation profiles with observations of the GC M10, and constrain the mass of a possible central IMBH in this cluster. We find that, although M10 does not seem to possess a cuspy surface density profile, the presence of an IMBH with a mass up to 0.75% of the total cluster mass, corresponding to about 600 Msun, cannot be excluded. This is also in agreement with the surface brightness profile, although we find it to be less constraining, as it is dominated by the light of giants, causing it to fluctuate significantly. We also find that the mass-segregation profile cannot be used to discriminate between models with and without IMBH. The reason is that M10 is not yet dynamically evolved enough for the quenching of mass segregation to take effect. Finally, detecting a velocity dispersion cusp in clusters with central densities as low as in M10 is extremely challenging, and has to rely on only 20-40 bright stars. It is only when stars with masses down to 0.3 Msun are included that the velocity cusp is sampled close enough to the IMBH for a significant increase above the core velocity dispersion to become detectable.

We study the production of Hot Jupiters (HJs) in stellar binaries. We show that the "eccentric Kozai-Lidov" (EKL) mechanism can play a key role in the dynamical evolution of a star-planet-star triple system. We run a large set of Monte Carlo simulations including the secular evolution of the orbits, general relativistic precession, and tides, and we determine the semi-major axis, eccentricity, inclination and spin-orbit angle distributions of the HJs that are produced. We explore the effect of different tidal friction parameters on the results. We find that the efficiency of forming HJs when taking the EKL mechanism into account is higher then previously estimated. Accounting for the frequency of stellar binaries, we find that this production mechanism can account for about 30% of the observed HJ population. Current observations of spin-orbit angles are consistent with this mechanism producing \sim 30% of all HJs, and up to 100% of the misaligned systems. Based on the properties of binaries without a HJ in our simulations, we predict the existence of many Jupiter-like planets with moderately eccentric and inclined orbits and semi-major axes of several AU.

A promising mechanism to form intermediate-mass black holes (IMBHs) is the runaway merger in dense star clusters, where main-sequence stars collide and form a very massive star (VMS), which then collapses to a black hole. In this paper we study the effects of primordial mass segregation and the importance of the stellar initial mass function (IMF) on the runaway growth of VMSs using a dynamical Monte Carlo code for N-body systems with N as high as 10^6 stars. Our code now includes an explicit treatment of all stellar collisions. We place special emphasis on the possibility of top-heavy IMFs, as observed in some very young massive clusters. We find that both primordial mass segregation and the shape of the IMF affect the rate of core collapse of star clusters and thus the time of the runaway. When we include primordial mass segregation we generally see a decrease in core collapse time (tcc). Moreover, primordial mass segregation increases the average mass in the core, thus reducing the central relaxation time, which also decreases tcc. The final mass of the VMS formed is always close to \sim 10^-3 of the total cluster mass, in agreement with the previous studies and is reminiscent of the observed correlation between the central black hole mass and the bulge mass of the galaxies. As the degree of primordial mass segregation is increased, the mass of the VMS increases at most by a factor of 3. Flatter IMFs generally increase the average mass in the whole cluster, which increases tcc. For the range of IMFs investigated in this paper, this increase in tcc is to some degree balanced by stellar collisions, which accelerate core collapse. Thus there is no significant change in tcc for the somewhat flatter global IMFs observed in very young massive clusters.

We review the current status of studies of the coalescence of binary neutron star systems. We begin with a discussion of the formation channels of merging binaries and we discuss the most recent theoretical predictions for merger rates. Next, we turn to the quasi-equilibrium formalisms that are used to study binaries prior to the merger phase and to generate initial data for fully dynamical simulations. The quasi-equilibrium approximation has played a key role in developing our understanding of the physics of binary coalescence and, in particular, of the orbital instability processes that can drive binaries to merger at the end of their lifetimes. We then turn to the numerical techniques used in dynamical simulations, including relativistic formalisms, (magneto-)hydrodynamics, gravitational-wave extraction techniques, and nuclear microphysics treatments. This is followed by a summary of the simulations performed across the field to date, including the most recent results from both fully relativistic and microphysically detailed simulations. Finally, we discuss the likely directions for the field as we transition from the first to the second generation of gravitational-wave interferometers and while supercomputers reach the petascale frontier.

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.

Motivated by suggestions that binaries with almost equal-mass components ("twins") play an important role in the formation of double neutron stars and may be rather abundant among binaries, we study the stability of synchronized close and contact binaries with identical components in circular orbits. In particular, we investigate the dependency of the innermost stable circular orbit on the core mass, and we study the coalescence of the binary that occurs at smaller separations. For twin binaries composed of convective main-sequence stars, subgiants, or giants with low mass cores (M_c <~0.15M, where M is the mass of a component), a secular instability is reached during the contact phase, accompanied by a dynamical mass transfer instability at the same or at a slightly smaller orbital separation. Binaries that come inside this instability limit transfer mass gradually from one component to the other and then coalesce quickly as mass is lost through the outer Lagrangian points. For twin giant binaries with moderate to massive cores (M_c >~0.15M), we find that stable contact configurations exist at all separations down to the Roche limit, when mass shedding through the outer Lagrangian points triggers a coalescence of the envelopes and leaves the cores orbiting in a central tight binary. In addition to the formation of binary neutron stars, we also discuss the implications of our results for the production of planetary nebulae with double degenerate central binaries.

About 25 per cent of `hot Jupiters' (extrasolar Jovian-mass planets with close-in orbits) are actually orbiting counter to the spin direction of the star. Perturbations from a distant binary star companion can produce high inclinations, but cannot explain orbits that are retrograde with respect to the total angular momentum of the system. Such orbits in a stellar context can be produced through secular (that is, long term) perturbations in hierarchical triple-star systems. Here we report a similar analysis of planetary bodies, including both octupole-order effects and tidal friction, and find that we can produce hot Jupiters in orbits that are retrograde with respect to the total angular momentum. With distant stellar mass perturbers, such an outcome is not possible. With planetary perturbers, the inner orbit's angular momentum component parallel to the total angular momentum need not be constant. In fact, as we show here, it can even change sign, leading to a retrograde orbit. A brief excursion to very high eccentricity during the chaotic evolution of the inner orbit allows planet-star tidal interactions to rapidly circularize that orbit, decoupling the planets and forming a retrograde hot Jupiter.

We model a circular mass-transferring binary system to calculate the exchange of angular momentum between stellar spins and the orbit due to direct impact of the mass transfer stream onto the surface of the accretor. We simulate mass transfer by calculating the ballistic motion of a point mass ejected from the $L_1$ point of the donor star, conserving the total linear and angular momentum of the system, and treating the stars as uniform density spheres with main sequence radii determined by their masses. We show that, contrary to previous assumptions in the literature, direct impact does not always act as a sink of orbital angular momentum and may in fact increase it by facilitating the transfer of angular momentum from the spin of the donor to the orbit. Here, we show an example of the exchange of angular momentum, as well as a measure of the orbital angular momentum changes for a variety of binary star systems with main sequence components.

Recent observations have revealed two new classes of planetary orbits. Rossiter- Mclaughlin (RM) measurements have revealed hot Jupiters in high-obliquity orbits. In addition, direct-imaging has discovered giant planets at large (~ 100 AU) separations via direct-imaging technique. Simple-minded disk-migration scenarios are inconsistent with the high-inclination (and even retrograde) orbits as seen in recent RM measurements. Furthermore, forming giant planets at large semi-major axis (a) may be challenging in the core-accretion paradigm. We perform many N-body simulations to explore the two above-mentioned orbital architectures. Planet-planet scattering in a multi-planet system can naturally excite orbital inclinations. Planets can also get scattered to large distances. Large-a planetary orbits created from planet-planet scattering are expected to have high eccentricities (e). Theoretical models predict that the observed long-period planets, such as Fomalhaut-b have moderate e \approx 0.3. Interestingly, these are also in systems with disks. We find that if a massive-enough outer disk is present, a scattered planet may be circularized at large a via dynamical friction from the disk and repeated scattering of the disk particles.

Recent discoveries of several transiting planets with clearly non-zero eccentricities and some large inclinations started changing the simple picture of close-in planets having circular and well-aligned orbits. Two major scenarios to form such planets are planet migration in a disk, and planet--planet interactions combined with tidal dissipation. The former scenario can naturally produce a circular and low-obliquity orbit, while the latter implicitly assumes an initially highly eccentric and possibly high-obliquity orbit, which are then circularized and aligned via tidal dissipation. We investigate the tidal evolution of transiting planets on eccentric orbits. We show that the current and future orbital evolution of these systems is likely dominated by tidal dissipation, and not by a more distant companion. Although most of these close-in planets experience orbital decay all the way to the Roche limit, there are two characteristic evolution paths for them, depending on the relative efficiency of tidal dissipation inside the star and the planet. We point out that the current observations may be consistent with one of them. Our results suggest that at least some of the close-in planets with non-zero orbital eccentricity may have been formed by tidally circularizing an initially eccentric orbit. We also find that even when the stellar spin-orbit misalignment is observed to be small at present, some systems could have had a highly misaligned orbit in the past. Finally, we also re-examine the recent claim by Levrard et. al., who found that all orbital and spin parameters evolve on a similar timescale to orbital decay.

The rapid circularization and synchronization of the stellar components in an eccentric binary system at the onset of Roche lobe overflow (RLO) is a fundamental assumption common to all binary stellar evolution and population synthesis codes, even though the validity of this assumption is questionable both theoretically and observationally. Here we calculate the evolution of the orbital elements of an eccentric binary through the direct three-body integration of a massive particle ejected through the inner Lagrangian point of the donor star at periastron. The trajectory of this particle leads to three possible outcomes: direct accretion (DA) onto the companion star within a single orbit, self-accretion (SA) back onto the donor star within a single orbit, or a quasi-periodic orbit around the companion star. We calculate the secular evolution of the binary orbit in the first two cases and conclude that DA can increase or decrease the orbital semi-major axis and eccentricity, while SA always decreases the orbital both orbital elements. In cases where mass overflow contributes to circularizing the orbit, circularization can set in on timescales as short as a few per cent of the mass transfer timescale. In cases where mass overflow increases the eccentricity, the orbital evolution is governed by competition between mass overflow and tidal torques. In the absence of tidal torques, mass overflow resulting in DI can lead to substantially subsynchronously rotating donor stars. Contrary to common assumptions, DI furthermore does not always provide a strong sink of orbital angular momentum in close mass-transferring binaries; in fact we instead find that a significant part can be returned to the orbit during the particle orbit. The formulation presented here can be combined with stellar and binary evolution codes to generate a better picture of the evolution of eccentric, RLO binary star systems.

The discovery of over 400 extrasolar planets allows us to statistically test our understanding of formation and dynamics of planetary systems via numerical simulations. Traditional N-body simulations of multiple-planet systems without gas disks have successfully reproduced the eccentricity (e) distribution of the observed systems, by assuming that the planetary systems are relatively closely packed when the gas disk dissipates, so that they become dynamically unstable within the stellar lifetime. However, such studies cannot explain the small semi-major axes (a) of extrasolar planetary systems, if planets are formed, as the standard planet formation theory suggests, beyond the ice line. In this paper, we numerically study the evolution of three-planet systems in dissipating gas disks, and constrain the initial conditions that reproduce the observed semi-major axis and eccentricity distributions simultaneously. We adopt the initial conditions that are motivated by the standard planet formation theory, and self-consistently simulate the disk evolution, and planet migration by using a hybrid N-body and 1D gas disk code. We also take account of eccentricity damping, and investigate the effect of saturation of corotation resonances on the evolution of planetary systems. We find that the semi-major axis distribution is largely determined in a gas disk, while the eccentricity distribution is determined after the disk dissipation. We also find that there may be an optimum disk mass which leads to the observed a-e distribution. Our simulations generate a larger fraction of planetary systems trapped in mean-motion resonances (MMRs) than the observations, indicating that the disk's perturbation to the planetary orbits may be important to explain the observed rate of MMRs. We also find much lower occurrence of planets on retrograde orbits than the current observations of close-in planets suggest.

The long-term habitability of Earth-like planets requires low orbital eccentricities. A secular perturbation from a distant stellar companion is a very important mechanism in exciting planetary eccentricities, as many of the extrasolar planetary systems are associated with stellar companions. Although the orbital evolution of an Earth-like planet in a stellar binary is well understood, the effect of a binary perturbation to a more realistic system containing additional gas giant planets has been very little studied. Here we provide analytic criteria confirmed by a large ensemble of numerical integrations that identify the initial orbital parameters leading to eccentric orbits. We show that an extra-solar earth is likely to experience a broad range of orbital evolution dictated by the location of a gas-giant planet, necessitating more focused studies on the effect of eccentricity on the potential for life.

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.

We study the dynamical evolution of the young star cluster Arches and its dependence on the assumed initial stellar mass function (IMF). We perform many direct $N$-body simulations with various initial conditions and two different choices of IMFs. One is a standard Kroupa IMF without any mass segregation. The other is a radially dependent IMF, as presently observed in the Arches. We find that it is unlikely for the Arches to have attained the observed degree of mass segregation at its current age starting from a standard non-segregated Kroupa IMF. We also study the possibility of a collisional runaway developing in the first $\sim 2-3 \rm{Myr}$ of dynamical evolution. We find that the evolution of this cluster is dramatically different depending on the choice of IMF: if a primordially mass segregated IMF is chosen, a collisional runaway should always occur between $2-3 \rm{Myr}$ for a broad range of initial concentrations. In contrast, for a standard Kroupa IMF no collisional runaway is predicted. We argue that if Arches was created with a mass segregated IMF similar to what is observed today then at the current cluster age a very unusual, high-mass star should be created. However, whether a collisional runaway leads to the formation of an intermediate-mass black hole (IMBH) depends strongly on the mass loss rate via winds from massive stars. Growth of stellar mass through collisions can be quenched by strong wind mass loss. In that case, the inter-cluster as well as intra-cluster medium are expected to have a significant Helium enrichment which may be observed via Helium recombination lines. The excess amount of gas lost in winds may also be observed via X-ray observations as diffused X-ray sources.

Planets have been observed in tight binary systems with separations less than 20 AU. A likely formation scenario for such systems involves a dynamical capture, after which high relative inclinations are likely and may lead to Kozai oscillations. We numerically investigate the fate of an initially coplanar double-planet system in a class of binaries with separation ranging between $12 - 20 $ AU. Dynamical integrations of representative four-body systems are performed, each including a hot Jupiter and a second planet on a wider orbit. We find that, although such systems can remain stable at low relative inclinations ($\lesssim 40^\circ$), high relative inclinations are likely to lead to instabilities. This can be avoided if the planets are placed in a \emph{Kozai-stable zone} within which mutual gravitational perturbations can suppress the Kozai mechanism. We investigate the possibility of inducing Kozai oscillations in the inner orbit by a weak coupling mechanism between the planets in which the coplanarity is broken due to a differential nodal precession. Propagating perturbations from the stellar companion through a planetary system in this manner can have dramatic effects on the dynamical evolution of planetary systems, especially in tight binaries and can offer a reasonable explanation for eccentricity trends among planets observed in binary systems. We find that inducing such oscillations into the orbit of a hot Jupiter is more likely in tight binaries and an upper limit can be set on the binary separation above which these oscillations are not observed.

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.

We investigate the secular evolution of the orbital semi-major axis and eccentricity due to mass transfer in eccentric binaries, allowing for both mass and angular momentum loss from the system. Adopting a delta function mass transfer rate at the periastron of the binary orbit, we find that, depending on the initial binary properties at the onset of mass transfer, the orbital semi-major axis and eccentricity can either increase or decrease at a rate linearly proportional to the magnitude of the mass transfer rate at periastron. The range of initial binary mass ratios and eccentricities that leads to increasing orbital semi-major axes and eccentricities broadens with increasing degrees of mass loss from the system and narrows with increasing orbital angular momentum loss from the binary. Comparison with tidal evolution timescales shows that the usual assumption of rapid circularization at the onset of mass transfer in eccentric binaries is not justified, irrespective of the degree of systemic mass and angular momentum loss. This work extends our previous results for conservative mass transfer in eccentric binaries and can be incorporated into binary evolution and population synthesis codes to model non-conservative mass transfer in eccentric binaries.

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.

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.

Mean-motion resonances (MMRs) are likely to play an important role both during and after the lifetime of a protostellar gas disk. We study the dynamical evolution and stability of planetary systems containing two giant planets on circular orbits near a 2:1 resonance and closer. We find that by having the outer planet migrate inward, the two planets can capture into either the 2:1, 5:3, or 3:2 MMR. We use direct numerical integrations of ~1000 systems in which the planets are initially locked into one of these resonances and allowed to evolve for up to ~10^7 yr. We find that the final eccentricity distribution in systems which ultimately become unstable gives a good fit to observed exoplanets. Next, we integrate ~500 two-planet systems in which the outer planet is driven to continuously migrate inward, resonantly capturing the inner; the systems are evolved until either instability sets in or the planets reach the star. We find that although the 5:3 resonance rapidly becomes unstable under migration, the 2:1 and 3:2 are very stable. Thus the lack of observed exoplanets in resonances closer than 2:1, if it continues to hold up, may be a primordial signature of the planet formation process.

Strong tidal interaction with the central star can circularize the orbits of close-in planets. With the standard tidal quality factor Q of our solar system, estimated circularization times for close-in extrasolar planets are typically shorter than the ages of the host stars. While most extrasolar planets with orbital radii a < 0.1 AU indeed have circular orbits, some close-in planets with substantial orbital eccentricities have recently been discovered. This new class of eccentric close-in planets implies that either their tidal Q factor is considerably higher, or circularization is prevented by an external perturbation. Here we constrain the tidal Q factor for transiting extrasolar planets by comparing their circularization times with accurately determined stellar ages. Using estimated secular perturbation timescales, we also provide constraints on the properties of hypothetical second planets exterior to the known ones.

The ensemble of now more than 250 discovered planetary systems displays a wide range of masses, orbits and, in multiple systems, dynamical interactions. These represent the end point of a complex sequence of events, wherein an entire protostellar disk converts itself into a small number of planetary bodies. Here, we present self-consistent numerical simulations of this process, which produce results in agreement with some of the key trends observed in the properties of the exoplanets. Analogs to our own solar system do not appear to be common, originating from disks near the boundary between barren and (giant) planet-forming.

Many recent observational studies have concluded that planetary systems commonly exist in multiple-star systems. At least ~20%, and presumably a larger fraction of the known extrasolar planetary systems are associated with one or more stellar companions. These stellar companions normally exist at large distances from the planetary systems (typical projected binary separations are on the orders 100-10000AU) and are often faint (ranging from F to T spectral types). Yet, secular cyclic angular momentum exchange with these distant stellar companions can significantly alter the orbital configuration of the planets around the primaries. One of the most interesting and fairly common outcomes seen in numerical simulations is the opening of a large mutual inclination angle between the planetary orbits, forced by differential nodal precessions caused by the binary companion. The growth of the mutual inclination angle between planetary orbits induces additional large-amplitude eccentricity oscillations of the inner planet due to the quadrupole gravitational perturbation by the outer planet. This eccentricity oscillation may eventually result in the orbital decay of the inner planet through tidal friction, as previously proposed as Kozai migration or Kozai cycles with tidal friction (KCTF). This orbital decay mechanism induced by the binary perturbation and subsequent tidal dissipation may stand as an alternative formation channel for close-in extrasolar planets.

(Abridged) In planetary systems with two or more giant planets, dynamical instabilities can lead to collisions or ejections through strong planet--planet scattering. Previous studies for simple initial configurations with two equal-mass planets revealed some discrepancies between the results of numerical simulations and the observed orbital elements of extrasolar planets. Here, we show that simulations with two unequal mass planets starting on nearly circular orbits predict a reduced frequency of collisions and a broader range of final eccentricities. The two-planet scattering model can easily reproduce the observed eccentricities with a plausible distribution of planet mass ratios. Further, the two-planet scattering model predicts a maximum eccentricity of about 0.8, independent of the distribution of planet mass ratios, provided that both planets are initially place on nearly circular orbits. This compares favorably with current observations and will be tested by future planet discoveries. The combination of planet--planet scattering and tidal circularization may be able to explain the existence of some giant planets with very short period orbits. Orbital migration due to planet scattering could play an important role in explaining the increased rate of giant planets with very short period orbits. We also re-examine and discuss various possible correlations between eccentricities and other properties of observed extrasolar planets. We demonstrate that the observed distribution of planet masses, orbital periods, and eccentricities can provide constraints for models of planet formation and evolution.

Recent Spitzer observations of the globular cluster M15 detected dust associated with its intracluster medium. Surprisingly, these observations imply that the dust must be very short-lived compared to the time since the last Galactic plane crossing of the cluster.Here we propose a simple mechanism to explain this short lifetime. We argue that the kinetic energy of the material ejected during a stellar collision may be sufficient to remove the gas and dust entirely from a cluster, or to remove the gas as a wind, in addition to partially destroying the dust. By calculating the rate of stellar collisions using an N-body model for the cluster, we find remarkable agreement between the average time between collisions and the inferred dust lifetime in this cluster, suggesting a possible close relation between the two phenomena. Furthermore, we also obtain the birthrate of blue stragglers formed through collisions in M15. By comparing with the observed number of blue stragglers, we derive an upper limit for their average lifetime which turns out to be consistent with recent model calculations, thereby lending further support to our model.

Observations in the past decade have revealed extrasolar planets with a wide range of orbital semimajor axes and eccentricities. Based on the present understanding of planet formation via core accretion and oligarchic growth, we expect that giant planets often form in closely packed configurations. While the protoplanets are embedded in a protoplanetary gas disk, dissipation can prevent eccentricity growth and suppress instabilities from becoming manifest. However, once the disk dissipates, eccentricities can grow rapidly, leading to close encounters between planets. Strong planet--planet gravitational scattering could produce both high eccentricities and, after tidal circularization, very short-period planets, as observed in the exoplanet population. We present new results for this scenario based on extensive dynamical integrations of systems containing three giant planets, both with and without residual gas disks. We assign the initial planetary masses and orbits in a realistic manner following the core accretion model of planet formation. We show that, with realistic initial conditions, planet--planet scattering can reproduce quite well the observed eccentricity distribution. Our results also make testable predictions for the orbital inclinations of short-period giant planets formed via strong planet scattering followed by tidal circularization.

Many recent observational studies have concluded that planetary systems commonly exist in multiple-star systems. At least ~20% of the known extrasolar planetary systems are associated with one or more stellar companions. The orbits of stellar binaries hosting planetary systems are typically wider than 100 AU and often highly inclined with respect to the planetary orbits. The effect of secular perturbations from such an inclined binary orbit on a coupled system of planets, however, is little understood theoretically. In this paper we investigate various dynamical classes of double-planet systems in binaries through numerical integrations and we provide an analytic framework based on secular perturbation theories. Differential nodal precession of the planets is the key property that separates two distinct dynamical classes of multiple planets in binaries: (1) dynamically-rigid systems in which the orbital planes of planets precess in concert as if they were embedded in a rigid disk, and (2) weakly-coupled systems in which the mutual inclination angle between initially coplanar planets grows to large values on secular timescales. In the latter case, the quadrupole perturbation from the outer planet induces additional Kozai cycles and causes the orbital eccentricity of the inner planet to oscillate with large amplitudes. The cyclic angular momentum transfer from a stellar companion propagating inward through planets can significantly alter the orbital properties of the inner planet on shorter timescales. This perturbation propagation mechanism may offer important constraints on the presence of additional planets in known single-planet systems in binaries.

Mergers of double neutron stars are considered the most likely progenitors for short gamma-ray bursts. Indeed such a merger can produce a black hole with a transient accreting torus of nuclear matter (Lee & Ramirez-Ruiz 2007, Oechslin & Janka 2006), and the conversion of a fraction of the torus mass-energy to radiation can power a gamma-ray burst (Nakar 2006). Using available binary pulsar observations supported by our extensive evolutionary calculations of double neutron star formation, we demonstrate that the fraction of mergers that can form a black hole -- torus system depends very sensitively on the (largely unknown) maximum neutron star mass. We show that the available observations and models put a very stringent constraint on this maximum mass under the assumption that a black hole formation is required to produce a short gamma-ray burst in a double neutron star merger. Specifically, we find that the maximum neutron star mass must be within 2 - 2.5 Msun. Moreover, a single unambiguous measurement of a neutron star mass above 2.5 Msun would exclude a black hole -- torus central engine model of short gamma-ray bursts in double neutron star mergers. Such an observation would also indicate that if in fact short gamma-ray bursts are connected to neutron star mergers, the gamma-ray burst engine is best explained by the lesser known model invoking a highly magnetized massive neutron star (e.g., Usov 1992; Kluzniak & Ruderman 1998; Dai et al. 2006; Metzger, Quataert & Thompson 2007).

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.

We use a 3-D relativistic SPH (Smoothed Particle Hydrodynamics) code to study mergers of black hole -- neutron star (BH--NS) binary systems with low mass ratios, adopting $M_{NS}/M_{BH} \simeq 0.1$ as a representative case. The outcome of such mergers depends sensitively on both the magnitude of the BH spin and its obliquity (i.e., the inclination of the binary orbit with respect to the equatorial plane of the BH). In particular, only systems with sufficiently high BH spin parameter $a$ and sufficiently low orbital inclinations allow any NS matter to escape or to form a long-lived disk outside the BH horizon after disruption. Mergers of binaries with orbital inclinations above $\sim60^o$ lead to complete prompt accretion of the entire NS by the BH, even for the case of an extreme Kerr BH. We find that the formation of a significant disk or torus of NS material around the BH always requires a near-maximal BH spin and a low initial inclination of the NS orbit just prior to merger.

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.

Using a population synthesis approach, we compute the total merger rate in the local Universe for double neutron stars, double black holes, and black hole -- neutron star binaries. These compact binaries are the prime source candidates for gravitational-wave detection by LIGO and VIRGO. We account for mergers originating from field populations and from dense stellar clusters. For both populations we use the same treatment of stellar evolution. Our results indicate that the merger rates of double neutron stars and black hole -- neutron star binaries are strongly dominated by field populations, while merging black hole binaries are formed much more effectively in dense stellar clusters. The overall merger rate of double compact objects depends sensitively on the (largely unknown) initial mass fraction contained in dense clusters f_cl. For f_cl < 0.0001, the Advanced LIGO detection rate will be dominated by field populations of double neutron star mergers, with a small but significant number of detections ~20 yr^-1. However for a higher mass fraction in clusters, f_cl > 0.001, the detection rate will be dominated by numerous mergers of double black holes originating from dense clusters, and it will be considerably higher, ~25 - 300 yr^-1. In addition, we show that, once mergers of double black holes are detected, it is easy to differentiate between systems formed in the field and in dense clusters, since the chirp mass distributions are strikingly different. Finally, we point out that there may exist a population of merging black hole binaries in intergalactic space.

Dynamical interactions that occur between objects in dense stellar systems are particularly important for the question of formation of X-ray binaries. We present results of numerical simulations of 70 globular clusters with different dynamical properties and a total stellar mass of 2*10^7 Msun. We find that in order to retain enough neutron stars to match observations we must assume that NSs can be formed via electron-capture supernovae. Our simulations explain the observed dependence of the number of LMXBs on ``collision number'' as well as the large scatter observed between different globular clusters. For millisecond pulsars, we obtain good agreement between our models and the numbers and characteristics of observed pulsars in the clusters Terzan 5 and 47 Tuc

In this contribution we discuss how neutron stars are produced and retained in globular clusters, outlining the most important dynamical channels and evolutionary events that affect thepopulation of mass-transferring binaries with neutron stars and result in the formation of recycled pulsars. We confirm the importance of electron-capture supernovae in globular clusters as the major supplier of retained neutron stars.By comparing the observed millisecond pulsar population and the results obtained from simulations, we discuss several constraints on the evolution of mass-transferring systems.In particular, we find that in our cluster model the following mass-gaining events create populations of MSPs that do not match the observations (with respect to binary periods and companion masses or the number of produced systems) and therefore likely do not lead to NSs spun up to millisecond periods: (i) accretion during a common envelope event with a NS formed through accretion-induced collapse, and (ii) mass transfer from a WD donor. By restricting ourselves to the evolutionary and dynamical paths that most likely lead to neutron star recycling, we obtain good agreement between our models and the numbers and characteristics of observed millisecond pulsars in the clusters Terzan 5 and 47 Tuc.

The orbital eccentricity of a single planet around a component of a stellar binary system with a sufficiently large mutual inclination angle is known to oscillate on a secular timescale through the Kozai mechanism. We have investigated the effects of the Kozai mechanism on double-planet systems in binaries. The evolutionary sequence of a pair of planets under the influence of a binary companion is fairly complex. Various dynamical outcomes are seen in numerical simulations. One interesting outcome is the rigid rotation of the planetary orbits in which the planetary orbital planes secularly precess in concert, while the orbital eccentricities oscillate synchronously. In such cases the outer planet acts as a propagator of the perturbation from the binary companion to the inner planet and drives the inner planetary orbit to precess at a rate faster than what is predicted by the Kozai mechanism.

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.

Binary black holes can form efficiently in dense young stellar clusters, such as the progenitors of globular clusters, via a combination of gravitational segregation and cluster evaporation. We use simple analytic arguments supported by detailed N-body simulations to determine how frequently black holes born in a single stellar cluster should form binaries, be ejected from the cluster, and merge through the emission of gravitational radiation. We then convolve this ``transfer function'' relating cluster formation to black hole mergers with (i) the distribution of observed cluster masses and (ii) the star formation history of the universe, assuming that a significant fraction gcl of star formation occurs in clusters and that a significant fraction gcand of clusters undergo this segregation and evaporation process. We predict future ground--based gravitational wave (GW) detectors could observe ~500 (gcl/0.5) (gcand/0.1) double black hole mergers per year, and the presently operating LIGO interferometer would have a chance (50%) at detecting a merger during its first full year of science data. More realistically, advanced LIGO and similar next-generation gravitational wave observatories provide unique opportunities to constrain otherwise inaccessible properties of clusters formed in the early universe.

We investigate the secular evolution of the orbital semi-major axis and eccentricity due to mass transfer in eccentric binaries, assuming conservation of total system mass and orbital angular momentum. Assuming a delta function mass transfer rate centered at periastron, we find rates of secular change of the orbital semi-major axis and eccentricity which are linearly proportional to the magnitude of the mass transfer rate at periastron. The rates can be positive as well as negative, so that the semi-major axis and eccentricity can increase as well as decrease in time. Adopting a delta-function mass-transfer rate of $10^{-9} M_\sun {\rm yr}^{-1}$ at periastron yields orbital evolution timescales ranging from a few Myr to a Hubble time or more, depending on the binary mass ratio and orbital eccentricity. Comparison with orbital evolution timescales due to dissipative tides furthermore shows that tides cannot, in all cases, circularize the orbit rapidly enough to justify the often adopted assumption of instantaneous circularization at the onset of mass transfer. The formalism presented can be incorporated in binary evolution and population synthesis codes to create a self-consistent treatment of mass transfer in eccentric binaries.

We show that interaction with a gas disk may produce young planetary systems with closely-spaced orbits, stabilized by mean-motion resonances between neighbors. On longer timescales, after the gas is gone, interaction with a remnant planetesimal disk tends to pull these configurations apart, eventually inducing dynamical instability. We show that this can lead to a variety of outcomes; some cases resemble the Solar System, while others end up with high-eccentricity orbits reminiscent of the observed exoplanets. A similar mechanism has been previously suggested as the cause of the lunar Late Heavy Bombardment. Thus, it may be that a large-scale dynamical instability, with more or less cataclysmic results, is an evolutionary step common to many planetary systems, including our own.

We present a comprehensive description of the population synthesis code StarTrack. The original code has been significantly modified and updated. Special emphasis is placed here on processes leading to the formation and further evolution of compact objects (white dwarfs, neutron stars, and black holes). Both single and binary star populations are considered. The code now incorporates detailed calculations of all mass-transfer phases, a full implementation of orbital evolution due to tides, as well as the most recent estimates of magnetic braking. This updated version of StarTrack can be used for a wide variety of problems, with relevance to many current and planned observatories, e.g., studies of X-ray binaries (Chandra, XMM-Newton), gravitational radiation sources (LIGO, LISA), and gamma-ray burst progenitors (HETE-II, Swift). The code has already been used in studies of Galactic and extra-galactic X-ray binary populations, black holes in young star clusters, Type Ia supernova progenitors, and double compact object populations. Here we describe in detail the input physics, we present the code calibration and tests, and we outline our current studies in the context of X-ray binary populations.

(Abridged) The discovery of extrasolar planetary systems revealed an unexpected diversity of planetary systems that has revolutionized planet formation theory. A strong program of theoretical research is essential to maximize both the discovery potential and the scientific returns of future observational programs, so as to achieve a deeper understanding of the formation and evolution of planetary systems. We outline three broad categories of theoretical research: detailed studies of specific planetary systems, testing planet formation models by comparing their predictions to the observed exoplanet population, and detailed modeling of specific physical processes. We describe how such theoretical research plays an important role in analyzing observations for a wide range detection methods and contributes to understanding the Earth's place in the universe and the potential for Earth-like life beyond our solar system. In this white paper, we suggest how to maximize the scientific return of future exoplanet observations. Our recommendations include a strong theory program, support for multiple observational programs that will study a diverse set of planets and stars, significant observing time devoted to follow-up observations, and healthy collaboration between observers and theorists.

Double black hole binaries are among the most important sources of gravitational radiation for ground-based detectors such as LIGO or VIRGO. Even if formed with lower efficiency than double neutron star binaries, they could dominate the predicted detection rates, since black holes are more massive than neutron stars and therefore could be detected at greater distances. Here we discuss an evolutionary process that can very significantly limit the formation of close double black hole binaries: the vast majority of their potential progenitors undergo a common envelope (CE) phase while the donor, one of the massive binary components, is evolving through the Hertzsprung gap. Our latest theoretical understanding of the CE process suggests that this will probably lead to a merger, inhibiting double black hole formation. Barring uncertainties in the physics of CE evolution, we use population synthesis calculations, and find that the corresponding reduction in the merger rate of double black holes formed in galactic fields is so great (by ~500) that their contribution to inspiral detection rates for ground-based detectors could become relatively small (1 in 10) compared to double neutron star binaries. A similar process also reduces the merger rates for double neutron stars, by factor of ~5, eliminating most of the previously predicted ultracompact NS-NS systems. Our predicted detection rates for Advanced LIGO are now much lower for double black holes (~2/yr), but are still quite high for double neutron stars (~20/yr). If double black holes were found to be dominant in the detected inspiral signals, this could indicate that they mainly originate from dense star clusters (not included here) or that our theoretical understanding of the CE phase requires significant revision.

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.

We derive detailed theoretical models for 1074 nearby stars from the SPOCS (Spectroscopic Properties of Cool Stars) Catalog. The California and Carnegie Planet Search has obtained high-quality echelle spectra of over 1000 nearby stars taken with the Hamilton spectrograph at Lick Observatory, the HIRES spectrograph at Keck, and UCLES at the Anglo Australian Observatory. A uniform analysis of the high-resolution spectra has yielded precise stellar parameters, enabling systematic error analyses and accurate theoretical stellar modeling. We have created a large database of theoretical stellar evolution tracks using the Yale Stellar Evolution Code (YREC) to match the observed parameters of the SPOCS stars. Our very dense grids of evolutionary tracks eliminate the need for interpolation between stellar evolutionary tracks and allow precise determinations of physical stellar parameters (mass, age, radius, size and mass of the convective zone, etc.). Combining our stellar models with the observed stellar atmospheric parameters and uncertainties, we compute the likelihood for each set of stellar model parameters separated by uniform time steps along the stellar evolutionary tracks. The computed likelihoods are used for a Bayesian analysis to derive posterior probability distribution functions for the physical stellar parameters of interest. We provide a catalog of physical parameters for 1074 stars that are based on a uniform set of high quality spectral observations, a uniform spectral reduction procedure, and a uniform set of stellar evolutionary models. We explore this catalog for various possible correlations between stellar and planetary properties, which may help constrain the formation and dynamical histories of other planetary systems.

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.

We have calculated the first dynamical evolutions of merging black hole-neutron star binaries that treat the combined spacetime in a nonperturbative general relativistic framework. Using the conformal flatness approximation, we have studied how the location of the tidal disruption radius with respect to the the black hole horizon and innermost stable circular orbit (ISCO) affects the qualitative evolution of the system. Based on simple arguments, we show that for a binary mass ratio q>~0.24, tidal disruption occurs outside the ISCO, while the opposite is true for q<~0.24. When tidal disruption occurs sufficiently far outside the ISCO, mass is transferred unstably from the neutron star to the black hole, resulting in the complete disruption of the neutron star. When tidal disruption occurs slightly within the ISCO, we find that some of the mass forms an extremely hot disk around the black hole. The resulting configurations in this case are excellent candidates for the progenitors of short-hard gamma ray bursts.

Using an updated population synthesis code we study the formation and evolution of black holes (BHs) in young star clusters following a massive starburst. This study continues and improves on the initial work described by Belczynski, Sadowski & Rasio (2004). In our new calculations we account for the possible ejections of BHs and their progenitors from clusters because of natal kicks imparted by supernovae and recoil following binary disruptions. The results indicate that the properties of both retained BHs in clusters and ejected BHs (forming a field population) depend sensitively on the depth of the cluster potential. In particular, most BHs ejected from binaries are also ejected from clusters with central escape speeds Vesc < 100 km/s. Conversely, most BHs remaining in binaries are retained by clusters with Vesc > 50 km/s. BHs from single star evolution are also affected significantly: about half of the BHs originating from primordial single stars are ejected from clusters with Vesc < 50 km/s. Our results lay a foundation for theoretical studies of the formation of BH X-ray binaries in and around star clusters, including possible ultra-luminous sources, as well as merging BH--BH binaries detectable with future gravitational-wave observatories.

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.

We present a new study of the collisional runaway scenario to form an intermediate-mass black hole (IMBH, MBH > 100 Msun) at the centre of a young, compact stellar cluster. The first phase is the formation of a very dense central core of massive stars (Mstar =~ 30-120 Msun) through mass segregation and gravothermal collapse. Previous work established the conditions for this to happen before the massive stars evolve off the main sequence (MS). In this and a companion paper, we investigate the next stage by implementing direct collisions between stars. Using a Monte Carlo stellar dynamics code, we follow the core collapse and subsequent collisional phase in more than 100 models with varying cluster mass, size, and initial concentration. Collisions are treated either as ideal, ``sticky-sphere'' mergers or using realistic prescriptions derived from 3-D hydrodynamics computations. In all cases for which the core collapse happens in less than the MS lifetime of massive stars (~3 Myr), we obtain the growth of a single very massive star (VMS, Mstar =~ 400-4000 Msun) through a runaway sequence of mergers. Mass loss from collisions, even for velocity dispersions as high as sigma1D ~ 1000 km/s, does not prevent the runaway. The region of cluster parameter space leading to runaway is even more extended than predicted in previous work because, in clusters with sigma1D > 300 km/s, collisions accelerate (and, in extreme cases, drive) core collapse. Although the VMS grows rapidly to > 1000 Msun in models exhibiting runaway, we cannot predict accurately its final mass. This is because the termination of the runaway process must eventually be determined by a complex interplay between stellar dynamics, hydrodynamics, and the stellar evolution of the VMS. [abridged]

We present the methods and preparatory work for our study of the collisional runaway scenario to form a very massive star (VMS, M > 400 Msun) at the centre of a young, compact stellar cluster. In the first phase of the process, a very dense central core of massive stars (M =~ 30-120 Msun) forms through mass segregation and gravothermal collapse. This leads to a collisional stage, likely to result in the formation of a VMS (itself a possible progenitor for an intermediate-mass black hole) through a runaway sequence of mergers between the massive stars. In this paper we present the runaway scenario in a general astrophysical context. We then explain the numerical method used to investigate it. Our approach is based on a Monte Carlo code to simulate the stellar dynamics of spherical star clusters using a very large number of particles (a few 1e5 to several 1e6). Finally, we report on test computations carried out to ensure that our implementation of the important physics is sound. In a second paper, we present results from more than 100 cluster simulations realized to determine the conditions leading to the collisional formation of a VMS and the characteristics of the runaway sequences.

Many of the known extrasolar planets are ``hot Jupiters,'' giant planets with orbital periods of just a few days. We use the observed distribution of hot Jupiters to constrain the location of its inner edge in the mass--period diagram. If we assume a slope corresponding to the classical Roche limit, then we find that the edge corresponds to a separation close to_twice_ the Roche limit, as expected if the planets started on highly eccentric orbits that were later circularized. In contrast, any migration scenario would predict an inner edge right at the Roche limit, which applies to planets approaching on nearly circular orbits. However, the current sample of hot Jupiters is not sufficient to provide a precise constraint simultaneously on both the location and slope of the inner edge.

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.

(abridged) We report the results of SPH calculations of parabolic collisions between a subgiant or slightly evolved red-giant star and a neutron star (NS). Such collisions are likely to form ultracompact X-ray binaries (UCXBs) observed today in old globular clusters. In particular, we compute collisions of a 1.4 Msun NS with realistically modelled parent stars of initial masses 0.8 and 0.9 Msun, each at three different evolutionary stages (corresponding to three different radii R). The distance of closest approach for the initial orbit varies from 0.04 R (nearly head-on) to 1.3 R (grazing). These collisions lead to the formation of a tight binary, composed of the NS and the subgiant or red-giant core, embedded in an extremely diffuse common envelope (CE) typically of mass ~0.1 to 0.3 Msun. Our calculations follow the binary for many hundreds of orbits, ensuring that the orbital parameters we determine at the end of the calculations are close to final. Some of the fluid initially in the envelope of the (sub)giant, from 0.003 to 0.023 Msun in the cases we considered, is left bound to the NS. The eccentricities of the resulting binaries range from about 0.2 for our most grazing collision to about 0.9 for the nearly head-on cases. In almost all the cases we consider, gravitational radiation alone will cause sufficiently fast orbital decay to form a UCXB within a Hubble time, and often on a much shorter timescale. Our hydrodynamics code implements the recent SPH equations of motion derived with a variational approach by Springel & Hernquist and by Monaghan. Numerical noise is reduced by enforcing an analytic constraint equation that relates the smoothing lengths and densities of SPH particles. We present tests of these new methods to help demonstrate their improved accuracy.

We calculate the first dynamical evolutions of merging black hole-neutron star binaries that construct the combined black hole-neutron star spacetime in a general relativistic framework. We treat the metric in the conformal flatness approximation, and assume that the black hole mass is sufficiently large compared to that of the neutron star so that the black hole remains fixed in space. Using a spheroidal spectral methods solver, we solve the resulting field equations for a neutron star orbiting a Schwarzschild black hole. The matter is evolved using a relativistic, Lagrangian, smoothed particle hydrodynamics (SPH) treatment. We take as our initial data recent quasiequilibrium models for synchronized neutron star polytropes generated as solutions of the conformal thin-sandwich (CTS) decomposition of the Einstein field equations. We are able to construct from these models relaxed SPH configurations whose profiles show good agreement with CTS solutions. Our adiabatic evolution calculations for neutron stars with low compactness show that mass transfer, when it begins while the neutron star orbit is still outside the innermost stable circular orbit, is more unstable than is typically predicted by analytical formalisms. This dynamical mass loss is found to be the driving force in determining the subsequent evolution of the binary orbit and the neutron star, which typically disrupts completely within a few orbital periods. The majority of the mass transferred onto the black hole is accreted promptly; a significant fraction (~30%) of the mass is shed outward as well, some of which will become gravitationally unbound and ejected completely from the system. The remaining portion forms an accretion disk around the black hole, and could provide the energy source for short-duration gamma ray bursts.

The most puzzling property of the extrasolar planets discovered by recent radial velocity surveys is their high orbital eccentricities, which are very difficult to explain within our current theoretical paradigm for planet formation. Current data reveal that at least 25% of these planets, including some with particularly high eccentricities, are orbiting a component of a binary star system. The presence of a distant companion can cause significant secular perturbations in the orbit of a planet. At high relative inclinations, large-amplitude, periodic eccentricity perturbations can occur. These are known as "Kozai cycles" and their amplitude is purely dependent on the relative orbital inclination. Assuming that every planet host star also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, we determine the resulting eccentricity distribution of planets and compare it to observations? We find that perturbations from a binary companion always appear to produce an excess of planets with both very high (e>0.6) and very low (e<0.1) eccentricities. The paucity of near-circular orbits in the observed sample implies that at least one additional mechanism must be increasing eccentricities. On the other hand, the overproduction of very high eccentricities observed in our models could be combined with plausible circularization mechanisms (e.g., friction from residual gas) to create more planets with intermediate eccentricities (e=0.1-0.6).

Many of the known extrasolar planets are ``hot Jupiters,'' giant planets with orbital periods of just a few days. We use the observed distribution of hot Jupiters to constrain the location of the ``inner edge'' and planet migration theory. If we assume the location of the inner edge is proportional to the Roche limit, then we find that this edge is located near twice the Roche limit, as expected if the planets were circularized from a highly eccentric orbit. If confirmed, this result would place significant limits on migration via slow inspiral. However, if we relax our assumption for the slope of the inner edge, then the current sample of hot Jupiters is not sufficient to provide a precise constraint on both the location and power law index of the inner edge.

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.

One of the most remarkable properties of extrasolar planets is their high orbital eccentricities. Observations have shown that at least 20% of these planets, including some with particularly high eccentricities, are orbiting a component of a wide binary star system. The presence of a distant binary companion can cause significant secular perturbations to the orbit of a planet. In particular, at high relative inclinations, a planet can undergo a large-amplitude eccentricity oscillation. This so-called "Kozai mechanism" is effective at a very long range, and its amplitude is purely dependent on the relative orbital inclination. In this paper, we address the following simple question: assuming that every host star with a detected giant planet also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, how well could secular perturbations reproduce the observed eccentricity distribution of planets? Our calculations show that the Kozai mechanism consistently produces an excess of planets with very high (e >0.6) and very low (e < 0.1) eccentricities. The paucity of near-circular orbits in the observed sample cannot be explained solely by the Kozai mechanism, because, even with high enough inclinations, the Kozai mechanism often fails to produce significant eccentricity perturbations when there are other competing sources of orbital perturbations on secular timescales, such as general relativity. On the other hand, the Kozai mechanism can produce many highly eccentric orbits. Indeed the overproduction of high eccentricities observed in our models could be combined with plausible circularizing mechanisms (e.g., friction from residual gas) to create more intermediate eccentricities (e=0.1-0.6).

We present results from dynamical Monte Carlo simulations of dense star clusters near the Galactic center. These clusters sink toward the center of the Galaxy by dynamical friction. During their inspiral, they may undergo core collapse and form an intermediate-mass black hole through runaway collisions. Such a cluster can then reach within a parsec of the Galactic center before it completely disrupts, releasing many young stars in this region. This scenario provides a natural explanation for the presence of the young stars observed near the Galactic center. Here we determine the initial conditions for this scenario to work and we derive the mass distribution of cluster stars as a function of distance from the Galactic center. For reasonable initial conditions, we find that clusters massive enough for rapid inspiral would include more massive stars (m >~ 30 Msun) than currently observed in the inspiral region. We point out several possible explanations for this apparent discrepancy.

We review recent progress in understanding the hydrodynamics of compact binary mergers using relativistic smoothed particle hydrodynamics (SPH) codes. Recent results are discussed for both double neutron stars and black hole - neutron star binaries. The gravitational wave signals from these mergers contain a lot of information about the masses and radii of the stars and the equation of state of the neutron star fluid. For black hole - neutron star mergers, the final fate of the fluid and the gravitational wave emission depend strongly on the mass ratio and the black hole spin.

Doppler spectroscopy has detected 136 planets around nearby stars. A major puzzle is why their orbits are highly eccentric, while all planets in our Solar System are on nearly circular orbits, as expected if they formed by accretion processes in a protostellar disk. Several mechanisms have been proposed to generate large eccentricities after planet formation, but so far there has been little observational evidence to support any particular one. Here we report that the current orbital configuration of the three giant planets around Upsilon Andromedae (ups And) provides evidence for a close dynamical interaction with another planet, now lost from the system. The planets started on nearly circular orbits, but chaotic evolution caused the outer planet (ups And d) to be perturbed suddenly into a higher-eccentricity orbit. The coupled evolution of the system then causes slow periodic variations in the eccentricity of the middle planet (ups And c). Indeed, we show that ups And c periodically returns to a very nearly circular state every 9000 years. Our analysis shows that strong planet-planet scattering, one of several mechanisms previously discussed for increasing orbital eccentricities, must have operated in this system.

Bright, ultracompact X-ray binaries observed in dense star clusters, such as Galactic globular clusters, must have formed relatively recently, since their lifetimes as persistent bright sources are short (e.g., ~10^8 yr above 10^36 erg/s for a 1.4 Msun neutron star accreting from a degenerate helium companion with an initial mass of ~0.2 Msun). Therefore, we can use the present conditions in a cluster core to study possible dynamical formation processes for these sources. Here we show that direct physical collisions between neutron stars and red giants can provide a sufficient formation rate to explain the observed numbers of bright sources. These collisions produce tight, eccentric neutron star -- white dwarf binaries that decay to contact by gravitational radiation on timescales ~10^6-10^10 yr, usually shorter and often much shorter than the cluster age.

More than 100 radio pulsars have been detected in 24 globular clusters. The largest observed samples are in Terzan 5 and 47 Tucanae, which together contain 45 pulsars. Accurate timing solutions, including positions in the cluster, are known for many of these pulsars. Here we provide an observational overview of some properties of pulsars in globular clusters, as well as properties of the globular clusters with detected pulsars. The many recent detections also provide a new opportunity to re-examine theoretically the formation and evolution of recycled pulsars in globular clusters. Our brief review considers the most important dynamical interaction and binary evolution processes: collisions, exchange interactions, mass transfer, and common-envelope phases.

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.

The population of compact binaries in dense stellar systems is affected strongly by frequent dynamical interactions between stars and their interplay with the stellar evolution. In this contribution, we consider these effects on binaries with a white dwarf accretor, in particular cataclysmic variables and AM CVns. We examine which processes can successfully lead to the creation of such X-ray binaries. Using numerical simulations, we identify predominant formation channels and predict the expected numbers of detectable systems. We discuss also why the distribution of cataclysmic variables has a weaker dependence upon the cluster density than the distribution of quiescent low-mass X-ray binaries and why dwarf nova outbursts may not occur among globular cluster cataclysmic variables.

We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be stripped from it, but in no case do we find enough energy transferred to the planet to lead to complete disruption. In light of the three new extrasolar planets discovered with periods shorter than two days, we argue that the shortest-period cases observed in the period-mass relation may be explained by a model whereby planets undergo strong tidal encounters with stars, after either being scattered by dynamical interactions into highly eccentric orbits, or tidally captured from nearly parabolic orbits. Although this scenario does provide a natural explanation for the edge found for planets at twice the Roche limit, it does not explain how such planets will survive the inevitable expansion that results from energy injection during tidal circularization.

Using Monte Carlo codes, we follow the collisional evolution of clusters in a variety of scenarios. We consider the conditions under which a cluster of main sequence stars may undergo rapid core collapse due to mass segregation, thus entering a phase of runaway collisions, forming a very massive star (VMS, M >= 1000 Msun) through repeated collisions between single stars. Although collisional mass loss is accounted for realistically, we find that a VMS forms even in proto-galactic nuclei models with a high velocity dispersion (many 100 km/s). Such a VMS may be a progenitor for an intermediate-mass black hole (M >= 100 Msun). In contrast, in galactic nuclei hosting a central massive black hole, collisions are found to be disruptive. The stars which are subject to collisions are progressively ground down by high-velocity collisions and a merger sequence appears impossible.

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 consider spherical stellar clusters with a broad mass function and a relaxation time short enough so that the segregation of massive stars toward the centre occurs before they have time to evolve off the main sequence. The relaxational and collisional dynamics of model clusters is followed with a high-resolution Monte Carlo code. Stellar collisions are treated in a realistic way, through the use of a outcome of a very large set of SPH simulations. We find that, even in proto-galactic nuclei models with high velocity dispersions, run-away growth of a very massive star (VMS, M>100 M_sun) occurs in all cases when the core collapse time is shorter than the MS life time of massive stars, i.e. 3 Myrs. The VMS is a likely progenitor for an intermediate-mass or massive black hole (IMBH/MBH).

We present theoretical models of black hole (BH) populations in young stellar environments, such as starbursts and young star clusters. Using a population synthesis approach we compute the formation rates and characteristic properties of single and binary BHs for various representative ages and choices of parameters. We find that most of the BHs (typically 80% for an initial 50% binary fraction) are single, but with many originating from primordial binaries (which either merged into a single massive star or were disrupted following a supernova explosion). A smaller but significant fraction (typically 20%) of the BHs remain in binary systems. Main-sequence stars are the most frequent BH companions, but massive BH--BH binaries are the next most numerous group. The most massive BHs found in our simulations reach 80 Msun, and are formed through mergers of massive binary components. If formed in a dense star cluster such a massive stellar BH may become the seed for growth to a more massive (``intermediate-mass'') BH. Although we do not include dynamical interactions, our results provide realistic initial conditions for N-body simulations of dense star clusters (e.g., globular clusters) including primordial BHs.

We present the first results from our new general relativistic, Lagrangian hydrodynamics code, which treats gravity in the conformally flat (CF) limit. The evolution of fluid configurations is described using smoothed particle hydrodynamics (SPH), and the elliptic field equations of the CF formalism are solved using spectral methodsin spherical coordinates. The code was tested on models for which the CF limit is exact, finding good agreement with the classical Oppenheimer-Volkov solution for a relativistic static spherical star as well as the exact semi-analytic solution for a collapsing spherical dust cloud. By computing the evolution of quasi-equilibrium neutron star binary configurations in the absence of gravitational radiation backreaction, we have confirmed that these configurations can remain dynamically stable all the way to the development of a cusp. With an approximate treatment of radiation reaction, we have calculated the complete merger of an irrotational binary configuration from the innermost point on an equilibrium sequence through merger and remnant formation and ringdown, finding good agreement withprevious relativistic calculations. In particular, we find that mass loss is highly suppressed by relativistic effects, but that, for a reasonably stiff neutron star equation of state, the remnant is initially stable against gravitational collapse because of its strong differential rotation. The gravity wave signal derived from our numerical calculation has an energy spectrum which matches extremely well with estimates based solely on quasi-equilibrium results, deviating from the Newtonian power-law form at frequencies below 1 kHz, i.e., within the reach of advanced interferometric detectors.

We study the early dynamical evolution of young, dense star clusters using Monte Carlo simulations for systems with up to N~10^7 stars. Rapid mass segregation of massive main-sequence stars and the development of the Spitzer instability can drive these systems to core collapse in a small fraction of the initial half-mass relaxation time. If the core collapse time is less than the lifetime of the massive stars, all stars in the collapsing core may then undergo a runaway collision process leading to the formation of a massive black hole. Here we study in detail the first step in this process, up to the occurrence of core collapse. We have performed about 100 simulations for clusters with a wide variety of initial conditions, varying systematically the cluster density profile, stellar IMF, and number of stars. We also considered the effects of initial mass segregation and stellar evolution mass loss. Our results show that, for clusters with a moderate initial central concentration and any realistic IMF, the ratio of core collapse time to initial half-mass relaxation time is typically ~0.1, in agreement with the value previously found by direct N-body simulations for much smaller systems. Models with even higher central concentration initially, or with initial mass segregation (from star formation) have even shorter core-collapse times. Remarkably, we find that, for all realistic initial conditions, the mass of the collapsing core is always close to ~10^-3 of the total cluster mass, very similar to the observed correlation between central black hole mass and total cluster mass in a variety of environments. We discuss the implications of our results for the formation of intermediate-mass black holes in globular clusters and super star clusters, ultraluminous X-ray sources, and seed black holes in proto-galactic nuclei.

We have designed and tested a new relativistic Lagrangian hydrodynamics code, which treats gravity in the conformally flat approximation to general relativity. We have tested the resulting code extensively, finding that it performs well for calculations of equilibrium single-star models, collapsing relativistic dust clouds, and quasi-circular orbits of equilibrium solutions. By adding in a radiation reaction treatment, we compute the full evolution of a coalescing binary neutron star system. We find that the amount of mass ejected from the system, much less than a percent, is greatly reduced by the inclusion of relativistic gravitation. The gravity wave energy spectrum shows a clear divergence away from the Newtonian point-mass form, consistent with the form derived from relativistic quasi-equilibrium fluid sequences.

In dense stellar systems the frequent dynamical interactions between stars play a crucial role in the formation and evolution of compact binaries. We study these processes using a novel approach combining a state-of-the-art binary population synthesis code with a simple treatment of dynamical interactions in dense star cluster cores. Here we focus on the dynamical and evolutionary processes leading to the formation of compact binaries containing white dwarfs in dense globular clusters. We demonstrate that dynamics can increase by factors ~ 2 - 100 the production rates of interesting binaries such as cataclysmic variables, ``nonflickerers'' (He white dwarfs with a heavier dark companion), merging white dwarf binaries with total masses above the Chandrasekhar limit, and white dwarf binaries emitting gravitational waves in the LISA band.

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.

We consider the implications of the presence of ~1 stellar-mass black hole (BH) at the center of a dense globular cluster. We show that BH X-ray binaries formed through exchange interactions are likely to have extremely low duty cycles (<0.001), consistent with the absence of BH X-ray transients in Galactic globular clusters. In contrast, we find that BH X-ray binaries formed through tidal capture would be persistent, bright X-ray sources. Given the absence of any such source and the very high interaction rates, we conclude that tidal capture of a main-sequence star by a BH most likely leads to the complete disruption of the star.

We have calculated the evolution of 60 model binary systems consisting of helium stars in the mass range of M_He= 2.5-6Msun with a 1.4Msun neutron star companion to investigate the formation of double neutron star systems.Orbital periods ranging from 0.09 to 2 days are considered, corresponding to Roche lobe overflow starting from the helium main sequence to after the ignition of carbon burning in the core. We have also examined the evolution into a common envelope phase via secular instability, delayed dynamical instability, and the consequence of matter filling the neutron star's Roche lobe. The survival of some close He-star neutron-star binaries through the last mass transfer episode (either dynamically stable or unstable mass transfer phase) leads to the formation of extremely short-period double neutron star systems (with P<~0.1days). In addition, we find that systems throughout the entire calculated mass range can evolve into a common envelope phase, depending on the orbital period at the onset of mass transfer. The critical orbital period below which common envelope evolution occurs generally increases with M_He. In addition, a common envelope phase may occur during a short time for systems characterized by orbital periods of 0.1-0.5 days at low He-star masses (~ 2.6-3.3Msun). The existence of a short-period population of double neutron stars increases the predicted detection rate of inspiral events by ground-based gravitational-wave detectors and impacts their merger location in host galaxies and their possible role as gamma-ray burst progenitors. We use a set of population synthesis calculations and investigate the implications of the mass-transfer results for the orbital properties of DNS populations.

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 review possible dynamical formation processes for central massive black holes in dense star clusters. We focus on the early dynamical evolution of young clusters containing a few thousand to a few million stars. One natural formation path for a central seed black hole in these systems involves the development of the Spitzer instability, through which the most massive stars can drive the cluster to core collapse in a very short time. The sudden increase in the core density then leads to a runaway collision process and the formation of a very massive merger remnant, which must then collapse to a black hole. Alternatively, if the most massive stars end their lives before core collapse, a central cluster of stellar-mass black holes is formed. This cluster will likely evaporate before reaching the highly relativistic state necessary to drive a runaway merger process through gravitational radiation, thereby avoiding the formation of a central massive black hole. We summarize the conditions under which these different paths will be followed, and present the results of recent numerical simulations demonstrating the process of rapid core collapse and runaway collisions between massive stars.

Close to 70 radio pulsars have now been discovered in 23 globular clusters, with a record 22 pulsars observed in 47 Tuc alone. Accurate timing solutions, including positions in the cluster, are known for many of the pulsars. These recent observations provide a unique opportunity to re-examine theoretically the formation and evolution of recycled pulsars in globular clusters. This brief review focuses on dynamical exchange interactions between neutron stars and primordial binaries, through which neutron stars can acquire intermediate-mass binary companions. The later evolution of these intermediate-mass binaries leads naturally to the two main types of binary millisecond pulsars observed in clusters.

Instabilities and strong dynamical interactions between multiple giant planets have been proposed as a possible explanation for the surprising orbital properties of extrasolar planetary systems. In particular, dynamical instabilities seem to provide a natural mechanism for producing the highly eccentric orbits seen in many systems. Previously, we performed numerical integrations for the dynamical evolution of planetary systems containing two giant planets of equal masses initially in nearly circular orbits very close to the dynamical stability limit. We found the ratio of collisions to ejections in these simulations was greater than the ratio of circular orbits to eccentric orbits among the known extrasolar planets. Further, the mean eccentricity of the planets remaining after an ejection was larger than the mean eccentricity of the known extrasolar planets. Recently, we have performed additional integrations, generalizing to consider two planets of unequal masses. Our new simulations reveal that the two-planet scattering model can produce a distribution of eccentricities consistent with the observed eccentricity distribution for plausible mass distributions. Additionally, this model predicts a maximum eccentricity of about 0.8, in agreement with observations. Early results from simulations of three equal-mass planets also reveal a reduced frequency of collisions and a broad range of final eccentricities for the retained inner planet.

Coalescing binary neutron stars (NS) are expected to be an important source of gravitational waves (GW) detectable by laser interferometers. We present here a simple method for determining the compactness ratio M/R of NS based on the observed deviation of the GW energy spectrum from point-mass behavior at the end of an inspiral event. Our method is based on the properties of quasi-equilibrium binary NS sequences and does not require the computation of the full GW signal h(t). Combined with the measurement of the NS masses from the GW signal during inspiral, the determination of M/R will allow very strong constraints to be placed on the equation of state of nuclear matter at high densities.

We examine the theoretical implications of a population of low-mass helium-core white dwarfs in globular clusters. In particular, we focus on the observed population in the core of NGC 6397, where several low-mass white dwarf canditates have been identified as ``non-flickerers'' by Cool and collaborators. Age and mass estimates from cooling models, combined with dynamical and evolutionary considerations, lead us to infer that the dark binary companions are C/O white dwarfs rather than neutron stars. Furthermore, we find that the progenitor binaries very likely underwent an exchange interaction within the last 1.e9 years. We examine the prospects for detecting a similar population in other globular clusters, with particular attention to the case of 47 Tuc.

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.

Using our new post-Newtonian (PN) smoothed particle hydrodynamics (SPH) code, we have studied numerically the mergers of neutron star binaries with irrotational initial configurations. Here we describe a new method for constructing numerically accurate initial conditions for irrotational binary systems with circular orbits in PN gravity. We then compute the 3D hydrodynamic evolution of these systems until the two stars have completely merged, and we determine the corresponding GW signals. We present results for systems with different binary mass ratios, and for neutron stars represented by polytropes with $\Gamma=2$ or $\Gamma=3$. Compared to mergers of corotating binaries, we find that irrotational binary mergers produce similar peak GW luminosities, but they shed almost no mass at all to large distances. The dependence of the GW signal on numerical resolution for calculations performed with N>10^5 SPH particles is extremely weak, and we find excellent agreement between runs utilizing N=10^5 and N=10^6 SPH particles (the largest SPH calculation ever performed to study such irrotational binary mergers). We also compute GW energy spectra based on all calculations reported here and in our previous works. We find that PN effects lead to clearly identifiable features in the GW energy spectrum of binary neutron star mergers, which may yield important information about the nuclear equation of state at extreme densities.

Collisions of main sequence stars occur frequently in dense star clusters. In open and globular clusters, these collisions produce merger remnants that may be observed as blue stragglers. Detailed theoretical models of this process require lengthy hydrodynamic computations in three dimensions. However, a less computationally expensive approach, which we present here, is to approximate the merger process (including shock heating, hydrodynamic mixing, mass ejection, and angular momentum transfer) with simple algorithms based on conservation laws and a basic qualitative understanding of the hydrodynamics. These algorithms have been fine tuned through comparisons with the results of our previous hydrodynamic simulations. We find that the thermodynamic and chemical composition profiles of our simple models agree very well with those from recent SPH (smoothed particle hydrodynamics) calculations of stellar collisions, and the subsequent stellar evolution of our simple models also matches closely that of the more accurate hydrodynamic models. Our algorithms have been implemented in an easy to use software package, which we are making publicly available (see http://vassun.vassar.edu/~lombardi/mmas/). This software could be used in combination with realistic dynamical simulations of star clusters that must take into account stellar collisions.

We investigate the orbits of compact binary systems during the final inspiral period before coalescence by integrating numerically the second-order post-Newtonian equations of motion. We include spin-orbit and spin-spin coupling terms, which, according to a recent study by Levin [J. Levin, Phys. Rev. Lett. 84, 3515 (2000)], may cause the orbits to become chaotic. To examine this claim, we study the divergence of initially nearby phase-space trajectories and attempt to measure the Lyapunov exponent gamma. Even for systems with maximally spinning objects and large spin-orbit misalignment angles, we find no chaotic behavior. For all the systems we consider, we can place a strict lower limit on the divergence time t_L=1/gamma that is many times greater than the typical inspiral time, suggesting that chaos should not adversely affect the detection of inspiral events by upcoming gravitational-wave detectors.

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 report here on two complementary population synthesis studies which relate directly to the formation and evolution of neutron star binaries in globular clusters. In the first, we compute the probability of retaining neutron stars in globular clusters, and quantitatively confirm the idea that the retention fraction for neutron stars born in binary systems is greatly enhanced over those born in isolated stars. However, the retention fraction may well be insufficient to explain the current population of neutron star binaries. In the second study, we follow a large population of primordial binaries and neutron stars throughout the lifetime of a globular cluster whose properties may be similar to 47 Tuc. We directly compute all 3-body interactions among binary systems, neutron stars, and isolated field stars throughout the history of the cluster. The evolution of certain types of neutron star binaries is followed up to the current epoch. The numbers of close, recycled, binary radio pulsars are evaluated and compared with the results of radio observations.

We study the dynamical evolution of globular clusters using our new 2-D Monte Carlo code, and we calculate the lifetimes of clusters in the Galactic environment. We include the effects of a mass spectrum, mass loss in the Galactic tidal field, and stellar evolution. We consider initial King models containing N = 10^5 - 3x10^5 stars, and follow the evolution up to core collapse, or disruption, whichever occurs first. We find that the lifetimes of our models are significantly longer than those obtained using 1-D Fokker-Planck (F-P) methods. We also find that our results are in very good agreement with recent 2-D F-P calculations, for a wide range of initial conditions. Our results show that the direct mass loss due to stellar evolution can significantly accelerate the mass loss through the tidal boundary, causing most clusters with a low initial central concentration (Wo <~ 3) to disrupt quickly in the Galactic tidal field. Only clusters born with high initial central concentrations (Wo >~ 7) or steep initial mass functions are likely to survive to the present and undergo core collapse. We also study the orbital characteristics of escaping stars, and find that the velocity distribution of escaping stars in collapsing clusters looks significantly different from the distribution in disrupting clusters. We calculate the lifetime of a cluster on an eccentric orbit in the Galaxy, such that it fills its Roche lobe only at perigalacticon. We find that such an orbit can extend the lifetime by at most a factor of a few compared to a circular orbit in which the cluster fills its Roche lobe at all times.

The final burst of gravitational radiation emitted by coalescing binary neutron stars carries direct information about the neutron star fluid, and, in particular, about the equation of state of nuclear matter at extreme densities. The final merger may also be accompanied by a detectable electromagnetic signal, such as a gamma-ray burst. In this paper, we summarize the results of theoretical work done over the past decade that has led to a detailed understanding of this hydrodynamic merger process for two neutron stars, and we discuss the prospects for the detection and physical interpretation of the gravity wave signals by ground-based interferometers such as LIGO. We also present results from our latest post-Newtonian SPH calculations of binary neutron star coalescence, using up to 10^6 SPH particles to compute with higher spatial resolution than ever before the merger of an initially irrotational system. We discuss the detectability of our calculated gravity wave signals based on power spectra.

We review recent 3D hydrodynamic calculations of stellar collisions using the Smoothed Particle Hydrodynamics (SPH) method, and we discuss the implications of the results for the formation and evolution of blue stragglers in globular clusters. We also discuss the construction of simple analytic models for merger remnants, approximating the mass loss, shock heating, mixing, and angular momentum transport during a collision with simple algorithms that can be calibrated using our 3D hydrodynamic results. The thermodynamic and chemical composition profiles predicted by these simple models are compared with those from our SPH simulations, demonstrating that our new models provide accurate representations of true collisional merger remnants.

Instabilities and strong dynamical interactions between several giant planets have been proposed as a possible explanation for the surprising orbital properties of extrasolar planetary systems. In particular, dynamical instabilities would seem to provide a natural mechanism for producing the highly eccentric orbits seen in many systems. Here we present results from numerical integrations for the dynamical evolution of planetary systems containing two identical giant planets in nearly circular orbits very close to the dynamical stability limit. We determine the statistical properties of the three main types of systems resulting from the development of an instability: systems containing one planet, following either a collision between the two initial planets, or the ejection of one of them to infinity, and systems containing two planets in a new, quasi-stable configuration. We discuss the implications of our results for the formation and evolution of observed extrasolar planetary systems. We conclude that the distributions of eccentricities and semimajor axes for observed systems cannot be explained easily by invoking dynamical interactions between two planets initially on circular orbits. While highly eccentric orbits can be produced naturally by these interactions, collisions between the two planets, which occur frequently in the range of observed semimajor axes, would result in many more nearly circular orbits than in the observed sample.

The radio millisecond pulsar PSR B1620-26 is part of an extraordinary triple star system in the globular cluster M4. The inner companion to the neutron star is thought to be a white dwarf of mass m1 ~0.3 Msun in an orbit of period ~0.5 year. The nature and orbital characteristics of the second, more distant companion, have remained a mystery for many years. A theoretical analysis of the latest available radio pulsar timing data is presented here, allowing us to determine approximately the mass and orbital parameters of the second companion. Remarkably, the current best-fit parameters correspond to a second companion of **planetary mass**, with m2 ~0.007 Msun, in an orbit of eccentricity e2 ~0.45 and with a large semimajor axis a2 ~60 AU. The short dynamical lifetime of this very wide triple in M4 suggests that large numbers of such planets must be present in globular clusters. We also address the question of the anomalously high eccentricity of the inner binary pulsar. While this eccentricity could have been induced during the same dynamical interaction that created the triple, we find that it could also naturally arise from long-term secular perturbation effects in the triple, combining the general relativistic precession of the inner orbit with the Newtonian gravitational perturbation by the outer planet.

Using our new Post-Newtonian SPH (smoothed particle hydrodynamics) code, we study the final coalescence and merging of neutron star (NS) binaries. We vary the stiffness of the equation of state (EOS) as well as the initial binary mass ratio and stellar spins. Results are compared to those of Newtonian calculations, with and without the inclusion of the gravitational radiation reaction. We find a much steeper decrease in the gravity wave peak strain and luminosity with decreasing mass ratio than would be predicted by simple point-mass formulae. For NS with softer EOS (which we model as simple $\Gamma=2$ polytropes) we find a stronger gravity wave emission, with a different morphology than for stiffer EOS (modeled as $\Gamma=3$ polytropes as in our previous work). We also calculate the coalescence of NS binaries with an irrotational initial condition, and find that the gravity wave signal is relatively suppressed compared to the synchronized case, but shows a very significant second peak of emission. Mass shedding is also greatly reduced, and occurs via a different mechanism than in the synchronized case. We discuss the implications of our results for gravity wave astronomy with laser interferometers such as LIGO, and for theoretical models of gamma-ray bursts (GRBs) based on NS mergers.

Globular clusters house a population of compact binaries that will be interesting gravitational wave sources for LISA. We provide estimates for the numbers of sources of several categories and discuss the sensitivity of LISA to detecting these sources. The estimated total number of detectable sources ranges from about 10 to about 1000 with gravitational wave frequencies above 1 mHz. These sources are typically undetectable by any other means and thus offer an opportunity for doing true gravitational-wave astronomy. The detection of these sources would provide information about both binary star evolution and the dynamics of globular clusters.

We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitted much higher numerical resolution of the hydrodynamics. We also present evolutionary tracks for the resulting collision products, which emerge as rapidly rotating blue stragglers. The rotating collision products are brighter, bluer and remain on the main sequence longer than their non-rotating counterparts. In addition, they retain their rapid rotation rates throughout their main sequence lifetime. Rotationally-induced mixing strongly affects the evolution of the collision products, although it is not sufficient to mix the entire star. We discuss the implications of these results for studies of blue straggler populations in clusters. This work shows that off-axis collision products cannot become blue stragglers unless they lose a large fraction of their initial angular momentum. The mechanism for this loss is not apparent, although some possibilities are discussed.

We have developed a new parallel supercomputer code based on Henon's Monte Carlo method for simulating the dynamical evolution of globular clusters. This new code allows us to calculate the evolution of a cluster containing a realistic number of stars (N ~ 10^5 - 10^6) in about a day of computing time. The discrete, star-by-star representation of the cluster in the simulation allows us to treat naturally a number of important processes, including single and binary star evolution, all dynamical interactions of single stars and primordial binaries, and tidal interactions with the Galaxy.

We present the first results from our Post-Newtonian (PN) Smoothed Particle Hydrodynamics (SPH) code, which has been used to study the coalescence of binary neutron star (NS) systems. The Lagrangian particle-based code incorporates consistently all lowest-order (1PN) relativistic effects, as well as gravitational radiation reaction, the lowest-order dissipative term in general relativity. We test our code on sequences of single NS models of varying compactness, and we discuss ways to make PN simulations more relevant to realistic NS models. We also present a PN SPH relaxation procedure for constructing equilibrium models of synchronized binaries, and we use these equilibrium models as initial conditions for our dynamical calculations of binary coalescence. Though unphysical, since tidal synchronization is not expected in NS binaries, these initial conditions allow us to compare our PN work with previous Newtonian results. We compare calculations with and without 1PN effects, for NS with stiff equations of state, modeled as polytropes with $\Gamma=3$. We find that 1PN effects can play a major role in the coalescence, accelerating the final inspiral and causing a significant misalignment in the binary just prior to final merging. In addition, the character of the gravitational wave signal is altered dramatically, showing strong modulation of the exponentially decaying waveform near the end of the merger. We also discuss briefly the implications of our results for models of gamma-ray bursts at cosmological distances.

We present the results of Monte Carlo simulations for the dynamical evolution of star clusters containing two stellar populations with individual masses m1 and m2 > m1, and total masses M1 and M2 < M1. We use both King and Plummer model initial conditions and we perform simulations for a wide range of individual and total mass ratios, m2/m1 and M2/M1. We ignore the effects of binaries, stellar evolution, and the galactic tidal field. The simulations use N = 10^5 stars and follow the evolution of the clusters until core collapse. We find that the departure from energy equipartition in the core follows approximately the theoretical predictions of Spitzer (1969) and Lightman & Fall (1978), and we suggest a more exact condition that is based on our results. We find good agreement with previous results obtained by other methods regarding several important features of the evolution, including the pre-collapse distribution of heavier stars, the time scale on which equipartition is approached, and the extent to which core collapse is accelerated by a small subpopulation of heavier stars. We briefly discuss the possible implications of our results for the dynamical evolution of primordial black holes and neutron stars in globular clusters.

Coalescing compact binaries are thought to be involved in a wide variety of astrophysical phenomena. In particular, they are important sources of gravitational radiation for both ground-based and space-based laser-interferometer detectors, and they may be sources of supernova explosions or gamma-ray bursts. Mergers of two white dwarfs may produce neutron stars with peculiar properties, including perhaps millisecond radio pulsars sometimes accompanied by planets (as observed in PSR 1257+12). According to a widely held belief, the coalescence of two neutron stars should produce a rapidly rotating black hole surrounded by an accretion disk or torus, but this is by no means certain. This review paper focuses on the final hydrodynamic coalescence and merger of double neutron stars and double white dwarfs, and addresses the question of the nature of the final merger products.

We present a new dynamical scenario for the formation of short-period binary millisecond pulsars in globular clusters. Our work is motivated by the recent observations of 20 radio pulsars in 47 Tuc. In a dense cluster such as 47 Tuc, most neutron stars acquire binary companions through exchange interactions with primordial binaries. The resulting systems have semimajor axes in the range \~0.1-1 AU and neutron star companion masses ~1-3 Msun. For many of these systems we find that, when the companion evolves off the main sequence and fills its Roche lobe, the subsequent mass transfer is dynamically unstable. This leads to a common envelope phase and the formation of short-period neutron star - white dwarf binaries. For a significant fraction of these binaries, the decay of the orbit due to gravitational radiation will be followed by a period of stable mass transfer driven by a combination of gravitational radiation and tidal heating of the companion. The properties of the resulting short-period binaries match well those of observed binary pulsars in 47 Tuc.

Smoothed Particle Hydrodynamics is a multidimensional Lagrangian method of numerical hydrodynamics that has been used to tackle a wide variety of problems in astrophysics. Here we develop the basic equations of the SPH scheme, and we discuss some of its numerical properties and limitations. As an illustration of typical astrophysical applications, we discuss recent calculations of stellar interactions, including collisions between main sequence stars and the coalescence of compact binaries.

We derive octupole-level secular perturbation equations for hierarchical triple systems, using classical Hamiltonian perturbation techniques. By extending previous work done to leading (quadrupole) order to octupole level (i.e., including terms of order alpha^3, where alpha = a_1/a_2 < 1 is the ratio of semimajor axes) we obtain expressions that are applicable to a much wider range of parameters. For triple systems containing a close inner binary, we also discuss the possible interaction between the classical Newtonian perturbations and the general relativistic precession of the inner orbit. In some cases we show that this interaction can lead to resonances and a significant increase in the maximum amplitude of eccentricity perturbations. We establish the validity of our analytic expressions by providing detailed comparisons with the results of direct numerical integrations of the three-body problem obtained for a large number of representative cases. We also discuss applications of the theory in the context of several observed triple systems of current interest, including the millisecond pulsar PSR B1620-26 in M4, the giant planet in 16 Cygni, and the protostellar binary TMR-1. [short version]

We present a new parallel supercomputer implementation of the Monte-Carlo method for simulating the dynamical evolution of globular star clusters. Our method is based on a modified version of Henon's Monte-Carlo algorithm for solving the Fokker-Planck equation. Our code allows us to follow the evolution of a cluster containing up to 5x10^5 stars to core collapse in < 40 hours of computing time. In this paper we present the results of test calculations for clusters with equal-mass stars, starting from both Plummer and King model initial conditions. We consider isolated as well as tidally truncated clusters. Our results are compared to those obtained from approximate, self-similar analytic solutions, from direct numerical integrations of the Fokker-Planck equation, and from direct N-body integrations performed on a GRAPE-4 special-purpose computer with N=16384. In all cases we find excellent agreement with other methods, establishing our new code as a robust tool for the numerical study of globular cluster dynamics using a realistic number of stars.

Coalescing compact binaries with neutron star (NS) or black hole (BH) components are important sources of gravitational waves for the laser-interferometer detectors currently under construction, and may also be sources of gamma-ray bursts at cosmological distances. This paper focuses on the final hydrodynamic coalescence and merger of NS--NS binaries, and addresses the question of whether black hole formation is the inevitable final fate of these systems.

Particle methods play an important role in the study of a wide variety of astrophysical fluid dynamics problems. The different methods currently in use are all variants of the so-called Smoothed Particle Hydrodynamics (SPH) scheme introduced by Lucy and Gingold & Monaghan more than twenty years ago. This paper presents a complete introduction to SPH in its modern form, and discusses some of the main numerical properties of the scheme. In particular, the convergence properties of SPH are studied, as a function of the number of particles N and the number of interacting neighbors N_n, using a simple analysis based on sound waves. It is shown that consistency of SPH (i.e., convergence towards a physical solution) requires both N -> infinity and N_n -> infinity, with the smoothing length h -> 0, i.e., N_n/N -> 0.

We present a new theoretical analysis of the PSR B1620-26 triple system in the globular cluster M4, based on the latest radio pulsar timing data, which now include measurements of five time derivatives of the pulse frequency. These data allow us to determine the mass and orbital parameters of the second companion completely (up to the usual unknown orbital inclination angle $i_2$). The current best-fit parameters correspond to a second companion of planetary mass, $m_2 \sin i_2 \simeq 7\times10^{-3} M_\odot$, in an orbit of eccentricity $e_2\simeq 0.45$ and semimajor axis $a_2\simeq 60 $AU. Using numerical scattering experiments, we study a possible formation scenario for the triple system, which involves a dynamical exchange interaction between the binary pulsar and a primordial star-planet system. The current orbital parameters of the triple are consistent with such a dynamical origin, and suggest that the separation of the parent star-planet system was very large, $\go 50 $AU. We also examine the possible origin of the anomalously high eccentricity of the inner binary pulsar. While this eccentricity could have been induced during the same dynamical interaction that created the triple, we find that it could equally well arise from long-term secular perturbation effects in the triple, combining the general relativistic precession of the inner orbit with the Newtonian gravitational perturbation of the planet. The detection of a planet in this system may be taken as evidence that large numbers of extrasolar planetary systems, not unlike those discovered recently in the solar neighborhood, also exist in old star clusters.

Twenty millisecond radio pulsars have now been observed in the globular cluster 47 Tuc. This is by far the largest sample of radio pulsars known in any globular cluster. These recent observations provide a unique opportunity to re-examine theoretically the formation and evolution of recycled pulsars in globular clusters.

Coalescing compact binaries with neutron star or black hole components provide the most promising sources of gravitational radiation for detection by the LIGO/VIRGO/GEO/TAMA laser interferometers now under construction. This fact has motivated several different theoretical studies of the inspiral and hydrodynamic merging of compact binaries. Analytic analyses of the inspiral waveforms have been performed in the Post-Newtonian approximation. Analytic and numerical treatments of the coalescence waveforms from binary neutron stars have been performed using Newtonian hydrodynamics and the quadrupole radiation approximation. Numerical simulations of coalescing black hole and neutron star binaries are also underway in full general relativity. Recent results from each of these approaches will be described and their virtues and limitations summarized.

Using the Yale stellar evolution code, we have calculated theoretical models for nearby stars with planetary-mass companions in short-period nearly circular orbits: 51 Pegasi, Tau Bootis, Upsilon Andromedae, Rho Cancri, and Rho Coronae Borealis. We present tables listing key stellar parameters such as mass, radius, age, and size of the convective envelope as a function of the observable parameters (luminosity, effective temperature, and metallicity), as well as the unknown helium fraction. For each star we construct best models based on recently published spectroscopic data and the present understanding of galactic chemical evolution. We discuss our results in the context of planet formation theory, and, in particular, tidal dissipation effects and stellar metallicity enhancements.

We have performed a series of systematic tests to evaluate the effects of spurious transport in three-dimensional smoothed particle hydrodynamics (SPH) calculations. Our tests investigate (i) particle diffusion, (ii) shock heating, (iii) numerical viscosity, and (iv) angular momentum transport. The results are useful for quantifying the accuracy of the SPH scheme, especially for problems where shear flows or shocks are present, as well as for problems where true hydrodynamic mixing is relevant. We examine the different forms of artificial viscosity (AV) which have been proposed by Monaghan, by Hernquist & Katz, and by Balsara. For each form, our tests suggest a single set of values for the AV parameters $\alpha$ and $\beta$ (coefficients of the linear and quadratic terms) which are appropriate in a large number of situations. We also discuss how these parameters should be adjusted depending on the goals of the particular application. We find that both the Hernquist & Katz and Balsara forms introduce relatively small amounts of numerical viscosity. Furthermore, both Monaghan's and Balsara's AV do well at treating shocks and at limiting the amount of spurious mixing. For these reasons, we endorse the Balsara AV for use in a broad range of applications.

Using an energy variational method, we calculate quasi-equilibrium configurations of binary neutron stars modeled as compressible triaxial ellipsoids obeying a polytropic equation of state. Our energy functional includes terms both for the internal hydrodynamics of the stars and for the external orbital motion. We add the leading post-Newtonian (PN) corrections to the internal and gravitational energies of the stars, and adopt hybrid orbital terms which are fully relativistic in the test-mass limit and always accurate to PN order. The total energy functional is varied to find quasi-equilibrium sequences for both corotating and irrotational binaries in circular orbits. We examine how the orbital frequency at the innermost stable circular orbit depends on the polytropic index n and the compactness parameter GM/Rc^2. We find that, for a given GM/Rc^2, the innermost stable circular orbit along an irrotational sequence is about 17% larger than the innermost secularly stable circular orbit along the corotating sequence when n=0.5, and 20% larger when n=1. We also examine the dependence of the maximum neutron star mass on the orbital frequency and find that, if PN tidal effects can be neglected, the maximum equilibrium mass increases as the orbital separation decreases.

We explore the evolution of collisionally merged stars in the blue straggler region of the HR diagram. The starting models for our stellar evolution calculations are the results of the smoothed particle hydrodynamics (SPH) simulations of parabolic collisions between main sequence stars performed by Lombardi, Rasio and Shapiro (1996). Since SPH and stellar evolution codes employ different and often contradictory approximations, it is necessary to treat the evolution of these products carefully. The mixture and disparity of the relevant timescales (hydrodynamic, thermal relaxation and nuclear burning) and of the important physical assumptions between the codes makes the combined analysis of the problem challenging, especially during the initial thermal relaxation of the star. In particular, the treatment of convection is important, and semiconvection must be modeled in some detail. The products of seven head-on collisions are evolved through their initial thermal relaxation, and then through the main sequence phase to the base of the giant branch. Their evolutionary tracks are presented. In contrast to the assumptions in previous work, these collision products do not develop substantial convective regions during their thermal relaxation, and therefore are not mixed significantly after the collision.

We present a general method for determining the masses and orbital parameters of binary millisecond pulsars with long orbital periods (P_orb >> 1 yr), using timing data in the form of pulse frequency derivatives. We apply our method to analyze the properties of the second companion in the PSR B1620-26 triple system. We use the latest timing data for this system to constrain the mass and orbital parameters of the second companion. We find that all possible solutions have a mass m_2 in the range 2.4 10^-4 M_sun <= m_2 sin i_2 <= 1.2 10^-2 M_sun, i.e., almost certainly excluding a second companion of stellar mass and suggesting instead that the system contains a planet or a brown dwarf. Using Monte-Carlo realizations of the triple configuration in three dimensions we find the most probable value of m_2 to be 0.010(5) M_sun, corresponding to a distance of 38(6) AU from the center of mass of the inner binary (the errors indicate 80% confidence intervals). We also apply our method to analyze the planetary system around PSR B1257+12, where a distant, giant planet may be present in addition to the three well-established Earth-mass planets. We find that the simplest interpretation of the frequency derivatives implies the presence of a fourth planet with a mass of ~100 M_earth in a circular orbit of radius ~40 AU.

Previous timing data for PSR~B1620$-$26 were consistent with a second companion mass $m_2$ anywhere in the range $\sim10^{-3}-1\,M_\odot$, i.e., from a Jupiter-type planet to a star. We present the latest timing parameters for the system, including a significant change in the projected semi-major axis of the inner binary, a marginal detection of the fourth time derivative of the pulse frequency, and the pulsar proper motion (which is in agreement with published values for the proper motion of M4), and use them to further constrain the mass $m_2$ and the orbital parameters. Using the observed value of $\stackrel{\ldots.}{f}$, we obtain a one-parameter family of solutions, all with $m_2 \lo 10^{-2}\,M_\odot$, i.e., excluding stellar masses. Varying $\stackrel{\ldots.}{f}$ within its formal $1\sigma$ error bar does not affect the mass range significantly. However, if we vary $\stackrel{\ldots.}{f}$ within a $4\sigma$ error bar, we find that stellar-mass solutions are still possible. We also calculate the predicted rate of change of the projected semi-major axis of the inner binary and show that it agrees with the measured value.

We derive analytical expressions for the change in the orbital eccentricity of a binary following a distant encounter with a third star on a hyperbolic or parabolic orbit. To establish the accuracy of these expressions, we present detailed comparisons with the results of direct numerical integrations of the equations of motion for the three bodies. We treat with particular care the difficult case of a binary with zero initial eccentricity. In this case, we show that the eccentricity $\delta e$ induced by the encounter declines in general as a power-law, $\delta e\propto (a/\rp)^{5/2}$, where $a$ is the binary semi-major axis and $\rp$ is the periastron distance of the encounter. This power-law arises from the octupole-level secular perturbation of the binary. In contrast, non-secular quadrupole-level perturbations induce an eccentricity change that declines exponentially with $\rp$. These non-secular effects can become dominant at sufficiently small $\rp$, for a sufficiently high relative velocity, or for a sufficiently massive perturber. We also derive cross sections for eccentricity change and compare our results with those of previous studies based on numerical scattering experiments. Our results have important implications for a number of astrophysical problems including, in particular, the evolution of binary millisecond pulsars in globular clusters.

The 4.2-day orbit of the newly discovered planet around 51~Pegasi is formally unstable to tidal dissipation. However, the orbital decay time in this system is longer than the main-sequence lifetime of the central star. Given our best current understanding of tidal interactions, a planet of Jupiter's mass around a solar-like star could have dynamically survived in an orbit with a period as short as $\sim10\,$hr. Since radial velocities increase with decreasing period, we would expect to find those planets close to the tidal limit first and, unless this is a very unusual system, we would expect to find many more. We also consider the tidal stability of planets around more evolved stars and we re-examine in particular the question of whether the Earth can dynamically survive the red-giant phase in the evolution of the Sun.

Equilibrium fluid configurations for close binary systems can become {\em globally unstable\/}. Instabilities arise from the strong tidal interaction between the two components, which tends to make the effective two-body potential governing the orbital motion steeper than $1/r$. As a result, a circular orbit can become unstable to small radial perturbations. The instability can be either secular or dynamical. In both cases it leads to the coalescence and merging of the two components on a timescale generally much shorter than the lifetime of a stable system. Binary coalescence is relevant to a wide variety of astrophysical problems including the formation of blue stragglers, Type~Ia supernovae, the evolution of massive X-ray binaries, and the detection of gravitational waves.

The hydrodynamics of collisions and mergers of main-sequence stars is discussed in the light of recent 3-D calculations using the smoothed particle hydrodynamics (SPH) method. Theoretical models for the formation of blue stragglers are reviewed in the context of recent comparisons between the observed properties of blue stragglers in dense globular clusters and the predictions of those models.

We report the results of new SPH calculations of parabolic collisions between two main-sequence stars in a globular cluster. Such collisions are directly relevant to the formation of blue stragglers. In particular, we consider parent stars of mass $M/M_{TO}=0.2,0.5,0.75$, and $1$, where $M_{TO}$ is the cluster turnoff mass (typically about $0.8\,M_\odot$). Low-mass stars (with $M=0.2 M_{TO}$ or $0.5 M_{TO}$) are assumed to be fully convective and are therefore modeled as $n=1.5$ polytropes. Stars at the turnoff (with $M= M_{TO}$) are assumed to be mostly radiative and are modeled as $n=3$ polytropes. Intermediate-mass stars (with $M=0.75 M_{TO}$) are modeled as composite polytropes consisting of a radiative core with polytropic index $n=3$ and a convective envelope with $n=1.5$. We focus our study on the question of hydrodynamic mixing of helium and hydrogen, which plays a crucial role in determining the observable characteristics of blue stragglers. In all cases we find that there is negligible hydrodynamic mixing of helium into the outer envelope of the merger remnant. The amount of hydrogen mixed into the core of the merger depends strongly on the entropy profiles of the two colliding stars. In contrast to what has been done in previous studies, blue stragglers formed by direct stellar collisions should not necessarily be assumed to have initially homogeneous composition profiles.

We discuss the formation of a common envelope system following dynamically unstable mass transfer in a close binary, and the subsequent dynamical evolution and final fate of the envelope. We base our discussion on new three-dimensional SPH calculations that we have performed for a close binary system containing a $4\,M_\odot$ red giant with a $0.7\,M_\odot$ main-sequence star companion. The initial parameters are chosen to model the formation of a system resembling V~471~Tau, a typical progenitor of a cataclysmic variable binary. In our highest-resolution calculation, we find evidence for a corotating region of gas around the central binary. This is in agreement with the theoretical model proposed by Meyer \& Meyer-Hofmeister (1979) for the evolution of common envelope systems, in which this central corotating region is coupled to the envelope through viscous angular momentum transport only. We also find evidence that the envelope is convectively unstable, in which case the viscous dissipation time could be as short as $\sim100$ dynamical times, leading to rapid ejection of the envelope. For V~471~Tau, our results, and the observed parameters of the system, are entirely consistent with rapid envelope ejection on a timescale $\sim1\,$yr and an efficiency parameter $\alpha_{CE}\simeq1$.

The hydrodynamics of stellar collisions and mergers is discussed in the context of blue straggler formation. Emphasis is placed on the very important question of hydrodynamic mixing during the merger process. Recent results of three-dimensional hydrodynamic calculations suggest that the merger remnants produced by stellar collisions are typically {\em not\/} well-mixed. However, comparisons between the observed colors and numbers of blue stragglers in dense clusters and the predictions of theoretical calculations for their stellar evolution appear to require that the initial blue-straggler models be close to chemically homogeneous. The resolution of this apparent conflict is likely to involve the development of convection or other thermal instabilitites that can provide efficient mixing during the contraction of the merger remnant to a thermal equilibrium state on the main sequence.

Low-mass binary millisecond pulsars (LMBPs) are born with very small orbital eccentricities, typically of order $e_i\sim10^{-6}$--$10^{-3}$. In globular clusters, however, higher eccentricities $e_f\gg e_i$ can be induced by dynamical interactions with passing stars. Here we show that the cross section for this process is much larger than previously estimated. This is because, even for initially circular binaries, the induced eccentricity $e_f$ for an encounter with pericenter separation $r_{\rm p}$ beyond a few times the binary semi-major axis $a$ declines only as a power-law, $e_f\propto(r_{\rm p}/a)^ {-5/2}$, and {\it not\/} as an exponential. We find that {\it all\/} currently known \lmbps in clusters were probably affected by interactions, with their current eccentricities typically an order of magnitude or more greater than at birth.

We report the results of new SPH calculations of parabolic collisions between main-sequence (MS) stars. The stars are assumed to be close to the MS turn-off point in a globular cluster and are therefore modeled as $n=3$, $\Gamma=5/3$ polytropes. We find that the high degree of central mass concentration in these stars has a profound effect on the hydrodynamics. In particular, very little hydrodynamic mixing occurs between the dense, helium-rich inner cores and the outer envelopes. As a result, and in contrast to what has been assumed in previous studies, blue stragglers formed by direct stellar collisions are not necessarily expected to have anomalously high helium abundances in their envelopes or to have their cores replenished with fresh hydrogen fuel.

The minimum mass ratio for tidal stability of a contact binary containing two unevolved main-sequence stars is calculated to be $\qmin\simeq0.09$ in the case of a mostly radiative primary, and is higher if an appreciable fraction of the mass lies in a convective envelope. At least one observed system, AW UMa, has a mass ratio just below this value ($q=0.075$), implying, if the system is stable, that the primary must be slightly evolved and must have a very shallow convective envelope. Contact binaries with mass ratios significantly below that of AW UMa should not be observed, since they are tidally unstable and quickly merge into a single, rapidly rotating object, on a timescale $\sim10^3$--$10^4\,$yr.

In addition to their possible relevance to gamma-ray bursts, coalescing binary neutron stars have long been recognized as important sources of gravitational radiation that should become detectable with the new generation of laser interferometers such as LIGO. Hydrodynamics plays an essential role near the end of the coalescence when the two stars finally merge together into a single object. The shape of the corresponding burst of gravitational waves provides a direct probe into the interior structure of a neutron star and the nuclear equation of state. The interpretation of the gravitational waveform data will require detailed theoretical models of the complicated three-dimensional hydrodynamic processes involved. Here we review the results of our recent work on this problem, using both approximate quasi-analytic methods and large-scale numerical hydrodynamics calculations on supercomputers. We also discuss briefly the coalescence of white-dwarf binaries, which are also associated with a variety of interesting astrophysical phenomena.

The hierarchical triple system containing the millisecond pulsar PSR B1620-26 in M4 is the first triple star system ever detected in a globular cluster. Such systems should form in globular clusters as a result of dynamical interactions between binaries. We propose that the triple system containing PSR B1620-26 formed through an exchange interaction between a wide primordial binary and a {\it pre-existing\/} binary millisecond pulsar. This scenario would have the advantage of reconciling the $\sim10^9\,$yr timing age of the pulsar with the much shorter lifetime of the triple system in the core of M4.

We have searched deep optical images of the globular cluster M4 for a possible optical counterpart for the third body in the PSR B1620-26 system. We identify as a possible optical counterpart a star located within $0.3''$ of the nominal pulsar position with $V\approx 20$. This magnitude is consistent with a main sequence cluster member with $M\approx 0.45M_\odot$. However, this star may be blended with fainter stars, and the probability of a chance superposition is non-negligible.

We present a new numerical study of the equilibrium and stability properties of close binary systems using the smoothed-particle hydrodynamics (SPH) technique. We adopt a simple polytropic equation of state $p=K\rho^\gam$ with $\gam=5/3$ and $K=\,$constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio $q\lo0.4$ we find that the (synchronized) Roche limit configuration is secularly unstable. Degenerate binary configurations remain hydrodynamically stable all the way to the Roche limit for all mass ratios $q\ne1$. Dynamically unstable mass transfer can also lead to the rapid coalescence of a binary system, but the details of the hydrodynamic evolution are quite different. We discuss the implications of our results for double white-dwarf and W Ursae Majoris systems.

A simple analytic calculation is presented for the secular evolution of the eccentricities in a hierarchical triple system such as the one containing the millisecond pulsar PSR 1620-26 in M4. If the second companion of PSR 1620-26 is of stellar mass ($m_2\ga0.1\,M_\odot$), an eccentricity as large as that observed today for the inner binary ($e_1\simeq0.03$) could very well have been induced by secular perturbations. In contrast, such a large eccentricity cannot be induced by a second companion of planetary mass.

Recent observations by the ROSAT X-ray satellite of the Vela supernova remnant have revealed, in addition to the previously identified compact nebula, a nearly circular emitting region with a radius of about 4 degrees. The Vela pulsar is slightly off the center of this circular region, consistent with its measured proper motion of about 100km/s and an age of about 10^4 yr. The emitting region is bounded by the main supernova shock. Just outside the shock, the X-ray image reveals several well-defined V-shaped features extending radially outwards. These features are most likely wakes produced by objects moving supersonically through the outside medium. The shapes and orientations of the wakes suggest that these objects have been ejected from the center of the supernova explosion. Their present positions indicate that they have been moving with a mean velocity of a few thousand km/s. We show that pre-existing objects, such as planets in orbit around the progenitor star, could not have been accelerated to sufficiently high velocities or would have been destroyed. Instead, we propose that the observed objects are fragments ejected during the formation of the neutron star. Fragmentation during gravitational collapse is a natural consequence of both convective and rotational instabilities.

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

The unusually large eccentricity ($e_1=0.025$) of the low-mass binary millisecond pulsar PSR B1620-26 can be explained naturally as arising from the secular perturbation of a second, more distant companion. Such a triple configuration has been proposed recently as the most likely cause of the anomalous second period derivative of the pulsar. The current timing data are consistent with a second companion mass $m_2$ as low as $\sim10^{-3}\,M_\odot$, i.e., comparable to that of Jupiter. However, {\em if\/} the eccentricity is indeed produced by secular perturbations, then the second companion must be another star, most likely of mass $m_2\lo1M_\odot$ and in a very eccentric ($e_2\go0.5$) orbit of period $P_2\sim10^2$--$10^3\,$yr. A second companion of planetary mass cannot induce the observed eccentricity. Independent of the mass of the second companion, small changes in the binary pulsar's orbit should become detectable with just a few additional years of timing data. This detection would provide direct confirmation of the triple nature of the system, and an accurate measurement of the effects would place important new constraints on the orbital parameters.

We have performed a series of three-dimensional hydrodynamic calculations of binary coalescence using the smoothed particle hydrodynamics (SPH) method. The initial conditions are exact polytropic equilibrium configurations with $\gam > 5/3$, on the verge of dynamical instability. We calculate the emission of gravitational radiation in the quadrupole approximation. The fully nonlinear development of the instability is followed until a new equilibrium configuration is reached. We find that the properties of this final configuration depend sensitively on both the compressibility and mass ratio. An {\em axisymmetric} merged configuration is always produced when $\gam\lo2.3$. As a consequence, the emission of gravitational radiation shuts off abruptly right after the onset of dynamical instability. In contrast, {\em triaxial\/} merged configurations are obtained when $\gam\go2.3$, and the system continues to emit gravitational waves after the final coalescence. Systems with mass ratios $q\ne1$ typically become dynamically unstable before the onset of mass transfer. Stable mass transfer from one neutron star to another in a close binary is therefore probably ruled out. The maximum amplitude $h_{max}$ and peak luminosity $L_{max}$ of the gravitational waves emitted during the final coalescence are nearly independent of $\gam$, but depend very sensitively on the mass ratio $q$.

The formation of quasar black holes during the hydrodynamic collapse of protogalactic gas clouds is discussed. The dissipational collapse and long- term dynamical evolution of these systems is analysed using three-dimensional numerical simulations. The calculations focus on the final collapse stages of the inner baryonic component and therefore ignore the presence of dark matter. Two types of initial conditions are considered: uniformly rotating spherical clouds, and irrotational ellipsoidal clouds. In both cases the clouds are initially cold, homogeneous, and not far from rotational support ($T/|W| \approx 0.1$). Although the details of the dynamical evolution depend sensitively on the initial conditions, the qualitative features of the final configurations do not. Most of the gas is found to fragment into small dense clumps, that eventually make up a spheroidal component resembling a galactic bulge. About 5\% of the initial mass remains in the form of a smooth disk of gas supported by rotation in the gravitational potential well of the outer spheroid. If a central seed black hole of mass $\gsim10^6M_\odot$ forms, it can grow by steady accretion from the disk and reach a typical quasar black hole mass $\sim10^8M_\odot$ in less than $5\times 10^8\,$yr. In the absence of a sufficiently massive seed, dynamical instabilities in a strongly self- gravitating inner region of the disk will inhibit steady accretion of gas and may prevent the immediate formation of a quasar.

We present a new analytic study of the equilibrium and stability properties of close binary systems containing polytropic components. Our method is based on the use of ellipsoidal trial functions in an energy variational principle. We consider both synchronized and nonsynchronized systems, constructing the compressible generalizations of the classical Darwin and Darwin-Riemann configurations. Our method can be applied to a wide variety of binary models where the stellar masses, radii, spins, entropies, and polytropic indices are all allowed to vary over wide ranges and independently for each component. We find that both secular and dynamical instabilities can develop before a Roche limit or contact is reached along a sequence of models with decreasing binary separation. High incompressibility always makes a given binary system more susceptible to these instabilities, but the dependence on the mass ratio is more complicated. As simple applications, we construct models of double degenerate systems and of low-mass main-sequence-star binaries. We also discuss the orbital evolution of close binary systems under the combined influence of fluid viscosity and secular angular momentum losses from processes like gravitational radiation. We show that the existence of global fluid instabilities can have a profound effect on the terminal evolution of coalescing binaries. The validity of our analytic solutions is examined by means of detailed comparisons with the results of recent numerical fluid calculations in three dimensions.

We study the importance of hydrodynamic effects on the evolution of coalescing binary neutron stars. Using an approximate energy functional constructed from equilibrium solutions for polytropic binary configurations, we incorporate hydrodynamic effects into the calculation of the orbital decay driven by gravitational wave emission. In particular, we follow the transition between the quasi-static, secular decay of the orbit at large separation and the rapid dynamical evolution of configurations approaching contact. We show that a purely Newtonian hydrodynamic instability can significantly accelerate the coalescence at small separation. Such an instability occurs in all close binary configurations containing sufficiently incompressible stars. Calculations are performed for various neutron star masses, radii, spins, and effective polytropic indices. Typically, we find that the radial infall velocity just prior to contact is about 10\% of the tangential orbital velocity. Post-Newtonian effects can move the stability limit to a larger binary separation, and may induce an even larger radial velocity. We also consider the possibility of mass transfer from one neutron star to the other. We show that stable mass transfer is impossible except when the mass of one of the components is very small (less than about 0.4 solar mass) and the viscosity is high enough to maintain corotation.