We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the signal in noisy detector data with theoretical templates derived from some combination of post-Newtonian (PN) approximants, effective one-body (EOB) models and %analytic fits to numerical relativity (NR) simulations. The accuracy of these templates is limited by errors in the NR simulations, by the approximate nature of the PN/EOB waveforms, and by the hybridization procedure used to combine them. In this paper, we estimate the impact of these errors by constructing and comparing a set of PN-NR hybrid waveforms, for the first time with NR waveforms from two different codes, namely, SpEC and SACRA, for such systems. We then attempt to recover the parameters of the binary using two non-precessing template approximants. We find that systematic errors are too large for tidal effects to be accurately characterized for any realistic NS equation of state model. We conclude that NSBH waveform models must be significantly improved if they are to be useful for the extraction of NS equation of state information or even for distinguishing NSBH systems from binary black holes.

Each of the potential signals from a black hole-neutron star merger should contain an imprint of the neutron star equation of state: gravitational waves via its effect on tidal disruption, the kilonova via its effect on the ejecta, and the gamma ray burst via its effect on the remnant disk. These effects have been studied by numerical simulations and quantified by semi-analytic formulae. However, most of the simulations on which these formulae are based use equations of state without finite temperature and composition-dependent nuclear physics. In this paper, we simulate black hole-neutron star mergers varying both the neutron star mass and the equation of state, using three finite-temperature nuclear models of varying stiffness. Our simulations largely vindicate formulae for ejecta properties but do not find the expected dependence of disk mass on neutron star compaction. We track the early evolution of the accretion disk, largely driven by shocking and fallback inflow, and do find notable equation of state effects on the structure of this early-time, neutrino-bright disk.

Neutron star-black hole binaries are among the strongest sources of gravitational waves detectable by current observatories. They can also power bright electromagnetic signals (gamma-ray bursts, kilonovae), and may be a significant source of production of r-process nuclei. A misalignment of the black hole spin with respect to the orbital angular momentum leads to precession of that spin and of the orbital plane, and has a significant effect on the properties of the post-merger remnant and of the material ejected by the merger. We present a first set of simulations of precessing neutron star-black hole mergers using a hot, composition dependent, nuclear-theory based equation of state (DD2). We show that the mass of the remnant and of the dynamical ejecta are broadly consistent with the result of simulations using simpler equations of state, while differences arise when considering the dynamics of the merger and the velocity of the ejecta. We show that the latter can easily be understood from assumptions about the composition of low-density, cold material in the different equations of state, and propose an updated estimate for the ejecta velocity which takes those effects into account. We also present an updated mesh-refinement algorithm which allows us to improve the numerical resolution used to evolve neutron star-black hole mergers.