Journal of Chemical and Engineering Data, Vol.64, No.9, 3701-3717, 2019
Histogram-Free Reweighting with Grand Canonical Monte Carlo: Post-simulation Optimization of Non-bonded Potentials for Phase Equilibria
Histogram reweighting (HR) is a standard approach for converting grand canonical Monte Carlo (GCMG) simulation output into vapor-liquid coexistence properties (saturated liquid density, rho(sat)(liq), saturated vapor density, rho(sat)(vap), saturated vapor pressures, P-vap(sat), and enthalpy of vaporization, Delta H-v). We demonstrate that a histogram-free reweighting approach, namely, the Multistate Bennett Acceptance Ratio (MBAR), is similar to the traditional HR method for computing, rho(sat)(liq), rho(sat)(vap), P-vap(sat), and Delta H-v. The primary advantage of MBAR is the ability to predict phase equilibria properties for an arbitrary force field parameter set that has not been simulated directly. Thus, MBAR can greatly reduce the number of GCMC simulations that are required to parameter ize a force field with phase equilibria data. Four different applications of GCMC-MBAR are presented in this study. First, we validate that GCMC-MBAR and GCMC-HR yield statistically indistinguishable results for rho(sat)(liq), rho(sat)(vap), P-vap(sat), and Delta H-v in a limiting test case. Second, we utilize GCMC-MBAR to optimize an individualized (compound-specific) parameter (psi) for 8 branched alkanes and 11 alkynes using the Mie Potentials for Phase Equilibria (MiPPE) force field. Third, we predict rho(sat)(liq), rho(sat)(vap), P-vap(sat), and Delta H-v for force field j by simulating force field i, where i and j are common force fields from the literature. In addition, we provide guidelines for determining the reliability of GCMC-MBAR predicted values. Fourth, we develop and apply a post-simulation optimization scheme to obtain new MiPPE non-bonded parameters for cyclohexane (epsilon(CH2), sigma(CH2), and lambda(CH2)).