We discuss the fermion sign problem and, by examining a very general Hubbard-Stratonovich transformation, argue that the sign problem cannot be solved with such methods. We propose a different kind of transformation which, while not solving the sign problem, shows more detailed information about the system. With our transformation it is trivial to tell which auxiliary field configurations give a positive sign and which give a negative sign. We then discuss briefly various properties of this transformation and construct an algorithm which with one simulation gives results for a whole range of particle densities and Hubbard U values, positive and negative. Our approach is in excellent agreement with exact calculations.
Auxiliary field quantum Monte Carlo simulations of interacting fermions require sampling over a Hubbard-Stratonovich field
h introduced to decouple the interactions. The weight for a given configuration involves the products of the determinant of matrices ${M}_{\ensuremath{\sigma}}(h)$, where $\ensuremath{\sigma}$ labels the species, and hence is typically not positive definite. Indeed, the average sign $\ensuremath{\langle}\mathcal{S}\ensuremath{\rangle}$ of the determinants goes to zero exponentially with increasing spatial size and decreasing temperature for most Hamiltonians of interest. This statement, however, does not explicitly separate two possible origins for the vanishing of $\ensuremath{\langle}\mathcal{S}\ensuremath{\rangle}$. Does $\ensuremath{\langle}\mathcal{S}\ensuremath{\rangle}\ensuremath{\rightarrow}0$ because randomly chosen field configurations have
det[M(h)]<0, or does the sign problem arise because the specific subset of configurations chosen by the weighting function have a greater preponderance of negative values? In the latter case, the process of weighting the configurations with
∣det[M(h)]∣ might steer the simulation to a region of configuration space of
h where positive and negative determinants are equally likely, even though randomly chosen
h would preferentially have determinants with a single dominant sign. In this paper, we address the relative importance of these two mechanisms for the vanishing of $\ensuremath{\langle}\mathcal{S}\ensuremath{\rangle}$ in quantum simulations.
The auxiliary-field quantum Monte Carlo method is reviewed. The Hubbard-Stratonovich transformation converts an interacting Hamiltonian into a non-interacting Hamiltonian in a time-dependent stochastic field, allowing calculation of the resulting functional integral by Monte Carlo methods. The method is presented in a sufficiently general form to be applicable to any Hamiltonian with one- and two-body terms, with special reference to the Heisenberg model and one- and many-band Hubbard models. Many physical correlation functions can be related to correlation functions of the auxiliary field; general results are given here. Issues relating to the choice of auxiliary fields are addressed; operator product identities change the relative dimensionalities of the attractive and repulsive parts of the interaction. Frequently the integrand is not positive-definite, rendering numerical evaluation unstable. If the auxiliary field violates time-reversal invariance, the integrand is complex and this sign problem becomes a phase problem. The origin of this sign or phase is examined from a number of geometrical and other viewpoints and illustrated by simple examples: the phase problem by the spin (1/2) Heisenberg model, and the sign problem by the attractive SU(N) Hubbard model on a triangular molecule with negative hopping integrals. In the latter case, widely studied in the Jahn-Teller literature, the sign is due neither to fermions nor spin, but to frustration. This system is used to illustrate a number of suggested interpretations of the sign problem.
Auxiliary field Quantum Monte Carlo simulations of interacting fermions require sampling over a Hubbard-Stratonovich field
x introduced to decouple the interactions. The weight for a given configuration involves the products of the determinant of matrices
Mσ(x), where
σ labels the species, and hence is typically not positive definite. Indeed, the average sign
⟨S⟩ of the determinants goes to zero exponentially with increasing spatial size and decreasing temperature for most Hamiltonians of interest. This statement, however, does not explicitly separate two possible origins for the vanishing of
⟨S⟩. Does
⟨S⟩→0 because {\it randomly} chosen field configurations have
det(M(x))<0, or does the `sign problem' arise as a result of {\it importance sampling} in the simulations? In the latter case, the process of weighting the configurations with
∣det(M(x))∣ might steer the simulation to a region of configuration space of
x where positive and negative determinants are equally likely, even though randomly chosen
x would preferentially have determinants with a single dominant sign. In this paper we address the relative importance of these two mechanisms for the vanishing of
⟨S⟩ in quantum simulations.
The fermion sign problem, the biggest obstacle in quantum Monte Carlo calculations, is completely solved in this paper. Here, we find a strategy, in which the contribution from those negative-weighted paths is thoroughly cancelled or replaced by some positive-weighted paths. The crucial point lies on the Feynman path integral formula proposed in our group, which allows us to deeply analyze the Boltzmann weight of each path. Through mathematical proof, we demonstrate that physical quantities can be exactly calculated within a specific kind of paths, which have positive the Boltzmann weight. With this finding, a new Monte Carlo method is proposed, in which the fermion sign problem is absent. As an example, the current method is applied to the two-dimensional Hubbard model, and the results do manifest the correctness.
The minus sign problem, as a problem of the non-positive-definite measure of the path integrals, is examined in the Monte Carlo simulation of the Hubbard model. A new correction method in terms of the ratio of two measured quantities is proposed to calculate physical properties. The momentum distribution and the spin correlation as well as the energy in the ground state are tested against the exactly diagonalized results. The results show that this ratio correction method is a promising prescription to reduce the minus sign problem.
We consider a class of Hubbard-Stratonovich transformations suitable for treating Hubbard interactions in the context of quantum Monte Carlo simulations. A tunable parameter p allows us to continuously vary from a discrete Ising auxiliary field (p=∞) to a compact auxiliary field that couples to electrons sinusoidally (p=0). In tests on the single-band square and triangular Hubbard models, we find that the severity of the sign problem decreases systematically with increasing p. Selecting p finite, however, enables continuous sampling methods such as the Langevin or Hamiltonian Monte Carlo methods. We explore the tradeoffs between various simulation methods through numerical benchmarks.