Anderson Transition in Disordered Graphene

We use the regularized kernel polynomial method (RKPM) to numerically study the effect disorder on a single layer of graphene. This accurate numerical method enables us to study very large lattices with millions of sites, and hence is almost free of finite size errors. Within this approach, both weak and strong disorder regimes are handled on the same footing. We study the tight-binding model with on-site disorder, on the honeycomb lattice. We find that in the weak disorder regime, the Dirac fermions remain extended and their velocities decrease as the disorder strength is increased. However, if the disorder is strong enough, there will be a {\em mobility edge} separating {\em localized states around the Fermi point}, from the remaining extended states. This is in contrast to the scaling theory of localization which predicts that all states are localized in two-dimensions (2D).Comment: 4 page