Best way to solve Schrodinger's wave equation numerically.


by Slide rule
Tags: maple, numerical pde, pde, pde wave, quantum mechanics
Slide rule
Slide rule is offline
#1
Feb18-13, 06:58 PM
P: 3
I have been trying to research the best way to solve the Schrodinger wave equation numerically so that I can plot and animate it in Maple. I'd also like to animate as it is affected by a potential. I have been trying for weeks to do this and I don't feel any closer than when I started. I have looked at finite difference method but I get so far and don't know what to do next.

Any help would be greatly appreciated.

The sort of thing I'm looking for is like in this presentation on youtube, especially at 11s leading onto something like the animation at 13s.

http://m.youtube.com/watch?hl=en-GB&client=mv-google&gl=GB&v=Xj9PdeY64rA&fulldescription=1

Thanks you very much
Phys.Org News Partner Physics news on Phys.org
Physicists design quantum switches which can be activated by single photons
'Dressed' laser aimed at clouds may be key to inducing rain, lightning
Higher-order nonlinear optical processes observed using the SACLA X-ray free-electron laser
DrClaude
DrClaude is offline
#2
Feb19-13, 04:34 AM
Sci Advisor
PF Gold
DrClaude's Avatar
P: 1,111
Can you give more detail? What is the Hamiltonian? In how many dimensions?

And there is no link to the YouTube video.
Slide rule
Slide rule is offline
#3
Feb19-13, 04:49 AM
P: 3
Hi, sorry for not posting the link, I was very tired when making this post, definitely an oversight on my part. I will be able to post the link in just over an hour.

With regards to dimensions in would only be in 1 dimension along x. Regarding the hamiltonian I am trying to solve the equation as
i x hbar x diff(psi, t) = -(hbar^2)/2m x diff(psi, x$2) + V(x) x psi
where psi = psi(x, t).

Thank you for replying.

DrClaude
DrClaude is offline
#4
Feb19-13, 05:43 AM
Sci Advisor
PF Gold
DrClaude's Avatar
P: 1,111

Best way to solve Schrodinger's wave equation numerically.


Since your potential is time independent, the fastest way to solve your problem is to first find the eigenfunctions of the Hamiltonian
[tex]
H \phi_i = E_i \phi_i
[/tex]
To do this, discretize space and write the Hamiltonian as a matrix. As you said, you can use a finite difference approximation for the momentum operator.

Once you have the [itex]\phi_i[/itex], find the initial coefficients of your wave function in this basis,
[tex]
\psi(x,t=0) = \sum_i c_i \phi_i(x)
[/tex]
by calculating
[tex]
c_i = \int \phi_i^*(x) \psi(x,t=0) dx
[/tex]

Then, the wave function at any time time is simply given by
[tex]
\psi(x,t) = \sum_i c_i \phi_i(x) \exp(-i E_i t / \hbar)
[/tex]
By advancing [itex]t[/itex] and refreshing the plot, you will get your animation. I have no idea how to do this in Maple
Bill_K
Bill_K is offline
#5
Feb19-13, 05:17 PM
Sci Advisor
Thanks
Bill_K's Avatar
P: 3,853
I have been trying to research the best way to solve the Schrodinger wave equation numerically so that I can plot and animate it in Maple. I'd also like to animate as it is affected by a potential. I have been trying for weeks to do this and I don't feel any closer than when I started. I have looked at finite difference method but I get so far and don't know what to do next.
Slide rule, Why not ask this question over on the Diff Eq forum.


Register to reply

Related Discussions
To solve an equation numerically using mathematica. Math & Science Software 6
who has the programma used to solve the coupled Dirac equation numerically? Quantum Physics 0
how to solve this nonlinear differential equation numerically? Differential Equations 10
solve numerically the self-consistent Bogoliubov-De Gennes equation Atomic, Solid State, Comp. Physics 0
solve complex equation numerically Programming & Computer Science 7