Split operator method
You can also download a PDF version of this article.
The split operator method (SOM) is a method to avoid manually solving Schrödinger's equation and obtaining an approximate solution for an arbitrary potential of the form . Later on we shall see what split-operator does mean. For simplicity, our discussion will be regarding the one-dimensional problem. Given the standard form of the time-dependent Schrödinger's equation
its formal solution is as follows:
Let us consider the following hamiltonian function
Since , there is no reason to prefer one of the two following options to set the order in which the hamiltonian momentum-dependent and position-dependent part act over the system ket:
In fact, none of the two would be correct.
The SOM consists on somehow setting up an approximation to the exponential. This approximation is
Thus, we have splitted the hamiltonian operator in a certain way to symmetrize the exponentials. This approximation is valid for a sufficiently small time step.
Now we are going to switch to momentum representation to see how this works for the wave packet. We do not choose position representation because that would require an additional step:
As we stated before, the time step must be sufficiently small. Apparently, we cannot solve Schrödinger's equation for an arbitrary value of . To solve this problem we are going to make the following substitution:
Now we can make as small as we need to match convergence requirements.
Given the closure, eigenvalues and dot product relations for the position and momentum basis sets
and keeping in mind that , it is true that
Hence, we can, in a straightforward manner, apply the closure relation twice, once in position and once in momentum, to get to the following result:
which is a double Fourier transform, once forward and once backward. To obtain we just need to transform again (the additional step I was referring to before). Thus, under the SOM approximation and given that we know the state of the system at a certain time, we can solve Schrödinger's equation by simply using Fourier transforms whatever the shape of the potential and obtain the wave function later.
I made the next video with data obtained from a SOM routine I programmed. Both the wave packet and potential barrier are gaussian. In the left side you can see position representation and in the right one you have momentum representation: