Approximate Solutions to the Binary Black Hole Initial Value Problem

We present approximate analytical solutions to the Hamiltonian and momentum constraint equations, corresponding to systems composed of two black holes with arbitrary linear and angular momentum. The analytical nature of these initial data solutions makes them easier to implement in numerical evolutions than the traditional numerical approach of solving the elliptic equations derived from the Einstein constraints.Comment: 2 pages, to appear in the 9th Marcel Grossman Meeting Proceeding

Binary black hole mergers: large kicks for generic spin orientations

We present results from several simulations of equal mass black holes with spin. The spin magnitudes are $S/m^2=0.8$ in all cases, but we vary the spin orientations arbitrarily, in and outside the orbital plane. We find that in all but one case the final merged black hole acquires a kick of more than 1000 km/s, indicating that kicks of this magnitude are likely to be generic and should be expected for mergers with general spin orientations. The maximum kick velocity we find is 2500 km/s and occurs for initial spins which are anti-aligned in the initial orbital plane.Comment: 5 pages, 3 figures, to be submitted to PR

The final mass and spin of black hole mergers

We consider black holes resulting from binary black hole mergers. By fitting to numerical results we construct analytic formulas that predict the mass and spin of the final black hole. Our formulas are valid for arbitrary initial spins and mass ratios and agree well with available numerical simulations. We use our spin formula in the context of two common merger scenarios for supermassive galactic black holes. We consider the case of isotropically distributed initial spin orientations (when no surrounding matter is present) and also the case when matter closely aligns the spins with the orbital angular momentum. The spin magnitude of black holes resulting from successive generations of mergers (with symmetric mass ratio $\eta$) has a mean of $1.73\eta + 0.28$ in the isotropic case and 0.94 for the closely aligned case.Comment: Corrected typos, shortened introductio

Numerical Models of Spin-Orbital Coupling in Neutron Star Binaries

We present a new numerical scheme for solving the initial value problem for quasiequilibrium binary neutron stars allowing for arbitrary spins. We construct sequences of circular-orbit binaries of varying separation, keeping the rest mass and circulation constant along each sequence. The spin angular frequency of the stars is shown to vary along the sequence, a result that can be derived analytically in the PPN limit. This spin effect, in addition to leaving an imprint on the gravitational waveform emitted during binary inspiral, is measurable in the electromagnetic signal if one of the stars is a pulsar visible from Earth.Comment: 4 pages, 3 figures. Submitted to the Proceedings of the "X Marcel Grossmann Meeting on General Relativity" in Rio de Janeiro, Brazil, July 20-26 (2003

Solving the Initial Value Problem of two Black Holes

We solve the elliptic equations associated with the Hamiltonian and momentum constraints, corresponding to a system composed of two black holes with arbitrary linear and angular momentum. These new solutions are based on a Kerr-Schild spacetime slicing which provides more physically realistic solutions than the initial data based on conformally flat metric/maximal slicing methods. The singularity/inner boundary problems are circumvented by a new technique that allows the use of an elliptic solver on a Cartesian grid where no points are excised, simplifying enormously the numerical problem.Comment: 4 pages, 3 figures. Minor corrections, some points clarified, and one reference added. To appear in Phys. Rev. Let

Binary black holes on a budget: Simulations using workstations

Binary black hole simulations have traditionally been computationally very expensive: current simulations are performed in supercomputers involving dozens if not hundreds of processors, thus systematic studies of the parameter space of binary black hole encounters still seem prohibitive with current technology. Here we show how the multi-layered refinement level code BAM can be used on dual processor workstations to simulate certain binary black hole systems. BAM, based on the moving punctures method, provides grid structures composed of boxes of increasing resolution near the center of the grid. In the case of binaries, the highest resolution boxes are placed around each black hole and they track them in their orbits until the final merger when a single set of levels surrounds the black hole remnant. This is particularly useful when simulating spinning black holes since the gravitational fields gradients are larger. We present simulations of binaries with equal mass black holes with spins parallel to the binary axis and intrinsic magnitude of S/m^2= 0.75. Our results compare favorably to those of previous simulations of this particular system. We show that the moving punctures method produces stable simulations at maximum spatial resolutions up to M/160 and for durations of up to the equivalent of 20 orbital periods.Comment: 20 pages, 8 figures. Final version, to appear in a special issue of Class. Quantum Grav. based on the New Frontiers in Numerical Relativity Conference, Golm, July 200