Search for: All records

The post-Newtonian formalism plays an integral role in the models used to extract information from gravitational wave data, but models that incorporate this formalism are inherently approximations. Disagreement between an approximate model and nature will produce mismodeling biases in the parameters inferred from data, introducing systematic error. We here carry out a proof-of- principle study of such systematic error by considering signals produced by quasi-circular, inspiraling black hole binaries through an injection and recovery campaign. In particular, we study how un- known, but calibrated, higher-order post-Newtonian corrections to the gravitational wave phase impact systematic error in recovered parameters. As a first study, we produce injected data of non-spinning binaries as detected by a current, second-generation network of ground-based observatories and recover them with models of varying PN order in the phase. We find that the truncation of higher order (>3.5) post-Newtonian corrections to the phase can produce significant systematic error even at signal-to-noise ratios of current detector networks. We propose a method to mitigate systematic error by marginalizing over our ignorance in the waveform through the inclusion of higher-order post-Newtonian coefficients as new model parameters. We show that this method can reduce systematic error greatly at the cost of increasing statistical error. 

Abstract Hundreds of millions of supermassive black hole binaries are expected to contribute to the gravitational-wave signal in the nanohertz frequency band. Their signal is often approximated either as an isotropic Gaussian stochastic background with a power-law spectrum or as an individual source corresponding to the brightest binary. In reality, the signal is best described as a combination of a stochastic background and a few of the brightest binaries modeled individually. We present a method that uses this approach to efficiently create realistic pulsar timing array data sets using synthetic catalogs of binaries based on the Illustris cosmological hydrodynamic simulation. We explore three different properties of such realistic backgrounds that could help distinguish them from those formed in the early universe: (i) their characteristic strain spectrum, (ii) their statistical isotropy, and (iii) the variance of their spatial correlations. We also investigate how the presence of confusion noise from a stochastic background affects detection prospects of individual binaries. We calculate signal-to-noise ratios of the brightest binaries in different realizations for a simulated pulsar timing array based on the NANOGrav 12.5 yr data set extended to a time span of 15 yr. We find that ∼6% of the realizations produce systems with signal-to-noise ratios larger than 5, suggesting that individual systems might soon be detected (the fraction increases to ∼41% at 20 yr). These can be taken as a pessimistic prediction for the upcoming NANOGrav 15 yr data set, since it does not include the effect of potentially improved timing solutions and newly added pulsars.