Algorithm for the simulation of many-electron systems at low temperatures

Abstract
We present a new algorithm for the numerical simulation of many-electron systems, which is particularly effective for the study of two-dimensional models at low temperatures. We present results for the two-dimensional Hubbard model on lattices as large as 8×8 and inverse temperatures as large as β=10. As the temperature is lowered, we find a continuous rise in the d-wave pairing susceptibility, indicating the possibility of a phase transition to a superconducting state.