Search

Kagome materials manifest rich physical properties due to the emergence of abundant electronic phases. Here, we carry out a high-throughput first-principles study of the kagome 1:6:6 family MT$_6$Z$_6$ materials in space group 191, focusing on their phonon instability and electronic flat bands. Different MT$_6$Z$_6$ kagome candidates reveal a remarkable variety of kagome flat bands ranging from unfilled, partially filled, to fully filled. Notably, the Mn/Fe-166 compounds exhibit partially filled flat bands with a pronounced sharp peak in the density of states near the Fermi level, leading to magnetic orders that polarize the bands and stabilize the otherwise unstable phonon. When the flat bands are located away from the Fermi level, we find a large number of phonon instabilities, which can be classified into three types, based on the phonon dispersion and vibrational modes. Type-I instabilities involve the in-plane distortion of kagome nets, while type-II and type-III present out-of-plane distortion of trigonal M and Z atoms. We take MgNi$_6$Ge$_6$ and HfNi$_6$In$_6$ as examples to illustrate the possible CDW structures derived from the emergent type-I and type-II instabilities. The type-I instability in MgNi$_6$Ge$_6$ suggests a nematic phase transition, governed by the local twisting of kagome nets. The type-II instability in HfNi$_6$In$_6$ may result in a hexagonal-to-orthorhombic transition, offering insight into the formation of MT$_6$Z$_6$ in other space groups. Additionally, the predicted ScNb$_6$Sn$_6$ is analyzed as an example of the type-III instability. Our predictions suggest a vast kagome family with rich properties induced by the flat bands, possible CDW transitions, and their interplay with magnetism., Comment: 14 pages, 7 figures, with 1000 pages of additional supplemental materials

In this work, we propose a mechanism of inducing stable nodal superconductivity for multiple-flat-band systems. This mechanism is based on the degenerate flat band nature, so that not only intra-eigen-band pairing, but also inter-eigen-band pairing has a significant influence on the superconductivity properties. Based on the Bogoliubov de Gennes formalism of the heavy fermion model for the twisted bilayer graphene as an example, we show that although the nodal nematic d-wave pairing has higher energy around the nodal points compared to the chiral d-wave pairing in the momentum space, the inter-eigen-band pairing can lower the energy in other momenta away from the nodes, so that the nematic d-wave pairing is energetically favored for its overall condensation energy. This type of mechanism is particular to multiple-flat-band systems with no Fermi surfaces., Comment: 22 pages, 9 figures

Electronic ordering is prevalent in correlated systems, which commonly exhibit competing interactions. Here, we use x-ray diffraction to show that the charge density wave transition temperature of FeGe increases with pressure and evolves towards a $\sqrt{3}\times\sqrt{3}$ periodic lattice modulation, $\mathbf{q}$$^*$=$\left(\frac{1}{3}\ \frac{1}{3}\ \frac{1}{2}\right)$. In the pressure interval between 4$<$$p$$<$12 GPa both orders coexist and the spatial extent of the $\sqrt{3}\times\sqrt{3}$ order at high pressure becomes nearly long-range, $\sim$30 unit cells, while the correlation length of the 2$\times$2 phase remains shorter-ranged. The $\sqrt{3}\times\sqrt{3}$ phase is the ground state above 15 GPa, consistent with harmonic DFT calculations that predict a dimerization induced $\sqrt{3}\times\sqrt{3}$ order without phonon softening. The pressure dependence of the integrated intensities of $\mathbf{q}$$_\mathrm{CDW}=\left(\frac{1}{2}\ 0\ \frac{1}{2}\right)$ and $\mathbf{q}$$^*$ indicates a competition between the 2$\times$2 and $\sqrt{3}\times\sqrt{3}$ and demonstrates that the ground state of FeGe is characterized by a rich landscape of metastable/fragile phases. We discuss possible scenarios based on an order-disorder transformation and the formation of Friedel oscillations.

We conducted a high-throughput search for topological magnetic materials on 522 new, experimentally reported commensurate magnetic structures from MAGNDATA, doubling the number of available materials on the Topological Magnetic Materials database. This brings up to date the previous studies which had become incomplete due to the discovery of new materials. For each material, we performed first-principle electronic calculations and diagnosed the topology as a function of the Hubbard U parameter. Our high-throughput calculation led us to the prediction of 250 experimentally relevant topologically non-trivial materials, which represent 47.89% of the newly analyzed materials. We present five remarkable examples of these materials, each showcasing a different topological phase: Mn${}_2$AlB${}_2$ (BCSID 1.508), which exhibits a nodal line semimetal to topological insulator transition as a function of SOC, CaMnSi (BCSID 0.599), a narrow gap axion insulator, UAsS (BCSID 0.594) a 5f-orbital Weyl semimetal, CsMnF${}_4$ (BCSID 0.327), a material presenting a new type of quasi-symmetry protected closed nodal surface and FeCr${}_2$S${}_4$ (BCSID 0.613), a symmetry-enforced semimetal with double Weyls and spin-polarised surface states., Comment: 249 pages, Full topologies and band structures are provided on the Topological Magnetic Materials Database https://www.topologicalquantumchemistry.fr/magnetic/

The intrinsic frustrated nature of a kagome lattice is amenable to the realization of exotic phases of matter, such as quantum spin liquids or spin ices, and more recently the multiple-$\mathrm{\textbf{q}}$ charge density waves (CDW) in the kagome metals. Despite intense efforts to understand the mechanism driving the electronic modulations, its origin is still unknown and hindered by competing interactions and intertwined orders. Here, we identify a dimerization-driven 2D hexagonal charge-diffuse precursor in the antiferromagnetic kagome metal FeGe and demonstrate that the fraction of dimerized/undimerized states is the relevant order parameter of the multiple-$\mathrm{\textbf{q}}$ CDW of a continuous phase transition. The pretransitional charge fluctuations with propagation vector $\mathrm{\textbf{q}=\textbf{q}_M}$ at T$_{\mathrm{CDW}}$$<$T$<$T$^*$(125 K) are anisotropic, hence holding a quasi-long-range bond-orientational order. The broken translational symmetry emerges from the anisotropic diffuse precursor, akin to the Ising scenario of antiferromagnetic triangular lattices. The temperature and momentum dependence of the critical scattering show parallels to the stacked hexatic $\mathrm{B}$-phases reported in liquid crystals and transient states of CDWs and highlight the key role of the topological defect-mediated melting of the CDW in FeGe.

Magic-angle twisted bilayer graphene (MATBG) hosts a multitude of strongly correlated states at partial fillings of its flat bands. In a magnetic field, these flat bands further evolve into a unique Hofstadter spectrum renormalized by strong Coulomb interactions. Here, we study the interacting Hofstadter states spontaneously formed within the topological magnetic subbands of an ultraclean MATBG device, notably including symmetry-broken Chern insulator (SBCI) states and fractional quantum Hall (FQH) states. The observed SBCI states form a cascade with their Chern numbers mimicking the main sequence correlated Chern insulators. The FQH states in MATBG form in Jain sequence; however, they disappear at high magnetic field, distinct from conventional FQH states which strengthen with increasing magnetic field. We reveal a unique magnetic field-driven phase transition from composite fermion phases to a dissipative Fermi liquid. Our theoretical analysis of the magnetic subbands hosting FQH states predicts non uniform quantum geometric properties far from the lowest Landau level. This points towards a more natural interpretation of these FQH states as in-field fractional Chern insulators of the magnetic subbands.

The fractional Chern insulators (FCIs) observed in pentalayer rhombohedral graphene/hexagonal boron nitride superlattices have a unique origin contrary to theoretical expectations: their non-interacting band structure is gapless, unlike standard FCIs and the Landau level. Hartree-Fock (HF) calculations at filling $\nu=1$ yield a gapped ground state with Chern number 1 through band mixing, identifying a possible parent state. However, many-body calculations restricted to the occupied HF band predispose the system towards FCIs and are essentially uncontrolled. In this work, we use unbiased multi-band exact diagonalization (ED) to allow fluctuations into the gapless bands for two normal-ordering schemes. In the "charge neutrality" scheme, the weak moir\'e potential leads to theoretical proposals based on Wigner crystal-like states. However, we find that FCIs seen in 1-band ED calculations are destroyed by band mixing, becoming gapless as fluctuations are included. In the "average" scheme, the Coulomb interaction with the periodic valence charge background sets up a stronger moir\'e potential. On small systems, FCIs at $\nu=1/3$ are destroyed in multi-band calculations, while those at $\nu=2/3$ are initially strengthened. However we do not converge to a stable FCI at $\nu=2/3$ even on the largest accessible systems. These findings question prior results obtained within projection to a single HF band. They suggest that current models do not support FCIs with correlation length small enough to be converged in accessible, unbiased ED calculations, or do not support FCIs at all.

Using comprehensive exact diagonalization calculations on $\theta \approx 3.7 ^{\circ}$ twisted bilayer MoTe$_2$ ($t$MoTe$_2$), as well as idealized Landau level models also relevant for lower $\theta$, we extract general principles for engineering fractional topological insulators (FTIs) in realistic situations. First, in a Landau level setup at $\nu=1/3+1/3$, we investigate what features of the interaction destroy an FTI. For both pseudopotential interactions and realistic screened Coulomb interactions, we find that sufficient suppression of the short-range repulsion is needed for stabilizing an FTI. We then study $\theta \approx 3.7 ^{\circ}$ $t$MoTe$_2$ with realistic band-mixing and anisotropic non-local dielectric screening. Our finite-size calculations only find an FTI phase at $\nu=-4/3$ in the presence of a significant additional short-range attraction $g$ that acts to counter the Coulomb repulsion at short distances. We discuss how further finite-size drifts, dielectric engineering, Landau level character, and band-mixing effects may reduce the required value of $g$ closer towards the experimentally relevant conditions of $t$MoTe$_2$. Projective calculations into the $n=1$ Landau level, which resembles the second valence band of $\theta\simeq 2.1^\circ$ $t$MoTe$_2$, do not yield FTIs for any $g$, suggesting that FTIs at low-angle $t$MoTe$_2$ for $\nu=-8/3$ and $-10/3$ may be unlikely. While our study highlights the challenges, at least for the fillings considered, to obtaining an FTI with transport plateaus, even in large-angle $t$MoTe$_2$ where fractional Chern insulators are experimentally established, we also provide potential sample-engineering routes to improve the stability of FTI phases., Comment: 6+36 pages

The recent observation of the fractional quantum anomalous Hall effect in moir\'e fractional Chern insulators (FCI) provides opportunities for investigating zero magnetic field anyons. So far, both experimental and theoretical results suggest that filling > 1/3 FCI states in the first Chern band share features with those of the lowest Landau level (LL). To create the possibility of realizing non-Abelian anyons, one route is to engineer higher flat Chern bands that mimic higher LLs. Here, we investigate the interaction, topology, and ferromagnetism of the second moir\'e miniband in twisted MoTe2 bilayer (tMoTe2). Around filling factor v = -3, i.e., half-filling of the second miniband, we uncover spontaneous ferromagnetism and an incipient Chern insulator state. By measuring the anomalous Hall effect as a function of twist angle, we find that the Chern numbers (C) of the top two moir\'e flat bands have opposite sign (C = -+1) at twist angles above 3.1{\deg} but the same sign (C = -1) around 2.6{\deg}. This observation is consistent with the recently predicted twist-angle dependent band topology, resulting from the competition between moir\'e ferroelectricity and piezoelectricity. As we increase the magnetic field, only the small twist-angle device (2.6{\deg}) experiences a topological phase transition with an emergent C = -2 state. This is attributed to a Zeeman field-induced band crossing between opposite valleys, with the determined C = -1 for the top two bands. Our results lay a firm foundation for understanding the higher flat Chern bands, which is essential for the prediction or discovery of non-Abelian FCIs., Comment: 24 pages, 4 figures

Ingap states are commonly observed in semiconductors and are often well characterized by a hydrogenic model within the effective mass approximation. However, when impurities are strong, they significantly perturb all momentum eigenstates, leading to deep-level bound states that reveal the global properties of the unperturbed band structure. In this work, we discover that the topology of band wavefunctions can impose zeros in the impurity-projected Green's function within topological gaps. These zeros can be interpreted as spectral attractors, defining the energy at which ingap states are pinned in the presence of infinitely strong local impurities. Their pinning energy is found by minimizing the level repulsion of band eigenstates onto the ingap state. We refer to these states as ring states, marked by a mixed band character and a node at the impurity site, guaranteeing their orthogonality to the bare impurity eigenstates and a weak energy dependence on the impurity strength. We show that the inability to construct symmetric and exponentially localized Wannier functions ensures topological protection of ring states. Linking ring states together, the edge or surface modes can be recovered for any topologically protected phase. Therefore, ring states can also be viewed as building blocks of boundary modes, offering a framework to understand bulk-boundary correspondence., Comment: 12 pages, 5 figures + 20 pages supplementary material

Although the strongly interacting flat bands in twisted bilayer graphene (TBG) have been approached using the minimal Bistritzer-MacDonald (BM) Hamiltonian, there is mounting evidence that strain and lattice relaxation are essential in correctly determining the order of the correlated insulator groundstates. These effects can be incorporated in an enhanced continuum model by introducing additional terms computed from the relaxation profile. To develop an analytical and physical understanding of these effects, we include strain and relaxation in the topological heavy fermion (HF) model of TBG. We find that strain and relaxation are very well captured in first order perturbation theory by projection onto the fully symmetric HF Hilbert space, and remarkably do not alter the interacting terms in the periodic Anderson model. Their effects are fully incorporated in the single-particle HF Hamiltonian, and can be reproduced in a minimal model with only 4 symmetry-breaking terms. Our results demonstrate that the heavy fermion framework of TBG is an efficient and robust representation of the perturbations encountered in experiment.

We propose a general principle for the low-energy theory of narrow bands with concentrated Berry curvature and Fubini-Study metric in the form of a map to Anderson-"+" models composed of heavy fermions hybridizing and interacting with semi-metallic modes. This map resolves the obstruction preventing topological bands from being realized in a local Hamiltonian acting on the low-energy degrees of freedom. The concentrated quantum geometry is reproduced through band inversion with a dispersive semi-metal, leaving a nearly flat, trivial band which becomes the heavy fermion. This representation is natural when the narrow band is not energetically isolated on the scale of the interaction and an enlarged Hilbert space is inescapable, but also provides analytical insight into the projected-interaction limit. First exemplified in twisted bilayer graphene (TBG), we extend it to (1) the twisted checkerboard, which we find has a chiral symmetric stable anomaly that forbids a lattice realization at all energies, and (2) the Lieb lattice with gapless flat bands, where we show the heavy fermions can be obtained by minimizing a Euclidean instanton action to saturate its BPS bound. The heavy fermion approach is widely applicable and physically transparent: heavy electrons carry the strong correlations and dispersive electrons carry the topology. This simple picture unifies the dichotomous phenomena observed in TBG and points to connections between moir\'e and stoichiometric materials.

We study the interacting transport properties of twisted bilayer graphene (TBG) using the topological heavy-fermion (THF) model. In the THF model, TBG comprises localized, correlated $f$-electrons and itinerant, dispersive $c$-electrons. We focus on the Seebeck coefficient, which quantifies the voltage difference arising from a temperature gradient. We find that the TBG's Seebeck coefficient shows unconventional (strongly-interacting) traits: negative values with sawtooth oscillations at positive fillings, contrasting typical band-theory expectations. This behavior is naturally attributed to the presence of heavy (correlated, short-lived $f$-electrons) and light (dispersive, long-lived $c$-electrons) electronic bands. Their longer lifetime and stronger dispersion lead to a dominant transport contribution from the $c$-electrons. At positive integer fillings, the correlated TBG insulators feature $c$- ($f$-)electron bands on the electron (hole) doping side, leading to an overall negative Seebeck coefficient. Additionally, sawtooth oscillations occur around each integer filling due to gap openings. Our results highlight the essential importance of electron correlations in understanding the transport properties of TBG and, in particular, of the lifetime asymmetry between the two fermionic species (naturally captured by the THF model). Our findings are corroborated by new experiments in both twisted bilayer and trilayer graphene, and show the natural presence of strongly-correlated heavy and light carriers in the system., Comment: 5+106 pages, 3+36 figures, 6 tables. See also the accompanying experimental papers arXiv:2402.11749 and arXiv:2402.12296. The accompanying Ref. 252 will be posted soon

Ever since the initial experimental observation of correlated insulators and superconductivity in the flat Dirac bands of magic angle twisted bilayer graphene, a search for the microscopic description that explains its strong electronic interactions has begun. While the seemingly disagreeing electronic transport and scanning tunneling microscopy experiments suggest a dichotomy between local and extended electronic orbitals, definitive experimental evidence merging the two patterns together has been much sought after. Here, we report on the local photothermoelectric measurements in the flat electronic bands of twisted symmetric trilayer graphene (TSTG). We use a cryogenic scanning near-field optical microscope with an oscillating atomic force microscopy (AFM) tip irradiated by the infrared photons to create a nanoscopic hot spot in the planar samples, which generates a photocurrent that we probe globally. We observe a breakdown of the non-interacting Mott formalism at low temperatures (10K), signaling the importance of the electronic interactions. Our measurements reveal an overall negative offset of the Seebeck coefficient and significant peaks of the local photovoltage values at all positive integer fillings of the TSTG's moir\'e superlattice, further indicating a substantial deviation from the classical two-band semiconductor Seebeck response. We explain these observations using the interacting topological heavy-fermion model. In addition, our data reveal a spatial variation of the relative interaction strength dependent on the measured local twist angle (1.2{\deg} - 1.6{\deg}). Our findings provide experimental evidence of heavy fermion behaviour in the topological flat bands of moir\'e graphene and epitomize an avenue to apply local thermoelectric measurements to other strongly correlated materials in the disorder-free limit.

It has been recently postulated, that the strongly correlated flat bands of magicangle twisted bilayer graphene (MATBG) can host coexisting heavy and light carriers. While transport and spectroscopic measurements have shown hints of this behavior, a more direct experimental proof is still lacking. Here, we explore the thermoelectric response of MATBG through the photo-thermoelectric (PTE) effect in gate-defined MATBG pn-junctions. At low temperatures, we observe sign-preserving, fillingdependent oscillations of the Seebeck coefficient at non-zero integer fillings of the moir\'e lattice, which suggest the preponderance of one carrier type despite tuning the Fermi level from hole to electron doping of the correlated insulator. Furthermore, at higher temperatures, the thermoelectric response provides distinct evidence of the strong electron correlations in the unordered, normal state. We show that our observations are naturally accounted for by the interplay of light and long-lived and heavy and short-lived electron bands near the Fermi level at non-zero integer fillings. Our observations firmly establish the electron and hole asymmetry of the correlated gaps in MATBG, and shows excellent qualitative agreement with the recently developed topological heavy fermion model (THF).

We introduce the concept of "quantum geometric nesting'' (QGN) to characterize the idealized ordering tendencies of certain flat-band systems implicit in the geometric structure of the flat-band subspace. Perfect QGN implies the existence of an infinite class of local interactions that can be explicitly constructed and give rise to solvable ground states with various forms of possible fermion bi-linear order, including flavor ferromagnetism, density waves, and superconductivity. For the ideal Hamiltonians constructed in this way, we show that certain aspects of the low-energy spectrum can also be exactly computed including, in the superconducting case, the phase stiffness. Examples of perfect QGN include flat bands with certain symmetries (e.g. chiral or time-reversal), and non-symmetry-related cases exemplified with an engineered model for pair-density-wave. Extending this approach, we obtain exact superconducting ground states with nontrivial pairing symmetry.

We investigate in detail the $\nu=+1$ displacement-field-tuned interacting phase diagram of $L=3,4,5,6,7$ layer rhombohedral graphene aligned to hBN (R$L$G/hBN). Our calculations account for the 3D nature of the Coulomb interaction, the inequivalent stacking orientations $\xi=0,1$, the effects of the filled valence bands, and the choice of `interaction scheme' for specifying the many-body Hamiltonian. We show that the latter has a dramatic impact on the Hartree-Fock phase boundaries and the properties of the phases, including for pentalayers (R5G/hBN) with large displacement field $D$ where recent experiments observed a Chern insulator at $\nu=+1$ and fractional Chern insulators for $\nu<1$. In this large $D$ regime, the low-energy conduction bands are polarized away from the aligned hBN layer, and are hence well-described by the folded bands of moir\'eless rhombohedral graphene at the non-interacting level. Despite this, the filled valence bands develop moir\'e-periodic charge density variations which can generate an effective moir\'e potential, thereby explicitly breaking the approximate continuous translation symmetry in the conduction bands, and leading to contrasting electronic topology in the ground state for the two stacking arrangements. Within time-dependent Hartree-Fock theory, we further characterize the strength of the moir\'e pinning potential in the Chern insulator phase by computing the low-energy $\mathbf{q}=0$ collective mode spectrum, where we identify competing gapped pseudophonon and valley magnon excitations. Our results emphasize the importance of careful examination of both the single-particle and interaction model for a proper understanding of the correlated phases in R$L$G/hBN.

Flat band materials such as the kagome metals or moir\'e superlattice systems are of intense current interest. Flat bands can result from the electron motion on numerous (special) lattices and usually exhibit topological properties. Their reduced bandwidth proportionally enhances the effect of Coulomb interaction, even when the absolute magnitude of the latter is relatively small. Seemingly unrelated to these cases is the large family of strongly correlated electron systems, which includes the heavy fermion compounds, cuprate and pnictide superconductors. In addition to itinerant electrons from large, strongly overlapping orbitals, they frequently contain electrons from more localized orbitals, which are subject to a large Coulomb interaction. The question then arises as to what commonality in the physical properties and microscopic physics, if any, exists between the two broad categories of materials? A rapidly increasing body of strikingly similar phenomena across the different platforms -- from electronic localization-delocalization transitions to strange metal behavior and unconventional superconductivity -- suggests that similar underlying principles could be at play. Indeed, it has recently been suggested that flat band physics can be understood in terms of Kondo physics. Inversely, the concept of electronic topology from lattice symmetry, which is fundamental in flat band systems, is enriching the field of strongly correlated electron systems where correlation-driven topological phases are increasingly being investigated. Here we elucidate this connection, survey the new opportunities for cross-fertilization in understanding across the platforms, and assess the prospect for new insights that may be gained into both the correlation physics and its intersection with electronic topology., Comment: 35 pages, 6 figures, perspective paper

The experimental discovery of fractional Chern insulators (FCIs) in rhombohedral pentalayer graphene twisted on hexagonal boron nitride (hBN) has preceded theoretical prediction. Supported by large-scale first principles relaxation calculations at the experimental twist angle of $0.77^\circ$, we obtain an accurate continuum model of $n=3,4,5,6,7$ layer rhombohedral graphene-hBN moir\'e systems. Focusing on the pentalayer case, we analytically explain the robust $|C|=0,5$ Chern numbers seen in the low-energy single-particle bands and their flattening with displacement field, making use of a minimal two-flavor continuum Hamiltonian derived from the full model. We then predict nonzero valley Chern numbers at the $\nu = -4,0$ insulators observed in experiment. Our analysis makes clear the importance of displacement field and the moir\'e potential in producing localized "heavy fermion" charge density in the top valence band, in addition to the nearly free conduction band. Lastly, we study doubly aligned devices as additional platforms for moir\'e FCIs with higher Chern number bands., Comment: Second paper in the moir\'e FCI series

The electronic structure and interactions of kagome materials such as 1:1 (FeGe class) and 1:6:6 (MgFe$_6$Ge$_6$ class) are complicated and involve many orbitals and bands at the Fermi level. Current theoretical models treat the systems in an $s$-orbital kagome representation, unsuited and incorrect both quantitatively and qualitatively to the material realities. In this work, we lay the basis of a faithful framework of the electronic model for this large class of materials. We show that the complicated ``spaghetti" of electronic bands near the Fermi level can be decomposed into three groups of $d$-Fe orbitals coupled to specific Ge orbitals. Such decomposition allows for a clear analytical understanding (leading to different results than the simple $s$-orbital kagome models) of the flat bands in the system based on the $S$-matrix formalism of generalized bipartite lattices. Our three minimal Hamiltonians can reproduce the quasi-flat bands, van Hove singularities, topology, and Dirac points close to the Fermi level, which we prove by extensive ab initio studies. We also obtain the interacting Hamiltonian of $d$ orbitals in FeGe using the constraint random phase approximation (cRPA) method. We then use this as a fundamental ``LEGO"-like building block for a large family of 1:6:6 kagome materials, which can be obtained by doubling and perturbing the FeGe Hamiltonian. We applied the model to its kagome siblings FeSn and CoSn, and also MgFe$_6$Ge$_6$. Our work serves as the first complete framework for the study of the interacting phase diagram of kagome compounds., Comment: 5+50 pages, 3+16 figures, previously submitted. This is the second paper of a series on kagome materials. See also the first paper: arXiv:2305.15469

Recent experiments observed fractional Chern insulators (FCI) in twisted bilayer MoTe$_2$ at zero magnetic field, yet even the single-particle model of this material is controversial, leading to unreliable predictions of the experimental phase diagram as discussed in [Yu et al., 2023]. In this light, we revisit the single-particle model of twisted bilayer MoTe$_2$. Utilizing large-scale density functional theory, we calculate the band structure of twisted AA-stacked bilayer MoTe$_2$ at various twist angles relevant to experiment. We find that a band inversion occurs near $4.41^\circ$ between the second and third bands. Our ab initio band structure is in qualitative agreement with [Wang et al., 2023], but shows important differences in the remote bands and in the $\Gamma$ valley. We incorporate two higher harmonic terms into the continuum model to capture the highest 3 valence bands per valley. We confirm that the two highest valence bands per valley have opposite Chern numbers with $|C|=1$ for small angles, and also use our model to predict a variety of Chern states in the remote bands accessible by displacement field. Finally, we perform DFT calculations and build models for the AB stacking configuration. Our work serves as a foundation for accurate determination of the correlated phases in twisted bilayer MoTe$_2$.

Fractionally filled Chern bands with strong interactions may give rise to fractional Chern insulator (FCI) states, the zero-field analogue of the fractional quantum Hall effect. Recent experiments have demonstrated the existence of FCIs in twisted bilayer MoTe$_2$ without external magnetic fields -- most robust at $\nu=-2/3$ -- as well as Chern insulators (CIs) at $\nu=-1$. Although the appearance of both of these states is theoretically natural in an interacting topological system, experiments repeatedly observe nonmagnetic states (lacking FCIs) at $\nu=-1/3$ and $-4/3$, a puzzling result which has not been fully theoretically explained. In this work, we perform Hartree-Fock and exact diagonalization calculations to test whether the standard MoTe$_2$ moir\'e model with the (greatly varying) parameter values available in the literature can reproduce the non-magnetic states at $\nu=-1/3$ and $-4/3$ in unison with the FCI at $\nu=-2/3$ and CI state at $\nu = -1$. We focus on the experimentally relevant twist angles and, crucially, include remote bands. We find that the parameters proposed in [Wang et al. (2023)] can nearly capture the experimental phenomena at $\nu=-1/3,-2/3,-1,-4/3$ simultaneously, though the predicted ground states at $\nu=-1/3$ are still mostly fully-spin-polarized and a larger dielectric constant $\epsilon>10$ than is typical of hexagonal boron nitride (h-BN) substrate $\epsilon\sim 6$ is required. Our results show the importance of remote bands in identifying the competing magnetic orders and lay the groundwork for further study of the realistic phase diagram., Comment: 16+32 pages, 10+28 figures, 2+1 tables

In magic angle twisted bilayer graphene, transport, thermodynamic and spectroscopic experiments pinpoint at a competition between distinct low-energy states with and without electronic order. We use Dynamical Mean Field Theory (DMFT) on the topological heavy Fermion (THF) model of twisted bilayer graphene to investigate the emergence of electronic correlations and long-range order in the absence of strain. We contrast moment formation, Kondo screening and ordering on a temperature basis and explain the nature of emergent correlated states based on three central phenomena: (i) the formation of local spin and valley isospin moments around 100K, (ii) the ordering of the local isospin moments around 10K preempting Kondo screening, and (iii) a cascadic redistribution of charge between localized and delocalized electronic states upon doping. At integer fillings, we find that low energy spectral weight is depleted in the symmetric phase, while we find insulating states with gaps enhanced by exchange coupling in the zero-strain ordered phases. Doping away from integer filling results in distinct metallic states: a "bad metal" above the ordering temperature, where scattering off the disordered local moments suppresses electronic coherence, and a "good metal" in the ordered states with coherence of quasiparticles facilitated by isospin order. This finding reveals coherence from order as the microscopic mechanism behind the Pomeranchuk effect observed experimentally. Upon doping, there is a periodic charge reshuffling between localized and delocalized electronic orbitals leading to cascades of doping-induced Lifshitz transitions, local spectral weight redistributions and periodic variations of the electronic compressibility. Our findings provide a unified understanding of the most puzzling aspects of scanning tunneling spectroscopy, transport, and compressibility experiments., Comment: 22 pages, 7 figures

In this work, we propose that Majorana zero modes can be realized at the corners of the two-dimensional unconventional insulator. We demonstrate that 1T-PtSe2 is a symmetry indicator-free (SI-free) unconventional insulator, originating from orbital hybridization between Pt $d$ and Se $p_{x,y}$ states. The kind of SI-free unconventionality has no symmetry eigenvalue indication. Instead, it is diagnosed directly by the Wannier charge centers by using the one-dimensional Wilson loop method. The obstructed edge states exhibit strong anisotropy and large Rashba splitting. By introducing superconducting proximity and an external magnetic field, the Majorana corner modes can be obtained in the 1T-PtSe2 monolayer. In the end, we construct a two-Bernevig-Hughes-Zhang model with anisotropy to capture the Majorana physics.

In a series of recent reports, doped lead apatite (LK-99) has been proposed as a candidate ambient temperature and pressure superconductor. However, from both an experimental and theoretical perspective, these claims are largely unsubstantiated. To this end, our synthesis and subsequent analysis of an LK-99 sample reveals a multiphase material that does not exhibit high-temperature superconductivity. We study the structure of this phase with single-crystal X-ray diffraction (SXRD) and find a structure consistent with doped $\text{Pb}_{10}(\text{PO}_4)_6(\text{OH})_2$. However, the material is transparent which rules out a superconducting nature. From ab initio defect formation energy calculations, we find that the material likely hosts $\text{OH}^-$ anions, rather than divalent $\text{O}^{2-}$ anions, within the hexagonal channels and that Cu substitution is highly thermodynamically disfavored. Phonon spectra on the equilibrium structures reveal numerous unstable phonon modes. Together, these calculations suggest it is doubtful that Cu enters the structure in meaningful concentrations, despite initial attempts to model LK-99 in this way. However for the sake of completeness, we perform ab initio calculations of the topology, quantum geometry, and Wannier function localization in the Cu-dominated flat bands of four separate doped structures. In all cases, we find they are atomically localized by irreps, Wilson loops, and the Fubini-Study metric. It is unlikely that such bands can support strong superfluidity, and instead are susceptible to ferromagnetism (or out-of-plane antiferromagnetism) at low temperatures, which we find in ab initio studies. In sum, $\text{Pb}_{9}\text{Cu}(\text{PO}_4)_6(\text{OH})_2$ could more likely be a magnet, rather than an ambient temperature and pressure superconductor., Comment: 39 pages including appendices. Updated defect calculations and energy-dispersive X-ray spectroscopy data

Transitions between distinct obstructed atomic insulators (OAIs) protected by crystalline symmetries, where electrons form molecular orbitals centering away from the atom positions, must go through an intermediate metallic phase. In this work, we find that the intermediate metals will become a scale-invariant critical metal phase (CMP) under certain types of quenched disorder that respect the magnetic crystalline symmetries on average. We explicitly construct models respecting average $C_{2z}T$, $m$, and $C_{4z}T$ and show their scale-invariance under chemical potential disorder by the finite-size scaling method. Conventional theories, such as weak anti-localization and topological phase transition, cannot explain the underlying mechanism. A quantitative mapping between lattice and network models shows that the CMP can be understood through a semi-classical percolation problem. Ultimately, we systematically classify all the OAI transitions protected by (magnetic) groups $Pm$, $P2'$, $P4'$, and $P6'$ with and without spin-orbit coupling, most of which can support CMP.

This is our reply to "Comment on 'Nontrivial Quantum Geometry and the Strength of Electron-Phonon Coupling', arXiv:2305.02340, J. Yu, C. J. Ciccarino, R. Bianco, I. Errea, P. Narang, B. A. Bernevig" by Prof. Pickett, which focuses on the MgB$_2$ part of our work. We show that the entirety of the criticism in Prof. Pickett's comment pertaining to our work (arXiv:2305.02340) is invalid., Comment: 3+3 pages, 1 figure

Kagome Materials with flat bands exhibit wildly different physical properties depending on symmetry group, and electron number. For the case of ScV$_6$Sn$_6$ in space group 191, we investigate the existence of a charge density wave (CDW) at vector $\bar{K}=(\frac{1}{3},\frac{1}{3},\frac{1}{3})$ and its relationship with the phonon behavior. The experimental findings reveal a $\sim$95K CDW without nesting/peaks in the electron susceptibility at $\bar{K}$. Notably, ScV$_6$Sn$_6$ exhibits a collapsed phonon mode at $H=(\frac{1}{3},\frac{1}{3},\frac{1}{2})$ and an imaginary flat phonon band in the vicinity of $H$. The soft phonon is attributed to triangular Sn ($Sn^T$) mirror-even vibrations along the $z$-direction. We develop a simple force constant model to describe the entire soft phonon dispersion. By employing a new (Gaussian) approximation of the hopping parameter, we demonstrate the renormalization of the phonon frequency and the consequent collapse of the $H$ phonon. Additionally, we propose an effective model with two order parameters (OPs) to explain the appearance of the CDW at $\bar{K}$, which competes with the collapsed phonon at $H$. Through comparisons with experimental data, we show that the $H$ OP undergoes a second-order phase transition while exhibiting substantial fluctuations, ultimately inducing the first-order transition of the $\bar{K}$ OP. Furthermore, we extend our analysis to the similar compound YV$_6$Sn$_6$, which lacks a CDW phase, attributing this difference to the participation of the heavier Y atom in the out-of-plane phonon behavior. Our comprehensive study not only elucidates the CDW in ScV$_6$Sn$_6$ but also presents a significant advancement in modeling complex electronic systems, fostering collaborations between ab-initio simulations and analytical approaches., Comment: 5+108 pages, 3+42 figures, previously submitted. See also the related experimental study arXiv:2304.09173

The recently introduced topological heavy fermion model (THFM) provides a means for interpreting the low-energy electronic degrees of freedom of the magic angle twisted bilayer graphene as hybridization amidst highly dispersing topological conduction and weakly dispersing localized heavy fermions. In order to understand the Landau quantization of the ensuing electronic spectrum, a generalization of THFM to include the magnetic field B is desired, but currently missing. Here we provide a systematic derivation of the THFM in B and solve the resulting model to obtain the interacting Hofstadter spectra for single particle charged excitations. While naive minimal substitution within THFM fails to correctly account for the total number of magnetic subbands within the narrow band i.e. its total Chern number, our method -- based on projecting the light and heavy fermions onto the irreducible representations of the magnetic translation group -- reproduces the correct total Chern number. Analytical results presented here offer an intuitive understanding of the nature of the (strongly interacting) Hofstadter bands., Comment: 16+65 pages, 5+21 figures. The manuscript layout is restructured and the references are modified. The derivation of effective models is added to the supplementary notes

The coupling of electrons to phonons (electron-phonon coupling) is crucial for the existence of various phases of matter, in particular superconductivity and density waves. Here, we devise a theory that incorporates the quantum geometry of the electron bands into the electron-phonon coupling, demonstrating the crucial contributions of the Fubini-Study metric or its orbital selective version to the dimensionless electron-phonon coupling constant. We apply the theory to two materials, graphene and MgB$_2$ where the geometric contributions account for approximately 50\% and 90\% of the total electron-phonon coupling constant, respectively. The quantum geometric contributions in the two systems are further bounded from below by topological contributions. Our results suggest that the nontrivial electron band geometry/topology might favor superconductivity with relatively high critical temperature., Comment: submitted on 04/14/2023; close to the published version

We studied electronic band structure and topological property of a topological insulator thin film under a moir\'e superlattice potential to search for two-dimensional (2D) $\mathbb Z_2$ non-trivial isolated mini-bands. To model this system, we assume the Fermi energy inside the bulk band gap and thus consider an effective model Hamiltonian with only two surface states that are located at the top and bottom surfaces and strongly hybridized with each other. The moir\'e potential is generated by another layer of 2D insulating materials on top of topological insulator films. In this model, the lowest conduction (highest valence) mini-bands can be $\mathbb Z_2$ non-trivial when the minima (maxima) of the moir\'e potential approximately forms a hexagonal lattice with six-fold rotation symmetry. For the nontrivial conduction mini-band cases, the two lowest Kramers' pairs of conduction mini-bands both have nontrivial $\mathbb Z_2$ invariant in presence of inversion, while applying external gate voltages to break inversion leads to only the lowest Kramers' pair of mini-bands to be topologically non-trivial. The Coulomb interaction can drive the lowest conduction Kramers' mini-bands into the quantum anomalous Hall state when they are half-filled, which is further stabilized by breaking inversion symmetry. We propose the monolayer Sb$_{2}$ on top of Sb$_2$Te$_3$ thin films to realize our model based on results from the first principles calculations.

The long range electronic modulations recently discovered in the geometrically frustrated kagome lattice have opened new avenues to explore the effect of correlations in materials with topological electron flat bands. The observation of the lattice response to the emergent new phases of matter, a soft phonon mode, has remained elusive and the microscopic origin of charge density waves (CDWs) is still unknown. Here, we show, for the first time, a complete melting of the ScV$_ 6$Sn$_ 6$ (166) kagome lattice. The low energy phonon with propagation vector $\frac{1}{3} \frac{1}{3} \frac{1}{2}$ collapses at 98 K, without the emergence of long-range charge order, which sets in with a propagation vector $\frac{1}{3} \frac{1}{3} \frac{1}{3}$. The CDW is driven (but locks at a different vector) by the softening of an overdamped phonon flat plane at k$_z$=$\pi$. We observe broad phonon anomalies in momentum space, pointing to (1) the existence of approximately flat phonon bands which gain some dispersion due to electron renormalization, and (2) the effects of the momentum dependent electron-phonon interaction in the CDW formation. Ab initio and analytical calculations corroborate the experimental findings to indicate that the weak leading order phonon instability is located at the wave vector $\frac{1}{3} \frac{1}{3} \frac{1}{2}$ of a rather flat collapsed mode. We analytically compute the phonon frequency renormalization from high temperatures to the soft mode, and relate it to a peak in the orbital-resolved susceptibility, obtaining an excellent match with both ab initio and experimental results, and explaining the origin of the approximately flat phonon dispersion. Our data report the first example of the collapse of a softening of a flat phonon plane and promote the 166 compounds of the kagome family as primary candidates to explore correlated flat phonon-topological flat electron physics., Comment: 10 pages, 4 figures

We develop an analytic theory to describe the interaction between electrons and K-phonons and study its influence on superconductivity in the bare bands of twisted bilayer graphene (TBG). We find that, due to symmetry and the two-center approximation, only one optical K-phonon (~ 160meV) of graphene is responsible for inter-valley electron-phonon interaction. This phonon has recently been found in angular-resolved photo-emission spectroscopy to be responsible for replicas of the TBG flat bands. By projecting the interaction to the TBG flat bands, we perform the full symmetry analysis of phonon-mediated attractive interaction and pairing channels in the Chern basis, and show that several channels are guaranteed to have gapless order parameters. From the linearized gap equations, we find that the highest Tc pairing induced by this phonon is a singlet gapped s-wave inter-Chern-band order parameter, followed closely by a gapless nematic d-wave intra-Chern-band order parameter. We justify these results analytically, using the topological heavy fermion mapping of TBG which has allowed us to obtain an analytic form of phonon-mediated attractive interaction and to analytically solve the linearized and T=0 gap equations. For the intra-Chern-band channel, the nematic state with nodes is shown to be stabilized in the chiral flat band limit. While the flat band Coulomb interaction can be screened sufficiently enough - around Van-Hove singularities - to allow for electron-phonon based superconductivity, it is unlikely that this effect can be maintained in the lower density of states excitation bands around the correlated insulator states., Comment: 72 pages, 20 figures

The unusual properties of superconductivity in magic-angle twisted bilayer graphene (MATBG) have sparked enormous research interest. However, despite the dedication of intensive experimental efforts and the proposal of several possible pairing mechanisms, the origin of its superconductivity remains elusive. Here, using angle-resolved photoemission spectroscopy with micrometer spatial resolution, we discover replicas of the flat bands in superconducting MATBG unaligned with its hexagonal boron nitride (hBN) substrate, which are absent in non-superconducting MATBG aligned with the hBN substrate. Crucially, the replicas are evenly spaced in energy, separated by 150 +- 15 meV, signalling the strong coupling of electrons in MATBG to a bosonic mode of this energy. By comparing our observations to simulations, the formation of replicas is attributed to the presence of strong inter-valley electron-phonon coupling to a K-point phonon mode. In total, the observation of these replica flat bands and the corresponding phonon mode in MATBG could provide important information for understanding the origin and the unusual properties of its superconducting phase., Comment: 17 pages, 4 figures

The recently-introduced Topological Heavy Fermion (THF) model [Phys. Rev. Lett. 129, 047601] of twisted bilayer graphene (TBG) aims to reconcile the quantum-dot-like electronic structure of the latter observed by scanning tunneling microscopy, with its electron delocalization seen in transport measurements. The THF model achieves this by coupling localized (heavy) fermions with anomalous conduction electrons. Originally, the parameters of the THF model were obtained numerically from the Bistritzer-Macdonald (BM) model of TBG [Phys. Rev. Lett. 129, 047601]. In this work, we derive analytical expressions for the THF model parameters as a function of the twist angle, the ratio between the tunneling amplitudes at the $AA$ and $AB$ regions ($w_0 / w_1$), and the screening length of the interaction potential. By numerically computing the THF model parameters across an extensive experimentally-relevant parameter space, we show that the resulting approximations are remarkably good, i.e. within the 30% relative error for almost the entire parameter space. At the single-particle level, the THF model accurately captures the energy spectrum of the BM model over a large phase space of angles and tunneling amplitude ratios. When interactions are included, we also show that the THF description of TBG is good around the magic angle for realistic values of the tunneling amplitude ratios ($0.6 \leq w_0/w_1 \leq 1.0$), for which the hybridization between the localized and conduction fermions $\gamma$ is smaller than the onsite repulsion of the heavy fermions $U_1$ (i.e. $|\gamma| < U_1$)., Comment: 18+144 pages and 6+57 figures. Published version

Interactions among electrons create novel many-body quantum phases of matter with wavefunctions that often reflect electronic correlation effects, broken symmetries, and novel collective excitations. A wide range of quantum phases has been discovered in MATBG, including correlated insulating, unconventional superconducting, and magnetic topological phases. The lack of microscopic information, including precise knowledge of possible broken symmetries, has thus far hampered our understanding of MATBG's correlated phases. Here we use high-resolution scanning tunneling microscopy to directly probe the wavefunctions of the correlated phases in MATBG. The squares of the wavefunctions of gapped phases, including those of the correlated insulators, pseudogap, and superconducting phases, show distinct patterns of broken symmetry with a $\sqrt{3}$ x $\sqrt{3}$ super-periodicity on the graphene atomic lattice that has a complex spatial dependence on the moir\'e superlattice scale. We introduce a symmetry-based analysis to describe our measurements of the wavefunctions of MATBG's correlated phases with a set of complex-valued local order parameters. For the correlated insulators in MATBG, at fillings of $\nu$ = $\pm$ 2 electrons per moir\'e unit cell relative to charge neutrality, we compare the observed quantum textures to those expected for proposed theoretical ground states. In typical MATBG devices, the textures of correlated insulators' wavefunctions closely match those of the theoretically proposed IKS order, while in ultra-low-strain samples our data has local symmetries like those of a T-IVC phase. We also study the wavefunction of MATBG's superconducting state, revealing strong signatures of intervalley coherence that can only be distinguished from those of the insulator with our phase-sensitive measurements., Comment: 5 figures

Electronic structure calculations indicate that the Sr2FeSbO6 double perovskite has a flat-band set just above the Fermi level that includes contributions from ordinary sub-bands with weak kinetic electron hopping plus a flat sub-band that can be attributed to the lattice geometry and orbital interference. To place the Fermi energy in that flat band, electron doped samples with formulas Sr2-xLaxFeSbO6 (0 < x < 0.3) were synthesized and their magnetism and ambient temperature crystal structures determined by high-resolution synchrotron X-ray powder diffraction. All materials appear to display an antiferromagnetic-like maximum in the magnetic susceptibility, but the dominant spin coupling evolves from antiferromagnetic to ferromagnetic on electron doping. Which of the three sub-bands or combinations is responsible for the behavior has not been determined., Comment: 31 pages, 8 figures

We use the topological heavy fermion (THF) model and its Kondo Lattice (KL) formulation to study the symmetric Kondo state in twisted bilayer graphene. Via a large-N approximation, we find a symmetric Kondo (SK) state in KL mode at fillings $\nu=0,\pm 1,\pm 2$. 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 $\nu=0,\pm 1, \pm 2, \pm 3$, we also find a similar symmetric state. We study the stability of the symmetric state by comparing its energy with the ordered states and find the ordered states to have lower energy. However, moving away from integer fillings by doping holes to the light bands, we find the energy difference is reduced, which suggests the loss of ordering and a tendency towards Kondo screening. In order to include many-body effects beyond the mean-field approximation, we perform dynamical mean-field theory (DMFT) calculations on the THF model. We find the spin susceptibility follows a Curie behavior at $\nu=0, \pm 1,\pm 2$ down to $\sim 2\text{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 non-integer filling $\nu=\pm 0.5,\pm 0.8,\pm 1.2$ DMFT 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 $\nu=0,\pm 1, \pm 2$. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at non-integer fillings or by applying strain., Comment: 40 pages, 15 figures

We apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle $(\theta=1.05^{\circ})$ twisted bilayer graphene (MATBLG) (Phys. Rev. Lett. 129, 047601 (2022)), 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 $\nu=0, \pm 1, \pm 2$ are controlled by the heavy $f$-electrons, while $\nu = \pm 3$ is at the border of a phase transition between two $f$-electron fillings. For $\nu=0, \pm 1, \pm 2$, we then calculate the 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 twisted bilayer graphene., Comment: 120 pages, 8 figures

Recently, Ref. [1] reformulated magic-angle twisted bilayer graphene (MATBG) as a topological heavy fermion problem, and used this reformulation to provide a deeper understanding for the correlated phases at integer fillings. In this work, we generalize this heavy-fermion paradigm to magic-angle twisted symmetric trilayer graphene (MATSTG), and propose a low-energy $f-c-d$ model that reformulates MATSTG as heavy localized $f$ modes coupled to itinerant topological semimetalic $c$ modes and itinerant Dirac $d$ modes. Our $f-c-d$ model well reproduces the single-particle band structure of MATSTG at low energies for displacement field $\mathcal{E}\in[0,300]$meV. By performing Hartree-Fock calculations with the $f-c-d$ model for $\nu=0,-1,-2$ electrons per Moir\'e unit cell, we reproduce all the correlated ground states obtain from the previous numerical Hartree-Fock calculations with the Bistritzer-MacDonald-type (BM-type) model, and we find additional new correlated ground states at high displacement field. Based on the numerical results, we propose a simple rule for the ground states at high displacement fields by using the $f-c-d$ model, and provide analytical derivation for the rule at charge neutrality. We also provide analytical symmetry arguments for the (nearly-)degenerate energies of the high-$\mathcal{E}$ ground states at all the integer fillings of interest, and make experimental predictions of which charge-neutral states are stabilized in magnetic fields. Our $f-c-d$ model provides a new perspective for understanding the correlated phenomena in MATSTG, suggesting that the heavy fermion paradigm of Ref. [1] should be the generic underpinning of correlated physics in multilayer moire graphene structures., Comment: Close to published version

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 (in part) 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., Comment: 5 + 36 pages

Phonons play a crucial role in many properties of solid state systems, such as thermal and electrical conductivity, neutron scattering and associated effects or superconductivity. Hence, it is expected that topological phonons will also lead to rich and unconventional physics and the search of materials hosting topological phonons becomes a priority in the field. In electronic crystalline materials, a large part of the topological properties of Bloch states can be indicated by their symmetry eigenvalues in reciprocal space. This has been adapted to the high-throughput calculations of topological materials, and more than half of the stoichiometric materials on the databases are found to be topological insulators or semi-metals. Based on the existing phonon materials databases, here we have performed the first catalogue of topological phonon bands for more than ten thousand three-dimensional crystalline materials. Using topological quantum chemistry, we calculate the band representations, compatibility relations, and band topologies of each isolated set of phonon bands for the materials in the phonon databases. We have also calculated the real space invariants for all the topologically trivial bands and classified them as atomic and obstructed atomic bands. In particular, surface phonon modes (dispersion) are calculated on different cleavage planes for all the materials. Remarkably, we select more than one thousand "ideal" non-trivial phonon materials to fascinate the future experimental studies. All the data-sets obtained in the the high-throughput calculations are used to build a Topological Phonon Database., Comment: 8+535 pages, 187 figures, 21 tables. The Topological Phonon Database is available at https://www.topologicalquantumchemistry.com/topophonons or https://www.topologicalquantumchemistry.fr/topophonons

Models with Hilbert space fragmentation are characterized by (exponentially) many dynamically disconnected subspaces, not associated with conventional symmetries but captured by nontrivial Krylov subspaces. These subspaces usually exhibit a whole range of thermalization properties, from chaotic to integrable, to quantum many-body scars. However, so far, they have not been properly defined, nor can they be easily found, given a Hamiltonian. In this work, we consider Hamiltonians that have integer representations, a common feature of many (most) celebrated models in condensed matter. We show the equivalence of the integer characteristic polynomial factorization and the existence of Krylov subspaces generated from integer vectors. Considering the pair-hopping model, we illustrate how the factorization property can be used as a method to unveil Hilbert space fragmentation. We discuss the generalization over other rings of integers, for example those based on the cyclotomic field which are relevant when working in a given ($\ne 0, \pi$) momentum sector., Comment: 5 pages