Numerical stabilization of entanglement computation in auxiliary field quantum Monte Carlo simulations of interacting many-fermion systems

Abstract

In the absence of a fermion sign problem, auxiliary field (or determinantal) quantum Monte Carlo (DQMC) approaches have long been the numerical method of choice for unbiased, large-scale simulations of interacting many-fermion systems. More recently, the conceptual scope of this approach has been expanded by introducing ingenious schemes to compute entanglement entropies within its framework. On a practical level, these approaches however suffer from a variety of numerical instabilities that have largely impeded their applicability. Here we report on a number of algorithmic advances to overcome many of these numerical instabilities and significantly improve the calculation of entanglement measures in the zero-temperature projective DQMC approach, ultimately allowing to reach similar system sizes as for the computation of conventional observables. We demonstrate the applicability of this improved DQMC approach by providing an entanglement perspective on the quantum phase transition from a magnetically ordered Mott insulator to a band insulator in the bilayer square lattice Hubbard model at half filling.

References in corpus (19)

• Entanglement entropy of fermions in any dimension and the Widom conjecture
• Generalized Directed Loop Method for Quantum Monte Carlo Simulations
• Fermionic quantum criticality in honeycomb and $\pi$-flux Hubbard models: Finite-size scaling of renormalization-group-invariant observables from quantum Monte Carlo
• Sign-problem-free quantum Monte Carlo of the onset of antiferromagnetism in metals
• Fermionic Quantum Critical Point of Spinless Fermions on a Honeycomb Lattice
• Quantum Monte Carlo Study of an Interaction-Driven Band Insulator to Metal Transition
• Band insulator to Mott insulator transition in a bilayer Hubbard model
• Entanglement entropy scaling in the bilayer Heisenberg spin system
• Phase diagram of the Kane-Mele-Coulomb model
• Stable Quantum Monte Carlo Simulations for Entanglement Spectra of Interacting Fermions
• Breakup of the Fermi surface near the Mott transition in low-dimensional systems
• Breakdown of Luttinger's theorem in two-orbital Mott insulators
• Magnetic and Transport Properties of a Coupled Hubbard Bilayer with Electron and Hole Doping
• Renyi Entropies of Interacting Fermions from Determinantal Quantum Monte Carlo Simulations
• A path integral Monte Carlo method for Rényi entanglement entropies
• Non-local updates for quantum Monte Carlo simulations
• Ground state phase diagram of the half-filled bilayer Hubbard model
• Entanglement properties of the antiferromagnetic-singlet transition in the Hubbard model on bilayer square lattices
• Kane-Mele-Hubbard model on the $\pi$-flux honeycomb lattice

Cited by in corpus (1)

• Benchmark study of an auxiliary-field quantum Monte Carlo technique for the Hubbard model with shifted-discrete Hubbard-Stratonovich transformations