本文链接：个人站 | 简书 | CSDN
版权声明：除特别声明外，本博客文章均采用 BY-NC-SA 许可协议。转载请注明出处。
前几天整理电脑的时候发现了本科上量子力学讨论班时做的一个 Slide，觉得挺有意思的。花了点时间整理成这篇博客。
一维谐振子
一个质量为 m m m 的粒子，在一维势场 V ( x ) = 1 2 m ω 2 x 2 V(x) = \dfrac12m\omega^2x^2 V(x)=21mω2x2 中运动。其哈密顿算符为
H ^ = p ^ 2 2 m + 1 2 m ω 2 x ^ 2 \hat H = \frac{\hat p^2}{2m} + \frac12m\omega^2\hat x^2 H^=2mp^2+21mω2x^2
其中 x ^ \hat x x^ 为位置算符， p ^ = − i ℏ d d x \hat p = -i\hbar\dfrac{\mathrm d}{\mathrm dx} p^=−iℏdxd 为动量算符。我们需要求解该体系的定态 Schrödinger 方程：
H ^ ∣ ψ > = E ∣ ψ > \hat H\left|\psi\right> = E\left|\psi\right> H^∣ψ⟩=E∣ψ⟩
一维谐振子是除了氢原子之外，为数不多的可以解析求解的体系。那么我们为什么要费劲求它的数值解呢？正因为绝大多数的量子体系都无法解析求解，数值方法才显得尤为重要。
有限差分法
回忆一下泰勒公式
f ( a + h ) = f ( a ) + f ′ ( a ) 1 ! h + f ′ ′ ( a ) 2 ! h 2 + o ( h 3 ) f(a+h) = f(a) + \frac{f'(a)}{1!}h+\frac{f''(a)}{2!}h^2+o(h^3) f(a+h)=f(a)+1!f′(a)h+2!f′′(a)h2+o(h3)
令 h = − h h=-h h=−h，有
f ( a − h ) = f ( a ) − f ′ ( a ) 1 ! h + f ′ ′ ( a ) 2 ! h 2 + o ( h 3 ) f(a-h) = f(a) - \frac{f'(a)}{1!}h+\frac{f''(a)}{2!}h^2+o(h^3) f(a−h)=f(a)−1!f′(a)h+2!f′′(a)h2+o(h3)
两式相加，可得
f ′ ′ ( a ) = f ( a − h ) + f ( a + h ) − 2 f ( a ) h 2 + o ( h 3 ) ≈ f ( a − h ) + f ( a + h ) − 2 f ( a ) h 2 \begin{aligned} f''(a) &= \frac{f(a-h)+f(a+h)-2f(a)}{h^2} + o(h^3)\\ &\approx \frac{f(a-h)+f(a+h)-2f(a)}{h^2} \end{aligned} f′′(a)=h2f(a−h)+f(a+h)−2f(a)+o(h3)≈h2f(a−h)+f(a+h)−2f(a)
将 ψ ( x ) \psi(x) ψ(x) 在 x ∈ [ − r , r ] x\in[-r, r] x∈[−r,r] 区间离散化为
ϕ i ≡ ψ ( x i ) = ψ ( i Δ x − r ) , i = 0 , 1 , 2 , ⋯ , N \phi_i \equiv \psi(x_i) = \psi(i\Delta x - r),\quad i=0, 1, 2, \cdots, N ϕi≡ψ(xi)=ψ(iΔx−r),i=0,1,2,⋯,N
其中 N = 2 r / Δ x N=2r/\Delta x N=2r/Δx，则 Schrödinger 方程差分化为
− ℏ 2 2 m ϕ i − 1 + ϕ i + 1 − 2 ϕ i Δ x 2 + 1 2 m ω 2 x i 2 ϕ i = E ϕ i -\frac{\hbar^2}{2m}\frac{\phi_{i-1}+\phi_{i+1}-2\phi_i}{\Delta x^2}+\frac12m\omega^2x_i^2\phi_i = E\phi_i −2mℏ2Δx2ϕi−1+ϕi+1−2ϕi+21mω2
最低0.47元/天 解锁文章
4054
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
