Time-Dependent-Perturbation

• Time-dependent Problem
  • Interaction Picture
  • Example: Nuclear Magnetic Resonance
• Time-dependet Perturbation
  • Dyson Series
  • Transition Probability
  • Constant Perturbation
    • Zero Order
    • First Order
    • Second Order
  • Harmonic Perturbation
  • Expotential Decay Perturbation-Energy Shift and Decay Width
    • First Order
    • Second Order
  • Lamb Shift

Time-dependent Problem

从一阶微扰可以知道：

$$|n^1\rangle=\sum_k c_k|k\rangle$$

随时间的演化：

$$|n^1(t)\rangle=\sum_k c_ke^{-\frac{iE_kt}{\hbar}}|k\rangle$$

当我们考虑和哈密顿量和时间相关的时候

$$|n^1(t)\rangle=\sum_k c_k(t)e^{-\frac{iE_kt}{\hbar}}|k\rangle$$

我们的问题转化为求$c_k(t)$满足的微分方程。在求解之前，不妨先采用相互作用绘景。

Interaction Picture

从相互作用绘景可以精确（而不是使用微扰）求解含时的演化问题，所以我们用任意态矢量代替上文的一阶态矢量修正：

$$|\alpha(t)\rangle=\sum_n c_n(t)e^{-\frac{iE_nt}{\hbar}}|n\rangle$$

定义相互作用绘景下的态矢量：

$$|\alpha_I(t)\rangle=e^{\frac{i\hat H_0t}{\hbar}}|\alpha_S(t)\rangle$$

角标$S$表示薛定谔绘景，$I$表示相互作用绘景。相互作用绘景的具体内容已经在 对称性和守恒律 中介绍过了。从相互作用绘景的薛定谔方程出发：

$$\hat{H}=\hat{H}_0+\hat{V}_I(t),\quad i\hbar\frac{d}{dt}|\alpha_I(t)\rangle=\hat V_I(t)|\alpha_I(t)\rangle$$

所以

$$i\hbar\frac{d}{dt}\langle n|\alpha_I(t)\rangle=\sum_m\langle n|\hat V_I(t)|m\rangle\langle m|\alpha_I(t)\rangle$$

其中：

$$\langle n|\hat V_I(t)|m\rangle=\langle n|e^{\frac{i\hat H_0t}{\hbar}}\hat V_I e^{-\frac{i\hat H_0t}{\hbar}}|m\rangle=\langle n|\hat V_I|m\rangle e^{i(E_n-E_m)t/\hbar}$$

也就是说：

$$i\hbar\frac{d}{dt}c_n(t)=\sum_mV_{nm} c_m(t)e^{i(E_n-E_m)t/\hbar}$$

用矩阵表示：

$$i\hbar\frac{d}{dt}\begin{pmatrix}c_1(t)\\c_2(t)\\\vdots\end{pmatrix}=\begin{pmatrix}V_{11}&V_{12}e^{i\omega_{12}t}&\cdots\\V_{21}e^{i\omega_{21}t}&V_{22}&\cdots\\\vdots&\vdots&\ddots\end{pmatrix}\begin{pmatrix}c_1(t)\\c_2(t)\\\vdots\end{pmatrix}$$

其中

$$\omega_{nm}=\frac{E_n-E_m}{\hbar}=-\omega_{mn}$$

这个和 Density-Matrix 中密度矩阵的演化很像，不过密度矩阵的演化是：

$$\frac{d}{dt}\rho_{mn}=\frac{1}{i\hbar}(E_m-E_n)\rho_{mn}\Rightarrow \rho_{mn}(t)=\rho_{mn}(0)e^{-i\omega_{mn}t}$$

多了个负号。

Example: Nuclear Magnetic Resonance

一个最简单的相互作用绘景应用例子是核磁共振，一个二能级体系：

$$\hat H_0=\begin{pmatrix}E_1&0\\0&E_2\end{pmatrix}$$ $$\hat V=\begin{pmatrix}0&Ve^{i\omega t}\\Ve^{-i\omega t}&0\end{pmatrix}$$

初始条件：$c_1(0)=1,c_2(0)=0$，求解$c_1(t),c_2(t)$。

矩阵形式：

$$i\hbar\frac{d}{dt}\begin{pmatrix}c_1(t)\\c_2(t)\end{pmatrix}=\begin{pmatrix}0&Ve^{i(\omega+\omega_{12}) t}\\Ve^{-i(\omega+\omega_{12})t}&0\end{pmatrix}\begin{pmatrix}c_1(t)\\c_2(t)\end{pmatrix}$$

记

$$\Omega=\sqrt{\frac{V^2}{\hbar^2}+\frac{(\omega+\omega_{12})^2}{4}}$$

解得：

$$c_2(t)=-\frac{\gamma}{2\Omega\hbar}e^{i\Omega t}e^{i\frac{\omega+\omega_{12}}{2}t}+\frac{\gamma}{2\Omega\hbar}e^{-i\Omega t}e^{i\frac{\omega+\omega_{12}}{2}t}=\frac{i\gamma}{\Omega\hbar}\sin{[\Omega t]}e^{i\frac{\omega+\omega_{12}}{2}t}$$ $$|c_2(t)|^2=\frac{\gamma^2}{\hbar^2\Omega^2}\sin^2{[\Omega t]}$$ $$|c_1(t)|^2=1-|c_2(t)|^2=1-\frac{\gamma^2}{\hbar^2\Omega^2}\sin^2{[\Omega t]}$$

这就是Rabi振荡。

Time-dependet Perturbation

Dyson Series

处理含时微扰的方法与上面处理非含时微扰的方法有所不同。我们可以通过Dyson级数来求解。

假设系数的近似解为：

$$c_n(t)=c_n^{(0)}(t)+c_n^{(1)}(t)+c_n^{(2)}(t)+\cdots$$

知道$c_n^{(0)}(t)=c_n(0)$，是容易求得$c_n^{(1)}(t)$的，将其作为初值重新代入，我们可以不断迭代得到更高阶的近似解。

相互作用中关于时间演化算符的微分方程为：

$$i\hbar\frac{d}{dt}U_I(t,t_0)=\hat V_I(t)U_I(t,t_0)$$

初始条件：

$$U_I(t_0,t_0)=1$$

得到积分方程：

$$U_I(t,t_0)=1-\frac{i}{\hbar}\int_{t_0}^t\hat V_I(t')U_I(t',t_0)dt'$$

继续迭代：

$$\begin{aligned} U_I(t,t_0)&=1-\frac{i}{\hbar}\int_{t_0}^t\hat V_I(t')[1-\frac{i}{\hbar}\int_{t_0}^{t'}\hat V_I(t'')U_I(t'',t_0)dt'']dt'\\ &=1-\frac{i}{\hbar}\int_{t_0}^t\hat V_I(t')dt'-\frac{1}{\hbar^2}\int_{t_0}^t\int_{t_0}^{t'}\hat V_I(t')\hat V_I(t'')U_I(t'',t_0)dt''dt' \end{aligned}$$

这其中的第n阶项为：

$$U_I(t,t_0) =\sum_{n=0}^\infty(-\frac{i}{\hbar})^n\int_{t_0}^t\cdots\int_{t_0}^{t_{n-1}}\hat V_I(t_1)\cdots\hat V_I(t_n)dt_1\cdots dt_n$$

Transition Probability

考虑跃迁的概率：

$$c_n(t)=\langle n|U_I(t,t_0)|i\rangle$$

把$U_I(t,t_0)$代入：

$$c_n(t)=c_n^{(0)}(t)+c_n^{(1)}(t)+c_n^{(2)}(t)+\cdots$$

其中：

$$c_n^{(0)}(t)=\langle n|i\rangle$$ $$c_n^{(1)}(t)=\frac{-i}{\hbar}\int_{t_0}^t\langle n|\hat V_I(t')|i\rangle dt'=\frac{-i}{\hbar}\int_{t_0}^te^{i\omega_{ni}t'}V_{ni}(t')dt'$$ $$c_n^{(2)}(t)=\frac{-1}{\hbar^2}\sum_m\int_{t_0}^t\int_{t_0}^{t'}\langle n|\hat V_I(t')|m\rangle\langle m|\hat V_I(t'')|i\rangle dt''dt'$$

所以跃迁概率为：

$$P_{ni}(t)=|\sum_{n=0}^\infty c_n(t)|^2$$

Constant Perturbation

考虑一个常数微扰：

$$V(t)=\begin{cases}0&t<0\\V&t\geq 0\end{cases}$$

假定在$t=0$的时候，体系处于态$|i\rangle$：

Zero Order

$$c_n^{(0)}(t)=\langle n|i\rangle$$

First Order

$$\begin{aligned}c_n^{(1)}(t)&=\frac{-i}{\hbar}\int_{0}^t\langle n|\hat V_I(t')|i\rangle dt'\\&=\frac{-i}{\hbar}\int_{0}^te^{i\omega_{ni}t'}V_{ni}dt'\\&=\frac{V_{ni}}{E_n-E_i}(1-e^{i\omega_{ni}t})\end{aligned}$$

或者：

$$|c_n^{(1)}(t)|^2=\frac{4|V_{ni}|^2}{(E_n-E_i)^2}\sin^2[\frac{\omega_{ni}t}{2}]$$

对这个式子可以进行更多的讨论：

$$|c_n^{(1)}(t)|^2=\frac{4|V_{ni}|^2}{(E_n-E_i)^2}\sin^2[\frac{\omega_{ni}t}{2}]=\frac{|V_{ni}|^2}{\hbar^2}t^2\frac{\sin^2[\frac{\omega_{ni}t}{2}]}{(\frac{\omega_{ni}t}{2})^2}$$

这意味着，当微扰刚刚打开的时候，能量变化可以很大，跃迁的选择有很多；但微扰打开的时候长了以后，能量变化的区间很窄。

当末态具有不同而相近的能级时，我们可以计算总的跃迁概率：

$$\sum_n|c_n^{(1)}(t)|^2=\int dE \rho(E)\frac{|V_{ni}|^2}{\hbar^2}t^2\frac{\sin^2[\frac{\omega_{ni}t}{2}]}{(\frac{\omega_{ni}t}{2})^2}=\frac{2\pi t}{\hbar}\overline{|V_{ni}^2}|\rho(E)|_{E_n=E_i},t\rightarrow \infty$$

所以对于长时间，总的跃迁概率正比于时间。

定义单位时间的跃迁概率——跃迁速率：

$$W_{ni}=\frac{2\pi}{\hbar}\overline{|V_{ni}^2|}\rho(E)|_{E_n=E_i}$$

该式被称为费米黄金规则。

Second Order

$$\begin{aligned}c_n^{(2)}(t)&=\frac{-1}{\hbar^2}\sum_m\int_{0}^t\int_{0}^{t'}\langle n|\hat V_I(t')|m\rangle\langle m|\hat V_I(t'')|i\rangle dt''dt'\\ &=\frac{-1}{\hbar^2}\sum_m\int_{0}^t\int_{0}^{t'}e^{i\omega_{nm}t'}V_{nm}e^{-i\omega_{mi}t''}V_{mi}dt''dt'\\ &=\frac{-1}{\hbar^2}\sum_mV_{nm}V_{mi}\int_{0}^t\int_{0}^{t'}e^{i\omega_{nm}t'}e^{-i\omega_{mi}t''}dt''dt'\\ &=\frac{-1}{\hbar^2}\sum_mV_{nm}V_{mi}\frac{1}{i\omega_{mi}}\int_{0}^t e^{i\omega_{nm}t'}(e^{i\omega_{mi}t'}-1)dt'\\ &=\frac{i}{\hbar}\sum_m\frac{V_{nm}V_{mi}}{E_m-E_i}\int_{0}^t (e^{i\omega_{ni}t'}-e^{i\omega_{nm}t'})dt'\\\end{aligned}$$

合并一阶修正，跃迁速率：

$$W_{ni}=\frac{2\pi}{\hbar}\overline{|V_{ni}+\sum_m\frac{V_{nm}V_{mi}}{E_i-E_m}|^2}\rho(E)|_{E_n=E_i}$$

Harmonic Perturbation

考虑一个简谐微扰：

$$V(t)=Ve^{i\omega t}+V^\dagger e^{-i\omega t}$$

同样的方法：

$$c_n^{(1)}(t)=\frac{-i}{\hbar}\int_{0}^t\langle n|\hat V_I(t')|i\rangle dt'=\frac{1}{\hbar}[V_{ni}\frac{1-e^{i(\omega_{ni}+\omega)t}}{\omega_{ni}+\omega}+V_{ni}^\dagger\frac{1-e^{i(\omega_{ni}-\omega)t}}{\omega_{ni}-\omega}]$$

相比于常数微扰，主要区别在于：

$$\omega_{ni}\rightarrow \omega_{ni}\pm \omega$$

这里的$\hbar\omega$让人不禁想到了光子，对应的两种跃迁就是受激发射和吸收：

$$Emission:W_{ni}=\frac{2\pi}{\hbar}\overline{|V_{ni}|^2}\rho(E)|_{E_n=E_i-\hbar\omega}$$ $$Absorption:W_{ni}=\frac{2\pi}{\hbar}\overline{|V_{ni}^\dagger|^2}\rho(E)|_{E_n=E_i+\hbar\omega}$$

由：

$$|V_{ni}|^2=|V_{ni}^\dagger|^2$$

可以知道吸收和发射的概率是相等的，这叫做细致平衡原理。

Expotential Decay Perturbation-Energy Shift and Decay Width

考虑一个指数衰减的微扰：

$$V(t)=Ve^{-\eta t}$$

同样的方法：

$$c_n^{(1)}(t)=\frac{-i}{\hbar}\int_{0}^t\langle n|\hat V_I(t')|i\rangle dt'=\frac{-i}{\hbar}V_{ni}\frac{e^{\eta t+i\omega_{ni}t}}{\eta+i\omega_{ni}}$$

跃迁概率：

$$|c_n^{(1)}(t)|^2=\frac{|V_{ni}|^2}{\hbar^2}\frac{e^{2\eta t}}{\eta^2+\omega_{ni}^2}$$

跃迁速率：

$$W_{ni}=\frac{2|V_{ni}|^2}{\hbar^2}\frac{\eta e^{2\eta t}}{\eta^2+\omega_{ni}^2}$$

现在要求$t\rightarrow \infty$或者说，$\eta\rightarrow 0^+$，这时候可以得到：

$$\lim_{\eta\rightarrow 0^+}\frac{\eta}{\eta^2+\omega_{ni}^2}=\pi\delta(\omega_{ni})$$ $$W_{ni}=\frac{2\pi}{\hbar}|V_{ni}|^2\delta(E_n-E_i)$$

这同样满足费米黄金规则。

First Order

$$c_i^{(1)}(t)=\frac{-i}{\hbar}\lim_{t_0\rightarrow \infty}\int_{t_0}^t\langle i|\hat V_I(t')|i\rangle dt'=\frac{-i}{\hbar\eta}V_{ii}e^{\eta t}$$

Second Order

$$c_i^{(2)}(t)=\frac{-1}{\hbar^2}|V_{ii}|^2\frac{e^{2\eta t}}{2\eta^2}+\frac{-i}{\hbar}\sum_{m\neq i}\frac{|V_{mi}^2|e^{2\eta t}}{2\eta(E_i-E_m+i\hbar \eta)}$$

合并起来，跃迁概率为：

$$c_i(t)=1-\frac{i}{\hbar\eta}V_{ii}e^{\eta t}-\frac{1}{\hbar^2}|V_{ii}|^2\frac{e^{2\eta t}}{2\eta^2}+\frac{-i}{\hbar}\sum_{m\neq i}\frac{|V_{mi}^2|e^{2\eta t}}{2\eta(E_i-E_m+i\hbar \eta)}$$ $$\frac{\dot{c_i}}{c_i}=\frac{-i}{\hbar\eta}V_{ii}+\frac{-i}{\hbar}\sum_{m\neq i}\frac{|V_{mi}^2|}{E_i-E_m+i\hbar \eta}$$

解得：

$$c_i(t)=e^{-i\Delta_i t/\hbar},E_i\rightarrow E_i+\Delta_i$$

其中：

$$\Delta_i^{(1)}=V_{ii}$$ $$\Re\Delta_i^{(2)}=\Pr.\sum_{m\neq i}\frac{|V_{mi}|^2}{E_i-E_m}$$ $$\Im\Delta_i^{(2)}=-\pi\sum_{m\neq i}|V_{mi}|^2\delta(E_i-E_m)=-\frac{\hbar}{2}\sum_{m\neq i}W_{im}$$

定义：

$$\Gamma_i=-2\Im\Delta_i$$

那么：

$$|c_i(t)|^2=e^{-\Gamma_i t/\hbar}$$

这意味着$1/\Gamma_i$是衰减时间，$\Gamma_i$是对应的能级展宽。

Lamb Shift

经典中的真空无光子涨落，自然不会出现电子光子耦合，但是在量子力学中，真空中的光子涨落会导致电子的能级发生变化，这就是Lamb Shift。

写出微扰项：

$$\hat \delta H=-\frac{e}{m}\vec{A}\vec{p}+\frac{e^2}{2m}\vec{A}^2$$