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

arXiv:cond-mat/0611550 · Physical Review B

Quantum Monte Carlo scheme for frustrated Heisenberg antiferromagnets

Jacek Wojtkiewicz

Published May 15, 2007• 2 citations
arXivHTMLPDF

Abstract

When one tries to simulate quantum spin systems by the Monte Carlo method, often the ``minus-sign problem'' is encountered. In such a case, an application of probabilistic methods is not possible. In this paper, the method has been proposed which allows to avoid the minus-sign problem for certain class of frustrated Heisenberg models. The systems where this method is applicable are, for instance, the pyrochlore lattice and the ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model. The method works in singlet sector. It relies on the expression of wave functions in dimer (pseudo)basis and writing down the Hamiltonian as a sum over plaquettes. In such a formulation, matrix elements of the exponent of Hamiltonian are positive.

DOI
10.1103/physrevb.75.174421

Related papers

Green Function Monte Carlo with Stochastic Reconfiguration

Sandro Sorella

A new method for the stabilization of the sign problem in the Green function Monte Carlo technique is proposed. The method is devised for real lattice Hamiltonians and is based on an iterative ``stochastic reconfiguration'' scheme which introduces some bias but allows a stable simulation with constant sign. The systematic reduction of this bias is possible in principle. The method is applied to the frustrated ${J}_{1}{\ensuremath{-}J}_{2}$ Heisenberg model, and tested against exact diagonalization data. Evidence of a finite spin gap for ${J}_{2}{/J}_{1}>\ensuremath{\sim}0.4$ is found in the thermodynamic limit.

arXiv:cond-mat/9803107•Physical Review Letters•May 18, 1998•267 citations

Comment on “Quantum Monte Carlo scheme for frustrated Heisenberg antiferromagnets”

K. S. D. Beach, Matthieu Mambrini, Fabien Alet

Quantum Monte Carlo methods are sophisticated numerical techniques for simulating interacting quantum systems. In some cases, however, they suffer from the notorious ``sign problem'' and become too inefficient to be useful. A recent publication [J. Wojtkiewicz, Phys. Rev. B 75, 174421 (2007)] claims to have solved the sign problem for a certain class of frustrated quantum spin systems through the use of a bipartite valence bond basis. We show in this Comment that the apparent positivity of the path integral is due to a misconception about the resolution of the identity operator in this basis and that consequently the sign problem remains a severe obstacle for the simulation of frustrated quantum magnets.

arXiv:0801.4498•Physical Review B•Apr 14, 2008•2 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

Sign-Problem-Free Monte Carlo Simulation of Certain Frustrated Quantum Magnets

Fabien Alet, Kedar Damle, Sumiran Pujari

We introduce a quantum Monte Carlo (QMC) method for efficient sign-problem-free simulations of a broad class of frustrated S=1/2 antiferromagnets using the basis of spin eigenstates of clusters to avoid the severe sign problem faced by other QMC methods. We demonstrate the utility of the method in several cases with competing exchange interactions and flag important limitations as well as possible extensions of the method.

arXiv:1511.01586•Physical Review Letters•Nov 4, 2016•55 citations

Reduction of the sign problem near T=0 in quantum Monte Carlo simulations

Jonathan D’Emidio, Stefan Weßel, Frédéric Mila

Building on a recent investigation of the Shastry-Sutherland model [S. Wessel et al., Phys. Rev. B 98, 174432 (2018)], we develop a general strategy to eliminate the Monte Carlo sign problem near the zero-temperature limit in frustrated quantum spin models. If the Hamiltonian of interest and the sign-problem-free Hamiltonian, obtained by making all off-diagonal elements negative in a given basis, have the same ground state and this state is a member of the computational basis, then the average sign returns to one as the temperature goes to zero. We illustrate this technique by studying the triangular and kagome lattice Heisenberg antiferrromagnet in a magnetic field above saturation, as well as the Heisenberg antiferromagnet on a modified Husimi cactus in the dimer basis. We also provide detailed Appendices on using linear programming techniques to automatically generate efficient directed loop updates in quantum Monte Carlo simulations.

arXiv:2003.13356•Physical review. B./Physical review. B•Aug 20, 2020•8 citations

Incorporation of density-matrix wave functions in Monte Carlo simulations: Application to the frustrated Heisenberg model

M. S. L. du Croo de Jongh, J. M. J. van Leeuwen, Wim van Saarloos

We combine the density-matrix renormalization-group (DMRG) technique with Green function Monte Carlo (GFMC) simulations using a special representation of the DMRG wave function. As a test case we apply the method to the two-dimensional frustrated Heisenberg antiferromagnet. By supplementing the branching in GFMC simulations with stochastic reconfiguration we get a stable simulation with a small variance also in the region where the fluctuations due to the minus sign problem are maximal. The sensitivity of the results to the choice of the guiding wave function is extensively investigated. In agreement with earlier calculations it is found from the DMRG wave function that for small ratios of the next-nearest-to-nearest neighbor coupling strength the system orders as a N\'eel-type antiferromagnet and for large ratios as a columnar antiferromagnet. The spin stiffness suggests an intermediate regime without magnetic long-range order. The energy curve indicates that the columnar phase is separated from the intermediate phase by a first-order transition. The combination of the DMRG and GFMC techniques allows us to substantiate this picture by calculating also the spin correlations in the system. We observe a pattern of spin correlations in the intermediate regime which is in between dimerlike and plaquette-type ordering, states that have recently been suggested. It is a state with strong dimerization in one direction and weaker dimerization in the perpendicular direction and thus it lacks the square symmetry of the plaquette state.

arXiv:cond-mat/0002116•Physical review. B, Condensed matter•Dec 1, 2000•45 citations
When one tries to simulate quantum spin systems by the Monte Carlo method, often the ``minus-sign problem'' is encountered. In such a case, an application of probabilistic methods is not possible. In this paper, the method has been proposed which allows to avoid the minus-sign problem for certain class of frustrated Heisenberg models. The systems where this method is applicable are, for instance, the pyrochlore lattice and the J1\ensuremath−J2{J}_{1}\ensuremath{-}{J}_{2}J1​\ensuremath−J2​ Heisenberg model. The method works in singlet sector. It relies on the expression of wave functions in dimer (pseudo)basis and writing down the Hamiltonian as a sum over plaquettes. In such a formulation, matrix elements of the exponent of Hamiltonian are positive.
A new method for the stabilization of the sign problem in the Green function Monte Carlo technique is proposed. The method is devised for real lattice Hamiltonians and is based on an iterative ``stochastic reconfiguration'' scheme which introduces some bias but allows a stable simulation with constant sign. The systematic reduction of this bias is possible in principle. The method is applied to the frustrated ${J}_{1}{\ensuremath{-}J}_{2}$ Heisenberg model, and tested against exact diagonalization data. Evidence of a finite spin gap for ${J}_{2}{/J}_{1}>\ensuremath{\sim}0.4$ is found in the thermodynamic limit.
Quantum Monte Carlo methods are sophisticated numerical techniques for simulating interacting quantum systems. In some cases, however, they suffer from the notorious ``sign problem'' and become too inefficient to be useful. A recent publication [J. Wojtkiewicz, Phys. Rev. B 75, 174421 (2007)] claims to have solved the sign problem for a certain class of frustrated quantum spin systems through the use of a bipartite valence bond basis. We show in this Comment that the apparent positivity of the path integral is due to a misconception about the resolution of the identity operator in this basis and that consequently the sign problem remains a severe obstacle for the simulation of frustrated quantum magnets.
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.
Building on a recent investigation of the Shastry-Sutherland model [S. Wessel et al., Phys. Rev. B 98, 174432 (2018)], we develop a general strategy to eliminate the Monte Carlo sign problem near the zero-temperature limit in frustrated quantum spin models. If the Hamiltonian of interest and the sign-problem-free Hamiltonian, obtained by making all off-diagonal elements negative in a given basis, have the same ground state and this state is a member of the computational basis, then the average sign returns to one as the temperature goes to zero. We illustrate this technique by studying the triangular and kagome lattice Heisenberg antiferrromagnet in a magnetic field above saturation, as well as the Heisenberg antiferromagnet on a modified Husimi cactus in the dimer basis. We also provide detailed Appendices on using linear programming techniques to automatically generate efficient directed loop updates in quantum Monte Carlo simulations.
We combine the density-matrix renormalization-group (DMRG) technique with Green function Monte Carlo (GFMC) simulations using a special representation of the DMRG wave function. As a test case we apply the method to the two-dimensional frustrated Heisenberg antiferromagnet. By supplementing the branching in GFMC simulations with stochastic reconfiguration we get a stable simulation with a small variance also in the region where the fluctuations due to the minus sign problem are maximal. The sensitivity of the results to the choice of the guiding wave function is extensively investigated. In agreement with earlier calculations it is found from the DMRG wave function that for small ratios of the next-nearest-to-nearest neighbor coupling strength the system orders as a N\'eel-type antiferromagnet and for large ratios as a columnar antiferromagnet. The spin stiffness suggests an intermediate regime without magnetic long-range order. The energy curve indicates that the columnar phase is separated from the intermediate phase by a first-order transition. The combination of the DMRG and GFMC techniques allows us to substantiate this picture by calculating also the spin correlations in the system. We observe a pattern of spin correlations in the intermediate regime which is in between dimerlike and plaquette-type ordering, states that have recently been suggested. It is a state with strong dimerization in one direction and weaker dimerization in the perpendicular direction and thus it lacks the square symmetry of the plaquette state.