氢原子低能级斯塔克效应的可视化.docVIP

目 录 1 技术要求 2 2 基本原理 2 2．1 斯塔克效应 2 2．2 微扰矩阵的简化处理 3 2．2．1 氢原子状态的一级修正的能级分裂 4 2．2．2 氢原子状态的一级修正的能级分裂 4 2．2．3 氢原子状态的一级修正的能级分裂 5 3 建立模型描述 7 4 源程序代码 8 4．1 主量子数n的回调函数。 8 4．2 外电场E的回调函数 8 4．3 能级表达式的计算 8 4．4 矩阵元的计算 9 4．5 矩阵的行列式及其显示 14 4．6 重置程序 18 4．7 绘图程序 18 4．8 能级分裂值ΔE的计算 19 5 调试过程及结论 19 5．1 调试过程 19 5．2 结论 20 6 心得体会 21 7 参考文献 21 氢原子低能级Stark效应的可视化 1 技术要求 设计一个用户界面，只需输入主量子数，即可给出分裂能级的解析同时画出能级分裂图。要求： （1）有用户任意输入量子数的界面； （2）根据量子力学中对氢原子的求解，设计出各个模块的参数（例如氢原子径向矩阵元的计算）； （3）用Matlab来进行模拟。 2 基本原理 2.1 斯塔克效应 在外电场中，不计电子自旋和相对论效应，氢原子的哈密顿算符为 （1） 其中是无氢原子的外电场时氢原子的哈密顿量，且 （2） 是电子在外电场中的势能，且 （3） 式中表示均匀外电场强度的大小，按惯例将的方向取为Z轴正方向。 的本证函数和本证值分别为 （4） （表示玻尔半径） （5） 根据简并微扰论确定氢原子激发态一级斯塔克分裂能级的久期方程为 （6） 式中与久期方程（5）对应的解，即微扰矩阵的本征值，就是氢原子激发态斯塔克分裂能级的一级修正值。 2.2微扰矩阵的简化处理 为了便于程序化计算，我们对微扰矩阵作如下简化处理。 首先，由球谐函数的性质和球谐函数的归一化性质，有 （7） 可求得 （8） 所以只有满足，时，才不为零。又因为为厄米算符，其积分结果为实数。又因为 （9） 所以只须计算或即可得到不为零的微扰矩阵元。在这里取即，即不为零的微扰矩阵元为 （10） 再者，由式（8）可知，在非零微扰矩阵元中，与对应的微扰矩阵元和间有如下关系 （11） 综合以上性质，我们可以将阶的微扰矩阵按照值分解为个低阶（阶到）矩阵 ，从而使得 计算过程得到简化。例如，当时，我们将原来的阶矩阵分解为一个阶（）、两个阶（）和两个阶（）矩阵。而与所对应的矩阵是相同的，因此我们总共只需要求解个低阶矩阵的本征值。 2.2.1氢原子状态的一级修正的能级分裂 当时，未受微扰的状态是度简并的，对应的简并波函数为 （12） 将式（11）中的各函数代入式（8）得 此时微扰矩阵简化为一个阶（）、两个阶（）的低阶矩阵： 和。对应的本征值为和，能级分裂图如图。 图 能级分裂图 2.2.2氢原子状态的一级修正的能级分裂 当时，未受微扰的状态是度简并的，对应的简并波函数为 （13） 将式（12）的各函数代入式（8）得： 其他矩阵元为零，此时微扰矩阵简化为一个阶（）、两个阶（）、两个阶 的低阶矩阵： 、和 对应的本征值为、和，能级分裂图如图。 图 能级分裂图 2.2.3氢原子状态的一级修正的能级分裂 当时，未受微扰的状态是度简并的，对应的简并波函数为 （14）