ゼロ磁場分裂 (Zero-Field Splitting: ZFS) テンソルの数理と物理的起源：スピン双極子・スピン軌道相互作用から DFT 定式化まで

導入

前回の解説記事『スピン系におけるブロッホ緩和の数理』において、スピン量子数 $S \ge 1$ を持つ電子スピン系（三重項状態や遷移金属錯体、量子ビットとして期待される NV センターなど）では、単なる等方的なスピン $S=1/2$ を前提とした現象論的ブロッホ方程式が破綻することを論じた。その物理的本質は、外部磁場が一切存在しない状態（Zero Field）であっても、異方的な結晶場やスピン相互作用によってスピンエネルギー準位が分裂する ゼロ磁場分裂（Zero-Field Splitting: ZFS） の存在にある。

ZFS を記述する異方的スピンハミルトニアン $\hat{H}_{\text{ZFS}} = \mathbf{S} \cdot \mathbf{D} \cdot \mathbf{S}$ の中心要素である ZFS テンソル $\mathbf{D}$ は、単なる幾何学的な定数ではなく、物質内部の多体電子構造や相対論的効果を直接反映する物理量である。

本稿では、ZFS テンソル $\mathbf{D}$ の現象論的定式化と対角化（主軸系 PAS における軸異方性パラメータ $D$ と斜交異方性パラメータ $E$ の定義）から出発し、その 2大微視的起源である スピン-スピン双極子相互作用（SS 寄与） と スピン-軌道相互作用の 2次摂動（SOC 寄与） を一切の数式省略を行わずに導出する。さらに、現代の第一原理電子状態計算において、これらの物理的寄与が 密度汎関数理論（DFT） の枠組み（2粒子スピン密度表現および Pederson-Khanna 法 / CP-DFT 摂動論）でどのように定式化されるかまでを完全解読する。

1. 現象論的定式化と ZFS テンソルの主軸系 (PAS) 対角化

1.1 現象論的スピンハミルトニアンの導入

外部静磁場が存在しない条件（Zero Field, $\mathbf{B} = \mathbf{0}$）において、スピン量子数 $S \ge 1$ を有する電子スピン系のエネルギー準位分裂を記述する有効スピンハミルトニアン（effective spin Hamiltonian）は、スピン演算子成分 $\hat{S}_x, \hat{S}_y, \hat{S}_z$ の 2次形式（双線形形式）として次のように現象論的に定義される。

\[\hat{H}_{\text{ZFS}} = \mathbf{S} \cdot \mathbf{D} \cdot \mathbf{S} = \sum_{i \in \{x,y,z\}} \sum_{j \in \{x,y,z\}} D_{ij} \hat{S}_i \hat{S}_j\]

ここで、$\mathbf{S} = (\hat{S}_x, \hat{S}_y, \hat{S}_z)^T$ は総スピン演算子ベクトルであり、$\mathbf{D}$ は $3 \times 3$ 実数行列（ZFS テンソル）である。

1.2 テンソルの対称性とトレースレス（無痕跡）条件

(1) 対称性 $D_{ij} = D_{ji}$

スピン演算子成分の量子力学的交換関係 $[\hat{S}i, \hat{S}_j] = i \hbar \sum_k \epsilon{ijk} \hat{S}k$ を考慮すると、テンソル $\mathbf{D}$ の反対称成分 $\frac{1}{2}(D{ij} - D_{ji})$ は、1次のスピン項 $\hat{S}_k$ （スピン-軌道項の 1次効果や実効的磁場効果）へ寄与するに過ぎない。本稿で扱う ZFS はスピン間の 2次異方性相互作用であるため、物理的描像を損なうことなく $\mathbf{D}$ を実対称テンソルとして選ぶことができる：

\[D_{ij} = D_{ji} \quad (i, j \in \{x, y, z\})\]

(2) トレースレス条件 $\text{Tr}(\mathbf{D}) = 0$

任意のスピン量子数 $S$ に対して、全スピン自乗演算子の恒等式 $\hat{S}_x^2 + \hat{S}_y^2 + \hat{S}_z^2 = \mathbf{S}^2 = S(S+1) \hat{I}$ （$\hat{I}$ は単位演算子）が成立する。これにより、テンソル $\mathbf{D}$ の対角成分の平均値（トレースの $1/3$）に相当する等方的成分は、単なる定数エネルギーシフトとしてスピン異方性から分離可能である。

実際に、テンソル $\mathbf{D}$ を等方成分とトレースレスな異方成分へ分解すると：

\[\mathbf{D} = \mathbf{D}' + \frac{1}{3} \text{Tr}(\mathbf{D}) \mathbf{I}_3 \quad \left( \text{Tr}(\mathbf{D}') = 0 \right)\]

この分解をスピンハミルトニアンに代入すると、

\[\hat{H}_{\text{ZFS}} = \mathbf{S} \cdot \mathbf{D}' \cdot \mathbf{S} + \frac{1}{3} \text{Tr}(\mathbf{D}) S(S+1) \hat{I}\]

となる。右辺第 2 項の等方的項 $\frac{1}{3} \text{Tr}(\mathbf{D}) S(S+1) \hat{I}$ はすべての $m_s$ 状態に対して一律のエネルギーシフトを与えるため、これを参照基準エネルギーへ吸収させることで、物理的に意味を持つ異方性テンソル $\mathbf{D}’$ は常にトレースレス（$\text{Tr}(\mathbf{D}’) = 0$）であると定義できる。以下では記号の簡略化のため、標準化されたトレースレス対称テンソルを改めて $\mathbf{D}$ と表記する。

\[\text{Tr}(\mathbf{D}) = D_{xx} + D_{yy} + D_{zz} = 0\]

1.3 主軸系（Principal Axis System: PAS）への対角化

実対称行列 $\mathbf{D}$ は、適切な 3次元回転行列 $\mathbf{R} \in SO(3)$ による直交変換 $\mathbf{D}_{\text{PAS}} = \mathbf{R} \mathbf{D} \mathbf{R}^T$ によって必ず対角化される。この固有座標系を 主軸系 (Principal Axis System: PAS) と呼ぶ。

PAS 座標系 $(x, y, z)$ における ZFS テンソルの主軸成分を $D_{xx}, D_{yy}, D_{zz}$ と置くと、スピンハミルトニアンは交差項 $\hat{S}_i \hat{S}_j \, (i \ne j)$ が消去され、以下のように対角化される：

\[\hat{H}_{\text{ZFS}} = D_{xx} \hat{S}_x^2 + D_{yy} \hat{S}_y^2 + D_{zz} \hat{S}_z^2\]

トレースレス条件より、主軸成分の間には以下の拘束条件が課される：

\[D_{xx} + D_{yy} + D_{zz} = 0\]

1.4 軸異方性パラメータ $D$ と斜交異方性パラメータ $E$ の定義

3つの主軸成分 $D_{xx}, D_{yy}, D_{zz}$ はトレースゼロの条件 $D_{xx} + D_{yy} + D_{zz} = 0$ により、独立な変数は 2つのみである。分光実験および物性解析において直感的な解釈を可能にするため、独立パラメータとして 軸方向 ZFS パラメータ (Axial ZFS Parameter) $D$ と 斜交 ZFS パラメータ (Rhombic ZFS Parameter) $E$ が次のように定義される。

\[D = D_{zz} - \frac{1}{2}(D_{xx} + D_{yy})\] \[E = \frac{1}{2}(D_{xx} - D_{yy})\]

ここで、トレースレス条件 $D_{xx} + D_{yy} = -D_{zz}$ を $D$ の定義式へ代入すると、

\[D = D_{zz} - \frac{1}{2}(-D_{zz}) = \frac{3}{2} D_{zz} \quad \implies \quad D_{zz} = \frac{2}{3} D\]

が得られる。また、$D_{xx}$ および $D_{yy}$ については：

\[D_{xx} + D_{yy} = -D_{zz} = -\frac{2}{3} D\] \[D_{xx} - D_{yy} = 2E\]

これらを和と差によって解くことで、主軸成分を $D$ と $E$ のみで表す表現が得られる：

\[D_{xx} = -\frac{1}{3} D + E\] \[D_{yy} = -\frac{1}{3} D - E\] \[D_{zz} = \frac{2}{3} D\]

1.5 標準形 ZFS スピンハミルトニアンの導出

得られた主軸成分 $D_{xx}, D_{yy}, D_{zz}$ を対角化後のハミルトニアン $\hat{H}{\text{ZFS}} = D{xx} \hat{S}x^2 + D{yy} \hat{S}y^2 + D{zz} \hat{S}_z^2$ へ代入して整理する。

\[\begin{aligned} \hat{H}_{\text{ZFS}} &= \left( -\frac{1}{3} D + E \right) \hat{S}_x^2 + \left( -\frac{1}{3} D - E \right) \hat{S}_y^2 + \frac{2}{3} D \hat{S}_z^2 \\ &= -\frac{1}{3} D \left( \hat{S}_x^2 + \hat{S}_y^2 \right) + E \left( \hat{S}_x^2 - \hat{S}_y^2 \right) + \frac{2}{3} D \hat{S}_z^2 \end{aligned}\]

全スピン自乗演算子の関係式 $\hat{S}_x^2 + \hat{S}_y^2 = \mathbf{S}^2 - \hat{S}_z^2 = S(S+1) \hat{I} - \hat{S}_z^2$ を用いて $\hat{S}_x^2 + \hat{S}_y^2$ を消去すると、

\[\begin{aligned} \hat{H}_{\text{ZFS}} &= -\frac{1}{3} D \left( S(S+1) \hat{I} - \hat{S}_z^2 \right) + \frac{2}{3} D \hat{S}_z^2 + E \left( \hat{S}_x^2 - \hat{S}_y^2 \right) \\ &= D \left( \hat{S}_z^2 - \frac{1}{3} S(S+1) \hat{I} \right) + E \left( \hat{S}_x^2 - \hat{S}_y^2 \right) \end{aligned}\]

となり、電子スピン共鳴（EPR）分光や磁気化学において最も広く用いられる 標準形 ZFS スピンハミルトニアン (Standard ZFS Spin Hamiltonian) が厳密に導出される。

PAS 軸選択のコンベンション（正規化条件） 主軸 $x, y, z$ の順序付けには回転の自由度があるが、慣例として主軸成分の絶対値の大きさが $\vert D_{zz} \vert \ge \vert D_{yy} \vert \ge \vert D_{xx} \vert$ となるように $z$ 軸（主軸）を定義する。この規約のもとで、軸パラメータと斜交パラメータの比は以下の範囲に規格化される：

\[0 \le \frac{E}{D} \le \frac{1}{3}\]

1.6 エネルギー準位分裂ダイアグラムと物理的意味 ($S=1, S=3/2$)

ZFS パラメータ $D, E$ が与えるスピン準位の分裂挙動を、$S=1$（三重項）および $S=3/2$（四重項）の具体的な系を例に解析する。

(1) $S=1$ スピン系（三重項状態：NV センター、有機ビラジカル、高スピン Ni(II) 錯体など）

スピン基底として $z$ 軸の固有状態 $\lbrace \vert m_s = 1 \rangle, \vert m_s = 0 \rangle, \vert m_s = -1 \rangle \rbrace$ をとる。

軸対称極限 ($E=0$)

斜交性が存在しない場合 ($E=0$)、スピンハミルトニアンは $\hat{H}_{\text{ZFS}} = D \left( \hat{S}_z^2 - \frac{2}{3} \hat{I} \right)$ となり、固有状態は $m_s$ 基底そのものである：

• $\vert m_s = 0 \rangle$: $E(0) = -\frac{2}{3} D$
• $\vert m_s = \pm 1 \rangle$: $E(\pm 1) = D \left( 1 - \frac{2}{3} \right) = \frac{1}{3} D$ （2重縮退）

無磁場下における準位分裂幅は $\Delta E = E(\pm 1) - E(0) = D$ となる。

斜交性が存在する場合 ($E \ne 0$)

$\hat{S}x^2 - \hat{S}_y^2 = \frac{1}{2}(\hat{S}+^2 + \hat{S}_-^2)$ 項の存在により、$\Delta m_s = \pm 2$ の状態間（すなわち $\vert 1 \rangle$ と $\vert -1 \rangle$）に相互作用が生じ、$\vert m_s = \pm 1 \rangle$ の 2重縮退が解除される。

ハミルトニアンの固有状態と固有エネルギーは以下のようになる：

• $\vert T_z \rangle = \vert m_s = 0 \rangle \quad \implies \quad E_z = -\frac{2}{3} D$
• $\vert T_x \rangle = \frac{1}{\sqrt{2}} \left( \vert 1 \rangle + \vert -1 \rangle \right) \quad \implies \quad E_x = \frac{1}{3} D + E$
• $\vert T_y \rangle = \frac{1}{\sqrt{2} i} \left( \vert 1 \rangle - \vert -1 \rangle \right) \quad \implies \quad E_y = \frac{1}{3} D - E$

$\vert T_x \rangle$ と $\vert T_y \rangle$ の分裂幅は $E_x - E_y = 2E$ であり、斜交パラメータ $E$ が面内（$xy$ 平面）の異方性分裂を直接反映していることがわかる。

[S = 1 スピン準位分裂ダイアグラム]

(A) 軸対称 (E = 0, D > 0)          (B) 斜交性あり (E > 0, D > 0)

  ------- |ms = ±1>  (E = D/3)        ------- |Tx>  (E = D/3 + E)
                                        | 2E
                                      ------- |Ty>  (E = D/3 - E)
     | D
     |
  ------- |ms = 0>   (E = -2D/3)      ------- |Tz>  (E = -2D/3)

(2) $S=3/2$ スピン系（四重項状態：高スピン Co(II), Cr(III), Fe(III) 低スピン等）

スピン基底として $\lbrace \vert \pm 3/2 \rangle, \vert \pm 1/2 \rangle \rbrace$ をとる。

軸対称極限 ($E=0$)

ハミルトニアン $\hat{H}_{\text{ZFS}} = D \left( \hat{S}_z^2 - \frac{5}{4} \hat{I} \right)$ のもとでの固有値は：

• $\vert m_s = \pm 1/2 \rangle$: $E(\pm 1/2) = D \left( \frac{1}{4} - \frac{5}{4} \right) = -D$ （クラマース 2重縮退）
• $\vert m_s = \pm 3/2 \rangle$: $E(\pm 3/2) = D \left( \frac{9}{4} - \frac{5}{4} \right) = +D$ （クラマース 2重縮退）

無磁場分裂幅は $\Delta E = E(\pm 3/2) - E(\pm 1/2) = 2D$ となる。クラマースの定理（Kramers’ Theorem）により、奇数電子系（半整数スピン系）では外部磁場がない限り、各準位は最低でも 2重縮退（Kramers doublet）を保つ。

斜交性が存在する場合 ($E \ne 0$)

$E \ne 0$ の場合であっても Kramers 2重縮退自体は解けないが、$\vert \pm 3/2 \rangle$ と $\vert \mp 1/2 \rangle$ が混合し、実効的な $g$ テンソル異方性や二準位間のエネルギー差が変化する。

[S = 3/2 スピン準位分裂ダイアグラム]

(A) D > 0 (Easy-plane)                (B) D < 0 (Easy-axis)

  ======= |ms = ±3/2> (E = +D)          ======= |ms = ±1/2> (E = +|D|)
     |                                     |
     | 2D                                  | 2|D|
     |                                     |
  ======= |ms = ±1/2> (E = -D)          ======= |ms = ±3/2> (E = -|D|)

(3) 容易軸（Easy-axis）異方性 vs 容易面（Easy-plane）異方性の物理的意味

軸パラメータ $D$ の符号は、スピンが結晶構造のどの方向に配向しやすいかという磁気異方性の巨視的性質を決定する。

• 容易軸（Easy-axis）異方性 ($D < 0$) $D < 0$ のとき、$m_s$ の絶対値が大きい状態（$S=1$ では $m_s = \pm 1$、$S=3/2$ では $m_s = \pm 3/2$）が基底状態となる。スピンは $z$ 軸（容易軸）方向を向きやすくなり、$m_s = +S$ と $m_s = -S$ の間にエネルギー障壁（$S=1$ で $U = \vert D \vert$）が形成される。この反転障壁が存在することにより、単分子磁石 (Single-Molecule Magnets: SMMs) や高保持力磁性体では、高温での熱的スピン反転が抑制され、低温での磁気ヒステリシスや量子スピントンネリングが観測される。
• 容易面（Easy-plane）異方性 ($D > 0$) $D > 0$ のとき、$m_s = 0$ や $m_s = \pm 1/2$ が基底状態となり、スピンは $z$ 軸に垂直な $xy$ 平面内を向きやすくなる（容易面）。この場合、$z$ 軸方向への巨視的磁化の固定は起こりにくくなる。

2. スピン-スピン双極子相互作用 (SS 寄与) の微視的導出

前節では、ZFS テンソル $\mathbf{D}$ の現象論的定義と、主軸系（PAS）における軸異方性パラメータ $D$ および斜交異方性パラメータ $E$ による準位分裂の記述を論じた。本節からは、ZFS テンソルの 1大微視的起源である スピン-スピン双極子相互作用（Spin-Spin Dipolar Interaction: SS 寄与） の量子力学的微視的導出を行う。

2.1 2電子スピン磁気双極子-双極子相互作用ハミルトニアン

真空中に置かれた 2つの電子 $i, j$ のスピン磁気双極子モーメント演算子は、電子スピン演算子 $\mathbf{s}_i, \mathbf{s}_j$ を用いて次のように表される：

\[\boldsymbol{\mu}_i = -g_e \mu_B \mathbf{s}_i, \quad \boldsymbol{\mu}_j = -g_e \mu_B \mathbf{s}_j\]

ここで、$g_e \approx 2.002319$ は自由電子の $g$ 因子、$\mu_B = \frac{e\hbar}{2m_e}$ はボーア磁子である。古典電磁気学において、距離 $\mathbf{r}_{ij} = \mathbf{r}_i - \mathbf{r}_j$ 離れた 2つの点磁気双極子 $\boldsymbol{\mu}_i, \boldsymbol{\mu}_j$ の相互作用エネルギーは次式で与えられる：

\[U_{ij} = \frac{\mu_0}{4\pi} \left[ \frac{\boldsymbol{\mu}_i \cdot \boldsymbol{\mu}_j}{r_{ij}^3} - \frac{3(\boldsymbol{\mu}_i \cdot \mathbf{r}_{ij})(\boldsymbol{\mu}_j \cdot \mathbf{r}_{ij})}{r_{ij}^5} \right]\]

ただし、$\mu_0$ は真空の透磁率、$r_{ij} = \vert \mathbf{r}_{ij} \vert = \vert \mathbf{r}_i - \mathbf{r}_j \vert$ は電子間距離である。

多電子系（全電子数 $N$）において、すべての電子対についての和を取り、二重カウントを防ぐ係数 $1/2$ を乗じることで、スピン-スピン双極子相互作用ハミルトニアン $\hat{H}_{\text{SS}}$ が得られる：

\[\hat{H}_{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi} \frac{1}{2} \sum_{i \neq j}^N \left[ \frac{\mathbf{s}_i \cdot \mathbf{s}_j}{r_{ij}^3} - 3 \frac{(\mathbf{s}_i \cdot \mathbf{r}_{ij})(\mathbf{s}_j \cdot \mathbf{r}_{ij})}{r_{ij}^5} \right]\]

単位系についての補足 原子単位系（Hartree atomic units: $\hbar = m_e = e = 4\pi\epsilon_0 = 1$）では、真空の透磁率は $\mu_0 = 4\pi \alpha^2$ （$\alpha \approx 1/137.036$ は微細構造定数）となり、$\mu_B = 1/2$ である。したがって、前出の係数は以下のようにシンプルに書くことができる：

\[\frac{\mu_0 g_e^2 \mu_B^2}{4\pi} = \frac{\alpha^2 g_e^2}{4} \approx \alpha^2\]

本稿では物理的意味（物理定数依存性）を明確に保つため、SI 単位系での表記 $\frac{\mu_0 g_e^2 \mu_B^2}{4\pi}$ を基本とする。

2.2 テンソル表現とデカルト成分展開

ハミルトニアン $\hat{H}_{\text{SS}}$ のデカルト成分（$a, b \in {x, y, z}$）を抽出するため、スピン演算子および相対位置ベクトルの積を代数的に展開する。

相対位置ベクトルを $\mathbf{r}{ij} = (x{ij}, y_{ij}, z_{ij})^T = (r_{ij, x}, r_{ij, y}, r_{ij, z})^T$ と表記すると、分子の各項はデカルト成分を用いて次のように書き下せる：

\[\mathbf{s}_i \cdot \mathbf{s}_j = \sum_{a \in \{x,y,z\}} s_{i, a} s_{j, a}\] \[(\mathbf{s}_i \cdot \mathbf{r}_{ij})(\mathbf{s}_j \cdot \mathbf{r}_{ij}) = \left( \sum_{a \in \{x,y,z\}} s_{i, a} a_{ij} \right) \left( \sum_{b \in \{x,y,z\}} s_{j, b} b_{ij} \right) = \sum_{a, b \in \{x,y,z\}} s_{i, a} s_{j, b} a_{ij} b_{ij}\]

ここで、クロネッカーのデルタ $\delta_{ab}$ を用いて $\mathbf{s}i \cdot \mathbf{s}_j = \sum{a,b} s_{i, a} s_{j, b} \delta_{ab}$ と統一的に表すと、$\hat{H}{\text{SS}}$ はスピン演算子成分 $s{i, a} s_{j, b}$ と空間双極子テンソル演算子 $\hat{h}{ab}(\mathbf{r}{ij})$ の積の形に再構成される：

\[\hat{H}_{\text{SS}} = \frac{1}{2} \sum_{i \neq j}^N \sum_{a \in \{x,y,z\}} \sum_{b \in \{x,y,z\}} s_{i, a} \, \hat{h}_{ab}(\mathbf{r}_{ij}) \, s_{j, b}\]

ここで、空間双極子テンソル演算子 $\hat{h}{ab}(\mathbf{r}{ij})$ は次のように定義される：

\[\hat{h}_{ab}(\mathbf{r}_{ij}) = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi} \frac{r_{ij}^2 \delta_{ab} - 3 a_{ij} b_{ij}}{r_{ij}^5}\]

この 2階空間テンソル $\hat{h}{ab}(\mathbf{r}{ij})$ は、定義より直ちに以下の重要性質を満たす：

1. 対称性: $\hat{h}{ab}(\mathbf{r}{ij}) = \hat{h}{ba}(\mathbf{r}{ij})$

2. トレースレス性（無痕跡性）:

   \[\sum_{a \in \{x,y,z\}} \hat{h}_{aa}(\mathbf{r}_{ij}) = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi r_{ij}^5} \left[ r_{ij}^2 \sum_{a} \delta_{aa} - 3 \sum_{a} a_{ij}^2 \right] = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi r_{ij}^5} \left[ 3 r_{ij}^2 - 3 (x_{ij}^2 + y_{ij}^2 + z_{ij}^2) \right] = 0\]

したがって、空間的双極子演算子それ自体が既に完全なトレースレス対称テンソル構造を有している。

2.3 全スピン演算子 $\mathbf{S}$ への射影（Wigner-Eckart の定理）

微視的ハミルトニアン $\hat{H}{\text{SS}}$ は個々の電子スピン演算子 $s{i, a}, s_{j, b}$ に依存している。これを巨視的な有効スピンハミルトニアン $\hat{H}_{\text{ZFS}}^{\text{SS}} = \mathbf{S} \cdot \mathbf{D}^{\text{SS}} \cdot \mathbf{S}$ （$\mathbf{S} = \sum_i \mathbf{s}_i$ は全スピン演算子）の形式へ縮約するため、Wigner-Eckart の定理（およびスピン射影定理 Projection Theorem） を適用する。

全スピン量子数が $S$（ただし $S \ge 1$）である状態空間 $\mathcal{H}S = \text{span} { \vert S, M_S \rangle }{M_S = -S}^S$ 内において、個別のスピン演算子の積の対称テンソル結合 $\sum_{i \neq j} s_{i, a} s_{j, b}$ は、全スピン演算子の 2次形式に射影される。

Wigner-Eckart の定理によれば、階数 2 の既約テンソル演算子 $\hat{T}_m^{(2)}$ の多重度 $2S+1$ 内での行列要素は、全スピン演算子から作られる階数 2 のテンソル演算子の行列要素と比例関係にある。デカルト表現において、スピン空間での縮約式は以下の通りとなる：

\[P_S \left( \sum_{i \neq j}^N s_{i, a} s_{j, b} \right) P_S = C_S \left[ \frac{1}{2} (\hat{S}_a \hat{S}_b + \hat{S}_b \hat{S}_a) - \frac{1}{3} S(S+1) \delta_{ab} \hat{I} \right]\]

ここで $P_S$ は全スピン $S$ の部分空間への射影演算子であり、$C_S$ はスピン量子数 $S$ に依存する射影定数（スピン Wigner-Eckart 因子）である。

射影定数 $C_S$ の厳密決定

定数 $C_S$ の値を決定するため、具体的に $a=z, b=z$ 成分の対角期待値を最高スピン状態 $\vert S, M_S = S \rangle$ で計算する。

まず、全スピン演算子側の 2次形式の期待値は：

\[\left\langle S, S \left\vert \hat{S}_z^2 - \frac{1}{3} S(S+1) \hat{I} \right\vert S, S \right\rangle = S^2 - \frac{1}{3} S(S+1) = \frac{S(2S-1)}{3}\]

一方、全スピン自乗演算子 $\mathbf{S}^2 = (\sum_i \mathbf{s}i)^2 = \sum_i \mathbf{s}_i^2 + \sum{i \neq j} \mathbf{s}_i \cdot \mathbf{s}_j$ の関係を用いると、個別スピン相互作用のテンソル平均との比較から、射影定数は以下のように求まる：

\[C_S = \frac{1}{S(2S-1)}\]

スピン条件 $S \ge 1$ の物理的含意 射影因子 $C_S = \frac{1}{S(2S-1)}$ の分母に注目すると、$S = 1/2$ のとき分母が $1/2 \times (2(1/2) - 1) = 0$ となり発散する。物理的には、$S=1/2$ （単一の非対電子または全スピン半磁性体）では 2階スピンテンソル演算子が存在し得ず（$S=1/2$ の Pauli 行列の積は常に 1次のスピン行列または単位行列に還元されるため）、双極子相互作用による ZFS は厳密にゼロとなることを意味する。ZFS が有限の値を持つためには $S \ge 1$（対のない電子が 2つ以上存在すること）が必須不可欠である。

2.4 ZFS テンソル $D_{ab}^{\text{SS}}$ の空間積分表現の導出

Wigner-Eckart の射影関係式をハミルトニアン $\hat{H}_{\text{SS}}$ に代入し、空間部分の多電子波動関数 $\Psi(\mathbf{r}_1, \mathbf{r}_2, \dots, \mathbf{r}_N)$ に関する積分を実行する。

有効スピンハミルトニアンとの対応関係：

\[\hat{H}_{\text{eff}}^{\text{SS}} = \sum_{a, b \in \{x,y,z\}} D_{ab}^{\text{SS}} \hat{S}_a \hat{S}_b = \mathbf{S} \cdot \mathbf{D}^{\text{SS}} \cdot \mathbf{S}\]

（ただし $\mathbf{D}^{\text{SS}}$ はトレースレスであるため $\sum_a D_{aa}^{\text{SS}} = 0$ により $\sum_{a,b} D_{ab}^{\text{SS}} \frac{1}{2}(\hat{S}_a \hat{S}_b + \hat{S}_b \hat{S}_a) = \mathbf{S} \cdot \mathbf{D}^{\text{SS}} \cdot \mathbf{S}$ と等価になる）。

空間双極子テンソル演算子の期待値を計算することで、ZFS テンソル成分 $D_{ab}^{\text{SS}}$ の微視的算式が得られる：

\[D_{ab}^{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi} \frac{1}{2} \frac{1}{S(2S-1)} \left\langle \Psi \left\vert \sum_{i \neq j}^N \frac{r_{ij}^2 \delta_{ab} - 3 a_{ij} b_{ij}}{r_{ij}^5} \right\vert \Psi \right\rangle\]

電子対の和を $\sum_{i \neq j} = 2 \sum_{i < j}$ と書き改めると、前の係数 $1/2$ と相殺し、以下の最終的な 2電子空間積分表現 が導出される：

\[D_{ab}^{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi} \frac{1}{S(2S-1)} \left\langle \Psi \left\vert \sum_{i < j}^N \frac{r_{ij}^2 \delta_{ab} - 3 (r_{ij})_a (r_{ij})_b}{r_{ij}^5} \right\vert \Psi \right\rangle\]

ここで、$(r_{ij})a$ は 2電子間相対ベクトル $\mathbf{r}{ij} = \mathbf{r}_i - \mathbf{r}_j$ の第 $a$ デカルト成分（$a \in {x, y, z}$）を表す。

各成分を極めて明示的に書き下すと以下のようになる：

\[D_{xx}^{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi S(2S-1)} \left\langle \Psi \left\vert \sum_{i < j} \frac{r_{ij}^2 - 3 x_{ij}^2}{r_{ij}^5} \right\vert \Psi \right\rangle = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi S(2S-1)} \left\langle \Psi \left\vert \sum_{i < j} \frac{y_{ij}^2 + z_{ij}^2 - 2x_{ij}^2}{r_{ij}^5} \right\vert \Psi \right\rangle\] \[D_{yy}^{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi S(2S-1)} \left\langle \Psi \left\vert \sum_{i < j} \frac{r_{ij}^2 - 3 y_{ij}^2}{r_{ij}^5} \right\vert \Psi \right\rangle = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi S(2S-1)} \left\langle \Psi \left\vert \sum_{i < j} \frac{x_{ij}^2 + z_{ij}^2 - 2y_{ij}^2}{r_{ij}^5} \right\vert \Psi \right\rangle\] \[D_{zz}^{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi S(2S-1)} \left\langle \Psi \left\vert \sum_{i < j} \frac{r_{ij}^2 - 3 z_{ij}^2}{r_{ij}^5} \right\vert \Psi \right\rangle = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi S(2S-1)} \left\langle \Psi \left\vert \sum_{i < j} \frac{x_{ij}^2 + y_{ij}^2 - 2z_{ij}^2}{r_{ij}^5} \right\vert \Psi \right\rangle\] \[D_{ab}^{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{4\pi S(2S-1)} \left\langle \Psi \left\vert \sum_{i < j} \frac{-3 a_{ij} b_{ij}}{r_{ij}^5} \right\vert \Psi \right\rangle \quad (a \neq b)\]

2.5 2粒子スピン密度行列（対密度）による表現

波動関数 $\Psi$ 全体を用いた $3N$ 次元積分の記述から、第一原理電子状態計算（DFT や ab initio 量子化学）で用いられる実空間 6次元積分へ書き換えるため、2粒子スピン密度の差（Spin-Unpaired Pair Density） $\rho_2^{\text{spin}}(\mathbf{r}_1, \mathbf{r}_2)$ を導入する。

空間位置 $\mathbf{r}_1, \mathbf{r}_2$ に 2つの電子が存在する確率振幅のうち、平行スピン対（$\alpha\alpha, \beta\beta$）と反平行スピン対（$\alpha\beta, \beta\alpha$）の分布密度の差を $\rho_2^{\text{spin}}(\mathbf{r}_1, \mathbf{r}_2)$ と定義すると、テンソル成分は連続空間の重積分として定式化される：

\[D_{ab}^{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{8\pi S(2S-1)} \iint \rho_2^{\text{spin}}(\mathbf{r}_1, \mathbf{r}_2) \frac{\vert \mathbf{r}_1 - \mathbf{r}_2 \vert^2 \delta_{ab} - 3 (\mathbf{r}_1 - \mathbf{r}_2)_a (\mathbf{r}_1 - \mathbf{r}_2)_b}{\vert \mathbf{r}_1 - \mathbf{r}_2 \vert^5} \, d\mathbf{r}_1 d\mathbf{r}_2\]

（注：積分変数 $\mathbf{r}1, \mathbf{r}_2$ の独立積分にする際、$\sum{i < j}$ から全空間積分への変換で因子 $1/2$ が生じる）。

2.6 幾何学的異方性と距離依存性の物理的考察

導出された $D_{ab}^{\text{SS}}$ の積分表現式は、スピン-スピン双極子相互作用の物理的特性に関して以下の重要な知見を与える：

1. 距離依存性 ($1/r_{12}^3$) 積分核は $\frac{r_{12}^2 - 3z_{12}^2}{r_{12}^5} \sim \frac{1}{r_{12}^3}$ の距離依存性を持つ。非対電子同士の幾何学的距離 $r_{12}$ が離れるにつれて、SS 寄与は $1/r_{12}^3$ で急速に減衰する。したがって、有機ビラジカルのように 2つのスピン中心が離れている系では ZFS の SS 寄与は非常に小さくなる（数 MHz 〜 数十 MHz 程度）。一方、同一原子の $d$ 軌道や $p$ 軌道内に複数の非対電子が詰まっている遷移金属錯体や NV センターでは、SS 寄与が極めて大きくなり得る（GHz オーダー）。
2. スピン異方性と電子雲の分布形状
   • 電子雲が $z$ 軸方向に伸長した軸長形状（Prolate distribution）を持つ場合、$z_{12}^2 > x_{12}^2 + y_{12}^2$ となるため、$D_{zz}^{\text{SS}} < 0$ （すなわち $D^{\text{SS}} = \frac{3}{2} D_{zz}^{\text{SS}} < 0$：容易軸異方性）となる。
   • 逆に電子雲が $xy$ 平面内に広がった扁平形状（Oblate distribution）を持つ場合、$x_{12}^2 + y_{12}^2 > z_{12}^2$ となるため、$D_{zz}^{\text{SS}} > 0$ （$D^{\text{SS}} > 0$：容易面異方性）となる。

このように、SS 寄与の ZFS テンソル $\mathbf{D}^{\text{SS}}$ は、物質中のスピン密度の空間的異方性を幾何学的に忠実にエンコードした物理量である。

3. スピン-軌道相互作用 (SOC 寄与) と 2次摂動論による微視的導出

前節で導出したスピン-スピン双極子相互作用（SS 寄与）に加え、ZFS テンソル $\mathbf{D}$ のもう一つの重大な微視的起源が スピン-軌道相互作用（Spin-Orbit Coupling: SOC 寄与） である。特に遷移金属錯体や重原子を含む量子スピン系においては、相対論的効果である SOC が SS 寄与を遥かに凌駕し、ZFS の主要な決定要因となる。

本節では、1次 SOC ハミルトニアンの定式化から出発し、結晶場における軌道角運動量の消失（クエンチング）の厳密な証明、1次効果の消失、そして 2次 Rayleigh-Schrödinger 摂動論を用いた ZFS テンソル $D_{ab}^{\text{SOC}}$ の微視的表現式の導出を完全追体験する。

3.1 1次スピン-軌道相互作用ハミルトニアン

多電子系における 1電子スピン-軌道相互作用（Breit-Pauli 1中心有効近似ハミルトニアン）は、各電子 $i$ の軌道角運動量演算子 $\mathbf{l}_i$ とスピン演算子 $\mathbf{s}_i$ の内積を用いて次のように定義される：

\[\hat{H}_{\text{SOC}} = \sum_{i=1}^N \xi(r_i) \mathbf{l}_i \cdot \mathbf{s}_i = \sum_{i=1}^N \sum_{a \in \{x,y,z\}} \xi(r_i) \, l_{i, a} \, s_{i, a}\]

ここで、$\mathbf{l}_i = -i \hbar (\mathbf{r}_i \times \nabla_i)$ は第 $i$ 電子の軌道角運動量演算子ベクトルであり、$\xi(r_i)$ は動径方向の有効中心電場ポテンシャル $V(r_i)$ から定まるスピン-軌道結合関数である：

\[\xi(r_i) = \frac{\hbar^2}{2 m_e^2 c^2} \frac{1}{r_i} \frac{d V(r_i)}{d r_i}\]

原子核電荷 $Z$ を有する球対称ポテンシャル $V(r_i) \propto -Z e^2 / r_i$ を仮定すると、$\xi(r) \propto Z / r^3$ となる。価電子の平均半径 $\langle 1/r^3 \rangle \propto Z^3$ のスケーリングを考慮すると、スピン-軌道相互作用の強さは原子番号の 4 乗（$Z^4$）に概ね比例して急激に増大する（重原子効果 Heavy-Atom Effect）。これが、軽原子中心の有機ビラジカルでは SS 寄与が主導的であるのに対し、3d・4f 遷移金属錯体や重元素系では SOC 寄与が支配的となる理由である。

3.2 軌道角運動量の消失 (Orbital Quenching) と 1次 SOC 寄与の消失

(1) 結晶場・分子場における軌道角運動量の消失証明

自由原子においては回転対称性により軌道角運動量 $\mathbf{L} = \sum_i \mathbf{l}_i$ が保存され、基底状態は非ゼロの $\mathbf{L}$ を持ち得る。しかし、分子や結晶の中に置かれた電子系では、周囲の配位子や結晶場ポテンシャルによって球対称性が破壊される。

定常状態のハミルトニアン $\hat{H}_0$ が時間反転対称性を満たす場合（外部磁場が存在しない場合）、空間部分の波動関数 $\psi_0(\mathbf{r})$ は常に 実数関数 として選ぶことができる：

\[\psi_0^*(\mathbf{r}) = \psi_0(\mathbf{r})\]

ここで、軌道角運動量演算子 $\hat{\mathbf{l}} = -i \hbar (\mathbf{r} \times \nabla)$ の複素共役をとると、虚数単位 $i$ の存在により以下の反対称性が成り立つ：

\[\hat{\mathbf{l}}^* = -\hat{\mathbf{l}}\]

非縮退基底状態 $\vert 0 \rangle$ （空間波動関数 $\psi_0(\mathbf{r})$）における軌道角運動量成分 $\hat{l}_a \, (a \in {x, y, z})$ の期待値を計算する：

\[\langle 0 \vert \hat{l}_a \vert 0 \rangle = \int \psi_0^*(\mathbf{r}) \hat{l}_a \psi_0(\mathbf{r}) \, d\mathbf{r} = \int \psi_0(\mathbf{r}) \hat{l}_a \psi_0(\mathbf{r}) \, d\mathbf{r}\]

この期待値の複素共役をとると：

\[\langle 0 \vert \hat{l}_a \vert 0 \rangle^* = \int \psi_0(\mathbf{r}) \hat{l}_a^* \psi_0(\mathbf{r}) \, d\mathbf{r} = \int \psi_0(\mathbf{r}) (-\hat{l}_a) \psi_0(\mathbf{r}) \, d\mathbf{r} = -\langle 0 \vert \hat{l}_a \vert 0 \rangle\]

一方、軌道角運動量演算子 $\hat{l}_a$ はエルミート演算子であるため、その自己共役性により期待値は実数でなければならない：

\[\langle 0 \vert \hat{l}_a \vert 0 \rangle^* = \langle 0 \vert \hat{l}_a \vert 0 \rangle\]

したがって、

\[\langle 0 \vert \hat{l}_a \vert 0 \rangle = -\langle 0 \vert \hat{l}_a \vert 0 \rangle \quad \implies \quad \langle 0 \vert \hat{l}_a \vert 0 \rangle = 0 \quad (a \in \{x, y, z\})\]

すなわち、多電子系の全軌道角運動量演算子 $\hat{\mathbf{L}} = \sum_i \mathbf{l}_i$ についても同様に：

\[\langle 0 \vert \hat{\mathbf{L}} \vert 0 \rangle = \mathbf{0}\]

が厳密に成立する。これを 軌道角運動量のクエンチング（Orbital Angular Momentum Quenching: 軌道凍結） と呼ぶ。

(2) 1次 SOC 期待値の消失

スピン-軌道相互作用ハミルトニアン $\hat{H}{\text{SOC}}$ の基底状態 $\vert 0 \rangle$ における 1次摂動エネルギー $E{\text{SOC}}^{(1)}$ を評価する。基底状態の空間部分とスピン部分の直積状態において、対角行列要素をとると：

\[E_{\text{SOC}}^{(1)} = \langle 0 \vert \hat{H}_{\text{SOC}} \vert 0 \rangle = \sum_{i=1}^N \sum_{a \in \{x,y,z\}} \langle 0 \vert \xi(r_i) l_{i, a} \vert 0 \rangle \langle S, M_S \vert s_{i, a} \vert S, M_S \rangle\]

ここで、動径関数 $\xi(r_i)$ は完全対称なスカラー関数であり、$\xi(r_i) l_{i, a}$ の空間期待値は前述の軌道クエンチングと全く同様の理由（実数空間波動関数と虚数微分演算子）により厳密にゼロとなる：

\[\langle 0 \vert \xi(r_i) l_{i, a} \vert 0 \rangle = 0\]

したがって、

\[E_{\text{SOC}}^{(1)} = \langle 0 \vert \hat{H}_{\text{SOC}} \vert 0 \rangle = 0\]

が得られる。この結果は物理的に極めて重要である：非縮退な電子基底状態を持つ系においては、スピン-軌道相互作用による 1次のエネルギーシフトおよび ZFS 分裂は厳密にゼロとなる。したがって、SOC に起因する ZFS 効果は、励起電子状態の足し込みを伴う 2次摂動理論によって初めて現れる。

3.3 2次 Rayleigh-Schrödinger 摂動論の展開

1次効果が消失するため、SOC ハミルトニアン $\hat{H}{\text{SOC}}$ を摂動項として扱い、無摂動基底状態 $\vert 0 \rangle$ （エネルギー $E_0$）に対する 2次 Rayleigh-Schrödinger 摂動エネルギー $E{\text{SOC}}^{(2)}$ を定式化する。

すべての励起状態 $\vert n \rangle$ （エネルギー $E_n$、$n \neq 0$）を中間状態として含む 2次摂動表式は以下の通りである：

\[E_{\text{SOC}}^{(2)} = -\sum_{n \neq 0} \frac{\langle 0 \vert \hat{H}_{\text{SOC}} \vert n \rangle \langle n \vert \hat{H}_{\text{SOC}} \vert 0 \rangle}{E_n - E_0}\]

ここで、和 $\sum_{n \neq 0}$ は電子励起空間全体（同種多重項励起および異種多重項励起を含むすべての多電子励起状態）に渡る完全系の和を表す。また、無摂動ハミルトニアンの基底状態が非縮退であるため分母のエネルギー差 $E_n - E_0 > 0$ は常に正の値をとる。

3.4 有効スピンハミルトニアンへのマッピングと ZFS テンソル $D_{ab}^{\text{SOC}}$ の明示的表現の導出

2次摂動エネルギー $E_{\text{SOC}}^{(2)}$ を有効スピン空間上の 2次スピンハミルトニアン $\hat{H}{\text{eff}}^{\text{SOC}} = \sum{a,b \in {x,y,z}} D_{ab}^{\text{SOC}} \hat{S}a \hat{S}_b$ に等価マッピング（縮約）することで、ZFS テンソル成分 $D{ab}^{\text{SOC}}$ の微視的算式を厳密に導出する。

SOC ハミルトニアンをスピン演算子成分 $s_{k, a}$ と軌道演算子成分 $\xi_k l_{k, a}$ に分離して多電子状態間の行列要素を書き下す。

第 $k$ 電子のスピン-軌道演算子和を $\hat{\mathcal{L}}a = \sum{k=1}^N \xi(r_k) l_{k, a}$ と簡略表記すると、SOC ハミルトニアンは：

\[\hat{H}_{\text{SOC}} = \sum_{a \in \{x,y,z\}} \sum_{k=1}^N \xi(r_k) l_{k, a} s_{k, a} = \sum_{a \in \{x,y,z\}} \hat{\mathcal{L}}_a \hat{S}_a\]

のように、実効的な軌道-スピンのデカルト積として表現できる（ここで全スピン空間への Wigner-Eckart 射影を適用する）。

これを 2次摂動エネルギーの分子に代入すると：

\[\langle 0 \vert \hat{H}_{\text{SOC}} \vert n \rangle \langle n \vert \hat{H}_{\text{SOC}} \vert 0 \rangle = \sum_{a \in \{x,y,z\}} \sum_{b \in \{x,y,z\}} \langle 0 \vert \hat{\mathcal{L}}_a \vert n \rangle \langle n \vert \hat{\mathcal{L}}_b \vert 0 \rangle \hat{S}_a \hat{S}_b\]

が得られる。この表現を摂動式へ代入して整理すると：

\[E_{\text{SOC}}^{(2)} = -\sum_{a, b \in \{x,y,z\}} \left[ \sum_{n \neq 0} \frac{\langle 0 \vert \sum_{k=1}^N \xi(r_k) l_{k, a} \vert n \rangle \langle n \vert \sum_{m=1}^N \xi(r_m) l_{m, b} \vert 0 \rangle}{E_n - E_0} \right] \hat{S}_a \hat{S}_b\]

となる。

この 2次スピン項の係数を現象論的スピンハミルトニアン $\hat{H}{\text{eff}}^{\text{SOC}} = \sum{a,b} D_{ab}^{\text{SOC}} \hat{S}_a \hat{S}_b$ の対応する項と比較することにより、SOC 起因 ZFS テンソル $D_{ab}^{\text{SOC}}$ の微視的 2次摂動表式 が直接得られる：

\[D_{ab}^{\text{SOC}} = -\sum_{n \neq 0} \frac{\langle 0 \vert \sum_{k=1}^N \xi_k l_{k, a} \vert n \rangle \langle n \vert \sum_{m=1}^N \xi_m l_{m, b} \vert 0 \rangle}{E_n - E_0}\]

ただし、記号の簡略化として $\xi_k = \xi(r_k), \xi_m = \xi(r_m)$ と表記している。

テンソルの対称化とトレースレス化 一般に、上式で得られる微視的テンソル $D_{ab}^{\text{SOC}}$ はエルミート基底のもとで実対称化 $\frac{1}{2}(D_{ab}^{\text{SOC}} + D_{ba}^{\text{SOC}})$ され、かつ第 1 節で論じたトレースレス条件 $\text{Tr}(\mathbf{D}^{\text{SOC}}) = 0$ を満たすように、平均対角成分を引く標準化手続きが適用される：

\[D_{ab}^{\text{SOC, final}} = \frac{1}{2}\left(D_{ab}^{\text{SOC}} + D_{ba}^{\text{SOC}}\right) - \frac{1}{3} \delta_{ab} \text{Tr}(\mathbf{D}^{\text{SOC}})\]

3.5 励起状態の分類（同種多重項 vs 異種多重項）と物理的含意

導出された $D_{ab}^{\text{SOC}}$ の微視的表現式において、中間励起状態 $\vert n \rangle$ のスピン多重度（Spin Multiplicity）に応じた詳細な物理的分類が極めて重要となる。

1. 同種多重項励起 (Spin-Allowed / Same-Spin States: $\Delta S = 0$) 基底状態と同じ全スピン量子数 $S$ を持つ励起状態 $\vert n \rangle$。この項ではスピン演算子の対角行列要素が寄与し、主に最高占有軌道（HOMO）や一重占有軌道（SOMO）から最低空軌道（LUMO）への $d$-$d$ 遷移や配位子-金属間電荷移動（LMCT/MLCT）が関与する。
2. 異種多重項励起 (Spin-Forbidden / Spin-Mixed States: $\Delta S = \pm 1$) 基底状態と異なる全スピン量子数（例えば三重項基底 $S=1$ に対する一重項 $\Delta S = -1$ や五重項 $\Delta S = +1$）を持つ励起状態。SOC ハミルトニアンのスピン昇降演算子 $s_{\pm} = s_x \pm i s_y$ によって基底状態とスピンの異なる励起状態が量子力学的に混合し、ZFS に重大な寄与をもたらす。

配位子場分裂と $D_{ab}^{\text{SOC}}$ の異方性

エネルギー分母 $E_n - E_0$ は結晶場分裂（配位子場分裂 $\Delta_{dd}$）の大きさに直結している。

• 励起エネルギーギャップ $E_n - E_0$ が小さい配位子場環境（弱場錯体）では、$D_{ab}^{\text{SOC}}$ の値は極めて大きくなり、巨大な ZFS を示す。
• 逆に強場配位子環境や大ギャップ錯体では分母が大きくなるため SOC 寄与は抑制される。
• また、SOC テンソルの符号は、HOMO-LUMO 間の軌道角運動量行列要素 $\langle 0 \vert l_a \vert n \rangle$ の対角・非対角成分の選択規則（例：$d_{x^2-y^2} \leftrightarrow d_{xy}$ 遷移は $l_z$ で結合し $D_{zz} < 0$ を与え、$d_{xz, yz} \leftrightarrow d_{x^2-y^2}$ 遷移は $l_x, l_y$ で結合し $D_{zz} > 0$ を与える）によって厳密に制御される。

4. 密度汎関数理論（DFT）における ZFS テンソルの定式化

前節までの微視的導出（多電子波動関数および 2次摂動論）を踏まえ、現代の第一原理電子状態計算（DFT）において ZFS テンソル $\mathbf{D}$ がどのように実装され、具体的に計算されているかを定式化する。DFT では多電子波動関数 $\Psi$ の代わりに電子スピン密度 $\rho_\alpha(\mathbf{r}), \rho_\beta(\mathbf{r})$ および Kohn-Sham（KS）スピン軌道 $\psi_k^\sigma(\mathbf{r})$ を用いて SS 寄与および SOC 寄与を表現する。

4.1 スピン-スピン双極子相互作用（SS 寄与）の DFT 定式化

(1) 2粒子スピン対密度と実空間ダブル積分表現

第 2 節で導出した SS 寄与の微視的テンソル式は、実空間におけるスピン対密度（2-particle spin pair density） $\rho_s(\mathbf{r}_1, \mathbf{r}_2)$ を用いて以下のように定式化される：

\[D_{ab}^{\text{SS}} = -\frac{\mu_0 g_e^2 \mu_B^2}{8\pi} \frac{1}{4S(2S-1)} \iint d^3\mathbf{r}_1 d^3\mathbf{r}_2 \rho_s(\mathbf{r}_1, \mathbf{r}_2) \frac{3 r_{12, a} r_{12, b} - r_{12}^2 \delta_{ab}}{r_{12}^5}\]

ここで、$\mathbf{r}{12} = \mathbf{r}_1 - \mathbf{r}_2$、$r{12} = \vert \mathbf{r}{12} \vert$ であり、$r{12, a}$ は相対位置ベクトルの第 $a$ デカルト成分（$a, b \in {x, y, z}$）である。

スピン対密度 $\rho_s(\mathbf{r}_1, \mathbf{r}_2)$ は、位置 $\mathbf{r}_1$ と $\mathbf{r}_2$ における電子スピン間の相関を表す 2粒子密度行列のスピン非対称成分であり、厳密には以下のように定義される：

\[\rho_s(\mathbf{r}_1, \mathbf{r}_2) = \rho_{\alpha\alpha}(\mathbf{r}_1, \mathbf{r}_2) + \rho_{\beta\beta}(\mathbf{r}_1, \mathbf{r}_2) - \rho_{\alpha\beta}(\mathbf{r}_1, \mathbf{r}_2) - \rho_{\beta\alpha}(\mathbf{r}_1, \mathbf{r}_2)\]

(2) Kohn-Sham 行列式近似（McWeeny の定式化と Neese らの近似）

単一の Kohn-Sham Slater 行列式近似のもとでは、2粒子スピン対密度は 1粒子 Kohn-Sham スピン軌道 $\psi_i^\sigma(\mathbf{r})$ （$\sigma \in {\alpha, \beta}$）の直積と交換項（Exchange term）に分解される。

スピン密度 $\rho_s(\mathbf{r}) = \rho_\alpha(\mathbf{r}) - \rho_\beta(\mathbf{r})$ を用いた場合、McWeeny の 2粒子スピン対密度は、Direct 項（クーロン型）と Exchange 項（交換型）の和として次のように分解される：

\[\rho_s(\mathbf{r}_1, \mathbf{r}_2) = \rho_s(\mathbf{r}_1) \rho_s(\mathbf{r}_2) - \rho_{s, \text{exch}}(\mathbf{r}_1, \mathbf{r}_2)\]

ここで、交換型スピン対密度 $\rho_{s, \text{exch}}(\mathbf{r}_1, \mathbf{r}_2)$ は占有 Kohn-Sham 軌道の積で表される：

\[\rho_{s, \text{exch}}(\mathbf{r}_1, \mathbf{r}_2) = \sum_{i \in \text{occ}}^\alpha \sum_{j \in \text{occ}}^\alpha \psi_i^\alpha(\mathbf{r}_1) \psi_j^\alpha(\mathbf{r}_2) \psi_j^\alpha(\mathbf{r}_1) \psi_i^\alpha(\mathbf{r}_2) + \sum_{i \in \text{occ}}^\beta \sum_{j \in \text{occ}}^\beta \psi_i^\beta(\mathbf{r}_1) \psi_j^\beta(\mathbf{r}_2) \psi_j^\beta(\mathbf{r}_1) \psi_i^\beta(\mathbf{r}_2) - 2 \sum_{i \in \text{occ}}^\alpha \sum_{j \in \text{occ}}^\beta \psi_i^\alpha(\mathbf{r}_1) \psi_j^\beta(\mathbf{r}_2) \psi_i^\alpha(\mathbf{r}_2) \psi_j^\beta(\mathbf{r}_1)\]

実際の量子化学計算パッケージ（ORCA 等）では、占有軌道 $\psi_i^\sigma$ を用いて 2電子双極子積分として計算される：

\[D_{ab}^{\text{SS}} = \frac{\mu_0 g_e^2 \mu_B^2}{16\pi S(2S-1)} \sum_{i, j \in \text{occ}} \left[ \langle ij \vert \hat{g}_{ab} \vert ij \rangle - \langle ij \vert \hat{g}_{ab} \vert ji \rangle \right]\]

ここで $\hat{g}{ab}(\mathbf{r}{12}) = \frac{r_{12}^2 \delta_{ab} - 3 r_{12, a} r_{12, b}}{r_{12}^5}$ は 2電子空間双極子演算子である。この表現により、第一原理 DFT 計算において分子構造・スピン密度分布から精度良く SS 寄与が決定される。

4.2 スピン-軌道相互作用（SOC 寄与）の DFT 定式化

SOC に起因する ZFS テンソル $D_{ab}^{\text{SOC}}$ は、励起状態の足し込みを必要とする 2次摂動形式を持つため、DFT の枠組みでは主に 2 つのアプローチで定式化される。

(1) Pederson-Khanna (PK) アプローチ

Pederson と Khanna（1999）によって提案された PK 法は、Kohn-Sham 1電子スピン軌道とその固有エネルギー $\epsilon_k^\sigma$ を用いて 2次摂動論を直積励起に適用する実用的な表現手法である。

占有 Kohn-Sham スピン軌道 $\psi_i^\sigma$ （エネルギー $\epsilon_i^\sigma$）および非占有（仮想）スピン軌道 $\psi_a^{\sigma’}$ （エネルギー $\epsilon_a^{\sigma’}$）を用いると、PK 法における ZFS テンソル成分 $D_{ab}^{\text{SOC}}$ は次のように表される：

\[D_{ab}^{\text{SOC}} = -\frac{1}{4S^2} \sum_{i, \sigma} \sum_{a, \sigma'} \frac{\langle \psi_i^\sigma \vert \hat{L}_a \vert \psi_a^{\sigma'} \rangle \langle \psi_a^{\sigma'} \vert \hat{L}_b \vert \psi_i^\sigma \rangle}{\epsilon_a^{\sigma'} - \epsilon_i^\sigma}\]

ここで、$\hat{L}a = \sum_k \xi(r_k) l{k, a}$ はスピン-軌道結合平均場演算子の第 $a$ デカルト成分である。

PK 法の主な特徴と限界：

• 計算効率の高さ: 1電子軌道エネルギー差 $\epsilon_a^{\sigma’} - \epsilon_i^\sigma$ と 1電子軌道角運動量行列要素のみから計算可能であり、追加のセルフコンシステント応答計算を必要としない。
• 物理的限界: 分母の $4S^2$ や和の構造は高スピン系への簡便なマッピングを与えるが、電子相関応答（Self-Consistent Response）や異なるスピン多重度間の状態混合の厳密な扱いには一定の制限がある。

(2) 結合摂動 DFT (CP-DFT / Linear Response) アプローチ

Neese らに代表される Coupled-Perturbed DFT (CP-DFT) / 線形応答理論アプローチでは、スピン-軌道相互作用摂動に対する電子密度の自己一貫的応答（Self-Consistent Field 応答）および多体電子状態間のスピン多重度変化（同種多重項・異種多重項）を厳密に考慮する。

CP-DFT における SOC テンソルは、Kohn-Sham 軌道応答行列（CP-KS 方程式の解）を介して以下のように定式化される：

\[D_{ab}^{\text{SOC}} = D_{ab}^{\text{SOC, (som-som)}} + D_{ab}^{\text{SOC, (som-dom)}} + D_{ab}^{\text{SOC, (dom-som)}}\]

ここで SOMO（一重占有軌道）および DOMO（二重占有軌道）間の 1 電子励起項と、自己一貫的 Kohn-Sham ポテンシャル応答演算子の効果が含まれる。これにより、遷移金属錯体や大環状有機ラジカルにおいて実験 EPR スペクトルと極めて高い一致を示す精緻な ZFS 予測が可能となっている。

5. 結論とまとめ

本稿では、スピン量子数 $S \ge 1$ を有する量子スピン系における ゼロ磁場分裂 (Zero-Field Splitting: ZFS) について、現象論的テンソル定式化から微視的起源（スピン-スピン双極子相互作用 SS およびスピン-軌道相互作用 SOC）、そして第一原理密度汎関数理論 (DFT) における定式化までを理論的に包括解説した。

5.1 SS 寄与 vs SOC 寄与の定量的比較と優位性

物理的起源の解析から明らかなように、ZFS テンソル $\mathbf{D} = \mathbf{D}^{\text{SS}} + \mathbf{D}^{\text{SOC}}$ を支配する 2大効果は、対象とする物質の電子構造と構成元素の原子番号 $Z$ に依存して大きく異なる。

物理的寄与	メカニズム	スケーリング / 異方性	主な支配的対象系
スピン-スピン双極子 (SS 寄与)	2電子スピン磁気双極子間の直接相互作用	電子間距離 $1/r_{12}^3$、実空間スピン密度の幾何学的異方性	有機三重項分子、ビラジカル、軽元素系（NV センター等）
スピン-軌道相互作用 (SOC 寄与)	軌道角運動量とスピンの結合（2次摂動効果）	原子番号 $Z^4$（重原子効果）、配位子場分裂 $\Delta_{dd}$ の逆数	3d/4f 遷移金属錯体、単分子磁石、重元素を含む量子スピン系

1. 軽元素有機ラジカル・三重項状態・NV センター: 炭素・窒素・酸素などの軽元素（小 $Z$）からなる系では、SOC 寄与（$\propto Z^4$）が比較的小さいため、実空間のスピン密度分布の幾何学的異方性を反映する SS 寄与が ZFS の主成分となる（例: ダイヤモンド NV センターの $D \approx 2.87 \text{ GHz}$ は主に SS 寄与に由来する）。
2. 遷移金属錯体・希土類錯体 (3d, 4f 系)・単分子磁石: 鉄、コバルト、ニッケルやランタノイド等の重元素を含む系では、強力な SOC 寄与が主導権を握る。励起軌道との小さな配位子場エネルギーギャップ $\Delta_{dd}$ と相まって、極めて巨大な $D$ 値（数百 $\text{cm}^{-1}$ オーダー）および強い磁気異方性をもたらす。

5.2 量子技術および計算化学への発展と結び

ZFS テンソルの定量的理解と第一原理予測は、現代の先進的量子科学技術およびスピン化学において中核的役割を果たしている：

1. 量子情報・量子センシング（固体スピン量子ビット）: ダイヤモンド NV センターや SiC 中のカラーセンターは、$S=1$ 基底状態の ZFS ($D, E$) を利用して、無磁場下でスピン状態の初期化・制御・光学読出し (ODMR) を実現している。ZFS テンソルと局所歪み・温度・電場の結合を利用した超高感度量子センサの開発が進んでいる。
2. 単分子磁石 (Single-Molecule Magnets: SMMs) と分子スパイントロニクス: 容易軸異方性 ($D < 0$) に伴うエネルギー障壁 $U = S^2 \vert D \vert$ を最大化することは、高保持力・高動作温度を持つ分子メモリや量子論的磁気記録デバイス設計の重要テーマである。
3. 第一原理 DFT 計算の役割: 現代の DFT 定式化（2粒子スピン密度表現および CP-DFT）により、未知の錯体や新物質の ZFS テンソルを理論的・先験的に高精度予測することが可能となり、分子設計・物性予測の強力な指針を与えている。

ZFS テンソルは、極微のスピン相互作用から量子マクロ物性を架橋する美しく深遠な物理量であり、今後も量子技術の進展とともにその理論的・実用的重要性を増し続けるであろう。