双态系统是简单的量子系统
氨分子翻转、氢分子离子、磁场中1/2自旋粒子等都是有实际意义的双态系统。
一般量子系统不会只有两个能级,不过如果两个能级靠得很近,而离开其它能级又较远,在特定问题中它们与其它能级间的跃迁无需考虑,这种情况也可以视为双态系统。

能级离散系统中薛定谔方程的矩阵形式

含时薛定谔方程

波函数 ψ(x)\psi(x) 可表示为 ∣ψ(t)⟩|\psi(t)\rang,称为态矢量。取一套完备正交归一函数 φj(j=1,2,3⋯ )\varphi_{j}(j=1,2,3\cdots ),表示为 ∣j⟩|j\rang,称为基矢。用基矢展开态矢量,有

∣ψ(t)⟩=∑jCj∣j⟩|\psi(t)\rang=\sum_{j}C_j|j\rang \quad

其中 Cj=⟨j∣ψ(t)⟩C_j=\lang j|\psi(t)\rang,CjC_j 称为展开系数,表示态矢量 ∣ψ(t)⟩|\psi(t)\rang 分解为 ∣j⟩|j\rang 的概率幅。

用态矢量表示薛定谔方程,有

iℏ∂∂t∣ψ(t)⟩=H^∣ψ(t)⟩i\hbar \frac{\partial }{\partial t}|\psi(t)\rang=\hat{H}|\psi(t)\rang

由

∣ψ(t)⟩=∑jCj∣j⟩=∑j∣j⟩Cj=∑j∣j⟩⟨j∣ψ(t)⟩|\psi(t)\rang=\sum_{j}C_j|j\rang = \sum_{j}|j\rang C_j = \sum_{j}|j\rang \lang j|\psi(t)\rang

代入薛定谔方程,得到

iℏ∂∂t∣ψ(t)⟩=∑jH^∣j⟩⟨j∣ψ(t)⟩→左右同乘 ⟨i∣iℏ∂∂t⟨i∣ψ(t)⟩=∑j⟨i∣H^∣j⟩⟨j∣ψ(t)⟩i\hbar \frac{\partial }{\partial t}|\psi(t)\rang = \sum_{j}\hat{H}|j\rang \lang j|\psi(t)\rang \xrightarrow{\text{左右同乘 }\lang i|} i\hbar \frac{\partial }{\partial t}\lang i|\psi(t)\rang = \sum_{j}\lang i|\hat{H}|j \rang \lang j|\psi(t) \rang

⇒iℏ∂∂tCi=∑jHijCj\Rightarrow i\hbar \frac{\partial }{\partial t}C_i=\sum_{j}H_{ij}C_j

式中矩阵元 Hij=⟨i∣H^∣j⟩=∫φi∗H^φjdx\displaystyle H_{ij}=\lang i|\hat{H}|j\rang=\int \varphi^{*}_{i}\hat{H}\varphi_{j}\mathrm{d}x,同时上式为薛定谔方程在离散能级情况下的形式。

对于双态系统,态空间是二维的,所以只有两个基。薛定谔方程的矩阵形式为:

iℏ∂∂t(C1C2)=(H11H12H21H22)(C1C2)i\hbar \frac{\partial }{\partial t} \begin{pmatrix} C_1 \\ C_2 \end{pmatrix} = \begin{pmatrix} H_{11}&H_{12}\\ H_{21}&H_{22} \end{pmatrix} \begin{pmatrix} C_1\\ C_2 \end{pmatrix}

此处哈密顿算符变成哈密顿矩阵。

定态薛定谔方程

(H11H12H21H22)(C1′C2′)=E(C1′C2′)\begin{pmatrix} H_{11}&H_{12}\\ H_{21}&H_{22} \end{pmatrix} \begin{pmatrix} C_1'\\ C_2' \end{pmatrix} = E \begin{pmatrix} C_1'\\ C_2' \end{pmatrix}

其中 EE 为本征值,(C1′C2′)\begin{pmatrix} C_1'\\C_2' \end{pmatrix} 为本征矢。

若取哈密顿算符的本征态为基,记为 ∣β⟩(β=I,II)|\beta\rang(\beta=\mathrm{I},\mathrm{II}),有

H^∣I⟩=EI∣I⟩H^∣II⟩=EII∣II⟩\begin{gather*} \hat{H}|\mathrm{I}\rang =E_{\mathrm{I}}|\mathrm{I}\rang\\ \hat{H}|\mathrm{II}\rang = E_{\mathrm{II}}|\mathrm{II}\rang \end{gather*}

矩阵元 Hαβ=⟨α∣H^∣β⟩=EαδαβH_{\alpha \beta}=\lang \alpha|\hat{H}|\beta\rang=E_{\alpha}\delta_{\alpha\beta},可以看出此时哈密顿矩阵是对角的。

(HI IHI IIHII IHII II)=(EI00EII)\begin{pmatrix} H_{\mathrm{I\ I}}&H_{\mathrm{I\ II}}\\ H_{\mathrm{II\ I}}&H_{\mathrm{II\ II}} \end{pmatrix} = \begin{pmatrix} E_{\mathrm{I}} &0\\ 0&E_{\mathrm{II}} \end{pmatrix}

此时定态薛定谔方程的形式为

(EI00EII)(C1′C2′)=E(C1′C2′)\begin{pmatrix} E_{\mathrm{I}} &0\\ 0&E_{\mathrm{II}} \end{pmatrix} \begin{pmatrix} C_1'\\C_2' \end{pmatrix} = E \begin{pmatrix} C_1'\\C_2' \end{pmatrix}

可以得到能量本征值为 Eα(α=I,II)E_{\alpha}(\alpha=\mathrm{I,II}),EIE_{\mathrm{I}} 对应本征矢 (CI′CII′)=(10)=∣I⟩\begin{pmatrix} C_{\mathrm{I}}'\\ C_{\mathrm{II}}' \end{pmatrix} =\begin{pmatrix} 1\\0 \end{pmatrix} =|\mathrm{I}\rang,EIIE_{\mathrm{II}} 对应本征矢 (CI′CII′)=(01)=∣II⟩\begin{pmatrix} C_{\mathrm{I}}'\\ C_{\mathrm{II}}' \end{pmatrix} =\begin{pmatrix} 0\\1 \end{pmatrix} =|\mathrm{II}\rang

氨分子

对于一个氨分子,其N原子有一个镜像对称位置,对应的量子态分别为 ∣1⟩|1\rang 和 ∣2⟩|2\rang,其他自由度不计。若将 ∣1⟩|1\rang 和 ∣2⟩|2\rang 取为基,那么任意时刻分子的状态 ∣ψ⟩=∣1⟩⟨1∣ψ⟩+∣2⟩⟨2∣ψ⟩≡∣1⟩C1+∣2⟩C2|\psi\rang =|1\rang \lang 1|\psi\rang+|2\rang \lang 2|\psi\rang\equiv |1\rang C_1+|2\rang C_2

假设情况

首先考虑一种特殊情况,即假如系统一开始处于态 ∣1⟩|1\rang,以后不再有机会进入态 ∣2⟩|2\rang,反之亦然。那么考虑到非对角元会引起两个态之间的跃迁,则 H12=0,H21=0H_{12}=0,H_{21}=0。由薛定谔方程的矩阵形式

iℏ∂∂t(C1C2)=(H11H12H21H22)(C1C2)i\hbar \frac{\partial }{\partial t} \begin{pmatrix} C_1 \\ C_2 \end{pmatrix} = \begin{pmatrix} H_{11}&H_{12}\\ H_{21}&H_{22} \end{pmatrix} \begin{pmatrix} C_1\\ C_2 \end{pmatrix}

可以得到

iℏ∂C1∂t=H11C1⇒C1(t)=C1(0)e−iH11t/ℏiℏ∂C2∂t=H22C2⇒C2(t)=C2(0)e−iH22t/ℏ\begin{gather*} i\hbar \frac{\partial C_1}{\partial t}=H_{11}C_1 \Rightarrow C_1(t)=C_1(0)e^{-iH_{11}t /\hbar}\\ \quad \\ i\hbar \frac{\partial C_2}{\partial t}=H_{22}C_2 \Rightarrow C_2(t)=C_2(0)e^{-iH_{22}t /\hbar} \end{gather*}

C1(t)C_1(t) 和 C2(t)C_2(t) 为氨分子处于 E1=H11,E2=H22E_1=H_{11},E_2=H_{22} 定态的概率幅。再由对称性,H11=H22=E0H_{11}=H_{22}=E_0。

真实情况

但这并非氨分子的实际行为。氨分子中的N原子可以通过隧道效应在 ∣1⟩,∣2⟩|1\rang,|2\rang 态之间穿越。此时则假设矩阵元 H12=H21=A≠0(A<0)H_{12}=H_{21}=A\neq 0(A<0),则薛定谔方程为

iℏ∂∂t(C1C2)=(E0AAE0)(C1C2)i\hbar \frac{\partial }{\partial t} \begin{pmatrix} C_1 \\ C_2 \end{pmatrix} = \begin{pmatrix} E_0&A\\ A&E_0 \end{pmatrix} \begin{pmatrix} C_1\\ C_2 \end{pmatrix}

通过一些数学上的技巧,可以解得

C1(t)=a2e−i(E0+A)t/ℏ+b2e−i(E0−A)t/ℏC2(t)=a2e−i(E0+A)t/ℏ−b2e−i(E0−A)t/ℏ\begin{gather*} C_1(t)=\frac{a}{2}e^{-i(E_0+A)t /\hbar}+\frac{b}{2}e^{-i(E_0-A)t/\hbar}\\ \quad\\ C_2(t)=\frac{a}{2}e^{-i(E_0+A)t /\hbar}-\frac{b}{2}e^{-i(E_0-A)t/\hbar}\\ \end{gather*}

式中 a,ba,b 为待定常量。而氨分子波函数 ∣ψ⟩=C1∣1⟩+C2∣2⟩|\psi\rang =C_1|1\rang+C_2|2\rang。

若 b=0b=0,则两项具有相同的频率 ω=(E0+A)/ℏ\omega=(E_0+A) /\hbar,具有相同的概率幅 C1=C2=a2e−i(E0+A)t/ℏ\displaystyle C_1=C_2=\frac{a}{2}e^{-i(E_0+A)t /\hbar}。那么此时分子具有确定的能量 EI=E0+AE_{\mathrm{I}}=E_0+A,分子处于定态。定态的态矢量为

∣ψI(t)⟩=12(∣1⟩+∣2⟩)e−i(E0+A)t/ℏ|\psi_{\mathrm{I}}(t)\rang =\frac{1}{\sqrt{2}}(|1\rang +|2\rang)e^{-i(E_0+A)t /\hbar}

态矢量通常又用与时间无关的态矢量来表示

∣I⟩=12(∣1⟩+∣2⟩)|\mathrm{I}\rang =\frac{1}{\sqrt{2}}(|1\rang+|2\rang)

该态矢量也记为 ∣+⟩|+\rang。同时可以看出,这就是 H^\hat{H} 的本征态, H^∣I⟩=EI∣I⟩\hat{H}|\mathrm{I}\rang=E_{\mathrm{I}}|\mathrm{I}\rang。

若 a=0a=0,和上方类似的,我们可以得到 C1=−C2=b2e−i(E0−A)t/ℏC_1=-C_2=\displaystyle \frac{b}{2}e^{-i(E_0-A)t /\hbar},表示N原子“在上”或者“在下”具有相反的概率幅。分子处于定态态矢量为

∣ψII(t)⟩=12(∣1⟩−∣2⟩)e−i(E0−A)t/ℏ|\psi_{\mathrm{II}}(t)\rang =\frac{1}{\sqrt{2}}(|1\rang -|2\rang)e^{-i(E_0-A)t /\hbar}

∣II⟩=12(∣1⟩−∣2⟩)|\mathrm{II}\rang =\frac{1}{\sqrt{2}}(|1\rang-|2\rang)

该态矢量记为 ∣−⟩|-\rang,也是 H^\hat{H} 的本征态,H^∣II⟩=EII∣II⟩\hat{H}|\mathrm{II}\rang=E_{\mathrm{II}}|\mathrm{II}\rang

如果从定态薛定谔方程讨论

(E0AAE0)(x1x2)=(x1x2)\begin{pmatrix} E_0&A\\ A &E_0 \end{pmatrix} \begin{pmatrix} x_1\\ x_2 \end{pmatrix} = \begin{pmatrix} x_1\\ x_2 \end{pmatrix}

同样可以求得能量本征值与本征矢。

讨论

如果 t=0t=0 时刻,分子处于态 ∣1⟩|1\rang,则 C1(0)=1,C2(0)=0⇒a=b=1C_1(0)=1,C_2(0)=0 \Rightarrow a=b=1,解得

C1(t)=e−iE0t/ℏcos⁡AtℏC2(t)=−e−iE0t/ℏsin⁡Atℏ\begin{gather*} C_1(t)=e^{-iE_0t /\hbar}\cos \frac{At}{\hbar}\\ \quad \\ C_2(t)=-e^{-iE_0t /\hbar}\sin \frac{At}{\hbar} \end{gather*}

此时系统能量不确定。同时波函数为两个振动的叠加,会出现拍的现象。在 tt 时刻,分子处于态 ∣1⟩,∣2⟩|1\rang,|2\rang 的概率分别为

P1=∣C1(t)∣2=cos⁡2AtℏP2=∣C2(t)∣2=sin⁡2AtℏP_1=\left\vert C_1(t) \right\vert ^{2}=\cos ^{2}\frac{At}{\hbar}\quad P_2=\left\vert C_2(t) \right\vert ^{2}=\sin ^{2}\frac{At}{\hbar}

量子与经典的重要区别:统一的能级分裂为二,分裂为两个定态 (E0→E0±A)(E_0\rightarrow E_0\pm A)。这就是能量的翻转分裂。同时从实验数据可知,AA 非常小,因此能极差也非常小。激发氨分子能级翻转仅需要微波提供的能量。