Metallic quantum criticality is among the central themes in the understanding of correlated electronic systems, and converging results between analytical and numerical approaches are still under review. In this work, we develop a state-of-the-art large-scale quantum Monte Carlo simulation technique and systematically investigate the itinerant quantum critical point on a 2D square lattice with antiferromagnetic spin fluctuations at wavevector Q=(π,π)—a problem that resembles the Fermi surface setup and low-energy antiferromagnetic fluctuations in high-Tc cuprates and other critical metals, which might be relevant to their non–Fermi-liquid behaviors. System sizes of 60×60×320 (L×L×Lτ) are comfortably accessed, and the quantum critical scaling behaviors are revealed with unprecedented high precision. We found that the antiferromagnetic spin fluctuations introduce effective interactions among fermions and the fermions in return render the bare bosonic critical point into a different universality, different from both the bare Ising universality class and the Hertz–Mills–Moriya RPA prediction. At the quantum critical point, a finite anomalous dimension η∼0.125 is observed in the bosonic propagator, and fermions at hotspots evolve into a non-Fermi liquid. In the antiferromagnetically ordered metallic phase, fermion pockets are observed as the energy gap opens up at the hotspots. These results bridge the recent theoretical and numerical developments in metallic quantum criticality and can serve as the stepping stone toward final understanding of the 2D correlated fermions interacting with gapless critical excitations.

This review summarizes recent developments in the study of fermionic quantum criticality, focusing on new progress in numerical methodologies, especially quantum Monte Carlo methods, and insights that emerged from recently large-scale numerical simulations. Quantum critical phenomena in fermionic systems have attracted decades of extensive research efforts, partially lured by their exotic properties and potential technology applications and partially awaked by the profound and universal fundamental principles that govern these quantum critical systems. Due to the complex and non-perturbative nature, these systems belong to the most difficult and challenging problems in the study of modern condensed matter physics, and many important fundamental problems remain open. Recently, new developments in model design and algorithm improvements enabled unbiased large-scale numerical solutions to be achieved in the close vicinity of these quantum critical points, which paves a new pathway towards achieving controlled conclusions through combined efforts of theoretical and numerical studies, as well as possible theoretical guidance for experiments in heavy-fermion compounds, Cu-based and Fe-based superconductors, ultra-cold fermionic atomic gas, twisted graphene layers, etc., where signatures of fermionic quantum criticality exist. CONTENTS

We report an implementation of the momentum space quantum Monte Carlo (QMC) method on the interaction model for the twisted bilayer graphene (TBG). The long-range Coulomb repulsion is treated exactly with the flat bands, spin and valley degrees of freedom of electrons taking into account. We prove the absence of the minus sign problem for QMC simulation when either the two valleys or the two spin degrees of freedom are considered. By taking the realistic parameters of the twist angle and interlayer tunnelings into the simulation, we benchmark the QMC data with the exact band gap obtained at the chiral limit, to reveal the insulating ground states at the charge neutrality point (CNP). Then, with the exact Green’s functions from QMC, we perform stochastic analytic continuation to obtain the first set of single-particle spectral function for the TBG model at CNP. Our momentum space QMC scheme therefore offers the controlled computation pathway for systematic investigation of the electronic states in realistic TBG model at various electron fillings.