The penalty method for random walks with uncertain energies
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 a 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.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- May 1999
- DOI:
- 10.1063/1.478034
- arXiv:
- arXiv:physics/9812035
- Bibcode:
- 1999JChPh.110.9812C
- Keywords:
-
- 05.40.-a;
- 02.50.-r;
- Fluctuation phenomena random processes noise and Brownian motion;
- Probability theory stochastic processes and statistics;
- Physics - Computational Physics;
- Physics - Chemical Physics
- E-Print:
- 17 pages, 7 figures, accepted to Journal of Chemical Physics