Quantom Mechanics Week 4

QM Week 4

2.20

a

∫ψ∗ψdx=∫A2e−2a∣x∣dx=1  ⟹  A=a\int \psi^* \psi dx = \int A^2 e^{-2a|x|} dx = 1\implies A = \sqrt a

Ψ(x,0)=ae−a∣x∣\Psi(x,0) = \sqrt a e^{-a|x|}

b

ϕ(k)=12π∫Ψ(x,0)e−ikx dx=2aa2π(a2+k2)\phi(k) = \frac{1}{\sqrt{2\pi}} \int \Psi(x,0) e^{-ikx}\, dx= \frac{2a\sqrt a }{\sqrt{2\pi}(a^2+k^2)}

c

Ψ(x,t)=12π∫ϕ(k)ei(kx−ℏk22mt) dk=aaπ∫1a2+k2ei(kx−ℏk22mt) dk\begin{aligned} \Psi(x,t) &= \frac{1}{\sqrt{2\pi}}\int \phi(k)e^{i(kx - \frac{\hbar k^2}{2m}t)} \, dk \\ & = \frac{a\sqrt a }{\pi} \int\frac{1}{a^2+k^2} e^{i(kx-\frac{\hbar k^2}{2m}t)} \, dk \\ \end{aligned}

d

lim⁡a→∞ϕ(k)=2πa≃0\lim_{a\to \infty} \phi(k) = \sqrt\frac{2}{\pi a}\simeq 0

lim⁡a→0ϕ(k)=∞\lim_{a\to0}\phi(k) = \infty

2.21

a

1=∫Ψ∗Ψdx=A2∫e−2ax2 dx=A2π2a  ⟹  A=(2aπ)1/41 = \int \Psi^* \Psi dx = A^2 \int e^{-2ax^2} \, dx = A^2\sqrt\frac{\pi}{2a}\\ \implies A = \left(\frac{2a}{\pi}\right)^{1/4}

b

ϕ(k)=12π∫Ψ(x,0)e−ikx dx=12π∫Ae−a∣x∣e−ikx dx\phi(k) = \frac{1}{\sqrt{2\pi}}\int \Psi(x,0) e^{-ikx}\, dx= \frac{1}{\sqrt{2\pi}}\int A e^{-a|x|} e^{-ikx} \, dx\\

Ψ(x,t)=12π∫ϕ(k)ei(kx−ℏk22mt) dk=12π∫Ae−ax2ei(kx−ℏk22mt) dk=Ae−ax22π∫ei(kx−ℏk22mt) dk=\begin{aligned} \Psi(x,t) &= \frac{1}{\sqrt{2\pi}}\int \phi(k)e^{i(kx - \frac{\hbar k^2}{2m}t)} \, dk \\ & = \frac{1}{\sqrt{2\pi}} \int Ae^{-ax^2} e^{i(kx-\frac{\hbar k^2}{2m}t)} \, dk \\ & = \frac{Ae^{-ax^2}}{\sqrt{2\pi}} \int e^{i(kx-\frac{\hbar k^2}{2m}t)} \, dk \\ & = \end{aligned}

拼尽全力无法积分

c

if

Ψ(x,t)=(2aπ)1/41γe−ax2/γ2γ=1+(2iℏat/m)\Psi(x,t) = \left(\frac{2a}\pi\right)^{1/4}\frac{1}\gamma e^{-ax^2/\gamma^2}\qquad \gamma = \sqrt{1+(2i\hbar a t/m)}

∣Ψ∣2=2πωe−2ω2x2|\Psi|^2 = \sqrt \frac{2}{\pi}\omega e^{-2\omega^2 x ^2}

d

⟨x⟩=∫x∣Ψ∣2dx=0\langle x \rangle = \int x|\Psi|^2 dx = 0

⟨p⟩=md⟨x⟩dt=0\langle p \rangle = m \frac{d\langle x \rangle}{dt} = 0

⟨x2⟩=∫x2∣Ψ∣2dx=14ω2\langle x^2 \rangle = \int x^2|\Psi|^2 dx = \frac{1}{4\omega^2}

⟨p2⟩=aℏ2\langle p^2 \rangle = a\hbar^2

σx=⟨x2⟩−⟨x⟩2=12ω\sigma_x = \sqrt{\langle x^2 \rangle - \langle x \rangle^2} = \frac{1}{2\omega}

σp=⟨p2⟩−⟨p⟩2=aℏ2\sigma_p = \sqrt{\langle p^2 \rangle - \langle p \rangle^2} = \sqrt{a\hbar^2}

σxσp=aℏ2ω=ℏ1+(2ℏat/m)22≥ℏ2\sigma_x\sigma_p = \frac{\sqrt a\hbar}{2\omega} = \frac{\hbar \sqrt{1+(2\hbar a t /m )^2}}{2}\geq \frac \hbar 2

2.24

⟨x⟩=⟨p⟩=0\langle x \rangle = \langle p \rangle = 0

⟨x2⟩=2∫mαℏ2x2e−2mα/ℏ2xdx=ℏ42(mα)2\langle x^2 \rangle = 2 \int \frac{m\alpha}{\hbar^2} x^2 e^{-2m\alpha/\hbar^2 x} dx = \frac{\hbar^4}{2(m\alpha)^2}

σxσp=⟨x2⟩⟨p2⟩=2ℏ2≥ℏ2\sigma_x\sigma_p = \sqrt{\langle x^2 \rangle\langle p^2\rangle} = \frac{\sqrt 2\hbar}2\geq \frac \hbar 2

2.25

ψ1(x)=mαℏe−mα∣x∣/ℏ2=Ce−mα∣x∣/ℏ2\psi_1(x) = \frac{\sqrt{m\alpha}}\hbar e^{-m\alpha |x| /\hbar^2} = C e^{-m\alpha |x| /\hbar^2}

ψ2(x)=Aeikx+Be−ikx={eikx+Re−ikxx≤0Teikxx>0\psi_2(x) =A e^{ikx} + B e^{-ikx} = \begin{cases} e^{ikx} + Re^{-ikx} &x \leq 0 \\T e^{ikx} & x > 0\end{cases}

1C⟨ψ1,ψ2⟩=1C∫ψ1∗ψ2 dx=∫−∞0(eikx+Re−ikx)e−mα∣x∣/ℏ2 dx+∫0∞Teikxe−mα∣x∣/ℏ2 dx=−1ik+mα/ℏ2+Rik−mα/ℏ2+Tik−mα/ℏ21=R+T=0\begin{aligned}\frac{1}C\langle \psi_1 ,\psi_2 \rangle &= \frac{1}C\int \psi_1^* \psi_2 \,dx\\ &=\int_{-\infty}^0 ( e^{ikx} + Re^{-ikx}) e^{-m\alpha |x| /\hbar^2}\, dx + \int_0^\infty Te^{ikx}e^{-m\alpha |x| /\hbar^2}\, dx\\ &=\frac{-1}{ik + m\alpha /\hbar^2} + \frac{R}{ik - m\alpha /\hbar^2} + \frac{T}{ik - m\alpha /\hbar^2}\\ 1= R+T\quad &\\ & = 0 \end{aligned}

2.27

a

Case 1, where the wave function is even:

ψ(x)={Bekxx<−aC(ekx+e−kx)−a<x<aBe−kxx>a\psi(x) = \begin{cases} Be^{kx} & x < -a\\ C (e^{kx} + e^{-kx}) & -a<x<a\\ B e^{-kx} & x > a\\ \end{cases}

Boundary condition (continuity):

Be−ak=C(e−ak+eak)Be^{-ak} = C(e^{-ak}+e^{ak})

Schorodinger Equation:

−ℏ22md2ψdx2+α[δ(x+a)+δ(x−a)]ψ=Eψ-\frac{\hbar^2}{2m} \frac{d^2\psi}{dx^2} + \alpha[\delta(x+a)+\delta(x-a)]\psi = E\psi

take the integral on [−a−ε,−a+ε][-a-\varepsilon,-a+\varepsilon] and [a−ε,a+ε][a-\varepsilon,a+\varepsilon] , and let ε→0\varepsilon \to 0 :

−ℏ22mdψdx+α∫δ(x+a)dx+α∫δ(x−a)dx=E∫ψdx-\frac{\hbar^2}{2m} \frac{d\psi}{dx} + \alpha\int \delta(x+a) dx + \alpha \int \delta(x-a) dx = E\int \psi dx

Δ(dψdx)∣−a=2mℏ2lim⁡ε→0∫−a−ε−a+εαδ(x+a)ψ(x) dx=−2mαℏ2ψ(−a)\Delta \left(\frac{d\psi}{dx}\right)\Bigg|_{-a} = \frac{2m}{\hbar^2} \lim_{\varepsilon\to 0} \int_{-a-\varepsilon}^{-a+\varepsilon} \alpha \delta(x+a)\psi(x)\, dx = -\frac{2m\alpha}{\hbar^2}\psi(-a)

The same applies at aa :

Δ(dψdx)∣a=−2mαℏ2ψ(a)\Delta \left(\frac{d\psi}{dx}\right)\Bigg|_{a} = -\frac{2m\alpha}{\hbar^2}\psi(a)

According to the form of ψ\psi :

dψdx={Bkekxx<−aC(kekx−ke−kx)−a<x<a−Bke−kxx>a\frac{d\psi}{dx} = \begin{cases} Bke^{kx} & x < -a\\ C(ke^{kx} -ke^{-kx})& -a<x<a\\ -Bke^{-kx} & x > a\\ \end{cases}

dψdx∣a−=C(keak−ke−ak)dψdx∣a+=−Bke−ak∴Δ(dψdx)∣a=−Bke−ak−C(keak−ke−ak)\begin{aligned} \frac{d\psi}{dx}\Bigg|_{a^-} &= C(ke^{ak} -ke^{-ak})\\ \frac{d\psi}{dx}\Bigg|_{a^+} &= -Bke^{-ak}\\ \therefore \Delta \left(\frac{d\psi}{dx}\right)\Bigg|_{a} &= -Bke^{-ak} - C(ke^{ak} -ke^{-ak}) \end{aligned}

The same appies at −a-a :

Δ(dψdx)∣−a=−Bke−ak+C(ke−ak−keak)\Delta \left(\frac{d\psi}{dx}\right)\Bigg|_{-a} = -Bke^{-ak} + C(ke^{-ak} -ke^{ak})\\

Solve the equation:

{−2mαℏ2ψ(a)=−Bke−ak−C(keak−ke−ak)−2mαℏ2ψ(−a)=−Bke−ak+C(ke−ak−keak)\begin{cases}-\frac{2m\alpha}{\hbar^2}\psi(a) = -Bke^{-ak} - C(ke^{ak} -ke^{-ak})\\ -\frac{2m\alpha}{\hbar^2}\psi(-a) = -Bke^{-ak} + C(ke^{-ak} -ke^{ak}) \end{cases}

  ⟹  −2mαℏ2C(eak+e−ak)=−2Ckeakmαℏ2(e2ak+1)=ke2ak#1\begin{aligned}\implies \frac{-2m\alpha}{\hbar^2} C(e^{ak}+e^{-ak}) &= -2Cke^{ak}\\ \frac{m\alpha}{\hbar^2}( e^{2ak} +1 ) &= k e^{2ak}\qquad \# 1 \end{aligned}

E=−ℏ2k22mE = -\frac{\hbar^2 k^2}{2m}

Case 2, where the wave function is odd:

ψ(x)={Bekxx<−aC(ekx−e−kx)−a<x<a−Be−kxx>a\psi(x) = \begin{cases} Be^{kx} & x < -a\\ C (e^{kx} - e^{-kx}) & -a<x<a\\ -B e^{-kx} & x > a\\ \end{cases}

Boundary condition (continuity):

Be−ak=C(e−ak−eak)Be^{-ak} = C(e^{-ak}-e^{ak})

Repeat the steps in case 1, we get:

mαℏ2(e2ak−1)=ke2ak#2\frac{m\alpha}{\hbar^2}( e^{2ak} -1 ) = k e^{2ak}\qquad \# 2

E=−ℏ2k22mE = -\frac{\hbar^2 k^2}{2m}

The number of bound states are the number of solutions to equations #1 and #2 accordingly.

b

  1. When α=ℏ2/ma\alpha = \hbar^2/ma :
    For case 1:

ak=1+e−2ak,ak≃1.109ak = 1+ e^{-2ak},\quad ak\simeq 1.109

For case 2:

ak=1−e−2ak,ak≃0.797ak = 1- e^{-2ak},\quad ak \simeq 0.797

  1. When α=ℏ2/4ma\alpha = \hbar^2/4ma :
    For case 1:

2ak=1+e−2ak,ak≃0.6402ak = 1+ e^{-2ak},\quad ak\simeq 0.640

For case 2:

2ak=1−e−2ak,ak=0not possible2ak = 1- e^{-2ak},\quad ak = 0\quad \text{not possible}

The according energy states can be caculated with E=−ℏ2k2/2mE = -\hbar^2k^2/2m .

c

  1. When a→0a\to 0 , there could only be one bound state, which takes the even wave function.

1+e−2ak=ℏ2kmα  ⟹  lim⁡a→0k=2mαℏ2  ⟹  E=−2mα2ℏ2 1 + e^{-2ak} = \frac{\hbar^2k}{m\alpha}\\ \implies \lim_{a\to 0} k = \frac{2m\alpha}{\hbar^2}\implies E = -\frac{2m\alpha^2}{\hbar^2}

The energy is the same to where V(x)=−2αδ(x)V(x) =- 2\alpha\delta(x) , as if the two heaps of the delta functions are now merged into one.

  1. When a→∞a\to \infty, there are two bound states.
    For case 1:

1+e−2ak=ℏ2kmα  ⟹  lim⁡a→∞k=mαℏ2  ⟹  E=−mα22ℏ21 +e^{-2ak} = \frac{\hbar^2k}{m\alpha}\\ \implies \lim_{a\to \infty} k = \frac{m\alpha}{\hbar^2}\implies E = -\frac{m\alpha^2}{2\hbar^2}

For case 2:

1−e−2ak=ℏ2kmα  ⟹  lim⁡a→∞k=mαℏ2  ⟹  E=−mα22ℏ21- e^{-2ak} = \frac{\hbar^2k}{m\alpha}\\ \implies \lim_{a\to \infty} k = \frac{m\alpha}{\hbar^2}\implies E = -\frac{m\alpha^2}{2\hbar^2}

2.28

Consider the particle traveling from left to right, we have the wave function:

ψ(x)={Aeikx+Be−ikxx<−aCeikx+De−ikx−a<x<aFeikxx>a\psi(x) = \begin{cases} Ae^{ikx} + Be^{-ikx} & x < -a\\ C e^{ikx} + De^{-ikx} & -a<x<a\\ F e^{ikx} & x > a \end{cases}

Boundary conditions(continuity):

Ae−iak+Beiak=Ce−iak+DeiakCeiak+De−iak=Feiak\begin{aligned} Ae^{-iak} + Be^{iak} &= C e^{-iak} + De^{iak}\\ Ce^{iak} + De^{-iak} &= Fe^{iak}\end{aligned}

let β=e−2iak\beta = e^{-2iak}

βA+B=βC+DC+βD=F\begin{align}\beta A+B &= \beta C+D\\C+\beta D &= F\end{align}

Δψ′∣−a=−mαℏ2ψ(−a)\Delta \psi'|_{-a} = - \frac{m\alpha}{\hbar^2}\psi(-a) :

(Cik−Aik)e−ika−(Dik−Bik)eika=−mαℏ2(Ae−ika+Beika)(Cik - Aik)e^{-ika} - (Dik -Bik)e^{ika} = -\frac{m\alpha}{\hbar^2}(Ae^{-ika}+Be^{ika})

β(C−A)−D−B=imαkℏ2(βA+B)\begin{align} \beta (C-A) - D -B& = \frac{im\alpha}{k\hbar^2}(\beta A+ B)\\\end{align}

Δψ′∣a=−mαℏ2ψ(a)\Delta \psi'|_{a} = - \frac{m\alpha}{\hbar^2}\psi(a) :

(F−C)ikeika+Dike−ika=−mαℏ2Fika(F-C)ike^{ika} + Dike^{-ika} = -\frac{m\alpha}{\hbar^2} F^{ika}

F−C+βD=imαkℏ2F\begin{align}F-C +\beta D = \frac{im\alpha}{k\hbar^2} F\end{align}

let γ=imαkℏ2\gamma = \frac{im\alpha}{k\hbar^2}

According to equation (1) to (4) and assume A=1A = 1 for computational ease:

[1−β−1001β−1−1−γβ−100−1β1−γ][BCDF]=[−β0(γ+1)β0]\begin{bmatrix} 1 &-\beta &-1 &0\\ 0 &1 &\beta &-1\\ -1-\gamma & \beta &-1 &0\\ 0 &-1 &\beta &1-\gamma \end{bmatrix} \begin{bmatrix} B\\ C\\ D\\ F\end{bmatrix} =\begin{bmatrix} -\beta\\ 0\\ (\gamma+1)\beta\\ 0 \end{bmatrix}

T=∣F∣2∣A∣2=∣F∣2T = \frac{|F|^2}{|A|^2} = |F|^2

这周作业是我抄答案都懒得抄完的level. 最后两题对着答案写的都写不明白 , 用 mathematica 倒腾半天才算出来. 🥲🥲🥲


Quantom Mechanics Week 4
https://blog.jacklit.com/2025/03/17/qm_hw4/
作者
Jack H
发布于
2025年3月17日
许可协议