Schröder C, Steinhauser O
Department of Computational Biological Chemistry, University of Vienna, A-1090 Vienna, Austria.
J Chem Phys. 2008 Jun 14;128(22):224503. doi: 10.1063/1.2929848.
The vast majority of molecular dynamics simulations are based on nonpolarizable force fields with fixed partial charges for all atoms. The traditional way to obtain these charges are quantum-mechanical calculations performed prior to simulation. Unfortunately, the set of the partial charges heavily relies on the method and the basis set used. Therefore, investigations of the influence of charge variation on simulation data are necessary in order to validate various charge sets. This paper elucidates the consequences of different charge sets on the structure and dynamics of the ionic liquid: 1-ethyl-3-methyl-imidazolium dicyanoamide. The structural features seem to be more or less independent of the partial charge set pointing to a dominance of shape force as modeled by Lennard-Jones parameters. This can be seen in the radial distribution and orientational correlation functions. The role of electrostatic forces comes in when studying dynamical properties. Here, significant deviations between different charge sets can be observed. Overall, dynamics seems to be governed by viscosity. In fact, all dynamical parameters presented in this work can be converted from one charge set to another by viscosity scaling.
绝大多数分子动力学模拟基于对所有原子具有固定部分电荷的非极化力场。获得这些电荷的传统方法是在模拟之前进行量子力学计算。不幸的是,部分电荷集严重依赖于所使用的方法和基组。因此,为了验证各种电荷集,有必要研究电荷变化对模拟数据的影响。本文阐明了不同电荷集对离子液体1-乙基-3-甲基咪唑二氰胺的结构和动力学的影响。结构特征似乎或多或少与部分电荷集无关,这表明由 Lennard-Jones 参数建模的形状力占主导地位。这可以在径向分布函数和取向相关函数中看出。在研究动力学性质时,静电力的作用就显现出来了。在这里,可以观察到不同电荷集之间存在显著偏差。总体而言,动力学似乎受粘度支配。事实上,本文中呈现的所有动力学参数都可以通过粘度缩放从一个电荷集转换为另一个电荷集。