Jellium at finite temperature

Abstract

We adopt the fixed node restricted path integral Monte Carlo method within the ‘Worm algorithm’ to simulate Wigner's Jellium model at finite, non zero, temperatures using free-particle nodes of the density matrix. The new element is that we incorporate the Worm algorithm paradigm of Prokof'ev and Svistunov in the grand canonical ensemble in order to more efficiently handle the fermionic exchanges. We present results for the structure and thermodynamic properties of the ideal Fermi gas and three points for the interacting electron gas. We treat explicitly the case of the partially polarized electron gas.