  We discuss the non-zero baryon number formulation of QCD in the quenched
limit at finite temperature. This describes the thermodynamics of gluons in the
background of static quark sources. Although a sign problem remains in this
theory, our simulation results show that it can be handled quite well
numerically. The transition region gets shifted to smaller temperatures and the
transition region broadens with increasing baryon number. Although the action
is in our formulation explicitly Z(3) symmetric the Polyakov loop expectation
value becomes non-zero already in the low temperature phase and the heavy quark
potential gets screened at non-vanishing number density already this phase.
