Understanding correlation effects in topological phases of matter is at the forefront of current research in condensed matter physics. Here we try to clarify some subtleties in studying topological behaviors of interacting Weyl semimetals. It is well known that there exist two topological invariants defined to identify their topological character. One is the many-body Chern number, which can be directly linked to the Hall conductivity and thus to the two-particle correlations. The other is the topological index constructed from the single-particle Green's functions. Because the information on Green's functions is achieved more easily than the many-body wave functions, usually only the topological index is employed in the literature. However, the approach based on the single-particle Green's function can break down in the strongly correlated phase. For illustration, an exactly solvable two-orbital model with momentum-local two-body interactions is discussed, in which both topological invariants can be calculated analytically. We find that the topological index calculated from the Green's function formalism can be nonzero even for a nontopological strongly correlated phase with a vanishing many-body Chern number. In addition, we stress that the physical surface states implied by nonzero many-body Chern numbers should be the edge modes of particle-hole collective excitations, rather than those of a quasiparticle nature derived from the Green's function formalism. Our observations thus demonstrate the limitation of the validity of the Green's function formalism in the investigations of interacting topological materials.
Understanding correlation effects in topological phases of matter is at the forefront of current research in condensed matter physics. Here we try to clarify some subtleties in studying topological behaviors of interacting Weyl semimetals. It is well known that there exist two topological invariants defined to identify their topological character. One is the many-body Chern number, which can be directly linked to the Hall conductivity and thus to the two-particle correlations. The other is the topological index constructed from the single-particle Green's functions. Because the information on Green's functions is achieved more easily than the many-body wave functions, usually only the topological index is employed in the literature. However, the approach based on the single-particle Green's function can break down in the strongly correlated phase. For illustration, an exactly solvable two-orbital model with momentum-local two-body interactions is discussed, in which both topological invariants can be calculated analytically. We find that the topological index calculated from the Green's function formalism can be nonzero even for a nontopological strongly correlated phase with a vanishing many-body Chern number. In addition, we stress that the physical surface states implied by nonzero many-body Chern numbers should be the edge modes of particle-hole collective excitations, rather than those of a quasiparticle nature derived from the Green's function formalism. Our observations thus demonstrate the limitation of the validity of the Green's function formalism in the investigations of interacting topological materials.
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.
Luttinger's theorem constrains the particle density of interacting fermions through global properties of the single-particle Green's function, and its violation signals a breakdown of the identification between the quantized Hall response and the Green-function-based Ishikawa-Matsuyama invariant. This phenomenon becomes especially compelling in strongly correlated topological phases, such as fractional Chern insulators, where fractionalized quasiparticles lack an adiabatic connection to electrons, raising the question of how Green's-function-based topological invariants manifest in such phases. Using exact diagonalization of the fermionic Harper-Hofstadter-Hubbard model, we compute bulk single-particle Green's functions deep inside a fractional Chern insulating phase and directly evaluate the Luttinger count, its possible correction (the Luttinger integral), and the Ishikawa-Matsuyama invariant N3[G]. We demonstrate a clear violation of Luttinger's theorem and show that the fractional nature of the many-body Chern number is encoded in the Středa response of the Luttinger integral, while the integer invariant N3[G] arises from the Středa response of the Luttinger count. We also analytically prove that N3[G] is fully determined by the Luttinger count together with the Chern number of the occupied Bloch band, upon neglecting Bloch-band mixing. Finally, we propose an experimental protocol to extract all Green-function-based topological invariants from local density-of-states measurements, experimentally accessible in fractional quantum Hall systems.
Understanding how topology survives in strongly correlated systems remains a central challenge, as most topological diagnostics rely on non-interacting band structures. Here we present a framework to characterize interacting topological phases within an effective single-particle description derived from the single-particle Green's function. Using the Su-Schrieffer-Heeger model with Hatsugai-Kohmoto interactions as an analytically tractable example, we construct the one-body reduced density matrix from the Green's function and use it to define an effective winding number together with quantum volume, a measurement of state geometry. These quantities allow us to distinguish three insulating phases including correlated Mott states directly from single-particle observables. Our results show that interacting topology can be interpreted in terms of the spectral weight distribution of single-particle excitations, providing an intuitive and computationally accessible route to diagnose topological phases in correlated systems. This approach is compatible with modern many-body simulation techniques and opens a pathway toward the identification of interacting topological materials.
Weyl semimetals are three-dimensional, topologically protected, gapless phases which show exotic phenomena such as Fermi arc surface states or negative magnetoresistance. It is an open question whether interparticle interactions can turn the topological semimetal into a topologically nontrivial Mott-insulating phase. We investigate an experimentally motivated model for Weyl physics of cold atoms in optical lattices, with the main focus on interaction effects and topological properties, by means of dynamical mean-field theory. We characterize topological phases by numerically evaluating the Chern number via the Ishikawa-Matsuyama formula for interacting phases. Within our studies, we find that the Chern numbers become trivial when interactions lead to insulating behavior. For a deeper understanding of the Weyl-semimetal-to-Mott-insulator topological phase transition, we evaluate the topological properties of quasiparticle bands as well as so-called blind bands. Finally, we consider a system with an open boundary along one spatial direction in order to study correlation effects of surface states. In a narrow regime close to the topological phase transition, we find a correlation-induced state in which the surface becomes metallic while the bulk is semimetallic.
Topological phases set themselves apart from other phases since they cannot be understood in terms of the usual Landau theory of phase transitions. This fact, which is a consequence of the property that topological phase transitions can occur without breaking symmetries, is reflected in the complicated form of topological order parameters. While the mathematical classification of phases through homotopy theory is known, an intuition for the relation between phase transitions and changes to the physical system is largely inhibited by the general complexity. In this thesis we aim to get back some of this intuition by studying the properties of the Chern number (a topological order parameter) in two scenarios. First, we investigate the effect of electronic correlations on topological phases in the Green's function formalism. By developing a statistical method that averages over all possible solutions of the manybody problem, we extract general statements about the shape of the phase diagram and investigate the stability of topological phases with respect to interactions. In addition, we find that in many topological models the local approximation, which is part of many standard methods for solving the manybody lattice model, is able to produce qualitatively correct phase transitions at low to intermediate correlations. We then extend the statistical method to study the effect of the lattice, where we evaluate possible applications of standard machine learning techniques against our information theoretical approach. We define a measure for the information about particular topological phases encoded in individual lattice parameters, which allows us to construct a qualitative phase diagram that gives a more intuitive understanding of the topological phase. Finally, we discuss possible applications of our method that could facilitate the discovery of new materials with topological properties.
The aim of this series of two papers is to discuss topological invariants for interacting topological insulators (TIs). In the first paper (I), we provide a paradigm of efficient numerical evaluation scheme for topological invariants, in which we demystify the procedures and techniques employed in calculating Z2 invariant and spin Chern number via zero-frequency single-particle Green's function in quantum Monte Carlo (QMC) simulations. Here we introduce an interpolation process to overcome the ubiquitous finite-size effect, so that the calculated spin Chern number shows ideally quantized values. We also show that making use of symmetry properties of the underlying systems can greatly reduce the computational effort. To demonstrate the effectiveness of our numerical evaluation scheme, especially the interpolation process, for calculating topological invariants, we apply it on two independent two-dimensional models of interacting topological insulators. In the subsequent paper (II), we apply the scheme developed here to wider classes of models of interacting topological insulators, for which certain limitation of constructing topological invariant via single-particle Green's functions will be presented.