フーリエ変換NMRのソースを表示

'''フーリエ変換NMR'''（フーリエへんかんNMR、'''FT-NMR'''）とは、[[静磁場]]中のサンプルに[[パルス]][[磁場]]を与え、その後観察されるインパルス応答である[[自由誘導減衰]] (FID) を[[フーリエ変換]]することで[[核磁気共鳴]] (NMR) の吸収スペクトルを得る手法である。

== 概要 ==
[[File:FT-NMR.jpg|thumb|400px|フーリエ変換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" を出願した<ref>{{US patent|2561490}}</ref>。[[1954年]]に[[久保亮五]]、[[冨田和久 (物理学者)|冨田和久]]らにより[[線形応答理論]]に基づいた'''フーリエ変換NMR'''の基礎理論が提唱された。[[1956年]]にRussell H. Varianがフーリエ変換NMRの概念に関して記述した "Gyromagnetic resonance methods and apparatus" を出願した<ref>{{US patent|3287629}}</ref>。[[1957年]]にフッ化カルシウムを用いてフーリエ変換NMRがはじめて測定された。[[1964年]]にヴァリアン社のパルスNMRの先駆者の一人である[[リヒャルト・R・エルンスト]]とWeston A. Andersonによってフーリエ変換NMRが開発された<ref>{{Cite book|洋書|url=|language=English|title=NMR and EPR Spectroscopy|author=Yong Zhou|publisher=Elsevier|date=2013-09-03|ISBN=9781483226699}}</ref><ref>{{Cite journal|url=|language=English|title=Application of Fourier transfom spectroscopy to magnetic resonance|author=リヒャルト・ローベルト・エルンスト, Weston A. Anderson|publisher=Review of Scientific Instruments|date=1966|pages=93-102|volume=37 |isbn=}}</ref>。エルンストは1968年に帰国して[[チューリヒ工科大学]]で1971年に{{仮リンク|ジャン・ジェーネル|en|Jean Jeener}} (Jean Jeener)が発表した[[二次元NMR]]の着想を基に二次元フーリエ変換分光法を開発して、フーリ変換NMRと多次元NMRの開発における業績<ref>{{US patent|4045723}} "Two-dimensional gyromagnetic resonance spectroscopy"</ref><ref>{{US patent|4070611}} "Gyromagnetic resonance Fourier transfom zeugmatography"</ref><ref>{{US patent|4134058}} "Selective detection of multiple quantum transitions in nuclear magnetic resonance"</ref>で1991年の[[ノーベル化学賞]]を受賞した。

== 理論 ==
以下では、磁気共鳴に関する議論を総括した[[久保亮五]]と[[冨田和久 (物理学者)|冨田和久]]の論文<ref>{{Cite journal|url=http://jpsj.ipap.jp/link?JPSJ/9/888/|language=English|title=A General Theory of Magnetic Resonance Absorption|author=[[久保亮五|Ryogo Kubo]]|author2=[[冨田和久 (物理学者)|Kazuhisa Tomita]]|journal=[http://jpsj.ipap.jp/index.html Journal of the Physical Society of Japan]|volume=1954|issue=9|pages=888-919|publisher=[[日本物理学会]]|date=1954-6-26|doi=10.1143/JPSJ.9.888}}<!--J. Phys. Soc. Jpn. 9 (1954) pp. 888-919--></ref>による理論的展開について説明する。彼らは階段型の磁場に対する応答を考えた。

===動磁化率を用いた磁化の記述===
z軸方向の静磁場<math>\mathbf{B}_0</math>に垂直にx軸方向に周期的に振動する[[RF磁場]]<math>\mathbf{B}_{1x}(t)=\mathbf{B}_{1x}\cos\omega t</math>に対する定常的な[[磁化]]の応答を考える。

<math>\mathbf{B}_{1x}(t)</math>と磁化の応答の間には[[線形応答関係]]が成立するものと仮定する。すると<math>\mathbf{B}_{1x}(t)</math>によって定常的に誘起される磁気モーメント<math>M_{ind}(t)</math>の各成分は<math>B_{1x}</math>に比例し、比例定数は<math>\mathbf{B}_{1x}(t)</math>と同位相の成分と90°遅れた成分で定義できるから、次式のように書ける。
:<math>\mathbf{M}_{ind}(t)=(\chi'\cos\omega t+\chi''\sin\omega t)\mathbf{B}_{1x}V</math>
この<math>\chi', \chi''</math>を'''動磁化率'''と呼ぶ。動磁化率<math>\chi', \chi''</math>は2階の[[テンソル]]量である。ここで複素磁化率<math>\chi^*=\chi'-i\chi''</math>を導入すると、次式のような簡単な形に書ける。
:<math>\mathbf{M}_{ind}(t)=\mathrm{Re}\{\chi^*\mathbf{B}_{1x}e^{i\omega t}\}V</math>

====xx成分====
NMRでは、[[振動磁場]]<math>\mathbf{B}_{1x}(t)</math>によって誘起される磁化<math>\mathbf{M}_{ind}(t)</math>のx成分
:<math>M_x(t)=(\chi_{xx}'\cos\omega t + \chi_{xx}''\sin \omega t )B_{1x}V \qquad(1)</math>
がとくに重要である。<math>\chi_{xx}', \chi_{xx}''</math>はただの数である。一般に問題になるのはxx成分だけなので以後、添字xxは省略することにする。

===緩和関数を用いた磁化の記述===
[[緩和関数]]<math>\phi(t-t_0)</math>を
:<math>\Delta \mathbf{M} \equiv \mathbf{M}(t)-\mathbf{M}_0=\phi(t-t_0)\mathbf{B}_1x</math>
:<math>\phi(\infty)=0, \quad \mathbf{M}(\infty)=\mathbf{M}_0</math>
を満たすテンソル量と定義する。すると<math>\mathbf{B}_{1x}</math>が<math>-\infty</math>からずっとかけられている場合の時刻tでの磁化の期待値は、足しあわせの原理により次のように書ける。
:<math>\Delta \mathbf{M} =\phi(0)\mathbf{B}_{1x}(t)
-\int_{-\infty}^{t}\phi(t-t')\frac{d\mathbf{B}_{1x}(t')}{dt'}dt' </math>

====xx成分====
磁化のx成分は次式のように書ける。
:<math>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)</math>

===動磁化率と緩和関数の関係===
(1)式と(2)式を比較すると、動磁化率のxx成分と緩和関数のxx成分についての関係式が得られる。
:<math>\chi'=\frac{\phi{xx}(0)}{V}-\frac{\omega}{V} \int_{0}^{\infty} \phi_{xx}(t)\sin \omega t dt \qquad(3)</math>
:<math>\chi''=\frac{\omega}{V} \int_{0}^{\infty} \phi_{xx}(t)\cos\omega t dt</math>
したがって緩和関数がわかると、その[[フーリエ正弦変換|正弦]]あるいは[[フーリエ余弦変換|余弦フーリエ変換]]によって動磁化率<math>\chi', \chi''</math>が求まることになる。

===量子統計力学を用いた緩和関数の決定===
{{main|緩和関数}}
緩和関数が求まれば動磁化率が求まる。動磁化率が求まればエネルギー吸収速度や吸収係数を求めることが出来る(後述)。緩和関数の具体的な中身を知るには量子統計力学が必要である。つまり[[シュレディンガー方程式]](または[[フォン・ノイマン方程式]])を解かなければならない。

<math>t=-\infty\sim t_0</math>の間の状態は熱平衡状態であったとする．このときの状態は次の[[密度行列]]で記述される．
:<math>\rho'=\frac{e^{-\beta(\hat{H}-\mathbf{B}_{1x}\mathbf{M})}}{Tr[e^{-\beta(\hat{H}-\mathbf{B}_{1x}\mathbf{M})}]}</math>
この時間発展はフォン・ノイマン方程式で記述されるので，時刻t<math>t(\gg t_0)</math>での磁化の[[期待値]]は次式のように書ける。
:<math>\mathbf{M}(t)=Tr[\rho'(t)\mathbf{M}]=Tr[\rho'\mathbf{M}(t-t_0)]</math>
ここで<math>\mathbf{M}(t-t_0)</math>は[[ハイゼンベルク描像]]での磁化である．

<math>\hat{H}\gg\mathbf{B}_{1x}\mathbf{M}</math>と近似することで次式を得る。
:<math>\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</math> 
緩和関数の定義式と比較すると、動磁化率を定量化するために必要な緩和関数が次式で得られる。
:<math>\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</math>

====xx成分====
動磁化率のxx成分を定義するために必要な<math>\phi</math>のxx成分は、<math>\mathbf{M}_0</math>のx成分はゼロであることから次式のようになる。
:<math>\phi_{xx}(\tau)=\int_{0}^{\beta} \langle M_x M_x(\tau+i\hbar\lambda) \rangle d\lambda </math>
ただし、<math>\mathbf{B}_{1x}</math>の無い場合の[[オブザーバブル]]の期待値を<math>\langle\quad\rangle</math>と表記した。

====高温近似====
これは[[高温近似]]<math>\beta \to 0</math>が成立する場合は、次のような簡単な形に書ける。
:<math>\phi_{xx}(\tau)=\beta \langle M_x M_x (\tau) \rangle \qquad(4)</math>

===量子統計力学を用いた動磁化率の決定===
高温近似が成立する場合、(3)式と(4)式より動磁化率<math>\chi', \chi''</math>のxx成分は
:<math>\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 t dt\right] </math>
:<math>\chi'' = \frac{\omega}{kTV}\int_{0}^{\infty} \langle M_x M_x(t) \rangle \cos \omega t dt</math>

===相関関数の導入===
ここで<math>M_x</math>の'''[[相関関数]]'''<math>G(t)</math>
:<math>G(t)=\langle M_x(\tau+t)M_x(\tau)\rangle</math>
を，時間差tの偶関数として定義する．すると動磁化率の虚部は次のように簡単に書ける。
:<math>\chi''(\omega)=\frac{\omega}{2kTV}\int_{-\infty}^{\infty}G(t)e^{-i\omega t}dt</math>

===エネルギー吸収速度・吸収係数===
z軸方向の静磁場<math>\mathbf{B}_0</math>に加えて振動磁場<math>\mathbf{B}_{1x}(t)</math>をかけた後、十分に時間が経過した後での定常状態におけるx方向の誘起磁化<math>M_x(t)</math>の時間変動より得られる<math>\chi', \chi''</math>より、エネルギーの吸収速度<math>Q(\omega)</math>や吸収係数<math>A(\omega)</math>（[[エネルギー散逸]]、またはパワーロス）は振動地場の周波数<math>\omega</math>の関数として次式で与えられる。
:<math>Q(\omega)=\frac{\omega\chi'}{2}'(\omega)B_{1x}^2</math>
:<math>A(\omega)=\frac{\omega}{2\pi}\chi''(\omega)</math>
よって吸収スペクトルは以下のように相関関数のフーリエ変換として表せる。
:<math>A(\omega)=\frac{\omega^2}{4\pi kTV}\int_{-\infty}^{\infty}G(t)e^{-i\omega t}dt</math>
逆に相関関数は，動磁化率の虚部<math>\chi''</math>または吸収スペクトル<math>A(\omega)</math>の逆フーリエ変換として表せる。
:<math>G(t)=\frac{kTV}{\pi}\int_{-\infty}^{\infty}\frac{\chi''(\omega)}{\omega}e^{i\omega t}d\omega</math>
:<math>G(t)=2kTV\int_{-\infty}^{\infty}\frac{A(\omega)}{\omega^2}e^{i\omega t}d\omega</math>
<math>G(t)</math>を偶関数として定義したので，<math>\chi''</math>は奇関数，<math>A(\omega)</math>は偶関数である。

== 注意点 ==
FT-NMRの理論は、線形系に対する一般原理「周波数応答関数は系のインパルス応答関数とフーリエ変換で結ばれる」を基にしている。スピン系では、単一のパルス入力についてであればこのことは成立するが、一般の多重パルス、多重共鳴法では系の線形性は成り立たない。通常、周波数スペクトルと呼んでいるものは、弱い正弦波入力にたいする応答である。非線形な系でも、入力が十分に弱ければ、入力に対して二次以上の依存性を示す項（非線形項）の存在を無視することができて、線形系としての取り扱いが可能である。しかし強いラジオ波パルスを入力した場合には、その応答は一般に非線形となる。

== 出典 ==
{{脚注ヘルプ}}
{{reflist}}

==参考文献==
{{参照方法|section=1|date=2018年1月}}
*{{Cite book|洋書|url=|language=English|title=Pulse and Fourier Transform NMR Introduction to Theory and Methods|author=[[Thomas C. Farrar]], [[Edwin D. Becker]]|publisher=[[Academic Press]]|date=1971-9-11|ISBN=978-0-122-49650-9}}
**{{Cite book|和書|url=http://www3.ocn.ne.jp/~yoshioka/single_titles/ISBN4-8427-0181-1.html|language=日本語|title=パルスおよびフーリエ変換NMR 理論および方法への入門|author1=Thomas C. Farrar|authorlink1=Thomas C. Farrar|author2=Edwin D. Becker|authorlink2=Edwin D. Becker|translator=[[赤坂一之]], [[井元敏明]]|work=現代科学|publisher=[[吉岡書店]]|date=1979-1-1|ISBN=978-4-842-70181-3}}
*{{Cite book|和書|url=http://www.kyoritsu-pub.co.jp/bookdetail/9784320042124|language=日本語|title=核磁気共鳴の基礎と原理|author=北丸竜三|authorlink=北丸竜三|publisher=[[共立出版]]|date=1987-9-1|ISBN=978-4-320-04212-4}}
*{{Cite book|洋書|url=http://ukcatalogue.oup.com/product/9780198555674.do|language=English|title=Modern NMR Spectroscopy A Guide for Chemists|author=[[Jeremy K. M. Sanders]], [[Brian K. Hunter]]|publisher=[[オックスフォード大学出版局|Oxford University Press]]|date=1993-3-25|ISBN=978-0-19-855567-4}}
**{{Cite book|和書|url=|language=日本語|title=NMRガイドブック|author1=Jeremy K. M. Sanders|authorlink1=Jeremy K. M. Sanders|author2=Brian K. Hunter|authorlink2=Brian K. Hunter|translator=[[坂口潮]]|publisher=[[廣川書店]]|date=1992-1-1|ISBN=978-4-567-25450-2}}
<!--*{{Cite book|洋書|url=http://global.oup.com/ushe/product/modern-nmr-spectroscopy-9780198555674?tab=|language=English|title=Modern NMR Spectroscopy A Guide for Chemists|author=[[Jeremy K. M. Sanders]], [[Brian K. Hunter]]|publisher=[[Oxford University Press]]|date=1993-3-25|ISBN=978-0-19-855567-4}}-->
*{{Cite book|洋書|url=|language=English|title=A Handbook of Nuclear Magnetic Resonance|author=[[Ray Freeman]]|publisher=[[Longman Publishing]]|date=1997-4-1|ISBN=978-0-582-25184-7}}
**{{Cite book|和書|url=http://www.kyoritsu-pub.co.jp/bookdetail/9784320042858|language=日本語|title=NMRハンドブック|author=Ray Freeman|authorlink=Ray Freeman|translator=[[坂口潮]], [[嶋田一夫]], [[荒田洋治]]|publisher=[[共立出版]]|date=1992-3-10|ISBN=978-4-320-04285-8}}

== 関連項目 ==
*[[核磁気共鳴]]
*[[二次元NMR]]
*[[電子スピン共鳴]](ESR)
*[[固体核磁気共鳴]](固体NMR)
*[[線形応答理論]]

{{DEFAULTSORT:ふうりえへんかんえぬえむあある}}

[[Category:核磁気共鳴]]
[[Category:量子力学]]