一维势垒贯穿时透射系数的计算与MATLAB分析
- 格式:doc
- 大小:15.01 KB
- 文档页数:4
势垒贯穿与应用 势垒贯穿设一个质量为m 的粒子,沿x 轴正方向运动,其势能为: U(x)=0 x<0 和x>a U(x)=U 0 0≤x ≤a这种势能分布称为一维势垒。
粒子在 x < 0 区域里,若其能量小于势垒高度,经典物理来看是不能越过势垒达到 x > a 的区域。
在量子力学中,情况又如果呢?为讨论方便,我们把整个空间分成三个区域: 在各个区域的波函数分别表示为ψ1 ψ2 ψ3三个区间的薛定谔方程简化为:求出解的形式是)(),0(),0(a x a x x ≥I ∏≤≤∏≤I ),()(212122x E dx x d m ϕϕ=- 0≤x ),()()(22202222x E x U dxx d m ϕϕϕ=+- ax ≤≤0),()(232322x E dxx d m ϕϕ=- a x ≥222 mEk =2021)(2 E U m k -=,0)()(12212≤=+x x k dxx d ϕϕa x x k dxx d ≤≤=-0,0)()(221222ϕϕa x x k dxx d ≥=+,0)()(32232ϕϕikxikx e A Ae -'+=ψ1x ik Be 12+=ψikx Ce =ψ3O(1)E>U 0按照经典力学观点,在E>U 0情况下,粒子应畅通无阻地全部通过势垒,而不会在势垒壁上发生反射而在微观粒子的情形,却会发生反射。
(2)E<U 0从解薛定谔方程的结果来看,在势垒内部存在波函数ψ。
即在势垒内部找出粒子的概率不为零,同时,在x>a 区域也存在波函数,所以粒子还可能穿过势垒进入x>a 区域粒子在总能量E 小于势垒高度时仍能贯穿势垒的现象称为隧道效应定义粒子穿过势垒的贯穿系数是:透射波的概率密度与入射波概率密度的比值。
势垒高度U 0越低、势垒宽a 度越小,则粒子穿过势垒的概率就越大。
隧道效应是经典力学所无法解释的由于电子的隧道效应,金属中的电子并不完全局限于表面边界之内,电子密度并不在表面边界处突变为零,而是在表面以外呈指数形式衰减,衰减长度约为1nm只要将原子线度的极细探针以及被研究物质的表面作为两个电极,当样品与针尖的距离非常接近时,它们的表面电子云就可能重叠若在样品与针尖之间加一微小电压U b 电子就会穿过电极间的势垒形成隧道电流。
一维势垒中的透射系数利用传递矩阵方法研究了粒子在一维势垒中运动时的粒子的透射系数,主要研究的是在一个方势垒两个方势垒中透射系数,对以上的透射系数的总结,推出了对于任意势垒中透射系数, 并讨论了透射系数、反射系数与势垒宽度的关系.一维方势垒势垒模型在方势垒中,遇到的问题和 值得注意的地方。
在求方势垒波 函数中,首先要知道这是一个什 么样问题,满足什么样的方程, 方程可以写成什么样的形式,在 求解方程中,波函数的形式应该怎样需要怎样的分段,分段的过程中,特别要强调的边界条件问题。
并且验证了概率流密度。
在量子力学中,粒子在势垒附近发生的现象是不一样的,能量E 大于势垒高度0u 的粒子在势垒中有一部分发生反射,而能量小于0u 的粒子也会有部分穿过势垒,这在经典力学中是不会发生的。
下面讨论的是一维散射(即在非束缚态下问题,在无穷远处波函数不趋于零)。
重点讨论的是粒子通过势垒的透射和反射,重点在于求出波函数,这就必须求解薛定谔方程,由于)(x U 是与时间无关的,此处是定态薛定谔方程。
定态薛定谔方程通式:ψψψE U m=+∇-222h 在量子力学里, 必须知道波函数ψ, 因此必须要解薛定谔方程t i U x m ∂∂=+∂∂-ψψψh h 2222一维散射问题是一个非束缚态问题(()U x 与时间无关, 而E 是正的).因此令t Ei ex t x h-=)(),(ψψ由此得到ψψψE U dx d m =+-2222h按照势能()U x 的形式, 方程(2)一般需要分成几个部分求解.将上式改写成如下形式0222=+ψψk dxd⎩⎨⎧><<<=.,0,0;0,)(0a x x a x u x U 先讨论0u E >的情形粒子满足薛定谔方程分解为三个区域:⎪⎪⎪⎩⎪⎪⎪⎨⎧>=-<<=+-<=-a x x E x dx d m a x x E x u x dx d m x x E x dx d m ),()(20),()()(20),()(233222220222211222ψψψψψψψh h h (1) ⎪⎪⎪⎩⎪⎪⎪⎨⎧>=+<<=-+<=+a x x mEx dx d a x x u E x dx d x x mEx dxd ,0)(2)(0,0)()()(0,0)(2)(323222022212122ψψψψψψh h特征方程02=++q pr r 的两个根21,r r方程 0=+'+''qy y p y 的通解两个不相等的实根21r r ≠ x r x r e C e C y 2121+= 两个相等的实根21r r = x r e x C C y 1)(21+= 一对共轭复根βαi r ±=2,1)sin cos (21x C x C e y x ββα+=注: 0=+''qy y 的通解:特征方程02=+q r ,当0<q 时,通解xq xq eC e C y ---+=21,当0>q 时,通解xq ixq ie C e C y -+=21方程(1)的解可以表示为:⎪⎪⎪⎩⎪⎪⎪⎨⎧>+=<<+=<+=-----a x de te x a x ce be x x re ae x x mEi x mE i x u E m i x u E m i x mE i x mE i ,)(0,)(0,)(223)(2)(2222100h h hh hh ψψψ (2)定态波函数321,,ψψψ再分别乘上一个含时间的因子Et i eh-,可以看到式子(2)的三式,第一项是左向右传播的平面波,第二项是由右向左传播的平面波,即入射波和反射波。
一维高斯积分点坐标和权重系数计算在数值分析和数值计算中,一维高斯积分是一种常用的数值积分方法。
它通过一组特定的积分点坐标和权重系数来近似计算给定函数的定积分,具有高精度和快速收敛的特点。
在实际工程和科学计算中,一维高斯积分常常被广泛应用于有限元分析、求解常微分方程和偏微分方程等领域。
掌握一维高斯积分点坐标和权重系数的计算方法及其在matlab中的实现是非常重要的。
一维高斯积分的基本思想是将定积分转化为对一组特定积分点的加权求和,从而实现对定积分的近似计算。
一维高斯积分点坐标和权重系数的选择对积分精度和收敛速度有着重要影响。
常用的一维高斯积分点坐标和权重系数的计算方法有多种,其中比较常用的是Legendre多项式求解法、Jacobi矩阵方法和递推关系法等。
这些方法在matlab 中都能够很方便地实现,下面将逐一介绍它们的原理和具体实现。
1. Legendre多项式求解法Legendre多项式是一种重要的正交多项式,在计算高斯积分点和权重系数时可以利用Legendre多项式的求解方法。
通过对Legendre多项式的根和权重的计算,可以得到一维高斯积分点坐标和权重系数的近似值。
在matlab中,可以使用`legendre`函数来计算Legendre多项式的根和权重,然后进一步得到一维高斯积分点和权重系数。
2. Jacobi矩阵方法Jacobi矩阵方法是另一种计算一维高斯积分点和权重系数的常用方法。
它通过Jacobi矩阵的特征值和特征向量来计算高斯积分点和权重系数,具有较高的计算精度和稳定性。
在matlab中,可以使用`jacobi`函数来实现Jacobi矩阵方法的计算,并进一步得到一维高斯积分点坐标和权重系数。
3. 递推关系法递推关系法是一种简单而有效的计算一维高斯积分点和权重系数的方法。
它通过递推关系来逐步计算高斯积分点和权重系数,具有较快的计算速度和较高的计算精度。
在matlab中,可以编写简单的递推算法来实现一维高斯积分点坐标和权重系数的计算。
一维势垒贯穿时透射系数的计算与MATLAB分析
Abstract:The paper firstly derives transmission coefficient for one-dimensional square barrier and one-dimensional general potential, and then uses Matlab to accurate transmission coefficients. This paper also discusses the influence of the energy of incident particles, the width and height of potential to the transmission coefficient. Then this paper discusses the influence of the width of potential to calculational error in a kind of approximate algorithm for calculating transmission coefficient in many textbooks.
Key Words:Potential barrier penetration;Transmission coefficient; Matlab
在量子力学中,对于纳米量级的薄势垒层[1],如果向它射入能量比势垒高度低的粒子,在势垒后侧能够发现虽然很小但并不为零的粒子概率密度,这种穿透势垒的现象称为量子隧道效应[2-3]。
在理论和实践中计算出粒子贯穿势垒时的透射系数都具有重要意义。
量子力学是描述微观粒子运动的理论,其计算常复杂繁琐,但随着计算机技术的发展,使得许多繁琐计算都迎t刃而解。
以往很多作者对单方势垒
和多方势垒条件下的透射系数进行了研究,该文对一维一般势垒的透射系数进行理论推导,然后用Matlab程序对透射系数进行定量计算。
1 一维势垒贯穿概率模型的构建
1.1 一维方势垒
3 粒子贯穿一维势垒程序运行示例及讨论
3.1 粒子贯穿方势垒程序运行示例及讨论
利用2.2节中程序源代码,输入相关参数和相关势垒函数便可快速得到电子贯穿一维势垒时透射系数和反射系数。
表1为相关运行示例。
从表1中可以看到,透射系数D和反射系数R之和S为1;从A.1组中知E为0时透射系数为0,反射系数为1,即没有粒子去贯穿势垒,认为全部粒子都被反射;从B组1号和C组1~6号看到电子能量从小于势垒能量到等于势垒能量,再到大于势垒能量对一定高度和宽度的势垒进行贯穿时透射系数逐渐增大,反射系数逐渐减小;从D组1~3号到B 组1号可知用一定能量的电子去贯穿一定宽度的势垒时,随着势垒高度越大,透射系数逐渐减小,反射系数逐渐增大;从B组1~4号可知,电子入射能量一定,势垒高度一定,势垒宽度渐增,透射系数递减,反射系数递增,并且相当敏感。
由于透射系数与反射系数之和为1,当知道其一时便知
其二,接下来将用B组数据与大多数教材中的计算隧道效应透射系数的一种近似算法进行比较。
得d1=54.13%,d2=47.14%,d3=28.27%,d4=15.75%。
从d的变化可看到随着a递增,D1与B组中越接近。
B组中数据是没有经过近似处理的,属于准确值,上述的近似处理算法在势垒宽度a越小时误差越大。
3.2 粒子贯穿一般势垒程序运行示例及讨论
上面讨论了电子贯穿方势垒透射系数,接下来讨论电子贯穿一般势垒程序运行示例。
从表3中E组1号和D组1-3号可得随着一维线性势垒高度递增,透射系数递减,反射系数递增,从E.1和F组可看出,相同能量的电子贯穿相同高度势垒,势垒越厚实,透射系数越小,反射系数越大。
至此定量计算了电子贯穿一维方势垒和一维一般势垒的一些透射系数和反射系数,并定性讨论了电子入射能量与势垒宽度及高度对透射系数与反射系数的影响。
文中程序只是对电子贯穿一维势垒透射系数反射系数的计算。
如果要定量的计算别的粒子贯穿一维势垒的透射系数与反射系数,只须将程序中u换为相应粒子质量,便可快速得到结果。
4 结语
该文对粒子贯穿一维方势垒及一维一般势垒的透射系数进行理论推导,再用Matlab程序对电子贯穿一般势垒时
的透射系数进行数值计算。
从计算程序来看只需在MATLAB 软件环境下输入相关参数和势垒函数,便可快速得到粒子贯穿一维势垒时的透射系数。
该文利用所编程序,定量计算了电子贯穿一维势垒时透射系数,并定性的讨论了电子入射能量与势垒宽度及高度对透射系数的影响,而且把所得数据与周世勋所著《量子力学教程》中透射系数数据进行比较,讨论了众多教材中一种计算隧道效应透射系数的近似算法中势垒宽度对透射系数误差的影响。