Carmen Rubio-Verdú, Simon Turkel, Yuan Song, Lennart Klebl, Rhine Samajdar, Mathias S. Scheurer, Jörn W. F. Venderbos, Kenji Watanabe, Takashi Taniguchi, Héctor Ochoa, Lede Xian, Dante M. Kennes, Rafael M. Fernandes, Ángel Rubio, Abhay N. Pasupathy, Department of Energy (US), Air Force Office of Scientific Research (US), European Commission, European Research Council, Flatiron Institute, German Research Foundation, Max Planck Society, National Science Foundation (US), RWTH Aachen University, Japan Society for the Promotion of Science, Ministry of Education, Culture, Sports, Science and Technology (Japan), and Simons Foundation
Graphene moiré superlattices display electronic flat bands. At integer fillings of these flat bands, energy gaps due to strong electron–electron interactions are generally observed. However, the presence of other correlation-driven phases in twisted graphitic systems at non-integer fillings is unclear. Here, we report the existence of three-fold rotational (C3) symmetry breaking in twisted double bilayer graphene. Using spectroscopic imaging over large and uniform areas to characterize the direction and degree of C3 symmetry breaking, we find it to be prominent only at energies corresponding to the flat bands and nearly absent in the remote bands. We demonstrate that the magnitude of the rotational symmetry breaking does not depend on the degree of the heterostrain or the displacement field, being instead a manifestation of an interaction-driven electronic nematic phase. We show that the nematic phase is a primary order that arises from the normal metal state over a wide range of doping away from charge neutrality. Our modelling suggests that the nematic instability is not associated with the local scale of the graphene lattice, but is an emergent phenomenon at the scale of the moiré lattice., S.T. and A.N.P. acknowledge funding from Programmable Quantum Materials, an Energy Frontier Research Center funded by the US Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), under award no. DE-SC0019443. STM equipment support (A.N.P.) and 2D sample synthesis (Y.S.) were provided by the Air Force Office of Scientific Research via grant no. FA9550-16-1-0601. C.R.-V. acknowledges funding from the European Union Horizon 2020 Research and Innovation Programme under the Marie Skłodowska-Curie grant agreement no. 844271. A.R. acknowledges funding by the European Research Council (ERC-2015-AdG-694097), Grupos Consolidados (IT1249-19) and the Flatiron Institute, a division of the Simons Foundation. L.K., D.M.K. and A.R. acknowledge funding by the Deutsche Forschungsgemeinschaft (DFG) under Germany’s Excellence Strategy-Cluster of Excellence Matter and Light for Quantum Computing (ML4Q) EXC 2004/1-390534769 and Advanced Imaging of Matter (AIM) EXC 2056−390715994 and funding by the Deutsche Forschungsgemeinschaft (DFG) under RTG 1995, within the Priority Program SPP 2244 ‘2DMP’ and GRK 2247. A.R. acknowledges support by the Max Planck Institute-New York City Center for Non-Equilibrium Quantum Phenomena. H.O. is supported by the NSF MRSEC programme grant no. DMR-1420634. Tight-binding and fRG simulations were performed with computing resources granted by RWTH Aachen University under projects rwth0496 and rwth0589. R.S. and M.S.S. acknowledge support from the National Science Foundation under grant no. DMR-2002850. R.M.F. was supported by the DOE-BES under award no. DE-SC0020045. K.W. and T.T. acknowledge support from the Elemental Strategy Initiative conducted by the MEXT, Japan (grant no. JPMXP0112101001), JSPS KAKENHI (grant no. JP20H00354) and the CREST (grant no. JPMJCR15F3) JST.