An algorithm for calculating two-loop propagator type Feynman diagrams with arbitrary masses and external momentum is proposed. Recurrence relations allowing to express any scalar integral in terms of basic integrals are given. A minimal set consisting of 15 essentially two-loop and 15 products of one-loop basic integrals is found. Tensor integrals and integrals with irreducible numerators are represented as a combination of scalar ones with a higher space-time dimension which are reduced to the basic set by using the generalized recurrence relations proposed in [1].

Recently a powerful duality between color and kinematics has been proposed for integrands of scattering amplitudes in quite general gauge theories. In this paper the duality proposal is extended to the more general class of gauge theory observables formed by form factors. After a discussion of the general setup the existence of the duality is verified in twoand three-loop examples in four dimensional maximally supersymmetric Yang-Mills theory which involve the stress energy tensor multiplet. In these cases the duality reproduces known results in a particularly transparent and uniform way. As a non-trivial application we obtain a very simple form of the integrand of the four-loop two-point (Sudakov) form factor which passes a large set of unitarity cut checks.

New types of relationships between Feynman integrals are presented. It is shown that Feynman integrals satisfy functional equations connecting integrals with different values of scalar invariants and masses. A method is proposed for obtaining such relations. The derivation of functional equations for one-loop propagator- and vertex - type integrals is given. It is shown that a propagator - type integral can be written as a sum of two integrals with modified scalar invariants and one propagator massless. The vertex - type integral can be written as a sum over vertex integrals with all but one propagator massless and one external momenta squared equal to zero. It is demonstrated that the functional equations can be used for the analytic continuation of Feynman integrals to different kinematic domains.Comment: LaTeX, 10 pages, several typos are corrected. To be published in Phys.Lett.

A recurrence relation between equal mass two-loop sunrise diagrams differing in dimensionality by 2 is derived and it's solution in terms of Gauss' 2F1 and Appell's F_2 hypergeometric functions is presented. For arbitrary space-time dimension d the imaginary part of the diagram on the cut is found to be the 2F1 hypergeometric function with argument proportional to the maximum of the Kibble cubic form. The analytic expression for the threshold value of the diagram in terms of the hypergeometric function 3F2 of argument -1/3 is given.Comment: 10 page

We present a new method for the momentum expansion of Feynman integrals with arbitrary masses and any number of loops and external momenta. By using the parametric representation we derive a generating function for the coefficients of the small momentum expansion of an arbitrary diagram. The method is applicable for the expansion w.r.t. all or a subset of external momenta. The coefficients of the expansion are obtained by applying a differential operator to a given integral with shifted value of the space-time dimension d and the expansion momenta set equal to zero. Integrals with changed d are evaluated by using the generalized recurrence relations proposed in [1]. We show how the method works for one-and two-loop integrals. It is also illustrated that our method is simpler and more efficient than others.