EMUS-QMC: Elective Momentum Ultra-Size Quantum Monte Carlo Method

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 ultra-size 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 manuscript, 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\times 48\times 32$ ($L\times L\times \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.

9 pages, 4 figures

References in corpus (18)

• Fermi-liquid instabilities at magnetic quantum phase transitions
• Quantum spin-liquid emerging in two-dimensional correlated Dirac fermions
• Quantum phase transitions of metals in two spatial dimensions: II. Spin density wave order
• Are non-Fermi-liquids stable to Cooper pairing?
• Self-Learning Monte Carlo Method
• Superconductivity in the vicinity of antiferromagnetic order in CrAs
• Critical fermi surfaces and non-fermi liquid metals
• Pressure induced superconductivity on the border of magnetic order in MnP
• Four-loop critical exponents for the Gross-Neveu-Yukawa models
• Sign-problem-free quantum Monte Carlo of the onset of antiferromagnetism in metals
• Interplay of quantum and thermal fluctuations in a frustrated magnet
• Instabilities near the onset of spin density wave order in metals
• Critical behavior of Dirac fermions from perturbative renormalization
• A simple fermionic model of deconfined phases and phase transitions
• Quantum critical properties of a metallic spin density wave transition
• Symmetry Enforced Self-Learning Monte Carlo Method Applied to the Holstein Model
• Anisotropic Non-Fermi Liquids
• Effect of nematic order on the low-energy spin fluctuations in detwinned BaFe${}_{1.935}$Ni${}_{0.065}$As${}_2$

Cited by in corpus (13)

• Quantum Monte Carlo simulation of the chiral Heisenberg Gross-Neveu-Yukawa phase transition with a single Dirac cone
• Momentum space quantum Monte Carlo on twisted bilayer Graphene
• Itinerant Quantum Critical Point with Fermion Pockets and Hot Spots
• Designer Monte Carlo Simulation for Gross-Neveu Transition
• Revealing Fermionic Quantum Criticality from New Monte Carlo Techniques
• Dynamical properties of collective excitations in twisted bilayer Graphene
• Identification of non-Fermi liquid fermionic self-energy from quantum Monte Carlo data
• Nematic Quantum Criticality in Dirac Systems
• Emus live on the Gross-Neveu-Yukawa archipelago
• Fermi-liquid ground state of interacting Dirac fermions in two dimensions
• A Sport and a Pastime: Model Design and Computation in Quantum Many-Body Systems
• Network-Initialized Monte Carlo Based on Generative Neural Networks
• Quantifying the fragility of unprotected quadratic band crossing points