The matching coefficients for the four-quark operators in NRQCD (NRQED) are
calculated at one loop using dimensional regularization for ultraviolet an
d infrared divergences. The matching for the electromagnetic current follow
s easily from our results. Both the unequal and equal mass cases are consid
ered. The role played by the Coulomb infrared singularities is explained in
detail. [S0556-2821(98)03621-2].