Opplestrup Tomas, Bulatov Vasily V, Gilmer George H, Kalos Malvin H, Sadigh Babak
Lawrence Livermore National Laboratory, University of California, Livermore, California 94551, USA.
Phys Rev Lett. 2006 Dec 8;97(23):230602. doi: 10.1103/PhysRevLett.97.230602. Epub 2006 Dec 4.
We present a novel Monte Carlo algorithm for N diffusing finite particles that react on collisions. Using the theory of first-passage processes and time dependent Green's functions, we break the difficult N-body problem into independent single- and two-body propagations circumventing numerous diffusion hops used in standard Monte Carlo simulations. The new algorithm is exact, extremely efficient, and applicable to many important physical situations in arbitrary integer dimensions.
我们提出了一种用于N个扩散有限粒子在碰撞时发生反应的新型蒙特卡罗算法。利用首达过程理论和含时格林函数,我们将困难的N体问题分解为独立的单体和两体传播,避免了标准蒙特卡罗模拟中使用的大量扩散跳跃。新算法精确、极其高效,适用于任意整数维中的许多重要物理情形。