# Topological Landau Lattice

###### Abstract

The concept of topological fermions, including Weyl and Dirac fermions, stems from the quantum Hall state induced by a magnetic field, but the definitions and classifications of topological fermions are formulated without using magnetic field. It is unclear whether and how the topological information of topological fermions can be probed once their eigen spectrum is completely rebuilt by a strong magnetic field. In this work, we provide an answer via mapping Landau levels (bands) of topological fermions in dimensions to the spectrum of a -dimensional lattice model. The resultant “Landau lattice” may correspond to a topological insulator, and its topological property can be determined by real-space topological invariants. Accordingly, each zero-energy Landau level (band) inherits the topological stability from the corresponding topological boundary state of the Landau lattice. The theory is demonstrated in detail by transforming 2D Dirac fermions under magnetic fields to the Su-Schrieffer-Heeger models in class AIII, and 3D Weyl fermions to the Chern insulators in class A.

Introduction– The practice of defining phases in quantum matter by their topological properties traces its root back to the quantum Hall effects of electron gases in strong magnetic fields Klitzing et al. (1980); Thouless et al. (1982a). In the past decades, its generalization in the absence of magnetic field led to the discoveries of new phases of matter such as the topological insulators Hasan and Kane (2010a); Qi and Zhang (2011a); Shen (2012) and semimetals. In particular, classifying and modeling topological massless fermions with various symmetries Volovik (2003); Hořava (2005); Zhao and Wang (2013); Shiozaki and Sato (2014); Chiu and Schnyder (2014); Chiu et al. (2016); Zhao et al. (2016); Bradlyn et al. (2016) lie at the frontiers of condensed-matter physics. A number of quantum materials have been proposed and experimentally verified Wan et al. (2011); Xu et al. (2015a); Lv et al. (2015); Lu et al. (2015); Murakami (2007); Wang et al. (2012); Liu et al. (2014); Xu et al. (2015b); Burkov et al. (2011); Bian et al. (2016); Schoop et al. (2016); Neupane et al. (2016); Wu et al. (2016); Wang et al. (2016); Bradlyn et al. (2016); Chen et al. (2017); Yan et al. (2017), among them the topological Weyl Wan et al. (2011); Xu et al. (2015a); Lv et al. (2015); Lu et al. (2015), Dirac Murakami (2007); Wang et al. (2012); Liu et al. (2014); Xu et al. (2015b), and nodal-line Burkov et al. (2011); Bian et al. (2016); Schoop et al. (2016); Neupane et al. (2016); Wu et al. (2016) semimetals. When exposed to a magnetic field, the topological massless fermions form novel Landau levels (bands) that are essentially different from ordinary Fermi liquids, leading to a number of exotic properties such as the Berry phase corrected half-integer Hall conductance Zheng and Ando (2002); Gusynin and Sharapov (2005); Novoselov et al. (2005); Zhang et al. (2005); Xu et al. (2014); Yoshimi et al. (2015); Zhang et al. (2015a), chiral anomaly Nielsen and Ninomiya (1983), charge density wave Alicea and Balents (2009); Wang and Zhang (2013) and Tomonaga-Luttinger liquid Zhang and Nagaosa (2017) with high Landau degeneracy, and Weyl node annihilation Zhang et al. (2015b). Although one may naturally expect that these phenomena originate from the topological properties of the field-free fermions, a critical question here is whether and how a correspondence can be established between them, knowing that the eigen spectrum of the massless fermions can be completely re-structured in the presence of magnetic fields and the topological information hidden.

In this work, we answer this question by providing a general theory for defining and characterizing the topological properties of the Landau levels, whereby the Landau levels of the topological massless fermions are mapped to a semi-infinite lattice model, which we refer to as a Landau lattice. The Landau lattice can be extended to be infinite, and the field-hidden topological information of massless fermions can be extracted by studying topological properties of the extended Landau lattice. Because the infinite Landau lattice has no translational symmetry, we employ the topological invariants formulated in real space Kitaev (2006); Hastings and Loring (2010, 2011); Prodan and Schulz-Baldes (2016), which has been developed in the context of strong topological insulators and superconductors Hasan and Kane (2010b); Qi and Zhang (2011b); Chiu et al. (2016); Schnyder et al. (2008); Kitaev (2009); Ryu et al. (2010). Moreover, we discover that the Landau levels (bands) crossing zero energy are topologically protected by the topological configuration of the states with high Landau indices and are stable against perturbations. As the low and high Landau levels (bands) correspond to the boundary and bulk states in the Landau lattice, respectively, this is a natural consequence of the bulk-boundary correspondence Hatsugai (1993); Volovik (2003); Kitaev (2006); Zhao and Wang (2014); Prodan and Schulz-Baldes (2016), which dictates that the topologically protected boundary states are determined by a bulk topological invariant. The low-energy effective theory for the zero-energy low Landau levels (bands) is then further established by studying the topological invariants of the bulk states in the lattice model.

Systems in magnetic fields as tight-binding models– We start with a 2D system in a perpendicular uniform magnetic field . The gauge covariant momentum operators are and with , which may be linearly recombined to form the boson creator and annihilator of a harmonic oscillator,

(1) |

where and . The operators act on the Hilbert space spanned by the orthonormal particle-number basis with integers . The creator adds a particle and the annihilator destructs one, which are implemented by

(2) |

The eigenvalue of the particle number is then the Landau index .

We then construct the real-space component of the Hilbert space for a tight-binding model on a 1D semi-infinite lattice. The lattice sites are labelled as , each assigned with a quantum state , as illustrated in Fig. 1(a). As a 1D lattice, the forward and backward translators, denoted by and , respectively, act on the Hilbert space as

(3) | |||

(4) |

Another useful operator for the lattice is the position operator which reads off the position of a state,

(5) |

It is noteworthy that and commute for every state except the end state of the semi-infinite lattice , namely

(6) |

We then identify the real-space Hilbert space of the semi-infinite 1D lattice with that of the harmonic oscillator. In this identification, the position operator is just the Landau index operator , and the translators can be expressed as , . Inversely, with ,

(7) |

Before proceeding, we note that topological configurations shall be defined for the entirety of Landau levels, namely in the quotient space by the degeneracy of each Landau level. This should be distinguished from the quantum hall effects, where the topology is defined in the degenerate subspaces.

In general, a coarse-grained Hamiltonian for topological fermions in a condensed-matter system may be cast into the form , where are a set of matrices, Dirac matrices for instance, accounting for the internal degrees of freedom for topological fermions, and is the effective momentum in a -dimensional momentum space. Such a Hamiltonian may describe low-energy excitations arising from band crossings in a band structure. If a field is applied normal to the - plane, the corresponding gauge potentials satisfy , and the momenta and are gauged to be and , respectively. The recombined operators and are mapped to be forward and backward translators and in a 1D lattice, whose sites are labelled by . By such construction, a D model of topological fermions under a constant field is transformed to a D semi-infinite lattice with the D boundary perpendicular to , which we refer to as a Landau lattice of the topological fermions. More generally, for , if constant fields are applied, we obtain a D lattice with boundaries perpendicular to the independent dimensions via the transformation. For physical dimensions and , a magnetic field reduces the dimensionality to and , respectively.

We then extend the semi-infinite Landau lattice to be infinite by adding states , as illustrated in Fig. 1(a). On the infinite Landau lattice, the two translators are denoted by “removing hats”, and accordingly act as

(8) |

with . The extended operators are invertible and mutually independent, satisfying the relations

(9) |

We now have a D infinite lattice model—the extended Landau lattice—derived from D topological fermions under a uniform magnetic field. The infinite lattice model may be a topological insulator or (semi)metal that hosts topological gapless boundary modes when a D boundary is opened perpendicular to the direction. Inversely, these topological boundary modes correspond to the Landau levels (bands) of fermions that are of small Landau indices, and are located at (go across) zero energy. Since the boundary modes of the semi-infinite lattice are topologically protected by the bulk topological invariants, these corresponding Landau levels of small indices are topologically protected by nontrivial topological configurations encoded deep into the Landau levels of large indices. Similar to the robustness of the topological boundary modes, these corresponding topological Landau levels (bands) with small Landau indices are also stable against (symmetry-preserving) disorder and perturbations. It is noteworthy that the resultant infinite lattice model has no translational symmetry along the -direction due to the position-dependent mappings in Eq. (7). Thus, in contrast to the common practice of using topological invariants in terms of momentum-space Hamiltonian or Berry connection in the presence of translational symmetry, it is required here to employ topological invariants directly formulated in real space,

2D Dirac fermions and the SSH models– In the following, we proceed to illustrate and demonstrate the somewhat abstract general theory by elementary examples, which are interesting in their own right. Let us first consider 2D Dirac fermions with the Hamiltonian

(10) |

where and . Here are the Pauli matrices. Coupling Dirac fermions to the magnetic field is realized by substitutions and , which lead to with . The corresponding Landau lattice can be derived using Eq. (7) as , which inherits the chiral symmetry and is therefore in class AIII. The chiral symmetry charges a state with positive or negative chirality, i.e., . The obtained lattice model has an energy gap in the bulk, and one can determine whether it has topologically protected zero end-states by extending it to an infinite lattice. Note that the extension is not unique, and the inversion-symmetric one is given as

(11) |

Let us recall some basics of topological insulators in class AIII. The Hamiltonian can be written in a general form

(12) |

with the chiral symmetry , and the topological invariant for the classification can be given in real space as

(13) |

where is the number of sites Kitaev (2006); Not (), which when sent to infinity quantizes the topological invariant to an integer. A positive (negative) corresponds to (-) zero-modes with positive (negative) chirality concentrated at the end of the semi-infinite lattice. Note that the topological invariant in real space requires no translation invariance and thus serves well for our purpose.

For the case of Eq. (11) with , the topological invariant can then be directly computed by Eq. (13). As shown in Fig. 2, the topological invariant converges to with increasing . Alternatively, one may reduce this case to the simplest topologically nontrivial Su-Schrieffer-Heeger (SSH) model Su et al. (1979, 1980) with by simplifying the topological invariant to , or through a continuous deformation with varying from to , in course of which the bulk gap is always open and the topological invariant unchanged.

According to the bulk-boundary correspondence, there is a topologically-protected zero mode residing at the end of the semi-infinite lattice. Translating this back to the field-coupled Dirac fermions, this implies that there exists a Landau level composed of states with small Landau indices. We find the zero-energy state to be , which has positive chirality, in accord with the bulk topological invariant. The rest of the Landau levels can be interpreted as bulk states, each of which is a superposition of states from nearest neighbors on the lattice. In other words, they correspond to the bulk states of the SSH model with the pairings of the bulk states, and , where , and the Hamiltonian is block-diagonalized in the subspace of each pairing, as illustrated in Fig. 1. Further diagonalizing the Hamiltonian in each block solves the spectrum of the eigenstates with .

Now we elaborate a significant difference in phase transition between the Landau lattice model in class AIII and the translation-invariant SSH model. Let us recall the generic SSH model with a symmetry-preserving on-site term, , which corresponds to . For , the model is in a trivial phase without end states since the on-site term is dominant. On the other hand, an arbitrarily large cannot demote the Landau lattice with to a trivial phase. As shown by numerical calculations in Fig. 2, always approaches for any finite when the lattice length is sufficiently large. More precisely, as seen from Fig. 2(a), there exists a critical site number , exceeding which the topological invariant jumps from to . Numerical fitting of the results in Fig. 2(b) shows that is approximately a quadratic function of , whose explicit form reads . We emphasize that the criticality here corresponds to no phase transition in the ordinary sense, since the notion of a phase is defined under the thermodynamic limit . The fact that the Landau lattice is always topologically nontrivial can be intuitively understood as follows: The nearest hopping term responsible for the topology will always dominate with increasing magnitude when going deep enough into the bulk. What actually happens with increasing is that the concentration of the end state is moved towards the bulk, and its spatial distribution is smoothed [see the Supplemental Material (SM) Sup () for details]. But it is inescapable to be bounded to the end for a finite on the semi-infinite lattice, since always approaches as . Accordingly, under an arbitrarily large perturbation , there is always a zero-energy Landau level that is a superposition dominated by states with small Landau indices.

Recalling that 1D topological insulators in class AIII have an integer classification, we can construct the corresponding Landau levels of topological fermions with an arbitrary number of topologically protected zero-modes mainly composed of small indices, which can be found in SM Sup ().

Weyl fermions and Chern insulators– We now consider a 3D example in class A without any symmetry, namely the Weyl-fermion Hamiltonian, , whose infinite Landau lattice is derived as

(14) |

where the -direction is expressed in real space without translational symmetry, while the -direction remains in momentum space with translational symmetry. This is a 2D lattice with an energy gap, since is invertible, and anticommutes with . The model can be continuously deformed to be without closing the gap, so that translational symmetry along the -direction is acquired, entirely parallel to what we did for Eq. (11). Now the model may be a Chern insulator Haldane (1988); Thouless et al. (1982b), and the Chern number can be calculated using translational symmetries, although a formula of the Chern number in real space is also available Kitaev (2006). In momentum space, the Hamiltonian is , and the Chern number is just the unit winding number of the vector field Volovik (2003), after normalized to be , as a mapping from to . It can be computed explicitly as . The positive (negative) unit Chern number implies a right (left)-moving edge band across zero energy for any 1D edge of a semi-infinite system. In this particular case where all sites left to are cut off, making an edge for the -dimension, the unit Chern number corresponds to a chiral band parametrized by , which are gapless edge modes for the semi-infinite lattice model, or equivalently a chiral Landau band dominated by states with small indices (see the SM Sup () for a graphic illustration). The existence of the chiral Landau band is topologically protected against weak disorder and perturbations. Chern insulators of higher Chern numbers correspond to Landau bands of multiply-charged Weyl fermions, which is elaborated in the SM similar to the previous case of symmetry class AIII Sup ().

In a two-band minimal model of Weyl semimetals, for instance , there are two Weyl points each with left and right chirality, conforming with the Nielsen-Ninomiya no-go theorem Nielsen and Ninomiya (1981); Zhao and Wang (2016). By applying a uniform magnetic field along the -direction, according to the discussions above, the right (left) Weyl point leads to a right (left) massless mover, which is concentrated at the smallest Landau index. In the virtual Landau band structure as graphically illustrated in the SM Sup (), the left and right moving modes should be connected smoothly, leading to a band parametrized by as a reminiscence of the Tomonaga-Luttinger liquid when all negative states are filled by electrons. Again such topological properties are insensitive to weak disorder and perturbations.

Summary - In summary, we bridge the Landau levels (bands) of -dimensional topological fermions to a D Landau lattice with . It is shown that the topological boundary modes of the Landau lattice, resulting from the bulk topological invariant, correspond to topological low Landau levels (bands) crossing zero energy. Thus, these zero modes among Landau levels (bands) are protected by topological configurations deep into the gapped states with large Landau indices, which are therefore robust against (symmetry-preserving) weak disorder and perturbations.

Acknowledgements - This work was supported by the National Key R & D Program (Grant No. 2016YFA0301700), National Natural Science Foun- dation of China (Grant No. 11574127), and Guangdong Innovative and Entrepreneurial Research Team Program (Grant No. 2016ZT06D348).

## References

- Klitzing et al. (1980) K. v. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett. 45, 494 (1980), URL http://link.aps.org/doi/10.1103/PhysRevLett.45.494.
- Thouless et al. (1982a) D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Phys. Rev. Lett. 49, 405 (1982a), URL http://link.aps.org/doi/10.1103/PhysRevLett.49.405.
- Hasan and Kane (2010a) M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010a), URL http://link.aps.org/doi/10.1103/RevModPhys.82.3045.
- Qi and Zhang (2011a) X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011a), URL http://link.aps.org/doi/10.1103/RevModPhys.83.1057.
- Shen (2012) S.-Q. Shen, Topological Insulators (Springer-Verlag, Berlin Heidelberg, 2012).
- Volovik (2003) G. E. Volovik, Universe in a helium droplet (Oxford University Press, Oxford UK, 2003), ISBN 0521670535.
- Hořava (2005) P. Hořava, Phys. Rev. Lett. 95, 016405 (2005), URL https://link.aps.org/doi/10.1103/PhysRevLett.95.016405.
- Zhao and Wang (2013) Y. X. Zhao and Z. D. Wang, Phys. Rev. Lett. 110, 240404 (2013), URL http://link.aps.org/doi/10.1103/PhysRevLett.110.240404.
- Shiozaki and Sato (2014) K. Shiozaki and M. Sato, Phys. Rev. B 90, 165114 (2014), URL http://link.aps.org/doi/10.1103/PhysRevB.90.165114.
- Chiu and Schnyder (2014) C.-K. Chiu and A. P. Schnyder, Phys. Rev. B 90, 205136 (2014), URL http://link.aps.org/doi/10.1103/PhysRevB.90.205136.
- Chiu et al. (2016) C.-K. Chiu, J. C. Y. Teo, A. P. Schnyder, and S. Ryu, Rev. Mod. Phys. 88, 035005 (2016), URL http://link.aps.org/doi/10.1103/RevModPhys.88.035005.
- Zhao et al. (2016) Y. X. Zhao, A. P. Schnyder, and Z. D. Wang, Phys. Rev. Lett. 116, 156402 (2016), URL http://link.aps.org/doi/10.1103/PhysRevLett.116.156402.
- Bradlyn et al. (2016) B. Bradlyn, J. Cano, Z. Wang, M. G. Vergniory, C. Felser, R. J. Cava, and B. A. Bernevig, Science 353 (2016), URL http://science.sciencemag.org/content/early/2016/07/20/science.aaf5037.
- Wan et al. (2011) X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, Phys. Rev. B 83, 205101 (2011), URL http://link.aps.org/doi/10.1103/PhysRevB.83.205101.
- Xu et al. (2015a) S. Y. Xu, I. Belopolski, N. Alidoust, M. Neupane, G. Bian, C. L. Zhang, R. Sankar, G. Q. Chang, Z. J. Yuan, C. C. Lee, et al., Science 349, 613 (2015a), URL http://www.sciencemag.org/content/349/6248/613.abstract.
- Lv et al. (2015) B. Q. Lv, H. M. Weng, B. B. Fu, X. P. Wang, H. Miao, J. Ma, P. Richard, X. C. Huang, L. X. Zhao, G. F. Chen, et al., Phys. Rev. X 5, 031013 (2015), URL http://link.aps.org/doi/10.1103/PhysRevX.5.031013.
- Lu et al. (2015) L. Lu, Z. Wang, D. Ye, L. Ran, L. Fu, J. D. Joannopoulos, and M. Soljačić, Science 349, 622 (2015), ISSN 0036-8075.
- Murakami (2007) S. Murakami, New. J. Phys. 9, 356 (2007).
- Wang et al. (2012) Z. Wang, Y. Sun, X. Q. Chen, C. Franchini, G. Xu, H. Weng, X. Dai, and Z. Fang, Phys. Rev. B 85, 195320 (2012), URL http://link.aps.org/doi/10.1103/PhysRevB.85.195320.
- Liu et al. (2014) Z. K. Liu, B. Zhou, Y. Zhang, Z. J. Wang, H. M. Weng, D. Prabhakaran, S. K. Mo, Z. X. Shen, Z. Fang, X. Dai, et al., Science 343, 864 (2014), URL http://www.sciencemag.org/content/343/6173/864.abstract.
- Xu et al. (2015b) S. Y. Xu, C. Liu, S. K. Kushwaha, R. Sankar, J. W. Krizan, I. Belopolski, M. Neupane, G. Bian, N. Alidoust, T. R. Chang, et al., Science 347, 294 (2015b), URL http://www.sciencemag.org/content/347/6219/294.abstract.
- Burkov et al. (2011) A. A. Burkov, M. D. Hook, and L. Balents, Phys. Rev. B 84, 235126 (2011), URL https://link.aps.org/doi/10.1103/PhysRevB.84.235126.
- Bian et al. (2016) G. Bian, T.-R. Chang, R. Sankar, S.-Y. Xu, H. Zheng, T. Neupert, C.-K. Chiu, S.-M. Huang, G. Chang, I. Belopolski, et al., 7, 10556 (2016), article, URL http://dx.doi.org/10.1038/ncomms10556.
- Schoop et al. (2016) L. M. Schoop, M. N. Ali, C. Straßer, A. Topp, A. Varykhalov, D. Marchenko, V. Duppel, S. S. P. Parkin, B. V. Lotsch, and C. R. Ast, Nature Commun. 7, 11696 (2016), URL http://www.nature.com/ncomms/2016/160531/ncomms11696/full/ncomms11696.html.
- Neupane et al. (2016) M. Neupane, I. Belopolski, M. M. Hosen, D. S. Sanchez, R. Sankar, M. Szlawska, S.-Y. Xu, K. Dimitri, N. Dhakal, P. Maldonado, et al., Phys. Rev. B 93, 201104 (2016), URL https://journals.aps.org/prb/abstract/10.1103/PhysRevB.93.201104.
- Wu et al. (2016) Y. Wu, L.-L. Wang, E. Mun, D. D. Johnson, D. Mou, L. Huang, Y. Lee, S. L. Bud’ko, P. C. Canfield, and A. Kaminski, Nature Phys. 12, 667 (2016), URL http://www.nature.com/nphys/journal/vaop/ncurrent/full/nphys3712.html?WT.feed_name=subjects_electronic-properties-and-materials.
- Wang et al. (2016) Z. Wang, A. Alexandradinata, R. J. Cava, and B. A. Bernevig, Nature 532, 189 (2016), ISSN 0028-0836, article, URL http://dx.doi.org/10.1038/nature17410.
- Chen et al. (2017) W. Chen, H.-Z. Lu, and J.-M. Hou, arXiv:1703.10886 (2017), URL http://arxiv.org/abs/1703.10886.
- Yan et al. (2017) Z. Yan, R. Bi, H. Shen, L. Lu, S.-C. Zhang, and Z. Wang, arXiv:1704.00655 (2017), URL http://arxiv.org/abs/1704.00655.
- Zheng and Ando (2002) Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002), URL https://link.aps.org/doi/10.1103/PhysRevB.65.245420.
- Gusynin and Sharapov (2005) V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 146801 (2005), URL https://link.aps.org/doi/10.1103/PhysRevLett.95.146801.
- Novoselov et al. (2005) K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005), ISSN 0028-0836, URL http://dx.doi.org/10.1038/nature04233.
- Zhang et al. (2005) Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 438, 201 (2005), URL http://www.nature.com/nature/journal/v438/n7065/abs/nature04235.html.
- Xu et al. (2014) Y. Xu, I. Miotkowski, C. Liu, J. Tian, H. Nam, N. Alidoust, J. Hu, C.-K. Shih, M. Z. Hasan, and Y. P. Chen, Nature Phys. 10, 956 (2014), URL http://www.nature.com/nphys/journal/v10/n12/full/nphys3140.html.
- Yoshimi et al. (2015) R. Yoshimi, K. Yasuda, A. Tsukazaki, K. S. Takahashi, N. Nagaosa, M. Kawasaki, and Y. Tokura, Nature Commun. 6 (2015), URL http://www.nature.com/articles/ncomms9530.
- Zhang et al. (2015a) S. B. Zhang, H. Z. Lu, and S. Q. Shen, Scientific reports 5, 13277 (2015a), URL http://www.nature.com/articles/srep13277.
- Nielsen and Ninomiya (1983) H. B. Nielsen and M. Ninomiya, Physics Letters B 130, 389 (1983), ISSN 0370-2693, URL http://www.sciencedirect.com/science/article/pii/0370269383915290.
- Alicea and Balents (2009) J. Alicea and L. Balents, Phys. Rev. B 79, 241101 (2009), URL https://link.aps.org/doi/10.1103/PhysRevB.79.241101.
- Wang and Zhang (2013) Z. Wang and S.-C. Zhang, Phys. Rev. B 87, 161107 (2013), URL https://link.aps.org/doi/10.1103/PhysRevB.87.161107.
- Zhang and Nagaosa (2017) X.-X. Zhang and N. Nagaosa, Phys. Rev. B 95, 205143 (2017), URL https://link.aps.org/doi/10.1103/PhysRevB.95.205143.
- Zhang et al. (2015b) C. Zhang, Z. Lin, C. Guo, S.-Y. Xu, C.-C. Lee, H. Lu, S.-M. Huang, G. Chang, C.-H. Hsu, H. Lin, et al., arXiv:1507.06301 (2015b), URL http://arxiv.org/abs/1507.06301.
- Kitaev (2006) A. Kitaev, Annals of Physics 321, 2 (2006), ISSN 0003-4916, january Special Issue, URL http://www.sciencedirect.com/science/article/pii/S0003491605002381.
- Hastings and Loring (2010) M. B. Hastings and T. A. Loring, Journal of Mathematical Physics 51, 015214 (2010), URL http://dx.doi.org/10.1063/1.3274817.
- Hastings and Loring (2011) M. B. Hastings and T. A. Loring, Annals of Physics 326, 1699 (2011), ISSN 0003-4916, july 2011 Special Issue, URL http://www.sciencedirect.com/science/article/pii/S0003491610002277.
- Prodan and Schulz-Baldes (2016) E. Prodan and H. Schulz-Baldes, Bulk and Boundary Invariants for Complex Topological Insulators: From K-Theory to Physics (Springer International Publishing, Switzerland, 2016).
- Hasan and Kane (2010b) M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010b), URL http://link.aps.org/doi/10.1103/RevModPhys.82.3045.
- Qi and Zhang (2011b) X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011b), URL http://link.aps.org/doi/10.1103/RevModPhys.83.1057.
- Schnyder et al. (2008) A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, Phys. Rev. B 78, 195125 (2008), URL http://link.aps.org/doi/10.1103/PhysRevB.78.195125.
- Kitaev (2009) A. Kitaev, AIP Conference Proceedings 1134, 22 (2009), URL http://aip.scitation.org/doi/abs/10.1063/1.3149495.
- Ryu et al. (2010) S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. Ludwig, New Journal of Physics 12, 065010 (2010), URL http://stacks.iop.org/1367-2630/12/i=6/a=065010.
- Hatsugai (1993) Y. Hatsugai, Phys. Rev. Lett. 71, 3697 (1993), URL https://link.aps.org/doi/10.1103/PhysRevLett.71.3697.
- Zhao and Wang (2014) Y. X. Zhao and Z. D. Wang, Phys. Rev. B 89, 075111 (2014), URL http://link.aps.org/doi/10.1103/PhysRevB.89.075111.
- (53) In numerical computing, edges with width of the order of tunnelling distance are cut off from the matrices in the trace to avoid boundary effects.
- Su et al. (1979) W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. Lett. 42, 1698 (1979), URL https://link.aps.org/doi/10.1103/PhysRevLett.42.1698.
- Su et al. (1980) W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. B 22, 2099 (1980), URL https://link.aps.org/doi/10.1103/PhysRevB.22.2099.
- (56) See Supplemental Materials for….
- Haldane (1988) F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988), URL http://link.aps.org/doi/10.1103/PhysRevLett.61.2015.
- Thouless et al. (1982b) D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Phys. Rev. Lett. 49, 405 (1982b), URL http://link.aps.org/doi/10.1103/PhysRevLett.49.405.
- Nielsen and Ninomiya (1981) H. Nielsen and M. Ninomiya, Nucl. Phys. B 185, 20 (1981), URL http://www.sciencedirect.com/science/article/pii/0550321381903618.
- Zhao and Wang (2016) Y. X. Zhao and Z. D. Wang, Phys. Rev. Lett. 116, 016401 (2016), URL http://link.aps.org/doi/10.1103/PhysRevLett.116.016401.
- Volovik (2011) G. E. Volovik, arXiv:1111.4627 (2011), URL https://arxiv.org/abs/1111.4627.
- Lu and Shen (2017) H.-Z. Lu and S.-Q. Shen, Front. Phys. 12, 127201 (2017), URL http://link.springer.com/article/10.1007%2Fs11467-016-0609-y.

Supplemental Material

## Appendix S1 2D Dirac fermions with muptile topological charges

Let us consider 2D massless Dirac fermions with quadratic band touching Volovik (2011),

(S15) |

The model has mirror reflection symmetries through the and -axis, corresponding to the operators and , respectively, with () being the inversion of the -coordinate (-coordinate), and therefore has inversion symmetry as the combination of and , . The model may appear in solid-state systems that satisfy these crystal symmetries. In a -direction perpendicular magnetic field, the Hamiltonian becomes

(S16) |

with . The corresponding infinite lattice model is

(S17) |

which gives . The model has an energy gap since is invertible with the inverse,

(S18) |

and can be continuously deformed to be via with , during which the gap is opened and the invertibility of is preserved. Now we calculate the topological invariant of in momentum space by using the acquired translation invariance, which turns out to be . Thus there exist two zero-energy end-states with positive chirality for the semi-infinite lattice with sites on the left of the 0th site being cut off. Accordingly, there are two-fold degenerate zero-energy Landau levels with positive chirality and small Landau indices.

This model can be generalized to be for any positive integer , which corresponds to a topological invariant , implying that the Landau levels with positive chirality at zero energy are -fold degenerate besides the Landau degeneracy. To construct models for , we just need to exchange and , which leads to . Correspondingly for this model there exist zero-energy end-states or equivalently zero-energy Landau levels with small indices, which have negative chirality.

## Appendix S2 Spatial Distributions of Zero-Energy Landau Levels corresponding to generic SSH models

The zero-energy wave function of , or the zero-energy end-state of the equivalent semi-infinite Landau lattice, , can be exactly solved as

(S19) |

The zero-energy state has been normalized for any , and has positive chirality for the chiral symmetry , which are consistent with the fact that the topological invariant is equal to for any . The probability at position is given by , which is plotted with different ’s in Fig. S3. It is observed that the concentration of the wave pocket is moved towards the bulk and its shape is smoothened, when is increased, although it is always bounded to the end for any .

## Appendix S3 Topological Landau bands of Weyl fermions and Semimetals

As shown in Fig. S4(a), in a magnetic field has a right-handed chiral band which crosses zero energy and has the smallest Landau index, or equivalently the corresponding Landau lattice is a semi-infinite Chern insulator that has a right-handed chiral edge band. The Landau-band structure for is just the mirror image of Fig. S4(a) through the -axis.

For the minimal model of Weyl semimetals Lu and Shen (2017), , the Landau-band structure is illustrated in Fig. S4(b). The low-energy expansion of consists of and at and , respectively, and thus the lowest Landau band equivalent to the edge states of the corresponding semi-infinite Landau lattice is obtained by connecting a left-handed chiral band crossing and a righ-handed one crossing . The other bands in Fig. S4(b) correspond to the gapped bulk bands of the semi-infinite Landau lattice.

The construction of a Weyl point with a Chern number of can be achieved by simply adding to , i.e., . The corresponding semi-infinite Landau lattice is a Chern insulator with a Chern number of , and has () flavors of right (left)-handed chiral bands on the edge when (), or equivelently () flavors of right (left)-handed chiral Landau bands with small indices crossing zero energy.