matlab画本征函数系,双势垒结构能量本征值和本征波函数系数的MATLAB计算

云光技术 第 45卷 第 2期 双势垒结构能量本征值和本征波函数系数的 MATLAB计算 王忆锋 (昆明物 理研 究所，云南 昆明 650223) 摘要：晶体中的电子状态可以用波函数来描述。从计算角度来看，必须知道能量 本征值以及波函数的系数才能确定波函数的具体形态。在计算一维方势阱束缚态 能量本征值时，基于边界连续性条件，可将相关波函数展开为一个以矩阵形式描 述的线性方程组。根据能量本征值必须满足该方程组的系数行列式等于零的要 求，在能量区间内逐点扫描，即可确定相应的能量本征值，由此进一步确定波函 数的系数。该方法不需要花费较多精力编程，具有概念简单、使用方便、实用性 强等特点。 关键词：半导体；能量本征值；薛定谔方程；波函数系数 中图分类号 ：0471．1 文献标识码 ：A 0 引言 晶体 中的 电子状态可 以用波 函数来描 述 。从计算角度来看 ，必须知道 能量本征 值 以及波函数 的系数才能确定波 函数 的具 体形态 。本文 以双势垒结构为例，介绍 了 利用 MATLAB计算能量本征值及本征波函 数系数的方法 。 1 势能函数坐标系的选取与波函数 界面连续性条件 一 维 方 势 阱可 以视 为一 个双 势 垒结 构。有关计算途径或计算结果 的表达 形式 会因势能函数坐标系 的不同取法而异。被 限制在势阱 中运动 的粒子所处的状态称 为 束缚态。由于一维束缚态是非简并的，故 在一个对称势中，每一个束缚态必须是偶 宇称 的，或者是奇宇称的L1 ；两者对应的是 数学中偶 函数触)= ～ )和奇函数 )=一 触 )的概念 。换言之 ，波函数可 以是偶函数， 也可 以是奇 函数 ，前提是首先构成一个满 足 )一 -x)关系的对称势场 。例如在 图 1(a)中，势阱的势能函数 )是关于纵轴对 称 的，因此可 以利用宇称性质；而在 图 l(b) 选择的坐标系下 ，势 阱中的势能函数 ) 则不是关于纵轴对称 的。它们 的数学表达 式可 以写为： f <0) ( )={0 (0≤x 口) I > ) f 口／2) f11 一 般情况下，界面连续性条件为： I ( )= +l(X ) { ’ ( )一(／9tj+1( ) 【 + f21 式(2)是能量本征值及本征 函数计算 的基础 。 2013笠 云光技术 第 45卷 第 2期 ) + I I ：2 - a／2 0 a／2 (a) 一 一 一 一 ) + 2 0 a (b) 一 一 一 一 图 1 一维方形势阱坐标轴 的 2种取法 2 称性 质求解 波 函数 cI勺台皂量 求出 解M 为 ATLAB的 dso 可以 本征值 ～ 一～ 为了利用宇称性质 ，坐标系取图 1(a) ( ) sin( )+ c。 ( ) (6) 的形式 。这时两个势垒区中的薛定谔方程 考虑到波函数的有界性，各区域 的波 为： 函数有下列形式： d2( , o 一 砰 =。 ： 2=x／2mb(E-Vo) {仍(x)=Clexp(klx) 乏 荔 式中：h 为普朗克常数；mb为势垒中的粒 子有效质 量 ；E 为入射粒 子能量 。利用 MATLAB的 dsolve()命令，可以求出其通 解为： ( )=cl exp(~x)+CE exp(一k,x) (4) 式中：cl、c2为待定系数。 势阱中的波函数满足下列薛定谔方程 2 +砖 ㈣ ： — 2~x／2 — mwE h 式 中：势阱中的波函数仍(x)适于利用 宇称 性质。对于偶态波函数，有 ∽ ： (一 ， 相应地，C3：0