Algebraic-Numerical Evaluation of Feynman Diagrams: Two-Loop Self-Energies

A recently proposed scheme for numerical evaluation of Feynman diagrams is extended to cover all two-loop two-point functions with arbitrary internal and external masses. The adopted algorithm is a modification of the one proposed by F. V. Tkachov and it is based on the so-called generalized Bernstein functional relation. On-shell derivatives of self-energies are also considered and their infrared properties analyzed to prove that the method which is aimed to a numerical evaluation of massive diagrams can handle the infrared problem within the scheme of dimensional regularization. Particular care is devoted to study the general massive diagrams around their leading and non-leading Landau singularities.