Recently exploratory studies were performed on the possibility of constraining the neutron star equation of state (EOS) using signals from coalescing binary neutron stars, or neutron star-black hole systems, as they will be seen in upcoming advanced
gravitational wave detectors such as Advanced LIGO and Advanced Virgo. In particular, it was estimated to what extent the combined information from multiple detections would enable one to distinguish between different equations of state through hypothesis ranking or parameter estimation. Under the assumption of zero neutron star spins both in signals and in template waveforms and considering tidal effects to 1 post-Newtonian (1PN) order, it was found that O(20) sources would suffice to distinguish between a hard, moderate, and soft equation of state. Here we revisit these results, this time including neutron star tidal effects to the highest order currently known, termination of gravitational waveforms at the contact frequency, neutron star spins, and the resulting quadrupole-monopole interaction. We also take the masses of neutron stars in simulated sources to be distributed according to a relatively strongly peaked Gaussian, as hinted at by observations, but without assuming that the data analyst will necessarily have accurate knowledge of this distribution for use as a mass prior. We find that especially the effect of the latter is dramatic, necessitating many more detections to distinguish between different EOS and causing systematic biases in parameter estimation, on top of biases due to imperfect understanding of the signal model pointed out in earlier work. This would get mitigated if reliable prior information about the mass distribution could be folded into the analyses.
The direct detection of gravitational waves with upcoming second-generation gravitational wave detectors such as Advanced LIGO and Virgo will allow us to probe the genuinely strong-field dynamics of general relativity (GR) for the first time. We pres
ent a data analysis pipeline called TIGER (Test Infrastructure for GEneral Relativity), which is designed to utilize detections of compact binary coalescences to test GR in this regime. TIGER is a model-independent test of GR itself, in that it is not necessary to compare with any specific alternative theory. It performs Bayesian inference on two hypotheses: the GR hypothesis $mathcal{H}_{rm GR}$, and $mathcal{H}_{rm modGR}$, which states that one or more of the post-Newtonian coefficients in the waveform are not as predicted by GR. By the use of multiple sub-hypotheses of $mathcal{H}_{rm modGR}$, in each of which a different number of parameterized deformations of the GR phase are allowed, an arbitrarily large number of testing parameters can be used without having to worry about a model being insufficiently parsimonious if the true number of extra parameters is in fact small. TIGER is well-suited to the regime where most sources have low signal-to-noise ratios, again through the use of these sub-hypotheses. Information from multiple sources can trivially be combined, leading to a stronger test. We focus on binary neutron star coalescences, for which sufficiently accurate waveform models are available that can be generated fast enough on a computer to be fit for use in Bayesian inference. We show that the pipeline is robust against a number of fundamental, astrophysical, and instrumental effects, such as differences between waveform approximants, a limited number of post-Newtonian phase contributions being known, the effects of neutron star spins and tidal deformability on the orbital motion, and instrumental calibration errors.
