Histogram-Free Reweighting with Grand Canonical Monte Carlo: Post-simulation Optimization of Non-bonded Potentials for Phase Equilibria
- Richard A. Messerly*
- Mohammad Soroush BarhaghiMohammad Soroush BarhaghiDepartment of Chemical Engineering and Materials Science, Wayne State University, Detroit, Michigan 48202, United StatesMore by Mohammad Soroush Barhaghi
- Jeffrey J. Potoff
- , and
- Michael R. Shirts
Histogram reweighting (HR) is a standard approach for converting grand canonical Monte Carlo (GCMC) simulation output into vapor–liquid coexistence properties (saturated liquid density, ρliqsat, saturated vapor density, ρvapsat, saturated vapor pressures, Pvapsat, and enthalpy of vaporization, ΔHv). We demonstrate that a histogram-free reweighting approach, namely, the Multistate Bennett Acceptance Ratio (MBAR), is similar to the traditional HR method for computing ρliqsat, ρvapsat, Pvapsat, and ΔHv. 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 parameterize 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 ρliqsat, ρvapsat, Pvapsat, and ΔHv in a limiting test case. Second, we utilize GCMC-MBAR to optimize an individualized (compound-specific) parameter (ψ) for 8 branched alkanes and 11 alkynes using the Mie Potentials for Phase Equilibria (MiPPE) force field. Third, we predict ρliqsat, ρvapsat, Pvapsat, and ΔHv 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 (ϵCH2, σCH2, and λCH2).
This article is cited by 4 publications.
- Marina P. Oliveira, Philippe H. Hünenberger. Influence of the Lennard-Jones Combination Rules on the Simulated Properties of Organic Liquids at Optimal Force-Field Parametrization. Journal of Chemical Theory and Computation 2023, 19
, 2048-2063. https://doi.org/10.1021/acs.jctc.2c01170
- Yan M. H. Gonçalves, Bruno A. C. Horta. gmak: A Parameter-Space Mapping Strategy for Force-Field Calibration. Journal of Chemical Theory and Computation 2023, 19
, 605-618. https://doi.org/10.1021/acs.jctc.2c00955
- Marina P. Oliveira, Maurice Andrey, Salomé R. Rieder, Leyla Kern, David F. Hahn, Sereina Riniker, Bruno A. C. Horta, Philippe H. Hünenberger. Systematic Optimization of a Fragment-Based Force Field against Experimental Pure-Liquid Properties Considering Large Compound Families: Application to Saturated Haloalkanes. Journal of Chemical Theory and Computation 2020, 16
, 7525-7555. https://doi.org/10.1021/acs.jctc.0c00683
- Marina P. Oliveira, Philippe H. Hünenberger. Systematic optimization of a fragment-based force field against experimental pure-liquid properties considering large compound families: application to oxygen and nitrogen compounds. Physical Chemistry Chemical Physics 2021, 23
, 17774-17793. https://doi.org/10.1039/D1CP02001C