  The contributions to the coefficient functions of the quark and the mixed
quark-gluon condensate to mesonic correlators are calculated for the first time
to all orders in the quark masses, and to lowest order in the strong coupling
constant. Existing results on the coefficient functions of the unit operator
and the gluon condensate are reviewed. The proper factorization of short- and
long-distance contributions in the operator product expansion is discussed in
detail. It is found that to accomplish this task rigorously the operator
product expansion has to be performed in terms of non-normal-ordered
condensates. The resulting coefficient functions are improved with the help of
the renormalization group. The scale invariant combination of dimension 5
operators, including mixing with the mass operator, which is needed for the
renormalization group improvement, is calculated in the leading order.
