On the Test Particle Monte-Carlo method to solve the steady state Boltzmann equation, the congruity of its results with experiments and its potential for shared memory parallelism

2021 ◽  
pp. 110590
Author(s):  
Maxime Rondeau ◽  
R. Arès
2019 ◽  
Vol 25 (4) ◽  
pp. 329-340 ◽  
Author(s):  
Preston Hamlin ◽  
W. John Thrasher ◽  
Walid Keyrouz ◽  
Michael Mascagni

Abstract One method of computing the electrostatic energy of a biomolecule in a solution uses a continuum representation of the solution via the Poisson–Boltzmann equation. This can be solved in many ways, and we consider a Monte Carlo method of our design that combines the Walk-on-Spheres and Walk-on-Subdomains algorithms. In the course of examining the Monte Carlo implementation of this method, an issue was discovered in the Walk-on-Subdomains portion of the algorithm which caused the algorithm to sometimes take an abnormally long time to complete. As the problem occurs when a walker repeatedly oscillates between two subdomains, it is something that could cause a large increase in runtime for any method that used a similar algorithm. This issue is described in detail and a potential solution is examined.


Sign in / Sign up

Export Citation Format

Share Document