Abstract:One bottleneck of quantum Monte Carlo (QMC) simulation of strongly correlated electron systems lies at the scaling relation of computational complexity with respect to the system sizes. For generic lattice models of interacting fermions, the best methodology at hand still scales with beta N-3 where beta is the inverse temperature and N is the system size. Such scaling behavior has greatly hampered the accessibility of the universal infrared (IR) physics of many interesting correlated electron models at (2+1)D, let alone (3+1)D. To reduce the computational complexity, we develop a new QMC method with inhomogeneous momentum-space mesh, dubbed elective-momentum ultrasize quantum Monte Carlo (EQMC) method. Instead of treating all fermionic excitations on an equal footing as in conventional QMC methods, by converting the fermion determinant into the momentum space, our method focuses on fermion modes that are directly associated with low-energy (IR) physics in the vicinity of the so-called hot spots, while other fermion modes irrelevant for universal properties are ignored. As shown in the paper, for any cutoff-independent quantities, e.g., scaling exponents, this method can achieve the same level of accuracy with orders of magnitude increase in computational efficiency. We demonstrate this method with a model of antiferromagnetic itinerant quantum critical point, realized via coupling itinerant fermions with a frustrated transverse-field Ising model on a triangle lattice. The system size of 48 x 48 x 32 (L x L x beta, almost 3 times of previous investigations) are comfortably accessed with EQMC. With much larger system sizes, the scaling exponents are unveiled with unprecedentedly high accuracy, and this result sheds new light on the open debate about the nature and the universality class of itinerant quantum critical points.

An integral algorithm of exponential observables for interacting fermions in quantum Monte Carlo simulation

High-efficiency quantum Monte Carlo algorithm for extracting entanglement entropy in interacting fermion systems

Incremental SWAP Operator for Entanglement Entropy: Application for Exponential Observables in Quantum Monte Carlo Simulation

Measuring Rényi entanglement entropy with high efficiency and precision in quantum Monte Carlo simulations

Computing entanglement entropy of interacting fermions with quantum Monte Carlo: Why we failed and how to get it right

Stable computation of entanglement entropy for 2D interacting fermion systems

Stable computation of entanglement entropy for two-dimensional interacting fermion systems

The Rényi Entanglement Entropy of a General Quantum Dimer Model at the RK Point: a Highly Efficient Algorithm.

Elective-momentum Ultrasize Quantum Monte Carlo Method

Non-equilibrium quantum Monte Carlo algorithm for stabilizer Renyi entropy in spin systems

Self-learning quantum Monte Carlo method in interacting fermion systems

Renyi Entanglement Entropy of Interacting Fermions Calculated Using Continuous-Time Quantum Monte Carlo Method

Renyi Entanglement Entropy of Interacting Fermions Calculated Using the Continuous-Time Quantum Monte Carlo Method

Resummation-based Quantum Monte Carlo for Entanglement Entropy Computation

Determinant Monte Carlo algorithm for dynamical quantities in fermionic quantum systems

Determinant Monte Carlo Algorithms for Dynamical Quantities in Fermionic Systems

Tracking the variation of entanglement Rényi negativity: an efficient quantum Monte Carlo method

Boosting Determinant Quantum Monte Carlo with Submatrix Updates: Unveiling the Phase Diagram of the 3D Hubbard Model

Post-measurement Quantum Monte Carlo

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices