  We introduce a general Monte Carlo scheme for achieving atomistic simulations
with monoelectronic Hamiltonians including the thermalization of both nuclear
and electronic degrees of freedom. The kinetic Monte Carlo algorithm is used to
obtain the exact occupation numbers of the electronic levels at canonical
equilibrium, and comparison is made with Fermi-Dirac statistics in infinite and
finite systems. The effects of a nonzero electronic temperature on the
thermodynamic properties of liquid silver and sodium clusters are presented.
