# Flat bands and perfect metal in trilayer moiré graphene

###### Abstract

We investigate the electronic structure of a twisted multilayer graphene system forming a moiré pattern. We consider small twist angles separating the graphene sheets and develop a low-energy theory to describe the coupling of Dirac Bloch states close to the K point in each individual plane. Extending beyond the bilayer case, we show that, when the ratio of the consecutive twist angles is rational, a periodicity emerges in quasimomentum space with moiré Bloch bands even when the system does not exhibit a crystalline lattice structure in real space. For a trilayer geometry, we find flat bands in the spectrum at certain rotation angles. Performing a symmetry analysis of the band model for the trilayer, we prove that the system is a perfect metal in the sense that it is gapless at all energies. This striking result originates from the three Dirac cones which can only gap in pairs and produce bands with an infinite connectivity. The full gapless property is protected by an emergent particle-hole symmetry valid at sufficiently small angles.

Two parallel layers of graphene twisted by a small angle exhibit a moiré pattern Pong and Durkan (2005), with a lattice periodicity much larger than the graphene unit cell, and non-trivial electronic properties Bistritzer and MacDonald (2010), such as band-flattening at certain magic angles Bistritzer and MacDonald (2011); Trambly de LaissardiÃ¨re et al. (2010); San-Jose et al. (2012); Lopes dos Santos et al. (2012); Fang and Kaxiras (2016); Tarnopolsky et al. (2018). Flat bands present a reduced kinetic energy, thereby artificially boosting electron correlations Kim et al. (2017). The recent discovery of correlated insulating phases at half-filling and possibly unconventional superconductivity Cao et al. (2018a, b); Carr et al. (2018); Yankowitz et al. (2018) have unveiled bilayer moiré graphene as a tunable device for exploring novel correlated states, at zero or finite magnetic field, spurring intense theoretical work in this direction Yuan and Fu (2018); Wu et al. (2018); Padhi et al. (2018); Po et al. (2018a); Isobe et al. (2018); Dodaro et al. (2018); Baskaran (2018); Huang et al. (2018); You and Vishwanath (2018); Wu et al. (2018); Roy and Juricic (2018); Xu and Balents (2018); Zhang et al. (2018); Koshino et al. (2018); Kang and Vafek (2018); Ray and Das (2018); Guinea and Walet (2018); Wu et al. (2018); Lian et al. (2018); Kang and Vafek (2018); Lian et al. (2018). On the other hand, moiré bands were also investigated for their topological properties San-Jose and Prada (2013); de Gail et al. (2011); Po et al. (2018b); Song et al. (2018); Ahn et al. (2018); Zou et al. (2018); Liu et al. (2018) and topological phase transitions were identified close to the magic angles Song et al. (2018); Hejazi et al. (2018).

In view of the great wealth of correlation and topological phenomena occuring with moiré bilayer graphene, it is desirable to extend studies to multilayer, and specifically trilayer geometries, in which moiré patterns also appear for small rotation angles. Flat bands in bilayer result from an interplay between the K (or K’) Dirac points in each layer and the situation with three or more Dirac points has yet to be explored. At low energy and close to half-filling, the band structure is formed only from the electron states of the Dirac cones in each layer Lopes dos Santos et al. (2007). A moiré band theory describing such as state, that does not require a crystalline lattice, is built in Ref. Bistritzer and MacDonald (2011). In this paper, we extend this theory to multilayer graphene and discuss in depth the symmetry and topology for three layers. We find magic angles of vanishing Dirac velocities; they are not related to a complete flattening of the spectrum, rather by a flattening along certain symmetry lines. We also characterize the different moiré bands by the irreducible representations they generate at the high-symmetry points and lines. Based on the compatibility between these representations Bradlyn et al. (2017); Vergniory et al. (2017), we are able to prove the remarkable result that all bands are connected such that no subset of bands can be energetically isolated from the others. The most obvious consequence is that the system remains metallic at arbitrary energy. This property requires particle-hole symmetry - neglecting the band curvature in the vicinity of the original K points.

Bloch band coupling. We detail the derivation of the band structure of the moiré pattern in twisted multilayer graphene. When the twisting angle is small, a moiré pattern is formed by the interference of lattices between the different layers. Restricting the analysis to Dirac fields near the K points of each layer Lopes dos Santos et al. (2007); Bistritzer and MacDonald (2011); Fang and Kaxiras (2016) (see also Ref. Amorim and Castro (2018)) and assuming a local short-range tunnel amplitude between atoms in consecutive planes, one derives the following Hamiltonian

(1) |

where are hopping energies between the neighboring layers and . The first term in Eq. (1) represents the Dirac cones in each layer and is a small momentum deviation from the point for the layer . We have introduced the matrices

(2) |

associated with the three symmetric momentum directions , , . the rotation around with the twist angle separating the layers and and with angle . For small twist angles, is perpendicular to so that all the are parallel for fixed (see below).

The magnitudes of the vectors depend on the twist angles. The second term in Eq. (1) couples momenta and in layers and which generates a lattice in momentum space for each pair of consecutive layers. To make the calculation tractable and maintain an emergent periodicity in momentum space regardless of whether the multilayer system itself is crystalline, one has to assume that the ratios of twist angles are rational numbers. This periodicity results in moiré bands which we compute numerically and classify according to their irreducible representations at the high-symmetry points and lines. On the contrary, if the twist angles were incommensurate, or the vectors are not parallel for fixed , then successive applications of the Hamiltonian (1) would reach arbitrary momentum and the moiré periodicity would be absent.

The Hamiltonian (1) reduces to the model of Ref. Bistritzer and MacDonald (2011) for bilayer and we henceforth focus on the trilayer geometry. We take the rotation angle as a reference and introduce the moiré magnitude . We rescale all momenta by and the Hamiltonian as . Fixing the direction of along , we use the complex notation

(3) |

and for all , where and are coprime integers. The Hamiltonian (1) can then be written as

(4) |

where we assume a uniform tunnel amplitude and introduce the dimensionless coupling between Dirac cones. The vectors form a -space lattice, see Fig. 1b, where each site is associated to a specific layer.

Equal twist angles. We first consider the most symmetric case of evenly rotated planes where . A representative set of moiré spectra obtained from Eq. (4) with different values of the coupling is displayed in Fig. 2a-d.

The moiré bands exhibit a rich structure. The first remarkable feature is that all bands are connected: it is impossible to isolate a set of bands which are detached from the rest. We provide below a formal proof for this statement based on irreducible representations at symmetric points and lines, and afterwards extend it to arbitrary and .

Three Dirac cones are attached at zero energy to the points , and (see Fig. 1b) as is varied. We display the corresponding Dirac velocities of the cones at and (with velocity linked by symmetry to that of ) in Fig. 2e and find a set of magic angles - in analogy with the bilayer case - where one of these velocities vanishes. The difference with the bilayer case is that these magic angles are not associated with a flattening of the whole spectrum which would be at odds with the fully connected band structure. However, we do see a flattening of part of the spectrum close to magic angles: on the line (for the second levels) when , close to when , for the first two magic angles. This flattening along one-dimensional directions in space open interesting perspectives for the realization of exotic correlated many-body physics.

Symmetries The moiré reciprocal lattice vectors , generate the whole lattice in Fig. 1b. The different layers (colors) are coupled by the vectors. Bloch periodicity takes the form

(5) |

The spectrum is thus invariant upon shifting the origin of by a combination of and .

The moiré lattice also transforms into itself by the action of a rotation around . In the Hamiltonian language, the corresponding operator is given by . We note that if is in a given layer so is , such that the three layers are not mixed by the rotation. The symmetry operates a reflexion accross the x-axis going through the point (red in Fig. 1b). As such, sends lattice sites of layer to and viceversa but keeps layer invariant. The corresponding symmetry operator is . The two symmetry operators induce the Hamiltonian transformation

(6) |

which leave the spectrum invariant. The antiunitary symmetry acts locally on the moiré lattice. It takes complex conjugation and reverses the pseudo-spin direction. It is represented by the operator , which squares to and commutes with the Hamiltonian . It also commutes with the spatial symmetries and .

The moiré model also possesses a unitary particle-hole (p-h) symmetry. The original Dirac Hamiltonian of the single layer graphene sheet has a unitary particle-hole symmetry due to the absence of terms in the Hamiltonian. Since our model is based on this low-energy expansion, it retains a similar symmetry with the operator

(7) |

where is for belonging to the lower and top layers and for the middle layer. With this, we have . Importantly, one checks that commutes with all other symmetry operators, , and , and satisfies

(8) |

This is to be contrasted with the p-h operator identified Song et al. (2018) for moiré bilayer which has: (i) , and (ii) anticommutes with , i.e. , instead of the commutation found for .

1 | 1 | 2 | 1 | 1 | 1 | 2 | 1 | 1 | ||||

1 | 1 | -1 | 1 | -1 | 1 | -1 | 1 | -1 | ||||

1 | -1 | 0 | 1 | -1 |

Based on the generators discussed so far, the symmetry group of the moiré lattice is the magnetic space group called (#177.151 in the BNS notation Gallego et al. (2012)). Although the same group describes moiré bilayer, the physics is different here. It indeed originates from Dirac cones, instead of two, and the extra particle-hole symmetry is essentially different. The high-symmetry points and their little co-groups are (), () and (). The symmetries on the high-symmetry lines are (, ) and (). The classification of the different irreducible representations at the symmetric points and lines are given in Table 1. At , , and on the line , each energy or band in Fig. 2 is characterized by a certain representation determined from the character, i.e. from the eigenvalues of the operators and restricted to this (possibly degenerate) energy.

First proof of all-connected bands. We now prove by contradiction that all bands are connected such that there is no gap in the spectrum at any energy. We assume a subspace of isolated bands between the energies and . By p-h symmetry, a symmetric set of bands exists in the energy window and, consequently, the bands between and must be disconnected from all other bands. We focus on these bands and investigate their transformation property under . remains a symmetry along the line connecting and such that their total character must coincide at each end, or .

We call the multiplicity of the representation in the set of bands, and the multiplicity of the representation in the set of bands at . Hence we have the equations: , , and which leads to

(9) |

For , we have three Dirac cones around , and and a gapped spectrum at . At , the zero energy subspace is doubly degenerate and the character of is simply given by corresponding to the irreducible representation as indicated in Table 1. We thus have , when restricted to zero energy and . Increasing away from zero, the symmetry is maintained and commutes with and such that any band associated to a given representation collapsing at (or departing from) zero energy, at either or , must move with its energy symmetric p-h partner associated to the same representation. As a result, the multiplicities and can only change by units of two in such processes. The same argument extends to non-zero energies where each band with energy and representation (or ) has a p-h partner with energy and the same representation. Since and are identified at and such that the interval is symmetric, we finally obtain by continuity with that the multiplicity must be odd while and are both even integers. It contradicts Eq. (9), thus completing the proof.

General proof for unequal twist angles. Our discussion was so far restricted to the symmetric configuration of equal and . The and symmetries are broken when and are different whereas and are maintained. However, a remnant of p-h symmetry still exists in the form of a mirror symmetry with respect to the plane orthogonal to the layer and crossing and , leaving each layer invariant. It corresponds to the operator with , acting as

(10) |

which associates pairs of mirror-symmetric momenta with opposite energies. Extending the arguments of Ref. Ahn et al. (2018), we show that a set of isolated bands with symmetry cannot accommodate an odd winding number , corresponding for example to an odd number of Dirac cones. The derivation is explicit in Ref. Song et al. (2018); Ahn et al. (2018) for two bands and a vanishing total Berry phase (Wilson loop) with , where the Euler class is an integer topological invariant. Adding more bands, windings around singularities can change sign but keep a definite parity while the parity of defines a topological invariant, the Stiefel-Whitney class . Then, the relation simply enforces that the winding number must be an even integer. More intuitively, we note that Dirac points are monopoles attaching Dirac strings. They can annihilate in pairs when of opposite signs or form a topological isolated band by combining pairs of same sign Song et al. (2018); Po et al. (2018b), but in all cases they need to pair to form an isolated set of bands.

We now show by contradiction that all bands are connected by gapless points in our trilayer moiré model for arbitrary and . As already discussed above, we can assume, without loss of generality, a set of disconnected bands symmetric around zero energy. P-h symmetry implies that all band crossings at non-zero energy come in pair such that the analysis of the parity of can be restricted to zero-energy modes. A single Dirac cone is protected by and is pinned at zero energy by p-h sP. p-h further protects the parity of for zero modes as is varied. By continuity with the case where we have three Dirac cones and , we finally obtain that is odd, in contradiction with , which completes our proof that all bands must be connected.

In summary, we showed that trilayer twisted graphene exhibits band flattening along symmetry lines and close to magic angles. We also proved, by compatibility of band representations for evenly twisted planes or by counting an odd number of Dirac points protected by , that the system is always a metal with an infinite connectivity, an unprecedented feature in standard materials Bradlyn et al. (2017); Vergniory et al. (2018); Zhang et al. (2018); Tang et al. (2018). This property relies on p-h symmetry emerging for small twisting angles. We checked that this condition is practically realized already for angles close to the first magic angle ^{1}^{1}1See Supplemental Material with the numerical investigation of p-h symmetry breaking, more details on the symmetries and the positions of the Dirac points at .. Since it originates from the underlying three Dirac cones, we conjecture that the property of infinite band connectivity will appear in many other configuration such as multilayer moiré graphene with an odd number of planes.

We would like to thank B. Estienne and Biao Lian for fruitful discussions. BB is supported by the Department of Energy Grant No. de- sc0016239, the National Science Foundation EAGER Grant No. noaawd1004957, Simons Investigator Grants No. ONRN00014-14-1-0330, No. ARO MURI W911NF- 12-1-0461, and No. NSF-MRSEC DMR- 1420541, the Packard Foundation, the Schmidt Fund for Innovative Research.

## References

- Pong and Durkan (2005) Wing-Tat Pong and Colm Durkan, “A review and outlook for an anomaly of scanning tunnelling microscopy (stm): superlattices on graphite,” Journal of Physics D: Applied Physics 38, R329 (2005).
- Bistritzer and MacDonald (2010) R. Bistritzer and A. H. MacDonald, “Transport between twisted graphene layers,” Phys. Rev. B 81, 245412 (2010).
- Bistritzer and MacDonald (2011) Rafi Bistritzer and Allan H MacDonald, “Moiré bands in twisted double-layer graphene,” Proceedings of the National Academy of Sciences 108, 12233–12237 (2011).
- Trambly de LaissardiÃ¨re et al. (2010) G. Trambly de LaissardiÃ¨re, D. Mayou, and L. Magaud, “Localization of dirac electrons in rotated graphene bilayers,” Nano Letters 10, 804–808 (2010).
- San-Jose et al. (2012) P. San-Jose, J. González, and F. Guinea, “Non-abelian gauge potentials in graphene bilayers,” Phys. Rev. Lett. 108, 216802 (2012).
- Lopes dos Santos et al. (2012) J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Castro Neto, “Continuum model of the twisted graphene bilayer,” Phys. Rev. B 86, 155449 (2012).
- Fang and Kaxiras (2016) Shiang Fang and Efthimios Kaxiras, ‘‘Electronic structure theory of weakly interacting bilayers,” Phys. Rev. B 93, 235153 (2016).
- Tarnopolsky et al. (2018) Grigory Tarnopolsky, Alex J Kruchkov, and Ashvin Vishwanath, “Origin of magic angles in twisted bilayer graphene,” arXiv preprint arXiv:1808.05250 (2018).
- Kim et al. (2017) Kyounghwan Kim, Ashley DaSilva, Shengqiang Huang, Babak Fallahazad, Stefano Larentis, Takashi Taniguchi, Kenji Watanabe, Brian J. LeRoy, Allan H. MacDonald, and Emanuel Tutuc, “Tunable moiré bands and strong correlations in small-twist-angle bilayer graphene,” Proceedings of the National Academy of Sciences 114, 3364–3369 (2017).
- Cao et al. (2018a) Yuan Cao, Valla Fatemi, Ahmet Demir, Shiang Fang, Spencer L Tomarken, Jason Y Luo, Javier D Sanchez-Yamagishi, Kenji Watanabe, Takashi Taniguchi, Efthimios Kaxiras, et al., “Correlated insulator behaviour at half-filling in magic-angle graphene superlattices,” Nature 556, 80 (2018a).
- Cao et al. (2018b) Yuan Cao, Valla Fatemi, Shiang Fang, Kenji Watanabe, Takashi Taniguchi, Efthimios Kaxiras, and Pablo Jarillo-Herrero, “Unconventional superconductivity in magic-angle graphene superlattices,” Nature 556, 43 (2018b).
- Carr et al. (2018) Stephen Carr, Shiang Fang, Pablo Jarillo-Herrero, and Efthimios Kaxiras, “Pressure dependence of the magic twist angle in graphene superlattices,” Phys. Rev. B 98, 085144 (2018).
- Yankowitz et al. (2018) Matthew Yankowitz, Shaowen Chen, Hryhoriy Polshyn, K Watanabe, T Taniguchi, David Graf, Andrea F Young, and Cory R Dean, ‘‘Tuning superconductivity in twisted bilayer graphene,” arXiv preprint arXiv:1808.07865 (2018).
- Yuan and Fu (2018) Noah F. Q. Yuan and Liang Fu, “Model for the metal-insulator transition in graphene superlattices and beyond,” Phys. Rev. B 98, 045103 (2018).
- Wu et al. (2018) Fengcheng Wu, A. H. MacDonald, and Ivar Martin, “Theory of phonon-mediated superconductivity in twisted bilayer graphene,” Phys. Rev. Lett. 121, 257001 (2018).
- Padhi et al. (2018) Bikash Padhi, Chandan Setty, and Philip W. Phillips, “Doped twisted bilayer graphene near magic angles: Proximity to wigner crystallization, not mott insulation,” Nano Letters 18, 6175–6180 (2018), pMID: 30185049, https://doi.org/10.1021/acs.nanolett.8b02033 .
- Po et al. (2018a) Hoi Chun Po, Liujun Zou, Ashvin Vishwanath, and T. Senthil, “Origin of mott insulating behavior and superconductivity in twisted bilayer graphene,” Phys. Rev. X 8, 031089 (2018a).
- Isobe et al. (2018) Hiroki Isobe, Noah F. Q. Yuan, and Liang Fu, “Unconventional superconductivity and density waves in twisted bilayer graphene,” Phys. Rev. X 8, 041041 (2018).
- Dodaro et al. (2018) J. F. Dodaro, S. A. Kivelson, Y. Schattner, X. Q. Sun, and C. Wang, “Phases of a phenomenological model of twisted bilayer graphene,” Phys. Rev. B 98, 075154 (2018).
- Baskaran (2018) G. Baskaran, “Theory of Emergent Josephson Lattice in Neutral Twisted Bilayer Graphene (Moi’re is Different),” ArXiv e-prints (2018), arXiv:1804.00627 [cond-mat.supr-con] .
- Huang et al. (2018) T. Huang, L. Zhang, and T. Ma, “Antiferromagnetically ordered Mott insulator and superconductivity in twisted bilayer graphene: A quantum Monte carlo study,” ArXiv e-prints (2018), arXiv:1804.06096 [cond-mat.supr-con] .
- You and Vishwanath (2018) Y.-Z. You and A. Vishwanath, “Superconductivity from Valley Fluctuations and Approximate SO(4) Symmetry in a Weak Coupling Theory of Twisted Bilayer Graphene,” ArXiv e-prints (2018), arXiv:1805.06867 [cond-mat.str-el] .
- Wu et al. (2018) X.-C. Wu, K. A. Pawlak, C.-M. Jian, and C. Xu, “Emergent Superconductivity in the weak Mott insulator phase of bilayer Graphene Moir’e Superlattice,” ArXiv e-prints (2018), arXiv:1805.06906 .
- Roy and Juricic (2018) B. Roy and V. Juricic, “Unconventional superconductivity in nearly flat bands in twisted bilayer graphene,” ArXiv e-prints (2018), arXiv:1803.11190 [cond-mat.mes-hall] .
- Xu and Balents (2018) Cenke Xu and Leon Balents, “Topological superconductivity in twisted multilayer graphene,” Phys. Rev. Lett. 121, 087001 (2018).
- Zhang et al. (2018) Y.-H. Zhang, D. Mao, Y. Cao, P. Jarillo-Herrero, and T. Senthil, ‘‘Moir’e Superlattice with Nearly Flat Chern Bands: Platform for (Fractional) Quantum Anomalous Hall Effects and Unconventional Superconductivity,” ArXiv e-prints (2018), arXiv:1805.08232 [cond-mat.str-el] .
- Koshino et al. (2018) Mikito Koshino, Noah F. Q. Yuan, Takashi Koretsune, Masayuki Ochi, Kazuhiko Kuroki, and Liang Fu, “Maximally localized wannier orbitals and the extended hubbard model for twisted bilayer graphene,” Phys. Rev. X 8, 031087 (2018).
- Kang and Vafek (2018) Jian Kang and Oskar Vafek, ‘‘Symmetry, maximally localized wannier states, and a low-energy model for twisted bilayer graphene narrow bands,” Phys. Rev. X 8, 031088 (2018).
- Ray and Das (2018) Sujay Ray and Tanmoy Das, “Wannier pairs in the superconducting twisted bilayer graphene and related systems,” arXiv preprint arXiv:1804.09674 (2018).
- Guinea and Walet (2018) Francisco Guinea and Niels R Walet, “Electrostatic effects and band distortions in twisted graphene bilayers,” arXiv preprint arXiv:1806.05990 (2018).
- Lian et al. (2018) Biao Lian, Zhijun Wang, and B Andrei Bernevig, “Twisted bilayer graphene: A phonon driven superconductor,” arXiv preprint arXiv:1807.04382 (2018).
- Kang and Vafek (2018) Jian Kang and Oskar Vafek, “Strong coupling phases of partially filled twisted bilayer graphene narrow bands,” arXiv e-prints , arXiv:1810.08642 (2018), arXiv:1810.08642 [cond-mat.str-el] .
- Lian et al. (2018) Biao Lian, Fang Xie, and B. Andrei Bernevig, ‘‘The Landau Level of Fragile Topology,” arXiv e-prints , arXiv:1811.11786 (2018), arXiv:1811.11786 [cond-mat.mes-hall] .
- San-Jose and Prada (2013) Pablo San-Jose and Elsa Prada, “Helical networks in twisted bilayer graphene under interlayer bias,” Phys. Rev. B 88, 121408 (2013).
- de Gail et al. (2011) R. de Gail, M. O. Goerbig, F. Guinea, G. Montambaux, and A. H. Castro Neto, “Topologically protected zero modes in twisted bilayer graphene,” Phys. Rev. B 84, 045436 (2011).
- Po et al. (2018b) Hoi Chun Po, Liujun Zou, T. Senthil, and Ashvin Vishwanath, “Faithful Tight-binding Models and Fragile Topology of Magic-angle Bilayer Graphene,” arXiv:1808.02482 [cond-mat] (2018b), arXiv: 1808.02482.
- Song et al. (2018) Zhida Song, Zhijun Wang, Wujun Shi, Gang Li, Chen Fang, and B Andrei Bernevig, “All” magic angles” are” stable” topological,” arXiv preprint arXiv:1807.10676 (2018).
- Ahn et al. (2018) Junyeong Ahn, Sungjoon Park, and Bohm-Jung Yang, “Failure of nielsen-ninomiya theorem and fragile topology in two-dimensional systems with space-time inversion symmetry: application to twisted bilayer graphene at magic angle,” arXiv preprint arXiv:1808.05375 (2018).
- Zou et al. (2018) Liujun Zou, Hoi Chun Po, Ashvin Vishwanath, and T. Senthil, ‘‘Band structure of twisted bilayer graphene: Emergent symmetries, commensurate approximants, and wannier obstructions,” Phys. Rev. B 98, 085435 (2018).
- Liu et al. (2018) Jianpeng Liu, Junwei Liu, and Xi Dai, “A complete picture for the band topology in twisted bilayer graphene,” arXiv preprint arXiv:1810.03103 (2018).
- Hejazi et al. (2018) Kasra Hejazi, Chunxiao Liu, Hassan Shapourian, Xiao Chen, and Leon Balents, ‘‘Multiple topological transitions in twisted bilayer graphene near the first magic angle,” arXiv preprint arXiv:1808.01568 (2018).
- Lopes dos Santos et al. (2007) J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Castro Neto, “Graphene bilayer with a twist: Electronic structure,” Phys. Rev. Lett. 99, 256802 (2007).
- Bradlyn et al. (2017) Barry Bradlyn, L. Elcoro, Jennifer Cano, M. G. Vergniory, Zhijun Wang, C. Felser, M. I. Aroyo, and B. Andrei Bernevig, “Topological quantum chemistry,” Nature 547, 298–305 (2017).
- Vergniory et al. (2017) M. G. Vergniory, L. Elcoro, Zhijun Wang, Jennifer Cano, C. Felser, M. I. Aroyo, B. Andrei Bernevig, and Barry Bradlyn, “Graph theory data for topological quantum chemistry,” Phys. Rev. E 96, 023310 (2017).
- Amorim and Castro (2018) B Amorim and Eduardo V Castro, “Electronic spectral properties of incommensurate twisted trilayer graphene,” arXiv preprint arXiv:1807.11909 (2018).
- Gallego et al. (2012) Samuel V Gallego, Emre S Tasci, G Flor, J Manuel Perez-Mato, and Mois I Aroyo, “Magnetic symmetry in the bilbao crystallographic server: a computer program to provide systematic absences of magnetic neutron diffraction,” Journal of Applied Crystallography 45, 1236–1247 (2012).
- Vergniory et al. (2018) M. G. Vergniory, L. Elcoro, C. Felser, B. A. Bernevig, and Z. Wang, “The (high quality) topological materials in the world,” arXiv preprint arXiv:1807.10271 (2018).
- Zhang et al. (2018) Tiantian Zhang, Yi Jiang, Zhida Song, He Huang, Yuqing He, Zhong Fang, Hongming Weng, and Chen Fang, ‘‘Catalogue of topological electronic materials,” arXiv preprint arXiv:1807.08756 (2018).
- Tang et al. (2018) Feng Tang, Hoi Chun Po, Ashvin Vishwanath, and Xiangang Wan, “Towards ideal topological materials: Comprehensive database searches using symmetry indicators,” arXiv preprint arXiv:1807.09744 (2018).
- (50) See Supplemental Material with the numerical investigation of p-h symmetry breaking, more details on the symmetries and the positions of the Dirac points at .
- Cao et al. (2016) Y. Cao, J. Y. Luo, V. Fatemi, S. Fang, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, “Superlattice-induced insulating states and valley-protected orbits in twisted bilayer graphene,” Phys. Rev. Lett. 117, 116804 (2016).

Supplementary information for: Flat bands and perfect metal in trilayer moiré graphene

## Appendix S-I Breaking particle-hole symmetry

The particle-hole (p-h) symmetry is important in proving that all bands in the spectrum of trilayer twisted graphene are connected. It gives a band landscape symmetric with respect to zero energy and implies that features, such as Dirac points and other gapless crossings, or symmetries occurring at finite energy always come in pairs. It maintains globally the gapleass structure of the spectrum. The other essential symmetry is the combined spatial-time symmetry - this symmetry locally maintains the gapless structure of the band spectrum. Breaking directly open gaps at Dirac points locally and hence isolates finite sets of bands. In contrast to that, we expect the gapless band structure to keep its integrity if the terms breaking p-h symmetry are not too strong.

In a real trilayer graphene system, the rotation angles between the planes are finite and there are corrections to the emergent Eq. (1) of the main text. There are two main effects that break p-h symmetry: (i) corrections to the Dirac cones approximation close to the points in each layer, (ii) the three K points are not aligned at finite rotation angle. We note that (i) is the first step towards the complete band structure of each graphene layer with both the K and K’ points. The vicinities of K and K’ are indeed fully decoupled in Eq. (1) of the main text. (ii) can be simply taken into account by choosing a different basis in each layer such that the vectors are parallel. This results into the Hamiltonian

(S-11) |

where is the in-plane rotation of angle , reducing to Eq. (1) of the main text for vanishingly small twist angles. Eq. (S-11) explicitely breaks p-h symmetry. The numerical resolution of Eq. (S-11) is displayed in Fig. S-3 for the representative case and for several values of . Inspecting the states close to the Fermi energy along symmetry lines, we find band crossings between all states despite the breaking of p-h symmetry.

## Appendix S-II Symmetries for arbitrary and

We first review the particular case . The effective single-valley model of Eq. (1) (see main text) is obtained by keeping only the states close to the K point in each layer. It is decoupled from its time-reversed counterpart built with the states around the K’ points. Hence, Eq. (1) does not respect time-reversal symmetry. For , the symmetry group of this one-valley model is the magnetic space group , the same as for twisted bilayer graphene, characterized by the generators

(S-12) |

The band spectrum is invariant with respect to these three symmetries. The antiunitary symmetry is a combination of time-reversal and symmetry acting locally on the moiré lattice. It commutes with the spatial symmetries and . Whereas maps each layer into itself, exchanges the layers and and leaves the layer invariant.

Besides the magnetic space group , the model also exhibits a particle-hole (p-h) symmetry inherited from the original p-h symmetry in each graphene layer. The corresponding operator - see Eq. (7) and (8) in the main text -

(S-13) |

commutes with the generators (S-12) of the magnetic space group, and squares to in contrast with the p-h operator in the bilayer case which squares to .

The rotation symmetry and the antiunitary symmetry are valid also for arbitrary and , i.e. when the rotation angles are different but still commensurate. Interestingly, if or is different from , the symmetries and are individually broken but their product remains a symmetry of the one-valley Hamiltonian (4) (see main text). We thus introduce the mirror p-h symmetry operator

(S-14) |

obtained by combining the product with the rotation symmetry . is for belonging to the bottom and top layers and for the middle layer. The invariant plane of the mirror symmetry crosses orthogonally the graphene layers along lines. These lines make an angle of with respect to the axis and go through the and points. By construction, each layer is thus invariant upon the mirror symmetry . The action on the Hamiltonian (4) is

(S-15) |

corresponding to a p-h symmetry. Eq. (S-15) implies a symmetric band spectrum around zero energy where each state has a mirror-symmetric partner with opposite energy.

## Appendix S-III Positions of the Dirac cones for

As discussed in the main text, the full connectivity of the moiré band model requires essentially the symmetry to protect Dirac cones, or at least to have them gapped by pairs, the mirror p-h symmetry to reduce the parity of Dirac points to their parity in the zero energy manifold, and the continuity with respect to the case where the positions of Dirac points are easily determined. As illustrated in Fig. S-4, we find two configurations at : either (i) the model has three Dirac points at distinct positions in the Brillouin zone, or (ii) two Dirac points are on top of each other and separated from the third one. The case (i) occurs for whereas (ii) occurs for .

## Appendix S-IV Densities of electrons close to the Fermi energy

Twisted bilayer graphene, with a small angle, exhibits a local density of states for the (almost) flat bands that is well-localized to the AA regions of the Moiré pattern Trambly de LaissardiÃ¨re et al. (2010); Fang and Kaxiras (2016); Cao et al. (2016), forming a triangular lattice. This spatial localization is understood at zero energy from the absence of tunneling between AA and AB, BA regions Lopes dos Santos et al. (2012). In the trilayer geometry, it is not possible to isolate a finite set of bands as we prove in the main text. However, we can integrate the local density of states (LDOS) over a narrow energy range close to zero energy similarily to what has been done in the bilayer case Lopes dos Santos et al. (2012). The result is shown in Fig. S-5 for . In the middle layer, localization around the region occurs more dramatically when bands are nearly dispersionless close, see . Although the whole pattern forms a triangular lattice, the localization around the region is never really strong for the lower and top layers.