第二章 波函数和薛定谔方程b
第二章 波函数和薛定谔方程
§2.1 学习指导
本章主要介绍微观粒子运动状态的描述方法、演化规律以及由此带来的新特点,并以一维情况作例子进行具体说明。
根据实验,微观粒子具有波粒二象性。经典波一般用振幅A (r , t ) 与位相ϕ(r , t ) 来描述,
v v v i ϕ(r
它们可以统一写为ψ(r , t ) =A (r , t ) e , t ) ,在量子力学中沿用坐标与时间的复值函数
v v
v v
ψ(r , t ) 来描述微观粒子的运动状态,称为波函数。经典情况下,模方|ψ(r , t ) |2表示波的
强度;量子情况下,|ψ(r , t ) |2表示粒子出现的概率密度,因此需要把波函数归一化。
波函数随时间的变化由薛定谔方程确定。按照波函数的演化形式,粒子运动可以分为定态和非定态。在定态中,粒子的概率密度不随时间变化。按照定态波函数的空间形式,粒子运动可以分为束缚态和非束缚态。在束缚态中,粒子的能量取离散值,形成能级,可以很好地说明原子光谱。散射态是典型的非束缚态,可以用来描述粒子之间的碰撞,解释微观粒子的隧道贯穿现象。
真实的物理空间是三维的,但是当系统具有某些对称性时,可以约化为一维问题,例如中心势场中粒子的径向运动。近来,实验中也制备出了某些类型的一维量子力学系统。一维薛定谔方程容易求解,便于初学者理解量子力学的基本概念、熟悉常用方法和领会核心思想。
本章的主要知识点有 1. 微观粒子运动状态的描述 1)波函数
波函数ψ(r , t ) 是描述微观粒子状态的复值函数,波函数需要满足的标准条件为单值性、连续性和有界性。实际体系波函数满足平方可积条件,即2)波函数的意义
波函数的模方
v
v
⎰⎰⎰
v 2
ψ(r , t ) d τ=N 2
v v 2
w (r , t ) =ψ(r , t ) (2-1)
给出t 时刻粒子出现在位置r 邻域单位体积内的概率,即概率密度。
因此,标准的波函数应该是归一化的,即满足归一化条件
v
⎰⎰⎰
v 2
ψ(r , t ) d τ=1 (2-2)
未归一化的波函数可以通过乘以一个归一化因子来实现归一化。 3)波函数的性质
波函数ψ(r , t ) 满足叠加原理,如果ψi (r , t ), i =1,2, L 为微观粒子的可能状态,则
ψ(r , t ) =
也是一个可能的状态。 2. 微观状态的演化 1)薛定谔方程
状态ψ(r , t ) 随时间演化满足薛定谔方程 i h
v v
v
∑
i
i
v
c ψi (r , t ) , i c ∈ C (2-3)
v
∂v v
ψ(r , t ) =ˆH ψ(r , t ) (2-4) ∂t
2
v ˆ=-h ∇2+U (r , t ) (2-5) 其中 H
2m
称为哈密顿算符,U (r , t ) 是势能。若已知初始状态ψ(r ,0) ,由薛定谔方程可求出任意时刻t 的状态ψ(r , t ) 。 2)连续性方程
由薛定谔方程可以推出连续性方程
v v
v
v ∂w
+∇J ⋅=0 (2-6) ∂t v i h
(ψ*∇ψ-ψ∇ψ*) (2-7) 其中J =-2m
称为概率流密度,即沿着给定方向单位时间通过单位截面的概率。连续性方程是概率守恒定律的定域表现。 3)定态薛定谔方程
ˆ不显含时间,即势场U 不含t 时,薛定谔方程可以分离变量,得到 若体系的哈密顿H
v v -h Et 定态波函数解 ψE (r , t ) =ψE (r ) e (2-8)
其中E 为能量本征值,ψE (r ) 为对应的本征函数,满足定态薛定谔方程
i
v
h 22v v v v
∇ψE (r ) +U (r ) ψE (r ) =E ψE (r ) (2-9) -2m
3. 一维束缚定态问题 1)问题的描述
一维束缚定态问题由下面的方程和边界条件组成
⎧h 2d 2ψ(x )
+U (x ) ψ(x ) =E ψ(x ) ⎪-2
(2-10) 2m dx ⎨
x →∞⎪ψ(x ) →0⎩
其中束缚态能量满足条件E
束缚定态中的能量取值不连续,形成能级。同一能级只对应一个本征函数,无简并现象。第n 个能级E n , n ∈N 对应的本征函数ψn (x ) 有n 个内部零点(不包括边界)。
束缚态本征函数ψn (x ) 可以归一化,归一化后的本征函数满足正交归一性
⎰
∞
-∞
*
ψm (x ) ψn (x ) dx =δm , n (2-11)
本征函数集合具有完备性,任何平方可积函数ψ(x ) 都可以展开为归一化本征函数的线性组合,即
ψ(x ) =∑n c n ψn (x ) (2-12)
其中展开系数为
*c n =⎰ψn (x ) ψ(x ) dx (2-13)
-∞∞
3)典型实例:一维简谐振子
一般的解析势阱在其极小值附近都可以近似为简谐振子势,其标准形式为
U (x ) =kx =m ωx (2-14) 在上述势场中,粒子作束缚运动,能级为
1
2
1
22
E n =(n +1) h ω, n ∈N (2-15) 对应的本征函数为
ψn (x ) =N n e
1
-α2x 22
H n (αx ) (2-16)
其中H n (x
) 为厄密多项式,参数α,归一化系数
N n =简谐振子的本征函数满足递推关系
x ψn (x ) =
1
α
(x ) +n -1n +1(x )]
(2-17)
d
n (x ) =α(x ) -n -1dx
4. 一维散射问题 1)问题的描述
n +1(x )]
以能量E >U (±∞) 自左边向势场U (x ) 入射的粒子满足下面的方程和边界条件
⎧h 2d 2ψ(x )
+U (x ) ψ(x ) =E ψ(x ) ⎪-2
(2-18)
⎨2m dx
x →-∞x →∞⎪ψ(x ) −−−→Ae ikx +A ' e -ikx , ψ(x ) −−−→Ce ik ' x ⎩
其中k =
k =
2)问题的意义
在上面的问题中,入射波的概率流密度为J =h k |A |2/m ,反射波的概率流密度为
J R =-h k |A '|2/m ,透射波的概率流密度为J D =h k '|C |2/m 。由此得到反射系数R 和
透射系数D 分别为
J R |A ' |2J D k ' |C |2
(2-19) R =||=, D =||=⋅22
J |A |J k |A |
3)典型实例:粒子对方势垒的透射
能量为E 的粒子入射到一个宽度为a ,高度为U 0的方形势垒
⎧U 0, 00, x a ⎩
反射系数和透射系数分别为
22(k 12-k 2) sin 2k 2a 4k 12k 2
,D =2 (2-21) R =2
2222
(k 1-k 2) sin 2k 2a +4k 12k 2(k 1-k 2) sin 2k 2a +4k 12k 2
其中
k 1=
k 2=
§2.2 习题分析与求解
2.1 证明在定态中, 概率流密度与时间无关. 【题意分析】
v v -Et
已知条件:粒子处于定态,波函数为ψ(r , t ) =ψ(r ) e h (2.1-1) v
待证问题:概率流密度J 与时间无关;
相互联系:概率流密度与波函数之间具有关系(2-7) 【求解过程】
将定态波函数的一般形式(2.1-1)式代入概率流密度公式(2-7),得到
i
i i i i
v v -Et Et i h v -h Et v -h Et *v *v h h J (r , t ) =[ψ(r ) e ∇ψ(r ) e -ψ(r ) e ∇ψ(r ) e ]
2m (2.1-2)
v i h v v v v v
=[ψ(r ) ∇ψ*(r ) -ψ*(r ) ∇ψ(r )]=J (r ) 2m
容易看出,由上式得出的结果与时间无关。 又解:
定态波函数满足关系i h
∂t
v v *v *v ∂ψ(r , t ) =E ψ(r , t ) 和i h ψ(r , t ) =-E ψ(r , t ) ,因此有 v
=i h [∇ψ*+ψ∇*-*∇ψ-ψ*∇]2m (2.1-3)
1=[E ψ∇ψ*-E ψ∇ψ*+E ψ*∇ψ-E ψ*∇ψ]=02m
【物理讨论】
不能简单地由定态中概率密度w =|ψ(r ) |2不随时间变化,就推断概率流密度也不随时间变化,粒子流绕z 轴对称均匀加速转动就是一个相反的例子。定态中概率流密度不随时间变化有更深刻的原因。按(2-7)式,概率流密度可以变形为
v
v i h i h w ψh w
J =-(ψ*ψ∇ln ψ-ψψ*∇ln ψ*) =-∇ln *=∇ϕ (2.1-4)
2m 2m ψm
v v
其中ϕ=Arg ψ为波函数ψ的位相,即幅角。将上式与经典的粒子流密度J =wv 比较,量
子力学中的h ∇ϕ/m 对应于经典运动的速度。
(2.1-4)式表明概率流密度完全由波函数的位相决定,与波函数的模无关。在定态的情况下,波函数的位相ϕ=Arg ψ(r ) -Et /h ,因此∇ϕ=∇Arg ψ(r ) 与时间无关。这表明在定态中,概率流动的速度是稳定的。
也不能由定态中概率密度w =|ψ(r ) |2不随时间变化,就推断概率流密度为零,定向传播的平面波就是一个相反的例子。 2.2 由下列两定态波函数计算概率流密度:
(1) ψ1=
v v
v
1ikr 1
e , (2) ψ2=e -ikr r r
从所得结果说明ψ1表示向外传播的球面波, ψ2表示向内(即向原点) 传播的球面波. 【题意分析】
已知条件:粒子处于定态,定态波函数分别为ψ1=ψ1(r ) 和ψ2=ψ2(r ) ;
v v
待求问题:对应的概率流密度J (r ) ;
相互联系:概率流密度与定态波函数之间满足关系式
v v i h v v v v J (r ) =[ψ(r ) ∇ψ*(r ) -ψ*(r ) ∇ψ(r )] (2.2-1)
2m
【求解过程】
将定态波函数ψ1代入(2.2-1)式,利用梯度算符在球坐标中的表示形式(附录A ),得到
v v i h v v v v J 1(r ) =[ψ1(r ) ∇ψ1*(r ) -ψ1*(r ) ∇ψ1(r )]
2m
i h 1ikr v ∂1-ikr 1-ikr v ∂1ikr =[e e r (e ) -e e r (e )] (2.2-2) 2m r ∂r r r ∂r r i h 1ik 1ik v h k v =[(-2-) -(-2+)]e r =2e r 2mr r r r r mr
同理可得
v v i h h k v v v *v *v J 2(r ) =[ψ2(r ) ∇ψ2(r ) -ψ2(r ) ∇ψ2(r )]=-2e r (2.2-3)
2m mr
上式也可以通过在定态波函数的表达式ψ1=e ikr /r 中作变换k →-k 直接得到。 又解:
根据已知条件,定态波函数ψ1的模为u 1=1/r ,位相为ϕ1=kr ,代入简化后的概率流密度公式(2.1-4)中,立即得到
v v v h u 12(r ) h k v v
J 1(r ) =∇ϕ1(r ) =e r (2.2-4)
m mr 2
同理可计算出J 2(r ) 。 【物理讨论】
本题中,概率流密度与角度变量θ和ϕ无关,具有球对称性。J 1(r ) 与径向单位向量e r 同方向,表示向外传播的球面波;J 2(r ) 与径向单位向量e r 反方向,表示向内传播的球面波。
对于状态ψ1,单位时间通过球面r =a 向外传出的概率分别为
v v
v
v
v
v
I (a ) =乙⎰⎰
r =a
v v
J 1(r ) dS =
⎰⎰
r =a
h k v v h k 4πh k 2
e dS =⋅4πa = (2.2-5) 2r 2
mr ma m
这个概率值不随球面半径的大小变化,说明进入任意球壳层中的概率与流出的概率总是相等的,即任意球壳层中的概率不变。然而,对于半径任意小的球面,总是有概率向外流出,这表明在原点处有一个强度为4πh k /m 的概率源。同理,状态ψ2中在原点处有一个强度为
-4πh k /m 的概率源(即概率汇)。
2.3 一粒子在一维势场
⎧∞, x a
U (x ) =⎨
⎩0, 0≤x ≤a
【题意分析】
中运动,求粒子的能级和对应的波函数.
已知条件:粒子处于一维无限深方势阱U (x ) 中运动; 待求问题:粒子的能量本征值E n 和定态波函数ψn (x , t ) ; 相互联系:定态波函数ψn (x , t ) =ψn (x ) e 足定态薛定谔方程
-iE n t /h
,其空间部分ψn (x ) 和能量本征值E n 满
h 2d 2-(x ) +U (x ) ψ(x ) =E ψ(x ) (2.3-1) 2m dx 2
【求解过程】
因为势场U (x ) 是分段函数,本征函数ψ(x ) 也应分段考虑。在x a 区间内,
U (x ) =∞,而能量为有限值,定态薛定谔方程要求ψ(x ) =0;在x a 区间内,U (x ) =0,(2.3-1)成为
ψ''(x ) +k 2ψ(x ) =0, 0其中k (2.3-3)
方程(2.3-2)的通解为
ψ(x ) =A sin kx +B cos kx (2.3-4)
由波函数在x =0与x =a 处的连续性条件得到
ψ(0)=ψ(a ) =0 (2.3-5)
将通解(2.3-4)代入条件(2.3-5),有
ψ(0)=A sin 0+B cos0=0
ψ(a ) =A sin ka +B cos ka =0
由上面第一式得到B =0,代入第二式后解出
k =n π/a , n =1,2, L (2.3-6)
将(2.3-6)式代入(2.3-3)式,得到能量本征值
E =h 2k 2/(2m ) =h 2n 2π2/(2ma 2), n =1,2, L (2.3-7)
将(2.3-6)式代入 (2.3-4)式中,得到本征函数
⎧A sin(n πx /a ), 0(2.3-8) ψn (x ) =⎨
0, x a ⎩
其中常数A 由归一化条件确定,即
∞
a
⎰
-∞
|ψ(x ) |2dx =⎰|ϕn (x ) |2dx =1
由此得到归一化系数A =
-iE n t /h
定态波函数为ψn (x , t ) =ψn (x ) e 又解:
。
对于宽度为2b 的对称一维无限深方势阱
⎧0, |x |≤b
(2.3-9) U (x ) =⎨
⎩∞, |x |>b
在该势场中运动粒子的能级为(见参考文献【1】§2.6节)
E n (b ) =
对应的本征函数为
π2h 2n 2
8mb 2
, n =1, 2, L (2.3-10)
n π(x +b ), |x |≤b (2.3-11) ϕn (x , b ) =2b
0, |x |>b ⎩
本题中研究的是宽度为a 的一维无限深方势阱,如果在上述势阱中取2b = a ,并把势阱的位置向x 轴正向平移b ,就成为本题的情况。由于能量本征值的大小与势阱的位置无关,因此在能级表达式(2.3-10)中取2b = a ,就得到本题情况下的能级
E n =E (a ) =
1
n π2h 2n 2
2ma 2
, n =1, 2, L (2.3-12)
将表达式(2.3-11)中的本征函数向x 轴正向平移b ,再取2b = a ,就得到本题情况下的本征函数
n πx , |x -1a |≤1a a ψn (x ) =ϕn (x -2a , 2a ) = (2.3-13)
0, |x -⎩2a |>2a
三解:
注意到本题中势阱恰好是宽度为2a 的对称一维无限深方势阱的一半,对称一维无限深方势阱中的奇宇称本征函数
ϕ2n (x ) =A sin
n πn π(x +a ) =Ae in πsin x , n =1, 2, L (2.3-14) a a
恰好满足本题中势阱内本征函数所要求的方程和边界条件
⎧ψ''(x ) +k 2ψ(x ) =0, 0(2.3-15) ⎨
⎩ψ(0)=ψ(a ) =0
因此本题的解为
ψn (x ) =Ae in πsin
n π
x , 0将上面的本征函数代入定态薛定谔方程(2.3-2)式后,立刻能量本征值为
d 2ψn (x ) π2h 2n 2h 2
(2.3-17) E n =-=22
2m ψn (x ) dx 2ma
【物理讨论】
第一解中在得到系数B =0后,另一个系数A ≠0,否则波函数在全空间中都等于零,对应的粒子不存在。同样,量子数n 也不能等于零,因此粒子的基态为n = 1。
由量子化条件(1-2),得到
Ñ⎰pdq =2由此解出
=nh , n ∈N (2.3-18)
h 2n 2π2h 2n 2
E n ==, n =0,1,2, L (2.3-19)
8ma 22ma 2
两者相差1个能级。
此外,按照量子力学的结果,粒子的概率分布为
π2
w n (x ) =|ψn (x ) |2=sin 2n x , 0但是按照经典力学,在(x , x +dx ) 区间内找到粒子的概率w (x ) dx 与该粒子在此区间内逗留的时间成正比,即w (x ) dx =dt /T ,其中T 为粒子运动的周期。考虑到在一个周期中,粒子两次经过同一个位置,于是得到概率密度为
w (x ) =
2dt 2
= (2.3-21) Tdx Tv
在本题中,T =2a /v ,于是得到w (x ) =1/a ,为一个与位置和能量都无关的常数。 2.4 证明宽度为2a 的一维对称无限深方势阱中本征函数ψ
n (x ) 的归一化因子是A ' = 【题意分析】
n π
⎧⎪A 'sin (x +a ), |x |
已知条件:本征函数为ψn (x ) =⎨ (2.4-1)
0, |x |≥a ⎪⎩
待求问题:归一化因子A ' ; 相互联系:归一化条件
⎰
∞
-∞
|ψ(x ) |2dx =1 (2.4-2)
【求解过程】
将本征函数(2.4-1)式代入归一化条件,得到
⎰
∞
-∞
'2n d x =⎰|A |s i 2a -
2
a
n π
2a
+x (
2
a ) d =x '|A =|a
1
由此得到归一化因子A '=又解:
i ϕ
。 从物理上看,一维无限深方势阱中运动粒子的能级和归一化因子完全由势阱的宽度决定,与势阱在x 轴上所处的位置无关。由上题,势阱宽度为a
因此当势阱宽度为2a
= 【物理讨论】
一般来说,归一化因子可以是复数,其幅角部分称为相因子。本题中归一化因子的模
相因子为e ,其中幅角ϕ可以取任意实数。相因子中幅角的取值既不影响系统的概率密度和概率流密度,也不影响能量、动量等物理量,没有可以观察的物理效应。因此,为了简单起见,通常忽略归一化因子中的相因子,即把幅角取为零。这样,本题中的归一化因
子成为A ' =1/
在本题中,所有本征函数的归一化因子都相同,这是一个非常特殊的情况。在一般情况下,归一化因子随量子数n 的变化而改变。 2.5 求一维谐振子处在第一激发态时概率最大的位置。 【题意分析】
已知条件:一维谐振子第一激发态的波函数
1-αx ψ(x ) =2αxe 2 (2.5-1)
22
i ϕ
待求问题:求x m ,使得w (x m ) =max x ∈R w (x ) ; 相互联系:w (x ) =|ψ(x ) |。
2
【求解过程】
由定义,概率密度为
w (x ) =|ψ(x ) |=
2
3
x 2e -α
22
x
(2.5-2)
满足束缚态条件w (±∞) =0。由极值条件
w '(x ) =
3
(2x -2αx ) e
23
-α2x 2
=0 (2.5-3)
x =±1α/。容易验证,在驻点处有w "(0)>0,得到概率密度函数的驻点为x =0,
w "(±1/α)
又解:
在物理学中常常考察物理量的相对变化率y '(x ) /y (x ) ,它往往比变化率y '(x ) 更有意义。相对变化率可以表示成对数导数的形式,即y '(x ) /y (x ) =[lny (x )]',计算很方便。从数学的角度看,对数函数是单调增函数,能够保持函数的增减性,不改变极值点的性质和位置。在本题中我们也可以用概率密度的对数导数来进行分析,即计算
[lnw (x )]'=[ln
3
2ln x -α2x 2]'=
2
-2α2x =0 (2.5-3) x
立刻得到x =±1/α。 【物理讨论】
数学上可以证明,在一维势场中运动的束缚态粒子,基态在区域内部没有零点,第n 个激发态有n 个零点(不包括边界上的两个零点),这个结论称为零点定理。概率密度函数为波函数的模方,因此也有n 个零点,加上边界上的2个零点,一共有n+2个零点,即最小值。连续函数的两个相邻最小值之间有一个极大值,因此第n 个能量本征态的概率密度具有n+1个极大值。由下题可知,一维束缚态的本征函数具有确定的宇称,对应的概率密度为偶函数,其极大值的分布关于y 轴对称。在本题的情况下,n = 1,因此概率密度有2个极大值,对称地分布在原点的两侧。
2.6 在一维势场中运动的粒子,势能对原点对称,证明粒子的定态波函数具有确定的宇称. 【题意分析】
已知条件:粒子在一维势场U (x ) 中运动,U (-x ) =U (x ) 。 待证问题:本征函数具有确定的宇称,即ψ(-x ) =±ψ(x ) 相互联系:本征函数满足定态薛定谔方程
h 2d 2
ψ(x ) +U (x ) ψ(x ) =E ψ(x ) (2.6-1) -
2m dx 2
【求解过程】
对定态薛定谔方程进行空间反演,即把x 换成-x ,得到
h 2d 2
ψ(-x ) +U (x ) ψ(-x ) =E ψ(-x ) (2.6-2) -2
2m dx
计算中已经利用了势能的对称性。与(2.6-1)式比较后,可以看出ψ(-x ) 也是定态薛定谔方程的一个解,即也是本征函数,与ψ(x ) 对应于同一个本征值E 。 如果ψ(-x ) 与ψ(x ) 线性相关,即
ψ(-x ) =c ψ(x ) (2.6-3) 在上式中再把x 换成-x ,得到
ψ(x ) =c ψ(-x ) =c 2ψ(x ) (2.6-4)
2
由于波函数不能恒等于零,上式给出了条件c =1,即c =±1。代入到(2.6-3)式后,得到
( ) (2.6-5) ψ(-x ) =±ψx
如果ψ(-x ) 与ψ(x ) 线性无关,则可以线性组合成两个新的独立本征函数
ψe (x ) =ψ(x ) +ψ(-x )
(2.6-6)
ψo (x ) =ψ(x ) -ψ(-x )
容易验证
ψe (-x ) =+ψe (x )
(2.6-7)
ψo (-x ) =-ψo (x )
由此说明了粒子的定态波函数具有确定的宇称。 【物理讨论】
在一维束缚态的情况下,波函数不存在简并现象,这时ψ(-x ) 与ψ(x ) 描写了同一个能量本征态,它们之间线性相关,故ψ(x ) 只能是奇函数或偶函数。由于第n 个能级的波函数有n 个零点,因此在对称势阱中,与偶数能级对应的波函数是偶函数;与奇数能级对应的波函数是奇函数。
在一维散射态的情况下,波函数存在简并现象,这时ψ(-x ) 与ψ(x ) 可以描写同一个能
量的两个不同的状态。它们能够组合成奇本征函数和偶本征函数,分别有确定的宇称。 2.7 一粒子在一维势阱
⎧U >0, |x |>a ,
U (x ) =⎨0
|x |≤a ⎩0,
中运动,求束缚态(0
已知条件:粒子处于一维有限深对称方势阱U (x ) 中运动; 待求问题:粒子的束缚态能级E n ; 相互联系:定态薛定谔方程(1-9)式 【求解过程】
因为势场U (x ) 是不连续的分段函数,我们将待求的本征函数ψ(x ) 也表示为分段形式
⎧ϕ-(x ), x
⎪
ψ(x ) =⎨ϕ(x ), |x |≤a (2.7-1)
⎪ϕ(x ), x >a ⎩+
将(2.7-1)式代入定态薛定谔方程(1-9)中,得到
⎧h 2d 2
⎪-2m dx 2ϕ-(x ) +U 0ϕ-(x ) =E ϕ-(x ), x
ϕ(x ) =E ϕ(x ), |x |≤a (2.7-2) ⎨-2
2m dx ⎪
⎪h 2d 2
ϕ(x ) +U 0ϕ+(x ) =E ϕ+(x ), x >a ⎪-2+
2m dx ⎩
由于束缚态能量满足条件0
k =κ= (2.7-3)
方程(2.7-2)可以简化为
⎧ϕ-''(x ) -κ2ϕ-(x ) =0, x
⎨ϕ''(x ) +k ϕ(x ) =0, |x |≤a (2.7-4) ⎪ϕ''(x ) -κ2ϕ(x ) =0, x >a
+⎩+
由此解出
⎧ϕ-=Ae κx +Be -κx , x
⎪
⎨ϕ=C sin kx +D cos kx , |x |≤a (2.7-5) ⎪ϕ=Fe κx +Ge -κx , x >a ⎩+
考虑到束缚态波函数满足无穷远边界条件ψ(±∞) =0,上式中的系数B 和F 都必须等于零。在x =-a 处,波函数满足连接条件ϕ-(x ) =ϕ(x ) 和ϕ-'(x ) =ϕ'(x ) ;在x =a 处,满足连接条件ϕ+(x ) =ϕ(x ) 和ϕ+'(x ) =ϕ'(x ) ,于是得到四个关系式
Ae -κa =-C sin ka +D cos ka
κAe -κa =kC cos ka +kD sin ka
Ge
-κa
=C sin ka +D cos ka
(2.7-6)
-κGe -κa =kC cos ka -kD sin ka
上式可以化为矩阵形式
⎛e -κa -κa κe 0 0⎝e -κa κe -κa
00
⎫⎛A ⎫⎪ ⎪
-k cos ka -k sin ka 0⎪ C ⎪
=0 (2.7-7) -κa ⎪ ⎪D -sin ka -cos ka e
⎪ ⎪
-k cos ka k sin ka -κe -κa ⎪⎭⎝G ⎭sin ka
-cos ka
sin ka -cos ka -k cos ka -k sin ka -sin ka -k cos ka
-cos ka k sin ka
e 00
-κa
由于波函数不能等于零,上面的方程应该有非零解,这要求系数行列式等于零,即
=0 (2.7-8)
-κe -κa
经过仔细的计算,得到
(κ2-k 2)sin ka cos ka +k κ(cos2ka -sin 2ka ) =0 (2.7-9)
上式可以简化为
k 2-κ2=2k κcot 2ka (2.7-10)
这就是束缚态能量所必须满足的条件,将关系(2.7-3)代入(2.7-10)式后即可确定能量本征值。 又解:
方程(2.7-4)的解又可以表示为
⎧ϕ-=Ae κx +Be -κx , x
⎪
⎨ϕ=C sin(kx +δ), |x |≤a (2.7-11) ⎪ϕ=Fe κx +Ge -κx , x >a ⎩+
无穷远边界条件要求系数B 和F 都必须等于零。在x =a 处,波函数满足连接条件
ϕ+(x ) =ϕ(x ) 和ϕ+'(x ) =ϕ'(x ) ,当ϕ(a ) ≠0时,这两个条件可以归结为对数导数连接条件
(lnϕ+)' =(lnϕ)' ;同理,在x =-a 处,满足连接条件(lnϕ+)' =(lnϕ)' ,于是得到二个关
系式
-κ=k cot(ka +δ)
(2.7-12)
κ=k cos(-ka +δ)
上式为束缚态能量所必须满足的条件。 三解:
由于势阱具有对称性,因此束缚态本征函数具有确定的宇称。这样,我们不需要在全空间对薛定谔方程求解,只要考虑x >0时的情况。这时,方程(2.7-4)可以简化为
2
⎧⎪ϕ''(x ) +k ϕ(x ) =0, 0
(2.7-13) ⎨2
⎪⎩ϕ+''(x ) -κϕ+(x ) =0, x >a
在奇宇称的情况下,ϕ(0)=0;结合无穷远条件ϕ(∞) =0,得到解函数
⎧ϕ=C sin kx , 0
(2.7-14) ⎨-κx
x >a ⎩ϕ+=Ge ,
由x =a 处的对数导数连接条件(lnϕ+)' =(lnϕ)' ,得到确定奇宇称能级的关系式
-κ=k cot ka (2.7-15)
在偶宇称的情况下,ϕ'(0)=0;结合无穷远条件ϕ(∞) =0,得到解函数
⎧ϕ=D cos kx , 0
(2.7-16) ⎨-κx
ϕ=Ge , x >a ⎩+
由x =a 处的连接条件(lnϕ+)' =(lnϕ)' ,得到确定偶宇称能级的关系式
-κ=-k tan ka (2.7-17)
【物理讨论】
在第二解(2.7-12)式中,利用正切函数的和角公式cot(x +y ) =去参数δ,得到
cot x cot y -1
,可以消
cot x +cot y
κ2-k 2
cot 2ka =cot[(ka +δ) +(ka -δ)]= (2.7-18)
-2k κ
这正是第一解中的结果。
将第一解中得到的束缚态能量条件(2.7-9)因式分解为
cos ka sin ka (κ-k tan ka )(κ+k cot ka ) =0 (2.7-19)
可见与第三解的结果(2.7-15)和(2.7-17)两式等价,但是第三解的物理意义更明显。
为了便于具体考察束缚态能级的性质,我们定义无量纲变量u =
,
ξ=ka ,于是有
η=κa = (2.7-20)
(2.7-17)式成为
η=ξtan ξ (2.7-21)
联立(2.7-20)与(2.7-21)
式就可以确定ξ=
的数值,这可以通过作图法得到。例如取
u =10,利用Mathematica 命令
Plot[{Sqrt[u^2-ξ^2],ξ Tan[ξ]},{ξ,0,u}] 得到图形
图中出现了4个交点,对应4个能级。对偶宇称情况,当n π
x
⎪U , 0≤x
U (x ) =⎨
⎪-U 1, a ≤x ≤b ⎪b
求束缚态的能级所满足的方程. 【题意分析】
已知条件:粒子处于一维势阱U (x ) 中运动;
待求问题:粒子的束缚态能量本征值E 所满足的方程; 相互联系:定态薛定谔方程(2-9) 。 【求解过程】
因为势场U (x ) 是不连续的分段函数,我们将本征函数ψ(x ) 也表示为分段形式。考虑到在区间x
⎧ϕ0(x ), 0≤x
ψ(x ) =⎨ϕ1(x ), a ≤x ≤b (2.8-1)
⎪ϕ(x ), b
将(2.8-1)式代入定态薛定谔方程(2-9)中,得到
⎧h 2d 2
⎪-2m dx 2ϕ0(x ) +U 0ϕ0(x ) =E ϕ0(x ), 0≤x
ϕ(x ) -U 1ϕ1(x ) =E ϕ1(x ), a ≤x ≤b (2.8-2) ⎨-21
2m dx ⎪
⎪h 2d 2
ϕ(x ) =E ϕ2(x ), b
2m dx ⎩
由于束缚态能量E
κ0=k =κ2= (2.8-3)
方程(2.8-2)可以简化为
⎧ϕ0''(x ) -κ2ϕ0(x ) =0, 0≤x
⎨ϕ1''(x ) +k ϕ1(x ) =0, a ≤x ≤b (2.8-4) ⎪2ϕ''(x ) +κϕ2(x ) =0, b
由此解出
⎧ϕ0=Ae κ0x +Be -κ0x , 0≤x
⎨ϕ1=C sin kx +D cos kx , a ≤x ≤b (2.8-5) ⎪κ2x -κ2x
b
考虑到无穷远边界条件ψ(∞) =0,上式中的系数F 必须等于零;在x =0处,波函数的连续性要求ϕ0(x ) =0,得到系数关系A +B =0。而在x =a 处,波函数满足连接条件
ϕ0(x ) =ϕ1(x ) 和ϕ0'(x ) =ϕ1'(x ) ;在x =b 处,满足连接条件ϕ2(x ) =ϕ1(x ) 和
ϕ2'(x ) =ϕ1'(x ) ,又得到四个关系式
Ae κ0a +Be -κ0a =C sin ka +D cos ka
κ0Ae κa -κ0Be -κa =kC cos ka -kD sin ka
Ge
-κ2b
=C sin kb +D cos kb
(2.8-6)
-κ2Ge -κ2b =kC cos kb -kD sin kb
利用系数关系B =-A ,上式可以化为
-sin ka -cos ka 0⎫⎛A ⎫⎛2sinh κ0a
⎪⎪2κcosh κa -k cos ka k sin ka 000 ⎪C ⎪=0 (2.8-7) 0-sin kb -cos kb e -κ2b ⎪D ⎪ ⎪-κ2b ⎪0-k cos kb k sin kb -κe ⎝2⎭⎝G ⎭
由于波函数不能等于零,这要求系数行列式等于零,经过耐心仔细的计算,得到
tan k (b -a ) =
k κ0coth κ0a +k κ2
(2.8-8)
k 2-κ0κ2coth κ0a
上式为束缚态能量所必须满足的条件,与(2.8-4)式联立后即可确定能量本征值。 二解:
方程(2.8-4)的解又可以写为下列形式
⎧ϕ0=A sinh κ0x +B cosh κ0x , 0≤x
a ≤x ≤b (2.8-9) ⎨ϕ1=C sin(kx +δ),
⎪κ2x -κ2x
b
考虑到无穷远边界条件,系数F =0;在x =0处,由连续性条件得到系数B =0。利用波函数连接条件的对数导数形式,即在x =a 处,波函数满足连接条件(lnϕ0)' =(lnϕ1)' ;在x =b 处,满足连接条件(lnϕ1)' =(lnϕ2)' ,得到二个关系式
κ0coth(κ0a ) =k cot(ka +δ), k cot(kb +δ) =-κ2 (2.8-10)
上式结合关系(2.8-3)可以确定能量本征值。 【物理讨论】
从表面看,两种解法的结果不同。利用差角公式tan(x -y ) =去第二解(2.8-11)式中的参数δ,得到
tan x -tan y
,可以消
1+tan x tan y
tan k (b -a ) =tan[(kb +δ) -(ka +δ)]=-
这正是第一解所得到的公式(2.8-8)。
k κ0+k κ2tanh(κ0a )
(2.8-11) 2
κ0κ2-k tanh(κ0a )
§2.3 扩展练习
E2-1 设ψ1(r , t ) 和ψ2(r , t ) 是体系的两个可能的运动状态,现由ψ1和ψ2构成以下波函数
c ψ+c ψ22
,ψb =c 1ψ1+c 2ψ2,ψc =e ψa =c 1ψ1+c 2ψ2
11
22
v v
其中c 1, c 2为常数。问以上这些波函数能否描述体系的状态?为什么?
【提示】ψb 是ψ1, ψ2的线性叠加态,满足叠加原理,ψa 和ψc 所表示的状态不满足态叠加原理。
E2.2 求线性谐振子处于第n 个激发态时,在经典界限外被发现的概率。
【提示】经典力学要求粒子的动能T =E n -U (x ) ≥0,由此得到经典运动范围为
E n -m ω
2
22
x ≥
0,即|x |≤x m =
。在经典界限内的概率为⎰
x m
-x m
|ψn |2dx 。
E2.3粒子在一维势场中作束缚运动,其基态波函数为ψ=A cosh -μkx , μ>0,求对应的基态能量和势能U (x ) 。
【提示】由定态薛定谔方程-2ψ=E ψ,得到 m ψ'' +U (x )
2ψ'' 2
U (x ) =+E =[μ(μ+1) k 2tanh 2kx -μk 2]+E (E2.3-1) 2m ψ2m
2
取无穷远为势能零点,即
2
U (∞) =[μ(μ+1) k 2-μk 2]+E =0 (E2.3-2) 2m
得到
22E =-μ2 (E2.3-3) 代回(E2.3-1)式,得到势能为
2222μ(μ+1) 2
U (x ) =[μ(μ+1) k 2tanh 2kx -μk 2]-μ2=- (E2.3-4)
2m 2m 2m cosh 2kx
E2.4 粒子在宽度为π的一维无限深势阱中运动,在t =0时刻的波函数为ψ(x ,0) =A sin 3x ,求状态随时间的演化规律。
h
2n 2, n ∈Z +,对应的本征函数为ψn =
【提示】能量本征值为E n =nx ,因此2m ψ(x ,0) =4A (3sinx -sin 3x ) =4ψ1-
ψ3) 。由归一化条件得到A =
ψ(x ,0) =ψ(x , t ) x -sin3x ) =
ψ1-ψ3) 。由定态波函数的性质和叠加原理,得到
ψ1e -iE 1t /h -ψ3e -iE 3t /h )
⎧0, x
E2.5 能量为E 的粒子从左边向势垒U (x ) =⎨U 1,0≤x ≤a 运动,求透射系数。
⎪U , a
【答案】当E U 2时,透射系数
D =
4k 2/k 1
(E2.5-1)
(1+k 2/k ) 2cos 2k 1a +(k 1/k +k 2/k 1) 2sin 2k 1a
2
其中k 2=2mE /h 2, k 12=2m (E -U 1)/h 2, k 2=2m (E -U 2)/h 2.
E2.6 设粒子处于二维无限深势阱V (x , y ) =⎨
⎧00,求粒子能量和相应
其它区域⎩∞
的本征态。如a =b ,试讨论前5条能级简并情况。 【提示】用分离变量法,得到能量本征值和对应的本征函数
(x )
E n , l =E n +E l (y ) , ψn , l (x , y ) =ψn (x ) φl (y ), n , l ∈Z + (E2.6-1)
其中
E
(x ) n
n πx (y ) π2h 2l 2l πy
(E2.6-2) =, ψ(x ) =; E =, φ(y ) =n l l 22
2ma a 2mb b
E 1=
π2h 2n 2
当a =b 时,E 1,1=2E 1, E 1,2=E 2,1=5E 1, E 2,2=8E 1, E 1,3=E 3,1=10E 1,
π2h 2
2ma 2
。
E2.7质量为m 的粒子在一维势场U (x ) =U 0tan 2(kx ) 中运动,分别就U 0很大和很小两种情况,估算粒子的前几个能级的能量E n ,并与严格解比较。
【提示】当U 0很小时,势场变为宽度为π/k 的一维无限深势阱,能级为
h 2π2n 2h 2k 2n 2E n ; =, n =1,2, L (E2.7-1) 22
2m π/k 2m
当U 0很大时,低能级中的粒子只能在原点附近作微振动,这时势场可以近似表示为
22
U (x ) ≈U 0k 2x 2=m ωx (E2.7-2)
2
其中ω=k /m ,V 0=2mU 0/h 。由此得到
2
E n ≈h ω(n =0,1, L (E2.7-3) 2+n ),
严格解为
h 2k 21
E n =(n +4λn -2λ) , λ=n =1,2, L (E2.7-4) 4, 2m 2