Magnetic and superconducting correlation in monolayer and twisted bilayer graphene

Using exact quantum Monte Carlo method, we identify the phase diagram of the half filled, the lightly doped and heavily doped graphene. At half filling, the system is driven to a Mott insulator with antiferromagnetic long range order by increasing interaction, and a transition from a d+id pairing to a p+ip pairing is revealed, depends on the next-nearest hoping and the electronic fillings. We also examine the recent novel electronic states seen in magic-angle graphene superlattices. We reveal that an antiferromagnetically ordered Mott insulator emerges beyond a critical Uc at half filling, and with a small doping, the pairing with d+id symmetry dominates over other pairings at low temperature. The effective d+id pairing interaction strongly increase as the on-site Coulomb interaction increases, indicating that the superconductivity is driven by electron-electron correlation. Our non-biased numerical results demonstrate that the twisted bilayer graphene share the similar superconducting mechanism of high temperature superconductors, which is a new and idea platform for further investigating the strongly correlated phenomena.