Constraining dark-sector effects using gravitational waves from compact binary inspirals

Gravitational wave observations are a powerful tool to constrain fundamental physics. This work considers dark matter that carries charge under a dark abelian massive vector field. If such dark matter is bound inside coalescing neutron stars, the presence of the new force will modify the total energy of the binary, and the emission of dark radiation modes will impact the rate of inspiral. At leading order, the dark matter corrections introduce a Yukawa term to the potential energy and a dipole radiation mode. We calculate modifications to the binary inspiral waveform to first post-Newtonian order, incorporate these corrections into a state-of-the-art waveform model for neutron star binaries, and use this model in a Bayesian parameter estimation analysis on real data collected by current, second-generation ground-based detectors. Through this study, we place the first robust constraints on the Yukawa interaction strength and dipole emission parameter. We find that the strongest constraints are obtained from analysis of the GW170817 event, for which the Yukawa interaction strength is constrained to be less than $O(10^{-2})$ and the dipole emission parameter is less than $O(10^{-4})$ at 90% credibility.