Abstract
We 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.
Original language | English (US) |
---|---|
Pages (from-to) | 318-332 |
Number of pages | 15 |
Journal | Journal of Computational Physics |
Volume | 257 |
Issue number | PA |
DOIs | |
State | Published - Jan 2014 |
Bibliographical note
KAUST Repository Item: Exported on 2020-10-01Acknowledged KAUST grant number(s): KUK-I1-007-43
Acknowledgements: The 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.
ASJC Scopus subject areas
- Physics and Astronomy (miscellaneous)
- Computer Science Applications