BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Zhi-Da Song

Zhi-Da Song contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mes-hallcond-mat.str-elcond-mat.mtrl-scicond-mat.supr-con

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
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

11 published item(s)

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)].

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gapless Criterion for Crystals from Effective Axion Field

Gapless criteria that can efficiently determine whether a crystal is gapless or not are particularly useful for identifying topological semimetals. In this work, we propose a sufficient gapless criterion for three-dimensional non-interacting crystals, based on the simplified expressions for the bulk average value of the static axion field. The brief logic is that two different simplified expressions give the same value in an insulator, and thus the gapless phase can be detected by the mismatch of them. We demonstrate the effectiveness of the gapless criterion in the magnetic systems with space groups 26 and 13, where mirror, glide, and inversion symmetries provide the simplified expressions. In particular, the gapless criterion can identify gapless phases that are missed by the symmetry representation approach, as illustrated by space group 26. Our proposal serves as a guiding principle for future discovery of topological semimetals.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access