Employing the self-learning quantum Monte Carlo algorithm, we investigate the frustrated transverse-field triangle-lattice Ising model coupled to a Fermi surface. Without fermions, the spin degrees of freedom undergoes a second-order quantum phase transition between paramagnetic and clock-ordered phases. This quantum critical point (QCP) has an emergent U(1) symmetry and thus belongs to the (2+1)D XY universality class. In the presence of fermions, spin fluctuations introduce effective interactions among fermions and distort the bare Fermi surface towards an interacting one with hot spots and Fermi pockets. Near the QCP, non-Fermi-liquid behavior are observed at the hot spots, and the QCP is rendered into a different universality with Hertz-Millis type exponents. The detailed properties of this QCP and possibly related experimental systems are also discussed.

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

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 βN 3 where β 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 × 48 × 32 (L × L × β, 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.