| Popis: |
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. |
