  By decomposing the distribution functions and color field to regular and
fluctuation parts, the solution of the semi-classical kinetic equations of
quark-gluon plasma is analyzed. Through expanding the kinetic equations of the
fluctuation parts to third order, the nonlinear permittivity including the
self-interaction of gauge field is obtained and a rough numerical estimate is
given out for the important $\vk =0$ modes of the pure gluon plasma.
