Researcher profile

B. Andrei Bernevig

B. Andrei Bernevig contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
41works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2022arXiv

All Topological Bands of All Nonmagnetic Stoichiometric Materials

Topological Quantum Chemistry and Symmetry-Based Indicators have facilitated large-scale searches for materials with topological properties at the Fermi energy ($E_{F}$). We report the completion of a publicly accessible catalog of stable and fragile topology in all of the bands both at and away from $E_{F}$ in the 96,196 processable entries in the Inorganic Crystal Structure Database. Our calculations represent the completion of the symmetry-indicated band topology of known nonmagnetic materials, and enable the discovery of repeat-topological and supertopological materials, which include rhombohedral bismuth and Bi$_2$Mg$_3$. We find that 52.65% of all materials are topological at $E_{F}$, roughly $2/3$ of bands across all materials exhibit symmetry-indicated stable topology, and that 87.99% of all materials contain at least one stable or fragile topological band.

preprint2022arXiv

Catalogue of Flat-Band Stoichiometric Materials

Topological electronic flatten bands near or at the Fermi level are a promising avenue towards unconventional superconductivity and correlated insulating states. However, the related experiments are mostly limited to the engineered materials, such as moire systems. Here we present a catalogue of all the three-dimensional stoichiometric materials with flat bands around the Fermi level that exist in nature. We consider 55,206 materials from the Inorganic Crystal Structure Database catalogued using the Topological Quantum Chemistry website which provides their structural parameters, space group (SG), band structure, density of states and topological characterization. We combine several direct signatures and properties of band flatness to a high-throughput analysis of all crystal structures. In particular, we identify materials hosting line-graph or bipartite sublattices - either in two or three dimensions - likely leading to flat bands. From this trove of information, we create the Materials Flatband Database website, a powerful search engine for future theoretical and experimental studies. We use it to extract a curated list of 2,379 materials, with among them 345 promising candidates, potentially hosting flat bands whose charge centers are not strongly localized on the atomic sites. We showcase five representative materials: KAg[CN]2 in SG 163 $(P\bar{3}1c)$, Pb2Sb2O7 in SG 227 $(Fd\bar{3}m)$, Rb2CaH4 in SG 139 $(I4/mmm)$, Ca2NCl in SG 166 $(R\bar{3}m)$ and WO3 in SG 221 $(Pm\bar{3}m)$. We provide a theoretical explanation for the origin of their flat bands close to the Fermi energy using the $S$-matrix method introduced in a parallel work [Calugaru et al., Nature Physics 18, 185 (2022)].

preprint2022arXiv

Evidence for a Monolayer Excitonic Insulator

The interplay between topology and correlations can generate a variety of quantum phases, many of which remain to be explored. Recent advances have identified monolayer WTe2 as a promising material for doing so in a highly tunable fashion. The ground state of this two-dimensional (2D) crystal can be electrostatically tuned from a quantum spin Hall insulator (QSHI) to a superconductor. However, much remains unknown about the gap-opening mechanism of the insulating state. Here we report evidence that the QSHI is also an excitonic insulator (EI), arising from the spontaneous formation of electron-hole bound states (excitons). We reveal the presence of an intrinsic insulating state at the charge neutrality point (CNP) in clean samples and confirm the correlated nature of this charge-neutral insulator by tunneling spectroscopy. We provide evidence against alternative scenarios of a band insulator or a localized insulator and support the existence of an EI phase in the clean limit. These observations lay the foundation for understanding a new class of correlated insulators with nontrivial topology and identify monolayer WTe2 as a promising candidate for exploring quantum phases of ground-state excitons.

preprint2022arXiv

Exact Quantum Scars in the Chiral Non-Linear Luttinger Liquid

While the chiral linear Luttinger liquid is integrable via bosonization, its non-linear counterpart does not admit for an analytic solution. In this work, we find a sub-extensive number of exact eigenstates for a large family of density-density interaction terms. These states are embedded in a continuum of strongly-correlated excited states. The real-space entanglement entropy of some exact states scales logarithmically with system size, while that of others has volume-law scaling. We introduce momentum-space entanglement as an unambiguous differentiator between these exact states and the remaining excited states. With regard to momentum space, the exact states behave as bona fide quantum many body scars: they exhibit identically zero momentum-space entanglement, while typical eigenstates behave thermally. We corroborate this finding by a level statistics analysis. Furthermore, we detail the general formalism for systematically finding all interaction terms and associated exact states, and present a number of infinite exact state sequences extending to arbitrarily high energies. Unlike many previous examples of quantum many body scars, the exact states uncovered here do not lie at equidistant energies, and do not follow from a special operator algebra. Instead, they are uniquely enabled by the interplay of Fermi statistics and chirality.

preprint2022arXiv

Magnetic Bloch Theorem and Reentrant Flat Bands in Twisted Bilayer Graphene at $2π$ Flux

Bloch's theorem is the centerpiece of topological band theory, which itself has defined an era of quantum materials research. However, Bloch's theorem is broken by a perpendicular magnetic field, making it difficult to study topological systems in strong flux. For the first time, moiré materials have made this problem experimentally relevant, and its solution is the focus of this work. We construct gauge-invariant irreps of the magnetic translation group at $2π$ flux on infinite boundary conditions, allowing us to give analytical expressions in terms of the Siegel theta function for the magnetic Bloch Hamiltonian, non-Abelian Wilson loop, and many-body form factors. We illustrate our formalism using a simple square lattice model and the Bistritzer-MacDonald Hamiltonian of twisted bilayer graphene, obtaining reentrant ground states at $2π$ flux under the Coulomb interaction.

preprint2022arXiv

Many-Body Superconductivity in Topological Flat Bands

In a flat band superconductor, bosonic excitations can disperse while unpaired electrons are immobile. To study this strongly interacting system, we construct a family of multi-band Hubbard models with exact eta-pairing ground states in all space groups. We analytically compute their many-body excitations and find that the Cooper pair bound states and density excitations obey an effective single-particle Hamiltonian written in terms of the non-interacting wavefunctions. These bound states possess a unique zero-energy excitation whose quadratic dispersion is determined by the minimal quantum metric. The rest of the bound state spectrum is classified by topological quantum chemistry, which we use to identify Cooper pairs with Weyl nodes, higher angular momentum pairing, and fragile topology. We also add electron kinetic energy as a perturbation to show that the strongest pairing occurs at half filling and not at the highest density of states. This is similar in spirit to the superconductivity observed in twisted bilayer graphene.

preprint2022arXiv

MATBG as Topological Heavy Fermion: I. Exact Mapping and Correlated Insulators

Magic-angle ($θ=1.05^\circ$) twisted bilayer graphene (MATBG) has shown two seemingly contradictory characters: the localization and quantum-dot-like behavior in STM experiments, and delocalization in transport experiments. We construct a model, which naturally captures the two aspects, from the Bistritzer-MacDonald (BM) model in a first principle spirit. A set of local flat-band orbitals ($f$) centered at the AA-stacking regions are responsible to the localization. A set of extended topological conduction bands ($c$), which are at small energetic separation from the local orbitals, are responsible to the delocalization and transport. The topological flat bands of the BM model appear as a result of the hybridization of $f$- and $c$-electrons. This model then provides a new perspective for the strong correlation physics, which is now described as strongly correlated $f$-electrons coupled to nearly free topological semimetallic $c$-electrons - we hence name our model as the topological heavy fermion model. Using this model, we obtain the U(4) and U(4)$\times$U(4) symmetries as well as the correlated insulator phases and their energies. Simple rules for the ground states and their Chern numbers are derived. Moreover, features such as the large dispersion of the charge $\pm1$ excitations and the minima of the charge gap at the $Γ_M$ point can now, for the first time, be understood both qualitatively and quantitatively in a simple physical picture. Our mapping opens the prospect of using heavy-fermion physics machinery to the superconducting physics of MATBG.

preprint2022arXiv

Progress and prospects in magnetic topological materials

Magnetic topological materials represent a class of compounds whose properties are strongly influenced by the topology of the electronic wavefunctions coupled with the magnetic spin configuration. Such materials can support chiral electronic channels of perfect conduction, and can be used for an array of applications from information storage and control to dissipationless spin and charge transport. Here, we review the theoretical and experimental progress achieved in the field of magnetic topological materials beginning with the theoretical prediction of the Quantum Anomalous Hall Effect without Landau levels, and leading to the recent discoveries of magnetic Weyl semimetals and antiferromagnetic topological insulators. We outline the recent theoretical progress that resulted in the tabulation, for the first time, of all magnetic symmetry group representations and topology. We describe several experiments realizing Chern insulators, Weyl and Dirac magnetic semimetals, and an array of axionic and higher-order topological phases of matter as well as survey future perspectives.

preprint2022arXiv

Quantum Many-Body Scars and Hilbert Space Fragmentation: A Review of Exact Results

The discovery of Quantum Many-Body Scars (QMBS) both in Rydberg atom simulators and in the Affleck-Kennedy-Lieb-Tasaki (AKLT) spin-1 chain model, have shown that a weak violation of ergodicity can still lead to rich experimental and theoretical physics. In this review, we provide a pedagogical introduction to and an overview of the exact results on weak ergodicity breaking via QMBS in isolated quantum systems with the help of simple examples such as the fermionic Hubbard model. We also discuss various mechanisms and unifying formalisms that have been proposed to encompass the plethora of systems exhibiting QMBS. We cover examples of equally-spaced towers that lead to exact revivals for particular initial states, as well as isolated examples of QMBS. Finally, we review Hilbert Space Fragmentation, a related phenomenon where systems exhibit a richer variety of ergodic and non-ergodic behaviors, and discuss its connections to QMBS.

preprint2022arXiv

Quasi-symmetry protected topology in a semi-metal

The crystal symmetry of a material dictates the type of topological band structures it may host, and therefore symmetry is the guiding principle to find topological materials. Here we introduce an alternative guiding principle, which we call 'quasi-symmetry'. This is the situation where a Hamiltonian has an exact symmetry at lower-order that is broken by higher-order perturbation terms. This enforces finite but parametrically small gaps at some low-symmetry points in momentum space. Untethered from the restraints of symmetry, quasi-symmetries eliminate the need for fine-tuning as they enforce that sources of large Berry curvature will occur at arbitrary chemical potentials. We demonstrate that a quasi-symmetry in the semi-metal CoSi stabilizes gaps below 2 meV over a large near-degenerate plane that can be measured in the quantum oscillation spectrum. The application of in-plane strain breaks the crystal symmetry and gaps the degenerate point, observable by new magnetic breakdown orbits. The quasi-symmetry, however, does not depend on spatial symmetries and hence transmission remains fully coherent. These results demonstrate a class of topological materials with increased resilience to perturbations such as strain-induced crystalline symmetry breaking, which may lead to robust topological applications as well as unexpected topology beyond the usual space group classifications.

preprint2022arXiv

Superfluid weight bounds from symmetry and quantum geometry in flat bands

Flat-band superconductivity has theoretically demonstrated the importance of band topology to correlated phases. In two dimensions, the superfluid weight, which determines the critical temperature through the Berezinksii-Kosterlitz-Thouless criteria, is bounded by the Fubini-Study metric at zero temperature. We show this bound is nonzero within flat bands whose Wannier centers are obstructed from the atoms - even when they have identically zero Berry curvature. Next, we derive general lower bounds for the superfluid weight in terms of momentum space irreps in all 2D space groups, extending the reach of topological quantum chemistry to superconducting states. We find that the bounds can be naturally expressed using the formalism of real space invariants (RSIs) that highlight the separation between electronic and atomic degrees of freedom. Finally, using exact Monte Carlo simulations on a model with perfectly flat bands and strictly local obstructed Wannier functions, we find that an attractive Hubbard interaction results in superconductivity as predicted by the RSI bound beyond mean-field. Hence, a nonzero superfluid weight constitutes a nontrivial bulk property that distinguishes obstructed bands from trivial bands in the presence of interactions.

preprint2022arXiv

Trions in Twisted Bilayer Graphene

The strong coupling phase diagram of magic angle twisted bilayer graphene (TBG) predicts a series of exact one particle charge $\pm 1$ gapped excitations on top of the integer-filled ferromagnetic ground-states. Finite-size exact diagonalization studies showed that these are the lowest charge $\pm 1$ excitations in the system (for $10$nm screening length), with the exception of charge $+1$ at filling $-1$ in the chiral limit. In the current paper we show that this "trion bound state", a $3$-particle, charge $1$ excitation of the insulating ferromagnetic ground-state of the projected Hamiltonian of TBG is the lowest charge $+1$ overall excitation at $ν=-1$, and also for some large ($\approx 20$nm) screening lengths at $ν=-2$ in the chiral limit and with very small binding energy. At other fillings, we show that trion bound states do exist, but only for momentum ranges that do not cover the entire moiré Brillouin zone. The trion bound states (at different momenta) exist for finite parameter range $w_0/w_1$ but they all disappear in the continuum far below the realistic values of $w_0/w_1= 0.8$. We find the conditions for the existence of the trion bound state, a good variational wavefunction for it, and investigate its behavior for different screening lengths, at all integer fillings, on both the electron and hole sides.

preprint2022arXiv

Twisted bilayer graphene I. Matrix elements, approximations, perturbation theory and a $k\cdot p$ 2-Band model

We investigate the Twisted Bilayer Graphene (TBG) model to obtain an analytic understanding of its energetics and wavefunctions needed for many-body calculations. We provide an approximation scheme which first elucidates why the BM $K_M$-point centered calculation containing only $4$ plane-waves provides a good analytical value for the first magic angle. The approximation scheme also elucidates why most many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a $Γ_M$-point centered model of 6 plane-waves, we analytically understand the small $Γ_M$-point gap between the active and passive bands in the isotropic limit $w_0=w_1$. Furthermore, we analytically calculate the group velocities of passive bands in the isotropic limit, and show that they are \emph{almost} doubly degenerate, while no symmetry forces them to be. Furthermore, away from $Γ_M$ and $K_M$ points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite it is defined only by vanishing $K_M$-point Dirac velocity. We derive analytically a connected "magic manifold" $w_1=2\sqrt{1+w_0^2}-\sqrt{2+3w_0^2}$, on which the bands remain extremely flat as $w_0$ is tuned between the isotropic ($w_0=w_1$) and chiral ($w_0=0$) limits. We analytically show why going away from the isotropic limit by making $w_0$ less (but not larger) than $w_1$ increases the $Γ_M$- point gap between active and passive bands. Finally, perturbatively, we provide an analytic $Γ_M$ point $k\cdot p$ $2$-band model that reproduces the TBG band structure and eigenstates in a certain $w_0,w_1$ parameter range. Further refinement of this model suggests a possible faithful $2$-band $Γ_M$ point $k\cdot p$ model in the full $w_0, w_1$ parameter range.

preprint2022arXiv

Twisted bilayer graphene III. Interacting Hamiltonian and exact symmetries

We derive the explicit Hamiltonian of twisted bilayer graphene (TBG) with Coulomb interaction projected into the flat bands, and study the symmetries of the Hamiltonian. First, we show that all projected TBG Hamiltonians can be written as Positive Semidefinite Hamiltonian, the first example of which was found in [PRL 122, 246401]. We then prove that the interacting TBG Hamiltonian exhibits an exact U(4) symmetry in the exactly flat band (nonchiral-flat) limit. We further define, besides a first chiral limit where the AA stacking hopping is zero, a new second chiral limit where the AB/BA stacking hopping is zero. In the first chiral-flat limit (or second chiral-flat limit) with exactly flat bands, the TBG is enhanced to have an exact U(4)$\times$U(4) symmetry, whose generators are different between the two chiral limits. While in the first chiral limit and in the non-chiral case these symmetries have been found in [PRX 10, 031034] for the $8$ lowest bands, we here prove that they are valid for projection into any $8 n_\text{max}$ particle-hole symmetric TBG bands, with $n_\text{max}>1$ being the practical case for small twist angles $<1^\circ$. Furthermore, in the first or second chiral-nonflat limit without flat bands, an exact U(4) symmetry still remains. We also elucidate the link between the U(4) symmetry presented here and the similar but different U(4) of [PRL 122, 246401]. Furthermore, we show that our projected Hamiltonian can be viewed as the normal-ordered Coulomb interaction plus a Hartree-Fock term from passive bands, and exhibits a many-body particle-hole symmetry which renders the physics symmetric around charge neutrality. We also provide an efficient parameterization of the interacting Hamiltonian. The existence of two chiral limits, with an enlarged symmetry suggests a possible duality of the model yet undiscovered.

preprint2022arXiv

Twisted Bilayer Graphene IV. Exact Insulator Ground States and Phase Diagram

We derive the exact insulator ground states of the projected Hamiltonian of magic-angle twisted bilayer graphene (TBG) flat bands with Coulomb interactions in various limits, and study the perturbations away from these limits. We define the (first) chiral limit where the AA stacking hopping is zero, and a flat limit with exactly flat bands. In the chiral-flat limit, the TBG Hamiltonian has a U(4)$\times$U(4) symmetry, and we find that the exact ground states at integer filling $-4\le ν\le 4$ relative to charge neutrality are Chern insulators of Chern numbers $ν_C=4-|ν|,2-|ν|,\cdots,|ν|-4$, all of which are degenerate. This confirms recent experiments where Chern insulators are found to be competitive low-energy states of TBG. When the chiral-flat limit is reduced to the nonchiral-flat limit which has a U(4) symmetry, we find $ν=0,\pm2$ has exact ground states of Chern number $0$, while $ν=\pm1,\pm3$ has perturbative ground states of Chern number $ν_C=\pm1$, which are U(4) ferromagnetic. In the chiral-nonflat limit with a different U(4) symmetry, different Chern number states are degenerate up to second order perturbations. In the realistic nonchiral-nonflat case, we find that the perturbative insulator states with Chern number $ν_C=0$ ($0<|ν_C|<4-|ν|$) at integer fillings $ν$ are fully (partially) intervalley coherent, while the insulator states with Chern number $|ν_C|=4-|ν|$ are valley polarized. However, for $0<|ν_C|\le4-|ν|$, the fully intervalley coherent states are highly competitive (0.005meV/electron higher). At nonzero magnetic field $|B|>0$, a first-order phase transition for $ν=\pm1,\pm2$ from Chern number $ν_C=\text{sgn}(νB)(2-|ν|)$ to $ν_C=\text{sgn}(νB)(4-|ν|)$ is expected, which agrees with recent experimental observations. Lastly, the TBG Hamiltonian reduces into an extended Hubbard model in the stabilizer code limit.

preprint2022arXiv

Twisted Bilayer Graphene V: Exact Analytic Many-Body Excitations in Twisted Bilayer Graphene Coulomb Hamiltonians: Charge Gap, Goldstone Modes and Absence of Cooper Pairing

We find exact analytic expressions for the energies and wavefunctions of the charged and neutral excitations above the exact ground states (at rational filling per unit cell) of projected Coulomb Hamiltonians in twisted bilayer graphene. Our exact expressions are valid for any form of the Coulomb interaction and any form of $AA$ and $AB/BA$ tunneling. The single charge excitation energy is a convolution of the Coulomb potential with a quantum geometric tensor of the TBG bands. The neutral excitations are (high-symmetry group) magnons, and their dispersion is analytically calculated in terms of the form factors of the active bands in TBG. The two-charge excitation energy and wavefunctions are also obtained, and a sufficient condition on the graphene eigenstates for obtaining a Cooper-pair from Coulomb interactions is obtained. For the actual TBG bands at the first magic angle, we can analytically show that the Cooper pair binding energy is zero in all such projected Coulomb models, implying that either phonons and/or non-zero kinetic energy are needed for superconductivity. Since the [Phys. Rev. Lett. 122, 246401] showed that the kinetic energy bounds on the superexchange energy are less $10^{-3}$ in Coulomb units, the phonon mechanism becomes then very likely. If nonetheless the superconductivity is due to kinetic terms which render the bands non-flat, one prediction of our theory is that the highest $T_c$ would not occur at the highest DOS.

preprint2022arXiv

Twisted bilayer graphene. II. Stable symmetry anomaly

We show that the entire continuous model of twisted bilayer graphene (TBG) (and not just the two active bands) with particle-hole symmetry is anomalous and hence incompatible with a lattice model. Previous works, e.g., [Phys. Rev. Lett. 123, 036401], [Phys. Rev. X 9, 021013], [Phys. Rev. B 99, 195455], and others [1-4] found that the two flat bands in TBG possess a fragile topology protected by the $C_{2z}T$ symmetry. [Phys. Rev. Lett. 123, 036401] also pointed out an approximate particle-hole symmetry ($\mathcal{P}$) in the continuous model of TBG. In this work, we numerically confirm that $\mathcal{P}$ is indeed a good approximation for TBG and show that the fragile topology of the two flat bands is enhanced to a $\mathcal{P}$-protected stable topology. This stable topology implies $4l+2$ ($l\in\mathbb{N}$) Dirac points between the middle two bands. The $\mathcal{P}$-protected stable topology is robust against arbitrary gap closings between the middle two bands the other bands. We further show that, remarkably, this $\mathcal{P}$-protected stable topology, as well as the corresponding $4l + 2$ Dirac points, cannot be realized in lattice models that preserve both $C_{2z}T$ and $\mathcal{P}$ symmetries. In other words, the continuous model of TBG is anomalous and cannot be realized on lattices. Two other topology related topics, with consequences for the interacting TBG problem, i.e., the choice of Chern band basis in the two flat bands and the perfect metal phase of TBG in the so-called second chiral limit, are also discussed.

preprint2022arXiv

Twisted bilayer graphene.VI. An Exact Diagonalization Study of Twisted Bilayer Graphene at Non-Zero Integer Fillings

Using exact diagonalization, we study the projected Hamiltonian with Coulomb interaction in the 8 flat bands of first magic angle twisted bilayer graphene. Employing the U(4) (U(4)$\times$U(4)) symmetries in the nonchiral (chiral) flat band limit, we reduced the Hilbert space to an extent which allows for study around $ν=\pm 3,\pm2,\pm1$ fillings. In the first chiral limit $w_0/w_1=0$ where $w_0$ ($w_1$) is the $AA$ ($AB$) stacking hopping, we find that the ground-states at these fillings are extremely well-described by Slater determinants in a so-called Chern basis, and the exactly solvable charge $\pm1$ excitations found in [arXiv:2009.14200] are the lowest charge excitations up to system sizes $8\times8$ (for restricted Hilbert space) in the chiral-flat limit. We also find that the Flat Metric Condition (FMC) used in [arXiv:2009.11301,2009.11872,2009.12376,2009.13530,2009.14200] for obtaining a series of exact ground-states and excitations holds in a large parameter space. For $ν=-3$, the ground state is the spin and valley polarized Chern insulator with $ν_C=\pm1$ at $w_0/w_1\lesssim0.9$ (0.3) with (without) FMC. At $ν=-2$, we can only numerically access the valley polarized sector, and we find a spin ferromagnetic phase when $w_0/w_1\gtrsim0.5t$ where $t\in[0,1]$ is the factor of rescaling of the actual TBG bandwidth, and a spin singlet phase otherwise, confirming the perturbative calculation [arXiv:2009.13530]. The analytic FMC ground state is, however, predicted in the intervalley coherent sector which we cannot access [arXiv:2009.13530]. For $ν=-3$ with/without FMC, when $w_0/w_1$ is large, the finite-size gap $Δ$ to the neutral excitations vanishes, leading to phase transitions. Further analysis of the ground state momentum sectors at $ν=-3$ suggests a competition among (nematic) metal, momentum $M_M$ ($π$) stripe and $K_M$-CDW orders at large $w_0/w_1$.

preprint2021arXiv

An $O(N)$ $Ab~initio$ Calculation Scheme for Large-Scale Moiré Structures

We present a two-step method specifically tailored for band structure calculation of the small-angle moiré-pattern materials which contain tens of thousands of atoms in a unit cell. In the first step, the self-consistent field calculation for ground state is performed with $O(N)$ Krylov subspace method implemented in OpenMX. Secondly, the crystal momentum dependent Bloch Hamiltonian and overlap matrix are constructed from the results obtained in the first step and only a small number of eigenvalues near the Fermi energy are solved with shift-invert and Lanczos techniques. By systematically tuning two key parameters, the cutoff radius for electron hopping interaction and the dimension of Krylov subspace, we obtained the band structures for both rigid and corrugated twisted bilayer graphene structures at the first magic angle ($θ=1.08^\circ$) and other three larger ones with satisfied accuracy on affordable costs. The band structures are in good agreement with those from tight binding models, continuum models, plane-wave pseudo-potential based $ab~initio$ calculations, and the experimental observations. This efficient two-step method is to play a crucial role in other twisted two-dimensional materials, where the band structures are much more complex than graphene and the effective model is hard to be constructed.

preprint2021arXiv

Cascades between light and heavy fermions in the normal state of magic angle twisted bilayer graphene

We present a framework for understanding the recently observed cascade transitions and the Landau level degeneracies at every integer filling of twisted bilayer graphene. The Coulomb interaction projected onto narrow bands causes the charged excitations at an integer filling to disperse, forming new bands. If the excitation moves the filling away from the charge neutrality point, then it has a band minimum at the moire Brillouin zone center with a small mass that compares well with the experiment; if towards the charge neutrality point, then it has a much larger mass and a higher degeneracy. At a non-zero density away from an integer filling the excitations interact. The system on the small mass side has a large bandwidth and forms a Fermi liquid. On the large mass side the bandwidth is narrow, the compressibility is negative and the Fermi liquid is likely unstable. This explains the observed sawtooth features in compressibility, the Landau fans pointing away from charge neutrality as well as their degeneracies. By providing a description of the charge itineracy in the normal state this framework sets the stage for superconductivity at lower temperatures.

preprint2021arXiv

General Construction and Topological Classification of All Magnetic and Non-Magnetic Flat Bands

Exotic phases of matter emerge from the interplay between strong electron interactions and non-trivial topology. Owing to their lack of dispersion at the single-particle level, systems harboring flat bands are excellent testbeds for strongly interacting physics, with twisted bilayer graphene serving as a prime example. On the other hand, existing theoretical models for obtaining flat bands in crystalline materials, such as the line-graph formalism, are often too restrictive for real-life material realizations. Here we present a generic technique for constructing perfectly flat bands from bipartite crystalline lattices. Our prescription encapsulates and generalizes the various flat band models in the literature, being applicable to systems with any orbital content, with or without spin-orbit coupling. Using Topological Quantum Chemistry, we build a complete topological classification in terms of symmetry eigenvalues of all the gapped and gapless flat bands, for all 1651 Magnetic Space Groups. In addition, we derive criteria for the existence of symmetry-protected band touching points between the flat and dispersive bands, and we identify the gapped flat bands as prime candidates for fragile topological phases. Finally, we show that the set of all (gapped and gapless) perfectly flat bands is finitely generated and construct the corresponding bases for all 1651 Shubnikov Space Groups.

preprint2021arXiv

Obstructed surface states as the origin of catalytic activity in inorganic heterogeneous catalysts

The discovery of new catalysts that are efficient, sustainable, and low-cost is a major research endeavor for many industrial chemical processes. This requires an understanding and determination of the catalytic origins for the given catalysts, which still remains a challenge. Here we describe a novel method to identify new catalysts based on searching for crystalline symmetry-protected obstructed atomic insulators (OAIs) that have metallic surface states on otherwise semiconducting or insulating compounds. The Wannier charge centers in OAIs are pinned by symmetries at some empty Wyckoff positions so that surfaces that accommodate these sites are guaranteed to have metallic obstructed surface states (OSSs). Beyond the well-studied 2H-MoS2, we further verified our theory on the catalysts, 2H-MoTe2, and 1T&#39;-MoTe2, whose catalytic active sites are consistent with our calculations of obstructed Wannier charge centers (OWCCs) and OSSs. In addition, we have predicted the location of catalytic active sites and confirmed these predictions by exploring the hydrogen evolution reaction on NiPS3 bulk single crystals, which we find to be one of the most promising new catalysts with high activity and, moreover, of low cost. Most importantly, we successfully identified several high-efficient catalysts just by considering the number of OWCCs and the crystal symmetry of the OAIs. Using the real space invariant theory and high-throughput computational methods applied to a database of 34013 topologically trivial insulators, we have identified 1788 unique OAIs (3383 ICSDs), of which 465 are potential high-quality catalysts for heterogeneous reactions. The Miller indices of the active surfaces are also obtained. Our new methodology will facilitate and accelerate the discovery of new catalysts for a wide range of heterogeneous redox reactions, where sustainability, toxicity, and cost must be considered.

preprint2021arXiv

Reentrant Correlated Insulators in Twisted Bilayer Graphene at 25T ($2π$ Flux)

Twisted bilayer graphene (TBG) is remarkable for its topological flat bands, which drive strongly-interacting physics at integer fillings, and its simple theoretical description facilitated by the Bistritzer-MacDonald Hamiltonian, a continuum model coupling two Dirac fermions. Due to the large moiré unit cell, TBG offers the unprecedented opportunity to observe reentrant Hofstadter phases in laboratory-strength magnetic fields near $25$T. This Letter is devoted to magic angle TBG at $2π$ flux where the magnetic translation group commutes. We use a newly developed gauge-invariant formalism to determine the exact single-particle band structure and topology. We find that the characteristic TBG flat bands reemerge at $2π$ flux, but, due to the magnetic field breaking $C_{2z} \mathcal{T}$, they split and acquire Chern number $\pm1$. We show that reentrant correlated insulating states appear at $2π$ flux driven by the Coulomb interaction at integer fillings, and we predict the characteristic Landau fans from their excitation spectrum. We conjecture that superconductivity can also be re-entrant at $2π$ flux.

preprint2021arXiv

Topological Materials Discovery from Crystal Symmetry

Topological materials discovery has evolved at a rapid pace over the past 15 years following the identification of the first nonmagnetic topological insulators (TIs), topological crystalline insulators (TCIs), and 3D topological semimetals (TSMs). Most recently, through complete analyses of symmetry-allowed band structures - including the theory of Topological Quantum Chemistry (TQC) - researchers have determined crystal-symmetry-enhanced Wilson-loop and complete symmetry-based indicators for nonmagnetic topological phases, leading to the discovery of higher-order TCIs and TSMs. The recent application of TQC and related methods to high-throughput materials discovery has revealed that over half of all of the known stoichiometric, solid-state, nonmagnetic materials are topological at the Fermi level, over 85% of the known stoichiometric materials host energetically isolated topological bands, and that just under $2/3$ of the energetically isolated bands in known materials carry the stable topology of a TI or TCI. In this Review, we survey topological electronic materials discovery in nonmagnetic crystalline solids from the prediction of the first 2D and 3D TIs to the recently introduced methods that have facilitated large-scale searches for topological materials. We also discuss future venues for the identification and manipulation of solid-state topological phases, including charge-density-wave compounds, magnetic materials, and 2D few-layer devices.

preprint2021arXiv

Topology invisible to eigenvalues in obstructed atomic insulators

We consider the extent to which symmetry eigenvalues reveal the topological character of bands. Specifically, we compare distinct atomic limit phases (band representations) that share the same irreducible representations (irreps) at all points in the Brillouin zone and, therefore, appear equivalent in a classification based on eigenvalues. We derive examples where such &#34;irrep-equivalent&#34; phases can be distinguished by a quantized Berry phase or generalization thereof. These examples constitute a generalization of the Su-Schrieffer-Heeger chain: neither phase is topological, in the sense that localized Wannier functions exist, yet there is a topological obstruction between them. We refer to two phases as &#34;Berry obstructed atomic limits&#34; if they have the same irreps, but differ by Berry phases. This is a distinct notion from eigenvalue obstructed atomic limits, which differ in their symmetry irreps at some point in the Brillouin zone. We compute exhaustive lists of elementary band representations that are irrep-equivalent, in all space groups, with and without time-reversal symmetry and spin-orbit coupling, and use group theory to derive a set of necessary conditions for irrep-equivalence. Finally, we conjecture, and in some cases prove, that irrep-equivalent elementary band representations that are not equivalent can be distinguished by a topological invariant.

preprint2020arXiv

Application of the induction procedure and the Smith Decomposition in the calculation and topological classification of electronic band structures in the 230 space groups

The electronic properties in a solid depend on the specific form of the wave-functions that represent the electronic states in the Brillouin zone. Since the discovery of topological insulators, much attention has been paid to the restrictions that the symmetry imposes on the electronic band structures. In this work we apply two different approaches to characterize all types of bands in a solid by the analysis of the symmetry eigenvalues: the induction procedure and the Smith Decomposition method. The symmetry eigenvalues or irreps of any electronic band in a given space group can be expressed as the superposition of the eigenvalues of a relatively small number of building units (the \emph{basic} bands). These basic bands in all the space groups are obtained following a group-subgroup chain starting from P1. Once the whole set of basic bands are known in a space group, all other types of bands (trivial, strong topological or fragile topological) can be easily derived. In particular, we confirm previous calculations of the fragile root bands in all the space groups. Furthermore, we define an automorphism group of equivalences of the electronic bands which allows to define minimum subsets of, for instance, independent basic or fragile root bands.

preprint2020arXiv

Cascade of electronic transitions in magic-angle twisted bilayer graphene

Magic-angle twisted bilayer graphene (MATBG) exhibits a rich variety of electronic states, including correlated insulators, superconductors, and topological phases. Understanding the microscopic mechanisms responsible for these phases requires determining the interplay between electron-electron interactions and quantum degeneracy due to spin and valley degrees of freedom. Signatures of strong electron-electron correlations have been observed at partial fillings of the flat electronic bands in recent spectroscopic measurements. Transport experiments have shown changes in the Landau level degeneracy at fillings corresponding to an integer number of electrons per moiré unit cell. However, the interplay between interaction effects and the degeneracy of the system is currently unclear. Using high-resolution scanning tunneling microscopy (STM), we observed a cascade of transitions in the spectroscopic properties of MATBG as a function of electron filling. We find distinct changes in the chemical potential and a rearrangement of the low-energy excitations at each integer filling of the moiré flat bands. These spectroscopic features are a direct consequence of Coulomb interactions, which split the degenerate flat bands into Hubbard sub-bands. We find these interactions, the strength of which we can extract experimentally, to be surprisingly sensitive to the presence of a perpendicular magnetic field, which strongly modifies the spectroscopic transitions. The cascade of transitions we report here characterizes the correlated high-temperature parent phase from which various insulating and superconducting ground-state phases emerge at low temperatures in MATBG.

preprint2020arXiv

Crystalline symmetry-protected non-trivial topology in prototype compound BaAl$_4$

The BaAl$_4$ prototype crystal structure is the most populous of all structure types, and is the building block for a diverse set of sub-structures including the famous ThCr$_2$Si$_2$ family that hosts high-temperature superconductivity and numerous magnetic and strongly correlated electron systems. The MA$_4$ family of materials (M=Sr, Ba, Eu; A=Al, Ga, In) themselves present an intriguing set of ground states including charge and spin orders, but have largely been considered as uninteresting metals. Using electronic structure calculations, symmetry analysis and topological quantum chemistry techniques, we predict the exemplary compound BaAl$_4$ to harbor a three-dimensional Dirac spectrum with non-trivial topology and possible nodal lines crossing the Brillouin zone, wherein one pair of semi-Dirac points with linear dispersion along the $k_z$ direction and quadratic dispersion along the $k_x/k_y$ direction resides on the rotational axis with $C_{4v}$ point group symmetry. Electrical transport measurements reveal the presence of an extremely large, unsaturating positive magnetoresistance in BaAl$_4$ despite an uncompensated band structure, and quantum oscillations and angle-resolved photoemission spectroscopy measurements confirm the predicted multiband semimetal structure with pockets of Dirac holes and a Van Hove singularity (VHS) remarkably consistent with the theoretical prediction. We thus present BaAl$_4$ as a new topological semimetal, casting its prototype status into a new role as building block for a vast array of new topological materials.

preprint2020arXiv

Detection of Topological Materials with Machine Learning

Databases compiled using ab-initio and symmetry-based calculations now contain tens of thousands of topological insulators and topological semimetals. This makes the application of modern machine learning methods to topological materials possible. Using gradient boosted trees, we show how to construct a machine learning model which can predict the topology of a given existent material with an accuracy of 90%. Such predictions are orders of magnitude faster than actual ab-initio calculations. Through extensive testing of different models we determine which properties help detect topological materials. We identify the sources of our model&#39;s errors and we discuss approaches to overcome them.

preprint2020arXiv

Eta-Pairing in Hubbard Models: From Spectrum Generating Algebras to Quantum Many-Body Scars

We revisit the $η$-pairing states in Hubbard models and explore their connections to quantum many-body scars to discover a universal scars mechanism. $η$-pairing occurs due to an algebraic structure known as a Spectrum Generating Algebra (SGA), giving rise to equally spaced towers of eigenstates in the spectrum. We generalize the original $η$-pairing construction and show that several Hubbard-like models on arbitrary graphs exhibit SGAs, including ones with disorder and spin-orbit coupling. We further define a Restricted Spectrum Generating Algebra (RSGA) and give examples of perturbations to the Hubbard-like models that preserve an equally spaced tower of the original model as eigenstates. The states of the surviving tower exhibit a sub-thermal entanglement entropy, and we analytically obtain parameter regimes for which they lie in the bulk of the spectrum, showing that they are exact quantum many-body scars. The RSGA framework also explains the equally spaced towers of eigenstates in several well-known models of quantum scars, including the AKLT model.

preprint2020arXiv

Fragile Phases As Affine Monoids: Classification and Material Examples

Topological phases in electronic structures contain a new type of topology, called fragile, which can arise, for example, when an Elementary Band Representation (Atomic Limit Band) splits into a particular set of bands. We obtain, for the first time, a complete classification of the fragile topological phases which can be diagnosed by symmetry eigenvalues, to find an incredibly rich structure which far surpasses that of stable/strong topological states. We find and enumerate all hundreds of thousands of different fragile topological phases diagnosed by symmetry eigenvalues (available at http://www.cryst.ehu.es/html/doc/FragileRoots.pdf), and link the mathematical structure of these phases to that of Affine Monoids in mathematics. Furthermore, we predict and calculate, for the first time, (hundred of realistic) materials where fragile topological bands appear, and show-case the very best ones.

preprint2020arXiv

Fragile topology and flat-band superconductivity in the strong-coupling regime

In flat bands, superconductivity can lead to surprising transport effects. The superfluid &#34;mobility&#34;, in the form of the superfluid weight $D_s$, does not draw from the curvature of the band but has a purely band-geometric origin. In a mean-field description, a non-zero Chern number or fragile topology sets a lower bound for $D_s$, which, via the Berezinskii-Kosterlitz-Thouless mechanism, might explain the relatively high superconducting transition temperature measured in magic-angle twisted bilayer graphene (MATBG). For fragile topology, relevant for the bilayer system, the fate of this bound for finite temperature and beyond the mean-field approximation remained, however, unclear. Here, we use numerically exact Monte Carlo simulations to study an attractive Hubbard model in flat bands with topological properties akin to those of MATBG. We find a superconducting phase transition with a critical temperature that scales linearly with the interaction strength. We then investigate the robustness of the superconducting state to the addition of trivial bands that may or may not trivialize the fragile topology. Our results substantiate the validity of the topological bound beyond the mean-field regime and further stress the importance of fragile topology for flat-band superconductivity.

preprint2020arXiv

Large Classes of Quantum Scarred Hamiltonians from Matrix Product States

Motivated by the existence of exact many-body quantum scars in the AKLT chain, we explore the connection between Matrix Product State (MPS) wavefunctions and many-body quantum scarred Hamiltonians. We provide a method to systematically search for and construct parent Hamiltonians with towers of exact eigenstates composed of quasiparticles on top of an MPS wavefunction. These exact eigenstates have low entanglement in spite of being in the middle of the spectrum, thus violating the strong Eigenstate Thermalization Hypothesis (ETH). Using our approach, we recover the AKLT chain starting from the MPS of its ground state, and we derive the most general nearest-neighbor Hamiltonian that shares the AKLT quasiparticle tower of exact eigenstates. We further apply this formalism to other simple MPS wavefunctions, and derive new families of Hamiltonians that exhibit AKLT-like quantum scars. As a consequence, we also construct a scar-preserving deformation that connects the AKLT chain to the integrable spin-1 pure biquadratic model. Finally, we also derive other families of Hamiltonians that exhibit new types of exact quantum scars, including a $U(1)$-invariant perturbed Potts model.

preprint2020arXiv

Observation of backscattering induced by magnetism in a topological edge state

The boundary modes of topological insulators are protected by the symmetries of the non-trivial bulk electronic states. Unless these symmetries are broken, they can give rise to novel phenomena, such as the quantum spin Hall effect in one-dimensional (1D) topological edge states, where quasiparticle backscattering is suppressed by time-reversal symmetry (TRS). Here, we investigate the properties of the 1D topological edge state of bismuth in the absence of TRS, where backscattering is predicted to occur. Using spectroscopic imaging and spin-polarized measurements with a scanning tunneling microscope, we compared quasiparticle interference (QPI) occurring in the edge state of a pristine bismuth bilayer with that occurring in the edge state of a bilayer, which is terminated by ferromagnetic iron clusters that break TRS. Our experiments on the decorated bilayer edge reveal an additional QPI branch, which can be associated with spin-flip scattering across the Brioullin zone center between time-reversal band partners. The observed QPI characteristics exactly match with theoretical expectations for a topological edge state, having one Kramer&#39;s pair of bands. Together, our results provide further evidence for the non-trivial nature of bismuth, and, in particular, demonstrate backscattering inside a helical topological edge state induced by broken TRS through local magnetism.

preprint2020arXiv

Spin-Orbit-Induced Topological Flat Bands in Line and Split Graphs of Bipartite Lattices

Topological flat bands, such as the band in twisted bilayer graphene, are becoming a promising platform to study topics such as correlation physics, superconductivity, and transport. In this work, we introduce a generic approach to construct two-dimensional (2D) topological quasi-flat bands from line graphs and split graphs of bipartite lattices. A line graph or split graph of a bipartite lattice exhibits a set of flat bands and a set of dispersive bands. The flat band connects to the dispersive bands through a degenerate state at some momentum. We find that, with spin-orbit coupling (SOC), the flat band becomes quasi-flat and gapped from the dispersive bands. By studying a series of specific line graphs and split graphs of bipartite lattices, we find that (i) if the flat band (without SOC) has inversion or $C_2$ symmetry and is non-degenerate, then the resulting quasi-flat band must be topologically nontrivial, and (ii) if the flat band (without SOC) is degenerate, then there exists an SOC potential such that the resulting quasi-flat band is topologically nontrivial. This generic mechanism serves as a paradigm for finding topological quasi-flat bands in 2D crystalline materials and meta-materials.

preprint2020arXiv

The Landau Level of Fragile Topology

We study the Hofstadter butterfly and Landau levels of the twisted bilayer graphene (TBG). We show that the nontrivial fragile topology of the lowest two bands near the charge neutral point makes their Hofstadter butterfly generically connected with higher bands, closing the gap between the first and second conduction (valence) bands at a certain magnetic flux per unit cell. We also develop a momentum space method for calculating the TBG Hofstadter butterfly, from which we identify three phases where the Hofstadter butterflies of the lowest two bands and the higher bands are connected in different ways. We show this leads to a crossing between the $ν=4$ Landau fan from the charge neutral point and the zero field band gap at one flux per Moiré unit cell, which corresponds to a magnetic field $25θ^2$T (twist angle $θ$ in degrees). This provides an experimentally testable feature of the fragile topology. In general, we expect it to be a generic feature that the Hofstadter butterfly of topological bands are connected with the Hofstadter spectra of other bands. We further show the TBG band theory with Zeeman splitting being the most sizable splitting could result in Landau fans at the charge neutral point and half fillings near the magic angle, and we predict their variations under an in-plane magnetic field.

preprint2019arXiv

Fermionic tensor networks for higher order topological insulators from charge pumping

We apply the charge pumping argument to fermionic tensor network representations of d-dimensional topological insulators (TIs) to obtain tensor network states for (d+1)-dimensional TIs. We exemplify the method by constructing a two-dimensional projected entangled pair states (PEPS)for a Chern insulator starting from a matrix product state (MPS) in d=1 describing pumping in the Su-Schrieffer-Heeger (SSH) model. In extending the argument to second-order TIs, we build a three-dimensional TNS for a chiral hinge TI from a PEPS in d=2 for the obstructed atomic insulator (OAI) of the quadrupole model. The (d+1)-dimensional TNSs obtained in this way have a constant bond dimension inherited from the d-dimensional TNSs in all but one spatial direction, making them candidates for numerical applications. From the d-dimensional models, we identify gapped next-nearest neighbour Hamiltonians interpolating between the trivial and OAI phases of the fully dimerized SSH and quadrupole models, whose ground states are given by an MPS and a PEPS with a constant bond dimension equal to 2, respectively.

preprint2019arXiv

Glide-resolved photoemission spectroscopy: Measuring topological invariants in nonsymmorphic space groups

The two classes of 3D, time-reversal-invariant insulators are known to subdivide into four classes in the presence of glide symmetry. Here, we extend this classification of insulators to include glide-symmetric Weyl metals, and find a finer $Z_4\oplus Z$ classification. We further elucidate the smoking-gun experimental signature of each class in the photoemission spectroscopy of surface states. Measuring the $Z_4$ topological invariant by photoemission relies on identifying the glide representation of the initial Bloch state before photo-excitation - we show how this is accomplished with relativistic selection rules, combined with standard spectroscopic techniques to resolve both momentum and spin. Our method relies on a novel spin-momentum locking that is characteristic of all glide-symmetric solids (inclusive of insulators and metals in trivial and topological categories). As an orthogonal application, given a glide-symmetric solid with an ideally symmetric surface, we may utilize this spin-momentum locking to generate a source of FULLY spin-polarized photoelectrons, which have diverse applications in solid-state spectroscopy. Our ab-initio calculations predict Ba$_2$Pb, stressed Na$_3$Bi, and KHgSb to realize all three, nontrivial insulating phases in the $Z_4$ classification.

preprint2019arXiv

Real Space Invariants: Twisted Bulk-Boundary Correspondence of Fragile Topology

In this paper, we propose a new type of bulk-boundary correspondence as a generic approach to theoretically and experimentally detect fragile topological states. When the fragile phase can be written as a difference of a trivial atomic insulator and the so-called obstructed atomic insulator, the gap between the fragile phase and other bands must close under a specific novel twist of the boundary condition of the system. We explicitly work out all the twisted boundary conditions (TBC) that can detect all the 2D fragile phases implied by symmetry eigenvalues in all wallpaper groups. We develop the concept of real space invariants - local good quantum numbers in real space - which fully characterize the eigenvalue fragile phases. We show that the number of unavoidable level crossings under the twisted boundary condition is completely determined by the real space invariants. Possible realizations of the TBC of the fragile band in metamaterial systems are discussed.

preprint2019arXiv

Topology-Bounded Superfluid Weight In Twisted Bilayer Graphene

While regular flat bands are good for enhancing the density of states and hence the gap, they are detrimental to the superfluid weight. We show that the predicted nontrivial topology of the two lowest flat bands of twisted bilayer graphene plays an important role in the enhancement of the superfluid weight and hence of superconductivity. We derive the superfluid weight (phase stiffness) of the TBLG superconducting flat bands with a uniform pairing, and show that it can be expressed as an integral of the Fubini-Study metric of the flat bands. This mirrors results already obtained for nonzero Chern number bands even though the TBLG flat bands have zero Chern number. We further show the metric integral is lower bounded by the topological $C_{2z}T$ Wilson loop winding number of the TBLG flat bands, which renders the superfluid weight has a topological lower bound proportional to the pairing gap. In contrast, trivial flat bands have a zero superfluid weight. The superfluid weight is crucial in determining the BKT transition temperature of the superconductor. Based on the transition temperature measured in TBLG experiments, we estimate the topological contribution of the superfluid weight in TBLG.

preprint2018arXiv

Entanglement Entropy From Tensor Network States for Stabilizer Codes

In this paper, we present the construction of tensor network states (TNS) for some of the degenerate ground states of 3D stabilizer codes. We then use the TNS formalism to obtain the entanglement spectrum and entropy of these ground-states for some special cuts. In particular, we work out the examples of the 3D toric code, the X-cube model and the Haah code. The latter two models belong to the category of &#34;fracton&#34; models proposed recently, while the first one belongs to the conventional topological phases. We mention the cases for which the entanglement entropy and spectrum can be calculated exactly: for these, the constructed TNS is the singular value decomposition (SVD) of the ground states with respect to particular entanglement cuts. Apart from the area law, the entanglement entropies also have constant and linear corrections for the fracton models, while the entanglement entropies for the toric code models only have constant corrections. For the cuts we consider, the entanglement spectra of these three models are completely flat. We also conjecture that the negative linear correction to the area law is a signature of extensive ground state degeneracy. Moreover, the transfer matrices of these TNS can be constructed. We show that the transfer matrices are projectors whose eigenvalues are either 1 or 0. The number of nonzero eigenvalues is tightly related to the ground state degeneracy.