Free energies from dynamic weighted histogram analysis using unbiased Markov state model

E Rosta, G Hummer - Journal of chemical theory and computation, 2015 - ACS Publications
Journal of chemical theory and computation, 2015ACS Publications
The weighted histogram analysis method (WHAM) is widely used to obtain accurate free
energies from biased molecular simulations. However, WHAM free energies can exhibit
significant errors if some of the biasing windows are not fully equilibrated. To account for the
lack of full equilibration, we develop the dynamic histogram analysis method (DHAM).
DHAM uses a global Markov state model to obtain the free energy along the reaction
coordinate. A maximum likelihood estimate of the Markov transition matrix is constructed by …
The weighted histogram analysis method (WHAM) is widely used to obtain accurate free energies from biased molecular simulations. However, WHAM free energies can exhibit significant errors if some of the biasing windows are not fully equilibrated. To account for the lack of full equilibration, we develop the dynamic histogram analysis method (DHAM). DHAM uses a global Markov state model to obtain the free energy along the reaction coordinate. A maximum likelihood estimate of the Markov transition matrix is constructed by joint unbiasing of the transition counts from multiple umbrella-sampling simulations along discretized reaction coordinates. The free energy profile is the stationary distribution of the resulting Markov matrix. For this matrix, we derive an explicit approximation that does not require the usual iterative solution of WHAM. We apply DHAM to model systems, a chemical reaction in water treated using quantum-mechanics/molecular-mechanics (QM/MM) simulations, and the Na+ ion passage through the membrane-embedded ion channel GLIC. We find that DHAM gives accurate free energies even in cases where WHAM fails. In addition, DHAM provides kinetic information, which we here use to assess the extent of convergence in each of the simulation windows. DHAM may also prove useful in the construction of Markov state models from biased simulations in phase-space regions with otherwise low population.
ACS Publications
以上显示的是最相近的搜索结果。 查看全部搜索结果