蒙特卡罗电子双缝衍射计算机模拟.docVIP

基于蒙特卡罗的电子双缝衍射的计算机模拟 电子双缝衍射是微观粒子具有波动性的重要证明实验，开始是作为假想实验而提出的， 1988年才由Tonomura等人做出了该实验[1]。鉴于一般教学仪器不具备进行上述实验的条件，本文依据电子衍射的几率密度函数，运用蒙特卡罗随机模拟方法，借助计算机的数据可视化技术、绘图技术，构建电子双缝衍射的动态随机过程，清晰地演示出电子衍射的全过程。 1电子双缝衍射的随机模拟 电子衍射是大量电子随机运动的统计结果，可运用蒙特卡罗方法对其进行模拟。运用蒙特卡罗方法来处理本文的问题，首先必须根据要处理问题的规律，建立一个概率模型，然后进行随机抽样试验，从而得出一组按已知分布的随机数序列。最后依据这一随机数序列，借助计算机程序设计语言或图形软件，就可实现电子双缝衍射动态随机过程的模拟。 概率模型的构建 图1是电子双缝衍射的示意图。衍射屏位于x′o′y′平面，观测屏位于xoy平面，缝S1和S2的宽度均为a，两缝中心间距为(a+b)，衍射屏到观测屏的距离为D。设动量为p、能量为E的自由电子沿Z轴正方向入射到双缝，其波函数为: (1) 设在t=0时刻，波前到达双缝处(z=0)，由式(1)可知，在双缝处波函数为一常数。为简单起见，设在衍射屏上的波函数为: (2) 根据费曼的路径积分理论，电子经过双缝而在t时刻到达观测屏上P点的衍射波函数为: (3) 式中，。所以，电子经过双缝在观测屏上P点出现的几率密度为: (4) 设D (a+b)，则有: (5) 代入式(3.4)有: (6) 其中，。 数据的生成 根据上述电子在观测屏上出现的几率密度函数进行随机抽样，便得到按此几率密度函数分布的随机数序列。在蒙特卡罗方法中，有多种方法可实现按已知分布的随机抽样。根据本文要处理的问题的特点，采用Von Neumann的舍选法。 舍选法的具体做法是：(1)计算机在一定的范围内随机地选取观测屏一坐标点，并计算。其中，是几率密度函数在点的值，是几率密度函数的最大值；(2)计算机产生一个0至1之间均匀分布的随机数M；(3)将H与M进行比较，若，则选取该点，若，则舍去该点，并重复(1)至(3)，重新选择坐标点。通过以上方法，便可以得到一个按分布的随机数序列。这个随机数序列便是我们实现电子衍射随机运动过程的依据。 电子双缝衍射动态随机运动过程的实现 利用上述生成的按分布的随机数序列，借助计算机程序设计语言或图形软件，就可以实现电子双缝衍射动态随机过程的模拟。本文采用编写Matlab程序的方法来实现。在编写Matlab程序时，用Matlab的rand函数来产生均匀分布的随机数，用plot函数来实现电子衍射点的显示。 2电子双缝衍射的Matlab程序 取，加速电子的电压为。根据可算出入射电子的波长。 从式(6)可看出，电子在观测屏上出现的最大几率密度为，即，所以： (7) 根据以上讨论编写的电子双缝衍的Matlab程序如附录. 3电子双缝衍射图 运行上述程序便会得到电子双缝衍射动态随机过程的演示，图2、图3、图4分别是电子数为100(cx1)、1000(cx2)、2000(cx3)时的衍射图。从中可看出，当电子数目较少时(如图2)，衍射特征不明显，电子表现出的是粒子性；但随着电子数目的增多(如图3、4)，衍射特征逐渐明显，电子表现出了波动性。将衍射结果与根据量子力学原理得出的衍射几率密度分布函数进行比较。可以看出，衍射图中所反映出的电子点的分布，与几率密度分布函数曲线所反映的几率密度大小分布完全一致。图5(cx4)电子数为20 000时的衍射图样与衍射几率密度分布函数曲线对比图（图5中markerSize取2）。 附录: CX1: clear h=6.62559e-34; m=9.10908e-31; e=1.6021e-19; U=1000; a=2e-7; b=1e-6; D=0.25; rd=h/sqrt(2*m*e*U); %计算电子波长rd A=(pi*a)/(rd*D); B=pi*(a+b)/(rd*D); %计算A和B axis([-5e-5,5e-5,-4e-5,4e-5]) %设置坐标轴标度范围 set(gca,color,[0.122,0.012,0.62]) %设置坐标面背景颜色 title(电子双缝衍射动态随机过程演示,fontsize,16,color,k) %设置标题 i=1; while i=100 %用循环结构控制电子的数目 x=5e-5*(2*rand-1); %产生5e-5至-5e-5之间均匀分布的随机数赋给