“…A multitude of different implementations exist, depending on the spatial discretization and the time propagation, each having their strengths and weaknesses. Moreover, the FDTD method has also been applied to various related equations such as the non-linear Schrödinger equation [20][21][22][23][24][25][26], the Gross-Pitaevskii equations [27], and the stationary Schrödinger equation by performing a Wick rotation [28][29][30][31].…”