Abstract
The present paper is dedicated to the rapid solution of boundary integral equations arising from solvation continuum models. We apply a fully discrete wavelet Galerkin scheme for the computation of the apparent surface charge on van der Waals or solvent accessible surfaces. The molecular surface is described in parametric form as a set of four-sided spherical patches. Each patch is exactly represented as rational Bézier surface. Numerical results are reported to illustrate the approach.
Similar content being viewed by others
References
Beylkin G, Coifman R, Rokhlin V (1991) The fast wavelet transform and numerical algorithms. Comm Pure Appl Math 44: 141–183
Brunnett G (1995) Geometric design with trimmed surfaces. Comput Suppl 10: 101–115
Cancés E, Mennucci B (1998) New applications of integral equations methods for solvation continuum models: ionic solutions and liquid crystals. J Math Chem 23: 309–326
Cohen A, Daubechies I, Feauveau J-C (1992) Biorthogonal bases of compactly supported wavelets. Pure Appl Math 45: 485–560
Connolly M (1983) Analytical molecular surface calculation. J Appl Cryst 16: 548–558
Connolly M (2006) Molecular surface guide. http://connolly.best.vwh.net
Dahmen W, Harbrecht H, Schneider R (2006) Compression techniques for boundary integral equations—optimal complexity estimates. SIAM J Numer Anal 43: 2251–2271
Dahmen W, Harbrecht H, Schneider R (2007) Adaptive methods for boundary integral equations: complexity and convergence estimates. Math Comput 76: 1243–1274
Dahmen W, Kunoth A (1992) Multilevel preconditioning. Numer Math 63: 315–344
DALTON, a molecular electronic structure program, Release 2.0 (2005). http://www.kjemi.uio.no/software/dalton/dalton.html
Dietz R, Hoschek J, Jüttler B (1995) Rational patches on quadric surfaces. Comput Aided Des 27: 27–40
Farin G (1999) Curves and surfaces for CAGD: a practical guide. Academic Press
Frediani L, Cammi R, Corni S, Tomasi J (2004) A polarizable continuum model for molecules at diffuse interfaces. J Chem Phys 120: 3893–3907
Greengard L, Rokhlin V (1987) A fast algorithm for particle simulation. J Comput Phys 73: 325–348
Hackbusch W, Khoromskij BN (2000) A sparse \({\mathcal{H}}\) -matrix arithmetic. II: application to multi-dimensional problems. Computing 64: 21–47
Hackbusch W, Nowak ZP (1989) On the fast matrix multiplication in the boundary element method by panel clustering. Numer Math 54: 463–491
Harbrecht H, Randrianarivony M (2007) From computer aided design to wavelet BEM. Bericht 07-18. Berichtsreihe des Mathematischen Seminars, Christian-Albrechts-University of Kiel, Germany (to appear in Comput Vis Sci)
Harbrecht H, Schneider R (2004) Biorthogonal wavelet bases for the boundary element method. Math Nachr 269–270: 167–188
Harbrecht H, Schneider R (2002) Wavelet galerkin schemes for boundary integral equations—implementation and quadrature. SIAM J Sci Comput 27: 1347–1370
Hochmuth R (1998) A-posteriori estimates and adaptive schemes for transmission problems. J Integral Equ Appl 10: 1–50
Hoschek J, Lasser D (1993) Fundamentals of computer aided geometric design. A.K. Peters, Wellesley
Hu Q, Wang G (2007) Necessary and sufficient conditions for rational quartic representation of conic sections. J Comput Appl Math 203: 190–208
Korte B, Vygen J (2006) Combinatorial optimization: theory and algorithms, 3rd edn. Springer, Berlin
Laug P, Borouchaki H (2002) Molecular surface modeling and meshing. Eng Comput 18: 199–210
Miertuš S, Scrocco E, Tomasi J (1981) Electrostatic interaction of a solute with a continuum. A direct utilization of ab initio molecular potentials for the prevision of solvent effects. Chem Phys 55: 117–129
Pomelli CS (2007) Cavity surfaces and their discretization. In: Mennucci B (ed), Cammi R (Co-Editor) Continuum solvation models in chemical physics: from theory to applications. Wiley, pp 49–63
Piegl L, Tiller W (1995) The NURBS book. Springer, Berlin
Pomelli C (2004) A tessellationless integration grid for the polarizable continuum model reaction field. J Comp Chem 25: 1532–1541
Randrianarivony M (2006) Geometric processing of CAD data and meshes as input of integral equation solvers. PhD thesis, Technische Universität Chemnitz
Randrianarivony M (2008) Quadrilateral decomposition by two-ear property resulting in CAD segmentation. In: Proc. Fifth international conference on computational geometry, pp 179–185
Sanner M, Olsen A, Spehner J (1996) Reduced surface: an efficient way to compute molecular surfaces. Biopolymers 38: 305–320
Schneider R (1998) Multiskalen- und Wavelet-Matrixkompression: Analysisbasierte Methoden zur Lösung großer vollbesetzter Gleichungssysteme. Teubner, Stuttgart
Steinbach O (2008) Numerical approximation methods for elliptic boundary value problems. Finite and boundary elements. Springer, New York
Tomasi J, Mennucci B, Cammi R (2005) Quantum mechanical continuum solvation models. Chem Rev 105: 2999–3094
Wang G-J, Wang G-Z (1992) The rational cubic Bézier representation of conics. Comput Aided Geom Des 9: 447–455
Weijo V, Randrianarivony M, Harbrecht H, Frediani L (2009) Wavelet formulation of the Polarizable Continuum Model. INS-Preprint 0906, Institute for Numerical Simulation, University of Bonn (submitted)
Whitley D (1998) Van der Waals graphs and molecular shape. J Math Chem 23: 377–397
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by W. Hackbusch.
Rights and permissions
About this article
Cite this article
Harbrecht, H., Randrianarivony, M. Wavelet BEM on molecular surfaces: parametrization and implementation. Computing 86, 1–22 (2009). https://doi.org/10.1007/s00607-009-0050-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00607-009-0050-y