Search

Due to the presence of flat Chern bands, moir\'e transition metal dichalcogenide (TMD) bilayers are a platform to realize strongly correlated topological phases of fermions such as fractional Chern insulators. TMDs are also known to host long-lived excitons, which inherit the topology of the underlying Chern bands. For the particular example of MoTe$_2$/WSe$_2$ heterobilayers we perform a time-dependent Hartree-Fock calculation to identify a regime in the phase diagram where the excitons themselves form a topological flat band. This paves a way towards realizing strongly correlated states of bosons in moir\'e TMDs., Comment: 7 pages, 4 figures (4 pages in supplement)

Electrons in a quantum Hall fluid can bind with an integer number of vortices to form composite fermions and composite bosons. We show that the quantum Hall bilayer at filling $\nu=\frac{1}{4}+\frac{3}{4}$ with interlayer separation $d$ can be well-described in terms of these composite particles. At small $d$ the system can be understood as interlayer paired electrons and holes, whereas at large $d$ the system is best understood in terms of composite fermions with four vortices attached to each electron. By computing the overlaps of trial wavefunctions with the ground state from exact diagonalization, we find that as $d$ increases, the number of vortices that attach to each electron increases. We also construct trial states for two types of excitation, the Goldstone mode and a meron excitation. These two trial states have good overlaps with the lowest excited states in the exact diagonalization spectrum for small and intermediate $d$ respectively., Comment: 6 pages, 2 figures (3 pages, 4 figures in supplement)

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

We introduce and study new interacting topological states that arise in time-reversal symmetric bands with an underlying obstruction to forming localized states. If the $U(1)$ valley symmetry linked to independent charge conservation in each time-reversal sector is spontaneously broken, the corresponding `excitonic' order parameter is forced to form a topologically non-trivial texture across the Brillouin zone. We show that the resulting phase, which we dub a textured exciton insulator, cannot be given a local-moment description due to a form of delicate topology. Using toy models of bands with Chern or Euler obstructions to localization we construct explicit examples of the Chern or Euler texture insulators (CTIs or ETIs) they support, and demonstrate that these are generically competitive ground states at intermediate coupling. We construct field theories that capture the response properties of these new states. Finally, we identify the incommensurate Kekul\'e spiral phase observed in magic-angle bi- and trilayer graphene as a concrete realization of an ETI., Comment: 27+17 pages

We explore the phase diagrams of moir\'e materials in search of a new class of intervalley-coherent correlated insulating state: the Chern texture insulator (CTI). This phase of matter, proposed in a companion paper, breaks valley $U(1)$ symmetry in a nontrivial fashion wherein the valley order parameter is forced to texture in momentum space as a consequence of band topology. Using detailed Hartree-Fock studies, we establish that the CTI emerges as an energetically competitive intermediate-coupling ground state in several moir\'e systems which lack a twofold rotation symmetry that forbids the single-particle topology essential to the formation of the CTI valley texture., Comment: 12 + 2 pages

Graphene moir\'e systems are ideal environments for investigating complex phase diagrams and gaining fundamental insights into the mechanisms underlying exotic states of matter, as they permit controlled manipulation of electronic properties. Magic-angle twisted trilayer graphene (MATTG) has emerged as a key platform to explore moir\'e superconductivity, owing to the robustness of its superconducting order and the displacement-field tunability of its energy bands. Recent measurements strongly suggest that superconductivity in MATTG is unconventional. Here, we report the first direct observation of double-dome superconductivity in MATTG. The temperature, magnetic field, and bias current dependence of the superconductivity of doped holes collectively show that it is significantly suppressed near moir\'e filling $\nu^* = -2.6$, leading to a double dome in the phase diagram within a finite window of the displacement field. The temperature dependence of the normal-state resistance and the $I-V$ curves straddling $\nu^*$ are suggestive of a phase transition and the potentially distinct nature of superconductivity in the two domes. Hartree-Fock calculations incorporating mild strain yield an incommensurate Kekul\'e spiral state whose effective spin polarization peaks in the regime where superconductivity is suppressed in experiments. This allows us to draw conclusions about the normal state as well as the unconventional nature of the superconducting order parameter.

We study the phase diagram of magic-angle twisted symmetric trilayer graphene in the presence of uniaxial heterostrain and interlayer displacement field. For experimentally reasonable strain, our mean-field analysis finds robust Kekul\'e spiral order whose doping-dependent ordering vector is incommensurate with the moir\'e superlattice, consistent with recent scanning tunneling microscopy experiments, and paralleling the behaviour of closely-related twisted bilayer graphene (TBG) systems. Strikingly, we identify a new possibility absent in TBG: the existence of $\textit{commensurate}$ Kekul\'e spiral order even at zero strain for experimentally realistic values of the interlayer potential in a trilayer. Our studies also reveal a complex pattern of charge transfer between weakly- and strongly-dispersive bands in strained trilayer samples as the density is tuned by electrostatic gating, that can be understood intuitively in terms of the `cascades' in the compressibility of magic-angle TBG., Comment: 5 + 7 pages; we are informed of an erroneous convention in the preprint version of an experimental paper by private communication from the authors, conclusion unchanged

It is a matter of current debate whether the gate-tunable superconductivity in twisted bilayer graphene is phonon-mediated or arises from electron-electron interactions. The recent observation of the strong coupling of electrons to so-called $K$-phonon modes in angle-resolved photoemission spectroscopy experiments has resuscitated early proposals that $K$-phonons drive superconductivity. We show that the bandwidth-enhancing effect of interactions drastically weakens both the intrinsic susceptibility towards pairing as well as the screening of Coulomb repulsion that is essential for the phonon attraction to dominate at low temperature. This rules out purely $K$-phonon-mediated superconductivity with the observed transition temperature of $\sim 1$ K. We conclude that the unflattening of bands by Coulomb interactions challenges any purely phonon-driven pairing mechanism, and must be addressed by a successful theory of superconductivity in moir\'e graphene., Comment: 8 pages, 2 figures (+16 pages, 8 figures supplement)

Twisted bilayer graphene is an exciting platform for exploring correlated quantum phases, extremely tunable with respect to both the single-particle bands and the interaction profile of electrons. Here, we investigate the phase diagram of twisted bilayer graphene as described by an extended Hubbard model on the honeycomb lattice with two fermionic orbitals (valleys) per site. Besides the special extended {\it cluster interaction} $Q$, we incorporate the effect of gating through an onsite Hubbard-interaction $U$. Within Quantum Monte Carlo (QMC), we find valence-bond-solid, N\'eel-valley antiferromagnetic or charge-density wave phases. Further, we elucidate the competition of these phases by noticing that the cluster interaction induces an exotic constraint on the Hilbert space, which we dub {\it the cluster rule}, in analogy to the famous pyrochlore spin-ice rule. Formulating the perturbative Hamiltonian by projecting into the cluster-rule manifold, we perform exact diagonalization and construct the fixed-point states of the observed phases. Finally, we compute the local electron density patterns as signatures distinguishing these phases, which could be observed with scanning tunneling microscopy. Our work capitalizes on the notion of cluster constraints in the extended Hubbard model of twisted bilayer graphene, and suggests a scheme towards realization of several symmetry-breaking insulating phases in a twisted-bilayer graphene sheet., Comment: 10 pages, 6 figures

Localized or propagating Majorana boundary modes are the key feature of topological superconductors. They are rare in naturally-occurring compounds, but the tailored manipulation of quantum matter offers opportunities for their realization. Specifically, lattices of Yu-Shiba-Rusinov bound states $-$ Shiba lattices $-$ that arise when magnetic adatoms are placed on the surface of a conventional superconductor can be used to create topological bands within the superconducting gap of the substrate. Here, using scanning tunnelling microscopy to create and probe adatom lattices with single atom precision we reveal two signatures consistent with the realization of two types of mirror symmetry protected topological superconductors. The first has edge modes as well as higher-order corner states, and the second has symmetry-protected bulk nodal points. In principle, their topological character and boundary modes should be protected by the spatial symmetries of the adatom lattice. Our results highlight the potential of Shiba lattices as a platform to design the topology and sample geometry of 2D superconductors., Comment: Published in Nature Physics (https://doi.org/10.1038/s41567-023-02104-5)

Despite much experimental and theoretical work, the nature of the charge order in the kagome metals belonging to the family of materials AV$_3$Sb$_5$ (A=Cs,Rb,K) remains controversial. A crucial ingredient for the identification of the ordering in these materials is their response to external perturbations, such as strain or magnetic fields. To this end, we provide a comprehensive symmetry classification of the possible charge orders in kagome materials with a $2\times2$ increase of the unit cell. Motivated by the experimental reports of time-reversal-symmetry breaking and rotational anisotropy, we consider the interdependence of flux and bond orders. Deriving the relevant Landau free energy for possible orders, we study the effect of symmetry-breaking perturbations such as strain and magnetic fields. Our results, thus, provide a roadmap for future tests of these intricate orders., Comment: 15 pages (+7 pages appendix), 4 figures (+6 figures appendix)

The kagome metals AV$_3$Sb$_5$ (A=K,Rb,Cs) present an ideal sandbox to study the interrelation between multiple coexisting correlated phases such as charge order and superconductivity. So far, no consensus on the microscopic nature of these states has been reached as the proposals struggle to explain all their exotic physical properties. Among these, field-switchable electric magneto-chiral anisotropy (eMChA) in CsV$_3$Sb$_5$ provides intriguing evidence for a rewindable electronic chirality, yet the other family members have not been likewise investigated. Here, we present a comparative study of magneto-chiral transport between CsV$_3$Sb$_5$ and KV$_3$Sb$_5$. Despite their similar electronic structure, KV$_3$Sb$_5$ displays negligible eMChA, if any, and with no field switchability. This is in stark contrast to the non-saturating eMChA in CsV$_3$Sb$_5$ even in high fields up to 35 T. In light of their similar band structures, the stark difference in eMChA suggests its origin in the correlated states. Clearly, the V kagome nets alone are not sufficient to describe the physics and the interactions with their environment are crucial in determining the nature of their low-temperature state.

Spontaneously broken symmetries are at the heart of many phenomena of quantum matter and physics more generally. However, determining the exact symmetries broken can be challenging due to imperfections such as strain, in particular when multiple electronic orders form complex interactions. This is exemplified by charge order in some kagome systems, which are speculated to show nematicity and flux order from orbital currents. We fabricated highly symmetric samples of a member of this family, CsV$_3$Sb$_5$, and measured their transport properties. We find the absence of measurable anisotropy at any temperature in the unperturbed material, however, a striking in-plane transport anisotropy appears when either weak magnetic fields or strains are present. A symmetry analysis indicates that a perpendicular magnetic field can indeed lead to in-plane anisotropy by inducing a flux order coexisting with more conventional bond order. Our results provide a unifying picture for the controversial charge order in kagome metals and highlight the need for microscopic materials control in the identification of broken symmetries.

Recent scanning tunneling microscopy experiments [K.P. Nuckolls et al., arXiv:2303.00024] have revealed the ubiquity of Kekul\'e charge-density wave order in magic-angle twisted bilayer graphene. Most samples are moderately strained and show `incommensurate Kekul\'e spiral' (IKS) order involving a graphene-scale charge density distortion uniaxially modulated on the scale of the moir\'e superlattice, in accord with theoretical predictions. However, ultra-low strain samples instead show graphene-scale Kekul\'e charge order that is uniform on the moir\'e scale. This order, especially prominent near filling factor $\nu=-2$, is unanticipated by theory which predicts a time-reversal breaking Kekul\'e current order at low strain. We show that including the coupling of moir\'e electrons to graphene-scale optical zone-corner (ZC) phonons stabilizes a uniform Kekul\'e charge ordered state at $|\nu|=2$ with a quantized topological (spin or anomalous Hall) response. Our work clarifies how this phonon-driven selection of electronic order emerges in the strong-coupling regime of moir\'e graphene., Comment: 5+4 pages

Bilayer quantum Hall states have been shown to be described by a BCS-paired state of composite fermions. However, finding a qualitatively accurate model state valid across all values of the bilayer separation is challenging. Here, we introduce two variational wavefunctions, each with a single variational parameter, which can be thought of as a proxy for the BCS order parameter. Studying systems of up to 9+9 electrons in a spherical geometry using Monte Carlo methods, we show that the ground state can be accurately described by these single-parameter variational states. In addition, for the first time we provide a numerically exact wavefunction for the Halperin-111 state in terms of composite fermions., Comment: 7 pages + 3 figures main text (7 pages + 4 figures supplement)

The DARPA Subterranean Challenge was designed for competitors to develop and deploy teams of autonomous robots to explore difficult unknown underground environments. Categorised in to human-made tunnels, underground urban infrastructure and natural caves, each of these subdomains had many challenging elements for robot perception, locomotion, navigation and autonomy. These included degraded wireless communication, poor visibility due to smoke, narrow passages and doorways, clutter, uneven ground, slippery and loose terrain, stairs, ledges, overhangs, dripping water, and dynamic obstacles that move to block paths among others. In the Final Event of this challenge held in September 2021, the course consisted of all three subdomains. The task was for the robot team to perform a scavenger hunt for a number of pre-defined artefacts within a limited time frame. Only one human supervisor was allowed to communicate with the robots once they were in the course. Points were scored when accurate detections and their locations were communicated back to the scoring server. A total of 8 teams competed in the finals held at the Mega Cavern in Louisville, KY, USA. This article describes the systems deployed by Team CSIRO Data61 that tied for the top score and won second place at the event.

A striking series of experiments have observed superconductivity in Bernal-stacked bilayer graphene (BBG) when the energy bands are flattened by applying an electrical displacement field. Intriguingly, superconductivity manifests only at non-zero magnetic fields, or when spin-orbit coupling is induced in BBG by coupling to a substrate. We present detailed functional renormalization group and random-phase approximation calculations that provide a unified explanation for the superconducting mechanism in both cases. Both calculations yield a purely electronic $p$-wave instability of the Kohn-Luttinger (KL) type. The latter can be enhanced either by magnetic fields or Ising spin-orbit coupling, naturally explaining the behaviour seen in experiments., Comment: 5 pages, 3 figures main text (+2 pages refs+10 pages, 5 figures Supplementary)

Topological quantum chemistry (TQC) is a successful framework for identifying (noninteracting) topological materials. Based on the symmetry eigenvalues of Bloch eigenstates at maximal momenta, which are attainable from first principles calculations, a band structure can either be classified as an atomic limit, in other words adiabatically connected to independent electronic orbitals on the respective crystal lattice, or it is topological. For interacting systems, there is no single-particle band structure and hence, the TQC machinery grinds to a halt. We develop a framework analogous to TQC, but employing $n$-particle Green's function to classify interacting systems. Fundamentally, we define a class of interacting reference states that generalize the notion of atomic limits, which we call Mott atomic limits, and are symmetry protected topological states. Our formalism allows to fully classify these reference states (with $n=2$), which can themselves represent symmetry protected topological states. We present a comprehensive classification of such states in one-dimension and provide numerical results on model systems. With this, we establish Mott atomic limit states as a generalization of the atomic limits to interacting systems.

Step edges of topological crystalline insulators can be viewed as predecessors of higher-order topology, as they embody one-dimensional edge channels embedded in an effective three-dimensional electronic vacuum emanating from the topological crystalline insulator. Using scanning tunneling microscopy and spectroscopy we investigate the behaviour of such edge channels in Pb$_{1-x}$Sn$_{x}$Se under doping. Once the energy position of the step edge is brought close to the Fermi level, we observe the opening of a correlation gap. The experimental results are rationalized in terms of interaction effects which are enhanced since the electronic density is collapsed to a one-dimensional channel. This constitutes a unique system to study how topology and many-body electronic effects intertwine, which we model theoretically through a Hartree-Fock analysis., Comment: 8 pages, 5 figures (main text) and 5 pages, 6 figures (Supplementary Material)

Standard superconductors display a ubiquitous discontinuous jump in the electronic specific heat at the critical superconducting transition temperature. In a growing class of unconventional superconductors, however, a second order parameter component may get stabilized and produce a second heat capacity jump at a lower temperature, typically associated with the spontaneous breaking of time-reversal symmetry. The splitting of the two specific heat discontinuities can be controlled by external perturbations such as chemical substitution, hydrostatic pressure, or uniaxial strain. We develop a theoretical quantitative multi-band framework to determine the ratio of the heat capacity jumps, given the band structure and the order parameter momentum structure. We discuss the conditions of the gap profile which determine the amplitude of the second jump. We apply our formalism to the case of Sr$_2$RuO$_4$, and using the gap functions from a microscopic random phase approximation calculation, we show that recently-proposed accidentally degenerate order parameters may exhibit a strongly suppressed second heat capacity jump. We discuss the origin of this result and consider also the role of spatial inhomogeneity on the specific heat. Our results provide a possible explanation of why a second heat capacity jump has so far evaded experimental detection in Sr$_2$RuO$_4$., Comment: 13+8 pages, 8+8 figures

We study the excitations that emerge upon doping the translationally-invariant correlated insulating states in magic angle twisted bilayer graphene at various integer filling factors $\nu$. We identify parameter regimes where these are associated with skyrmion textures in the spin or pseudospin degrees of freedom, and explore both short-distance pairing effects and the formation of long-range ordered skyrmion crystals. We perform a comprehensive analysis of the pseudospin skyrmions that emerge upon doping insulators at even $\nu$, delineating the regime in parameter space where these are the lowest-energy charged excitations by means of self-consistent Hartree-Fock calculations on the interacting Bistritzer-MacDonald model. We explicitly demonstrate the purely electron-mediated pairing of skyrmions, a key ingredient behind a recent proposal of skyrmion superconductivity. Building upon this, we construct hopping models to extract the effective masses of paired skyrmions, and discuss our findings and their implications for skyrmion superconductivity in relation to experiments, focusing on the dome-shaped dependence of the transition temperature on the twist angle. We also investigate the properties of spin skyrmions about the quantized anomalous Hall insulator at $\nu=+3$. In both cases, we demonstrate the formation of robust spin/pseudospin skyrmion crystals upon doping to a finite density away from integer filling., Comment: 17pp, 12 figures, 4pp references (main text) + 6pp (supplement)

We investigate the full doping and strain-dependent phase diagram of the normal state of magic-angle twisted bilayer graphene (TBG). Using comprehensive Hartree-Fock calculations, we show that at temperatures where superconductivity is absent the global phase structure can be understood based on the competition and coexistence between three types of intertwined orders: a fully symmetric phase, spatially uniform flavor-symmetry-breaking states, and an incommensurate Kekul\'e spiral (IKS) order. For small strain, the IKS phase, recently proposed as a candidate order at all non-zero integer fillings of the moir\'e unit cell, is found to be ubiquitous for non-integer doping as well. We demonstrate that the corresponding electronic compressibility and Fermi surface structure are consistent with the 'cascade' physics and Landau fans observed experimentally. This suggests a unified picture of the phase diagram of TBG in terms of IKS order., Comment: 4+3 pages main text, 6 pages supplement (2 figures main text, 10 figures supplement), accepted version

We introduce a new variational wavefunction for a quantum Hall bilayer at total filling $\nu = 1$, which is based on $s$-wave BCS pairing between composite-fermion electrons in one layer and composite-fermion holes in the other. We compute the overlap of the optimized trial function with the ground state from exact diagonalization calculations of up to 14 electrons in a spherical geometry, and we find excellent agreement over the entire range of values of the ratio between the layer separation and the magnetic length. The trial wavefunction naturally allows for charge imbalance between the layers and provides important insights into how the physics at large interlayer separations crosses over to that at small separations in a fashion analogous to the BEC-BCS crossover., Comment: 18 pages, 8 figures

We study magic angle graphene in the presence of both strain and particle-hole symmetry breaking due to non-local inter-layer tunneling. We perform a self-consistent Hartree-Fock study that incorporates these effects alongside realistic interaction and substrate potentials, and explore a comprehensive set of competing orders including those that break translational symmetry at arbitrary wavevectors. We find that at all non-zero integer fillings very small strains, comparable to those measured in scanning tunneling experiments, stabilize a fundamentally new type of time-reversal symmetric and spatially non-uniform order. This order, which we dub the 'incommensurate Kekul\'e spiral' (IKS) order, spontaneously breaks both the emergent valley-charge conservation and moir\'e translation symmetries, but preserves a modified translation symmetry $\hat{T}'$ -- which simultaneously shifts the spatial coordinates and rotates the $U(1)$ angle which characterizes the spontaneous inter-valley coherence. We discuss the phenomenological and microscopic properties of this order. We argue that our findings are consistent with all experimental observations reported so far, suggesting a unified explanation of the global phase diagram in terms of the IKS order., Comment: 18pp + 6pp references, 8 figures (main text); 16pp, 22 figures (supplement); v3: published version

Heterogeneous teams of robots, leveraging a balance between autonomy and human interaction, bring powerful capabilities to the problem of exploring dangerous, unstructured subterranean environments. Here we describe the solution developed by Team CSIRO Data61, consisting of CSIRO, Emesent and Georgia Tech, during the DARPA Subterranean Challenge. These presented systems were fielded in the Tunnel Circuit in August 2019, the Urban Circuit in February 2020, and in our own Cave event, conducted in September 2020. A unique capability of the fielded team is the homogeneous sensing of the platforms utilised, which is leveraged to obtain a decentralised multi-agent SLAM solution on each platform (both ground agents and UAVs) using peer-to-peer communications. This enabled a shift in focus from constructing a pervasive communications network to relying on multi-agent autonomy, motivated by experiences in early circuit events. These experiences also showed the surprising capability of rugged tracked platforms for challenging terrain, which in turn led to the heterogeneous team structure based on a BIA5 OzBot Titan ground robot and an Emesent Hovermap UAV, supplemented by smaller tracked or legged ground robots. The ground agents use a common CatPack perception module, which allowed reuse of the perception and autonomy stack across all ground agents with minimal adaptation.

The long-standing quest to determine the superconducting order of Sr$_2$RuO$_4$ (SRO) has received renewed attention after recent nuclear magnetic resonance (NMR) Knight shift experiments have cast doubt on the possibility of spin-triplet pairing in the superconducting state. As a putative solution, encompassing a body of experiments conducted over the years, a $d+ig$-wave order parameter caused by an accidental near-degeneracy has been suggested [S. A. Kivelson et al., npj Quantum Materials $\bf{5}$, 43 (2020)]. Here we develop a general Ginzburg--Landau theory for multiband superconductors. We apply the theory to SRO and predict the relative size of the order parameter components. The heat capacity jump expected at the onset of the second order parameter component is found to be above the current threshold deduced by the experimental absence of a second jump. Our results tightly restrict theories of $d+ig$ order, and other candidates caused by a near-degeneracy, in SRO. We discuss possible solutions to the problem., Comment: 13 pages, 8 figures

We consider magic-angle twisted bilayer graphene (TBG) at filling $\nu=+3$, where experiments have observed a robust quantized anomalous Hall effect. This has been attributed to the formation of a valley- and spin-polarized Chern insulating ground state that spontaneously breaks time-reversal symmetry, and is stabilized by a hexagonal boron nitride (hBN) substrate. We identify three different types of domain wall, and study their properties and energetic selection mechanisms via theoretical arguments and Hartree-Fock calculations adapted to deal with inhomogeneous moir\'e systems. We comment on the implications of these results for transport and scanning probe experiments., Comment: 19 pages, 13 figures

We apply the semi-classical quantum Boltzmann formalism for the computation of transport properties to multilayer graphene. We compute the electrical conductivity as well as the thermal conductivity and thermopower for Bernal-stacked multilayers with an even number of layers. We show that the window for hydrodynamic transport in multilayer graphene is similar to the case of bilayer graphene. We introduce a simple hydrodynamic model which we dub the multi-fluid model and which can be used to reproduce the results for the electrical conductivity and thermopower from the quantum Boltzmann equation., Comment: 15 pages, 8 figures