博 士 （ 理 学 ） 朝 日 信 義

学位論文題名NrvI
RandrvIolec
ularDynamicsStudyofLi
quidand
SupercriticalOrganic Solvents（液体 ・超臨界有機溶媒のNMRおよび分子動力学による研究）

学 位論文内容の要旨
超臨界流体とは臨界点を超えた流体のことをいう。臨界点近くの流体は、気体と液体の中間の性質を示し、その構造は均一ではなく、液体のように密度の高い部分と気体のように密度の小さい部分の混合であり、密度揺らぎの大きい状態であると考えられる。超臨界流体ではマクロな物性を、温度や圧カの変化によって気体に近い性質から液体に近い性質まで幅広く連続的に変化させることができる。その特徴は、特異な溶解能や反応の活性化を制御できることであり、近年、溶媒として非常に注目されている。
液体の温度と圧カを上げていくと、臨界点を越えたところで超臨界流体となるので、超臨界流体の構造と動的な性質および溶質溶媒相互作用を調べることは、超臨界状態のみならず液体状態の研究においても重要な知見をもたらすと考えられる。
測定に用いた有機溶媒は、ベンゼン(C6H6)、塩化メチレン(CH2CI2)、ヘキサフルオロベンゼン(C6F6)、アセトニトリル(CH3CD、アセ卜ン((CH3)
2CO)
、メタノール(CH30H）である。これらの溶媒においては、それぞれ分子間相互作用や極性が異なっており、それぞれの分子運動の大きさやその性質、構造などを調べることは非常に興味がある。また、溶質溶媒相互作用を調べるために、ベンゼンーヘキサフルオロ
ベンゼン系とベンゼンーメタノール系などのニ成分系についても実験を試みた。本論文では、これら超臨界 流体 の物 性の研究を核磁気共鳴実験(NMR)と分子動力学計算(MD)をもちぃて行った。 NMRによって、密度、自己拡散係数、スピンー格子 緩和時間、化学シフトなどの物理量を得た。MDからは、自己拡 散係数 、動 径分布 関数 、相関 時間 などを 計算 し、実 験と 比較し 考察 した。
本論文は8章で構成される。
第 1章は 序論として、超臨界流 体の特徴やこれまでの基礎的研究とその応用例、NMRや分子動力学シミュレーションから得られる物理量の意味と超臨界流体の物性研究のためにそれらの手段を用いる意義 、 そ し て 本 研 究 を 行 う に 至 っ た 動 機 と そ の 目 的 お よ び 本 論 文 の 位 置 づ け に つい て述 べた 。
第 2章で は、まずNMRから得られるスピン―格子緩和のメカ ニズムである磁気双極子―双極子緩和とス ピン 回転 緩和 の緩 和 速度 を簡 単に 導出した。次にN
MRに よる自己拡散係数測定の原理として、
Carr‑Purcell‑
Meiboom‑GiIl(CPMG)
法について説明した。
第 3章で は、液体から超臨界域 を含むベンゼンについて、プ口トンのNMR測定を行った結果にっいて考察した。気液共存下における気体と液体の密度は、一次元投影法によって測定し、臨界指数をp=o．38土0.07と得た。これは、分子性流体の値として知られている理論値1/3に測定誤差の範囲内で一致している。 自己拡 散係数の 測定は CPMG法を用い て行っ た。従来 、NMRによる自 己拡散係 数の測 定にはスピンエコー法またはパルス磁場勾配法が用いられていたが、超臨界流体は一般にスピン―格子緩和時間が長いために多大の測定時間を要する。本研究においては、超臨界流体が大きい拡散係数を持っことに注目し、CPMG法を初めて超臨界流体に適用した。この方法には、時間を短縮する以外にも利点があり、得られた信頼性の高いデータは、超臨界域においてChapman‑Enskogの式から予想される値に一致した。
スピン―格子緩和時間は、温度の上昇とともにゆるやかに増加し、臨界温度付近からゆるやかに減少する結果となった。この結果から、得られた密度と拡散係数のデータを用いて、分子間緩和と分子内緩和の緩和速度を分離した。Bloembergen‑
Purcell‑Pound(BPP)
理論とHubbardの式から、分子内緩和速度の温度依存 性が極 大を示す470Kにおいて、分子内双極子一双極子相互作用とスピン回転相互作用の相関時間を求めた。得られた相関時間とスビン―格子緩和時間の密度依存性のデータから、超臨界域においてはスピン回転相互作用が緩和機構として支配的であると結論した。
第4章では、水素結合性液体であるメタノールにっいて、室温から超臨界域を含む温度範囲で行った
プロト ンNMR実験と分 子動力 学計算に っいて考 察した 。メタノ ールにっいては、NMRによって自己拡散係数と化学シフトの結果を報告した。自己拡散係数は、ベンゼンと同様の密度依存性を示し、超臨界域にお いては Chapman‑Enskogの式に一致した。このことから、本研究における温度、密度範囲においては、水素結合の効果は拡散係数に影響を与えないと結論した。化学シフトの結果から、温度の増加とともに水素結合によるメタノールクラスターは減少し、密度の滅少とともにクラスタ一数が減少することを示した。また、超臨界域におけるクラスター分布の計算を熟力学モデルを用いて行った。この結果
から、575K、49.
9kgm‑3の高温、低密度の状態においても、メタノールクラスターが存在していることがわか った。 MD計算には 、Jorgensenの3サ イ卜ポ テンシャ ルを用いた。MD計算によって、自己拡散係数、動径分布関数、メタノールクラスターのサイズ分布などを求めた。自己拡散係数の計算結果は、測定した全ての温度、密度範囲において実測値と非常によく一致した。このことから、この温度、密度範囲において、Jorgensenのポテンシヤルが妥当であると結論した。MDから得られたクラスター分布は、熱 力学 モデ ル によ っ て得 ら れた クラ ス ター 分布と 非常によい一致 を示した。
第5章では 、その他 の有機 純溶媒に ついてNMR
実験 を行った 結果について考察した。全ての流体における自己拡散係数の結果は、ベンゼンとメタノールで得られたのと同様の温度依存性を示した。これ
ら全ての超臨界流体において、本研究で測定した温度と密度範囲では、希薄気体に対するC
hapman‑Enskogの式か ら予想 される拡 散係数 と一致し た。ヘキ サフル オ口ベン ゼンとアセ卜ニ卜リルについては、
Lennard‑Jonesポテンシャルパラメータを、Lennard‑Jones
流体に対する対応状態の法則から導いた。スピン 一格 子緩 和 時間 の 結果 は 、ベ ンゼ ン と同 様の温 度、密度依存性 を示した。
第6章では 、ベンゼ ンーヘ キサフル オロベン ゼン系 のNMR実 験を行った結果を報告した。混合比は1:1．1であり、1Hと19Fにっいてそれぞれスピンー格子緩和時間と自己拡散係数を測定した。拡散係数の結果は、両成分ともに近い値となった。スピン―格子緩和時間は、ベンゼンとへキサフルオロベンゼンのそれぞれの純流体で得られた温度依存性と同様の結果を示した。
第7章では、重ベンゼンーメタノール系のプ口卜ン化学シフ卜を測定した結果を考察した。化学シフ
卜の温度依存性は純メタノール系で得られたのと同様の結果となった。このことから、水素結合クラス
タ ーの 温 度、 密 度依 存 性は メ タノ ー ル純 溶 媒の 結 果 と同 様 であ る と結 論され た。
第8章において、本研究における総括を行った。

学位論文審査の要旨
主 査 教 授 中 村 義 男
副 査 教 授 井 川 駿 一
副 査 教 授 稲 辺保
学位論文題名NrvIR and N/Iolecular DynamlCSStudyofLiquidand
S
uper
Cri
tiCalorgani
CSOlVentS（液体・超臨界有機溶媒のNMRおよび分子動力学による研究）

近年、 超臨界領域の流体は、その 特異な溶解能や反応活性によって、溶媒として非常に注目されている。
超臨界流 体では気体に近い密度から 液体に近い密度まで連続的に 変化させて、そのマクロな物性を調べることが可 能である。液体状態から超 臨界状態に至る間の流体の構 造や動的性質を測定することにより、それぞれの 液体を構成する分子間の相 互作用についての、基本的な 知見をることができる。著者は有機溶媒としてべ ンゼン、塩化メチレン、ヘ キサフルオロベンゼン、アセ トニトリル、アセトン、メタノールを選んでいる 。これらの溶媒においては 、極性がそれぞれ異なってお り、とくにメタノールでは水素結合が存在するた め、ここで選んだ溶媒の構 造や動的性質を比較検討する ことは極めて興味深い。さらに著者はべンゼンー ヘキサフルオロベンゼン系 、重ベンゼンーメタノール系 などの二成分系についても実験を試みている。
著 者 の 行 っ た 研 究 は 、 核 磁 気 共 鳴 (NMR)の 測 定、 およ び 分子 動力 学(MD)計算 によ る、 ミク ロ なレ
ベ ルで の分 子 間の 相互 作用 の 解明 であ る。 NMRによ って、試料の密度、自己拡 散係数、スピンー格子緩和 時間 、化 学 シフ トな どの 測定を行っている。MDからは、自己拡散係数、動径 分布関数、相関時間などを計算し て実験と比較している。
本研 究に お ける 高温 ・高 圧 NMR測 定に 関 して は、 著者独自の工夫が各所に盛 り込まれている。下川らが開発し た一次元投影法による密度 測定を、現実の系に適用しその有用性を示すことができた。また従来、
横 緩 和 速 度 の 測 定 に 用 い ら れ て い る Carr‑Purcell‑Meiboom‑Gill（ CPMG） 法 が 、 超 臨 界 流 体 の 自
己 拡散 係数 の 測定 に極 めて 有用であることを初め て示した。このCPMG法の採用 により、測定時間が大幅に短縮さ れ、精度の良い自己拡散係 数の測定が可能になったので ある。
著者は 先ず、無極性溶媒であるべ ンゼンについて報告している 。一次元投影法で気ー液共存下の液体とー23―気体の 密度を 測定し 、臨界 指数の 値を決 めてい る。こ れらより 、高温 高圧下 におい ては、 NMRの一次元投 影法 に よ る 密度 測 定 が 極め て 有 用 であ る こ とを示 した。 またCPMG法 によっ て求め た自己 拡散係数 Dは 超臨 界 領 域 で大 き な 値 を示 し 、 密 度に 強 く 依 存し た 。 超 臨界 領 域 で はChapmanー Enskogの 理論 から予想 できる 値と良 く一致 した。 また液 体を含 む広い 温度範囲 でMDの 計算結果と実測値は一致した。プロトン の緩和 速度の 測定か らは、超臨界領域ではスピン回転相互作用が緩和機構として支配的であることを明か にした 。
次に 水素結 合をもつ 純メタノールに関する研究結果を報告している。ここではプロトンの化学シフトの温度・ 密度依 存性の 測定結 果を熟力学的モデルを用いて解析し、水素結合によるクラスターの形成について調ぺ ている 。その 結果、 超臨界領域においても二量体を主とするクラスターが存在することを明かにした。こ のこと を、別 に行っ たMD計算 からも 確かめ ている 。また メタノ ール分子の自己拡散係数の実測値とMD計 算の結 果を比較 し、液 体状態 では分 子の極 性すな わち水 素結合 が、拡散係数に強く影響することを 示し た 。 こ れに 対 し て 、超 臨 界 領 域で は 実 測 値は べ ン ゼ ンの 場 合 と 同様 に Chapmanー Enskogの 理論と良 く一致 し、分 子の極 性、二量体の形成が分子のダイナミックスには、事実上影響していないことを示した 。さら にMD計算 より超 臨界領 域では 分子の 拡散運 動が分 子間相 互作用ボテンシヤルの詳細には依存しな いこと を示し た。
その 他の溶 媒として ジクロロメタン、ヘキサフルオロペンゼン、アセトニト1Jル、アセトンの自己拡散係数と スピン ー格子 緩和時 間を測定している。その結果は、ベンゼンおよびメタノールで得られた結果と同 様の 傾 向 を 示し て い る 。す な わ ち 、超 臨 界 領 域の 自 己 拡 散係 数 の 実 測値 は 、 Chapmanー Enskogの理論式 から予 想され る値と 良くー致する。また臨界点近傍より、スピン回転相互作用が緩和の機構の中で支配的 となる 。さら にべン ゼンー ヘキサ フルオ ロベン ゼンニ成 分系の プロト ンとFー 19の自己 拡散係数とスピ ンー格 子緩和 時間、 および重ベンゼンーメタノール系のプロトンの化学シフトの測定結果を報告して いる 。

以 上を要 約する と、著 者は超 臨界領 域を含む 広い温 度範囲 で、い くっか の有機 溶媒の 分子運動をNMRと 分子動 力学シ ミュレーションから調ベ、分子性液体の分子間の相互作用について重要な知見を得て、この 分野の 研究に 大きく 貢献して いる。
よ っ て 、 著 者は 北 海 遭 大学 博 士 （ 理学 ） の 学 位を 受 け る のに 充 分 の 資格 が あ る もの と 認 め る。

― 24―