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:2012.12283 · Physical review. B./Physical review. B

Quantum Monte Carlo simulation of generalized Kitaev models

Toshihiro Sato, Fakher F. Assaad

Published Aug 10, 2021• 25 citations
arXivHTMLPDF

Abstract

Frustrated spin systems generically suffer from the negative sign problem inherent in Monte Carlo methods. Since the severity of this problem is formulation dependent, optimization strategies can be put forward. We introduce a phase pinning approach in the realm of the auxiliary field quantum Monte Carlo algorithm. If we can find an antiunitary operator that commutes with the one-body Hamiltonian coupled to the auxiliary field, then the phase of the action is pinned to 0 and $\ensuremath{\pi}$. For generalized Kitaev models, we can successfully apply this strategy and observe a remarkable improvement of the average sign. We use this method to study the thermodynamical and dynamical properties of the Kitaev-Heisenberg model down to temperatures corresponding to half of the exchange coupling constant. Our dynamical data reveal finite temperature properties of ordered and spin-liquid phases inherent in this model.

DOI
10.1103/physrevb.104.l081106

Related papers

Negative sign free formulations of generalized Kitaev models with higher symmetries

Toshihiro Sato, Fakher F. Assaad

We provide a negative sign free formulation of the auxiliary field quantum Monte Carlo algorithm for generalized Kitaev models with higher symmetries. Our formulation is based on the Abrikosov fermion representation of the spin-1/2 degree of freedom and the phase pinning approach [T. Sato and F. F. Assaad, Phys. Rev. B 104, L081106 (2021)]. Enhancing the number of fermion flavors or orbitals from 1 to $N$ allows one to generalize the inherent ${Z}_{2}$ global symmetry to ${Z}_{2}\ifmmode\times\else\texttimes\fi{}\mathrm{S}\mathrm{U}{(N)}_{o}$. Using this general approach, we study the ${Z}_{2}\ifmmode\times\else\texttimes\fi{}\mathrm{S}\mathrm{U}{(2)}_{o}$ Kitaev-Heisenberg model reflecting the competition between the isotropic Heisenberg exchange and Kitaev-type bond-directional exchange interactions. We show that the symmetry enhancement provides a path to escape frustration and that the spin liquid phases in the original ${Z}_{2}$-symmetric model are not present in this model. Nevertheless, the ground-state phase diagram is extremely rich and has points with higher global and local continuous symmetries as well as deconfined quantum critical points.

arXiv:2204.12859•Physical review. B./Physical review. B•Oct 4, 2022•0 citations

Quantum Monte Carlo method on asymptotic Lefschetz thimbles for quantum spin systems: An application to the Kitaev model in a magnetic field

Petr A. Mishchenko, Yasuyuki Kato, Yukitoshi Motome

The quantum Monte Carlo method on asymptotic Lefschetz thimbles is a numerical algorithm devised specifically for alleviation of the sign problem appearing in the simulations of quantum many-body systems. In this method, the sign problem is alleviated by shifting the integration domain for the auxiliary fields, appearing, for example, in the conventional determinant quantum Monte Carlo method, from real space to an appropriate manifold in complex space. Here, we extend this method to quantum spin models with generic two-spin interactions, by using the Hubbard-Stratonovich transformation to decouple the exchange interactions and the Popov-Fedotov transformation to map the quantum spins to complex fermions. As a demonstration, we apply the method to the Kitaev model in a magnetic field whose ground state is predicted to deliver a topological quantum spin liquid with non-Abelian anyonic excitations. To illustrate how the sign problem is alleviated in this method, we visualize the asymptotic Lefschetz thimbles in complex space, together with the saddle points and the zeros of the fermion determinant. We benchmark our method in the low-temperature region in a magnetic field and show that the sign of the action is recovered considerably and unbiased numerical results are obtained with sufficient precision.

arXiv:2106.07937•Physical review. D/Physical review. D.•Oct 20, 2021•8 citations

Sign-problem-free Bosonic Auxiliary Field Quantum Monte Carlo

Yin-Kai Yu, Zhong-Chao Wei, Zi-Xiang Li

The sign problem poses a fundamental barrier to numerically exact quantum Monte Carlo simulations, severely hindering our understanding of quantum many-body physics across disciplines ranging from condensed matter to particle and nuclear physics. Geometrical frustration makes this barrier especially acute for interacting bosons. Here, we develop a bosonic auxiliary-field quantum Monte Carlo (BAFQMC) method and rigorously prove that time-reversal symmetry and reflection positivity each provide a general sufficient condition for sign-problem-free simulations. Crucially, under either criterion, BAFQMC evades the sign problem even on frustrated lattices where conventional world-line and stochastic-series-expansion methods fail. Benchmarking our approach on a two-flavor Bose-Hubbard model on the triangular lattice yields excellent quantitative agreement with finite-temperature exact diagonalization. By establishing unifying principles to solve the sign problem in BAFQMC, our work substantially expands the reach of unbiased, numerically exact simulations for frustrated bosonic matter.

arXiv:2610.31541•Oct 22, 2026•0 citations

Field-induced phase transitions in the Kitaev-Heisenberg model: A sign-problem-free quantum Monte Carlo study and possible application to α − RuCl 3

Xiangbin Zou, Shuo Liu, Wenan Guo +1 more

The frustrated magnet $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$ is one of the prime candidates for realizing a Kitaev quantum spin liquid (QSL). However, the existence of a field-induced intermediate QSL phase in this material remains under debate. Here, we employ sign-free numerically exact quantum Monte Carlo simulations to investigate the Kitaev-Heisenberg (KH) model on the honeycomb lattice with $K=\ensuremath{-}2J$ under an applied magnetic field along the $z$ direction. Our findings reveal that the system undergoes a direct quantum phase transition from a zigzag magnetically ordered phase to a spin-polarized phase at zero temperature, which belongs to the 3D XY universality class. At finite temperatures, a Berezinskii-Kosterlitz-Thouless transition line separates the spin-polarized phase from a quasi-long-range ordered state, eventually terminating at the quantum critical point. Our results convincingly show that there is no intermediate QSL phase in the KH model with a $z$-direction magnetic field, which we believe will shed important light on understanding experimental observations in $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$.

arXiv:2501.12437•Physical review. B./Physical review. B•Sep 8, 2025•2 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

Variational and auxiliary field Monte Carlo for the Hubbard and Hubbard-Holstein models: an accurate finite-size scaling and a "sign problem" solution

Seher Karakuzu

Understanding the properties of strongly correlated materials is a tough work since most of the times analytical tools cannot be applied due to the complicated nature of the problem. Here we apply quantum Monte Carlo techniques in order to tackle with one of the most important and unknown phenomena of many materials: the mechanism of superconductivity. We study the Hubbard-Holstein model which is one of the most simple theoretical models including strong correlation and electron-phonon coupling. By using twist-averaging boundary conditions we eliminate finite-size errors in the Monte Carlo simulations. With this useful tool we first solve the Hubbard-Holstein model by variational Monte Carlo by reporting the phase diagram of the model for different phonon frequencies at half filling. Then we investigate the phase separation away from half filling. Finally we try to attack the problem with an essentially unbiased method based on the auxiliary field quantum Monte Carlo technique and accelerated Langevin dynamics. By curing the sign problem via Cauchy integration in the complex plane where the auxiliary fields have been defined, we report the effect of the electron-phonon coupling on some observables such as magnetization, spin and charge structure factors. We show preliminary results that are already meaningful to understand the nature of the transition between magnetism and charge order in the model at half filling.

•Oct 5, 2018•0 citations
Frustrated spin systems generically suffer from the negative sign problem inherent in Monte Carlo methods. Since the severity of this problem is formulation dependent, optimization strategies can be put forward. We introduce a phase pinning approach in the realm of the auxiliary field quantum Monte Carlo algorithm. If we can find an antiunitary operator that commutes with the one-body Hamiltonian coupled to the auxiliary field, then the phase of the action is pinned to 0 and \ensuremathπ\ensuremath{\pi}\ensuremathπ. For generalized Kitaev models, we can successfully apply this strategy and observe a remarkable improvement of the average sign. We use this method to study the thermodynamical and dynamical properties of the Kitaev-Heisenberg model down to temperatures corresponding to half of the exchange coupling constant. Our dynamical data reveal finite temperature properties of ordered and spin-liquid phases inherent in this model.
We provide a negative sign free formulation of the auxiliary field quantum Monte Carlo algorithm for generalized Kitaev models with higher symmetries. Our formulation is based on the Abrikosov fermion representation of the spin-1/2 degree of freedom and the phase pinning approach [T. Sato and F. F. Assaad, Phys. Rev. B 104, L081106 (2021)]. Enhancing the number of fermion flavors or orbitals from 1 to NNN allows one to generalize the inherent Z2{Z}_{2}Z2​ global symmetry to ${Z}_{2}\ifmmode\times\else\texttimes\fi{}\mathrm{S}\mathrm{U}{(N)}_{o}$. Using this general approach, we study the ${Z}_{2}\ifmmode\times\else\texttimes\fi{}\mathrm{S}\mathrm{U}{(2)}_{o}$ Kitaev-Heisenberg model reflecting the competition between the isotropic Heisenberg exchange and Kitaev-type bond-directional exchange interactions. We show that the symmetry enhancement provides a path to escape frustration and that the spin liquid phases in the original Z2{Z}_{2}Z2​-symmetric model are not present in this model. Nevertheless, the ground-state phase diagram is extremely rich and has points with higher global and local continuous symmetries as well as deconfined quantum critical points.
The quantum Monte Carlo method on asymptotic Lefschetz thimbles is a numerical algorithm devised specifically for alleviation of the sign problem appearing in the simulations of quantum many-body systems. In this method, the sign problem is alleviated by shifting the integration domain for the auxiliary fields, appearing, for example, in the conventional determinant quantum Monte Carlo method, from real space to an appropriate manifold in complex space. Here, we extend this method to quantum spin models with generic two-spin interactions, by using the Hubbard-Stratonovich transformation to decouple the exchange interactions and the Popov-Fedotov transformation to map the quantum spins to complex fermions. As a demonstration, we apply the method to the Kitaev model in a magnetic field whose ground state is predicted to deliver a topological quantum spin liquid with non-Abelian anyonic excitations. To illustrate how the sign problem is alleviated in this method, we visualize the asymptotic Lefschetz thimbles in complex space, together with the saddle points and the zeros of the fermion determinant. We benchmark our method in the low-temperature region in a magnetic field and show that the sign of the action is recovered considerably and unbiased numerical results are obtained with sufficient precision.
The sign problem poses a fundamental barrier to numerically exact quantum Monte Carlo simulations, severely hindering our understanding of quantum many-body physics across disciplines ranging from condensed matter to particle and nuclear physics. Geometrical frustration makes this barrier especially acute for interacting bosons. Here, we develop a bosonic auxiliary-field quantum Monte Carlo (BAFQMC) method and rigorously prove that time-reversal symmetry and reflection positivity each provide a general sufficient condition for sign-problem-free simulations. Crucially, under either criterion, BAFQMC evades the sign problem even on frustrated lattices where conventional world-line and stochastic-series-expansion methods fail. Benchmarking our approach on a two-flavor Bose-Hubbard model on the triangular lattice yields excellent quantitative agreement with finite-temperature exact diagonalization. By establishing unifying principles to solve the sign problem in BAFQMC, our work substantially expands the reach of unbiased, numerically exact simulations for frustrated bosonic matter.
The frustrated magnet $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$ is one of the prime candidates for realizing a Kitaev quantum spin liquid (QSL). However, the existence of a field-induced intermediate QSL phase in this material remains under debate. Here, we employ sign-free numerically exact quantum Monte Carlo simulations to investigate the Kitaev-Heisenberg (KH) model on the honeycomb lattice with $K=\ensuremath{-}2J$ under an applied magnetic field along the zzz direction. Our findings reveal that the system undergoes a direct quantum phase transition from a zigzag magnetically ordered phase to a spin-polarized phase at zero temperature, which belongs to the 3D XY universality class. At finite temperatures, a Berezinskii-Kosterlitz-Thouless transition line separates the spin-polarized phase from a quasi-long-range ordered state, eventually terminating at the quantum critical point. Our results convincingly show that there is no intermediate QSL phase in the KH model with a zzz-direction magnetic field, which we believe will shed important light on understanding experimental observations in $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{RuCl}}_{3}$.
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.
Understanding the properties of strongly correlated materials is a tough work since most of the times analytical tools cannot be applied due to the complicated nature of the problem. Here we apply quantum Monte Carlo techniques in order to tackle with one of the most important and unknown phenomena of many materials: the mechanism of superconductivity. We study the Hubbard-Holstein model which is one of the most simple theoretical models including strong correlation and electron-phonon coupling. By using twist-averaging boundary conditions we eliminate finite-size errors in the Monte Carlo simulations. With this useful tool we first solve the Hubbard-Holstein model by variational Monte Carlo by reporting the phase diagram of the model for different phonon frequencies at half filling. Then we investigate the phase separation away from half filling. Finally we try to attack the problem with an essentially unbiased method based on the auxiliary field quantum Monte Carlo technique and accelerated Langevin dynamics. By curing the sign problem via Cauchy integration in the complex plane where the auxiliary fields have been defined, we report the effect of the electron-phonon coupling on some observables such as magnetization, spin and charge structure factors. We show preliminary results that are already meaningful to understand the nature of the transition between magnetism and charge order in the model at half filling.