Leveraging NRTidalv3 to develop gravitational waveform models with higher-order modes for binary neutron star systems

Accurate and reliable gravitational waveform models are crucial in determining the properties of compact binary mergers. In particular, next-generation gravitational-wave detectors will require more accurate waveforms to avoid biases in the analysis. In this work, we extend the recent model to account for higher-mode corrections in the tidal phase contributions for binary neutron star systems. The higher-mode, multipolar model is then attached to several binary-black-hole baselines, such as the phenomenological and models, and the effective-one-body-based model . We test the performance and validity of the newly developed models by comparing them with numerical-relativity simulations and other tidal models. Finally, we employ them in parameter estimation analyses on simulated signals from both comparable-mass and high-mass-ratio systems, as well as on the gravitational-wave event GW170817, for which we find consistent results with respect to previous analyses.