Improved diffusion monte carlo

Martin Hairer, Jonathan Weare

Research output: Contribution to journalArticlepeer-review

Abstract

We propose a modification, based on the RESTART (repetitive simulation trials after reaching thresholds) and DPR (dynamics probability redistribution) rare event simulation algorithms, of the standard diffusion Monte Carlo (DMC) algorithm. The new algorithm has a lower variance per workload, regardless of the regime considered. In particular, it makes it feasible to use DMC in situations where the "naïve" generalization of the standard algorithm would be impractical due to an exponential explosion of its variance. We numerically demonstrate the effectiveness of the new algorithm on a standard rare event simulation problem (probability of an unlikely transition in a Lennard-Jones cluster), as well as a high-frequency data assimilation problem.

Original languageEnglish (US)
Pages (from-to)1995-2021
Number of pages27
JournalCommunications on Pure and Applied Mathematics
Volume67
Issue number12
DOIs
StatePublished - Dec 1 2014

ASJC Scopus subject areas

  • General Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Improved diffusion monte carlo'. Together they form a unique fingerprint.

Cite this