We calculate paer Monte Carlo evaluation on an 8 3 × 3 lattice the energy density ϵ G of the gluon sector of QCD, including virtual quark loops up to the fourth power in the hopping parameter expansion. For light quarks of one flavour, we observe at T/ Λ L 95 ± 10 a rapid variation of ϵ G in T, accompanied by strong fluctuations from iteration to iteration. as clear signal of the deconfinement transition.