  We propose a beyond mean field approach to evaluate Yang-Mills thermodynamics
from the partition function with n-body gluon contribution, in the presence of
a uniform background Polyakov field. Using a path integral based formalism, we
obtain, unlike the previous mean field studies within this model framework,
physically consistent results with good agreement to the lattice data
throughout the temperature range.
