TY - JOUR
T1 - Robust estimation of diffusion-optimized ensembles for enhanced sampling
AU - Tian, Pengfei
AU - Jónsson, Sigurður Ægir
AU - Ferkinghoff-Borg, Jesper
AU - Krivov, Sergei V .
AU - Lindorff-Larsen, Kresten
AU - Irbäck, Anders
AU - Boomsma, Wouter Krogh
PY - 2014/2/11
Y1 - 2014/2/11
N2 - The multicanonical, or flat-histogram, method is a common technique to improve the sampling efficiency of molecular simulations. The idea is that free-energy barriers in a simulation can be removed by simulating from a distribution where all values of a reaction coordinate are equally likely, and subsequently reweight the obtained statistics to recover the Boltzmann distribution at the temperature of interest. While this method has been successful in practice, the choice of a flat distribution is not necessarily optimal. Recently, it was proposed that additional performance gains could be obtained by taking the position-dependent diffusion coefficient into account, thus placing greater emphasis on regions diffusing slowly. Although some promising examples of applications of this approach exist, the practical usefulness of the method has been hindered by the difficulty in obtaining sufficiently accurate estimates of the diffusion coefficient. Here, we present a simple, yet robust solution to this problem. Compared to current state-of-the-art procedures, the new estimation method requires an order of magnitude fewer data to obtain reliable estimates, thus broadening the potential scope in which this technique can be applied in practice.
AB - The multicanonical, or flat-histogram, method is a common technique to improve the sampling efficiency of molecular simulations. The idea is that free-energy barriers in a simulation can be removed by simulating from a distribution where all values of a reaction coordinate are equally likely, and subsequently reweight the obtained statistics to recover the Boltzmann distribution at the temperature of interest. While this method has been successful in practice, the choice of a flat distribution is not necessarily optimal. Recently, it was proposed that additional performance gains could be obtained by taking the position-dependent diffusion coefficient into account, thus placing greater emphasis on regions diffusing slowly. Although some promising examples of applications of this approach exist, the practical usefulness of the method has been hindered by the difficulty in obtaining sufficiently accurate estimates of the diffusion coefficient. Here, we present a simple, yet robust solution to this problem. Compared to current state-of-the-art procedures, the new estimation method requires an order of magnitude fewer data to obtain reliable estimates, thus broadening the potential scope in which this technique can be applied in practice.
KW - Faculty of Science
U2 - 10.1021/ct400844x
DO - 10.1021/ct400844x
M3 - Journal article
SN - 1549-9618
VL - 10
SP - 543
EP - 553
JO - Journal of Chemical Theory and Computation
JF - Journal of Chemical Theory and Computation
IS - 2
ER -