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

Abstract

The modified Poisson-Boltzmann (MPB) equations are often used to describe the equilibrium particle distribution of ionic systems. In this paper, we propose a fast algorithm to solve the 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 requirement of solving a high-dimensional partial differential equation. Our algorithm combines the selected inversion with hierarchical interpolative factorization for the self Green’s function, building upon our previous result of two dimensions. This approach yields an algorithm with a complexity of $O(N\log N)$ by strategically leveraging the locality and low-rank characteristics of the corresponding operators. Additionally, the theoretical O(N) complexity is obtained by applying cubic edge skeletonization at each level for thorough dimensionality reduction. Extensive numerical results are conducted to demonstrate the accuracy and efficiency of the proposed algorithm for problems in three dimensions.

Keywords

Selected inversion / Hierarchical interpolative factorization / Linear scaling / Self Green’s function / Modified Poisson-Boltzmann (MPB) equations / 65N22 / 65F50 / 82B21

References

[1] Bazant MZ, Storey BD, Kornyshev AA. Double layer in ionic liquids: overscreening versus crowding. Phys. Rev. Lett.. 2011, 106(4): 046102

[2] Boroudjerdi H, Kim Y-W, Naji A, Netz RR, Schlagberger X, Serr A. Statics and dynamics of strongly charged soft matter. Phys. Rep.. 2005, 416: 129-199

[3] Borukhov I, Andelman D, Orland H. Steric effects in electrolytes: a modified Poisson-Boltzmann equation. Phys. Rev. Lett.. 1998, 79(3): 435-438

[4] Brandt D. Multi-level adaptive solutions to boundary-value problems. Math. Comput.. 1977, 138: 333-390

[5] Chapman DL. A contribution to the theory of electrocapillarity. Phil. Mag.. 1913, 25: 475-481

[6] Cheng H, Gimbutas Z, Martinsson P, Rokhlin V. On the compression of low rank matrices. SIAM J. Sci. Comput.. 2005, 26(4): 1389-1404

[7] Corry B, Kuyucak S, Chung SH. Dielectric self-energy in Poisson-Boltzmann and Poisson-Nernst-Planck models of ion channels. Biophys. J.. 2003, 84(6): 3594-3606

[8] Daiguji H, Yang P, Majumdar A. Ion transport in nanofluidic channels. Nano Lett.. 2004, 4(1): 137-142

[9] Duff IS, Reid JK. The multifrontal solution of indefinite sparse symmetric linear. ACM Trans. Math. Softw.. 1983, 9(3): 302-325

[10] George A. Nested dissection of a regular finite element mesh. SIAM J. Numer. Anal.. 1973, 10(2): 345-363

[11] Gillman A, Martinsson PG. A direct solver with ${O(N)}$ complexity for variable coefficient elliptic PDEs discretized via a high-order composite spectral collocation method. SIAM J. Sci. Comput.. 2013, 36(4): 2023-2046

[12] Gillman A, Martinsson PG. An ${O(N)}$ algorithm for constructing the solution operator to 2D elliptic boundary value problems in the absence of body loads. Adv. Comput. Math.. 2014, 40(4): 773-796

[13] Gouy G. Constitution of the electric charge at the surface of an electrolyte. J. Phys.. 1910, 9: 457-468

[14] Grasedyck L, Kriemann R, Borne SL. Domain-decomposition based ${{{\mathcal{H} }}}$-LU preconditioners. Numer. Math.. 2009, 112(4): 565-600

[15] Ho KL, Ying L. Hierarchical interpolative factorization for elliptic operators: differential equations. Comm. Pure Appl. Math.. 2015, 69(8): 1415-1451

[16] Ho KL, Ying L. Hierarchical interpolative factorization for elliptic operators: integral equations. Comm. Pure Appl. Math.. 2016, 69(7): 1314-1353

[17] Ji L, Liu P, Xu Z, Zhou S. Asymptotic analysis on dielectric boundary effects of modified Poisson-Nernst-Planck equations. SIAM J. Appl. Math.. 2018, 78: 1802-1822

[18] Kenneth LH. FLAM: fast linear algebra in MATLAB—algorithms for hierarchical matrices. J. Open Source Softw.. 2020, 5: 1906

[19] Liljeström V, Seitsonen J, Kostiainen M. Electrostatic self-assembly of soft matter nanoparticle cocrystals with tunable lattice parameters. ACS Nano. 2015, 9(11): 11278-11285

[20] Lin L, Lu J, Ying L, Car R, E W. Fast algorithm for extracting the diagonal of the inverse matrix with application to the electronic structure analysis of metallic systems. Commun. Math. Sci.. 2009, 7(3): 755-777

[21] Lin L, Yang C, Lu J, Ying L, E W. A fast parallel algorithm for selected inversion of structured sparse matrices with application to 2D electronic structure calculations. SIAM J. Sci. Comput.. 2011, 33(3): 1329-1351

[22] Lin L, Yang C, Meza JC, Lu J, Y L, E W. SelInv–an algorithm for selected inversion of a sparse symmetric matrix. ACM Trans. Math. Softw.. 2011, 37(4): 40:1-40:19

[23] Liu C, Wang C, Wise S, Yue X, Zhou S. A positivity-preserving, energy stable and convergent numerical scheme for the Poisson-Nernst-Planck system. Math. Comput.. 2021, 90: 2071-2106

[24] Liu H, Wang Z. A free energy satisfying finite difference method for Poisson-Nernst-Planck equations. J. Comput. Phys.. 2014, 268(2): 363-376

[25] Liu J-L, Eisenberg RS. Molecular mean-field theory of ionic solutions: a Poisson-Nernst-Planck-Bikerman model. Entropy. 2020, 22(5): 550

[26] Liu JWH. The multifrontal method for sparse matrix solution: theory and practice. SIAM Rev.. 1992, 34(1): 82-109

[27] Liu P, Ji X, Xu Z. Modified Poisson-Nernst-Planck model with accurate Coulomb correlation in variable media. SIAM J. Appl. Math.. 2018, 78: 226-245

[28] Ma M, Xu Z. Self-consistent field model for strong electrostatic correlations and inhomogeneous dielectric media. J. Chem. Phys.. 2014, 141(24): 244903

[29] Ma M, Xu Z, Zhang L. Modified Poisson-Nernst-Planck model with Coulomb and hard-sphere correlations. SIAM J. Appl. Math.. 2021, 81: 1645-1667

[30] Netz RR, Orland H. Beyond Poisson-Boltzmann: fluctuation effects and correlation functions. Eur. Phys. J. E. 2000, 1(2): 203-214

[31] Netz RR, Orland H. Variational charge renormalization in charged systems. Eur. Phys. J. E. 2003, 11(3): 301-311

[32] Podgornik R. Electrostatic correlation forces between surfaces with surface specific ionic interactions. J. Chem. Phys.. 1989, 91: 5840-5849

[33] Schmitz PG, Ying L. A fast direct solver for elliptic problems on general meshes in 2D. J. Comput. Phys.. 2012, 231(4): 1314-1338

[34] Schoch RB, Han J, Renaud P. Transport phenomena in nanofluidics. Rev. Mod. Phys.. 2008, 80: 839-883

[35] Tu Y, Pang Q, Yang H, Xu Z. Linear-scaling selected inversion based on hierarchical interpolative factorization for self Green’s function for modified Poisson-Boltzmann equation in two dimensions. J. Comput. Phys.. 2022, 461: 110893

[36] Wang Z-G. Fluctuation in electrolyte solutions: the self energy. Phys. Rev. E. 2010, 81: 021501

[37] Xia J, Chandrasekaran S, Gu M, Li X. Superfast multifrontal method for large structured linear systems of equations. SIAM J. Matrix Anal. Appl.. 2009, 31(3): 1382-1411

[38] Xia J, Xi Y, Cauley S, Balakrishnan V. Fast sparse selected inversion. SIAM J. Matrix Anal. Appl.. 2015, 36(3): 1283-1314

[39] Xu Z, Maggs AC. Solving fluctuation-enhanced Poisson-Boltzmann equations. J. Comput. Phys.. 2014, 36(3): 310-322

[40] Xu Z, Ma M, Liu P. Self-energy-modified Poisson-Nernst-Planck equations: WKB approximation and finite-difference approaches. Phys. Rev. E. 2014, 90(1): 013307

Funding

National Natural Science Foundation of China(12071288 and 12325113)

Science and Technology Innovation Plan Of Shanghai Science and Technology Commission(21JC1403700)

National Science Foundation(DMS-2244988 and DMS-2206333)

Office of Naval Research(N00014-23-1-2007)

RIGHTS & PERMISSIONS

Shanghai University