We perform an extensive exact diagonalization study of interaction driven insulators in spin-and valleypolarized moiré flat bands of twisted bilayer graphene aligned with its hexagonal boron nitride substrate. In addition to previously reported fractional Chern insulator phases, we provide compelling evidence for competing charge-density-wave phases at multiple fractional fillings of a realistic single-band model. A thorough analysis at different interlayer hopping parameters, motivated by experimental variability, and the role of kinetic energy at various Coulomb interaction strengths highlight the competition between these phases. The interplay of the single-particle and the interaction induced hole dispersion with the inherent Berry curvature of the Chern bands is intuitively understood to be the driving mechanism for the ground-state selection. The resulting phase diagram features remarkable agreement with experimental findings in a related moiré heterostructure and affirms the relevance of our results beyond the scope of graphene based materials.

Twisted double- and mono-bilayer graphene are graphene-based moiré materials hosting strongly correlated fermions in a gate-tunable conduction band with a topologically non-trivial character. Using unbiased exact diagonalization complemented by unrestricted Hartree-Fock calculations, we find that the strong electron-electron interactions lead to a non-coplanar magnetic state, which has the same symmetries as the tetrahedral antiferromagnet on the triangular lattice and can be thought of as a skyrmion lattice commensurate with the moiré scale, competing with a set of ferromagnetic, topological charge density waves featuring an approximate emergent O(3) symmetry, “rotating” the different charge density wave states into each other. Direct comparison with exact diagonalization reveals that the ordered phases are accurately described within the unrestricted Hartree-Fock approximation. Exhibiting a finite charge gap and Chern number |C|=1, the formation of charge density wave order which is intimately connected to a skyrmion lattice phase is consistent with recent experiments on these systems.