Lefschetz thimbles decomposition for the Hubbard model on the hexagonal lattice

Abstract

We propose a framework to study the properties of the Lefschetz thimbles decomposition for lattice fermion models approaching the thermodynamic limit. The proposed set of algorithms includes the Schur complement solver and the exact computation of the derivatives of the fermion determinant. It allows us to solve the gradient flow (GF) equations taking into account the fermion determinant exactly, with high performance. We can find both real and complex saddle points and describe the structure of the Lefschetz thimbles decomposition for large enough lattices to extrapolate the results to the thermodynamic limit. The algorithms are described for a general lattice fermion model, with emphasis on two types of lattice discretizations for relativistic fermions (staggered and Wilson), as well as on interacting tight-binding models for condensed matter systems. We apply these algorithms to the Hubbard model on a hexagonal lattice, dealing with lattice volumes as large as 12x12 with 256 steps in Euclidean time, in order to capture the properties of the thimbles decomposition as the thermodynamic, low-temperature, and continuum limits are approached. The complexity of the thimbles decomposition appears to be dependent on the form of the Hubbard-Stratonovich (HS) transformation. In particular, the evidence is provided for the existence of an optimal regime for the hexagonal lattice Hubbard model, with a reduced number of thimbles being important in the overall sum. We have performed quantum Monte Carlo (QMC) simulations using GF to deform the integration contour into the complex plane and demonstrated the agreement with exact diagonalization on small volumes (8 sites in space) if optimal setup for HS transformation is used. The residual sign problem is compared with the state-of-the-art BSS-QMC. We show that the average sign can be kept substantially higher using the Lefschetz thimbles approach.

Text substantially updated, algorithms for QCD fermions added. Accepted for publication in Phys. Rev. D

References in corpus (44)

• Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations
• Nearly Perfect Fluidity: From Cold Atomic Gases to Hot Quark Gluon Plasmas
• High density QCD on a Lefschetz thimble?
• Absence of a Spin Liquid Phase in the Hubbard Model on the Honeycomb Lattice
• Pinning the order: the nature of quantum criticality in the Hubbard model on honeycomb lattice
• Simulating full QCD at nonzero density using the complex Langevin equation
• Monte-Carlo study of the semimetal-insulator phase transition in monolayer graphene with realistic inter-electron interaction potential
• Monte Carlo simulations on the Lefschetz thimble: taming the sign problem
• Hybrid Monte Carlo on Lefschetz Thimbles -- A study of the residual sign problem
• Controlling complex Langevin dynamics at finite density
• Introductory lectures on lattice QCD at nonzero baryon number
• Strange metallicity in the doped Hubbard model
• A Monte Carlo algorithm for simulating fermions on Lefschetz thimbles
• Lefschetz-thimble analysis of the sign problem in one-site fermion model
• Structure of Lefschetz thimbles in simple fermionic systems
• Complex Langevin dynamics and zeroes of the fermion determinant
• Revisiting the Hybrid Quantum Monte Carlo Method for Hubbard and Electron-Phonon Models
• Phase diagram of the Kane-Mele-Coulomb model
• Lefschetz thimble structure in one-dimensional lattice Thirring model at finite density
• Monte-Carlo simulation of the tight-binding model of graphene with partially screened Coulomb interactions
• The ALF (Algorithms for Lattice Fermions) project release 1.0. Documentation for the auxiliary field quantum Monte Carlo code
• Complex Langevin and boundary terms
• Monte Carlo calculations of the finite density Thirring model
• Fermions at Finite Density in (2+1)d with Sign-Optimized Manifolds
• Thimble regularization at work: from toy models to chiral random matrix theories
• Hybrid-Monte-Carlo study of competing order in the extended fermionic Hubbard model on the hexagonal lattice
• Lefschetz-thimble techniques for path integral of zero-dimensional $O(n)$ sigma models
• Applying the tempered Lefschetz thimble method to the Hubbard model away from half-filling
• Complex Langevin Simulation of a Random Matrix Model at Nonzero Chemical Potential
• One-dimensional QCD in thimble regularization
• Wilson Fermion Determinant in Lattice QCD
• Finite-temperature equation of state of polarized fermions at unitarity
• Finite Density $QE{D}_{1+1}$ Near Lefschetz Thimbles
• Complex Langevin analysis of the spontaneous symmetry breaking in dimensionally reduced super Yang-Mills models
• Gauge cooling for the singular-drift problem in the complex Langevin method --a test in Random Matrix Theory for finite density QCD
• Pairing and chiral spin density wave instabilities on the honeycomb lattice: a comparative quantum Monte Carlo study
• Numerical evidence of conformal phase transition in graphene with long-range interactions
• Spin Polarized Non-Relativistic Fermions in 1+1 Dimensions
• Avoiding Ergodicity Problems in Lattice Discretizations of the Hubbard Model
• Reweighting Lefschetz Thimbles
• Taming the sign problem of the finite density Hubbard model via Lefschetz thimbles
• Schur complement solver for Quantum Monte-Carlo simulations of strongly interacting fermions
• Path integral representation for the Hubbard model with reduced number of Lefschetz thimbles
• A Hybrid-Monte-Carlo study of monolayer graphene with partially screened Coulomb interactions at finite spin density

Cited by in corpus (36)

• Machine Learning in Nuclear Physics
• Quantum Critical Points and the Sign Problem
• Applying the tempered Lefschetz thimble method to the Hubbard model away from half-filling
• Spinfoam on Lefschetz Thimble: Markov Chain Monte-Carlo Computation of Lorentzian Spinfoam Propagator
• Path integral contour deformations for observables in $SU(N)$ gauge theory
• Path integral contour deformations for noisy observables
• Worldvolume approach to the tempered Lefschetz thimble method
• Leveraging Machine Learning to Alleviate Hubbard Model Sign Problems
• Implementation of the HMC algorithm on the tempered Lefschetz thimble method
• The Antiferromagnetic Character of the Quantum Phase Transition in the Hubbard Model on the Honeycomb Lattice
• Fermion sign bounds theory in quantum Monte Carlo simulation
• Backpropagating Hybrid Monte Carlo algorithm for fast Lefschetz thimble calculations
• Sign Problem in Quantum Monte Carlo Simulation
• Mitigating the Hubbard Sign Problem with Complex-Valued Neural Networks
• Real-time lattice gauge theory actions: unitarity, convergence, and path integral contour deformations
• Simulating both parity sectors of the Hubbard Model with Tensor Networks
• Mitigating the fermion sign problem by automatic differentiation
• Deconvolving the components of the sign problem
• A density of states approach to the hexagonal Hubbard model at finite density
• Simulating Yang-Mills theories with a complex coupling
• Fermionic Sign Problem Minimization by Constant Path Integral Contour Shifts
• On the Infinite Variance Problem in Fermion Models
• Lefschetz Thimble Quantum Monte Carlo for Spin Systems
• Mitigating spikes in fermion Monte Carlo methods by reshuffling measurements
• Settling an old story: solution of the Thirring model in thimble regularization
• Instanton gas approach to the Hubbard model
• Lefschetz Thimbles and Quantum Phases in Zero-Dimensional Bosonic Models
• Single Particle Spectrum of Doped ${\mathrm{C}}_{20}{\mathrm{H}}_{12}$-Perylene
• Link auxiliary field method in the Extended Hubbard model
• A quantum Monte Carlo method on asymptotic Lefschetz thimbles for quantum spin systems: An application to the Kitaev model in a magnetic field
• Bridging the gap between numerics and experiment in free standing graphene
• Is it worth the effort to find Lefschetz thimbles? Integration contours with numerically optimal signal-to-noise ratios in simple fermionic toy models
• Stochastic and Tensor Network simulations of the Hubbard Model
• Reducing the Sign Problem with simple Contour Deformation
• The Hubbard Model on the Honeycomb Lattice with Hybrid Monte Carlo
• Negative sign free formulations of generalized Kitaev models with higher symmetries