Goal-oriented adaptivity and multilevel preconditioning for the poisson-boltzmann equation

Journal of Scientific Computing

Volumn 52, Issue 1, 2012, Pages 202-225

  Aksoylu, Burak ^a,b   Bond, Stephen D ^c   Cyr, Eric C ^c   Holst, Michael ^d  

^a TOBB UNIVERSITY OF ECONOMICS AND TECHNOLOGY   (Turkey)

^b LOUISIANA STATE UNIVERSITY   (United States)

^c SANDIA NATIONAL LABORATORIES   (United States)

^d UNIVERSITY OF CALIFORNIA   (United States)

Author keywords

Adaptive finite element methods; Electrostatics; Goal oriented a posteriori error estimation; Multilevel preconditioning; Poisson Boltzmann equation; Solvation free energy

Indexed keywords

ADAPTIVE FINITE ELEMENT METHODS; MULTILEVEL PRECONDITIONING; POISSON-BOLTZMANN EQUATIONS; POSTERIORI ERROR ESTIMATION; SOLVATION FREE ENERGIES;

ADAPTIVE ALGORITHMS; ELECTROSTATICS; FINITE ELEMENT METHOD; FREE ENERGY; SOLVATION;

POISSON EQUATION;

EID: 84861455916     PISSN: 08857474     EISSN: None     Source Type: Journal    
DOI: 10.1007/s10915-011-9539-6     Document Type: Article

References (86)

1. McQuarrie, D.A.: Statistical Mechanics. Harper and Row, New York (1973)
2. Tanford, C.: Physical Chemistry of Macromolecules. Wiley, New York (1961)
3. N. Baker D. Sept S. Joseph M. Holst J.A. McCammon 2001 Electrostatics of nanosystems: application to microtubules and the ribosome Proc. Natl. Acad. Sci. USA 98 10037 10041 10.1073/pnas.181342398 (Pubitemid 32802969)
4. B.R. Brooks C.L. Brooks III A.D.Jr. Mackerell L. Nilsson R.J. Petrella B. Roux Y. Won G. Archontis C. Bartels S. Boresch A. Caflisch L. Caves Q. Cui A.R. Dinner M. Feig S. Fischer J. Gao M. Hodoscek W. Im K. Kuczera T. Lazaridis J. Ma V. Ovchinnikov E. Paci R.W. Pastor C.B. Post J.Z. Pu M. Schaefer B. Tidor R.M. Venable H.L. Woodcock X. Wu W. Yang D.M. York M. Karplus 2009 CHARMM: The biomolecular simulation program J. Comput. Chem. 30 10 1545 1614 10.1002/jcc.21287
5. W. Rocchia S. Sridharan A. Nicholls E. Alexov A. Chiabrera B. Honig 2002 Rapid grid-based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objects J. Comput. Chem. 23 1 128 137 10.1002/jcc.1161 (Pubitemid 34063137)
6. J.D. Madura J.M. Briggs R.C. Wade M.E. Davis B.A. Luty A. Ilin J. Antosiewicz M.K. Gilson B. Bagheri L.R. Scott J.A. McCammon 1995 Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian Dynamics program Comput. Phys. Comm. 91 1-3 57 95 10.1016/0010-4655(95)00043-F
7. M. Holst N. Baker F. Wang 2000 Adaptive multilevel finite element solution of the Poisson-Boltzmann equation I: algorithms and examples J. Comput. Chem. 21 1319 1342 10.1002/1096-987X(20001130)21:15<1319::AID-JCC1>3.0. CO;2-8
8. N. Baker M. Holst F. Wang 2000 Adaptive multilevel finite element solution of the Poisson-Boltzmann equation II: refinement at solvent accessible surfaces in biomolecular systems J. Comput. Chem. 21 1343 1352 10.1002/1096-987X(20001130)21:15<1343::AID-JCC2>3.0.CO;2-K
9. L. Chen M. Holst J. Xu 2007 The finite element approximation of the nonlinear Poisson-Boltzmann equation SIAM J. Numer. Anal. 45 6 2298 2320 2361891 10.1137/060675514
10. Bangerth, W., Rannacher, R.: Adaptive Finite Element Methods for Differential Equations. Birkhauser, Boston (2003)
11. M. Holst J.A. McCammon Z. Yu Y.C. Zhou Y. Zhu 2012 Adaptive finite element modeling techniques for the Poisson-Boltzmann equation Commun. Comput. Phys. 11 179 214 2841952
12. J.H. Chaudhry S.D. Bond L.N. Olson 2011 Finite element approximation to a finite-size modified Poisson-Boltzmann equation J. Sci. Comput. 47 3 347 364 2793588 10.1007/s10915-010-9441-7
13. Baker, N.A., Bashford, D., Case, D.A.: Implicit solvent electrostatics in biomolecular simulation. In: Leimkuhler, B., Chipot, C., Elber, R., Laaksonen, A., Mark, A., Schlick, T., Schutte, C., Skeel, R. (eds.) New Algorithms for Macromolecular Simulation. Lecture Notes in Computational Science and Engineering, vol. 49, pp. 263-295. Springer, Berlin (2006)
14. M.K. Gilson M.E. Davis B.A. Luty J.A. McCammon 1993 Computationn of electrostatic forces on solvated molecules using the Poisson-Boltzmann equation J. Phys. Chem. 97 3591 3600 10.1021/j100116a025
15. Z. Zhou P. Payne M. Vasquez N. Kuhn M. Levitt 1996 Finite-difference solution of the Poisson-Boltzmann equation: complete elimination of self-energy J. Comput. Chem. 17 1344 1351 10.1002/(SICI)1096-987X(199608)17:11<1344::AID- JCC7>3.0.CO;2-M (Pubitemid 126576751)
16. L. Greengard V. Rokhlin 1987 A fast algorithm for particle simulations J. Comput. Phys. 73 2 325 348 918448 10.1016/0021-9991(87)90140-9
17. R.D. Skeel I. Tezcan D.J. Hardy 2002 Multiple grid methods for classical molecular dynamics J. Comput. Chem. 23 6 673 684 10.1002/jcc.10072 (Pubitemid 34312089)
18. Hardy, D.J.: Multilevel summation for the fast evaluation of forces for the simulation of biomolecules. Ph.D. thesis, Department of Computer Science, University of Illinois at Urbana-Champaign (2006)
19. U. Essmann L. Perera M.L. Berkowitz T. Darden H. Lee L.G. Pedersen 1995 A smooth particle mesh Ewald method J. Chem. Phys. 103 19 8577 8593 10.1063/1.470117
20. M. Holst G. Nagy G. Tsogtgerel 2009 Rough solutions of the Einstein constraints on closed manifolds without near-CMC conditions Commun. Math. Phys. 288 2 547 613 2500992 10.1007/s00220-009-0743-2
21. F. Fogolari P. Zuccato G. Esposito P. Viglino 1999 Biomolecular electrostatics with the linearized Poisson-Boltzmann equation Biophys. J. 76 1 1 16 10.1016/S0006-3495(99)77173-0 (Pubitemid 29202433)
22. B. Lu Y.C. Zhou G.A. Huber S.D. Bond M.J. Holst J.A. McCammon 2007 Electrodiffusion: A continuum modeling framework for biomolecular systems with realistic spatiotemporal resolution J. Chem. Phys. 127 13 135102 10.1063/1.2775933 (17 pages) (Pubitemid 47531645)
23. A.H. Boschitsch M.O. Fenley 2004 Hybrid boundary element and finite difference method for solving the nonlinear Poisson-Boltzmann equation J. Comput. Chem. 25 7 935 955 10.1002/jcc.20000
24. B. Lu Y. Zhou M. Holst J.A. McCammon 2008 Recent progress in numerical methods for the Poisson-Boltzmann equation in biophysical applications Commun. Comput. Phys. 3 5 973 1009
25. W. Im D. Beglov B. Roux 1998 Continuum solvation model: computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equation Comput. Phys. Commun. 111 59 75 10.1016/S0010-4655(98)00016-2 (Pubitemid 128400465)
26. J.A. Wagoner N.A. Baker 2006 Assessing implicit models for nonpolar mean solvation forces: The importance of dispersion and volume terms Proc. Natl. Acad. Sci. USA 103 22 8331 8336 10.1073/pnas.0600118103 (Pubitemid 43838455)
27. S. Yu W. Geng G.W. Wei 2007 Treatment of geometric singularities in implicit solvent models J. Chem. Phys. 126 244108 10.1063/1.2743020 (13 pages) (Pubitemid 47141404)
28. Strang, G., Fix, G.J.: An Analysis of the Finite Element Method. Prentice Hall, New York (1973)
29. Brenner, S.C., Scott, L.R.: The Mathematical Theory of Finite Element Methods, 2nd edn. Springer, Berlin (2002)
30. Braess, D.: Finite elements: Theory, Fast Solvers, and Applications in Solid Mechanics, 2nd edn. Cambridge University Press, Cambridge (2001)
31. Ern, A., Guermond, J.-L.: Theory and Practice of Finite Elements. Applied Mathematical Sciences, vol. 159, 1st edn. Springer, Berlin (2004)
32. S.C. Eisenstat H.F. Walker 1994 Globally convergent inexact Newton methods SIAM J. Optim. 4 2 393 422 1273766 10.1137/0804022
33. R.E. Bank D.J. Rose 1982 Analysis of a multilevel iterative method for nonlinear finite element equations Math. Comput. 39 160 453 465 669639 10.1090/S0025-5718-1982-0669639-X
34. Hackbusch, W.: Multi-grid Methods and Applications. Springer, Berlin (1985)
35. R.E. Bank T.F. Dupont 1981 An optimal order process for solving finite element equations Math. Comput. 36 153 35 51 595040 10.1090/S0025-5718-1981- 0595040-2
36. J. Xu 1992 Iterative methods by space decomposition and subspace correction SIAM Rev. 34 4 581 613 1193013 10.1137/1034116
37. Holst, M.: The Poisson-Boltzmann equation: Analysis and multilevel numerical solution (Monograph based on the Ph.D. Thesis: Multilevel Methods for the Poisson-Boltzmann Equation). Tech. rep., Applied Mathematics and CRPC, California Institute of Technology (1994)
38. M. Holst F. Saied 1993 Multigrid solution of the Poisson-Boltzmann equation J. Comput. Chem. 14 1 105 113 10.1002/jcc.540140114
39. M. Holst F. Saied 1995 Numerical solution of the nonlinear Poisson-Boltzmann equation: Developing more robust and efficient methods J. Comput. Chem. 16 3 337 364 10.1002/jcc.540160308
40. Bank, R.E., Xu, J.: The hierarchical basis multigrid method and incomplete LU decomposition. In: Keyes, D., Xu, J. (eds.) Seventh International Symposium on Domain Decomposition Methods for Partial Differential Equations, pp. 163-173. AMS, Providence (1994)
41. R.E. Bank J. Xu 1996 An algorithm for coarsening unstructured meshes Numer. Math. 73 1 36 1379277 10.1007/s002110050181
42. A. Brandt 1986 Algebraic multigrid theory: The symmetric case Appl. Math. Comput. 19 23 56 849831 10.1016/0096-3003(86)90095-0
43. Chan, T.F., Smith, B., Zou, J.: Overlapping Schwarz methods on unstructured meshes using non-matching coarse grids. Tech. Rep. CAM 94-8, Department of Mathematics, UCLA (1994)
44. Chan, T.F., Go, S., Zikatanov, L.: Lecture notes on multilevel methods for elliptic problems on unstructured meshes. Tech. rep., Dept. of Mathematics, UCLA (1997)
45. Brandt, A., McCormick, S., Ruge, J.: Algebraic multigrid (AMG) for sparse matrix equations. In: Evans, D.J. (ed.) Sparsity and Its Applications. Cambridge Univ. Press, Cambridge (1984)
46. Ruge, J.W., Stüben, K.: Algebraic multigrid (AMG). In: McCormick, S.F. (ed.) Multigrid Methods. Frontiers in Applied Mathematics, vol. 3, pp. 73-130. Philadelphia, SIAM (1987)
47. Vanek, P., Mandel, J., Brezina, M.: Algebraic multigrid on unstructured meshes. Tech. Rep. UCD/CCM 34, Center for Computational Mathematics, University of Colorado at Denver (1994)
48. Vanek, P., Mandel, J., Brezina, M.: Algebraic multigrid by smoothed aggregation for second and fourth order elliptic problems. Tech. Rep. UCD/CCM 36, Center for Computational Mathematics, University of Colorado at Denver (1995)
49. W. Dörfler 1996 A convergent adaptive algorithm for Poisson's equation SIAM J. Numer. Anal. 33 3 1106 1124 1393904 10.1137/0733054 (Pubitemid 126423840)
50. M.C. Rivara 1984 Design and data structure of fully adaptive, multigrid, finite-element software ACM Trans. Math. Softw. 10 3 242 264 791990 10.1145/1271.1274 (Pubitemid 15421453)
51. Ainsworth, M., Oden, J.: A Posteriori Error Estimation in Finite Element Analysis. Wiley, New York (2000)
52. Cyr, E.C.: Numerical methods for computing the free-energy of coarse-grained molecular systems. Ph.D. thesis, University of Illinois at Urbana-Champaign (2008)
53. S.D. Bond J.H. Chaudhry E.C. Cyr L.N. Olson 2010 A first-order systems least-squares finite element method for the Poisson-Boltzmann equation J. Comput. Chem. 31 8 1625 1635
54. Chaudhry, J.H., Bond, S.D., Olson, L.N.: A weighted adaptive least-squares finite element method for the Poisson-Boltzmann equation (2011, submitted)
55. R. Becker R. Rannacher 2001 An optimal control approach to a posteriori error estimation in finite element methods Acta Numer. 10 1 102 2009692 10.1017/S0962492901000010 (Pubitemid 33305809)
56. J.T. Oden S. Prudhomme 2001 Goal-oriented error estimation and adaptivity for the finite element method Comput. Math. Appl. 41 735 756 1822600 10.1016/S0898-1221(00)00317-5 (Pubitemid 32874156)
57. S. Prudhomme J.T. Oden 1999 On goal-oriented error estimation for elliptic problems: application to the control of pointwise errors Comput. Methods Appl. Mech. Eng. 176 1-4 313 331 1665351 10.1016/S0045-7825(98)00343-0
58. Bastian, P.: Locally refined solution of unsymmetric and nonlinear problems. In: Proc. of the 8th GAMM Seminar. Notes on Numerical Fluid Mechanics, vol. 46, pp. 12-21. Vieweg, Wiesbaden (1993)
59. Bastian, P., Wittum, G.: On robust and adaptive multigrid methods. In: Wesseling, P., Hemker, P. (eds.) Proc. of the 4th European Multigrid Conference. Birkhäuser, Basel (1994)
60. J.H. Bramble J.E. Pasciak J. Wang J. Xu 1991 Convergence estimates for product iterative methods with applications to domain decomposition Math. Comput. 57 1 21 1090464 10.1090/S0025-5718-1991-1090464-8
61. M.C. Rivara 1984 Algorithms for refining triangular grids for adaptive and multigrid techniques Int. J. Numer. Methods Eng. 20 4 745 756 739618 10.1002/nme.1620200412 (Pubitemid 14569502)
62. J.H. Bramble J.E. Pasciak J. Xu 1990 Parallel multilevel preconditioners Math. Comput. 55 191 1 22 1023042 10.1090/S0025-5718-1990-1023042-6
63. Briggs, W.L., Henson, V.E., McCormick, S.F.: A Multigrid Tutorial, 2nd edn. Philadelphia, SIAM (2000)
64. Trottenberg, U., Oosterlee, C., Schüller, A.: Multigrid. Academic Press, London (2001)
65. I. Yavneh 2006 Why multigrid methods are so efficient Comput. Sci. Eng. 8 6 12 22 10.1109/MCSE.2006.125 (Pubitemid 44711773)
66. P. Bastian W. Hackbusch G. Wittum 1998 Additive and multiplicative multi-grid - a comparison Computing 60 345 364 1634835 10.1007/BF02684380 (Pubitemid 128509005)
67. Wittum, G.: Multi-grid methods - an introduction. In: Hergert, W., Ernst, A., Däne, M. (eds.) Computational Materials Science. Lect. Notes Phys., vol. 642, pp. 283-311. Springer, Berlin (2004)
68. R.E. Bank 1996 Hierarchical bases and the finite element method Acta Numer. 5 1 43 1624587 10.1017/S0962492900002610
69. R.E. Bank T.F. Dupont H. Yserentant 1988 The hierarchical basis multigrid method Numer. Math. 52 427 458 932709 10.1007/BF01462238
70. H. Yserentant 1993 Old and new convergence proofs for multigrid methods Acta Numer. 2 285 326 1224685 10.1017/S0962492900002385
71. B. Aksoylu M. Holst 2006 Optimality of multilevel preconditioners for local mesh refinement in three dimensions SIAM J. Numer. Anal. 44 3 1005 1025 2231853 10.1137/S0036142902406119 (Pubitemid 47135912)
72. B. Aksoylu S. Bond M. Holst 2003 An odyssey into local refinement and multilevel preconditioning III: Implementation and numerical experiments SIAM J. Sci. Comput. 25 2 478 498 2058071 10.1137/S1064827502407676
73. Aksoylu, B., Bond, S., Holst, M.: Implementation and theoretical aspects of the BPX preconditioner in the three dimensional local mesh refinement setting. Tech. rep., UT-Austin ICES Report 04-50 (2004)
74. Aksoylu, B., Holst, M.: An odyssey into local refinement and multilevel preconditioning I: Optimality of the BPX preconditioner. Tech. rep., UT-Austin ICES Report 05-03 (2005)
75. Aksoylu, B., Holst, M.: An odyssey into local refinement and multilevel preconditioning II: Stabilizing hierarchical basis methods. Tech. rep., UT-Austin ICES Report 05-04 (2005)
76. B. Aksoylu A. Khodakovsky P. Schröder 2005 Multilevel solvers for unstructured surface meshes SIAM J. Sci. Comput. 26 4 1146 1165 2143479 10.1137/S1064827503430138 (Pubitemid 41205807)
77. F. Bornemann B. Erdmann R. Kornhuber 1993 Adaptive multilevel methods in three space dimensions Int. J. Numer. Methods Eng. 36 3187 3203 1236370 10.1002/nme.1620361808
78. F. Bornemann H. Yserentant 1993 A basic norm equivalence for the theory of multilevel methods Numer. Math. 64 455 476 1213412 10.1007/BF01388699
79. J.H. Bramble J.E. Pasciak 1993 New estimates for multilevel algorithms including the V-cycle Math. Comput. 60 202 447 471 1176705
80. W. Dahmen A. Kunoth 1992 Multilevel preconditioning Numer. Math. 63 315 344 1186345 10.1007/BF01385864
81. M.H. le Du P. Marchot P.E. Bougis J.C. Fontecilla-Camps 1992 1.9-Å resolution structure of fasciculin 1, an anti-acetylcholinesterase toxin from green mamba snake venom J. Biol. Chem. 267 31 22122 22130
82. M. Holst 2001 Adaptive numerical treatment of elliptic systems on manifolds Adv. Comput. Math. 15 1-4 139 191 1887732 10.1023/A:1014246117321 (Pubitemid 33634300)
83. T.J. Dolinsky J.E. Nielsen J.A. McCammon N.A. Baker 2004 PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations Nucleic Acids Res. 32 W665 W667 10.1093/nar/gkh381 (Pubitemid 38997419)
84. T. Hayashi M.E. Martone Z. Yu A. Thor M. Doi M. Holst M.H. Ellisman M. Hoshijima 2009 Three-dimensional reconstruction reveals new details of membrane systems for calcium signaling in the heart J. Cell Sci. 122 7 1005 1013 10.1242/jcs.028175
85. Z. Yu M. Holst Y. Cheng J.A. McCammon 2008 Feature-preserving adaptive mesh generation for molecular shape modeling and simulation J. Mol. Graph. Model. 26 1370 1380 10.1016/j.jmgm.2008.01.007
86. Z. Yu M. Holst J.A. McCammon 2008 High-fidelity geometric modeling for biomedical applications Finite Elem. Anal. Des. 44 11 715 723 2441978 10.1016/j.finel.2008.03.004