量子圍栏的理论计算及计算机模拟.docVIP

原子物理读书报告 ——量子围栏的理论计算及计算机模拟 原子物理上到第三章量子力学导论的时候，老师在课件中展示了一幅由STM扫描所得到的量子围栏的图片，课后我查阅了有关书藉，并进行了相应的思考，发现可以用我们现有的理论基础，就可以定性地讨论这个问题，并用Matlab软件进行了计算机模拟。 理论推导 理论推导可将铁原子的围栏视为无限深的势阱，就可以用类似于求解一维无限深势阱的方法来求解二维无限深势阱的问题。 表面电子的波函数在铁原子处反射，形成图中所示的驻波，类似于一维的问题，可视围栏处为波节，即该处的波函数（概率幅）为0，问题就化为以围栏为边界条件的定态薛定谔方程。 可见，方程化为亥姆霍兹方程，在极坐标系下用分离变量法求解该方程。 计算机模拟 方程解的化简： 参考其它书藉知，在电子表面态的费米面上主要有 | n,l = |5,0 、|4,2 、|2,7三个态，即总的概率函数可由这三个态的波函数模方适当地组合起来，即可以很好地将实验曲线很好地拟合出来。而这三个波函数的模方是按以下条件归一化的： 对于这三个贝塞尔函数中的 ，即l阶贝塞尔函数的第n个零点（原点不算），通过传统的方面求这些 比较困难，在这个问题中可以用Matlab求得它的近似值： 取x为0到20，y为x的l阶（以0阶为例）贝塞尔函数，作图，观察零点的大体位置： x=0:0.01:20; y=besselj(0,x); plot(x,y) 结果如图： 可见，0阶贝塞尔函数的第5个零点在x=14到16之间，下面求x50 x50=fzero(besselj(0,x),[14,16]) x50 = 14.9309 计算该贝塞尔函数值矩阵(备用)： r=0:0.01:1; u50=(besselj(0,x50.*r)).^2; 计算系数a50,由 可得系数，由于积分比较复杂，可用Matlab得出数值积分（取R=1）： a50=1/(2*pi*trapz(r,u50.*r)) a50 = 7.4642 同理，用同样的办法可得： x42 = 14.7960 a42 = 7.4617 x27 = 14.8213 a27 = 8.3890 于是，可得到三个函数模方以各自系数为权的叠加，即为从一维观察的情况： u=a50.*u50+a42.*u42+a27.*u27; plot(r,u) 量子围栏中波函数的三维再现： 三维情况中，取rx,ry为两坐标，范围分别为（-1，1），则r= ,于是： rx=-1:0.01:1; ry=rx; [x,y]=meshgrid(rx,ry); u50=(besselj(0,x50.*sqrt(x.^2+y.^2))).^2; u42=(besselj(2,x42.*sqrt(x.^2+y.^2))).^2; u27=(besselj(7,x27.*sqrt(x.^2+y.^2))).^2; u=a50.*u50+a42.*u42+a27.*u27; surf(x,y,u) 所得的结果如下图： 计算机模拟图 实际照片 分析 在这个问题中，我近似认为铁原子处的无限深的势阱，而事实并非无限深，所以计算有一些误差，而且由于势阱不是无限深，从而会有隧道效应的存在，实际图中围栏外部的波函数图像可能就是隧道效应的缘故。 关于为什么各态的波函数归一化条件为 ，还是不清楚。 参考书藉 《原子物理学》，杨福家著，高等教育出版社2000年版。 《数学物理方法》，沈 著，同济大学出版社2002年版。 《量子物理》，赵凯华，罗蔚茵著，高等教育出版社2001年版。 《MATLAB在大学物理在的应用》，鲍鸿吉著，同济大学物理系。 《理论力学计算机模拟》，彭芳麟等著，清华大学出版社2002年版。  《量子物理》，第139页到141页，高等教育出版社。 应用物理 013249 金赞 5 专业：应用物理 学号：013249 姓名：金赞