Topology-Bounded Superfluid Weight in Twisted Bilayer Graphene

Auto-resolving the atomic structure at van der Waals interfaces using a generative model

The high-resolution visualization of atomic structures is significant for understanding the relationship between the microscopic configurations and macroscopic properties of materials. However, a rapid, accurate, and robust approach to automatically resolve complex patterns in atomic-resolution microscopy remains difficult to implement. Here, we present a Trident strategy-enhanced disentangled representation learning method (a generative model), which utilizes a few unlabelled experimental images with abundant low-cost simulated images to generate a large corpus of annotated simulation data that closely resembles experimental results, producing a high-quality large-volume training dataset. A structural inference model is then trained via a residual neural network which can directly deduce the interlayer slip and rotation of diversified and complicated stacking patterns at van der Waals (vdW) interfaces with picometer-scale accuracy across various materials (e.g. MoS2, WS2, ReS2, ReSe2, and 1 T'-MoTe2) with different layer numbers (bilayer and trilayers), demonstrating robustness to defects, imaging quality, and surface contaminations. The framework can also identify pattern transition interfaces, quantify subtle motif variations, and discriminate moiré patterns that are difficult to distinguish in frequency domains. Finally, the high-throughput processing ability of our method provides insights into a vdW epitaxy mode where various thermodynamically favorable slip stackings can coexist.

Berry Curvature Dipole and Nonlinear Hall Effect in Type-II Semi-Dirac Systems

The Berry curvature dipole (BCD) and the resulting nonlinear Hall effect have been investigated in various time-reversal (TR) invariant but inversion-breaking materials, where the primary mechanisms are typically attributed to low-energy Dirac models with tilt, Fermi surface (FS) warping, or semi-Dirac dispersion with quadratic momentum dependence in one direction. This study proposes that a nonzero BCD arises in a special type-II semi-Dirac model formed by the merging of three conventional Dirac points in the absence of TR symmetry. The BCD in this model initially increases, then decreases, as the chemical potential varies or as the Dirac cones merge. This non-monotonic behavior is strongly linked to the evolution of the FS and the distribution of the BCD density across different parameter regimes. Detailed comparison with other models further reveals distinct characteristics of this model. Additionally, it is shown that a pair of TR counterparts can exhibit a net BCD, suggesting potential applications in TR-invariant materials. As a concrete example, a pronounced BCD is demonstrated in bilayer graphene, where type-II semi-Dirac dispersion is achieved by fine-tuning interlayer sliding. These findings provide insights into BCD behavior in semi-Dirac materials and establish a foundation for exploring the nonlinear Hall effect in related systems.

Quantum geometry in condensed matter

One of the most celebrated accomplishments of modern physics is the description of fundamental principles of nature in the language of geometry. As the motion of celestial bodies is governed by the geometry of spacetime, the motion of electrons in condensed matter can be characterized by the geometry of the Hilbert space of their wave functions. Such quantum geometry, comprising Berry curvature and the quantum metric, can thus exert profound influences on various properties of materials. The dipoles of both Berry curvature and the quantum metric produce nonlinear transport. The quantum metric plays an important role in flat-band superconductors by enhancing the transition temperature. The uniformly distributed momentum-space quantum geometry stabilizes the fractional Chern insulators and results in the fractional quantum anomalous Hall effect. Here we review in detail quantum geometry in condensed matter, paying close attention to its effects on nonlinear transport, superconductivity and topological properties. Possible future research directions in this field are also envisaged.

Superfluid weight cross-over and critical temperature enhancement in singular flat bands

Nonanalytic Bloch eigenstates at isolated band degeneracy points exhibit singular behavior in the quantum metric. Here, a description of superfluid weight for zero-energy flat bands in proximity to other high-energy bands is presented, where they together form a singular band gap system. When the singular band gap closes, the geometric and conventional contributions to the superfluid weight as a function of the superconducting gap exhibit different cross-over behaviors. The scaling behavior of superfluid weight with the band gap is studied in detail, and the effect on the Berezinskii-Kosterlitz-Thouless transition temperature is explored. It is found that tuning the singular band gap provides a unique mechanism for enhancing the supercurrent and critical temperature of two-dimensional superconductors.

Superfluid stiffness of magic-angle twisted bilayer graphene

The physics of superconductivity in magic-angle twisted bilayer graphene (MATBG) is a topic of keen interest in moiré systems research, and it may provide an insight into the pairing mechanism of other strongly correlated materials such as high-critical-temperature superconductors. Here we use d.c. transport and microwave circuit quantum electrodynamics to directly measure the superfluid stiffness of superconducting MATBG through its kinetic inductance. We find the superfluid stiffness to be much larger than expected from conventional Fermi liquid theory. Rather, it is comparable to theoretical predictions1 and recent experimental indications2 of quantum geometric effects that are dominant at the magic angle. The temperature dependence of the superfluid stiffness follows a power law, which contraindicates an isotropic Bardeen-Cooper-Schrieffer (BCS) model. Instead, the extracted power-law exponents indicate an anisotropic superconducting gap, whether interpreted in the Fermi liquid framework or by considering the quantum geometry of flat-band superconductivity. Moreover, a quadratic dependence of the superfluid stiffness on both d.c. and microwave current is observed, which is consistent with the Ginzburg-Landau theory. Taken together, our findings show that MATBG is an unconventional superconductor with an anisotropic gap and strongly suggest a connection between quantum geometry, superfluid stiffness and unconventional superconductivity in MATBG. The combined d.c.-microwave measurement platform used here is applicable to the investigation of other atomically thin superconductors.

Corner Charge Fluctuation as an Observable for Quantum Geometry and Entanglement in Two-Dimensional Insulators

Measuring bipartite fluctuations of a conserved charge, such as the particle number, is a powerful approach to understanding quantum systems. When the measured region has sharp corners, the bipartite fluctuation receives an additional contribution known to exhibit a universal angle dependence in 2D isotropic and uniform systems. Here we establish that, for generic lattice systems of interacting particles, the corner charge fluctuation is directly related to quantum geometry. We first provide a practical scheme to isolate the corner contribution on lattices and analytically prove that its angle dependence in the "small-angle limit" measures exclusively the integrated quantum metric. A model of a compact obstructed atomic insulator is introduced to illustrate this effect analytically, while numerical verification for various Chern insulator models further demonstrate the experimental relevance of the corner charge fluctuation in a finite-size quantum simulator as a probe of quantum geometry. Last but not least, for free fermions, we unveil an intimate connection between quantum geometry and quantum information through the lens of corner entanglement entropies.

Three-Dimensional Multiorbital Flat Band Models and Materials

Flat band (FB) systems are essential for uncovering exotic quantum phenomena associated with strong electron correlations. Here we present a systematic theoretical framework for constructing multiorbital FB models and identifying feasible material candidates. This framework integrates group theory and crystallography into a symmetry-adapted tight-binding model incorporating lattice, site, and orbital degrees of freedom. Using this approach, we unveil a novel three-dimensional (3D) multiorbital FB model in the face-centered-cubic lattice, distinct from well-known single-orbital Lieb and kagome models. Critically, we identify numerous high-quality binary materials with ultraclean 3D FBs near the Fermi level. Furthermore, we explore diverse orbital bases within this model and extend our analysis to other cubic lattices with different space groups, broadening the scope for realizing 3D multiorbital FB systems. Our findings provide a foundational platform for exploring correlated physics in multiorbital FB systems and guiding the discovery of new quantum materials.

Generation of Isolated Flat Bands with Tunable Numbers through Moiré Engineering

In contrast to the Dirac-cone materials in which the low-energy spectrum features a pseudospin-1/2 structure, Lieb and Dice lattices both host triply degenerate low-energy excitations. Here, we discuss moiré structures involving twisted bilayers of these lattices, which are shown to exhibit a tunable number of isolated flat bands near the Fermi level due to the bipartite nature of their structures. These flat bands remain isolated from the high-energy bands even in the presence of small higher-order terms and chiral-symmetry-breaking interlayer tunneling. At small twist angles, many isolated flat bands can be generated with a notable Berry curvature, which could provide a geometric contribution to the superfluid weight under a BCS-type pairing potential. Remarkably, the emergence of isolated flat bands is insensitive to the twist angle, so that fine-tuning of the twist angle in an experimental setup would not be required. Our study suggests a promising new and efficient avenue for exploring and engineering flat bands based on the twisted bilayer Lieb and Dice lattices.

Energy spectrum theory of incommensurate systems

Because of the lack of translational symmetry, calculating the energy spectrum of an incommensurate system has always been a theoretical challenge. Here, we propose a natural approach to generalize energy band theory to incommensurate systems without reliance on the commensurate approximation, thus providing a comprehensive energy spectrum theory of incommensurate systems. Except for a truncation-dependent weighting factor, the formulae of this theory are formally almost identical to that of Bloch electrons, making it particularly suitable for complex incommensurate structures. To illustrate the application of this theory, we give three typical examples: one-dimensional bichromatic and trichromatic incommensurate potential models, as well as a moiré quasicrystal. Our theory establishes a fundamental framework for understanding incommensurate systems.

Quantized Integrated Shift Effect in Multigap Topological Phases

We show that certain three-dimensional multigap topological insulators can host quantized integrated shift photoconductivities due to bulk invariants that are defined under reality conditions imposed by additional symmetries. We recast the quantization in terms of the integrated torsion tensor and the non-Abelian Berry connection constituting Chern-Simons forms. Physically, we recognize that the topological quantization emerges purely from virtual transitions contributing to the optical response. Our findings provide another quantized electromagnetic dc response due to the nontrivial band topology, beyond the quantum anomalous Hall effect of Chern insulators and quantized circular photogalvanic effect found in Weyl semimetals.

Molecular Pairing in Twisted Bilayer Graphene Superconductivity

We propose a theory for how the weak phonon-mediated interaction (J_{A}=1-4  meV) wins over the prohibitive Coulomb repulsion (U=30-60  meV) and leads to a superconductor in magic-angle twisted bilayer graphene (MATBG). We find the pairing mechanism akin to that in the A_{3}C_{60} family of molecular superconductors: Each AA stacking region of MATBG resembles a C_{60} molecule, in that optical phonons can dynamically lift the degeneracy of the moiré orbitals, in analogy to the dynamical Jahn-Teller effect. Such induced J_{A} has the form of an intervalley anti-Hund's coupling and is less suppressed than U by the Kondo screening near a Mott insulator. Additionally, we also considered an intraorbital Hund's coupling J_{H} that originates from the on-site repulsion of a carbon atom. Under a reasonable approximation of the realistic model, we prove that the renormalized local interaction between quasiparticles has a pairing (negative) channel in a doped correlated insulator at ν=±(2+δν), albeit the bare interaction is positive definite. The proof is nonperturbative and based on exact asymptotic behaviors of the vertex function imposed by Ward identities. Existence of an optimal U for superconductivity is predicted. In a large area of the parameter space of J_{A}, J_{H}, the ground state is found to have a nematic d-wave singlet pairing, which, however, can lead to a p-wave-like nodal structure due to the Berry's phase on Fermi surfaces (or Euler obstruction).

Emergent phases in graphene flat bands

Electronic correlations in two-dimensional materials play a crucial role in stabilising emergent phases of matter. The realisation of correlation-driven phenomena in graphene has remained a longstanding goal, primarily due to the absence of strong electron-electron interactions within its low-energy bands. In this context, magic-angle twisted bilayer graphene has recently emerged as a novel platform featuring correlated phases favoured by the low-energy flat bands of the underlying moiré superlattice. Notably, the observation of correlated insulators and superconductivity, and the interplay between these phases have garnered significant attention. A wealth of correlated phases with unprecedented tunability was discovered subsequently, including orbital ferromagnetism, Chern insulators, strange metallicity, density waves, and nematicity. However, a comprehensive understanding of these closely competing phases remains elusive. The ability to controllably twist and stack multiple graphene layers has enabled the creation of a whole new family of moiré superlattices with myriad properties. Here, we review the progress and development achieved so far, encompassing the rich phase diagrams offered by these graphene-based moiré systems. Additionally, we discuss multiple phases recently observed in non-moiré multilayer graphene systems. Finally, we outline future opportunities and challenges for the exploration of hidden phases in this new generation of moiré materials.

Enhancing shift current response via virtual multiband transitions

Materials exhibiting a significant shift current response could potentially outperform conventional solar cell materials. The myriad of factors governing shift-current response, however, poses significant challenges in finding such strong shift-current materials. Here we propose a general design principle that exploits inter-orbital mixing to excite virtual multiband transitions in materials with multiple flat bands to achieve an enhanced shift current response. We further relate this design principle to maximizing Wannier function spread as expressed through the formalism of quantum geometry. We demonstrate the viability of our design using a 1D stacked Rice-Mele model. Furthermore, we consider a concrete material realization - alternating angle twisted multilayer graphene (TMG) - a natural platform to experimentally realize such an effect. We identify a set of twist angles at which the shift current response is maximized via virtual transitions for each multilayer graphene and highlight the importance of TMG as a promising material to achieve an enhanced shift current response at terahertz frequencies. Our proposed mechanism also applies to other 2D systems and can serve as a guiding principle for designing multiband systems that exhibit an enhanced shift current response.

Geometric Stiffness in Interlayer Exciton Condensates

Recent experiments have confirmed the presence of interlayer excitons in the ground state of transition metal dichalcogenide bilayers. The interlayer excitons are expected to show remarkable transport properties when they undergo Bose condensation. In this Letter, we demonstrate that quantum geometry of Bloch wave functions plays an important role in the phase stiffness of the interlayer exciton condensate. Notably, we identify a geometric contribution that amplifies the stiffness, leading to the formation of a robust condensate with an increased Berezinskii-Kosterlitz-Thouless temperature. Our results have direct implications for the ongoing experimental efforts on interlayer excitons in materials that have nontrivial quantum geometry. We provide estimates for the geometric contribution in transition metal dichalcogenide bilayers through a realistic continuum model with gated Coulomb interaction, and find that the substantially increased stiffness may allow an interlayer exciton condensate to be realized at amenable experimental conditions.

Engineering Flat Bands in Twisted-Bilayer Graphene away from the Magic Angle with Chiral Optical Cavities

Twisted bilayer graphene (TBG) is a recently discovered two-dimensional superlattice structure which exhibits strongly correlated quantum many-body physics, including strange metallic behavior and unconventional superconductivity. Most of TBG exotic properties are connected to the emergence of a pair of isolated and topological flat electronic bands at the so-called magic angle, θ≈1.05°, which are nevertheless very fragile. In this work, we show that, by employing chiral optical cavities, the topological flat bands can be stabilized away from the magic angle in an interval of approximately 0.8°<θ<1.3°. As highlighted by a simplified theoretical model, time reversal symmetry breaking (TRSB), induced by the chiral nature of the cavity, plays a fundamental role in flattening the isolated bands and gapping out the rest of the spectrum. Additionally, TRSB suppresses the Berry curvature and induces a topological phase transition, with a gap closing at the Γ point, towards a band structure with two isolated flat bands with Chern number equal to 0. The efficiency of the cavity is discussed as a function of the twisting angle, the light-matter coupling and the optical cavity characteristic frequency. Our results demonstrate the possibility of engineering flat bands in TBG using optical devices, extending the onset of strongly correlated topological electronic phases in moiré superlattices to a wider range in the twisting angle.

Interacting topological quantum chemistry in 2D with many-body real space invariants

The topological phases of non-interacting fermions have been classified by their symmetries, culminating in a modern electronic band theory where wavefunction topology can be obtained from momentum space. Recently, Real Space Invariants (RSIs) have provided a spatially local description of the global momentum space indices. The present work generalizes this real space classification to interacting 2D states. We construct many-body local RSIs as the quantum numbers of a set of symmetry operators on open boundaries, but which are independent of the choice of boundary. Using the U(1) particle number, they yield many-body fragile topological indices, which we use to identify which single-particle fragile states are many-body topological or trivial at weak coupling. To this end, we construct an exactly solvable Hamiltonian with single-particle fragile topology that is adiabatically connected to a trivial state through strong coupling. We then define global many-body RSIs on periodic boundary conditions. They reduce to Chern numbers in the band theory limit, but also identify strongly correlated stable topological phases with no single-particle counterpart. Finally, we show that the many-body local RSIs appear as quantized coefficients of Wen-Zee terms in the topological quantum field theory describing the phase.

Spin-Triplet Superconductivity from Quantum-Geometry-Induced Ferromagnetic Fluctuation

We show that quantum geometry induces ferromagnetic fluctuation resulting in spin-triplet superconductivity. The criterion for ferromagnetic fluctuation is clarified by analyzing contributions from the effective mass and quantum geometry. When the non-Kramers band degeneracy is present near the Fermi surface, the Fubini-Study quantum metric strongly favors ferromagnetic fluctuation. Solving the linearized gap equation with the effective interaction obtained by the random phase approximation, we show that the spin-triplet superconductivity is mediated by quantum-geometry-induced ferromagnetic fluctuation.

Ginzburg-Landau Theory of Flat-Band Superconductors with Quantum Metric

Recent experimental studies unveiled highly unconventional phenomena in the superconducting twisted bilayer graphene (TBG) with ultraflat bands, which cannot be described by the conventional BCS theory. For example, given the small Fermi velocity of the flat bands, the superconducting coherence length predicted by BCS theory is more than 20 times shorter than the measured values. A new theory is needed to understand many of the unconventional properties of flat-band superconductors. In this Letter, we establish a Ginzburg-Landau (GL) theory from a microscopic flat-band Hamiltonian. The GL theory shows how the properties of the physical quantities such as the critical temperature, superconducting coherence length, upper critical field, and superfluid density are governed by the quantum metric of the Bloch states. One key conclusion is that the superconducting coherence length is not determined by the Fermi velocity but by the size of the optimally localized Wannier functions which are limited by the quantum metric. Applying the theory to TBG, we calculated the superconducting coherence length and the upper critical fields. The results match the experimental ones well without fine-tuning of parameters. The established GL theory provides a new and general theoretical framework for understanding flat-band superconductors with the quantum metric.

Symmetric Kondo Lattice States in Doped Strained Twisted Bilayer Graphene

We use the topological heavy fermion (THF) model and its Kondo lattice (KL) formulation to study the possibility of a symmetric Kondo (SK) state in twisted bilayer graphene. Via a large-N approximation, we find a SK state in the KL model at fillings ν=0,±1,±2 where a KL model can be constructed. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at ν=0,±1,±2,±3 we also find a similar symmetric state that is adiabatically connected to the symmetric Kondo state. We study the stability of the symmetric state by comparing its energy with the ordered (symmetry-breaking) states found in [H. Hu et al., Phys. Rev. Lett. 131, 026502 (2023).PRLTAO0031-900710.1103/PhysRevLett.131.026502, Z.-D. Song and B. A. Bernevig, Phys. Rev. Lett. 129, 047601 (2022).PRLTAO0031-900710.1103/PhysRevLett.129.047601] and find the ordered states to have lower energy at ν=0,±1,±2. However, moving away from integer fillings by doping the light bands, our mean-field calculations find the energy difference between the ordered state and the symmetric state to be reduced, which suggests the loss of ordering and a tendency toward Kondo screening. In order to include many-body effects beyond the mean-field approximation, we also performed dynamical mean-field theory calculations on the THF model in the nonordered phase. The spin susceptibility follows a Curie behavior at ν=0,±1,±2 down to ∼2  K where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At noninteger filling ν=±0.5,±0.8,±1.2 dynamical mean-field theory shows deviations from a 1/T susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Finally, we study the effect of a C_{3z}-rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks C_{3z} symmetry) can be stabilized at sufficiently large strain at ν=0,±1,±2. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at noninteger fillings or by applying strain.

Quantum Optics Measurement Scheme for Quantum Geometry and Topological Invariants

We show how a quantum optical measurement scheme based on heterodyne detection can be used to explore geometrical and topological properties of condensed matter systems. Considering a 2D material placed in a cavity with a coupling to the environment, we compute correlation functions of the photons exiting the cavity and relate them to the hybrid light-matter state within the cavity. Different polarizations of the intracavity field give access to all components of the quantum geometric tensor on contours in the Brillouin zone defined by the transition energy. Combining recent results based on the metric-curvature correspondence with the measured quantum metric allows us to characterize the topological phase of the material. Moreover, in systems where S_{z} is a good quantum number, the procedure also allows us to extract the spin Chern number. As an interesting application, we consider a minimal model for twisted bilayer graphene at the magic angle, and discuss the feasibility of extracting the Euler number.

Fubini-Study metric and topological properties of flat band electronic states: the case of an atomic chain withs - porbitals

The topological properties of the flat band states of a one-electron Hamiltonian that describes a chain of atoms withs - porbitals are explored. This model is mapped onto a Kitaev-Creutz type model, providing a useful framework to understand the topology through a nontrivial winding number and the geometry introduced by theFubini-Study (FS)metric. This metric allows us to distinguish between pure states of systems with the same topology and thus provides a suitable tool for obtaining the fingerprint of flat bands. Moreover, it provides an appealing geometrical picture for describing flat bands as it can be associated with a local conformal transformation over circles in a complex plane. In addition, the presented model allows us to relate the topology with the formation of compact localized states and pseudo-Bogoliubov modes. Also, the properties of the squared Hamiltonian are investigated in order to provide a better understanding of the localization properties and the spectrum. The presented model is equivalent to two coupled SSH chains under a change of basis.

Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund's Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory

We apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)PRLTAO0031-900710.1103/PhysRevLett.129.047601], to obtain its Kondo lattice limit. In this limit localized f electrons on a triangular lattice interact with topological conduction c electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0,±1,±2 are controlled by the heavy f electrons, while ν=±3 is at the border of a phase transition between two f-electron fillings. For ν=0,±1,±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the f moments in the full model and analytically prove the SU(4) Hund's rule for the ground state which maintains that two f electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite "bond-operators" method. We then compute the excitation spectrum of the f moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the f moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.

Pair Density Waves from Local Band Geometry

A band-projection formalism is developed for calculating the superfluid weight in two-dimensional multiorbital superconductors with an orbital-dependent pairing. It is discovered that, in this case, the band geometric superfluid stiffness tensor can be locally nonpositive definite in some regions of the Brillouin zone. When these regions are large enough or include nodal singularities, the total superfluid weight becomes nonpositive definite due to pairing fluctuations, resulting in the transition of a BCS state to a pair density wave (PDW). This geometric BCS-PDW transition is studied in the context of two-orbital superconductors, and proof of the existence of a geometric BCS-PDW transition in a generic topological flat band is established.

Superconductivity, Charge Density Wave, and Supersolidity in Flat Bands with a Tunable Quantum Metric

Predicting the fate of an interacting system in the limit where the electronic bandwidth is quenched is often highly nontrivial. The complex interplay between interactions and quantum fluctuations driven by the band geometry can drive competition between various ground states, such as charge density wave order and superconductivity. In this work, we study an electronic model of topologically trivial flat bands with a continuously tunable Fubini-Study metric in the presence of on-site attraction and nearest-neighbor repulsion, using numerically exact quantum Monte Carlo simulations. By varying the electron filling and the minimal spatial extent of the localized flat-band Wannier wave functions, we obtain a number of intertwined orders. These include a phase with coexisting charge density wave order and superconductivity, i.e., a supersolid. In spite of the nonperturbative nature of the problem, we identify an analytically tractable limit associated with a "small" spatial extent of the Wannier functions and derive a low-energy effective Hamiltonian that can well describe our numerical results. We also provide unambiguous evidence for the violation of any putative lower bound on the zero-temperature superfluid stiffness in geometrically nontrivial flat bands.

Suppression of Nonequilibrium Quasiparticle Transport in Flat-Band Superconductors

We study nonequilibrium transport through a superconducting flat-band lattice in a two-terminal setup with the Schwinger-Keldysh method. We find that quasiparticle transport is suppressed and coherent pair transport dominates. For superconducting leads, the ac supercurrent overcomes the dc current, which relies on multiple Andreev reflections. With normal-normal and normal-superconducting leads, the Andreev reflection and normal currents vanish. Flat-band superconductivity is, thus, promising not only for high critical temperatures, but also for suppressing unwanted quasiparticle processes.

Evidence for Dirac flat band superconductivity enabled by quantum geometry

In a flat band superconductor, the charge carriers' group velocity v_F is extremely slow. Superconductivity therein is particularly intriguing, being related to the long-standing mysteries of high-temperature superconductors¹ and heavy-fermion systems². Yet the emergence of superconductivity in flat bands would appear paradoxical, as a small v_F in the conventional Bardeen-Cooper-Schrieffer theory implies vanishing coherence length, superfluid stiffness and critical current. Here, using twisted bilayer graphene^3-7, we explore the profound effect of vanishingly small velocity in a superconducting Dirac flat band system^8-13. Using Schwinger-limited non-linear transport studies^14,15, we demonstrate an extremely slow normal state drift velocity v_n ≈ 1,000 m s^-1 for filling fraction ν between -1/2 and -3/4 of the moiré superlattice. In the superconducting state, the same velocity limit constitutes a new limiting mechanism for the critical current, analogous to a relativistic superfluid¹⁶. Importantly, our measurement of superfluid stiffness, which controls the superconductor's electrodynamic response, shows that it is not dominated by the kinetic energy but instead by the interaction-driven superconducting gap, consistent with recent theories on a quantum geometric contribution^8-12. We find evidence for small Cooper pairs, characteristic of the Bardeen-Cooper-Schrieffer to Bose-Einstein condensation crossover^17-19, with an unprecedented ratio of the superconducting transition temperature to the Fermi temperature exceeding unity and discuss how this arises for ultra-strong coupling superconductivity in ultra-flat Dirac bands.

Topological Flat Bands in Self-Complementary Plasmonic Metasurfaces

Photonics can be confined in real space with dispersion vanishing in the momentum space due to destructive interference. In this Letter, we report the experimental realization of flat bands with nontrivial topology in a self-complementary plasmonic metasurface. The band diagram and compact localized states are measured. In these nontrivial band gaps, we observe the topological edge states by near-field measurements. Furthermore, we propose a digitalized metasurface by loading controllable diodes with C_{3} symmetry in every unit cell. By pumping a digital signal into the metasurface, we investigate the interaction between incident waves and the dynamic metasurface. Experimental results indicate that compact localized states in the nontrivial flat band could enhance the wave-matter interactions to convert more incident waves to time-modulated harmonic photonics. Although our experiments are conducted in the microwave regime, extending the related concepts into the optical plasmonic systems is feasible. Our findings pave an avenue toward planar integrated photonic devices with nontrivial flat bands and exotic transmission phenomena.

Magic-Angle Twisted Bilayer Graphene as a Topological Heavy Fermion Problem

Magic-angle (θ=1.05°) 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 semimetallic 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 c electrons-we hence name our model as the topological heavy fermion model. Using this model, we obtain the U(4) and U(4)×U(4) symmetries of Refs. [1-5] 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 ±1 excitations [2,6,7], 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.

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 implementation 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, which represent the completion of the symmetry-indicated band topology of known nonmagnetic materials, have enabled the discovery of repeat-topological and supertopological materials, including rhombohedral bismuth and Bi₂Mg₃. We find that 52.65% of all materials are topological at E_F, roughly two-thirds of bands across all materials exhibit symmetry-indicated stable topology, and 87.99% of all materials contain at least one stable or fragile topological band.

Catalogue of flat-band stoichiometric materials

Topological electronic flattened bands near or at the Fermi level are a promising route towards unconventional superconductivity and correlated insulating states. However, the related experiments are mostly limited to engineered materials, such as moiré systems^1-3. Here we present a catalogue of the naturally occuring three-dimensional stoichiometric materials with flat bands around the Fermi level. We consider 55,206 materials from the Inorganic Crystal Structure Database catalogued using the Topological Quantum Chemistry website^4,5, which provides their structural parameters, space group, band structure, density of states and topological characterization. We combine several direct signatures and properties of band flatness with a high-throughput analysis of all crystal structures. In particular, we identify materials hosting line-graph or bipartite sublattices-in either two or three dimensions-that probably lead 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 the database to extract a curated list of 2,379 high-quality flat-band materials, from which we identify 345 promising candidates that potentially host flat bands with charge centres that are not strongly localized on the atomic sites. We showcase five representative materials and 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⁶.

Diagnosing topological phase transitions in 1D superconductors using Berry singularity markers

In this work I demonstrate how to characterize topological phase transitions in BDI symmetry class superconductors (SCs) in 1D, using the recently introduced approach of Berry singularity markers (BSMs). In particular, I apply the BSM method to the celebrated Kitaev chain model, as well as to a variant of it, which contains both nearest and next nearest neighbor equal spin pairings. Depending on the situation, I identify pairs of external fields which can detect the topological charges of the Berry singularities which are responsible for the various topological phase transitions. These pairs of fields consist of either a flux knob which controls the supercurrent flow through the SC, or, strain, combined with a field which can tune the chemical potential of the system. Employing the present BSM approach appears to be within experimental reach for topological nanowire hybrids.

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, obstructed bands are distinguished from trivial bands in the presence of interactions by the nonzero lower bound imposed on their superfluid weight.

Reproducibility in the fabrication and physics of moiré materials

Overlaying two atomic layers with a slight lattice mismatch or at a small rotation angle creates a moiré superlattice, which has properties that are markedly modified from (and at times entirely absent in) the 'parent' materials. Such moiré materials have progressed the study and engineering of strongly correlated phenomena and topological systems in reduced dimensions. The fundamental understanding of the electronic phases, such as superconductivity, requires a precise control of the challenging fabrication process, involving the rotational alignment of two atomically thin layers with an angular precision below 0.1 degrees. Here we review the essential properties of moiré materials and discuss their fabrication and physics from a reproducibility perspective.

Superfluid density, Josephson relation and pairing fluctuations in a multi-component fermion superfluid

In this work, a Josephson relation is generalized to a multi-component fermion superfluid. Superfluid density is expressed through a two-particle Green function for pairing states. When the system has only one gapless collective excitation mode, the Josephson relation is simplified, which is given in terms of the superfluid order parameters and the trace of two-particle normal Green function. In addition, it is found that the matrix elements of two-particle Green function is directly related to the matrix elements of the pairing fluctuations of superfluid order parameters. Furthermore, in the presence of inversion symmetry, the superfluid density is given in terms of the pairing fluctuation matrix. The results of the superfluid density in Haldane model show that the generalized Josephson relation can be also applied to a multi-band fermion superfluid in lattice.

Competing Zero-Field Chern Insulators in Superconducting Twisted Bilayer Graphene

The discovery of magic angle twisted bilayer graphene has unveiled a rich variety of superconducting, magnetic, and topologically nontrivial phases. Here, we show that the zero-field states at odd integer filling factors in h-BN nonaligned devices are consistent with symmetry broken Chern insulators, as is evidenced by the observation of the anomalous Hall effect near moiré cell filling factor ν=+1. The corresponding Chern insulator has a Chern number C=±1 and a relatively high Curie temperature of T_{c}≈4.5  K. In a perpendicular magnetic field above B>0.5  T we observe a transition of the ν=+1 Chern insulator from Chern number C=±1 to C=3, characterized by a quantized Hall plateau with R_{yx}=h/3e^{2}. These observations demonstrate that interaction-induced symmetry breaking leads to zero-field ground states that include almost degenerate and closely competing Chern insulators, and that states with larger Chern numbers couple most strongly to the B field. In addition, the device reveals strong superconducting phases with critical temperatures of up to T_{c}≈3.5  K. By providing the first demonstration of a system that allows gate-induced transitions between magnetic and superconducting phases, our observations mark a major milestone in the creation of a new generation of quantum electronics.

Quantum Geometry and Flat Band Bose-Einstein Condensation

We study the properties of a weakly interacting Bose-Einstein condensate (BEC) in a flat band lattice system by using the multiband Bogoliubov theory and discover fundamental connections to the underlying quantum geometry. In a flat band, the speed of sound and the quantum depletion of the condensate are dictated by the quantum geometry, and a finite quantum distance between the condensed and other states guarantees stability of the BEC. Our results reveal that a suitable quantum geometry allows one to reach the strong quantum correlation regime even with weak interactions.

Realization of topological Mott insulator in a twisted bilayer graphene lattice model

Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG system. Employing the large-scale density matrix renormalization group on the lattice model containing the projected Coulomb interactions only, we identify a first-order quantum phase transition between the insulating stripe phase and the quantum anomalous Hall state with the Chern number of ±1. Our results not only shed light on the mechanism of the quantum anomalous Hall state discovered at three-quarters filling, but also provide an example of the topological Mott insulator, i.e., the quantum anomalous Hall state in the strong coupling limit.

Optical spectral weight, phase stiffness, and T c bounds for trivial and topological flat band superconductors

We present exact results that give insight into how interactions lead to transport and superconductivity in a flat band where the electrons have no kinetic energy. We obtain bounds for the optical spectral weight for flat-band superconductors that lead to upper bounds for the superfluid stiffness and the two-dimensional (2D) [Formula: see text] We focus on on-site attraction [Formula: see text] on the Lieb lattice with trivial flat bands and on the π-flux model with topological flat bands. For trivial flat bands, the low-energy optical spectral weight [Formula: see text] with [Formula: see text] , where n is the flat-band density and Ω is the Marzari-Vanderbilt spread of the Wannier functions (WFs). We also obtain a lower bound involving the quantum metric. For topological flat bands, with an obstruction to localized WFs respecting all symmetries, we again obtain an upper bound for [Formula: see text] linear in [Formula: see text] We discuss the insights obtained from our bounds by comparing them with mean-field and quantum Monte Carlo results.

Coulomb interaction, phonons, and superconductivity in twisted bilayer graphene

The polarizability of twisted bilayer graphene, due to the combined effect of electron-hole pairs, plasmons, and acoustic phonons, is analyzed. The screened Coulomb interaction allows for the formation of Cooper pairs and superconductivity in a significant range of twist angles and fillings. The tendency toward superconductivity is enhanced by the coupling between longitudinal phonons and electron-hole pairs. Scattering processes involving large momentum transfers, Umklapp processes, play a crucial role in the formation of Cooper pairs. The magnitude of the superconducting gap changes among the different pockets of the Fermi surface.

Superfluid density and collective modes of fermion superfluid in dice lattice

The superfluid properties of attractive Hubbard model in dice lattice are investigated. It is found that three superfluid order parameters increase as the interaction increases. When the filling factor falls into the flat band, due to the infinite large density of states, the resultant superfluid order parameters are proportional to interaction strength, which is in striking contrast with the exponentially small counterparts in usual superfluid (or superconductor). When the interaction is weak, and the filling factor is near the bottom of the lowest band (or the top of highest band), the superfluid density is determined by the effective mass of the lowest (or highest) single-particle band. When the interaction is strong and filling factor is small, the superfluid density is inversely proportional to interaction strength, which is related to effective mass of tightly bound pairs. In the strong interaction limit and finite filling, the asymptotic behaviors of superfluid density can be captured by a parabolic function of filling factor. Furthermore, when the filling is in flat band, the superfluid density shows a logarithmic singularity as the interaction approaches zero. In addition, there exist three undamped collective modes for strong interactions. The lowest excitation is gapless phonon, which is characterized by the total density oscillations. The two others are gapped Leggett modes, which correspond relative density fluctuations between sublattices. The collective modes are also reflected in the two-particle spectral functions by sharp peaks. Furthermore, it is found that the two-particle spectral functions satisfy an exact sum-rule, which is directly related to the filling factor (or density of particle). The sum-rule of the spectral functions may be useful to distinguish between the hole-doped and particle-doped superfluid (superconductor) in experiments.

Theory of optical responses in clean multi-band superconductors

Electromagnetic responses in superconductors provide valuable information on the pairing symmetry as well as physical quantities such as the superfluid density. However, at the superconducting gap energy scale, optical excitations of the Bogoliugov quasiparticles are forbidden in conventional Bardeen-Cooper-Schrieffer superconductors when momentum is conserved. Accordingly, far-infrared optical responses have been understood in the framework of a dirty-limit theory by Mattis and Bardeen for over 60 years. Here we show, by investigating the selection rules imposed by particle-hole symmetry and unitary symmetries, that intrinsic momentum-conserving optical excitations can occur in clean multi-band superconductors when one of the following three conditions is satisfied: (i) inversion symmetry breaking, (ii) symmetry protection of the Bogoliubov Fermi surfaces, or (iii) simply finite spin-orbit coupling with unbroken time reversal and inversion symmetries. This result indicates that clean-limit optical responses are common beyond the straightforward case of broken inversion symmetry. We apply our theory to optical responses in FeSe, a clean multi-band superconductor with inversion symmetry and significant spin-orbit coupling. This result paves the way for studying clean-limit superconductors through optical measurements.

Fragile Topology and Flat-Band Superconductivity in the Strong-Coupling Regime

In flat bands, superconductivity can lead to surprising transport effects. The superfluid "mobility", 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 nonzero 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 numerically use 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. Then, we 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.

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 Letter, we introduce a generic approach to construct two-dimensional (2D) topological quasiflat 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 quasiflat 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 nondegenerate, then the resulting quasiflat band must be topologically nontrivial, and (ii) if the flat band (without SOC) is degenerate, then there exists a SOC potential such that the resulting quasiflat band is topologically nontrivial. This generic mechanism serves as a paradigm for finding topological quasiflat bands in 2D crystalline materials and metamaterials.

Graphene bilayers with a twist

Near a magic twist angle, bilayer graphene transforms from a weakly correlated Fermi liquid to a strongly correlated two-dimensional electron system with properties that are extraordinarily sensitive to carrier density and to controllable environmental factors such as the proximity of nearby gates and twist-angle variation. Among other phenomena, magic-angle twisted bilayer graphene hosts superconductivity, interaction-induced insulating states, magnetism, electronic nematicity, linear-in-temperature low-temperature resistivity and quantized anomalous Hall states. We highlight some key research results in this field, point to important questions that remain open and comment on the place of magic-angle twisted bilayer graphene in the strongly correlated quantum matter world.

Current-Induced Reversal of Anomalous Hall Conductance in Twisted Bilayer Graphene

It is observed experimentally that the sign of the Hall resistance can be flipped by a dc electric current in the twisted bilayer graphene (TBG) at 3/4 filling of the fourfold degenerate conduction flat bands. The experiment implies a switching of the valley polarization (VP) and topology in TBG. Here we present a theory on the current-induced switching of VP and topology. The presence of current in the bulk causes the redistribution of electron occupation in bands near the Fermi energy, which then deforms and shifts the band dispersion due to the Coulomb interaction. Above a critical current, the original occupied and empty bands can be swapped, resulting in the switching of VP and topology.