基于GUI类氢离子中电子概率密度的可视化设计

2019-03-22高志华

物理与工程 2019年1期

李轩袁琦高志华，2

(东北师范大学 1物理学院, 2物理学国家级实验教学示范中心，吉林长春 130024)

量子力学作为四大力学之一，在物理类本科学习阶段非常重要。但因其涉及的诸多物理问题与经典理论相差较大，相关理论的推导过程亦较为繁杂，故量子力学对于初学者来说接受起来往往较为困难。随着现代化教学手段的发展，基于Matlab等一些软件平台对一些物理问题进行可视化模拟被得到广泛研究[1-5]。本文利用Matlab软件中的GUI功能对氢原子或类氢离子中电子的概率密度分布问题进行可视化设计，利用该程序输入核电荷数、主量子数、角量子数和磁量子数等就可以直接得到相应的概率密度分布图像，方便快捷。

1 类氢离子中电子的概率密度

类氢离子中电子的定态波函数为

ψn lm(r，θ，φ)=Rn l(r)Yl m(θ，φ)

(1)

其中，Rn l(r)为径向分布函数；Yl m(θ，φ)为角向分布函数；n、l、m分别是主量子数、角量子数和磁量子数，每一组量子数(n，l，m)对应一组波函数。

类氢离子的径向分布函数和角向分布函数分别为

利用径向波函数，可以求得电子在(r，r+dr)球壳中的概率为

(4)

同理，可以求得电子在(θ，φ)方向的立体角dΩ中的概率为

wl m(θ，φ)=|Yl m(θ，φ)|2dΩ

(5)

2 基于GUI类氢离子中电子概率密度的可视化设计

2.1 类氢离子中电子概率密度的径向角向图像可视化设计

图1 类氢离子中电子概率密度分布可视化流程图

图1给出基于GUI得到类氢离子中电子概率密度分布的流程图。由图可知，概率密度分布界面包括了设置可编辑文本框、单功能按钮、绘图等控件属性及编写相应的回调函数(call-back)。考虑可编辑文本框输入的量子数之间必须满足l=0，1，…，n-1；m=-l，-l+1，…，l；在编写回调函数时，利用如图2所示的错误对话框以确保所输入的量子数准确有效。

图2 错误对话框提示

在量子力学教材和一些文献中，容易找到低阶量子态氢原子或类氢离子中电子的概率密度分布图像。为验证程序的正确性，下面以氦离子为例，取Z=2，给出部分量子数下的电子概率密度分布图像。同时考虑角向概率密度分布为三维图像，利用max函数和text语句得到角向分布的最大概率密度值，如图3所示。

图3 氦离子(Z=2)中电子概率密度分布运行图示(a) n=11，l=10，m=1时,He+中电子概率密度分布图像；(b) n=11，l=10，m=2时,He+中电子概率密度分布图像；(c) n=11，l=10，m=3时,He+中电子概率密度分布图像；(d) n=11，l=10，m=4时,He+中电子概率密度分布图像；(e) n=11，l=10，m=5时,He+中电子概率密度分布图像；(f) n=11，l=10，m=6时,He+中电子概率密度分布图像；(g) n=11，l=10，m=7时,He+中电子概率密度分布图像；(h) n=11，l=10，m=8时,He+中电子概率密度分布图像；(i) n=11，l=10，m=9时,He+中电子概率密度分布图像；(j) n=11，l=10，m=10时,He+中电子概率密度分布图像