UA OPTI570 量子力学29 摄动理论简介
UA OPTI570 量子力學(xué)29 攝動理論簡介
這一篇介紹一些簡單的perturbation theory基礎(chǔ),攝動理論(perturbation theory)的思想是用微擾方法將復(fù)雜的量子系統(tǒng)用簡單的理想化模型近似,以此獲得復(fù)雜量子系統(tǒng)的一些性質(zhì)。
攝動方法描述 已知系統(tǒng)的特征方程
H0∣?ni?=En0∣?ni?,i∈{1,?,gn}H_0|\phi_n^i \rangle = E_n^0 |\phi_n^i \rangle,i \in \{1,\cdots,g_n\}H0?∣?ni??=En0?∣?ni??,i∈{1,?,gn?}
其中gng_ngn?是第nnn個本征態(tài)的degeneracy index,
- gn=1g_n=1gn?=1: non-degenerate case,特征方程簡寫為H0∣?n?=En0∣?n?H_0|\phi_n \rangle=E_n^0|\phi_n \rangleH0?∣?n??=En0?∣?n??
- gn>1g_n>1gn?>1: degenerate case,能量本征值En0E_n^0En0?對應(yīng)gng_ngn?個本征態(tài)
引入另一個哈密頓量
H(λ)=H0+W=H0+λW^,λ∈R,∣λ∣<<1H(\lambda)=H_0+W = H_0 + \lambda \hat W, \lambda \in \mathbb{R},|\lambda|<<1H(λ)=H0?+W=H0?+λW^,λ∈R,∣λ∣<<1
稱λW^\lambda \hat WλW^是H0H_0H0?的一個perturbation。這個哈密頓量的特征方程為
H(λ)∣ψn,j?=En,j(λ)∣ψn,j?H(\lambda)|\psi_{n,j} \rangle = E_{n,j}(\lambda)|\psi_{n,j} \rangleH(λ)∣ψn,j??=En,j?(λ)∣ψn,j??
假設(shè)這個H(λ)H(\lambda)H(λ)就是我們要研究的更復(fù)雜的量子系統(tǒng)的Hamitonian,那么攝動理論要解決的問題就是怎么找一個合適的H0H_0H0?,能讓我們用哈密頓量為H0H_0H0?的系統(tǒng)近似這個更復(fù)雜的量子系統(tǒng)。
2-D Q.H.O.攝動
考慮用2-D Q.H.O.作為更簡單的量子系統(tǒng),也就是
H0=P22m+12mw02(X2+Y2)H_0=\frac{P^2}{2m}+\frac{1}{2}mw_0^2(X^2+Y^2)H0?=2mP2?+21?mw02?(X2+Y2)
我們給它加上一個perturbation,比如
W=?w(X/σ)4W=\hbar w(X/\sigma)^4W=?w(X/σ)4
其中σ\sigmaσ是quantum length scale,σ=?mw\sigma=\sqrt{\frac{\hbar}{mw}}σ=mw???,w<<w0w<<w_0w<<w0?,引入λ=w/w0\lambda=w/w_0λ=w/w0?,顯然λ<<1\lambda<<1λ<<1,并且
W=λ?w0(X/σ)4?λW^W=\lambda \hbar w_0 (X/\sigma)^4 \triangleq \lambda \hat WW=λ?w0?(X/σ)4?λW^
先考慮比較簡單的non-degenerate case:H=H0+WH=H_0+WH=H0?+W,它的特征方程為
H∣ψn?=En∣ψn?H|\psi_n \rangle = E_n |\psi_n \rangleH∣ψn??=En?∣ψn??
不加證明地給出下面兩個結(jié)果
En≈En0+λ??n∣W^∣?n?+λ2∑p≠n∑i=1gp∣??pi∣W^∣?n?∣2En0?Ep0∣ψn?≈∣?n?+λ∑p≠n∑i=1gp??pi∣W^∣?n?En0?Ep0∣?pi?E_n\approx E_n^0+\lambda \langle \phi_n | \hat W | \phi_n \rangle + \lambda^2 \sum_{p \ne n}\sum_{i=1}^{g_p}\frac{|\langle \phi_p^i|\hat W | \phi_n \rangle|^2}{E_n^0-E_p^0} \\ |\psi_n \rangle \approx |\phi_n \rangle + \lambda \sum_{p \ne n}\sum_{i=1}^{g_p}\frac{\langle \phi_p^i|\hat W|\phi_n \rangle}{E_n^0-E_p^0}|\phi_p^i \rangleEn?≈En0?+λ??n?∣W^∣?n??+λ2p?=n∑?i=1∑gp??En0??Ep0?∣??pi?∣W^∣?n??∣2?∣ψn??≈∣?n??+λp?=n∑?i=1∑gp??En0??Ep0???pi?∣W^∣?n???∣?pi??
例:non-generate stationary perturbation的簡單算例
假設(shè)
?n(x)=2Lsin?(nπx/L)En0=n2?2π22mL2=n2A,n=1,2,3,?\phi_n(x)=\sqrt{\frac{2}{L}}\sin(n\pi x/L) \\ E_n^0 =\frac{n^2 \hbar^2 \pi^2}{2mL^2}=n^2A,n=1,2,3,\cdots?n?(x)=L2??sin(nπx/L)En0?=2mL2n2?2π2?=n2A,n=1,2,3,?
其中AAA表示基態(tài)能量本征值E10E_1^0E10?,perturbation為W=?X/L,?<<AW=\epsilon X/L,\epsilon<<AW=?X/L,?<<A,引入λ=?/A\lambda=\epsilon/Aλ=?/A,則W=λA^=λ(AX/L)W=\lambda \hat A=\lambda(AX/L)W=λA^=λ(AX/L)
問題1:WWW如何改變基態(tài)能量本征值?
E1≈E10+λ??1∣W^∣?1?+λ2∑p≠n∣??p∣W^∣?1?∣2E10?Ep0=A+λA??1∣X/L∣?1?+λ2A2∑p=2+∞∣??p∣X/L∣?1?∣2A?p2A\begin{aligned}E_1 & \approx E_1^0+\lambda \langle \phi_1|\hat W |\phi_1 \rangle+\lambda^2 \sum_{p \ne n} \frac{|\langle \phi_p|\hat W|\phi_1 \rangle|^2}{E_1^0-E_p^0} \\ & = A+\lambda A \langle \phi_1 |X/L|\phi_1 \rangle +\lambda^2A^2 \sum_{p=2}^{+\infty} \frac{|\langle \phi_p|X/L|\phi_1 \rangle|^2}{A-p^2A} \end{aligned}E1??≈E10?+λ??1?∣W^∣?1??+λ2p?=n∑?E10??Ep0?∣??p?∣W^∣?1??∣2?=A+λA??1?∣X/L∣?1??+λ2A2p=2∑+∞?A?p2A∣??p?∣X/L∣?1??∣2??
其中
??p∣X/L∣?1?=∫0L2Lsin?(πx/L)sin?(pπx/L)xLdx={1/2,p=1?8pπ2(p2?1)2,p≥2,peven0,p≥3,podd\begin{aligned} \langle \phi_p|X/L|\phi_1 \rangle &= \int_0^L \frac{2}{L}\sin(\pi x/L)\sin(p \pi x/L) \frac{x}{L}dx \\ & = \begin{cases}1/2,p=1 \\ \frac{-8p}{\pi^2(p^2-1)^2},p \ge 2,p\ even \\ 0,p \ge 3,p \ odd\end{cases}\end{aligned}??p?∣X/L∣?1???=∫0L?L2?sin(πx/L)sin(pπx/L)Lx?dx=??????1/2,p=1π2(p2?1)2?8p?,p≥2,p?even0,p≥3,p?odd??
綜上可得,
E1≈A+A2λ?λ2A∑p=2,peven1p2?164p2π4(p2?1)4E_1 \approx A +\frac{A}{2}\lambda - \lambda^2 A \sum_{p=2,p\ even} \frac{1}{p^2-1}\frac{64p^2}{\pi^4(p^2-1)^4}E1?≈A+2A?λ?λ2Ap=2,p?even∑?p2?11?π4(p2?1)464p2?
問題2:將攝動后的基態(tài)寫為原來的本征態(tài)的疊加。
∣ψ1?≈∣?1?+λA∑p=2+∞??p∣X/L∣?1?A?Ap2∣?p?=∣?1?+8λπ2(227∣?2?+43375∣?3?+642875∣?4?+?)\begin{aligned} |\psi_1 \rangle & \approx |\phi_1 \rangle + \lambda A \sum_{p=2}^{+\infty}\frac{\langle \phi_p|X/L|\phi_1 \rangle}{A-Ap^2}|\phi_p \rangle \\ & =|\phi_1 \rangle + \frac{8 \lambda}{\pi^2} \left( \frac{2}{27}|\phi_2 \rangle + \frac{4}{3375}|\phi_3 \rangle + \frac{6}{42875}|\phi_4 \rangle + \cdots \right) \end{aligned}∣ψ1???≈∣?1??+λAp=2∑+∞?A?Ap2??p?∣X/L∣?1???∣?p??=∣?1??+π28λ?(272?∣?2??+33754?∣?3??+428756?∣?4??+?)?
可以發(fā)現(xiàn)攝動雖然會導(dǎo)致基本偏移,但與更高階的能量本征態(tài)之間的線性系數(shù)下降得非常快,當(dāng)λ\lambdaλ本身就很小時,可以粗略地認(rèn)為∣ψ1?=∣?1?|\psi_1 \rangle = |\phi_1 \rangle∣ψ1??=∣?1??,或者只用基態(tài)和下一級本征態(tài)的疊加:
∣ψ1?≈∣?1?+16λ27π2∣?2?|\psi_1 \rangle \approx |\phi_1 \rangle + \frac{16 \lambda}{27 \pi^2}|\phi_2 \rangle∣ψ1??≈∣?1??+27π216λ?∣?2??
總結(jié)
以上是生活随笔為你收集整理的UA OPTI570 量子力学29 摄动理论简介的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA OPTI501 电磁波 Loren
- 下一篇: UA OPTI501 电磁波 Loren