分子動力学計算による密度ゆらぎの分子論とその実験的検証

Research Project

Project/Area Number	06214206
Research Category	Grant-in-Aid for Scientific Research on Priority Areas
Allocation Type	Single-year Grants
Research Institution	Tokyo Institute of Technology
Principal Investigator	岡崎進東京工業大学, 大学院総合理工学研究所, 助手 (70194339)
Co-Investigator(Kenkyū-buntansha)	岡田勲東京工業大学, 大学院総合理工学研究所, 教授 (60011582)
Project Period (FY)	1994
Project Status	Completed (Fiscal Year 1994)
Budget Amount *help	¥1,700,000 (Direct Cost: ¥1,700,000) Fiscal Year 1994: ¥1,700,000 (Direct Cost: ¥1,700,000)
Keywords	超臨界流体 / 密度ゆらぎ / 中性子散乱 / 分子動力学シミュレーション
Research Abstract	PVT曲線や溶質の溶解度等の諸物性と密接な関係を有していると考えられ、超臨界流体を特徴づける最も重要な物理的性質のひとつである臨界点近傍における大きな密度ゆらぎの現象の解明に焦点を絞り、分子動力学計算と中性子小角散乱実験の双方から系統的な研究を進め、これまで全く未知であったその密度ゆらぎに対するミクロスコピックな分子論的描像に関するいくつかの新たな知見を得た。用いた試料は、二酸化炭素及びキセノンである。高エネ研の中性子源を用い、中性子小角散乱実験はWINK分光器を用い、また広角散乱実験はHIT分光器を用いて行った。また、密度ゆらぎの検討のためのMD計算は、超臨界キセノン分子約2万個の系に対し、LJポテンシャルでの臨界点より約15K高い温度で行った。中性子散乱実験に関しては、二つの分光器を合わせ用いることにより、0.018A^<-1>から30A^<-1>まで3桁半に及ぶ広いQ領域での構造因子S(Q)の測定に成功した。ここにおいて、密度ゆらぎの大きな超臨界流体の構造に関する測定手法のひとつが確立されたということができる。この構造因子S(Q)のフーリエ変換として、あるいはMD計算からの軌跡に基づいた直接計算から求めた動径分布関数g(r)は、いずれも1からの系統的偏りを示し、ここにおいて密度ゆらぎを反映した実験量を示し得たということになる。特に、MD計算からの分子の投影図の直接観察、またボロノイ多面体解析、またその他の様々な統計力学関数から流体中における分子クラスターの生成が確認され、さらにはその大きさの分布を構成分子数や占有体積の分布により定量化した。そこでは、分布関数の関数形が、フラクショナルな指数関数となるということを見いだし、前年度までに行った、確率論モデルにおける取り扱いとの良好な対応を示した。