Rapid Computation of Thermodynamic Properties Over Multidimensional Nonbonded Parameter Spaces using Adaptive Multistate Reweighting

TL;DR

This paper introduces an efficient method combining multistate reweighting and basis functions to rapidly compute thermodynamic properties across extensive multidimensional parameter spaces, significantly reducing computational costs.

Contribution

It presents a novel adaptive sampling and linear basis function approach that enables fast, accurate thermodynamic property estimation without extensive simulations for each parameter set.

Findings

Reduced computational cost from thousands of CPU years to days

Accurate estimation of solvation free energies and thermodynamic properties

Effective detection of poor configuration space overlap and convergence

Abstract

We show how thermodynamic properties of molecular models can be computed over a large, multidimensional parameter space by combining multistate reweighting analysis with a linear basis function approach. This approach reduces the computational cost to estimate thermodynamic properties from molecular simulations for over 130,000 tested parameter combinations from over a thousand CPU years to tens of CPU days. This speed increase is achieved primarily by computing the potential energy as a linear combination of basis functions, computed from either modified simulation code or as the difference of energy between two reference states, which can be done without any simulation code modification. The thermodynamic properties are then estimated with the Multistate Bennett Acceptance Ratio (MBAR) as a function of multiple model parameters without the need to define a priori how the states are connected by a pathway. Instead, we adaptively sample a set of points in parameter space to create mutual configuration space overlap. The existence of regions of poor configuration space overlap are detected by analyzing the eigenvalues of the sampled states' overlap matrix. The configuration space overlap to sampled states is monitored alongside the mean and maximum uncertainty to determine convergence, as neither the uncertainty or the configuration space overlap alone is a sufficient metric of convergence. This adaptive sampling scheme is demonstrated by estimating with high precision the solvation free energies of charged particles of Lennard-Jones plus Coulomb functional form. We also compute entropy, enthalpy, and radial distribution functions of unsampled parameter combinations using only the data from these sampled states and use the free energies estimates to examine the deviation of simulations from the Born approximation to the solvation free energy.