方势垒

                     

贡献者: addis

预备知识 无限深势阱,一维散射态的正交归一化

   本文采用原子单位制。我们要解一维定态薛定谔方程(参考式 9 )为

(1)12md2ψdx2+Vψ=Eψ .
令方势垒长度为 2l,关于原点对称。势能函数为
(2)V(x)={0(x<l)V0(lx<l)0(lx) .
一些教材把势垒放在区间 [0,a],这时只需令 a=2l,把解出的波函数向右平移 l 即可。以下讨论两种常见的解,一种类似于 sin(kx),cos(kx) 另一种类似于 exp(ikx),exp(ikx)

1. 奇偶函数解

   对称势能的好处是存在奇和偶的实函数解,且它们自动正交。

EV0 的情况

   令

(3)k=2mE ,b=2m(EV0) .
令对称解和反对称的波函数分别为
(4)ψke(x)={A1cos(bx)(0xl)C1cos(kx)+D1sin(kx)(l<x)ψke(x)(x<0) ,
(5)ψko(x)={B2sin(bx)(0xl)C2cos(kx)+D2sin(kx)(l<x)ψko(x)(x<0) ,
其中 l<x 的部分也可以分别记为
(6)Gisin(kx+ϕi)(i=1,2) ,
(7)Gi=Ci2+Di2ϕi=Arctan(Ci,Di)(i=1,2) .

   在 x=l 处匹配波函数和一阶导数,解得

(8){C1A1=cos(kl)cos(bl)+bksin(kl)sin(bl)D1A1=bkcos(kl)sin(bl)+sin(kl)cos(bl) ,
(9){C2B2=cos(kl)sin(bl)bksin(kl)cos(bl)D2B2=bkcos(kl)cos(bl)+sin(kl)sin(bl) .
对波函数归一化使 ψk(x)ψk(x)=δ(kk),令无穷远处振幅为 1/π定理 1 ),得
(10)A1=1πcos2(bl)+(b2/k2)sin2(bl) ,B2=1πsin2(bl)+(b2/k2)cos2(bl) .

0<E<V0 的情况

   令 κ=2m(V0E)

(11)ψke(x)={A1cosh(κx)(0xl)C1cos(kx)+D1sin(kx)(l<x) ,
(12)ψko(x)={A2sinh(κx)(0xl)C2cos(kx)+D2sin(kx)(l<x) .
x=l 处匹配波函数和一阶导数,分别得
(13){A1cosh(κl)=C1cos(kl)+D1sin(kl)κA1sinh(κl)=kC1sin(kl)+kD1cos(kl) ,
(14){A2sinh(κl)=C2cos(kl)+D2sin(kl)κA2cosh(κl)=kC2sin(kl)+kD2cos(kl) .
解得
(15){C1A1=cosh(κl)cos(kl)κksinh(κl)sin(kl)D1A1=cosh(κl)sin(kl)+κksinh(κl)cos(kl) ,
(16){C2A2=sinh(κl)cos(kl)κkcosh(κl)sin(kl)D2A2=sinh(κl)sin(kl)+κkcosh(κl)cos(kl) .
对波函数归一化,同样令无穷远处振幅为 1/π
(17)A1=1πcosh2(κl)+(κ/k)2sinh2(kl) ,
(18)A2=1πsinh2(κl)+(κ/k)2cosh2(kl) .

2. 行波解

   令左、中、右三个区间 exp 项的系数分别为 A,BC,DE,F。左边右入射和出射,右边只有出射没有反射

(19)ψk,1(x)={Aexp(ikx)+Bexp(ikx)(x<L)Cexp(ibx)+Dexp(ibx)(LxL)Fexp(ikx)(x<L) ,
(20)ψk,2(x)=ψk,1(x) ,
要归一化只需满足 A=1/2π式 19 )。

E>V0 的情况

   当 E>V0 时,系数解为

(21){A=1/2πB=N(k2b2)sin(2bl)exp(ikl)C=Nik(k+b)exp(ibl)D=Nik(bk)exp(ibl)E=N2ikbexp(ikl) ,
其中
(22)N=12π(k2+b2)sin(2bl)2ikbcos(2bl)(k2b2)2sin2(2bl)+4k2b2exp(ikl) .
可以验证 ψk,1,ψk,2 正交的条件(式 25 Re[BE]=0

   透射率,反射率分别为

(23)T=|EA|2=4k2b2(k2b2)2sin2(2bl)+4k2b2 ,
(24)R=|BA|2=(k2b2)sin2(2bl)(k2b2)2sin2(2bl)+4k2b2 .

0<E<V0 的情况

   当 0<E<V0 时,系数解为

(25){A=1/2πB=N(κ2+k2)sinh(2κl)exp(ikl)C=Nk(iκk)exp(κl)D=Nk(iκ+k)exp(κl)E=N2ikκexp(ikl) ,
(26)N=12π(k2κ2)sinh(2κl)2ikκcosh(2κl)(k2+κ2)2sinh2(2κl)+4k2κ2exp(ikl) ,
(27)T=|EA|2=4k2κ2(k2+κ2)2sinh2(2κl)+4k2κ2 ,
(28)R=|BA|2=(κ2+k2)sinh2(2κl)(k2+κ2)2sinh2(2κl)+4k2κ2 .

  

未完成:未完成:透射率反射率画图

投影系数

   任意平方可积的波函数可以用正交归一化的行波展开为

(29)ψ(x)=0+C1(k)ψk,1(x)+C2(k)ψk,2(x)dk ,
ψ(x) 只在 x<l 不为零,则有
(30)C1(k)=ψk,1|ψ=ψ~(k)+2πBψ~(k) ,C2(k)=ψk,2|ψ=2πEψ~(k) ,
其中 ψ~(k)ψ(x)傅里叶变换
(31)ψ~(k)=12π+ψ(x)eikxdx .
具体例子见 “方势垒散射数值计算”,我们将把一个高斯波包展开为散射态的积分。


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利