• Support PF! Buy your school textbooks, materials and every day products Here!

Split operator with FFT

  • Thread starter coelurus
  • Start date
1
0
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:

http://coelurus.thorntwig.se/data/pics/splitop.jpg" [Broken]
http://coelurus.thorntwig.se/data/tjo.m" [Broken]

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 by a moderator:

Related Threads for: Split operator with FFT

  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
7
Views
5K
  • Last Post
Replies
4
Views
1K
  • Last Post
Replies
2
Views
905
  • Last Post
Replies
2
Views
3K
  • Last Post
Replies
4
Views
7K
  • Last Post
Replies
2
Views
3K
Top