Accelerating the convergence and reducing the variance of path integral calculations of quantum mechanical free energies by using local reference Potentials

Steven L. Mielke, Donald G. Truhlar

Research output: Contribution to journalArticlepeer-review

6 Scopus citations

Abstract

We present two new methods to accelerate the convergence of Feynman path integral calculations of thermodynamic partition functions. The first enhancement uses information from instantaneous normal mode (INM) calculations to decrease the number of discretized points necessary to represent the Feynman paths and is denoted the local generalized Pitzer-Gwinn (LGPG) scheme. The second enhancement, denoted harmonically guided variance reduction (HGVR), reduces the variance in Monte Carlo (MC) calculations by exploiting the correlation between the sampling error associated with the sum over paths at a particular centroid location for the accurate potential and for the INM approximation of a model potential, the latter of which can be exactly calculated. The LGPG scheme can reduce the number of quadrature points required along the paths by nearly an order of magnitude, and the HGVR scheme can reduce the number of MC samples needed to achieve a target accuracy by more than an order of magnitude. Numerical calculations are presented for H 2O 2, a very anharmonic system where torsional motion is important, and H 2O, a system more amenable to harmonic reference treatment.

Original languageEnglish (US)
Pages (from-to)1589-1596
Number of pages8
JournalJournal of Chemical Theory and Computation
Volume8
Issue number5
DOIs
StatePublished - May 8 2012

Fingerprint

Dive into the research topics of 'Accelerating the convergence and reducing the variance of path integral calculations of quantum mechanical free energies by using local reference Potentials'. Together they form a unique fingerprint.

Cite this