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:1704.01973 · SciPost Physics

Efficient Quantum Monte Carlo simulations of highly frustrated magnets: the frustrated spin-1/2 ladder

Stefan Weßel, Bruce Normand, Frédéric Mila, A. Honecker

Published Jul 18, 2017• 43 citations
arXivHTMLPDF

Abstract

Quantum Monte Carlo simulations provide one of the more powerful and versatile numerical approaches to condensed matter systems. However, their application to frustrated quantum spin models, in all relevant temperature regimes, is hamstrung by the infamous “sign problem.” Here we exploit the fact that the sign problem is basis-dependent. Recent studies have shown that passing to a dimer (two-site) basis eliminates the sign problem completely for a fully frustrated spin model on the two-leg ladder. We generalize this result to all partially frustrated two-leg spin-1/2 ladders, meaning those where the diagonal and leg couplings take any antiferromagnetic values. We find that, although the sign problem does reappear, it remains remarkably mild throughout the entire phase diagram. We explain this result and apply it to perform efficient quantum Monte Carlo simulations of frustrated ladders, obtaining accurate results for thermodynamic quantities such as the magnetic specific heat and susceptibility of ladders up to L = 200 L = 200 rungs (400 spins 1/2) and down to very low temperatures.

DOI
10.21468/scipostphys.3.1.005

Related papers

Thermodynamic properties of highly frustrated quantum spin ladders: Influence of many-particle bound states

A. Honecker, Stefan Weßel, René John Kerkdyk +3 more

Quantum antiferromagnets have proven to be some of the cleanest realizations available for theoretical, numerical, and experimental studies of quantum fluctuation effects. At finite temperatures, however, the additional effects of thermal fluctuations in the restricted phase space of a low-dimensional system have received much less attention, particularly the situation in frustrated quantum magnets, where the excitations may be complex collective (bound or even fractionalized) modes. We investigate this problem by studying the thermodynamic properties of the frustrated two-leg $S=\frac{1}{2}$ spin ladder, with particular emphasis on the fully frustrated case. We present numerical results for the magnetic specific heat and susceptibility, obtained from exact diagonalization and quantum Monte Carlo studies, which we show can be rendered free of the sign problem even in a strongly frustrated system and which allow us to reach unprecedented sizes of $L=200$ ladder rungs. We find that frustration effects cause an unconventional evolution of the thermodynamic response across the full parameter regime of the model. However, close to the first-order transition they cause a highly anomalous reduction in temperature scales with no concomitant changes in the gap; the specific heat shows a very narrow peak at very low energies and the susceptibility rises abruptly at extremely low temperatures. Unusually, the two quantities have different gaps over an extended region of the parameter space. We demonstrate that these results reflect the presence of large numbers of multiparticle bound-state excitations, whose energies fall below the one-triplon gap in the transition region.

arXiv:1511.01501•Physical review. B./Physical review. B•Feb 8, 2016•48 citations

Quantum Monte Carlo simulations in the trimer basis: first-order transitions and thermal critical points in frustrated trilayer magnets

Lukas Weber, A. Honecker, Bruce Normand +3 more

The phase diagrams of highly frustrated quantum spin systems can exhibit first-order quantum phase transitions and thermal critical points even in the absence of any long-ranged magnetic order. However, all unbiased numerical techniques for investigating frustrated quantum magnets face significant challenges, and for generic quantum Monte Carlo methods the challenge is the sign problem. Here we report on a general quantum Monte Carlo approach with a loop-update scheme that operates in any basis, and we show that, with an appropriate choice of basis, it allows us to study a frustrated model of coupled spin-1/2 trimers: simulations of the trilayer Heisenberg antiferromagnet in the spin-trimer basis are sign-problem-free when the intertrimer couplings are fully frustrated. This model features a first-order quantum phase transition, from which a line of first-order transitions emerges at finite temperatures and terminates in a thermal critical point. The trimer unit cell hosts an internal degree of freedom that can be controlled to induce an extensive entropy jump at the quantum transition, which alters the shape of the first-order line. We explore the consequences for the thermal properties in the vicinity of the critical point, which include profound changes in the lines of maxima defined by the specific heat. Our findings reveal trimer quantum magnets as fundamental systems capturing in full the complex thermal physics of the strongly frustrated regime.

arXiv:2105.05271•SciPost Physics•Feb 8, 2022•25 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

Quantum Monte Carlo Calculation for Spin Systems in Two-Leg Ladders with Alternating Heisenberg-Bond Chains

Takafumi Momose, Toshihiko Suzuki, Yuhei Natsume

Thermodynamic properties of spin=1/2 two-leg ladders with alternating chains including ferromagnetic and antiferromagnetic interactions are discussed on the basis of the numerical results obtained by quantum Monte Carlo simulations. As well as the role of exchange interaction of rungs, the essential role of the bond alternation in the chain on ladder systems is investigated in connection with the appearance of characteristic properties in magnetic susceptibilities and heat capacities. Furthermore, quantum effects on the ground-state energy are obtained. On the basis of calculated results, the relationship between singlet dimer phases and Haldane states is discussed.

Journal of the Physical Society of Japan•Jun 15, 1997•0 citations

Thermodynamic properties of the Shastry-Sutherland model from quantum Monte Carlo simulations

Stefan Weßel, Ido Niesen, Jonas Stapmanns +4 more

We investigate the minus-sign problem that afflicts quantum Monte Carlo (QMC) simulations of frustrated quantum spin systems, focusing on spin $S=1/2$, two spatial dimensions, and the extended Shastry-Sutherland model. We show that formulating the Hamiltonian in the diagonal dimer basis leads to a sign problem that becomes negligible at low temperatures for small and intermediate values of the ratio of the inter- and intradimer couplings. This is a consequence of the fact that the product state of dimer singlets is the exact ground state both of the extended Shastry-Sutherland model and of a corresponding ``sign-problem-free'' model, obtained by changing the signs of all positive off-diagonal matrix elements in the dimer basis. By exploiting this insight, we map the sign problem throughout the extended parameter space from the Shastry-Sutherland to the fully frustrated bilayer model and compare it with the phase diagram computed by tensor-network methods. We use QMC to compute with high accuracy the temperature dependence of the magnetic specific heat and susceptibility of the Shastry-Sutherland model for large systems up to a coupling ratio of 0.526(1) and down to zero temperature. For larger coupling ratios, our QMC results assist us in benchmarking the evolution of the thermodynamic properties by systematic comparison with exact diagonalization calculations and interpolated high-temperature series expansions.

arXiv:1808.02043•Physical review. B./Physical review. B•Nov 27, 2018•47 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
Quantum Monte Carlo simulations provide one of the more powerful and versatile numerical approaches to condensed matter systems. However, their application to frustrated quantum spin models, in all relevant temperature regimes, is hamstrung by the infamous “sign problem.” Here we exploit the fact that the sign problem is basis-dependent. Recent studies have shown that passing to a dimer (two-site) basis eliminates the sign problem completely for a fully frustrated spin model on the two-leg ladder. We generalize this result to all partially frustrated two-leg spin-1/2 ladders, meaning those where the diagonal and leg couplings take any antiferromagnetic values. We find that, although the sign problem does reappear, it remains remarkably mild throughout the entire phase diagram. We explain this result and apply it to perform efficient quantum Monte Carlo simulations of frustrated ladders, obtaining accurate results for thermodynamic quantities such as the magnetic specific heat and susceptibility of ladders up to L = 200 L = 200 rungs (400 spins 1/2) and down to very low temperatures.
Quantum antiferromagnets have proven to be some of the cleanest realizations available for theoretical, numerical, and experimental studies of quantum fluctuation effects. At finite temperatures, however, the additional effects of thermal fluctuations in the restricted phase space of a low-dimensional system have received much less attention, particularly the situation in frustrated quantum magnets, where the excitations may be complex collective (bound or even fractionalized) modes. We investigate this problem by studying the thermodynamic properties of the frustrated two-leg S=12S=\frac{1}{2}S=21​ spin ladder, with particular emphasis on the fully frustrated case. We present numerical results for the magnetic specific heat and susceptibility, obtained from exact diagonalization and quantum Monte Carlo studies, which we show can be rendered free of the sign problem even in a strongly frustrated system and which allow us to reach unprecedented sizes of L=200L=200L=200 ladder rungs. We find that frustration effects cause an unconventional evolution of the thermodynamic response across the full parameter regime of the model. However, close to the first-order transition they cause a highly anomalous reduction in temperature scales with no concomitant changes in the gap; the specific heat shows a very narrow peak at very low energies and the susceptibility rises abruptly at extremely low temperatures. Unusually, the two quantities have different gaps over an extended region of the parameter space. We demonstrate that these results reflect the presence of large numbers of multiparticle bound-state excitations, whose energies fall below the one-triplon gap in the transition region.
The phase diagrams of highly frustrated quantum spin systems can exhibit first-order quantum phase transitions and thermal critical points even in the absence of any long-ranged magnetic order. However, all unbiased numerical techniques for investigating frustrated quantum magnets face significant challenges, and for generic quantum Monte Carlo methods the challenge is the sign problem. Here we report on a general quantum Monte Carlo approach with a loop-update scheme that operates in any basis, and we show that, with an appropriate choice of basis, it allows us to study a frustrated model of coupled spin-1/2 trimers: simulations of the trilayer Heisenberg antiferromagnet in the spin-trimer basis are sign-problem-free when the intertrimer couplings are fully frustrated. This model features a first-order quantum phase transition, from which a line of first-order transitions emerges at finite temperatures and terminates in a thermal critical point. The trimer unit cell hosts an internal degree of freedom that can be controlled to induce an extensive entropy jump at the quantum transition, which alters the shape of the first-order line. We explore the consequences for the thermal properties in the vicinity of the critical point, which include profound changes in the lines of maxima defined by the specific heat. Our findings reveal trimer quantum magnets as fundamental systems capturing in full the complex thermal physics of the strongly frustrated regime.
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.
Thermodynamic properties of spin=1/2 two-leg ladders with alternating chains including ferromagnetic and antiferromagnetic interactions are discussed on the basis of the numerical results obtained by quantum Monte Carlo simulations. As well as the role of exchange interaction of rungs, the essential role of the bond alternation in the chain on ladder systems is investigated in connection with the appearance of characteristic properties in magnetic susceptibilities and heat capacities. Furthermore, quantum effects on the ground-state energy are obtained. On the basis of calculated results, the relationship between singlet dimer phases and Haldane states is discussed.
We investigate the minus-sign problem that afflicts quantum Monte Carlo (QMC) simulations of frustrated quantum spin systems, focusing on spin S=1/2S=1/2S=1/2, two spatial dimensions, and the extended Shastry-Sutherland model. We show that formulating the Hamiltonian in the diagonal dimer basis leads to a sign problem that becomes negligible at low temperatures for small and intermediate values of the ratio of the inter- and intradimer couplings. This is a consequence of the fact that the product state of dimer singlets is the exact ground state both of the extended Shastry-Sutherland model and of a corresponding ``sign-problem-free'' model, obtained by changing the signs of all positive off-diagonal matrix elements in the dimer basis. By exploiting this insight, we map the sign problem throughout the extended parameter space from the Shastry-Sutherland to the fully frustrated bilayer model and compare it with the phase diagram computed by tensor-network methods. We use QMC to compute with high accuracy the temperature dependence of the magnetic specific heat and susceptibility of the Shastry-Sutherland model for large systems up to a coupling ratio of 0.526(1) and down to zero temperature. For larger coupling ratios, our QMC results assist us in benchmarking the evolution of the thermodynamic properties by systematic comparison with exact diagonalization calculations and interpolated high-temperature series expansions.
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.