Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

A Monte Carlo study of the minus sign problem in the t-J model using an Intel iPSC/860 hypercube

Mitchel D. Kovarik, Ted Barnes

Published Jan 1, 1993• 0 citations
OpenAlex

Abstract

We describe a Monte Carlo simulation of the 2-dimensional t-J model on an Intel iPSC/860 hypercube. The problem studied is the determination of the dispersion relation of a dynamical hole in the t-J model of the high temperature superconductors. Since this problem involves the motion of many fermions in more than one spatial dimensions, it is representative of the class of systems that suffer from the ``minus sign problem`` of dynamical fermions which has made Monte Carlo simulation very difficult. We demonstrate that for small values of the hole hopping parameter one can extract the entire hole dispersion relation using the GRW Monte Carlo algorithm, which is a simulation of the Euclidean time Schroedinger equation, and present results on 4 {times} 4 and 6 {times} 6 lattices. We demonstrate that a qualitative picture at higher hopping parameters may be found by extrapolating weak hopping results where the minus sign problem is less severe. Generalization to physical hopping parameter values will only require use of an improved trial wavefunction for importance sampling.

Related papers

Simulating the Euclidean time Schroedinger equations using an Intel iPSC/860 hypercube: Application to the t-J model of high-{Tc} superconductivity

Mitchel D. Kovarik, Ted Barnes

We describe a Monte Carlo simulation of a dynamical fermion problem in two spatial dimensions on an Intel iPSC/860 hypercube. The problem studied is the determination of the dispersion relation of a dynamical hole in the t-J model of the high temperature superconductors. Since this problem involves the motion of many fermions in more than one spatial dimensions, it is representative of the class of systems that suffer from the ``minus sign problem`` of dynamical fermions which has made Monte Carlo simulation very difficult. We demonstrate that for small values of the hole hopping parameter one can extract the entire hole dispersion relation using the GRW Monte Carlo algorithm, which is a simulation of the Euclidean time Schroedinger equation, and present results on 4 {times} 4 and 6 {times} 6 lattices. Generalization to physical hopping parameter values wig only require use of an improved trial wavefunction for importance sampling.

•Oct 1, 1993•0 citations

Monte Carlo diagonalization of many-body problems: Application to fermion systems

Hans A. De Raedt, Wolfgang von der Linden

Existing quantum Monte Carlo algorithms suffer the so-called minus-sign problem. We propose a stochastic algorithm to compute ground-state properties of quantum lattice models that is free of minus-sign problems. Illustrative results for the two-dimensional Hubbard model are presented.

Physical review. B, Condensed matter•Apr 15, 1992•31 citations

MONTE CARLO DIAGONALIZATION OF VERY LARGE MATRICES: APPLICATION TO FERMION SYSTEMS

Hans A. De Raedt, Wolfgang von der Linden

All known Quantum-Monte-Carlo algorithms for fermions suffer from the so-called “minus-sign-problem” which is detrimental to the application of these simulation methods to fermion systems at very low temperatures and/or of very many lattice sites. We identify the origin of this fundamental problem, demonstrate that it is a very general feature, not necessarily related to the presence of fermionic degrees of freedom. We describe an novel algorithm which does not suffer from the minus-sign problem. Illustrative results for the two-dimensional Hubbard model are presented

International Journal of Modern Physics C•Feb 1, 1992•6 citations

Minus Sign Problem in the Monte Carlo Simulation of Lattice Fermion Systems

Nobuo Furukawa, Masatoshi Imada

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.

Journal of the Physical Society of Japan•Mar 15, 1991•32 citations

Perturbative solution of fermionic sign problem in quantum Monte Carlo computations

Sergei Iskakov, Mikhail I. Katsnelson, Alexander I. Lichtenstein

Abstract We have developed a strong-coupling perturbation scheme for a general doped Hubbard model around a particle-hole-symmetric reference system, which is free from the fermionic sign problem. Our approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in both continuous and discrete time versions for large periodic clusters in a fermionic bath. By considering the first-order perturbation in the shift of the chemical potential and the second-neighbor hopping, we are able to obtain an accurate electronic spectral function for a range of parameters that correspond to optimally doped cuprate systems at temperatures of up to T = 0.1 t , which are challenging to access using straightforward lattice QMC calculations. We also discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system with the interaction parameter U equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter $${t}^{{\prime} }$$ t ′ for high-temperature superconducting cuprates.

arXiv:2303.01607•npj Computational Materials•Feb 15, 2024•6 citations

Sign problem in Monte Carlo simulations of frustrated quantum spin systems

Patrik Henelius, Anders W. Sandvik

We discuss the sign problem arising in Monte Carlo simulations of frustrated quantum spin systems. We show that for a class of ``semifrustrated'' systems [Heisenberg models with ferromagnetic couplings ${J}_{z}(r)<0$ along the z axis and antiferromagnetic couplings ${J}_{\mathrm{xy}}(r)=\ensuremath{-}{J}_{z}(r)$ in the $\mathrm{xy}$ plane, for arbitrary distances $r],$ the sign problem present for algorithms operating in the z basis can be solved within a recent ``operator-loop'' formulation of the stochastic series expansion method [a cluster algorithm for sampling the diagonal matrix elements of the power series expansion of $\mathrm{exp}(\ensuremath{-}\ensuremath{\beta}H)$ to all orders]. The solution relies on the identification of operator loops which change the configuration sign when updated (``merons'') and is similar to the meron-cluster algorithm recently proposed by Chandrasekharan and Wiese for solving the sign problem for a class of fermion models [Phys. Rev. Lett. 83, 3116 (1999)]. Some important expectation values, e.g., the internal energy, can be evaluated in the subspace with no merons, where the weight function is positive definite. Calculations of other expectation values require sampling of configurations with only a small number of merons (typically zero or two), with an accompanying sign problem which is not serious. We also discuss problems which arise in applying the meron concept to more general quantum spin models with frustrated interactions.

arXiv:cond-mat/0001351•Physical review. B, Condensed matter•Jul 1, 2000•134 citations
We describe a Monte Carlo simulation of the 2-dimensional t-J model on an Intel iPSC/860 hypercube. The problem studied is the determination of the dispersion relation of a dynamical hole in the t-J model of the high temperature superconductors. Since this problem involves the motion of many fermions in more than one spatial dimensions, it is representative of the class of systems that suffer from the minus sign problem of dynamical fermions which has made Monte Carlo simulation very difficult. We demonstrate that for small values of the hole hopping parameter one can extract the entire hole dispersion relation using the GRW Monte Carlo algorithm, which is a simulation of the Euclidean time Schroedinger equation, and present results on 4 {times} 4 and 6 {times} 6 lattices. We demonstrate that a qualitative picture at higher hopping parameters may be found by extrapolating weak hopping results where the minus sign problem is less severe. Generalization to physical hopping parameter values will only require use of an improved trial wavefunction for importance sampling.
We describe a Monte Carlo simulation of a dynamical fermion problem in two spatial dimensions on an Intel iPSC/860 hypercube. The problem studied is the determination of the dispersion relation of a dynamical hole in the t-J model of the high temperature superconductors. Since this problem involves the motion of many fermions in more than one spatial dimensions, it is representative of the class of systems that suffer from the ``minus sign problem`` of dynamical fermions which has made Monte Carlo simulation very difficult. We demonstrate that for small values of the hole hopping parameter one can extract the entire hole dispersion relation using the GRW Monte Carlo algorithm, which is a simulation of the Euclidean time Schroedinger equation, and present results on 4 {times} 4 and 6 {times} 6 lattices. Generalization to physical hopping parameter values wig only require use of an improved trial wavefunction for importance sampling.
All known Quantum-Monte-Carlo algorithms for fermions suffer from the so-called “minus-sign-problem” which is detrimental to the application of these simulation methods to fermion systems at very low temperatures and/or of very many lattice sites. We identify the origin of this fundamental problem, demonstrate that it is a very general feature, not necessarily related to the presence of fermionic degrees of freedom. We describe an novel algorithm which does not suffer from the minus-sign problem. Illustrative results for the two-dimensional Hubbard model are presented
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.
Abstract We have developed a strong-coupling perturbation scheme for a general doped Hubbard model around a particle-hole-symmetric reference system, which is free from the fermionic sign problem. Our approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in both continuous and discrete time versions for large periodic clusters in a fermionic bath. By considering the first-order perturbation in the shift of the chemical potential and the second-neighbor hopping, we are able to obtain an accurate electronic spectral function for a range of parameters that correspond to optimally doped cuprate systems at temperatures of up to T = 0.1 t , which are challenging to access using straightforward lattice QMC calculations. We also discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system with the interaction parameter U equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter t′{t}^{{\prime} }t′ t ′ for high-temperature superconducting cuprates.
We discuss the sign problem arising in Monte Carlo simulations of frustrated quantum spin systems. We show that for a class of ``semifrustrated'' systems [Heisenberg models with ferromagnetic couplings Jz(r)<0{J}_{z}(r)<0Jz​(r)<0 along the z axis and antiferromagnetic couplings ${J}_{\mathrm{xy}}(r)=\ensuremath{-}{J}_{z}(r)$ in the xy\mathrm{xy}xy plane, for arbitrary distances r],r],r], the sign problem present for algorithms operating in the z basis can be solved within a recent ``operator-loop'' formulation of the stochastic series expansion method [a cluster algorithm for sampling the diagonal matrix elements of the power series expansion of $\mathrm{exp}(\ensuremath{-}\ensuremath{\beta}H)$ to all orders]. The solution relies on the identification of operator loops which change the configuration sign when updated (``merons'') and is similar to the meron-cluster algorithm recently proposed by Chandrasekharan and Wiese for solving the sign problem for a class of fermion models [Phys. Rev. Lett. 83, 3116 (1999)]. Some important expectation values, e.g., the internal energy, can be evaluated in the subspace with no merons, where the weight function is positive definite. Calculations of other expectation values require sampling of configurations with only a small number of merons (typically zero or two), with an accompanying sign problem which is not serious. We also discuss problems which arise in applying the meron concept to more general quantum spin models with frustrated interactions.