Binary Neutron Star Mergers: Multi-Messenger Systematics and Prospects with Next-Generation Facilities

Abstract: Multi-messenger astronomy was galvanized by the detection of gravitational waves (GWs) from the binary neutron star (BNS) merger GW170817 and electromagnetic (EM) emission from the subsequent kilonova and short gamma ray burst. Maximizing multi-messenger constraints on these systems requires combining models of the progenitors and products of BNS mergers within a single framework in anticipation of future GW and EM detectors. We demonstrate how combining models of different aspects of the BNS progenitor-merger-remnant system into a full multi-messenger modeling pipeline reveals insight into modeling challenges that will need to be addressed in the coming decade of multi-messenger astronomy. Motivated by GW170817, this combined model relates the progenitor astrophysics of a BNS population with their GW observability and localizability, kilonova light curves, gamma-ray burst afterglow flux, and kilonova remnant evolution. We find joint correlations between the GW and EM observables that depend on a complicated interplay between modeling assumptions and theoretical uncertainties. Next generation detectors will generically provide multi-messenger constraints on BNSs with median GW network signal-to-noise ratio $\approx 10$, median of the 90$^{\rm th}$ percentile sky area $\approx 10$ sq. deg., and kilonova $i$-band apparent magnitudes ranging from $\approx 33$ to $23$. We find no more than 4\% of the BNSs are simultaneously detectable by a network of two Cosmic Explorers and one Einstein Telescope and by the Roman and Rubin telescopes for 55 and 180 sec. exposures in a $K$-like band and $i$ and $g$ bands, respectively. We discuss key modeling assumptions that will be transmitted as multi-messenger systematics in the analysis of future datasets, such as uncertain astrophysical processes of binary populations and the nuclear equation of state for tests of nuclear physics and dark matter.