Using exact diagonalization and quantum Monte Carlo calculations we investigate the effects of disorder on the phase diagram of both noninteracting and interacting models of two-dimensional topological insulators. In the fermion sign problem-free interacting models we study, electron-electron interactions are described by an on-site repulsive Hubbard interaction and disorder is included via the one-body hopping operators. In both the noninteracting and interacting models we make use of recent advances in highly accurate real-space numerical evaluation of topological invariants to compute phase boundaries and in the noninteracting models determine critical exponents of the transitions. We find different models exhibit distinct stability conditions of the topological phase with respect to interactions and disorder. We provide a general analytical theory that accurately predicts these trends.
Using exact diagonalization and quantum Monte Carlo calculations we investigate the effects of disorder on the phase diagram of both noninteracting and interacting models of two-dimensional topological insulators. In the fermion sign problem-free interacting models we study, electron-electron interactions are described by an on-site repulsive Hubbard interaction and disorder is included via the one-body hopping operators. In both the noninteracting and interacting models we make use of recent advances in highly accurate real-space numerical evaluation of topological invariants to compute phase boundaries and in the noninteracting models determine critical exponents of the transitions. We find different models exhibit distinct stability conditions of the topological phase with respect to interactions and disorder. We provide a general analytical theory that accurately predicts these trends.
We study the boundary criticality in 2D interacting topological insulators. Using the determinant quantum Monte Carlo method, we present a nonperturbative study of the boundary quantum phase diagram in the Kane-Mele-Hubbard-Rashba model. Our results reveal rich boundary critical phenomena at the quantum phase transition between a topological insulator and an antiferromagnetic insulator, encompassing ordinary, special, and extraordinary transitions. Combining analytical derivation of the boundary theory with unbiased numerically exact quantum Monte Carlo simulations, we demonstrate that the presence of topological edge states enriches the ordinary transition that renders a continuous boundary scaling dimension and, more intriguingly, leads to a special transition of the Berezinskii-Kosterlitz-Thouless type. Our work establishes a framework for the nonperturbative study of boundary criticality in two-dimensional topological systems with strong electron correlations.
We theoretically study topological phase transitions in four generalized versions of the Kane-Mele-Hubbard model with up to $2\ifmmode\times\else\texttimes\fi{}{18}^{2}$ sites. All models are free of the fermion-sign problem allowing numerically exact quantum Monte Carlo (QMC) calculations to be performed to extremely low temperatures. We numerically compute the Z2 invariant and spin Chern number ${C}_{\ensuremath{\sigma}}$ directly from the zero-frequency single-particle Green's functions, and study the topological phase transitions driven by the tight-binding parameters at different on-site interaction strengths. The Z2 invariant and spin Chern number, which are complementary to each another, characterize the topological phases and identify the critical points of topological phase transitions. Although the numerically determined phase boundaries are nearly identical for different system sizes, we find strong system-size dependence of the spin Chern number, where quantized values are only expected upon approaching the thermodynamic limit. For the Hubbard models we considered, the QMC results show that correlation effects lead to shifts in the phase boundaries relative to those in the noninteracting limit, without any spontaneously symmetry breaking. The interaction-induced shift is nonperturbative in the interactions and cannot be captured within a ``simple'' self-consistent calculation either, such as Hartree-Fock. Furthermore, our QMC calculations suggest that quantum fluctuations from interactions stabilize topological phases in systems where the one-body terms preserve the D3 symmetry of the lattice, and destabilize topological phases when the one-body terms break the D3 symmetry.
We study the effect of disorder and doping on the metal-insulator transition in a repulsive Hubbard model on a square lattice using the determinant quantum Monte Carlo method. First, with the aim of making our results reliable, we compute the sign problem with various parameters such as temperature, disorder, on-site interactions, and lattice size. We show that in the presence of randomness in the hopping elements, the metal-insulator transition occurs and the critical disorder strength differs at different fillings. We also demonstrate that doping is a driving force behind the metal-insulator transition.
The topological Anderson and Mott insulators are two phases that have so far been separately and widely explored beyond topological band insulators. Here we combine the two seemingly different topological phases into a system of spin-1/2 interacting fermionic atoms in a disordered optical lattice. We find that the topological Anderson and Mott insulators in the noninteracting and clean limits can be adiabatically connected without gap closing in the phase diagram of our model. Lying between the two phases, we uncover a disordered correlated topological insulator, which is induced from a trivial band insulator by the combination of disorder and interaction, as the generalization of topological Anderson insulators to the many-body interacting regime. The phase diagram is determined by computing various topological properties and confirmed by unsupervised and automated machine learning. We develop an approach to provide a unified and clear description of topological phase transitions driven by interaction and disorder. The topological phases can be detected from disorder-/interaction-induced edge excitations and charge pumping in optical lattices.
Topological phase transitions in free fermion systems can be characterized by the closing of single-particle gap and the change in topological invariants. However, in the presence of electronic interactions, topological phase transitions can be more complicated. In paper I of this series [Phys. Rev. B 93, 195163 (2016)], we have proposed an efficient scheme to evaluate the topological invariants based on the single-particle Green's function formalism. Here, in paper II, we demonstrate several interaction-driven topological phase transitions (TPTs) in two-dimensional (2D) interacting topological insulators (TIs) via large-scale quantum Monte Carlo (QMC) simulations, based on the scheme of evaluating topological invariants presented in paper I. Across these transitions, the defining symmetries of the TIs have been neither explicitly nor spontaneously broken. In the first two models, the topological invariants calculated from the Green's function formalism succeed in characterizing the topologically distinct phases and identifying interaction-driven TPTs. However, in the other two models, we find that the single-particle gap does not close and the topological invariants constructed from the single-particle Green's function acquire no change across the TPTs. Unexpected breakdown of the Green's function formalism in constructing the topological invariants is thus discovered. We thence classify the topological phase transitions in interacting TIs into two categories in practical computation: Those that have noninteracting correspondence can be characterized successfully by the topological invariants constructed from the Green's functions, while for the others that do not have noninteracting correspondence, the Green's function formalism experiences a breakdown, but more interesting and exciting phenomena, such as emergent collective critical modes at the transition, arise. Discussion on the success and breakdown of topological invariants constructed from the Green's function formalism in the context of symmetry protected topological (SPT) states is presented.
We theoretically study topological phase transitions in two generalized versions of the Kane-Mele-Hubbard model on a finite portion of the honeycomb lattice with up to 2×182 sites. Both models are free of the fermion-sign problem allowing numerically exact quantum Monte Carlo calculations to be performed in the low-temperature regime. We numerically compute the Z2 invariant and the spin-Chern number directly from the zero-frequency single-particle Greens function for different values of on-site interaction and tune tight-binding parameters to drive topological phase transitions. While phase boundaries are nearly identical for different system sizes, we find a rather strong system-size dependence of the spin Chern number, with the expected quantized value only slowly approached for the largest system sizes. Moreover, we find that quantum fluctuation effects from interactions can act both to stabilize and destabilize topological phases, depending on the details of the model.