AI 正在绞尽脑汁想思路 ING···
AI 摘要
DeepSeek & Kimi

§ 一维谐振子

图1-1 弹簧振子

力:F(x)=−kxF(x)=-kx
势函数:V(x)=12kx2=12kmmx2=12ω02mx2V(x)=\frac{1}{2}kx^2=\frac{1}{2}\frac{k}{m}mx^2=\frac{1}{2}\omega_{0}^{2}mx^2
弹簧振子的固有角频率:ω02=km\omega_{0}^{2}=\frac{k}{m}\,

Schrödinger 解法:

一维谐振子势函数:V(x)=12ω02mx2V(x)=\frac{1}{2}\omega_{0}^{2}mx^2
相应 Schrödinger 方程:

iℏ∂∂tψ(x,t)=−ℏ22m∂2∂x2ψ(x,t)+12ω02mx2ψ(x,t)i\hbar\frac{\partial}{\partial t}\psi(x,t)=-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}\psi(x,t)+\frac{1}{2}\omega_{0}^{2}mx^2\psi(x,t)

使用分离变量法,把 tt 分离出来

令 ψ(x,t)=ϕ(x)e−iωt\psi\left(x,t\right)=\phi\left(x\right)e^{-i\omega t} ,ω\omega 为待定系数
Schrödinger 方程改写为

iℏ(−iω)ϕ(x)e−iωt=[−ℏ22md2ϕ(x)dx2+12mω02x2ϕ(x)]e−iωti\hbar\left(-i\omega\right)\phi\left(x\right)e^{-i\omega t}=\left[-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2\phi\left(x\right)}{\mathrm{d}x^2}+\frac{1}{2}m\omega_{0}^{2}x^2\phi\left(x\right)\right]e^{-i\omega t}

ℏωϕ(x)=Eϕ(x)=−ℏ22md2ϕ(x)dx2+12mω02x2ϕ(x)\hbar\omega\phi\left(x\right)=E\phi\left(x\right)=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2\phi\left(x\right)}{\mathrm{d}x^2}+\frac{1}{2}m\omega_{0}^{2}x^2\phi\left(x\right)

化为解微分方程的定态问题

d2ϕ(x)dx2+2mℏ2(E−12mω02x2)ϕ(x)=0\frac{\mathrm{d}^2\phi\left(x\right)}{\mathrm{d}x^2}+\frac{2m}{\hbar^2}\left(E-\frac{1}{2}m\omega_{0}^{2}x^2\right)\phi\left(x\right)=0

奇异点分析:

当 ∣x∣→∞\left|x\right|\rightarrow\infty 时,V(x)→∞V\left(x\right)\rightarrow\infty ,粒子为束缚态。lim⁡∣x∣→∞φ(x)=0\underset{\left|x\right|\rightarrow\infty}{\lim}\varphi\left(x\right)=0

−ℏ22md2ϕ(x)dx2+12mω02x2ϕ(x)≅0-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2\phi\left(x\right)}{\mathrm{d}x^2}+\frac{1}{2}m\omega_{0}^{2}x^2\phi\left(x\right)\cong 0

令方程解为 ϕ~=eλx\tilde{\phi}=e^{\lambda x} ,则 ϕ′=λeλx\phi^{\prime}=\lambda e^{\lambda x} ,ϕ′′=λ2eλx\phi^{\prime\prime}=\lambda^2e^{\lambda x}\,

−ℏ22mλ2eλx+12mω02x2eλx=0-\frac{\hbar^2}{2m}\lambda^2\cancel{e^{\lambda x}}+\frac{1}{2}m\omega_{0}^{2}x^2\cancel{e^{\lambda x}}=0

则

λ=±mω0ℏx\lambda=\pm\frac{m\omega_0}{\hbar}x

ϕ~=exp⁡(±mω02ℏx2)\tilde{\phi}=\exp\left(\pm\frac{m\omega_0}{2\hbar}x^2\right)

dϕ~(x)dx=exp⁡(−mω02ℏx2)(−mω0ℏx)\frac{d\tilde{\phi}\left(x\right)}{dx}=\exp\left(-\frac{m\omega_0}{2\hbar}x^2\right)\left(-\frac{m\omega_0}{\hbar}x\right)

d2ϕ~(x)dx2=exp⁡(−mω02ℏx2)(−mω0ℏx)2−mω0ℏexp⁡(−mω02ℏx2)\frac{d^2\tilde{\phi}\left(x\right)}{dx^2}=\exp\left(-\frac{m\omega_0}{2\hbar}x^2\right)\left(-\frac{m\omega_0}{\hbar}x\right)^2-\frac{m\omega_0}{\hbar}\exp\left(-\frac{m\omega_0}{2\hbar}x^2\right)

代入渐近方程:

−ℏ22md2ϕ~(x)dx2+12mω02x2ϕ~(x)=[−ℏ22m(mω0ℏx)2−mω0ℏ(−ℏ22m)]exp⁡(−mω02ℏx2)+12mω02x2exp⁡(−mω02ℏx2)≅(−12mω02x2+12mω02x2)exp⁡(−mω02ℏx2)=0\begin{aligned} -\frac{\hbar^2}{2m}\frac{d^2\tilde{\phi}\left(x\right)}{dx^2}+\frac{1}{2}m\omega_{0}^{2}x^2\tilde{\phi}\left(x\right)&= \left[-\frac{\hbar^2}{2m}\left(\frac{m\omega_0}{\hbar}x\right)^2-\frac{m\omega_0}{\hbar}\left(-\frac{\hbar^2}{2m}\right)\right]\exp\left(-\frac{m\omega_0}{2\hbar}x^2\right)+\frac{1}{2}m\omega_{0}^{2}x^2\exp\left(-\frac{m\omega_0}{2\hbar}x^2\right)\\ &\cong\left(-\frac{1}{2}m\omega_{0}^{2}x^2+\frac{1}{2}m\omega_{0}^{2}x^2\right)\exp\left(-\frac{m\omega_0}{2\hbar}x^2\right)=0 \end{aligned}

推得

−ℏ22md2ϕ~(x)dx2+12mω02x2ϕ~(x)≅0-\frac{\hbar^2}{2m}\frac{d^2\tilde{\phi}\left(x\right)}{dx^2}+\frac{1}{2}m\omega_{0}^{2}x^2\tilde{\phi}\left(x\right)\cong 0

解得

ϕ~1=exp⁡(−mω02ℏx2),ϕ~2=exp⁡(mω02ℏx2)(x→∞,是发散的,舍弃)\tilde{\phi}_1=\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\text{,}\tilde{\phi}_2=\exp\left(\frac{m\omega_{0}}{2\hbar}x^2\right)\text{(}x\rightarrow\infty\text{,是发散的,舍弃)}

我们要求波函数是平方可积的,所以要求 x→∞x\rightarrow\infty ,波函数一定趋近于 00
令

ϕ(x)=ϕ~(x)χ(x)=exp⁡(−mω02ℏx2)χ(x)\phi\left(x\right)=\tilde{\phi}\left(x\right)\chi\left(x\right)=\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi\left(x\right)

ϕ′(x)=−mω0ℏxexp⁡(−mω02ℏx2)χ(x)+exp⁡(−mω02ℏx2)χ′(x)\phi^{\prime}\left(x\right)=-\frac{m\omega_{0}}{\hbar}x\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi\left(x\right)+\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi^{\prime}\left(x\right)

ϕ′′(x)=[(mω0ℏx)2−mω0ℏ]exp⁡(−mω02ℏx2)χ(x)−2mω0ℏxexp⁡(−mω02ℏx2)χ′(x)+exp⁡(−mω02ℏx2)χ′′(x)\phi^{\prime\prime}\left(x\right)=\left[\left(\frac{m\omega_{0}}{\hbar}x\right)^2-\frac{m\omega_{0}}{\hbar}\right]\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi\left(x\right)-\frac{2m\omega_{0}}{\hbar}x\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi^{\prime}\left(x\right)+\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi^{\prime\prime}\left(x\right)

代入

−ℏ22mϕ′′(x)+12mω02x2ϕ(x)=Eϕ(x)-\frac{\hbar^2}{2m}\phi^{\prime\prime}\left(x\right)+\frac{1}{2}m\omega_{0}^{2}x^2\phi\left(x\right)=E\phi\left(x\right)

−ℏ22m[(mω0ℏx)2−mω0ℏ]exp⁡(−mω02ℏx2)χ(x)+12mω02x2exp⁡(−mω02ℏx2)χ(x)−2mω0ℏxexp⁡(−mω02ℏx2)χ′(x)=Eexp⁡(−mω02ℏx2)χ(x)+exp⁡(−mω02ℏx2)χ′′(x)-\frac{\hbar^2}{2m}\left[\left(\frac{m\omega_{0}}{\hbar}x\right)^2-\frac{m\omega_{0}}{\hbar}\right]\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi\left(x\right)+\frac{1}{2}m\omega_{0}^{2}x^2\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi\left(x\right)-\frac{2m\omega_{0}}{\hbar}x\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi^{\prime}\left(x\right)=E\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi\left(x\right)+\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi^{\prime\prime}\left(x\right)

−ℏ22mχ′′(x)+ℏω0xχ′(x)+12ℏω0χ(x)−12mω02x2χ(x)+12mω02x2χ(x)=Eχ(x)-\frac{\hbar^2}{2m}\chi^{\prime\prime}\left(x\right)+\hbar\omega_{0}x\chi^{\prime}\left(x\right)+\frac{1}{2}\hbar\omega_{0}\chi\left(x\right)-\frac{1}{2}m\omega_{0}^{2}x^2\chi\left(x\right)+\frac{1}{2}m\omega_{0}^{2}x^2\chi\left(x\right)=E\chi\left(x\right)

ℏ22mχ′′(x)−ℏω0xχ′(x)+(E−12ℏω)χ(x)=0(福克斯Fuchs微分方程)\frac{\hbar^2}{2m}\chi^{\prime\prime}\left(x\right)-\hbar\omega_{0}x\chi^{\prime}\left(x\right)+\left(E-\frac{1}{2}\hbar\omega\right)\chi\left(x\right)=0\text{(福克斯}Fuchs\text{微分方程)}

奇异点分析做完,把奇异点扣除,由于方程是二阶的,所以一定有两个线性无关的解在零点处解析(可在零点展成无穷级数)

χ(x)=C0+C1x+C2x2+⋯χ′(x)=C1+2C2x+3C3x2+⋯χ′′(x)=2C2+6C3x+⋯\begin{aligned} \chi\left(x\right)&=\mathrm{C}_0+\mathrm{C}_1x+\mathrm{C}_2x^2+\cdots\\ \chi^{\prime}\left(x\right)&=\mathrm{C}_1+2\mathrm{C}_2x+3\mathrm{C}_3x^2+\cdots\\ \chi^{\prime\prime}\left(x\right)&=2\mathrm{C}_2+6\mathrm{C}_3x+\cdots \end{aligned}

ℏ22mχ′′−ℏω0xχ′+(E−12ℏω)χ=0\frac{\hbar^2}{2m}\chi^{\prime\prime}-\hbar\omega_{0}x\chi^{\prime}+\left(E-\frac{1}{2}\hbar\omega\right)\chi=0

ℏ22m(2C2+6C3x+⋯)−ℏω0x(C1+2C2x+3C3x2+⋯)+(E−12ℏω)(C0+C1x+C2x2+⋯)=0\frac{\hbar^2}{2m}\left(2\mathrm{C}_2+6\mathrm{C}_3x+\cdots\right)-\hbar\omega_0x\left(\mathrm{C}_1+2\mathrm{C}_2x+3\mathrm{C}_3x^2+\cdots\right)+\left(E-\frac{1}{2}\hbar\omega\right)\left(\mathrm{C}_0+\mathrm{C}_1x+\mathrm{C}_2x^2+\cdots\right)=0

展开每项

ℏ22m⋅2C2+ℏ22m⋅6C3x+⋯−ℏω0xC1−ℏω0x⋅2C2x+⋯+(E−12ℏω)C0+(E−12ℏω)C1x+(E−12ℏω)C2x2+⋯=0\begin{aligned} \quad\qquad\frac{\hbar^2}{2m}\cdot 2\mathrm{C}_2&+\frac{\hbar^2}{2m}\cdot 6\mathrm{C}_3x+\cdots\\ &-\hbar\omega _0x\mathrm{C}_1-\hbar\omega_{0}x\cdot 2\mathrm{C}_2x+\cdots\\ +\left(E-\frac{1}{2}\hbar\omega\right)\mathrm{C}_0&+\left(E-\frac{1}{2}\hbar\omega\right)\mathrm{C}_1x+\left(E-\frac{1}{2}\hbar\omega\right)\mathrm{C}_2x^2+\cdots=0 \end{aligned}

xx 的一次幂、二次幂、三次幂作为函数是线性无关的,若其线性关系等于 00 ,则所乘系数必须等于 00 。

x0=ℏ22m⋅2C2+(E−12ℏω)C0=0x1=ℏ22m⋅6C3−ℏω0xC1+(E−12ℏω)C1=0x2=−ℏω0x⋅2C2+(E−12ℏω)C2=0⋮\begin{aligned} x^0&=\frac{\hbar^2}{2m}\cdot 2\mathrm{C}_2+\left(E-\frac{1}{2}\hbar\omega\right)\mathrm{C}_0=0\\ x^1&=\frac{\hbar^2}{2m}\cdot 6\mathrm{C}_3-\hbar\omega_{0}x\mathrm{C}_1+\left(E-\frac{1}{2}\hbar\omega\right)\mathrm{C}_1=0\\ x^2&=-\hbar\omega_{0}x\cdot 2\mathrm{C}_2+\left(E-\frac{1}{2}\hbar\omega\right)\mathrm{C}_2=0\\ \end{aligned} \vdots

C0\mathrm{C}_0 和 C1\mathrm{C}_1 不确定(有一定自由度),一旦给定,其余系数也都能给定。

  1. 取 C0=1\mathrm{C}_0=1 ,C1=0\mathrm{C}_1=0 ,则 C1=C3=C5=⋯=0\mathrm{C}_1=\mathrm{C}_3=\mathrm{C}_5=\cdots =0\,

ℏ22m⋅2C2=−(E−12ℏω)C0=−(E−12ℏω)⇓C2=−mℏ(E−12ℏω),C4≠0,C6≠0\begin{aligned} \frac{\hbar^2}{2m}\cdot 2\mathrm{C}_2=-\left(E-\frac{1}{2}\hbar\omega\right)&\mathrm{C}_0=-\left(E-\frac{1}{2}\hbar\omega\right)\\ \Downarrow\quad&\\ \mathrm{C}_2=-\frac{m}{\hbar}\left(E-\frac{1}{2}\hbar\omega\right)&\text{,}\mathrm{C}_4\ne 0\text{,}\mathrm{C}_6\ne 0 \end{aligned}

χ偶(x)=C0+C2x2+C4x4+⋯\chi_{\text{偶}}\left(x\right)=\mathrm{C}_0+\mathrm{C}_2x^2+\mathrm{C}_4x^4+\cdots

当∣x∣→∞,χ偶(x)=exp⁡(−mω0ℏx2)\text{当}\left|x\right|\rightarrow\infty\text{,}\chi_{\text{偶}}\left(x\right)=\exp\left(-\frac{m\omega_{0}}{\hbar}x^2\right)

ϕ(x)=exp⁡(−mω02ℏx2)χ偶(x)∣x∣→∞→{exp⁡(−mω02ℏx2)exp⁡(mω0ℏx2)exp⁡(mω02ℏx2)(发散需截断)\phi\left(x\right)=\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi_{\text{偶}}\left(x\right) \underrightarrow{\left|x\right|\rightarrow\infty} \begin{cases} \exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\exp\left(\frac{m\omega_{0}}{\hbar}x^2\right)\\ \exp\left(\frac{m\omega_{0}}{2\hbar}x^2\right)\text{(发散需截断)} \end{cases}

令 E−12ℏω=0E-\frac{1}{2}\hbar\omega=0 ,E0=12ℏωE_0=\frac{1}{2}\hbar\omega ,则 C2=0,C4=0,C6=0,⋯\mathrm{C}_2=0,\mathrm{C}_4=0,\mathrm{C}_6=0,\cdots

χ0(x)=C0+C2x2+C4x4+⋯=C0=1\chi_0\left(x\right)=\mathrm{C}_0+\mathrm{C}_2x^2+\mathrm{C}_4x^4+\cdots =\mathrm{C}_0=1

ϕ(x)=A0exp⁡(−mω02ℏx2)χ0(x)=A0exp⁡(−mω02ℏx2),A0为归一化系数。\phi\left(x\right)=A_{0}\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\chi_{0}\left(x\right)=A_0\exp\left(-\frac{m\omega_{0}}{2\hbar}x^2\right)\text{,}A_{0}\text{为归一化系数。}

  1. 取 C0=0\mathrm{C}_0=0 ,C1=1\mathrm{C}_1=1 则 C0=C2=C4=⋯=0\mathrm{C}_0=\mathrm{C}_2=\mathrm{C}_4=\cdots=0

3ℏ2mC3=−(E−32ℏω)E=32ℏω(需截断否则发散)\frac{3\hbar^2}{m}\mathrm{C}_3=-\left(E-\frac{3}{2}\hbar\omega\right)\quad\quad E=\frac{3}{2}\hbar\omega\text{(需截断否则发散)}

C3=0,C5=0,⋯\mathrm{C}_3=0,\mathrm{C}_5=0,\cdots

χ1(x)=x\chi_1\left(x\right)=x

结论:当 EE 取值为 En=(n+12)ℏω0E_n=\left(n+\frac{1}{2}\right)\hbar\omega_0 ,n=0,1,2,⋯n=0,1,2,\cdots 时,χ(x)\chi\left(x\right) 的级数解会在某一阶截断,成为厄密多项式 Hn(y)H_n\left(y\right)(为实函数),χn(x)=Hn(y)=Hn(mω0ℏx)\chi_n\left(x\right)=H_n\left(y\right)=H_n\left(\sqrt{\frac{m\omega_{0}}{\hbar}}x\right) 。

ϕn(x)=Anexp⁡(−mω02ℏx2)Hn(mω0ℏx),An=[mω0ℏπ2nn!]12\phi_n\left(x\right)=A_n\exp\left(-\frac{m\omega_0}{2\hbar}x^2\right)H_n\left(\sqrt{\frac{m\omega_0}{\hbar}}x\right)\text{,}A_n=\left[\frac{\sqrt{\frac{m\omega_0}{\hbar}}}{\sqrt{\pi}2^nn!}\right]^{\frac{1}{2}}

引入归一化系数 AnA_n ,使得波函数归一化,即

ϕn(x)=Anexp⁡(−mω02ℏx2)Hn(mω0ℏx),An=[mω0ℏπ2nn!]12\phi_n\left(x\right)=A_n\exp\left(-\frac{m\omega_0}{2\hbar}x^2\right)H_n\left(\sqrt{\frac{m\omega_0}{\hbar}}x\right)\text{,}A_n=\left[\frac{\sqrt{\frac{m\omega_0}{\hbar}}}{\sqrt{\pi}2^nn!}\right]^{\frac{1}{2}}

∫−∞∞ϕn∗(x)ϕn(x)dx=∫−∞∞∣ϕn(x)∣2dx=1\int_{-\infty}^{\infty}{\phi_{n}^{*}\left(x\right)\phi_{n}\left(x\right)\mathrm{d}x=}\int_{-\infty}^{\infty}{\left|\phi_{n}\left(x\right)\right|^2\mathrm{d}x=}1

(ϕm,ϕn)=∫−∞∞ϕm∗(x)ϕn(x)dx=0(若m≠n)\left(\phi_{m},\phi_{n}\right)=\int_{-\infty}^{\infty}{\phi_{m}^{*}\left(x\right)\phi_{n}\left(x\right)\mathrm{d}x=}0\quad\left(\text{若}m\ne n\right)

∫−∞∞ϕm∗(x)ϕn(x)dx=δmn\int_{-\infty}^{\infty}{\phi_{m}^{*}\left(x\right)\phi_{n}\left(x\right)\mathrm{d}x=}\delta_{mn}

引入记号 α=mω0ℏ\alpha=\sqrt{\frac{m\omega_0}{\hbar}}

ϕ0(x)=απ1/4e−α2x22,E0=12ℏω0ϕ1(x)=2απ1/4αxe−α2x22,E1=32ℏω0ϕ2(x)=α/2π1/4(2α2x2−1)e−α2x22,E2=52ℏω0}En=(n+12)ℏω0\begin{rcases} \phi_0\left(x\right)=\frac{\sqrt{\alpha}}{\pi^{1/4}}e^{-\frac{\alpha^2x^2}{2}}, & E_0=\frac{1}{2}\hbar\omega_{0}\\ \phi_1\left(x\right)=\frac{\sqrt{2\alpha}}{\pi^{1/4}}\alpha xe^{-\frac{\alpha^2x^2}{2}}, & E_1=\frac{3}{2}\hbar\omega_{0}\\ \phi_2\left(x\right)=\frac{\sqrt{\alpha /2}}{\pi^{1/4}}\left(2\alpha^2x^2-1\right)e^{-\frac{\alpha ^2x^2}{2}}, & E_2=\frac{5}{2}\hbar\omega_{0} \end{rcases} E_n=\left(n+\frac{1}{2}\right)\hbar\omega_{0}

E0=12ℏω0E_0=\frac{1}{2}\hbar\omega_{0} 为零点能。

估算:

x∼0⇒xˉ∼(Δx)2‾,Px∼0⇒Px‾∼(ΔPx)2‾x\sim 0\Rightarrow\bar{x}\sim\sqrt{\overline{\left(\varDelta x\right)^2}},\quad P_x\sim 0\Rightarrow\overline{P_x}\sim\sqrt{\overline{\left(\varDelta P_x\right)^2}}

xˉ⋅Px‾∼(Δx)2‾(ΔPx)2‾∼14ℏ2=12ℏ\bar{x}\cdot\overline{P_x}\sim\sqrt{\overline{\left(\varDelta x\right)^2}\overline{\left(\varDelta P_x\right)^2}}\sim\sqrt{\frac{1}{4}\hbar^2}=\frac{1}{2}\hbar

⇒xˉ=ℏ21Px‾\Rightarrow\bar{x}=\frac{\hbar}{2}\frac{1}{\overline{P_x}}

E≅P2‾2m+12mω02xˉ2=P2‾2m+12mω02ℏ241P2‾E\cong\frac{\overline{P^2}}{2m}+\frac{1}{2}m\omega_{0}^{2}\bar{x}^2=\frac{\overline{P^2}}{2m}+\frac{1}{2}m\omega_{0}^{2}\frac{\hbar ^2}{4}\frac{1}{\overline{P^2}}

令 t=P2‾t=\overline{P^2} ,

dE(t)dt=12m+18mω02ℏ2(−1t2)=0\frac{dE\left(t\right)}{dt}=\frac{1}{2m}+\frac{1}{8}m\omega_{0}^{2}\hbar^2\left(-\frac{1}{t^2}\right)=0

12m=−18mω02ℏ2(1t2)\frac{1}{2m}=-\frac{1}{8}m\omega_{0}^{2}\hbar^2\left(\frac{1}{t^2}\right)

⇒t2=(12mω0ℏ)2⇒tmin⁡=12mω0ℏ\Rightarrow t^2=\left(\frac{1}{2}m\omega_{0}\hbar\right)^2\Rightarrow t_{\min}=\frac{1}{2}m\omega_{0}\hbar

Emin⁡=E(tmin⁡)=tmin⁡2m+12mω02ℏ24tmin⁡=12mmω0ℏ2+12mω02ℏ24mω0ℏ2=14ℏω0+14ℏω0=12ℏω0\begin{aligned} E_{\min}&=E\left(t_{\min}\right)\\ &=\frac{t_{\min}}{2m}+\frac{1}{2}m\omega_{0}^{2}\frac{\hbar^2}{4t_{\min}}\\ &=\frac{1}{2m}\frac{m\omega_{0}\hbar}{2}+\frac{1}{2}m\omega_{0}^{2}\frac{\hbar^2}{4\frac{m\omega_{0}\hbar}{2}}\\ &=\frac{1}{4}\hbar\omega_{0}+\frac{1}{4}\hbar\omega_{0}\\ &=\frac{1}{2}\hbar\omega_{0} \end{aligned}

§ 一维定态体系的一些基本性质

定理一

设势能函数 V(x)V(x) 为实函数。若 ψ(x)\psi(x) 是本征方程的解,则 ψ∗(x)\psi^{*}(x) 也是一个解,对应本征值也是 EE(也适用于二维)。

证明:

若V∗(x)=V∗(x),即V∗(x)为实函数\text{若}V^{*}(x)=V^{*}(x), \text{即}V^{*}(x)\text{为实函数}

Eψ(x)=−ℏ22md2dx2ψ(x)+V(x)ψ(x)E\psi(x)=-\frac{\hbar^{2}}{2m}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}\psi(x)+V(x)\psi(x)

[Eψ(x)]∗=Eψ∗(x)=[−ℏ22md2dx2ψ(x)+V(x)ψ(x)]∗=−ℏ22md2dx2ψ∗(x)+V(x)ψ∗(x)[E\psi(x)]^{*}=E\psi^{*}(x)=[-\frac{\hbar^{2}}{2m}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}\psi(x)+V(x)\psi(x)]^{*}=-\frac{\hbar^{2}}{2m}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}\psi^{*}(x)+V(x)\psi^{*}(x)

推论

若方程的本征值 EE 是非简并的,则相应的本征函数可以取作实函数。
非简并就是方程只有一个线性无关的解,一定有 ψ∗(x)=Cψ(x)\psi^{*}(x)=C\psi(x) ,

∫∣ψ∗(x)∣2dx=∫∣Cψ(x)∣2dx=∣C∣2∫∣ψ(x)∣2dx\int|\psi^{*}(x)|^{2}\mathrm{d}x=\int|C\psi(x)|^{2}\mathrm{d}x=|C|^{2}\int|\psi(x)|^{2}\mathrm{d}x

假设 ∫∣ψ(x)∣2dx=1\int|\psi(x)|^{2}\mathrm{d}x=1 ,因 ∫∣ψ∗(x)∣2dx=1\int|\psi^{*}(x)|^{2}\mathrm{d}x=1 ,则 ∣C∣2=1|C|^{2}=1 ,C=eiαC=e^{i\alpha} ,ψ∗=ψ∣eiα∣2\psi^{*}=\psi|e^{i\alpha}|^2 。


定理二

设势能函数 V(x)V(x) 为实函数,则对应于任何方程的本征值 EE ,总可以找到方程的一组实函数解。凡是属于本征值 EE 的任何一个解,都可以表示成这组实函数解的线性叠加。

假设 ψ(x)\psi(x),ψ∗(x)\psi^{*}(x) 线性无关,即 C1ψ(x)+C2ψ∗(x)=0C_1\psi(x)+C_2\psi^{*}(x)=0 ,则只有 C1=0,C2=0C_1=0,\;C_2=0
定义 φ(x)=ψ(x)+ψ∗(x),χ(x)=−i[ψ(x)−ψ∗(x)]\varphi(x)=\psi(x)+\psi^{*}(x),\;\chi(x)=-i\left[\psi(x)-\psi^{*}(x)\right](两函数都是实的,取复共轭等于其自身)
假设 D1φ+D2χ=0D_1\varphi+D_2\chi=0

D1[ψ(x)+ψ∗(x)]+D2{−i[ψ(x)−ψ∗(x)]}=0D1ψ+D1ψ∗−iD2ψ+iD2ψ∗=0(D1−iD2)ψ+(D1+iD2)ψ∗=0\begin{aligned} D_1\left[\psi(x)+\psi^{*}(x)\right]+D_2\{-i\left[\psi(x)-\psi^{*}(x)\right]\}&=0\\ D_1\psi+D_1\psi^{*}-iD_2\psi+iD_2\psi^{*}&=0\\ \left(D_1-iD_2\right)\psi+\left(D_1+iD_2\right)\psi^{*}&=0 \end{aligned}

因为 ψ\psi,ψ∗\psi^{*} 线性无关,则

{D1−iD2=0D1+iD2=0⇒{D1=0D2=0\begin{cases} D_1-iD_2=0\\ D_1+iD_2=0 \end{cases} \Rightarrow \begin{cases} D_1=0\\ D_2=0 \end{cases}


定理三

设势能函数 V(x)V(x) 具有空间反射对称性( V(−x)=V(x)V(-x)=V(x) )。若 ψ(x)\psi(x) 是方程的一个解,则 ψ(−x)\psi(-x) 亦是该方程的一个解。

证明:令 x′=−xx^{\prime}=-x ,

Eψ(−x)=−ℏ22md2dx2ψ(−x)+V(x)ψ(−x)=−ℏ22md2dx2ψ(x′)+V(−x′)ψ(x′)⇓V(x)具有空间反射对称性。=−ℏ22md2dx2ψ(x′)+V(x′)ψ(x′)=Eψ(x′)\begin{aligned} E\psi(-x)&=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi(-x)+V(x)\psi(-x)\\ &=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi(x^{\prime})+V(-x^{\prime})\psi(x^{\prime})\\ &\Downarrow V(x)\text{具有空间反射对称性。}\\ &=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi(x^{\prime})+V(x^{\prime})\psi(x^{\prime})\\ &=E\psi(x^{\prime}) \end{aligned}

当 EE 为非简并时,有 P^ψn(x)=ψn(−x)=Cψ(x)\hat{P}\psi_n(x)=\psi_n(-x)=\mathrm{C}\psi(x)
经过两次反射后,P^(P^ψ(x))=P^(ψ(−x))=P^(Cψ(x))=C2ψ(x)\hat{P}\left(\hat{P}\psi(x)\right)=\hat{P}\left(\psi(-x)\right)=\hat{P}\left(\mathrm{C}\psi(x)\right)=\mathrm{C}^2\psi(x)
因此 C2=1\mathrm{C}^2=1 或 C=±1\mathrm{C}=\pm 1 ,C\mathrm{C} 为波函数的宇称值,P^\hat{P} 称为++宇称算符++{.wavy}。

若 P^ψn(x)=ψn(−x)\hat{P}\psi_n(x)=\psi_{n}(-x) ,则称 ψn(x)\psi_{n}(x) 具有偶宇称;若 P^ψn(x)=−ψn(x)\hat{P}\psi_{n}(x)=-\psi_{n}(x) ,则称 ψn(x)\psi_{n}(x) 具有奇宇称。
对于一维谐振子,由于它的每一条能级都是非简并的,因此相应的本征函数都具有确定的宇称。En=(n+12)ℏω0E_n=(n+\frac{1}{2})\hbar\omega_{0} ,P^φn(x)=φn(−x)=(−1)nφn(x)\hat{P}\varphi_n(x)=\varphi_n(-x)=(-1)^{n}\varphi_n(x) 。


定理四

设势能函数 V(x)V(x) 具有空间反射对称性。则对应于任何一个本征值 EE ,我们总可以找到方程的一组解,它们的每一个都具有确定的宇称。而且,任何对应于本征值 EE 的解 ψ(x)\psi(x) 都可以按照它们来展开。

证明:ψ(−x)=Cψ(x)\psi(-x)=\mathrm{C}\psi(x),ψ(x)\psi(x)和ψ(−x)\psi(-x) 是线性无关的,即 C1ψ(x)+C2ψ(−x)=0(C1=0,C2=0)\mathrm{C}_1\psi(x)+\mathrm{C}_2\psi(-x)=0\;(\mathrm{C}_1=0,\mathrm{C}_2=0)
定义

ϕ(x)=ψ(x)+ψ(−x)ϕ(−x)=ψ(−x)+ψ(−(−x))=ψ(−x)+ψ(x)=ϕ(x)\begin{aligned} \phi(x)&=\psi(x)+\psi(-x)\\ \phi(-x)&=\psi(-x)+\psi(-(-x))=\psi(-x)+\psi(x)=\phi(x) \end{aligned}

ϕ(x)\phi(x)和ϕ(−x)\phi(-x) 也是线性无关的,上面看出 ϕ(x)\phi(x) 具有偶宇称。
而在定理二的证明中,χ(−x)=ψ(−x)−ψ(−(−x))=ψ(−x)−ψ(x)=−χ(x)\chi(-x)=\psi(-x)-\psi(-(-x))=\psi(-x)-\psi(x)=-\chi(x),即χ(x)\chi(x)具有奇宇称。


定理五

对于一维空间的 Schrödinger 方程,若 ψ1(x)\psi_{1}(x) 和 ψ2(x)\psi_{2}(x) 是对应于同一本征值 EE 的两个解,则有 ψ1(x)ψ2′(x)−ψ1′(x)ψ2(x)=常数\psi_{1}(x)\psi_{2}^{\prime}(x)-\psi_{1}^{\prime}(x)\psi_{2}(x)=\text{常数}。

证明:

{Eψ1(x)=−ℏ22md2dx2ψ1(x)+V(x)ψ1(x)(1)Eψ2(x)=−ℏ22md2dx2ψ2(x)+V(x)ψ2(x)(2)\begin{cases} E\psi_{1}(x)=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_{1}(x)+V(x)\psi_{1}(x) & \text{(1)}\\ E\psi_{2}(x)=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_{2}(x)+V(x)\psi_{2}(x) & \text{(2)} \end{cases}

(1)×ψ2(x)−(2)×ψ1(x)\text{(1)}\times\psi_{2}(x)-\text{(2)}\times\psi_{1}(x)得

−ℏ22m[ψ2(x)d2dx2ψ1(x)−ψ1(x)d2dx2ψ2(x)]=0-\frac{\hbar^2}{2m}\left[\psi_{2}(x)\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_{1}(x)-\psi_{1}(x)\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_{2}(x)\right]=0

从方程中得

ddx[ψ1(x)ψ2′(x)−ψ1′(x)ψ2(x)]=0\frac{\mathrm{d}}{\mathrm{d}x}\left[\psi_{1}(x)\psi_{2}^{\prime}(x)-\psi_{1}^{\prime}(x)\psi_{2}(x)\right]=0

因此

ψ1(x)ψ2′(x)−ψ1′(x)ψ2(x)=常数(仅一维空间)\psi_{1}(x)\psi_{2}^{\prime}(x)-\psi_{1}^{\prime}(x)\psi_{2}(x)=\text{常数(仅一维空间)}


定理六

假设粒子在一个一维空间中规则(regular)(即没有奇异点)的势函数 V(x)V(x) 中运动,则其束缚态必为非简并的。

证明:势函数 V(x)V(x) 为规则的 ⇒\Rightarrow Schrödinger 方程的任何一个解 ψ(x)\psi(x) 及其一阶导数 ψ′(x)\psi^{\prime}(x) 在空间的每一点都连续。

V(x)没有太多奇异性(regular)⇒{(1)不包含δ函数(2)V(x)有有限个阶跃V(x)\text{没有太多奇异性(regular)}\Rightarrow\begin{cases} \text{(1)不包含}\delta\text{函数}\\ \text{(2)}V(x)\text{有有限个阶跃} \end{cases}

假设 EnE_n 是粒子的一个束缚态的本征值且简并,则可以找到至少两个本征函数 ψ1(x)\psi_{1}(x) 和 ψ2(x)\psi_{2}(x) 与它对应。
另一方面根据定理五,有

ψ1(x)ψ2′(x)−ψ1′(x)ψ2(x)=常数\psi_{1}(x)\psi_{2}^{\prime}(x)-\psi_{1}^{\prime}(x)\psi_{2}(x)=\text{常数}

对于束缚态,有渐进关系

lim⁡∣x∣→∞ψ1(x)=lim⁡∣x∣→∞ψ2(x)=0\lim_{|x|\to\infty}\psi_{1}(x)=\lim_{|x|\to\infty}\psi_{2}(x)=0

则常数为零,方程改写为

ψ1′(x)ψ1(x)=ψ2′(x)ψ2(x)1ψ1dψ1dx=1ψ2dψ2dxddxln⁡ψ1=ddxln⁡ψ2ln⁡ψ1=ln⁡ψ2+C或ψ1(x)=eCψ2(x)=Aψ2\begin{aligned} \frac{\psi_{1}^{\prime}(x)}{\psi_{1}(x)}&=\frac{\psi_{2}^{\prime}(x)}{\psi_{2}(x)}\\ \frac{1}{\psi_{1}}\frac{\mathrm{d}\psi_{1}}{\mathrm{d}x}&=\frac{1}{\psi_{2}}\frac{\mathrm{d}\psi_{2}}{\mathrm{d}x}\\ \frac{\mathrm{d}}{\mathrm{d}x}\ln\psi_{1}&=\frac{\mathrm{d}}{\mathrm{d}x}\ln\psi_{2}\\ \ln\psi_{1}&=\ln\psi_{2}+C\\ \text{或}\psi_{1}(x)&=e^{C}\psi_{2}(x)=A\psi_{2} \end{aligned}

即 ψ1(x)\psi_{1}(x) 与 ψ2(x)\psi_{2}(x) 线性相关,这与假设不符,因此得证。

总结:

{(1)定理一、二:势函数为实函数的情况;(2)定理三、四:势函数具有空间反射对称性的情况;(3)定理五、六:仅适用于一维空间,高维空间不成立。\begin{cases} \text{(1)定理一、二:势函数为实函数的情况;}\\ \text{(2)定理三、四:势函数具有空间反射对称性的情况;}\\ \text{(3)定理五、六:仅适用于一维空间,高维空间不成立。} \end{cases}


§ 一维方势阱

在量子力学中,一般给定势能函数较复杂,实际中可以近似于方势阱,其结果在定性上与真实出入不大。

图3-1 一维谐振子势函数

V(x)=12kx2V(x)=\frac{1}{2}kx^{2}

图3-2 一维谐振子势函数

V(x)=ch⁡x=ex+e−x2V(x)=\ch{x}=\frac{e^{x}+e^{-x}}{2}

§ 无穷深一维方势阱

对于一维谐振子、双曲余弦函数的近似。
图3-3

V(x)={+∞,x⩽00,0⩽x⩽a−∞,x⩾aV(x)= \begin{cases} +\infty, & x\leqslant 0\\ 0, & 0\leqslant x\leqslant a\\ -\infty, & x\geqslant a \end{cases}

阱内:

Eψ(x)=−ℏ22md2dx2ψ(x)E\psi(x)=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi(x)

阱外:令 ψ(x=0)=0\psi(x=0)=0 ,ψ(x=a)=0\psi(x=a)=0

Eψ(x)+ℏ22md2dx2ψ(x)=0d2ψ(x)dx2+2mEℏ2ψ(x)=0\begin{aligned} E\psi(x)+\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi(x)&=0\\ \frac{\mathrm{d}^2\psi(x)}{\mathrm{d}x^2}+\frac{2mE}{\hbar^2}\psi(x)&=0 \end{aligned}

令解 ψ=eλx\psi=e^{\lambda x} ,ψ′=λeλx\psi^{\prime}=\lambda e^{\lambda x} ,ψ′′=λ2eλx\psi^{\prime\prime}=\lambda^{2} e^{\lambda x}\,

λ2eλx+2mEℏ2eλx=0λ=−2mEℏ2\begin{aligned} \lambda^{2} e^{\lambda x}&+\frac{2mE}{\hbar^2}e^{\lambda x}=0\\ \lambda&=-\frac{2mE}{\hbar^2} \end{aligned}

解得:

ψ1=exp⁡(−i2mEℏ2x),ψ2=exp⁡(i2mEℏ2x)\psi_{1}=\exp\left(-i\sqrt{\frac{2mE}{\hbar^2}}x\right)\text{,}\; \psi_{2}=\exp\left(i\sqrt{\frac{2mE}{\hbar^2}}x\right)

因为 ψ1(x=0)=1\psi_{1}(x=0)=1 ,ψ1(x=a)=exp⁡(−i2mEℏ2x)≠0\psi_{1}(x=a)=\exp\left(-i\sqrt{\frac{2mE}{\hbar^2}}x\right)\neq 0
ψ1\psi_{1} 和 ψ2\psi_{2} 线性无关,取其线性组合(原因:Schrödinger 方程是线性的)。

ψ=C1ψ1+C2ψ2=C1exp⁡(−i2mEℏ2x)+C2exp⁡(i2mEℏ2x)\psi=C_{1}\psi_{1}+C_{2}\psi_{2}=C_{1}\exp\left(-i\sqrt{\frac{2mE}{\hbar^2}}x\right)+C_{2}\exp\left(i\sqrt{\frac{2mE}{\hbar^2}}x\right)

  1. 由条件 ψ(x=0)=0\psi(x=0)=0 得,C1e0+C2e0=C1+C2⟹C1=−C2C_{1}e^{0}+C_{2}e^{0}=C_{1}+C_{2}\implies C_{1}=-C_{2}\,

ψ(x)=C1ψ1+(−C1)ψ2=C1(ψ1−ψ2)=C1[exp⁡(−i2mEℏ2x)−exp⁡(i2mEℏ2x)]=−2iC1[exp⁡(−i2mEℏ2x)−exp⁡(i2mEℏ2x)−2i]=(−2iC1)sin⁡2mEℏ2x=Dsin⁡2mEℏ2x\begin{aligned} \psi(x)&=C_{1}\psi_{1}+(-C_{1})\psi_{2}\\ &=C_{1}(\psi_{1}-\psi_{2})\\ &=C_{1}\left[\exp\left(-i\sqrt{\frac{2mE}{\hbar^2}}x\right)-\exp\left(i\sqrt{\frac{2mE}{\hbar^2}}x\right)\right]\\ &=-2iC_{1}\left[\frac{\exp\left(-i\sqrt{\frac{2mE}{\hbar^2}}x\right)-\exp\left(i\sqrt{\frac{2mE}{\hbar^2}}x\right)}{-2i}\right]\\ &=(-2iC_{1})\sin\sqrt{\frac{2mE}{\hbar^2}}x\\ &=D\sin\sqrt{\frac{2mE}{\hbar^2}}x \end{aligned}

sin⁡x=eix−e−ix2icos⁡x=eix+e−ix2\begin{aligned} \sin{x}&=\frac{e^{ix}-e^{-ix}}{2i}\\ \cos{x}&=\frac{e^{ix}+e^{-ix}}{2} \end{aligned}

  1. 由条件 ψ(x=a)=0\psi(x=a)=0 得,

ψ(x=a)=Dsin⁡2mEℏ2a\psi(x=a)=D\sin\sqrt{\frac{2mE}{\hbar^2}}a

令 2mEℏ2a=nπ,n=1,2,3,⋯\sqrt{\frac{2mE}{\hbar^2}}a=n\pi,\;n=1,2,3,\cdots

En=n2π2ℏ22ma2,ψn(x)=Dnsin⁡nπaxE_{n}=\frac{n^2\pi^2\hbar^2}{2ma^2},\quad \psi_{n}(x)=D_{n}\sin\frac{n\pi}{a}x

∫0n∣ψn(x)∣2dx=1∫0n∣Dn∣2sin⁡2nπaxdx=1=∣Dn∣2∫0nsin⁡2nπaxdx=∣Dn∣2∫0n12(1−cos⁡2nπax)dx=∣Dn∣2a2\begin{aligned} \int_{0}^{n}\left|\psi_{n}(x)\right|^2\mathrm{d}x&=1\\ \int_{0}^{n}\left|D_{n}\right|^2\sin^2\frac{n\pi}{a}x\mathrm{d}x&=1\\ &=\left|D_{n}\right|^2\int_{0}^{n}\sin^2\frac{n\pi}{a}x\mathrm{d}x\\ &=\left|D_{n}\right|^2\int_{0}^{n}\frac{1}{2}\left(1-\cos\frac{2n\pi}{a}x\right)\mathrm{d}x\\ &=\left|D_{n}\right|^2\frac{a}{2} \end{aligned}

∣Dn∣=2a\left|D_{n}\right|=\sqrt{\frac{2}{a}} ,取 Dn=2aD_{n}=\sqrt{\frac{2}{a}} ,
最终得到,

ψn(x)=2asin⁡nπax\boxed{\psi_{n}(x)=\sqrt{\frac{2}{a}}\sin\frac{n\pi}{a}x}

§ 有限深对称方势阱

核力相互作用势的近似。

图3-4

将核力相互作用分开需要 2∼8MeV2\sim 8MeV 的能量,强力是短程力 a=10−13cma=10^{-13}cm( 11 费米),分子的结合能 V0=10eVV_{0}=10eV 。

从物理学角度看,所谓爱情只不过是头脑中的分子作用。
“我爱你有多深(不超过 10−8cm10^{-8}cm ),我爱你有几分(不多于 10eV10eV )”

图3-5

V(x)={0,−a2⩽x⩽a2V0,∣x∣⩾a2V(x)= \begin{cases} 0,\;-\frac{a}{2}\leqslant x\leqslant \frac{a}{2}\\ V_{0},\qquad\;|x|\geqslant \frac{a}{2} \end{cases}

{1 区:−ℏ22md2dx2ψ1+V0ψ1=Eψ12 区:−ℏ22md2dx2ψ2=Eψ23 区:−ℏ22md2dx2ψ3+V0ψ3=Eψ3\begin{cases} \text{1 区:} & -\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_1+V_0\psi_1=E\psi_1\\[0.3em] \text{2 区:} & -\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_2=E\psi_2 \\[0.3em] \text{3 区:} & -\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_3+V_0\psi_3=E\psi_3 \end{cases}

2 区:

ψ2=C1exp⁡(i2mEℏ2x)+C2exp⁡(−i2mEℏ2x)=C1eikx+C2e−ikx,k=2mEℏ2\begin{aligned} \psi_2&=\mathrm{C}_1\exp{\left(i\sqrt{\frac{2mE}{\hbar^2}}x\right)}+\mathrm{C}_2\exp{\left(-i\sqrt{\frac{2mE}{\hbar^2}}x\right)}\\ &=\mathrm{C}_1e^{ikx}+\mathrm{C}_2e^{-ikx},\quad k=\sqrt{\frac{2mE}{\hbar^2}} \end{aligned}

1 区:

d2dx2ψ1(x)+2mℏ2(E−V0)⏟⩽0ψ1(x)=0\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_1(x)+\frac{2m}{\hbar^2}\underbrace{(E-V_0)}_{\leqslant 0}\psi_1(x)=0

d2dx2ψ1(x)−2mℏ2(V0−E)ψ1(x)=0\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_1(x)-\frac{2m}{\hbar^2}(V_0-E)\psi_1(x)=0

令解 ψ1=eλx,ψ1′=λeλx,ψ1′′=λ2eλx\psi_1=e^{\lambda x},\;\psi_1^{\prime}=\lambda e^{\lambda x},\;\psi_1^{\prime\prime}=\lambda^2 e^{\lambda x}\,

λ2eλx−2mℏ2(V0−E)eλx=0\lambda^2 e^{\lambda x}-\frac{2m}{\hbar^2}(V_0-E)e^{\lambda x}=0

λ2=2mℏ2(V0−E)⩾0\lambda^2=\frac{2m}{\hbar^2}(V_0-E)\geqslant 0

λ1=2mℏ(V0−E)=β,λ2=−2mℏ(V0−E)=−β(两个线性无关的解)\lambda_1=\sqrt{\frac{2m}{\hbar}(V_0-E)}=\beta,\quad \lambda_2=-\sqrt{\frac{2m}{\hbar}(V_0-E)}=-\beta\quad\text{(两个线性无关的解)}

ψ1(a)(x)=Beβx,ψ1(b)(x)=Ae−βx→x→−∞∞发散舍弃\psi_1^{(a)}(x)=\mathrm{B}e^{\beta x},\quad \psi_1^{(b)}(x)=\mathrm{A}e^{-\beta x}\xrightarrow{x\to -\infty}\infty\text{发散舍弃}

方程整理为:

{ψ1=Beβxψ2=C1eikx+C2e−ikxψ3=Ae−βx\begin{cases} \psi_1=\mathrm{B}e^{\beta x} \\ \psi_2=\mathrm{C}_1e^{ikx}+\mathrm{C}_2e^{-ikx} \\ \psi_3=\mathrm{A}e^{-\beta x} \end{cases}

对于阱右边界,

{ψ2(x=a2)=ψ3(x=a2)ψ2′(x=a2)=ψ3′(x=a2)(波函数在边界连续)\begin{cases} \psi_2(x=\frac{a}{2})=\psi_3(x=\frac{a}{2}) \\[0.5em] \psi_2^{\prime}(x=\frac{a}{2})=\psi_3^{\prime}(x=\frac{a}{2}) \end{cases} \text{(波函数在边界连续)}

{C1eika2+C2e−ika2=Ae−βa2(1)C1(ik)eika2−C2(ik)e−ika2=−Aβe−βa2(2)\begin{cases} \mathrm{C}_1e^{ik\frac{a}{2}}+\mathrm{C}_2e^{-ik\frac{a}{2}}=\mathrm{A}e^{-\beta\frac{a}{2}} & \text{(1)}\\[0.5em] \mathrm{C}_1(ik)e^{ik\frac{a}{2}}-\mathrm{C}_2(ik)e^{-ik\frac{a}{2}}=-\mathrm{A}\beta e^{-\beta\frac{a}{2}} & \text{(2)} \end{cases}

(2)(1):ik(C1eika2−C2e−ika2)C1eika2+C2e−ika2=−Aβe−βa2Ae−βa2=−β\frac{\text{(2)}}{\text{(1)}}\text{:} \frac{ik\left(\mathrm{C}_1e^{ik\frac{a}{2}}-\mathrm{C}_2e^{-ik\frac{a}{2}}\right)}{\mathrm{C}_1e^{ik\frac{a}{2}}+\mathrm{C}_2e^{-ik\frac{a}{2}}} =\frac{-\mathrm{A}\beta e^{-\beta\frac{a}{2}}}{\mathrm{A}e^{-\beta\frac{a}{2}}} =-\beta

对于阱左边界,

{ψ1(x=−a2)=ψ2(x=−a2)ψ1′(x=−a2)=ψ2′(x=−a2)\begin{cases} \psi_1(x=-\frac{a}{2})=\psi_2(x=-\frac{a}{2}) \\[0.5em] \psi_1^{\prime}(x=-\frac{a}{2})=\psi_2^{\prime}(x=-\frac{a}{2}) \end{cases}

同理

β=ikC1e−ika2−C2eika2C1e−ika2+C2eika2\beta=ik\frac{\mathrm{C}_1e^{-ik\frac{a}{2}}-\mathrm{C}_2e^{ik\frac{a}{2}}}{\mathrm{C}_1e^{-ik\frac{a}{2}}+\mathrm{C}_2e^{ik\frac{a}{2}}}

比较两式得,C1=±C2\mathrm{C}_1=\pm\mathrm{C}_2\,

  1. 当 C1=C2\mathrm{C}_1=\mathrm{C}_2 时,

    ψ2=C1eikx+C2e−ikx=C1(eikx+e−ikx)\psi_2=\mathrm{C}_1e^{ikx}+\mathrm{C}_2e^{-ikx}=\mathrm{C}_1\left(e^{ikx}+e^{-ikx}\right)

    ψ(x)=ψ(−x),为偶宇称\psi(x)=\psi(-x),\;\text{为偶宇称}

    ψ1(−x)=ψ3(x)ψ3(−x)=ψ1(x)}⟹A=B\begin{rcases} \psi_1(-x)=\psi_3(x) \\ \psi_3(-x)=\psi_1(x) \end{rcases}\implies \mathrm{A}=\mathrm{B}

    {ψ1=Aeβxψ2=C1(eikx+e−ikx)ψ3=Ae−βx\begin{cases} \psi_1=\mathrm{A}e^{\beta x} \\ \psi_2=\mathrm{C}_1\left(e^{ikx}+e^{-ikx}\right) \\ \psi_3=\mathrm{A}e^{-\beta x} \end{cases}

    β=ikC1e−ika2−C2eika2C1eika2+C2e−ika2=2i⋅ik2e−ika2−eika22ieika2+e−ika22=−k−sin⁡ka2cos⁡ka2=ktg⁡ka2\begin{aligned} \beta&=ik\frac{\mathrm{C}_1e^{-ik\frac{a}{2}}-\mathrm{C}_2e^{ik\frac{a}{2}}}{\mathrm{C}_1e^{ik\frac{a}{2}}+\mathrm{C}_2e^{-ik\frac{a}{2}}} \\ &=\frac{2i\cdot ik}{2}\frac{\frac{e^{-ik\frac{a}{2}}-e^{ik\frac{a}{2}}}{2i}}{\frac{e^{ik\frac{a}{2}}+e^{-ik\frac{a}{2}}}{2}} \\ &=-k\frac{-\sin\frac{ka}{2}}{\cos\frac{ka}{2}} \\ &=k\tg\frac{ka}{2} \end{aligned}

    βa2=ka2tg⁡ka2\beta\frac{a}{2}=k\frac{a}{2}\tg\frac{ka}{2}

    令 ξ=ka2,η=βa2,k=2mEℏ2,β=2m(V0−E)ℏ2\xi=\frac{ka}{2},\;\eta=\frac{\beta a}{2},\;k=\sqrt{\frac{2mE}{\hbar^2}},\;\beta=\sqrt{\frac{2m(V_0-E)}{\hbar^2}}\,

    η=ξtg⁡ξ,ξ2+η2=(ka2)2+(βa2)2=a24(k2+β2)=a24[2mEℏ2+2m(V0−E)ℏ2]=a242mV0ℏ2=ma2V02ℏ2\begin{aligned} \eta=\xi\tg\xi,\quad \xi^2+\eta^2&=\left(\frac{ka}{2}\right)^2+\left(\frac{\beta a}{2}\right)^2 \\ &=\frac{a^2}{4}\left(k^2+\beta^2\right) \\ &=\frac{a^2}{4}\left[\frac{2mE}{\hbar^2}+\frac{2m(V_0-E)}{\hbar^2}\right] \\ &=\frac{a^2}{4}\frac{2mV_0}{\hbar^2} \\ &=\frac{ma^2V_0}{2\hbar^2} \end{aligned}

    图3-6

    圆 ξ2+η2=ma2V02ℏ2\xi^2+\eta^2=\frac{ma^2V_0}{2\hbar^2} 与正切函数 η=ξtg⁡ξ\eta=\xi\tg\xi 交点的 ξ\xi 值为解,阱深有限,圆的半径也有限,则有限个解,有限个束缚态。

    {ψ1=Aeβxψ2=C1(eikx+e−ikx)=2C1cos⁡kx=Dcos⁡kxψ3=Ae−βx\begin{cases} \psi_1=\mathrm{A}e^{\beta x} \\ \psi_2=\mathrm{C}_1\left(e^{ikx}+e^{-ikx}\right)=2\mathrm{C}_1\cos{kx}=\mathrm{D}\cos{kx} \\ \psi_3=\mathrm{A}e^{-\beta x} \end{cases}

    A,D可定出{ψ2(x=a2)=ψ3(x=a2)⟹∫−∞∞∣ψ∣2dx=∫−∞−a2∣ψ1∣2dx+∫−a2a2∣ψ2∣2dx+∫a2∞∣ψ3∣2dx=∫−∞−a2A2e2βxdx+∫−a2a2D2cos⁡2kxdx+∫a2∞A2e−2βxdx\mathrm{A},\;\mathrm{D}\text{可定出} \begin{cases} \psi_2(x=\frac{a}{2})=\psi_3(x=\frac{a}{2})\implies \\[2em] \begin{aligned} \int_{-\infty}^{\infty}|\psi|^2\,\mathrm{d}x&=\int_{-\infty}^{-\frac{a}{2}}|\psi_1|^2\,\mathrm{d}x+\int_{-\frac{a}{2}}^{\frac{a}{2}}|\psi_2|^2\,\mathrm{d}x+\int_{\frac{a}{2}}^{\infty}|\psi_3|^2\,\mathrm{d}x \\ &=\int_{-\infty}^{-\frac{a}{2}}\mathrm{A}^{2}e^{2\beta x}\,\mathrm{d}x+\int_{-\frac{a}{2}}^{\frac{a}{2}}\mathrm{D}^{2}\cos^{2}{kx}\,\mathrm{d}x+\int_{\frac{a}{2}}^{\infty}\mathrm{A}^{2}e^{-2\beta x}\,\mathrm{d}x \end{aligned} \end{cases}

  2. 当 C1=−C2\mathrm{C}_1=-\mathrm{C}_2 时,

    −β=ikC1eika2−C2e−ika2C1e−ika2+C2eika2-\beta=ik\frac{\mathrm{C}_1e^{ik\frac{a}{2}}-\mathrm{C}_2e^{-ik\frac{a}{2}}}{\mathrm{C}_1e^{-ik\frac{a}{2}}+\mathrm{C}_2e^{ik\frac{a}{2}}}

    β=ikC1e−ika2−C2eika2C1e−ika2+C2eika2\beta=ik\frac{\mathrm{C}_1e^{-ik\frac{a}{2}}-\mathrm{C}_2e^{ik\frac{a}{2}}}{\mathrm{C}_1e^{-ik\frac{a}{2}}+\mathrm{C}_2e^{ik\frac{a}{2}}}

    β=−kctg⁡ka2\beta=-k\ctg{\frac{ka}{2}}

    βa2=−ka2ctg⁡ka2⟹{η=−ξctg⁡ξξ2+η2=mV0a22ℏ2\frac{\beta a}{2}=-\frac{ka}{2}\ctg{\frac{ka}{2}}\implies \begin{cases} \eta=-\xi\ctg\xi \\[0.5em] \xi^2+\eta^2=\frac{mV_0a^2}{2\hbar^2} \end{cases}

    图3-7

    ψ2=C1(eikx−e−ikx)=Dsin⁡kx\psi_2=\mathrm{C}_1\left(e^{ikx}-e^{-ikx}\right)=\mathrm{D}\sin{kx}

    ψ1=−Aeβx,ψ3=Ae−βx\psi_1=-\mathrm{A}e^{\beta x},\;\psi_3=\mathrm{A}e^{-\beta x}

    {ψ1(−x)=ψ3(x)ψ3(−x)=ψ1(x)⟹A=−B\begin{cases} \psi_1(-x)=\psi_3(x) \\ \psi_3(-x)=\psi_1(x) \end{cases}\implies\mathrm{A}=-\mathrm{B}

    ψ(x)=Dsin⁡kx,(−a2⩽x⩽a2)\psi(x)=\mathrm{D}\sin{kx},\;\left(-\frac{a}{2}\leqslant x \leqslant\frac{a}{2}\right)

    对于方势阱,无论 V0V_0 数值多小,都会存在一个束缚态,其宇称是偶的(仅一维)。

    对于处在束缚态的粒子,区域 ∣x∣⩾a2|x|\geqslant\frac{a}{2} 叫++经典禁绝区++{.dot}(经典力学中,粒子不可能出现在该区域),但在量子力学中,波函数 ψ(x)∽e±βx\psi(x)\backsim e^{\pm\beta x} ,β=2m(V0−E)ℏ2\beta=\frac{2m(V_0-E)}{\hbar^2} ,粒子在 x0x_0 点出现几率 ρ(x)dx=∣ψ(x0)∣2dx∽e±2βx0dx\rho(x)\,\mathrm{d}x=|\psi(x_0)|^2\,\mathrm{d}x\backsim e^{\pm 2\beta x_0}\,\mathrm{d}x 不为零,尽管它随 ∣x∣|x| 的增大而指数级衰减。这一现象称为++隧道穿透++{.dot}(隧道效应)。

图1-1 弹簧振子

经典力学中:

E0=P22m+12mω02x2=12ω02xmax2⟹xmax2=E012mω02E_0=\frac{P^2}{2m}+\frac{1}{2}m\omega_{0}^{2}x^2=\frac{1}{2}\omega_{0}^{2}x_{\text{max}}^{2}\implies x_{\text{max}}^2=\frac{E_0}{\frac{1}{2}m\omega_{0}^{2}}

α=mω0ℏ,α−1=ℏmω0=12ω0ℏ12mω02=E012mω2=xmax=x禁绝\alpha=\sqrt{\frac{m\omega_0}{\hbar}},\quad\alpha^{-1}=\sqrt{\frac{\hbar}{m\omega_0}}=\sqrt{\frac{\frac{1}{2}\omega_0\hbar}{\frac{1}{2}m\omega_{0}^{2}}}=\sqrt{\frac{E_0}{\frac{1}{2}m\omega^2}}=x_{\text{max}}=x_{\text{禁绝}}

§ 方势垒的反射与穿透(贯穿)

图3-8

{1 区:Eψ1=−ℏ22md2dx2ψ1+V(x)ψ1=−ℏ22md2dx2ψ12 区:Eψ2=−ℏ22md2dx2ψ2+V0ψ23 区:Eψ3=−ℏ22md2dx2ψ3\begin{cases} \text{1 区:}E\psi_1=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_1+V(x)\psi_1=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_1 \\[0.5em] \text{2 区:}E\psi_2=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_2+V_0\psi_2 \\[0.5em] \text{3 区:}E\psi_3=-\frac{\hbar^2}{2m}\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_3 \end{cases}

1 区:

d2dx2ψ1+2mEℏ2ψ1=0\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_1+\frac{2mE}{\hbar^2}\psi_1=0

k2=2mEℏ2k^2=\frac{2mE}{\hbar^2}

ψ1(a)=eikx,ψ1(b)=e−ikx\psi_1^{(a)}=e^{ikx},\quad\psi_1^{(b)}=e^{-ikx}

行波有向前的波和反射的波:

ψ1=C1eikx+C2e−ikx=eikx+Re−ikx\psi_1=\mathrm{C}_1e^{ikx}+\mathrm{C}_2e^{-ikx}=e^{ikx}+\mathrm{R}e^{-ikx}

2 区:

d2dx2ψ2−2mEℏ2(V0−E)ψ2=0\frac{\mathrm{d}^2}{\mathrm{d}x^2}\psi_2-\frac{2mE}{\hbar^2}(V_0-E)\psi_2=0

ψ2(a)=eβx,ψ2(b)=e−βx\psi_2^{(a)}=e^{\beta x},\quad\psi_2^{(b)}=e^{-\beta x}

ψ2=Aeβx+Be−βx\psi_2=\mathrm{A}e^{\beta x}+\mathrm{B}e^{-\beta x}

3 区:

ψ3(a)=eikx,ψ3(b)=e−ikx\psi_3^{(a)}=e^{ikx},\quad\psi_3^{(b)}=e^{-ikx}

ψ1(a)=eikx(无反射波)\psi_1^{(a)}=e^{ikx}\text{(无反射波)}

ψ入射=eikx\psi_{\text{入射}}=e^{ikx}

jin=ℏ2mi[e−ikxddxeikx−eikxddxe−ikx]=ℏ2mi[ik−(−ik)]=ℏ2mi⋅2ik=ℏkm\begin{aligned} j_{in}&=\frac{\hbar}{2mi}\left[e^{-ikx}\frac{\mathrm{d}}{\mathrm{d}x}e^{ikx}-e^{ikx}\frac{\mathrm{d}}{\mathrm{d}x}e^{-ikx}\right] \\ &=\frac{\hbar}{2mi}\left[ik-(-ik)\right] \\ &=\frac{\hbar}{2mi}\cdot 2ik \\ &=\frac{\hbar k}{m} \end{aligned}

ℏkm=Pm=v入射(入射波速度,仅是量纲分析,速率在量子力学中无意义)\frac{\hbar k}{m}=\frac{P}{m}=v_{\text{入射}}\text{(入射波速度,仅是量纲分析,速率在量子力学中无意义)}

几率流密度jin=ℏ2mi[ψ∗∇ψ−ψ∇ψ∗]=ℏ2mi[ψin∗∇ψin−ψin∇ψin∗],∂ρ∂t+∇⋅j⃗=0\text{几率流密度}j_{in}=\frac{\hbar}{2mi}\left[\psi^{*}\nabla\psi-\psi\nabla\psi^{*}\right]=\frac{\hbar}{2mi}\left[\psi_{in}^{*}\nabla\psi_{in}-\psi_{in}\nabla\psi_{in}^{*}\right],\quad\frac{\partial\rho}{\partial t}+\nabla\cdot\vec{j}=0

反射系数R=j反射jin=∣R∣2vinvin=∣R∣2透射系数T=j透射jin=∣S∣2vinvin=∣S∣2\begin{aligned} \text{反射系数}\mathcal{R}=\frac{j_{\text{反射}}}{j_{in}}=|\mathrm{R}|^{2}\frac{v_{in}}{v_{in}}=|\mathrm{R}|^{2} \\ \text{透射系数}\mathcal{T}=\frac{j_{\text{透射}}}{j_{in}}=|\mathrm{S}|^{2}\frac{v_{in}}{v_{in}}=|\mathrm{S}|^{2} \end{aligned}

  1. x=0x=0 :

    {ψ1(x=0)=ψ2(x=0)ψ1′(x=0)=ψ2′(x=0)⟹{1+R=A+B(1)ik(1+R)=β(A−B)(2)\begin{cases} \psi_1(x=0)=\psi_2(x=0) \\ \psi_1^{\prime}(x=0)=\psi_2^{\prime}(x=0) \end{cases} \implies \begin{cases} 1+\mathrm{R}=\mathrm{A}+\mathrm{B} & \text{(1)} \\ ik\left(1+\mathrm{R}\right)=\beta\left(\mathrm{A}-\mathrm{B}\right) & \text{(2)} \end{cases}

    (1)+(2){A=12[(1+ikβ)+R(1−ikβ)]B=12[(1−ikβ)+R(1+ikβ)]\text{(1)+(2)} \begin{cases} \mathrm{A}=\frac{1}{2}\left[\left(1+\frac{ik}{\beta}\right)+\mathrm{R}\left(1-\frac{ik}{\beta}\right)\right] \\[1em] \mathrm{B}=\frac{1}{2}\left[\left(1-\frac{ik}{\beta}\right)+\mathrm{R}\left(1+\frac{ik}{\beta}\right)\right] \end{cases}

  2. x=ax=a :

    {ψ2(x=a)=ψ3(x=a)ψ2′(x=a)=ψ3′(x=a)⟹{Aeβa+Be−βa=Seika(3)β(Aeβa−Be−βa)=ikSeika(4)\begin{cases} \psi_2(x=a)=\psi_3(x=a) \\ \psi_2^{\prime}(x=a)=\psi_3^{\prime}(x=a) \end{cases} \implies \begin{cases} \mathrm{A}e^{\beta a}+\mathrm{B}e^{-\beta a}=\mathrm{S}e^{ika} & \text{(3)} \\ \beta\left(\mathrm{A}e^{\beta a}-\mathrm{B}e^{-\beta a}\right)=ik\mathrm{S}e^{ika} & \text{(4)} \end{cases}

    (3)+(4){A=S2(1+ikβ)eika−βaB=S2(1−ikβ)eika+βa\text{(3)+(4)} \begin{cases} \mathrm{A}=\frac{\mathrm{S}}{2}\left(1+\frac{ik}{\beta}\right)e^{ika-\beta a} \\[1em] \mathrm{B}=\frac{\mathrm{S}}{2}\left(1-\frac{ik}{\beta}\right)e^{ika+\beta a} \end{cases}

整理得:

{12[(1+ikβ)+R(1−ikβ)]=S2(1+ikβ)eika−βa(5)12[(1−ikβ)+R(1+ikβ)]=S2(1−ikβ)eika+βa(6)\begin{cases} \frac{1}{2}\left[\left(1+\frac{ik}{\beta}\right)+\mathrm{R}\left(1-\frac{ik}{\beta}\right)\right]=\frac{\mathrm{S}}{2}\left(1+\frac{ik}{\beta}\right)e^{ika-\beta a} & \text{(5)} \\[1em] \frac{1}{2}\left[\left(1-\frac{ik}{\beta}\right)+\mathrm{R}\left(1+\frac{ik}{\beta}\right)\right]=\frac{\mathrm{S}}{2}\left(1-\frac{ik}{\beta}\right)e^{ika+\beta a} & \text{(6)} \end{cases}

(5)×(1+ikβ)−(6)×(1−ikβ)\text{(5)}\times\left(1+\frac{ik}{\beta}\right)-\text{(6)}\times\left(1-\frac{ik}{\beta}\right) :

(1+ikβ)2−(1−ikβ)2=Seika[(1+ikβ)2eβa−(1−ikβ)eβa]\left(1+\frac{ik}{\beta}\right)^2-\left(1-\frac{ik}{\beta}\right)^2=\mathrm{S}e^{ika}\left[\left(1+\frac{ik}{\beta}\right)^2e^{\beta a}-\left(1-\frac{ik}{\beta}\right)e^{\beta a}\right]

双曲余弦函数:ch⁡x=ex+e−x2\ch{x}=\frac{e^x+e^{-x}}{2}\,
双曲正弦函数:sh⁡x=ex−e−x2\sh{x}=\frac{e^x-e^{-x}}{2}\,
关系:ch⁡2x+sh⁡2x=1\ch^2{x}+\sh^2{x}=1\,

Seika=4ikβ[−2+2(kβ)2]sh⁡βa+4i(kβ)ch⁡βa=−2ikβ[1−(kβ)2]sh⁡βa−2i(kβ)ch⁡βa\begin{aligned} \mathrm{S}e^{ika}&=\frac{\frac{4ik}{\beta}}{\left[-2+2\left(\frac{k}{\beta}\right)^2\right]\sh{\beta a}+4i\left(\frac{k}{\beta}\right)\ch{\beta a}} \\ &=\frac{-\frac{2ik}{\beta}}{\left[1-\left(\frac{k}{\beta}\right)^2\right]\sh{\beta a}-2i\left(\frac{k}{\beta}\right)\ch{\beta a}} \end{aligned}

T=∣S∣2=4k2β2(β2−k2)sh⁡2βa+4k2β2ch⁡2βa=4k2β2(β2−k2)2sh⁡2βa+4k2β2−4k2β2sh⁡2βach⁡2x=1−sh⁡2x→=4k2β2(β2+k2)2sh⁡2βa+4k2β2=11+(β2+k2)24k2β2sh⁡2βa\begin{aligned} \mathcal{T}=|S|^2&=\frac{4k^2\beta^2}{\left(\beta^2-k^2\right)\sh^2{\beta a}+4k^2\beta^2\ch^2{\beta a}} \\ &=\frac{4k^2\beta^2}{\left(\beta^2-k^2\right)^2\sh^2{\beta a}+4k^2\beta^2-4k^2\beta^2\sh^2{\beta a}} \\ &\underrightarrow{\ch^2{x}=1-\sh^2{x}} \\ &=\frac{4k^2\beta^2}{\left(\beta^2+k^2\right)^2\sh^2{\beta a}+4k^2\beta^2} \\ &=\frac{1}{1+\frac{\left(\beta^2+k^2\right)^2}{4k^2\beta^2}\sh^2{\beta a}} \end{aligned}

k2=2mEℏ2,β=2m(V0−E)ℏ2k^2=\frac{2mE}{\hbar^2},\quad\beta=\frac{2m(V_0-E)}{\hbar^2}

β2+k2=2mV0ℏ2\beta^2+k^2=\frac{2mV_0}{\hbar^2}

(β2+k2)24k2β2=4m2V02ℏ442mℏ2E⋅2mℏ2(V0−E)=4m2V02ℏ416m2E(V0−E)ℏ4=V024E(V0−E)=14EV0(1−EV0)\frac{\left(\beta^2+k^2\right)^2}{4k^2\beta^2}=\frac{\frac{4m^2V_0^2}{\hbar^4}}{4\frac{2m}{\hbar^2}E\cdot\frac{2m}{\hbar^2}(V_0-E)}=\frac{\frac{4m^2V_0^2}{\hbar^4}}{\frac{16m^2E(V_0-E)}{\hbar^4}}=\frac{V_0^2}{4E(V_0-E)}=\frac{1}{4\frac{E}{V_0}\left(1-\frac{E}{V_0}\right)}

T=∣S∣2=[1+sh⁡2βa4EV0(1−EV0)]−1\mathcal{T}=|\mathrm{S}|^2=\left[1+\frac{\sh^{2}{\beta a}}{4\frac{E}{V_0}\left(1-\frac{E}{V_0}\right)}\right]^{-1}

R=∣R∣2=(k2+β2)sh⁡2βa(k2+β2)sh⁡2βa+4k2β2\mathcal{R}=|\mathrm{R}|^2=\frac{\left(k^2+\beta^2\right)\sh^2{\beta a}}{\left(k^2+\beta^2\right)\sh^2{\beta a}+4k^2\beta^2}

R+T=∣R∣2+∣S∣2=1\mathcal{R}+\mathcal{T}=|\mathrm{R}|^2+|\mathrm{S}|^2=1

当 βa≫1\beta a\gg 1 ,也就是说 2m(V0−E)a2ℏ2≫1\sqrt{\frac{2m\left(V_0-E\right)a^2}{\hbar^2}}\gg 1 时,近似地有 sh⁡βa=eβa−e−βa2≅12eβa\sh{\beta a}=\frac{e^{\beta a}-e^{-\beta a}}{2}\cong \frac{1}{2}e^{\beta a} ,代入 T\mathcal{T} 后,有

T≅[1+14EV0(1−EV0)⋅14e2βa]−1=[14EV0(1−EV0)⋅14e2βa]−1=16E(V0−E)V02e−2βa=16E(V0−E)V02exp⁡(−2aℏ2m(V0−E))\mathcal{T}\cong\left[1+\frac{1}{4\frac{E}{V_0}\left(1-\frac{E}{V_0}\right)}\cdot\frac{1}{4}e^{2\beta a}\right]^{-1}=\left[\frac{1}{4\frac{E}{V_0}\left(1-\frac{E}{V_0}\right)}\cdot\frac{1}{4}e^{2\beta a}\right]^{-1}=\frac{16E(V_0-E)}{V_0^2}e^{-2\beta a}=\frac{16E(V_0-E)}{V_0^2}\exp\left(-\frac{2a}{\hbar}\sqrt{2m(V_0-E)}\right)

卡西米尔1948年提出,“真空之力”

卡西米尔效应

两块中性金属板之间充满电磁场形成的驻波。

(非严格证明)

E总=12ℏω0(1)+12ℏω0(2)+⋯=∞E_{\text{总}}=\frac{1}{2}\hbar\omega_{0}^{(1)}+\frac{1}{2}\hbar\omega_{0}^{(2)}+\cdots=\infty

左移金属板,ω\omega 改变

E~总=12ℏω~0(1)+12ℏω~0(2)+⋯=∞\tilde{E}_{\text{总}}=\frac{1}{2}\hbar\tilde{\omega}_{0}^{(1)}+\frac{1}{2}\hbar\tilde{\omega}_{0}^{(2)}+\cdots=\infty

E~总−E总=ΔE\tilde{E}_{\text{总}}-E_{\text{总}}=\Delta E

卡西米尔力 FF 为负值,也就是向右的吸引力。

范德华力(范德瓦尔斯力,又称分子间相互作用力)[1]
两个中性原子(或分子):

伦敦指出,两个原子在不断振动,会有一个瞬时的偶极矩。范德华力是卡西米尔效应在分子尺度的宏观表现,卡西米尔效应是范德华力的微观效应。常见是蜘蛛、壁虎爬墙的力。

壁虎
卢卡斯

§ 参考资料


  1. Condensed Matter Field Theory(凝聚态场论) ↩︎

更新于 阅读次数 次

请我喝[茶]~( ̄▽ ̄)~*

星旅人 微信支付

微信支付

星旅人 支付宝

支付宝


在提问之前,你应该学会如何提问