Research ArticlePHYSICS

Numerical observation of emergent spacetime supersymmetry at quantum criticality

See allHide authors and affiliations

Science Advances  02 Nov 2018:
Vol. 4, no. 11, eaau1463
DOI: 10.1126/sciadv.aau1463


No definitive evidence of spacetime supersymmetry (SUSY) that transmutes fermions into bosons and vice versa has been revealed in nature so far. Moreover, the question of whether spacetime SUSY in 2 + 1 and higher dimensions can emerge in generic lattice microscopic models remains open. Here, we introduce a lattice realization of a single Dirac fermion in 2 + 1 dimensions with attractive interactions that preserves both time-reversal and chiral symmetries. By performing sign problem–free determinant quantum Monte Carlo simulations, we show that an interacting single Dirac fermion in 2 + 1 dimensions features a superconducting quantum critical point (QCP). We demonstrate that the Embedded Image spacetime SUSY in 2 + 1 dimensions emerges at the superconducting QCP by showing that the fermions and bosons have identical anomalous dimensions 1/3, a hallmark of the emergent SUSY. We further show some experimental signatures that may be measured to test such emergent SUSY in candidate systems.


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 (114), 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 Embedded Image 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 (1720), 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 (2532).

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 Embedded Image 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 Embedded Image 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 Embedded Image 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 Embedded Image SUSY, which can be measured in experiments such as scanning tunneling microscopy (STM) to test the predicted SUSY.

Fig. 1 The phase diagram of a single-flavor Dirac fermion with interactions.

From sign problem–free start-of-the-art QMC simulations, we show that the Embedded Image spacetime SUSY emerges at the SC QCP. The QMC calculation is performed at zero temperature, and the phase diagram at finite temperature is schematic.


The single Dirac fermion model

To regularize a single Dirac fermion on the square lattice while preserving both time-reversal and chiral symmetries, we consider the following single-particle HamiltonianEmbedded Image(1)where Embedded Image creates an electron at site i with spin polarization σ = ↑/↓, and tij = tR (R = rirj) is the hopping amplitude given byEmbedded Image(2)with Embedded Image, Rα = 1, ⋯, Lα − 1, and Lα denoting the number of sites along the α = x, y direction of the square lattice. This type of regularized single Dirac fermion on the lattice was dubbed “SLAC” fermions (35). Note that the feature of hopping only along either the x or y direction is not essential, and appropriate hopping along other directions can be added without qualitatively changing the main physics discussed below. We now show that the above lattice model (or SLAC fermion) satisfies all of the requirements expected for a single Dirac fermion for all practical reasons.

It is straightforward to perform a Fourier transform from real to to momentum space and obtainEmbedded Image(3)where ψp = (cp, cp)T, with cpσ annihilating a fermion with momentum p = (px, py) and spin σ, vF is the Fermi velocity (we set vF = 1 hereafter), and σa denotes Pauli matrices. Note that the momentum eigenvalues p run over the first Brillouin zone and are quantized as Embedded Image for periodic boundary conditions. It is clear that the lattice model has a single Dirac point at p = 0 (namely, Γ point) with a linear dispersion all the way to the edges of the first Brillouin zone, as shown in Fig. 2A. Moreover, it can be easily verified that the model does not vary under both time-reversal and chiral symmetries. Note that our model does not directly contradict the fermion-doubling theorem because the hopping here is not local. The hopping amplitudes decay as 1/r at long distance.

Fig. 2 The band structure of lattice fermions with a single Dirac cone.

(A) Energy dispersion for a massless lattice Dirac Hamiltonian. The energy dispersion is perfectly linear and given by Embedded Image. (B) Berry flux of the same model with a finite mass m = 1. Although the total Berry phase vanishes, we obtain θB ≃ π upon excluding the boundaries of the Brillouin zone. (C) We obtain a nearly gapless left-moving (right-moving) boundary state near x = 0 (x = L/2) by creating a domain wall in the mass term. The results are obtained for |m| = 1, and Lx = Ly = L = 41. The gap at ky = 0 is nonzero and equal to 0.0483. (D) Scaling of the edge-state gap with system size for |m| = 0.1. This plot implies that ΔE ∝ 1/L when m → 0.

Besides linear dispersion around the single Dirac point, the lattice model above also exhibits most of the other physical properties expected for Dirac fermions such as π Berry phase around the Dirac point and chiral edge states along mass domain walls. By considering the mass term in the lattice model, namely, Embedded Image, it is straightforward to verify that the lattice model gives a Berry phase, which is sgn(m)π for the whole Brillouin zone excluding its boundaries, as shown in Fig. 2B. However, the total Berry phase vanishes due to the − sgn(m)π contribution of the Brillouin zone boundaries. Although this observation seemingly implies the absence of protected zero modes and gapless edge states |m| > 0 according to the Atiyah-Singer’s index theorem, we shall show below that the edge states along domain walls are nearly gapless with a tiny gap that vanishes as 1/L, where L is the distance between two domain walls.

We now explicitly consider a domain wall for the mass term along the x direction and periodic boundary condition along the y direction. The local mass term has the following profile: m(x < Lx/2) = m0 and m(xLx/2) = − m0, where m0 is a finite constant. In Fig. 2C, the energy eigenvalues are plotted against ky. Two nearly gapless modes with opposite chiralities appear because of the presence of two domain walls. A direct examination of the single-particle wave functions reveals that the chiral (anti-chiral) branch of edge states is localized around Embedded Image (x = Lx). For finite Lx, because of the nonzero value of the direct hopping between the two domain walls, the edge states exhibit a tiny gap that decays to zero algebraically (∝ 1/Lx), as shown in Fig. 2D. The emergence of the chiral modes along boundaries implies that Embedded Image. Recall that m > 0 and m < 0 regions are time-reversal partners and thus must have opposite Chern numbers, namely, Embedded Image. We thus obtain Embedded Image. Therefore, for all practical reasons, the effective Chern number of the above hopping model can be considered as sgn(m)/2, similar to the surface of 3D topological insulators.

SC quantum criticality

Having shown that the regularized lattice model exhibits almost all physical properties of a single Dirac fermion, we are ready to consider interactions in such a system with the following HamiltonianEmbedded Image(4)where U denotes the strength of onsite Hubbard interactions and Embedded Image. With the onsite Hubbard interactions, the model still respects the particle-hole symmetry such that the system stays at half-filling. When the Hubbard interaction is attractive, namely, U < 0, this model is sign problem free in QMC (the details of QMC are discussed in the supplementary materials). Consequently, the interacting effects can be investigated by large-scale numerically reliable QMC simulation. Here, we use projector QMC in the Majorana representation to study the ground-state properties as well as the nature of quantum phase transitions of the model in Eq. 4 with attractive Hubbard interaction. (Previously, quantum criticality was only studied by QMC in 2D Dirac semimetals with even number of Dirac cones.)

It is expected that singlet SC pairing can be generated when the attractive Hubbard interaction is sufficiently strong. To study the quantum phase transition into the putative SC phase, we calculate the structure factor of onsite singlet pairing on a system with size Embedded Image, where Δj = cjcj. The SC long-range order can be extracted through finite-size scaling (FSS) Embedded Image. Besides, we also measure the quasiparticle excitation gap from time-dependent Green’s function. From the state-of-the-art QMC simulations (shown in the supplementary materials), we show that the SC order parameter is finite and the single-particle gap is opened when the Hubbard interaction exceeds a critical value.

To accurately identify the QCP, we evaluate the renormalization group (RG)-invariant quantity Binder ratio that is independent of the system sizes at the critical point. The Binder ratio is defined as Embedded Image, where Embedded Image and Embedded Image for a system with N = L2 sites. The quantum phase transition point is identified as the crossing point of the Binder ratio for different system sizes L. The results for the Binder ratio, as shown in Fig. 3A, convincingly demonstrate that there is a quantum phase transition from the Dirac semimetal phase to the SC phase occurring at U = Uc ≈ − 0.83 (in unit of the bandwidth). In the SC phase, our QMC calculations show evidence of expected Goldstone modes and Higgs bosons.

Fig. 3 The QMC results of SC quantum criticality in a single Dirac fermion.

(A) The Binder ratio results show that SC phase transition occurs at Uc ≈ 0.83. (B) From the structure factor of the SC order parameter plotted versus L, we obtain the boson anomalous dimension ηb = 0.32 ± 0.02. (C) From the fermion correlation function at largest separation Embedded Image plotted versus linear system size L, we obtain fermion anomalous dimension ηf = 0.34 ± 0.05. (D) By using data collapse analysis of SC structure factors near U = Uc for L = 11, ⋯, 21, we obtain the transition point Uc = 0.827 and critical exponents ηb = 0.32 ± 0.02 and ν = 0.87 ± 0.05.

Emergent 2 + 1D spacetime SUSY

At the SC QCP U = Uc, the system features a single Dirac fermion mode as well as a single complex boson (here, the complex boson is the SC order parameter fluctuation). It was argued from the perturbative renormalization group analysis in 4 − ε spacetime dimensions that a 2 + 1D Embedded Image SUSY might emerge by setting ε = 1 (57). However, it is not known a priori whether such spacetime SUSY can emerge in a microscopic model at the QCP and nonperturbative methods such as QMC are needed to address this unambiguously.

If the 2 + 1D Embedded Image SUSY emerges at the SC QCP, the anomalous dimensions of fermions and bosons at the QCP should be identical and are equal to 1/3, namely, Embedded Image. The equivalence of fermion and boson anomalous dimensions is a hallmark of SUSY. To verify whether the SC QCP in our model features an emergent spacetime SUSY, we study the critical properties of this quantum phase transition systematically through FSS analysis (the details of FSS are shown in the supplementary materials). The anomalous dimensions of the boson and fermion can be extracted via the correlation functions Embedded Image and Embedded Image according to the definition of anomalous dimensions. Here, Embedded Image is the largest separation between two sites in the system.

The bosonic and fermionic correlation functions at the QCP are shown in Fig. 3 (B and C). The anomalous dimensions of the boson and fermion are equal to each other within error bar: ηb = 0.32 ± 0.02 and ηf = 0.34 ± 0.05. Moreover, the values of the bosonic and fermionic anomalous dimensions obtained from QMC are consistent with the exact result of 1/3 associated with the 2 + 1D Embedded Image SUSY. These results provide convincing evidence that the SC QCP in our regularized lattice model features emergent SUSY. In the supplementary materials, we present the QMC results of bosonic and fermionic correlation functions away from the QCP. It is clear that the boson and fermion anomalous dimensions in the disordered phase away from the QCP are not equal to each other and deviate from 1/3, which indicates that the consistency of the anomalous dimensions between our QMC result and the exact result associated with the Embedded Image SUSY theory only occurs at the SC QCP. The consistency between the anomalous dimension of the bosons in the model and the one in the SUSY theory is further supported by the results of the data collapse analysis, as shown in Fig. 3D. In addition, from the data collapse analysis, we extract the correlation-length critical exponent ν = 0.87 ± 0.05, which is consistent with the nearly exact result of 0.917 obtained from the conformal bootstrap calculation of the Embedded Image SUSY theory in 2 + 1D (34) and with the results from RG calculations (36, 37). This again indicates that the SC QCP in the interacting quantum model features the emergent spacetime SUSY.

Experimental signatures

The SC QCP of a single Dirac fermion is potentially observable in realistic materials, such as the surface of 3D topological insulators. There are various experimental ways to check the putative emergent SUSY at the SC QCP. For instance, at the QCP, the zero-temperature optical conductivity Embedded Image, where K is a constant known exactly due to the emergent SUSY (38). Moreover, the SUSY dictates that the local density of states (LDOS) ρ(ω) of electrons satisfies the scaling law Embedded Image for ω ≪ 1, which can be measured by STM in experiments. In our model, the LDOS can be calculated by evaluating imaginary-time single-particle Green’s function Gf(τ) = 〈ci(0)ci(τ)〉 in QMC simulations and then performing analytical continuation (see the supplementary materials for details). As shown in Fig. 4A, Gf(τ) obtained from QMC simulations at the SC QCP behaves as Embedded Image with the exponent α = 2.34 ± 0.02, which is consistent with the one in the Embedded Image SUSY. Moreover, by analytical continuation, we obtain the LDOS at SC QCP: ρ(ω) ∝ |ω|a with a ≈ 1.37 ± 0.07, as shown in Fig. 4B. This scaling of LDOS is consistent with the exact result of 4/3 given by the Embedded Image SUSY within error bar. The LDOS can be measured by STM measurements to experimentally test the emergent SUSY.

Fig. 4 The QMC results of LDOS and unequal-time single-particle Green’s function at SC QCP.

(A) From the slope of linear fitting in the ln-ln plot of the Green’s function versus imaginary time τ, we obtain Embedded Image, with α = 2.34 ± 0.02. (B) From the slope of linear fitting in the ln-ln plot of LDOS versus frequency, we obtain the LDOS ρ(ω) ∝ |ω|a, with a = 1.37 ± 0.07.


The emergent SUSY observed at the SC QCP in the 2D microscopic model above suggests that the microscopic model can capture all essential physics of a single Dirac cone in 2 + 1D. In particular, it may be used to investigate novel properties of a single Dirac cone on the surface of 3D interacting topological insulators, such as non-Abelian Majorana zero modes at magnetic vortex cones when the single Dirac fermion is SC (39). Moreover, it has been recently argued that it is possible to gap out the single Dirac cone surface states of 3D topological insulators without breaking any symmetry through strong interactions, and the resulting exotic gapped ground state exhibits nontrivial topological order (1720). One of the approaches to justifying the quantum phase transition involves disordering the time-reversal symmetric Fu-Kane state via the multiple-vortex proliferation mechanism. We think that adding the interaction Embedded Image with V > 0 in Eq. 4 can destroy the SC phase coherence. It would be interesting to study in the future whether such regularized models of a single Dirac fermion can realize the nontrivial surface topological order.

Our work also motivates further studies of other types of 2 + 1D spacetime SUSY in microscopic models by nonperturbative methods. For instance, a sign problem–free microscopic lattice model similar to the one in the present work may be constructed for a single 2 + 1D massless Majorana fermion that can emerge on the surface of 3 + 1D topological superconductors. Strong short-range interactions can gap out Majorana fermions by breaking the time-reversal symmetry, and its QCP may realize an Embedded Image SUSY (5).

We now discuss distinctions between our work and previous simulations of the supersymmetric Wess-Zumino model [see, e.g., (40)]. Previous works mainly investigated lattice regularization of supersymmetric field theories. In contrast, we consider the microscopic lattice model that can mimic the surface of 3D topological insulators with attractive interactions and then investigate whether spacetime SUSY emerges at the SC quantum criticality. Moreover, lattice models studied in previous works are all sign problematic. In contrast, the microscopic lattice model proposed in the present work is sign problem free such that we are able to perform large-scale QMC simulation to reliably verify whether spacetime SUSY emerges at long distance and low energy.

Our unbiased and numerically reliable simulations of the microscopic quantum model of a single Dirac cone have provided convincing evidence of emergent Embedded Image spacetime SUSY in 2 + 1 dimensions at the SC QCP. The results presented here can lend concrete support to potentially realize emergent spacetime SUSY in quantum materials such as the surface of 3D topological insulators, e.g., Bi2Se3. If realized experimentally, it will not only shed light on the intriguing interplay between topology and symmetry but also provide a promising arena to explore SUSY as well as its spontaneous breaking.


Quantum Monte Carlo

We used projector QMC to investigate the ground-state properties of the model of single Dirac fermion described by the Hamiltonian in Eq. 1 with attractive Hubbard interaction. In the projector QMC, the expectation value of an observable O in the ground state can be evaluated as Embedded Image, where ψ0 is the true ground-state wave function and |ψT〉 is a trial wave function that should have a finite overlap with the true ground-state wave function. Note that Θ is a projection parameter in the simulation. Although Θ → ∞ is needed to reach the exact ground state, in numerical calculations a sufficient large Θ works for practical purposes of obtaining physical quantities with required accuracy. Because of the absence of sign problem, we can perform large-scale QMC simulations with large system sizes and sufficiently large Θ values. In our QMC simulation, we used periodic boundary condition on the square lattice L × L with largest L = 21. The imaginary-time projection parameter is 2Θ = 60/t for most systems in the calculation. In the calculation of single-particle gap, the systems with large sizes were computed using 2Θ = 70/t. We have checked that all the results stayed nearly the same when larger Θ values were used, which ensured desired convergence to the limit of Θ → ∞.

Because the computational complexity of the determinant QMC algorithm scales in cubic power with the system’s volume, it is highly challenging to simulate a fermionic lattice model with a very large size by QMC. Our sign problem–free QMC simulation studied the fermionic model with the largest size accessible within a reasonable computation time. Note that, although numerical results obtained from simulating fermionic models with finite system sizes could have a certain degree of uncertainty, the numerical evidences of emergent SUSY provided by our state-of-the-art QMC simulations are already convincing. It is partly because the raw data obtained from determinant QMC typically have a small statistical error. Moreover, the results of critical exponents extracted from different FSS methods are consistent with each other within error bars, which also indicate that the critical exponents obtained by our QMC simulations are reliable.


Supplementary material for this article is available at

Supplementary Text

Fig. S1. The QMC results of superconductivity structure factors and single-particle gaps.

Fig. S2. The scaling behaviors of superconductivity structure factors and fermion correlation away from the QCP.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We would like to thank S. Kivelson and S. Shenker for helpful discussions. Funding: This work was supported, in part, by the MOST of China under grant nos. 2016YFA0301001 and 2018YFA0305604 (H.Y.) and by the NSFC under grant no. 11474175 (Z.-X.L. and H.Y.). A.V. was funded by the Gordon and Betty Moore Foundation’s EPiQS Initiative through grant GBMF4302, and C.B.M. acknowledges support from the DOE Office of Basic Energy Sciences under grant no. DE-AC02-76SF00515 and from the Alexander von Humboldt Foundation. Author contributions: Z.-X.L., A.V., C.B.M., and H.Y. performed research, analyzed data, and wrote the paper. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article