Rahbari Ahmadreza, Hens Remco, Moultos Othonas A, Dubbeldam David, Vlugt Thijs J H
Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands.
Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands.
J Chem Theory Comput. 2020 Mar 10;16(3):1757-1767. doi: 10.1021/acs.jctc.9b01097. Epub 2020 Feb 12.
We introduce an alternative method to perform free energy calculations for mixtures at multiple temperatures and pressures from a single simulation, by combining umbrella sampling and the continuous fractional component Monte Carlo method. One can perform a simulation of a mixture at a certain pressure and temperature and accurately compute the chemical potential at other pressures and temperatures close to the simulation conditions. This method has the following advantages: (1) Accurate estimates of the chemical potential as a function of pressure and temperature are obtained from a single state simulation without additional postprocessing. This can potentially reduce the number of simulations of a system for free energy calculations for a specific temperature and/or pressure range. (2) Partial molar volumes and enthalpies are obtained directly from the estimated chemical potentials. We tested our method for a Lennard-Jones system, aqueous mixtures of methanol at = 298 K and = 1 bar, and a mixture of ammonia, nitrogen, and hydrogen at = 573 K and = 800 bar. For pure methanol ( = 410 molecules), we observed that the estimated chemical potentials from umbrella sampling are in excellent agreement with the reference values obtained from independent simulations, for Δ = ±15 K and Δ = 100 bar (with respect to the simulated system). For larger systems, this range becomes smaller because the relative fluctuations of energy and volume become smaller. Without sufficient overlap, the performance of the method will become poor especially for nonlinear variations of the chemical potential.
我们引入一种替代方法,通过结合伞形采样和连续分数组分蒙特卡罗方法,从单个模拟中对多种温度和压力下的混合物进行自由能计算。人们可以在特定压力和温度下对混合物进行模拟,并准确计算接近模拟条件的其他压力和温度下的化学势。该方法具有以下优点:(1)从单个状态模拟中可获得作为压力和温度函数的化学势的准确估计值,无需额外的后处理。这有可能减少针对特定温度和/或压力范围进行自由能计算时系统的模拟次数。(2)偏摩尔体积和焓可直接从估计的化学势中获得。我们对一个 Lennard-Jones 系统、298 K 和 1 bar 下甲醇的水混合物以及 573 K 和 800 bar 下氨、氮和氢的混合物测试了我们的方法。对于纯甲醇((N = 410) 个分子),我们观察到,对于 (\Delta T = \pm15) K 和 (\Delta P = 100) bar(相对于模拟系统),伞形采样估计的化学势与独立模拟获得的参考值高度吻合。对于更大的系统,这个范围会变小,因为能量和体积的相对波动变小。如果没有足够的重叠,该方法的性能会变差,特别是对于化学势的非线性变化。