一般に任意のエルミート演算子 は、以下のように(規格直交化 )された演算子 の完全系 { F ^ i } {\displaystyle \{{\hat {F}}_{i}\}} で展開できる。
T r [ F ^ i F ^ j ] = c δ i j {\displaystyle Tr[{\hat {F}}_{i}{\hat {F}}_{j}]=c\delta _{ij}} よって密度演算子も、これらを用いて展開できる。
ρ ^ = ∑ i ρ i F ^ i {\displaystyle {\hat {\rho }}=\sum _{i}\rho _{i}{\hat {F}}_{i}} そこで基底 演算子 { F ^ i } {\displaystyle \{{\hat {F}}_{i}\}} として、各スピン角運動量演算子 の直積 (直積演算子 )を採用したのが直積演算子法 である。
直積演算子法を適用するには次の3つの前提がある。
弱く結合したスピン系を扱う (高温近似 )が成り立つ 緩和 の影響を無視する1スピン系の例 1スピン系では21 ×21 =4個の密度行列から成る。
E / 2 = 1 / 2 ( 1 0 0 1 ) {\displaystyle E/2=1/2{\begin{pmatrix}1&0\\0&1\end{pmatrix}}} :単位演算子 I x = 1 / 2 ( 0 1 1 0 ) {\displaystyle I_{x}=1/2{\begin{pmatrix}0&1\\1&0\end{pmatrix}}} : I {\displaystyle I} スピンの磁化のx成分 I y = i / 2 ( 0 1 − 1 0 ) {\displaystyle I_{y}=i/2{\begin{pmatrix}0&1\\-1&0\end{pmatrix}}} : I {\displaystyle I} スピンの磁化のy成分 I z = 1 / 2 ( 1 0 0 − 1 ) {\displaystyle I_{z}=1/2{\begin{pmatrix}1&0\\0&-1\end{pmatrix}}} : I {\displaystyle I} スピンの磁化のz成分弱く結合した2スピン系の例 弱く結合した2スピン系(スピンIとスピンS)では直積演算子は22 ×22 =16個ある。ここでは直積 A × B {\displaystyle A\times B} を単に A B {\displaystyle AB} と記す。
E / 2 = 1 / 2 ( 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ) {\displaystyle 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 × E S = 1 / 2 ( 0 1 1 0 ) × ( 1 0 0 1 ) = 1 / 2 ( 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0 ) {\displaystyle 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 ( 0 0 − 1 0 0 0 0 − 1 1 0 0 0 0 1 0 0 ) {\displaystyle 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 ( 1 0 0 0 0 1 0 0 0 0 − 1 0 0 0 0 − 1 ) {\displaystyle 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 ( 0 1 0 0 1 0 0 0 0 0 0 1 0 0 1 0 ) {\displaystyle 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 ( 0 − 1 0 0 1 0 0 0 0 0 0 − 1 0 0 1 0 ) {\displaystyle 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 ( 1 0 0 0 0 − 1 0 0 0 0 1 0 0 0 0 − 1 ) {\displaystyle S_{z}=1/2{\begin{pmatrix}1&0&0&0\\0&-1&0&0\\0&0&1&0\\0&0&0&-1\end{pmatrix}}} 2 I x S z = 2 × 1 / 2 ( 0 1 1 0 ) × 1 / 2 ( 1 0 0 − 1 ) = 1 / 2 ( 0 0 1 0 0 0 0 − 1 1 0 0 0 0 − 1 0 0 ) {\displaystyle 2I_{x}S_{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 {\displaystyle I} スピンの磁化( I {\displaystyle I} スピンの磁化のx成分が S {\displaystyle S} スピンの縦方向の2つの可能な状態に対応して2つの反位相成分に分裂することを表し、ベクトルモデル で表すことは可能である。) 2 I y S z = i / 2 ( 0 0 − 1 0 0 0 0 1 1 0 0 0 0 − 1 0 0 ) {\displaystyle 2I_{y}S_{z}=i/2{\begin{pmatrix}0&0&-1&0\\0&0&0&1\\1&0&0&0\\0&-1&0&0\end{pmatrix}}} :反位相の I {\displaystyle I} スピンの磁化 2 I z S x = 1 / 2 ( 0 1 0 0 1 0 0 0 0 0 0 − 1 0 0 − 1 0 ) {\displaystyle 2I_{z}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 {\displaystyle S} スピンの磁化 2 I z S y = i / 2 ( 0 − 1 0 0 1 0 0 0 0 0 0 1 0 0 − 1 0 ) {\displaystyle 2I_{z}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 {\displaystyle S} スピンの磁化 2 I z S z = 1 / 2 ( 1 0 0 0 0 − 1 0 0 0 0 − 1 0 0 0 0 1 ) {\displaystyle 2I_{z}S_{z}=1/2{\begin{pmatrix}1&0&0&0\\0&-1&0&0\\0&0&-1&0\\0&0&0&1\end{pmatrix}}} :縦2-スピンオーダー 2 I x S x = 1 / 2 ( 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 0 ) {\displaystyle 2I_{x}S_{x}=1/2{\begin{pmatrix}0&0&0&1\\0&0&1&0\\0&1&0&0\\1&0&0&0\end{pmatrix}}} :2-スピンコヒーレンス 2 I y S y = 1 / 2 ( 0 0 0 − 1 0 0 1 0 0 1 0 0 − 1 0 0 0 ) {\displaystyle 2I_{y}S_{y}=1/2{\begin{pmatrix}0&0&0&-1\\0&0&1&0\\0&1&0&0\\-1&0&0&0\end{pmatrix}}} :2-スピンコヒーレンス 2 I x S y = i / 2 ( 0 0 0 − 1 0 0 1 0 0 − 1 0 0 1 0 0 0 ) {\displaystyle 2I_{x}S_{y}=i/2{\begin{pmatrix}0&0&0&-1\\0&0&1&0\\0&-1&0&0\\1&0&0&0\end{pmatrix}}} :2-スピンコヒーレンス 2 I y S x = i / 2 ( 0 0 0 − 1 0 0 − 1 0 0 1 0 0 1 0 0 0 ) {\displaystyle 2I_{y}S_{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 {\displaystyle I_{z}} で表される I {\displaystyle I} スピンの磁化のz成分に励起パルス β ( x ) {\displaystyle \beta (x)} が働いた場合、以下のように時間変化する。
I z → β I x I z c o s β − I y s i n β {\displaystyle I_{z}{\xrightarrow {\quad \beta I_{x}\quad }}I_{z}cos\beta -I_{y}sin\beta } 例として90度パルスの場合は、以下のようになる。
I z → ( π / 2 ) I x − I y {\displaystyle I_{z}{\xrightarrow {\quad (\pi /2)I_{x}\quad }}-I_{y}} S z → ( π / 2 ) S x − S y {\displaystyle S_{z}{\xrightarrow {\quad (\pi /2)S_{x}\quad }}-S_{y}} また13 C や15 Nの感度向上のために異種核ISにおける分極移動 を用いるINEPT では、Iスピンの磁化を回転座標系 の±x軸にそって範囲層にそろえて、Iスピンに+y軸に沿って90度パルスを与える。この回転は次のように書くことができる。
2 I x S z → ( π / 2 ) I y − 2 I z S z {\displaystyle 2I_{x}S_{z}{\xrightarrow {\quad (\pi /2)I_{y}\quad }}-2I_{z}S_{z}} これによって縦2-スピンオーダーが作り出される。
化学シフト 化学シフトによるスピンの時間変化はつぎのように表すことができる。
I y → ( π δ I t ) I z I y c o s ( 2 π δ I t ) − I x s i n ( 2 π δ I t ) {\displaystyle I_{y}{\xrightarrow {\quad (\pi \delta _{I}t)I_{z}\quad }}I_{y}cos(2\pi \delta _{I}t)-I_{x}sin(2\pi \delta _{I}t)} ここで δ I {\displaystyle \delta _{I}} は、トランスミッター周波数から測ったIスピンのシフトである。Sスピンの化学シフトも同様に扱うことができる。
異種核スピンの場合は、個々のスピンに別の回転系が設定されるから、それぞれの化学シフトはそれぞれの回転系におけるトランスミッター周波数からの相対値として表現する。
スピン結合(Jカップリング) 直積演算子法では、スピン結合はつぎのように表される。
I y → ( 2 π J I S τ ) 2 I z S z I y c o s ( 2 π J I S τ ) − 2 I x S z s i n ( 2 π J I S τ ) {\displaystyle I_{y}{\xrightarrow {\quad (2\pi J_{IS}\tau )2I_{z}S_{z}\quad }}I_{y}cos(2\pi J_{IS}\tau )-2I_{x}S_{z}sin(2\pi J_{IS}\tau )} もしも τ = 1 / ( 2 J I S ) {\displaystyle \tau =1/(2J_{IS})} となるように選ぶと、
I y → ( π / 2 ) 2 I z S z − 2 I x S z {\displaystyle I_{y}{\xrightarrow {\quad (\pi /2)2I_{z}S_{z}\quad }}-2I_{x}S_{z}} さらに次の自由歳差運動の期間を考えると、
− 2 I x S z → ( π J I S τ ) 2 I z S z − 2 I x S z c o s ( π J I S τ ) − I y s i n ( π J I S τ ) {\displaystyle -2I_{x}S_{z}{\xrightarrow {\quad (\pi J_{IS}\tau )2I_{z}S_{z}\quad }}-2I_{x}S_{z}cos(\pi J_{IS}\tau )-I_{y}sin(\pi J_{IS}\tau )} となる。もしこの時間を τ = 1 / ( 2 J I S ) {\displaystyle \tau =1/(2J_{IS})} となるように設定すると、2つのベクトルは再び-y軸にそって配向する。
− 2 I x S z → ( π / 2 ) 2 I z S z − I y {\displaystyle -2I_{x}S_{z}{\xrightarrow {\quad (\pi /2)2I_{z}S_{z}\quad }}-I_{y}}