Ground state superconducting pair correlations in twisted bilayer graphene

Motivated by the recent novel electronic features extracted from the magic-angle graphene superlattices, we studied the ground state superconducting pairing correlations within the Hubbard model on a twisted bilayer honeycomb lattice. Using Constrained-Path Quantum Monte Carlo method, we found that the [Formula: see text] pairing correlation dominates over other pairing patterns among various electron fillings and interaction strengths, and the effective pairing interaction was enhanced as the on-site Coulomb interaction increased. We further examined the effect of the nearest neighbor interaction [Formula: see text], and the effective pairing interaction with [Formula: see text] pairing symmetry was also enhanced by either a repulsive or attractive interaction. Our intensive numerical results confirm the interaction driven superconductivity with a dominant [Formula: see text] pairing symmetry in twisted bilayer graphene.