Improving the NRTidal model for binary neutron star systems

Accurate and fast gravitational waveform (GW) models are essential to extract information about the properties of compact binary systems that generate GWs. Building on previous work, we present an extension of the NRTidal model for binary neutron star (BNS) waveforms. The upgrades are (i) a new closed-form expression for the tidal contribution to the GW phase which includes further analytical knowledge and is calibrated to more accurate numerical relativity data than previously available; (ii) a tidal correction to the GW amplitude; and (iii) an extension of the spin-sector incorporating equation-of-state-dependent finite size effects at quadrupolar and octupolar order; these appear in the spin-spin tail terms and cubic-in-spin terms, both at 3.5 PN. We add the new description to the precessing binary black hole waveform model IMRPhenomPv2 to obtain a frequency-domain precessing binary neutron star model. In addition, we extend the SEOBNRv4_ROM and IMRPhenomD aligned-spin binary black hole waveform models with the improved tidal phase corrections. Focusing on the new IMRPhenomPv2_NRTidalv2 approximant, we test the model by comparing with numerical relativity waveforms as well as hybrid waveforms combining tidal effective-one-body and numerical relativity data. We also check consistency against a tidal effective-one-body model across large regions of the BNS parameter space.