9.化学交換のある時のスペクトル。
A成分の横磁化(FID)は
dMA+
b g
t dt=iωAMA+b g
t (54) に従う。ただしωAはAの核スピンのLarmor周波数。AとBの化学種が交換に より入れ代わるとすると各々の確率PAおよびPBはd dt
P t P t
W W
W W
P t P t
A B
A B
b g b g b g
L b g NM O
QP
=L
− −NM O
QP L NM O
QP
(55)にしたがって時間変化する。ただし、Wは A と B の間の交換速度。また(55)で は左辺をゼロと置いてやると平衡状態のAとBの確率がもとまる。これらは等 しい:PAeq. = PBeq.= 12。 核磁化の大きさは各化学種の存在確率に比例するので(54) 式と(55)式を組み合わせて
d dt
M t M t
i W W
W i W
M t M t
A B
A
B
A B +
+
+ +
L NM O
QP
=L
− −NM O
QP L
NM O
b g QP
b g b g
ω
b g
ω (56)
をえる。テキスト(8)式のようにスペクトルは FID を FT したものの実部で与え られる。すなわち
( ) { ( ) ( ) }
( ) ( )
{
21 0}
Re Re
A B
i t
A B
I I I
dt M t M t e ω
π
ω + ω + ω
∞ + + −
= +
⎡ ⎤
=
∫
⎣ + ⎦ (57)となる。(56)の両辺にe−i tω をかけて積分しIA+
( )
ω とIB+( )
ω を与える式を導出せよ。そして(57)式を使ってスペクトルを計算せよ。ただし時間の長いところで横磁化 が ゼ ロ に 減 衰 す る 条 件 : MA+
b g
∞ = MB+b g
∞ =0 お よ び 横 磁 化 の 初 期 値 がMA+
b g
0 = PAeq.= MB+b g
0 = PBeq.= 12で与えられるという条件を用いる。10.量子力学によるスピン(S=1/2)の記述。
スピン角運動量演算子(プランク定数を単位として測る。)の3つの成分は Sz= ±12の状態 Sz;± を基底関数にとると次式で与えられる。
; ; ; ;
; ; ; ;
S S S S S S S
S S S S S S
z
z z z z z z
z z z z z z
= + + + −
− + − −
L NMM O
11.密度行列 密度行列を
ρ= =
L
N MM O
Q PP
+ + ∗ + − ∗
− + ∗ − − ∗
Ψ Ψ c c c c
c c c c
z z z z
z z z z
c h c h
c h c h
(60)で定義する。
(a)Sα,α =x y z, , の期待値は
Ψ Sα Ψ =Tr S
n s
αρ (61)で与えられることを示せ。ただしTr A Ann
n
n s
=∑
は行列の対角成分の和を意味す る。(b)一般に二つの行列に対してTr A B
n s
+ =0となる時、これらの行列は直行し ているという。ただしA+は行列Aのエルミート共役で複素共役をとり転置行列 にするという意味である。(A nm Amn+ ∗
=
b g
)単位行列EとSα,α = x y z, , が互いに 直行することを示せ。(c)cz+ =cosθおよびcz− =sinθeiφと置いた時の密度行列を
ρ=C E0 +C Sx x+C Sy y+C Sz z (62) のように4つの直行行列の重ね合わせで表し係数を求めよ。密度行列で記述す ると磁化がどの方向を向いているのかが直感でわかるので有利である。
12.系の時間発展
波動の時間依存性は一般にϕ≈exp
b g
−i tω ϕ0のように表される。Planck-Eisteinの関 係(E= ω;Eは波動関数のエネルギー)使うとω=E である。量子力学では エネルギーをハミルトニアンと呼ばれる行列(またはオペレーター)Hで表す。ただしここではHˆ をNMRの理論でよくやるようにE =ωの単位で示す。波動 関数の時間発展は
Ψ
b g
t =expn s b g
−iHt Ψ 0 (63)あるいは
∂Ψ
b g
t ∂ = −t iH Ψb g
t (64) で与えられる。U tb g
=expn s
−iHt を時間発展演算子という。ここに現れた行列の 指数関数とはスカラー量と同様にexp !
U t iHt
n iHt
n
b g
=n s
− =n s
− n=
∑
∞ 1 0(65)
で与えられる。行列の計算でスカラー量と異なる点は積の交換が成立しない点
(AB ≠ BA )である。
(a)化学シフト相互作用のハミルトニアンはHCS= Δω0Sz、回転系でみたX方 向のラジオ波とスピンの相互作用ハミルトニアンはHX-rf =ω1Sxで与えられる。
exp
UCS
b g
t =n
−iH tCSs
とUX-rfb g
t =expn
−iHX-rfts
の行列を求めよ。ヒント:次の関係をもちい、奇数次の項と偶数次の項を別々にまとめよ。
1 4
2 1
4 1
4
ˆ ˆ 0 , , ,
Sα = E=⎡⎢0 ⎤⎥ α =x y z
⎣ ⎦ (66)
(b)密度行列の時間発展は
ρ
b g b g b g
t = Ψ t Ψ t =U Ψb g b g
0 Ψ 0 U−1=Uρb g
0U−1 (67) で与えられる。系の初期状態がρb g
0 =Sxの時、化学シフトにより時間推進させる とcos cos
UCS
b g b g b g
t ρ 0UCS t 1 Sxb
ω0tg
Syb
ω0tg
− = Δ + Δ (68)
となることを示せ。またρ
b g
0 =Szに対しUX-rfb g
t による時間推進を計算せよ。13.スピン-スピン結合
スピン-スピン結合のハミルトニアンは
HJ =2πJS S 1z 2z (69)
で与えられる.ここでSkzはk番目のスピンのz成分である。
(a) S1z,S2z; ,± ± を基底にとりHJを4 4× の行列で書き表せ。同じ基底を用いて S1xの行列を書き表せ。
(b)U t S U tJ
b g
1x Jb g
−1 =S1xcosb g
πJt +2S S1y 2zsinb g
πJt が成立することを示せ。[4,9-11]14.J結合のある時の1次元スペクトルの計算。
(90°)y
上の図の実験のFIDは
S t+
b g
=Tr So d
1x+S2x+iS1y+iS2yi b gd
U t S1x+S2xi b g
U t −1t
(70)で与えられる。位相敏感検波を行っているので観測量はS1x+S2x+i S
d
1y+S2yi
であり、また90度yパルスの直後を初期状態と考えればρ
b g
0 =S1x +S2xである。3つ のハミルトニアンHJ、HCS1 = Δω01S1z、HCS2 = Δω02S2zが互いに可換であることを考 慮すると、全ハミルトニアンHˆ =HˆJ +HˆCS1 +HˆCS2 に対する時間発展演算子は U tb g b g b g b g
=U t UJ CS1 t UCS2 t のように書け、それぞれのハミルトニアンに対応する時 間発展演算子を別々に掛けていけば良いことがわかる。(a)問題12、13の結果を利用して(70)式を計算せよ。
(b)(70)式の右辺に横緩和の効果を取り入れるためにexp
b
−t T2g
を掛けてフーリェ変換せよ。
(c)初期状態がρ