フーリエ変換NMR (フーリエへんかんNMR、FT-NMR )とは、静磁場 中のサンプルにパルス 磁場 を与え、その後観察されるインパルス応答である自由誘導減衰 (FID) をフーリエ変換 することで核磁気共鳴 (NMR) の吸収スペクトルを得る手法である。
概要
フーリエ変換NMRの概要
連続波(CW)法では様々な周波数 の磁場 を掃引しながら系に与える。一方でフーリエ変換NMRではパルス磁場を系に与える。この2つの入力、つまり「多数の周波数成分の和」と「パルス」は、互いにフーリエ変換の関係になっている。また系が線形応答 するときは、入力の間の数学的関係(つまりフーリエ変換)が、そのまま出力の数学的関係に伝達される。つまりCW法で得られる応答と、パルス法で得られる応答との間にもフーリエ変換の関係が成り立っている。このことからパルス磁場を与えたときのインパルス応答をフーリエ変換することにより、CW法で得られる吸収曲線と分散曲線と同じものが得られる。
線形応答理論 によればインパルス応答関数 のフーリエ変換は周波数応答関数 を与える。周波数応答関数はある周波数の電磁波が吸収される程度を表す関数であるから、これはNMRスペクトルに他ならない。それゆえにインパルス(パルス状の電磁波)を試料に当ててすべての核を一斉に励起し、その結果生じる磁化ベクトルの変化、すなわち自由誘導減衰 (Free Induction Decay, FID) を測定し、これをフーリエ変換したものは虚部が分散スペクトル 、実部が吸収スペクトル (エネルギー散逸 、パワーロス)になっている。またこの2つのスペクトルの間にはクラマース・クローニッヒの関係 が成立する。
NMRのほとんどの応用は、完全なNMRスペクトル、つまり周波数の関数としてのNMRシグナルの強度を含む。単純な連続波 (CW) 法よりも効率的にNMRスペクトルを得るための初期の試みでは、2つ以上の周波数で同時に標的に光を当てる手法が使われていた。NMRにおける革命は、高周波の短パルスが使われ始めた時に起こった。簡単に言えば、任意の「キャリア」周波数の矩形パルスは、さまざまな範囲の周波数を「含んで」おり、励起の幅(バンド幅)はパルスの持続時間と反比例する。近似方形波のフーリエ変換は、主周波数の隣接領域における全ての周波数の寄与を含んでいる。NMR周波数の範囲が制限されていることによって、全NMRスペクトルを励起するための短い(ミリ秒からマイクロ秒)高周波パルスを使うことが比較的容易となっている。
こういったパルスを一連の核 スピン に印加すると、全ての単一量子 NMR遷移 が同時に励起される。総磁化ベクトルの観点からは、(外部磁場に沿って並んだ)平衡位置から磁化ベクトルが傾くことに対応する。平衡から外れた磁化ベクトルは、スピンのNMR周波数における外部磁場ベクトルに対して歳差運動する。この周期的に振動する磁化ベクトルはすぐ近くの検出コイルに電流 を誘導 し、NMR周波数における電気シグナルの周期的な振動を作る。このシグナルは自由誘導減衰 (FID) として知られており、全ての励起スピンからのNMR応答のベクトル和を含んでいる。周波数領域のNMRスペクトル(NMR吸収強度 vs. NMR周波数)を得るためには、この時間領域シグナル(強度 vs. 時間)をフーリエ変換しなければならない。幸運なことに、フーリエ変換NMRの開発は、デジタルコンピュータやデジタル高速フーリエ変換 の開発と同時期に起こった。フーリエ法は多くの分光法 の種類に適応することができる。
歴史
1948年 にRussell H. Varianがヴァリアン・アソシエイツ を設立して自由誘導減衰 信号の検出に関して記述した "Method and means for correlating nuclear properties of atoms and magnetic fields" を出願した[ 1] 。1954年 に久保亮五 、冨田和久 らにより線形応答理論 に基づいたフーリエ変換NMR の基礎理論が提唱された。1956年 にRussell H. Varianがフーリエ変換NMRの概念に関して記述した "Gyromagnetic resonance methods and apparatus" を出願した[ 2] 。1957年 にフッ化カルシウムを用いてフーリエ変換NMRがはじめて測定された。1964年 にヴァリアン社のパルスNMRの先駆者の一人であるリヒャルト・R・エルンスト とWeston A. Andersonによってフーリエ変換NMRが開発された[ 3] [ 4] 。エルンストは1968年に帰国してチューリヒ工科大学 で1971年にジャン・ジェーネル (英語版 ) (Jean Jeener)が発表した二次元NMR の着想を基に二次元フーリエ変換分光法を開発して、フーリ変換NMRと多次元NMRの開発における業績[ 5] [ 6] [ 7] で1991年のノーベル化学賞 を受賞した。
理論
以下では、磁気共鳴に関する議論を総括した久保亮五 と冨田和久 の論文[ 8] による理論的展開について説明する。彼らは階段型の磁場に対する応答を考えた。
動磁化率を用いた磁化の記述
z軸方向の静磁場
B
0
{\displaystyle \mathbf {B} _{0}}
に垂直にx軸方向に周期的に振動するRF磁場
B
1
x
(
t
)
=
B
1
x
cos
-->
ω ω -->
t
{\displaystyle \mathbf {B} _{1x}(t)=\mathbf {B} _{1x}\cos \omega t}
に対する定常的な磁化 の応答を考える。
B
1
x
(
t
)
{\displaystyle \mathbf {B} _{1x}(t)}
と磁化の応答の間には線形応答関係 が成立するものと仮定する。すると
B
1
x
(
t
)
{\displaystyle \mathbf {B} _{1x}(t)}
によって定常的に誘起される磁気モーメント
M
i
n
d
(
t
)
{\displaystyle M_{ind}(t)}
の各成分は
B
1
x
{\displaystyle B_{1x}}
に比例し、比例定数は
B
1
x
(
t
)
{\displaystyle \mathbf {B} _{1x}(t)}
と同位相の成分と90°遅れた成分で定義できるから、次式のように書ける。
M
i
n
d
(
t
)
=
(
χ χ -->
′
cos
-->
ω ω -->
t
+
χ χ -->
″
sin
-->
ω ω -->
t
)
B
1
x
V
{\displaystyle \mathbf {M} _{ind}(t)=(\chi '\cos \omega t+\chi ''\sin \omega t)\mathbf {B} _{1x}V}
この
χ χ -->
′
,
χ χ -->
″
{\displaystyle \chi ',\chi ''}
を動磁化率 と呼ぶ。動磁化率
χ χ -->
′
,
χ χ -->
″
{\displaystyle \chi ',\chi ''}
は2階のテンソル 量である。ここで複素磁化率
χ χ -->
∗ ∗ -->
=
χ χ -->
′
− − -->
i
χ χ -->
″
{\displaystyle \chi ^{*}=\chi '-i\chi ''}
を導入すると、次式のような簡単な形に書ける。
M
i
n
d
(
t
)
=
R
e
{
χ χ -->
∗ ∗ -->
B
1
x
e
i
ω ω -->
t
}
V
{\displaystyle \mathbf {M} _{ind}(t)=\mathrm {Re} \{\chi ^{*}\mathbf {B} _{1x}e^{i\omega t}\}V}
xx成分
NMRでは、振動磁場
B
1
x
(
t
)
{\displaystyle \mathbf {B} _{1x}(t)}
によって誘起される磁化
M
i
n
d
(
t
)
{\displaystyle \mathbf {M} _{ind}(t)}
のx成分
M
x
(
t
)
=
(
χ χ -->
x
x
′
cos
-->
ω ω -->
t
+
χ χ -->
x
x
″
sin
-->
ω ω -->
t
)
B
1
x
V
(
1
)
{\displaystyle M_{x}(t)=(\chi _{xx}'\cos \omega t+\chi _{xx}''\sin \omega t)B_{1x}V\qquad (1)}
がとくに重要である。
χ χ -->
x
x
′
,
χ χ -->
x
x
″
{\displaystyle \chi _{xx}',\chi _{xx}''}
はただの数である。一般に問題になるのはxx成分だけなので以後、添字xxは省略することにする。
緩和関数を用いた磁化の記述
緩和関数
ϕ ϕ -->
(
t
− − -->
t
0
)
{\displaystyle \phi (t-t_{0})}
を
Δ Δ -->
M
≡ ≡ -->
M
(
t
)
− − -->
M
0
=
ϕ ϕ -->
(
t
− − -->
t
0
)
B
1
x
{\displaystyle \Delta \mathbf {M} \equiv \mathbf {M} (t)-\mathbf {M} _{0}=\phi (t-t_{0})\mathbf {B} _{1}x}
ϕ ϕ -->
(
∞ ∞ -->
)
=
0
,
M
(
∞ ∞ -->
)
=
M
0
{\displaystyle \phi (\infty )=0,\quad \mathbf {M} (\infty )=\mathbf {M} _{0}}
を満たすテンソル量と定義する。すると
B
1
x
{\displaystyle \mathbf {B} _{1x}}
が
− − -->
∞ ∞ -->
{\displaystyle -\infty }
からずっとかけられている場合の時刻tでの磁化の期待値は、足しあわせの原理により次のように書ける。
Δ Δ -->
M
=
ϕ ϕ -->
(
0
)
B
1
x
(
t
)
− − -->
∫ ∫ -->
− − -->
∞ ∞ -->
t
ϕ ϕ -->
(
t
− − -->
t
′
)
d
B
1
x
(
t
′
)
d
t
′
d
t
′
{\displaystyle \Delta \mathbf {M} =\phi (0)\mathbf {B} _{1x}(t)-\int _{-\infty }^{t}\phi (t-t'){\frac {d\mathbf {B} _{1x}(t')}{dt'}}dt'}
xx成分
磁化のx成分は次式のように書ける。
M
x
(
t
)
=
{
[
ϕ ϕ -->
x
x
(
0
)
− − -->
∫ ∫ -->
0
∞ ∞ -->
ω ω -->
ϕ ϕ -->
x
x
(
t
)
(
τ τ -->
)
sin
-->
ω ω -->
τ τ -->
d
τ τ -->
]
cos
-->
ω ω -->
t
+
[
∫ ∫ -->
0
∞ ∞ -->
ω ω -->
ϕ ϕ -->
x
x
(
t
)
(
τ τ -->
)
cos
-->
ω ω -->
τ τ -->
d
τ τ -->
]
sin
-->
ω ω -->
t
}
B
1
x
(
2
)
{\displaystyle M_{x}(t)=\left\{\left[\phi _{xx}(0)-\int _{0}^{\infty }\omega \phi _{xx}(t)(\tau )\sin \omega \tau d\tau \right]\cos \omega t+\left[\int _{0}^{\infty }\omega \phi _{xx}(t)(\tau )\cos \omega \tau d\tau \right]\sin \omega t\right\}B_{1x}\qquad (2)}
動磁化率と緩和関数の関係
(1)式と(2)式を比較すると、動磁化率のxx成分と緩和関数のxx成分についての関係式が得られる。
χ χ -->
′
=
ϕ ϕ -->
x
x
(
0
)
V
− − -->
ω ω -->
V
∫ ∫ -->
0
∞ ∞ -->
ϕ ϕ -->
x
x
(
t
)
sin
-->
ω ω -->
t
d
t
(
3
)
{\displaystyle \chi '={\frac {\phi {xx}(0)}{V}}-{\frac {\omega }{V}}\int _{0}^{\infty }\phi _{xx}(t)\sin \omega tdt\qquad (3)}
χ χ -->
″
=
ω ω -->
V
∫ ∫ -->
0
∞ ∞ -->
ϕ ϕ -->
x
x
(
t
)
cos
-->
ω ω -->
t
d
t
{\displaystyle \chi ''={\frac {\omega }{V}}\int _{0}^{\infty }\phi _{xx}(t)\cos \omega tdt}
したがって緩和関数がわかると、その正弦 あるいは余弦フーリエ変換 によって動磁化率
χ χ -->
′
,
χ χ -->
″
{\displaystyle \chi ',\chi ''}
が求まることになる。
量子統計力学を用いた緩和関数の決定
緩和関数が求まれば動磁化率が求まる。動磁化率が求まればエネルギー吸収速度や吸収係数を求めることが出来る(後述)。緩和関数の具体的な中身を知るには量子統計力学が必要である。つまりシュレディンガー方程式 (またはフォン・ノイマン方程式 )を解かなければならない。
t
=
− − -->
∞ ∞ -->
∼ ∼ -->
t
0
{\displaystyle t=-\infty \sim t_{0}}
の間の状態は熱平衡状態であったとする.このときの状態は次の密度行列 で記述される.
ρ ρ -->
′
=
e
− − -->
β β -->
(
H
^ ^ -->
− − -->
B
1
x
M
)
T
r
[
e
− − -->
β β -->
(
H
^ ^ -->
− − -->
B
1
x
M
)
]
{\displaystyle \rho '={\frac {e^{-\beta ({\hat {H}}-\mathbf {B} _{1x}\mathbf {M} )}}{Tr[e^{-\beta ({\hat {H}}-\mathbf {B} _{1x}\mathbf {M} )}]}}}
この時間発展はフォン・ノイマン方程式で記述されるので,時刻t
t
(
≫ ≫ -->
t
0
)
{\displaystyle t(\gg t_{0})}
での磁化の期待値 は次式のように書ける。
M
(
t
)
=
T
r
[
ρ ρ -->
′
(
t
)
M
]
=
T
r
[
ρ ρ -->
′
M
(
t
− − -->
t
0
)
]
{\displaystyle \mathbf {M} (t)=Tr[\rho '(t)\mathbf {M} ]=Tr[\rho '\mathbf {M} (t-t_{0})]}
ここで
M
(
t
− − -->
t
0
)
{\displaystyle \mathbf {M} (t-t_{0})}
はハイゼンベルク描像 での磁化である.
H
^ ^ -->
≫ ≫ -->
B
1
x
M
{\displaystyle {\hat {H}}\gg \mathbf {B} _{1x}\mathbf {M} }
と近似することで次式を得る。
M
(
t
)
=
M
0
+
B
1
x
T
r
[
ρ ρ -->
0
∫ ∫ -->
0
β β -->
e
λ λ -->
H
^ ^ -->
M
e
− − -->
λ λ -->
H
^ ^ -->
M
(
t
− − -->
t
0
)
d
λ λ -->
]
− − -->
β β -->
B
1
x
M
0
M
0
{\displaystyle \mathbf {M} (t)=\mathbf {M} _{0}+\mathbf {B} _{1x}Tr\left[\rho _{0}\int _{0}^{\beta }e^{\lambda {\hat {H}}}\mathbf {M} e^{-\lambda {\hat {H}}}\mathbf {M} (t-t_{0})\,d\lambda \right]-\beta \mathbf {B} _{1x}\mathbf {M} _{0}\mathbf {M} _{0}}
緩和関数の定義式と比較すると、動磁化率を定量化するために必要な緩和関数が次式で得られる。
ϕ ϕ -->
(
τ τ -->
)
=
∫ ∫ -->
0
β β -->
T
r
[
ρ ρ -->
0
M
M
(
τ τ -->
+
i
ℏ ℏ -->
λ λ -->
)
]
d
λ λ -->
− − -->
β β -->
M
0
M
0
{\displaystyle \phi (\tau )=\int _{0}^{\beta }Tr[\rho _{0}\mathbf {M} \mathbf {M} (\tau +i\hbar \lambda )]d\lambda -\beta \mathbf {M} _{0}\mathbf {M} _{0}}
xx成分
動磁化率のxx成分を定義するために必要な
ϕ ϕ -->
{\displaystyle \phi }
のxx成分は、
M
0
{\displaystyle \mathbf {M} _{0}}
のx成分はゼロであることから次式のようになる。
ϕ ϕ -->
x
x
(
τ τ -->
)
=
∫ ∫ -->
0
β β -->
⟨ ⟨ -->
M
x
M
x
(
τ τ -->
+
i
ℏ ℏ -->
λ λ -->
)
⟩ ⟩ -->
d
λ λ -->
{\displaystyle \phi _{xx}(\tau )=\int _{0}^{\beta }\langle M_{x}M_{x}(\tau +i\hbar \lambda )\rangle d\lambda }
ただし、
B
1
x
{\displaystyle \mathbf {B} _{1x}}
の無い場合のオブザーバブル の期待値を
⟨ ⟨ -->
⟩ ⟩ -->
{\displaystyle \langle \quad \rangle }
と表記した。
高温近似
これは高温近似
β β -->
→ → -->
0
{\displaystyle \beta \to 0}
が成立する場合は、次のような簡単な形に書ける。
ϕ ϕ -->
x
x
(
τ τ -->
)
=
β β -->
⟨ ⟨ -->
M
x
M
x
(
τ τ -->
)
⟩ ⟩ -->
(
4
)
{\displaystyle \phi _{xx}(\tau )=\beta \langle M_{x}M_{x}(\tau )\rangle \qquad (4)}
量子統計力学を用いた動磁化率の決定
高温近似が成立する場合、(3)式と(4)式より動磁化率
χ χ -->
′
,
χ χ -->
″
{\displaystyle \chi ',\chi ''}
のxx成分は
χ χ -->
′
=
1
k
T
V
[
⟨ ⟨ -->
M
x
M
x
(
t
)
⟩ ⟩ -->
− − -->
ω ω -->
∫ ∫ -->
0
∞ ∞ -->
⟨ ⟨ -->
M
x
M
x
(
t
)
⟩ ⟩ -->
sin
-->
ω ω -->
t
d
t
]
{\displaystyle \chi '={\frac {1}{kTV}}\left[\langle M_{x}M_{x}(t)\rangle -\omega \int _{0}^{\infty }\langle M_{x}M_{x}(t)\rangle \sin \omega tdt\right]}
χ χ -->
″
=
ω ω -->
k
T
V
∫ ∫ -->
0
∞ ∞ -->
⟨ ⟨ -->
M
x
M
x
(
t
)
⟩ ⟩ -->
cos
-->
ω ω -->
t
d
t
{\displaystyle \chi ''={\frac {\omega }{kTV}}\int _{0}^{\infty }\langle M_{x}M_{x}(t)\rangle \cos \omega tdt}
相関関数の導入
ここで
M
x
{\displaystyle M_{x}}
の相関関数
G
(
t
)
{\displaystyle G(t)}
G
(
t
)
=
⟨ ⟨ -->
M
x
(
τ τ -->
+
t
)
M
x
(
τ τ -->
)
⟩ ⟩ -->
{\displaystyle G(t)=\langle M_{x}(\tau +t)M_{x}(\tau )\rangle }
を,時間差tの偶関数として定義する.すると動磁化率の虚部は次のように簡単に書ける。
χ χ -->
″
(
ω ω -->
)
=
ω ω -->
2
k
T
V
∫ ∫ -->
− − -->
∞ ∞ -->
∞ ∞ -->
G
(
t
)
e
− − -->
i
ω ω -->
t
d
t
{\displaystyle \chi ''(\omega )={\frac {\omega }{2kTV}}\int _{-\infty }^{\infty }G(t)e^{-i\omega t}dt}
エネルギー吸収速度・吸収係数
z軸方向の静磁場
B
0
{\displaystyle \mathbf {B} _{0}}
に加えて振動磁場
B
1
x
(
t
)
{\displaystyle \mathbf {B} _{1x}(t)}
をかけた後、十分に時間が経過した後での定常状態におけるx方向の誘起磁化
M
x
(
t
)
{\displaystyle M_{x}(t)}
の時間変動より得られる
χ χ -->
′
,
χ χ -->
″
{\displaystyle \chi ',\chi ''}
より、エネルギーの吸収速度
Q
(
ω ω -->
)
{\displaystyle Q(\omega )}
や吸収係数
A
(
ω ω -->
)
{\displaystyle A(\omega )}
(エネルギー散逸 、またはパワーロス)は振動地場の周波数
ω ω -->
{\displaystyle \omega }
の関数として次式で与えられる。
Q
(
ω ω -->
)
=
ω ω -->
χ χ -->
′
2
′
(
ω ω -->
)
B
1
x
2
{\displaystyle Q(\omega )={\frac {\omega \chi '}{2}}'(\omega )B_{1x}^{2}}
A
(
ω ω -->
)
=
ω ω -->
2
π π -->
χ χ -->
″
(
ω ω -->
)
{\displaystyle A(\omega )={\frac {\omega }{2\pi }}\chi ''(\omega )}
よって吸収スペクトルは以下のように相関関数のフーリエ変換として表せる。
A
(
ω ω -->
)
=
ω ω -->
2
4
π π -->
k
T
V
∫ ∫ -->
− − -->
∞ ∞ -->
∞ ∞ -->
G
(
t
)
e
− − -->
i
ω ω -->
t
d
t
{\displaystyle A(\omega )={\frac {\omega ^{2}}{4\pi kTV}}\int _{-\infty }^{\infty }G(t)e^{-i\omega t}dt}
逆に相関関数は,動磁化率の虚部
χ χ -->
″
{\displaystyle \chi ''}
または吸収スペクトル
A
(
ω ω -->
)
{\displaystyle A(\omega )}
の逆フーリエ変換として表せる。
G
(
t
)
=
k
T
V
π π -->
∫ ∫ -->
− − -->
∞ ∞ -->
∞ ∞ -->
χ χ -->
″
(
ω ω -->
)
ω ω -->
e
i
ω ω -->
t
d
ω ω -->
{\displaystyle G(t)={\frac {kTV}{\pi }}\int _{-\infty }^{\infty }{\frac {\chi ''(\omega )}{\omega }}e^{i\omega t}d\omega }
G
(
t
)
=
2
k
T
V
∫ ∫ -->
− − -->
∞ ∞ -->
∞ ∞ -->
A
(
ω ω -->
)
ω ω -->
2
e
i
ω ω -->
t
d
ω ω -->
{\displaystyle G(t)=2kTV\int _{-\infty }^{\infty }{\frac {A(\omega )}{\omega ^{2}}}e^{i\omega t}d\omega }
G
(
t
)
{\displaystyle G(t)}
を偶関数として定義したので,
χ χ -->
″
{\displaystyle \chi ''}
は奇関数,
A
(
ω ω -->
)
{\displaystyle A(\omega )}
は偶関数である。
注意点
FT-NMRの理論は、線形系に対する一般原理「周波数応答関数は系のインパルス応答関数とフーリエ変換で結ばれる」を基にしている。スピン系では、単一のパルス入力についてであればこのことは成立するが、一般の多重パルス、多重共鳴法では系の線形性は成り立たない。通常、周波数スペクトルと呼んでいるものは、弱い正弦波入力にたいする応答である。非線形な系でも、入力が十分に弱ければ、入力に対して二次以上の依存性を示す項(非線形項)の存在を無視することができて、線形系としての取り扱いが可能である。しかし強いラジオ波パルスを入力した場合には、その応答は一般に非線形となる。
出典
^ アメリカ合衆国特許第 2,561,490号
^ アメリカ合衆国特許第 3,287,629号
^ Yong Zhou (2013-09-03) (English). NMR and EPR Spectroscopy . Elsevier. ISBN 9781483226699
^ リヒャルト・ローベルト・エルンスト, Weston A. Anderson (1966) (English). Application of Fourier transfom spectroscopy to magnetic resonance . 37 . Review of Scientific Instruments. pp. 93-102.
^ アメリカ合衆国特許第 4,045,723号 "Two-dimensional gyromagnetic resonance spectroscopy"
^ アメリカ合衆国特許第 4,070,611号 "Gyromagnetic resonance Fourier transfom zeugmatography"
^ アメリカ合衆国特許第 4,134,058号 "Selective detection of multiple quantum transitions in nuclear magnetic resonance"
^ Ryogo Kubo ; Kazuhisa Tomita (1954-6-26). “A General Theory of Magnetic Resonance Absorption” (English). Journal of the Physical Society of Japan (日本物理学会 ) 1954 (9): 888-919. doi :10.1143/JPSJ.9.888 . http://jpsj.ipap.jp/link?JPSJ/9/888/ .
参考文献
関連項目