1
$\begingroup$

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

  • 0
    Not sure what your problem is?2011-01-30

1 Answers 1

2

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.

  • 0
    @chutsu: Your FFT algorithm is a discretization of the original problem. The time dependence is not really very explicit anymore, but it is there. First as mjqxxxx points out, there is the time interval $dt$, but more importantly, the wave function at time-step $t+dt$ is given as a function of the wave function at time step $t$ with some extra operations performed on it. So, if you want the full time development, you have to iterate and repeat the same procedure over again but with $\psi(t+dt)$ in the input to obtain $\psi(t+2dt)$ in the output, etc...2011-02-03