Abstract Interaction-driven topological phases significantly enrich the class of topological materials and thus are of great importance. Here, we study the phase diagram of interacting spinless fermions filling the two-dimensional checkerboard lattice with a quadratic band touching (QBT) point. By developing new diagnosis based on the state-of-the-art density-matrix renormalization group and exact diagonalization, we determine accurate quantum phase diagram for such a system at half-filling with three distinct phases. For weak nearest-neighboring interactions, we demonstrate the instability of the QBT towards an interaction-driven spontaneous quantum anomalous Hall (QAH) effect. For strong interactions, the system breaks the rotational symmetry realizing a nematic charge-density-wave (CDW) phase. Interestingly, for intermediate interactions we discover a symmetry-broken bond-ordered critical phase sandwiched in between the QAH and CDW phases, which splits the QBT into two Dirac points driven by interaction. Instead of the direct transition between QAH and CDW phases, our identification of an intermediate phase sheds new light on the theoretical understanding of the interaction-driven phases in QBT systems.
The interaction-driven quantum anomalous Hall (QAH) insulator has been sought for a long time in a Dirac semimetal with linear band touching points at the Fermi level. By combining exact diagonalization, density matrix renormalization group, and analytical methods, we study a spinless fermion system on the checkerboard lattice with twofold rotational symmetry, which realizes two Dirac band touching points in the absence of interaction. At weak coupling, the Dirac semimetal is stable. At a finite density-density repulsive interaction, we analyze possible symmetry broken states, and find that an QAH state is stabilized when the interaction strength exceeds the energy scale controlling the separation between the Dirac points. Through numerical simulations, we verify the existence of the QAH phase with spontaneous time-reversal symmetry breaking and quantized Chern number
C=1.
We study the quantum phases driven by interaction in a semimetal with a quadratic band touching at the Fermi level. By combining the density matrix renormalization group (DMRG), analytical power expanded Gibbs potential method, and the weak coupling renormalization group, we study a spinless fermion system on a checkerboard lattice at half-filling, which has a quadratic band touching in the absence of interaction. In the presence of strong nearest-neighbor (
V1) and next-nearest-neighbor (
V2) interactions, we identify a site nematic insulator phase, a stripe insulator phase, and a phase separation region, in agreement with the phase diagram obtained analytically in the strong coupling limit (i.e., in the absence of fermion hopping). In the intermediate interaction regime we establish a quantum anomalous Hall phase in the DMRG as evidenced by the spontaneous time-reversal symmetry breaking and the appearance of a quantized Chern number
C=1. For weak interaction we utilize the power expanded Gibbs potential method that treats
V1 and
V2 on equal footing, as well as the weak coupling renormalization group. Our analytical results reveal that not only the repulsive
V1 interaction, but also the
V2 interaction (both repulsive and attractive), can drive the quantum anomalous Hall phase. We also determine the phase boundary in the ${V}_{1}\text{\ensuremath{-}}{V}_{2}$ plane that separates the semimetal from the quantum anomalous Hall state. Finally, we show that the nematic semimetal, which was proposed for $|{V}_{2}|\ensuremath{\gg}{V}_{1}$ at weak coupling in a previous study, is absent, and the quantum anomalous Hall state is the only weak coupling instability of the spinless quadratic band touching semimetal.
We propose a general scheme for diagnosing interaction-driven topological phases in the weak interaction regime using exact diagonalization (ED). The scheme comprises the analysis of eigenvalues of the point-group operators for the many-body eigenstates and the correlation functions for physical observables to extract the symmetries of the order parameters and the topological numbers of the underlying ground states at the thermodynamic limit from a relatively small size system afforded by ED. As a concrete example, we investigate the interaction effects on the half-filled spinless fermions on the checkerboard lattice with a quadratic band crossing point. Numerical results support the existence of a spontaneous quantum anomalous Hall phase purely driven by a nearest-neighbor weak repulsive interaction, separated from a nematic Mott insulator phase at strong repulsive interaction by a first-order phase transition.
Electronic systems with quadratic band touchings, commonly found in two- and three-dimensional materials such as Bernal-stacked bilayer graphene, kagome metals, HgTe, and pyrochlore iridates, have attracted significant interest concerning the role of interactions in shaping their electronic properties. However, even in the simplest model of spinless fermions on a two-dimensional checkerboard lattice, the quantum phase diagram as a function of nearest-neighbor interaction remains under debate. We employ constrained-path quantum Monte Carlo simulations (CP-QMC) simulations to investigate the problem using a two-dimensional torus geometry. We cross-validate our results on small lattices by comparing them with density-matrix renormalization group calculations, finding quantitative agreement. In particular, we implement an improved optimization scheme within the CP-QMC simulations, enabling the identification of a bond-nematic Dirac semimetal phase that was found in tensor-network studies on cylindrical geometries, but remains inaccessible to Hartree-Fock mean-field methods. The CP-QMC approach makes it possible to establish the emergence of this phase in a geometry that preserves lattice rotational symmetry and permits extrapolation to the thermodynamic limit. Our results show that the quantum phase diagram of spinless fermions on the checkerboard lattice with nearest-neighbor repulsion features three interaction-induced phases at half filling: a quantum anomalous Hall insulator at weak coupling, a bond-nematic Dirac semimetal at intermediate coupling, and a site-nematic insulator at strong coupling.
Quadratic band touching (QBT) points are widely observed in two- and three-dimensional (2D and 3D) materials, including bilayer graphene and Luttinger semimetals, and attract significant attention from theory to experiment. However, even in its simplest form, the 2D checkerboard lattice QBT model, the phase diagram characterized by temperature and interaction strength, still remains unknown beyond the weak-coupling regime. Intense debates persist regarding the existence of various interaction-driven insulating states in this system. To address these uncertainties, we employ thermal tensor network simulations, specifically exponential tensor renormalization group and tangent space tensor renormalization group, along with density matrix renormalization group calculations to provide a comprehensive finite-temperature phase diagram for this model and shed light on previous ambiguities. Notably, our findings reveal the emergence of a robust bond-nematic Dirac semimetal phase with distinct thermodynamic properties that set it part from the nematic insulating state and other symmetry-broken states. This previously overlooked feature is found to be ubiquitous in interacting QBT systems. We also discuss the implications of these results for experimental systems such as bilayer graphene and iridate compounds.
We study the quantum phases of spinless fermions at one-third filling on a kagome lattice featuring a quadratic band touching Fermi point. In the presence of weak first and second nearest-neighbor repulsive interactions (
V1 and
V2), we demonstrate an interaction driven quantum anomalous Hall effect by employing exact diagonalization and density-matrix renormalization group methods. The time-reversal symmetry is broken spontaneously by forming loop currents that exhibit long-range correlation. Quantized Hall conductance corresponding to a Chern number of $\ifmmode\pm\else\textpm\fi{}1$ is obtained by measuring the pumped charge through inserting flux in a cylinder geometry. We find that the energy gap, which topologically protects the emerging ground states, can be enhanced remarkably by a moderate
V2<V1 via calculating the spectrum and charge excitation gaps, which highlights the experimentally feasible scheme of realizing the interaction driven topological phase by spatially decaying interactions on topologically trivial lattice models.