1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Split operator with FFT

  1. Jan 8, 2007 #1
    Hi! This is my first post on these forums.

    I'm having some problems with the split operator + FFT algorithm to solve the Schrödinger equation in time. A real Gauss curve in a zero potential environment should simply flatten out, but I get two peaks as the following Matlab run shows:

    Little pic.
    Matlab code.

    The real problem I assume is that I have no experience at all in using FFT for numerical simulations similar to this method and I can't find any good references on it either.
    Any hints, ideas or pointers would be much appreciated!

    EDIT: I found a "solution", but I am not sure how to interpret it yet. I had a look in the WavePacket package and it seems that one has to shift both position and momentum space (in Matlab, one would use "fftshift") before and after each FFT. If anybody knows a rigid answer to that I would be overjoyed :) So yes, it works, but I'd like to see why...
    Last edited: Jan 8, 2007
  2. jcsd
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Can you help with the solution or looking for help too?

Similar Discussions: Split operator with FFT
  1. Heat engine operations (Replies: 0)