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.

Improvements to the Stochastic Series Expansion method for the $JQ_2$ model with a magnetic field

High temperature series expansions of S = 1/2 Heisenberg spin models: algorithm to include the magnetic field with optimized complexity

Quantum Monte Carlo simulations of highly frustrated magnets in a cluster basis: The two-dimensional Shastry-Sutherland model

Elective-momentum Ultrasize Quantum Monte Carlo Method

Non-perturbative Many-Body Treatment of Molecular Magnets

Magnetic Behaviors of Classical Spin Model on the Shastry–Sutherland Lattice: Monte Carlo Simulation

Rapidly mixing loop representation quantum Monte Carlo for Heisenberg models on star-like bipartite graphs

Two-dimensional coherent spectrum of high-spin models via a quantum computing approach

Post-measurement Quantum Monte Carlo

Magnon, doublon and quarton excitations in 2D S=1/2 trimerized Heisenberg models

Loop ensembles in Stochastic Series Expansion of Two-Dimensional Heisenberg Antiferromagnets

Hierarchical equations of motion approach for accurate characterization of spin excitations in quantum impurity systems

Loop Algorithm for Quantum Transverse Ising Model in a Longitudinal Field

Quantum Monte Carlo study of long-range transverse-field Ising models on the triangular lattice

Monte Carlo Based Techniques for Quantum Magnets with Long-Range Interactions

Competition between spin excitation and Kondo correlation in magnetic molecular junctions: Theoretical insight from first-principles-based simulations

A quantum Monte Carlo algorithm for arbitrary spin-1/2 Hamiltonians

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

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

Quantum channel locally interacting with environment