Spacetime supersymmetry (SUSY) was originally proposed as a fundamental symmetry of nature more than four decades ago, but no experimental evidence of SUSY in particle physics has been confirmed. Recently, it has been theoretically argued that SUSY can also spontaneously emerge in certain condensed matter systems (1–14), e.g., near the superconducting (SC) quantum critical point (QCP) of an interacting single-flavor Dirac fermions in 2 + 1–dimensional (2 + 1D) systems (5, 6). However, verification of this fascinating SUSY of a single Dirac fermion in microscopic lattice models in 2 + 1D by nonperturbative and unbiased approaches is still lacking and is thus highly desired.
Dirac fermions are essential ingredients of modern physics that can appear as either elementary particles such as electrons and positrons or emergent quasiparticles, e.g., massless Dirac fermions in graphene and on the surface of 3D topological insulators (15, 16). For a single flavor of massless interacting Dirac fermion in 2 + 1 dimensions, there are numerous interesting phenomena and theoretical predictions, from emergent spacetime SUSY at the SC QCP (5, 6) to the surface topological order (17–20), as well as fermion dualities (21). Although a single Dirac cone can occur on the surface of 3D topological insulators, studying such interacting problems in 2 + 1D microscopic models has been highly challenging due to the notorious no-go theorem of fermion doubling. According to this theorem, it is impossible to realize a single Dirac fermion in local lattice models in two spatial dimensions while maintaining time-reversal and chiral symmetries. Usual lattice regularization of a single-flavor Dirac fermion violates some of those symmetry requirements such that existing approaches cannot reveal many fascinating features associated with a single Dirac fermion.
In this study, we investigate a novel 2D lattice model of spin-1/2 fermions that features a single Dirac point at Γ, with perfectly linear dispersion and quantized π Berry phase around Γ, and preserves both time-reversal and chiral symmetries. Fermions in this model can hop along either the x or y direction, with hopping amplitudes that decay in power law at long distances. At half-filling, namely, when the Fermi level is exactly at the neutral point of single Dirac cone, sufficiently strong attractions between fermions should induce superconductivity in the system. If the lattice regularization can capture low-energy physics of a single Dirac cone, spacetime SUSY could emerge at the SC QCP. Consequently, it is highly desired to investigate universal properties of this putative SC quantum phase transition by a reliable and nonperturbative method like quantum Monte Carlo (QMC) (22) without encountering the fermion sign problem (23, 24). However, QMC methods are sign problem free only for limited classes of interacting models (25–32).
Our lattice model of a single Dirac cone with onsite Hubbard attractive interaction U is sign problem free, which allows us to study the emergent behaviors of the SC quantum phase transition in a numerically reliable way. From the state-of-the-art QMC simulations, we provide convincing evidence that the spacetime SUSY emerges at the SC QCP, as shown schematically in Fig. 1. First, the fermions and order parameter bosons have identical anomalous dimensions that are consistent with the exact value of 1/3 (33) associated with the SUSY. Moreover, we obtain the correlation-length exponent ν = 0.87 ± 0.05, which is consistent with the nearly exact result of 0.917 obtained from conformal bootstrap calculations (34) of the SUSY in 2 + 1 dimensions. Moreover, our QMC calculations show that the local electronic density of states ρ(ω) at ω ≪ 1 behaves like ρ(ω) ∝ ωa with the exponent a = 1.37 ± 0.07, close to the exact value of 4/3 associated with the SUSY, which can be measured in experiments such as scanning tunneling microscopy (STM) to test the predicted SUSY.