Local Coarse-grained Approach to Fermion Sign Problem in Path-Integral Monte Carlo Simulations

Siu A. Chin

Ethan Ethan, Jonathan L. DuBois, David M. Ceperley

Deyan Sun

An approximate treatment of exchange in finite-temperature path integral Monte Carlo simulations for fermions has been proposed. In this method, some of the fine details of density matrix due to permutations have been smoothed over or averaged out by using the coarse-grained approximation. The practical usefulness of the method is tested for interacting fermions in a three dimensional harmonic well. The results show that, the present method not only reduces the sign fluctuation of the density matrix, but also avoid the fermion system collapsing into boson system at low temperatures. The method is substantiated to be exact when applied to free particles.

Siu A. Chin

It is widely known that there is no sign problem in path integral Monte Carlo (PIMC) simulations of fermions in one dimension. As far as the author is aware, there is no direct proof of this in the literature. This work shows that the sign of the N-fermion antisymmetric free propagator is given by the product of all possible pairs of particle separations, or relative displacements. For a nonvanishing closed-loop product of such propagators, as required by PIMC, all relative displacements from adjacent propagators are paired into perfect squares, and therefore the loop product must be positive, but only in one dimension. By comparison, permutation sampling, which does not evaluate the determinant of the antisymmetric propagator exactly, remains plagued by a low-level sign problem, even in one dimension.

Saad A. Khairallah, Erik W. Draeger, John B. Shumway

Chi H. Mak, Reinhold Egger, H. Weber-Gottschick

A general algorithm toward the solution of the fermion sign problem in finite-temperature quantum Monte Carlo simulations has been formulated for discretized fermion path integrals with nearest-neighbor interactions in the Trotter direction. This multilevel approach systematically implements a simple blocking strategy in a recursive manner to synthesize the sign cancellations among different fermionic paths throughout the whole configuration space. The practical usefulness of the method is demonstrated for interacting electrons in a quantum dot.