  In this paper we consider the matching coefficients up to two loops between
Quantum Chromodynamics (QCD) and Non-Relativistic QCD (NRQCD) for the vector,
axial-vector, scalar and pseudo-scalar currents. The structure of the effective
theory is discussed and analytical results are presented. Particular emphasis
is put on the singlet diagrams.
