A convergent 2D finite-difference scheme for the Dirac-Poisson system and the simulation of graphene
KAUST DepartmentComputer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
Applied Mathematics and Computational Science Program
MetadataShow full item record
AbstractWe present a convergent finite-difference scheme of second order in both space and time for the 2D electromagnetic Dirac equation. We apply this method in the self-consistent Dirac-Poisson system to the simulation of graphene. The model is justified for low energies, where the particles have wave vectors sufficiently close to the Dirac points. In particular, we demonstrate that our method can be used to calculate solutions of the Dirac-Poisson system where potentials act as beam splitters or Veselago lenses. © 2013 Elsevier Inc.
SponsorsThe authors acknowledge support from King Abdullah University of Science and Technology (KAUST) Award Number KUK-I1-007-43 and from the WWTF (Vienna Science and Technology Fund) Project Number MA09-028.
JournalJournal of Computational Physics