Binary Black Hole Mergers from Planet-like Migrations

If supermassive black holes (BHs) are generically present in galaxy centers, and if galaxies are built up through hierarchical merging, BH binaries are at least temporary features of most galactic bulges. Observations suggest, however, that binary BHs are rare, pointing towards a binary lifetime far shorter than the Hubble time. We show that, regardless of the detailed mechanism, all stellar-dynamical processes are insufficient to reduce significantly the orbital separation once orbital velocities in the binary exceed the virial velocity of the system. We propose that a massive gas disk surrounding a BH binary can effect its merger rapidly, in a scenario analogous to the orbital decay of super-jovian planets due to a proto-planetary disk. As in the case of planets, gas accretion onto the secondary (here a supermassive BH) is integrally connected with its inward migration. Such accretion would give rise to quasar activity. BH binary mergers could therefore be responsible for many or most quasars.Comment: 8 pages, submitted to ApJ Letter

Pseudoparticle Multipole Method: A Simple Method to Implement High-Accuracy Treecode

In this letter we describe the pseudoparticle multipole method (P2M2), a new method to express multipole expansion by a distribution of pseudoparticles. We can use this distribution of particles to calculate high order terms in both the Barnes-Hut treecode and FMM. The primary advantage of P2M2 is that it works on GRAPE. GRAPE is a special-purpose hardware for the calculation of gravitational force between particles. Although the treecode has been implemented on GRAPE, we could handle terms only up to dipole, since GRAPE can calculate forces from point-mass particles only. Thus the calculation cost grows quickly when high accuracy is required. With P2M2, the multipole expansion is expressed by particles, and thus GRAPE can calculate high order terms. Using P2M2, we implemented an arbitrary-order treecode on GRAPE-4. Timing result shows GRAPE-4 accelerates the calculation by a factor between 10 (for low accuracy) to 150 (for high accuracy). Even on general-purpose programmable computers, our method offers the advantage that the mathematical formulae and therefore the actual program is much simpler than that of the direct implementation of multipole expansion.Comment: 6 pages, 4 figures, latex, submitted to ApJ Letter

Formation of Protoplanets from Massive Planetesimals in Binary Systems

More than half of stars reside in binary or multiple star systems and many planets have been found in binary systems. From theoretical point of view, however, whether or not the planetary formation proceeds in a binary system is a very complex problem, because secular perturbation from the companion star can easily stir up the eccentricity of the planetesimals and cause high-velocity, destructive collisions between planetesimals. Early stage of planetary formation process in binary systems has been studied by restricted three-body approach with gas drag and it is commonly accepted that accretion of planetesimals can proceed due to orbital phasing by gas drag. However, the gas drag becomes less effective as the planetesimals become massive. Therefore it is still uncertain whether the collision velocity remains small and planetary accretion can proceed, once the planetesimals become massive. We performed {\it N}-body simulations of planetary formation in binary systems starting from massive planetesimals whose size is about 100-500 km. We found that the eccentricity vectors of planetesimals quickly converge to the forced eccentricity due to the coupling of the perturbation of the companion and the mutual interaction of planetesimals if the initial disk model is sufficiently wide in radial distribution. This convergence decreases the collision velocity and as a result accretion can proceed much in the same way as in isolated systems. The basic processes of the planetary formation, such as runaway growth and oligarchic growth and final configuration of the protoplanets are essentially the same in binary systems and single star systems, at least in the late stage where the effect of gas drag is small.Comment: 26pages, 11 figures. ApJ accepte

The influence of gas expulsion and initial mass-segregation on the stellar mass-function of globular star clusters

Recently de Marchi, Paresce & Pulone (2007) studied a sample of twenty globular clusters and found that all clusters with high concentrations have steep stellar mass-functions while clusters with low concentration have comparatively shallow mass-functions. No globular clusters were found with a flat mass-function and high concentration. This seems curious since more concentrated star clusters are believed to be dynamically more evolved and should have lost more low-mass stars via evaporation, which would result in a shallower mass-function in the low-mass part. We show that this effect can be explained by residual-gas expulsion from initially mass-segregated star clusters, and is enhanced further through unresolved binaries. If gas expulsion is the correct mechanism to produce the observed trend, then observation of these parameters would allow to constrain cluster starting conditions such as star formation efficiency and the time-scale of gas expulsion.Comment: accepted for publication in MNRAS, 10 pages, 6 figure

The PCI Interface for GRAPE Systems: PCI-HIB

We developed a PCI interface for GRAPE systems. GRAPE(GRAvity piPE) is a special-purpose computer for gravitational N-body simulations. A GRAPE system consists of GRAPE processor boards and a host computer. GRAPE processors perform the calculation of gravitational forces between particles. The host computer performs the rest of calculations. The newest of GRAPE machines, the GRAPE-4, achieved the peak performance of 1.08 Tflops. The GRAPE-4 system uses TURBOChannel for the interface to the host, which limits the selection of the host computer. The TURBOChannel bus is not supported by any of recent workstations. We developed a new host interface board which adopts the PCI bus instead of the TURBOChannel. PCI is an I/O bus standard developed by Intel. It has fairly high peak transfer speed, and is available on wide range of computers, from PCs to supercomputers. Thus, the new interface allows us to connect GRAPE-4 to a wide variety of host computers. In test runs with a Barnes-Hut treecode, we found that the performance of new system with PCI interface is 40% better than that of the original system.Comment: 15 pages, 10 Postscript figures, 3 tables, Latex, submitted to Publications of the Astronomical Society of Japan. corrected figure 2 which contained non standard font

Monte Carlo Simulations of Globular Cluster Evolution - II. Mass Spectra, Stellar Evolution and Lifetimes in the Galaxy

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

Brownian Motion of Black Holes in Dense Nuclei

We evaluate the Brownian motion of a massive particle ("black hole") at the center of a galaxy using N-body simulations. Our galaxy models have power-law central density cusps like those observed at the centers of elliptical galaxies. The simulations show that the black hole achieves a steady-state kinetic energy that is substantially different than would be predicted based on the properties of the galaxy model in the absence of the black hole. The reason appears to be that the black hole responds to stars whose velocities have themselves been raised by the presence of the black hole. Over a wide range of density slopes and black hole masses, the black hole's mean kinetic energy is equal to what would be predicted under the assumption that it is in energy equipartition with stars lying within a distance ~r_h/2 from it, where r_h is the black hole's influence radius. The dependence of the Brownian velocity on black hole mass is approximately ~ 1/M^{1/(3-gamma)} with gamma the power-law index of the stellar density profile, rho~1/r^gamma. This is less steep than the 1/M dependence predicted in a model where the effect of the black hole on the stellar velocities is ignored. The influence of a stellar mass spectrum on the black hole's Brownian motion is also evaluated and found to be consistent with predictions from Chandrasekhar's theory. We use these results to derive a probability function for the mass of the Milky Way black hole based on a measurement of its proper motion velocity. Interesting constraints on M will require a velocity resolution exceeding 0.5 km/s.Comment: 11 pages, uses emulateapj.st

Evolution of Galactic Nuclei. I. orbital evolution of IMBH

Resent observations and theoretical interpretations suggest that IMBHs (intermediate-mass black hole) are formed in the centers of young and compact star clusters born close to the center of their parent galaxy. Such a star cluster would sink toward the center of the galaxy, and at the same time stars are stripped out of the cluster by the tidal field of the parent galaxy. We investigated the orbital evolution of the IMBH, after its parent cluster is completely disrupted by the tidal field of the parent galaxy, by means of large-scale N-body simulations. We constructed a model of the central region of our galaxy, with an SMBH (supermassive black hole) and Bahcall-Wolf stellar cusp, and placed an IMBH in a circular orbit of radius 0.086pc. The IMBH sinks toward the SMBH through dynamical friction, but dynamical friction becomes ineffective when the IMBH reached the radius inside which the initial stellar mass is comparable to the IMBH mass. This is because the IMBH kicks out the stars. This behavior is essentially the same as the loss-cone depletion observed in simulations of massive SMBH binaries. After the evolution through dynamical friction stalled, the eccentricity of the orbit of the IMBH goes up, resulting in the strong reduction in the merging timescale through gravitational wave radiation. Our result indicates that the IMBHs formed close to the galactic center can merge with the central SMBH in short time. The number of merging events detectable with DECIGO is estimated to be around 50 per year. Event rate for LISA would be similar or less, depending on the growth mode of IMBHs.Comment: 12 pages, 24 figures, submitted to Ap