In the post-Minkowskian limit approximation, we study gravitational wave solutions for general fourthorder theories of gravity. Specifically, we consider a Lagrangian with a generic function of curvature invariants f (R, R αβ R αβ , R αβγ δ R αβγ δ ). It is well known that when dealing with General Relativity such an approach provides massless spin-two waves as propagating degree of freedom of the gravitational field while this theory implies other additional propagating modes in the gravity spectra. We show that, in general, fourth order gravity, besides the standard massless graviton is characterized by two further massive modes with a finite-distance interaction. We find out the most general gravitational wave solutions in terms of Green functions in vacuum and in presence of matter sources. If an electromagnetic source is chosen, only the modes induced by R αβ R αβ are present, otherwise, for any f (R) gravity model, we have the complete analogy with tensor modes of General Relativity. Polarizations and helicity states are classified in the hypothesis of plane wave. B S. Capozziello capozziello@na.infn.it A. Stabile