Subtraction of infrared divergences in light-quark QCD sum rules

In QCD sum rules for light-quark systems, infrared (IR) divergences can appear in the Wilson coefficients of certain condensates. These divergences manifest explicitly in the coordinate-space expressions of the light-quark propagators. We propose an improved method to eliminate these IR divergences at the propagator level and present a subtraction formula that implements this procedure. Compared to the existing methods that rely on the mixing between quark and gluon condensates of the same dimension to eliminate IR divergences, this method is more intuitive and easier to apply in practical QCD sum rule calculations.