哈密顿算符的群

哈密顿算符的对称群

一个体系可用哈密顿量 H^\hat{H} 来描述,体系的对称性就是由哈密顿量的对称性决定,那么我们就要研究哈密顿量在群下的变换情况。

在量子力学中,算符变换为

A^经变换 P^A^=P^A^P^1\hat{A}\xrightarrow{\text{经变换 }\hat{P}}\hat{A}^{\prime}=\hat{P}\hat{A}\hat{P}^{-1}

A^\hat{A} 为体系的哈密顿算符 H^\hat{H},则变换后哈密顿量为

H^=P^H^P^1\hat{H}^{\prime}=\hat{P}\hat{H}\hat{P}^{-1}

若体系哈密顿量在变换 P^\hat{P} 下不变,即 H^=P^H^P^1=H^\hat{H}^{\prime}=\hat{P}\hat{H}\hat{P}^{-1}=\hat{H},那么我们称 H^\hat{H} 具有 P^\hat{P} 操作的对称性,等价于

P^H^=H^P^[H^,P^]=0\hat{P}\hat{H}=\hat{H}\hat{P}\quad \Rightarrow \quad \begin{bmatrix} \hat{H},\hat{P} \end{bmatrix}=0

也就是当 H^\hat{H}P^\hat{P} 对易时,说明 H^\hat{H} 具有 P^\hat{P} 的对称性。

所有使体系哈密顿量 H^\hat{H} 不变的变换 P^\hat{P} 构成一个群

G={P^P^H^P^1=H^}G = \{ \hat{P} \mid \hat{P}\hat{H}\hat{P}^{-1} = \hat{H} \}

称为哈密顿算符的群或薛定谔方程的群,或称 GGH^\hat{H} 的对称群,显然 H^\hat{H} 与该群中任一变换对易 [H^,P^]=0P^G[\hat{H}, \hat{P}] = 0, \quad \forall \hat{P} \in G

可以证明,若哈密顿量是 r\boldsymbol{r} 的函数 H^(r)\hat{H}(\boldsymbol{r}) 且变换为空间变换 P^R\hat{P}_R 时有

P^RH^(r)P^R1=H^(R1r)\hat{P}_R \hat{H}(\boldsymbol{r}) \hat{P}_R^{-1} = \hat{H}(R^{-1}\boldsymbol{r})

H^(r)\hat{H}(\boldsymbol{r}) 变换为 H^(R1r)\hat{H}(R^{-1}\boldsymbol{r}),形式上和函数的变换一样。如果哈密顿量在该变换下不变,则要求 H^(R1r)=H^(r)\hat{H}(R^{-1}\boldsymbol{r}) = \hat{H}(\boldsymbol{r})

定理一: H^\hat{H} 的具有同一能量本征值的所有线性无关的简并波函数构成其对称群 GG 的一个表示的基

EnE_nH^\hat{H} 的某个简并能级,且该能级的简并波函数共有 dd 个,记为 ψi(i=1,2,,d)\psi_i (i = 1,2,\cdots,d)

H^ψi=Enψi\hat{H}\psi_i = E_n\psi_i

任取 P^G\hat{P} \in GH^P^=P^H^\hat{H}\hat{P} = \hat{P}\hat{H},于是

H^(P^ψi)=P^H^ψi=P^Enψi=En(P^ψi)\hat{H}(\hat{P}\psi_i) = \hat{P}\hat{H}\psi_i = \hat{P}E_n\psi_i = E_n(\hat{P}\psi_i)

说明 P^ψi\hat{P}\psi_i 依然是 H^\hat{H} 的属于本征值 EnE_n 的本征函数,又 {ψi}\{\psi_i\} 彼此线性无关,那么 P^ψi\hat{P}\psi_i 必然可用 {ψi}\{\psi_i\} 线性展开为

P^ψi=j=1dψjDji(P^)\hat{P}\psi_i = \sum_{j=1}^d \psi_j D_{ji}(\hat{P})

于是,其线性组合系数构成的矩阵 D(P^)D(\hat{P}) 形成了群 GG 在以 {ψ1,ψ2,,ψd}\{\psi_1, \psi_2, \cdots, \psi_d\} 为基时的一个 dd 维表示。一般来说,这个表示 DD 可能是可约的,设约化后为

T1D(g)T=papDp(g)T^{-1}D(g)T = \sum_p \bigoplus a_p D^p(g)

相应不可约表示的基 {ψiαpp=1,,r;i=1,,dp;α=1,ap}\{\psi_{i\alpha}^p | p = 1,\cdots,r; i = 1,\cdots,d_p; \alpha = 1,\cdots a_p\} 则是原来 dd 个波函数 {ψi}\{\psi_i\} 经线性变换 TT 而得到,仍是相应于 EnE_n 的本征态

P^ψiαp=jψjαpDjip(P^)a1d1+a2d2++ardr=d\hat{P}\psi_{i\alpha}^p = \sum_j \psi_{j\alpha}^p D_{ji}^p(\hat{P}) \quad \quad \begin{matrix} a_1 d_1 + a_2 d_2 + \\ \cdots + a_r d_r = d \end{matrix}

形式上将对应不可约表示 DpD^p 的能级标记为 εnpα(α=1,,ap)\varepsilon_{np}^\alpha \quad (\alpha = 1,\cdots,a_p),约化后的基对应本征值 εn11,εn12,,εn1a1p=1,,εnr1,εnr2,,εnrarp=r\underbrace{\varepsilon_{n1}^1, \varepsilon_{n1}^2, \cdots, \varepsilon_{n1}^{a_1}}_{p=1} , \cdots , \underbrace{\varepsilon_{nr}^1, \varepsilon_{nr}^2, \cdots, \varepsilon_{nr}^{a_r}}_{p=r},只不过这里都等于 EnE_n 罢了: εnpα=En\varepsilon_{np}^\alpha = E_n

将属于同一不可约表示 DpD^pdpd_p 重简并能级 εnpα\varepsilon_{np}^\alpha (对于固定的 ppα\alpha) 称为正常简并,也称必然简并;而对应不同 p,αp, \alpha 的简并称为偶然简并


  • 也就是需要同一个且出现的顺序也相同的不可约表示的简并态称为必然简并,否则称为偶然简并。
  • 必然简并是受对称群保护的。每一个不可约表示都可以看作是对角的单元,无法对不可约表示进一步约化,所以不可约表示对应的基都是简并的。偶然简并是由群对称性之外的其他因素引起的。
  • 可约表示对应的简并都是偶然简并。

定理二:H^\hat{H} 的对称群为 GG,若无偶然简并,则 H^\hat{H} 的一组线性无关本征态荷载群 GG 的一个不可约表示的充要条件是它们属于同一能量本征值且是一个最大线性无关组

定理三:H^\hat{H} 的对称群为 GG,函数集 {ϕnp}\{\phi_n^p\}{ψmq}\{\psi_m^q\} 分别是群 GG 的第 ppqq 个不可约么正表示的基,那么算符 H^\hat{H} 的矩阵元满足下式

(ϕ_np,H^ψ_mq)=δ_pqδ_nm(ϕ_kp,H^ψ_kp)(\phi\_n^p, \hat{H}\psi\_m^q) = \delta\_{pq}\delta\_{nm}(\phi\_k^p, \hat{H}\psi\_k^p)

其中 (ϕkp,H^ψkp)(\phi_k^p, \hat{H}\psi_k^p) 是与 kk 无关的常数

  • ϕnp{\phi_n^p}ψmq{\psi_m^q} 不是 H^\hat{H} 的本征函数,而是群 GG 的不可约表示的基函数。
  • {H^ψmq}\{\hat{H}\psi_m^q\} 同样是 GG 的第 qq 个不可约表示的基

哈密顿量的简化求解

有了上面三个定理后,我们就可以利用 H^\hat{H} 的对称群 GG 简化对 H^\hat{H} 的求解。

设有 H^\hat{H} 的一组完备基 {ϕii=1,,n}\{|\phi_i\rangle | i = 1, \cdots, n\} 且正交归一 ϕiϕj=δij\langle\phi_i|\phi_j\rangle = \delta_{ij}求解 H^\hat{H} 本征值 EE 及本征态 ψ|\psi\rangle:

H^ψ=Eψ\hat{H}|\psi\rangle = E|\psi\rangle

ψ=jcjϕj|\psi\rangle = \sum_j c_j |\phi_j\rangle,得到 H^jcjϕj=Ejcjϕj\hat{H} \sum_j c_j |\phi_j\rangle = E \sum_j c_j |\phi_j\rangle。对两边同时作用 ϕi\langle\phi_i|,得到

jϕiH^ϕjHijcj=Ejϕiϕjδijcj\sum_j \underbrace{\langle\phi_i|\hat{H}|\phi_j\rangle}_{H_{ij}} c_j = E \sum_j \underbrace{\langle\phi_i|\phi_j\rangle}_{\delta_{ij}} c_j

也就是

jHijcj=EciHc=Ec\sum_j H_{ij}c_j = E c_i \quad \Rightarrow \quad Hc = Ec

得到矩阵形式

(H11H12H1nH21H22H2nHn1Hn2Hnn)(c1c2cn)=E(c1c2cn)\begin{pmatrix} H_{11} & H_{12} & \cdots & H_{1n} \\[0.5em] H_{21} & H_{22} & \cdots & H_{2n} \\[0.5em] \vdots & \vdots & \ddots & \vdots \\[0.5em] H_{n1} & H_{n2} & \cdots & H_{nn} \end{pmatrix} \begin{pmatrix} c_1 \\[0.5em] c_2 \\[0.5em] \vdots \\[0.5em] c_n \end{pmatrix} = E \begin{pmatrix} c_1 \\[0.5em] c_2 \\[0.5em] \vdots \\[0.5em] c_n \end{pmatrix}

可以看到,在给定基后,求解算符 H^\hat{H} 的本征值问题就转化为求解 n×nn \times n 矩阵 HH 的本征值问题,但是当 nn 很大时难于求解,我们可以利用 H^\hat{H} 的对称群 GG 来简化求解过程:将矩阵 HH 化为块对角形式,此时若要求解 H^\hat{H} 的本征值问题,只需分别求解各个对角块的本征值问题即可,大大简化了计算量。


对已有的基重新线性组合并按群 GG 的不可约表示分类(可利用投影算符)为 对称化的波函数 {ϕiαp},p=1,,r; i=1,,dp; α=1,,ap\{|\phi_{i\alpha}^p\rangle\},\quad p=1,\cdots,r; \ i=1,\cdots,d_p; \ \alpha=1,\cdots,a_p。因为线性组合不改变基的个数,那么有

a1d1+a2d2++ardr=na_1 d_1 + a_2 d_2 + \cdots + a_r d_r = n

根据正交归一性 ϕiαpϕjβq=δpqδijδαβ\langle\phi_{i\alpha}^p|\phi_{j\beta}^q\rangle = \delta_{pq}\delta_{ij}\delta_{\alpha\beta},那么在 ϕiαp\ket{\phi_{i\alpha}^p}ϕjβq|\phi_{j\beta}^q\rangle 基下,记算符 H^\hat{H} 的矩阵元为 Hiα,jβp,q=ϕiαpH^ϕjβqH_{i\alpha,j\beta}^{p,q} = \langle\phi_{i\alpha}^p|\hat{H}|\phi_{j\beta}^q\rangle,那么哈密顿量的矩阵方程变为

q,β,jHiα,jβp,qCjβq=ECiαp\sum_{q,\beta,j} H_{i\alpha,j\beta}^{p,q} C_{j\beta}^{q} = E C_{i\alpha}^{p}

根据定理三有

Hiα,jβp,qϕiαpH^ϕjβq=δpqδijHαβ(p)H_{i\alpha,j\beta}^{p,q} \equiv \langle\phi_{i\alpha}^p|\hat{H}|\phi_{j\beta}^q\rangle = \delta_{pq} \delta_{ij} H_{\alpha\beta}^{(p)}

此时矩阵已经块对角化了,指标 p,qp,q 对应不同的块(不可约表示),而每个块内又划分出了不同的块区域,i,ji,j 对应同一块内的不同子块,α,β\alpha,\beta 对应同一子块内的矩阵元(同一不可约表示的不同出现顺序)。块内的子块划分如下

H(p=1)=(i=1i=2(H11H12H21H22)α,β00(H11H12H21H22)α,β)H^{(p=1)} = \begin{pmatrix} \begin{array}{c|c} \mathbf{i=1} & \mathbf{i=2} \\ \hline \begin{pmatrix} \textcolor{red}{H_{11}} & \textcolor{red}{H_{12}} \\ \textcolor{red}{H_{21}} & \textcolor{red}{H_{22}} \end{pmatrix}_{\alpha,\beta} & \mathbf{0} \\ \hline \mathbf{0} & \begin{pmatrix} \textcolor{red}{H_{11}} & \textcolor{red}{H_{12}} \\ \textcolor{red}{H_{21}} & \textcolor{red}{H_{22}} \end{pmatrix}_{\alpha,\beta} \end{array} \end{pmatrix}

其中 Hαβ(p)H_{\alpha\beta}^{(p)} 对应 Hiα,iβp,pH_{i\alpha,i\beta}^{p,p}。因此我们只需求块内小矩阵的本征值为题即可,代入化简得

βHαβ(p)Ciβp=ECiαp\sum_{\beta} H_{\alpha\beta}^{(p)} C_{i\beta}^{p} = E C_{i\alpha}^{p}

每个小矩阵都是 ap×apa_p \times a_p 维的,所以只需对不同的 pp 求解一个 ap×apa_p \times a_p 的矩阵 H(p)[Hαβ(p)]H^{(p)} \equiv [H_{\alpha\beta}^{(p)}](α,β=1,2,,ap)(\alpha,\beta = 1,2,\cdots,a_p) 即可,且每个 ii 都对应同一个 H(p)H^{(p)} 矩阵。

其具体矩阵表示如下

哈密顿量矩阵块对角化

每个小矩阵的结构如下

哈密顿量矩阵子块

小矩阵块重复的次数取决于表示的维数 dpd_p(简并重数),大小取决于表示出现的次数 apa_p。由于相同 pp 的小矩阵块是相同的,所以总共只需求解 rr 个小矩阵的本征值问题,每个小矩阵的维数为 ap×apa_p \times a_pn×nn \times n 的矩阵本征值问题化为 rrap×apa_p \times a_p 矩阵的本征值问题,大大简化了计算量。

我们找到了一组对称性适配基底,使得巨大的哈密顿矩阵崩解成了许多互不相干的小块。这样就不需要对角化整个大矩阵,只需要分别对角化每一个小块 Hαβ(p)H_{\alpha\beta}^{(p)} 即可,这就是对称群的简化作用。

微扰引起的能级劈裂

设体系 H^0\hat{H}_0 具有群 G0G_0 的对称性,其本征值和本征态按 G0G_0 的不可约表示分类为 εαp\varepsilon_\alpha^pψiαp\psi_{i\alpha}^p (属于第 α\alpha 次出现的第 pp 个不可约表示的第 ii 列基):

H^0ψiαp=εαpψiαp\hat{H}_0 \psi_{i\alpha}^p = \varepsilon_\alpha^p \psi_{i\alpha}^p

其中 {ψiαpi=1,,dp}\{\psi_{i\alpha}^p | i = 1, \cdots, d_p\} 对应同一能量本征值 εαp\varepsilon_\alpha^p,该本征值具有 dpd_p 重简并。

现在引入一个微扰项 H^\hat{H}',使得体系的哈密顿量变为

H^=H^0+H^\hat{H} = \hat{H}_0 + \hat{H}'

体系哈密顿量仍保持 G0G_0 的对称性。若微扰 H^\hat{H}' 同样具有 G0G_0 的对称性,那么仍旧可以应用定理三

ψiαpH^ψjβq=δpqδijΔαβp(Δαβp=ψiαpH^ψiβp)\langle\psi_{i\alpha}^p|\hat{H}'|\psi_{j\beta}^q\rangle = \delta_{pq}\delta_{ij}\Delta_{\alpha\beta}^p \quad \quad (\Delta_{\alpha\beta}^p = \langle\psi_{i\alpha}^p|\hat{H}'|\psi_{i\beta}^p\rangle)

则本征能量的改变 Δε\Delta\varepsilonΔ(p)=[Δαβp]\Delta^{(p)} = [\Delta_{\alpha\beta}^p] 有关。

现假设 H^0\hat{H}_0 存在偶然简并,有如下 99 重简并态ψiαt(i=1,2,3;α=1,2)\psi_{i\alpha}^t (i = 1,2,3; \alpha = 1,2)ψiαs(i=1;α=1,2,3)\psi_{i\alpha}^s (i = 1; \alpha = 1,2,3),也就是说 dt=3,at=2d_t = 3, a_t = 2ds=1,as=3d_s = 1, a_s = 3。设其简并能量为 ε0\varepsilon_0,则有 ε0=ε1t=ε2t=ε1s=ε2s=ε3s\varepsilon_0 = \varepsilon_1^t = \varepsilon_2^t = \varepsilon_1^s = \varepsilon_2^s = \varepsilon_3^s

{ψiαp}\{\psi_{i\alpha}^p\} (遍历所有 p,α,ip,\alpha,i) 为基构造 H^\hat{H} 的哈密顿矩阵为,得到块对角后的 99 个小矩阵

H=(H1)Hiα,jβp,qψiαpH^ψjβqH = \begin{pmatrix} \ddots & & \\ & \boxed{H_1} & \\ & & \ddots \end{pmatrix} \quad \quad H_{i\alpha,j\beta}^{p,q} \equiv \langle\psi_{i\alpha}^p|\hat{H}|\psi_{j\beta}^q\rangle

其中 H^0\hat{H}_0 在九重简并的条件下得到相同的能量 ε0\varepsilon_0,所以其对应 ε0\varepsilon_0 倍的单位矩阵,所以上式中的 H1H_1 可写为

H1=ε0I+(Δ(t)Δ(t)Δ(t)Δ(s))H_1 = \varepsilon_0 I + \begin{pmatrix} \boxed{\Delta^{(t)}} & & & \\ & \boxed{\Delta^{(t)}} & & \\ & & \boxed{\Delta^{(t)}} & \\ & & & \boxed{\Delta^{(s)}} \end{pmatrix}

其中 Δ(t)\Delta^{(t)}2×22 \times 2 矩阵, 重复出现 33 次,有 22 个本征值 Δε1t,Δε2t\Delta\varepsilon_1^t, \Delta\varepsilon_2^t,每个本征值都是三重简并的
Δ(s)\Delta^{(s)}3×33 \times 3 矩阵, 只出现 11 次,有 33 个本征值 Δε1s,Δε2s,Δε3s\Delta\varepsilon_1^s, \Delta\varepsilon_2^s, \Delta\varepsilon_3^s,每个本征值一重简并(不简并)。

因此,微扰引起的能级劈裂结果为

对称微扰引起的能级劈裂

可见,若微扰 H^\hat{H}' 同样具有 G0G_0 的对称性,则原来正常简并的态间依然正常简并,也就是处于块矩阵内的各个小矩阵的相同位置的本征值依然简并;原来偶然简并 (不同 ppα\alpha) 的态可能发生劈裂。

H^\hat{H}' 的对称群为 G1G_1 且是 G0G_0 的子群,说明微扰项使对称性降低,原来 dpd_p 重简并的波函数荷载 G0G_0 的不可约表示 DpD^p,这些波函数同样可荷载子群 G1G_1 的表示,即 DpD^pG1G_1 中的分导表示(父群的表示用于子群时称为分导表示)。不过这个表示可能是可约的了,按 G1G_1 的不可约表示约化后为

Dp(G1)=a1D11a2D12arD1rD^p(G_1) = a_1 D_1^1 \oplus a_2 D_1^2 \oplus \cdots \oplus a_r D_1^r

D1p (p=1,2,,r)D_1^p \ (p = 1,2,\cdots,r)G1G_1 的不可约表示。若无偶然简并,原 dpd_p重简并在微扰下将劈裂为 a1+a2++ar=npa_1 + a_2 + \cdots + a_r = n_p 个不同的能级。

一般情况下,对称性降低将会使原来简并的能级发生劈裂。因为父群的不可约表示对子群来说一般是可约的,而可约表示对应的简并态都是偶然简并,所以对称性降低后,原来简并态将不再受对称群保护,从而可能发生劈裂。


现考虑一个晶体场劈裂(把孤立原子放入晶体场前后的能级劈裂),即考虑 O(3)\text{O}(3) 的不可约表示向具体的晶体点群的分导表示。

一个球对称孤立原子具有 O(3)\text{O}(3) 对称性,由于 O(3)\text{O}(3) 的不可约表示可由 SO(3)\text{SO}(3) 的不可约表示 D(l)D^{(l)} 得到,若忽略其自旋,则处于某壳层 ll 轨道内的电子态是 2l+12l + 1 重简并的,例如 ss 轨道: l=0l=0,一重简并;pp 轨道: l=1l=1,三重简并;dd 轨道: l=2l=2,五重简并。

现在把该原子放入具有某一点群对称性的晶体场中,例如 O\text{O} 点群,对称性降低: O(3)O\text{O}(3) \rightarrow \text{O},相当于加了一个具有 O\text{O} 群对称性的微扰势 VV,那么原来简并的能级就可能劈裂,那么如何劈裂呢?这就需要对 O(3)\text{O}(3) 群在 O\text{O} 群中的分导表示进行约化!而 O(3)\text{O}(3) 的不可约表示可由 SO(3)\text{SO}(3) 的不可约表示 D(l)D^{(l)} 构造出来,虽然 D(l)D^{(l)} 的矩阵元很复杂,但是约化表示时只需其特征标即可

O(3)\text{O}(3) 群可由 SO(3)\text{SO}(3) 群与反演 ii 生成

O(3)=SO(3)+iSO(3)\text{O}(3) = \text{SO}(3) +i \cdot \text{SO}(3)

ii 与 所有 RSO(3)R \in \text{SO}(3) 对易,则 ii 的表示矩阵与所有 D(l)D^{(l)} 对易,所以其表示为单位矩阵常数倍 D(i)=cID(i) = cI,又 D(i)2=D(e)=ID(i)^2 = D(e) = I,所以 c2=1 c=±1c^2 = 1 \Rightarrow \ c = \pm 1,于是 D(l)±(iR)=D(l)(i)D(l)(R)=±D(l)(R)D^{(l)\pm}(iR) =D^{(l)}(i)D^{(l)}(R) = \pm D^{(l)}(R),但是 对于 RR 来说仍有 D(l)±(R)=D(l)(R)D^{(l)\pm}(R) = D^{(l)}(R)

SO(3)\text{SO}(3) 的不可约表示 D(l)D^{(l)} (l=0,1,2,)(l = 0,1,2,\cdots) 可构造出 O(3)\text{O}(3) 的两个不可约表示 D(l)+D^{(l)+}D(l)D^{(l)-},具体是哪个根据基函数决定。

其基函数由角向和径向两部分组成

ψ(r,θ,ϕ)=f(r)Yml(θ,ϕ)\psi(r, \theta, \phi) = f(r)Y_m^l(\theta, \phi)

其中 Yml(θ,ϕ)Y_m^l(\theta, \phi) 就是 SO(3)\text{SO}(3) 群的表示。 dd 轨道 (l=2)(l=2) 的球谐函数可线性组合成 {dz2dx2y2,dxy,dxz,dyz}\{d_{z^2},d_{x^2-y^2}, d_{xy}, d_{xz}, d_{yz}\} 这五个常用的 dd 轨道

dz2=Y02=14r215π13(3z2r2)d_{z^2} = Y_0^2 = \frac{1}{4r^2}\sqrt{\frac{15}{\pi}} \textcolor{red}{\sqrt{\frac{1}{3}}(3z^2 - r^2)}

dx2y2=12(Y22+Y22)=14r215π(x2y2)d_{x^2-y^2} = \frac{1}{\sqrt{2}}(Y_2^2 + Y_{-2}^2) = \frac{1}{4r^2}\sqrt{\frac{15}{\pi}}\textcolor{red}{(x^2 - y^2)}

dxy=i2(Y22Y22)=14r215π2xyd_{xy} = \frac{-i}{\sqrt{2}}(Y_2^2 - Y_{-2}^2) = \frac{1}{4r^2}\sqrt{\frac{15}{\pi}}\textcolor{red}{2xy}

dxz=12(Y12Y12)=14r215π2xzd_{xz} = \frac{-1}{\sqrt{2}}(Y_1^2 - Y_{-1}^2) = \frac{1}{4r^2}\sqrt{\frac{15}{\pi}}\textcolor{red}{2xz}

dyz=i2(Y12+Y12)=14r215π2yzd_{yz} = \frac{i}{\sqrt{2}}(Y_1^2 + Y_{-1}^2) = \frac{1}{4r^2}\sqrt{\frac{15}{\pi}}\textcolor{red}{2yz}

而这五个 dd 轨道可以构成 O(3)\text{O}(3) 表示的基,由于他们都是二次型的,因此在空间反演下保持不变,所以能荷载与 D(2)+D^{(2)+} 等价的不可约表示,那么就具有相同特征标,D(l)+D^{(l)+} 的特征标计算式为

χ(l)+(α)=sin(l+12)αsinα2\chi^{(l)+}(\alpha)=\frac{\sin\left(l+\dfrac{1}{2}\right)\alpha}{\sin\dfrac{\alpha}{2}}

三个 pp 轨道荷载 D(1)D^{(1)-},不同于 dd 轨道的是,pp 轨道是线性型的,在空间反演下变号,因此荷载与 D(1)D^{(1)-} 等价的不可约表示。

计算可得

χ(1)(c2)=1,χ(1)(c3)=0,χ(1)(c4)=1χ(2)+(c2)=1,χ(2)+(c3)=1,χ(2)+(c4)=1\begin{aligned} \chi^{(1)-}(c_2) &= -1, & \chi^{(1)-}(c_3) &= 0, & \chi^{(1)-}(c_4) &= 1 \\[1em] \chi^{(2)+}(c_2) &= 1, & \chi^{(2)+}(c_3) &= -1, & \chi^{(2)+}(c_4) &= -1 \end{aligned}

即便对于 χ(l)(α)\chi^{(l)-}(\alpha),由于我们计算的都是正当旋转,因此特征标仍然不反号。

Oe8c33c426c26c4A111111A211111E21200T130111T230111χ(1)30111χ(2)+51111\begin{array}{c|ccccc} O & e & 8c_3 & 3c_4^2 & 6c_2 & 6c_4 \\[0.5em] \hline A_1 & 1 & 1 & 1 & 1 & 1 \\[0.5em] A_2 & 1 & 1 & 1 & -1 & -1 \\[0.5em] E & 2 & -1 & 2 & 0 & 0 \\[0.5em] T_1 & 3 & 0 & -1 & -1 & 1 \\[0.5em] T_2 & 3 & 0 & -1 & 1 & -1 \\[0.5em] \hline \chi^{(1)-} & 3 & 0 & -1 & -1 & 1 \\[0.5em] \chi^{(2)+} & 5 & -1 & 1 & 1 & -1 \end{array}

pp 轨道的劈裂情况:D(1)(O)=T1D^{(1)-}(O) = T_1。以 33pp 轨道为基的三维表示在 O 中不劈裂,仍是一个三维表示 T1T_1

dd 轨道的劈裂情况:D(2)+(O)=ET2D^{(2)+}(O) = E \oplus T_2。以 55dd 轨道为基的五维表示在 OO 中劈裂为两组,一组属二维不可约表示 EE,另一组属三维不可约表示 T2T_2

晶体场劈裂

矩阵元定理与选择定则

矩阵元定理

量子体系在含时微扰 HH' 的作用下从一个态到另一个态的跃迁速率可用“费米黄金规则”来计算

Tif=2πfHi2ρf=2πHfi2ρfT_{i \to f} = \frac{2\pi}{\hbar} |\langle f|H'|i\rangle|^2 \rho_f = \frac{2\pi}{\hbar} |H'_{fi}|^2 \rho_f

其中 TifT_{i \to f} 是跃迁速率 (ii 初态, ff 末态)Hfi2,|H'_{fi}|^2 是跃迁矩阵元,ρf\rho_f 是末态态密度。若 Hfi=0H'_{fi} = 0Tif=0T_{i \to f} = 0,初态 i|i\rangle 到末态 f|f\rangle 的跃迁是禁戒的仅当 Hfi0H'_{fi} \neq 0 时的跃迁才是允许的。当然这里我们忽略了 ρf\rho_f 的影响,因为其与对称性无关。

根据微扰 HH' 的对称性可以判断跃迁矩阵元 HfiH'_{fi} 是否为 00,那么这就可以判断哪些跃迁是禁戒的,哪些跃迁是允许的,这就是选择定则

设未被微扰体系为 H^0\hat{H}_0,其对称群为 GGH^0\hat{H}_0 的本征态按 GG 的不可约表示分类为 {ψiαp},p=1,,r; i=1,,dp; α=1,,ap\{\psi_{i\alpha}^p\},\quad p=1,\cdots,r; \ i=1,\cdots,d_p; \ \alpha=1,\cdots,a_p,那么从 ψiαp\psi_{i\alpha}^p 态到 ψjβq\psi_{j\beta}^q 态的跃迁矩阵元为

ψjβqHψiαp\langle\psi_{j\beta}^q | H' | \psi_{i\alpha}^p\rangle

矩阵元定理: 如果态 HψiαpH'|\psi_{i\alpha}^p\rangle 中并不包含依群 GG 的第 qq 个不可约表示的第 jj 列基变换的 ψjβq|\psi_{j\beta}^q\rangle 态的分量的话,则矩阵元

ψjβqHψiαp=0\langle\psi_{j\beta}^q | H' | \psi_{i\alpha}^p\rangle = 0

它们之间的跃迁是禁戒的

但是微扰 HH' 有时是一组量比如 {Hx,Hy,Hz}\{H'_x, H'_y, H'_z\},记为 {Hμ}\{H'_\mu\},我们要求不仅初末态能荷载群的表示,还要求微扰算符 {Hμ}\{H'_\mu\} 也能荷载群 GG 的一个表示,设这个表示为 DD' (如仅一个则为一维表示)

ΓgHμΓg1=νHνDνμ(g)\Gamma_g H'_\mu \Gamma_g^{-1} = \sum_\nu H'_\nu D'_{\nu\mu}(g)

那么就有

ΓgHμψiαp=(ΓgHμΓg1)Γgψiαp=νHνDνμ(g)jψjαpDjip(g)=νjHνψjαpDνμ(g)Djip(g)=ν,jHνψjαp[D(g)Dp(g)]νj,μi\begin{aligned} \Gamma_g H'_\mu \psi_{i\alpha}^p &= (\Gamma_g H'_\mu \Gamma_g^{-1}) \Gamma_g \psi_{i\alpha}^p = \sum_\nu H'_\nu D'_{\nu\mu}(g) \sum_j \psi_{j\alpha}^p D_{ji}^p(g) \\[1em] &= \sum_\nu \sum_j H'_\nu \psi_{j\alpha}^p D'_{\nu\mu}(g) D_{ji}^p(g) \\[1em] &= \sum_{\nu,j} H'_\nu \psi_{j\alpha}^p [D'(g) \otimes D^p(g)]_{\nu j, \mu i} \end{aligned}

这表示函数集 {Hμψiαp}\{H'_\mu \psi_{i\alpha}^p\} (μ,i\mu, i 变化) 形成了 群 GG 的直积表示 DDpD' \otimes D^p 的基。若要判断跃迁矩阵元 ψjβqHμψiαp\langle\psi_{j\beta}^q | H'_\mu | \psi_{i\alpha}^p\rangle 是否为 00,只需对表示 DDpD' \otimes D^p 进行约化,应用矩阵元定理,判断其中是否含有不可约表示 DqD^q 即可

DDp=qaqDqaq=1GgGχq(g)χp(g)χ(g)D' \otimes D^p = \sum_q \bigoplus a_q D^q \\[1em] a_q = \dfrac{1}{|G|} \sum_{g \in G} \chi^q(g)^* \chi^p(g) \chi'(g)

aq=0a_q = 0,则 HμψiαpH'_\mu \psi_{i\alpha}^p 中不包含属 DqD^q 的所有态 {ψjβq}\{\psi_{j\beta}^q\} 的分量,应用基函数的正交性定理,跃迁积分 =0= 0,跃迁禁戒

aq0a_q \neq 0,只能确定 HμψiαpH'_\mu \psi_{i\alpha}^p 中包含属 DqD^q 的某些态的分量,还需进一步确定 HμψiαpH'_\mu \psi_{i\alpha}^p 中是否包含 DqD^q 的第 jj 列基ψjβq\psi_{j\beta}^q 的分量,如有则跃迁允许,否则跃迁禁戒

电偶极跃迁的选择定则

电偶极矩算符 P=er=(ex,ey,ez)\boldsymbol{P} = e\boldsymbol{r} = (ex, ey, ez) 是一个矢量算符。在光与物质相互作用中,当光场较弱时,微扰项 H^\hat{H}' 可近似正比于电偶极矩,P\boldsymbol{P} 的三个分量 Px,Py,PzP_x, P_y, P_z 代表三个方向上的偏振,分别对应于 H^\hat{H}' 的三个方向上的分量。

电偶极跃迁C4v

C4vC_{4v} 群为例,前两列为不同幂次的基函数。表示约化为 D=A1ED' = A_1 \oplus E,不可约表示与各分量的对应关系为 A1:Pz,E:Px,PyA_1: P_z, \quad E: P_x, P_y

  1. A1Di=Di (i=1,,5)\textcolor{red}{A_1} \otimes D^i = D^i \ (i = 1, \cdots, 5),此时微扰表示为 A1A_1,初态表示为 DiD^iDDpD' \otimes D^p 只有 DiD^i 分量。即 zz 方向偏振光 Pz\textcolor{red}{P_z} 作用下,只有自身与自身 DiDiD^i \leftrightarrow D^i 之间跃迁是允许的,其他跃迁是禁戒的

  2. EDi=E (i=1,2,3,4), EE=A1A2B1B2\textcolor{blue}{E} \otimes D^i = E \ (i = 1,2,3,4), \ E \otimes E = A_1 \oplus A_2 \oplus B_1 \oplus B_2。也就是说,在 Px,Py\textcolor{blue}{P_x, P_y} 作用下, E{A1,A2,B1,B2}E \leftrightarrow \{A_1, A_2, B_1, B_2\} 间跃迁是允许的,其他跃迁则是禁戒的

  • 直积表示的特征标就是特征标直接相乘

磁偶极跃迁的选择定则

此时微扰项 H^\hat{H}' 可近似正比于磁偶极矩算符

m=12r×J(r)dV\boldsymbol{m} = \frac{1}{2} \int \boldsymbol{r} \times \boldsymbol{J}(\boldsymbol{r}) dV

这是一个轴矢量算符变换,其性质如同轴矢量 (Rx,Ry,Rz)(R_x, R_y, R_z)。我们仍以 C4vC_{4v} 群为例。微扰项表示约化为 D=A2ED' = \textcolor{red}{A_2} \oplus \textcolor{blue}{E}, 不可约表示与各分量的对应关系为 A2:mz,E:mx,myA_2: \textcolor{red}{m_z}, \quad E: \textcolor{blue}{m_x, m_y}

磁偶极跃迁C4v
  1. A2A1=A2, A2A2=A1, A2B1=B2, A2B2=B1, A2E=E\textcolor{red}{A_2} \otimes A_1 = A_2, \ \textcolor{red}{A_2} \otimes A_2 = A_1, \ \textcolor{red}{A_2} \otimes B_1 = B_2, \ \textcolor{red}{A_2} \otimes B_2 = B_1, \ \textcolor{red}{A_2} \otimes E = E。这表明,在 mz\textcolor{red}{m_z} 作用下, A1A2A_1 \leftrightarrow A_2B1B2B_1 \leftrightarrow B_2EEE \leftrightarrow E 间跃迁是允许的,其他禁戒

  2. EDi=E (i=1,2,3,4), EE=A1A2B1B2\textcolor{blue}{E} \otimes D^i = E \ (i = 1,2,3,4), \ \textcolor{blue}{E} \otimes E = A_1 \oplus A_2 \oplus B_1 \oplus B_2。这表示,在 mx,my\textcolor{blue}{m_x, m_y} 作用下, E{A1,A2,B1,B2}E \leftrightarrow \{A_1, A_2, B_1, B_2\} 间跃迁是允许的,其他跃迁则是禁戒的,这与电偶极矩跃迁相同。

中心反演对称体系的跃迁

不可约表示分为奇宇称 {Dui}\{D_u^i\} 和偶宇称 {Dgi}\{D_g^i\} 两类,基函数满足奇宇称: ψ(r)=ψ(r)\psi(-\boldsymbol{r}) = -\psi(\boldsymbol{r}) \quad 偶宇称: ψ(r)=ψ(r)\psi(-\boldsymbol{r}) = \psi(\boldsymbol{r})

若微扰算符 HH' 具有奇宇称,比如电偶极跃迁(正比于位置矢量),若初末态为相同宇称,则跃迁矩阵元为

==0,==0\langle\text{奇}|\text{奇}|\text{奇}\rangle = \langle\text{奇}|\text{偶}\rangle = 0, \quad \langle\text{偶}|\text{奇}|\text{偶}\rangle = \langle\text{偶}|\text{奇}\rangle = 0

相同宇称间跃迁禁戒,不同宇称间跃迁允许。

若微扰算符 HH' 具有偶宇称,比如磁偶极跃迁(轴矢量空间反演不变号),若初末态为不同宇称,则跃迁矩阵元为

==0,==0\langle\text{奇}|\text{偶}|\text{偶}\rangle = \langle\text{奇}|\text{偶}\rangle = 0, \quad \langle\text{偶}|\text{偶}|\text{奇}\rangle = \langle\text{偶}|\text{奇}\rangle = 0

不同宇称间跃迁禁戒,相同宇称间跃迁允许。

原子振动对跃迁的影响

有时候我们不仅需要考虑电子的波函数,还需要考虑声子的波函数。

我们知道 dd 态电子波函数都是偶宇称 (x2y2,z2,xy,xz,yz)(\propto x^2 - y^2, z^2, xy, xz, yz),那么在中心反演对称体系中 ddd \leftrightarrow d 电子态间的电偶极跃迁就应该是禁戒的。然而实验上却能观察到这种跃迁,比如在过渡金属离子中。原因在于实际中存在热振动的参与,即原子振动与电子态的耦合

ψ=ψaψe \psi = \psi_a\psi_e \

其中 ψa\psi_a 为原子振动的波函数。应判断 ψaψeHψaψe\langle\psi'_a\psi'_e|H'|\psi_a\psi_e\rangle 是否为 0,而不是仅判断 ψeHψe\langle\psi'_e|H'|\psi_e\rangle。若存在 ψa\psi_a 为偶宇称而 ψa\psi'_a 为奇宇称,虽然电子部分是禁戒的

dPid==0\langle d'|P_i|d\rangle = \langle\text{偶}|\text{奇}|\text{偶}\rangle = 0

但考虑声子部分后仍有

ψadPiψad=奇偶偶偶=\langle\psi'_a d'|P_i|\psi_a d\rangle = \langle\text{奇偶}|\text{奇}|\text{偶偶}\rangle = \langle\text{奇}|\text{奇}|\text{偶}\rangle

跃迁矩阵元可以不为 00

红外吸收和拉曼跃迁

电磁辐射不太强时与物质的相互作用可近似为电偶极微扰。红外吸收便是如此, HP=(Px,Py,Pz)(x,y,z)H' \propto \boldsymbol{P} = (P_x, P_y, P_z) \propto (x, y, z),具有中心反演体系中仅不同宇称态间存在红外吸收。

拉曼效应则是一束近单色光入射到物质上与其发生作用,散射光的频率发生变化的效应

ωsc=ωin±ωba,ωba=(EbEa)/\omega_{sc} = \omega_{in} \pm \omega_{ba}, \quad \omega_{ba} = (E_b - E_a)/\hbar

详细研究表明,仅当初态 a|a\rangle,末态 b|b\rangle 及某些中间态 i|i\rangle 间存在偶极矩矩阵元时拉曼跃迁才可能发生

briirabrra0\langle b|\boldsymbol{r}|i\rangle\langle i|\boldsymbol{r}|a\rangle \sim \langle b|\boldsymbol{r}\boldsymbol{r}|a\rangle \neq 0

拉曼跃迁算符 HRrirjH'_R \propto r_i r_j,因而具有偶宇称,中心反演对称体系中仅允许相同宇称态间发生拉曼跃迁。

这表明,中心反演对称体系中,红外吸收仅允许不同宇称态间发生跃迁,而拉曼跃迁仅允许相同宇称态间发生跃迁。

排斥定理:对于中心反演对称的体系,相同振动频率不可能同时是拉曼激活和红外激活的。

不可约张量算符和 Wigner-Eckart 定理

矢量算符

如果一个算符有三个分量 {A^1,A^2,A^3}\{\hat{A}_1, \hat{A}_2, \hat{A}_3\},且在旋转 RR 的作用下如下变换

P^RA^iP^R1=j=13A^jRjiP^RA^P^R1=R1A^\hat{P}_R \hat{A}_i \hat{P}_R^{-1} = \sum_{j=1}^3 \hat{A}_j R_{ji} \quad \text{或} \quad \hat{P}_R \hat{\boldsymbol{A}} \hat{P}_R^{-1} = R^{-1} \hat{\boldsymbol{A}}

则称该算符为矢量算符,记为 A^=(A^1,A^2,A^3)T\hat{\boldsymbol{A}} = (\hat{A}_1, \hat{A}_2, \hat{A}_3)^T,因为它与位置算符 r^\hat{\boldsymbol{r}} 有着相同的变换性质,或说它们构成了旋转 RR 自身表示 D(R)=RD(R) = R 的基。

算符为行矢量和列矢量时的表达式分别如下

P^R(A^1,A^2,A^3)P^R1=(A^1,A^2,A^3)RR 正交:RT=R1P^R(A^1A^2A^3)P^R1=R1(A^1A^2A^3)\hat{P}_R (\hat{A}_1, \hat{A}_2, \hat{A}_3) \hat{P}_R^{-1} = (\hat{A}_1, \hat{A}_2, \hat{A}_3)R \quad \xLeftrightarrow{R \text{ 正交:} R^T = R^{-1}} \quad \hat{P}_R \begin{pmatrix} \hat{A}_1 \\[0.5em] \hat{A}_2 \\[0.5em] \hat{A}_3 \end{pmatrix} \hat{P}_R^{-1} = R^{-1} \begin{pmatrix} \hat{A}_1 \\[0.5em] \hat{A}_2 \\[0.5em] \hat{A}_3 \end{pmatrix}

算符 P^R\hat{P}_R 不受转置的作用

我们还只是考虑在正当旋转下的变换,若考虑空间反演 ii,则可以对矢量算符进一步区分。在空间反演下变号,真矢量算符:比如位置算符 r^\hat{\boldsymbol{r}},动量算符 p^\hat{\boldsymbol{p}};在空间反演下不变号,赝/轴矢量算符:角动量算符 L^\hat{\boldsymbol{L}},自旋算符 S^\hat{\boldsymbol{S}}

需要注意的是:矢量算符的变换与矢量的变换不同!

V^RV^=P^RV^P^R1=R1V^vRv=Rv\hat{\boldsymbol{V}} \xrightarrow{R} \hat{\boldsymbol{V}}' = \hat{P}_R \hat{\boldsymbol{V}} \hat{P}_R^{-1} = R^{-1}\hat{\boldsymbol{V}} \quad \quad \quad \boldsymbol{v} \xrightarrow{R} \boldsymbol{v}' = R\boldsymbol{v}

算符的每一分量被看成是位置 r\boldsymbol{r} 的函数,函数的旋转引入 R1R^{-1},以 r^i\hat{r}_i 算符为例,其作用在任意位置函数上为: r^iψ(r)=riψ(r)fi(r)\hat{r}_i\psi(\boldsymbol{r}) = r_i\psi(\boldsymbol{r}) \equiv f_i(\boldsymbol{r})。对这个体系做 RR 旋转,一方面有 fi(r)P^Rfi(r)f_i(\boldsymbol{r}) \rightarrow \hat{P}_R f_i(\boldsymbol{r})

P^Rfi(r)=fi(R1r)=(R1r)iψ(R1r)=(R1r^)iψ(R1r)\hat{P}_R f_i(\boldsymbol{r}) = f_i(R^{-1}\boldsymbol{r}) = (R^{-1}\boldsymbol{r})_i \psi(R^{-1}\boldsymbol{r}) = \textcolor{red}{(R^{-1}\hat{\boldsymbol{r}})_i} \psi(R^{-1}\boldsymbol{r})

另一方面有

P^Rfi(r)=P^Rr^iψ(r)=P^Rr^iP^R1P^Rψ(r)=(P^Rr^iP^R1)ψ(R1r)\hat{P}_R f_i(\boldsymbol{r}) = \hat{P}_R \hat{r}_i \psi(\boldsymbol{r}) = \hat{P}_R \hat{r}_i \hat{P}_R^{-1} \hat{P}_R \psi(\boldsymbol{r}) = \textcolor{red}{(\hat{P}_R \hat{r}_i \hat{P}_R^{-1})} \psi(R^{-1}\boldsymbol{r})

ψ(r)\psi(\boldsymbol{r}) 的任意性有

P^Rr^iP^R1=(R1r^)iP^Rr^P^R1=R1r^\hat{P}_R \hat{r}_i \hat{P}_R^{-1} = (R^{-1}\hat{\boldsymbol{r}})_i \quad \Rightarrow \quad \textcolor{red}{\hat{P}_R \hat{\boldsymbol{r}} \hat{P}_R^{-1} = R^{-1}\hat{\boldsymbol{r}}}

注意其中 r^\hat{\boldsymbol{r}} 的含义是其三个分量算符组成的列矩阵 r^(r^1r^2r^3)\hat{\boldsymbol{r}} \equiv \begin{pmatrix} \hat{r}_1 \\ \hat{r}_2 \\ \hat{r}_3 \end{pmatrix}。因矢量是一阶张量,故矢量算符也是一阶张量算符。有一阶张量算符自然也有更高阶的张量算符

张量算符

如果一个算符有 3r3^r 个分量 {A^i1i2iri1,,ir=1,2,3}\{\hat{A}_{i_1 i_2 \cdots i_r} | i_1, \cdots, i_r = 1, 2, 3\},且在旋转 RR 的作用下如下变换

P^RA^i1irP^R1=j1jr=13A^j1jrRj1i1Rj2i2Rjrir=j1jr=13A^j1jr(RRr 个)j1jr,i1ir\begin{aligned} \hat{P}_R \hat{A}_{i_1 \cdots i_r} \hat{P}_R^{-1} &= \sum_{j_1 \cdots j_r=1}^3 \hat{A}_{j_1 \cdots j_r} R_{j_1 i_1} R_{j_2 i_2} \cdots R_{j_r i_r} \\ &= \sum_{j_1 \cdots j_r=1}^3 \hat{A}_{j_1 \cdots j_r} (\underbrace{R \otimes \cdots \otimes R}_{r \text{ 个}})_{j_1 \cdots j_r, i_1 \cdots i_r} \end{aligned}

则称该算符为 rr 阶(或秩)张量算符

这说明 rr 阶张量算符荷载着旋转群的直积表示: RRr 个\underbrace{R \otimes \cdots \otimes R}_{r \text{ 个}}

我们同样要注意 rr 阶张量算符与 rr 阶张量 TT 的定义的区别:

Ti1irRTi1ir=j1jr=13Ri1j1Ri2j2RirjrTj1jrT_{i_1 \cdots i_r} \quad \xrightarrow{R}\quad T'_{i_1 \cdots i_r} = \sum_{j_1 \cdots j_r=1}^3 R_{i_1 j_1} R_{i_2 j_2} \cdots R_{i_r j_r} T_{j_1 \cdots j_r}

矢量算符的右边是对行指标求和,而张量的右边是对列指标求和。


设有二阶张量算符 {T^ij}\{\hat{T}_{ij}\},以它为基可以荷载 SO(3)\text{SO}(3) 群的一个 99 维表示 RRR \otimes R,其与 D(1)D(1)D^{(1)} \otimes D^{(1)} 等价,但是一个可约表示,在约化之前需要计算其特征标

χ(1)(1)(α)=[sin3α2/sinα2]2=3+4cosα+2cos2α\chi^{(1)\otimes(1)}(\alpha) = \left[\sin\frac{3\alpha}{2} / \sin\frac{\alpha}{2}\right]^2 = 3 + 4\cos\alpha + 2\cos2\alpha

经过计算可以得到

χ(1)(1)(α)=χ(0)(α)+χ(1)(α)+χ(2)(α)\chi^{(1)\otimes(1)}(\alpha) = \chi^{(0)}(\alpha) + \chi^{(1)}(\alpha) + \chi^{(2)}(\alpha)

D(1)D(1)=D(0)D(1)D(2)\Rightarrow \quad D^{(1)} \otimes D^{(1)} = D^{(0)} \oplus D^{(1)} \oplus D^{(2)}

可以对张量算符 {T^ij}\{\hat{T}_{ij}\} 重新线性组合,使得它们荷载的是 SO(3)\text{SO}(3) 的不可约表示,新组合得到的 99 个基为

A^=13(T^11+T^22+T^33)B^k=12i,jεijkT^ijC^ij=12(T^ij+T^ji)A^δij\begin{aligned} \hat{A} &= \frac{1}{3}(\hat{T}_{11} + \hat{T}_{22} + \hat{T}_{33}) \\[1em] \hat{B}_k &= \frac{1}{2}\sum_{i,j} \varepsilon_{ijk}\hat{T}_{ij} \\[1em] \hat{C}_{ij} &= \frac{1}{2}(\hat{T}_{ij} + \hat{T}_{ji}) - \hat{A}\delta_{ij} \end{aligned}

那么张量算符的 99 个分量 {T^ij}\{\hat{T}_{ij}\} 可以拆解为

T^=(T^11T^12T^13T^21T^22T^23T^31T^32T^33)=(A^000A^000A^)+(0B^3B^2B^30B^1B^2B^10)+(C^11C^12C^13C^12C^22C^23C^13C^23(C^11+C^22))\hat{T} = \begin{pmatrix} \hat{T}_{11} & \hat{T}_{12} & \hat{T}_{13} \\[0.5em] \hat{T}_{21} & \hat{T}_{22} & \hat{T}_{23} \\[0.5em] \hat{T}_{31} & \hat{T}_{32} & \hat{T}_{33} \end{pmatrix} = \begin{pmatrix} \hat{A} & 0 & 0 \\[0.5em] 0 & \hat{A} & 0 \\[0.5em] 0 & 0 & \hat{A} \end{pmatrix} + \begin{pmatrix} 0 & \hat{B}_3 & -\hat{B}_2 \\[0.5em] -\hat{B}_3 & 0 & \hat{B}_1 \\[0.5em] \hat{B}_2 & -\hat{B}_1 & 0 \end{pmatrix} + \begin{pmatrix} \hat{C}_{11} & \hat{C}_{12} & \hat{C}_{13} \\[0.5em] \hat{C}_{12} & \hat{C}_{22} & \hat{C}_{23} \\[0.5em] \hat{C}_{13} & \hat{C}_{23} & -(\hat{C}_{11} + \hat{C}_{22}) \end{pmatrix}

A^\hat{A} 是单位矩阵常数倍,荷载恒等表示 D(0)D^{(0)}{B^k}\{\hat{B}_k\} 是反对称矩阵,荷载 RR 表示,与 D(1)D^{(1)} 等价。{C^ij}\{\hat{C}_{ij}\} 是无迹矩阵,荷载一个与 D(2)D^{(2)} 等价的表示。

原来可约的直积表示 RRR \otimes R 变换的 99 个分量按变换性质的不同,由简单到复杂被分成了三类不可约表示的基:

  • 一个 A^\hat{A} 按恒等表示变换,即不变
  • 三个 B^k\hat{B}_kRR 变换,如同一个矢量
  • 五个 C^ij\hat{C}_{ij} 按一个五维不可约表示变换,是一个无迹对称张量

这些按不可约表示变换的算符便是不可约张量算符。但以上只是我们凑出来的不可约表示的算符基,不可约张量算符可通过 SO(3)\text{SO}(3) 群的旋转严格定义

定义一组 2l+12l+1 个算符 {T^il}, i=l,l1,,l\{\hat{T}_i^l\}, \ i = l, l-1, \cdots, -l,如果它们在旋转 RR 下依 SO(3) 群的不可约表示 D(l)D^{(l)} 变换的话,即

P^RT^ilP^R1=j=llT^jlDji(l)(R)\hat{P}_R \hat{T}_i^l \hat{P}_R^{-1} = \sum_{j=l}^{-l} \hat{T}_j^l D_{ji}^{(l)}(R)

这组算符被称为 ll 阶的(球面)不可约张量算符,可记为 T^l={T^il}\hat{T}^l = \{\hat{T}_i^l\}

例:普通矢量算符 {V^x,V^y,V^z}\{\hat{V}_x, \hat{V}_y, \hat{V}_z\} 荷载的表示为 RR,只能称为矢量算符而不能称为(球面)不可约张量算符,因为其不可约表示不是 D(l)D^{(l)}。其写成 11 阶不可约张量算符为:

T^±11=12V^±=12(V^x±iV^y)T^01=V^z\hat{T}_{\pm 1}^1 = \mp \frac{1}{\sqrt{2}}\hat{V}_{\pm} = \mp \frac{1}{\sqrt{2}} (\hat{V}_x \pm i\hat{V}_y) \quad \quad \hat{T}_0^1 = \hat{V}_z

更一般的,可定义一组依群 GG 的某个不可约表示 DpD^p 变换的算符 {J^μpμ=1,,dp}\{\hat{J}_\mu^p | \mu = 1, \cdots, d_p\}不可约张量算符,满足

Γ^gJ^μpΓ^g1=ν=1dpJ^νpDνμp(g)\hat{\Gamma}_g \hat{J}_\mu^p \hat{\Gamma}_g^{-1} = \sum_{\nu=1}^{d_p} \hat{J}_\nu^p D_{\nu\mu}^p(g)

荷载了某个群的不可约表示就是不可约张量算符,荷载了可约表示(例如 RRR \otimes R)的则只是张量算符

若一组波函数 {ϕiq}\{\phi_i^q\} 是群 GG 的不可约表示 DqD^q 的基,且不可约张量算符 {J^jp}\{\hat{J}_j^p\} 也能荷载群 GG 的不可约表示 DpD^p,那么 {J^jp}\{\hat{J}_j^p\} 作用在这些态上后依然可荷载 GG 的一个表示

Γ^gJ^jpϕiq=(Γ^gJ^jpΓ^g1)Γ^gϕiq=jJ^jpDjjp(g)iϕiqDiiq(g)=j,iJ^jpϕiqDjjp(g)Diiq(g)=j,i(J^jpϕiq)[Dp(g)Dq(g)]ji,ji\begin{aligned} \hat{\Gamma}_g \hat{J}_j^p \phi_i^q &= (\hat{\Gamma}_g \hat{J}_j^p \hat{\Gamma}_g^{-1}) \hat{\Gamma}_g \phi_i^q = \sum_{j'} \hat{J}_{j'}^p D_{j'j}^p(g) \sum_{i'} \phi_{i'}^q D_{i'i}^q(g) \\ &= \sum_{j',i'} \hat{J}_{j'}^p \phi_{i'}^q D_{j'j}^p(g) D_{i'i}^q(g) = \sum_{j',i'} (\hat{J}_{j'}^p \phi_{i'}^q) [D^p(g) \otimes D^q(g)]_{j'i', ji} \end{aligned}

这表明 {J^jpϕiq}\{\hat{J}_j^p \phi_i^q\} 荷载了群 GG 的直积表示 DpDqD^p \otimes D^q,它一般是可约的。设约化后的基为: {ψlαtt=1,,r;l=1,,dt;α=1,,at}\{\psi_{l\alpha}^t | t=1,\cdots,r; l=1,\cdots,d_t; \alpha=1,\cdots,a_t\}

约化前后两组基之间的变换关系正好就是我们之前讲到的 CG\text{CG} 系数

ψlαt=j,ipqjitlαJ^jpϕiqJ^jpϕiq=t,l,αpqjitlαψlαt\psi_{l\alpha}^t = \sum_{j,i} \langle \begin{matrix} p & q \\ j & i \end{matrix} | \begin{matrix} t \\ l \end{matrix} \alpha \rangle \hat{J}_j^p \phi_i^q \quad \quad \hat{J}_j^p \phi_i^q = \sum_{t,l,\alpha} \langle \begin{matrix} p & q \\ j & i \end{matrix} | \begin{matrix} t \\ l \end{matrix} \alpha \rangle^* \psi_{l\alpha}^t

CG\text{CG} 系数就是同一个群的 两个不可约表示 的 直积表示 约化到一个 不可约表示 时基函数的变换矩阵的矩阵元

然后来看 J^jp\hat{J}_j^p 作为微扰时,从初态 ϕiq\phi_i^q 到末态 ϕks\phi_k^s 的跃迁矩阵元

ϕksJ^jpϕiq=(ϕks,J^jpϕiq)\langle \phi_k^s | \hat{J}_j^p | \phi_i^q \rangle = (\phi_k^s, \hat{J}_j^p \phi_i^q)

借助 CG\text{CG} 系数