Premium
Infinite Parallel Plates Algorithm
Author(s) -
Kim Mingyo,
Hwang ChiOk
Publication year - 2020
Publication title -
advanced theory and simulations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.068
H-Index - 17
ISSN - 2513-0390
DOI - 10.1002/adts.202000014
Subject(s) - random walk , monte carlo method , diffusion , algorithm , isomorphism (crystallography) , jump , series (stratigraphy) , jump diffusion , monte carlo integration , sampling (signal processing) , markov chain monte carlo , mathematics , computer science , hybrid monte carlo , physics , paleontology , statistics , chemistry , quantum mechanics , biology , crystal structure , thermodynamics , crystallography , filter (signal processing) , computer vision
In a diffusion Monte Carlo simulation, when diffusion starts between two parallel planes, the “walks‐on‐spheres” (WOS) algorithm is usually used. In this case, an absorption ε‐layer is used to terminate the diffusion. Here, based on isomorphism between the electrostatic Poisson problem and the corresponding diffusion motion expectation of the first passage, a new parallel‐plates algorithm is developed without an absorption layer, which induces another diffusion error in addition to the intrinsic Monte Carlo error. As there is not an analytic closed solution for the induced charge density on the parallel plates by a charge at the center between them, a series solution is used, combined with the acceptance–rejection sampling method. Using this algorithm, even though a series solution is used, it is shown that an exact sampling can be performed, which means that the random walk jump can be performed to the parallel planes without using the detailed WOS random walk jumps. It is verified that the proposed parallel plates algorithm is significantly more efficient than the current WOS algorithm.