Fast and Accurate Multidimensional Free Energy Integration

Enhanced sampling and free energy calculation algorithms of the thermodynamic integration family (such as the adaptive biasing force (ABF) method) are not based on the direct computation of a free energy surface but rather of its gradient. Integrating the free energy surface is nontrivial in dimensi...

Full description

Saved in:
Bibliographic Details
Published in:Journal of chemical theory and computation 2021-11, Vol.17 (11), p.6789-6798
Main Author: Hénin, Jérôme
Format: Article
Language:English
Subjects:
Algorithms
Biological Physics
Boundary conditions
Computation
Convergence
Energy
Free energy
Physics
Poisson equation
Portability
Sampling methods
Simulation
Statistical Mechanics
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Enhanced sampling and free energy calculation algorithms of the thermodynamic integration family (such as the adaptive biasing force (ABF) method) are not based on the direct computation of a free energy surface but rather of its gradient. Integrating the free energy surface is nontrivial in dimensions higher than one. Here, the author introduces a flexible, portable implementation of a Poisson equation formalism to integrate free energy surfaces from estimated gradients in dimensions 2 and 3 using any combination of periodic and nonperiodic (Neumann) boundary conditions. The algorithm is implemented in portable C++ and provided as a standalone tool that can be used to integrate multidimensional gradient fields estimated on a grid using any algorithm, such as umbrella integration as a post-treatment of umbrella sampling simulations. It is also included in the implementation of ABF (and its extended-system variant eABF) in the Collective Variables Module, enabling the seamless computation of multidimensional free energy surfaces within ABF and eABF simulations. A Python-based analysis toolchain is provided to easily plot and analyze multidimensional ABF simulation results, including metrics to assess their convergence. The Poisson integration algorithm can also be used to perform Helmholtz decomposition of noisy gradient estimates on the fly, resulting in an efficient implementation of the projected ABF (pABF) method proposed by Leliévre and co-workers. In numerical tests, pABF is found to lead to faster convergence with respect to ABF in simple cases of low intrinsic dimension but seems detrimental to convergence in a more realistic case involving degenerate coordinates and hidden barriers due to slower exploration. This suggests that variance reduction schemes do not always yield convergence improvements when applied to enhanced sampling methods.
ISSN:1549-9618
1549-9626
DOI:10.1021/acs.jctc.1c00593