張子珍，楊成全，康占成

（山西大同大學(xué)物理與電子科學(xué)學(xué)院,山西大同 037009）

1905年,A.Einstein用Plank的量子假設(shè)去解決光電效應(yīng)問題時,提出了光量子的概念。De Broglie仔細分析了光的微粒說及波動說發(fā)展的歷史,提出了物質(zhì)波的假說。1926年,Born把微觀粒子的波動性與粒子性統(tǒng)一起來,提出的幾率波概念。1927年,Heisenberg提出測不準關(guān)系,將經(jīng)典粒子的概念對微觀世界的適用程度作了最集中和最形象的概括。微觀粒子的量子態(tài)用波函數(shù)來描述,波函數(shù)隨時間變化的規(guī)律就是Schrodinger方程,Schrodinger方程是量子力學(xué)中最基本的方程,其地位與Newton方程在經(jīng)典力學(xué)中的地位相當(dāng)。在不同的勢場作用下,求出Schrodinger方程的解析解比較困難,對有些勢場(如wood-Sxaon勢,一維模型勢等)要嚴格求出其解析解幾乎不可能。為此,本文借助計算機編程,以有限深方勢阱為例來求Schrodinger方程的數(shù)值解,并將結(jié)果與解析解進行比較,為求解Schrodinger方程提供一種簡捷的方法。

1 Schrodinger方程的引入

與具有一定能量E及動量p→的粒子相聯(lián)系的波是平面單色波,其表達式是

將上式對時間及坐標求導(dǎo),并利用自由粒子的能量動量關(guān)系 E＝p2／2m,可得

如果粒子在勢場中運動,經(jīng)典粒子的能量動量關(guān)系式是

作用于波函數(shù)上,即得

(5)式就是定態(tài)Schrodinger方程。

其中E是分離變量時引入的常數(shù)。方程(4)的解是

2 定態(tài)Schrodinger方程的解

2.1 有限深方勢阱中的解析解

有限深方勢阱[1]的表達式為:

在勢阱內(nèi),Schrodinger方程為

奇宇稱態(tài)

(10)和(11)是超越方程組,用圖解法近似計算。在計算中取V0=1 000 MeV,a=4.0 Fm,m取質(zhì)子的靜止質(zhì)量,mc2=938 MeV?？勺鞒雠c方程組(10)和(11)對應(yīng)的圖 1和圖2。在圖1中,實線代表虛線代表 ξtgξ＝η。在圖 2 中,實線代表,虛線代表－ξctgξ＝η。

從兩曲線的交點坐標可算出,當(dāng)V0=1 000 MeV,a=4.0 Fm時,各束縛態(tài)能級分別是E1=10.99 MeV,E2=44.51 MeV,E3=99.96 MeV,E4=176.85 MeV,E5=276.73 MeV。

圖1 偶宇稱態(tài)

圖2 奇宇稱態(tài)

2.2 有限深方勢阱中的數(shù)值解

2.2.1 坐標空間的實穩(wěn)定方法

實穩(wěn)定方法是通過基展開的方法進行求解。它是基于這樣的事實:在分立能區(qū),各本征態(tài) (束縛態(tài))的能量不會隨基空間的維數(shù)而改變;而在連續(xù)能區(qū),本征態(tài)的能量多數(shù)會隨著基空間維數(shù)的增加而逐漸降低,這些態(tài)對應(yīng)于散射態(tài);還有一些本征態(tài)能量隨基空間維數(shù)的增加會出現(xiàn)一個平臺,這些態(tài)對應(yīng)于共振態(tài)[2]。實穩(wěn)定方法也可以在坐標空間(離散的空間格點上)求解。通過改變坐標空間的大小,可以得到本征能量隨格點數(shù)的變化關(guān)系。與基展開的方法類似,能量不隨坐標空間格點數(shù)而改變的態(tài)就對應(yīng)于束縛態(tài),能量隨坐標空間格點數(shù)的增加而減小的態(tài)就對應(yīng)于散射態(tài);而隨坐標空間格點數(shù)的增加能量出現(xiàn)一個穩(wěn)定的平臺,再增加格點數(shù),能量會逐漸減小,這樣的態(tài)就對應(yīng)共振態(tài)[3]。

2.2.2 數(shù)值求解過程

1)給定初始能量、兩端初始波函數(shù)及其一階導(dǎo)數(shù)值。

2)通過Runge-Kutta方法[4]向match點推近。

3)利用match點波函數(shù)是否光滑對接及節(jié)點數(shù)判斷初始能量是否合適。

4)若節(jié)點數(shù)符合,波函數(shù)光滑對接,歸一化波函數(shù)之后輸出計算結(jié)果。

5)若節(jié)點數(shù)不符,則改變能量初值,重新進行計算,直到求得本征能量[5]。

2.2.3 數(shù)值計算結(jié)果

阱寬固定,阱深不斷加深,數(shù)值計算得出各束縛態(tài)的能級如表1。阱深為V0=1 000 MeV,阱寬為a=4.0 fm時數(shù)值計算得到的前3級波函數(shù)如圖3。

阱深固定,阱寬不斷增加,數(shù)值計算得到的各束縛態(tài)能級如表2。

從表2中可以看出,當(dāng)勢阱的深度與寬度比較接近時,幾乎找不到束縛態(tài)。而固定阱寬,當(dāng)阱深加深時能級逐漸升高。當(dāng)阱深加深到V0=1.0×106MeV時,它的各能級高度與無限深方勢阱的各能級高度就相當(dāng)接近了?？梢酝茢?阱寬不變,繼續(xù)加深阱深,各能級與無限深勢阱的極限值將會相當(dāng)吻合。

從表2中發(fā)現(xiàn),固定阱深,寬度越小,能級越高,寬度越大,能級越低。而且各能級的能量值與阱寬a2基本上符合反比規(guī)律。從表1和表2中均可發(fā)現(xiàn),當(dāng)阱寬和阱深固定時,各能級的能量仍保持與量子數(shù)n2成正比的規(guī)律。

表1 阱寬,a=4.0 Fm,阱深由10～107MeV 變化時各束縛態(tài)能級 單位:MeV

圖3 阱深為V0=1000 MeV,阱寬為a=4.0 Fm時的前3級波函數(shù)圖

表2 阱深為V0=106MeV,阱寬由1.0～100.0 Fm變化時各束縛態(tài)能級 單位:MeV

3 兩種解的比較

為說明數(shù)值計算方法的有效性,下面將V0=1 000 MeV,a=4.0 Fm,得出的解析解與數(shù)值解進行比較。表3中第2行為解析求解得到的能級,第3行為數(shù)值計算得到的能級。從表3中的數(shù)據(jù)可以看出,數(shù)值解得到的能級與解析解得出的結(jié)果基本上是一致的。數(shù)值求解不僅計算準確,更重要的是對于那些比較復(fù)雜的勢 (如wood-Sxaon勢,一維模型勢等),根本無法求出其解析解,但用數(shù)值求解時,只要在程序中把勢的表達式換一下,方程就可迎刃而解。因此,數(shù)值計算為求解Schrodinger方程提供一種簡捷而有效的方法。

表3 兩種解的比較 單位:MeV

[1]曾謹言.量子力學(xué)(卷Ⅰ)[M].北京:科學(xué)出版社,1995:101-105.

[2]Hazi A U,Taylor H S.Stabilization Method of Calculating Resonance Energies:Model Problem[J].Phys Rev A,1970,1:1109.

[3]張力,周善貴,孟杰,等.單粒子共振態(tài)的實穩(wěn)定方法研究[J].物理學(xué)報,2007,56(7):3839-3844.

[4]聶鐵軍.數(shù)值計算方法[M].北京:高等教育出版社,1990:247-256.

[5]譚浩強,田淑清.FORTRAN語言程序設(shè)計[M].北京:高等教育出版社,1986.