skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A whole-path importance-sampling scheme for Feynman path integral calculations of absolute partition functions and free energies

Journal Article · · Journal of Chemical Physics
DOI:https://doi.org/10.1063/1.4939869· OSTI ID:1469169

Using Feynman path integrals, a molecular partition function can be written as a double integral with the inner integral involving all closed paths centered at a given molecular configuration, and the outer integral involving all possible molecular configurations. In previous work employing Monte Carlo methods to evaluate such partition functions, we presented schemes for importance sampling and stratification in the molecular configurations that constitute the path centroids, but we relied on free-particle paths for sampling the path integrals. At low temperatures, the path sampling is expensive because the paths can travel far from the centroid configuration. In this paper, we now present a scheme for importance sampling of whole Feynman paths based on harmonic information from an instantaneous normal mode calculation at the centroid configuration, which we refer to as harmonically guided whole-path importance sampling (WPIS). We obtain paths conforming to our chosen importance function by rejection sampling from a distribution of free-particle paths. Sample calculations on CH4 demonstrate that at a temperature of 200 K, about 99.9% of the free-particle paths can be rejected without integration, and at 300 K, about 98% can be rejected. Finally, we also show that it is typically possible to reduce the overhead associated with the WPIS scheme by sampling the paths using a significantly lower-order path discretization than that which is needed to converge the partition function.

Research Organization:
Univ. of Minnesota, Minneapolis, MN (United States)
Sponsoring Organization:
USDOE Office of Science (SC), Basic Energy Sciences (BES)
Grant/Contract Number:
FG02-86ER13579
OSTI ID:
1469169
Alternate ID(s):
OSTI ID: 1235492
Journal Information:
Journal of Chemical Physics, Vol. 144, Issue 3; ISSN 0021-9606
Publisher:
American Institute of Physics (AIP)Copyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 1 work
Citation information provided by
Web of Science

References (22)

Set of Co-Ordinate Systems Which Diagonalize the Kinetic Energy of Relative Motion journal February 1959
Erratum: “Calculation of converged rovibrational energies and partition function for methane using vibrational–rotational configuration interaction” [J. Chem. Phys. 121, 2071 (2004)] journal September 2011
Improved methods for Feynman path integral calculations and their application to calculate converged vibrational–rotational partition functions, free energies, enthalpies, entropies, and heat capacities for methane journal January 2015
Algorithm 806: SPRNG: a scalable library for pseudorandom number generation journal September 2000
Improved Methods for Feynman Path Integral Calculations of Vibrational−Rotational Free Energies and Application to Isotopic Fractionation of Hydrated Chloride Ions journal April 2009
A simple model for the treatment of imaginary frequencies in chemical reaction rates and molecular liquids journal August 2009
Displaced-points path integral method for including quantum effects in the Monte Carlo evaluation of free energies journal July 2001
Path integrals in the theory of condensed helium journal April 1995
Classical trajectory studies of the reaction CH 4 +H→CH 3 +H 2 journal April 1995
Approximate path integral methods for partition functions journal March 1993
Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 journal April 1996
Convenient and accurate discretized path integral methods for equilibrium quantum mechanical calculations journal August 1981
Calculation of converged rovibrational energies and partition function for methane using vibrational–rotational configuration interaction journal August 2004
Equation of State Calculations by Fast Computing Machines journal June 1953
Accelerating the Convergence and Reducing the Variance of Path Integral Calculations of Quantum Mechanical Free Energies by Using Local Reference Potentials journal April 2012
Efficient methods for including quantum effects in Monte Carlo calculations of large systems: Extension of the displaced points path integral method and other effective potential methods to calculate properties and distributions journal January 2013
The Ziggurat Method for Generating Random Variables journal January 2000
Monte Carlo Statistical Methods book January 2004
Monte Carlo sampling methods using Markov chains and their applications journal April 1970
Instantaneous normal mode analysis as a probe of cluster dynamics journal July 1990
Accurate vibrational-rotational partition functions and standard-state free energy values for H2O2 from Monte Carlo path-integral calculations journal September 2004
A ‘path-by-path’ monotone extrapolation sequence for Feynman path integral calculations of quantum mechanical free energies journal September 2003