Ben Farr

and 16 more

INTRODUCTION As we enter the advanced-detector era of ground-based gravitational-wave (GW) astronomy, it is critical that we understand the abilities and limitations of the analyses we are prepared to conduct. Of the many predicted sources of GWs, binary neutron-star (BNS) coalescences are paramount; their progenitors have been directly observed , and the advanced detectors will be sensitive to their GW emission up to ∼400 Mpc away . When analyzing a GW signal from a circularized compact binary merger, strong degeneracies exist between parameters describing the binary (e.g., distance and inclination). To properly estimate any particular parameter(s) of interest, the marginal distribution is estimated by integrating the joint posterior probability density function (PDF) over all other parameters. In this work, we sample the posterior PDF using software implemented in the LALINFERENCE library . Specifically we use results from LALINFERNCE_NEST , a nest sampling algorithm , and LALINFERENCE_MCMC , a Markov-chain Monte Carlo algorithm \citep[chapter 12]{Gregory2005}. Previous studies of BNS signals have largely assessed parameter constraints assuming negligible neutron-star (NS) spin, restricting models to nine parameters. This simplification has largely been due to computational constraints, but the slow spin of NSs in short-period BNS systems observed to date \citep[e.g.,][]{Mandel_2010} has also been used as justification. However, proper characterization of compact binary sources _must_ account for the possibility of non-negligible spin; otherwise parameter estimates will be biased . This bias can potentially lead to incorrect conclusions about source properties and even misidentification of source classes. Numerous studies have looked at the BNS parameter estimation abilities of ground-based GW detectors such as the Advanced Laser Interferometer Gravitational-Wave Observatory \citep[aLIGO;][]{Aasi_2015} and Advanced Virgo \citep[AdV;][]{Acernese_2014} detectors. assessed localization abilities on a simulated non-spinning BNS population. looked at several potential advanced-detector networks and quantified the parameter-estimation abilities of each network for a signal from a fiducial BNS with non-spinning NSs. demonstrated the ability to characterize signals from non-spinning BNS sources with waveform models for spinning sources using Bayesian stochastic samplers in the LALINFERENCE library . used approximate methods to quantify the degeneracy between spin and mass estimates, assuming the compact objects’ spins are aligned with the orbital angular momentum of the binary \citep[but see][]{Haster_2015}. simulated a collection of loud signals from non-spinning BNS sources in several mass bins and quantified parameter estimation capabilities in the advanced-detector era using non-spinning models. introduced precession from spin–orbit coupling and found that the additional richness encoded in the waveform could reduce the mass–spin degeneracy, helping BNSs to be distinguished from NS–black hole (BH) binaries. conducted a similar analysis of a large catalog of sources and found that it is difficult to infer the presence of a mass gap between NSs and BHs , although, this may still be possible using a population of a few tens of detections . Finally, and the follow-on represent an (almost) complete end-to-end simulation of BNS detection and characterization during the first 1–2 years of the advanced-detector era. These studies simulated GWs from an astrophysically motivated BNS population, then detected and characterized sources using the search and follow-up tools that are used for LIGO–Virgo data analysis . The final stage of the analysis missing from these studies is the computationally expensive characterization of sources while accounting for the compact objects’ spins and their degeneracies with other parameters. The present work is the final step of BNS characterization for the simulations using waveforms that account for the effects of NS spin. We begin with a brief introduction to the source catalog used for this study and in section [sec:sources]. Then, in section [sec:spin] we describe the results of parameter estimation from a full analysis that includes spin. In section [sec:mass] we look at mass estimates in more detail and spin-magnitude estimates in section [sec:spin-magnitudes]. In section [sec:extrinsic] we consider the estimation of extrinsic parameters: sky position (section [sec:sky]) and distance (section [sec:distance]), which we do not expect to be significantly affected by the inclusion of spin in the analysis templates. We summarize our findings in section [sec:conclusions]. A comparison of computational costs for spinning and non-spinning parameter estimation is given in appendix [ap:CPU].