# The topological Anderson insulator phase in the Kane-Mele model

###### Abstract

It has been proposed that adding disorder to a topologically trivial mercury telluride/cadmium telluride (HgTe/CdTe) quantum well can induce a transition to a topologically nontrivial state. The resulting state was termed topological Anderson insulator and was found in computer simulations of the Bernevig-Hughes-Zhang model. Here, we show that the topological Anderson insulator is a more universal phenomenon and also appears in the Kane-Mele model of topological insulators on a honeycomb lattice. We numerically investigate the interplay of the relevant parameters, and establish the parameter range in which the topological Anderson insulator exists. A staggered sublattice potential turns out to be a necessary condition for the transition to the topological Anderson insulator. For weak enough disorder, a calculation based on the lowest-order Born approximation reproduces quantitatively the numerical data. Our results thus considerably increase the number of candidate materials for the topological Anderson insulator phase.

###### pacs:

73.20.Fz,03.65.Vf,73.43.NqTopological insulators (TIs) are novel materials which have raised a great deal of interest over the past decade.hasan10 (); qi11 () One of their distinguishing features is the existence of conducting boundary states together with an insulating bulk. The boundary states are protected by time-reversal symmetry (TRS) and exist both in two-dimensional (2D) and three-dimensional (3D) TIs. In 2D TIs, the boundary states lead to an edge conductance of one conductance quantum per edge for chemical potentials inside the bulk band gap.bernevig06 (); kane05 (); kane05b ()

It is a challenging task to find candidate materials for TIs. So far, only a limited number of materials are known. The most prominent 2D TIs are HgTe/(Hg,Cd)Te quantum wells (HgTeQWs) koenig07 () and InAs/GaSb heterostructures,knez11 (); li15 () whereas 3D TIs were found for instance in BiSb. hsieh08 () The fact that their metallic surface states emerge due to a topological property of the bulk band structure means that they are robust to weak disorder. However, one expects that a large amount of disorder should ultimately localize the surface states and render them insulating.

All the more surprising, it was predicted that the opposite transition can happen in certain parameter ranges: adding strong disorder can convert a trivial insulator without edge states into a topological insulator with perfectly conducting edge states. Materials that exhibit this new state have been termed topological Anderson insulators (TAIs).

This effect was first theoretically predicted based on the lattice version of the Bernevig-Hughes-Zhang (BHZ) model for HgTeQWs in the presence of Anderson disorder.Li09 (); jiang09 () For Anderson disorder, a random on-site potential, uniformly distributed in an energy window of width , is assigned to each lattice site of a tight-binding model. From Anderson’s theory of localization anderson58 () one expects that a system with finite conductance without disorder undergoes a transition to a system with localized states and suppressed conductance as the disorder is increased beyond a certain threshold value. The behavior of TAIs instead is quite different. A TAI is an ordinary band insulator in the clean limit. Above a critical disorder strength , an interesting topological state appears, in which the material features a quantized conductance. For even stronger , above the disorder strength at which the states of the conduction and valence band localize, it was proposed that tunneling across the bulk becomes possible,chen12 () probably enabled by percolating states,girschik15 () and the conductance is again suppressed.

The disorder-induced transition can be understood by a renormalization of the model parameters. The BHZ model with disorder and band mass can be approximated by an effective model of a clean system and renormalized mass . Using an effective-medium theory and the self-consistent Born approximation (SCBA), it was shown that for certain model parameters, can become negative even if the bare mass is positive.groth09 () As a consequence, the effective model becomes that of a TI and features edge states with a quantized conductance of .prodan11 ()

Furthermore, TAIs have been predicted in several related systems, for instance in a honeycomb lattice described by the time-reversal-symmetry breaking Haldane model,xing11 () a modified Dirac model,xing11 () the BHZ model with non-conserving spin-orbit coupling,yamakage11 () as well as in 3D topological insulators.guo10 () Moreover, similar transitions from a topologically trivial to a topologically nontrivial phase have been found to be generated by periodically varying potentials fu14 () or phonons.garate13 () In contrast to on-site Anderson disorder, certain kinds of bond disorder cannot produce a TAI as they lead only to a positive correction to .song12 (); lv13 () So far, however, the TAI phase was not found in the Kane-Mele model on a honeycomb lattice, describing for example graphene or proposed TIs such as silicene, germanene and stanene.aufray10 (); kara12 (); davila14 (); zhu15 () First indications to this phase were already found, showing that the Kane-Mele model without a staggered sublattice potential hosts extended bulk states even for large disorder strengths.prodan11b ()

In this paper we show the existence of TAIs in the Kane-Mele model by means of tight-binding calculations. The interplay between the parameters characterizing intrinsic spin-orbit coupling (SOC) , extrinsic Rashba SOC , and a staggered sublattice potential turns out to be crucial for the visibility of TAIs, and we calculate the parameter ranges in which TAIs can be observed. We find analytically that to lowest order in , the parameters and are not renormalized with increasing disorder strength, in contrast to . However, a new effective hopping is generated due to the disorder, which is related but not identical to . Although is not a crucial ingredient for the existence of TAIs, it significantly alters the physics of topological insulators in various ways orth13 (); rod15 () and, as we will show below, strongly affects the TAI state.

Even though recently first signs of a TAI phase may have been found experimentally in evanescently coupled waveguides,stutzer15 () there has been no experimental evidence so far for the existence of the TAI phase in fermionic systems. The main difficulty is the requirement of a rather large and specific amount of disorder, which is difficult to control in the topological insulators currently investigated, where the 2D TI layer is buried inside a semiconductor structure. In contrast, producing and controlling disorder in 2D materials described by the Kane-Mele model could be much easier. Disorder in 2D materials with honeycomb structure can be produced by randomly placed adatoms weeks11 (); hua12 () or a judicious choice of substrate material.ando06 (); ishigami07 (); fratini08 (); varlet15 () Moreover, a sizeable staggered sublattice potential can be generated via a suitable substrate material.nevius15 () Other means of engineering disorder were proposed in periodically driven systems.titum15 (); yang15 () Finally, honeycomb structures with the SOC necessary to produce a topological phase have already been realized using ultracold atoms in optical lattices,jotzu14 () in which disorder can in principle be engineered.

## Results

### Setup

The basis of our calculations is the Kane-Mele model kane05b () given by the following Hamiltonian on a tight-binding honeycomb lattice

(1) |

which has been supplemented by an on-site Anderson disorder term with disorder strength and uniformly distributed random variables . The summations over the lattice sites and include all nearest neighbors and next-nearest neighbors, respectively. The operators , are creation and annihilation operators for the site of the lattice. The parameters , and are the nearest-neighbor hopping strength, intrinsic SOC, and Rashba SOC, respectively. If the next-nearest neighbor hopping from site to site corresponds to a right-turn on the honeycomb lattice, then , otherwise . Furthermore, is the vector of Pauli matrices for the spin degree of freedom, and is the unit vector between sites and . The Wannier states at the two basis atoms of the honeycomb lattice are separated in energy by twice the staggered sublattice potential , with for the sublattice and for the sublattice. The lattice constant is .

The band structure of this model depends strongly on the parameter set , , and . In the clean limit and for , the system will be a topological insulator for and a trivial insulator otherwise.kane05b () The tight-binding lattice and examples for the band structure in the clean limit are displayed in Fig. 1. For , the system will be a topological insulator if and a metal or semimetal otherwise. For finite and the situation is more complex and a topological transition appears for values within these two boundaries.

### Numerical solution

For , we find a TAI phase for parameters close to the topological transition at . Changing this ratio corresponds to changing the band mass in the case of the BHZ model. Figure 2 shows the conductance for different values of . We find that for the system is a topological insulator. For , i.e., in the clean case, the conduction and valence bands are separated by a red region with a quantized conductance of . Remarkably, with increasing disorder strength, the states in the conduction and valence bands localize, but the helical edge states that are responsible for the conductance of exist for an even larger energy window. The conductances and the vanishing error bars for the two distinct energy values , in the lower row of Fig. 2 show that the conductance quantization, and with it the topological nature of the system, persist for the vast majority of microscopic disorder configurations. Interestingly, for , the system is a trivial insulator at . The trivial gap closes however, and at the system develops a topologically non-trivial gap and edge states. This can be seen from the quantized conductance. Finally, for , the closing of the trivial gap and re-opening of the topological gap happens at a disorder strength which is strong enough to destabilize the emergent topological phase. Features of the conductance quantization can still be seen, but this behavior is not that robust anymore. As no averaging is done in the upper row of Fig. 2, and a new disorder configuration is taken for every data point, destabilization of the topological phase can be seen by red and white speckles in the figure.

We find that no TAI exists without staggered sublattice potential (). If both and are finite, the TAI phase is in general less pronounced, see Fig. 3. The plot on the right shows the closing of a trivial gap and emergence of a topological phase at .

Furthermore, we observe that the simultaneous presence of intrinsic and Rashba SOC (both and ) destroys the particle-hole symmetry in the spectrum. In the absence of Rashba SOC, the symmetry operator , which acts on the lattice operators as and for the sublattices and , leaves the (disorder-free) Hamiltonian invariant. can be viewed as particle-hole conjugation combined with spatial inversion, and the inversion is needed to leave the staggered sublattice potential term invariant.

### Lowest-order Born approximation

In the self-consistent Born approximation, the self-energy for a finite disorder strength is given by the following integral equation bruus04 (); groth09 ()

(2) |

where is the Fourier transform of in the clean limit.kane05b () The coefficient originates from the second moment of the uniform distribution function of the disorder amplitudes, and is the chemical potential. The integration is over the full first Brillouin zone. We use the lowest-order Born approximation, which means setting on the right-hand side of the equation.

After a low-energy expansion of , the integral can be evaluated analytically groth09 () for . This requires keeping the terms up to second order in wherever this is the leading -dependent order. The evaluation yields the renormalized staggered sublattice potential

(3) |

For a certain set of parameters, the logarithm can be negative and is reduced compared to . Moreover, we find that is not renormalized to order . Therefore, it is possible to obtain . The system thus makes a transition from a trivial insulator to a topological insulator with increasing .

For a more quantitative treatment, we evaluate the integral for the full Hamiltonian numerically. The self-energy is then written as a linear combination of several independent contributions

(4) |

with and . Here, denote the Pauli matrices for the sublattice index. This leads to the following equations for the renormalized quantities

(5) |

whereas and remain unrenormalized to lowest order in . Surprisingly, a new coupling is created by the disorder. This coupling has the matrix structure , which is similar but not identical to the one for Rashba SOC. Expressing this new term in the lattice coordinates of Eq. (Setup) reveals that it corresponds to a Rashba-type nearest-neighbor hopping term which is asymmetric and appears only for bonds that are parallel to the unit vector ,

(6) |

where stands for summations over strictly vertical bonds only. Furthermore, we find to lowest order in that for .

For , the upper and lower edge of the gap are at the energies . This is the case for both topological and trivial insulators. Extrapolation of these equations to finite leads to the conditions . The solid black lines in Fig. 2 are the two solutions to these equations and describe the closing and reopening of the gap qualitatively for small .

For finite and therefore finite , there is no analytical expression of the gap energy. In this case, we read off the positions of the gap edges from band structure calculations for several values of and . An interpolation leads to two functions for the upper and lower band edge in the clean system. Replacing the unperturbed by the renormalized parameters yields two equations

(7) |

The solutions of these equations are indicated by the solid black lines in Fig. 3. Hence, these results agree with the numerical data for small without any fitting parameter. Deviations appear for larger , when the lowest-order Born approximation is not applicable.

### Phase diagram

Figure 4 shows a phase diagram as a function of and based on the tight-binding simulations. The dark color marks the regions for which a critical disorder strength exists above which the system is a TAI (blue for , red for ). The TAI phase is located along the boundary separating trivial from topological insulators in the clean case. Towards larger , the TAI region becomes narrower and eventually vanishes above a critical . Figure 5 shows the critical disorder strength as a function of for a fixed value of .

In Figs. 2 and 3 rather large values of the parameters and were chosen to better visualize the effect. The TAI phenomenon scales down also to smaller values of the parameters, as the red region in the Fig. 4 indicates, but the TAI phase becomes less pronounced in the conductance plots and is harder to identify. Material parameters for stanene for example are , xu13 () and .liu11c () We suspect that disorder, e.g., originating from missing or dislocated atoms, can reach disorder strengths in the eV range.

### Alternative disorder models

Anderson disorder is a special model for disorder which is not necessarily representative for all TI materials. To better understand the effect of the disorder model, we briefly remark on the following disorder Hamiltonian

(8) |

In contrast to the Anderson disorder model, where a random potential is assigned to every lattice site, here the distribution function for is such that only a fraction of the sites are affected by disorder. Denoting the total number of sites by , we assume on sites, on sites, and on the remaining sites. The disorder amplitude is constant. Because , the normalization factor in Eq. (8) ensures that the mean squared disorder strength is equal to the Anderson disorder case for .

For general , the prefactor in Eq. (2) is thus replaced by . The lowest-order Born approximation for the disorder model (8) therefore predicts that a reduced disorder density can be exactly compensated by an increased amplitude . For large enough , this is indeed confirmed in the tight-binding simulations.

However, because a single impurity () cannot destroy the topological phase, it is clear that the TAI phase should eventually vanish for at arbitrary . Nevertheless, we find numerical evidence for the TAI phase at surprisingly low impurity densities. A TAI region remains pronounced for densities as low as .

## Discussion

In conclusion, we have shown that the topological Anderson insulator is a significantly more universal phenomenon than previously thought. Using a combination of an analytical approach and tight-binding simulations, we have established that the topological Anderson insulator appears in the Kane-Mele model that describes potential topological insulators such as silicene, germanene, and stanene and that can also be realized in optical lattices. We have observed a transition from a trivially insulating phase to a topological phase at a finite disorder strength and have mapped out the phase diagram as a function of the staggered sublattice potential () and the Rashba spin-orbit coupling (). The new Anderson insulator exists at the boundary between trivial and topological insulators for small and finite , but not at the boundary between a semimetal and a topological insulator for small and finite . Since the Kane-Mele model on a honeycomb lattice describes a wide class of candidate materials for topological insulators, we hope that our work will trigger experimental efforts to confirm the existence of the topological Anderson insulator.

## Methods

The numerical simulations were done with the tight-binding Hamiltonian (Setup) on a honeycomb lattice with rectangular shape of width and length using the Kwant code.groth14 () A smaller version of the sample is shown in Fig. 1. Both the upper and lower edge are taken to be of zigzag type. At the left and right edges two semi-infinite, metallic leads of width are attached. The leads are also modeled by a honeycomb lattice with only nearest-neighbor hopping and a finite on-site energy of to bring them into the metallic regime.

## Acknowledgments

CPO, TS, and CB acknowledge financial support by the Swiss SNF and the NCCR Quantum Science and Technology. TLS acknowledges support by National Research Fund, Luxembourg (ATTRACT 7556175).

## Author contributions

CPO and TS performed the numerical simulations with input from CB and TLS. All authors contributed to writing the manuscript.

## Additional information

Competing financial interests: The authors declare no competing financial interests.

## References

- (1) Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82, 3045 (2010).
- (2) Qi, X. & Zhang, S. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057 (2011).
- (3) Bernevig, B. A., Hughes, T. L. & Zhang, S.-C. Quantum spin Hall effect and topological phase transition in HgTe quantum wells. Science 314, 1757 (2006).
- (4) Kane, C. L. & Mele, E. J. Quantum spin Hall effect in graphene. Phys. Rev. Lett. 95, 226801 (2005).
- (5) Kane, C. L. & Mele, E. J. topological order and the quantum spin Hall effect. Phys. Rev. Lett. 95, 146802 (2005).
- (6) König, M. et al. Quantum spin Hall insulator state in HgTe quantum wells. Science 318, 766 (2007).
- (7) Knez, I., Du, R.-R. & Sullivan, G. Evidence for helical edge modes in inverted InAs/GaSb quantum wells. Phys. Rev. Lett. 107, 136603 (2011).
- (8) Li, T. et al. Observation of a helical Luttinger liquid in InAs/GaSb quantum spin Hall edges. Phys. Rev. Lett. 115, 136804 (2015).
- (9) Hsieh, D. et al. A topological Dirac insulator in a quantum spin Hall phase. Nature 452, 970 (2008).
- (10) Li, J., Chu, R.-L., Jain, J. K. & Shen, S.-Q. Topological Anderson insulator. Phys. Rev. Lett. 102, 136806 (2009).
- (11) Jiang, H., Wang, L., Sun, Q.-f. & Xie, X. C. Numerical study of the topological Anderson insulator in HgTe/CdTe quantum wells. Phys. Rev. B 80, 165316 (2009).
- (12) Anderson, P. W. Absence of diffusion in certain random lattices. Phys. Rev. 109, 1492 (1958).
- (13) Chen, L., Liu, Q., Lin, X., Zhang, X. & Jiang, X. Disorder dependence of helical edge states in HgTe/CdTe quantum wells. New J. Phys. 14, 043028 (2012).
- (14) Girschik, A., Libisch, F. & Rotter, S. Percolating states in the topological Anderson insulator. Phys. Rev. B 91, 214204 (2015).
- (15) Groth, C. W., Wimmer, M., Akhmerov, A. R., Tworzydło, J. & Beenakker, C. W. J. Theory of the topological Anderson insulator. Phys. Rev. Lett. 103, 196805 (2009).
- (16) Prodan, E. Three-dimensional phase diagram of disordered HgTe/CdTe quantum spin-Hall wells. Phys. Rev. B 83, 195119 (2011).
- (17) Xing, Y., Zhang, L. & Wang, J. Topological Anderson insulator phenomena. Phys. Rev. B 84, 035110 (2011).
- (18) Yamakage, A., Nomura, K. & Imura, K.-I., Kuramoto, Y. Disorder-Induced Multiple Transition Involving Z2 Topological Insulator. J. Phys. Soc. Jpn. 80, 053703 (2011).
- (19) Guo, H.-M., Rosenberg, G., Refael, G. & Franz, M. Topological Anderson insulator in three dimensions. Phys. Rev. Lett. 105, 216601 (2010).
- (20) Fu, B., Zheng, H., Li, Q., Shi, Q. & Yang, J. Topological phase transition driven by a spatially periodic potential. Phys. Rev. B 90, 214502 (2014).
- (21) Garate, I. Phonon-induced topological transitions and crossovers in Dirac materials. Phys. Rev. Lett. 110, 046402 (2013).
- (22) Song, J., Liu, H., Jiang, H., Sun, Q.-f. & Xie, X. C. Dependence of topological Anderson insulator on the type of disorder. Phys. Rev. B 85, 195125 (2012).
- (23) Lv, S.-H., Song, J. & Li, Y.-X. Topological Anderson insulator induced by inter-cell hopping disorder. J. Appl. Phys. 114, 183710 (2013).
- (24) Aufray, B. et al. Graphene-like silicon nanoribbons on Ag(110): A possible formation of silicene. Appl. Phys. Lett. 96, 183102 (2010).
- (25) Kara, A. et al. A review on silicene - new candidate for electronics. Surf. Sci. Rep. 67, 1 (2012).
- (26) Dávila, M. E., Xian, L., Cahangirov, S., Rubio, A. & Lay, G. L. Germanene: a novel two-dimensional germanium allotrope akin to graphene and silicene. New J. Phys. 16, 095002 (2014).
- (27) Zhu, F.-f. et al. Epitaxial growth of two-dimensional stanene. Nat. Mater. 14, 1020 (2015).
- (28) Prodan, E. Disordered topological insulators: a non-commutative geometry perspective. J. Phys. A: Math. Theor. 44, 113001 (2011).
- (29) Orth, C. P., Strübi, G. & Schmidt, T. L. Point contacts and localization in generic helical liquids. Phys. Rev. B 88, 165315 (2013).
- (30) Rod, A., Schmidt, T. L. & Rachel, S. Spin texture of generic helical edge states. Phys. Rev. B 91, 245112 (2015).
- (31) Stützer, S. et al. Experimental realization of a topological Anderson insulator. Paper presented at CLEO: QELS Fundamental Science 2015, San Jose (CA), United States, 10-15 May 2015, doi:10.1364/CLEO_QELS.2015.FTh3D.2.
- (32) Weeks, C., Hu, J., Alicea, J., Franz, M. & Wu, R. Engineering a robust quantum spin Hall state in graphene via adatom deposition. Phys. Rev. X 1, 021001 (2011).
- (33) Jiang, H., Qiao, Z., Liu, H., Shi, J. & Niu, Q. Stabilizing topological phases in graphene via random adsorption. Phys. Rev. Lett. 109, 116803 (2012).
- (34) Ando, T. Screening effect and impurity scattering in monolayer graphene. J. Phys. Soc. Jpn. 75, 074716 (2006).
- (35) Ishigami, M., Chen, J. H., Cullen, W. G., Fuhrer, M. S. & Williams, E. D. Atomic structure of graphene on SiO2. Nano Letters 7, 1643–1648 (2007).
- (36) Fratini, S. & Guinea, F. Substrate-limited electron dynamics in graphene. Phys. Rev. B 77, 195415 (2008).
- (37) Varlet, A. et al. Tunable Fermi surface topology and Lifshitz transition in bilayer graphene. ArXiv e-prints: 1508.02922 (2015).
- (38) Nevius, M. S. et al. Semiconducting Graphene from Highly Ordered Substrate Interactions. Phys. Rev. Lett. 115, 136802 (2015).
- (39) Titum, P., Lindner, N. H., Rechtsman, M. C. & Refael, G. Disorder-induced Floquet topological insulators. Phys. Rev. Lett. 114, 056801 (2015).
- (40) Yang, Z. et al. Topological acoustics. Phys. Rev. Lett. 114, 114301 (2015).
- (41) Jotzu, G. et al. Experimental realization of the topological Haldane model with ultracold fermions. Nature 515, 237–240 (2014).
- (42) Groth, C. W., Wimmer, M., Akhmerov, A. R. & Waintal, X. Kwant: a software package for quantum transport. New J. Phys. 16, 063065 (2014).
- (43) Bruus, H. & Flensberg, K. Many-Body Quantum Theory in Condensed Matter Physics (Oxford University Press, 2004).
- (44) Xu, Y. et al. Large-gap quantum spin Hall insulators in tin films. Phys. Rev. Lett. 111, 136804 (2013).
- (45) Liu, C.-C., Jiang, H. & Yao, Y. Low-energy effective Hamiltonian involving spin-orbit coupling in silicene and two-dimensional germanium and tin. Phys. Rev. B 84, 195430 (2011).