No Arabic abstract
The first detections of black hole - neutron star mergers (GW200105 and GW200115) by the LIGO-Virgo-Kagra Collaboration mark a significant scientific breakthrough. The physical interpretation of pre- and post-merger signals requires careful cross-examination between observational and theoretical modelling results. Here we present the first set of black hole - neutron star simulations that were obtained with the numerical-relativity code BAM. Our initial data are constructed using the public LORENE spectral library which employs an excision of the black hole interior. BAM, in contrast, uses the moving-puncture gauge for the evolution. Therefore, we need to ``stuff the black hole interior with smooth initial data to evolve the binary system in time. This procedure introduces constraint violations such that the constraint damping properties of the evolution system are essential to increase the accuracy of the simulation and in particular to reduce spurious center-of-mass drifts. Within BAM we evolve the Z4c equations and we compare our gravitational-wave results with those of the SXS collaboration and results obtained with the SACRA code. While we find generally good agreement with the reference solutions and phase differences $lesssim 0.5$ rad at the moment of merger, the absence of a clean convergence order in our simulations does not allow for a proper error quantification. We finally present a set of different initial conditions to explore how the merger of black hole neutron star systems depends on the involved masses, spins, and equations of state.
Gravitational waves from binary neutron star (BNS) and black hole/neutron star (BHNS) inspirals are primary sources for detection by the Advanced Laser Interferometer Gravitational-Wave Observatory. The tidal forces acting on the neutron stars induce changes in the phase evolution of the gravitational waveform, and these changes can be used to constrain the nuclear equation of state. Current methods of generating BNS and BHNS waveforms rely on either computationally challenging full 3D hydrodynamical simulations or approximate analytic solutions. We introduce a new method for computing inspiral waveforms for BNS/BHNS systems by adding the post-Newtonian (PN) tidal effects to full numerical simulations of binary black holes (BBHs), effectively replacing the nontidal terms in the PN expansion with BBH results. Comparing a waveform generated with this method against a full hydrodynamical simulation of a BNS inspiral yields a phase difference of $<1$ radian over $sim 15$ orbits. The numerical phase accuracy required of BNS simulations to measure the accuracy of the method we present here is estimated as a function of the tidal deformability parameter ${lambda}$.
Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant $ell=m=2$ gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of $sim 10^{-4}$. The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% ($90^{text{th}}$ percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms .
Only numerical relativity simulations can capture the full complexities of binary black hole mergers. These simulations, however, are prohibitively expensive for direct data analysis applications such as parameter estimation. We present two new fast and accurate surrogate models for the outputs of these simulations: the first model, NRSur7dq4, predicts the gravitational waveform and the second model, RemnantModel, predicts the properties of the remnant black hole. These models extend previous 7-dimensional, non-eccentric precessing models to higher mass ratios, and have been trained against 1528 simulations with mass ratios $qleq4$ and spin magnitudes $chi_1,chi_2 leq 0.8$, with generic spin directions. The waveform model, NRSur7dq4, which begins about 20 orbits before merger, includes all $ell leq 4$ spin-weighted spherical harmonic modes, as well as the precession frame dynamics and spin evolution of the black holes. The final black hole model, RemnantModel, models the mass, spin, and recoil kick velocity of the remnant black hole. In their training parameter range, both models are shown to be more accurate than existing models by at least an order of magnitude, with errors comparable to the estimated errors in the numerical relativity simulations. We also show that the surrogate models work well even when extrapolated outside their training parameter space range, up to mass ratios $q=6$.
We present the first set of numerical relativity simulations of binary neutron mergers that include spin precession effects and are evolved with multiple resolutions. Our simulations employ consistent initial data in general relativity with different spin configurations and dimensionless spin magnitudes $sim 0.1$. They start at a gravitational-wave frequency of $sim392$~Hz and cover more than $1$ precession period and about 15 orbits up to merger. We discuss the spin precession dynamics by analyzing coordinate trajectories, quasi-local spin measurements, and energetics, by comparing spin aligned, antialigned, and irrotational configurations. Gravitational waveforms from different spin configuration are compared by calculating the mismatch between pairs of waveforms in the late inspiral. We find that precession effects are not distinguishable from nonprecessing configurations with aligned spins for approximately face-on binaries, while the latter are distinguishable from a nonspinning configurations. Spin precession effects are instead clearly visible for approximately edge-on binaries. For the parameters considered here, precession does not significantly affect the characteristic postmerger gravitational-wave frequencies nor the mass ejection. Our results pave the way for the modeling of spin precession effects in the gravitational waveform from binary neutron star events.
The RIT numerical relativity group is releasing a public catalog of black-hole-binary waveforms. The initial release of the catalog consists of 126 recent simulations that include precessing and non precessing systems with mass ratios $q=m_1/m_2$ in the range $1/6leq qleq1$. The catalog contains information about the initial data of the simulation, the waveforms extrapolated to infinity, as well as information about the peak luminosity and final remnant black hole properties. These waveforms can be used to independently interpret gravitational wave signals from laser interferometric detectors and