李明明，劉海英(濟(jì)南大學(xué)物理科學(xué)與技術(shù)學(xué)院，山東濟(jì)南　250022)

?

計(jì)算δ函數(shù)勢(shì)阱中平均動(dòng)能的方法

李明明，劉海英
(濟(jì)南大學(xué)物理科學(xué)與技術(shù)學(xué)院，山東濟(jì)南250022)

摘要:用量子力學(xué)基本公式F=∫ψ*F^ψdτ計(jì)算粒子在δ勢(shì)阱中運(yùn)動(dòng)的平均動(dòng)能時(shí)，結(jié)果為負(fù)，這是由求解δ勢(shì)阱中運(yùn)動(dòng)方程的特殊方法造成的。采用由解定態(tài)方程得到的定態(tài)能量減去平均勢(shì)能得到平均動(dòng)能的新方法，計(jì)算公式為T= E－∫ψ*Vψdτ。并研究了粒子在勢(shì)阱V( r) =－γδ( r－a)基態(tài)的平均動(dòng)能。

關(guān)鍵詞:δ勢(shì)阱;平均動(dòng)能;維里定理

在量子力學(xué)中，計(jì)算粒子的力學(xué)量F^平均值F—是一個(gè)核心問(wèn)題［1－7］，其計(jì)算公式［8］為:

式中:ψ是粒子的歸一化波函數(shù);ψ*是ψ的復(fù)共軛。

在微觀領(lǐng)域中，δ函數(shù)勢(shì)阱是一種反映微觀現(xiàn)象共同特征的常用的理想模型［9－12］。然而，當(dāng)粒子處于δ函數(shù)勢(shì)阱束縛定態(tài)時(shí)，用式( 1)算出的平均動(dòng)能是負(fù)的。這說(shuō)明式( 1)不能用來(lái)計(jì)算δ函數(shù)勢(shì)阱中粒子的動(dòng)能。鑒于此，本文提出了一種函數(shù)勢(shì)阱中粒子的平均動(dòng)能的計(jì)算方法，并討論了式( 1)為負(fù)的原因。

1　在勢(shì)阱V( x) =－γδ( x)中的定態(tài)能量和波函數(shù)

粒子波函數(shù)ψ( x)滿足的定態(tài)方程為:

式中:μ為粒子的質(zhì)量; E為束縛態(tài)能量;γ為代表δ勢(shì)阱強(qiáng)度的常數(shù);為普朗克常數(shù)。

式( 2)存在束縛態(tài)的條件是E＜0。令

式( 2)變?yōu)?/p>

求解式( 3)的方法是去掉勢(shì)能發(fā)散點(diǎn)x=0，再解方程。在不包含x=0的全部x空間，式( 3)變?yōu)?/p>

式( 4)的解為ψ=e±ax。式( 4)滿足束縛態(tài)條件ψ( +∞)及ψ在x=0連續(xù)的一般解為

由波函數(shù)ψ的歸一化條件∫+∞ψ( x)2dx = 1，得A=槡α，則式( 5)可寫為:

－∞

式( 6)中α決定粒子的能量E，由ψ'( x)在x=0處的不連續(xù)條件［13］

確定為

由此得束縛定態(tài)能量

將式( 6)代入式( 1)，可得:

平均動(dòng)能為負(fù)，這顯然不對(duì)。因此，式( 1)不能用來(lái)計(jì)算平均動(dòng)能，同樣不能用來(lái)計(jì)算平均能量，因?yàn)?/p>

然而，可以證明，式( 1)計(jì)算平均勢(shì)能是正確的:

根據(jù)維里定理［14］，在束縛定態(tài)下

如果V( r→)是r→的v次齊次函數(shù)，則r→·▽V( r→) = vV( r→)，式( 9)寫為:

對(duì)于V( x) =－γδ( x)，V(λx) =－γδ(λx) =λ－1γδ( x) =λ－1V( x)，V( x)是x的v=－1的齊次函數(shù)。于是，式( 10)改寫為:

將式( 8)代入式( 11)，得T =μγ2/( 22)。由此得到的能量為:

這個(gè)結(jié)果同式( 7)一致。這就證明了，用式( 1)來(lái)計(jì)算粒子的平均勢(shì)能是正確的。因此，正確計(jì)算平均動(dòng)能的方法是由解定態(tài)方程得到的定態(tài)能量E減去由公式( 1)的算出平均勢(shì)能V ，得:

由式( 12)算出平均動(dòng)能T=μγ2/( 22)，這是正確的結(jié)果。雖然，式( 11)也給出了正確的平均動(dòng)能，但是式( 11)不具有一般性，如果δ函數(shù)勢(shì)阱形式稍有改變，如改為三維勢(shì)V(→r) =－γδ(→r－→α)，式( 11)就不適用了。而式( 12)具有普遍性，因?yàn)樗鼇?lái)源于基本的能量公式H=T+V。

2　粒子在勢(shì)阱V( r) =－γδ( r－a)基態(tài)的平均動(dòng)能

在中心力場(chǎng)V( r) =－γδ( r－a)中，軌道角動(dòng)量是守恒量。令定態(tài)波函數(shù)

式中Ylm(θ，φ)為球諧函數(shù)。

將式( 13)代入定態(tài)方程

得徑向波函數(shù)u1( r)滿足的方程［15］

基態(tài)必定是l=0的態(tài)。令l=0的波函數(shù)u0( r)≡u(píng)( r)。在式( 14)中，令l=0得u( r)滿足的方程

在勢(shì)阱V( r) =－γδ( r－a)中，要形成束縛態(tài)，能量E＜0。令

式( 15)變?yōu)?

在r≠a的全空間，方程( 17)改寫為:

式( 17)的解為u( r) = e±kr。方程( 16)的一般解由ekr與e－kr的線性組合構(gòu)成，而滿足條件u( 0) = 0與u(∞) = 0的解為

由波函數(shù)的連續(xù)條件u1( a) = u2( a)得B=( e2ka－1) A，代入式( 18)，得

于是，l=0的基態(tài)波函數(shù)為

由波函數(shù)歸一化條件∫∞ψ( r)24πr2dr = 1，得ψ( r)中的歸一化常數(shù)

0

徑向波函數(shù)一階導(dǎo)數(shù)u'( r)在r=a處的不連續(xù)條件為:

波函數(shù)k中以及定態(tài)能量E決定。將式( 19)代入式( 21)，得

式( 22)是k或k中所含定態(tài)能量E滿足的超越方程。這個(gè)方程只能用作圖法或數(shù)值計(jì)算法求解。由作圖法得，方程( 22)有解的條件是勢(shì)阱強(qiáng)度γ＞2/( 2μa)。選擇勢(shì)阱強(qiáng)度γ=2/(μa)，用數(shù)值計(jì)算法計(jì)算值與定態(tài)能量E。令γ=2/(μa)，并令x=2ka，式( 22)變?yōu)?/p>

用數(shù)值計(jì)算法求得式( 23)的解為x=1．593 6。再由x=2ka=槡2μ E/2，算出

將式( 20)代入式( 1)得

計(jì)算出的平均動(dòng)能是負(fù)的。這又一次說(shuō)明公式( 1)是不適用的，必須用式( 12)計(jì)算平均動(dòng)能。式( 12)中的定態(tài)能量E已經(jīng)求出，E=－0．317 42/(μa2)，則平均勢(shì)能

3　討論

平均動(dòng)能取負(fù)值的原因是由δ勢(shì)場(chǎng)定態(tài)方程的求解方法決定的。以式( 13)為例。由于在此勢(shì)阱的束縛態(tài)能量E＜0，令E=－|E|。方程變?yōu)?/p>

式( 25)含有δ函數(shù)，這個(gè)方程是無(wú)法求解的。為了能求解這個(gè)方程，不得不去掉r=a點(diǎn)，在不包括r=a點(diǎn)的空間，方程變?yōu)檫@正是動(dòng)能算符T=－2▽2/( 2μ)的本征值T=－|E|的本征方程。這就是說(shuō)，要求解的波函數(shù)是動(dòng)能算符的本征值為－|E|(這也是束縛定態(tài)能量)的本征函數(shù)E。因此，由這個(gè)波函數(shù)算出的動(dòng)能平均值必然是負(fù)的。

在求解定態(tài)方程時(shí)，去掉了勢(shì)能發(fā)散點(diǎn)，因而去掉了δ函數(shù)勢(shì)，卻能求出δ函數(shù)勢(shì)中的定態(tài)能量。這是因?yàn)橛?jì)算中用到ψ'在勢(shì)能發(fā)散點(diǎn)r=a的不連續(xù)條件這個(gè)條件是由含有δ函數(shù)勢(shì)的定態(tài)方程推導(dǎo)出的。

參考文獻(xiàn):

［1］廖榮寶．波函數(shù)塌縮問(wèn)題及物理量期望值表達(dá)式的討論［J］．中國(guó)科技信息，2011( 21) : 41．

［2］馬廣清，肖鴻飛．計(jì)算力學(xué)量平均值的一種方法［J］．沈陽(yáng)大學(xué)學(xué)報(bào)，2005，17( 2) : 98－100．

［3］劉樹(shù)信．力學(xué)量的平均值［J］．煙臺(tái)師范學(xué)院學(xué)報(bào)(自然科學(xué)版)，1996，12 ( 2) : 153－155．

［4］樊西漢．量子力學(xué)平均值計(jì)算時(shí)應(yīng)注意的問(wèn)題［J］．臨沂師范學(xué)院學(xué)報(bào)，2001，23 ( 4) : 24－26．

［5］高廣君．量子力學(xué)中力學(xué)量的平均值［J］．齊齊哈爾師范學(xué)院學(xué)報(bào)(自然科學(xué)版)，1997，17( 4) : 28－31．

［6］王杰芳，李玉曉，賈瑜，等．量子物理學(xué)基本概念的教學(xué)研究［J］．洛陽(yáng)理工學(xué)院學(xué)報(bào)(自然科學(xué)版)，2010，20( 2) : 91－94．

［7］王慶領(lǐng)，張力．算符與量子力學(xué)中的力學(xué)量［J］．新疆師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，2002，21( 3) : 19－23．

［8］張永德．量子力學(xué)［M］．2版．北京:科學(xué)出版社，2008．

［9］張昆實(shí)．量子力學(xué)中的勢(shì)阱問(wèn)題解析［J］．高等函授學(xué)報(bào)(自然科學(xué)版)，2005，19( 6) : 17－19．

［10］林瓊桂．一維δ勢(shì)阱中的相對(duì)論粒子［J］．大學(xué)物理，2002，21( 2) : 15－19．

［11］何鴻斌．δ勢(shì)阱束縛態(tài)的一種求法［J］．遼寧師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，1989( 1) : 77－80．

［12］EUGEV Merzbacher．Quantum mechanics［M］．Second edition．New York: John Wiley＆Sons，Inc，1977: 82－85．

［13］陳鄂生，李明明．量子力學(xué)基礎(chǔ)教程［M］．6版．濟(jì)南:山東大學(xué)出版社，2014．

［14］錢伯初，曾謹(jǐn)言．量子力學(xué)習(xí)題精選與剖析［M］．3版．北京:科學(xué)出版社，2008．

［15］曾謹(jǐn)言．量子力學(xué)［M］．北京:科學(xué)出版社，1981．

(責(zé)任編輯:郎偉鋒)

Method for Calculating Average Value of Kinetic Energy in the Potential Well Function

LI Mingming，LIU Haiying
( School of Physics and Technology，University of Jinan，Jinan 250022，China)

Abstract:In quantum mechanics，the general formula of the average value for physical quantity F^is F= ∫ψ*F^ψdτ．However，it is not valid for the calculation of the kinetic energy of the particle in a δ-function potential well，since it gives a negative value for the kinetic energy，which is caused by the special method of solving the motion equation of the potential well．A new method to obtain kinetic energy is adopted by using the steady state energy obtained by the solution of the steady state equation minus the average potential energy，the function is T = E－∫ψ*Vψdτ．The average kinetic energy of particles in the ground state of the potential well V ( r) =－γδ( r－a) is studied．

Key words:δ－function potential well; average kinetic energy; Virial theorem

作者簡(jiǎn)介:李明明( 1965—)，女，山東壽光人，副教授，主要研究方向?yàn)槔碚撐锢恚?/p>

基金項(xiàng)目:濟(jì)南大學(xué)教研項(xiàng)目( W201402)

收稿日期:2015－05－10

DOI:10．3969/j．issn．1672－0032．2015．02．016

文章編號(hào):1672－0032( 2015) 02－0082－05

文獻(xiàn)標(biāo)志碼:A

中圖分類號(hào):O41