Abstract: We present the computation of the two-loop form factors for diphoton production in the quark annihilation channel. These quantities are relevant for the NNLO QCD corrections to diphoton production at LHC recently presented in [1]. The computation is performed retaining full dependence on the mass of the heavy quark in the loops. The master integrals are evaluated by means of differential equations which are solved exploiting the generalised power series technique.