Kondo breakdown is one of the most intriguing problems in strongly correlated electron systems, as it is rooted in many anomalous electron behaviors found in heavy-fermion materials. In Kondo lattice systems, Kondo breakdown can arise from either strong magnetic frustrations or critical fluctuations of collective modes. Here, we reveal another type of Kondo breakdown with a fully different origin in interacting topological Kondo insulators. By employing numerically exact quantum Monte Carlo simulations, we show that with open boundary conditions, Kondo screening is inexorably destroyed by interaction effects on edges or corners in these systems. We argue that the Kondo breakdown is enforced by the symmetries of the system, because the ground states are symmetry-protected Haldane phases.
Kondo breakdown is one of the most intriguing problems in strongly correlated electron systems, as it is rooted in many anomalous electron behaviors found in heavy-fermion materials. In Kondo lattice systems, Kondo breakdown can arise from either strong magnetic frustrations or critical fluctuations of collective modes. Here, we reveal another type of Kondo breakdown with a fully different origin in interacting topological Kondo insulators. By employing numerically exact quantum Monte Carlo simulations, we show that with open boundary conditions, Kondo screening is inexorably destroyed by interaction effects on edges or corners in these systems. We argue that the Kondo breakdown is enforced by the symmetries of the system, because the ground states are symmetry-protected Haldane phases.
Topological Kondo insulators are strongly correlated materials where itinerant electrons hybridize with localized spins, giving rise to a topologically nontrivial band structure. Here, we use nonperturbative bosonization and renormalization-group techniques to study theoretically a one-dimensional topological Kondo insulator, described as a Kondo-Heisenberg model, where the Heisenberg spin-1/2 chain is coupled to a Hubbard chain through a Kondo exchange interaction in the p-wave channel (i.e., a strongly correlated version of the prototypical Tamm-Schockley model). We derive and solve renormalization-group equations at two-loop order in the Kondo parameter, and find that, at half filling, the charge degrees of freedom in the Hubbard chain acquire a Mott gap, even in the case of a noninteracting conduction band (Hubbard parameter U=0). Furthermore, at low enough temperatures, the system maps onto a spin-1/2 ladder with local ferromagnetic interactions along the rungs, effectively locking the spin degrees of freedom into a spin-1 chain with frozen charge degrees of freedom. This structure behaves as a spin-1 Haldane chain, a prototypical interacting topological spin model, and features two magnetic spin-1/2 end states for chains with open boundary conditions. Our analysis allows us to derive an insightful connection between topological Kondo insulators in one spatial dimension and the well-known physics of the Haldane chain, showing that the ground state of the former is qualitatively different from the predictions of the naive mean-field theory.
We study, by means of the density-matrix renormalization group (DMRG) technique, the evolution of the ground state in a one-dimensional topological insulator, from the noninteracting to the strongly interacting limit, where the system can be mapped onto a topological Kondo-insulator model. We focus on a toy model Hamiltonian (i.e., the interacting ``sp-ladder'' model), which could be experimentally realized in optical lattices with higher orbitals loaded with ultracold fermionic atoms. Our goal is to shed light on the emergence of the strongly interacting ground state and its topological classification as the Hubbard U interaction parameter of the model is increased. Our numerical results show that the ground state can be generically classified as a symmetry-protected topological phase of the Haldane type, even in the noninteracting case U=0 where the system can be additionally classified as a time-reversal Z2-topological insulator, and evolves adiabatically between the noninteracting and strongly interacting limits.
We investigate the ground-state properties of a recently proposed model for a topological Kondo insulator in one dimension (i.e., the p-wave Kondo-Heisenberg lattice model) by means of the density-matrix renormalization-group method. The nonstandard Kondo interaction in this model is different from the usual (i.e., local) Kondo interaction in that the localized spins couple to the ``p-wave'' spin density of conduction electrons, inducing a topologically nontrivial insulating ground state. Based on the analysis of the charge- and spin-excitation gaps, the string order parameter, and the spin profile in the ground state, we show that, at half filling and low energies, the system is in the Haldane phase and hosts topologically protected spin-1/2 end states. Beyond its intrinsic interest as a useful ``toy model'' to understand the effects of strong correlations on topological insulators, we show that the p-wave Kondo-Heisenberg model could be experimentally implemented in $p\text{\ensuremath{-}}\mathrm{band}$ optical lattices loaded with ultracold Fermi gases.
The effect of a local Kondo coupling and Hubbard interaction on the topological phase of the one-dimensional topological Kondo insulator (TKI) is numerically investigated using the infinite matrix-product state density-matrix renormalization group algorithm. The ground state of the TKI is a symmetry-protected topological (SPT) phase protected by inversion symmetry. It is found that on its own, the Hubbard interaction that tends to force fermions into a one-charge per site order is insufficient to destroy the SPT phase. However, when the local Kondo Hamiltonian term that favors a topologically trivial ground state with a one-charge per site order is introduced, the Hubbard interaction assists in the destruction of the SPT phase. This topological phase transition occurs in the charge sector where the correlation length of the charge excitation diverges while the correlation length of the spin excitation remains finite. The critical exponents, central charge, and the phase diagram separating the SPT phase from the topologically trivial phase are presented.
Distinguishing the nontrivial symmetry-protected topological (SPT) phase from the trivial insulator phase in the presence of electron-electron interaction is an urgent question to the study of topological insulators, due to the fact that most of the topological indices defined for free electron systems are very likely unsuitable for interacting cases. In this work, we demonstrate that the strange correlator is a sensitive diagnosis to detect SPT states in interacting systems. Employing large-scale quantum Monte Carlo (QMC) simulations, we investigate the interaction-driven quantum phase transition in the Kane-Mele-Hubbard model. The transition from the quantum spin Hall insulator at weak interaction to an antiferromagnetic Mott insulator at strong interaction can be readily detected by the momentum space behavior of the strange correlator in single-particle, spin, and pairing sectors. The interaction effects on the symmetry-protected edge states in various sectors, i.e., the helical Luttinger liquid behavior, are well captured in the QMC measurements of strange correlators. Moreover, we demonstrate that the strange correlator is technically easier to implement in QMC and more robust in performance than other proposed numerical diagnoses for interacting topological states, as only static correlations are needed. The attempt in this work paves the way for using the strange correlator to study interaction-driven topological phase transitions in fermionic as well as bosonic systems.
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.