分子の充填によるモデリングと液体の分子動力学シミュレーション#

分子の充填機能#

Advance/NanoLaboにはユニットセルの空隙に任意の分子を充填する機能があります。この機能を用いることで、固液界面や溶液のモデルを容易に作成することができます。

分子の充填機能を用いたナノシミュレーション向けのモデリング例を3例ほど示します。1例目は水-メタノール混合系、2例目はアナターゼ二酸化チタン-水界面、3例目はグラフェン-水界面です。あらかじめ作成したスラブモデルなどに対して溶媒分子種とその濃度を指定するだけで、溶媒を含んだモデルを自動で作成することができます。

Tips

真空のユニットセルに分子を充填することで、気相または液相のみを含むモデル（本事例では水-メタノール混合系）を作成することができます。

溶媒の存在は反応機構や物性に大きな影響を与えるため、溶媒分子を考慮したモデルを用いることで、現実系により近い精緻なシミュレーションを行うことが可能になります。


(例1) 水-メタノール混合系	(例2) アナターゼ二酸化チタン-水界面	(例3) グラフェン-水界面

液体の分子動力学シミュレーション#

本事例では分子の充填機能を用いて水・メタノール・水-メタノール混合系（体積分率1:1）のモデルをそれぞれ作成し、常温常圧（300K・1bar）の条件下でNPTアンサンブルによる分子動力学シミュレーションを行いました。力場にはOPLS-AA¹を採用しています。シミュレーション時間は100psです。

以下に3つの系（水・メタノール・水-メタノール混合系）のモデル・温度変化・密度変化・アニメーションを示します。

水#


水のモデル	温度変化	密度変化

メタノール#


メタノールのモデル	温度変化	密度変化

水-メタノール混合系#


水-メタノール混合系のモデル	温度変化	密度変化

解析結果#

3つの系（水、メタノール、水-メタノール混合系）いずれにおいても、温度は300Kに収束しています。また、いずれの系においても、水分子とメタノール分子がユニットセルの中を激しく動き回っている様子が動画より分かります。

各系の密度を下表に示します。水およびメタノールの密度は実験値とよい一致を示しています。また、水-メタノール混合系の密度は0.89g/cm³となっており、各1成分系（水・メタノール）の中間的な値をとっています。

系	密度(計算値)² (g/cm³)	密度(実験値) (g/cm³)
水	1.00	0.9950 (25℃) ³
メタノール	0.75	0.7866 (25℃) ³
水-メタノール混合系	0.89	0.9272 (20℃) ⁴⁵