## Solvable lattice models for metals with Z2 topological order

Brin Verheijden 1 , Yuhao Zhao 2 , Matthias Punk 3,4*

- 1 Physics Department, Ludwig-Maximilians-University Munich, D-80333 M¨ unchen 2 Department of Physics, ETH Zurich, 8093 Z¨ urich, Switzerland 3 Physics Department, Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Ludwig-Maximilians-University Munich, D-80333 M¨ unchen 4 Munich Center for Quantum Science and Technology (MCQST), Schellingstr. 4, D-80799 M¨ unchen, Germany

* matthias.punk@lmu.de

October 23, 2019

## Abstract

We present quantum dimer models in two dimensions which realize metallic ground states with Z2 topological order. Our models are generalizations of a dimer model introduced in [PNAS 112, 9552-9557 (2015)] to provide an e GLYPH&lt;11&gt; ective description of unconventional metallic states in holedoped Mott insulators. We construct exact ground state wave functions in a specific parameter regime and show that the ground state realizes a fractionalized Fermi liquid. Due to the presence of Z2 topological order the Luttinger count is modified and the volume enclosed by the Fermi surface is proportional to the density of doped holes away from half filling. We also comment on possible applications to magic-angle twisted bilayer graphene.

| Contents                                                 | Contents      |
|----------------------------------------------------------|---------------|
| 1 Introduction                                           | 2             |
| 2 Triangular lattice dimer model                         | 3             |
| 3 Exact ground state solution                            | 6             |
| 4 Perturbing away from the exactly solvable line: Z2-FL* | 7             |
| 5 Application to Twisted Bilayer Graphene                | 10            |
| 6 Discussion and Conclusions                             | 12            |
| A Kagome lattice dimer model                             | 12            |
| References 14                                            | References 14 |

## 1 Introduction

Landau's Fermi liquid theory is one of the cornerstones of condensed matter physics and is remarkably successful in describing conventional metallic phases of interacting electrons. Despite its wide success, various strongly correlated electron materials show unconventional metallic behavior which does not fit into the Fermi liquid framework. One prime example are the cuprate highTc superconductors, which exhibit a distinct non-Fermi liquid or 'strange metal' phase around optimal doping [1, 2], as well as a metallic 'pseudogap' phase at low hole-doping featuring Fermi-liquid like transport properties, but an anomalously low charge carrier concentration [3-5]. Interestingly, phases with a similar non-Fermi liquid phenomenology have been recently observed in magic-angle twisted bilayer graphene [6,7] .

The theoretical description of unconventional metallic phases in dimensions d GLYPH&lt;21&gt; 2 largely focuses on two broad classes of non-Fermi liquids. The first are metals without well-defined electronic quasiparticle excitations and appear e.g. in the phenomenological theory of marginal Fermi liquids [8], in the vicinity of metallic quantum critical points [9], or have recently been discussed in the context of SYKmodels [10,11]. A second class of models with non-Fermi liquid phenomenology is based on the concepts of topological order and fractionalization, where electronic degrees of freedom fractionalize into partons, each carrying some of the quantum numbers of an electron [12-17]. One striking consequence of topological order in such models is the possibility to modify Luttinger's theorem, which states that the volume enclosed by the Fermi surface is proportional to the density of electrons in the conduction band for ordinary Fermi liquids [18]. By contrast, so-called fractionalized Fermi liquids, which have been introduced originally in the context of heavy Fermion systems [19], provide a prime example for metallic phases with a modified Luttinger count and feature a reconstructed, small Fermi surface in the absence of broken translational symmetries.

In this work we present exctly solvable, two-dimensional lattice models that exhibit metallic ground states with Z2 topological order. These models are defined on the non-bipartite triangular and kagome lattices and are generalizations of a quantum dimer model introduced in Ref. [20]. The latter is defined on the square lattice and has been argued to capture some of the unusual electronic properties of the pseudogap phase in underdoped cuprates and is itself a generalization of the wellknown Rokhsar-Kivelson (RK) model [21]. Its Hilbert space is spanned by hard-core configurations of bosonic spin-singlet dimers, as well as fermionic spin-1 2 dimers carrying charge / q = + e , both living on nearest neighbor bonds. Fermionic dimers represent a hole in a bonding orbital between two neighboring lattice sites and can be viewed as bound states of a spinon and a holon in a doped resonating valence bond liquid. The density of fermionic dimers is equal to the density of doped holes away from half filling and the model reduces to the RK model at half filling. Subsequent numerical studies of the square lattice model computed single-electron spectral functions and revealed the presence of an anti-nodal pseudogap as well as Fermi arc-like features in the electron spectral function [22]. Furthermore, exact ground-state wave functions were constucted for a specific parameter regime in Ref. [23].

Due to the fact that the model in Ref. [20] is defined on the bipartite square lattice and only allows for nearest neighbor dimers, a fractionalized Fermi liquid ground state without broken symmetries only appears by fine-tuning to the special RK-point, where the bosonic dimers form a U (1) spin liquid. This is not a stable phase of matter, however, as the U (1) spin liquid is considered to be confining at large length scales [24, 25]. By contrast, the analogous quantum dimer models on the non-bipartite triangular and kagome lattices constructed in this work feature a stable fractionalized Fermi liquid (FL*) phase and don't require fine-tuning. At half-filling, i.e. at a vanishing density of

Figure 1: Example of a hard-core dimer configuration spanning the Hilbert space of our triangular lattice model. White ellipses represent bosonic spin-singlet dimers, black ellipses are fermionic dimers. The parity of the number of bosonic and fermionic dimers crossing the red dashed line does not change under arbitrary local dimer rearrangements, as expected for a state with Z 2 topological order.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00103v2-figure-08d00688-60c7-46dd-89fe-d56d304fd3d9.jpg)

fermionic dimers, these models exhibit an extended Z2 spin liquid phase [26] and realize a metallic Z2-FL* ground state without broken symmetries upon doping with holes [17].

Lastly we also note that a metallic phase with an unusually low charge carrier density has been observed in magic-angle twisted bilayer graphene (TBG) [6]. It appears upon hole-doping the Mott-like insulating phase at half filling of the lower Moir´ mini-band close to charge neutrality, in remarkable e analogy with the pseudogap phase in underdoped cuprates. In this work we argue that the triangular lattice quantum dimer model presented here could provide a toy model for the description of this unconventional metallic phase and we point out specific signatures which can be checked in future experiments. In particular we argue that the Fermi surface in the Z2-FL* phase of the triangular lattice dimer model consists of small hole pockets centered at the M points of the Brillouin zone (i.e. the Moir´ e mini-Brillouin zone in the case of twisted bilayer graphene). Angle resolved photoemission (ARPES) experiments with a su GLYPH&lt;14&gt; ciently high momentum resolution or quasi-particle interference in scanning tunnneling microscopy (STM) experiments on TBG should be able to test this prediction.

The rest of the paper is outlined as follows: in Sec. 2 we introduce the two-species quantum dimer model on the triangular lattice and construct exact ground states at a specific line in parameter space in Sec. 3. In Sec. 4 we consider perturbations away from the exactly solvable line and compare the perturbative results with numerical exact diagonalization data. Moreover, we discuss the emergence of a Z2-FL* ground state and compute the single electron spectral function. Finally, in Sec. 5 we discuss possible applications of this model to magic-angle twisted bilayer graphene. In the appendix we briefely present a construction of exact ground states for analogous dimer models on the kagome lattice.

## 2 Triangular lattice dimer model

As in Ref. [20] the Hilbert space of our model is spanned by hard-core coverings of the triangular lattice with two kinds of dimers living on nearest neighbor bonds. A bosonic dimer which represents two electrons in a spin-singlet configuration as in the usual Rokhsar-Kivelson model, as well as a fermionic spin-1 2 dimer carrying electric charge / + e with respect to a bosonic dimer background. The latter represents an electron in a bonding orbital delocalized between two neighboring lattice sites and can also be viewed as a tightly bound spinon-holon pair in a doped resonating valence bond liquid.

Figure 2: The three elementary plaquettes per lattice site, labelled P 1 2 3 ; ; . The possible configurations of a dimer, indicated by the three ellipses, are labelled by the lattice site and the direction GLYPH&lt;17&gt; 2 f 1 2 3 ; ; g of the unit vector e GLYPH&lt;17&gt; .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00103v2-figure-928ef04d-835e-4e27-a160-38dd56d1965d.jpg)

An example of a dimer configuration is shown in Fig. 1.

Dynamics on this Hilbert space is generated by a Hamiltonian including various local dimer resonance terms which can change the dimer configuration on each elementary plaquette consisting of an up- and a down-facing triangle as shown in Fig. 2, as well as potential energy terms for two parallel dimers on a plaquette. The terms which we consider here are sketched in Fig. 3. We define the canonical bosonic (fermionic) creation and annihilation operators D j y ;GLYPH&lt;17&gt; and Dj ;GLYPH&lt;17&gt; ( F j y ;GLYPH&lt;17&gt;;GLYPH&lt;27&gt; and Fj ;GLYPH&lt;17&gt;;GLYPH&lt;27&gt; ) of a bosonic (fermionic) dimer emanating from lattice site j in one of the three directions e GLYPH&lt;17&gt; with GLYPH&lt;17&gt; 2 f 1 2 3 , as depicted in Fig. 2. For fermionic dimers the index ; ; g GLYPH&lt;27&gt; labels the two spin components. In order to make a connection to the electronic Hilbert space of a Hubbard- or t-J model, the dimer operators can be expressed in terms of electron creation and annihilation operators c y j and cj GLYPH&lt;27&gt; as

GLYPH&lt;27&gt;

$$D _ { j, \eta } ^ { \dagger } \ \sim \ \frac { 1 } { \sqrt { 2 } } \left ( c _ { j \uparrow } ^ { \dagger } c _ { j + \mathrm e _ { \eta } \downarrow } ^ { \dagger } - c _ { j \downarrow } ^ { \dagger } c _ { j + \mathrm e _ { \eta } \uparrow } ^ { \dagger } \right ) \,,$$

$$F ^ { \dagger } _ { j, \eta, \sigma } \ \sim \ \frac { 1 } { \sqrt { 2 } } \left ( c ^ { \dagger } _ { j \sigma } + c ^ { \dagger } _ { j + \mathbf e _ { \eta } \sigma } \right ) \,,$$

up to a phase factor which depends on a gauge choice [20]. Using these dimer operators the Hamiltonian takes the form

$$\begin{array} { r c l } \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{array} \end{array} \begin{array} { r c l } \quad \\ \quad \\ H & = & - J \sum _ { j } D _ { j, 1 } ^ { \dagger } D _ { j + \epsilon _ { 2 }, 1 } ^ { \dagger } D _ { j, 2 } D _ { j + \epsilon _ { 1 }, 2 } + \dots \\ & & + V \sum _ { j } D _ { j, 1 } ^ { \dagger } D _ { j + \epsilon _ { 2 }, 1 } ^ { \dagger } D _ { j + \epsilon _ { 2 }, 1 } ^ { \dagger } D _ { j, 1 } + \dots \\ & & - t _ { 1 } \sum _ { j, \sigma } F _ { j, 1, \sigma } ^ { \dagger } D _ { j + \epsilon _ { 2 }, 1 } ^ { \dagger } F _ { j + \epsilon _ { 2 }, 1, \sigma } D _ { j, 1 } + \dots \\ & & + v _ { 1 } \sum _ { j, \sigma } F _ { j, 1, \sigma } ^ { \dagger } D _ { j + \epsilon _ { 2 }, 1 } ^ { \dagger } D _ { j + \epsilon _ { 2 }, 1 } F _ { j, 1, \sigma } + \dots \\ & & - t _ { 2 } \sum _ { j, \sigma } F _ { j, 1, \sigma } ^ { \dagger } D _ { j + \epsilon _ { 2 }, 1 } ^ { \dagger } D _ { j, 2 } F _ { j + \epsilon _ { 1 }, 2, \sigma } + \dots \,, \\ & & & \end{array}$$

Figure 3: Definition of plaquette terms appearing in our Hamiltonian in Eq. (3) and their respective transition amplitudes. Equivalent terms on the other two elementary plaquettes indicated in Fig. 2 as well as symmetry related terms are not shown.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00103v2-figure-e47dd678-28a1-4052-943e-4688b77601be.jpg)

where the dots indicate hermitean conjugate terms and analogous terms defined on the other two elementary plaquettes shown in Fig. 2, as well as symmetry related t 1, t 2 and v 1 terms, where bosonic and fermionic dimers are interchanged and or rotated in the initial or final state. As already mentioned, / the density of fermionic dimers is equal to the density of doped holes away from the half-filled electron band. At half-filling, i.e. without any fermionic dimers, only the terms in the first two lines of Eq. (3) remain and the Hamiltonian reduces to the Rokhsar-Kivelson model, which is exactly solvable at the so-called RK-point J = V . At this special point the Hamiltonian can be written as a sum of projectors on each elementary plaquette and the ground state is an equal weight superposition of all bosonic dimer coverings. It is important to note that the RK-point in the triangular lattice model is part of an extended Z2 spin liquid phase [26]. This is in contrast to the analogous model on the square lattice, where the RK-point is a singular point in the phase diagram exhibiting a U (1) spin liquid ground state.

Since the Hamiltonian (3) is local, the Hilbert space splits into di GLYPH&lt;11&gt; erent topological sectors depending on the lattice topology, in precise analogy to the RK-model on the triangular lattice. Using periodic boundary conditions e.g. in x-direction, the parity associated with the number of dimers crossing an arbitrary closed path around the lattice in x-direction, as indicated by the dashed line in Fig. 1, is conserved under arbitrary local operations on the dimer Hilbert space, which thus splits into two topological superselection sectors. This Z2 topological order is an important property of metallic phases realized in this model and allows for a modification of the conventional Luttinger count [27]. As we will see below, the model in Eq. (3) exhibits a metallic ground state where the Fermi volume is proportional to the density of fermionic dimers, i.e. proportional to the density of holes ( p ) away from half-filling. This is in contrast to the conventional Luttinger count for an ordinary Fermi-liquid metal, where the Fermi volume is proportional to the total density of holes measured from the fully filled band (i.e. 1 + p ).

In the following we are going to construct exact ground states for the full model including fermionic dimers. Note that we did not include purely fermionic plaquette resonance terms in Eq. (3), as this model is intended to describe systems at small hole doping slightly below half filling, where the density of fermionic dimers is small and such terms are not expected to play a prominent role.

## 3 Exact ground state solution

We proceed by constructing the exact ground state of a single fermionic dimer interacting with a background of bosonic dimers at a specific line in parameter space. A generalization to an arbitrary number of fermions (also with spin) is then straightforward. Our notation and construction follows closely Ref. [23].

Exact ground states can be constructed for the special choice of parameters J = V and v 1 = t 2 = GLYPH&lt;0&gt; t 1, in which case the Hamiltonian reduces to a sum of projectors

$$H = H _ { \text{RK} } + v _ { 1 } \sum _ { j, \alpha } P _ { j, \eta } \,,$$

where H RK is the standard Rokhsar-Kivelson Hamiltonian (i.e. the first two lines of Eq. (3)) at the special RK-point J = V [21] and we defined the three di GLYPH&lt;11&gt; erent plaquette projectors Pj ;GLYPH&lt;17&gt; = j GLYPH&lt;30&gt; j ;GLYPH&lt;17&gt; ih GLYPH&lt;30&gt; j ;GLYPH&lt;17&gt; j with GLYPH&lt;17&gt; 2 f 1 2 3 ; ; g for every lattice site j and plaquette GLYPH&lt;17&gt; as defined in Fig. 2. Furthermore,

$$| \phi _ { j, 1 } \rangle \ = \ | \nmid \rangle + | \nmid \rangle - | \nmid \rangle - | \nmid \rangle$$

$$| \phi _ { j, 2 } \rangle \ = \ | \bigtriangledown \Delta \rangle + | \bigtriangledown \Delta \rangle - | \bigtriangledown \rangle - | \bigtriangledown \rangle$$

$$| \phi _ { j, 3 } \rangle \ = \ | \not \subset \Pi \rangle + | \not \subset \Pi \rangle - | \not \subset \rangle - | \not \subset \rangle \.$$

Here white ellipses denote bosonic spin-singlet dimers and black ellipses represent fermionic dimers. Since we focus on a single fermionic dimer here, the spin index is suppressed. Also note that the Hamiltonian in Eq. (4) is positive semi-definite and thus all states with zero energy are automatically ground states.

We define basis states with one fermionic dimer fixed at a position specified by the lattice site j and the direction GLYPH&lt;17&gt; 2 f 1 2 3 ; ; g (see Fig. 2 for a definition)

$$| ( j, \eta ) \rangle = F _ { j, \eta } ^ { \dagger } D _ { j, \eta } | R K \rangle$$

where j RK i is the usual Rokhsar-Kivelson wave function, i.e. the equal weight superposition of all hard-core coverings with bosonic spin-singlet dimers in a given topological sector. It is important to note that these basis states are ground states of the RK Hamiltonian H RK at the exactly solvable RK point J = V per construction, i.e. H RK ( j j ; GLYPH&lt;17&gt; ) i = 0, because applying the RK Hamiltonian to a plaquette with a fermionic dimer gives zero. Moreover, applying the plaquette projectors P ';GLYPH&lt;11&gt; to these basis states gives

$$P _ { \ell, 1 } | ( j, \eta ) \rangle \ = \ \left [ \delta _ { \eta, 1 } ( \delta _ { j, \ell } + \delta _ { j, \ell + e _ { 2 } } ) - \delta _ { \eta, 2 } ( \delta _ { j, \ell } + \delta _ { j, \ell + e _ { 1 } } ) \right ] | \phi _ { \ell, 1 } \rangle \,,$$

$$P _ { \ell, 2 } | ( j, \eta ) \rangle \ = \ \left [ \delta _ { \eta, 2 } ( \delta _ { j, \ell } + \delta _ { j, \ell + e _ { 2 } } ) - \delta _ { \eta, 3 } ( \delta _ { j, \ell } + \delta _ { j, \ell + e _ { 2 } } ) \right ] | \phi _ { \ell, 2 } \rangle \,,$$

$$P _ { \ell, 3 } | ( j, \eta ) \rangle \ = \ \left [ \delta _ { \eta, 1 } ( \delta _ { j, \ell } + \delta _ { j, \ell + e _ { 3 } } ) - \delta _ { \eta, 3 } ( \delta _ { j, \ell } + \delta _ { j, \ell + e _ { 1 } } ) \right ] | \phi _ { \ell, 3 } \rangle \,.$$

At the exactly solvable line the ground state turns out to be highly degenerate and di GLYPH&lt;11&gt; erent ground states j p i can be parametrized by a lattice momentum p . We make the ansatz

$$| \mathfrak { p } \rangle = \sum _ { j, \eta } e ^ { i \mathfrak { p } \cdot \mathbf R _ { j } } c _ { \eta } ( \mathfrak { p } ) | ( j, \eta ) \rangle \,,$$

where R j is the lattice position of site j and determine the complex coe GLYPH&lt;14&gt; cients c GLYPH&lt;17&gt; ( p ) by demanding that H j p i = 0. Since H RK j p i = 0 per construction, this immediately leads to the three equations

$$\sum _ { \ell } P _ { \ell, \eta } | { \mathbf p } \rangle = 0$$

for each GLYPH&lt;17&gt; 2 f 1 2 3 . Using Eqs. (9), (10), (11) these equations reduce to ; ; g

$$c _ { 1 } ( \mathbf p ) \left ( 1 + e ^ { i p _ { 2 } } \right ) - c _ { 2 } ( \mathbf p ) \left ( 1 + e ^ { i p _ { 1 } } \right ) \ = \ 0$$

$$c _ { 2 } ( \mathbf p ) \left ( 1 + e ^ { i p _ { 3 } } \right ) - c _ { 3 } ( \mathbf p ) \left ( 1 + e ^ { i p _ { 2 } } \right ) \ = \ 0$$

$$c _ { 1 } ( \mathbf p ) \left ( 1 + e ^ { i p _ { 3 } } \right ) - c _ { 3 } ( \mathbf p ) \left ( 1 + e ^ { i p _ { 1 } } \right ) \ = \ 0$$

where we defined p GLYPH&lt;17&gt; GLYPH&lt;17&gt; p GLYPH&lt;1&gt; e GLYPH&lt;17&gt; . Together with the normalization condition

$$\langle \mathfrak { p } | \mathfrak { p } \rangle = N Q _ { c } [ ( j, \eta ) ] \sum _ { \eta } | c _ { \eta } ( \mathfrak { p } ) | ^ { 2 } = 1$$

with N the number of lattice sites and Qc [( j ; GLYPH&lt;17&gt; )] = h ( j ; GLYPH&lt;17&gt; ) ( j j ; GLYPH&lt;17&gt; ) i = 1 6 as the average dimer density per = bond, these equations are solved by

$$c _ { \eta } ( \mathbf p ) = \sqrt { \frac { 6 } { N } } \, \frac { 1 + e ^ { i p _ { \eta } } } { \sqrt { \sum _ { \xi } | 1 + e ^ { i p _ { \zeta } } | ^ { 2 } } } \,.$$

The above construction can be straightforwardly generalized to an arbitrary number Nf of fermionic dimers, since any of the projectors in the Hamiltonian (4) applied to a plaquette with more than one fermionic dimer gives zero. Accordingly the states

$$| \mathbf p _ { 1 }, \dots \mathbf p _ { N _ { f } } \rangle = \sum _ { j _ { 1 }, \eta _ { 1 }, \dots j _ { N _ { f } }, \eta _ { N _ { f } } } e ^ { i \sum _ { j } \mathbf p _ { j }, \mathbf R _ { j } } c _ { \eta _ { 1 } } ( \mathbf p _ { 1 } ) \dots c _ { \eta _ { N _ { f } } } ( \mathbf p _ { N _ { f } } ) | ( j _ { 1 }, \eta _ { 1 } ), \dots, ( j _ { N _ { f } }, \eta _ { N _ { f } } ) \rangle \,,$$

with the coe GLYPH&lt;14&gt; cients c GLYPH&lt;17&gt; ( p ) given in Eq. (18) are exact, zero-energy ground states of the Hamiltonian (4) with Nf fermionic dimers. Also note that these states are anti-symmetric as required, since the basis states j ( j 1 ; GLYPH&lt;17&gt; 1) ; : : : ; ( j Nf ; GLYPH&lt;17&gt; Nf ) i = F j y 1 ;GLYPH&lt;17&gt; 1 D j 1 ;GLYPH&lt;17&gt; 1 : : : F jN f y ;GLYPH&lt;17&gt; N f D jN f ;GLYPH&lt;17&gt; N f j RK i are anti-symmetric under the exchange of two fermionic dimers.

## 4 Perturbing away from the exactly solvable line: Z2-FL*

The massive ground state degeneracy is immediately lifted away from the exactly solvable line. Within first order perturbation theory the ground states j p i obtain an energy shift " ( p ), which can be interpreted as a fermion dispersion. Accordingly, at a finite density of fermionic dimers the ground state corresponds to a Fermi sea with the lowest energy states j p i occupied up to the Fermi energy. This state realizes a fractionalized Fermi liquid with Z2 topological order ( Z2-FL*) and a small Fermi surface, the volume of which is determined by the density of fermionic dimers, i.e. the density of holes ( p ) away from half filling, rather than the total density of holes (1 + p ) measured from the completely filled band.

Here we compute the dispersion " ( p ) for small deviations of the amplitude t 1 ! t 1 + GLYPH&lt;14&gt; t 1 in the Hamiltonian (4) away from the exactly solvable line v 1 = t 2 = GLYPH&lt;0&gt; t 1 using first order perturbation theory. A straightforward computation gives

$$\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \\ \varepsilon ( \mathbf p ) \ & = \ \langle \mathbf p | \Delta H | \mathbf p \rangle \\ & = \ - \delta t _ { 1 } \sum _ { i, \eta, \cdot, j, \eta, \ell } e ^ { - i \mathbf p \cdot ( \mathbf R, - \mathbf R ) } c _ { \eta, } ^ { * } ( \mathbf p ) c _ { \eta, } ( \mathbf p ) \, \langle ( i, \eta _ { i } ) | F _ { \ell, 1, \sigma } ^ { \dagger } D _ { \ell + e _ { 2 }, 1, \sigma } ^ { \dagger } F _ { \ell + e _ { 2 }, 1, \sigma } D _ { \ell, 1 } + \dots | ( j, \eta _ { j } ) \rangle \\ & = \ - \frac { 2 \, \delta t _ { 1 } } { 3 } \frac { \epsilon _ { \mu \nu \lambda } ^ { 2 } ( 1 + \cos p _ { \mu } ) ( \cos p _ { \nu } + \cos p _ { \lambda } ) } { \sum _ { \zeta } | 1 + e ^ { i p _ { \zeta } } | ^ { 2 } }$$

Figure 4: Fermion dispersion relation " ( k ) as function of momentum along the high symmetry line GLYPH&lt;0&gt; GLYPH&lt;0&gt; M GLYPH&lt;0&gt; K GLYPH&lt;0&gt; GLYPH&lt;0&gt; in the vicinity of the exactly solvable line for GLYPH&lt;14&gt; t 1 = 0 01. Blue solid line: perturbative : result from Eq. (20), red line with circles: exact diagonalization results on a 4 GLYPH&lt;2&gt; 4 lattice, green line with squares: exact diagonalization results on a 6 GLYPH&lt;2&gt; 6 lattice.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00103v2-figure-71e773ea-0be4-4d0f-9d9b-a72a709c804f.jpg)

where GLYPH&lt;15&gt; GLYPH&lt;22&gt;GLYPH&lt;23&gt;GLYPH&lt;21&gt; is the antisymmetric tensor and sums over repeated greek indices GLYPH&lt;22&gt;; GLYPH&lt;23&gt; ; GLYPH&lt;21&gt; 2 f 1 2 3 ; ; g are implied. In the last line we used that h ( '; 1) ( j ' + e 2 ; 1) i = Qc [( '; 1) ; ( ' + e 2 ; 1)] = 1 18 is the clas-= sical dimer correlation function of two parallel dimers on a plaquette, which can be computed using standard Grassmann techniques [28].

In Fig. 4 we check the validity of the perturbative expansion by comparing Equ. (20) with numerical results. Shown is the dispersion along the high symmetry line GLYPH&lt;0&gt; GLYPH&lt;0&gt; M GLYPH&lt;0&gt; K GLYPH&lt;0&gt; GLYPH&lt;0&gt; for GLYPH&lt;14&gt; t 1 = GLYPH&lt;0&gt; 0 01. The : numerical results were obtained using Lanczos exact diagonalization of the Hamiltonian (3), where we computed the ground state energy of a single fermionic dimer with fixed momentum p interacting with a background of bosonic dimers on lattices of size 4 GLYPH&lt;2&gt; 4 and 6 GLYPH&lt;2&gt; 6 with twisted boundary conditions to increase the momentum resolution. The numerical results show reasonably good agreement with the perturbative result in Eq. (20) and the small discrepancies can be attributed to finite size e GLYPH&lt;11&gt; ects.

Fig. 5 shows a contour plot of the perturbative dispersion " ( p ) for GLYPH&lt;14&gt; t 1 = GLYPH&lt;0&gt; 0 25. : Note that the dispersion minima are located at the M points on the edges of the Brillouin zone. This feature persists away from the exactly solvable line and is a generic feature of the fermion dispersion for realistic parameter choices. Indeed, the dimer resonance amplitudes can be estimated by computing matrix elements of a simple nearest-neighbor tight binding Hamiltonian Ht = GLYPH&lt;0&gt; t P h i ; j i c c j y i of electrons on the triangular lattice of the form

$$t _ { 1 } \ \simeq \ - \langle \oint \nabla | H _ { t } | \nabla \nabla \rangle = - \frac { 3 t } { 4 }$$

$$t _ { 2 } \ \simeq \ - \langle \bigcap | H _ { t } | \bigcup \rangle = \frac { t } { 2 }$$

where t is the nearest-neighbor electron hopping amplitude. Note in particular that t 1 is negative and 50% larger in magnitude than t 2. This puts the parameters in the vicinity of the exactly solvable line, where t 1 = GLYPH&lt;0&gt; t 2. Accordingly we might suspect that the exact solution captures physical properties in the realistic parameter regime. In Fig. 6 we show a plot of the numerically obtained fermion dispersion for the choice of parameters J = V = 1, t 1 = GLYPH&lt;0&gt; 3, t 2 = 2 and v 1 = 2, corresponding to a nearest neighbor electron hopping amplitude t = 4 J . Again, the data was obtained using Lanczos exact diagonalization of the Hamiltonian (3) for a single fermionic dimer with fixed momentum on a lattice of size 4 GLYPH&lt;2&gt; 4 with twisted boundary conditions. Note that the minima of the fermion dispersion remain at the M points for this choice of parameters.

Figure 5: Plot of the perturbative result for the fermion dispersion relation " ( k ) from Eq. (20) for GLYPH&lt;14&gt; t 1 = GLYPH&lt;0&gt; 0 25. :

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00103v2-figure-16b94393-efbe-4528-8779-0326353e20d4.jpg)

Figure 6: Exact diagonalization result for the dispersion of a single fermionic dimer interacting with a background of bosonic dimers on a lattice of 4 GLYPH&lt;2&gt; 4 sites with twisted boundary conditions. Parameters: J = V = 1, t 1 = GLYPH&lt;0&gt; 3, t 2 = 2 and v 1 = 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00103v2-figure-570b72f4-ef69-4082-a29e-d8e4d6b02600.jpg)

It is also important to note that the Fermi surface of the fermionic dimers as obtained from the dispersion " ( p ) in Eq. (20) has a direct imprint on the electronic Fermi surface. The argument is analogous to Ref. [23] and here we focus on the main idea. Defining the operators

$$f _ { \mathfrak { p }, \sigma } ^ { \dagger } & = \sum _ { j, \eta } e ^ { i \mathfrak { p } \cdot \mathbf R _ { j } } c _ { \eta } ( \mathfrak { p } ) F _ { j, \eta, \sigma } ^ { \dagger } D _ { j, \eta }$$

and acting with them on the bosonic RK ground state j RK i generates the exact ground states in Equ. (19). Note that we've reintroduced the spin label here. The operators f y p ;GLYPH&lt;27&gt; and f p ;GLYPH&lt;27&gt; obey canonical anti-commutation relations GLYPH&lt;8&gt; f p ;GLYPH&lt;27&gt; ; f y q ;GLYPH&lt;27&gt; 0 GLYPH&lt;9&gt; = GLYPH&lt;14&gt; p q ; GLYPH&lt;14&gt; GLYPH&lt;27&gt; ;GLYPH&lt;27&gt; 0 in the thermodynamic limit on the Hilbert space spanned by the states in Equ. (19). This can be shown by computing GLYPH&lt;12&gt; GLYPH&lt;12&gt; GLYPH&lt;12&gt; GLYPH&lt;8&gt; f p ;GLYPH&lt;27&gt; ; f y q ;GLYPH&lt;27&gt; 0 GLYPH&lt;9&gt; j RK i GLYPH&lt;12&gt; GLYPH&lt;12&gt; GLYPH&lt;12&gt; 2 and realizing that it is proportional to the Fourier transform of the classical dimer correlation function. Within first order perturbation theory the Hamiltonian thus takes the form of a free fermion Hamiltonian H = P p ;GLYPH&lt;27&gt; " ( p ) f y p ;GLYPH&lt;27&gt; f p ;GLYPH&lt;27&gt; in the vicinity of the exactly solvable line.

Lastly, the operators f y p ;GLYPH&lt;27&gt; can be directly related to electron annihilation operators c p ;GLYPH&lt;27&gt; , which can be seen as follows. On the dimer Hilbert space the electron destruction operator cj ;GLYPH&lt;27&gt; on lattice site j can be uniquely written in terms of the dimer operators in Eqs. (1), (2) as [20]

$$c _ { j, \sigma } = \frac { \epsilon _ { \sigma \sigma ^ { \prime } } } { 2 } \sum _ { \eta } \left ( F _ { j, \eta, \sigma ^ { \prime } } ^ { \dagger } D _ { j, \eta } + F _ { j - \mathbf e _ { \eta }, \eta, \sigma ^ { \prime } } ^ { \dagger } D _ { j - \mathbf e _ { \eta }, \eta } \right ) \,,$$

which can be checked by computing matrix elements on both sides of the equation in the dimer Hilbert space. Fourier transforming this relation it is straightforward to show that

$$c _ { \mathbf p, \sigma } = \epsilon _ { \sigma, \sigma ^ { \prime } } \sqrt { Z _ { \mathbf p } } f ^ { \dagger } _ { - \mathbf p, \sigma ^ { \prime } }$$

where Z p = P GLYPH&lt;17&gt; j 1 + e GLYPH&lt;0&gt; ip GLYPH&lt;17&gt; j 2 = 24 is the electronic quasiparticle residue. In the vicinity of the exactly solvable line the electron spectral function thus takes the form

$$A _ { e } ( \mathbf p, \omega ) = Z _ { \mathbf p } \, \delta ( \omega - \varepsilon ( \mathbf p ) + \mu )$$

where we introduced the chemical potential GLYPH&lt;22&gt; to fix the average density of fermionic dimers. This demonstrates that the small Fermi surface of fermionic dimers is directly imprinted on the electronic Fermi surface and proves that the ground state is indeed a fractionalized Fermi liquid with Z2 topological order. It features a small Fermi surface with an enclosed volume proportional to the density of doped holes away from half filling in the absence of any broken symmetries.

## 5 Application to Twisted Bilayer Graphene

In this section we comment on a possible application of the previously introduced triangular lattice quantum dimer model in the context of magic-angle twisted bilayer graphene (TBG). Recent experiments on TBG showed correlated Mott-like insulating behavior at half filling of the lower Moir´ e mini-band near the charge neutrality point [6,29,30]. Moreover, TBG becomes superconducting upon doping the correlated insulator with electons or holes. Interestingly, Shubnikov - de Haas as well as Hall resistivity measurements above the superconducting transition temperature on the hole-doped side of the correlated insulator indicate an anomalously small charge carrier density, which is equal to the density of doped carriers measured from the correlated insulator at half filling [6]. This is in

remarkable analogy to Hall measurements in the pseudogap phase of underdoped cuprates [5]. One possible explanation of a small carrier density involves translational symmetry breaking orders, which enlarge the size of the unit cell and reconstruct the Fermi surface into small pockets. A di GLYPH&lt;11&gt; erent possibility is the presence of topological order, which can give rise to small pocket Fermi surfaces and a small carrier density as well, as discussed in this work. So far no symmetry breaking orders have been observed in TBG which could explain the small carrier density, consequently the Z2-FL* scenario discussed here could provide a viable explanation.

The microscopic details of TBG are rather complex and a lot of e GLYPH&lt;11&gt; ort has been put into developing realistic tight binding models [31-34]. In the following we argue that the triangular lattice dimer model studied here can be viewed as a very simplistic toy model to gain intuition about the unconventional metallic state near the Mott-like insulator at half filling of the lower Moir´ mini-band. Indeed, electric e charge in TBG is localized at the AA stacked regions of the two twisted graphene layers, which form a triangular lattice [35]. Even though symmetry constraints preclude the definition of a triangular lattice tight binding model for TBG and microscopic models have a lower rotational symmetry, it is expected that the six-fold rotational symmetry of the triangular lattice emerges in e GLYPH&lt;11&gt; ective low-energy theories. As such, triangular lattice models like the one discussed here should be able to capture e GLYPH&lt;11&gt; ects of local correlations on the low-energy properties of TBG [36].

An additional important complication is the two-fold valley degeneracy in TBG. In our simplistic picture this would lead to two degenerate orbitals per triangular lattice site due to the valley degree of freedom, i.e. a four fold degeneracy per lattice site in total, if spin is included. This fourfold degeneracy is indeed observed in the Landau level structure emanating from the charge neutrality point [6]. However, the Landau level degeneracy in the metallic state emanating from the Mott-like insulator at half filling is only two-fold. It has been argued that this could be due to the presence of so-called inter-valley coherent order, where the valley quantum number is lost [33]. The spin quantum number is likely to remain una GLYPH&lt;11&gt; ected, since the superconductor, which develops out of this unconventional metallic state, can be suppressed by a parallel magnetic field and thus spin-singlet pairing is probable. Using this scenario as a basis, our triangular quantum dimer model can be viewed as a simple toy-model for the remaining two-fold degenerate spin-states per triangular lattice site in a conjectured inter-valley coherent phase and o GLYPH&lt;11&gt; ers a possible explanation for the low-charge carrier density upon hole-doping the Mott-like insulator. As discussed in the main part of this manuscript, a prominent signature of the Z2-FL* state are the small hole-pockets centered at the M points of the Moir´ e mini-Brillouin zone in the absence of translational symmetry breaking orders. However, the Fermi-pockets at the three distinct M points should give rise to an additional three-fold degeneracy of the Landau-levels, which is not observed in experiments. This could be due to two reasons: either the Fermi surface undergoes an additional reconstruction in the presence of a strong magnetic field, or our simple estimate of the dimer hopping amplitudes from a single-band, nearest-neighbor triangular tightbinding model does not hold for TBG and the dispersion minimum is in fact at the GLYPH&lt;0&gt; point. In any case, the location of the small Fermi pockets in TBG could be determined in principle using either ARPES measurements with very high energy and momentum resolution, which is currently out of reach, or using quasi-particle interference in STM experiments. We also note here that a di GLYPH&lt;11&gt; erent route towards a description of metallic states with Z2 topological order in TBG is discussed in Ref. [36].

## 6 Discussion and Conclusions

We presented a generalized, two-species quantum dimer model on the triangular lattice which features a fractionalized Fermi liquid ground state with Z2 topological order (Z2-FL*). An exact ground state solution was constructed for a specific line in parameter space. At this exactly solvable line the ground state is highly degenerate and can be interpreted as a flat fermionic band. Using perturbation theory we computed the fermion dispersion away from the exactly solvable line and showed that the ground state is indeed a Z2-FL* with a modified Luttinger count, where the electronic Fermi surface encloses a volume proportional to the density of doped holes away from half filling. In particular we showed that the Fermi surface consists of small hole pockets in the vicinity of the M points on the edges of the Brillouin zone at low doping. This feature is quite robust for realistic parameter choices derived from a single-band tight-binding model of electrons on a triangular lattice and can be used as a signature of the Z2-FL* state on the triangular lattice.

Even though we focused on the unconventional metallic Z2-FL* state in this work, the complete phase diagram of this model also contains symmetry broken states that we did not consider here. In the undoped case, transitions to symmetry broken valence bond solid (VBS) phases can be described as confinement transitions of the e GLYPH&lt;11&gt; ective Z2 gauge theory for the resonating valence bond phase [37,38]. It is an interesting open problem how such confinement transitions are a GLYPH&lt;11&gt; ected by the presence of a finite density of fermionic dimers. Related questions on the square lattice have been studied recently using sign-problem free quantum Monte-Carlo [39, 40] and it would be interesting to extend this approach to the model studied here.

## Acknowledgements

We are grateful to Johannes Feldmeier and Sebastian Huber for helpful discussions.

Funding information This research was supported by the German Excellence Initiative via the Nano Initiative Munich (NIM) as well as by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy - EXC-2111 - 390814868.

## A Kagome lattice dimer model

The construction presented in Sec. 3 can be straightforwardly generalized to di GLYPH&lt;11&gt; erent lattice geometries, such as the kagome lattice, which we will briefly discuss in this appendix. A general dimer Hamiltonian again consists of various plaquette resonance terms, in analogy to Eq. 3 on the triangular lattice. We consider terms on the four di GLYPH&lt;11&gt; erent plaquettes shown in Fig. 7, namely the elementary hexagon, as well as the three diamond shaped loops of length eight. Instead of specifying the various resonance terms in the Hamiltonian in detail, we directly define an exactly solvable Hamiltonian in terms of plaquette projectors, as in Eq. (4). Again, we keep an RK-like term H RK, which ensures that the ground state without fermionic dimers is an equal weight superposition of all bosonic dimer coverings. Moreover, we define the four plaquette projectors Pj ;GLYPH&lt;17&gt; = j GLYPH&lt;30&gt; j ;GLYPH&lt;17&gt; ih GLYPH&lt;30&gt; j ;GLYPH&lt;17&gt; j involving a fermionic

Figure 7: Kagome lattice unit cell with independent bonds enumerated from one to six (left). We consider four di GLYPH&lt;11&gt; erent plaquette terms: the elementary hexagon (marked in blue), as well as the three diamond shaped loops of length 8 (marked by green dots, red stripes and blue bricks).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00103v2-figure-2eb86c43-eeea-4620-83f3-54f590f451fd.jpg)

dimer, where

$$| \phi _ { j, 1 } \rangle \ = \ | \nwarrow \psi \rangle \ + \ | \nwarrow \psi \rangle \ + \ | \nwarrow \psi \rangle \ + \ | \nwarrow \psi \rangle \ + \ | \nwarrow \psi \rangle$$

$$| \phi _ { j, 2 } \rangle \ = \ | \underline { \phi } \nearrow \rangle + | \underline { \phi } \nearrow \rangle + | \underline { \phi } \nearrow \rangle + | \underline { \phi } \nearrow \rangle + | \underline { \phi } \nearrow \rangle + | \underline { \phi } \nearrow \rangle \quad ( 2 8 )$$

and analogous terms j GLYPH&lt;30&gt; j ; 3 i and j GLYPH&lt;30&gt; j ; 4 i for the other two diamond shaped plaquettes. Note that we have six independent bonds per kagome unit cell (see Fig. 7) and thus we need to determine six coe GLYPH&lt;14&gt; cients c GLYPH&lt;17&gt; ( p ) with GLYPH&lt;17&gt; 2 f 1 ; : : : ; 6 g in our ansatz for the ground state wave function in Eq. (12). Consequently we would need five linearly independent plaquette projectors together with the normalization condition in order to uniquely specify the coe GLYPH&lt;14&gt; cients c GLYPH&lt;17&gt; ( p ). With our four plaquette operators the coe GLYPH&lt;14&gt; cients are underdetermined, but we can define exact ground states nonetheless.

Applying the four projectors to our ansatz (12) and requiring that the eigenenergy is zero leads to the four equations

$$c _ { 1 } ( \mathbf p ) + c _ { 2 } ( \mathbf p ) + c _ { 3 } ( \mathbf p ) + c _ { 4 } ( \mathbf p ) + c _ { 5 } ( \mathbf p ) + c _ { 6 } ( \mathbf p ) \ = \ 0 \, \quad ( 2 9 )$$

$$c _ { 1 } ( \mathbf p ) \left ( 1 + e ^ { i 2 p _ { 1 } } \right ) + c _ { 2 } ( \mathbf p ) \left ( 1 + e ^ { i 2 p _ { 3 } } \right ) + c _ { 4 } ( \mathbf p ) \left ( 1 + e ^ { - i 2 p _ { 1 } } \right ) + c _ { 5 } ( \mathbf p ) \left ( 1 + e ^ { - i 2 p _ { 3 } } \right ) \ = \ 0 \, \quad ( 3 0 )$$

$$c _ { 2 } ( \mathbf p ) \left ( 1 + e ^ { - i 2 p _ { 2 } } \right ) + c _ { 3 } ( \mathbf p ) \left ( 1 + e ^ { - i 2 p _ { 1 } } \right ) + c _ { 5 } ( \mathbf p ) \left ( 1 + e ^ { j 2 p _ { 2 } } \right ) + c _ { 6 } ( \mathbf p ) \left ( 1 + e ^ { j 2 p _ { 1 } } \right ) \ = \ 0 \,, \quad ( 3 1 )$$

$$c _ { 1 } ( \mathbf p ) \left ( 1 + e ^ { - i 2 p _ { 2 } } \right ) + c _ { 3 } ( \mathbf p ) \left ( 1 + e ^ { i 2 p _ { 3 } } \right ) + c _ { 4 } ( \mathbf p ) \left ( 1 + e ^ { i 2 p _ { 2 } } \right ) + c _ { 6 } ( \mathbf p ) \left ( 1 + e ^ { - i 2 p _ { 3 } } \right ) \ = \ 0 \. \quad ( 3 2 )$$

Solving these equations for c 2 3 5 6 ; ; ; together with the normalization condition Eq. (17), where the mean dimer density per bond on the kagome lattice is Qc [( i ; GLYPH&lt;17&gt; )] = 1 6, and setting = c 4 = GLYPH&lt;0&gt; c 1 gives a particularly simple solution. We finally obtain

$$c _ { 1 } ( \mathbf p ) \ = \ - c _ { 4 } ( \mathbf p ) \ = \ \sqrt { \frac { 6 } { N } } \frac { | \sin ( 2 p _ { 3 } ) | } { \sqrt { 3 - \cos ( 4 p _ { 1 } ) - \cos ( 4 p _ { 2 } ) - \cos ( 4 p _ { 3 } ) } } \,,$$

$$c _ { 2 } ( \mathbf p ) \ = \ - c _ { 5 } ( \mathbf p ) \ = \ - c _ { 1 } ( \mathbf p ) \, \frac { e ^ { i 2 ( p _ { 1 } - p _ { 3 } ) } \stackrel { \mu } { \left ( 1 - e ^ { - i 4 p _ { 1 } } \right ) } } { \stackrel { 1 - e ^ { - i 4 p _ { 3 } } \right ) } { \stackrel { \mu } { \left ( 1 - e ^ { - i 4 p _ { 3 } } \right ) } } } \,,$$

$$c _ { 3 } ( \mathbf p ) \ = \ - c _ { 6 } ( \mathbf p ) \ = \ c _ { 1 } ( \mathbf p ) \, \frac { e ^ { i 2 ( p _ { 2 } - p _ { 3 } ) } ( 1 - e ^ { - i 4 p _ { 2 } } ) } { 1 - e ^ { - i 4 p _ { 3 } } } \,.$$

Again, this solution can be straightforwardly extended to an arbitrary number of fermionic dimers.

## References

- [1] M. Gurvitch and A. T. Fiory, Resistivity of La 1 825 : Sr 0 175 : CuO 4 and YBa Cu O 2 3 7 to 1100 K: Absence of saturation and its implications , Phys. Rev. Lett. 59 , 1337 (1987), doi:10.1103 PhysRevLett.59.1337. /
- [2] A. Legros, S. Benhabib, W. Tabis, F. Lalibert´ e, M. Dion, M. Lizaire, B. Vignolle, D. Vignolles, H. Ra GLYPH&lt;11&gt; y, Z. Z. Li, P. Auban-Senzier, N. Doiron-Leyraud, P. Fournier, D. Colson, L. Taillefer, C. Proust, Universal T-linear resistivity and Planckian dissipation in overdoped cuprates , Nature Physics 15 , 142?147 (2019), doi:10.1038 s41567-018-0334-2. /
- [3] S. I. Mirzaei, D. Stricker, J. N. Hancock, C. Berthod, A. Georges, E. van Heumen, M. K. Chan, X. Zhao, Y. Li, M. Greven, N. Barisi´ c, D. van der Marel, Spectroscopic evidence for Fermi liquid-like energy and temperature dependence of the relaxation rate in the pseudogap phase of the cuprates , Proc. Natl. Acad. Sci. 110 , 5774-5778 (2013), doi:10.1073 pnas.1218846110. /
- [4] M. K. Chan, M. J. Veit, C. J. Dorow, Y. Ge, Y. Li, W. Tabis, Y. Tang, X. Zhao, N. Barisi´ c, and M. Greven, In-Plane Magnetoresistance Obeys Kohler's Rule in the Pseudogap Phase of Cuprate Superconductors , Phys. Rev. Lett. 113 , 177005 (2014), doi:10.1103 PhysRevLett.113.177005. /
- [5] S. Badoux, W. Tabis, F. Lalibert´ e, G. Grissonnanche, B. Vignolle, D. Vignolles, J. B´ eard, D. A. Bonn, W. N. Hardy, R. Liang, N. Doiron-Leyraud, L. Taillefer, C. Proust, Change of carrier density at the pseudogap critical point of a cuprate superconductor , Nature 531 , 210 (2016), doi:10.1038 nature16983. /
- [6] Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxiras, P. Jarillo-Herrero, Magicangle graphene superlattices: a new platform for unconventional superconductivity , Nature 556 , 43-50 (2018), doi:10.1038 nature26160. /
- [7] Y. Cao, D. Chowdhury, D. Rodan-Legrain, O. Rubies-Bigorda, K. Watanabe, T. Taniguchi, T. Senthil, P. Jarillo-Herrero, Strange metal in magic-angle graphene with near Planckian dissipation , http://arxiv.org/abs/1901.03710 .
- [8] C. M. Varma, P. B. Littlewood, S. Schmitt-Rink, E. Abrahams, and A. E. Ruckenstein, Phenomenology of the normal state of Cu-O high-temperature superconductors , Phys. Rev. Lett. 63 , 1996 (1989), doi:10.1103 PhysRevLett.63.1996. /
- [9] H. v. L¨hneysen, A. Rosch, M. Vojta, and P. W¨lfle, o o Fermi-liquid instabilities at magnetic quantum phase transitions , Rev. Mod. Phys. 79 , 1015 (2007), doi:10.1103 RevModPhys.79.1015. /
- [10] S. Sachdev and J. Ye, Gapless spin-fluid ground state in a random quantum Heisenberg magnet , Phys. Rev. Lett. 70 , 3339 (1993), doi:10.1103 PhysRevLett.70.3339. /
- [11] S. Sachdev, Bekenstein-Hawking Entropy and Strange Metals , Phys. Rev. X 5 , 041025 (2015), doi:10.1103 PhysRevX.5.041025. /
- [12] P. W. Anderson, The Resonating Valence Bond State in La CuO 2 4 and Superconductivity , Science 235 , 1196-1198 (1987), doi:10.1126 science.235.4793.1196. /
- [13] G. Kotliar, A. E. Ruckenstein, New Functional Integral Approach to Strongly Correlated Fermi Systems: The Gutzwiller Approximation as a Saddle Point , Phys. Rev. Lett. 57 , 1362 (1986), doi:10.1103 PhysRevLett.57.1362. /

- [14] T. Senthil, M. P. A. Fisher, Z 2 gauge theory of electron fractionalization in strongly correlated systems , Phys. Rev. B 62 , 7850 (2000), doi:10.1103 PhysRevB.62.7850. /
- [15] S. Florens, A. Georges, Slave-rotor mean-field theories of strongly correlated systems and the Mott transition in finite dimensions , Phys. Rev. B 70 , 035114 (2004), doi:10.1103 PhysRevB.70.035114. /
- [16] R. Nandkishore, M. A. Metlitski, T. Senthil, Orthogonal metals: The simplest non-Fermi liquids , Phys. Rev. B 86, 045128 (2012), doi:10.1103 PhysRevB.86.045128. /
- [17] S. Sachdev, Topological order, emergent gauge fields, and Fermi surface reconstruction , Reports on Progress in Physics 82 , 014001 (2019), doi:10.1088 1361-6633 aae110. / /
- [18] M. Oshikawa, Topological Approach to Luttinger's Theorem and the Fermi Surface of a Kondo Lattice , Phys. Rev. Lett. 84 , 3370 (2000), doi:10.1103 PhysRevLett.84.3370. /
- [19] T. Senthil, S. Sachdev, M. Vojta, Fractionalized Fermi Liquids , Phys. Rev. Lett. 90 , 216403 (2003), doi:10.1103 PhysRevLett.90.216403. /
- [20] M. Punk, A. Allais, S. Sachdev, Quantum dimer model for the pseudogap metal , Proc. Natl. Acad. Sci. 112 9552-9557 (2015), doi:10.1073 pnas.1512206112. /
- [21] D. S. Rokhsar, S. A. Kivelson, Superconductivity and the Quantum Hard-Core Dimer Gas , Phys. Rev. Lett. 61 , 2376 (1988), doi:10.1103 PhysRevLett.61.2376. /
- [22] S. Huber, J. Feldmeier, M. Punk, Electron spectral functions in a quantum dimer model for topological metals , Phys. Rev. B 97 , 075144 (2018), doi:110.1103 PhysRevB.97.075144. /
- [23] J. Feldmeier, S. Huber, M. Punk, Exact solution of a two-species quantum dimer model for pseudogap metals , Phys. Rev. Lett. 120 , 187001 (2018), doi:10.1103 PhysRevLett.120.187001. /
- [24] A. M. Polyakov, Quark confinement and topology of gauge theories , Nucl. Phys. B 120 , 429 (1977), doi:10.1016 0550-3213(77)90086-4. /
- [25] M. Hermele, T. Senthil, M. P. A. Fisher, P. A. Lee, N. Nagaosa, X.-G. Wen, Stability of U(1) spin liquids in two dimensions , Phys. Rev. B 70 , 214437 (2004), doi:10.1103 PhysRevB.70.214437. /
- [26] R. Moessner , S. L. Sondhi, Resonating Valence Bond Phase in the Triangular Lattice Quantum Dimer Model , Phys. Rev. Lett. 86 , 1881 (2001), doi:10.1103 PhysRevLett.86.1881. /
- [27] S. Sachdev, D. Chowdhury, Progress of Theoretical and Experimental Physics 12C102 (2016), doi:10.1093 ptep ptw110. / /
- [28] P. Fendley, R. Moessner, S. L. Sondhi, Classical dimers on the triangular lattice , Phys. Rev. B 66 , 214513 (2002), doi:10.1103 PhysRevB.66.214513. /
- [29] Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, R. C. Ashoori, P. Jarillo-Herrero, Correlated Insulator Behaviour at Half-Filling in Magic Angle Graphene Superlattices , Nature 556 , 80-84 (2018), doi:10.1038 nature26154. /

- [30] X. Lu, P. Stepanov, W. Yang, M. Xie, M. A. Aamir, I. Das, C. Urgell, K. Watanabe, T. Taniguchi, G. Zhang, A. Bachtold, A. H. MacDonald, D. K. Efetov, Superconductors, Orbital Magnets, and Correlated States in Magic Angle Bilayer Graphene , http://arxiv.org/abs/1901.03710 .
- [31] M. Koshino, N. F. Q. Yuan, T. Koretsune, M. Ochi, K. Kuroki, L. Fu, Maximally-localized Wannier orbitals and the extended Hubbard model for the twisted bilayer graphene , Phys. Rev. X 8 , 031087 (2018), doi:10.1103 PhysRevX.8.031087. /
- [32] J. Kang, O. Vafek, Symmetry, maximally localized Wannier states, and low energy model for the twisted bilayer graphene narrow bands , Phys. Rev. X 8 , 031088 (2018), doi:10.1103 PhysRevX.8.031088. /
- [33] H. C. Po, L. Zou, A. Vishwanath, T. Senthil, Origin of Mott insulating behavior and superconductivity in twisted bilayer graphene , Phys. Rev. X 8 , 031089 (2018), doi:10.1103 PhysRevX.8.031089. /
- [34] H. C. Po, L. Zou, T. Senthil, A. Vishwanath, Faithful Tight-binding Models and Fragile Topology of Magic-angle Bilayer Graphene , Phys. Rev. B 99 , 195455 (2019), doi:10.1103 PhysRevB.99.195455. /
- [35] A. Luican, G. Li, A. Reina, J. Kong, R. R. Nair, K. S. Novoselov, A. K. Geim, and E. Y. Andrei, Single-Layer Behavior and Its Breakdown in Twisted Graphene Layers , Phys. Rev. Lett. 106 , 126802 (2011), doi:10.1103 PhysRevLett.106.126802. /
- [36] A. Thomson, S. Chatterjee, S. Sachdev, and M. S. Scheurer, Triangular antiferromagnetism on the honeycomb lattice of twisted bilayer graphene , Phys. Rev. B 98 , 075109 (2018), doi:10.1103 PhysRevB.98.075109. /
- [37] F. Wegner, Duality in Generalized Ising Models and Phase Transitions without Local Order Parameter , J. Math. Phys. 12 , 2259 (1971), doi:10.1063 1.1665530. /
- [38] Y. Huh, M. Punk, S. Sachdev, Vison states and confinement transitions of Z2 spin liquids on the kagome lattice , Phys. Rev. B 94 , 094419 (2011), doi:10.1103 PhysRevB.94.094419. /
- [39] S. Gazit, M. Randeria, A. Vishwanath, Emergent Dirac fermions and broken symmetries in confined and deconfined phases of Z2 gauge theories , Nature Physics 13 , 484 (2017), doi:10.1038 nphys4028. /
- [40] S. Gazit, F. F. Assaad, S. Sachdev, A. Vishwanath, C. Wang, Confinement transition of Z2 gauge theories coupled to massless fermions: emergent QCD3 and SO(5) symmetry , Proc. Natl. Acad. Sci. 115 , E6987 (2018), doi:10.1073 pnas.1806338115. /
Brin Verheijden Yuhao Zhao Matthias Punk
Solvable lattice models for metals with Z2 topological order
We present quantum dimer models in two dimensions which realize metallic
ground states with Z2 topological order. Our models are generalizations of a
dimer model introduced in [PNAS 112,9552-9557 (2015)] to provide an effective
description of unconventional metallic states in hole-doped Mott insulators. We
construct exact ground state wave functions in a specific parameter regime and
show that the ground state realizes a fractionalized Fermi liquid. Due to the
presence of Z2 topological order the Luttinger count is modified and the volume
enclosed by the Fermi surface is proportional to the density of doped holes
away from half filling. We also comment on possible applications to magic-angle
twisted bilayer graphene.