by Konstantinidis, Simos and Kokkotas, Kostas D.
24 pages, 28 figures

We present a new C++ code for collisional N-body simulations of star clusters. The code uses the Hermite fourth-order scheme with block time steps, for advancing the particles in time, while the forces and neighboring particles are computed using the GRAPE-6 board. Special treatment is used for close encounters, binary and multiple sub-systems that either form dynamically or exist in the initial configuration. The structure of the code is modular and allows the appropriate treatment of more physical phenomena, such as stellar and binary evolution, stellar collisions and evolution of close black-hole binaries. Moreover, it can be easily modified so that the part of the code that uses GRAPE-6, could be replaced by another module that uses other accelerating-hardware like the Graphics Processing Units (GPUs). Appropriate choice of the free parameters give a good accuracy and speed for simulations of star clusters up to and beyond core collapse. Simulations of Plummer models consisting of equal-mass stars reached core collapse at t~17 half-mass relaxation times, which compares very well with existing results, while the cumulative relative error in the energy remained below 0.001. Also, comparisons with published results of other codes for the time of core collapse for different initial conditions, show excellent agreement. Simulations of King models with an initial mass-function, similar to those found in the literature, reached core collapse at t~0.17, which is slightly smaller than the expected result from previous works. Finally, the code accuracy becomes comparable and even better than the accuracy of existing codes, when a number of close binary systems is dynamically created in a simulation. This is due to the high accuracy of the method that is used for close binary and multiple sub-systems.

by Lubich, Christian and Walther, Benny and Bruegmann, Bernd
9 pages, 6 figures

We present a non-canonically symplectic integration scheme tailored to numerically computing the post-Newtonian motion of a spinning black-hole binary. Using a splitting approach we combine the flows of orbital and spin contributions. In the context of the splitting, it is possible to integrate the individual terms of the spin-orbit and spin-spin Hamiltonians analytically, exploiting the special structure of the underlying equations of motion. The outcome is a symplectic, time-reversible integrator, which can be raised to arbitrary order by composition. A fourth-order version is shown to give excellent behavior concerning error growth and conservation of energy and angular momentum in long-term simulations. Favorable properties of the integrator are retained in the presence of weak dissipative forces due to radiation damping in the full post-Newtonian equations.

by Manca, Gian Mario and Vallisneri, Michele
RevTeX4, 21 pages, 9 PDF figures

The efficient placement of signal templates in source-parameter space is a crucial requisite for exhaustive matched-filtering searches of modeled gravitational-wave sources. Unfortunately, the current placement algorithms based on regular parameter-space meshes are difficult to generalize beyond simple signal models with few parameters. Various authors have suggested that a general, flexible, yet efficient alternative can be found in randomized placement strategies such as random placement and stochastic placement, which enhances random placement by selectively rejecting templates that are too close to others. In this article we explore several theoretical and practical issues in randomized placement: the size and performance of the resulting template banks; the effects of parameter-space boundaries; the use of quasi-random (self avoiding) number sequences; most important, the implementation of these algorithms in curved signal manifolds with and without the use of a Riemannian signal metric, which may be difficult to obtain. Specifically, we show how the metric can be replaced with a discrete triangulation-based representation of local geometry. We argue that the broad class of randomized placement algorithms offers a promising answer to many search problems, but that the specific choice of a scheme and its implementation details will still need to be fine-tuned separately for each problem.

by Barack, Leor and Sago, Norichika
42 pages

We present a numerical code for calculating the local gravitational self-force acting on a pointlike particle in a generic (bound) geodesic orbit around a Schwarzschild black hole. The calculation is carried out in the Lorenz gauge: For a given geodesic orbit, we decompose the Lorenz-gauge metric perturbation equations (sourced by the delta-function particle) into tensorial harmonics, and solve for each harmonic using numerical evolution in the time domain (in 1+1 dimensions). The physical self-force along the orbit is then obtained via mode-sum regularization. The total self-force contains a dissipative piece as well as a conservative piece, and we describe a simple method for disentangling these two pieces in a time-domain framework. The dissipative component is responsible for the loss of orbital energy and angular momentum through gravitational radiation; as a test of our code we demonstrate that the work done by the dissipative component of the computed force is precisely balanced by the asymptotic fluxes of energy and angular momentum, which we extract independently from the wave-zone numerical solutions. The conservative piece of the self force does not affect the time-averaged rate of energy and angular-momentum loss, but it influences the evolution of the orbital phases; this piece is calculated here for the first time in eccentric strong-field orbits. As a first concrete application of our code we recently reported the value of the shift in the location and frequency of the innermost stable circular orbit due to the conservative self-force [Phys. Rev. Lett.\ {\bf 102}, 191101 (2009)]. Here we provide full details of this analysis, and discuss future applications.

by Lousto, Carlos O. and Nakano, Hiroyuki and Zlochower, Yosef and Campanelli, Manuela
4 pages, 5 figures, revtex4

We study black-hole binaries in the intermediate-mass-ratio regime 0.01 < q < 0.1 with a new technique that makes use of nonlinear numerical trajectories and efficient perturbative evolutions to compute waveforms at large radii for the leading and nonleading modes. As a proof-of-concept, we compute waveforms for q=1/10. We discuss applications of these techniques for LIGO/VIRGO data analysis and the possibility that our technique can be extended to produce accurate waveform templates from a modest number of fully-nonlinear numerical simulations.

by Mueller, Doreen and Bruegmann, Bernd
15 pages, submitted to CQG for NRDA 2009 conference proceedings

Moving puncture simulations of black hole binaries rely on a specific gauge choice that leads to approximately stationary coordinates near each black hole. Part of the shift condition is a damping parameter, which has to be properly chosen for stable evolutions. However, a constant damping parameter does not account for the difference in mass in unequal mass binaries. We introduce a position dependent shift damping that addresses this problem. Although the coordinates change, the changes in the extracted gravitational waves are small.

by Zenginoglu, Anil
9 pages, 8 figures

We study linear gravitational perturbations of Schwarzschild spacetime by solving numerically Regge-Wheeler-Zerilli equations in time domain using hyperboloidal surfaces and a compactifying radial coordinate. We stress the importance of including the asymptotic region in the computational domain in studies of gravitational radiation. The hyperboloidal approach should be helpful in a wide range of applications employing black hole perturbation theory.

by Pollney, Denis and Reisswig, Christian and Schnetter, Erik and Dorband, Nils and Diener, Peter

We present results from a new code for binary black hole evolutions using the moving-puncture approach, implementing finite differences in generalised coordinates, and allowing the spacetime to be covered with multiple communicating non-singular coordinate patches. Here we consider a regular Cartesian near zone, with adapted spherical grids covering the wave zone. The efficiencies resulting from the use of adapted coordinates allow us to maintain sufficient grid resolution to an artificial outer boundary location which is causally disconnected from the measurement. For the well-studied test-case of the inspiral of an equal-mass non-spinning binary (evolved for more than 8 orbits before merger), we determine the phase and amplitude to numerical accuracies better than 0.010% and 0.090% during inspiral, respectively, and 0.003% and 0.153% during merger. The waveforms, including the resolved higher harmonics, are convergent and can be consistently extrapolated to $latex r\to\infty$ throughout the simulation, including the merger and ringdown. Ringdown frequencies for these modes (to $latex (\ell,m)=(6,6)$) match perturbative calculations to within 0.01%, providing a strong confirmation that the remnant settles to a Kerr black hole with irreducible mass $latex M_{\rm irr} = 0.884355\pm20\times10^{-6}$ and spin $latex S_f/M_f^2 = 0.686923 \pm 10\times10^{-6}$

by Chung, Shin Kee and Wen, Linqing and Blair, David and Cannon, Kipp and Datta, Amitava
12 pages, 4 figures, submitted to CQG

We report a novel application of graphics processing units (GPUs) for the purpose of accelerating the search pipelines for gravitational waves from coalescing binaries of compact objects. A speed-up of 16 fold has been achieved compared with a single central processing unit (CPU). We show that substantial improvements are possible and discuss the reduction in CPU count required for the detection of inspiral sources afforded by the use of GPUs.

by Preto, Miguel and Saha, Prasenjit
21 pages, 5 figures, submitted to Astrophysical Journal

Stars near the Galactic center reach a few percent of light speed during pericenter passage, which makes post-Newtonian effects potentially detectable. We formulate the orbit equations in Hamiltonian form such that the $latex O(v^2/c^2)$ and $latex O(v^3/c^3)$ post-Newtonian effects of the Kerr metric appear as a simple generalization of the Kepler problem. A related perturbative Hamiltonian applies to photon paths. We then derive a symplectic integrator with adaptive time-steps, for fast and accurate numerical calculation of post-Newtonian effects. Using this integrator, we explore relativistic effects. Taking the star S2 as an example, we find that general relativity would contribute tenths of mas in astrometry and tens of $latex \rm km s^{-1}$ in kinematics. (For eventual comparison with observations, redshift and time-delay contributions from the gravitational field on light paths will need to be calculated, but we do attempt these in the present paper.) The contribution from stars, gas, and dark matter in the Galactic center region is still poorly constrained observationally, but current models suggest that the resulting Newtonian perturbation on the orbits could plausibly be of the same order as the relativistic effects for stars with semi-major axes $latex \gtrsim 0.01$ pc (or 250 mas). Nevertheless, the known and distinctive {\it time dependence} of the relativistic perturbations may make it possible to disentangle and extract both effects from observations.

by Lodato, G. and Nayakshin, S. and King, A. R. and Pringle, J. E.
11 pages, 10 figures, MNRAS accepted

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre–existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary’s gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework. The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high–photon–energy variability on the relatively long orbital timescale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary. We also find that the disc can drive the binary to merger on a reasonable timescale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e. $latex a_0 \lesssim 0.05$ pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is $latex \la 0.001$. Gas discs thus appear not to be an effective solution to the `last parsec’ problem for a significant class of mergers.

by Megevand, Miguel and Anderson, Matthew and Frank, Juhan and Hirschmann, Eric W. and Lehner, Luis and Liebling, Steven L. and Motl, Patrick M. and Neilsen, David
10 pages, 13 figures

The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk’s angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk’s angular momentum, the disk’s dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.

by Palenzuela, Carlos and Anderson, Matthew and Lehner, Luis and Liebling, Steven L. and Neilsen, David
4 pages, 4 figures

In addition to producing gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary’s dynamics induce a variability in possible electromagnetically induced emissions as well as a possible enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves.