期刊
CHEMICAL PHYSICS
卷 370, 期 1-3, 页码 294-305出版社
ELSEVIER
DOI: 10.1016/j.chemphys.2010.02.014
关键词
Isothermal-isobaric ensemble; Measure-preserving integrator; Liouville operator; Holonomic constraints; ROLL algorithm
资金
- NSF [CHE-0704036]
- DOE [FG05-08OR23334]
- Conacyt
The Liouville operator approach is employed to derive a new measure-preserving geometric integrator for molecular dynamics simulations in the isothermal-isobaric (NPT) ensemble. Recently, we introduced such a scheme for NPT simulations with isotropic cell fluctuations in the absence of holonomic constraints [M. E. Tuckerman et al., J. Phys. A 39 (2006) 5629]. Here, we extend this approach to include both fully flexible cell fluctuations and holonomic constraints via a new and simpler formulation of the ROLL algorithm of Martyna et al. [Martyna et al., Mol. Phys. 87 (1996) 1117]. The new algorithm improves on earlier schemes in that it possesses a simpler mathematical structure and rigorously preserves the phase space metric. The new algorithm is illustrated on two example systems, ice and liquid n-decane. (C) 2010 Elsevier B. V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据