Optimal Error Estimates of Compact Finite Difference Discretizations for the Schrodinger-Poisson System

We study compact finite difference methods for the Schrodinger-Poisson equation in a bounded domain and establish their optimal error estimates under proper regularity assumptions on wave functionψand external potential V(x).The CrankNicolson compact finite difference method and the semi-implicit compact finite difference method are both of order O(h^(4)+τ^(2))in discrete l^(2),H^(1) and l^(∞) norms with mesh size h and time step τ.For the errors of compact finite difference approximation to the second derivative and Poisson potential are nonlocal,thus besides the standard energy method and mathematical induction method,the key technique in analysis is to estimate the nonlocal approximation errors in discrete l^(∞) and H^(1) norm by discrete maximum principle of elliptic equation and properties of some related matrix.Also some useful inequalities are established in this paper.Finally,extensive numerical results are reported to support our error estimates of the numerical methods.