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.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
