Numerical-Relativity-Informed Effective-One-Body model for Black-Hole-Neutron-Star Mergers with Higher Modes and Spin Precession

Alejandra Gonzalez,Rossella Gamba,Matteo Breschi,Francesco Zappa,Gregorio Carullo,Sebastiano Bernuzzi,Alessandro Nagar
DOI: https://doi.org/10.1103/PhysRevD.107.084026
2022-12-08
Abstract:We present the first effective-one-body (EOB) model for generic-spins quasi-circular black-hole--neutron-star (BHNS) inspiral-merger-ringdown gravitational waveforms (GWs). Our model is based on a new numerical-relativity (NR) informed expression of the BH remnant and its ringdown. It reproduces the NR $(\ell,m)=(2,2)$ waveform with typical phase agreement of ${\lesssim0.5}\,$rad (${\lesssim1}\,$rad) to merger (ringdown). The maximum (minimum) mismatch between the $(2,2)$ and the NR data is 4% (0.6%). Higher modes (HMs) $(2,1)$, $(3,2)$, $(3,3)$, $(4,4)$ and $(5,5)$ are included and their mismatch with the available NR waveforms are up to (down to) a 60% (1%) depending on the inclination. Phase comparison with a 16 orbit precessing simulation shows differences within the NR uncertainties. We demonstrate the applicability of the model in GW parameter estimation by perfoming the first BHNS Bayesian analysis with HMs (and non-precessing spins) of the event GW190814, together with new $(2,2)$-mode analysis of GW200105 and GW200115. For the GW190814 study, the inclusion of HMs gives tighter parameter posteriors. The Bayes factors of our analyses on this event show decisive evidence for the presence of HMs, but no clear preference for a BHNS or a binary black hole (BBH) source. Similarly, we confirm GW200105 and GW200115 show no evidence for tidal effects.
General Relativity and Quantum Cosmology
What problem does this paper attempt to address?
This paper aims to solve the problem of gravitational waveform modeling in black - hole - neutron - star merger systems (BHNS). Specifically, the paper proposes a new Effective - One - Body (EOB) model for describing the gravitational waveforms (GWs) in the quasi - circular - orbit black - hole - neutron - star merger process with spin precession and high - order modes. This model is based on Numerical Relativity (NR) data and can accurately simulate the gravitational wave signals throughout the entire process from inspiral, merger to ringdown. ### Main problems and solutions 1. **Problems**: - How to construct a model that can accurately describe the gravitational waveform in the black - hole - neutron - star merger process? - How to include the effects of spin precession and high - order modes in the model? - How to ensure the robustness and accuracy of the model in different parameter spaces? 2. **Solutions**: - **Model design**: The paper proposes a new EOB model TEOBResumS, which combines Numerical Relativity data and is improved especially in describing the properties of the remnant black hole (remnant BH) and the characteristics of the ringdown stage. - **High - order modes**: The model includes high - order modes (Higher Modes, HMs), such as (2,1), (3,2), (3,3), (4,4) and (5,5), which are particularly important for describing unequal - mass binary systems. - **Spin precession**: The model takes into account the spin precession effect, which is achieved by introducing the spin - orbit interaction term. - **Data verification**: The accuracy and robustness of the model are verified by comparing with the existing Numerical Relativity data. The model performs well in different parameter spaces, especially in describing the post - merger ringdown stage. ### Specific implementation of the model 1. **Remnant black hole model**: - Use Numerical Relativity data fitting formulas to describe the mass \( M_{\bullet} \) and dimensionless spin \( a_{\bullet} \) of the remnant black hole. These formulas can accurately predict the properties of the remnant black hole under different symmetric mass ratios \( \nu \), initial black - hole spins \( a_{\text{BH}} \) and neutron - star tidal polarizabilities \( \Lambda \). - The formula form is as follows: \[ F: (\nu, a_{\text{BH}}, \Lambda) \rightarrow (M_{\bullet}, a_{\bullet}) \] 2. **Ringdown model**: - The ringdown emission must consider the effect of neutron - star tidal disruption, which will significantly suppress the excitation of quasi - normal modes (QNMs) in some parameter regions. - The ringdown model is based on the deformation of the ringdown model of the binary black - hole system in the EOB method, and uses the multiplicative hypothesis of multipole modes to describe the ringdown signal: \[ \bar{h}_{\ell m}(\tau) = e^{\sigma_1 \tau + i \phi_0} h(\tau) \equiv A_{\bar{h}}(\tau) e^{i \phi_{\bar{h}}(\tau)} \] - Among them, \( \sigma_1 = \alpha_{\ell m 1} + i \omega_{\ell m 1} \) is the complex frequency of the fundamental quasi - normal mode, and \( \tau = (t - t_0) / M_{\bullet} \) is the dimensionless time parameter. 3. **Inspiral - merger - ringdown waveform**: - The Inspiral - merger - ringdown (IMR) waveform construction method of TEOBResumS is similar to that of the binary black - hole system, but uses