Gravitational waves from eccentric intermediate-mass black hole binaries
If binary intermediate-mass black holes (IMBHs; with masses between 100 and ) form in dense stellar clusters, their inspiral will be detectable with the planned Laser Interferometer Space Antenna (LISA) out to several Gpc. Here we present a study of the dynamical evolution of such binaries using a combination of direct -body techniques (when the binaries are well separated) and three-body relativistic scattering experiments (when the binaries are tight enough that interactions with stars occur one at a time). We find that for reasonable IMBH masses there is only a mild effect on the structure of the surrounding cluster even though the binary binding energy can exceed the binding energy of the cluster. We demonstrate that, contrary to standard assumptions, the eccentricity in the LISA band can be in some cases as large as and that it induces a measurable phase difference from circular binaries in the last year before merger. We also show that, even though energy input from the binary decreases the density of the core and slows down interactions, the total time to coalescence is short enough (typically less than a hundred million years) that such mergers will be unique snapshots of clustered star formation.
Subject headings:Black hole physics, gravitational waves, stellar dynamics, methods: N-body simulations
The existence of intermediate-mass black holes (IMBHs; masses ) is not as certain as that of stellar-mass or supermassive black holes because there is as yet no conclusively established dynamical mass for any candidate, although there is strong circumstantial evidence for this mass range in several cases (see Miller & Colbert 2004 and references therein for a review). Mergers of IMBHs would, however, be strong sources of gravitational waves.
The best studied scenario is the runaway growth of a star in a young cluster via physical collisions among the most massive stars in the center, which have sunk through mass segregation (Portegies Zwart & McMillan 2000; Gürkan et al. 2004; Portegies Zwart et al. 2004; Freitag et al. 2006). Recently, Gürkan et al. (2006) addressed the same configuration but added a fraction of primordial binaries to the stellar system. Using a Monte-Carlo stellar-dynamics code, they found that not one but two very massive stars grow in rich clusters in which 10% or more of stars are in primordial hard binaries, suggesting the formation of two IMBHs. However, this result has not been confirmed yet using more accurate direct body simulations. Portegies Zwart et al. (2004) have a simulation with primordial binaries but they do not see this formation, though it is also currently unclear how different core concentrations will affect binary IMBH formation with a certain fraction of primordial binaries. It is also possible that wind losses may drive away mass more rapidly than it accretes through further collisions (see Belkus, van Bever, & Vanbeveren 2007), although this relies on uncertain extrapolations from the that is the top of their range (see their Table 2) to the masses observed in N-body simulations.
Fregeau et al. (2006) considered for the first time the possibility that such a binary could be observed thanks to the emission of gravitational waves in the coalescence phase and estimated that one can expect the Laser Interferometer Space Antenna (LISA) to detect tens of them depending on the distribution of cluster masses and densities. Amaro-Seoane & Freitag (2006) addressed the evolution of a binary of two IMBHs formed as the result of the collision of two independent stellar clusters and followed the parameters of the binary orbit down to the region in which it will emit gravitational waves in the Hz LISA domain. To do this, they combined direct-summation simulations with an analytical model to evolve the binary from a point in which it was hard.
Here we assume that an IMBH binary has been produced in a single dense stellar cluster, and study the subsequent sinking of the IMBHs and the evolution and properties of the binary when it forms. In § 2 we discuss our numerical method, which combines direct -body studies with three-body scattering integrations. In § 3 we discuss the astrophysical implications of our results.
2. Evolution of the IMBH pair: Numerical method
2.1. Direct body simulations
Direct body codes integrate all gravitational accelerations in a stellar system without supposing any special symmetries. They are thus the most general and robust tools for numerical analysis of stellar clusters (Aarseth 1999, 2003). The code we use, NBODY4, includes a variety of sophisticated approaches that improve speed and accuracy, including KS regularization (Kustaanheimo & Stiefel 1965), as well as triple (3-body subsystems), quad (4-body subsystems), and chain regularization (Aarseth 1999, 2003). It also does not make use of any softening, which would lead to unrealistic evolution of the orbital parameters of the binary of massive black holes. The disadvantage of this or any direct body code is the required computational time. However, our calculations are accelerated thanks to the special-purpose hardware GRAPE-6A single PCI cards of the AEI cluster Tuffstein used for the simulations. Each card has a peak performance of 130 Gflops (Fukushige et al. 2005), so that a single node is comparable to a cluster of 100 individual CPUs working in parallel.
|128002||1900 + 380||6||Kroupa||0.1|
Table 1 gives the initial conditions for the different simulations that we feature. We ran six cases with varying number densities and concentrations, of which two had a Kroupa (2001) mass function instead of single-mass stars. The IMBHs have equal mass except in simulation , which has a mass ratio of 5. In our simulations the individual time steps led to fractional energy errors that were always less than per N-body unit of time and, globally, the total energy error of the cluster (i.e. the accumulated error in the integration of all particles) is in the case of our fiducial model ().
2.2. Evolution of the binary: gravitational radiation versus dynamics
Our approach is to evolve the cluster up to Myrs using the direct code with a central bound binary. As we discuss below, we can see in Fig.(1) that the stellar cluster experiences a very moderate expansion during the hardening of the IMBHs. Once the IMBHs are hard enough relative to each other they can be treated as an isolated binary that interacts occasionally with a passing star. We note that although the subsequent evolution of the binary will not be identical to that from the N-body runs, due to the stochastic nature of the encounters, the general development is similar. Such interactions tend to increase the binding energy of the binary, hence shrinking its semi-major axis. The eccentricity is also changed, both by Newtonian three-body interactions (which can increase or decrease the eccentricity; see Sesana et al. 2007 for a recent treatment) and by gravitational radiation, which circularizes the binary. The combination of the two determines the eccentricity of the binary when it enters the frequency range of LISA.
In Figure 2 we show the inspiral of the binary for models A and C. The irregular lines correspond to the body simulations. We take the last point of these evolutions and the number density of field stars as input to relativistic scattering experiments which we performed following Gültekin et al. (2006). The equations of motion we use for the three-body encounters include relativistic precession to first post-Newtonian order, as well as radiation reaction caused by gravitational waves. Between encounters, we evolve the semimajor axis and eccentricity of the IMBH binary using the Peters quadrupolar formulae (Peters 1964). The stars that interact with the binary are sent with a velocity at infinity of =10 km s, typical of cluster velocity dispersions. The interaction time is drawn from an exponential distribution with a mean time of , where is the stellar number density (taken from the body simulations) and is the scattering cross section including gravitational focusing. The typical region in which the IMBH binary wanders is larger than its radius of influence, hence there is no loss cone as there is for supermassive black holes.
We ran four sets of 40 simulations, two sets starting at large separations with zero eccentricity that established agreement with the direct-summation body simulations and two sets at the endpoints of Models A and B. In none of these runs was the binary itself ejected from the cluster, as expected given its large mass. For the Model A endpoint run (binary mass , initial semimajor axis AU, and initial eccentricity ), the eccentricity was when the gravitational wave frequency (equal to twice the orbital frequency) was at the Hz low end of the LISA band; for the Model B endpoint run (binary mass , initial semimajor axis AU, and initial eccentricity ) we found . We show the envelope of the Model B endpoint runs in the left panel of Figure 3. We also did scattering experiments corresponding to the endpoint of Model C, which had a Kroupa mass function. The results are shown in the right panel of Figure 3. Compared to the single-mass runs we see considerably greater variance in the eccentricity as a function of semimajor axis, and although the range of eccentricities in the LISA band overlaps those in the single-mass runs there are also a number of cases in which the binary is nearly circular by the time the gravitational wave frequency reaches Hz. This could be a general feature of scattering interactions when there is a broad mass function, but we have not performed enough runs to determine this with confidence.
3. Discussion and Conclusions
We have addressed the inspiral of two massive black holes in a single young stellar cluster. Our three main results are: (1) the cluster itself experiences only mild structural changes as a result of the inspiral, (2) the coalescence takes a short enough time (typically 100 Myrs) that mergers occurs close to the time of formation of the cluster, and (3) there is a significant residual eccentricity by the time the binary enters the LISA band. We now discuss these conclusions in order.
The stability of clusters against IMBH mergers is consistent with analytic expectations even though the binding energy of the IMBH binary can exceed the total binding energy of the cluster by a significant factor. To see this, consider a circular IMBH binary of component masses and , with total mass and reduced mass . As shown by Quinlan (1996), a star with low speed at infinity that interacts with the binary will typically be ejected with a speed , where is the relative speed of the two objects. For equal masses , this is . Suppose now that the cluster has an escape speed . If then the star will be retained and share its kinetic energy with the cluster. Otherwise, the star will be ejected from the cluster without depositing significant energy, because the dynamical time of escape is much less than the relaxation time (which is the time required for the star to give up energy). The binding energy of the IMBH binary when will be . In comparison, if the cluster has a three-dimensional velocity dispersion and a mass , the binding energy of the cluster is . The ratio is then . Typically , so only if the larger black hole mass is could the release of energy unbind the cluster. We also note that subsequent to this point, the loss of mass from stars being thrown out would also soften the cluster. However, since typically interaction with of order the binary mass changes the semimajor axis by a factor of , just in stars will shrink the binary by enough of a factor to produce coalescence. Therefore, as verified by our numerical simulations, hardening of an IMBH binary has only a minor effect on the cluster.
For the time to merger, we note that hardening from large separations to a few hundred AU takes Myr, based on our simulations. Our three-body runs then indicate that the total time from that point to merger is virtually always less than 10 Myr, meaning that conservatively the total time from formation to merger is less than yr. This is significantly shorter than the age of the universe. One consequence of this is that if star formation in massive clusters was more common at redshift than it is now, and if binary IMBH formation was also thus more common, then LISA observations of IMBH mergers will serve as a unique snapshot of star formation as well as of cluster dynamics (see also Fregeau et al. 2006).
Figure (3) shows that the eccentricity of the binary will be in the range when the dominant gravitational wave frequency is of Hz. Consistent with Quinlan (1996) we find that the eccentricity does not undergo a random walk, but instead tends to higher eccentricities when the binary is hard but before gravitational radiation circularization is important. As discussed in section 4 of Amaro-Seoane & Freitag (2006), a residual eccentricity will induce a difference in the phase evolution of the second harmonic compared to a circular orbit, even if it is as small as 0.07, as Amaro-Seoane & Freitag (2006) found. In our case, if we use an eccentricity in equation (4) of Amaro-Seoane & Freitag (2006), we find that the accumulated phase shift if observations cover a time of at least
before merger, where is the redshift. If we set , then we have to cover a time days before merger. This means that if we are able to observe the system during that period of time before the final coalescence, we will recover enough information to determine that the orbit is not circular. On the other hand, if we use a residual eccentricity of 0.07, as in Amaro-Seoane & Freitag (2006), we would need 3-4 years of observation for a binary before merger.
In conclusion, if young massive clusters form binary IMBHs then they will be strong and moderately eccentric LISA sources that could serve as unique signposts of clustered star formation. The non-zero residual eccentricity has an impact on the detection of such sources, since it is generally assumed that an equal-mass massive binary will have a zero eccentricity when entering the LISA band. Our results show that in our scenario is non-negligible for certain cases -though for some other models it is very low but detectable-; notably, case and , which are the only models in which we have a mass fraction and thus, they are the more realistic ones. The process of formation must of course be studied carefully from the standpoints of stellar dynamics and merger product evolution, but if binary IMBHs can form then their mergers are promising sources for future LISA detections.
- Aarseth (1999) Aarseth, S. J. 1999, The Publications of the Astronomical Society of the Pacific, 111, 1333
- Aarseth (2003) Aarseth, S. J. 2003, Gravitational N-Body Simulations. ISBN 0521432723. Cambridge, UK: Cambridge University Press, November 2003.
- Amaro-Seoane & Freitag (2006) Amaro-Seoane, P., & Freitag, M. 2006, ApJ Lett., 653, L53
- Belkus, van Bever, & Vanbeveren (2007) Belkus, H., van Bever, J., & Vanbeveren, D. 2007, ApJ, 659, 1756
- Fregeau et al. (2006) Fregeau, J. M., Larson, S. L., Miller, M. C., O’Shaughnessy, R., & Rasio, F. A. 2006, ApJ, 646, L135
- Freitag et al. (2006) Freitag, M., Gürkan, M. A., & Rasio, F. A. 2006, MNRAS, 368, 141
- Fukushige et al. (2005) Fukushige, T., Makino, J., & Kawai, A. 2005, PASJ, 57, 1009
- Gültekin et al. (2006) Gültekin, K., Miller, M. C., & Hamilton, D. P. 2006, ApJ, 640, 156
- Gürkan et al. (2004) Gürkan, M. A., Freitag, M., & Rasio, F. A. 2004, ApJ, 604, 632
- Gürkan et al. (2006) Gürkan, M. A., Fregeau, J. M., & Rasio, F. A. 2006, ApJ, 640, L39
- King (1966) King, I. R. 1966, AJ, 71, 64
- Kroupa (2001) Kroupa, P. 2001, MNRAS, 322, 231
- Kustaanheimo & Stiefel (1965) Kustaanheimo, P. E., & Stiefel, E. L. 1965, J. Reine Angew. Math.
- Miller & Colbert (2004) Miller, M. C., & Colbert, E. J. M. 2004, International Journal of Modern Physics D, 13, 1
- Peters (1964) Peters, P. C. 1964, Physical Review, 136, 1224
- Portegies Zwart et al. (2004) Portegies Zwart, S. F., Baumgardt, H., Hut, P., Makino, J., & McMillan, S. L. W. 2004, Nat, 428, 724
- Portegies Zwart & McMillan (2000) Portegies Zwart, S. F., & McMillan, S. L. W. 2000, ApJ Lett., 528, L17
- Quinlan (1996) Quinlan, G. D. 1996, New Astronomy, 1, 255
- Sesana et al. (2007) Sesana, A., Haardt, F., & Madau, P. 2007, ArXiv e-prints, 710