Systematic effects from black hole-neutron star waveform model uncertainties on the neutron star equation of state

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

Systematic effects from black hole-neutron star waveform model uncertainties on the neutron star equation of state

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 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. As expected, these errors have negligible effect on detectability. Mass and spin estimates are moderately affected by systematic errors for near equal-mass binaries, while the recovered masses can be inaccurate at higher mass ratios. Large uncertainties are also found in the tidal deformability Λ , due to differences in PN base models used in hybridization, numerical relativity NR errors, and inherent limitations of the hybridization method. 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

Black hole-neutron star mergers using a survey of finite-temperature equations of state

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 semianalytic formulas. However, most of the simulations on which these formulas 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 formulas 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

Systematic effects from black hole-neutron star waveform model uncertainties on the neutron star equation of state

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 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. As expected, these errors have negligible effect on detectability. Mass and spin estimates are moderately affected by systematic errors for near equal-mass binaries, while the recovered masses can be inaccurate at higher mass ratios. Large uncertainties are also found in the tidal deformability Λ , due to differences in PN base models used in hybridization, numerical relativity NR errors, and inherent limitations of the hybridization method. 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

Survey of finite-temperature, nuclear-theory based Equations of State in Black Hole-Neutron Star Mergers

Neutron star-black hole binaries are one of the primary sources ofgravitational waves.These systems can also produce high-powered,bright electromagnetic counterparts via short-duration gammaray bursts and kilonovae, the latter of which is powered by the productionof heavy r-process elements.We consider systems where we assume the same initial black hole massand spin for all simulations, varying the equation of state of the neutronstar companion.We use three finite-temperature, composition-dependent, nuclear-theory basedequations of state (SFHo, DD2, FSU2.1) and assume neutron star massesin the range 1.2 - 1.4 M_\\odot.We show that ejecta masses mostly agree with predictions fit from simplerequations of state, although not all, while the ejecta velocities do agree with theupdated fitting-formula.We also determine that the dynamics of bound matter may be of particular future interest,where the bound tail material (fallback) and the early-stage circularizationof fluid near the horizon (protodisk) admixture is highly energetic andoptically bright in neutrinos.This distinction could be important in understanding the origins ofgamma ray bursts