李 軒 袁 琦 高志華，2

(東北師范大學(xué) 1物理學(xué)院, 2物理學(xué)國(guó)家級(jí)實(shí)驗(yàn)教學(xué)示范中心，吉林 長(zhǎng)春 130024)

量子力學(xué)作為四大力學(xué)之一，在物理類本科學(xué)習(xí)階段非常重要。但因其涉及的諸多物理問(wèn)題與經(jīng)典理論相差較大，相關(guān)理論的推導(dǎo)過(guò)程亦較為繁雜，故量子力學(xué)對(duì)于初學(xué)者來(lái)說(shuō)接受起來(lái)往往較為困難。隨著現(xiàn)代化教學(xué)手段的發(fā)展，基于Matlab等一些軟件平臺(tái)對(duì)一些物理問(wèn)題進(jìn)行可視化模擬被得到廣泛研究[1-5]。本文利用Matlab軟件中的GUI功能對(duì)氫原子或類氫離子中電子的概率密度分布問(wèn)題進(jìn)行可視化設(shè)計(jì)，利用該程序輸入核電荷數(shù)、主量子數(shù)、角量子數(shù)和磁量子數(shù)等就可以直接得到相應(yīng)的概率密度分布圖像，方便快捷。

1 類氫離子中電子的概率密度

類氫離子中電子的定態(tài)波函數(shù)為

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

(1)

其中，Rn l(r)為徑向分布函數(shù)；Yl m(θ，φ)為角向分布函數(shù)；n、l、m分別是主量子數(shù)、角量子數(shù)和磁量子數(shù)，每一組量子數(shù)(n，l，m)對(duì)應(yīng)一組波函數(shù)。

類氫離子的徑向分布函數(shù)和角向分布函數(shù)分別為

利用徑向波函數(shù)，可以求得電子在(r，r+dr)球殼中的概率為

(4)

同理，可以求得電子在(θ，φ)方向的立體角dΩ中的概率為

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

(5)

2 基于GUI類氫離子中電子概率密度的可視化設(shè)計(jì)

2.1 類氫離子中電子概率密度的徑向角向圖像可視化設(shè)計(jì)

圖1 類氫離子中電子概率密度分布可視化流程圖

圖1給出基于GUI得到類氫離子中電子概率密度分布的流程圖。由圖可知，概率密度分布界面包括了設(shè)置可編輯文本框、單功能按鈕、繪圖等控件屬性及編寫(xiě)相應(yīng)的回調(diào)函數(shù)(call-back)?？紤]可編輯文本框輸入的量子數(shù)之間必須滿足l=0，1，…，n-1；m=-l，-l+1，…，l；在編寫(xiě)回調(diào)函數(shù)時(shí)，利用如圖2所示的錯(cuò)誤對(duì)話框以確保所輸入的量子數(shù)準(zhǔn)確有效。

圖2 錯(cuò)誤對(duì)話框提示

在量子力學(xué)教材和一些文獻(xiàn)中，容易找到低階量子態(tài)氫原子或類氫離子中電子的概率密度分布圖像。為驗(yàn)證程序的正確性，下面以氦離子為例，取Z=2，給出部分量子數(shù)下的電子概率密度分布圖像。同時(shí)考慮角向概率密度分布為三維圖像，利用max函數(shù)和text語(yǔ)句得到角向分布的最大概率密度值，如圖3所示。

圖3 氦離子(Z=2)中電子概率密度分布運(yùn)行圖示(a) n=11，l=10，m=1時(shí),He+中電子概率密度分布圖像；(b) n=11，l=10，m=2時(shí),He+中電子概率密度分布圖像；(c) n=11，l=10，m=3時(shí),He+中電子概率密度分布圖像；(d) n=11，l=10，m=4時(shí),He+中電子概率密度分布圖像；(e) n=11，l=10，m=5時(shí),He+中電子概率密度分布圖像；(f) n=11，l=10，m=6時(shí),He+中電子概率密度分布圖像；(g) n=11，l=10，m=7時(shí),He+中電子概率密度分布圖像；(h) n=11，l=10，m=8時(shí),He+中電子概率密度分布圖像；(i) n=11，l=10，m=9時(shí),He+中電子概率密度分布圖像；(j) n=11，l=10，m=10時(shí),He+中電子概率密度分布圖像

對(duì)于高階量子態(tài)氫原子或類氫離子中電子的概率密度分布，由于波函數(shù)形式非常復(fù)雜，所以一般在量子力學(xué)教材中很少涉及。通過(guò)該可視化軟件可以得到任意一組量子數(shù)對(duì)應(yīng)的電子概率密度分布圖像。

圖4所示的氦離子中電子概率密度分布圖像由于核電荷數(shù)和量子數(shù)大小的增加，圖像變得更為復(fù)雜。徑向概率密度分布曲線存在14個(gè)極大值和13個(gè)節(jié)點(diǎn)，主極大距原子核距離較大，同時(shí)主極大的相對(duì)概率密度較低；角向概率密度分布圖像呈帶有圓盤(pán)的雙缽狀，與相關(guān)理論相符。

圖4 n=25，l=11，m=9時(shí)，He+(Z=2)中電子概率密度分布運(yùn)行圖示

2.2 主量子數(shù)一定時(shí)，類氫離子中電子徑向概率密度分布比較的可視化設(shè)計(jì)

為更好地研究在主量子數(shù)n一定時(shí)，氫原子或類氫離子中電子徑向概率密度隨角量子數(shù)l的分布規(guī)律，設(shè)計(jì)出圖5所示的可視化程序，通過(guò)該程序輸入核電荷數(shù)和主量子數(shù)即可得到不同角量子數(shù)下的徑向概率密度分布比較圖像。

圖5 n=3時(shí)，氫原子(Z=1)中電子概率密度分布運(yùn)行圖示

圖5說(shuō)明在主量子數(shù)n已知的情況下，角量子數(shù)可取值l=0，1，…，n-1；每條曲線有n-l個(gè)峰，n-l-1個(gè)節(jié)點(diǎn)。特別地，對(duì)于l=n-1對(duì)應(yīng)的“圓軌道”，峰對(duì)應(yīng)的橫坐標(biāo)為r=n2az，此時(shí)峰值為曲線的最大值，與已知理論一致，表明該可視化程序是正確的。

3 結(jié)語(yǔ)

利用Matlab制作人機(jī)交互界面，直接輸入或改變相應(yīng)參數(shù)就可得到所需氫原子或類氫離子中電子概率密度分布圖像，有助于量子力學(xué)等課程的多維度教學(xué)，為理解高階量子態(tài)氫原子或類氫離子中電子概率密度分布提供了一種途徑，對(duì)Matlab更好地在物理教學(xué)中發(fā)揮作用做出有益的嘗試。