Calculation of free energy landscapes: a histogram reweighted metadynamics approach

J Comput Chem. 2011 Jul 30;32(10):2084-96. doi: 10.1002/jcc.21790. Epub 2011 Apr 15.

Abstract

We present an efficient method for the calculation of free energy landscapes. Our approach involves a history-dependent bias potential, which is evaluated on a grid. The corresponding free energy landscape is constructed via a histogram reweighting procedure a posteriori. Because of the presence of the bias potential, it can be also used to accelerate rare events. In addition, the calculated free energy landscape is not restricted to the actual choice of collective variables and can in principle be extended to auxiliary variables of interest without further numerical effort. The applicability is shown for several examples. We present numerical results for the alanine dipeptide and the Met-Enkephalin in explicit solution to illustrate our approach. Furthermore, we derive an empirical formula that allows the prediction of the computational cost for the ordinary metadynamics variant in comparison with our approach, which is validated by a dimensionless representation.

Keywords: free energy landscape; grid calculations; metadynamics; molecular dynamics; reweighting.