Kekulé valence bond order in an extended Hubbard model on the honeycomb lattice with possible applications to twisted bilayer graphene

Abstract

Motivated by recent experimental findings of correlating insulator and superconductor behavior in twisted bilayer graphene, the authors study a cluster charge interaction model on a honeycomb lattice. The cluster charge is defined on the dual lattice of honeycomb to mimic charge centers forming a triangular lattice found in experiments. Based on unbiased sign-problem-free quantum Monte Carlo simulations, the authors determine the phase diagram for this model. It features a fermionic quantum critical point with chiral XY universality and a Kekul\'e valence bond solid phase at intermediate interaction, which may be related to the correlating insulator phase found in experiments. This is an example of the rich phases found in a simple model without artificial design.