Beyond the instanton gas approach: dominant thimbles approximation for the Hubbard model
read the original 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.
This paper has not been read by Pith yet.
Forward citations
Cited by 5 Pith papers
-
Enhancing the ergodicity of Worldvolume HMC via embedding generalized thimble HMC
Embedding generalized thimble HMC into worldvolume HMC improves ergodicity and phase-space exploration for sign-problem mitigation in 2D doped Hubbard model simulations, enabling larger lattices and controlled extrapolations.
-
Strain-Tuned Incommensurate Kekul\'e Spiral Order in Twisted Bilayer Graphene: a Quantum Many-Body Study
Combined QMC, ED, and HF calculations show a strain-driven transition from Kramers intervalley coherent to incommensurate Kekulé spiral order at ν = ±2 in twisted bilayer graphene.
-
Analyzing the two-dimensional doped Hubbard model with the Worldvolume HMC method
WV-HMC successfully simulates the doped 2D Hubbard model on 8x8 lattices at U/t=8 and T/t≈0.156 with controlled statistical errors.
-
Strain-Tuned Incommensurate Kekul\'e Spiral Order in Twisted Bilayer Graphene: a Quantum Many-Body Study
Quantum many-body calculations identify a strain-tuned transition from the Kramers intervalley coherent state to the incommensurate Kekulé spiral state at filling factors ±2 in twisted bilayer graphene.
-
Applying the Worldvolume Hybrid Monte Carlo method to the Hubbard model away from half filling
WV-HMC computes number and energy densities for the doped 2D Hubbard model on 6x6 and 8x8 lattices at U/t=8 and T/t≈0.156, showing effectiveness where standard DQMC fails.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.