Infinite Parallel Plates Algorithm. Issue 6 (5th May 2020)

• Record Type: Journal Article
• Authors: Kim, Mingyo; Hwang, Chi‐Ok
• Is Part Of: Advanced theory and simulations. Volume 3:Issue 6(2020)
• Page Start: n/a
• Page End: n/a
• Publication Date: 2020-05-05
• Subjects: Monte Carlo simulation -- parallel plates -- walks‐on‐spheres algorithm; Science -- Simulation methods -- Periodicals; Science -- Methodology -- Periodicals; Engineering -- Simulation methods -- Periodicals; Engineering -- Methodology -- Periodicals; 507.21
• DOI: 10.1002/adts.202000014
• Languages: English
• ISSNs: 2513-0390

Abstract

Abstract: 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. Abstract : A new parallel‐plates algorithm without a diffusion absorption layer is developed. 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. It is shown that an exact random walk sampling on the parallel planes is possible.