Monte Carlo Study of the Semimetal-Insulator Phase Transition in Monolayer Graphene with a Realistic Interelectron Interaction Potential

We report on the results of the first-principles numerical study of spontaneous breaking of chiral (sublattice) symmetry in suspended monolayer graphene due to electrostatic interaction, which takes into account the screening of Coulomb potential by electrons on sigma orbitals. In contrast to the results of previous numerical simulations with unscreened potential, we find that suspended graphene is in the conducting phase with unbroken chiral symmetry. This finding is in agreement with recent experimental results by the Manchester group [D. C. Elias et al., Nat. Phys. 7, 701 (2011); A. S. Mayorov et al., Nano Lett. 12, 4629 (2012)]. Further, by artificially increasing the interaction strength, we demonstrate that suspended graphene is quite close to the phase transition associated with spontaneous chiral symmetry breaking, which suggests that fluctuations of chirality and nonperturbative effects might still be quite important.