Exploring intermediate and massive black-hole binaries with the Einstein Telescope
We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope (ET) to enhance our astrophysical knowledge through detections of gravitational waves emitted by binaries including intermediate-mass and massive black holes. The design target for such instruments calls for improved sensitivity at low frequencies, specifically in the –Hz range. This will allow the detection of gravitational waves generated in binary systems containing black holes of intermediate mass, –. We primarily discuss two different source types — mergers between two intermediate mass black holes (IMBHs) of comparable mass, and intermediate-mass-ratio inspirals (IMRIs) of smaller compact objects with mass – into IMBHs. IMBHs may form via two channels: (i) in dark matter halos at high redshift through direct collapse or the collapse of very massive metal-poor Population III stars, or (ii) via runaway stellar collisions in globular clusters. In this paper, we will discuss both formation channels, and both classes of merger in each case. We review existing rate estimates where these exist in the literature, and provide some new calculations for the approximate numbers of events that will be seen by a detector like the Einstein Telescope. These results indicate that the ET may see a few to a few thousand comparable-mass IMBH mergers and as many as several hundred IMRI events per year. These observations will significantly enhance our understanding of galactic black-hole growth, of the existence and properties of IMBHs and of the astrophysics of globular clusters. We finish our review with a discussion of some more speculative sources of gravitational waves for the ET, including hypermassive white dwarfs and eccentric stellar-mass compact-object binaries.
The Einstein Telescope (ET), a proposed third-generation ground-based gravitational-wave (GW) detector discussed in greater detail elsewhere in this volume, will be able to probe GWs in a frequency range reaching down to Hz Hild:2008 (); Freise:2009 (). This is lower than the limit of 40 Hz available to current ground-based interferometric GW detectors such as LIGO, Virgo, and GEO-600 or the Hz limit that could be reached by their second generation LIGO (); Virgo (); GEO600 (). On the other hand, GWs in the range above Hz will not be accessible to the planned Laser Interferometer Space Antenna (LISA, Bender:1998 ()), which will have sensitivity in the mHz–Hz range. The frequency range determines the typical masses of coalescing binaries that could be detected by an interferometer; for example, the frequency of GWs emitted from the innermost stable circular orbit of a test particle around a Schwarzschild black hole of mass is . The Einstein Telescope will therefore probe sources with total masses of hundreds or a few thousand which are less likely to be detected by LISA or the current ground-based detectors. This places the ET in a position to make complementary observations to LISA and LIGO/Virgo/GEO-600 and to carry out unique searches for several very exciting source types, particularly those involving light seeds of massive black holes and intermediate-mass black holes.
There is a significant body of evidence that massive black holes (MBHs) are generically found in the centers of massive galaxies Ferrarese2005 (). These MBHs merge during mergers of their host galaxies, and such mergers therefore trace the history of structure formation in the universe. Gravitational waves emitted during the mergers of MBHs with masses in the – will be detectable by LISA; dozens of detections could be made during the LISA mission Sesana:2004 (). According to some predictions, these massive black holes grow from light seeds of through accretion and mergers Madau2001 (); VHM (); SVH:2007 (). The typical frequencies of gravitational radiation emitted during the mergers of such systems will fall in the – Hz range, however, and will only be accessible to GW detectors sensitive in that range. The Einstein Telescope may be able to detect tens of such sources, determining their masses to an accuracy of a few percent and the luminosity distances to Sesana:2009ET (); Gair:2009ET ().
Additionally, globular clusters may host intermediate-mass black holes (IMBHs) with masses in the – range (see MillerColbert:2004 (); Miller:2009 () for reviews). Intermediate-mass-ratio inspirals (IMRIs) of neutron stars or stellar-mass black holes into these IMBHs could be detected by the second generation of ground-based detectors Mandel:2007rates (); however, the Einstein Telescope should be able to detect far greater numbers of events, up to as many as several hundred per year, and these events will range to higher IMBH masses. If the binary fraction in a globular cluster is sufficiently high or two globular clusters hosting IMBHs merge, an IMBH-IMBH binary can form and then coalesce, emitting gravitational waves Fregeau:2006 (). The Einstein Telescope could detect as many as thousands of such events, although, given the present uncertainty about the very existence of IMBHs, all such estimates must be viewed with a great deal of caution.
The Einstein Telescope may also be able to detect a number of other, more speculative, sources. These include the inspirals of stellar-mass black holes into IMBHs that may reside at the centers of dwarf galaxies, although we do not expect a significant rate of detectable signals of this type. There is also the possibility that ET will detect orbiting white dwarfs near the upper end of their allowed mass range, or eccentric compact object binaries.
This paper is organized as follows. In Section 2, we discuss the methodology for the event-rate calculations. We describe the adopted detector and network models, the formalism for estimating the signal-to-noise ratio and the waveform families used in the analysis. We then discuss in detail several types of GW sources of particular relevance to the ET. In Section 3, we consider sources involving IMBHs in globular clusters, both intermediate-mass-ratio inspirals and IMBH-IMBH coalescences. In Section 4, we focus on light massive black holes. We discuss how mergers between galaxies at high redshift can lead to IMBH-IMBH binaries detectable by ET, provided the black hole seeds in the galaxies are light. We also describe how light massive black holes in the centres of dwarf galaxies could act as IMRI sources. In Section 5, we discuss several speculative sources, including hypermassive white dwarfs and eccentric binaries. We finish with a discussion, in Section 6, of some of the potential scientific implications of ET observations of these sources. Section 7 provides a brief summary.
2 Methodology for event-rate calculations
2.1 The Einstein Telescope configuration
The design target for the Einstein telescope is a 10km scale interferometer, with a factor of increase in sensitivity over Advanced LIGO, and improved sensitivity at low frequencies. The ET design also calls for the ability to measure polarisation at a single site, which requires at least two non-coaligned coplanar detectors at the site. The currently favoured configuration is a triangular facility, with 10km long arms, and containing three independent detectors with 60 opening angles, as this has lower infrastructure costs and slightly better sensitivity than two right angle detectors placed at 45 to one another Freise:2009 (). We refer to this triangular design as a “single ET”. In Figure 1 we show the target ET noise curve, labelled “ET baseline” Hild:2008 (). This noise curve is for a single right-angle interferometer with the ET design sensitivity. Unless otherwise stated, signal-to-noise ratios (SNRs) etc. will be quoted for this configuration. The sensitivities of one interferometer, two right-angle interferometers and a single ET are changed relative to this by factors of , and respectively. The “ET baseline” design has recently been superseded by the curve labelled “ET B” in Figure 1, but we have checked that this change does not significantly affect our results, since the noise curves are largely similar in the frequency range where massive systems accumulate most of their SNR. Figure 1 also shows an alternative ‘xylophone’ configuration for ET that was described in Hild:2009 (). Such a noise curve is realised by operating two detectors within the same vacuum system, one optimised for low-frequency sensitivity and the second for high frequency sensitivity. The net effect on the composite noise curve is an improved sensitivity near Hz. The ET design and noise curve has not yet been finalised, nor what fraction of time ET would spend in high-frequency and low-frequency operation under the xylophone configuration. We will see in Figure 4 in Section 4.1 that the xylophone mode is to be preferred for the detection of black hole binaries in the – range.
In Section 4.1 we will quote results for the parameter estimation accuracies that are achievable for mergers of light seeds of massive black holes detected by ET. As the events are short lived, parameter estimation requires the existence of a network containing multiple detectors. The results we will present are taken from Sesana:2009ET (); Gair:2009ET (), in which four third-generation network configurations were considered — (i) one ET at the geographic location of Virgo, plus a second right-angle 10km detector at the location of LIGO Hanford or Perth (Australia); (ii) as configuration (i) plus a third 10km detector at the location of LIGO Livingston; (iii) as configuration (i) but with the Hanford/Perth 10km detector replaced by a second ET; and (iv) three ETs, one at each of the three sites.
2.2 Signal-to-noise ratio
In Sections 3.1 and 3.2 we will quote signal-to-noise ratios for events detected with ET. In this section, we describe how these were calculated. The signal-to-noise ratio (SNR) for a waveform measured by a single detector with one-sided noise power spectral density is given by . Here is the noise-weighted inner product
where is the Fourier transform of the waveform and denotes the complex conjugate.
We define the horizon distance as the distance at which an optimally oriented, overhead source produces an SNR of . The actual gravitational-wave emission is not isotropic and nor is the response of the detector. We can define the average range as the radius of a sphere whose volume is equal to the true (non-spherical) volume in which inspiral sources can be detected with an SNR greater than . For uniformly distributed, randomly oriented sources, and without applying corrections for cosmological redshift, this range is related to the horizon distance by FinnChernoff:1993 (). In general, the relationship between and for sources at cosmological distances will depend on the cosmology, the redshift dependence of the source distribution and even on the intrinsic luminosity of the sources. In this paper, we will ignore these issues and use the usual correction factor, since a more rigorous calculation is not currently available. We expect that the error that arises from this approximation will be of the same order as the SNR uncertainties that arise from inaccuracies in the waveform model, which we will discuss in the next section, but emphasise that this must be verified by a proper calculation in the future.
For multiple detectors, the inner product in equation (1) should be replaced by sums of inner products over individual detectors. Thus, the existence of identical interferometers will increase the range by a factor of over the single-interferometer range if the network SNR threshold is fixed.
In this paper, we will focus on two different types of source — binaries consisting of two intermediate mass black holes of comparable mass; and intermediate-mass-ratio inspirals (IMRIs) of stellar mass compact objects (neutron stars or black holes) into IMBHs. We model the gravitational waves emitted by these systems in different ways. For comparable mass binary systems at the upper end of the detectable mass range, a significant amount of energy is radiated during the merger and ringdown phases and so it is important to include these in waveform models for signal-to-noise ratio calculation and parameter estimation. The recent advances in numerical relativity have allowed the construction of hybrid waveform models that include inspiral, merger and ringdown in a self-consistent way in a single template. There are two models currently available, both of which are for systems containing non-spinning black holes. The assumption of zero spin and the intrinsic waveform uncertainties lead to uncertainties in the SNRs at the level of a few tens of percent, which we will discuss later. However, the corresponding uncertainties in the event rates are much smaller than the typical uncertainties that arise from the astrophysics, since the lower limit on the number of intermediate-mass black holes in the Universe is zero.
The non-spinning “phenomenological” inspiral, merger, ringdown model (NSphenom) Ajith:2008 (); Ajith:2008rev (), is constructed by taking a simple ansatz for the waveform model, inspired by post-Newtonian theory, and fitting the coefficients in this model to numerical relativity simulations. The phenomenological model has recently been extended to spinning black holes with non-precessing spins Ajith:2009spin (). The “effective-one-body, numerical relativity” (EOBNR) model Buonanno:2007EOBNR () uses post-Newtonian expressions to model the inspiral radiation, which are matched onto fits to numerical relativity simulations for the merger radiation and then onto analytic expressions for the quasinormal mode ringdown radiation. The EOBNR waveforms are constructed to match perturbative results in the test-particle limit where the mass-ratio tends to zero. Denoting the mass of the most massive object in the binary by and the less massive object by , the NSphenom and EOBNR models give waveforms as a function of the total mass of the binary, , the reduced mass ratio, and the time of merger, . Including the detector response introduces six additional extrinsic parameters expressing the relative location and orientation of the source and detector: the distance, the two sky-location angles, two binary orientation angles, and the phase at some fiducial time, e.g., .
In Section 4.1 we quote results for the SNRs of comparable mass IMBH mergers detected by ET. These SNRs were computed in Sesana:2009ET () using the NSphenom model. ET SNRs for the same systems were recomputed in Gair:2009ET () using both models and the results were found to agree to , with the NSphenom predictions being higher for low mass ratios, , and the EOBNR predictions being higher for high mass ratios, . The two waveform families are constructed in different ways and have been matched to different numerical relativity simulations, so it is not surprising that the results they predict differ. The level of difference provides a guide to how much the SNRs computed using either waveform family will differ from the true SNRs of these systems that nature provides. This SNR uncertainty from the waveform model must be compared to other uncertainties, such as the omission of spin from the waveform model and uncertainties in the intrinsic astrophysical rates. The effect of spin on the average SNR will also be at the level Mandel:2007spin (), although it is more significant for the highest mass systems as it can bring otherwise undetectable sources into the frequency band of the detector. However, the net effect will be to increase the rate of detectable events so using non-spinning models can be thought of as conservative in that regard. A SNR uncertainty leads to an uncertainty in the event rate of a factor of , but typical uncertainties in the astrophysical rates are an order of magnitude. The SNR could be out by as much as a factor of before it would be comparable to these astrophysical uncertainties, so we consider the waveform uncertainties as fairly negligible. Recent comparisons with additional numerical data suggest that the errors in the EOBNR waveforms are significantly smaller than those in the NSphenom waveforms, at least for equal-mass sources Buonanno:2009EOBNR (). For this reason, we use the EOBNR model to compute the new results presented in this paper, specifically the SNRs for IMBH-IMBH mergers in globular clusters, which we will describe in Section 3.2.
The NSphenom and EOBNR waveforms have been matched to numerical relativity simulations, but only for mass ratios of 1:4 and higher. Computational requirements suggest it is unlikely that numerical simulations using current techniques will go beyond mass ratios of 1:10 in the near future, although this may be possible using innovative new approaches. Post-Newtonian theory also breaks down once the mass-ratio becomes too extreme. For very extreme mass ratios, –, gravitational waveforms can be computed using black hole perturbation theory Poisson:2004 (), in which the smaller object is regarded as a perturbing field of the background spacetime of the larger object and radiation reaction is described in terms of the ‘self-force’. Significant progress has been made over the past few years in self-force calculations, which has led to the calculation of the self-force for circular orbits in the Schwarzschild spacetime BarackSago:2007 (), including the shift in the location of the innermost-stable-circular-orbit (ISCO) that results from the action of this force BarackSago:2009 (). However, even at mass ratios of , the terms that are missing in the first-order self-force formalism are estimated to have a marginal effect on the phasing of waveforms for LISA sources Huerta:2009 (). The mass ratios for typical IMRI sources for ET lie somewhere between these extremes, being typically –. In this regime neither post-Newtonian nor perturbative waveforms will be adequate on their own to model the true waveforms MandelGair:2009 (). More research is needed to devise waveforms that are suitable for bridging this gap. In the meantime we must make do with the available waveforms, with the understanding that these are not completely accurate.
In Section 3.1, we will present ET event rates for IMRI sources. These were estimated using SNRs computed with EOBNR waveforms, since these should be accurate both for and in the limit . To check the validity of the results, we also computed SNRs using the NSphenom waveforms and using perturbative waveforms for circular and equatorial inspirals, as described and tabulated in FinnThorne:2000 (), and as used to estimate SNRs for LISA in Gair:2009 (). The latter waveform model includes only the inspiral phase, and so we compared those results to the inspiral contribution to the EOBNR and NSphenom SNRs. For the early portion of the inspiral, significantly before the last stable orbit, we found that the SNRs obtained with the three waveform families agreed to within ; this difference is of a comparable magnitude to the effect of omitting relativistic corrections in post-Newtonian inspirals. We find, however, that for more massive systems, when merger and ringdown contribute a significant portion of the SNR, estimates from the NSphenom and EOBNR waveforms differ very significantly, with NSphenom waveforms predicting SNRs that are greater by more than an order of magnitude. In Figure 2 we show how the SNRs computed using the two models vary as a function of for a fixed luminosity distance and a fixed redshifted total mass of . This figure clearly shows the large difference between the two models for small . Our theoretical expectation is that the energy emitted during the ringdown should scale as , so that the ringdown SNR should scale as for a fixed total mass. However, the SNR predicted by the ringdown portion of the NSphenom waveform scales roughly as , and is therefore significantly over-predicted for small . For this reason we used EOBNR waveforms, which exhibit the correct scaling, for estimating IMRI SNRs, with the understanding that there is a clear need for more careful and accurate modelling of IMRI radiation in the future. This is essential not least because detection of these systems will almost certainly rely on matched filtering, for which accurate template waveforms are needed.
3 IMBH sources in Globular Clusters
A particularly exciting possibility for GW astronomy is the observations of intermediate-mass black holes (IMBHs) in globular clusters. Unlike stellar-mass and supermassive black holes, IMBHs have not been identified dynamically, that is, through the observation of one or more orbiting objects and the application of Kepler’s law. Thus, searches for their existence have relied on circumstantial evidence (for reviews and recent results see MillerColbert:2004 (); Trenti:2006 (); Noyola:2008 (); MaccaroneServillat:2008 (); Farrell:2009 ()). This evidence includes observed fluxes of radiation that, if isotropic and sub-Eddington, would imply IMBH masses above in many cases; thermal peaks at a few tenths of a keV, which, with normal scalings, would imply masses of tens to thousands of solar masses; and variability features that, if identified with corresponding features observed in systems containing stellar-mass black holes, would also suggest IMBHs weighing in at hundreds of solar masses. Depending on the threshold one adopts for IMBH likelihood, candidates exist in one out of every tens to hundreds of galaxies. Of course, the actual number of IMBHs is likely to be far greater than this, in the same way that, although we have only about 25 strong candidates for stellar-mass black holes in the Galaxy, the total number is likely to be . The lack of definitive dynamical evidence for IMBHs means, however, that their existence is still inconclusive and alternate explanations for the observations have been proposed (see Hurley:2007 (); Colbert:2008 () for recent discussions).
One channel leading to IMBH formation is via the runaway collisions of massive stars on time scales too short to allow for stellar evolution, Myr PZwart:2004 (). Recent simulations of runaway collisions with the inclusion of stellar winds suggest that winds will prevent the growth of IMBHs in all but the most metal-poor environments Glebbeek:2009 (), although these simulations extrapolate wind rates from much less massive stars and the collision rates are so high that the collision products are likely to be extended bags of gas rather than relaxed stars. Alternatively, IMBHs could form through mergers of stellar-mass black holes in dense subclusters at the cores of globular clusters; however, recoil kicks may eject the products of such collisions from the host globulars OLeary:2006 (). In addition to runaway collisions and repeated stellar-mass BH collisions, other IMBH formation channels include IMBH growth through gas accretion early in the cluster history PZMcMillanGieles:2010 () and direct collapse of Population III stars (these IMBHs could also be seeds for massive black holes, which we discuss in Section 4.1). We refer the reader to MillerColbert:2004 (); Miller:2009 () for reviews.
Numerical simulations of globular clusters suggest that IMBHs could merge with numerous lower-mass compact objects (COs) during the lifetime of the cluster Taniguchi:2000 (); MillerHamilton:2002a (); MillerHamilton:2002b (); MouriTaniguchi:2002a (); MouriTaniguchi:2002b (); Gultekin:2004 (); Gultekin:2006 (); OLeary:2006 (); Mandel:2007rates (), through a combination of the emission of gravitational radiation, binary exchange processes, and secular evolution of hierarchical triple systems. For IMBH mass , the GWs generated during the inspiral of a stellar-mass object (black hole or neutron star, since a white dwarf or a main sequence star would be tidally disrupted) into an IMBH are potentially detectable by the Einstein Telescope. Ringdown radiation could even be detected from more massive IMBHs.
When the primordial binary fraction in a globular cluster is sufficiently high, , stellar collisions during binary scattering interactions may lead to the production of two IMBHs in a single cluster, according to Monte Carlo simulations carried out by Gurkan:2006 (). Since observations and numerical calculations suggest that clusters may be born with large binary fractions (e.g., Ivanova:2005 ()), the formation of two IMBHs may be generic in sufficiently dense and massive clusters. If this happens, the two IMBHs will exchange into a common binary which, after shrinking via dynamical friction and dynamical encounters with other stars, will merge through radiation reaction; all of these processes occur on a timescale of Myr Fregeau:2006 ().
Additionally, two IMBHs from different globular clusters may merge during the merger of their parent clusters Amaro:2006imbh (). N-body simulations suggest that, as the clusters merge, the IMBHs should form a binary with a peak eccentricity of Amaro:2006imbh (), although the residual eccentricity would be negligible by the time the frequency reaches Hz and the system enters the ET band. The binary should merge on a timescale of a few hundred million years through a combination of dynamical interactions with stars and gravitational-wave emission. If the probability of forming an IMBH in a cluster is high, and if clusters merge with a probability – as discussed in Amaro:2006imbh (), the rates of mergers of IMBHs originating in different host clusters could be competitive with the single cluster channel (see also pauIMBH ()).
In the following two subsections, we estimate the rates with which the Einstein Telescope could detect gravitational waves from these two globular cluster channels. In Section 3.1 we estimate the rates of detectable inspirals of stellar-mass compact-objects into IMBHs (IMRIs) and in Section 3.2 we estimate the rate of detectable coalescences of IMBH-IMBH binaries formed within a single globular cluster. These subsections contain detailed derivations which are presented here for the first time; readers who are not interested in the details of these calculations will find the results for IMRIs in Table 1 and the results for single-cluster IMBH-IMBH binary mergers in Equation (15).
3.1 Intermediate-mass-ratio inspirals into IMBHs
In an earlier work, a subset of the authors analyzed the possibility of detecting intermediate-mass-ratio inspirals (IMRIs) of compact objects into IMBHs with Advanced LIGO Mandel:2007rates (). In that case, it was found that binary tightening via three-body interactions was the dominant channel that led to IMRIs. The IMBH, as the most massive object in the cluster, readily switches into compact-object (CO) binaries. Once a sufficiently hard CO-IMBH binary is formed, the binary will be hardened rather than disrupted by three-body interactions with other stars in the cluster. Eventually, as the interacting stars take away energy from the binary, the binary will tighten to the point where radiation reaction from gravitational-wave emission becomes dominant and drives the binary to coalesce. For COs that are neutron-stars or black-holes, it is possible to compute the distance to which the gravitational waves can be detected and convert this into an estimated detection rate.
Here, we repeat that calculation for ET sources, with the following two major changes. First, we take advantage of the recent development of hybrid waveforms that describe all three phases of the coalescence – inspiral, merger, and ringdown – to compute the SNR from the full GW signal, rather than just the inspiral portion. We use the EOBNR waveforms Buonanno:2007EOBNR () for this calculation. As discussed in Section 2.3, these waveforms have not been tested for mass ratios below , but, unlike the NSphenom waveforms, they do appear to behave correctly in the extreme-mass-ratio limit, . We use these waveforms as there are no better IMRI models available at the present time, but emphasise that there will be some corresponding uncertainty in the results, at the level of tens of percent. Second, because ET has a low frequency cutoff ( Hz) that is lower than Advanced LIGO ( Hz), we consider inspirals into IMBHs along with inspirals into IMBHs. We note, however, that for higher IMBH masses the IMBH could dominate the dynamics in the center of the cluster and a cusp could be formed around the IMBH, possibly increasing the importance of the direct-capture scenario Hopman:private (); additional discussion of this possibility can be found in Section 2.3 of Mandel:2007rates ().
Approximately close interactions with stars of mass are required to reduce the semimajor axis of a CO-IMBH binary with IMBH mass by one -folding Quinlan:1996 (). Stars come within a distance equal to the semimajor axis separation, , of the binary at a rate
where the bracketed expression is the gravitationally focused cross-section, is the velocity dispersion, and is the number density of stars in a globular cluster, with fiducial values for core-collapsed globulars taken from PryorMeylan:1993 (). The last -folding time dominates the hardening rate, so the hardening time-scale is
Meanwhile, the gravitational-wave merger timescale for a binary of semimajor axis , eccentricity , reduced mass approximately equal to the CO mass , and total mass is Peters:1964 ()
where in the last equality we set as the eccentricity after the final three-body encounter, following Gultekin:2006 (). Minimizing the total merger time over , while setting , and to their fiducial values, allows us to compute the CO-IMBH coalescence rate per globular cluster, .
To compute the volume within which the Einstein Telescope can detect such IMRIs, we follow the procedure outlined in Section 2.2. We use EOBNR waveforms and ignore the spin of the IMBH, which we expect to be small, , after a significant number of minor mergers Mandel:2007spin (). We compute the range for a “single ET” configuration. The range is a function of the redshifted masses of the IMBH, , and the compact object, . After computing the range, we convert it into a redshift, , by inverting the following expression for the luminosity distance Hogg:1999 ():
Here, we implicitly assume a flat universe () and use , , km s Mpc; and Mpc. We assume that the typical source is located near the redshift that corresponds to the search range, and obtain the source-frame masses by dividing the redshifted masses by ; we use these source-frame masses to compute the merger timescale from Equations (3) and (4).
We additionally assume that of clusters form an IMBH and are sufficiently dense to be relevant to the rate calculation, and that globular clusters have a fixed comoving space density of PZwart:2000 (). For , this yields a density of for relevant clusters. We compute the comoving volume up to redshift by integrating the following expression for Hogg:1999 (), with the cosmological parameters defined above:
The rate of detectable events can then be estimated as , where the factor of is included to convert the coalescence time measured in the source frame to time measured in the observer frame.
Table 1 summarizes the rate predictions for four combinations of and . Although the lack of knowledge about IMBHs and their mass distributions makes it impossible to generate firm predictions, and even a lower limit of zero IMRIs is possible, it appears that ET may detect hundreds of compact-object IMRIs into IMBHs over three years of operation. If ET is operated in the xylophone configuration these rates would increase further.
3.2 IMBH-IMBH inspirals
In this section, we wish to estimate the rate at which the single cluster channel generates IMBH-IMBH binaries that are detectable with the Einstein Telescope. To do this, we follow the event rate calculation for LISA and Advanced LIGO described in Fregeau:2006 (). Once a pair of IMBHs is formed in a single cluster, they sink rapidly to the center where they form a binary and merge via three-body interactions with the stars in the cluster (see Fregeau:2006 (); Amaro:2010 () for more details). Therefore, the rate of IMBH binary mergers is just the rate at which pairs of IMBHs form in clusters. The rate of detectable coalescences is
Here is the time measured in our observer’s frame and is the time measured at the redshift of the merger; is the total mass of the coalescing IMBH-IMBH binary and is the mass ratio between the IMBHs; is the maximum redshift to which the ET could detect a merger between two IMBHs of total mass and mass ratio ; is the relation between local time and our observed time, and is the change of comoving volume with redshift, given by Eq. (6).
We make the following assumptions:
IMBH pairs form in a fraction of all globular clusters.
We neglect the delay between cluster formation and IMBH coalescence, since it is expected to be no more than a few tens of millions of years Fregeau:2006 ().
When an IMBH pair forms in a cluster, its total mass is a fixed fraction of the cluster mass, . This assumption is based on what is typically seen in simulations Gurkan:2004 (). As there are no current constraints on the mass ratio, we take it to be uniformly distributed between and . We restrict our attention to systems with a total mass in the IMBH range, which we adopt to be between and . This means we confine our attention to clusters with masses (note that the lower limit is different from that chosen in Fregeau:2006 () since here we set for IMBH sources). Thus,
The distribution of cluster masses scales as independently of redshift ZhangFall:1999 () and the total mass formed in all clusters in this mass range at a given redshift is a redshift-independent fraction of the total star formation rate per comoving volume:
which provides the normalization for :
Under these asssumptions, Eq. (7) predicts the rate of detectable coalescences per year as
Note that here is measured in solar masses and is measured in Mpc.
Rather than computing for all values of and , we rely on the following fitting formula for the average range as a function of the redshifted total mass , obtained by using EOBNR waveforms to model the coalescence (see Section 2.2):
where , for and , for .
Lensing of gravitational wave sources adds some uncertainty to this picture. Individual sources can be magnified or de-magnified by lensing, making them visible at a greater or lesser distance than predicted by the preceding formula. Flux conservation ensures that the expected magnification of a source is , meaning no net change in flux (this argument was first elucidated in weinberg76 ()), which leads us to expect that the total change in the event rate will be small. The magnification distribution that arises from lensing peaks at less than , i.e., a demagnification, and shows an exponential fall-off for large magnifications peacock82 (); wang02 (). The peak moves toward greater demagnification for sources at higher redshifts, but this is compensated by a longer tail toward very high magnifications. The amount of volume added to a flux-limited sample by these highly magnified lines of sight compensates for the smaller volume lost by each of the (greater number) of demagnified lines of sight. In peacock82 (), it was shown that weak lensing did not significantly change the number counts in a flux limited radio sample, but the number of events tended to be increased by the strong lensing tail. The total change in number counts was only a few percent. We can therefore ignore the effect of weak lensing on the number counts of gravitational wave sources, although it will have an impact on the precision to which distances can be measured with gravitational wave observations.
Ignoring lensing, we can compute by inverting Eq. (5). For a given choice of and , the maximum detectable redshift is then obtained by finding a self-consistent solution of
The integrals over and in Eq. (3.2) were evaluated for two specific values of . For , the total rate was found to be ; for , it was . The range varies smoothly with and so we can estimate the integral over to be the average of these two rates. This yields a final estimate of the total rate as
where we arbitrarily chose and as the reference values for these unknown parameters.
4 Sources in low-mass galaxies
4.1 Light seeds of MBHs at high redshifts
Supermassive black holes (SMBHs) weighing millions to billions of solar masses are nowadays believed to reside in most local galaxies (Ferrarese2005 () and references therein). The masses of today’s SMBHs exhibit clear correlations with the properties of their host galaxies (luminosity, mass, and stellar velocity dispersion), suggesting there is a single mechanism for assembling SMBHs and forming galaxies. The evidence therefore favours a co-evolution between galaxies and SMBHs.
In the currently favoured cold dark matter cosmology, galaxies today are expected to have been built up, via a series of mergers, from small-mass building blocks that condensed out at early cosmic times. A single big galaxy can be traced back to hundreds of smaller components with individual masses as low as M. Similarly, we expect the SMBHs found in galaxies today to have grown partially by accretion and partially by mergers following mergers between galaxies (e.g., VHM (); Malbon2007 ()), so that a single SMBH can be traced back to some number of ‘seed’ black holes at early times VLN2008 (). There are large uncertainties in this picture, however. Did seed black holes form efficiently in small galaxies (with shallow potential wells) at early times, or was their formation delayed until substantial galaxies with deeper potential wells had been formed? This is a key question, as the mass and the occupation number of the seeds ultimately dictates the occupation number of SMBHs in galactic centers.
The formation of SMBHs is far less well understood than that of their light stellar-mass counterparts. The ‘flow chart’ presented by Rees1978 () still stands as a guideline for the possible paths leading to the formation of SMBH seeds in galactic nuclei. One possibility is that the seeds of SMBHs were the remnants of the first generation of stars, formed out of zero-metallicity gas Madau2001 (). In a cold dark matter universe, structure builds up hierarchically, so the smaller clumps at the earliest cosmic times have shallower potential wells. Stars cannot form until the clumps are sufficiently big to provide a potential well deep enough to pull in gas that can cool radiatively and contract to make a protostar. This requires dark matter clumps – minihalos – of M at redshifts of . The first stars forming in these minihalos develop under very different conditions from present-day stars: there are no heavy elements (so that molecular hydrogen is the only effective coolant), no dust, and no magnetic fields. These conditions mean that these ‘Population III’ stars were likely very massive, having characteristic masses of the order of M (e.g., Bromm1999 (); Nakamura2001 (); Abel:2002 (); Yoshida2006 ()). This prediction relies on the absence of efficient cooling agents in the primordial metal–free gas. If Population III stars form with masses 40 M140 M or M M, they are predicted to collapse and form IMBHs directly with little mass loss Fryer2001 (), i.e., leaving behind seed IMBHs with masses . This is a plausible formation mechanism for the seeds upon which supermassive black holes are grown VHM (), although more massive black holes may have been formed after the epoch of the first stars in dark-matter halos with virial temperatures of K Bromm2003 (); Spaans2006 (); BVR2006 () via ‘direct collapse’, as described in the following.
Direct collapse models for MBH formation rely on the collapse of supermassive objects formed directly out of dense gas haehnelt1993 (); LoebRasio1994 (); BrommLoeb2003 (); Koushiappas2004 (); BVR2006 (); LN2006 (). The physical conditions (density, gas content) in the inner regions of mainly gaseous proto-galaxies make these loci natural candidates, because the very first proto-galaxies were, by definition, metal-free, or at the very least very metal-poor. Enriched halos have a more efficient cooling, which in turn favours fragmentation and star formation over the efficient collection of gas conducive to MBH formation. In a typical galaxy, however, the tidally induced angular momentum would still be enough to provide centrifugal support at a distance pc from the centre, and halt collapse, ultimately leading to the formation of a disk. Additional mechanisms inducing transport of angular momentum are needed to further condense the gas until conditions fostering MBH formation are achieved. An appealing route to efficient angular momentum shedding is by global dynamical instabilities, such as the “bars-within-bars” mechanism, that relies on global gravitational instability and dynamical infall Shlosman1989 (); BVR2006 (). Self-gravitating gas clouds become bar-unstable when the level of rotational support surpasses a certain threshold. A bar can transport angular momentum outward on a dynamical timescale via gravitational and hydrodynamical torques, allowing the radius to shrink. Provided that the gas is able to cool, this shrinkage leads to even greater instability, on shorter timescales, and the process cascades. This mechanism is a very attractive candidate for collecting gas in the centres of halos, because it works on a dynamical time and can operate over many decades of radius. It has also been proposed that gas accumulation in the central regions of protogalaxies can be described by local, rather than global, instabilities. During the assembly of a galaxy disc, the disc can become self-gravitating. As soon as the disc becomes massive enough to be marginally stable, it will develop structures that will redistribute angular momentum and mass through the disc, preventing the surface density from becoming too large and the disc from becoming too unstable. To evaluate the stability of the disc, the Toomre stability parameter formalism can be used Koushiappas2004 (); LN2006 (). The gas made available in the central compact region can then form a central massive object, for instance via the intermediate stage of a ‘supermassive’ star hoyle1963 (); Baumgarte1999 (), or a ‘quasistar’, an initially low-mass black hole rapidly accreting within a massive, radiation-pressure-supported envelope, BVR2006 (); Begelman2008 (). In both cases, the mass function of seeds is predicted to peak at VLN2008 ().
As described in Section 3.1, the formation of an IMBH as a result of dynamical interactions in dense stellar systems is a long-standing idea, which could also create intermediate mass MBH seeds. This process could have been very effective in the very first stellar clusters that formed in high-redshift proto-galaxies, when the Universe was not as metal-rich as now. Low metallicity favors the growth of a very massive star, the precursor of an IMBH remnant. The mass loss due to winds is significantly reduced in metal-poor stars, which greatly helps in increasing the mass of the final IMBH remnant (cf. yungelson2008 ()). The formation of stellar clusters and the possible evolution of the stellar systems up to IMBH formation are explored in Devecchi2009 (). Figure 3 shows three mass functions for three different MBH ‘seed’ scenarios: direct collapse BVR2006 (), runaway stellar mergers in high-redshift clusters, and Population III remnants VHM ().
It is uncertain how many MBH ‘seeds’ formed, and in which mass range. Equally uncertain is how these ‘seed’ black holes grew within their host minihalos. It is not obvious if efficient accretion onto these seeds could have taken place, at least early on, in the fragile environment that the shallow potential wells of minihalos represent Milos2009 (). It is likely that seed IMBHs can grow efficiently only if they are hosted in the most massive galaxies at these early cosmic epochs, while IMBHs in an ‘average’ galaxy could have experienced intermittent and inefficient accretion, thus leaving behind a population of underfed IMBHs with a mass range similar to that of the original seeds, .
The Einstein Telescope will be able to probe mergers between black-hole seeds at high redshift, and thus help to distinguish between these various channels for seed formation. An estimate for the ET event rate under the Pop III model can be computed using Monte-Carlo merger-tree realizations based on the extended Press-Schechter formalism PressSchechter:1974 (), as described in VHM (); Volonteri:2007 (). This was done in Sesana:2009ET (), and the results which we now quote are all taken from that work and the companion paper Gair:2009ET (). These papers considered four different models that were based on the same merger tree realisations (taken from VHM (); Volonteri:2007 ()), but differed in the initial mass distribution of seeds and in the prescription for accretion onto the seed black holes. In these scenarios, which were all based on having the light, Pop III, remnants described earlier in this section as the seeds for black hole formation, a single ET would detect – seed mergers, depending on the model. The detected mergers would be between black holes with total mass ranging from up to , where is the mass of the lightest seed black hole in the initial mass distribution. This minimum seed mass is rather uncertain and depends on the details of the model used, as discussed earlier. In the scenarios considered in Sesana:2009ET (), was either or . The detected events would be seen at redshifts –, although this could extend to for the lightest seed model, which had . If ET was operated in the xylophone configuration described in Section 2.1, the number of events seen would be increased to several tens, and these would be out to a redshift Gair:2009ET ().
Figure 4, reproduced from the data in Gair:2009ET (), shows how the number of events seen by ET over three years varies as a function of the signal-to-noise ratio required in a single 10km right-angle interferometer for detection. The SNR required in the network of detectors is likely to be , although this depends somewhat on data-analysis issues, and on the amount of source confusion in the data stream. A network SNR of corresponds to an SNR in the single right-angle detector of for a single ET, or SNRs of , , and for the network configurations (i) – (iv) described in Section 2.1. In Figure 4 we show results for two of the four light-seed models considered in Sesana:2009ET (), and for both the baseline and xylophone configurations of the detector. The rate for the baseline ET configuration is rather sensitive to the SNR that is ultimately required for a confident detection, but the xylophone configuration is more robust, as it has improved sensitivity at just the right frequency for systems with mass in the – range. The mergers seen by ET will be complementary to mergers between heavier black holes that will be seen by space-based detectors such as LISA, ALIA or DECIGO Gair:2009ET (). The combination of detectors will provide a nearly complete survey of mergers between galactic black holes, yielding important constraints on astrophysical models of galaxy formation and growth. The utility of observations of binary black hole systems with multiple detectors was also discussed in pauIMBH (), with specific reference to IMBH-IMBH binary mergers arising from the mergers of globular clusters containing central IMBHs. They found that, if LISA and ET were operating concurrently, the same IMBH binary could be detected by both detectors with a time separation of a few months. Even if observations of the same system are not made with different detectors, each detector will provide a measurement of the rate of black hole mergers in a different black hole mass range. These rate observations over the black hole mass spectrum will provide important constraints on models of black hole growth.
One important question is whether ET will be able to distinguish between black-hole mergers coming from this channel, and those described in Section 3.2 that arise in globular clusters. To provide constraints on merger histories, it is necessary to know that an observed event is associated with a galaxy merger. The masses and redshifts of events will provide some information, but more work is required to understand what observational signatures provide the best discriminating power. We would expect mergers between seed black holes to occur over a range of redshifts, with some events at redshifts . In the mechanism described in Section 3.2, the black-hole binaries form and merge very quickly, so this could also produce events over a range of redshifts. However, the distinction between these two formation channels becomes increasingly vague at high redshift, when galaxies are in the process of formation. What is important for the light-seed scenario is that black holes of low mass, , exist at high redshift. Therefore, being able to identify an event as being between two black holes at redshift would be an important constraint, regardless of how that seed had initially formed. A single ET cannot measure the six extrinsic parameters of a merger source on its own — at least one additional non-colocated detector will be required. Possible network configurations were discussed in Section 2.1. With one additional km detector at the location of LIGO Hanford, the ET network will be able to determine the luminosity distance of a source to an accuracy of . Adding a third km detector at the site of LIGO Livingston or upgrading the detectors to ETs improves this modestly to Sesana:2009ET (); Gair:2009ET (). There will also be an additional distance error due to weak lensing of the signal, but this will be considerably smaller (, see, for example, shapiro09 ()) than the intrinsic error from the gravitational wave observations. If we assume that the luminosity distance is converted into a redshift using the concordance cosmology at that time, the redshift error will be comparable to the distance error. Thus, an ET network should be able to say with confidence if an event is indeed occurring between two black holes at high redshift, .
4.2 MBHs in dwarf galaxies
There are two simple arguments that lead us to believe that – black holes might inhabit the nuclei of dwarf galaxies today. Firstly, the mass of SMBHs detected in neighboring galaxies scales with the bulge mass — or stellar velocity dispersion () — of their host galaxy Magorrian1998 (); Gebhardt2000 (); Ferrarese2000 (); Haring2004 (). The lowest-mass galaxies currently known have velocity dispersions – km s Walker2009 (). If we extrapolate the correlation to these values, we expect the putative IMBHs to have masses in the range of hundreds to thousands of solar masses.
Secondly, as SMBHs grow from lower-mass seeds, it is natural to expect that a leftover population of progenitor IMBHs should also exist in the present universe. Indeed, one of the best diagnostics of ‘seed’ formation mechanisms would be to measure the masses of IMBHs in dwarf galaxies. This can be understood in terms of the cosmological bias. The progenitors of massive galaxies have a high probability that the central SMBH is not “pristine”, that is, it has increased its mass by accretion, or it has experienced mergers and dynamical interactions. Any dependence of on the initial seed mass is largely erased. However, low-mass galaxies undergo a quieter merger history, and as a result, at low masses the BH occupation fraction and the distribution of BH masses still retain some “memory” of the original seed mass distribution. The signature of the efficiency of the formation of SMBH seeds will consequently be stronger in isolated dwarf galaxies VLN2008 ().
One hopes that the next generation of 25-30m optical/IR telescopes operating at their diffraction limit ( milliarcsec) can provide the first constraints on the presence of IMBHs in dwarf galaxies, but the detection of gravitational waves from a central IMBH in a dwarf galaxy undergoing a merger is possibly a more promising probe. Dwarf galaxies have a very quiet merger history, hence we do not expect many IMBH-IMBH mergers involving dwarf galaxies at the present epoch, or in the low–redshift universe. The seed black hole mergers discussed in Section 4.1 probe a separate population of mergers, between the progenitors of galaxies which are more massive today. However, gravitational waves may also be generated in dwarf galaxies by mergers between the central IMBH and stellar remnants in the centre of the dwarf. These are analogous to the globular-cluster IMRI sources described in Section 3.1.
We can derive an estimate of the event rate based on the expected number of dwarf galaxies which can possibly host IMBHs in the interesting mass range. Theoretical models of SMBH formation and evolution, where the seeds of MBHs are Population III remnants VHM (), can be used to look for the distribution of IMBHs in dwarf galaxies. Using the dynamical model of vanwass09 (), we estimate a number density of IMBHs, – Mpc.
When we calculate the event rate of BH-IMBH mergers in dwarf galaxies, we have to further correct for the fact that only a small fraction of these tiny satellites do indeed form stars Bovill2009 (). Based on Gnedin2006 (), we estimate that a fraction of dwarfs in the km s range formed stars (which will eventually leave behind stellar mass BHs that can merge with the central IMBH). The number density of IMBHs that can be ET sources is therefore –Mpc. This number density is about an order of magnitude lower than the number density of globular clusters used to normalise the rates in Section 3.1.
The capture mechanisms that lead to IMRIs in dwarf galaxies are likely to be the same as those that operate in globular clusters. The event rate for the binary-hardening mechanism scales with the stellar density, , as , while the other mechanisms, such as direct capture, should scale approximately with (this is the same that enters Eqs. (2)–(3)). The core stellar densities in nearby dwarf galaxies are typically much lower than in core-collapsed globular clusters, e.g., the estimate for Fornax is pc Mateo:1998 () and for Sagittarius is pc Majewski:2005 (), compared to pc for globulars PryorMeylan:1993 (). The IMRI rates for dwarf galaxies are thus likely to be orders of magnitude lower than those for globular clusters. Therefore, although it is not inconceivable that ET will detect events from dwarf galaxies, any events would be serendipitous. Moreover, there are no obvious characteristics which would allow an observer to distinguish between an event in a dwarf galaxy from one in a globular cluster based on the GW signature alone. Nonetheless, it might be possible to make qualitative statements about dwarf galaxy IMBH populations. For instance, if ET does not detect any mergers between seed black holes at high redshift, of the type described in Section 4.1, it is very likely that BH seeds were heavy and not light. This would suggest dwarf galaxies would not contain light leftover BH seeds, and consequently that all of the observed IMRIs are occurring in globular clusters. Similarly, if seed mergers are detected but the rate of IMRIs is low or zero, it might suggest that IMBH formation in globular clusters is inefficient and any observed IMRIs are in dwarf galaxies. More refined modelling and detailed calculations are needed to understand/prove the robustness of these expectations, especially in view of the small number of seed black hole merger events and dwarf galaxy IMRIs that are predicted. In summary, while the dwarf galaxy channel should not be ignored completely, it is very unlikely to be a significant contributor to ET events or science.
5 Speculative sources
In this section we discuss some more speculative sources that might be observed by a future low-frequency ground-based interferometer such as ET. We examine first the possibility of observing orbiting or rotating white dwarfs near the high end of allowed masses, then discuss how the observation of eccentric compact binaries could illuminate their dynamical origin.
5.1 Orbiting white dwarfs
A gravitationally bound object of average density has a maximum orbital, rotational, or acoustic frequency . For neutron stars this maximum is Hz. White dwarfs are much more extended objects, but near their maximum masses their densities are sufficient to reach Hz. For example, from the classic work HamadaSalpeter:1961 (), a magnesium white dwarf with maximum mass has a radius cm and therefore Hz. In the few-Hz range, therefore, one will potentially see gravitational waves from the most massive white dwarfs.
If we consider specifically such a white dwarf in a binary orbit with a neutron star, black hole, or another white dwarf, then the orbital frequency at the point of tidal disruption of the dwarf depends weakly on the mass of the companion. For example, suppose that the equilibrium mass and radius of the white dwarf are respectively and , and that the companion is a compact object of mass . When the orbital separation is , tidal stripping begins (WigginsLai:2000, , and others). The orbital frequency at this point is
The gravitational wave frequency is , implying a maximum frequency of Hz for comparable-mass objects such as a neutron star and a heavy white dwarf, and a maximum that is 70% of this if the companion is a much more massive object such as an IMBH.
We have relatively few candidates for massive white dwarfs, hence although there is a significant literature related to lower-frequency radiation from white dwarf binaries (e.g., Farmer:2003 () and many subsequent papers) their numbers are difficult to estimate (see VennesKawka:2008 () for a recent discussion). Models of the mass distribution suggest that perhaps % of white dwarfs have masses near (e.g., see figure 10 of Catalan:2008 ()). Our requirement that both white dwarfs have masses near the maximum means that the mass ratio is greater than 2/3, and thus there will be a merger instead of stable mass transfer (see Marsh:2004 ()). If we estimate that that there are double white dwarf systems in a galaxy like the Milky Way Nelemans:2001popsynth (), and that % of the massive ones have semimajor axes that allow merger by gravitational radiation within yr (corresponding to the % merger fraction from Nelemans:2001popsynth ()), then we expect massive white-dwarf binaries to merge at a rate per galaxy of yr. At the high end this is similar to the low end of NS-NS merger rate estimates 2007PhR…442…75K (). If the ET is sensitive to such mergers out to Mpc, which may be optimistic given their low GW frequencies, one event per few years could be detected. Detection of these events would indicate rather precisely the maximum average density of white dwarfs, and would thus be a mechanism for establishing their mass-radius relation near the maximum mass.
5.2 Rotating hypermassive white dwarfs
Another possibility, suggested to us by Ott:private (), is that two white dwarfs with more typical masses might merge in a binary and produce a hypermassive white dwarf that spins rapidly enough that it is deformed into an ellipsoid. This is a promising candidate to explain some fraction of Type Ia supernovae 2009ApJ…699.2026R ().
To evaluate this prospect we note that if a Newtonian perfect fluid (a good model for a white dwarf) rotates uniformly then above a certain critical angular momentum, , for a given mass, , the equilibrium configuration splits off from the axisymmetric Maclaurin spheroids (which emit no gravitational radiation) to the Jacobi ellipsoids. If the three axes of the ellipsoids are , then, according to (Chandrasekhar:1969, , section 39), the critical angular momentum is
where . If two white dwarfs both of mass and radius spiral slowly together, then their angular momentum at the point of contact is . Since the equilibrium radius of the hypermassive object is smaller than the radii of the original white dwarfs, the angular momentum is sufficient to produce an ellipsoidal figure. Again from (Chandrasekhar:1969, , section 39), the angular velocity of this configuration will be and hence the dominant gravitational wave frequency will be .
The amplitude of gravitational waves depends on the ellipticity , where indicates the moment of inertia along axis . Near the critical angular momentum, slight changes in produce large changes in , and of several tenths is possible. Gravitational waves remove rotational energy from the star, such that
where . As a result, the characteristic spindown time is
The sweep rate at 1 Hz is, therefore, . For an integration of s the frequency would stay in a single frequency bin of Hz. Since the spindown rate will remain constant for a much longer time a search for a simple linear drift may make practical integrations over weeks to months. This would partially offset the low expected amplitudes. For comparison, continuous wave searches in LIGO are routinely done for spindown times as low as years at frequencies of Hz LIGOS5cw (), so a search for spindown times of years at Hz is certainly feasible.
Type Ia supernovae are estimated to occur once per 1000 years in galaxies such as the Milky Way 2009ApJ…699.2026R (), so even if only 1–10% of SNe Ia are binary mergers, the overall astrophysical rate is competitive with double neutron star mergers. Even though the ET sensitivity to gravitational waves from these binary white dwarf mergers will be much lower than for double neutron star mergers, the detection of gravitational waves from any such event may provide a new view on these important supernovae.
5.3 Eccentric binaries
In the sensitivity bands of second-generation gravitational wave detectors such as Advanced LIGO and Advanced Virgo, most compact binaries will be very close to circular. (For a proposed scenario where this may not be true, see OLeary:2008 (); another possibility includes direct captures of compact objects by IMBHs as precursors to eccentric IMRIs in globular clusters, although this formation mechanism is uncommon relative to the one described in Section 3.1, which will produce circular IMRIs). This is because for moderate to high eccentricities, gravitational radiation essentially reduces the semimajor axis of a binary while keeping the pericenter fixed. Therefore, to have palpable eccentricity at a given frequency, the pericenter at formation or at the last dynamical interaction must be inside the radius of a circular orbit at that frequency. For example, a binary of two neutron stars must have a pericenter less than 700 km to be significantly eccentric at a gravitational wave frequency Hz. This is highly improbable for a field binary, and is even difficult to arrange for binary-single scattering in dense stellar environments.
Somewhat higher eccentricities can be obtained via the Kozai secular resonance Kozai:1962 (). As explored in the context of black holes by MillerHamilton:2002b (); Wen:2003 (), a binary-binary interaction can result in a stable hierarchical triple in some tens of percent of encounters. If the inner binary and the outer tertiary have orbital planes that are inclined significantly with respect to each other, then over many orbital periods the inclination and eccentricity of the inner binary change periodically, leading at points in the cycle to very small pericenters and thus potentially observable eccentricity after the gravitational-wave driven inspiral. The eccentricity at 40 Hz is almost always very small (below 0.1), but at 10 Hz there are a few orientations in which the eccentricity can be a few tenths Wen:2003 (). At still lower frequencies the eccentricity will be yet higher, because for low eccentricities, , .
The preceding discussion implies that detector sensitivity at low frequencies will be important to determine the origin of compact binaries. In-situ formation from a massive main-sequence binary is still highly unlikely to produce detectable eccentricities: in order to have eccentricity at 1 Hz, the pericenter distance would have to be km immediately after the second supernova. In contrast, dynamical effects such as the Kozai resonance are expected to produce eccentric orbits at a few Hz. As a result, observation of a few BH-BH or BH-NS inspirals at a few Hz will illuminate their formation processes in a way that is not as easy at higher frequencies. We note, however, that simulations such as those in OLeary:2006 () suggest that of the few per year to few tens per year of black hole mergers in globulars that are expected to be seen with Advanced LIGO, less than 10% are initiated by the Kozai process. The greater reach of the Einstein telescope will enhance the total numbers, but binaries with palpable eccentricity in the ET band are still expected to be a minority.
6 Scientific impact of ET observations
ET detections of any of the systems described in this paper will yield important science products, which we now discuss.
The very existence of BHs in the – range is uncertain, so a single robust detection of an IMBH by ET will be of huge significance. If ET detects any seed black hole mergers at high redshift, it will be strong evidence that black hole seeds were light, which will help discriminate between light and heavy seed scenarios for the growth of structure in the Universe. Observations of mergers between more massive black holes with LISA do not have the same discriminating power, as they cannot distinguish between MBHs that formed through direct collapse or the collapse of a massive Pop III star and those that formed through a sequence of mergers SVH:2007 (); Gair:2009ET (). A significant number of ET detections of seed black hole mergers may provide constraints on the mass distribution of black hole seeds, and their early accretion history.
Detection of a significant number of IMRIs with ET will indicate that IMBHs form readily in globular clusters (since the rate of IMRIs in dwarf galaxies is so low). The characteristics of the IMRI events will provide constraints on the astrophysics of dense stellar environments, and on the efficiency of capture processes operating within them.
If ET detects white dwarfs undergoing tidal disruption, it will provide important constraints on the physics of degenerate matter, including the maximum density and mass that white dwarfs can reach. Detections of rotating hypermassive white dwarfs would provide information about proposed channels leading to supernovae. Finally, the detection of a significant population of eccentric coalescing binaries will shed light on the efficiency of the processes that drive eccentricity growth in binaries, such as the Kozai mechanism.
6.2 Fundamental physics
ET IMRI sources can also be used for testing aspects of relativity theory, in particular verifying that the central object is indeed a black hole as described by the Kerr metric of general relativity. This has been explored extensively in the context of extreme-mass-ratio inspiral events detectable by LISA (see, for example, Amaro:2007 () and references therein). In the course of an inspiral, the orbit of the smaller object traces out the spacetime geometry of the large body and hence the emitted gravitational waves encode a map of the spacetime structure. One way to characterize this is in terms of the multipole moments of the spacetime. It was demonstrated by Ryan Ryan:1995 (), for nearly circular and nearly equatorial orbits, that successive multipole moments of an arbitrary spacetime are encoded at different orders in an expansion of the orbital precession frequencies as functions of the orbital frequency. Since these frequencies can be measured from the emitted gravitational waves, a multipole map of the spacetime can in principle be measured. Similar multipole measurements are also possible from observations of ringdown radiation following mergers Berti:2006 (). For a Kerr black hole, the mass, , and angular momentum, , determine all higher-order mass, , and current, , multipole moments of the spacetime:
Measuring just three multipole moments and finding them to be inconsistent with this formula is therefore enough to demonstrate that the central object is not a Kerr black hole.
For IMRIs, it has been shown that Advanced LIGO could measure an fractional deviation in the mass quadrupole moment, , for typical systems Brown:2007 (). Corresponding results have not yet been computed for ET. However, ET will improve this significantly for two reasons — (i) the SNR of a source at fixed distance will increase by a factor of or more; and (ii) ET will observe the sources at lower frequencies. The ability to measure multipole moments improves significantly with the number of gravitational-wave cycles observed. At the leading-order Newtonian approximation, a system has cycles remaining until plunge when the frequency is Hz, but this increases to for a frequency of Hz, for Hz and for a frequency of Hz FinnThorne:2000 (). ET should thus be able to carry out tests of the Kerr nature of the central object that are significantly better than those possible with Advanced LIGO. Further research is required to quantify the improvement that will be possible, and how this will compare to expected results from LISA EMRI events.
There are various uncertainties which will affect the scientific impact of ET measurements discussed above. One important consideration is how to distinguish between IMBH events that arise from seed black holes and those that arise from IMBHs formed in globular clusters. Using ET measurements to constrain hierarchical structure formation relies on identification of mergers as seed black hole mergers, but, as we have seen, there may also be IMBH binary mergers in globular clusters. The masses and redshifts of the events may provide a robust discriminator, but more work is needed to understand if this is indeed the case, or whether other characteristic features exist that can be exploited.
The eventual sensitivity that is achieved by ET also has bearing on these results. The speculative sources that were discussed in Section 5 rely on ET having sensitivity in the –Hz band, and low-frequency sensitivity should also improve the accuracy of tests of relativity using IMRIs. ET may only have sensitivity down to a frequency of Hz, which will impact all of this science and perhaps eliminate the possibility of detecting gravitational radiation from massive white dwarfs. This must be properly quantified in the future.
Finally, there are open questions regarding ET data-analysis. The ET data stream will be very source-rich, and so the identification of individual sources of different types in the presence of this confusion will be a challenging problem. For instance, neutron star binary systems will create a confusion background near Hz Regimbau:2009 (). The data-analysis challenges for ET will inevitably change the SNRs required for detection of individual sources and therefore the rate predictions, and the accuracy with which source parameters can be estimated. However, the rate uncertainties arising from the data analysis will most likely be small compared to the order-of-magnitude uncertainties that are present in the astrophysical rate predictions.
We have discussed gravitational waves generated by intermediate-mass black holes as possible sources for the Einstein Telescope. Intermediate-mass black holes may be formed via two alternative channels — (i) they may be formed in the early Universe if MBH seeds are light (seed IMBH); (ii) they may form in globular clusters via runaway collisions between stars (cluster IMBH). In both cases, there are two distinct types of system that might be sources of gravitational waves for ET — (a) mergers between binaries containing two IMBHs; (b) mergers of stellar remnants with IMBHs (IMRIs).
Mergers between seed IMBHs occur following galaxy mergers during the hierarchical assembly of structure. If MBH seeds are light, ET could detect a few to a few tens of seed black hole merger events over three years at redshifts as high as –. An ET network would, in addition, be able to determine the luminosity distance to these events to an accuracy , which is sufficient to say confidently that an event involves intermediate-mass black holes and is occurring at high redshift. IMRIs involving seed IMBHs could occur in dwarf galaxies, but the event rate is probably very low, which makes it unlikely that this will be a significant contributor to the ET event rate. If cluster IMBHs form readily, binary IMBHs in globular clusters might be detected by ET at a rate of per year. Core-collapsed globular clusters are also a more promising host for IMRIs detectable by ET and the IMRI event rate for ET could be as high as a few hundred per year. However, there are significant uncertainties, not least of which is whether IMBHs form at all in the stellar environments of globular clusters.
The improved sensitivity of ET at low frequency may also allow the detection of several speculative sources. High-mass white dwarfs can survive tidal disruption long enough to reach orbital frequencies Hz in binaries. Hypermassive white dwarfs formed by the mergers of normal white dwarfs in binaries could also be sources for gravitational waves at frequencies around Hz as they will be rapidly rotating and can support relatively significant ellipticities. ET could detect these two types of event at a rate of one per few years, but this number is extremely uncertain. Finally, dynamical processes such as the Kozai mechanism can excite sufficiently high eccentricities in BH-BH and BH-NS binaries, that there would be significant residual eccentricity when their orbital frequency is in the –Hz range that ET will probe. These systems would circularize before reaching orbital frequencies in the Advanced LIGO band. ET might detect several eccentric binaries per year, but this rate depends on the fraction of binaries with residual eccentricity and as yet unknown details of the ET data analysis.
ET detections of any of these sources would have significant impact on our understanding of various astrophysical processes, as well as being useful for fundamental physics.
JG’s work is supported by a Royal Society University Research Fellowship. IM and MV acknowledge support from NASA ATP Grant NNX07AH22G. MCM acknowledges NASA ATP grant NNX08AH29G.
- (1) B. P. Abbott et al., Einstein@Home search for periodic gravitational waves in early S5 LIGO data, Phys. Rev. D80 (2009), 042003.
- (2) T. Abel, G. L. Bryan, and M. L. Norman, The Formation of the First Star in the Universe, Science 295 (2002), 93–98.
- (3) F. Acernese et al., Status of Virgo, Classical and Quantum Gravity 25 (2008), no. 11, 114045–+.
- (4) P. Ajith, Gravitational-wave data analysis using binary black-hole waveforms, Classical and Quantum Gravity 25 (2008), no. 11, 114033–+.
- (5) 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–+, Erratum: ibid. 79 (2009), 129901.
- (6) P. Ajith, M. Hannam, S. Husa, Y. Chen, B. Bruegmann, N. Dorband, D. Mueller, F. Ohme, D. Pollney, C. Reisswig, L. Santamaria, and J. Seiler, ”Complete” gravitational waveforms for black-hole binaries with non-precessing spins, ArXiv e-prints (2009).
- (7) P. Amaro-Seoane, C. Eichhorn, E. K. Porter, and R. Spurzem, Binaries of massive black holes in rotating clusters: dynamics, gravitational waves, detection and the role of eccentricity, MNRAS 401 (2010), 2268–2284.
- (8) P. Amaro-Seoane and M. Freitag, Intermediate-Mass Black Holes in Colliding Clusters: Implications for Lower Frequency Gravitational-Wave Astronomy, Astrophysical Journal Letters 653 (2006), L53–L56.
- (9) P. Amaro-Seoane, J. R. Gair, M. Freitag, M. C. Miller, I. Mandel, C. J. Cutler, and S. Babak, Intermediate and extreme mass-ratio inspirals – astrophysics, science applications and detection using LISA, Classical and Quantum Gravity 24 (2007), 113–+.
- (10) P. Amaro-Seoane and L. Santamaria, Detection of IMBHs with ground-based gravitational wave observatories: A biography of a binary of black holes, from birth to death, ArXiv e-prints (2009), 0910.0254.
- (11) L. Barack and N. Sago, Gravitational self-force on a particle in circular orbit around a Schwarzschild black hole, Phys. Rev. D 75 (2007), no. 6, 064021–+.
- (12) , Gravitational Self-Force Correction to the Innermost Stable Circular Orbit of a Schwarzschild Black Hole, Physical Review Letters 102 (2009), no. 19, 191101–+.
- (13) T. W. Baumgarte and S. L. Shapiro, Evolution of Rotating Supermassive Stars to the Onset of Collapse, Astrophysical Journal 526 (1999), 941–952.
- (14) M. C. Begelman, E. M. Rossi, and P. J. Armitage, Quasi-stars: accreting black holes inside massive envelopes, MNRAS 387 (2008), 1649–1659.
- (15) M. C. Begelman, M. Volonteri, and M. J. Rees, Formation of supermassive black holes by direct collapse in pre-galactic haloes, MNRAS 370 (2006), 289–298.
- (16) P.L. Bender et al., Lisa pre-phase a report; second edition, Tech. Report MPQ233, 1998.
- (17) C. T. Berghea, K. A. Weaver, E. J. M. Colbert, and T. P. Roberts, Testing the Paradigm that Ultraluminous X-ray Sources as a Class Represent Accreting Intermediate-Mass Black Holes, ArXiv e-prints (2008), 0807.1547.
- (18) E. Berti, V. Cardoso, and C. M. Will, Gravitational-wave spectroscopy of massive black holes with the space interferometer LISA, Phys. Rev. D 73 (2006), no. 6, 064030–+.
- (19) M. S. Bovill and M. Ricotti, Pre-Reionization Fossils, Ultra-Faint Dwarfs, and the Missing Galactic Satellite Problem, Astrophysical Journal 693 (2009), 1859–1870.
- (20) V. Bromm, P. S. Coppi, and R. B. Larson, Forming the First Stars in the Universe: The Fragmentation of Primordial Gas, ApJL 527 (1999), L5–L8.
- (21) V. Bromm and A. Loeb, Formation of the First Supermassive Black Holes, ApJ 596 (2003), 34–46.
- (22) , Formation of the First Supermassive Black Holes, ApJ 596 (2003), 34–46.
- (23) 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–+.
- (24) 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–+.
- (25) A. Buonanno, Y. Pan, H. P. Pfeiffer, M. A. Scheel, L. T. Buchman, and L. E. Kidder, Effective-one-body waveforms calibrated to numerical relativity simulations: Coalescence of nonspinning, equal-mass black holes, Phys. Rev. D 79 (2009), no. 12, 124028–+.
- (26) S. Catalán, J. Isern, E. García-Berro, and I. Ribas, The initial-final mass relationship of white dwarfs revisited: effect on the luminosity function and mass distribution, MNRAS 387 (2008), 1693–1706.
- (27) S. Chandrasekhar, Ellipsoidal figures of equilibrium: The Silliman Foundation Lectures, Yale University Press, New Haven, 1969.
- (28) B. Devecchi and M. Volonteri, Formation of the First Nuclear Clusters and Massive Black Holes at High Redshift, Astrophysical Journal 694 (2009), 302–313.
- (29) A. J. Farmer and E. S. Phinney, The gravitational wave background from cosmological compact binaries, MNRAS 346 (2003), 1197–1214.
- (30) 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.
- (31) L. Ferrarese and H. Ford, Supermassive Black Holes in Galactic Nuclei: Past, Present and Future Research, Space Science Reviews 116 (2005), 523–624.
- (32) L. Ferrarese and D. Merritt, A Fundamental Relation between Supermassive Black Holes and Their Host Galaxies, ApJ 539 (2000), L9–L12.
- (33) L. S. Finn and D. F. Chernoff, Observing binary inspiral in gravitational radiation: One interferometer, Phys. Rev. D 47 (1993), 2198–2219.
- (34) Lee Samuel Finn and Kip S. Thorne, Gravitational waves from a compact star in a circular, inspiral orbit, in the equatorial plane of a massive, spinning black hole, as observed by lisa, Phys. Rev. D62 (2000), 124021.
- (35) 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.
- (36) 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–+.
- (37) C. L. Fryer, S. E. Woosley, and A. Heger, Pair-Instability Supernovae, Gravity Waves, and Gamma-Ray Transients, ApJ 550 (2001), 372–382.
- (38) J. R. Gair, Probing black holes at low redshift using LISA EMRI observations, Classical and Quantum Gravity 26 (2009), no. 9, 094034–+.
- (39) J. R Gair, I. Mandel, A. Sesana, and A. Vecchio, Probing seed black holes using future gravitational-wave detectors, ArXiv e-prints (2009), 0907.3292.
- (40) K. Gebhardt, R. Bender, G. Bower, A. Dressler, S. M. Faber, A. V. Filippenko, R. Green, C. Grillmair, L. C. Ho, J. Kormendy, T. R. Lauer, J. Magorrian, J. Pinkney, D. Richstone, and S. Tremaine, A Relationship between Nuclear Black Hole Mass and Galaxy Velocity Dispersion, ApJ 539 (2000), L13–L16.
- (41) E. Glebbeek, E. Gaburov, S. E. de Mink, O. R. Pols, and S. F. Portegies Zwart, The evolution of runaway stellar collision products, A&A 497 (2009), 255–264.
- (42) N. Y. Gnedin and A. V. Kravtsov, Fossils of Reionization in the Local Group, Astrophysical Journal 645 (2006), 1054–1061.
- (43) H. Grote and the LIGO Scientific Collaboration, The status of GEO 600, Classical and Quantum Gravity 25 (2008), no. 11, 114043–+.
- (44) K. Gültekin, M. C. Miller, and D. P. Hamilton, Growth of Intermediate-Mass Black Holes in Globular Clusters, Astrophysical Journal 616 (2004), 221–230.
- (45) , Three-Body Dynamics with Gravitational Wave Emission, Astrophysical Journal 640 (2006), 156–166.
- (46) M. A. Gürkan, J. M. Fregeau, and F. A. Rasio, Massive Black Hole Binaries from Collisional Runaways, Astrophysical Journal 640 (2006), L39–L42.
- (47) M. A. Gürkan, M. Freitag, and F. A. Rasio, Formation of Massive Black Holes in Dense Star Clusters. I. Mass Segregation and Core Collapse, Astrophysical Journal 604 (2004), 632–652.
- (48) M. G. Haehnelt and M. J. Rees, The formation of nuclei in newly formed galaxies and the evolution of the quasar population, MNRAS 263 (1993), 168–178.
- (49) T. Hamada and E. E. Salpeter, Models for Zero-Temperature Stars., Astrophysical Journal 134 (1961), 683–+.
- (50) N. Häring and H.-W. Rix, On the Black Hole Mass-Bulge Mass Relation, Astrophysical Journal 604 (2004), L89–L92.
- (51) S. Hild, S. Chelkowski, and A. Freise, Pushing towards the ET sensitivity using ’conventional’ technology, ArXiv e-prints (2008), 0810.0604.
- (52) S. Hild, S. Chelkowski, A. Freise, J. Franc, N. Morgado, R. Flaminio, and R. DeSalvo, A Xylophone Configuration for a third Generation Gravitational Wave Detector, ArXiv e-prints (2009), 0906.2655.
- (53) D. W. Hogg, Distance measures in cosmology, ArXiv Astrophysics e-prints (1999), astro-ph/9905116.
- (54) C. Hopman, private communication, 2007.
- (55) F. Hoyle and W. A. Fowler, On the nature of strong radio sources, MNRAS 125 (1963), 169–+.
- (56) E. A. Huerta and J. R. Gair, Influence of conservative corrections on parameter estimation for extreme-mass-ratio inspirals, Phys. Rev. D 79 (2009), no. 8, 084021–+.
- (57) J. R. Hurley, Ratios of star cluster core and half-mass radii: a cautionary note on intermediate-mass black holes in star clusters, MNRAS 379 (2007), 93–99.
- (58) N. Ivanova, K. Belczynski, J. M. Fregeau, and F. A. Rasio, The evolution of binary fractions in globular clusters, MNRAS 358 (2005), 572–584.
- (59) V. Kalogera, K. Belczynski, C. Kim, R. O’Shaughnessy, and B. Willems, Formation of double compact objects, Physics Reports 442 (2007), 75–108.
- (60) S. M. Koushiappas, J. S. Bullock, and A. Dekel, Massive black hole seeds from low angular momentum material, MNRAS 354 (2004), 292–304.
- (61) Y. Kozai, Secular perturbations of asteroids with high inclination and eccentricity, Astronomical Journal 67 (1962), 591–+.
- (62) G. Lodato and P. Natarajan, Supermassive black hole formation during the assembly of pre-galactic discs, MNRAS 371 (2006), 1813–1823.
- (63) A. Loeb and F. A. Rasio, Collapse of primordial gas clouds and the formation of quasar black holes, ApJ 432 (1994), 52–61.
- (64) T. J. Maccarone and M. Servillat, Radio observations of NGC 2808 and other globular clusters: constraints on intermediate-mass black holes, MNRAS 389 (2008), 379–384.
- (65) P. Madau and M. J. Rees, Massive Black Holes as Population III Remnants, ApJL 551 (2001), L27–L30.
- (66) J. Magorrian, S. Tremaine, D. Richstone, R. Bender, G. Bower, A. Dressler, S. M. Faber, K. Gebhardt, R. Green, C. Grillmair, J. Kormendy, and T. Lauer, The Demography of Massive Dark Objects in Galaxy Centers, AJ 115 (1998), 2285–2305.
- (67) S. R. Majewski, P. M. Frinchaboy, W. E. Kunkel, R. Link, R. R. Muñoz, J. C. Ostheimer, C. Palma, R. J. Patterson, and D. Geisler, Exploring Halo Substructure with Giant Stars. VI. Extended Distributions of Giant Stars around the Carina Dwarf Spheroidal Galaxy: How Reliable Are They?, Astronomical Journal 130 (2005), 2677–2700.
- (68) R. K. Malbon, C. M. Baugh, C. S. Frenk, and C. G. Lacey, Black hole growth in hierarchical galaxy formation, MNRAS 382 (2007), 1394–1414.
- (69) I. Mandel, Spin distribution following minor mergers and the effect of spin on the detection range for low-mass-ratio inspirals, ArXiv e-prints (2007), 0707.0711.
- (70) 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.
- (71) I. Mandel and J. R. Gair, Can we detect intermediate mass ratio inspirals?, Classical and Quantum Gravity 26 (2009), no. 9, 094036–+.
- (72) T. R. Marsh, G. Nelemans, and D. Steeghs, Mass transfer between double white dwarfs, MNRAS 350 (2004), 113–128.
- (73) Mario Mateo, Dwarf galaxies of the local group, Annu. Rev. Astron. Astrophys. 36 (1998), 435–506.
- (74) M. C. Miller, Intermediate-Mass Black Holes as LISA Sources, ArXiv e-prints (2008), 0812.3028.
- (75) M. C. Miller and E. J. M. Colbert, Intermediate-Mass Black Holes, International Journal of Modern Physics D 13 (2004), 1–64.
- (76) M. C. Miller and D. P. Hamilton, Four-Body Effects in Globular Cluster Black Hole Coalescence, Astrophysical Journal 576 (2002), 894–898.
- (77) , Production of intermediate-mass black holes in globular clusters, MNRAS 330 (2002), 232–240.
- (78) M. Milosavljević, S. M. Couch, and V. Bromm, Accretion Onto Intermediate-Mass Black Holes in Dense Protogalactic Clouds, Astrophysical Journal 696 (2009), L146–L149.
- (79) H. Mouri and Y. Taniguchi, Mass Segregation in Star Clusters: Analytic Estimation of the Timescale, Astrophysical Journal 580 (2002), 844–849.
- (80) , Runaway Merging of Black Holes: Analytical Constraint on the Timescale, Astrophysical Journal 566 (2002), L17–L20.
- (81) F. Nakamura and M. Umemura, On the Initial Mass Function of Population III Stars, ApJ 548 (2001), 19–32.
- (82) G. Nelemans, L. R. Yungelson, S. F. Portegies Zwart, and F. Verbunt, Population synthesis for double white dwarfs . I. Close detached systems, A&A 365 (2001), 491–507.
- (83) E. Noyola, K. Gebhardt, and M. Bergmann, Gemini and Hubble Space Telescope Evidence for an Intermediate-Mass Black Hole in Centauri, Astrophysical Journal 676 (2008), 1008–1015.
- (84) 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.
- (85) R. M. O’Leary, F. A. Rasio, J. M. Fregeau, N. Ivanova, and R. O’Shaughnessy, Binary Mergers and Growth of Black Holes in Dense Star Clusters, Astrophysical Journal 637 (2006), 937–951.
- (86) C. Ott, private communication, 2009.
- (87) J. A. Peacock, Gravitational lenses and. cosmological evolution, MNRAS 199 (1982), 987–1006.
- (88) P. C. Peters, Gravitational radiation and the motion of two point masses, Phys. Rev. 136 (1964), B1224.
- (89) E. Poisson, The Motion of Point Particles in Curved Spacetime, Living Reviews in Relativity 7 (2004), 6–+.
- (90) S. Portegies Zwart, S. McMillan, and M. Gieles, Young massive star clusters, ArXiv e-prints (2010), 1002.1961.
- (91) S. F. Portegies Zwart, H. Baumgardt, P. Hut, J. Makino, and S. L. W. McMillan, Formation of massive black holes through runaway collisions in dense young star clusters, Nature 428 (2004), 724–726.
- (92) S. F. Portegies Zwart and S. L. W. McMillan, Black hole mergers in the universe, Astrophysical Journal 528 (2000), L17.
- (93) W. H. Press and P. Schechter, Formation of Galaxies and Clusters of Galaxies by Self-Similar Gravitational Condensation, Astrophysical Journal 187 (1974), 425–438.
- (94) C. Pryor and G. Meylan, Velocity Dispersions for Galactic Globular Clusters, Structure and Dynamics of Globular Clusters (S. G. Djorgovski and G. Meylan, eds.), Astronomical Society of the Pacific Conference Series, vol. 50, January 1993, pp. 357–+.
- (95) G. D. Quinlan, The dynamical evolution of massive black hole binaries I. Hardening in a fixed stellar background, New Astronomy 1 (1996), 35–56.
- (96) M. J. Rees, Emission from the nuclei of nearby galaxies - Evidence for massive black holes, Structure and Properties of Nearby Galaxies (E. M. Berkhuijsen and R. Wielebinski, eds.), IAU Symposium, vol. 77, 1978, pp. 237–242.
- (97) T. Regimbau and S. A. Hughes, Gravitational-wave confusion background from cosmological compact binaries: Implications for future terrestrial detectors, Phys. Rev. D 79 (2009), no. 6, 062002–+.
- (98) A. J. Ruiter, K. Belczynski, and C. Fryer, Rates and Delay Times of Type Ia Supernovae, Astrophysical Journal 699 (2009), 2026–2036.
- (99) 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.
- (100) 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.
- (101) A. Sesana, F. Haardt, P. Madau, and M. Volonteri, Low-Frequency Gravitational Radiation from Coalescing Massive Black Hole Binaries in Hierarchical Cosmologies, Astrophysical Journal 611 (2004), 623–632.
- (102) A. Sesana, M. Volonteri, and F. Haardt, The imprint of massive black hole formation models on the LISA data stream, MNRAS 377 (2007), 1711–1716.
- (103) C. Shapiro, D. Bacon, M. Hendry, and B. Hoyle, Delensing Gravitational Wave Standard Sirens with Shear and Flexion Maps, ArXiv e-prints (2009).
- (104) I. Shlosman, J. Frank, and M. C. Begelman, Bars within bars - A mechanism for fuelling active galactic nuclei, Nature 338 (1989), 45–47.
- (105) D. Sigg and the LIGO Scientific Collaboration, Status of the LIGO detectors, Classical and Quantum Gravity 25 (2008), no. 11, 114041–+.
- (106) M. Spaans and J. Silk, Pregalactic Black Hole Formation with an Atomic Hydrogen Equation of State, Astrophysical Journal 652 (2006), 902–906.
- (107) C. C. Steidel, K. L. Adelberger, M. Giavalisco, M. Dickinson, and M. Pettini, Lyman-Break Galaxies at z¿4 and the Evolution of the Ultraviolet Luminosity Density at High Redshift, Astrophysical Journal 519 (1999), 1–17.
- (108) Y. Taniguchi, Y. Shioya, T. G. Tsuru, and S. Ikeuchi, Formation of Intermediate-Mass Black Holes in Circumnuclear Regions of Galaxies, PASJ 52 (2000), 533–537.
- (109) M. Trenti, Dynamical evidence for intermediate mass black holes in old globular clusters, ArXiv Astrophysics e-prints (2006), arXiv:astro-ph/0612040.
- (110) S. Van Wassenhove, M. Volonteri, M. G. Walker, and J. R. Gair, Massive black holes lurking in Milky Way satellites, ArXiv e-prints (2010).
- (111) S. Vennes and A. Kawka, On the empirical evidence for the existence of ultramassive white dwarfs, MNRAS 389 (2008), 1367–1374.
- (112) M. Volonteri, F. Haardt, and P. Madau, The Assembly and Merging History of Supermassive Black Holes in Hierarchical Models of Galaxy Formation, ApJ 582 (2003), 559–573.
- (113) M. Volonteri, G. Lodato, and P. Natarajan, The evolution of massive black hole seeds, MNRAS 383 (2008), 1079–1088.
- (114) M. Volonteri, R. Salvaterra, and F. Haardt, Constraints on the accretion history of massive black holes from faint X-ray counts, MNRAS 373 (2006), 121–127.
- (115) M. G. Walker, M. Mateo, and E. W. Olszewski, Stellar Velocities in the Carina, Fornax, Sculptor, and Sextans dSph Galaxies: Data From the Magellan/MMFS Survey, Astronomical Journal 137 (2009), 3100–3108.
- (116) Y. Wang, D. E. Holz, and D. Munshi, A Universal Probability Distribution Function for Weak-lensing Amplification, Astrophysical Journal 572 (2002), L15–L18.
- (117) S. Weinberg, Apparent luminosities in a locally inhomogeneous universe, Astrophysical Journal 208 (1976), L1–L3.
- (118) L. Wen, On the Eccentricity Distribution of Coalescing Black Hole Binaries Driven by the Kozai Mechanism in Globular Clusters, Astrophysical Journal 598 (2003), 419–430.
- (119) P. Wiggins and D. Lai, Tidal Interaction between a Fluid Star and a Kerr Black Hole in Circular Orbit, Astrophysical Journal 532 (2000), 530–539.
- (120) N. Yoshida, K. Omukai, L. Hernquist, and T. Abel, Formation of Primordial Stars in a CDM Universe, Astrophysical Journal 652 (2006), 6–25.
- (121) L. R. Yungelson, E. P. J. van den Heuvel, J. S. Vink, S. F. Portegies Zwart, and A. de Koter, On the evolution and fate of super-massive stars, A&A 477 (2008), 223–237.
- (122) Q. Zhang and S. M. Fall, The Mass Function of Young Star Clusters in the “Antennae” Galaxies, Astrophysical Journal 527 (1999), L81–L84.