Spatially Resolved Thermodynamic Integration: An Efficient Method To Compute Chemical Potentials of Dense Fluids
Spatially Resolved Thermodynamic Integration: An Efficient Method To Compute Chemical Potentials of Dense Fluids
复制标题
DOI：
10.1021/acs.jctc.8b00002
复制
发表时间：
2018-07-01
影响因子：
5.5
通讯作者：
Potestio, Raffaello
中科院分区：
文献类型：
--
作者：
Heidari, Maziar;Kremer, Kurt;Potestio, Raffaello
计算化学势的许多流行方法都依赖于将测试粒子插入目标系统。在液体和液体混合物的情况下，随着密度或浓度的增加，这一程序的难度增加，因此不可避免地要使用复杂的强化取样技术。在这项工作中，我们提出了一种替代策略，空间分辨热力学积分，或简称SPARTIAN。在这里，分子被描述为在模拟子区域中具有原子级分辨率，并且在更大的储层中作为理想气体粒子。所有的分子都可以在子域之间自由扩散，动态地调整它们的分辨率。为了在整个模拟箱中实施均匀的密度分布，计算、施加单分子外部势，并通过两种分辨率之间的化学势差来识别。由于储层被表示为理想的气浴，因此该差值正好等于目标系统的过量化学势。本方法超越了粒子插入方法的高密度/浓度限制，因为进入目标系统区域的理想气体分子自发地适应局部环境。理想气体表示法对模拟的计算成本的贡献可以忽略不计，从而允许人们以最小的费用利用大型储层。该方法已通过计算过量化学势纯Lennard-Jones液体和混合物，SPC和SPC/E液体水，和氯化钠水溶液进行了验证。报告的结果很好地再现了这些系统的文献数据。
Many popular methods for the calculation of chemical potentials rely on the insertion of test particles into the target system. In the case of liquids and liquid mixtures, this procedure increases in difficulty upon increasing density or concentration, and the use of sophisticated enhanced sampling techniques becomes inevitable. In this work, we propose an alternative strategy, spatially resolved thermodynamic integration, or SPARTIAN for short. Here, molecules are described with atomistic resolution in a simulation subregion and as ideal gas particles in a larger reservoir. All molecules are free to diffuse between subdomains adapting their resolution on the fly. To enforce a uniform density profile across the simulation box, a single-molecule external potential is computed, applied, and identified with the difference in chemical potential between the two resolutions. Since the reservoir is represented as an ideal gas bath, this difference exactly amounts to the excess chemical potential of the target system. The present approach surpasses the high density/concentration limitation of particle insertion methods because the ideal gas molecules entering the target system region spontaneously adapt to the local environment. The ideal gas representation contributes negligibly to the computational cost of the simulation, thus allowing one to make use of large reservoirs at minimal expenses. The method has been validated by computing excess chemical potentials for pure Lennard-Jones liquids and mixtures, SPC and SPC/E liquid water, and aqueous solutions of sodium chloride. The reported results well reproduce literature data for these systems.