Histogram-Free Reweighting with Grand Canonical Monte
Carlo: Post-simulation Optimization of Non-bonded Potentials for Phase
Equilibria
Version 6 2019-06-11, 19:43
Version 5 2019-06-10, 10:03
Version 4 2019-06-09, 12:29
Version 3 2019-06-08, 10:13
Version 2 2019-06-07, 20:44
Version 1 2019-04-15, 18:36
Posted on 2019-06-11 - 19:43
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).
CITE THIS COLLECTION
DataCite
3 Biotech
3D Printing in Medicine
3D Research
3D-Printed Materials and Systems
4OR
AAPG Bulletin
AAPS Open
AAPS PharmSciTech
Abhandlungen aus dem Mathematischen Seminar der Universität Hamburg
ABI Technik (German)
Academic Medicine
Academic Pediatrics
Academic Psychiatry
Academic Questions
Academy of Management Discoveries
Academy of Management Journal
Academy of Management Learning and Education
Academy of Management Perspectives
Academy of Management Proceedings
Academy of Management Review
Messerly, Richard A.; Barhaghi, Mohammad Soroush; Potoff, Jeffrey J.; Shirts, Michael R. (2019). Histogram-Free Reweighting with Grand Canonical Monte
Carlo: Post-simulation Optimization of Non-bonded Potentials for Phase
Equilibria. ACS Publications. Collection. https://doi.org/10.1021/acs.jced.8b01232
or
Select your citation style and then place your mouse over the citation text to select it.
SHARE
Usage metrics
Read the peer-reviewed publication
AUTHORS (4)
RM
Richard A. Messerly
MB
Mohammad Soroush Barhaghi
JP
Jeffrey J. Potoff
MS
Michael R. Shirts
KEYWORDS
Grand Canonical Monte CarloMiPPE non-bonded parametersforce field parameterGCMC-MBARhistogram-free reweighting approachGCMCHRMultistate Bennett Acceptance RatioPhase Equilibria Histogram reweightingcanonical Monte CarloΔ H vρ liqMBARP vappost-simulation optimization schemesimulating force field iλ CH 2force field jforce fieldρ vapGCMC-HRphase equilibria dataphase equilibria propertiesσ CH 2vapor