test

19 minute read

導入

核磁気共鳴（NMR）や電子スピン共鳴（ESR / EPR）、さらには量子ビットや固体スピン（NVセンター等）の制御において、スピン系の時間の進みに伴う緩和現象（Relaxation）を正確に理解することは極めて重要である。

1946年に Felix Bloch によって提唱された ブロッホ方程式（Bloch Equations） は、磁場中の磁化ベクトルの歳差運動と、熱平衡への復帰（縦緩和 $T_1$）および位相相関の減衰（横緩和 $T_2$）を記述する現象論的かつ直感的な理論モデルである。しかし、現代の電子スピン物性や量子スピン系、あるいは $S \ge 1$ を持つ電子スピンにおける ゼロ磁場分裂（Zero-Field Splitting: ZFS）テンソル や複雑なコヒーレンス制御を扱うためには、単なる古典的磁化ベクトルの枠組みを超え、量子力学的な 密度行列（Density Matrix） と微視的な多体相互作用のゆらぎに基づく Bloch-Wangsness-Redfield (BWR) 緩和理論 への拡張が不可欠となる。

本稿では、現象論的な 3次元 Bloch 方程式の導出および回転座標系（Rotating Frame）での厳密解法から出発し、密度行列定式化による Liouville-von Neumann 方程式、微視的相互作用表示での摂動展開とマルコフ近似による $T_1, T_2$ のスペクトル密度関数表示、そして将来的な ZFS テンソル（$\mathbf{S} \cdot \mathbf{D} \cdot \mathbf{S}$）への布石となるスピン異方性の数理まで、一切の数式省略を行わずに解説する。

1. 磁場中のスピン運動と現象論的 Bloch 方程式

1.1 磁気モーメントと古典的トルク方程式

個々のスピン粒子（電子や原子核）は、固有の角運動量（スピン角運動量） $\mathbf{S}$ に比例する磁気モーメント $\boldsymbol{\mu}$ を所有している：

\[\boldsymbol{\mu} = \gamma \mathbf{S}\]

ここで $\gamma$ は磁気回転比（Gyromagnetic ratio）であり、粒子の電荷 $q$ および質量 $m$、 $g$ 因子を用いて $\gamma = g \frac{q}{2m}$ （電子の場合は $q = -e < 0$ より $\gamma = -g_e \mu_B / \hbar < 0$）と表される。

体積 $V$ 内に存在する多数のスピン粒子の統計的平均（巨視的アンサンブル平均）として定義される体積磁化ベクトル $\mathbf{M}$ は次のように与えられる：

\[\mathbf{M} = \frac{1}{V} \sum_{i} \langle \boldsymbol{\mu}_i \rangle\]

外部から磁場 $\mathbf{B}$ を印加すると、磁気モーメント $\boldsymbol{\mu}$ にはトルク $\mathbf{T} = \boldsymbol{\mu} \times \mathbf{B}$ が作用する。ニュートンの運動方程式（あるいは角運動量定理 $d\mathbf{S}/dt = \mathbf{T}$）より、巨視的磁化ベクトル $\mathbf{M}$ の時間変化を表す運動方程式は以下の古典的トルク方程式となる：

\[\frac{d\mathbf{M}}{dt} = \gamma \mathbf{M} \times \mathbf{B}\]

静磁場 $\mathbf{B}_0 = (0, 0, B_0)^T$（$z$ 軸方向）が印加された場合、各成分の方程式は以下のように書き下される：

\[\begin{aligned} \frac{d M_x}{dt} &= \gamma M_y B_0 \\ \frac{d M_y}{dt} &= -\gamma M_x B_0 \\ \frac{d M_z}{dt} &= 0 \end{aligned}\]

この微分方程式の解は、$z$ 軸まわりの右ねじ（または左ねじ）方向の円運動、すなわちラーモア歳差運動（Larmor precession）を表す：

\[\begin{aligned} M_x(t) &= M_x(0) \cos(\omega_0 t) + M_y(0) \sin(\omega_0 t) \\ M_y(t) &= -M_x(0) \sin(\omega_0 t) + M_y(0) \cos(\omega_0 t) \\ M_z(t) &= M_z(0) \end{aligned}\]

ここで角振動数 $\omega_0 = -\gamma B_0$ はラーモア振動数（Larmor frequency）と呼ばれる。純粋なトルク方程式のもとでは、磁化ベクトルの大きさ $\vert \mathbf{M} \vert$ および $z$ 成分 $M_z$ は永久に保存され、横成分 $M_x, M_y$ は周期的に振動し続ける。

1.2 熱平衡状態における巨視的磁化ベクトル

熱力学的平衡状態にあるスピン系を考える。温度 $T$ の熱浴（格子系）と熱平衡にある静磁場 $\mathbf{B}_0 = (0, 0, B_0)^T$ 中のスピン集団では、ゼーマン相互作用ハミルトニアン $\hat{\mathcal{H}}_0 = -\boldsymbol{\mu} \cdot \mathbf{B}_0 = -\gamma B_0 \hat{S}_z = \hbar \omega_0 \hat{S}_z$ に従い、エネルギー準位が分裂する（Zeeman Splitting）。

ボルツマン分布に従うスピンの配向分布により、$z$ 軸方向には有限の熱平衡磁化 $M_0$ が発生する。例えばスピン $S=1/2$ の系における熱平衡磁化 $M_0$ は、全スピン密度を $N_0$ とするとキュリーの法則（あるいはブリュアン関数）により以下のように求められる：

\[M_0 = N_0 \gamma \hbar S B_{S}\!\left( \frac{\gamma \hbar S B_0}{k_B T} \right) \approx \frac{N_0 \gamma^2 \hbar^2 S(S+1)}{3 k_B T} B_0 \quad (\hbar \omega_0 \ll k_B T \text{ の高温度近似})\]

一方、$x, y$ 平面（横方向）においては、熱的ゆらぎによって各スピンの横成分の位相（Phase）が無秩序にランダム化しているため、アンサンブル平均をとると横磁化は完全に相殺される。したがって、熱平衡状態における磁化ベクトル $\mathbf{M}_0$ は次のように一意に定まる：

\[\mathbf{M}_0 = \begin{pmatrix} 0 \\ 0 \\ M_0 \end{pmatrix}\]

古典的トルク方程式 $\frac{d\mathbf{M}}{dt} = \gamma \mathbf{M} \times \mathbf{B}_0$ のみでは、ラジオ波（RF）パルスやマイクロ波パルスによって平衡から崩された磁化 $\mathbf{M}(t)$ が、時間経過とともに再び熱平衡状態 $\mathbf{M}_0$ へ戻るという不可逆な熱力学的緩和現象を説明することができない。

1.3 縦緩和時間 $T_1$（スピン–格子緩和）の物理的メカニズム

スピン系が外部磁場や電磁波パルスからエネルギーを吸収して励起された後、$z$ 成分磁化 $M_z$ が熱平衡値 $M_0$ へ復帰するプロセスを縦緩和（Longitudinal Relaxation）またはスピン–格子緩和（Spin-Lattice Relaxation）と呼ぶ。

この過程における物理的メカニズムの本質は以下の通りである：

  1. エネルギー交換: スピンの $z$ 成分 $M_z$ の変化は、ゼーマンエネルギー $E = -\mathbf{M} \cdot \mathbf{B}_0 = -M_z B_0$ の変化を伴う。したがって、$M_z$ が $M_0$ に戻るためには、スピン系と周囲の環境（結晶格子、溶媒分子、フォノン場など、総称して「格子 (Lattice)」と呼ぶ）との間でエネルギーをやり取りする必要がある。
  2. 熱浴のゆらぎ: 格子の熱運動や分子運動によって生じる局所微視的磁場のゆらぎのうち、ラーモア振動数 $\omega_0$ に等しい周波数成分がスピンの誘導放出・吸収を引き起こし、スピンの反転（Flip）が生じる。
  3. 現象論的減衰方程式: Felix Bloch は、この $M_z$ の熱平衡値 $M_0$ への復帰速度が、現在のずれ $(M_z - M_0)$ に比例すると仮定した：
\[\left( \frac{d M_z}{dt} \right)_{\text{relax}} = -\frac{M_z - M_0}{T_1}\]

ここで $T_1$ は縦緩和時間（Longitudinal relaxation time）と呼ばれる特徴的自発時間定数である。この微分方程式を単独で解くと、以下の指数関数的な回復挙動が得られる：

\[M_z(t) = M_0 + \left( M_z(0) - M_0 \right) e^{-t / T_1}\]

例えば、完全な反転パルス（$\pi$ パルス）によって $M_z(0) = -M_0$ とされた反転回復（Inversion Recovery）実験では、$M_z(t) = M_0 (1 - 2e^{-t/T_1})$ となり、時間 $t = T_1 \ln 2 \approx 0.693 T_1$ で $M_z = 0$ を通過して熱平衡 $M_0$ へ回復する。

1.4 横緩和時間 $T_2$（スピン–スピン緩和）の物理的メカニズム

一方、$x, y$ 平面内の横磁化 $M_x, M_y$ が減衰して $0$ （熱平衡値）へ収束するプロセスを横緩和（Transverse Relaxation）またはスピン–スピン緩和（Spin-Spin Relaxation）と呼ぶ。

この過程の物理的メカニズムは縦緩和とは本質的に異なる：

  1. 局所磁場の不均一性と位相拡散（Dephasing）: 試料内部の各スピンは、外部磁場 $\mathbf{B}_0$ に加えて、近傍の他のスピン双極子が作る局所磁場や化学環境のゆらぎによる微小な磁場不均一性 $\delta B_z(\mathbf{r})$ を経験している。これにより、場所やスピンごとに局所ラーモア振動数が $\omega(\mathbf{r}) = \omega_0 + \gamma \delta B_z(\mathbf{r})$ と僅かに異なる。
  2. 位相の一致の喪失: パルス照射直後は横平面内で揃っていたスピンの位相（Coherence）が、時間の経過とともに徐々にバラバラになり（位相拡散）、横磁化ベクトル和の絶対値 $\sqrt{M_x^2 + M_y^2}$ が失われる。
  3. エネルギー非保存プロセス（Pure Dephasing）: 近傍スピン同士のエネルギー交換（スピンフリップ・フロップ運動 $\hat{S}{1+} \hat{S}{2-} + \hat{S}{1-} \hat{S}{2+}$）や局所磁場ゆらぎによる位相拡散は、スピン系全体のゼーマンエネルギーを変更することなく（格子とのエネルギー交換なしに）横コヒーレンスのみを破壊する（Pure Dephasing, 時間定数 $T_2^*$ または $T_{\phi}$）。
  4. $T_1$ 寄与との関係: スピンが縦緩和を起こしてエネルギー順位間を遷移する（寿命が有限である）こと自体も、不確定性原理 $\Delta E \Delta t \sim \hbar$ を通じて横コヒーレンスの減衰を引き起こす。したがって、総和としての横緩和率 $1/T_2$ は、純粋位相拡散率 $1/T_{\phi}$ と縦緩和からの寄与 $1/(2T_1)$ の和として以下のように表される：
\[\frac{1}{T_2} = \frac{1}{2T_1} + \frac{1}{T_{\phi}}\]

これによって、一般に不等式 $T_2 \le 2 T_1$（多くの常磁性体や液体NMRでは $T_2 \le T_1$）が常に成立する。

Bloch は横磁化 $M_x, M_y$ の減衰速度を次のように現象論的にモデル化した：

\[\begin{aligned} \left( \frac{d M_x}{dt} \right)_{\text{relax}} &= -\frac{M_x}{T_2} \\ \left( \frac{d M_y}{dt} \right)_{\text{relax}} &= -\frac{M_y}{T_2} \end{aligned}\]

ここで $T_2$ は横緩和時間（Transverse relaxation time）である。

1.5 現象論的 3次元 Bloch 方程式の完全な記述

古典的歳差運動の項 $\gamma (\mathbf{M} \times \mathbf{B})$ と、現象論的緩和項（$T_1, T_2$）を重畳することにより、一般の磁場 $\mathbf{B}(t) = (B_x(t), B_y(t), B_z(t))^T$ のもとでの現象論的 3次元 Bloch 方程式（Phenomenological 3D Bloch Equations）の完全な連立常微分方程式系が得られる：

\[\frac{d M_x}{dt} = \gamma (\mathbf{M} \times \mathbf{B})_x - \frac{M_x}{T_2}\] \[\frac{d M_y}{dt} = \gamma (\mathbf{M} \times \mathbf{B})_y - \frac{M_y}{T_2}\] \[\frac{d M_z}{dt} = \gamma (\mathbf{M} \times \mathbf{B})_z - \frac{M_z - M_0}{T_1}\]

外積成分 $\mathbf{M} \times \mathbf{B} = (M_y B_z - M_z B_y, M_z B_x - M_x B_z, M_x B_y - M_y B_x)^T$ を各成分ごとに明示的に展開すると、以下の通りである：

\[\begin{aligned} \frac{d M_x}{dt} &= \gamma (M_y B_z - M_z B_y) - \frac{M_x}{T_2} \\ \frac{d M_y}{dt} &= \gamma (M_z B_x - M_x B_z) - \frac{M_y}{T_2} \\ \frac{d M_z}{dt} &= \gamma (M_x B_y - M_y B_x) - \frac{M_z - M_0}{T_1} \end{aligned}\]

ベクトル・行列形式で記述すると、次のようになる：

\[\frac{d}{dt} \begin{pmatrix} M_x \\ M_y \\ M_z \end{pmatrix} = \begin{pmatrix} -\frac{1}{T_2} & \gamma B_z & -\gamma B_y \\ -\gamma B_z & -\frac{1}{T_2} & \gamma B_x \\ \gamma B_y & -\gamma B_x & -\frac{1}{T_1} \end{pmatrix} \begin{pmatrix} M_x \\ M_y \\ M_z \end{pmatrix} + \begin{pmatrix} 0 \\ 0 \\ \frac{M_0}{T_1} \end{pmatrix}\]

この方程式系は、常磁性体や核磁気スピン集団における古典的磁化ベクトルの動的応答および緩和挙動の第一歩となる基礎方程式である。

2. 高周波（RF）磁場と回転座標系における Bloch 方程式の解

前節の古典的 Bloch 方程式は静磁場 $\mathbf{B}_0 = (0, 0, B_0)^T$ のみの挙動を扱ったが、実際の NMR / ESR やスピン制御実験では、横平面内に時間周期変動する高周波（RF）磁場 $\mathbf{B}_1(t)$ を印加してスピン状態の操作や共鳴吸収を観測する。

本節では、静止座標系（Laboratory Frame）から回転振動する回転座標系（Rotating Frame）への階層的ユニタリ座標変換を導入し、時間依存項を消去した有効磁場 $\mathbf{B}_{\text{eff}}$ とラビ振動、および定常状態における定常解（吸収・分散ローレンツ線形、RF 飽和因子）を厳密に導出する。

2.1 高周波磁場（RF磁場）の導入と全磁場

$z$ 軸方向に平行な静磁場 $\mathbf{B}_0 = (0, 0, B_0)^T$ に加え、$xy$ 平面内で角振動数 $\omega$ で円偏光回転する高周波磁場（RF 磁場） $\mathbf{B}_1(t)$ を以下のように与える：

\[\mathbf{B}_1(t) = \begin{pmatrix} B_1 \cos\omega t \\ -B_1 \sin\omega t \\ 0 \end{pmatrix} = B_1 (\cos\omega t \, \hat{\mathbf{x}} - \sin\omega t \, \hat{\mathbf{y}})\]

したがって、スピン系に作用する全磁場 $\mathbf{B}(t) = \mathbf{B}_0 + \mathbf{B}_1(t)$ は次のように表される：

\[\mathbf{B}(t) = \begin{pmatrix} B_1 \cos\omega t \\ -B_1 \sin\omega t \\ B_0 \end{pmatrix}\]

【注意】直線偏光磁場との関係 実験で用いられる高周波磁場は、多くの場合 $x$ 軸方向の直線偏光磁場 $2 B_1 \cos\omega t \, \hat{\mathbf{x}}$ である。これは時計回りおよび反時計回りに回転する 2 つの円偏光成分の和として表される： \(2 B_1 \cos\omega t \, \hat{\mathbf{x}} = B_1 (\cos\omega t \, \hat{\mathbf{x}} - \sin\omega t \, \hat{\mathbf{y}}) + B_1 (\cos\omega t \, \hat{\mathbf{x}} + \sin\omega t \, \hat{\mathbf{y}})\) ラーモア歳差運動と同方向に回転する成分（上記第1項）のみが強力に共鳴相互作用し、逆方向に回転する成分（第2項）は共鳴条件から大きく離れているため、高次な幾何学的位相シフト（Bloch-Siegert シフト）を除いて無視できる（回転波近似: Rotating Wave Approximation, RWA）。

2.2 回転座標系（Rotating Frame）への座標変換

磁場 $\mathbf{B}(t)$ が陽に時間依存するため、静止座標系（Laboratory Frame）における Bloch 方程式は時間依存係数を持つ複雑な連立微分方程式となる。これを解決するため、$z$ 軸のまわりに角振動数 $\omega$ で回転する標準的な直交座標系（回転座標系, Rotating Frame）を導入する。

静止座標系の単位基底ベクトル $(\hat{\mathbf{x}}, \hat{\mathbf{y}}, \hat{\mathbf{z}})$ に対し、回転座標系の単位基底ベクトル $(\hat{\mathbf{x}}’, \hat{\mathbf{y}}’, \hat{\mathbf{z}}’)$ は次のように定義される：

\[\begin{aligned} \hat{\mathbf{x}}' &= \cos\omega t \, \hat{\mathbf{x}} - \sin\omega t \, \hat{\mathbf{y}} \\ \hat{\mathbf{y}}' &= \sin\omega t \, \hat{\mathbf{x}} + \cos\omega t \, \hat{\mathbf{y}} \\ \hat{\mathbf{z}}' &= \hat{\mathbf{z}} \end{aligned}\]

静止座標系における磁化ベクトル $\mathbf{M}{\text{lab}} = (M_x, M_y, M_z)^T$ と回転座標系における磁化ベクトル $\mathbf{M}{\text{rot}} = (u, v, M_z)^T$ の変換関係は、$z$ 軸周りの回転行列 $\mathbf{R}_z(\omega t)$ を用いて以下のように記述される：

\[\mathbf{M}_{\text{rot}} = \mathbf{R}_z(\omega t) \mathbf{M}_{\text{lab}}\] \[\begin{pmatrix} u \\ v \\ M_z \end{pmatrix} = \begin{pmatrix} \cos\omega t & -\sin\omega t & 0 \\ \sin\omega t & \cos\omega t & 0 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} M_x \\ M_y \\ M_z \end{pmatrix}\]

逆に、静止座標系への逆変換は以下の通りである：

\[\mathbf{M}_{\text{lab}} = \mathbf{R}_z(-\omega t) \mathbf{M}_{\text{rot}} = \begin{pmatrix} \cos\omega t & \sin\omega t & 0 \\ -\sin\omega t & \cos\omega t & 0 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} u \\ v \\ M_z \end{pmatrix}\]

古典力学におけるベクトルの回転運動論より、静止座標系での時間微分 $(d\mathbf{M}/dt){\text{lab}}$ と回転座標系での時間微分 $(d\mathbf{M}/dt){\text{rot}}$ の間には、回転角速度ベクトル $\boldsymbol{\omega} = (0, 0, \omega)^T = \omega \hat{\mathbf{z}}$ を用いて次の関係式が成立する：

\[\left( \frac{d\mathbf{M}}{dt} \right)_{\text{lab}} = \left( \frac{d\mathbf{M}}{dt} \right)_{\text{rot}} + \boldsymbol{\omega} \times \mathbf{M}\]

これをトルク方程式 $(d\mathbf{M}/dt)_{\text{lab}} = \gamma \mathbf{M} \times \mathbf{B}$ に代入すると、回転座標系における運動方程式が得られる：

\[\left( \frac{d\mathbf{M}}{dt} \right)_{\text{rot}} = \gamma \mathbf{M} \times \mathbf{B} - \boldsymbol{\omega} \times \mathbf{M} = \gamma \mathbf{M} \times \left( \mathbf{B} + \frac{\boldsymbol{\omega}}{\gamma} \right)\]

2.3 有効磁場 $\mathbf{B}_{\text{eff}}$・ラーモア共鳴条件・ラビ振動数

回転座標系における巨視的運動方程式は、実効的なトルク方程式として以下のようにシンプルにまとめられる：

\[\left( \frac{d\mathbf{M}}{dt} \right)_{\text{rot}} = \gamma \mathbf{M}_{\text{rot}} \times \mathbf{B}_{\text{eff}}\]

ここで、有効磁場（Effective Magnetic Field） $\mathbf{B}_{\text{eff}}$ は、回転座標系での全磁場に回転に伴う仮想磁場（慣性磁場） $-\boldsymbol{\omega}/\gamma$ を加えたものである：

\[\mathbf{B}_{\text{eff}} = \mathbf{B}_{\text{rot}} + \frac{\boldsymbol{\omega}}{\gamma}\]

回転座標系において、高周波磁場 $\mathbf{B}1(t)$ は常に $x’$ 軸正方向に静止して見える（$\mathbf{B}{\text{rot}} = (B_1, 0, B_0)^T$）。したがって、有効磁場の各成分は時間変化を含まず完全に定数となる：

\[\mathbf{B}_{\text{eff}} = \begin{pmatrix} B_1 \\ 0 \\ B_0 - \frac{\omega}{\gamma} \end{pmatrix} = \begin{pmatrix} B_1 \\ 0 \\ \frac{\Delta\omega}{\gamma} \end{pmatrix}\]

ここで $\Delta\omega$ はオフ・レゾナンス（Off-resonance parameter / デチューニング）であり、静磁場による固有ラーモア振動数 $\omega_0 = -\gamma B_0$ と高周波磁場の照射振動数 $\omega$ の差分として定義される。

回転波近似（RWA）のもとで、ラーモア歳差運動と同方向に回転するフレームを取ることで、有効磁場の $z’$ 成分はデチューニング $\Delta\omega = \omega - \omega_0$ を用いて以下のように表される：

\[B_{z, \text{eff}} = \frac{\Delta\omega}{\gamma}\]

【共鳴条件とラビ振動】 照射角振動数 $\omega$ がスピン固有のラーモア振動数 $\omega_0$ に一致するとき（ラーモア共鳴条件: Larmor Resonance Condition $\omega = \omega_0$）、オフゾナンスは完全消滅し $\Delta\omega = 0$ となる。

このとき有効磁場の $z’$ 成分は $0$ となり、有効磁場は横軸 $x’$ 方向の成分のみとなる： \(\mathbf{B}_{\text{eff}} = \begin{pmatrix} B_1 \\ 0 \\ 0 \end{pmatrix}\) したがって、回転座標系において磁化ベクトル $\mathbf{M}$ は $x’$ 軸のまわりに角振動数 $\omega_1$ で倒れ込むような歳差運動（ラビ振動, Rabi Oscillation）を行う。この特徴的振動数 $\omega_1$ はラビ振動数（Rabi Frequency）と呼ばれる： \(\omega_1 = -\gamma B_1\)

2.4 回転座標系における Bloch 方程式と $u, v$ 成分

回転座標系における磁化ベクトルの成分 $(u, v, M_z)^T$ において、

  • $u$: RF 磁場 $\mathbf{B}_1(t)$ と同位相（In-phase）の横磁化成分（分散成分 / Dispersion component）
  • $v$: RF 磁場 $\mathbf{B}_1(t)$ と 90° 位相のずれた（Out-of-phase / Quadrature）横磁化成分（吸収成分 / Absorption component）

トルク項 $\gamma (\mathbf{M}{\text{rot}} \times \mathbf{B}{\text{eff}})$ の外積計算を実行する：

\[\mathbf{M}_{\text{rot}} \times \mathbf{B}_{\text{eff}} = \begin{pmatrix} u \\ v \\ M_z \end{pmatrix} \times \begin{pmatrix} B_1 \\ 0 \\ \frac{\Delta\omega}{\gamma} \end{pmatrix} = \begin{pmatrix} v \frac{\Delta\omega}{\gamma} \\ M_z B_1 - u \frac{\Delta\omega}{\gamma} \\ -v B_1 \end{pmatrix}\]

これに各成分の現象論的緩和項（$-u/T_2, -v/T_2, -(M_z - M_0)/T_1$）を加算すると、回転座標系における現象論的 Bloch 方程式が以下の自立系（時間非依存な連立微分方程式）として得られる：

\[\frac{du}{dt} = \Delta\omega v - \frac{u}{T_2}\] \[\frac{dv}{dt} = -\Delta\omega u + \gamma B_1 M_z - \frac{v}{T_2}\] \[\frac{d M_z}{dt} = -\gamma B_1 v - \frac{M_z - M_0}{T_1}\]

ラビ振動数 $\omega_1 = -\gamma B_1$（すなわち $\gamma B_1 = -\omega_1$）を用いると、上式は以下のように表現することもできる：

\[\begin{aligned} \frac{du}{dt} &= \Delta\omega v - \frac{u}{T_2} \\ \frac{dv}{dt} &= -\Delta\omega u - \omega_1 M_z - \frac{v}{T_2} \\ \frac{d M_z}{dt} &= \omega_1 v - \frac{M_z - M_0}{T_1} \end{aligned}\]

2.5 定常状態解（Steady-State Solutions）の厳密導出

連続波（Continuous Wave: CW）実験のように、一定振幅の RF 磁場を十分長時間照射し続けると、スピン系は励起と緩和が釣り合う定常状態（Steady State）に達する。

定常状態においては、時間の微係数がすべてゼロとなる：

\[\frac{du}{dt} = 0, \quad \frac{dv}{dt} = 0, \quad \frac{d M_z}{dt} = 0\]

これにより、以下の代数方程式系が得られる：

\[\begin{aligned} \Delta\omega v - \frac{u}{T_2} &= 0 \quad \implies \quad u = \Delta\omega T_2 v \qquad \text{--- (1)} \\ -\Delta\omega u + \gamma B_1 M_z - \frac{v}{T_2} &= 0 \qquad \text{--- (2)} \\ -\gamma B_1 v - \frac{M_z - M_0}{T_1} &= 0 \quad \implies \quad M_z = M_0 - \gamma B_1 T_1 v \qquad \text{--- (3)} \end{aligned}\]

式 (1) および式 (3) を式 (2) に代入すると、$v$ のみに関する代数方程式が得られる：

\[-\Delta\omega (\Delta\omega T_2 v) + \gamma B_1 (M_0 - \gamma B_1 T_1 v) - \frac{v}{T_2} = 0\]

$v$ について同類項をまとめると：

\[v \left[ \frac{1}{T_2} + (\Delta\omega)^2 T_2 + \gamma^2 B_1^2 T_1 \right] = \gamma B_1 M_0\]

両辺に $T_2$ を乗じる：

\[v \left[ 1 + (\Delta\omega)^2 T_2^2 + \gamma^2 B_1^2 T_1 T_2 \right] = \gamma B_1 M_0 T_2\]

したがって、吸収成分 $v(\Delta\omega)$ の定常解は以下の通り決定される：

\[v(\Delta\omega) = \frac{\gamma B_1 M_0 T_2}{1 + (\Delta\omega)^2 T_2^2 + \gamma^2 B_1^2 T_1 T_2}\]

この $v(\Delta\omega)$ を式 (1) および式 (3) に代入することにより、分散成分 $u(\Delta\omega)$ および縦磁化成分 $M_z(\Delta\omega)$ の定常解が厳密に求められる：

\[u(\Delta\omega) = \frac{\gamma B_1 M_0 \Delta\omega T_2^2}{1 + (\Delta\omega)^2 T_2^2 + \gamma^2 B_1^2 T_1 T_2}\] \[M_z(\Delta\omega) = M_0 \frac{1 + (\Delta\omega)^2 T_2^2}{1 + (\Delta\omega)^2 T_2^2 + \gamma^2 B_1^2 T_1 T_2}\]

2.6 ローレンツ線形（Lorentzian Lineshapes）と RF 電力飽和因子

導出された定常解 $u(\Delta\omega), v(\Delta\omega)$ は、スペクトルの基本線形であるローレンツ線形（Lorentzian Lineshape）および RF 磁場強度による非線形な飽和効果（Saturation Effect）を明確に表している。

ここで、無次元の RF 電力飽和因子（RF Power Saturation Factor） $S_{\text{sat}}$ を次のように定義する：

\[S_{\text{sat}} = \frac{1}{1 + \gamma^2 B_1^2 T_1 T_2}\]

飽和因子 $S_{\text{sat}}$ を用いると、吸収成分 $v(\Delta\omega)$ および分散成分 $u(\Delta\omega)$ は以下のように簡潔に表現される：

\[v(\Delta\omega) = \gamma B_1 M_0 T_2 \, \frac{S_{\text{sat}}}{1 + (\Delta\omega)^2 T_2^2 S_{\text{sat}}}\] \[u(\Delta\omega) = \gamma B_1 M_0 \Delta\omega T_2^2 \, \frac{S_{\text{sat}}}{1 + (\Delta\omega)^2 T_2^2 S_{\text{sat}}}\]

1. 非飽和極限（弱 RF 磁場近似: $\gamma^2 B_1^2 T_1 T_2 \ll 1$）

RF 磁場が十分に弱いとき（$S_{\text{sat}} \approx 1$）、定常解は以下の標準的なローレンツ関数に帰着する：

\[v(\omega) = \gamma B_1 M_0 \frac{T_2}{1 + (\omega - \omega_0)^2 T_2^2}\] \[u(\omega) = \gamma B_1 M_0 \frac{(\omega - \omega_0) T_2^2}{1 + (\omega - \omega_0)^2 T_2^2}\]
  • 吸収信号 $v(\omega)$: 共鳴点 $\omega = \omega_0$ で最大値 $\gamma B_1 M_0 T_2$ をとる対称な山型の曲線。半値全幅（Full Width at Half Maximum: FWHM） $\Delta\omega_{1/2}$ は以下で与えられる： \(\Delta\omega_{1/2} = \frac{2}{T_2}\) すなわち、スペクトルの線幅から横緩和時間 $T_2$ を直接決定することができる。
  • 分散信号 $u(\omega)$: 共鳴点 $\omega = \omega_0$ で $0$ を横切り、共鳴の前後で符号を変える反対称な曲線。

2. パワー飽和現象（Power Saturation）とパワー広がり（Power Broadening）

RF 磁場強度 $B_1$ を極めて大きくすると、$\gamma^2 B_1^2 T_1 T_2 \gg 1$ となり、$S_{\text{sat}} \to 0$ に近づく。

  1. ポピュレーションの飽和: 縦磁化定常値は $M_z(\omega_0) = \frac{M_0}{1 + \gamma^2 B_1^2 T_1 T_2} \to 0$ となる。これは、スピンの RF 吸収速度が格子への縦緩和速度 $1/T_1$ を上回り、エネルギーの上下準位占有数が均等化（ポピュレーション差が消滅）することを意味する。
  2. 吸収強度の飽和限界: 共鳴点での吸収極大値は $v(\omega_0) = \frac{\gamma B_1 M_0 T_2}{1 + \gamma^2 B_1^2 T_1 T_2}$ となる。低 $B_1$ では $B_1$ に比例して強度が増加するが、$B_1 \to \infty$ では逆に $v(\omega_0) \sim 1/B_1 \to 0$ へと減衰する。
  3. パワー広がり（Power Broadening）: 飽和状態における吸収スペクトルの FWHM は以下のように拡大する： \(\Delta\omega_{1/2, \text{sat}} = \frac{2}{T_2} \sqrt{1 + \gamma^2 B_1^2 T_1 T_2} = \frac{2}{T_2 \sqrt{S_{\text{sat}}}}\) これは、実効的な横緩和時間が $T_{2, \text{eff}} = T_2 \sqrt{S_{\text{sat}}}$ に短縮されたことと等価であり、強電磁場下における不可逆なスペクトル線幅拡張を表している。

3. 密度行列定式化と Liouville-von Neumann 方程式

前節までは、巨視的磁化ベクトル $\mathbf{M}$ の運動を記述する古典的な現象論的 Bloch 方程式を扱ってきた。しかし、多体スピン系の量子干渉（コヒーレンス）、ゼロ磁場分裂（ZFS）テンソル、あるいは微視的環境との結合による緩和を定量的に議論するためには、密度演算子（Density Operator / Density Matrix） を用いた量子統計力学的な定式化への展開が必要不可欠である。

本節では、純粋状態のアンサンブルから密度演算子を導出し、その時間発展を支配する Liouville-von Neumann 方程式を提示する。さらに、緩和を考慮しない場合の量子力学的期待値の時間変化が、古典的トルク方程式 $\frac{d\mathbf{M}}{dt} = \gamma \mathbf{M} \times \mathbf{B}$ を厳密に再現することを完全な代数的証明によって示す。

3.1 密度演算子と状態アンサンブル

量子スピン系（スピンアンサンブル）を考える。巨視的試料中に存在する多数のスピンは、必ずしも単一の波動関数で記述される純粋状態（Pure State）にあるとは限らず、各種の微視的状態 $\vert \psi_k(t) \rangle$ が統計的確率 $p_k$ （$\sum_k p_k = 1, p_k \ge 0$）で混合した混合状態（Mixed State）にある。

この状態アンサンブルを統一的に記述するため、密度演算子（Density Operator） $\hat{\rho}(t)$ を次のように定義する：

\[\hat{\rho}(t) = \sum_{k} p_k \vert \psi_k(t) \rangle \langle \psi_k(t) \vert\]

スピン系の任意の完全正規直交基底 ${ \vert a \rangle }$ （例えば $z$ 軸スピン固有状態 $\hat{S}z \vert a \rangle = m_a \vert a \rangle$）を用いて展開すると、密度演算子の行列要素（密度行列） $\rho{ab}(t) = \langle a \vert \hat{\rho}(t) \vert b \rangle$ を用いて以下のように表現される：

\[\hat{\rho}(t) = \sum_{a, b} \rho_{ab}(t) \vert a \rangle \langle b \vert\]

対角成分 $\rho_{aa}$ と非対角成分 $\rho_{ab}$ の物理的意味

密度行列の各成分は、以下のように物理的に完全に異なる概念を担っている：

  1. 対角成分 $\rho_{aa}$ （ポピュレーション / Population）: 対角要素 $\rho_{aa} = \langle a \vert \hat{\rho} \vert a \rangle = \sum_k p_k \vert \langle a \vert \psi_k \rangle \vert^2$ は、アンサンブルの中でスピンが基底状態 $\vert a \rangle$ に存在する確率密度（占有率 / 占有数分布）を表す。
    • $z$ 軸磁化成分 $\langle M_z \rangle$ やエネルギーレベルの熱的占有数分布は、この対角成分（ポピュレーションの差 $\rho_{aa} - \rho_{bb}$）によって直接決定される。
  2. 非対角成分 $\rho_{ab} \, (a \neq b)$ （コヒーレンス / Coherence）: 非対角要素 $\rho_{ab} = \langle a \vert \hat{\rho} \vert b \rangle = \sum_k p_k c_{k, a} c_{k, b}^*$ （ただし $\vert \psi_k \rangle = \sum_a c_{k,a} \vert a \rangle$）は、状態 $\vert a \rangle$ と $\vert b \rangle$ の間の量子力学的位相相関（Phase coherence）および相対的干渉項を表す。
    • 横磁化成分 $\langle M_x \rangle, \langle M_y \rangle$ やスピンの歳差運動の位相は、この非対角成分（コヒーレンス）によって担われる。
    • 各スピンの位相がランダム化すると、アンサンブル平均において非対角成分 $\rho_{ab}$ が相互に打ち消し合って $0$ に減衰する（Dephasing）。

密度演算子の必要十分条件と基本性質

物理的に妥当な密度演算子 $\hat{\rho}$ は、以下の 3 つの数学的性質を満たす必要がある：

  1. エルミート性（Hermiticity）: \(\hat{\rho}^\dagger = \hat{\rho} \quad \implies \quad \rho_{ba} = \rho_{ab}^*\) これにより、任意の物理量演算子 $\hat{A}$ の期待値 $\langle \hat{A} \rangle = \text{Tr}(\hat{\rho} \hat{A})$ が必ず実数となることが保証される。

  2. 規格化条件（Normalization / Trace is unity）: \(\text{Tr}(\hat{\rho}) = \sum_{a} \rho_{aa} = 1\) 全占有確率の総和は常に 1 である。

  3. 半正値性（Positive Semi-definiteness）: 任意の状態で $\langle \phi \vert \hat{\rho} \vert \phi \rangle \ge 0$ であり、純粋度（Purity） $\gamma_{\text{purity}}$ について以下の不等式が成立する： \(\text{Tr}(\hat{\rho}^2) \le 1\) 等号 $\text{Tr}(\hat{\rho}^2) = 1$ はスピン系が完全な純粋状態にある場合にのみ成立し、完全混合状態（無限大温度極限 $T \to \infty$）では $\hat{\rho} = \frac{1}{N} \hat{I}$ （$N$ は状態空間の次元）となり $\text{Tr}(\hat{\rho}^2) = 1/N$ となる。

3.2 Liouville-von Neumann 方程式の導出

アンサンブルを構成する各純粋状態 $\vert \psi_k(t) \rangle$ は、時間依存のシュレーディンガー方程式に従って時間発展する：

\[i\hbar \frac{\partial}{\partial t} \vert \psi_k(t) \rangle = \hat{H} \vert \psi_k(t) \rangle\]

その複素共役（ブラベクトルの時間発展方程式）は、ハミルトニアンのエルミート性 $\hat{H}^\dagger = \hat{H}$ より以下のようになる：

\[-i\hbar \frac{\partial}{\partial t} \langle \psi_k(t) \vert = \langle \psi_k(t) \vert \hat{H}\]

密度演算子 $\hat{\rho}(t) = \sum_k p_k \vert \psi_k(t) \rangle \langle \psi_k(t) \vert$ の時間微分に積の微分法則を適用する（占有確率 $p_k$ は時間非依存であると仮定）：

\[\frac{d\hat{\rho}}{dt} = \sum_{k} p_k \left[ \left( \frac{d}{dt} \vert \psi_k(t) \rangle \right) \langle \psi_k(t) \vert + \vert \psi_k(t) \rangle \left( \frac{d}{dt} \langle \psi_k(t) \vert \right) \right]\]

シュレーディンガー方程式およびその共役を代入する：

\[\frac{d\hat{\rho}}{dt} = \sum_{k} p_k \left[ \left( -\frac{i}{\hbar} \hat{H} \vert \psi_k(t) \rangle \right) \langle \psi_k(t) \vert + \vert \psi_k(t) \rangle \left( \frac{i}{\hbar} \langle \psi_k(t) \vert \hat{H} \right) \right]\] \[\frac{d\hat{\rho}}{dt} = -\frac{i}{\hbar} \hat{H} \left( \sum_{k} p_k \vert \psi_k(t) \rangle \langle \psi_k(t) \vert \right) + \frac{i}{\hbar} \left( \sum_{k} p_k \vert \psi_k(t) \rangle \langle \psi_k(t) \vert \right) \hat{H}\] \[\frac{d\hat{\rho}}{dt} = -\frac{i}{\hbar} \left( \hat{H} \hat{\rho}(t) - \hat{\rho}(t) \hat{H} \right)\]

交換子 $[\hat{H}, \hat{\rho}] \equiv \hat{H}\hat{\rho} - \hat{\rho}\hat{H}$ を定義することにより、量子統計力学における基礎方程式である Liouville-von Neumann 方程式（Liouville-von Neumann Equation） が完璧に導出される：

\[\frac{d\hat{\rho}}{dt} = -\frac{i}{\hbar} [\hat{H}, \hat{\rho}]\]

【Liouvillian スーパーオペレーターと形式解】 演算子に作用して新しい演算子を生成するスーパーオペレーター（Superoperator）として Liouvillian スーパーオペレーター $\hat{\mathcal{L}}$ を以下のように定義する： \(\hat{\mathcal{L}} \hat{\rho} \equiv \frac{1}{\hbar} [\hat{H}, \hat{\rho}]\) これを用いると、Liouville-von Neumann 方程式は $\frac{d\hat{\rho}}{dt} = -i \hat{\mathcal{L}} \hat{\rho}$ とコンパクトに書ける。 ハミルトニアン $\hat{H}$ が時間に依存しない場合、初期密度演算子 $\hat{\rho}(0)$ に対する形式解は次のように与えられる： \(\hat{\rho}(t) = e^{-\frac{i}{\hbar} \hat{H} t} \hat{\rho}(0) e^{\frac{i}{\hbar} \hat{H} t}\)

3.3 緩和のない量子運動方程式と古典トルク方程式の等価性証明

ここで、量子力学的な密度行列定式化（Liouville-von Neumann 方程式）から、第1節で導入した古典的トルク方程式 $\frac{d\mathbf{M}}{dt} = \gamma \mathbf{M} \times \mathbf{B}$ が厳密に再現されることを代数的に証明する。

スピン系に外部磁場 $\mathbf{B} = (B_x, B_y, B_z)^T$ が作用しているとき、スピン演算子ベクトルを $\mathbf{S} = (\hat{S}_x, \hat{S}_y, \hat{S}_z)^T$ とすると、ゼーマンハミルトニアンは以下のようになる：

\[\hat{H} = -\boldsymbol{\mu} \cdot \mathbf{B} = -\gamma \mathbf{B} \cdot \mathbf{S} = -\gamma (B_x \hat{S}_x + B_y \hat{S}_y + B_z \hat{S}_z)\]

巨視的磁化ベクトルの量子力学的アンサンブル期待値 $\langle \mathbf{M} \rangle = (\langle M_x \rangle, \langle M_y \rangle, \langle M_z \rangle)^T$ は、密度演算子 $\hat{\rho}$ とスピン演算子 $\mathbf{S}$ のトレースとして定義される：

\[\langle \mathbf{M} \rangle = \gamma \langle \mathbf{S} \rangle = \gamma \text{Tr}(\hat{\rho}(t) \mathbf{S})\]

期待値の時間変化率 $\frac{d}{dt} \langle \mathbf{M} \rangle$ を計算する：

\[\frac{d}{dt} \langle \mathbf{M} \rangle = \gamma \text{Tr}\left( \frac{d\hat{\rho}}{dt} \mathbf{S} \right)\]

Liouville-von Neumann 方程式 $\frac{d\hat{\rho}}{dt} = -\frac{i}{\hbar} [\hat{H}, \hat{\rho}]$ を代入する：

\[\frac{d}{dt} \langle \mathbf{M} \rangle = -\frac{i\gamma}{\hbar} \text{Tr}\left( [\hat{H}, \hat{\rho}] \mathbf{S} \right)\]

ここで、任意の演算子 $A, B, C$ に対するトレースの循環可換性（Cyclic property of trace）:

\[\text{Tr}([A, B] C) = \text{Tr}((AB - BA)C) = \text{Tr}(ABC - BAC) = \text{Tr}(CAB - ACB) = \text{Tr}(B[C, A]) = \text{Tr}(A[B, C])\]

より、次式が成り立つ：

\[\text{Tr}\left( [\hat{H}, \hat{\rho}] \mathbf{S} \right) = \text{Tr}\left( (\hat{H}\hat{\rho} - \hat{\rho}\hat{H}) \mathbf{S} \right) = \text{Tr}\left( \hat{\rho} \mathbf{S} \hat{H} - \hat{\rho} \hat{H} \mathbf{S} \right) = \text{Tr}\left( \hat{\rho} [\mathbf{S}, \hat{H}] \right)\]

したがって、時間微分は以下の形に変形される：

\[\frac{d}{dt} \langle \mathbf{M} \rangle = -\frac{i\gamma}{\hbar} \text{Tr}\left( \hat{\rho} [\mathbf{S}, \hat{H}] \right)\]

次に、スピン演算子ベクトル $\mathbf{S}$ とハミルトニアン $\hat{H} = -\gamma \mathbf{B} \cdot \mathbf{S}$ の交換子 $[\mathbf{S}, \hat{H}]$ の $i$ 番目の成分を計算する。スピン演算子の基本交換関係：

\[[\hat{S}_i, \hat{S}_j] = i\hbar \sum_{k} \epsilon_{ijk} \hat{S}_k\]

（ただし $\epsilon_{ijk}$ はレビ・チビタの完全反対称テンソル）を用いると、以下のように展開される：

\[[\hat{S}_i, \hat{H}] = \left[ \hat{S}_i, -\gamma \sum_{j} B_j \hat{S}_j \right] = -\gamma \sum_{j} B_j [\hat{S}_i, \hat{S}_j] = -\gamma \sum_{j} B_j \left( i\hbar \sum_{k} \epsilon_{ijk} \hat{S}_k \right)\] \[[\hat{S}_i, \hat{H}] = -i\hbar \gamma \sum_{j, k} \epsilon_{ijk} B_j \hat{S}_k = i\hbar \gamma \sum_{j, k} \epsilon_{ikj} \hat{S}_k B_j = i\hbar \gamma (\mathbf{S} \times \mathbf{B})_i\]

これを成分表示からベクトル表示に戻すと、スピン演算子とハミルトニアンの交換子は以下のようになる：

\[[\mathbf{S}, \hat{H}] = i\hbar \gamma (\mathbf{S} \times \mathbf{B})\]

この結果を $\frac{d}{dt} \langle \mathbf{M} \rangle$ に代入する：

\[\frac{d}{dt} \langle \mathbf{M} \rangle = -\frac{i\gamma}{\hbar} \text{Tr}\left( \hat{\rho} \cdot i\hbar \gamma (\mathbf{S} \times \mathbf{B}) \right)\] \[\frac{d}{dt} \langle \mathbf{M} \rangle = \gamma^2 \text{Tr}\left( \hat{\rho} (\mathbf{S} \times \mathbf{B}) \right) = \gamma \left( \gamma \text{Tr}(\hat{\rho} \mathbf{S}) \right) \times \mathbf{B}\]

$\langle \mathbf{M} \rangle = \gamma \text{Tr}(\hat{\rho} \mathbf{S})$ であるから、最終的に以下の式が得られる：

\[\frac{d\langle \mathbf{M} \rangle}{dt} = \gamma \langle \mathbf{M} \rangle \times \mathbf{B}\]

証明完了（Q.E.D.）

この完全な証明により、微視的な量子密度行列の Liouville-von Neumann 方程式は、緩和（環境との相互作用）が存在しない閉じた量子系において、古典的トルク方程式と物理的・数学的に完全に一致することが示された。

4. 微視的 Bloch-Wangsness-Redfield (BWR) 緩和理論

第3節で示した Liouville-von Neumann 方程式 $\frac{d\hat{\rho}}{dt} = -\frac{i}{\hbar}[\hat{H}, \hat{\rho}]$ は、スピン系が完全に孤立している（可逆なユニタリ時間発展を行う）ことを前提としている。しかし現実の物質系では、スピン系は周囲の巨視的熱浴（結晶格子、溶媒の熱運動、電子雲など）と絶えず相互作用しており、この相互作用が熱力学的不可逆な緩和現象（$T_1, T_2$）を引き起こす。

本節では、R. K. Wangsness, F. Bloch, A. Redfield によって確立された Bloch-Wangsness-Redfield (BWR) 緩和理論 を構築する。システム・リザーバハミルトニアンから出発し、相互作用表示での2次摂動展開とマルコフ近似を経て、縦緩和率 $1/T_1$ および横緩和率 $1/T_2$ を局所微視的磁場ゆらぎのスペクトル密度関数 $J(\omega)$ の言葉で厳密に定式化する。

4.1 システム・リザーバハミルトニアンとゆらぐ局所磁場

全量子系は、観測対象であるスピン系（System: $S$）と、自由度の極めて大きい環境熱浴（Reservoir / Lattice: $R$）の複合系から構成される。全ハミルトニアン $\hat{H}$ は以下のように与えられる：

\[\hat{H} = \hat{H}_0 + \hat{H}_1(t) = \hat{H}_S + \hat{H}_R + \hat{H}_1(t)\]

各項の意味は以下の通りである：

  1. 静的スピンハミルトニアン $\hat{H}_S$: 静磁場 $\mathbf{B}_0 = (0, 0, B_0)^T$ によるゼーマン相互作用： \(\hat{H}_S = \hbar \omega_0 \hat{S}_z \quad (\omega_0 = -\gamma B_0)\)

  2. 熱浴ハミルトニアン $\hat{H}_R$: 格子振動（フォノン）、溶媒分子のランダムブラウン運動などの環境自由度を記述するハミルトニアン。

  3. 微視的相互作用ハミルトニアン $\hat{H}_1(t)$: 熱浴の熱運動によってスピン位置に引き起こされる時間変化・確率変動する局所微視的磁場 $\mathbf{B}(t) = (B_x(t), B_y(t), B_z(t))^T$ とスピンの結合： \(\hat{H}_1(t) = -\gamma \mathbf{B}(t) \cdot \mathbf{S} = -\gamma \left( B_x(t) \hat{S}_x + B_y(t) \hat{S}_y + B_z(t) \hat{S}_z \right)\) ただし、局所ゆらぎ磁場の熱浴平均はゼロであるとする（平均磁場は $\hat{H}_S$ に含まれている）： \(\langle \mathbf{B}(t) \rangle_R = \text{Tr}_R (\hat{\rho}_R \mathbf{B}(t)) = 0\)

4.2 相互作用表示と摂動展開によるマスター方程式

スピン系の静的運動（主たるラーモア歳差運動）と微小な相互作用による緩和動力を分離するため、相互作用表示（Interaction Picture） へユニタリ変換を行う。

主ハミルトニアン $\hat{H}0 = \hat{H}_S + \hat{H}_R$ を用いて、全密度演算子 $\tilde{\rho}{\text{total}}(t)$ および相互作用ハミルトニアン $\tilde{H}_1(t)$ の相互作用表示を定義する：

\[\tilde{\rho}_{\text{total}}(t) = e^{\frac{i}{\hbar} \hat{H}_0 t} \hat{\rho}_{\text{total}}(t) e^{-\frac{i}{\hbar} \hat{H}_0 t}\] \[\tilde{H}_1(t) = e^{\frac{i}{\hbar} \hat{H}_0 t} \hat{H}_1(t) e^{-\frac{i}{\hbar} \hat{H}_0 t}\]

相互作用表示における全系の Liouville-von Neumann 方程式は以下の通りである：

\[\frac{d\tilde{\rho}_{\text{total}}(t)}{dt} = -\frac{i}{\hbar} [\tilde{H}_1(t), \tilde{\rho}_{\text{total}}(t)]\]

この微分方程式を時刻 $0$ から $t$ まで両辺積分すると、形式的厳密解が得られる：

\[\tilde{\rho}_{\text{total}}(t) = \tilde{\rho}_{\text{total}}(0) - \frac{i}{\hbar} \int_{0}^{t} dt' [\tilde{H}_1(t'), \tilde{\rho}_{\text{total}}(t')]\]

これを元の微分方程式の右辺の $\tilde{\rho}_{\text{total}}(t)$ に代入することで、逐次展開（ボルン展開）を行う：

\[\frac{d\tilde{\rho}_{\text{total}}(t)}{dt} = -\frac{i}{\hbar} [\tilde{H}_1(t), \tilde{\rho}_{\text{total}}(0)] - \frac{1}{\hbar^2} \int_{0}^{t} dt' \left[ \tilde{H}_1(t), [\tilde{H}_1(t'), \tilde{\rho}_{\text{total}}(t')] \right]\]

ここで、以下の 2 つの物理的近似を導入する：

  1. ボルン近似（Born Approximation / 弱結合近似）: スピン系と熱浴の結合 $\hat{H}1$ は十分弱く、熱浴の熱容量が巨大であるため、全系の密度演算子は常にスピン系の被約密度行列 $\tilde{\rho}(t) = \text{Tr}_R (\tilde{\rho}{\text{total}}(t))$ と熱浴の熱平衡密度行列 $\hat{\rho}_R = e^{-\hat{H}_R / k_B T} / \text{Tr}_R(e^{-\hat{H}_R / k_B T})$ の直積で近似できる： \(\tilde{\rho}_{\text{total}}(t) \approx \tilde{\rho}(t) \otimes \hat{\rho}_R\)

  2. 熱浴平均の消滅: 両辺の熱浴自由度についてトレース $\text{Tr}_R$ をとる。局所ゆらぎ磁場の平均がゼロ $\text{Tr}_R(\tilde{H}_1(t) \hat{\rho}_R) = 0$ であるため、展開第1項 $-\frac{i}{\hbar} \text{Tr}_R([\tilde{H}_1(t), \tilde{\rho}(0) \otimes \hat{\rho}_R])$ は完全に消失する。

したがって、スピン系の被約密度行列 $\tilde{\rho}(t)$ に関する 2次摂動マスター方程式（Second-Order Quantum Master Equation） が得られる：

\[\frac{d\tilde{\rho}(t)}{dt} = -\frac{1}{\hbar^2} \int_{0}^{t} dt' \, \text{Tr}_R \left( \left[ \tilde{H}_1(t), \left[ \tilde{H}_1(t'), \tilde{\rho}(t') \otimes \hat{\rho}_R \right] \right] \right)\]

4.3 マルコフ近似と Redfield 緩和マトリックス

得られた 2 次摂動マスター方程式は、右辺の積分に過去の履歴 $\tilde{\rho}(t’)$ が含まれているため、非マルコフ的（メモリー効果を持つ）微分積分方程式である。

マルコフ近似（Markov Approximation）の導入

物質系（液体や固体）における熱浴の揺らぎの記憶時間（相関時間 $\tau_c$） は極めて短く（通常 $\tau_c \sim 10^{-12} \text{ s} = 1 \text{ ps}$）、スピン系が緩和を起こす巨視的時間スケール（$T_1, T_2 \sim 10^{-3} \text{ – } 1 \text{ s}$）に比べて圧倒的に短い：

\[\tau_c \ll T_1, T_2\]

この時間スケールの強烈な分離に基づき、以下の 2 つの操作を行う：

  1. メモリー効果の消滅（Redfield 近似）: 時刻 $t’$ におけるスピンの状態 $\tilde{\rho}(t’)$ を、現在の時刻 $t$ における状態 $\tilde{\rho}(t)$ で置き換える（$\tilde{\rho}(t’) \to \tilde{\rho}(t)$）。
  2. 積分時間軸の無次元化と無限大拡張: 時間差 $\tau = t - t’$ を導入する。熱浴相関は $\tau > \tau_c$ で急激に $0$ に減衰するため、観測時間 $t \gg \tau_c$ において積分の積分上限を $t \to \infty$ へ拡張する。

これらのマルコフ近似を適用すると、マスター方程式は完全に現在の状態 $\tilde{\rho}(t)$ のみに依存するマルコフ的運動方程式に変換される：

\[\frac{d\tilde{\rho}(t)}{dt} = -\frac{1}{\hbar^2} \int_{0}^{\infty} d\tau \, \text{Tr}_R \left( \left[ \tilde{H}_1(t), \left[ \tilde{H}_1(t-\tau), \tilde{\rho}(t) \otimes \hat{\rho}_R \right] \right] \right)\]

Redfield 緩和マトリックス（Redfield Relaxation Matrix）

基底状態 ${ \vert a \rangle, \vert b \rangle, \vert c \rangle, \vert d \rangle }$ （エネルギー固有値 $E_a = \hbar \omega_a$）での行列要素を展開すると、Redfield 方程式（Redfield Equation） が得られる：

\[\frac{d\tilde{\rho}_{ab}(t)}{dt} = \sum_{c, d} R_{abcd} \, \tilde{\rho}_{cd}(t) \, e^{i(\omega_{ab} - \omega_{cd})t}\]

ここで $\omega_{ab} = \omega_a - \omega_b$ であり、$R_{abcd}$ は Redfield 緩和テンソル（Redfield Relaxation Superoperator Elements） と呼ばれる定数行列要素である。 さらに、$\vert \omega_{ab} - \omega_{cd} \vert \gg 1/T_{1,2}$ となる快速振動項を平均化して消去する近似（Secular Approximation / 永年近似）を行うことで、ポピュレーション間および異なるコヒーレンス間のデカップリングが達成される。

4.4 相関関数とスペクトル密度関数 $J(\omega)$ による $1/T_1, 1/T_2$ の微視的表現

ここから、微視的な局所磁場ゆらぎ $\mathbf{B}(t)$ の統計的性質から、巨視的緩和率 $1/T_1$ および $1/T_2$ の具体的な解析的表示を明示的に導出する。

局所磁場ゆらぎの時間相関関数とスペクトル密度関数

局所微視的磁場 $B_k(t) \, (k = x, y, z)$ の時間相関を表す 自己相関関数（Autocorrelation Function） $g_k(\tau)$ を以下のように定義する：

\[g_k(\tau) = \langle B_k(t) B_k(t-\tau) \rangle_R = \text{Tr}_R \left( \hat{\rho}_R B_k(t) B_k(t-\tau) \right)\]

熱浴の異方性がなく、異なる成分間のゆらぎが無相関（$\langle B_j(t) B_k(t’) \rangle_R = \delta_{jk} g_k(t-t’)$）であると仮定する。

ウィーナー・ヒンチンの定理（Wiener-Khinchin Theorem）に従い、時間相関関数 $g_k(\tau)$ のフーリエ変換として スペクトル密度関数（Spectral Density Function） $J_k(\omega)$ を次のように定義する：

\[J_k(\omega) = \int_{-\infty}^{\infty} \langle B_k(t) B_k(0) \rangle_R \, e^{-i\omega t} dt = \int_{-\infty}^{\infty} g_k(\tau) \, e^{-i\omega \tau} d\tau\]

【ローレンツ型スペクトル密度の典型例】 相関関数が特徴的相関時間 $\tau_c$ で指数関数的に減衰する標準的なモデル $g_k(\tau) = \langle B_k^2 \rangle e^{-\vert \tau \vert / \tau_c}$ を採用すると、スペクトル密度関数 $J_k(\omega)$ は以下のローレンツ型（Lorentzian form）となる： \(J_k(\omega) = \langle B_k^2 \rangle \frac{2 \tau_c}{1 + \omega^2 \tau_c^2}\)

$S=1/2$ スピン系における $1/T_1$ および $1/T_2$ の明示的導出

相互作用ハミルトニアン $\hat{H}1(t) = -\gamma \mathbf{B}(t) \cdot \mathbf{S}$ を、スピン昇降演算子 $\hat{S}\pm = \hat{S}_x \pm i \hat{S}_y$ を用いて書き直す：

\[\hat{H}_1(t) = -\gamma \left[ \frac{1}{2} B_+(t) \hat{S}_- + \frac{1}{2} B_-(t) \hat{S}_+ + B_z(t) \hat{S}_z \right]\]

ただし $B_\pm(t) = B_x(t) \pm i B_y(t)$ である。

相互作用表示への変換 $\tilde{S}\pm(t) = e^{i \omega_0 \hat{S}_z t} \hat{S}\pm e^{-i \omega_0 \hat{S}z t} = \hat{S}\pm e^{\pm i \omega_0 t}$ および $\tilde{S}_z(t) = \hat{S}_z$ を考慮し、2次摂動マスター方程式を実行する。

  1. 縦緩和率 $1/T_1$ （対角成分 $\rho_{aa}$ / 縦磁化 $\langle S_z \rangle$ の運動）: 縦磁化の減衰速度方程式 $\frac{d\langle S_z \rangle}{dt} = -\frac{1}{T_1} (\langle S_z \rangle - S_0)$ を計算すると、スピン反転演算子 $\hat{S}+, \hat{S}-$ に結合する横方向の磁場ゆらぎ $B_x(t), B_y(t)$ のみが寄与し、スピン固有のラーモア振動数 $\omega_0$ におけるスペクトル密度成分によって以下のように与えられる：
\[\frac{1}{T_1} = \gamma^2 \left( J_x(\omega_0) + J_y(\omega_0) \right)\]

熱浴が $xy$ 平面内で等方的な場合（$J_x = J_y = J$）、縦緩和率は次のように表現される：

\[\frac{1}{T_1} = 2 \gamma^2 J(\omega_0)\]
  1. 横緩和率 $1/T_2$ （非対角成分 $\rho_{ab}$ / 横磁化 $\langle S_+ \rangle$ の運動）: 横磁化（コヒーレンス）の減衰速度方程式 $\frac{d\langle S_+ \rangle}{dt} = -\frac{1}{T_2} \langle S_+ \rangle$ を計算すると、スピン反転に伴うエネルギー遷移の寄与（$\frac{1}{2T_1}$）に加えて、$z$ 軸方向の局所磁場ゆらぎ $B_z(t)$ による極限周波数 $\omega = 0$ の純粋位相散逸（Pure Dephasing）が寄与する：
\[\frac{1}{T_2} = \frac{1}{2 T_1} + \gamma^2 J_z(0)\]

$J_x(\omega_0), J_y(\omega_0)$ を用いて明示的に展開すると、以下の通りである：

\[\frac{1}{T_2} = \frac{\gamma^2}{2} \left( J_x(\omega_0) + J_y(\omega_0) \right) + \gamma^2 J_z(0)\]

微視的 BWR 緩和理論が解き明かす物理的描像

導出された微視的表現は、第1節で現象論的に示された緩和時間の本質を見事に説明している：

  • 縦緩和 $T_1$ の微視的描像: スピンが $z$ 軸方向の占有数を変える（エネルギー順位間を遷移する）ためには、ゼーマンエネルギーの差 $\hbar \omega_0$ に等しいエネルギー量子を熱浴とやり取りする必要がある。そのため、ラーモア振動数 $\omega = \omega_0$ に一致する局所磁場の熱ゆらぎ成分 $J(\omega_0)$ のみが $T_1$ 緩和を引き起こす。

  • 横緩和 $T_2$ の微視的描像: 横コヒーレンスの失調には、エネルギー遷移を伴うスピン反転プロセス（$1/(2T_1)$）だけでなく、エネルギーを変化させずに局所ラーモア振動数を場所ごとに変調させる $z$ 方向の静的・低周波磁場ゆらぎ $J_z(0)$ （Pure Dephasing）が加算される。 これにより、常に以下の基本不等式が自発的に導かれる： \(\frac{1}{T_2} \ge \frac{1}{2 T_1} \quad \implies \quad T_2 \le 2 T_1\)

5. 高次スピン（$S \ge 1$）におけるゼロ磁場分裂（ZFS）テンソルと Bloch 方程式の破綻

前節までの論議では、スピン $S = 1/2$ の系を中心として現象論的 Bloch 方程式の挙動から量子密度行列および BWR 緩和理論への展開を追ってきた。スピン $S = 1/2$ の系では、状態空間が 2 次元（$\vert \uparrow \rangle, \vert \downarrow \rangle$）であり、密度行列は単位行列 $\hat{I}$ と 3 つのパウリ行列 $\boldsymbol{\sigma} = (\hat{\sigma}_x, \hat{\sigma}_y, \hat{\sigma}_z)^T$ の線形結合として完全に展開できる。このため、スピン演算子の期待値 $\langle \mathbf{S} \rangle$ （3次元の磁化ベクトル $\mathbf{M}$）のみを追跡することで、系のすべての動力学と緩和を不足なく記述することができた。

しかし、電子スピン $S \ge 1$ （あるいは四重極モーメントを持つ核スピン $I \ge 1$）の系においては、静磁場が非存在（$\mathbf{B}_0 = 0$）の状態であっても、電子間のスピン–スピン双極子相互作用や異方的なスピン–軌道相互作用に由来する ゼロ磁場分裂（Zero-Field Splitting: ZFS） が発生する。これにより、等方的な $S = 1/2$ を前提とした現象論的 3次元 Bloch 方程式は根本的に破綻する。

本節では、二次スピン異方性ハミルトニアンの数理、ZFS による Bloch 方程式の破綻メカニズム、そして NV センターや単一分子磁性体（SMM）、ESR/EPR における物理的影響について解説する。

5.1 二次スピン異方性ハミルトニアンと ZFS パラメータ $D, E$

スピン $S \ge 1$ を持つ電子スピン系において、スピン空間の結晶場や微視的電子構造による異方性を記述する二次スピンハミルトニアン $\hat{H}_{\text{ZFS}}$ は、対称 2 階テンソル $\mathbf{D}$ （ZFS テンソル）を用いて一般に以下のように定義される：

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

ここで ZFS テンソル $\mathbf{D}$ は $3 \times 3$ の実対称行列（$D_{ij} = D_{ji}$）であり、スピン空間のトレースをハミルトニアンの定数シフトとして取り除くことでトレースレス（$\text{Tr}(\mathbf{D}) = D_{xx} + D_{yy} + D_{zz} = 0$）に選ぶことができる。

テンソル $\mathbf{D}$ の主軸系（Principal Axis System: PAS）に直交座標系をとると、非対角成分は消去され、ハミルトニアンは主軸成分 $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\]

ここで、トレースレス条件 $D_{xx} + D_{yy} + D_{zz} = 0$ を用いて、独立な 2 つのスピン異方性パラメータである 軸対称 ZFS パラメータ $D$ および 斜方（菱形）ZFS パラメータ $E$ を以下のように定義する：

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

この定義式から逆算すると、主軸成分は $D$ と $E$ を用いて以下のように表される：

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

これらを対角化ハミルトニアンに代入して整理する：

\[\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 \\ &= D \left( \hat{S}_z^2 - \frac{1}{3} (\hat{S}_x^2 + \hat{S}_y^2 + \hat{S}_z^2) \right) + E (\hat{S}_x^2 - \hat{S}_y^2) \end{aligned}\]

スピン角運動量の自乗演算子の性質 $\hat{S}_x^2 + \hat{S}_y^2 + \hat{S}_z^2 = \mathbf{S}^2 = S(S+1) \hat{I}$ を使用すると、標準的な ZFS スピンハミルトニアン の形式が得られる：

\[\hat{H}_{\text{ZFS}} = D \left( \hat{S}_z^2 - \frac{1}{3} S(S+1) \hat{I} \right) + E (\hat{S}_x^2 - \hat{S}_y^2)\]
  • $D$ 項（軸対称異方性）: $z$ 軸方向を特別視し、$m_s = 0$ 状態と $m_s = \pm 1$ 状態のエネルギーを分裂させる。
  • $E$ 項（斜方異方性）: $xy$ 平面内の非軸対称な歪みを表し、$m_s = +1$ と $m_s = -1$ の状態間を量子力学的に混合（クーリング）させる。

5.2 ゼロ磁場準位分裂による等方現象論的 Bloch 方程式の破綻

スピン $S \ge 1$ において ZFS テンソルが存在する場合、第1節で導入した 3次元現象論的 Bloch 方程式が成立しなくなる物理的・数学的理由は以下の 3 点に集約される：

1. 状態空間の次元拡張と多重極子（Quadrupole）コヒーレンス

スピン $S$ の Hilbert 空間の次元は $N = 2S + 1$ であり、密度行列 $\hat{\rho}$ の独立な実パラメータ数は $N^2 - 1$ 個存在する。

  • $S = 1/2$ ($N=2$): 独立パラメータは $2^2 - 1 = 3$ 個であり、これはまさに 3次元ベクトル $\mathbf{M} = (M_x, M_y, M_z)^T$ と 1 対 1 に対応する。
  • $S = 1$ ($N=3$): 独立パラメータは $3^2 - 1 = 8$ 個存在する。 3次元磁化ベクトル $\mathbf{M} \propto \langle \mathbf{S} \rangle$ は、スピン双極子演算子（$\hat{S}x, \hat{S}_y, \hat{S}_z$）の期待値である 3 つの成分しか追跡できない。残りの 5 つの成分は、四重極演算子（$\hat{S}_z^2$, $\hat{S}_x^2 - \hat{S}_y^2$, $\hat{S}_x \hat{S}_y + \hat{S}_y \hat{S}_x$ など）に対応する非線形スピン相関や双子状態間のコヒーレンス（$\Delta m_s = \pm 2$ 遷移など）である。 ZFS ハミルトニアン $\hat{H}{\text{ZFS}}$ に含まれる演算子は $\hat{S}_i \hat{S}_j$ の二次形式であるため、時間発展において双極子演算子 $\langle \mathbf{S} \rangle$ と四重極演算子 $\langle \hat{S}_i \hat{S}_j \rangle$ が強く結合し、3次元ベクトル $\mathbf{M}$ のみで閉じた常微分方程式系を構成することが不可能となる。

2. ゼロ磁場におけるエネルギーレベルの固有分裂

外力としての静磁場が存在しない条件（$\mathbf{B}0 = 0$）であっても、例えば $S=1, E=0$ の系では $\hat{H}{\text{ZFS}} = D (\hat{S}_z^2 - \frac{2}{3} \hat{I})$ により、$\vert m_s = 0 \rangle$ 状態（エネルギー $-2D/3$）と 2重縮退した $\vert m_s = \pm 1 \rangle$ 状態（エネルギー $+D/3$）の間に固有のエネルギー差 $\Delta E = D$ が生じる。 古典 Bloch 方程式では $\mathbf{B}_0 = 0$ のときトルク項 $\gamma \mathbf{M} \times \mathbf{B}$ が消滅し、単なる指数関数的衰退しか予想しないが、ZFS が存在する高次スピン系では磁場ゼロであっても角振動数 $\omega = D/\hbar$ の固有歳差運動や量子ビート（Quantum Beats）が観測される。

3. 異方的な非単一指数関数的緩和

ZFS を持つスピン系が熱浴のゆらぎを受けると、Redfield 緩和テンソル $R_{abcd}$ の要素は遷移する準位の対によって大きく異なる値をとる。その結果、縦磁化や横磁化の緩和曲線は単一の指数関数 $e^{-t/T_1}$ や $e^{-t/T_2}$ ではなく、複数の異なる減衰定数を持つ多重指数関数（Multiexponential Relaxation）となる。

5.3 NV センター・単一分子磁性体（SMM）・EPR における物理的影響

ZFS テンソルとスピン異方性は、現代の量子情報科学や電子スピン物性において極めて重要かつ中心的な役割を果たしている。

1. ダイヤモンド NV センター（$\text{NV}^-$ Center in Diamond）

ダイヤモンド中の窒素–空孔色中心（NVセンター）の基底状態は、電子スピン $S=1$ のトリプレット状態（三重項）である。

  • 結晶軸に沿った C3v 対称性により、強固な軸対称 ZFS パラメータ $D \approx 2.87 \text{ GHz}$ を有する。
  • ゼロ磁場で $\vert m_s = 0 \rangle$ と $\vert m_s = \pm 1 \rangle$ が微視的に大きく分裂しているため、光ポンピング（532 nm 緑色レーザー照射）によって基底状態 $\vert m_s = 0 \rangle$ への高効率なスピン偏極（Initial State Preparation）が可能となる。
  • 微小な歪みや局所電場による斜方パラメーター $E$ の存在は、$\vert m_s = \pm 1 \rangle$ の縮退を解いて基底状態を $\frac{1}{\sqrt{2}} (\vert +1 \rangle \pm \vert -1 \rangle)$ へと重畳させる。
  • 格子振動（フォノン）による ZFS パラメータの熱的ゆらぎ $D(t), E(t)$ は、室温における NV センターの純粋位相散逸（Pure Dephasing）および $T_1$ 緩和プロセス（Orbach 過程や Raman 過程）の主要な経路を開く。

2. 単一分子磁性体（Single-Molecule Magnets: SMMs）

$\text{Fe}8$ や $\text{Mn}{12}$ に代表される単一分子磁性体では、大スピン（$S = 10$ など）と負の大きな軸対称異方性（$D < 0$）が結合することにより、$m_s = \pm S$ 状態を最安定とする巨視的なエネルギー障壁 $U = \vert D \vert S^2$ が形成される。

  • この異方性障壁により、分子単体で低温において永久磁石のような磁化の熱的緩和の遅延（Slow Relaxation of Magnetization）を示す。
  • 斜方項 $E (\hat{S}_x^2 - \hat{S}_y^2)$ や横磁場は、エネルギー障壁を挟んだ $+m_s$ 状態と $-m_s$ 状態との間に量子トンネル遷移（Quantum Tunneling of Magnetization: QTM）を引き起こす。これにより、古典的な Bloch 型の熱励起緩和とは全く異なる、共鳴的トンネリングによる階段状の非平衡磁化緩和減衰が実現する。

3. 電子スピン共鳴（EPR / ESR）スペクトロスコピー

常磁性金属錯体やラジカル、固体内欠陥の EPR 観測において、ZFS テンソルはスペクトルの分裂や線形の決定因子となる。

  • 多結晶（粉末）試料や生体高分子試料では、ZFS テンソルの主軸に対する静磁場 $\mathbf{B}_0$ の配向角がランダムに分布するため、巨大な配向異方性線形（Powder Pattern）が形成される。
  • BWR 緩和理論を適用すると、配向角ごとに実効的な $1/T_1$ および $1/T_2$ が大幅に異なるため、スペクトルの部位ごとに異なるパワー飽和特性が観測される。

まとめと展望：古典 Bloch 緩和から量子密度行列、そして ZFS テンソル解析へ

本稿では、スピン物理学における緩和現象の階層的な数理体系を、現象論的アプローチから微視的量子統計アプローチに至るまで厳密に解説してきた。

  1. 現象論的 Bloch 方程式: 古典的トルク方程式 $\frac{d\mathbf{M}}{dt} = \gamma \mathbf{M} \times \mathbf{B}$ に現象論的パラメータ $T_1, T_2$ を導入し、回転座標系（Rotating Frame）へのユニタリ座標変換によって RF 照射下の定常解（吸収・分散ローレンツ線形、飽和因子 $S_{\text{sat}}$、パワー広がり）を厳密に導出した。

  2. 密度行列と Liouville-von Neumann 方程式: アンサンブルの統計的状態を表す密度演算子 $\hat{\rho}$ を導入し、その基礎方程式 $\frac{d\hat{\rho}}{dt} = -\frac{i}{\hbar}[\hat{H}, \hat{\rho}]$ から、孤立系における期待値 $\langle \mathbf{M} \rangle$ の運動が古典トルク方程式を厳密に再現することを数学的に証明した。

  3. Bloch-Wangsness-Redfield (BWR) 微視的緩和理論: スピン系と熱浴の弱結合近似（Born 近似）および時間相関の短時間減衰（Markov 近似）のもとで2次摂動マスター方程式を展開し、縦緩和率 $1/T_1 = 2\gamma^2 J(\omega_0)$ および横緩和率 $1/T_2 = \frac{1}{2T_1} + \gamma^2 J_z(0)$ を局所磁場ゆらぎのスペクトル密度関数 $J(\omega)$ の言葉で微視的に関係づけた。

  4. 高次スピンへの拡張と ZFS テンソルの架け橋: スピン $S \ge 1$ における二次スピン異方性ハミルトニアン $\hat{H}_{\text{ZFS}} = D (\hat{S}_z^2 - \frac{1}{3}S(S+1)\hat{I}) + E (\hat{S}_x^2 - \hat{S}_y^2)$ を提示し、四重極コヒーレンスの存在とゼロ磁場準位分裂によって 3次元 Bloch 方程式が破綻する物理的理由を解明した。

次回予告：ゼロ磁場分裂（ZFS）テンソルの幾何学と数値シミュレーション

本稿で架け橋を渡した ZFS テンソル $\mathbf{D}$ は、単にスピンハミルトニアンの対角成分にとどまらず、スピン空間におけるオイラー角回転変換（$D_{ij} \to \mathbf{R} \mathbf{D} \mathbf{R}^T$）、Spherical Tensor Operator（球面テンソル演算子 $T_{2m}$）による代数的展開、そして EasySpin 等の量子スピンシミュレータにおける計算基盤となっている。

次回の記事では、ZFS テンソル $\mathbf{D}$ のテンソル表現、3次元配向異方性の幾何学的可視化、オイラー角による座標変換、および高次スピン系のスピンハミルトニアン対角化とスペクトルシミュレーションの完全な数理を徹底的に掘り下げて解き明かしていく。

You May Also Enjoy