Gravitational wave constraints on Einstein-\ae{}ther theory with LIGO/Virgo data

Lorentz symmetry is a fundamental property of Einstein's theory of general relativity that one may wish to test with gravitational wave observations. Einstein-aether theory is a model that introduces Lorentz-symmetry breaking in the gravitational sector through an aether vector field, while still leading to second-order field equations. This well-posed theory passes particle physics constraints because it modifies directly only the gravitational sector, yet it predicts deviations in the inspiral and coalescence of compact objects. We here, for the first time, put this theory to the test by comparing its gravitational wave predictions directly against LIGO/Virgo gravitational wave data. We construct a waveform model for Einstein-aether theory, EA_IMRPhenomD_NRT, through modifications of the general relativity IMRPhenomD_NRTidalv2 model (used by the LIGO/VIRGO collaboration). This model constructs a reponse function that not only contains the transverse-traceless polarization, but also additional Einstein-aether (scalar and vectorial) polarizations simultaneously. We then use the many current constraints on the theory to construct non-trivial priors for the Einstein-aether coupling constants. After testing the waveform model, we conduct parameter estimation studies on two gravitational wave events: GW170817 and GW190425. We find that these data are not sufficiently informative to place constraints on the theory that are stronger than current bounds from binary pulsar, solar system and cosmological observations. This is because, although Einstein-aether modifications include additional polarizations and have been computed beyond leading post-Newtonian order, these modifications are dominated by (already-constrained) dipole effects. These difficulties make it unclear whether future gravitational wave observations will be able to improve on current constraints on Einstein-aether theory.