Molecular Dynamics by EXCEL

Abstract

EXCELによる分子動力学プログラムを開発した．採用した分子模型はレナードージョーンズ関数である．実行方法，結果の読み方，設定条件の入力の仕方などを示した．ワークシートから系の状態を規定する温度と数密度を与える．熱力学量の他にニ体相関関数，積算配位数，平均ニ乗変位，速度自己相関関数，自己拡散係数などが計算できる．プログラムはEXCEL worksheet として付録に示される．

Figures

Figure 1.

 The coordinate in the periodic boundary condition X (1)/σ vs. step.

Figure 2.

 The coordinate Xr (1)/σ vs. step.

Figure 3.

 The pair correlation function gr (r) and running coordination number rcn (r) vs. molecular distance r/σ.

Figure 4.

 The mean square displacement msd (t)/σ² vs. time t/τ.

Figure 5.

 The velocity autocorrelation function c (t)/(σ²/τ²) and accumulation number s (t) vs. time t/τ.

References

• [1]   I. Okada, E. Osawa, (ed),”Bunshi Simyureshon Nyumon”, Kaibundo (1989)
• [2]   P. W. Atkins, Physical Chemistry, Oxford Univ. Press, Oxford (1998).
• [3]   F. Cuadros, I. Cachadiña, W. Ahamuda, Mol. Eng., 6, 319 (1996).
• [4]   https://www.jstage.jst.go.jp/browse/jccj/-char/ja/
• [5]   M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids, Clarendon Press, Oxford (1992).
• [6]   Y. Kataoka, Y. Yamada, J. Comput. Chem. Jpn., 13, 130 (2014).
• [7]   K. Terasaka,”Kagakukei Gakuseinotameno Excel/VBA Nyumon”, Koronasha (2009).
• [8]   Y. Kataoka, Y. Yamada, Bull. Res. Cent. Comp. Multimedia Stud. Hosei Univ., (2015), (to be published).
• [9]   Y. Kataoka,”Bunshidorikigaku To Monte Carlo Hou", Koudansha, (1994).