volume | PIER Journals

2005-11-01

A Parallelized 3D Floating Random-Walk Algorithm for the Solution of the Nonlinear Poisson-Boltzmann Equation

Kausik Chatterjee,

Dr. Kausik Chatterjee

Department of Electrical & Computer Engineering

California State University, Fresno

USA

Homepage

Jonathan Poggie

Dr. Jonathan Poggie

Computational Sciences Center

Air Force Research Laboratory

USA

Homepage

Progress In Electromagnetics Research, Vol. 57, 237-252, 2006

doi:10.2528/PIER05072802

Abstract

This paper presents a new three-dimensional floating random-walk (FRW) algorithm for the solution of the Nonlinear Poisson-Boltzmann (NPB) equation. The FRW method has not been previously used in the numerical solution of the NPB equation (and other nonlinear equations) because of the non-availability of analytical expressions for volumetric Green's functions. In the past, numerical studies using the FRW method have examined only the linearized Poisson-Boltzmann equation, producing solutions that are only accurate for small values of the potential. No such linearization is required for this algorithm. An approximate expression for a volumetric Green's functions has been calculated with the help of a novel use of perturbation theory, and the resultant integral form has been incorporated within the FRW framework. The algorithm requires no discretization of either the volume or the surface of the problem domains, and hence the memory requirements are expected to be lower than approaches based on spatial discretization, such as finite-difference methods. Another advantage of this algorithm is that each random walk is independent, so that the computational procedure is inherently parallelizable and an almost linear increase in computational speed is expected with increase in the number of processors. We have recently published the preliminary results for benchmark problems in one and two dimensions. In this work, we present our results for benchmark problems in three dimensions and demonstrate excellent agreement between the FRW- and finite-difference based algorithms. We also present the results of parallelization of the newly developed FRW algorithm. The solution of the NPB equation has applications in diverse branches of science and engineering including (but not limited to) the modeling of plasma discharges, semiconductor device modeling and the modeling of biomolecular structures and dynamics.

Citation

Copy Export

Kausik Chatterjee, and Jonathan Poggie, "A Parallelized 3D Floating Random-Walk Algorithm for the Solution of the Nonlinear Poisson-Boltzmann Equation," Progress In Electromagnetics Research, Vol. 57, 237-252, 2006.
doi:10.2528/PIER05072802

References

1. Mitchner, M. and C. H. Kruger, Partial ly Ionized Gases, John Wiley & Sons, 1973.

2. Sze, S. M., Physics of Semiconductor Devices, 366-368, 2nd ed., 1999.

3. Davis, M. E. and J. A. McCammon, "Electrostatics in bio-molecular structure and dynamics," Chemical Reviews, Vol. 90, 509-521, 1990.
doi:10.1021/cr00101a005

4. Brown, G. M., Modern Mathematics for Engineers, E. F. Beckenbach (ed.), 1956.

5. Le Coz, Y. L., H. J. Greub, and R. B. Iverson, "Performance of random walk capacitance extractors for IC interconnects: a numerical study," Solid-State Electronics, Vol. 42, 581-588, 1998.
doi:10.1016/S0038-1101(97)00283-9

6. Le Coz, Y. L., R. B. Iverson, T. L. Sham, H. F. Tiersten, and M. S. Shepard, "Theory of a floating random walk algorithm for solving the steady-state heat equation in complex materially inhomogeneous rectilinear domains," Numerical Heat Transfer, Vol. 26, 353-366, 1994.

7. Haberman, R., Elementary Applied Partial Differential Equations, 3rd ed., 1998.

8. Hwang, C.-O. and M. Mascagni, "Efficient modified 'walk on spheres' algorithm for the linearized Poisson-Boltzmann equation," App. Phys. Lett., Vol. 78, No. 6, 787-789, 2001.
doi:10.1063/1.1345817

9. Mascagni, M. and N. A. Simonov, "Monte carlo methods for calculating the electrostatic energy of a molecule," Proceedings of the 2003 International Conference on Computational Science (ICCS 2003)..

10. Chatterjee, K. and J. Poggie, "A meshless stochastic algorithm for the solution of the nonlinear Poisson-Boltzmann equation in the context of plasma discharge modeling: 1D analytical benchmark," Proceedings of the 17th AIAA Computational Fluid Dynamics Conference, 6-9.

11. Chatterjee, K. and J. Poggie, "A two-dimensional floating random-walk algorithm for the solution of the nonlinear Poisson-Boltzmann equation: application to the modeling of plasma sheaths," Proceedings of the 3rd MIT Conference on Computational Fluid and Solid Mechanics, 14-17, 2005.

12. Sobol, I. M., A Primer for the Monte Carlo Method, CRC Press, 1994.

13. Curtiss, J. H., "Monte carlo methods for the iteration of linear operators," J. Math. and Phys., Vol. 32, 209-232, 1954.

14. Sabelfeld, K. K. and N. A. Simonov, Random Walks on Boundary for Solving PDEs, 119-120, 119-120, 1994.

15. Chatterjee, K., "Development of a floating random walk algorithm for solving Maxwell's equations in complex IC-interconnect structures," Rensselaer Polytechnic Institute, 48106-1346.

16. Chatterjee, K. and Y. L. Le Coz, "A floating random-walk algorithm based on iterative perturbation theory: solution of the 2D vector-potential maxwell-helmholtz equation," Applied Computational Electromagnetics Society Journal, Vol. 18, No. 1, 48-57, 2003.

17. Hammersley, J. M. and D. C. Handscomb, Monte Carlo Methods, 50-54, 50-54, 1964.