Accelerating the Convergence and Reducing the Variance of Path Integral Calculations of Quantum Mechanical Free Energies by Using Local Reference Potentials
| Autor: | Steven L. Mielke, Donald G. Truhlar |
|---|---|
| Rok vydání: | 2012 |
| Předmět: |
Mathematical optimization
Discretization Monte Carlo method Mathematical analysis Centroid Computer Science Applications Quadrature (mathematics) symbols.namesake Path integral formulation Convergence (routing) symbols Feynman diagram Variance reduction Physical and Theoretical Chemistry Mathematics |
| Zdroj: | Journal of Chemical Theory and Computation. 8:1589-1596 |
| ISSN: | 1549-9626 1549-9618 |
| DOI: | 10.1021/ct300098p |
| Popis: | 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 H2O2, a very anharmonic system where torsional motion is important, and H2O, a system more amenable to harmonic reference treatment. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|