Single-cone finite difference scheme for the (2+1)D Dirac von Neumann equation
An explicit finite difference scheme is presented for the von Neumann equation for (2+1)D Dirac fermions. It is founded upon a staggered space–time grid which ensures a single-cone energy dispersion and performs the time-derivative in one sweep using a three-step leap-frog procedure. It enables a space–time-resolved numerical treatment of the mixed-state dynamics of Dirac fermions within the effective single-particle density matrix formalism. Energy-momentum dispersion, stability and convergence properties are derived. Elementary numerical tests to demonstrate stability properties use parameters which pertain to topological insulator surface states. A method for the simulation of charge injection from an electric contact is presented and tested numerically. Potential extensions of the scheme to a Dirac–Lindblad equation, real-space–time Green's function formulations, and higher-order finite-difference schemes are discussed.
- OSTI ID:
- 22701624
- Journal Information:
- Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 348; ISSN JCTPAH; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics