The adaptive multilevel finite element solution of the Poisson-Boltzmann equation on massively parallel computers

By using new methods for the parallel solution of elliptic partial differential equations, the teraflops computing power of massively parallel computers can be leveraged to perform electrostatic calculations on large biological systems. This paper describes the adaptive multilevel finite element solution of the Poisson-Boltzmann equation for a microtubule on the NPACI Blue Horizon - a massively parallel IBM RS/6000 (R) SID with eight POWER3 SMP nodes. The microtubule system is 40 nm in length and 24 nm in diameter, consists of roughly 600000 atoms, and has a net charge of -1800 e. Poisson-Boltzmann calculations are performed for several processor configurations, and the algorithm used shows excellent parallel scaling.