For $D$-dimensional weakly interacting topological insulators in certain symmetry classes, the topological invariant can be calculated from a $D$- or $(D+1)$-dimensional integration over a certain curvature function that is expressed in terms of single-particle Green's functions. Based on the divergence of curvature function at the topological phase transition, we demonstrate how a renormalization group approach circumvents these integrations and reduces the necessary calculation to that for the Green's function alone, rendering a numerically efficient tool to identify topological phase transitions in a large parameter space. The method further unveils a number of statistical aspects related to the quantum criticality in weakly interacting topological insulators, including correlation function, critical exponents, and scaling laws, that can be used to characterize the topological phase transitions driven by either interacting or noninteracting parameters. We use 1D class BDI and 2D class A Dirac models with electron-electron and electron-phonon interactions to demonstrate these principles and find that interactions may change the critical exponents of the topological insulators.
For D-dimensional weakly interacting topological insulators in certain symmetry classes, the topological invariant can be calculated from a D- or (D+1)-dimensional integration over a certain curvature function that is expressed in terms of single-particle Green's functions. Based on the divergence of curvature function at the topological phase transition, we demonstrate how a renormalization group approach circumvents these integrations and reduces the necessary calculation to that for the Green's function alone, rendering a numerically efficient tool to identify topological phase transitions in a large parameter space. The method further unveils a number of statistical aspects related to the quantum criticality in weakly interacting topological insulators, including correlation function, critical exponents, and scaling laws, that can be used to characterize the topological phase transitions driven by either interacting or noninteracting parameters. We use 1D class BDI and 2D class A Dirac models with electron-electron and electron-phonon interactions to demonstrate these principles and find that interactions may change the critical exponents of the topological insulators.
Topological phase transitions occur in real materials as well as quantum engineered systems, all of which differ greatly in terms of dimensionality, symmetries, interactions, and driving, and hence require a variety of techniques and concepts to describe their topological properties. For instance, topology may be accessed from single-particle Bloch wave functions, Green's functions, or many-body wave functions. We demonstrate that despite this diversity, all topological phase transitions display a universal feature: namely, a divergence of the curvature function that composes the topological invariant at the critical point. This feature can be exploited via a renormalization-group–like methodology to describe topological phase transitions. This approach serves to extend notions of correlation function, critical exponents, scaling laws and universality classes used in Landau theory to characterize topological phase transitions in a unified manner.
Abstract In topological insulators and topological superconductors, the discrete jump of the topological invariant upon tuning a certain system parameter defines a topological phase transition. A unified framework is employed to address the quantum criticality of the topological phase transitions in one to three spatial dimensions, which simultaneously incorporates the symmetry classification, order of band crossing, m -fold rotational symmetry, correlation functions, critical exponents, scaling laws, and renormalization group approach. We first classify higher-order Dirac models according to the time-reversal, particle-hole, and chiral symmetries, and determine the even–oddness of the order of band crossing in each symmetry class. The even–oddness further constrains the rotational symmetry m permitted in a symmetry class. Expressing the topological invariant in terms of a momentum space integration over a curvature function, the order of band crossing determines the critical exponent of the curvature function, as well as that of the Wannier state correlation function introduced through the Fourier transform of the curvature function. The conservation of topological invariant further yields a scaling law between critical exponents. In addition, a renormalization group approach based on deforming the curvature function is demonstrated for all dimensions and symmetry classes. Through clarification of how the critical quantities, including the jump of the topological invariant and critical exponents, depend on the nonspatial and the rotational symmetry, our work introduces the notion of universality class into the description of topological phase transitions.
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.
A recently proposed curvature renormalization group scheme for topological phase transitions defines a generic ``curvature function'' as a function of the parameters of the theory and shows that topological phase transitions are signaled by the divergence of this function at certain parameter values, called critical points, in analogy with usual phase transitions. A renormalization group procedure was also introduced as a way of flowing away from the critical point toward a fixed point, where an appropriately defined correlation function goes to zero and topological quantum numbers characterizing the phase are easy to compute. In this paper, using two independent models, a model in the AIII symmetry class and a model in the BDI symmetry class, in one dimension as examples, we show that there are cases where the fixed-point curve and the critical-point curve appear to intersect, which turn out to be multicritical points, and focus on understanding its implications.
We present a Green's function formalism to investigate the topological properties of weakly interacting one-dimensional topological insulators, including the bulk-edge correspondence and the quantum criticality near topological phase transitions, and using the interacting Su-Schrieffer-Heeger model as an example. From the many-body spectral function, we find that closing of the bulk gap remains a defining feature even if the topological phase transition is driven by interactions. The existence of edge state in the presence of interactions can be captured by means of a T-matrix formalism combined with Dyson's equation, and the bulk-edge correspondence is shown to be satisfied even in the presence of interactions. The critical exponent of the edge state decay length is shown to be affiliated with the same universality class as the noninteracting limit.
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.