We study electron excitations at finite temperature by giving the electron positions and sizes initial random velocities, and propagating the resulting dynamics at constant energy and volume. We use as a test case a NaCl lattice of 64 electrons with = 2 bohr. At low temperatures, the electrons make small excursions about their equilibrium positions, but at higher temperatures they mix more freely (Figure 6.8).
We can measure the heat capacity of the electron liquid by plotting the total energy as a function of temperature (T = 50 to 500 K, Figure 6.10). In a classical solid, the heat capacity at low temperatures is given by the Dulong and Petit expression , since there is an equipartition of energy among all the degrees of freedom in the solid. In metals, the heat capacity at temperatures below the Fermi temperature ( = 140,000 K for = 2 bohr) scales as ; the heat capacity is much lower than would be expected from a classical solid, because only states near the Fermi level can be excited  (Figure 6.9).
The eFF electron gas has a heat capacity that matches classical, not quantum, statistics (Figure 6.10), suggesting that all of the available modes are being excited uniformly, a contention further supported by the spectrum of phonon excitations (Figure 6.11) derived by computing the Fourier transform of the velocity autocorrelation function .
Since we are simulating electrons as classical particles interacting via effective potentials, it is not surprising that we reproduce classical and not quantum statistics. Other researchers who have created quasiclassical models of the uniform electron gas have reproduced the correct Fermi-Dirac distribution of momenta using momentum dependent potentials, which spread out the particles in momentum phase space. For example, in 1987, Dorso and Randrup  applied a Pauli potential of the form
In the future we will try adding a momentum-dependent repulsion function to eFF to reproduce Fermi-Dirac distributions of momenta. However, we will need to proceed carefully to ensure (1) that energy remains conserved, and that there is no energy loss via hysteresis effects and (2) that we still obtain correct dynamics when electrons are associated with nuclei.