直積演算子
核磁気共鳴において直積演算子(プロダクト演算子、プロダクトオペレーター)法とは、複数のスピン系の状態を表す密度演算子\hat{\rho} \ を直積演算子で表すことで、スピン系の時間発展を記述する方法である。
一般的な量子統計力学による取り扱いでは密度行列を用いて状態を表さなければならないが、これは数個のスピン系でさえも大変複雑で、物理的イメージも分かりにくい。たとえば最も単純な2スピン系であっても、密度行列の行列要素は合計16個になる。そこでSørensenらによって考案されたのが直積演算子法である。
スペクトルの強度を問題にしないならば、行列要素の値そのものは問題ではなく、どのように時間発展するか分かればよい。
一般に任意のエルミート演算子は、以下のように規格直交化された演算子の完全系\{ \hat{F}_i \} で展開できる。
Tr[\hat{F}_i \hat{F}_j ] = c\delta_{ij}
よって密度演算子も、これらを用いて展開できる。
\hat{\rho} = \sum_i \rho_i \hat{F}_i
そこで基底演算子\{ \hat{F}_i \} として、各スピン角運動量演算子の直積(直積演算子)を採用したのが直積演算子法である。
直積演算子法を適用するには次の3つの前提がある。
1スピン系では21×21=4個の密度行列から成る。
E/2 = 1/2 \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}:単位演算子
I_x = 1/2 \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}:Iスピンの磁化のx成分
I_y = i/2 \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix}:Iスピンの磁化のy成分
I_z = 1/2 \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}:Iスピンの磁化のz成分
弱く結合した2スピン系(スピンIとスピンS)では直積演算子は22×22=16個ある。ここでは直積A\times Bを単にABと記す。
E/2 = 1/2\begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{pmatrix}
I_x = I_x \times E_S = 1/2 \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} \times \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} = 1/2\begin{pmatrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \end{pmatrix}
I_y = i/2\begin{pmatrix} 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & -1 \\ 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \end{pmatrix}
I_z = 1/2\begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & -1 \end{pmatrix}
S_x = 1/2\begin{pmatrix} 0 & 1 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 1 & 0 \end{pmatrix}
S_y = i/2\begin{pmatrix} 0 & -1 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & -1 \\ 0 & 0 & 1 & 0 \end{pmatrix}
S_z = 1/2\begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & -1 \end{pmatrix}
2I_xS_z = 2 \times 1/2 \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} \times 1/2\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} = 1/2\begin{pmatrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & -1 \\ 1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \end{pmatrix}
反位相のIスピンの磁化(Iスピンの磁化のx成分がSスピンの縦方向の2つの可能な状態に対応して2つの反位相成分に分裂することを表し、ベクトルモデルで表すことは可能である。)
2I_yS_z = i/2\begin{pmatrix} 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & 1 \\ 1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \end{pmatrix}:反位相のIスピンの磁化
2I_zS_x = 1/2\begin{pmatrix} 0 & 1 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & -1 \\ 0 & 0 & -1 & 0 \end{pmatrix}:反位相のSスピンの磁化
2I_zS_y = i/2\begin{pmatrix} 0 & -1 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & -1 & 0 \end{pmatrix}:反位相のSスピンの磁化
2I_zS_z = 1/2\begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & 1 \end{pmatrix}:縦2-スピンオーダー
2I_xS_x = 1/2\begin{pmatrix} 0 & 0 & 0 & 1 \\ 0 & 0 & 1 & 0 \\ 0 & 1 & 0 & 0 \\ 1 & 0 & 0 & 0 \end{pmatrix}:2-スピンコヒーレンス
2I_yS_y = 1/2\begin{pmatrix} 0 & 0 & 0 & -1 \\ 0 & 0 & 1 & 0 \\ 0 & 1 & 0 & 0 \\ -1 & 0 & 0 & 0 \end{pmatrix}:2-スピンコヒーレンス
2I_xS_y = i/2\begin{pmatrix} 0 & 0 & 0 & -1 \\ 0 & 0 & 1 & 0 \\ 0 & -1 & 0 & 0 \\ 1 & 0 & 0 & 0 \end{pmatrix}:2-スピンコヒーレンス
2I_yS_x = i/2\begin{pmatrix} 0 & 0 & 0 & -1 \\ 0 & 0 & -1 & 0 \\ 0 & 1 & 0 & 0 \\ 1 & 0 & 0 & 0 \end{pmatrix}:2-スピンコヒーレンス
縦2-スピンオーダーと2-スピンコヒーレンスはベクトルモデルで表現することができない量子的現象である。
液体のパルスNMRでは、次の3つの時間変化が重要である。
スピン結合(Jカップリング)に由来する磁化の回転
あるパルス系列の働きは、その中のパルスおよびパルスの間に起こる自由歳差運動の時間の順序に依存する。しかし、もし弱く結合したスピン系を考えているならば、この時間の間に化学シフトとスピン結合による時間展開の順番を自由に変えることができる(交換できる)。
時間変化は、密度演算子の時間発展を記述する式(フォン・ノイマン方程式)を解くことにより導ける。以下では、その結果だけを示す。
ラジオ波パルスを与えている間、化学シフトとスピン結合による回転は無いと考えて良い。ラジオ波パルスによる回転は、x-y平面上の軸(通常はx軸)のまわりに起きる。必要なら、傾いた有効磁場のまわりの回転を考えることができる。
まず第一に、I_zで表されるIスピンの磁化のz成分に励起パルス\beta(x)が働いた場合、以下のように時間変化する。
I_z \xrightarrow{\quad \beta I_x \quad} I_z cos \beta - I_y sin \beta
例として90度パルスの場合は、以下のようになる。
I_z \xrightarrow{\quad (\pi/2) I_x \quad} -I_y
S_z \xrightarrow{\quad (\pi/2) S_x \quad} -S_y
また13Cや15Nの感度向上のために異種核ISにおける分極移動を用いるINEPTでは、Iスピンの磁化を回転座標系の±x軸にそって範囲層にそろえて、Iスピンに+y軸に沿って90度パルスを与える。この回転は次のように書くことができる。
2I_xS_z \xrightarrow{\quad (\pi/2) I_y \quad} -2I_zS_z
これによって縦2-スピンオーダーが作り出される。
化学シフトによるスピンの時間変化はつぎのように表すことができる。
I_y \xrightarrow{\quad (\pi \delta_I t) I_z \quad} I_ycos(2\pi \delta_I t)-I_xsin
(2\pi \delta_I t)
ここで\delta_Iは、トランスミッター周波数から測ったIスピンのシフトである。Sスピンの化学シフトも同様に扱うことができる。
異種核スピンの場合は、個々のスピンに別の回転系が設定されるから、それぞれの化学シフトはそれぞれの回転系におけるトランスミッター周波数からの相対値として表現する。
直積演算子法では、スピン結合はつぎのように表される。
I_y \xrightarrow{\quad (2\pi J_{IS} \tau)2I_zS_z \quad} I_ycos(2\pi J_{IS} \tau) -2I_xS_zsin(2\pi J_{IS} \tau)
もしも\tau=1/(2J_{IS})となるように選ぶと、
I_y \xrightarrow{\quad (\pi/2)2I_zS_z \quad} -2I_xS_z
さらに次の自由歳差運動の期間を考えると、
-2I_xS_z \xrightarrow{\quad (\pi J_{IS} \tau) 2I_zS_z \quad}
-2I_xS_zcos(\pi J_{IS} \tau)-I_ysin(\pi J_{IS} \tau)
となる。もしこの時間を\tau=1/(2J_{IS})となるように設定すると、2つのベクトルは再び-y軸にそって配向する。
-2I_xS_z \xrightarrow{\quad (\pi/2)2I_zS_z \quad}
-I_y
== 脚注 ==