|
本研究では、過冷却液体、ガラス転移の動的密度凡関数方程式の数値計算を主な目的とした。
この方程式はそのままでは密度が負になり得るための数値的不安定性の困難がある。これを避けるには、われわれが開発したモンテカルロシミュレーションモデルへのマッピングが有力な武器になった。これまさに密度の時間相関関数の計算を行い、現在可能な比較的短時間では妥当な結果を得た。一方サンプル数を限って行った剛体球液体模型の長時間のシミュレーションでは以下の結果を得た。(1)液体密度が融点を超えると、多数の準安定状態が出現し、その間の状態の重なりを表わす量の分布が広がり始める。広がりは密度の増加と共に大きくなる。(2)ひとつの状態の時間変化を初期状態との重なりを追及することにより調べた。MCTのクロスオーバー密度を超えると1つの準安定状態から別のそれに移り変わる所謂ホッピングが観測された。このようにして動的密度汎関数方程式はMCTの結果を含むと同時に長時間での熱的活性化仮定を記述することが確かめられた。これらの業績は国際的にも一定の評価を受け、平成9、10年だけでもイタリー、スペイン、ドイツ、アメリカ、韓国、日本での国際会議で代表者が発表を行った。
|
