一维非对称方势阱中的粒子研讨.docVIP

一维非对称方势阱中的粒子 实验任务 求解实物粒子在一维非对称方势阱中的分布。 实验目的 1、加强对一维势阱理论知识的掌握。 2、掌握薛定谔方程及能量本征方程的求解。 3、理解微分方程数值解求解的思路。 理论基础 薛定谔方程 薛定谔方程描述了实物粒子在的波函数随着时间、空间变化的规律，是量子力学的基本假设之一。其形式如下： 其中，m是粒子质量， 是拉普拉斯算符，描述空间中的势场分布。显然，薛定谔方程为偏微分方程。 量子力学中，粒子的波函数没有任何物理意义，但是表示在空间点找到该粒子的概率。 能量本征方程 不含时薛定谔方程与时间无关，它预言波函数可以形成驻波，称为定态（在原子物理学里，又称为轨道，例如，原子轨道或分子轨道），假若能够计算出这些定态，分析出其量子行为，则解析含时薛定谔方程会变得更为简易。不含时薛定谔方程为描述定态的方程。 3、一维方势阱 一维方势阱是量子力学中最为简单的势场分布，研究粒子在其中的分布规律，有助于对基本原理的理解，同时又是其它复杂问题的基础。一般地，一维方势阱V(x)可以表示为： 其中，为常数。 一维方势阱中，能量本征方程的形式为： 显然，能量本征方程变成了二阶线性常微分方程。 对于阱外，以区域为例，能量本征方程为： 令，则，考虑到时，，那么，可以写成：。 同理，对于区域，令，则。 对于阱内，能量本征方程为： 同时由于有限深势阱中，波函数是连续的，可以得到势阱内能量本征方程的边界条件为： 这样一来，阱内的能量本征方程转化为了带有边界条件的二阶线性常微分方程的求解问题。这样的方程不一定有解析解，但是，可以利用一些数学工具进行变形，在MATLAB中求出其数值解。 4、微分方程数值解——有限差分 二阶常微分方程的及其边界条件，如下式： （1） 下面将应用差分方法来解决这个问题。差分方法的关键，在于恰当的选取差商逼近微分方程中的导数，我们知道，逼近一阶导数可用向前差商，也可用向后差商或者中心差商.中心差商是向前差商和向后差商的算术平均.为逼近二阶导数，一般用二阶差商——向前差商的向后差商（即向后差商的向前差商） (2) 设将积分区间划分为等分，步长，节点. 差商替代相应的导数，可将边值问题(1)离散化得到下面的公式： （3） 如果函数是非线性的，那么所归结出的差分方程也是非线性的，这时实际求解困难。 如果所给方程(1)是如下形式的线性方程： (4) 则差分方程(2)相应的形式为： (5) 其中的下标表示在节点的取值。 利用边界条件(3)消除式(5)中的和，整理得到关于的下列方程组： （6） 这样归结出的方程组是所谓的三对角形的，即： (7) 将带有边界条件的二阶线性常微分方程转化为线性方程组，线性方程组的解就是微分方程的数值解。 四、实验仿真过程 1、首先我们经过小组内的讨论提出了三种形式的非对称势阱，分别是不等高方势阱、斜方势阱、正弦势阱。 ①、不等高方势阱 ②、斜方势阱 ③、正弦势阱 2、根据我们提出的三种模式结合我们已有的理论知识以及查阅的资料，我们决定通过matlab软件对我们非对称势阱进行仿真 ①、不对称高方势阱 a、matlab程序 clear; V1=3.2e-19; V2=4.8e-19; E=1.6e-19; m=9.1e-31; hb=1.05e-34; a=2e-10; beta1=sqrt(2*m*(V1-E))/hb; beta2=sqrt(2*m*(V2-E))/hb; s=exp((-1)*beta1*a/2); t=exp((-1)*beta2*a/2); lb=(-1)*a/2;ub=a/2; N=101; h=(ub-lb)/(N-1); x=linspace((-1)*a/2,a/2,N); x(1)=[]; N=N-1; for i=1:N-1 q(i)=2*m*E/hb.^2; p(i)=0; r(i)=0; end for i=1:N-1 b1(i)=-2+h.^2*q(i); end for i=1:N-2 c1(i)=1+h*p(i)/2; end fo