Accuracy and efficiency in computing electrostatic potential for an ion channel model in layered dielectric/electrolyte media

Authors:
Huimin Lin;Huazhong Tang;Wei Cai
Affiliations:
HEDPS, CAPT & LMAM, School of Mathematical Sciences, Peking University, Beijing 100871, China;HEDPS, CAPT & LMAM, School of Mathematical Sciences, Peking University, Beijing 100871, China;Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223, USA
Venue:
Journal of Computational Physics
Year:
2014

Citing 3
Cited 0

Extending the fast multipole method to charges inside or outside a dielectric sphere

Journal of Computational Physics
Fast Analytical Methods for Macroscopic Electrostatic Models in Biomolecular Simulations

SIAM Review
Image Approximations to Electrostatic Potentials in Layered Electrolytes/Dielectrics and an Ion-Channel Model

Journal of Scientific Computing

Quantified Score

Hi-index	31.45

Visualization

Abstract

This paper will investigate the numerical accuracy and efficiency in computing the electrostatic potential for a finite-height cylinder, used in an explicit/implicit hybrid solvation model for ion channel and embedded in a layered dielectric/electrolyte medium representing a biological membrane and ionic solvents. A charge locating inside the cylinder cavity, where ion channel proteins and ions are given explicit atomistic representations, will be influenced by the polarization field of the surrounding implicit dielectric/electrolyte medium. Two numerical techniques, a specially designed boundary integral equation method and an image charge method, will be investigated and compared in terms of accuracy and efficiency for computing the electrostatic potential. The boundary integral equation method based on the three-dimensional layered Green@?s functions provides a highly accurate solution suitable for producing a benchmark reference solution, while the image charge method is found to give reasonable accuracy and highly efficient and viable to use the fast multipole method for interactions of a large number of charges in the atomistic region of the hybrid solvation model.