Scalar diagrammatic rules for Born amplitudes in QCD
  We show that all Born amplitudes in QCD can be calculated from scalar propagators and a set of three- and four-valent vertices. In particular, our approach includes amplitudes with any number of quark pairs. The quarks may be massless or massive. The proof of the formalism is given entirely within quantum field theory. 
Christian Schwinn and Stefan Weinzierl