Comparing generalized ensemble methods for sampling of systems with many degrees of freedom.
We compare two standard replica exchange methods using temperature and dielectric constant as the scaling variables for independent replicas against two new corresponding enhanced sampling methods based on non-equilibrium statistical cooling (temperature) or descreening (dielectric). We test the four methods on a rough 1D potential as well as for alanine dipeptide in water, for which their relatively small phase space allows for the ability to define quantitative convergence metrics. We show that both dielectric methods are inferior to the temperature enhanced sampling methods, and in turn show that temperature cool walking (TCW) systematically outperforms the standard temperature replica exchange (TREx) method. We extend our comparisons of the TCW and TREx methods to the 5 residue met-enkephalin peptide, in which we evaluate the Kullback-Leibler divergence metric to show that the rate of convergence between two independent trajectories is faster for TCW compared to TREx. Finally we apply the temperature methods to the 42 residue amyloid-β peptide in which we find non-negligible differences in the disordered ensemble using TCW compared to the standard TREx. All four methods have been made available as software through the OpenMM Omnia software consortium (http://www.omnia.md/).
Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA.
James
J
Lincoff
Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA.
Sukanya
S
Sasmal
Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA.
Teresa
T
Head-Gordon