The Finite Element Approximation of the Nonlinear Poisson-Boltzmann Equation

A widely used electrostatics model in the biomolecular modeling community, the nonlinear Poisson-Boltzmann equation, along with its finite element approximation, are analyzed in this paper. A regularized Poisson-Boltzmann equation is introduced as an auxiliary problem, making it possible to study the original nonlinear equation with delta distribution sources. A priori error estimates for the finite element approximation are obtained for the regularized Poisson-Boltzmann equation based on certain quasi-uniform grids in two and three dimensions. Adaptive finite element approximation through local refinement driven by an a posteriori error estimate is shown to converge. The Poisson-Boltzmann equation does not appear to have been previously studied in detail theoretically, and it is hoped that this paper will help provide molecular modelers with a better foundation for their analytical and computational work with the Poisson-Boltzmann equation. Note that this article apparently gives the first rigorous convergence result for a numerical discretization technique for the nonlinear Poisson-Boltzmann equation with delta distribution sources, and it also introduces the first provably convergent adaptive method for the equation. This last result is currently one of only a handful of existing convergence results of this type for nonlinear problems.

[1]  Long Chen,et al.  Convergence and optimality of adaptive mixed finite element methods , 2010, Math. Comput..

[2]  Joseph W. Jerome,et al.  CONSISTENCY OF SEMICONDUCTOR MODELING: AN EXISTENCE/STABILITY ANALYSIS FOR THE STATIONARY VAN ROOSBROECK SYSTEM* , 1985 .

[3]  R. Bank,et al.  Some Refinement Algorithms And Data Structures For Regular Local Mesh Refinement , 1983 .

[4]  Michael J. Holst,et al.  The adaptive multilevel finite element solution of the Poisson-Boltzmann equation on massively parallel computers , 2001, IBM J. Res. Dev..

[5]  Giuseppe Savaré,et al.  Regularity Results for Elliptic Equations in Lipschitz Domains , 1998 .

[6]  Shen Yifan,et al.  A mortar finite element approximation for the linear Poisson-Boltzmann equation , 2005 .

[7]  Michael J. Holst,et al.  Numerical solution of the nonlinear Poisson–Boltzmann equation: Developing more robust and efficient methods , 1995, J. Comput. Chem..

[8]  Carsten Carstensen,et al.  Convergence analysis of an adaptive nonconforming finite element method , 2006, Numerische Mathematik.

[9]  R. Verfürth.  A posteriori error estimates for nonlinear problems: finite element discretizations of elliptic equations , 1994 .

[10]  J. Milovich,et al.  Solution of the nonlinear Poisson-Boltzmann equation using pseudo-transient continuation and the finite element method. , 2002, Journal of colloid and interface science.

[11]  M. Rivara.  Mesh Refinement Processes Based on the Generalized Bisection of Simplices , 1984 .

[12]  J. Zou,et al.  Finite element methods and their convergence for elliptic and parabolic interface problems , 1998 .

[13]  James H. Bramble,et al.  A finite element method for interface problems in domains with smooth boundaries and interfaces , 1996, Adv. Comput. Math..

[14]  William F. Mitchell,et al.  A comparison of adaptive refinement techniques for elliptic problems , 1989, TOMS.

[15]  E. H ckel,et al.  Zur Theorie der Elektrolyte , 1924 .

[16]  Jinchao Xu.  Two-grid Discretization Techniques for Linear and Nonlinear PDEs , 1996 .

[17]  Ricardo H. Nochetto,et al.  An Adaptive Uzawa FEM for the Stokes Problem: Convergence without the Inf-Sup Condition , 2002, SIAM J. Numer. Anal..

[18]  Jinchao Xu,et al.  Iterative Methods by Space Decomposition and Subspace Correction , 1992, SIAM Rev..

[19]  W. Dörfler.  A convergent adaptive algorithm for Poisson's equation , 1996 .

[20]  Ludmil T. Zikatanov,et al.  A monotone finite element scheme for convection-diffusion equations , 1999, Math. Comput..

[21]  Michael Levitt,et al.  Finite‐difference solution of the Poisson–Boltzmann equation: Complete elimination of self‐energy , 1996, Journal of computational chemistry.

[22]  H. Orland,et al.  Steric Effects in Electrolytes: A Modified Poisson-Boltzmann Equation , 1997, cond-mat/9803258.

[23]  R. Skeel,et al.  Unified Multilevel Adaptive Finite Element Methods for Elliptic Problems , 1988 .

[24]  Michael J. Holst,et al.  An Odyssey into Local Refinement and Multilevel Preconditioning III: Implementation and Numerical Experiments , 2003, SIAM J. Sci. Comput..

[25]  James M. Briggs,et al.  Computation Unravels Mysteries of Molecular Biophysics , 1992 .

[26]  Michael J. Holst,et al.  Optimality of Multilevel Preconditioners for Local Mesh Refinement in Three Dimensions , 2006, SIAM J. Numer. Anal..

[27]  Rüdiger Verfürth,et al.  A posteriori error estimation and adaptive mesh-refinement techniques , 1994 .

[28]  Willy Dörfler,et al.  An adaptive finite element method for a linear elliptic equation with variable coefficients , 2000 .

[29]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[30]  Ricardo H. Nochetto,et al.  Convergence of Adaptive Finite Element Methods for General Second Order Linear Elliptic PDEs , 2005, SIAM J. Numer. Anal..

[31]  Michael Holst,et al.  The Poisson-Boltzmann equation: Analysis and multilevel numerical solution , 1994 .

[32]  Ivo Babuska,et al.  The finite element method for elliptic equations with discontinuous coefficients , 1970, Computing.

[33]  F. Bornemann,et al.  Adaptive multivlevel methods in three space dimensions , 1993 .

[34]  Wenbin Chen,et al.  A mortar finite element approximation for the linear Poisson-Boltzmann equation , 2005, Appl. Math. Comput..

[35]  Rob P. Stevenson,et al.  Optimality of a Standard Adaptive Finite Element Method , 2007, Found. Comput. Math..

[36]  E. G. Sewell,et al.  Automatic generation of triangulations for piecewise polynomial approximation , 1972 .

[37]  Rüdiger Verfürth,et al.  A posteriori error estimates for nonlinear problems , 1994 .

[38]  Gary Patterson,et al.  Physical Chemistry of Macromolecules , 2007 .

[39]  Jinchao Xu,et al.  The method of alternating projections and the method of subspace corrections in Hilbert space , 2002 .

[40]  Michael J. Holst,et al.  Adaptive multilevel finite element solution of the Poisson–Boltzmann equation I. Algorithms and examples , 2001 .

[41]  Carsten Carstensen.  Convergence of adaptive FEM for a class of degenerate convex minimization problems , 2007 .

[42]  Michael J. Holst,et al.  Adaptive multilevel finite element solution of the Poisson-Boltzmann equation II. Refinement at solvent-accessible surfaces in biomolecular systems , 2000, J. Comput. Chem..

[43]  Mariette Yvinec,et al.  Variational tetrahedral meshing , 2005, ACM Trans. Graph..

[44]  Rob P. Stevenson,et al.  An Optimal Adaptive Finite Element Method , 2004, SIAM J. Numer. Anal..

[45]  Nathan A. Baker,et al.  Adaptive multilevel finite element solution of the Poisson–Boltzmann equation I. Algorithms and examples , 2000 .

[46]  Jan-Markus Schwindt.  Physik , 1973, Universum ohne Dinge.

[47]  Dexuan Xie,et al.  A new minimization protocol for solving nonlinear Poisson–Boltzmann mortar finite element equation , 2007 .

[48]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[49]  W. Doerfler,et al.  A robust adaptive strategy for the nonlinear Poisson equation , 1995, Computing.

[50]  Joseph W. Jerome,et al.  L∞ stability of finite element approximations to elliptic gradient equations , 1990 .

[51]  Ricardo H. Nochetto,et al.  Data Oscillation and Convergence of Adaptive FEM , 2000, SIAM J. Numer. Anal..

[52]  Andreas Veeser,et al.  Convergent adaptive finite elements for the nonlinear Laplacian , 2002, Numerische Mathematik.

[53]  Ricardo H. Nochetto,et al.  Convergence of Adaptive Finite Element Methods , 2002, SIAM Rev..

[54]  Carsten Carstensen,et al.  Convergence of adaptive finite element methods in elastoplasticity , 2005 .

[55]  Ricardo H. Nochetto,et al.  Pointwise a posteriori error estimates for monotone semi-linear equations , 2006, Numerische Mathematik.

[56]  Wolfgang Dahmen,et al.  Adaptive Finite Element Methods with convergence rates , 2004, Numerische Mathematik.

[57]  James H. Bramble,et al.  The analysis of multigrid methods , 2000 .

[58]  María Cecilia Rivara,et al.  Design and data structure of fully adaptive, multigrid, finite-element software , 1984, ACM Trans. Math. Softw..