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:2407.09452

Beyond the instanton gas approach: dominant thimbles approximation for the Hubbard model

Maksim Ulybyshev, Fakher F. Assaad

Published Jul 12, 2024• 0 citations
arXivHTMLPDF

Abstract

To each complex saddle point of an action, one can attach a Lefschetz thimble on which the imaginary part of the action is constant. Cauchy theorem states that summation over a set of thimbles produces the exact result. This reorganization of the path integral, is an appealing starting point for various approximations: In the realm of auxiliary quantum Monte Carlo methods it provides a framework to alleviate the negative sign problem. Here, we suggest to constrain the integration to the \textit{dominant} thimbles: the thimbles attached to the saddle points with the largest statistical weight. For the Hubbard model, in a formulation where the the Hubbard Stratonovitch field couples to the charge, this provides a \textit{symmetry} consistent approximation to the physics of the Hubbard model: constraining the integration domain does not explicitly break a symmetry. We can test this approach for the Hubbard model at half-filling on a bipartite lattice. The paper builds on the previously developed instanton gas approach, where an exhaustive saddle point approximation was constructed. We present results, showing that the dominant thimbles approximation provides results that are in remarkable agreement with the exact results for various fermionic observables including spin and charge order parameters and single electron spectral functions. We discuss implications of our results for simulations away of half filling.

DOI
10.48550/arxiv.2407.09452

Related papers

Instanton gas approach to the Hubbard model

Maksim Ulybyshev, Christopher Winterowd, Fakher F. Assaad +1 more

In this paper, we consider a path integral formulation of the Hubbard model based on a Hubbard-Stratonovich transformation that couples the auxiliary field to the local electronic density. This decoupling is known to have a saddle-point structure that shows a remarkable regularity: The field configuration at each saddle point can be understood in terms of a set of elementary field configurations localized in space and imaginary time which we coin instantons. The interaction between instantons is short ranged. Here, we formulate a classical partition function for the instanton gas that has predictive power. For a given set of physical parameters, we can predict the distribution of instantons and show that the instanton number is sharply defined in the thermodynamic limit, thereby defining a unique dominant saddle point. Decoupling in the charge channel conserves SU(2) spin symmetry for each field configurations. Hence, the instanton approach provides an SU(2) spin-symmetric approximation to the Hubbard model. It fails, however, to capture the magnetic transition inherent to the Hubbard model on the honeycomb lattice despite being able to describe local moment formation. In fact, the instanton itself corresponds to local moment formation and concomitant short-ranged antiferromagnetic correlations. This aspect is also seen in the single particle spectral function that shows clear signs of the upper and lower Hubbard bands. Our instanton approach bears remarkable similarities to local dynamical approaches, such as dynamical mean-field theory, in the sense that it has the unique property of allowing for local moment formation without breaking the SU(2) spin symmetry. In contrast to local approaches, it captures short-ranged magnetic fluctuations. Furthermore, it also offers possibilities for systematic improvements by taking into account fluctuations around the dominant saddle point. Finally, we show that the saddle point structure depends upon the choice of lattice geometry. For the square lattice at half filling, the saddle-point structure reflects the itinerant to localized nature of the magnetism as a function of the coupling strength. The implications of our results for Lefschetz thimble approaches to alleviate the sign problem are also discussed.

arXiv:2207.06297•Physical review. B./Physical review. B•Jan 30, 2023•3 citations

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

Taming the sign problem of the finite density Hubbard model via Lefschetz thimbles

Maksim Ulybyshev, Christopher Winterowd, Savvas Zafeiropoulos

We study the sign problem in the Hubbard model on the hexagonal lattice away from half-filling using the Lefschetz thimbles method. We identify the saddle points, reduce their amount, and perform quantum Monte Carlo (QMC) simulations using the holomorphic gradient flow to deform the integration contour in complex space. Finally, the results are compared against exact diagonalization (ED). We show that the sign problem can be substantially weakened, even in the regime with low temperature and large chemical potential, where standard QMC techniques exhibit an exponential decay of the average sign.

arXiv:1906.02726•Jun 6, 2019•12 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 analysis of the sign problem in one-site fermion model

Yuya Tanizaki, Yoshimasa Hidaka, Tomoya Hayata

The Lefschetz-thimble approach to path integrals is applied to a one-site model of electrons, i.e., the one-site Hubbard model. Since the one-site Hubbard model shows a non-analytic behavior at the zero temperature and its path integral expression has the sign problem, this toy model is a good testing ground for an idea or a technique to attack the sign problem. Semiclassical analysis using complex saddle points unveils the significance of interference among multiple Lefschetz thimbles to reproduce the non-analytic behavior by using the path integral. If the number of Lefschetz thimbles is insufficient, we found not only large discrepancies from the exact result, but also thermodynamic instabilities. Analyzing such singular behaviors semiclassically, we propose a criterion to identify the necessary number of Lefschetz thimbles. We argue that this interference of multiple saddle points is a key issue to understand the sign problem of the finite-density quantum chromodynamics.

arXiv:1509.07146•New Journal of Physics•Mar 1, 2016•91 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
To each complex saddle point of an action, one can attach a Lefschetz thimble on which the imaginary part of the action is constant. Cauchy theorem states that summation over a set of thimbles produces the exact result. This reorganization of the path integral, is an appealing starting point for various approximations: In the realm of auxiliary quantum Monte Carlo methods it provides a framework to alleviate the negative sign problem. Here, we suggest to constrain the integration to the \textit{dominant} thimbles: the thimbles attached to the saddle points with the largest statistical weight. For the Hubbard model, in a formulation where the the Hubbard Stratonovitch field couples to the charge, this provides a \textit{symmetry} consistent approximation to the physics of the Hubbard model: constraining the integration domain does not explicitly break a symmetry. We can test this approach for the Hubbard model at half-filling on a bipartite lattice. The paper builds on the previously developed instanton gas approach, where an exhaustive saddle point approximation was constructed. We present results, showing that the dominant thimbles approximation provides results that are in remarkable agreement with the exact results for various fermionic observables including spin and charge order parameters and single electron spectral functions. We discuss implications of our results for simulations away of half filling.
In this paper, we consider a path integral formulation of the Hubbard model based on a Hubbard-Stratonovich transformation that couples the auxiliary field to the local electronic density. This decoupling is known to have a saddle-point structure that shows a remarkable regularity: The field configuration at each saddle point can be understood in terms of a set of elementary field configurations localized in space and imaginary time which we coin instantons. The interaction between instantons is short ranged. Here, we formulate a classical partition function for the instanton gas that has predictive power. For a given set of physical parameters, we can predict the distribution of instantons and show that the instanton number is sharply defined in the thermodynamic limit, thereby defining a unique dominant saddle point. Decoupling in the charge channel conserves SU(2) spin symmetry for each field configurations. Hence, the instanton approach provides an SU(2) spin-symmetric approximation to the Hubbard model. It fails, however, to capture the magnetic transition inherent to the Hubbard model on the honeycomb lattice despite being able to describe local moment formation. In fact, the instanton itself corresponds to local moment formation and concomitant short-ranged antiferromagnetic correlations. This aspect is also seen in the single particle spectral function that shows clear signs of the upper and lower Hubbard bands. Our instanton approach bears remarkable similarities to local dynamical approaches, such as dynamical mean-field theory, in the sense that it has the unique property of allowing for local moment formation without breaking the SU(2) spin symmetry. In contrast to local approaches, it captures short-ranged magnetic fluctuations. Furthermore, it also offers possibilities for systematic improvements by taking into account fluctuations around the dominant saddle point. Finally, we show that the saddle point structure depends upon the choice of lattice geometry. For the square lattice at half filling, the saddle-point structure reflects the itinerant to localized nature of the magnetism as a function of the coupling strength. The implications of our results for Lefschetz thimble approaches to alleviate the sign problem are also discussed.
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.
We study the sign problem in the Hubbard model on the hexagonal lattice away from half-filling using the Lefschetz thimbles method. We identify the saddle points, reduce their amount, and perform quantum Monte Carlo (QMC) simulations using the holomorphic gradient flow to deform the integration contour in complex space. Finally, the results are compared against exact diagonalization (ED). We show that the sign problem can be substantially weakened, even in the regime with low temperature and large chemical potential, where standard QMC techniques exhibit an exponential decay of the average sign.
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.
The Lefschetz-thimble approach to path integrals is applied to a one-site model of electrons, i.e., the one-site Hubbard model. Since the one-site Hubbard model shows a non-analytic behavior at the zero temperature and its path integral expression has the sign problem, this toy model is a good testing ground for an idea or a technique to attack the sign problem. Semiclassical analysis using complex saddle points unveils the significance of interference among multiple Lefschetz thimbles to reproduce the non-analytic behavior by using the path integral. If the number of Lefschetz thimbles is insufficient, we found not only large discrepancies from the exact result, but also thermodynamic instabilities. Analyzing such singular behaviors semiclassically, we propose a criterion to identify the necessary number of Lefschetz thimbles. We argue that this interference of multiple saddle points is a key issue to understand the sign problem of the finite-density quantum chromodynamics.
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.