Multistate reweighting and configuration mapping together accelerate the efficiency of thermodynamic calculations as a function of molecular geometry by orders of magnitude

J Chem Phys. 2013 Apr 21;138(15):154108. doi: 10.1063/1.4801332.

Abstract

We present an approach to calculate free energy and other thermodynamic property differences between molecules which have very little or no overlap in configuration space, but where a one-to-one mapping between the molecule geometries exists. The approach combines multistate reweighting with remapping of phase space between simulated states. We apply this method to calculate the free energy differences between non-overlapping, truncated harmonic oscillators, the free energy, enthalpy, and entropy differences between different parameterizations of rigid water, and differences in free energy of solvation between dipoles of different lengths. Previously difficult or impossible problems become either trivially easy or are improved in efficiency by two to five orders of magnitude.

Publication types

  • Research Support, U.S. Gov't, Non-P.H.S.

MeSH terms

  • Algorithms
  • Computer Simulation
  • Models, Chemical
  • Periodicity
  • Thermodynamics*
  • Water / chemistry*

Substances

  • Water