© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2106.07937 · Physical review. D/Physical review. D.

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

Published Oct 20, 2021• 8 citations

Abstract

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.

DOI: 10.1103/physrevd.104.074517

Related papers

Path integral representation for the Hubbard model with reduced number of Lefschetz thimbles

Maksim Ulybyshev, S. N. Valgushev

The concept of Lefschetz thimble decomposition is one of the most promising possible modifications of Quantum Monte Carlo (QMC) algorithms aimed at alleviating the sign problem which appears in many interesting physical situations, e.g. in the Hubbard model away from half filling. In this approach one utilizes the fact that the integral over real variables with an integrand containing a complex fluctuating phase is equivalent to the sum of integrals over special manifolds in complex space ("Lefschetz thimbles"), each of them having a fixed complex phase factor. Thus, the sign problem can be reduced if the resulting sum contains terms with only a few different phases. We explore the complexity of the sign problem for the few-site Hubbard model on square lattice combining a semi-analytical study of saddle points and thimbles in small lattices with several steps in Euclidean time with results of test QMC calculations. We check different variants of conventional Hubbard-Stratonovich transformation based on Gaussian integrals. On the basis of our analysis we reveal a regime with minimal number of relevant thimbles in the vicinity of half filling. In this regime we found only two relevant thimbles for the few-site lattices studied in the paper. There is also indirect evidence of the existence of this regime in more realistic systems with large number of Euclidean time slices. In addition, we derive a new non-Gaussian representation of the interaction term, where the number of relevant Lefschetz is also reduced in comparison with conventional Gaussian Hubbard-Stratonovich transformation.

arXiv:1712.02188•Dec 6, 2017•11 citations

Lefschetz thimble quantum Monte Carlo for spin systems

T. C. Mooney, Jacob Bringewatt, Neill C. Warrington +1 more

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Therefore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

arXiv:2110.10699•Physical review. B./Physical review. B•Dec 15, 2022•5 citations

Lefschetz Thimbles and Wigner Functions

Vladimir S. Filinov, Alexander S. Larkin

The major difficulty for the Feynman Path Integral Monte Carlo (PIMC) simulations of the quantum systems of particles is the so called “sign problem”, arising due to the fast oscillations of the path integral integrand depending on the complex-valued action. Our aim is to find universal techniques being able to solve this problem. The new method combines the basic ideas of the Metropolis and Hasting algorithms and is based on the Picard-Lefschetz theory and complex-valued version of Morse theory. The basic idea is to choose the Lefschetz thimbles as manifolds approaching the saddle point of the integrand. On this thimble the imaginary part of the complex-valued action remains constant. As a result the integrand on each thimble does not oscillate, so the “sign problem” disappears and the integral can be calculated much more effectively. The developed approach allows also finding saddle points in the complexified space of path integral integration. Some simple test calculations and comparisons with available analytical results have been carried out.

Journal of Applied Mathematics and Physics•Jan 1, 2020•1 citation

Tempered Lefschetz thimble method and its application to the Hubbard model away from half filling

Masafumi Fukuma, Nobuyuki Matsumoto, Naoya Umeda

The tempered Lefschetz thimble method (TLTM) is a parallel-tempering algorithm towards solving the numerical sign problem. It tames both the sign and ergodicity problems simultaneously by tempering the system with the flow time of continuous deformations of the integration region. In this article, after reviewing the basics of the TLTM, we explain a new algorithm within the TLTM that enables us to estimate the expectation values precisely with a criterion ensuring global equilibrium and the sufficiency of the sample size. To demonstrate the effectiveness of the algorithm, we apply the TLTM to the quantum Monte Carlo simulation of the Hubbard model away from half filling on a two-dimensional lattice of small size, and show that the obtained numerical results agree nicely with exact values.

arXiv:2001.01665•International Symposium on Lattice Field Theory (LATTICE)•Jan 4, 2020•4 citations

Lefschetz thimble Monte Carlo for many-body theories: A Hubbard model study

Abhishek Mukherjee, M. Cristoforetti

Recently, a method based on stochastic integration on the surfaces of steepest descent of the action was introduced to tackle the sign problem in quantum-field theories. We show how this method can be used in many-body theories to perform fully nonperturbative calculations of quantum corrections about mean-field solutions. We discuss an explicit algorithm for implementing our method, and present results for the repulsive Hubbard model away from half-filling at intermediate temperatures. Our results are consistent with those from other state-of-the-art calculations.

Physical Review B•Jul 25, 2014•49 citations

Lefschetz thimble Monte Carlo for many body theories: application to the repulsive Hubbard model away from half filling

Mukherjee, Abhishek, M. Cristoforetti

Recently, a new method, based on stochastic integration on the surfaces of steepest descent of the action, was introduced to tackle the sign problem in quantum field theories. We show how this method can be used in many body theories to perform fully non-perturbative calculations of quantum corrections about mean field solutions. We discuss an explicit algorithm for implementing our method, and present results for the repulsive Hubbard model away from half-filling at intermediate temperatures. Our results are consistent with those from the recent state of the art cluster dynamical mean field theory calculations.

arXiv:1403.5680•Mar 22, 2014•17 citations

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 concept of Lefschetz thimble decomposition is one of the most promising possible modifications of Quantum Monte Carlo (QMC) algorithms aimed at alleviating the sign problem which appears in many interesting physical situations, e.g. in the Hubbard model away from half filling. In this approach one utilizes the fact that the integral over real variables with an integrand containing a complex fluctuating phase is equivalent to the sum of integrals over special manifolds in complex space ("Lefschetz thimbles"), each of them having a fixed complex phase factor. Thus, the sign problem can be reduced if the resulting sum contains terms with only a few different phases. We explore the complexity of the sign problem for the few-site Hubbard model on square lattice combining a semi-analytical study of saddle points and thimbles in small lattices with several steps in Euclidean time with results of test QMC calculations. We check different variants of conventional Hubbard-Stratonovich transformation based on Gaussian integrals. On the basis of our analysis we reveal a regime with minimal number of relevant thimbles in the vicinity of half filling. In this regime we found only two relevant thimbles for the few-site lattices studied in the paper. There is also indirect evidence of the existence of this regime in more realistic systems with large number of Euclidean time slices. In addition, we derive a new non-Gaussian representation of the interaction term, where the number of relevant Lefschetz is also reduced in comparison with conventional Gaussian Hubbard-Stratonovich transformation.

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Therefore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

The major difficulty for the Feynman Path Integral Monte Carlo (PIMC) simulations of the quantum systems of particles is the so called “sign problem”, arising due to the fast oscillations of the path integral integrand depending on the complex-valued action. Our aim is to find universal techniques being able to solve this problem. The new method combines the basic ideas of the Metropolis and Hasting algorithms and is based on the Picard-Lefschetz theory and complex-valued version of Morse theory. The basic idea is to choose the Lefschetz thimbles as manifolds approaching the saddle point of the integrand. On this thimble the imaginary part of the complex-valued action remains constant. As a result the integrand on each thimble does not oscillate, so the “sign problem” disappears and the integral can be calculated much more effectively. The developed approach allows also finding saddle points in the complexified space of path integral integration. Some simple test calculations and comparisons with available analytical results have been carried out.

The tempered Lefschetz thimble method (TLTM) is a parallel-tempering algorithm towards solving the numerical sign problem. It tames both the sign and ergodicity problems simultaneously by tempering the system with the flow time of continuous deformations of the integration region. In this article, after reviewing the basics of the TLTM, we explain a new algorithm within the TLTM that enables us to estimate the expectation values precisely with a criterion ensuring global equilibrium and the sufficiency of the sample size. To demonstrate the effectiveness of the algorithm, we apply the TLTM to the quantum Monte Carlo simulation of the Hubbard model away from half filling on a two-dimensional lattice of small size, and show that the obtained numerical results agree nicely with exact values.

Recently, a method based on stochastic integration on the surfaces of steepest descent of the action was introduced to tackle the sign problem in quantum-field theories. We show how this method can be used in many-body theories to perform fully nonperturbative calculations of quantum corrections about mean-field solutions. We discuss an explicit algorithm for implementing our method, and present results for the repulsive Hubbard model away from half-filling at intermediate temperatures. Our results are consistent with those from other state-of-the-art calculations.

Recently, a new method, based on stochastic integration on the surfaces of steepest descent of the action, was introduced to tackle the sign problem in quantum field theories. We show how this method can be used in many body theories to perform fully non-perturbative calculations of quantum corrections about mean field solutions. We discuss an explicit algorithm for implementing our method, and present results for the repulsive Hubbard model away from half-filling at intermediate temperatures. Our results are consistent with those from the recent state of the art cluster dynamical mean field theory calculations.