Arxiv:2004.00638V3 [Cond-Mat.Str-El] 17 Mar 2021 Tron filling Νt = ±4 Per Moire´ Unit Cell
Total Page:16
File Type:pdf, Size:1020Kb
Charged Skyrmions and Topological Origin of Superconductivity in Magic Angle Graphene Eslam Khalaf,1 Shubhayu Chatterjee,2 Nick Bultinck,2, 3 Michael P. Zaletel,2 and Ashvin Vishwanath1, ∗ 1Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 2Department of Physics, University of California, Berkeley, CA 94720, USA 3Department of Physics, Ghent university, 9000 Gent, Belgium (Dated: March 18, 2021) Spontaneous symmetry breaking plays a pivotal role in stance takes place in quantum Hall ferromagnets, where topo- many areas of physics, engendering a variety of excitations logical textures of spin in the form of skyrmions acquire a from sound modes in solids to pions in nuclear physics. charge due to the Landau level topology and are found to be While these represent the linear response of the symme- the lowest energy charge excitations [8, 15, 17, 22, 31, 52, 65]. try breaking condensate, equally important are nonlin- On the other hand, finding situations where Cooper pair- ear configurations, solitons, which enjoy exceptional sta- ing occurs between charged topological textures, rather than bility when associated with a topological charge . A no- between electrons, is even harder to come by. For exam- table example is in the Skyrme model of nuclear forces, ple, in the aforementioned quantum Hall ferromagnets where where the neutron and proton are interpreted as topologi- charged topological textures have been experimentally estab- cal solitons of the pion condensate. We argue that similar lished, strong breaking of time-reversal symmetry makes su- models may describe a remarkable new material - magic perconductivity highly unlikely. In fact, obtaining robust su- angle graphene - obtained on twisting and stacking two perconductivity from topological textures typically requires sheets of graphene. When the relative angle between the simultaneously satisfying two conditions: (a) an unbroken sheets is near the magic angle of ∼ 1◦ degree - insulating time-reversal symmetry and (b) the existence of stable low en- behavior resulting from ordering of electrons is observed, ergy charge 2e topological textures. The latter can be achieved which gives way to superconductivity on changing the elec- if the fundamental defects have charge 2e or by pairing charge tron density. Here we propose a unifying description of e defects of the same electric charge. Here, we show that all both the order underlying the insulating state as well as these critera are satisfied in a simple model consisting of two the superconductivity. In our theory, while the symme- time-reversal-related quantum Hall (or flat Chern band) ferro- try breaking condensate leads to the ordered phase, topo- magnets with purely repulsive interactions, coupled via tun- logical solitons in the condensate - skyrmions - are shown neling. Such a model, we point out, captures the essential to be bosons that carry an electric charge 2e. Condensa- physics of magic angle graphene making it a promising can- tion of skyrmions leads to a superconductor whose pair- didate for superconductivity arising from the pairing of topo- ing strength, symmetry and other properties are inferred. logical textures, a mechanism that is fundamentally different More generally, we show how topological textures can mit- from the conventional electron-phonon mechanism. igate Coulomb repulsion to pair electrons and provide a new route to superconductivity. Our mechanism poten- tially applies to much wider class of systems but crucially Let us start with a brief review of magic angle twisted bi- invokes certain key ingredient such as inversion symmetry layer graphene (MATBG). Two sheets of graphene twisted present in magic angle graphene. We discuss how these relative to one another generate a moire´ pattern and corre- insights not only clarify why certain correlated moire ma- spondingly a reduced Brillouin zone in reciprocal space. Pre- terials do not superconduct, they also point to promising vious work demonstrated that the reconstruction of the elec- new platforms where robust superconductivity is antici- tronic bands by the moire´ lattice leads to minibands with ex- pated. tremely narrow bandwidth near charge neutrality at a magic angle ∼ 1o [5, 25, 29, 57], separated by a band gap from other bands. Recent experiments [8, 9, 30, 63] revealed dramatic I. INTRODUCTION new physics near the same magic angle. Band gaps arising from the moire´ potential are expected and observed at elec- Typically when charge is added to an insulator, electron or arXiv:2004.00638v3 [cond-mat.str-el] 17 Mar 2021 tron filling νT = ±4 per moire´ unit cell. In addition, insula- hole excitations are produced leading to charge conduction. tors at other integer filling including νT = ±2 [8, 30, 63] and Can solitons play the role of charge carriers? This unusual in some experiments also at νT = 0 [30, 55] and at certain scenario can be realized through electrically charged topologi- odd integer fillings have also been observed and are attributed cal textures which have been theoretically proposed in various to the effects of electron-electron interaction. Furthermore, contexts [1, 16, 18, 19, 23, 36, 56], although physical realiza- superconductivity has been repeatedly observed in twisted bi- tions are rare. In fact, the only experimentally established in- layer graphene, though its precise relation to the correlated insulating phase remains to be determined. While early ex- periments observed the superconductivity in the vicinity of the νT = −2 correlated insulator, subsequently a wider extent of ∗ [email protected] superconductivity has been observed [44, 55]. 2 II. QUANTUM HALL FERROMAGNETISM MODEL OF From valley/sublattice to Chern /pseudospin MAGIC ANGLE GRAPHENE (τ; σ) ! (γ; η) γ = (γx; γy; γz) = (σx; σyτz; σzτz) η = (η ; η ; η ) = (σ τ ; σ τ ; τ ) We begin by noting an important feature of the nearly flat x y z x x x y z bands near charge neutrality. Electrons in graphene carry both 0 Basis a spin (s ="; #) and a twofold valley degeneracy (τ = K; K ) 0 Valley τz = K=K which would suggest that filling each moire´ lattice site would Sublattice σ = A=B require four electrons. In reality, it takes twice as many elec- z Chern sector γ = σ τ = +=− trons to fill the nearly flat bands, since they consist of two con- z z z Pseudospin η = τ =" = # nected bands as shown in Figure 1. Conventionally, one dis- z z ps ps tinguishes the two bands by their kinetic energy but in view of the narrow bandwidth, other linear combinations may be Symmetries preferable. A different choice is the sublattice basis - which Symm (τ; σ) basis (γ; η) basis separates bands on the basis of their weight on the two sites T τxK γxηxK σ = A; B of the honeycomb lattice of monolayer graphene. C2 σxτx ηx iφτz iφηz This basis arises naturally in the ideal chiral limit of Ref. [24] UV (1) e e where the sublattice polarization is complete, but can also be defined at physical parameter values as explained in Ref. [7]. TABLE I: Definition of new variables - Chern Sector and Even for realistic parameters, the sublattice polarization is pseudospin, from the valley and sublattice degrees of substantial, so a reasonable approximation is to ignore sublat- freedom. Action of symmetries on the internal indices are tice off diagonal terms in the density operator [7]. This sublat- shown on the right, in both set of variables. Additionally, the tice polarized (SP) approximation will be assumed throughout independent pseudospin rotations in the Chern sectors are this paper. 1 generated by the ηP± where the projectors P± = 2 (1 ± γz) A remarkable feature of the sublattice basis is that individ- single out a specific Chern sector. ual bands in this basis carry Chern number. This allows for a mapping of MATBG to a pair of time-reversal related copies of quantum-Hall like systems. The presence of flavors then wavefunctions uγη;k(r), which can be conveniently encapsu- relates our problem to quantum Hall ferromagnetism, as ad- lated in the “form” factor expression for the density which can vocated in Ref. [7]. This mapping to a generalized spin-valley be written within the SP approximation as [7]: ferromagnet is broadly consistent with the observation of a X cascade of polarization transitions on varying electron den- n = F (k)eiγΦq (k)Ψy Ψ (2) sity [62, 66] and it acquires a particularly simple form in the q q γ;k γ;k+q k,γ=± vicinity of half-filling ν = ±2 when the spin degree of free- dom can be neglected. In this case, the low energy physics is Note that the form factor is diagonal in the Chern sectors dominated by four flat bands, each band being labelled by a with equal amplitude F and opposite phase Φ for opposite valley and sublattice index. Chern sectors. It follows that nq - and hence the interac- The valley-sublattice resolved flat bands are characterized tion - is symmetric under independent pseudo-spin rotations by Chern number C = ±1. Labelling sublattice (valley) by U+(2) × U−(2) in each C = ± sector. σ (τ), we have C = στ which is opposite for opposite valleys Remarkably, at half-filling the ground state of the interac- due time-reversal symmetry T and also opposite for opposite 1 P tion: HC = 2 q δnqVqδn−q can be determined exactly for sublattices within the same valley due to the combination of any purely repulsive interaction, Vq > 0 8q. In this case, two-fold rotation and time-reversal C2T . Thus, the flat bands HC is a sum of positive definite terms, so any state which is can be conveniently studied by introducing the pseudospin annihilated by each δnq jΩi = 0 is a ground state. This is T T spinors Ψ+ = (cKA; cK0B) and Ψ− = (cKB; cK0A) for satisfied by1 a manifold of C = 0 states, jΩi, obtained by the ± Chern sectors as shown in Fig.