Article thumbnail
Location of Repository

The Penalty Method for Random Walks with Uncertain Energies

By D. M. Ceperley and M. Dewing

Abstract

We generalize the Metropolis et al. random walk algorithm to the situation where the energy is noisy and can only be estimated. Two possible applications are for long range potentials and for mixed quantum-classical simulations. If the noise is normally distributed we are able to modify the acceptance probability by applying a penalty to the energy difference and thereby achieve exact sampling even with very strong noise. When one has to estimate the variance we have an approximate formula, good in the limit of large number of independent estimates. We argue that the penalty method is nearly optimal. We also adapt an existing method by Kennedy and Kuti and compare to the penalty method on a one dimensional double well.Comment: 17 pages, 7 figures, accepted to Journal of Chemical Physics; Corrected swap of Figures 2 and 3. Added Figure

Topics: Physics - Computational Physics, Physics - Chemical Physics
Publisher: 'AIP Publishing'
Year: 1999
DOI identifier: 10.1063/1.478034
OAI identifier: oai:arXiv.org:physics/9812035

Suggested articles


To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.