夏力钢
Numerical solution to Schrodinger Equation
2019-5-24 22:23
阅读:2654


Schrodinger equation plays an important role in Quantum Mechanics. 


{\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi (\mathbf {r} ,t)=\left[{\frac {-\hbar ^{2}}{2m}}\nabla ^{2}+V(\mathbf {r} ,t)\right]\Psi (\mathbf {r} ,t)}


Despite the mathematical difficulties for analytic solutions, we can always try to obtain the numberical solution. In the following paragraphs, I am going to show the numerical simulation of the tunneling effect and the double-slit inteference based on the time-dependent Schrodinger equation.


The reference is this one.

https://web.pa.msu.edu/people/duxbury/courses/phy480/SchrodingerDynamics.pdf


The python script can be found in my github.

https://github.com/xialigang/schrodinger



  1. Tunneling effect

    This effect has been used by Gamow to explain the alpha decay in nuclear physics. It is also the principle behind the scanning-tunneling microscope. It says that a particle can pass through a potential in quantum mechanics, which is forbidden in classic physics.

    Let us look at the simulation at different times. We start with a Gaussian wave package. The barrier is put at x=50. Clearly, we can see a part of the wave package penetrating the barrier at late times.

    Cs_psix_t0.pngCs_psix_t3000.png

    Cs_psix_t5000.pngCs_psix_t7000.png

  2. Double-slit interference experiment

    The two-slit apparatus is put at x=50. From the second plot, we can clearly see two slits. From the last plot, we can see the interference pattern on both sides of the two-slit barrier. This is because we have scattered (forward) and reflected (backward) waves.

    Cs_psix_t0.pngCs_psix_t4000.pngCs_psix_t7000.pngCs_psix_t11000.png


It is interesting absolutely, helping us to see how a wave envolve with time. You are welcome to use the codes in my github ( https://github.com/xialigang/schrodinger ). You only need to install PYTHON and ROOT, which are used in the field of experimental high-energy physics (and other fields ofcourse).


I would say the numerical simulation algorithm is very simple. You can try by youself (one day is enough). The thing I did not tell you is that I spend a lot time adjusting the parameters and spend a lot of CPU time to produce these beautiful plots.

转载本文请联系原作者获取授权,同时请注明本文来自夏力钢科学网博客。

链接地址:https://wap.sciencenet.cn/blog-3401446-1180986.html?mobile=1

收藏

分享到:

当前推荐数:0
推荐到博客首页
网友评论0 条评论
确定删除指定的回复吗?
确定删除本博文吗?