水、オクタン、ベンゼン、エタノール、メタノールの300Kにおける比誘電率を全原子モデルのMD計算により評価しました。水はSPC-FWモデルを使用し、その他の分子はJ-OCTAでのモデリングを行いました（力場はGAFFを用い、電荷はMO計算により求めた値を付与しました）。それぞれの分子で単一組成の系をそれぞれ作成しました。なお、粒子数が700個程度となるように分子数を決定し、境界は周期境界条件としました。 MD計算には弊社ソルバーである「VSOP」を使用しました。

  NPTアンサンブルによる300KのMD計算を実施して緩和状態とした後、さらに計算を実施して誘電率評価用の時系列データを取得しました。誘電率の算出には以下の式を用いています^[1]。

  MDにより推算された比誘電率の実験値との比較を図1および表1に示します。シミュレーションの結果は実験結果に近い値となりました。ただしオクタンおよびベンゼンのように分子内の電子の偏りが小さい物質では比誘電率が低く評価される傾向が見られます。