No Arabic abstract
We evolve a binary black hole system bearing a mass ratio of $q=m_1/m_2=2/3$ and individual spins of $S^z_1/m_1^2=0.95$ and $S^z_2/m_2^2=-0.95$ in a configuration where the large black hole has its spin antialigned with the orbital angular momentum, $L^z$, and the small black hole has its spin aligned with $L^z$. This configuration was chosen to measure the maximum recoil of the remnant black hole for nonprecessing binaries. We find that the remnant black hole recoils at 500km/s, the largest recorded value from numerical simulations for aligned spin configurations. The remnant mass, spin, and gravitational waveform peak luminosity and frequency also provide a valuable point in parameter space for source modeling.
We demonstrate that numerical relativity codes based on the moving punctures formalism are capable of evolving nearly maximally spinning black hole binaries. We compare a new evolution of an equal-mass, aligned-spin binary with dimensionless spin chi=0.99 using puncture-based data with recent simulations of the SXS Collaboration. We find that the overlap of our new waveform with the published results of the SXS Collaboration is larger than 0.999. To generate our new waveform, we use the recently introduced HiSpID puncture data, the CCZ4 evolution system, and a modified lapse condition that helps keep the horizon radii reasonably large.
The construction of constraint-satisfying initial data is an essential element for the numerical exploration of the dynamics of compact-object binaries. While several codes have been developed over the years to compute generic quasi-equilibrium configurations of binaries comprising either two black holes, or two neutron stars, or a black hole and a neutron star, these codes are often not publicly available or they provide only a limited capability in terms of mass ratios and spins of the components in the binary. We here present a new open-source collection of spectral elliptic solvers that are capable of exploring the major parameter space of binary black holes (BBHs), binary neutron stars (BNSs), and mixed binaries of black holes and neutron stars (BHNSs). Particularly important is the ability of the spectral-solver library to handle neutron stars that are either irrotational or with an intrinsic spin angular momentum that is parallel to the orbital one. By supporting both analytic and tabulated equations of state at zero or finite temperature, the new infrastructure is particularly geared towards allowing for the construction of BHNS and BNS binaries. For the latter, we show that the new solvers are able to reach the most extreme corners in the physically plausible space of parameters, including extreme mass ratios and spin asymmetries, thus representing the most extreme BNS computed to date. Through a systematic series of examples, we demonstrate that the solvers are able to construct quasi-equilibrium and eccentricity-reduced initial data for BBHs, BNSs, and BHNSs, achieving spectral convergence in all cases. Furthermore, using such initial data, we have carried out evolutions of these systems from the inspiral to after the merger, obtaining evolutions with eccentricities $lesssim 10^{-4}-10^{-3}$, and accurate gravitational waveforms.
The standard post-Newtonian approximation to gravitational waveforms, called T-approximants, from non-spinning black hole binaries are known not to be sufficiently accurate close to the last stable orbit of the system. A new approximation, called P-approximants, is believed to improve the accuracy of the waveforms rendering them applicable up to the last stable orbit. In this study we apply P-approximants to the case of a test-particle in equatorial orbit around a Kerr black hole parameterized by a spin parameter q that takes values between -1 and 1. In order to assess the performance of the two approximants we measure their effectualness (i.e. larger overlaps with the exact signal), and faithfulness (i.e. smaller biases while measuring the parameters of the signal) with the exact (numerical) waveforms. We find that in the case of prograde orbits, that is orbits whose angular momentum is in the same sense as the spin angular momentum of the black hole, T-approximant templates obtain an effectualness of ~ 0.99 for spins q < 0.75. For 0.75 < q < 0.95, the effectualness drops to about 0.82. The P-approximants achieve effectualness of > 0.99 for all spins up to q = 0.95. The bias in the estimation of parameters is much lower in the case of P-approximants than T-approximants. We find that P-approximants are both effectual and faithful and should be more effective than T-approximants as a detection template family when q>0. For q<0 both T- and P-approximants perform equally well so that either of them could be used as a detection template family.
We explore the benefits of adapted gauges to small mass ratio binary black hole evolutions in the moving puncture formulation. We find expressions that approximate the late time behavior of the lapse and shift, $(alpha_0,beta_0)$, and use them as initial values for their evolutions. We also use a position and black hole mass dependent damping term, $eta[vec{x}_1(t),vec{x}_2(t),m_1,m_2]$, in the shift evolution, rather than a constant or conformal-factor dependent choice. We have found that this substantially reduces noise generation at the start of the numerical integration and keeps the numerical grid stable around both black holes, allowing for more accuracy with lower resolutions. We test our choices for this gauge in detail in a case study of a binary with a 7:1 mass ratio, and then use 15:1 and 32:1 binaries for a convergence study. Finally, we apply our new gauge to a 64:1 binary and a 128:1 binary to well cover the comparable and small mass ratio regimes.
We construct an approximate metric that represents the spacetime of spinning binary black holes (BBH) approaching merger. We build the metric as an analytical superposition of two Kerr metrics in harmonic coordinates, where we transform each black hole term with time-dependent boosts describing an inspiral trajectory. The velocities and trajectories of the boost are obtained by solving the post-Newtonian (PN) equations of motion at 3.5 PN order. We analyze the spacetime scalars of the new metric and we show that it is an accurate approximation of Einsteins field equations in vacuum for a BBH system in the inspiral regime. Furthermore, to prove the effectiveness of our approach, we test the metric in the context of a 3D general relativistic magneto-hydrodynamical (GRMHD) simulation of accreting mini-disks around the black holes. We compare our results with a previous well-tested spacetime construction based on the asymptotic matching method. We conclude that our new spacetime is well-suited for long-term GRMHD simulations of spinning binary black holes on their way to the merger.