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.
| 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 15 2014 |
Fingerprint
Keywords
- Beam splitter
- Dirac equation
- Dirac-Poisson system
- Finite differences
- Graphene
- Veselago lens
ASJC Scopus subject areas
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- Physics and Astronomy(all)
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics
Cite this
Brinkman, D., Heitzinger, C., & Markowich, P. A. (2014). A convergent 2D finite-difference scheme for the Dirac-Poisson system and the simulation of graphene. Journal of Computational Physics, 257(PA), 318-332. https://doi.org/10.1016/j.jcp.2013.09.052