DeMille Robert C, Molinero Valeria
Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, USA.
J Chem Phys. 2009 Jul 21;131(3):034107. doi: 10.1063/1.3170982.
A coarse-grained model of NaCl in water is presented where the ions are modeled without charge to avoid computationally challenging electrostatics. A monatomic model of water [V. Molinero and E. B. Moore, J. Phys. Chem. B 113, 4008 (2009)] is used as the basis for this coarse-grain approach. The ability of Na(+) to disrupt the native tetrahedral arrangement of water molecules, and of Cl(-) to integrate within this organization, is preserved in this mW-ion model through parametrization focused on water's solvation of these ions. This model successfully reproduces the structural effect of ions on water, referenced to observations from experiments and atomistic molecular dynamics simulations, while using extremely short-ranged potentials. Without Coulomb interactions the model replicates details of the ion-water structure such as distinguishing contact and solvent-separated ion pairs and the free energy barriers between them. The approach of mimicking ionic effects with short-ranged interactions results in performance gains of two orders of magnitude compared to Ewald methods. Explored over a broad range of salt concentration, the model reproduces the solvation structure and trends of diffusion relative to atomistic simulations and experimental results. The functional form of the mW-ion model can be parametrized to represent other electrolytes. With increased computational efficiency and reliable structural fidelity, this model promises to be an asset for accessing significantly longer simulation time scales with an explicit solvent in a coarse-grained system involving, for example, polyelectrolytes such as proteins, nucleic acids, and fuel-cell membranes.
本文提出了一种水中氯化钠的粗粒度模型,其中离子建模时不带电荷,以避免计算上具有挑战性的静电作用。水的单原子模型[V. 莫利内罗和E. B. 摩尔,《物理化学杂志B》113, 4008 (2009)]被用作这种粗粒度方法的基础。在这个mW - 离子模型中,通过专注于水对这些离子的溶剂化作用的参数化,保留了Na⁺破坏水分子天然四面体排列的能力以及Cl⁻融入这种结构的能力。该模型成功地再现了离子对水的结构效应,参考了实验和原子分子动力学模拟的观察结果,同时使用了极短程的势。没有库仑相互作用时,该模型复制了离子 - 水结构的细节,比如区分接触离子对和溶剂分隔离子对以及它们之间的自由能垒。与埃瓦尔德方法相比,用短程相互作用模拟离子效应的方法在性能上提高了两个数量级。在广泛的盐浓度范围内进行探索,该模型相对于原子模拟和实验结果再现了溶剂化结构和扩散趋势。mW - 离子模型的函数形式可以进行参数化以表示其他电解质。随着计算效率的提高和可靠的结构保真度,该模型有望成为一种有用的工具,用于在涉及例如蛋白质、核酸和燃料电池膜等聚电解质的粗粒度系统中,使用显式溶剂访问显著更长的模拟时间尺度。