Analytic results on the massive three-loop form factors: quarkonic contributions