Abstract
Using a reaction operator approach, we derive the multiple-scattering induced gluon number distribution function to all orders in powers of opacity at finite temperature. The detailed balance effect is analyzed by taking into account both gluon emission and absorption in a thermal medium. We also calculate virtual corrections and show that the infrared divergence cancels out in the gluon distribution function at finite temperature.