Quantum wave packet propagation, how to use it in FFT? So I used the split step method on the Schrodinger equation and have produced the following equation:
$\Psi(x,t+dt)=F^{-1} \left\{ e^{-i\frac{\hbar^2k^2}{2m}\frac{dt}{\hbar}}F\left\{e^{-iV(x)\frac{dt}{\hbar}}\Psi(x,t) \right\}\right\}$
Which when scaled to dimensionless the time evolution step can be written as:
$\Psi(x,t+dt)=F^{-1}\left\{e^{-ik^2dt}F\left\{e^{-iV(x)dt}\Psi(x,t)\right\}\right\}$
F above represents the Fourier transform operator acting on the equations.
Now my problem is that I don't know how to implement the above in Fast Fourier Transform.
For example: If I were to take the exponential factor with V(x), do I multiply -iV(x) by dt?
Not sure I understand what is going on here....
Thanks
 A: So, given the dimensionless Hamiltonian
$$
\hat{\cal{H}} = \frac{1}{2}\hat{p}^2 + V(\hat{x}),
$$
the wave function evolves as
$$
\vert\Psi(t+dt)\rangle = e^{-i\hat{\cal{H}}dt}\vert\Psi(t)\rangle = e^{-\frac{1}{2}i\hat{p}^2 dt}e^{-iV(\hat{x})dt}e^{O(dt^2)}\vert\Psi(t)\rangle.
$$
The idea is to apply the position-space part of the evolution operator ($e^{-iV(\hat{x})dt}$) to the position-space wave function $\Psi(x, t)$, where it is just a multiplication by $e^{-iV(x)dt}$, and then to apply the momentum-space part of the evolution operator ($e^{-\frac{1}{2}i\hat{p}^2 dt}$) to the momentum-space wave function $\tilde{\Psi}(k, t+dt)$, where it is just a multiplication by $e^{-\frac{1}{2}ik^2 dt}$.  To evolve the wave function in this approximation (errors are introduced because terms of order $dt^2$ arising from the noncommutativity of $\hat{x}$ and $\hat{p}$ are being omitted), you would typically choose a small value for $dt$ and a discretization for $x$ and $k$, and then repeatedly apply the four operations shown in your expression above: multiply by the potential energy term, apply the FFT to convert to momentum space, multiply by the kinetic energy term, and apply the inverse FFT to convert back to position space.
