We find that clouds of optically-thin, pressure-confined gas are prone to fragmentation as they cool below $\sim10^6$ K. This fragmentation follows the lengthscale $\sim{c}_{\text{s}}\,t_{\text{cool}}$, ultimately reaching very small scales ($\sim{0.1} \text{pc}/n$) as they reach the temperature $\sim10^4$ K at which hydrogen recombines. While this lengthscale depends on the ambient pressure confining the clouds, we find that the column density through an individual fragment $N_{\text{cloudlet}}\sim10^{17} \text{cm}^{-3}$ is essentially independent of environment; this column density represents a characteristic scale for atomic gas at $10^4$ K. We therefore suggest that "clouds" of cold, atomic gas may in fact have the structure of a mist or a fog, composed of tiny fragments dispersed throughout the ambient medium. We show that this scale emerges in hydrodynamic simulations, and that the corresponding increase in the surface area may imply rapid entrainment of cold gas. We also apply it to a number of observational puzzles, including the large covering fraction of diffuse gas in galaxy halos, the broad line widths seen in quasar and AGN spectra, and the entrainment of cold gas in galactic winds. While our simulations make a number of assumptions and thus have associated uncertainties, we show that this characteristic scale is consistent with a number of observations, across a wide range of astrophysical environments. We discuss future steps for testing, improving, and extending our model.

Heartbeat stars (HB stars) are a class of eccentric binary stars with close periastron passages. The characteristic photometric HB signal evident in their light curves is produced by a combination of tidal distortion, heating, and Doppler boosting near orbital periastron. Many HB stars continue to oscillate after periastron and along the entire orbit, indicative of the tidal excitation of oscillation modes within one or both stars. These systems are among the most eccentric binaries known, and they constitute astrophysical laboratories for the study of tidal effects. We have undertaken a radial velocity (RV) monitoring campaign of Kepler HB stars in order to measure their orbits. We present our first results here, including a sample of 21 Kepler HB systems, where for 19 of them we obtained the Keplerian orbit and for 3 other systems we did not detect a statistically significant RV variability. Results presented here are based on 218 spectra obtained with the Keck/HIRES spectrograph during the 2015 Kepler observing season, and they have allowed us to obtain the largest sample of HB stars with orbits measured using a single instrument, which roughly doubles the number of HB stars with an RV measured orbit. The 19 systems measured here have orbital periods from 7 to 90 d and eccentricities from 0.2 to 0.9. We show that HB stars draw the upper envelope of the eccentricity - period distribution. Therefore, HB stars likely represent a population of stars currently undergoing high eccentricity migration via tidal orbital circularization, and they will allow for new tests of high eccentricity migration theories.

The dynamical formation of stellar-mass black hole-black hole binaries has long been a promising source of gravitational waves for the Laser Interferometer Gravitational-Wave Observatory (LIGO). Mass segregation, gravitational focusing, and multibody dynamical interactions naturally increase the interaction rate between the most massive black holes in dense stellar systems, eventually leading them to merge. We find that dynamical interactions, particularly three-body binary formation, enhance the merger rate of black hole binaries with total mass M_tot roughly as ~M_tot^beta, with beta >~ 4. We find that this relation holds mostly independently of the initial mass function, but the exact value depends on the degree of mass segregation. The detection rate of such massive black hole binaries is only further enhanced by LIGO's greater sensitivity to massive black hole binaries with M_tot <~ 80 solar masses. We find that for power-law BH mass functions dN/dM ~ M^-alpha with alpha <~ 2, LIGO is most likely to detect black hole binaries with a mass twice that of the maximum initial black hole mass and a mass ratio near one. Repeated mergers of black holes inside the cluster result in about ~5% of mergers being observed between two and three times the maximum initial black hole mass. Using these relations, one may be able to invert the observed distribution to the initial mass function with multiple detections of merging black hole binaries.

We measure the proper motion of the pulsar PSR J1745-2900 relative to the Galactic Center massive black hole, Sgr A*, using the Very Long Baseline Array (VLBA). The pulsar has a transverse velocity of 236 +/- 11 km s^-1 at position angle 22 +/- 2 deg East of North at a projected separation of 0.097 pc from Sgr A*. Given the unknown radial velocity, this transverse velocity measurement does not conclusively prove that the pulsar is bound to Sgr A*; however, the probability of chance alignment is very small. We do show that the velocity and position is consistent with a bound orbit originating in the clockwise disk of massive stars orbiting Sgr A* and a natal velocity kick of <~ 500 km s^-1. An origin among the isotropic stellar cluster is possible but less probable. If the pulsar remains radio-bright, multi-year astrometry of PSR J1745-2900 can detect its acceleration and determine the full three-dimensional orbit. We also demonstrate that PSR J1745-2900 exhibits the same angular broadening as Sgr A* over a wavelength range of 3.6 cm to 0.7 cm, further confirming that the two sources share the same interstellar scattering properties. Finally, we place the first limits on the presence of a wavelength-dependent shift in the position of Sgr A*, i.e., the core shift, one of the expected properties of optically-thick jet emission. Our results for PSR J1745-2900 support the hypothesis that Galactic Center pulsars will originate from the stellar disk and deepens the mystery regarding the small number of detected Galactic Center pulsars.

We study the dynamical evolution of the putative gas clouds G1 and G2 recently discovered in the Galactic center. Following earlier studies suggesting that these two clouds are part of a larger gas streamer, we combine their orbits into a single trajectory. Since the gas clouds experience a drag force from background gas, this trajectory is not exactly Keplerian. By assuming the G1 and G2 clouds trace this trajectory, we fit for the drag force they experience and thus extract information about the accretion flow at a distance of thousands of Schwarzschild radii from the black hole. This range of radii is important for theories of black hole accretion, but is currently unconstrained by observations. In this paper we extend our previous work by accounting for radial forces due to possible inflow or outflow of the background gas. Such radial forces drive precession in the orbital plane, allowing a slightly better fit to the G1 and G2 data. This precession delays the pericenter passage of G2 by 4-5 months relative to estimates derived from a Keplerian orbital fit; if it proves possible to identify the pericenter time observationally, this enables an immediate test of whether G1 and G2 are gas clouds part of a larger gas streamer. If G2 is indeed a gas cloud, its closest approach likely occurred in late summer 2014, after many of the observing campaigns monitoring G2's anticipated pericenter passage ended. We discuss how this affects interpretation of the G2 observations.

Gamma-ray observations have shown pulsars to be efficient converters of rotational energy into GeV photons and it is of wide-ranging interest to determine their contribution to the gamma-ray background. We arrive at flux predictions from both the young (<~ Myr) and millisecond (~Gyr) Galactic pulsar populations. We find that unresolved pulsars can yield both a significant fraction of the excess GeV gamma rays near the Galactic Center and an inverse Compton flux in the inner kpc similar to that inferred by Fermi. We compare models of the young pulsar population and millisecond pulsar population to constraints from gamma-ray and radio observations. Overall, we find that the young pulsars should outnumber millisecond pulsars as unassociated gamma-ray point sources in this region. The number of young radio pulsars discovered near the Galactic Center is in agreement with our model of the young pulsar population. Deeper radio observations at higher latitudes can constrain the total gamma-ray emission from both young and millisecond pulsars from the inner galaxy. While this is a step towards better understanding of pulsars, cosmic rays in the Milky Way, and searches for dark matter, we also discuss a few interesting puzzles that arise from the underlying physics of pulsar emission and evolution.

The redshifted 21 cm transition line of hydrogen tracks the thermal evolution of the neutral intergalactic medium (IGM) at "cosmic dawn," during the emergence of the first luminous astrophysical objects (~100 Myr after the Big Bang) but before these objects ionized the IGM (~400-800 Myr after the Big Bang). Because X-rays, in particular, are likely to be the chief energy courier for heating the IGM, measurements of the 21 cm signature can be used to infer knowledge about the first astrophysical X-ray sources. Using analytic arguments and a numerical population synthesis algorithm, we argue that the progenitors of supermassive black holes (SMBHs) should be the dominant source of hard astrophysical X-rays---and thus the primary driver of IGM heating and the 21 cm signature---at redshifts $z < 20$, if (i) they grow readily from the remnants of Population III stars and (ii) produce X-rays in quantities comparable to what is observed from active galactic nuclei and high-mass X-ray binaries. We show that models satisfying these assumptions dominate over contributions to IGM heating from stellar populations, and cause the 21 cm brightness temperature to rise at $z > 20$. An absence of such a signature in the forthcoming observational data would imply that SMBH formation occurred later (e.g. via so-called direct collapse scenarios), that it was not a common occurrence in early galaxies and protogalaxies, or that it produced far fewer X-rays than empirical trends at lower redshifts, either due to intrinsic dimness (radiative inefficiency) or Compton-thick obscuration close to the source.

The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form.

Studies of Fermi data indicate an excess of GeV gamma rays around the Galactic center (GC), possibly due to dark matter. We show that young gamma-ray pulsars can yield a similar signal. First, a high concentration of GC supernovae naturally leads to a population of kicked pulsars symmetric about the GC. Second, while very-young pulsars with soft spectra reside near the Galactic plane, pulsars with spectra that have hardened with age accumulate at larger angles. This combination, including unresolved foreground pulsars, traces the morphology and spectrum of the Excess.

We present three-dimensional magnetohydrodynamic simulations of magnetized gas clouds accelerated by hot winds. We initialize gas clouds with tangled internal magnetic fields and show that this field suppresses the disruption of the cloud: rather than mixing into the hot wind as found in hydrodynamic simulations, cloud fragments end up co-moving and in pressure equilibrium with their surroundings. We also show that a magnetic field in the hot wind enhances the drag force on the cloud by a factor ~(1+v_A^2/v_wind^2)$, where v_A is the Alfven speed in the wind and v_wind measures the relative speed between the cloud and the wind. We apply this result to gas clouds in several astrophysical contexts, including galaxy clusters, galactic winds, the Galactic center, and the outskirts of the Galactic halo. Our results can explain the prevalence of cool gas in galactic winds and galactic halos and how such cool gas survives in spite of its interaction with hot wind/halo gas. We also predict that drag forces can lead to a deviation from Keplerian orbits for the G2 cloud in the galactic center.

All gravitationally bound clusters expand, due to both gas loss from their most massive members and binary heating. All are eventually disrupted tidally, either by passing molecular clouds or the gravitational potential of their host galaxies. However, their interior evolution can follow two very different paths. Only clusters of sufficiently large initial population and size undergo the combined interior contraction and exterior expansion that leads eventually to core collapse. In all other systems, core collapse is frustrated by binary heating. These clusters globally expand for their entire lives, up to the point of tidal disruption. Using a suite of direct N-body calculations, we trace the "collapse line" in r_v-N space that separates these two paths. Here, r_v and N are the cluster's initial virial radius and population, respectively. For realistic starting radii, the dividing N-value is from 10^4 to over 10^5. We also show that there exists a minimum population, N_min, for core collapse. Clusters with N < N_min tidally disrupt before core collapse occurs. At the Sun's Galactocentric radius, R_G = 8.5 kpc, we find N_min >~ 300. The minimum population scales with Galactocentric radius as R_G^{-9/8}. The position of an observed cluster relative to the collapse line can be used to predict its future evolution. Using a small sample of open clusters, we find that most lie below the collapse line, and thus will never undergo core collapse. Most globular clusters, on the other hand, lie well above the line. In such a case, the cluster may or may not go through core collapse, depending on its initial size. We show how an accurate age determination can help settle this issue.

Pulsars orbiting the Galactic center black hole, Sgr A*, would be potential probes of its mass, distance and spin, and may even be used to test general relativity. Despite predictions of large populations of both ordinary and millisecond pulsars in the Galactic center, none have been detected within 25 pc by deep radio surveys. One explanation has been that hyperstrong temporal scattering prevents pulsar detections, but the recent discovery of radio pulsations from a highly magnetized neutron star (magnetar) within 0.1 pc shows that the temporal scattering is much weaker than predicted. We argue that an intrinsic deficit in the ordinary pulsar population is the most likely reason for the lack of detections to date: a "missing pulsar problem" in the Galactic center. In contrast, we show that the discovery of a single magnetar implies efficient magnetar formation in the region. If the massive stars in the central parsec form magnetars rather than ordinary pulsars, their short lifetimes could explain the missing pulsars. Efficient magnetar formation could be caused by strongly magnetized progenitors, or could be further evidence of a top-heavy initial mass function. Furthermore, current high-frequency surveys should already be able to detect bright millisecond pulsars, given the measured degree of temporal scattering.

We explore the tidal excitation of stellar modes in binary systems using Kepler observations of the remarkable eccentric binary KOI-54 (HD 187091; KIC 8112039), which displays strong ellipsoidal variation as well as a variety of linear and nonlinear pulsations. We report the amplitude and phase of over 120 harmonic and anharmonic pulsations in the system. We use pulsation phases to determine that the two largest-amplitude pulsations, the 90th and 91st harmonics, most likely correspond to axisymmetric m=0 modes in both stars, and thus cannot be responsible for resonance locks as had been recently proposed. We find evidence that the amplitude of at least one of these two pulsations is decreasing with a characteristic timescale of ~100 yr. We also use the pulsations' phases to confirm the onset of the traveling wave regime for harmonic pulsations with frequencies <~50 Omega_orbit, in agreement with theoretical expectations. We present evidence that many pulsations that are not harmonics of the orbital frequency correspond to modes undergoing simultaneous nonlinear coupling to multiple linearly driven parent modes. Since coupling among multiple modes can lower the threshold for nonlinear interactions, nonlinear phenomena may be easier to observe in highly eccentric systems, where broader arrays of driving frequencies are available. This may help to explain why the observed amplitudes of the linear pulsations are much smaller than the theoretical threshold for decay via three-mode coupling.

The centers of most known galaxies host supermassive black holes (SMBHs). In orbit around these black holes are a centrally-concentrated distribution of stars, both in single and in binary systems. Occasionally, these stars are perturbed onto orbits that bring them close to the SMBH. If the star is in a binary system, the three-body interaction with the SMBH can lead to large changes in orbital energy, depositing one of the two stars on a tightly-bound orbit, and its companion into a hyperbolic orbit that may escape the galaxy. In this Letter, we show that the disruption of solitary stars can also lead to large positive increases in orbital energy. The kick velocity depends on the amount of mass the star loses at pericenter, but not on the ratio of black hole to stellar mass, and are at most the star's own escape velocity. We find that these kicks are usually too small to result in the ejection of stars from the Milky Way, but can eject the stars from the black hole's sphere of influence, reducing their probability of being disrupted again. We estimate that ~10^5 stars, ~1% of all stars within 10 pc of the galactic center, are likely to have had mass removed by the central black hole through tidal interaction, and speculate that these "turbovelocity" stars will at first be redder, but eventually bluer, and always brighter than their unharrassed peers.

We perform a suite of cosmological simulations in the LCDM paradigm of the formation of the first structures in the Universe prior to astrophysical reheating and reionization (15<~z<200). These are the first simulations initialized in a manner that self consistently account for the impact of pressure on the rate of growth of modes, temperature fluctuations in the gas, and the dark matter-baryon supersonic velocity difference. Even with these improvements, these are still difficult times to simulate accurately as the Jeans length of the cold intergalactic gas must be resolved while also capturing a representative sample of the Universe. Our simulations support the finding of recent studies that the dark matter-baryon velocity difference has a surprisingly large impact on the accretion of gas onto the first star-forming minihalos (with masses of ~10^6 Msun). In fact, the halo gas is often significantly downwind of such halos and with lower densities, which delays the formation of the first stars in most locations in the Universe. We also show that dynamical friction plays an important role in the nonlinear evolution of the differential velocity, acting to erase this velocity difference quickly in overdense gas as well as sourcing visually-apparent bow shocks and Mach cones throughout the Universe. We use simulations with both the GADGET and Enzo cosmological codes to test the robustness of these conclusions. We find that particle coupling in GADGET between the gas and dark matter particles can result in spurious growth that mimics nonlinear growth in the matter power spectrum. In a companion paper, we use the simulations presented here to make detailed estimates for the impact of the dark matter--baryon velocity differential on redshifted 21cm radiation. The initial conditions generator used in this study CICsASS can be publicly downloaded.

Recently, Tseliakhovich and Hirata (2010) showed that during the cosmic Dark Ages the baryons were typically moving supersonically with respect to the dark matter with a spatially variable Mach number. Such supersonic motion may source shocks that heat the Universe. This motion may also suppress star formation in the first halos. Even a small amount of coupling of the 21cm signal to this motion has the potential to vastly enhance the 21cm brightness temperature fluctuations at 15<z<40 as well as to imprint acoustic oscillations in this signal. We present estimates for the size of this coupling, which we calibrate with a suite of cosmological simulations. Our simulations, discussed in detail in a companion paper, are initialized to self-consistently account for gas pressure and the dark matter-baryon relative velocity, v_bc (in contrast to prior simulations). We find that the supersonic velocity difference dramatically suppresses structure formation at 10-100 comoving kpc scales, it sources shocks throughout the Universe, and it impacts the accretion of gas onto the first star-forming minihalos (even for halo masses as large as ~10^7 Msun). However, we find that the v_bc-sourced temperature fluctuations can contribute only as much as ~10% of the fluctuations in the 21cm signal. We do find that v_bc could source an O(1) component in the power spectrum of the 21cm signal via the X-ray (but not ultraviolet) backgrounds produced once the first stars formed. In a scenario in which ~10^6 Msun minihalos reheated the Universe via their X-ray backgrounds, we find that the pre-reionization 21cm signal would be larger than previously anticipated and exhibit significant acoustic features. We show that structure formation shocks are unable to heat the Universe sufficiently to erase a strong 21cm absorption trough at z ~ 20 that is found in most models of the sky-averaged 21cm intensity.

During the formation of the Milky Way, > 100 central black holes (BHs) may have been ejected from their small host galaxies as a result of asymmetric gravitational wave emission. We previously showed that many of these BHs are surrounded by a compact cluster of stars that remained bound to the BH during the ejection process. In this paper, we perform long term N-body simulations of these star clusters to determine the distribution of stars in these clusters today. These numerical simulations, reconciled with our Fokker-Planck simulations, show that stellar density profile follows a power-law with slope ~ -2.15, and show that large angle scattering and tidal disruptions remove 20 - 90% of the stars by ~10^10 yr. We then analyze the photometric and spectroscopic properties of recoiled clusters accounting for the small number of stars in the clusters. We use our results to perform a systematic search for candidates in the Sloan Digital Sky Survey. We find no spectroscopic candidates, in agreement with our expectations from the completeness of the survey. Using generic photometric models of present day clusters we identify ~100 recoiling cluster candidates. Follow-up spectroscopy would be able to determine the nature of these candidates.

Stellar mass black holes (BHs) are expected to segregate and form a steep density cusp around supermassive black holes (SMBHs) in galactic nuclei. We follow the evolution of a multi-mass system of BHs and stars by numerically integrating the Fokker-Planck energy diffusion equations for a variety of BH mass distributions. We find that the BHs "self-segregate'', and that the rarest, most massive BHs dominate the scattering rate closest to the SMBH (< 0.1 pc). BH--BH binaries form out of gravitational wave emission during BH encounters. We find that the expected rate of BH coalescence events detectable by Advanced LIGO is ~1 - 100/yr, depending on the initial mass function of stars in galactic nuclei and the mass of the most massive BHs. We find that the actual merger rate is likely ~10 times larger than this due to the intrinsic scatter of stellar densities in many different galaxies. The BH binaries that form this way in galactic nuclei have significant eccentricities as they enter the LIGO band (90% with e > 0.9), and are therefore distinguishable from other binaries, which circularise before becoming detectable. We also show that eccentric mergers can be detected to larger distances and greater BH masses than circular mergers, up to ~ 700 solar masses. Future ground-based gravitational wave observatories will be able to constrain both the mass function of BHs and stars in galactic nuclei.

Gravitational wave emission by coalescing black holes (BHs) kicks the remnant BH with a typical velocity of hundreds of km/s. This velocity is sufficiently large to remove the remnant BH from a low-mass galaxy but is below the escape velocity from the Milky Way (MW) galaxy. If central BHs were common in the galactic building blocks that merged to make the MW, then numerous BHs that were kicked out of low-mass galaxies should be freely floating in the MW halo today. We use a large statistical sample of possible merger tree histories for the MW to estimate the expected number of recoiled BH remnants present in the MW halo today. We find that hundreds of BHs should remain bound to the MW halo after leaving their parent low-mass galaxies. Each BH carries a compact cluster of old stars that populated the core of its original host galaxy. Using the time-dependent Fokker-Planck equation, we find that present-day clusters are ~< 1 pc in size, and their central bright regions should be unresolved in most existing sky surveys. These compact systems are distinguishable from globular clusters by their internal (Keplerian) velocity dispersion greater than one hundred km/s and their high mass-to-light ratio owing to the central BH. An observational discovery of this relic population of star clusters in the MW halo, would constrain the formation history of the MW and the dynamics of BH mergers in the early Universe. A similar population should exist around other galaxies, and may potentially be detectable in M31 and M33.

Hypervelocity stars (HVSs) discovered in the Milky Way (MW) halo are thought to be ejected from near the massive black hole (MBH) at the galactic centre. In this paper we investigate the spatial and velocity distributions of the HVSs which are expected to be similarly produced in the Andromeda galaxy (M31). We consider three different HVS production mechanisms: (i) the disruption of stellar binaries by the galactocentric MBH; (ii) the ejection of stars by an in-spiraling intermediate mass black hole; and (iii) the scattering of stars off a cluster of stellar-mass black holes orbiting around the MBH. While the first two mechanisms would produce large numbers of HVSs in M31, we show that the third mechanism would not be effective in M31. We numerically calculate 1.2*10^6 trajectories of HVSs from M31 within a simple model of the Local Group and hence infer the current distribution of these stars. Gravitational focusing of the HVSs by the MW and the diffuse Local Group medium leads to high densities of low mass (~ solar mass) M31 HVSs near the MW. Within the virialized MW halo, we expect there to be of order 1000 HVSs for the first mechanism and a few hundred HVSs for the second mechanism; many of these stars should have distinctively large approach velocities (< -500 km/s). In addition, we predict ~5 hypervelocity RGB stars within the M31 halo which could be identified observationally. Future MW astrometric surveys or searches for distant giants could thus find HVSs from M31.

Stars within 0.1 pc of the supermassive black hole Sgr A* at the Galactic centre are expected to encounter a cluster of stellar-mass black holes (BHs) that have segregated to that region. Some of these stars will scatter off an orbiting BH and be kicked out of the Galactic centre with velocities up to ~2000 km/s. We calculate the resulting ejection rate of hypervelocity stars (HVSs) by this process under a variety of assumptions, and find it to be comparable to the tidal disruption rate of binary stars by Sgr A*, first discussed by Hills (1988). Under some conditions, this novel process is sufficient to account for all of the hypervelocity B-stars observed in the halo, and may dominate the production rate of all HVSs with lifetimes much less than the relaxation time-scale at a distance ~2 pc from Sgr A* (>~ 2 Gyr). Since HVSs are produced by at least two unavoidable processes, the statistics of HVSs could reveal bimodal velocity and mass distributions, and can constrain the distribution of BHs and stars in the innermost 0.1 pc around Sgr A*.

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