Hierarchical Interpolative Factorization for Self Green's Function in 3D Modified Poisson-Boltzmann Equations

Yihui Tu,Zhenli Xu,Haizhao Yang
2023-07-03
Abstract:The modified Poisson-Boltzmann (MPB) equations are often used to describe equilibrium particle distribution of ionic systems. In this paper, we propose a fast algorithm to solve MPB equations with the self Green's function as the self energy in three dimensions, where the solution of the self Green's function poses a computational bottleneck due to the need to solve a high-dimensional partial differential equation. Our algorithm combines the selected inversion with hierarchical interpolative factorization for the self Green's function by extending our previous result of two dimensions. This leads to an $O(N\log N)$ algorithm through the strategical utilization of locality and low-rank characteristics of the corresponding operators. Furthermore, the estimated $O(N)$ complexity is obtained by applying cubic edge skeletonization at each level for thorough dimensionality reduction. Extensive numerical results are performed to demonstrate the accuracy and efficiency of the proposed algorithm for problems in three dimensions.
Computational Physics
What problem does this paper attempt to address?