Studies on free energy calculations. I. Thermodynamic integration using a polynomial path

Haluk Resat, Mihaly Mezei

Research output: Contribution to journalArticlepeer-review

86 Scopus citations

Abstract

The difference in the free energy of hydration between the C7 and αR conformations of alanine dipeptide at infinite dilution is computed using Monte Carlo method in a canonical ensemble with the united-atom AMBER force field. The recently introduced polynomial path with different exponents for the calculation of liquid state free energies using thermodynamic integration (TI) is tested. This is achieved by separating the interaction energy between the solute and the solvent molecules into the constituent r-12, r-6 and Coulomb terms. The separate contribution of each term is calculated and the comparison shows that the shape of the TI integrand as a function of the coupling parameter is mainly determined by the Coulomb term. Analysis of the convergence characteristics shows that five-point Gaussian quadrature integration would be sufficient to obtain quantitatively reliable results when a polynomial path is employed.

Original languageEnglish
Pages (from-to)6052-6061
Number of pages10
JournalJournal of Chemical Physics
Volume99
Issue number8
DOIs
StatePublished - 1993
Externally publishedYes

Fingerprint

Dive into the research topics of 'Studies on free energy calculations. I. Thermodynamic integration using a polynomial path'. Together they form a unique fingerprint.

Cite this