江俊勤

(廣東第二師范學(xué)院 物理與信息工程系， 廣東 廣州 510303)

?

有限高量子圍欄中電子的能級和波函數(shù)

江俊勤

(廣東第二師范學(xué)院 物理與信息工程系， 廣東 廣州 510303)

摘要:用數(shù)值分析的方法對有限高量子圍欄內(nèi)電子的能級和波函數(shù)進(jìn)行全面的研究，糾正了現(xiàn)有文獻(xiàn)中的一些紕漏. 借助Mathematica強(qiáng)大的計算能力和杰出的數(shù)字繪圖功能，給出了能級圖和三維概率分布圖.

關(guān)鍵詞:有限高量子圍欄; 修正的貝塞爾函數(shù)；能級; 概率分布

1993年，克羅米(M. F. Crommie)等人使用掃描隧道顯微鏡(STM)將48個鐵原子吸附在銅表面上并排列成一個半徑為71.3 ?的圓環(huán)，鐵原子的平均間距為9.5 ?，這個原子圓圈能將銅表面上的電子圍住，所以這樣的原子圈被稱為“量子圍欄”或“原子圍欄”，該實(shí)驗(yàn)通過STM觀測到了圍欄內(nèi)電子概率的駐波圖.

這是人類第一次直接觀測到的電子駐波圖，對于物理學(xué)的理論和應(yīng)用都具有極其重要的意義，因此已有多位作者通過求解定態(tài)Schr?dinger方程來解釋圍欄內(nèi)電子的駐波現(xiàn)象，文獻(xiàn)[1]把問題簡化為無限高圍欄并計算了電子的能級和波函數(shù), 但是無限高圍欄顯然與實(shí)際量子圍欄存在差距；文獻(xiàn)[2]和[3]考慮了有限高圍欄，但由于圍欄的外面是無限的自由平面，電子不被圍欄束縛、其能量不能具體確定；文獻(xiàn)[4]討論了高度有限而厚度無限的圍欄模型，這樣的量子圍欄能把能量低于圍欄高度的電子束縛在其內(nèi)部，電子具有確定的能級和本征波函數(shù)，可惜的是文獻(xiàn)[4]沒有對電子的能級和狀態(tài)進(jìn)行全面的研究，而且計算過程存在多處紕漏.

本文借助Mathematica強(qiáng)大的計算能力和杰出的數(shù)字繪圖功能，對高度有限而厚度無限的圍欄內(nèi)的電子態(tài)進(jìn)行全面深入的研究：糾正了文獻(xiàn)[4]的一些紕漏、 討論能級數(shù)目與圍欄尺寸的關(guān)系、繪制出能級圖和三維概率分布圖.

1 定態(tài)Schr?dinger及其貝塞爾函數(shù)解

(1)

用分離變量法求解，令Ψ=ψ(r)Φ(φ)，代入(1)式，得

(2)

和

(3)

(2)式的解為

(4)

式中,m為角量子數(shù)(決定軌道角動量的大小)，由波函數(shù)單值性的要求Φm(φ+2π)=Φm(φ), 可得m=0,±1,±2,…； |m|的最大值由R和U0的大小決定，具體取多少，將在后面進(jìn)一步討論.

圓內(nèi)外的勢能不同，(3)式必須分兩個區(qū)域求解.

當(dāng)0≤r

(5)

式(5)是標(biāo)準(zhǔn)的貝塞爾方程，有兩個獨(dú)立的解Jm(ρ)和Ym(ρ)，分別為第一類貝塞爾函數(shù)和第二類貝塞爾函數(shù).

但是Ym(0)→-，不符合波函數(shù)有限性要求(ψ(0)=有限值)，所以符合波函數(shù)有限性要求的解為

ψ(ρ)=c1Jm(ρ)(0≤r

當(dāng)r≥R時, 若令β2=k2(U0-E)，ρ=βr，則

(6)

式(6)也有兩個獨(dú)立的解Im(ρ)和Km(ρ)，分別稱為第一類修正的貝塞爾函數(shù)和第二類修正的貝塞爾函數(shù).Im(+)→, 不符合波函數(shù)要求，所以符合波函數(shù)要求的解為ψ(ρ)=c2Km(ρ)(r≥R).

這是關(guān)于c1和c2的齊次線性方程組，有非零解的條件是系數(shù)行列式值等于0，

這是決定能量本征值的條件，求其數(shù)值解就可求得體系的能級. 第一個方程左邊接近于0，會影響計算精確度；第二個方程分母多處為零，會干擾正常計算；第三個方程最優(yōu)越，令

(7)

則體系的能級為f(E)=0的根. 體系的徑向波函數(shù)可寫成

(8)

幾點(diǎn)說明：

(iii) 由于J-m(ρ)=(-1)mJm(ρ),K-m(ρ)=Km(ρ)，從式(7)和式(8)可知-m和m所對應(yīng)的能級和徑向波函數(shù)是相同的，但是每一個能級有兩個角向波函數(shù)e±i|m|φ，所以能級是二度簡并的(m=0除外)；由于m的正負(fù)不影響能級和概率分布，所以在具體計算時可以只考慮m=0,1,2,….

2能級的數(shù)值結(jié)果

下面以Mathematica為工具，取圍欄的半徑為R=71.3 ?、圍欄的高度為U0=0.5 eV，對圍欄中電子的能級和波函數(shù)進(jìn)行具體計算，并繪制概率分布圖.

圖1　當(dāng)R=71.3 ?，U0=0.5 eV,m=0時f(E)的曲線圖

對應(yīng)于m=0，f(E)的曲線如圖1所示. 在[0, 0.05]內(nèi)E0初值太接近，可放大該區(qū)間再畫另一幅f(E)曲線圖(略)，從而確定E0的8個初值分別為0.004，0.022，0.052，0.09，0.16, 0.23, 0.32和0. 41(曲線在0.5處并未與橫軸相交).

圖2　當(dāng)R=71.3 ?，U0=0.5 eV時f(E)的曲線實(shí)線:m=10；長劃線:m=21；短劃線:m=22

對于m=1,2,…的情況，用同樣的步驟逐一繪制f(E)的曲線圖，結(jié)果發(fā)現(xiàn): 對于固定圍欄半徑R和高度U0，隨著m的增大f(E)=0的實(shí)根個數(shù)逐漸減少；對于R=71.3 ?和U0=0.5 eV，當(dāng)m=10時只有四個實(shí)根(E=0不是本問題的解，不計算在內(nèi))，當(dāng)m=18～21時都只有一個正實(shí)根(與橫軸正向只有一個交點(diǎn))，從m=22開始沒有實(shí)根，圖2給出了m=10，21和22時f(E)的曲線. 所以對應(yīng)于這樣的圍欄參數(shù)，角量子數(shù)的取值范圍為m=0,±1,±2,…,±21(對應(yīng)的初值不在這里一一列出).

表1 給出了m=0,1,2,…,21所對應(yīng)的能級個數(shù)(-m與m的能級數(shù)相同).

表1　有限高量子圍欄(R=71.3 ?和U0=0.5 eV)的角量子數(shù)及其能級個數(shù)N

所以當(dāng)R=71.3 ?和U0=0.5eV時，能級個數(shù)為90；圍欄半徑R或者高度U0越大，m的取值范圍就越大. 如果把圍欄高度近似簡化為無限高，則m的取值范圍就是無窮(一切整數(shù))，這就是文[1]的情形.

有了能級的大概值(找根的初值)，就可以尋找當(dāng)R=71.3 ?和U0=0.5eV時m=0,1,2的能級了.對應(yīng)于m=0的能級(以eV為單位，下同)分別為: 0.004 016 71, 0.021 154 6, 0.051 947 5, 0.096 326 5, 0.154 16, 0.225 217, 0.309 051和0.404 554；

圖3　前三個角量子數(shù)的能級總圖實(shí)線,長劃線，短劃線分別對應(yīng)于m=0,1,2

對應(yīng)于m=1的能級分別為: 0.010 196, 0.034 158 8, 0.071 757, 0.122 887, 0.187 376, 0.264 909, 0.354 822和0.454 964；

對應(yīng)于m=2的能級分別為:0.018 3131, 0.049 154 6, 0.093 557 3, 0.151 414, 0.222 503, 0.306 39和0.402 018.把m=0,1,2三種不同角動量的能級畫在一起，得如圖3所示的能級總圖.

對于m=0,1,2三種不同角量子數(shù)，能級是很有規(guī)律的：從低到高分別是m=0,1,2的三個最低能級(m=0的第一個能級為基態(tài)能)，接著是m=0,1,2的三個次低能級(m=0的次低能級略高于m=2的最低能級)，…，最后是m=0的第八個能級及m=1的第八個能級．

如果把m=3,4,…,21的能級也一起排列，情況就復(fù)雜了，具體情況有待進(jìn)一步計算.

3波函數(shù)與概率分布的數(shù)值結(jié)果

把已經(jīng)得到的能級逐個代入式(8)即可確定每個能級所對應(yīng)的波函數(shù)，借助Mathematica可以自動快速得到波函數(shù)，并繪制電子的概率分布圖(立體圖)，對于m=0，R=71.3 ?和U0=0.5eV，如圖4～圖11所示. 對于m=1，只給出第四個能級所對應(yīng)的電子概率分布圖，如圖12所示；對于m=2，只給出第三個能級所對應(yīng)的電子概率分布圖，如圖13所示.

以上結(jié)果表明：角動量為零 (m=0) 的量子態(tài)在圍欄中心(r=0)是概率主峰，角動量不為零 (m≠0) 的量子態(tài)在圍欄中心(r=0)是概率低谷(波函數(shù)的節(jié)點(diǎn)). 此外，不論角量子數(shù)m是多少，概率分布圖的波峰個數(shù)都與能級E(m,n)的序號n相一致.

最后，研究束縛態(tài)的最高能級(m=21, 即能量最接近0.5eV)的量子狀態(tài)，數(shù)值結(jié)果如圖14所示. 電子概率分布主要集中在圍欄邊緣.

圖14m=21的能級(唯一且最高的能級)所對應(yīng)的概率密度

參考文獻(xiàn)：

[1] 周佩瑤，王萊，陳宜生. 原子圍欄中金屬表面電子的態(tài)函數(shù)及能級[J]．工科物理，2000，10(4)：9-11.

[2] 戴淇, 蔣平. 量子圍欄的電子結(jié)構(gòu)(英文)[J].復(fù)旦學(xué)報(自然科學(xué)版), 2000, 39(2):132-134.

[3] 塔立, 單崳瓊, 周向仁, 等. 量子圍欄中的表面電子局域態(tài)密度及其模擬[J]．云南大學(xué)學(xué)報(自然科學(xué)版)，2009，31(S2)：411-416.

[4] 董建．Mathematica與大學(xué)物理計算 [M]．2版. 北京：清華大學(xué)出版社，2013：289-296.

收稿日期：2016-01-21

基金項(xiàng)目：廣東省高等學(xué)校專業(yè)綜合改革試點(diǎn)項(xiàng)目(2014)

作者簡介：江俊勤，男，廣東揭陽人，廣東第二師范學(xué)院物理與信息工程系教授.

中圖分類號：O562.1

文獻(xiàn)標(biāo)識碼：A

文章編號：2095-3798(2016)03-0057-06

Energy Level and Wave Function of Electron inside Finite High Quantum Corral

JIANG Jun-qin

(Department of Physics and Information Engineering, Guangdong University of Education,Guangzhou，Guangdong, 510303, P.R.China)

Abstract:The energy levels and wave functions of the electron in the finite high quantum corral are numerically studied. Some errors in the existing literature are corrected. The energy level diagrams and the three dimensional probability distribution diagrams are given by using MATHEMATICA.

Key words:finite high quantum corral; modified Bessel function; energy level; probability distribution