Given the recent discovery of gravitational waves there is a pressing need for high accuracy templates in order to, not only extract precise values of black hole masses and spins, but also to learn about the equation of state of neutron stars. While late stages of the inspiral necessitates numerical simulations, the early stages can be calculated analytically in an expansion in the relative velocity (the Post-Newtonian approximation). This situation is reminiscent of quarkonia as both General Relativity and QCD have non-linearities controlled by the relative velocity of the bound state constituents. In this talk I will discuss how state of the art predictions have been generated using NRGR (Non-Relativistic General Relativity) which draws on ideas from NRQCD. I will also show how modern scattering amplitude techniques can be used to derive higher order Post-Newtonian potentials, which via the double copy relation, generates solutions to Einsteins equations directly from QCD. Finally I will discuss how the remarkable claim that black holes have vanishing tidal Love numbers implies, from the field theory point of view, that black holes are finely tuned objects. Remarkably, this claim can be tested in gravitational wave experiments.