Search

Numerical relativity (NR) simulations of binary black hole (BBH) systems provide the most accurate gravitational wave predictions, but at a high computational cost—especially when the black holes have nearly extremal spins (i.e. spins near the theoretical upper limit) or very unequal masses. Recently, the technique of reduced order modeling has enabled the construction of 'surrogate models' trained on an existing set of NR waveforms. Surrogate models enable the rapid computation of the gravitational waves emitted by BBHs. Typically these models are used for interpolation to compute gravitational waveforms for BBHs with mass ratios and spins within the bounds of the training set. Because simulations with nearly extremal spins are so technically challenging, surrogate models almost always rely on training sets with only moderate spins. In this paper, we explore how well surrogate models can extrapolate to nearly extremal spins when the training set only includes moderate spins. For simplicity, we focus on one-dimensional surrogate models trained on NR simulations of BBHs with equal masses and equal, aligned spins. We assess the performance of the surrogate models at higher spin magnitudes by calculating the mismatches between extrapolated surrogate model waveforms and NR waveforms, by calculating the differences between extrapolated and NR measurements of the remnant black-hole mass, and by testing how the surrogate model improves as the training set extends to higher spins. We find that while extrapolation in this one-dimensional case is viable for current detector sensitivities, surrogate models for next-generation detectors should use training sets that extend to nearly extremal spins. [ABSTRACT FROM AUTHOR]

Black hole-neutron star (BHNS) binaries are important sources of gravitational waves for second-generation interferometers, and BHNS mergers are also a proposed engine for short, hard gamma-ray bursts. The behavior of both the spacetime (and thus the emitted gravitational waves) and the neutron-star matter in a BHNS merger depend strongly and nonlinearly on the black hole's spin. While there is a significant possibility that astrophysical black holes could have spins that are nearly extremal (i.e. near the theoretical maximum), to date fully relativistic simulations of BHNS binaries have included black-hole spins only up to S/M² = 0.9, which corresponds to the black hole having approximately half as much rotational energy as possible, given the black hole's mass. In this paper, we present a new simulation of a BHNS binary with a mass ratio q = 3 and black-hole spin S/M² = 0.97, the highest simulated to date. We find that the black hole's large spin leads to the most massive accretion disc and the largest tidal tail outflow of any fully relativistic BHNS simulations to date, even exceeding the results implied by extrapolating results from simulations with lower black-hole spin. The disc appears to be remarkably stable. We also find that the high black-hole spin persists until shortly before the time of merger; afterward, both merger and accretion spin down the black hole. [ABSTRACT FROM AUTHOR]

Motivated by the possibility of observing gravitational waves from merging black holes whose spins are nearly extremal (i.e. 1 in dimensionless units), we present numerical waveforms from simulations of merging black holes with the highest spins simulated to date: (1) a 25.5-orbit inspiral, merger and ringdown of two holes with equal masses and spins of magnitude 0.97 aligned with the orbital angular momentum; and (2) a previously reported 12.5-orbit inspiral, merger and ringdown of two holes with equal masses and spins of magnitude 0.95 anti-aligned with the orbital angular momentum. First, we consider the horizon mass and spin evolution of the new aligned-spin simulation. During the inspiral, the horizon area and spin evolve in remarkably close agreement with Alvi's analytic predictions, and the remnant hole's final spin agrees reasonably well with several analytic predictions. We also find that the total energy emitted by a real astrophysical system with these parameters--almost all of which is radiated during the time included in this simulation--would be 10.952% of the initial mass at infinite separation. Second, we consider the gravitational waveforms for both simulations. After estimating their uncertainties, we compare the waveforms to several post-Newtonian approximants, finding significant disagreement well before merger, although the phase of the TaylorT4 approximant happens to agree remarkably well with the numerical prediction in the aligned-spin case. We find that the post-Newtonian waveforms have sufficient uncertainty that hybridized waveforms will require far longer numerical simulations (in the absence of improved post-Newtonian waveforms) for accurate parameter estimation of low-mass binary systems. [ABSTRACT FROM AUTHOR]

Simulations of binary black hole systems using the Spectral Einstein Code (SpEC) are done on a computational domain that excises the regions inside the black holes. It is imperative that the excision boundaries are outflow boundaries with respect to the hyperbolic evolution equations used in the simulation. We employ a time-dependent mapping between the fixed computational frame and the inertial frame through which the black holes move. The time-dependent parameters of the mapping are adjusted throughout the simulation by a feedback control system in order to follow the motion of the black holes, to adjust the shape and size of the excision surfaces so that they remain outflow boundaries, and to prevent large distortions of the grid. We describe in detail the mappings and control systems that we use. We show how these techniques have been essential in the evolution of binary black hole systems with extreme configurations, such as large spin magnitudes and high mass ratios, especially during the merger, when apparent horizons are highly distorted and the computational domain becomes compressed. The techniques introduced here may be useful in other applications of partial differential equations that involve time-dependent mappings. [ABSTRACT FROM AUTHOR]