Yiwen Huang, Carl-Johan Haster, Salvatore Vitale, Vijay Varma, Francois Foucart, Sylvia Biscoveanu

Gravitational waves emitted by neutron star black hole mergers encode key properties of neutron stars -- such as their size, maximum mass and spins -- and black holes. However, the presence of matter and the high mass ratio makes generating long and accurate waveforms from these systems hard to do with numerical relativity, and not much is known about systematic uncertainties due to waveform modeling. We simulate gravitational waves from neutron star black hole mergers by hybridizing numerical relativity waveforms produced with the SpEC code with recent numerical relativity surrogate waveforms. These signals are analyzed using a range of available waveform families, and statistical and systematic errors are reported. We find that at a network signal-to-noise ratio (SNR) of 30, statistical uncertainties are usually larger than systematic offsets, while at SNR of 70 the two become comparable. All neutron stars in our simulations are non-spinning, but in no case we can constrain the neutron star spin to be smaller than $\sim0.4$ (90% credible interval). The individual black hole and neutron star masses, as well as the mass ratios, are typically measured very precisely, though not always accurately at high SNR. At a SNR of 30 the neutron star tidal deformability can only be bound from above, while for louder sources it can be measured and constrained away from zero. Waveform families whose late inspiral has been tuned specifically for neutron star black hole signals yield measurements which are in tension with those obtained using waveform families tuned against binary neutron stars, even for mass ratios that could be relevant for both binary neutron stars and neutron star black holes mergers.