A First-Order System Least-Squares Finite Element Method for the Poisson-Boltzmann Equation
A First-Order System Least-Squares Finite Element Method for the Poisson-Boltzmann Equation
复制标题
DOI:
10.1002/jcc.21446
复制
发表时间:
2010-06-01
影响因子:
3
通讯作者:
Olson, Luke N.
中科院分区:
文献类型:
--
作者:
Bond, Stephen D.;Chaudhry, Jehanzeb Hameed;Olson, Luke N.
The Poisson-Boltzmann equation is an important tool in modeling solvent in biomolecular systems. In this article. we focus on numerical approximations to the electrostatic potential expressed in the regularized linear Poisson-Boltzmann equation. We expose the flux directly through a first-order system form of the equation. Using this formulation, we propose a system that yields a tractable least-squares finite element formulation and establish theory to support this approach. The least-squares finite element approximation naturally provides an a posteriori error estimator and we present numerical evidence in support of the method. The computational results highlight optimality in the case of adaptive mesh refinement for a variety of molecular configurations. In particular, we show promising performance for the Born ion. Fasciculin I. methanol, and a dipole, which highlights robustness of our approach. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 31: 1625-1635, 2010