Macroscopically degenerate localized zero-energy states of quasicrystalline bilayer systems in strong coupling limit 1 1, 2, 3, Hyunsoo Ha and Bohm-Jung Yang ∗ 1Department of Physics and Astronomy, Seoul National University, Seoul 08826, Korea 2Center for Correlated Electron Systems, Institute for Basic Science (IBS), Seoul 08826, Korea 3Center for Theoretical Physics (CTP), Seoul National University, Seoul 08826, Korea (Dated: March 30, 2021) When two identical two-dimensional (2D) periodic lattices are stacked in parallel after rotating one layer by a certain angle relative to the other layer, the resulting bilayer system can lose lattice periodicity completely and become a 2D quasicrystal. Twisted bilayer graphene with 30◦ rotation is a representative example. We show that such quasicrystalline bilayer systems generally develop macroscopically degenerate localized zero-energy states (ZESs) in strong coupling limit where the interlayer couplings are overwhelmingly larger than the intralayer couplings. The emergent chiral symmetry in strong coupling limit and aperiodicity of bilayer quasicrystals guarantee the existence of the ZESs. The macroscopically degenerate ZESs are analogous to the flat bands of periodic systems, in that both are composed of localized eigenstates, which give divergent density of states. For monolayers, we consider the triangular, square, and honeycomb lattices, comprised of homogenous tiling of three possible planar regular polygons: the equilateral triangle, square, and regular hexagon. We construct a compact theoretical framework, which we call the quasiband model, that describes the low energy properties of bilayer quasicrystals and counts the number of ZESs using a subset of Bloch states of monolayers. We also propose a simple geometric scheme in real space which can show the spatial localization of ZESs and count their number. Our work clearly demonstrates that bilayer quasicrystals in strong coupling limit are an ideal playground to study the intriguing interplay of flat band physics and the aperiodicity of quasicrystals. Introduction.— Flat bands with divergent density of states (a) (b) (c) are ideal playgrounds to explore strong correlation physics [1– 7]. Recent discovery of twisted bilayer graphene (TBG) has provided a new platform for studying flat bands. In TBG, 휋/6 휋/4 휋/6 where two graphene layers are stacked with relative inplane rotation by an angle θ, the electronic properties change dra- matically depending on θ. Although TBG is generally not periodic due to incommensurate interlayer potentials, when θ is small (θ < 10◦), its low energy properties can be well described by Moire effective theory with translational invari- FIG. 1. Lattice structures of (a) twisted bilayer triangular lattice ance. Especially, when θ takes special values, called the (TBTL) quasicrystal with θ = 30◦, (b) twisted bilayer square lattice (TBSL) quasicrystal with θ = 45◦, (c) twisted bilayer honeycomb magic angles, the Dirac cones of graphene layers coupled by lattice quasicrystal with θ = 30◦. the Moire superlattice potential are strongly renormalized and turn into nearly flat bands [8–12], which induce intriguing strong correlation physics [13–15]. riodic interlayer hopping. To establish the general properties Another interesting limit of TBG is when θ = 30 exactly. ◦ of 2D quasicrystalline bilayers beyond the TBG with θ = 30 , In this case, the system completely loses translational symme- ◦ we take the triangular, square, and honeycomb lattices as a try [16, 17], but instead, develops a quasicrystalline order with monolayer structure, and compose the corresponding bilayer 12-fold rotation symmetry. Such an extrinsic quasicrystal of quasicrystals with the rotation angle θ = 30 , 45 , and 30 , TBG, composed of two periodic graphene layers [18, 19], was ◦ ◦ ◦ respectively [see Fig. 1]. We note that these monolayers are recently realized in experiment [20]. Angle-resolved photoe- made of homogenous tiling of three allowed planar regular mission spectroscopy has shown many replicas of Dirac cones polygons, i.e., the equilateral triangles, squares, and regular with 12-fold rotational symmetry in TBG quasicrystal, aris- hexagons. Surprisingly, we find that quasicrystalline bilayers arXiv:2103.08851v2 [cond-mat.mes-hall] 28 Mar 2021 ing from the interlayer coupling with the quasicrystalline or- in strong coupling limit generally host macroscopically de- der [20]. However, most of the previous studies are limited generate localized zero-energy states (ZESs). The chiral sym- to the weak coupling limit where the aperiodic interlayer hop- metry related to layer degrees of freedom, which emerges in ping is much smaller than the intralayer one, so that the low strong coupling limit, guarantees the presence of the ZESs. energy properties of the TBG quasicrystal are governed by We note that the macroscopically degenerate ZESs are similar massless Dirac fermions as in monolayer graphene [21, 22]. to the flat bands of periodic lattices, both of which induce di- In this Letter, we unveil the low energy properties of qua- vergent density of states (DOS). Hence the quasicrystalline bi- sicrystalline bilayer systems in strong coupling limit where layer in strong coupling limit provides a fascinating platform the intralayer hopping is negligibly small compared to the ape- to study the interplay of the flat energy dispersion and the qua- 2 sicrystalline order, which is distinct from the weak coupling (a) (b) physics dominated by masselss Dirac fermions. region Mixed state ۧ ۧ Twisted bilayer triangular lattice (TBTL) quasicrystal.— ● A ห푞 1, ห푞 2 B ห푞ۧ , ห푞 + 퐺 ۧ , ห푞ۧ ● We first consider the quasicrystal comprised of twisted bilayer 1 2 1 2 ۧ ۧ ۧ 푮ퟐ ● C ห푞 1, ห푞 2, ห푞 + 퐺1 2 triangular lattices (TBTLs) with θ = 30◦. For comparison 풒 + 푮ퟐ 풒 (c) with TBG quasicrystal, we construct a model where a pz or- .) bital is placed at each lattice site. The full Hamiltonian is a.u gap given by H = H1 + H2 + H12 where Hi=1,2 describes in- tralayer hoppings in the layer i while H indicates interlayer Zero energy states 12 D.O.S.( H ㅡ Layer1 ㅡ Layer2 ㅡ UBZ hoppings. In strong coupling limit, we neglect i=1,2 whose q ∈ layer1 1st BZ ∩ layer2 1st BZ Energy(eV) influence is discussed later. The resulting Hamiltonian is 풍 풍 (d) (e) 풕ퟐ 풕ퟏ 푩풕 푻ퟏ H = T (R1 R2) R1 R2 + h.c. (1) − − | ih | 풀풕 R1,R2 X 푨풕 푿풕 풁풕 where R ( R ) indicates the p orbital located at the po- 풍풕ퟑ | 1i | 2i z sition R (R ) in the layer 1 (layer 2). T (R R ) is the 풒풄 1 2 1 − 2 transfer integral between two orbitals at R1 and R2, respec- tively, generally taking the following form [40], DFR R e 2 R e 2 T (R) = V 1 · z + V · z , − ppπ − R ppσ R FIG. 2. (a) The first Brillouin zone (BZ) of the layer 1 and 2, and the " # union of BZs (UBZ) taken as the effective BZ (EBZ). EBZ is divided 0 ( R a0)/r? 0 ( R d0)/r? Vppπ = Vppπe − − ,Vppσ = Vppσe − − , (2) into three regions A (yellow), B (blue), C (orange). (b) Bloch states coupled in each region. (c) Density of states (DOS) from the full tight-binding model and the quasiband model. (d) Regions B and C where a0 is the inplane lattice constant, d0 denotes the vertical distance between two layers, R is the relative displacement illustrated in the extended BZ scheme. The dark filled region (DFR) where regions B and C overlap can couple ZESs in the EBZ. (e) vector between two atoms, R = R , e denotes the unit vec- | | z Real space geometric consideration for ZESs. Orange dots denote tor in the z-direction. Vppπ (Vppσ) indicates the transfer inte- the layer 2 where ZESs are confined. The small triangle with blue gral between nearest-neighbor pz orbitals forming π-bonding dots inside indicates the possible center positions of localized ZESs (σ-bonding). The length scale r? determines the range of hop- with maximum amplitude. ping transfer. For comparision to TBG quasicrystal, we as- sume a0 0.142nm, d0 0.335 nm, r? = r0 0.0453 nm, 0 ≈ 0≈ ≈ j Vppπ 2.6 eV, and Vppσ 0.48 eV. Here r0 indicates the vector of the th layer. Eq. (3) show that the interlayer Hamil- ≈ − ≈ tonian couples Bloch states k (q) and k (q) satisfying interaction range in graphene [21, 35, 36]. | 1 i1 | 2 i2 Interestingly, the strong coupling Hamiltonian in Eq. (1) the generalized Umklapp scattering condition, k1 + G1 = k + G q. The coupling strength is proportional to t(q) has chiral symmetry related to layer degrees represented by 2 2 ≡ | | the operator U = R1 R1 R2 R2 . As which decays exponentially as q increases, so that the low R1 | ih | − R2 | ih | | | H, U = 0, for an eigenstate ψ satisfying H ψ = E ψ , energy properties of the Hamiltonian can be described by the { } | i | i | i one can always find anotherP eigenstate UPψ with the energy momentum region with small q . This motivate us to intro- | i | | E so that the energy spectrum is particle-hole symmetric. duce the quasiband model, composed of the Bloch states with − Quasiband model.— Though translational invariance is lost small momentum within the effective BZ (EBZ). Here we take in the bilayer quasicrystal, its low energy properties can be the union of the first BZs (UBZ) of two monolayers as the successfully described by the Bloch states of monolayers EBZ, and neglect the states with the momentum outside the 1 ikj Rj kj j = e · Rj (j = 1, 2) where kj indicates EBZ. | i √N Rj | i Explicitly, the quasiband model HQB = q EBZ H(q) the momentum in the Brillouin zone (BZ) of the jth layer and ∈ P with H(q) = [t(q) k1(q) 12 k2(q) + h.c.] can be divided N is the number of lattice sites per layer. The subscript j under | i h | P ket (or bra) is the layer index of the Bloch basis.