Cluster expansion is a basic but one of the few universal techniques in mathematical statistical physics. It is usually applied in the high-temperature regime where one usually obtains a unique Gibbs measure and exponential decay of correlation function. In this talk, we demonstrate the cluster expansion in the thermodynamic limit of the Yukawa gas. For the 1-charge component case, a rudimentary expansion is available. For the 2-charge component case, a clever interpolation is introduced, which hints the usefulness of multi-scale analysis.