Finite-Difference Approximations and Local Algorithm for the Poisson and Poisson-Boltzmann Electrostatics

TL;DR

This paper develops finite-difference methods and local algorithms for solving Poisson and Poisson-Boltzmann electrostatics problems, proving convergence, uniqueness, and error bounds, with numerical validation.

Contribution

It introduces a new local algorithm with shift for variable coefficients and provides rigorous analysis of convergence, error estimates, and characterization of minimizers.

Findings

Unique minimizers correspond to solutions of Poisson and PB equations.

Local algorithms converge and satisfy energy minimization conditions.

Error estimates of order h^2 are established for the discretized solutions.

Abstract

We study finite-difference approximations of both Poisson and Poisson-Boltzmann (PB) electrostatic energy functionals for periodic structures constrained by Gauss' law and a class of local algorithms for minimizing the finite-difference discretization of such functionals. The variable of Poisson energy is the vector field of electric displacement and that for the PB energy consists of an electric displacement and ionic concentrations. The displacement is discretized at midpoints of edges of grid boxes while the concentrations are discretize at grid points. The local algorithm is an iteration over all the grid boxes that locally minimizes the energy on each grid box, keeping Gauss' law satisfied. We prove that the energy functionals admit unique minimizers that are solutions to the corresponding Poisson's and charge-conserved PB equation, respectively. Local equilibrium conditions are identified to characterize the finite-difference minimizers of the discretized energy functionals. These conditions are the curl free for the Poisson case and the discrete Boltzmann distributions for the PB case, respectively. Next, we obtain the uniform bound with respect to the grid size h and O(h2)-error estimates in maximum norm for the finite-difference minimizers. The local algorithms are detailed, and a new local algorithm with shift is proposed to treat the general case of a variable coefficient for the Poisson energy. We prove the convergence of all these local algorithms, using the characterization of the finite-difference minimizers. Finally, we present numerical tests to demonstrate the results of our analysis.