量子力学(七):其他常见的一维势

到目前为止我们已经介绍了两种典型的一维势:一维无限深方势阱以及一维谐振子。接下来我们介绍其他几种同样典型的一维势。

一维delta势阱

一维delta势,顾名思义,其势函数和delta函数有关,也即:

δ(x)={0,  x0,  x=0\delta(x)=\left\{\begin{aligned}&0,\;x\neq 0\\ &\infty ,\; x=0 \end{aligned} \right. δ(x)dx=1\int_{-\infty}^\infty \delta(x)dx=1

我们考虑这样一个delta势阱:

V(x)=aδ(x),  a>0V(x)=-a\delta(x),\;a>0

x<0x<0 时,我们有一维定态薛定谔方程:

2md2ψdx2=Eψ-{\hbar\over 2m}{d^2\psi\over dx^2}=E\psi

我们考虑该势阱下的一个束缚态解。束缚态的定义是:当x|x|\to \infty 时,ψ(x)2dx0|\psi(x)|^2dx\to 0 ,即在无穷远处找到粒子的概率为00 ,而这就要求E<V()E<V(\infty) 。在我们当下讨论的delta势中,则有

E<V()=0E<V(\infty)=0

那么我们令k=2mEk={\sqrt{-2mE}\over\hbar} ,得到关于xx 的微分方程:

k2ψ2=d2ψdx2k^2\psi^2={d^2\psi\over dx^2}

得到在x<0x<0 时,ψ(x)\psi(x) 的形式

ψ(x)=Aekx+Bekx\psi(x)=Ae^{kx}+Be^{-kx}

考虑到x,  ψ(x)0x\to -\infty,\;\psi(x)\to 0 ,就应该有

ψ(x)=Aekx\psi(x)=Ae^{kx}

类似地,在x>0x>0 时,我们就得到

ψ(x)=Cekx\psi(x)=Ce^{-kx}

根据波函数的连续性条件,我们得到A=CA=C 。接下来我们就考虑处理delta势,即x=0x=0 时的情况。设一个很小很小的正数ε>0\varepsilon>0 ,我们对一维定态薛定谔方程两边积分:

εε22md2ψdx2aδ(x)ψ(x)dx=εεEψ(x)dx(1)\int_{-\varepsilon}^\varepsilon -{\hbar^2\over 2m} {d^2\psi\over dx^2}-a\delta(x)\psi(x)dx=\int_{-\varepsilon}^\varepsilon E\psi(x)dx\tag{1}

其中 εεd2ψdx2dx=ψ(ε)ψ(ε)=ψ(0+)ψ(0)\int_{-\varepsilon}^\varepsilon{d^2\psi\over dx^2}dx=\psi'(\varepsilon)-\psi'(-\varepsilon)=\psi'(0^+)-\psi(0^-)

ψ(0+)=(Aekx)x=0=kA\psi'(0^+)=(Ae^{-kx})'|_{x=0}=-kA ψ(0)=(Aekx)x=0=kA\psi(0^-)=(Ae^{kx})' |_{x=0}=kA

εεd2ψdx2dx=ψ(0+)ψ(0)=2kA\int_{-\varepsilon}^\varepsilon{d^2\psi\over dx^2}dx=\psi'(0^+)-\psi(0^-)=-2kA εεaδ(x)ψ(x)dx=aψ(0)=aA\int_{-\varepsilon}^\varepsilon a\delta(x)\psi(x)dx=a\psi(0)=aA εεEψ(x)dx=0(2)\tag{2}\int_{-\varepsilon}^\varepsilon E\psi(x)dx=0

其中式(2)(2) 是考虑到EEψ(x)\psi(x) 均为有限量,当ε0\varepsilon \to 0 时,自然有式(2)(2) 成立。

全部带入式(1)(1) 我们得到

k2mAaA=0{k\hbar^2\over m}A-aA=0

得到k=ma2k={ma\over \hbar^2}

k=2mEk={\sqrt{-2mE}\over \hbar}

解得束缚态能量

E=m2a222E=-{m^2a^2\over 2\hbar^2}

归一化:

ψ(x)=20A2e2kxdx=2A2(12k)(e2kxxe2kxx=0)=A2k=1\begin{align*}\int_{-\infty}^\infty\psi(x)&=2\int_0^\infty |A|^2e^{-2kx}dx\\&=2|A|^2(-{1\over 2k})(e^{-2kx}|_{x\to\infty}-e^{-2kx}|_{x=0})\\&={|A|^2\over k}\\&=1\end{align*}

得到归一化常数 A=k=maA=\sqrt{k}={\sqrt{ma}\over\hbar}

由此可以得到一维delta势阱的束缚态解:

ψ(x)={maekx,  x>0maekx,  x<0\psi(x)=\left\{\begin{aligned}{\sqrt{ma}\over\hbar}e^{-kx},\;x>0\\ {\sqrt{ma}\over\hbar}e^{kx},\;x<0\end{aligned}\right.

其中k=ma2k={ma\over\hbar^2}

一维自由粒子

顾名思义,自由粒子就是不受任何势束缚的粒子,也即在全空间上均有:V(x)=0V(x)=0 在这种情况下,粒子的束缚态解是不存在的,因此,我们采用类似于平面波的方式来描述它。我们还是从一维定态薛定谔方程出发:

2md2ψdx2=Eψ-{\hbar\over 2m}{d^2\psi\over dx^2}=E\psi

在这里我们要讨论的情况和一维delta势有些许不同:在一维delta势中,我们通过考虑束缚态解,得到了E<0E<0 的限制条件;但在一维自由势中,如果E<0E<0 ,那么在全空间中均有E<VE<V ,而这样的粒子是不能在空间中存在的;而在一维delta势中存在特殊的点,也即x=0x=0 处,使得粒子的能量不至于在全空间中都小于粒子的势能;尽管在x=0x=0 之外的所有地方仍然是E<VE<V ,也就是粒子的总能量小于其势能的情况出现,但是通过量子力学效应,我们让这个粒子得以存在。量子力学,很神奇吧。

换言之,此时我们应该考虑的是E>0E>0 的情况;并且我们一开始就已经说过,在自由势中粒子的束缚态是不存在的,因为根本就没有势来束缚它,所以我们设k=2mEk={\sqrt{2mE}\over \hbar} ,就得到

d2ψdx2=k2ψ{d^2\psi\over dx^2}=-k^2\psi

解得

ψ(x)=Aeikx+Beikx\psi(x)=Ae^{ikx}+Be^{-ikx}

如果我们引入时间演化项φ(t)=eiEt\varphi(t)=e^{-i{E\over\hbar}t} ,我们就可以得到:

Ψ(x,t)=ψ(x)φ(t)=Aei(Et+kx)+Bei(Etkx)\Psi(x,t)=\psi(x)\varphi(t)=Ae^{i(-{E\over\hbar}t+kx)}+Be^{i(-{E\over\hbar}t-kx)}

应该注意到这是非常类似于平面简谐波的表达式形式。一般地,我们认为eikxe^{ikx} 项是从左到右传播的波,相应地eikxe^{-ikx} 项表示从右往左传播的波。

一维势垒穿透

我们现在考虑这样一种情况:粒子从无穷远处,自左向右地传播;而在x=0x=0 处存在一个阶跃势垒:

V(x)={0,x<0V0,x0V(x)=\left\{\begin{aligned}0,x<0\\V_0,x\geq 0\end{aligned}\right.

这个粒子波撞上势垒后,一部分反射回去,另一部分穿过势垒继续传播。方便起见,我们设入射波的强度为11 ,也即

ψ1(x)=eikx,  x<0\psi_1(x)=e^{ikx},\;x<0

其中k=2mEk={\sqrt{-2mE}\over \hbar}

把这个势和波函数画出来,大概会长这样:

这里画出的图例是波的能量E>V0E>V_0 时的情况,我们暂且先讨论这种情况。

我们设反射波的波函数:

ψ2(x)=Reikx,  x<0\psi_2(x)=Re^{-ikx},\;x<0

其中RR 为未定常数。

接下来我们考虑透射波的波函数形式。注意到在x<0x<0 处的波函数,我们可以直接套用一维自由粒子的解的形式ψ(x)=Aeikx+Beikx\psi(x)=Ae^{ikx}+Be^{-ikx} ,并且由规定的传播方向直接得到解的形式。但在x>0x>0 处,由于阶跃势的存在,我们不能直接套用,而是要重新求解薛定谔方程。

我们首先写出x>0x>0 处,透射波满足的薛定谔方程:

22md2ψ3dx2+V0ψ3(x)=Eψ3(x)-{\hbar^2\over 2m}{d^2\psi_3\over dx^2}+V_0\psi_3(x)=E\psi_3(x)

注意到E>V0E>V_0 ,我们不妨就设κ=2m(EV0)\kappa={\sqrt{-2m(E-V_0)}\over \hbar}

就得到ψ3=Seiκx,  x>0\psi_3=Se^{i\kappa x},\; x>0 。其中SS 为未定常数。

为计算未定常数RRSS ,我们需要考虑波函数的连续性条件(怎么又是你)。在x=0x=0 处,由

ψ(0+)=ψ(0)\psi(0^+)=\psi(0^-) ψ(0+)=ψ(0)\psi'(0^+)=\psi'(0^-)

得到

1+R=S1+R=S ikikR=iκSik-ikR=i\kappa S

可以解得

R=kκk+κR={k-\kappa\over k+\kappa} S=2kk+κS={2k\over k+\kappa}

一般称RR反射系数,记反射率R2=(kκk+κ)2|R|^2=({k-\kappa\over k+\kappa})^2透射率1R2=4kκ(k+κ)21-|R|^2={4k\kappa\over(k+\kappa)^2} 。非常需要注意的是S21R2|S|^2\neq 1-|R|^2 ,也即不能记透射率为S2|S|^2 ,因为这里的反射率和透射率是以入射波为基准的。因此我们解得这个势的波函数为:

ψ(x)={eikx+Reikx,  x<0Seiκx,  x>0\psi(x)=\left\{ \begin{align*}e^{ikx}+Re^{-ikx},\;x<0\\ Se^{i\kappa x},\;x>0 \end{align*} \right.

接下来我们再讨论E<V0E<V_0 的情况。同样的,在x<0x<0 时,我们仍然可以得到

ψ(x)=eikx+Reikx,  x<0\psi(x)=e^{ikx}+Re^{-ikx},\;x<0

但在x>0x>0 处,由于势的改变,我们得到解的形式有所不同:

ψ(x)=Seκx\psi(x)=Se^{-\kappa x}

其中κ=2m(EV0)\kappa={\sqrt{-2m(E-V_0)}\over \hbar}这是一个实的指数函数,而非一个平面波的形式,可以预见,在x>0x>0 时,随着xx 的增大,在xdxxdx 附近找到粒子的概率将会迅速衰减,并且在xx\to \infty 时,找到粒子的概率为00

继续带入波函数的连续性条件,在x=0x=0 处有:

1+R=S1+R=S ikikR=Sκik-ikR=-S\kappa

解得

R=ik+κikκR={ik+\kappa\over ik-\kappa} S=2ikikκS={2ik\over ik-\kappa}

同样的,反射率为R2|R|^2 ,透射率为1R21-|R|^2

  • Copyright: Copyright is owned by the author. For commercial reprints, please contact the author for authorization. For non-commercial reprints, please indicate the source.

扫一扫,分享到微信

微信分享二维码

请我喝杯咖啡吧~

支付宝
微信