Chebyshev-Galerkin algorithm for computing 3D Solvation
Contact the author: Please use for correspondence this email.
Submission date: 19. Dec. 2011
Download full preprint: PDF (381 kB)
To solve 3D integral equations of molecular liquids, we have developed a numer- ical scheme based on the Galerkin method. Using a domain decomposition of the interface region, we reduce the problem to calculations of approximating coefficients and the kernel matrix in spherical shell elements (SSE). Applying the linear trans- formation of coordinates for each SSE we result in calcualtions of the approximating coefficients in cubic volumes. Using the conventional triple Chebyshev series as a basis set, we derive formulas for calculations of the approximating coefficients and evaluate the computational costs of these operations. We have described the general properties of the Chebyshev-Galerkin matrix and derived analytical expressions for recursion calculations of the matrix elements. We have also outlined an iterative method for the solutions of the nonlinear equations obtained for the approximating coefficients, which is based on the direct inversion in the iterative space. Finally, we have estimated the total computational cost of the proposed scheme and compared it with current algorithms for computing 3D solvation problem. It was found the proposed scheme to be by 2-3 orders eective than the current algorithms based on the uniform FFT.