# Compact Binary Coalescences in the Band of Ground-based Gravitational-Wave Detectors

###### Abstract

As the ground-based gravitational-wave telescopes LIGO, Virgo, and GEO 600 approach the era of first detections, we review the current knowledge of the coalescence rates and the mass and spin distributions of merging neutron-star and black-hole binaries. We emphasize the bi-directional connection between gravitational-wave astronomy and conventional astrophysics. Astrophysical input will make possible informed decisions about optimal detector configurations and search techniques. Meanwhile, rate upper limits, detected merger rates, and the distribution of masses and spins measured by gravitational-wave searches will constrain astrophysical parameters through comparisons with astrophysical models. Future developments necessary to the success of gravitational-wave astronomy are discussed.

## 1 Introduction

We are likely only a few years away from the first direct detections of gravitational waves from coalescing compact-object (CO) binaries. The Laser Interferometer Gravitational-Wave Observatory (LIGO, http://www.ligo.caltech.edu, [82]) has completed two years of data taking at the designed sensitivity (the S5 run) in the fall of 2007. Simultaneously, the Virgo detector (http://www.virgo.infn.it, [7]) ended its VSR1 run. Other ground-based detectors in operation include GEO 600 [36]. By 2014-2015, the Advanced versions of the LIGO and Virgo detectors are expected to operate at sensitivities – times greater than Initial LIGO [83]. These instruments will be able to detect gravitational waves (GWs) from inspirals, mergers, and ringdowns associated with coalescences of compact binary systems composed of neutron stars (NSs) or black holes (BHs). Although estimates of the rates of such coalescences cover a rather wide range, Advanced LIGO may be able to detect tens or even hundreds of these events per year (see, e.g., O’Shaughnessy et al. [63]).

In this paper, we review the mechanisms of formation of coalescing NS and BH binaries and summarize the predictions for the merger rates of NS-NS, NS-BH, and BH-BH systems. For NS-NS systems these predictions can be made by extrapolating from the small sample of observed binary pulsars, while for other binary types they are based on population-synthesis models of isolated binary evolution. We discuss the potential importance of the dynamical formation of BH-BH binaries in dense stellar environments, such as globular and nuclear star clusters, and the possible coalescenses involving intermediate-mass black holes (IMBHs). We summarize the existing knowledge about the primary characteristics of NSs and BHs in binaries, namely, their masses and spins.

We focus on the tight connection between (i) astrophysical observations and modeling and (ii) gravitational-wave searches and data analysis. Input from astrophysics can inform decisions about which detector configurations should be chosen, which searches are the most likely to yield interesting results, and which waveform families should be used for detection. On the other hand, looking beyond the excitement of the first detection, the detected merger rate (or even non-trivial upper limits) will allow us to limit the space of allowed astrophysical models, thereby providing constraints on currently ill-determined astrophysical parameters. Even more precise constraints will come from matching the observed parameters of the binaries, particularly their masses and spins, with the model predictions.

## 2 Mechanisms of Formation and Rate Predictions

The standard formation channel for compact-object binaries in the field is isolated binary evolution. In this case, a primordial binary of main-sequence stars evolves through several stages of mass transfer and possibly common-envelope evolution, while the binary components age and eventually undergo supernovae, occasionally leaving behind a binary composed of neutron stars or black holes (see [41] for a review of this process). A circular binary will merge through radiation-reaction from gravitational-wave emission on a timescale of [70]

(1) |

where and are the component masses and is the orbital radius. In this section, we describe two different ways to evaluate the merger rates for CO binaries: directly from observations for NS-NS systems and via population-synthesis modeling for other binary types. We also discuss dynamical formation scenarios that could be significant for BH-BH binaries in dense stellar environments.

### 2.1 Binary pulsars: observational constraints and extrapolated NS-NS rates

Among double-CO binaries, NS-NS binaries are the only ones for which we have direct observational evidence. This evidence comes from observations of Galactic binary pulsars: double neutron-star systems where at least one of the NSs is a pulsar. Approximately ten such systems have been detected so far. By modeling the selection effects of the searches which detected these systems, it is possible to estimate the total rate of merging NS-NS binaries in the galaxy, whether or not they are detected as binary pulsars (the pioneering efforts to do this are described in [71] and [57]). Five of the detected systems are expected to merge within a Hubble time through radiation reaction from the emission of gravitational waves, and therefore contribute to extrapolated merger rate estimates. These systems are: B1913+16 (the Hulse-Taylor pulsar), B1534+12, J0737-3039 (a double pulsar), J1756-2251, and J1906+0746; there is also a sixth merging system, located in a globular cluster, which is thought to be atypical and is not generally included in merger rate calculations.

Kim et al. [43] developed the statistical framework necessary to reconstruct the Galactic NS-NS binary population from observations, while rigorously estimating the uncertainties which arise due to the small number of observed systems. Kalogera et al. [42] reported the extrapolated rates from the first three systems of the five listed above. They used a variety of models for the pulsar luminosity distribution, which plays a very significant role in determining the selection effects. The distribution is described by two variables: the minimum pulsar luminosity and the slope of the luminosity distribution power law . For the preferred pulsar luminosity distribution with mJy kpc and (see model 6 in Table 1 of [42]), the authors found a most likely rate of NS-NS mergers per Myr in the Galaxy; the fifth and 95th percentile statistical rate estimates were and per Myr, respectively.

Although such extrapolation has significantly fewer free parameters than the population-synthesis models described below, the numbers should nevertheless be used with caution. Different choices of the variables describing the pulsar luminosity distribution, which are all consistent with the observations, could shift the merger rates by an order of magnitude. The fifth percentile for the model yielding the most pessimistic merger rates suggests only merger per Myr in the Galaxy, while rates as high as per Myr are possible at the 95th percentile for the most optimistic model. There are several possibilities for systematic biases related to the imperfect knowledge of pulsar age and the beaming fraction [65]. The implicit assumption that the observed pulsars form a good representation of the overall Galactic NS-NS population is difficult to verify. Lastly, the small number of observed systems means that a single system can significantly change the rates predictions. For example, according to [44], the inclusion of the relatively short-lived system J1906+0746 could increase the rate estimates above by almost a factor of .

### 2.2 Population-synthesis models for isolated binaries

In the absence of direct observations of NS-BH and BH-BH binaries, the only way to predict the rates of binary mergers is through population-synthesis modeling. A number of different models has been applied over the past two decades (see, e.g., [35, 73, 41] and references therein). Generally, these models rely on a number of ill-constrained astrophysical parameters. It is therefore necessary to explore the full allowed parameter space in order to quantify the uncertainties in the merger rates predictions. We follow the procedure described in [67], since most other models do not attempt a thorough exploration of the full parameter space.

The StarTrack population-synthesis code used in [67] has seven free parameters that can significantly affect the model outcomes, as identified in [13]. These include: the power-law index in the binary mass ratio; 3 parameters used to describe the supernova kick velocity distribution; the strength of the massive stellar wind; common-envelope efficiency; and the fractional mass loss during non-conservative mass transfer. Flat priors are used on these parameters to sample the widest possible set of models. However, several additional constraints are then applied to disallow those models (i.e., those combinations of the seven input parameters) that disagree with observations. The five constraints applied in [67] come from the following observations: (i) merging Galactic NS-NS binaries (i.e., binary pulsars expected to merge in a Hubble time); (ii) wide NS-NS binaries (those not expected to merge in a Hubble time); (iii) white dwarf - NS binaries; (iv) observed Type II supernova rate; and (v) observed Type Ib/c supernova rate. Only the first three observations provide non-trivial constraints. The predicted Galactic rates for NS-NS coalescences range from to per Myr, with the most likely value at per Myr [67]. For NS-BH binaries, the range is from to per Myr, with the most likely value at per Myr [67]. And for BH-BH binaries, model predictions range from to per Myr, with a most likely value of per Myr [41].

In addition to the ranges stated above, which reflect the impact of different choices of astrophysical model parameters, further uncertainties exist. For instance, if hypercritical accretion during the common-envelope phase is disallowed, merger rates can drop significantly [95, 14]. Many potential BH-BH progenitors may enter the common-envelope phase while the donor star is evolving through the Hertzsprung gap; according to [14], such systems may merge directly, thus significantly reducing the rate of BH-BH mergers. If kicks are preferentially aligned with pre-supernova spin, NS-NS coalescence rates could decrease by up to a factor of five [74]. Additionally, uncertainties in the star-forming conditions, particularly metallicity, can lead to systematic errors in the rates predictions at least at the -% level [68]. Recent calculations using updated models of stellar winds and direct observations of the black-hole X-ray binary IC10 X-1 suggest that low-metallicity environments may be even more important [19, 12]. Meanwhile, all of the applied observational constraints suffer from the uncertainties and possible biases discussed in the previous subsection. Nonetheless, the default rate predictions for NS-BH and BH-BH merger from all of the population-synthesis models mentioned above, as well as those in [58] and [26], fall within the uncertainty range listed in the previous paragraph, suggesting that predictions generally converge within the uncertainties.

The rates quoted in Sections 2.1 and 2.2 are computed for our Galaxy. To convert them to rates for other galaxies, they can be scaled by the star formation rate, which is tracked by the blue-light luminosity for spiral galaxies like the Milky Way [45]. However, while this tracer is relatively accurate for nearby galaxies, it fails to account for mergers that are significantly delayed relative to star formation in older elliptical galaxies [25]. This omission is particularly significant for BH-BH mergers, and therefore the actual BH-BH merger rate in the volume to which ground-based detectors will be sensitive in the Advanced LIGO-Virgo era can be significantly higher than the contribution from spiral galaxies alone would indicate [63].

### 2.3 Dynamical formation scenarios

Dynamical effects can make a significant contribution to CO binary mergers in dense stellar environments, particularly globular star clusters and nuclear clusters at the centers of galaxies. There, two-body encounters can lead to direct captures, while three-body single-binary interaction and four-body binary-binary interactions can tighten hard binaries to a point where radiation-reaction from gravitational-wave emission can take over and drive the binary to coalescence.

There have been some suggestions regarding the contribution of dynamically-formed binaries in globular clusters to short hard gamma ray bursts (SGRBs) [39, 37], which are believed to be powered by mergers of NS-NS or NS-BH binaries. However, extrapolating from the rate of observed SGRBs, as attempted, for example, in [27], is highly uncertain, because of unknown selection effects such as SGRB beaming fractions [55], which may be different for field and cluster SGRBs [34]; additionally, binary mergers involving NSs are not the only possible formation scenario for SGRBs [94, 46]. In fact, most simulations indicate that dynamical effects are not important for binaries containing neutron stars, as heavier black holes will segregate to the dense cores and substitute in for neutron stars in binaries [79, 40]. In particular, Grindlay et al. [34] estimate the NS-NS merger rate from all globular clusters in the Galaxy to be 40 per Gyr, a factor of less than predictions for field mergers.

On the other hand, dynamical effects can make a very significant contribution to BH-BH merger rates, since black holes are the more massive and will sink to the central, dense regions of the clusters. O’Leary et al. [61] and Sadowski et al. [79] both examined BH-BH binary formation in globular clusters. Their estimates range over several orders of magnitude, from to BH-BH mergers per globular cluster per Myr, due to different assumptions about the fraction of star formation that occurs in clusters, the fraction of clusters that may form black-hole subclusters, and whether these BH subclusters decouple or remain in thermal equilibrium with other stars in the cluster. Similarly, Banerjee et al. [11] find that 3-body encounters in intermediate-age (a few Gyr) massive ( ) clusters lead to the formation of a few merging BH-BH binaries per cluster, some of which coalesce after ejection from the cluster. O’Leary et al. [60] argue that the presence of a massive black hole in galactic nuclear clusters can create a sufficiently dense cusp to allow as many as BH-BH binaries per Gyr to form in such a cluster through direct 2-body encounters. These values may be based on very optimistic assumptions regarding the fraction of BHs in galactic nuclei and the number density of low-mass galaxies, but if such BH binaries did form, they could be easily distinguished due to their significant eccentricity, which is not expected for other stellar-mass binaries in the band of ground-based detectors. Meanwhile, Miller and Lauburg [54] calculate that a few BH-BH mergers per Myr per galaxy should occur in nuclear clusters of small galaxies without massive black holes, primarily through three-body interactions.

If intermediate-mass black holes (IMBHs) with a mass around are generic in globular clusters (see [53, 52] for reviews on IMBHs, and [28] for an announcement of a recently discovered ultra-luminous X-ray source that represents a possible IMBH detection), they can contribute to binary merger rates. Since an IMBH will be the most massive object in the cluster, it will readily sink to the center and substitute into a binary with a compact-object companion. The binary will then harden through three-body interactions and eventually merge via an intermediate-mass-ratio inspiral on timescales of less than one billion years [50]. If the binary fraction in a young dense cluster is % and the core-collapse time is Myr, two IMBHs can form in a cluster, yielding an IMBH-IMBH merger [29]. Such a merger is also possible through the merger of two distinct globular clusters, each containing an IMBH [9]. While Advanced LIGO and Virgo will be sensitive to coalescences involving IMBHs in globular clusters, third-generation detectors like the Einstein Telescope, with improved low-frequency sensitivity [30], could detect high-redshift mergers of IMBHs formed from population III stars which may be the seeds of today’s massive black holes [81, 33, 32]. However, the present lack of knowledge about the existence and occupation fraction of IMBHs in globular clusters and about early structure formation makes very speculative all rates predictions for mergers involving IMBHs.

## 3 Expected characteristics of binaries

The astrophysical characteristics of merging binaries, namely, their mass and spin distributions, are of great interest to gravitational-wave astrophysicists. However, these are not well understood. As with merger rates, there are two ways to constrain the binary parameters: through direct observations, and through modeling. As discussed above, direct observations are limited to a few systems for NS-NS binaries, and non-existent for NS-BH and BH-BH binaries. And although BH mass and spin measurements exist from mass-transfering X-ray binaries, these do not necessarily have the same mass and spin distributions as BHs in compact-object binaries because of the particular evolutionary histories of double-CO binaries. Meanwhile, models are sensitive to assumptions and can differ significantly in their predictions. We discuss what is presently known about mass and spin distributions for NSs and BHs below.

### 3.1 NS and BH masses

The masses of pulsars in binaries have been estimated to date [47, 48]. Until a decade ago, most observed binary pulsar masses were extremely consistent with a “canonical” value of [85]. Since then, several significantly more massive pulsars have been found with much larger confirmed masses, such as J1903+0327 with a confirmed timing mass [24]. In globular clusters, there are four pulsars in binaries with weakly constrained orbital inclinations whose masses appear to be higher than “canonical”, from for a pulsar in Terzan 5 up to for J1748-2021B; see Figure 28 in [47].

More than 40 black-hole X-ray binaries have been detected to date. The masses of approximately 20 of them have been measured by observing the companion dynamics [76, 62]. These masses vary between and solar masses. However, because of the mass-transfer episodes involving common envelopes, which play a crucial role in the formation of NS-BH and BH-BH binaries, these masses do not necessarily reflect the typical BH masses in such binaries.

In the absence of direct evidence, masses can be predicted by population-synthesis models. Most studies in the literature focus on predictions of the chirp mass, , since this parameter determines the gravitational-wave amplitude. The chirp mass is of the total mass for for equal-mass binaries. It is typically found to be between and solar masses for NS-BH systems and between and solar masses for BH-BH systems [18, 66, 14, 67, 63]. However, mean chirp mass values vary by % between these predictions depending on assumptions about the cutoff NS mass, accretion rate, etc. The mass ratio is typically very poorly constrained in existing population-synthesis studies (see, e.g., Fig. 4 of [67]).

Meanwhile, a number of channels have been suggested through which significantly more massive black-hole binaries could form. Recent observations of IC10 X-1, a binary with a massive black-hole accreting from a Wolf-Rayet companion star, suggest that more massive BH-BH binaries can form through isolated binary evolution, with chirp masses [20]. Dynamical formation channels will favor relatively massive BH-BH binary components, as the most massive black holes are likely to sink to the centers of clusters through mass segregation and substitute into binaries during three-body encounters. Moreover, the BH merger products in such dense clusters can be reused if they are not ejected from the cluster due to recoil kicks, leading to higher-mass mergers in subsequent generations; chirp masses for BH-BH mergers in globular clusters can therefore range to [79]. Additionally, although stellar winds in high-metallicity environments may prevent the formation of massive black holes, mass loss through stellar winds would be much less significant in low-metallicity environments, allowing more massive black holes to form [12]. Finally, IMBHs discussed above could range to hundreds or thousands of solar masses, although the more massive IMBHs will not be detectable with ground-based detectors.

### 3.2 NS and BH spins

The observed pulsar spin periods and assumptions about NS spindown rates place the NS spin periods at birth in the range ms [47], which corresponds to dimensionless spins . Meanwhile, the fastest known recycled pulsar in a double NS binary has a spin of [23]. Thus, the spins of neutron stars are likely to be irrelevant for gravitational-wave detection and parameter estimation.

On the other hand, direct measurements of BH spins in X-ray binaries through quasi-periodic oscillations [86], continuum spectra of the inner edge of the accretion disk [56], or iron-line profiles [51] indicate that a significant fraction of BHs are rapidly spinning with , with some apparently very near to maximum spin, although a few could have low spins . However, all of three methods for spin measurement potentially suffer from model-dependent assumptions. Moreover, these spins could have been enhanced by accretion, and may therefore not be relevant for birth spins in NS-BH or BH-BH binaries.

Population-synthesis models indicate that a majority of both NS-BH and BH-BH systems undergo a common-envelope phase with hypercritical accretion, i.e., highly super-Eddington accretion with neutrino cooling [41] (however, see [14]). The details of hypercritical accretion, including the associated spin-up of the black hole, are poorly understood. Under the assumption that all of the angular momentum at the last stable orbit is transferred to the BH, which maximizes the BH spin, O’Shaughnessy et al. [64] found typical BH spins in access of , although spins tended to get smaller for more massive BHs, which are more difficult to spin up. However, using more conservative assumptions about accretion and spin up, Belczynski et al. [15] conclude that unless a BH has a significant birth spin, it is unlikely to be spinning above .

The processes governing spin orientations are not fully understood at present. For some isolated binaries, the spins may be preferentially aligned with each other and with the orbital angular momentum. This will happen if all three of the following conditions are satisfied: (i) if the binary was aligned before the supernovae (e.g., through mass transfer episodes), (ii) if the supernovae kicks are small enough that the post-kick orbital plane is not too different from the pre-kick one (this is particularly likely for the more massive BHs, which are likely to have lower kicks, or for binaries that already have very short orbital periods before the second supernova [20], and (iii) if the asymmetric supernovae themselves don’t change the spin direction relative to the spin of the progenitor star. The conditions may be slightly relaxed for the primary, as its spin can be re-aligned with the orbital plane through accretion and the Bardeen-Peterson effect. For dynamically formed binaries, there is no clear mechanism for preferential spin orientations, so the spin directions should be isotropic.

## 4 Gravitational-wave observations and astrophysics

### 4.1 Why do we need astrophysics to inform GW searches?

Astrophysical inputs can significantly aid searches for gravitational waves. Perhaps the most obvious such input has to do with the selection of the optimal configurations for advanced detectors. Advanced LIGO will have several possibly tunable parameters, including the laser power, the reflectivity of the signal recycling mirror and particularly the detuning of the signal recycling cavity. Choices of these parameters can dramatically change the shape of the noise spectra, as shown in the Advanced LIGO Reference Design [3]. For example, the figure in [1] shows six different possible configurations. Some of them trade off high-frequency sensitivity for improved sensitivity at low frequencies, or broad-band sensitivity for a lower noise floor in a narrow frequency band. Future detectors, such as the Einstein Telescope [2], currently have an even greater range of possible configurations. Predictions for the rates of different classes of compact-binary coalescences and the astrophysical significance of specific searches will be very important in selecting detector configurations, particularly in the era after the first detections.

Similarly, astrophysical input can aid decisions regarding the allocation of limited human and computational resources in data analysis: Should searches be expanded to the mass range of intermediate-mass black holes? Given scarce resources, does it make sense to develop a new search for eccentric binaries at the expense of possible minor improvements in searches for circular binaries? And what waveforms should be used for existing searches [5, 6, 80]? For instance, the relatively low masses predicted for most NS-BH and BH-BH mergers suggest that inspiral-only waveforms are sufficient for such searches, but templates that include merger and ringdown signals [22, 8] will be necessary to search for mergers involving IMBHs. Although current searches do not include spins because existing non-physical spinning template families increase the false alarm rates unacceptably [87], the likelihood of significant BH spins requires spinning templates to be used for parameter estimation [90]. And the computational difficulty of including spinning waveforms in the template banks being used for detection may be partially circumvented by ignoring those regions of the parameter space which are astrophysically unlikely, such as NS-BH binaries with very rapidly spinning () massive () black holes [69, 64].

### 4.2 What can GW searches teach us about astrophysics?

Once the first detections of gravitational waves are achieved, GW astronomy will allow us to utilize this new observational window to study the astrophysics of stellar and binary evolution and cluster dynamics.

Some individual detections will already carry very exciting astrophysical information. For instance, a single detection with an accurate BH mass measurement of will provide an unambiguous proof of the existence of an intermediate-mass black hole. Any measurement indicating a significant eccentricity in a coalescing binary would indicate that dynamical effects played an important role in the binary’s formation. On the other hand, statistics from multiple observations can be compared with models to constrain astrophysical parameters.

Perhaps the most obvious statistic is the observed rate of binary coalescences. As discussed in [68], detections will constrain the rate of mergers to within a factor , according to Poisson statistics. A Bayesian framework can be used to convert observed rates into additional constraints on the population-synthesis models discussed in Sec. 2.2, thereby making it possible to more accurately constrain the underlying astrophysical parameters. For a given choice of astrophysical model parameters , population synthesis codes coupled to information about galaxy distributions and detector sensitivity provide a distribution of the detectable event rate, . If an actual merger rate is measured for a given binary type, then the likelihood that a model with parameters fits the data is

(2) |

where includes both the statistical uncertainty in the measurement (equal to ), and the systematic uncertainty in the models. Then the constrained posterior probability density function (PDF) of astrophysical parameters is given by the Bayes’ rule,

(3) |

where the prior represents any a priori knowledge about the parameters, and is the evidence which is set by the overall normalization requirement.

Even in the absence of detections, stringent upper limits can significantly inform our understanding of the astrophysical model parameters. Again, a Bayesian approach can be used. However, for illustration purposes, we can consider a simple example, where the upper limit is employed as a firm cutoff. Consider, for example, the population-synthesis prediction for the BH-BH merger rate per Milky Way Equivalent Galaxy (MWEG), constrained by the extrapolated observations of merging binary pulsars [66]; this rate is plotted in blue in the left-hand panel of Figure 1. If observations placed a firm upper limit of 1 BH-BH merger per MWEG per Myr, models predicting higher rates are ruled out, resulting in the shifted distribution shown in red. The right-hand panel of Figure 1 shows the distribution of the dimensionless parameter describing the strength of the stellar wind. Its prior distribution is assumed to be uniform in . The blue histogram indicates that binary-pulsar observations do not constrain this astrophysical parameter, as the distribution remains flat (the drops near 0 and 1 are features of the plotting routine). However, as the red histogram indicates, the additional constraint from the putative upper limit on the BH-BH merger rate would preferentially rule out models with low values of the stellar-wind strength, thereby significantly informing our knowledge about this astrophysical parameter.

In addition to rates, the distributions of the measured parameters of the coalescing binaries, particularly masses and spins, can be matched to the model predictions. This approach requires several steps. First, for each detected signal, the binary parameters have to be measured. This involves measuring not only the best-fit parameters, but also estimating the accuracy of the measurement. It is best accomplished by sampling the posterior PDF in the parameter space with a Markov Chain Monte Carlo [77, 89, 90] or an alternative Bayesian sampling technique [91]; such techniques have been demonstrated to work for the full 15-dimensional parameter space of spinning LIGO waveforms [75]. Next, the PDFs from multiple detections have to be combined into a single statement about the parameter distribution in the underlying population of CO binaries [49], while taking into account selection biases, such as the mass dependence of the detection range [16]. At this stage, comparisons can be made to the distributions predicted by population-synthesis models.

The fully Bayesian machinery described above is essential to exploit the information available in most parameter measurements, as the low-amplitude detections likely in second-generation detectors will not tightly confine their permissible ranges. The chirp mass, however, is measured more precisely than the scatter expected between independent samples drawn from the (detected) mass distribution. In this exceptional case the number of samples, rather than individual measurements’ accuracy, limits our ability to reconstruct the intrinsic distribution. Assuming we do not attempt to exploit other degrees of freedom, the sampled chirp mass distribution can therefore be directly compared with the simulated distribution; see for example [21]. A preliminary self-consistent treatment of both the number of detections and their chirp mass distribution is shown in Figure 2.

Finally, a wealth of information can be gained from simultaneous gravitational and electromagnetic observations of binary coalescences, such as the observations of gravitational waves associated with short GRBs [4]. We refer the reader to [72, 84] and references therein for detailed discussions of the scientific promise of multi-messenger astronomy, which spans the range from exposing GRB progenitors to enabling precise cosmography [59].

## 5 Future developments

There is a great deal of ongoing work to improve the detection and parameter-estimation pipelines within the LIGO-Virgo collaboration. These efforts range from the inclusion of new waveform templates and detection statistics in the pipeline [80] to improved sampling techniques for parameter estimation [88, 92] to simulated runs on numerical-relativity injections as part of the NINJA project [10].

In order to take full advantage of the data that will soon come from gravitational-wave detectors, more work is also needed on the astrophysical side. Larger population-synthesis runs with improved models are needed to enable comparisons against observational results. Most simulations of dynamical-formation scenarios relevant to ground-based detectors have not yet quantitatively explored the uncertainties in the predicted rates. In general, detailed studies of all input uncertainties in the models are needed.

We need to prepare for the era of routine detections by developing a framework for converting observed LIGO-Virgo merger rates or upper limits for compact-binary sources into improved measurements of astrophysical parameters. In order to utilize the full wealth of information that will be available in the measured binary parameters, such as masses and spins, we must also be ready to combine the parameter PDFs from individual events into statements about the population distribution while carefully taking selection biases into account, and to compare these observed distributions against those predicted by astrophysical models.

Another direction of theoretical work relates to probing strong-field gravity, and perhaps even testing general relativity, with gravitational-wave observations. There is a growing body of work on LISA extreme-mass-ratio inspirals as probes of the spacetimes of massive compact bodies (see, e.g., [78, 31, 93]). However, it has been suggested that even ground-based detectors can be useful tools for studying possible deviations from the Kerr metric. For example, gravitational waves from intermediate-mass-ratio inspirals (IMRIs) of NSs and BHs into IMBHs may encode enough information to allow comparisons of the mass-quadrupole moment of the IMBH against its expected Kerr value, thereby measuring the “bumpiness” of the black hole [17]. Further development effort is needed to make such measurements part of the standard data analysis of detected signals. Finally, there are ongoing efforts to study the ability of ground-based detectors to probe the equation of state of neutron stars through observations of NS-BH coalescences [38].

## Acknowledgments

IM is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-0901985. IM acknowledges support by NSF grant PHY-0653321 to Northwestern University for part of this work. R.O. was supported by National Science Foundation awards PHY 06-53462. Both authors appreciate the support of the NRDA-2 conference organizers.

## References

## References

- [1] Advanced LIGO anticipated sensitivity curves, https://dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid=2974.
- [2] The Einstein Telescope Project, https://www.et-gw.eu/et/.
- [3] Advanced LIGO Reference Design, 2007, LIGO Tech. Rep. M060056, https://dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid=1507.
- [4] B. Abbott et al., Implications for the Origin of GRB 070201 from LIGO Observations, Astrophysical Journal 681 (2008), 1419–1430.
- [5] B. Abbott et al., Search for gravitational waves from low mass binary coalescences in the first year of LIGO’s S5 data, Phys. Rev. D 79 (2009), no. 12, 122001.
- [6] , Search for Gravitational Waves from Low Mass Compact Binary Coalescence in 186 Days of LIGO’s fifth Science Run, Phys. Rev. D 80 (2009), 047101.
- [7] F. Acernese et al., Status of Virgo, Classical and Quantum Gravity 25 (2008), no. 11, 114045.
- [8] P. Ajith, S. Babak, Y. Chen, M. Hewitson, B. Krishnan, A. M. Sintes, J. T. Whelan, B. Brügmann, P. Diener, N. Dorband, J. Gonzalez, M. Hannam, S. Husa, D. Pollney, L. Rezzolla, L. Santamaría, U. Sperhake, and J. Thornburg, Template bank for gravitational waveforms from coalescing binary black holes: Nonspinning binaries, Phys. Rev. D 77 (2008), no. 10, 104017.
- [9] P. Amaro-Seoane and M. Freitag, Intermediate-Mass Black Holes in Colliding Clusters: Implications for Lower Frequency Gravitational-Wave Astronomy, Astrophysical Journal 653 (2006), L53–L56.
- [10] B. Aylott et al., Testing gravitational-wave searches with numerical relativity waveforms: results from the first Numerical INJection Analysis (NINJA) project, Classical and Quantum Gravity 26 (2009), no. 16, 165008.
- [11] S. Banerjee, H. Baumgardt, and P. Kroupa, Stellar-mass black holes in star clusters: implications for gravitational wave radiation, MNRAS 402 (2010), 371–380.
- [12] K. Belczynski, T. Bulik, C. L. Fryer, A. Ruiter, J. S. Vink, and J. R. Hurley, On The Maximum Mass of Stellar Black Holes, ArXiv e-prints (2009), 0904.2784.
- [13] K. Belczynski, V. Kalogera, F. A. Rasio, R. E. Taam, A. Zezas, T. Bulik, T. J. Maccarone, and N. Ivanova, Compact Object Modeling with the StarTrack Population Synthesis Code, Astrophysical Journal Supplement 174 (2008), 223–260.
- [14] K. Belczynski, R. E. Taam, V. Kalogera, F. A. Rasio, and T. Bulik, On the Rarity of Double Black Hole Binaries: Consequences for Gravitational Wave Detection, Astrophysical Journal 662 (2007), 504–511.
- [15] K. Belczynski, R. E. Taam, E. Rantsiou, and M. van der Sluys, Black Hole Spin Evolution: Implications for Short-Hard Gamma-Ray Bursts and Gravitational Wave Detection, Astrophysical Journal 682 (2008), 474–486.
- [16] P. R. Brady and S. Fairhurst, Interpreting the results of searches for gravitational waves from coalescing binaries, Classical and Quantum Gravity 25 (2008), no. 10, 105002.
- [17] D. A. Brown, J. Brink, H. Fang, J. R. Gair, C. Li, G. Lovelace, I. Mandel, and K. S. Thorne, Prospects for Detection of Gravitational Waves from Intermediate-Mass-Ratio Inspirals, Physical Review Letters 99 (2007), no. 20, 201102.
- [18] T. Bulik and K. Belczyński, Constraints on the Binary Evolution from Chirp Mass Measurements, Astrophysical Journal 589 (2003), L37–L40.
- [19] T. Bulik and K. Belczynski, Observational evidence for stellar mass binary black holes and their coalescence rate, ArXiv e-prints (2009), 0907.5515.
- [20] T. Bulik, K. Belczynski, and A. Prestwich, IC10 X-1: the immediate progenitor of a double black hole binary, ArXiv e-prints (2008), 0803.3516.
- [21] T. Bulik, K. Belczyński, and B. Rudak, Astrophysical significance of the detection of coalescing binaries with gravitational waves, A&A 415 (2004), 407–414.
- [22] A. Buonanno, Y. Pan, J. G. Baker, J. Centrella, B. J. Kelly, S. T. McWilliams, and J. R. van Meter, Approaching faithful templates for nonspinning binary black holes using the effective-one-body approach, Phys. Rev. D 76 (2007), no. 10, 104049.
- [23] Marta Burgay et al., An increased estimate of the merger rate of double neutron stars from observations of a highly relativistic system, Nature. 426 (2003), 531–533.
- [24] D. J. Champion, S. M. Ransom, P. Lazarus, F. Camilo, C. Bassa, V. M. Kaspi, D. J. Nice, P. C. C. Freire, I. H. Stairs, J. van Leeuwen, B. W. Stappers, J. M. Cordes, J. W. T. Hessels, D. R. Lorimer, Z. Arzoumanian, D. C. Backer, N. D. R. Bhat, S. Chatterjee, I. Cognard, J. S. Deneva, C.-A. Faucher-Giguère, B. M. Gaensler, J. Han, F. A. Jenet, L. Kasian, V. I. Kondratiev, M. Kramer, J. Lazio, M. A. McLaughlin, A. Venkataraman, and W. Vlemmings, An Eccentric Binary Millisecond Pulsar in the Galactic Plane, Science 320 (2008), 1309–1312.
- [25] J.A de Freitas Pacheco, T Regimbau, Vincent S., and A. Spallicci, Expected coalescence rates of ns-ns binaries for laser beam interferometers, Int.J.Mod.Phys. D15 (2006), 235–250.
- [26] J. D. M. Dewi, P. Podsiadlowski, and A. Sena, Double-core evolution and the formation of neutron star binaries with compact companions, MNRAS 368 (2006), 1742–1748.
- [27] A. Dietz, Chance Estimations for Detecting Gravitational Waves with LIGO/Virgo Associated with Gamma Ray Bursts, ArXiv e-prints (2009), 0904.0347.
- [28] S. A. Farrell, N. A. Webb, D. Barret, O. Godet, and J. M. Rodrigues, An intermediate-mass black hole of over 500 solar masses in the galaxy ESO243-49, Nature 460 (2009), 73–75.
- [29] J. M. Fregeau, S. L. Larson, M. C. Miller, R. O’Shaughnessy, and F. A. Rasio, Observing IMBH-IMBH Binary Coalescences via Gravitational Radiation, Astrophysical Journal Letters 646 (2006), L135–L138.
- [30] A. Freise, S. Chelkowski, S. Hild, W. Del Pozzo, A. Perreca, and A. Vecchio, Triple Michelson interferometer for a third-generation gravitational wave detector, Classical and Quantum Gravity 26 (2009), no. 8, 085012.
- [31] J. R. Gair, C. Li, and I. Mandel, Observable properties of orbits in exact bumpy spacetimes, Phys. Rev. D 77 (2008), no. 2, 024035.
- [32] J. R. Gair, I. Mandel, M. C. Miller, and M. Volonteri, Exploring intermediate and massive black-hole binaries with the Einstein Telescope, ArXiv e-prints (2009), 0907.5450.
- [33] J. R. Gair, I. Mandel, A. Sesana, and A. Vecchio, Probing seed black holes using future gravitational-wave detectors, Classical and Quantum Gravity 26 (2009), no. 20, 204009.
- [34] J. Grindlay, S. Portegies Zwart, and S. McMillan, Short gamma-ray bursts from binary neutron star mergers in globular clusters, Nature Physics 2 (2006), 116–119.
- [35] L. P. Grishchuk, V. M. Lipunov, K. A. Postnov, M. E. Prokhorov, and B. S. Sathyaprakash, REVIEWS OF TOPICAL PROBLEMS: Gravitational wave astronomy: in anticipation of first sources to be detected, Physics Uspekhi 44 (2001), 1.
- [36] H. Grote and the LIGO Scientific Collaboration, The status of GEO 600, Classical and Quantum Gravity 25 (2008), no. 11, 114043.
- [37] D. Guetta and L. Stella, Short -ray bursts and gravitational waves from dynamically formed merging binaries, A&A 498 (2009), 329–333.
- [38] T. Hinderer, B. D. Lackey, R. N. Lang, and J. S. Read, Tidal deformability of neutron stars with realistic equations of state and their gravitational wave signatures in binary inspiral, ArXiv e-prints (2009), 0911.3535.
- [39] C. Hopman, D. Guetta, E. Waxman, and S. Portegies Zwart, The Redshift Distribution of Short Gamma-Ray Bursts from Dynamically Formed Neutron Star Binaries, Astrophysical Journal 643 (2006), L91–L94.
- [40] N. Ivanova, C. O. Heinke, F. A. Rasio, K. Belczynski, and J. M. Fregeau, Formation and evolution of compact binaries in globular clusters - II. Binaries with neutron stars, MNRAS 386 (2008), 553–576.
- [41] V. Kalogera, K. Belczynski, C. Kim, R. O’Shaughnessy, and B. Willems, Formation of double compact objects, Physics Reports 442 (2007), 75–108.
- [42] V. Kalogera, C. Kim, D. R. Lorimer, M. Burgay, N. D’Amico, A. Possenti, R. N. Manchester, A. G. Lyne, B. C. Joshi, M. A. McLaughlin, M. Kramer, J. M. Sarkissian, and F. Camilo, The cosmic coalescence rates for double neutron star binaries, Astrophysical Journal 601 (2004), L179–L182, Erratum-ibid. 614 (2004) L137.
- [43] C. Kim, V. Kalogera, and D. R. Lorimer, The Probability Distribution of Binary Pulsar Coalescence Rates. I. Double Neutron Star Systems in the Galactic Field, Astrophysical Journal 584 (2003), 985–995.
- [44] , Effect of PSR J0737-3039 on the DNS Merger Rate and Implications for GW Detection, A life with stars, 2006, astro-ph/0608280.
- [45] R. K. Kopparapu, C. R. Hanna, V. Kalogera, R. O’Shaughnessy, G. Gonzalez, P. R. Brady, and S. Fairhurst, Host Galaxies Catalog Used in LIGO Searches for Compact Binary Coalescence Events, Astrophysical Journal 675 (2008), no. 2, 1459–1467.
- [46] W. H. Lee, E. Ramirez-Ruiz, and G. van de Ven, Short gamma-ray bursts from tidal capture and collisions of compact stars in globular clusters, ArXiv e-prints (2009), 0909.2884.
- [47] D. R. Lorimer, Binary and Millisecond Pulsars, Living Reviews in Relativity 11 (2008), 8–+.
- [48] R. N. Manchester, G. B. Hobbs, A. Teoh, and M. Hobbs, The Australia Telescope National Facility Pulsar Catalogue, Astronomical Journal 129 (2005), 1993–2006.
- [49] I. Mandel, Parameter estimation on gravitational waves from multiple coalescing binaries, ArXiv e-prints (2009), 0912.5531.
- [50] I. Mandel, D. A. Brown, J. R. Gair, and M. C. Miller, Rates and Characteristics of Intermediate Mass Ratio Inspirals Detectable by Advanced LIGO, Astrophysical Journal 681 (2008), 1431–1447.
- [51] J. M. Miller, A. C. Fabian, M. A. Nowak, and W. H. G. Lewin, Relativistic Iron Lines in Galactic Black Holes: Recent Results and Lines in the ASCA Archive, The Tenth Marcel Grossmann Meeting. On recent developments in theoretical and experimental general relativity, gravitation and relativistic field theories (M. Novello, S. Perez Bergliaffa, & R. Ruffini, ed.), January 2005, p. 1296.
- [52] M. C. Miller, Intermediate-mass black holes as LISA sources, Classical and Quantum Gravity 26 (2009), no. 9, 094031.
- [53] M. C. Miller and E. J. M. Colbert, Intermediate-Mass Black Holes, International Journal of Modern Physics D 13 (2004), 1–64.
- [54] M. C. Miller and V. M. Lauburg, Mergers of Stellar-Mass Black Holes in Nuclear Star Clusters, Astrophysical Journal 692 (2009), 917–923.
- [55] Ehud Nakar, Short-hard gamma-ray bursts, Phys. Rept. 442 (2007), 166–236.
- [56] R. Narayan, J. E. McClintock, and R. Shafee, Estimating the Spins of Stellar-Mass Black Holes by Fitting Their Continuum Spectra, Astrophysics of Compact Objects (Y.-F. Yuan, X.-D. Li, & D. Lai, ed.), American Institute of Physics Conference Series, vol. 968, January 2008, pp. 265–272.
- [57] R. Narayan, T. Piran, and A. Shemi, Neutron star and black hole binaries in the Galaxy, Astrophysical Journal 379 (1991), L17–L20.
- [58] G. Nelemans, Galactic Binaries as Sources of Gravitational Waves, The Astrophysics of Gravitational Wave Sources (J. M. Centrella, ed.), American Institute of Physics Conference Series, vol. 686, October 2003, pp. 263–272.
- [59] S. Nissanke, S. A. Hughes, D. E. Holz, N. Dalal, and J. L. Sievers, Exploring short gamma-ray bursts as gravitational-wave standard sirens, ArXiv e-prints (2009), 0904.1017.
- [60] R. M. O’Leary, B. Kocsis, and A. Loeb, Gravitational waves from scattering of stellar-mass black holes in galactic nuclei, MNRAS 395 (2009), 2127–2146.
- [61] R. M. O’Leary, R. O’Shaughnessy, and F. A. Rasio, Dynamical interactions and the black-hole merger rate of the Universe, Phys. Rev. D 76 (2007), no. 6, 061504.
- [62] J. A. Orosz, J. E. McClintock, R. Narayan, C. D. Bailyn, J. D. Hartman, L. Macri, J. Liu, W. Pietsch, R. A. Remillard, A. Shporer, and T. Mazeh, A 15.65-solar-mass black hole in an eclipsing binary in the nearby spiral galaxy M 33, Nature 449 (2007), 872–875.
- [63] R. O’Shaughnessy, V. Kalogera, and K. Belczynski, Binary compact object coalescence rates: The role of elliptical galaxies, ArXiv e-prints (2009), 0908.3635.
- [64] R. O’Shaughnessy, J. Kaplan, V. Kalogera, and K. Belczynski, Bounds on Expected Black Hole Spins in Inspiraling Binaries, Astrophysical Journal 632 (2005), 1035–1041.
- [65] R. O’Shaughnessy and C. Kim, Pulsar Binary Birthrates with Spin-Opening Angle Correlations, ArXiv e-prints (2009), 0908.2049.
- [66] R. O’Shaughnessy, C. Kim, T. Fragos, V. Kalogera, and K. Belczynski, Constraining Population Synthesis Models via the Binary Neutron Star Population, Astrophysical Journal 633 (2005), 1076–1084.
- [67] R. O’Shaughnessy, C. Kim, V. Kalogera, and K. Belczynski, Constraining Population Synthesis Models via Empirical Binary Compact Object Merger and Supernova Rates, Astrophysical Journal 672 (2008), 479–488.
- [68] R. O’Shaughnessy, R. Kopparapu, and K. Belczynski, Impact of star formation inhomogeneities on merger rates and interpretation of LIGO results, ArXiv e-prints (2008), 0812.0591.
- [69] Y. Pan, A. Buonanno, Y. Chen, and M. Vallisneri, Physical template family for gravitational waves from precessing binaries of spinning compact objects: Application to single-spin binaries, Phys. Rev. D 69 (2004), no. 10, 104017.
- [70] P. C. Peters, Gravitational Radiation and the Motion of Two Point Masses, Physical Review 136 (1964), 1224–1232.
- [71] E. S. Phinney, The rate of neutron star binary mergers in the universe: Minimal predictions for gravity wave detector, Astrophysical Journal 380 (1991), L17.
- [72] E. S. Phinney, Finding and Using Electromagnetic Counterparts of Gravitational Wave Sources, Astro2010: The Astronomy and Astrophysics Decadal Survey, Astronomy, vol. 2010, 2009, p. 235.
- [73] K. Postnov and L. Yungelson, The Evolution of Compact Binary Star Systems, Living Rev. Rel. 6 (2007).
- [74] K. A. Postnov and A. G. Kuranov, Neutron star spin-kick velocity correlation effect on binary neutron star coalescence rates and spin-orbit misalignment of the components, MNRAS 384 (2008), 1393–1398.
- [75] V. Raymond, M. V. van der Sluys, I. Mandel, V. Kalogera, C. Röver, and N. Christensen, Degeneracies in sky localization determination from a spinning coalescing binary through gravitational wave observations: a Markov-chain Monte Carlo analysis for two detectors, Classical and Quantum Gravity 26 (2009), no. 11, 114007.
- [76] R. A. Remillard and J. E. McClintock, X-Ray Properties of Black-Hole Binaries, ARA&A 44 (2006), 49–92.
- [77] C. Röver, R. Meyer, and N. Christensen, Bayesian inference on compact binary inspiral gravitational radiation signals in interferometric data, Classical and Quantum Gravity 23 (2006), 4895–4906.
- [78] F. D. Ryan, Gravitational waves from the inspiral of a compact object into a massive, axisymmetric body with arbitrary multipole moments, Phys. Rev. D52 (1995), 5707–5718.
- [79] A. Sadowski, K. Belczynski, T. Bulik, N. Ivanova, F. A. Rasio, and R. O’Shaughnessy, The Total Merger Rate of Compact Object Binaries in the Local Universe, Astrophysical Journal 676 (2008), 1162–1169.
- [80] A. S. Sengupta, for the LIGO Scientific Collaboration, and the Virgo collaboration, LIGO-Virgo searches for gravitational waves from coalescing binaries: a status update, ArXiv e-prints (2009), 0911.2738.
- [81] A. Sesana, J. Gair, I. Mandel, and A. Vecchio, Observing Gravitational Waves from the First Generation of Black Holes, Astrophysical Journal 698 (2009), L129–L132.
- [82] D. Sigg and the LIGO Scientific Collaboration, Status of the LIGO detectors, Classical and Quantum Gravity 25 (2008), no. 11, 114041.
- [83] J. R. Smith for the LIGO Scientific Collaboration, The path to the enhanced and advanced LIGO gravitational-wave detectors, ArXiv e-prints (2009), 0902.0381.
- [84] M. Stamatikos, N. Gehrels, F. Halzen, P. Mészáros, and P. W. A. Roming, Multi-Messenger Astronomy with GRBs: A White Paper for the Astro2010 Decadal Survey, Astro2010: The Astronomy and Astrophysics Decadal Survey, ArXiv Astrophysics e-prints, vol. 2010, 2009, p. 284.
- [85] S. E. Thorsett and Deepto Chakrabarty, Neutron star mass measurements. i. radio pulsars, Astrophys. J. 512 (1999), 288.
- [86] G. Török, M. A. Abramowicz, W. Kluźniak, and Z. Stuchlík, The orbital resonance model for twin peak kHz quasi periodic oscillations in microquasars, A&A 436 (2005), 1–8.
- [87] C. van den Broeck, D. A. Brown, T. Cokelaer, I. Harry, G. Jones, B. S. Sathyaprakash, H. Tagoshi, and H. Takahashi, Template banks to search for compact binaries with spinning components in gravitational wave data, Phys. Rev. D 80 (2009), no. 2, 024009.
- [88] M. van der Sluys, I. Mandel, V. Raymond, V. Kalogera, C. Röver, and N. Christensen, Parameter estimation for signals from compact binary inspirals injected into LIGO data, Classical and Quantum Gravity 26 (2009), no. 20, 204010.
- [89] M. van der Sluys, V. Raymond, I. Mandel, C. Röver, N. Christensen, V. Kalogera, R. Meyer, and A. Vecchio, Parameter estimation of spinning binary inspirals using Markov chain Monte Carlo, Classical and Quantum Gravity 25 (2008), no. 18, 184011.
- [90] M. V. van der Sluys, C. Röver, A. Stroeer, V. Raymond, I. Mandel, N. Christensen, V. Kalogera, R. Meyer, and A. Vecchio, Gravitational-Wave Astronomy with Inspiral Signals of Spinning Compact-Object Binaries, Astrophysical Journal 688 (2008), L61–L64.
- [91] J. Veitch and A. Vecchio, Assigning confidence to inspiral gravitational wave candidates with Bayesian model selection, Classical and Quantum Gravity 25 (2008), no. 18, 184010.
- [92] , Bayesian coherent analysis of in-spiral gravitational wave signals with a detector network, ArXiv e-prints (2009), 0911.3820.
- [93] S. J. Vigeland and S. A. Hughes, Spacetime and orbits of bumpy black holes, ArXiv e-prints (2009), 0911.1756.
- [94] F. J. Virgili, B. Zhang, P. O’Brien, and E. Troja, Are all short-hard gamma-ray bursts produced from mergers of compact stellar objects?, ArXiv e-prints (2009), 0909.1850.
- [95] R. Voss and T. M. Tauris, Galactic distribution of merging neutron stars and black holes - prospects for short gamma-ray burst progenitors and LIGO/VIRGO, MNRAS 342 (2003), 1169–1184.