On master integrals for two loop Bhabha scattering
  All scalar master integrals (MIs) for massive 2-loop QED Bhabha scattering are identified. The 2- and 3-point MIs have been calculated in terms of harmonic polylogarithms with the differential equation method. The calculation of 4-point MIs is underway. We sketch some alternative methods which help to solve (mainly) singularities of some MIs. 
M. Czakon, J. Gluza, T. Riemann