QuantumFDTD -- A computational framework for the relativistic Schrödinger equation
We extend the publicly available quantumfdtd code. It was originally intended for solving the time-independent three-dimensional Schrödinger equation via the finite-difference time-domain (FDTD) method and for extracting the ground, first, and second excited states. We (a) include the case of the relativistic Schrödinger equation and (b) add two optimized FFT-based kinetic energy terms for the non-relativistic case. All the three new kinetic terms are computed using Fast Fourier Transform (FFT). We release the resulting code as version 3 of quantumfdtd. Finally, the code now supports arbitrary external file-based potentials and the option to project out distinct parity eigenstates from the solutions. Our goal is quark models used for phenomenological descriptions of QCD bound states, described by the three-dimensional Schrödinger equation. However, we target any field where solving either the non-relativistic or the relativistic three-dimensional Schrödinger equation is required.