KのNPT計算を行って平衡構造を得ます。emim-ntf2.logvmdで可視化して結果を確かめてみます。deform後の構造(周期境界で分子内の結合が切れているように見えます実際は切れていません)シミュレーションで得られた統計データをプロットします。-wで統計量を評価する範囲を指定します。エネルギーや密度が次のように得られます。0.0136792計算が終了すると最終ステップで得られた平衡構造のemim-ntf2_eq.dataが保存されてます。構造やダイナミクスの評価、粘性計算はemim-ntf2_eq.dataを読んで平衡構造からリスタートすることとします。平衡構造からリスタートして次のような動径分布関数と平均二乗変位を計算してみます。emim-ntf2-1.log動径分布関数をプロットします。プロットして、拡散係数を求めています。psのデータを使って、直線を求め拡散計算を計算します。悪いですが、次のように自己拡散計算が得られます。構造とダイナミクス評価と同様に、平衡構造からリスタートして粘性の計算を行います。nsの時間windowで相関関数を求め積分による粘性係数を計算します。統計が悪いですが粘性が次のように求まります。Pa·s講習会では計算時間が限られているため、統計の良い相関関数で粘性の計算ができません。ファイルを作成し、構造や物性の計算行います。する系のインプットファイルを作成し計算を行います。イルを作成し計算を行います。イルを作成してみます。(100LiCl-200NaCl-50NaBr.inと100LiCl-200NaCl-50NaBr.data)が生成されます。