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

Abstract

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.

Mitigating spikes in fermion Monte Carlo methods by reshuffling measurements

“…where (18) and the pair (α, β) parametrizes the surface of constant observable O in 3D space (x 1 ; x 2 ; x2 ).…”

“…Here, the zeros of the determinant form "domain walls" in configuration space such that the dimensionality the zero-manifolds is given V L τ − 1, where V L τ corresponds to the number of auxiliary fields. The connection to the power law in heavytailed distributions was considered in details in [18] and for the particular case of zeros of the determinant forming domain walls, ξ = 5/2. A detailed proof is given in Appendix B.…”

Schur complement solver for Quantum Monte-Carlo simulations of strongly interacting fermions

“…This form can be advantageous for preserving some symmetries of the original Hamiltonian (6) at the level of the discretized path integral [13]. However, since the Trotter decomposition anyway introduces a discretization error of order O ∆τ 2 in the partition function (4) and observables (5), one can also expand the exponential in (14) to the leading order in ∆τ :…”

“…In practice, we have found that also many elements of the non-sparse matrix (14) are numerically very small, of order 10 −5 and smaller, and can be set to zero without introducing any noticeable error in the results of Monte-Carlo simulations. This allows to use sparse linear algebra to speed up the algorithm even for the exponential representation (14).…”

Accelerating Hybrid Monte Carlo simulations of the Hubbard model on the hexagonal lattice

“…All our calculations here are performed with m s = 0 . (In [25] it was pointed out that certain fermion operator discretizations with m s = 0 suffer from ergodicity issues. We stress that our discretization in eq.…”