この記事は現在、プロジェクトメンバーによる査読中のものです。草稿段階ですので、内容・表現の正確さについて責任を負いかねます。 リンクを正しく張れていないところが存在しますのでご注意ください。 正式公開まで、いましばらくお待ちください。
交換相互作用と任意の方向を向いた二電子スピン

交換相互作用は, \bm{S}_1 \cdot \bm{S}_2 を 用いて, - \dfrac{J_{ab}}{2}(\hbar^2 + 4 \bm{S}_1 \cdot \bm{S}_2 ) と表せられます. この記事では,\bm{S}_1 \cdot \bm{S}_2 \neq \dfrac{\hbar^2}{4},-\dfrac{3 \hbar^2}{4} の時は, どうなるんだろう?という問いに答えます.

一電子の固有スピノール

我々は, 任意の方向を向いたスピンのxyz方向固有状態での展開 において, 一電子スピンの任意の方向を向いた状態の記述方法を学びました. それを復習しておきましょう.極座標での (\theta , \phi) 方向を向いたスピンは, その任意方向単位ベクトル

\bm{n}=\begin{pmatrix}\sin \theta \cos \phi \\ \sin \theta \sin \phi \\ \cos \theta \end{pmatrix} \tag{1}

として,また,パウリのスピン行列

\sigma_x = \begin{pmatrix} 0 &  1 \\ 1 &  0  \end{pmatrix}, \ \ \ \sigma_y = \begin{pmatrix} 0 & -i \\ i &  0  \end{pmatrix}, \ \ \ \sigma_z = \begin{pmatrix} 1 &  0 \\ 0 & -1  \end{pmatrix} \tag{2} \\\bm{\sigma} = (\sigma_x,\sigma_y,\sigma_z)\tag{3}

として,そして,スピノール

\chi = \begin{pmatrix}c_+ \\c_-\end{pmatrix} \tag{4}

として,次の固有値方程式を満たすものでした.

(\bm{n} \cdot \bm{\sigma}) \chi &= \left( \sum_{i=1}^3 n_i \sigma_i \right) \chi \\&= \begin{pmatrix} \cos \theta & \sin \theta e^{-i \phi} \\ \sin \theta e^{i \phi} & - \cos \theta  \end{pmatrix} \begin{pmatrix} c_+ \\ c_- \end{pmatrix} \\ &= \lambda  \begin{pmatrix} c_+ \\ c_- \end{pmatrix} \tag{5}

この固有値方程式は解が二つあり,それぞれ,正の方向と負の方向の解 \chi_+,\chi_- が出てきます. これらは固有値 \lambda=\pm1 を持ち,固有ベクトル(固有スピノール)を書くと,

\chi_+ &= \begin{pmatrix}  \cos \dfrac{\theta}{2} \\ \sin \dfrac{\theta}{2} e^{i \phi}  \end{pmatrix} \tag{6} \\\chi_- &= \begin{pmatrix} -\sin \dfrac{\theta}{2} \\ \cos \dfrac{\theta}{2} e^{i \phi}  \end{pmatrix} \tag{7}

となります.

二電子の交換相互作用

ここで,少し話が変わりますが,二電子のスピンと軌道の状態を記述する波動関数を確認しておきます. 隣接する二つの原子を考え,簡単のためそれぞれの電子軌道を \phi_a(\bm{r}),\phi_b(\bm{r}), とし, それらは局在していて,互いに直交しているとします.そして, \phi_a,\phi_b にそれぞれ一個ずつ電子が入っているとします. スピン三重項を t+,t0,t- で表し,一重項を s で表すと,

\Psi_{t+} &= \dfrac{1}{\sqrt{2}}(\phi_a(\bm{r}_1)\phi_b(\bm{r}_2)-\phi_b(\bm{r}_1)\phi_a(\bm{r}_2))\alpha(\sigma_1)\alpha(\sigma_2) \tag{8} \\\Psi_{t0} &= \dfrac{1}{2}(\phi_a(\bm{r}_1)\phi_b(\bm{r}_2)-\phi_b(\bm{r}_1)\phi_a(\bm{r}_2))(\alpha(\sigma_1)\beta(\sigma_2)+\beta(\sigma_1)\alpha(\sigma_2)) \tag{9} \\\Psi_{t-} &= \dfrac{1}{\sqrt{2}}(\phi_a(\bm{r}_1)\phi_b(\bm{r}_2)-\phi_b(\bm{r}_1)\phi_a(\bm{r}_2))\beta(\sigma_1)\beta(\sigma_2) \tag{10} \\\Psi_{s} &= \dfrac{1}{2}(\phi_a(\bm{r}_1)\phi_b(\bm{r}_2)+\phi_b(\bm{r}_1)\phi_a(\bm{r}_2))(\alpha(\sigma_1)\beta(\sigma_2)-\beta(\sigma_1)\alpha(\sigma_2)) \tag{11}

となります.ただし,ここでスピン座標 \sigma_1,\sigma_2 は前述したパウリ行列とは関係なくて,軌道 \phi_a,\phi_b にたいする座標 \bm{r}_1\bm{r}_2 と同じようにスピンに対する離散的な座標であり, \sigma_1,\sigma_21/2,-1/2 いずれかをとります.すこし紛らわしいですが,この表記が慣例となっているようです.二電子の波動関数は軌道とスピン合わせて全体として反対称ですから,三重項 (t) では,軌道部分が反対称でスピン部分が対称, 一重項 (s) では軌道部分が対称でスピン部分が反対称ですね.ここで,通常の原子の二電子ハミルトニアンは,

\mathcal{H}_tot =  \dfrac{-\hbar^2}{2m}\bm{\nabla}_1^2 + \dfrac{-\hbar^2}{2m}\bm{\nabla}_2^2 + \dfrac{1}{4 \pi \varepsilon_0}\dfrac{e^2}{r_{12}}

となります.ここで, \bm{\nabla}_1 は一番目の電子に作用することを意味します.二番目の電子についても同様です. さて,二電子スピンの三重項と一重項では,運動エネルギーは変わりません.これを示しましょう.例えば,電子1の運動エネルギーは,二電子とも上向き(αスピン)の時,ラプラシアン \triangle_1 は電子1のみに働くことを考えれば,

\langle \Psi_{t+} | \dfrac{p_1^2}{2m} | \Psi_{t+}\rangle &= \int \Psi^\ast_{t+} \dfrac{-\hbar^2 \triangle_1 }{2m} \Psi_{t+} d \tau \\&= \dfrac{1}{2} \int ( \phi_a(\bm{r}_1)^\ast \phi_b(\bm{r}_2)^\ast - \phi_b(\bm{r}_1)^\ast \phi_a(\bm{r}_2)^\ast ) \dfrac{-\hbar^2}{2m} \triangle_1 ( \phi_a(\bm{r}_1) \phi_b(\bm{r}_2) - \phi_b(\bm{r}_1) \phi_a(\bm{r}_2) ) d\bm{r}_1d\bm{r}_2 \langle \alpha_1 \alpha_2 | \alpha_1 \alpha_2 \rangle \\&= \dfrac{1}{2} \int \phi_a(\bm{r}_1)^\ast ( \triangle_1 \phi_a(\bm{r}_1) ) d\bm{r}_1 \int \phi_b(\bm{r}_2)^\ast \phi_b(\bm{r}_2) d\bm{r}_2 \\&- \dfrac{1}{2} \int \phi_a(\bm{r}_1)^\ast ( \triangle_1 \phi_b(\bm{r}_1) ) d\bm{r}_1 \int \phi_b(\bm{r}_2)^\ast \phi_a(\bm{r}_2) d\bm{r}_2 \\&- \dfrac{1}{2} \int \phi_b(\bm{r}_1)^\ast ( \triangle_1 \phi_a(\bm{r}_1) ) d\bm{r}_1 \int \phi_a(\bm{r}_2)^\ast \phi_b(\bm{r}_2) d\bm{r}_2 \\&= \dfrac{1}{2} \int \phi_b(\bm{r}_1)^\ast ( \triangle_1 \phi_b(\bm{r}_1) ) d\bm{r}_1 \int \phi_a(\bm{r}_2)^\ast \phi_a(\bm{r}_2) d\bm{r}_2 \\&= \dfrac{1}{2} \langle \phi_{a1} | \triangle_1 | \phi_{a1} \rangle +\dfrac{1}{2} \langle \phi_{b1} | \triangle_1 | \phi_{b1} \rangle

となります. ただし,

\langle \alpha_1 \alpha_2 | \alpha_1 \alpha_2 \rangle &= 1 \\\int \phi_b(\bm{r}_2)^\ast \phi_a(\bm{r}_2) d\bm{r}_2 &= 0 \\\int \phi_a(\bm{r}_2)^\ast \phi_b(\bm{r}_2) d\bm{r}_2 &= 0 \tag{12}

を使いました.よって,この式はスピンの依存性が無いことが分かります. しかし,そのスピンは,クーロン相互作用に顕著な違いをもって働きます.実際に,式 (8)(11) に対し,クーロン相互作用の期待値を取ると,三重項 (t) では, r_{12}=|\bm{r}_1-\bm{r}_2| として,

\langle \Psi_{t+} | \dfrac{e^2}{r_{12}} | \Psi_{t+}\rangle &= \int \Psi^\ast_{t+} \dfrac{e^2}{r_{12}} \Psi_{t+} d \tau \\&= \dfrac{1}{2} \int ( \phi_a(\bm{r}_1)^\ast \phi_b(\bm{r}_2)^\ast - \phi_b(\bm{r}_1)^\ast \phi_a(\bm{r}_2)^\ast ) \dfrac{e^2}{r_{12}} ( \phi_a(\bm{r}_1) \phi_b(\bm{r}_2) - \phi_b(\bm{r}_1) \phi_a(\bm{r}_2) ) d\bm{r}_1d\bm{r}_2 \langle \alpha_1 \alpha_2 | \alpha_1 \alpha_2 \rangle \\&= \dfrac{1}{2} \int \dfrac{e^2}{r_{12}} |\phi_a(\bm{r}_1)|^2 |\phi_b(\bm{r}_2)|^2 d \bm{r}_1 d \bm{r}_2 \\&+ \dfrac{1}{2} \int \dfrac{e^2}{r_{12}} |\phi_b(\bm{r}_1)|^2 |\phi_a(\bm{r}_2)|^2 d \bm{r}_1 d \bm{r}_2 \\&- \dfrac{1}{2} \int \dfrac{e^2}{r_{12}} \phi_a^\ast(\bm{r}_1) \phi_b(\bm{r}_1) \phi_b^\ast(\bm{r}_2) \phi_a(\bm{r}_2)  \bm{r}_1 d \bm{r}_2 \\&- \dfrac{1}{2} \int \dfrac{e^2}{r_{12}} \phi_b^\ast(\bm{r}_1) \phi_a(\bm{r}_1) \phi_a^\ast(\bm{r}_2) \phi_b(\bm{r}_2)  \bm{r}_1 d \bm{r}_2 \\&= K_{ab} - J_{ab} \tag{13}

となります.ただし,

K_{ab} &= \int \dfrac{e^2}{r_{12}} |\phi_a(\bm{r}_1)|^2 |\phi_b(\bm{r}_2)|^2 d \bm{r}_1 d \bm{r}_2 \tag{14} \\J_{ab} &= \int \dfrac{e^2}{r_{12}} \phi_a^\ast(\bm{r}_1) \phi_b(\bm{r}_1) \phi_b^\ast(\bm{r}_2) \phi_a(\bm{r}_2)  \bm{r}_1 d \bm{r}_2 \tag{15}

です. 一重項 (s) で同様の計算を行うと,

\langle \Psi_{s} | \dfrac{e^2}{r_{12}} | \Psi_{s}\rangle = K_{ab} + J_{ab} \tag{16}

となります. 2J_{ab} だけ,三重項と一重項ではエネルギーが異なります.ここでは 示しませんが, J_{ab} \geq 0 が示せます.参考文献として,下に斯波先生の本を挙げておきます.

そして,この結果を少々技巧的な手法によって,変形します.右に来る関数の二電子を交換する演算子 P_{ex} というものを考えるのです,これにより,

P_{ex}| \alpha_1 \alpha_2 \rangle &= | \alpha_1 \alpha_2 \rangle \tag{17} \\P_{ex} \dfrac{1}{\sqrt{2}} \left( | \alpha_1 \beta_2 \rangle + | \beta_1 \alpha_2 \rangle  \right) &= \dfrac{1}{\sqrt{2}} \left( | \alpha_1 \beta_2 \rangle + | \beta_1 \alpha_2 \rangle  \right) \tag{18} \\P_{ex}| \beta_1 \beta_2 \rangle &= | \beta_1 \beta_2 \rangle \tag{19} \\P_{ex} \dfrac{1}{\sqrt{2}} \left( | \alpha_1 \beta_2 \rangle - | \beta_1 \alpha_2 \rangle  \right) &= - \dfrac{1}{\sqrt{2}} \left( | \alpha_1 \beta_2 \rangle - | \beta_1 \alpha_2 \rangle  \right) \tag{20}

となります.すると,クーロンエネルギー E_{t},E_{s} は,一重項か三重項かを問わず, K_{ab}-J_{ab} P_{ex} をスピン部分で挟むことにより,

E_{t+} &= \langle \alpha_1 \alpha_2 | K_{ab}-J_{ab}P_{ex} | \alpha_1 \alpha_2 \rangle \\&= K_{ab}-J_{ab} \tag{21} \\E_{t0} &=\dfrac{1}{2} (\langle \alpha_1 \beta_2 | + \langle \beta_1 \alpha_2) | K_{ab}-J_{ab}P_{ex} | ( \alpha_1 \beta_2 \rangle + |\beta_1 \alpha_2 \rangle) \\&= K_{ab}-J_{ab} \tag{22} \\E_{t-} &= \langle \beta_1 \beta_2 | K_{ab}-J_{ab}P_{ex} | \beta_1 \beta_2 \rangle \\&= K_{ab}-J_{ab} \tag{23} \\E_{s} &=\dfrac{1}{2} (\langle \alpha_1 \beta_2 | - \langle \beta_1 \alpha_2) | K_{ab}-J_{ab}P_{ex} | ( \alpha_1 \beta_2 \rangle - |\beta_1 \alpha_2 \rangle) \\&=\dfrac{1}{2} (\langle \alpha_1 \beta_2 | - \langle \beta_1 \alpha_2) | K_{ab} | ( \alpha_2 \beta_1 \rangle - |\beta_2 \alpha_1 \rangle) \\&- \dfrac{1}{2} (\langle \alpha_1 \beta_2 | - \langle \beta_1 \alpha_2) | J_{ab} | ( \alpha_2 \beta_1 \rangle - |\beta_2 \alpha_1 \rangle) \\&= K_{ab}+J_{ab} \tag{24}

で求まります.おそらく,この交換の演算子 P_{ex} が交換相互作用と呼ばれるゆえんではないかと,(筆者は勝手に)思っています. この辺の事情に詳しい方がいらっしゃったら,どうぞ教えてください.

交換相互作用についてよく使われる形

ここで,交換相互作用と言ったときによく用いられる形について言及しておきます. いきなりですが,スピン射影演算子は次のように右辺のパウリ行列で表されます.

S_x&=\dfrac{\hbar}{2} \sigma_x \\&= \dfrac{\hbar}{2} \begin{pmatrix}  0 & 1 \\1 & 0\end{pmatrix} \tag{25} \\S_y&=\dfrac{\hbar}{2} \sigma_y \\&=\dfrac{\hbar}{2} \begin{pmatrix}  0 & -i \\i & 0\end{pmatrix} \tag{26} \\S_z&=\dfrac{\hbar}{2} \sigma_z \\&=\dfrac{\hbar}{2} \begin{pmatrix}  1 & 0 \\0 & -1\end{pmatrix} \tag{27}

これらに電子 1 ,電子 2 を組み合わせて,テンソル積を構成します. それは 4 \times 4 の行列と同一視され,第一列〜第四列,第一行〜第四行は | \alpha_1 \alpha_2 \rangle, | \alpha_1 \beta_2 \rangle, | \beta_1 \alpha_2 \rangle, | \beta_1 \beta_2 \rangle を表します.

つまり,演算子を A とすると,

A = \begin{pmatrix}\langle \alpha_1 \alpha_2 |A| \alpha_1 \alpha_2 \rangle & \langle \alpha_1 \alpha_2 |A| \alpha_1 \beta_2 \rangle & \langle \alpha_1 \alpha_2 |A| \beta_1 \alpha_2 \rangle & \langle \alpha_1 \alpha_2 |A| \beta_1 \beta_2 \rangle \\ \langle \alpha_1 \beta_2 |A| \alpha_1 \alpha_2 \rangle & \langle \alpha_1 \beta_2 |A| \alpha_1 \beta_2 \rangle & \langle \alpha_1 \beta_2 |A| \beta_1 \alpha_2 \rangle & \langle \alpha_1 \beta_2 |A| \beta_1 \beta_2 \rangle \\ \langle \beta_1 \alpha_2 |A| \alpha_1 \alpha_2 \rangle & \langle \beta_1 \alpha_2 |A| \alpha_1 \beta_2 \rangle & \langle \beta_1 \alpha_2 |A| \beta_1 \alpha_2\ \rangle & \langle \beta_1 \alpha_2 |A| \beta_1 \beta_2 \rangle \\ \langle \beta_1 \beta_2 |A| \alpha_1 \alpha_2 \rangle & \langle \beta_1 \beta_2 |A| \alpha_1 \beta_2 \rangle & \langle \beta_1 \beta_2 |A| \beta_1 \alpha_2 \rangle & \langle \beta_1 \beta_2 |A| \beta_1 \beta_2 \rangle\end{pmatrix}

ということです.さて,具体的に演算子を考えていきましょう!

S_{1x} \otimes 1_{2} &= \dfrac{\hbar}{2} \sigma_{1x} \otimes 1_{2} \\&= \dfrac{\hbar}{2} \begin{pmatrix}0 & 1_2 \\ 1_2 & 0\end{pmatrix}&= \dfrac{\hbar}{2} \begin{pmatrix}0 & 0 & 1 & 0 \\0 & 0 & 0 & 1 \\1 & 0 & 0 & 0 \\0 & 1 & 0 & 0 \end{pmatrix}\tag{28} \\S_{1y} \otimes 1_{2}  &= \dfrac{\hbar}{2} \sigma_{1y} \otimes 1_{2} \\&= \dfrac{\hbar}{2} \begin{pmatrix}0 & -i \cdot 1_2 \\ i \cdot 1_2 & 0\end{pmatrix}&= \dfrac{\hbar}{2} \begin{pmatrix}0 & 0 & -i & 0 \\0 & 0 & 0 & -i \\i & 0 & 0 & 0 \\0 & i & 0 & 0 \end{pmatrix}\tag{29} \\S_{1z} \otimes 1_{2}  &= \dfrac{\hbar}{2} \sigma_{1z} \otimes 1_{2} \\&= \dfrac{\hbar}{2} \begin{pmatrix}1 \cdot 1_2 & 0 \\ 0 & -1 \cdot 1_2\end{pmatrix}&= \dfrac{\hbar}{2} \begin{pmatrix}1 & 0 & 0 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & -1 & 0 \\0 & 0 & 0 & -1 \end{pmatrix}\tag{30}

この規則はつまり,前の \sigma_i (電子1に対するスピン演算子)で 4 \times 4 行列を大ざっぱに 2 \times 2 行列に分割します.この時,その要素の値がその部分にかかります.そして,後ろの電子2に対する演算子, 1_2 等を掛けるのです. 1_2 とは,二次の単位行列を表したものです.これで,テンソル積の前の演算子は電子1にのみ作用し,後ろの演算子は電子2にのみ作用することになります.電子2についても書いておくと,

1_{1} \otimes S_{2x} &= \dfrac{\hbar}{2} 1_{1} \otimes \sigma_{2x} \\&= \dfrac{\hbar}{2} \begin{pmatrix}\sigma_x & 0 \\ 0 & \sigma_x\end{pmatrix}&= \dfrac{\hbar}{2} \begin{pmatrix}0 & 1 & 0 & 0 \\1 & 0 & 0 & 0 \\0 & 0 & 0 & 1 \\0 & 0 & 1 & 0 \end{pmatrix}\tag{31} \\1_{1} \otimes S_{2y} &= \dfrac{\hbar}{2} 1_{1} \otimes \sigma_{2y} \\&= \dfrac{\hbar}{2} \begin{pmatrix}\sigma_y & 0 \\0 & \sigma_y\end{pmatrix}&= \dfrac{\hbar}{2} \begin{pmatrix}0 & -i & 0 & 0 \\i & 0 & 0 & 0 \\0 & 0 & 0 & -i \\0 & 0 & i & 0 \end{pmatrix}\tag{32} \\1_{1} \otimes S_{2z} &= \dfrac{\hbar}{2} 1_{1} \otimes \sigma_{2z} \\&= \dfrac{\hbar}{2} \begin{pmatrix}\sigma_z & 0 \\ 0 & \sigma_z\end{pmatrix}&= \dfrac{\hbar}{2} \begin{pmatrix}1 & 0 & 0 & 0 \\0 & -1 & 0 & 0 \\0 & 0 & 1 & 0 \\0 & 0 & 0 & -1 \end{pmatrix}\tag{33}

そうしたら,次もいいでしょうか.

S_+ &= S_x  + i S_y \\ &=\hbar \begin{pmatrix}0 & 1 \\0 & 0\end{pmatrix}\tag{34} \\S_- &= S_x  - i S_y \\&=\hbar \begin{pmatrix}0 & 0 \\1 & 0\end{pmatrix}\tag{35}

で定義すると,

S_{1+} \otimes S_{2-} &= \begin{pmatrix}0 & S_{-} \\ 0 & 0\end{pmatrix}&=\hbar^2\begin{pmatrix}0 & 0 & 0 & 0 \\0 & 0 & 1 & 0 \\0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 \end{pmatrix}\tag{36} \\S_{1-} \otimes S_{2+} &= \begin{pmatrix}0 & 0 \\ S_{+} & 0\end{pmatrix}&=\hbar^2\begin{pmatrix}0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & 0 & 0 \end{pmatrix}\tag{37}

長いですが,もう少し頑張ってください.

S_{1z} \otimes S_{2z} &= \dfrac{\hbar^2}{4} \begin{pmatrix}\sigma_z & 0 \\ 0 & -\sigma_z\end{pmatrix} \\&=\dfrac{\hbar^2}{4}\begin{pmatrix}1 & 0 & 0 & 0 \\0 & -1 & 0 & 0 \\0 & 0 & -1 & 0 \\0 & 0 & 0 & 1 \end{pmatrix} \tag{38}

最後に,スピン射影演算子の内積です.

\bm{S}_1 \cdot \bm{S}_2 = \sum_{i=x}^z S_{1i} \otimes S_{2i} \tag{39}

もしくは,

\bm{S}_1 \cdot \bm{S}_2 = \dfrac{1}{2} \left( S_{1+} \otimes S_{2-} + S_{1+} \otimes S_{2-} \right) + S_{1z} \otimes S_{2z} \tag{40}

となります.どちらの式を使っても同じく,

\bm{S}_1 \cdot \bm{S}_2 &= \dfrac{\hbar^2}{4} \begin{pmatrix}1 & 0 & 0 & 0 \\0 & -1 & 2 & 0 \\0 & 2 & -1 & 0 \\0 & 0 & 0 & 1 \end{pmatrix} \tag{41}

となります.

少し寄り道

(41) は先ほどの交換演算子を用いて,美しく表すことができます.

交換演算子のスピンに対する行列は,

P_{ex} = \begin{pmatrix}1 & 0 & 0 & 0 \\0 & 0 & 1 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & 0 & 1 \end{pmatrix}

と書けるので,

2 \hbar^2 P_{ex} &= \hbar^2 + 4 (\bm{S}_1 \cdot \bm{S}_2) \\&= \hbar^2 + \hbar^2 (\bm{\sigma}_1 \cdot \bm{\sigma}_2) \tag{42}

つまり,

2 \hbar^2 \begin{pmatrix}1 & 0 & 0 & 0 \\0 & 0 & 1 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & 0 & 1 \end{pmatrix} &=\hbar^2\begin{pmatrix}1 & 0 & 0 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & 1 & 0 \\0 & 0 & 0 & 1 \end{pmatrix}+\hbar^2\begin{pmatrix}1 & 0 & 0 & 0 \\0 & -1 & 2 & 0 \\0 & 2 & -1 & 0 \\0 & 0 & 0 & 1 \end{pmatrix} \tag{43}

と書くことができます.これにより,交換相互作用は次の様に書くことができます.

K_{ab} - J_{ab} P_{ex} = K_{ab} - \dfrac{J_{ab}}{2}(1+ \dfrac{ 4 \bm{S}_1 \cdot \bm{S}_2 }{\hbar^2}) \tag{44}

これはよく定数部分を除いて,交換相互作用のハミルトニアン \mathcal{H}_{ex} として,

\mathcal{H}_{ex} &= - 2 J \bm{S}_1 \cdot \bm{S}_2 \tag{45}

と書かれます.そしてこれは,スピンが三重項の時は, \bm{S}_1 \cdot \bm{S}_2 = \dfrac{\hbar^2}{4} 一重項の時は \bm{S}_1 \cdot \bm{S}_2 = \dfrac{- 3 \hbar^2}{4} となります.これは式 (41) を対角化することによって,分かります. (式 (41) では二電子スピンが | \alpha_1 \alpha_2 \rangle ,| \alpha_1 \beta_2 \rangle ,| \beta_1 \alpha_2 \rangle ,| \beta_1 \beta_2 \rangle の順に並んでいるのでした.)

いよいよ本題

さて,ここでこれから進む方向を確認しておきます. 電子スピン1が方向 (\theta_1,\phi_1) を持ち,電子スピン2が方向 (\theta_2,\phi_2) を持つ状態を考えます. 少しばかり複雑な計算をして,その状態が, |\alpha_1 \alpha_2 \rangle , |\alpha_1 \beta_2 \rangle , |\beta_1 \alpha_2 \rangle , |\beta_1 \beta_2 \rangle の成分をどれだけ持っているかを調べます.

電子スピン1が方向 (\theta_1,\phi_1) を持ち,電子スピン2が方向 (\theta_2,\phi_2) を持つ状態はどうあらわされるでしょうか?それは,スピノールのテンソル積を考えればいいのです.正の方向を向いたスピノールだけを考えると,

つまり,

\chi_{1+} = \begin{pmatrix} \cos \dfrac{\theta_1}{2} \\ \sin \dfrac{\theta_1}{2} e^{i \phi_1} \end{pmatrix} \tag{46} \chi_{2+} = \begin{pmatrix} \cos \dfrac{\theta_2}{2} \\ \sin \dfrac{\theta_2}{2} e^{i \phi_2} \end{pmatrix} \tag{47}

簡単のため,前者を \begin{pmatrix} A \\ B \end{pmatrix} ,後者を \begin{pmatrix} C \\ D \end{pmatrix} とすると, スピノールのテンソル積は,

\chi_{1+} \otimes \chi_{2+} &= \begin{pmatrix} AC \\ AD \\ BC \\ BD \end{pmatrix} \\&= \begin{pmatrix} \cos \dfrac{\theta_1}{2} \cos \dfrac{\theta_2}{2} \\\cos \dfrac{\theta_1}{2} \sin \dfrac{\theta_2}{2} e^{i \phi_2} \\\sin \dfrac{\theta_1}{2} \cos \dfrac{\theta_2}{2} e^{i \phi_1} \\\sin \dfrac{\theta_1}{2} \sin \dfrac{\theta_2}{2} e^{i(\phi_1+\phi_2)}\end{pmatrix} \\&\equiv \begin{pmatrix} p \\ q \\ r \\ s \end{pmatrix}\tag{48}

となります.これが,スピン1が (\theta_1, \phi_1) ,スピン2が (\theta_2, \phi_2) を向いた時の |\alpha_1 \alpha_2 \rangle 状態です.要するに,テンソル積とは演算子1はスピノール1だけに,演算子2はスピノール2だけに作用するということを, 行列(テンソル)で表したものなのです.ちなみにこの式は規格化されています.

他にも, |\alpha_1 \beta_2 \rangle|\beta_1 \alpha_2 \rangle|\beta_1 \beta_2 \rangle の状態があり,それぞれ,次に対応します.

\chi_{1+} \otimes \chi_{2-} &= \begin{pmatrix} -\cos \dfrac{\theta_1}{2} \sin \dfrac{\theta_2}{2} \\\cos \dfrac{\theta_1}{2} \cos \dfrac{\theta_2}{2} e^{i \phi_2} \\-\sin \dfrac{\theta_1}{2} \sin \dfrac{\theta_2}{2} e^{i \phi_1} \\\sin \dfrac{\theta_1}{2} \cos \dfrac{\theta_2}{2} e^{i(\phi_1+\phi_2)}\end{pmatrix} \tag{49} \chi_{1-} \otimes \chi_{2+} &= \begin{pmatrix} -\sin \dfrac{\theta_1}{2} \cos \dfrac{\theta_2}{2} \\-\sin \dfrac{\theta_1}{2} \sin \dfrac{\theta_2}{2} e^{i \phi_2} \\\cos \dfrac{\theta_1}{2} \cos \dfrac{\theta_2}{2} e^{i \phi_1} \\\cos \dfrac{\theta_1}{2} \sin \dfrac{\theta_2}{2} e^{i(\phi_1+\phi_2)}\end{pmatrix} \tag{50} \chi_{1-} \otimes \chi_{2-} &= \begin{pmatrix} \sin \dfrac{\theta_1}{2} \sin \dfrac{\theta_2}{2} \\-\sin \dfrac{\theta_1}{2} \cos \dfrac{\theta_2}{2} e^{i \phi_2} \\-\cos \dfrac{\theta_1}{2} \sin \dfrac{\theta_2}{2} e^{i \phi_1} \\\cos \dfrac{\theta_1}{2} \cos \dfrac{\theta_2}{2} e^{i(\phi_1+\phi_2)}\end{pmatrix} \tag{51}

この4つのスピノール は, \bm{n}_i = \begin{pmatrix} \sin \theta_i \cos \phi_i \\ \sin \theta_i \sin \phi_i \\ \cos \theta_i \end{pmatrix} とした時,

四次の行列の固有値を \lambda とする固有値方程式

\sum_{x}^z \left( \bm{n}_1 \cdot (\bm{S}_1 \otimes 1_2) + \bm{n}_2 \cdot (1_2 \otimes \bm{S}_2 ) \right) \chi_4 = \lambda \chi_4 \tag{52}

の固有状態となるようです.それぞれ,順に \lambda=2,0,0,-2 の固有値を持っています. これは式 (5) の一電子から二電子への拡張です.

ただ,私が具体例を計算したところ,それぞれの状態に直感的に納得できないところがありました. その辺の事情はまた別の機会に書こうと思います.

さあ,いよいよ一般の形の交換相互作用の値 (\bm{S}_1 \cdot \bm{S}_2)_{general} を求めてみましょう! 状態ベクトル(スピノール) \chi_4 は式 (48) の最後の形をしているものとしましょう. それには,次のようにすればいいのだと思います.

(\bm{S}_1 \cdot \bm{S}_2)_{general} &= \chi_4^\dagger \bm{S}_1 \cdot \bm{S}_2 \chi_4 \\&= \dfrac{\hbar^2}{4} \begin{pmatrix} p^\ast & q^\ast & r^\ast & s^\ast \end{pmatrix}\begin{pmatrix}1 & 0 & 0 & 0 \\0 & -1 & 2 & 0 \\0 & 2 & -1 & 0 \\0 & 0 & 0 & 1 \end{pmatrix}\begin{pmatrix} p \\ q \\ r \\ s \end{pmatrix}\tag{53}

三重項のときは, \begin{pmatrix} p \\ q \\ r \\ s \end{pmatrix} = \begin{pmatrix} 1 \\ 0 \\ 0 \\ 0 \end{pmatrix}, \dfrac{1}{\sqrt{2}}\begin{pmatrix} 0 \\ 1 \\ 1 \\ 0 \end{pmatrix}, \begin{pmatrix} 0 \\ 0 \\ 0 \\ 1 \end{pmatrix} であり,一重項のときは, \begin{pmatrix} p \\ q \\ r \\ s \end{pmatrix} = \dfrac{1}{\sqrt{2}}\begin{pmatrix} 0 \\ 1 \\ -1 \\ 0 \end{pmatrix} です.代入して確認すると,確かに (\bm{S}_1 \cdot \bm{S}_2)_{general} = \dfrac{\hbar^2}{4},\dfrac{\hbar^2}{4},\dfrac{\hbar^2}{4},\dfrac{-3 \hbar^2}{4} となります.一重項と三重項の拡張になっていることが分かるかと思います.今日はこの辺で,お疲れ様でした.