FFT-split-operator code for solving the Dirac equation in 2+1 dimensions

Published: 01-06-2008| Version 1 | DOI: 10.17632/43v3vvkwwf.1
Guido R. Mocken,
Christoph H. Keitel


Abstract The main part of the code presented in this work represents an implementation of the split-operator method [J.A. Fleck, J.R. Morris, M.D. Feit, Appl. Phys. 10 (1976) 129-160; R. Heather, Comput. Phys. Comm. 63 (1991) 446] for calculating the time-evolution of Dirac wave functions. It allows to study the dynamics of electronic Dirac wave packets under the influence of any number of laser pulses and its interaction with any number of charged ion potentials. The initial wave function can be eith... Title of program: Dirac++ or (abbreviated) d++ Catalogue Id: AEAS_v1_0 Nature of problem The relativistic time evolution of wave functions according to the Dirac equation is a challenging numerical task. Especially for an electron in the presence of high intensity laser beams and/or highly charged ions, this type of problem is of considerable interest to atomic physicists. Versions of this program held in the CPC repository in Mendeley Data AEAS_v1_0; Dirac++ or (abbreviated) d++; 10.1016/j.cpc.2008.01.042 This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2019)