对于一个Delta位势阱的散射 。往左与往右的行进波 的振幅与方向都分别表示于图内。用来计算透射系数 与反射系数 的行进波 都以红色表示。
在量子力学 里,Delta位势阱 是一个阱 内位势为负狄拉克Delta函数 ,阱外位势为0的位势阱。Delta位势阱问题 专门研讨,在这种位势的作用中,一个粒子的量子行为。这是一个常见的理论问题。假若,粒子的能量是正值的,我们想要知道的是,在被Delta位势垒散射 的状况下,粒子的反射系数 与透射系数 。假若,粒子的能量是负值的,这粒子会被束缚于Delta位势阱的阱内。这时,我们想要知道的是粒子的能量与束缚的量子态。
一个粒子独立于时间 的薛定谔方程 为
−
ℏ
2
2
m
d
2
ψ
(
x
)
d
x
2
+
V
(
x
)
ψ
(
x
)
=
E
ψ
(
x
)
{\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}\psi (x)}{dx^{2}}}+V(x)\psi (x)=E\psi (x)\,\!}
;
其中,
ℏ
{\displaystyle \hbar \,\!}
是约化普朗克常数 ,
m
{\displaystyle m\,\!}
是粒子质量,
x
{\displaystyle x\,\!}
是粒子位置,
E
{\displaystyle E\,\!}
是能量,
ψ
(
x
)
{\displaystyle \psi (x)\,\!}
是波函数 ,
V
(
x
)
{\displaystyle V(x)\,\!}
是位势,表达为
V
(
x
)
=
−
λ
δ
(
x
)
{\displaystyle V(x)=-\lambda \delta (x)\,\!}
;
其中,
δ
(
x
)
{\displaystyle \delta (x)\,\!}
是狄拉克Delta函数 ,
λ
{\displaystyle \lambda \,\!}
是狄拉克Delta函数的强度。
这位势阱将一维空间分为两个区域:
x
<
0
{\displaystyle x<0\,\!}
与
x
>
0
{\displaystyle x>0\,\!}
。在任何一个区域内,位势为常数,薛定谔方程的解答可以写为往右与往左传播的波函数的的叠加 (参阅自由粒子 ):
ψ
L
(
x
)
=
A
r
e
i
k
x
+
A
l
e
−
i
k
x
x
<
0
{\displaystyle \psi _{L}(x)=A_{r}e^{ikx}+A_{l}e^{-ikx}\quad x<0\,\!}
,
ψ
R
(
x
)
=
B
r
e
i
k
x
+
B
l
e
−
i
k
x
x
>
0
{\displaystyle \psi _{R}(x)=B_{r}e^{ikx}+B_{l}e^{-ikx}\quad x>0\,\!}
;
其中,
A
r
{\displaystyle A_{r}\,\!}
、
A
l
{\displaystyle A_{l}\,\!}
、
B
r
{\displaystyle B_{r}\,\!}
、
B
l
{\displaystyle B_{l}\,\!}
都是必须由边界条件 决定的常数,下标
r
{\displaystyle r\,\!}
与
l
{\displaystyle l\,\!}
分别标记波函数往右或往左的方向。
k
=
2
m
E
/
ℏ
2
{\displaystyle k={\sqrt {2mE/\hbar ^{2}}}\,\!}
是波数 。
当
E
>
0
{\displaystyle E>0\,\!}
时,
ψ
L
{\displaystyle \psi _{L}\,\!}
与
ψ
R
{\displaystyle \psi _{R}\,\!}
都是行进波 。可是,当
E
<
0
{\displaystyle E<0\,\!}
时,
ψ
L
{\displaystyle \psi _{L}\,\!}
与
ψ
R
{\displaystyle \psi _{R}\,\!}
都随着坐标
x
{\displaystyle x\,\!}
呈指数递减或指数递增。
在
x
=
0
{\displaystyle x=0\,\!}
处,边界条件是:
ψ
L
=
ψ
R
{\displaystyle \psi _{L}=\psi _{R}\,\!}
,
d
d
x
ψ
L
=
d
d
x
ψ
R
−
2
m
λ
ℏ
2
ψ
R
{\displaystyle {\frac {d}{dx}}\psi _{L}={\frac {d}{dx}}\psi _{R}-{\frac {2m\lambda }{\hbar ^{2}}}\psi _{R}\,\!}
。
特别注意第二个边界条件方程,波函数随位置的导数在
x
=
0
{\displaystyle x=0\,\!}
并不是连续的,在位势阱两边的差额有
2
λ
ℏ
2
ψ
R
{\displaystyle {\frac {2\lambda }{\hbar ^{2}}}\psi _{R}\,\!}
这么多。这方程的推导必须用到薛定谔方程。将薛定谔方程积分于
x
=
0
{\displaystyle x=0\,\!}
的一个非常小的邻域:
−
ℏ
2
2
m
∫
−
ϵ
ϵ
d
2
ψ
d
x
2
d
x
+
∫
−
ϵ
ϵ
V
(
x
)
ψ
d
x
=
E
∫
−
ϵ
ϵ
ψ
d
x
{\displaystyle -{\frac {\hbar ^{2}}{2m}}\int _{-\epsilon }^{\epsilon }{\frac {d^{2}\psi }{dx^{2}}}\,dx+\int _{-\epsilon }^{\epsilon }V(x)\psi \,dx=E\int _{-\epsilon }^{\epsilon }\psi \,dx\,\!}
;(1)
其中,
ϵ
{\displaystyle \epsilon \,\!}
是一个非常小的数值。
方程(1)右边的能量项目是
E
∫
−
ϵ
ϵ
ψ
d
x
≈
E
⋅
2
ϵ
⋅
ψ
(
0
)
{\displaystyle E\int _{-\epsilon }^{\epsilon }\psi \,dx\approx E\cdot 2\epsilon \cdot \psi (0)\,\!}
。(2)
当
ϵ
→
0
{\displaystyle \epsilon \to 0\,\!}
时,该项趋向于0。
方程(1)左边是
−
ℏ
2
2
m
(
d
ψ
R
d
x
|
ϵ
−
d
ψ
L
d
x
|
−
ϵ
)
+
λ
∫
−
ϵ
ϵ
δ
(
x
)
ψ
d
x
=
0
{\displaystyle -{\frac {\hbar ^{2}}{2m}}\left({\frac {d\psi _{R}}{dx}}{\bigg |}_{\epsilon }-{\frac {d\psi _{L}}{dx}}{\bigg |}_{-\epsilon }\right)+\lambda \int _{-\epsilon }^{\epsilon }\delta (x)\psi \,dx=0\,\!}
(3)
根据狄拉克Delta函数 的定义,
∫
−
ϵ
ϵ
δ
(
x
)
ψ
d
x
=
ψ
R
(
0
)
{\displaystyle \int _{-\epsilon }^{\epsilon }\delta (x)\psi \,dx=\psi _{R}(0)\,\!}
。(4)
而在
ϵ
→
0
{\displaystyle \epsilon \to 0\,\!}
的极限,
lim
ϵ
→
0
d
ψ
L
d
x
|
−
ϵ
=
d
ψ
L
d
x
|
0
{\displaystyle \lim _{\epsilon \to 0}{\frac {d\psi _{L}}{dx}}{\bigg |}_{-\epsilon }={\frac {d\psi _{L}}{dx}}{\bigg |}_{0}\,\!}
,(5)
lim
ϵ
→
0
d
ψ
R
d
x
|
ϵ
=
d
ψ
R
d
x
|
0
{\displaystyle \lim _{\epsilon \to 0}{\frac {d\psi _{R}}{dx}}{\bigg |}_{\epsilon }={\frac {d\psi _{R}}{dx}}{\bigg |}_{0}\,\!}
。(6)
将这些结果(4),(5),(6)代入方程(3),整理后,可以得到第二个边界条件方程:在
x
=
0
{\displaystyle x=0\,\!}
,
d
ψ
L
d
x
=
d
ψ
R
d
x
−
2
m
λ
ℏ
2
ψ
R
{\displaystyle {\frac {d\psi _{L}}{dx}}={\frac {d\psi _{R}}{dx}}-{\frac {2m\lambda }{\hbar ^{2}}}\psi _{R}\,\!}
。
从这两个边界条件方程。稍加运算,可以得到以下方程:
A
r
+
A
l
=
B
r
+
B
l
{\displaystyle A_{r}+A_{l}=B_{r}+B_{l}\,\!}
,
i
k
(
A
r
−
A
l
−
B
r
+
B
l
)
=
2
m
λ
ℏ
2
(
B
r
+
B
l
)
{\displaystyle ik(A_{r}-A_{l}-B_{r}+B_{l})={\frac {2m\lambda }{\hbar ^{2}}}(B_{r}+B_{l})\,\!}
。
一个Delta位势阱的反射系数
R
{\displaystyle R\,\!}
(用红线表示)与透射系数
T
{\displaystyle T\,\!}
(用绿线表示)随着能量
E
{\displaystyle E\,\!}
的变化。在这里,能量
E
>
0
{\displaystyle E>0\,\!}
。能量的单位是
λ
2
2
m
ℏ
2
{\displaystyle {\frac {\lambda ^{2}}{2m\hbar ^{2}}}\,\!}
。经典力学的答案用虚线表示,量子力学的答案用实线表示。
假若,能量是正值的,粒子可以自由的移动于位势阱外的两个半空间,
x
<
0
{\displaystyle x<0\,\!}
或
x
>
0
{\displaystyle x>0\,\!}
。在这里,粒子的量子行为主要是由Delta位势阱造成的散射 行为。称这粒子的量子态 为散射态 。设定粒子从左边入射。在Delta位势阱,粒子可能会被反射回去,或者会被透射过去。我们想要知道散射的反射系数 与透射系数 。设定
A
r
=
1
{\displaystyle A_{r}=1\,\!}
,
A
l
=
r
{\displaystyle A_{l}=r\,\!}
,
B
l
=
0
{\displaystyle B_{l}=0\,\!}
,
B
r
=
t
{\displaystyle B_{r}=t\,\!}
。求算反射的概率幅
r
{\displaystyle r\,\!}
与透射的概率幅
t
{\displaystyle t\,\!}
:
r
=
−
1
i
ℏ
2
k
m
λ
+
1
{\displaystyle r=-\ {\cfrac {1}{{\cfrac {i\hbar ^{2}k}{m\lambda }}+1}}\,\!}
,
t
=
1
−
i
m
λ
ℏ
2
k
+
1
{\displaystyle t={\cfrac {1}{-\ {\cfrac {im\lambda }{\hbar ^{2}k}}+1}}\,\!}
。
反射系数是
R
=
|
r
|
2
=
1
1
+
ℏ
4
k
2
m
2
λ
2
=
1
1
+
2
ℏ
2
E
m
λ
2
{\displaystyle R=|r|^{2}={\cfrac {1}{1+{\cfrac {\hbar ^{4}k^{2}}{m^{2}\lambda ^{2}}}}}={\cfrac {1}{1+{\cfrac {2\hbar ^{2}E}{m\lambda ^{2}}}}}\,\!}
。
这纯粹是一个量子力学的效应;在经典力学里,这是不可能发生的。
透射系数是
T
=
|
t
|
2
=
1
−
R
=
1
1
+
m
2
λ
2
ℏ
4
k
2
=
1
1
+
m
λ
2
2
ℏ
2
E
{\displaystyle T=|t|^{2}=1-R={\cfrac {1}{1+{\cfrac {m^{2}\lambda ^{2}}{\hbar ^{4}k^{2}}}}}={\cfrac {1}{1+{\cfrac {m\lambda ^{2}}{2\hbar ^{2}E}}}}\,\!}
。
由于模型的对称性,假若,粒子从右边入射,我们也会得到同样的答案。
很奇异地,给予同样的能量、质量、与狄拉克Delta函数的强度,Delta位势垒与Delta位势阱有同样的反射系数与透射系数。
Delta位势阱的束缚态,在任何一个位置,波函数都是连续的;可是,除了在
x
=
0
{\displaystyle x=0\,\!}
以外,在其它任何位置,波函数随位置的导数都是连续的。
每一个一维的吸引位势,都至少会存在着一个束缚态 (bound state )。由于
E
<
0
{\displaystyle E<0\,\!}
,波数变为复数。设定
κ
=
−
i
k
=
2
m
|
E
|
/
ℏ
2
{\displaystyle \kappa =-ik={\sqrt {2m|E|/\hbar ^{2}}}\,\!}
。前述的振荡的波函数
ψ
L
{\displaystyle \psi _{L}\,\!}
与
ψ
R
{\displaystyle \psi _{R}\,\!}
,现在却随着坐标
x
{\displaystyle x\,\!}
呈指数递减或指数递增。为了要符合物理的真实性,我们要求波函数不发散 于
x
→
±
∞
{\displaystyle x\to \pm \infty \,\!}
。那么,
A
r
{\displaystyle A_{r}\,\!}
与
B
l
{\displaystyle B_{l}\,\!}
必须被设定为0。波函数变为
ψ
L
(
x
)
=
A
l
e
κ
x
{\displaystyle \psi _{L}(x)=A_{l}e^{\kappa x}\,\!}
,
ψ
R
(
x
)
=
B
r
e
−
κ
x
{\displaystyle \psi _{R}(x)=B_{r}e^{-\kappa x}\,\!}
。
从边界条件与归一条件 ,可以得到
A
l
=
B
r
=
κ
{\displaystyle A_{l}=B_{r}={\sqrt {\kappa }}\,\!}
,
κ
=
m
λ
ℏ
2
{\displaystyle \kappa ={\frac {m\lambda }{\hbar ^{2}}}\,\!}
。
Delta位势阱只能有一个束缚态。束缚态的能量是
E
=
−
ℏ
2
κ
2
2
m
=
−
m
λ
2
2
ℏ
2
{\displaystyle E=-\ {\frac {\hbar ^{2}\kappa ^{2}}{2m}}=-\ {\frac {m\lambda ^{2}}{2\hbar ^{2}}}\,\!}
。
束缚态的波函数是
ψ
(
x
)
=
m
λ
ℏ
e
−
m
λ
∣
x
∣
/
ℏ
2
{\displaystyle \psi (x)={\frac {\sqrt {m\lambda }}{\hbar }}e^{-m\lambda \mid x\mid /\hbar ^{2}}\,\!}
。
Delta位势阱是有限深方形阱 的一个特别案例。在有限深位势阱的深度
V
0
→
∞
{\displaystyle V_{0}\to \infty \,\!}
与阱宽
L
→
0
{\displaystyle L\to 0\,\!}
的极限,同时保持
V
0
L
=
λ
{\displaystyle V_{0}L=\lambda \,\!}
,就可以从有限深位势阱的波函数,得到Delta位势阱的波函数。
当核间距R=2时,双势井狄拉克Delta函数模型中的对称与反对称的波函数
Delta函数模型其实是氢原子 的一维版本根据维度比例由 达德利·赫施巴赫 (“Dudley R. Herschbach”)[ 1] 团队所研发。此 delta函数模型以双井迪拉克Delta函数模型最有用,因其代表一维版的水分子离子。
双井迪拉克Delta函数模型是用以下薛定谔方程描述:
−
ℏ
2
2
m
d
2
ψ
d
x
2
(
x
)
+
V
(
x
)
ψ
(
x
)
=
E
ψ
(
x
)
{\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}\psi }{dx^{2}}}(x)+V(x)\psi (x)=E\psi (x)}
电势现为:
V
(
x
)
=
−
q
[
δ
(
x
−
R
2
)
+
λ
δ
(
x
+
R
2
)
]
{\displaystyle V(x)=-q\left[\delta (x-{\frac {R}{2}})+\lambda \delta (x+{\frac {R}{2}})\right]}
其中
0
<
R
<
∞
{\displaystyle 0<R<\infty }
是“核间”距离于迪拉克Delta函数(负)峰值位于
x
=
±
R
2
{\displaystyle x=\pm {\textstyle {\frac {R}{2}}}}
(图表中棕色所示)。记得此模型与其三维分子版本的关系,我们用原子单位制 且设
ℏ
=
m
=
1
{\displaystyle \hbar =m=1}
。此处
0
<
λ
<
1
{\displaystyle 0<\lambda <1}
为一可调参数。从单井的例子,可推论拟设 于此解为:
ψ
(
x
)
=
A
e
−
d
|
x
−
R
2
|
+
B
e
−
d
|
x
+
R
2
|
{\displaystyle \psi (x)~=~Ae^{-d\left|x-{\frac {R}{2}}\right|}+Be^{-d\left|x+{\frac {R}{2}}\right|}}
令波函数于Delta函数峰值相等可得行列式 :
|
q
−
d
q
e
−
d
R
q
λ
e
−
d
R
q
λ
−
d
|
=
0
,
E
=
−
d
2
2
.
{\displaystyle \left|{\begin{array}{cc}q-d&qe^{-dR}\\q\lambda e^{-dR}&q\lambda -d\end{array}}\right|=0~,\qquad E=-{\frac {d^{2}}{2}}~.}
因此,
d
{\displaystyle d}
是由伪二次式方程:
d
±
(
λ
)
=
1
2
q
(
λ
+
1
)
±
1
2
{
q
2
(
1
+
λ
)
2
−
4
λ
q
2
[
1
−
e
−
2
d
±
(
λ
)
R
]
}
1
/
2
{\displaystyle d_{\pm }(\lambda )~=~{\textstyle {\frac {1}{2}}}q(\lambda +1)\pm {\textstyle {\frac {1}{2}}}\left\{q^{2}(1+\lambda )^{2}-4\,\lambda q^{2}\lbrack 1-e^{-2d_{\pm }(\lambda )R}]\right\}^{1/2}}
它有两解
d
=
d
±
{\displaystyle d=d_{\pm }}
。若等价情况(对称单核),
λ
=
1
{\displaystyle \lambda =1}
则伪二次式化为:
d
±
=
q
[
1
±
e
−
d
±
R
]
{\displaystyle d_{\pm }=q[1\pm e^{-d_{\pm }R}]}
此“+”代表了对称于中点的波函数(图中红色)而
A
=
B
{\displaystyle A=B}
称为偶态。接着,“-”情况为反对称于中点的波函数其
A
=
−
B
{\displaystyle A=-B}
称为非偶态(图中绿色)。它们代表着三维
H
2
+
{\displaystyle H_{2}^{+}}
的两种最低能态之近似且有助于其分析。对称电价的特征能分析解为[ 2] :
d
±
=
q
+
W
(
±
q
R
e
−
q
R
)
/
R
{\displaystyle d_{\pm }=q~+~W(\pm qRe^{-qR})/R}
其中W是标准朗伯W函数 注意此最低能对应于对称解
d
+
{\displaystyle d_{+}}
。当非等电价,此为三维分子问题,其解为一般化Lambert W函数(见一般化朗伯W函数 章节与相关参考)。
^ D.R Herschbach , J.S. Avery, and O. Goscinski (eds.), Dimensional Scaling in Chemical Physics , Springer, (1992). [1] (页面存档备份 ,存于互联网档案馆 )
^ T.C. Scott, J.F. Babb, Alexander Dalgarno and John D. Morgan III, "The Calculation of Exchange Forces: General Results and Specific Models", J. Chem. Phys. , 99, pp. 2841-2854, (1993). [2]