In this paper, we study gravitational waves generated by binary systems within an extension of General Relativity which is described by the addition of quadratic in curvature tensor terms to the Einstein–Hilbert action. Treating quadratic gravity as an effective theory valid in the low energy/curvature regime, we argue that reliable calculations can be performed in the early inspiral phase, and furthermore, no flux of additional massive waves can be detected. We then compute the Yukawa suppressed dipole-like corrections to the post-Newtonian (PN) expansion of the standard waveform. By confronting these theoretical calculations with available experimental data, we constrain both unknown parameters of quadratic gravity to be $$0 \le \gamma \, \lesssim 5.9\times 10^{76}$$
0
≤
γ
≲
5.9
×
10
76
, and $$-\frac{\gamma }{4} \le \beta \, \lesssim 9.8\times 10^{75} - \frac{\gamma }{4}$$
-
γ
4
≤
β
≲
9.8
×
10
75
-
γ
4
.