English

# Item

ITEM ACTIONSEXPORT
Fast Solution of the Poisson-Boltzmann Equation with nonaffine Parametrized Boundary Conditions Using the Reduced Basis Method

Benner, P., Feng, L., Kweyu, C. M., & Stein, M. (in preparation). Fast Solution of the Poisson-Boltzmann Equation with nonaffine Parametrized Boundary Conditions Using the Reduced Basis Method.

Item is

show hide
Genre: Paper

### Files

show Files
hide Files
:
1705.08349.pdf (Preprint), 4MB
Name:
1705.08349.pdf
Description:
Visibility:
Public
MIME-Type / Checksum:
application/pdf / [MD5]
-
-

show

### Creators

show
hide
Creators:
Benner, Peter1, Author
Feng, Lihong1, Author
Kweyu, Cleophas M.1, 2, Author
Stein, Matthias3, Author
Affiliations:
1Computational Methods in Systems and Control Theory, Max Planck Institute for Dynamics of Complex Technical Systems, Max Planck Society, ou_1738141
2International Max Planck Research School (IMPRS), Max Planck Institute for Dynamics of Complex Technical Systems, Max Planck Society, DE, ou_1738143
3Molecular Simulations and Design, Max Planck Institute for Dynamics of Complex Technical Systems, Max Planck Society, ou_1738148

### Content

show
hide
Free keywords: Mathematics, Numerical Analysis, math.NA
Abstract: The Poisson-Boltzmann equation (PBE) is a nonlinear elliptic parametrized partial differential equation that arises in biomolecular modeling and is a fundamental tool for structural biology. It is used to calculate electrostatic potentials around an ensemble of fixed charges immersed in an ionic solution. Efficient numerical computation of the PBE yields a high number of degrees of freedom in the resultant algebraic system of equations, ranging from several hundred thousands to millions. Coupled with the fact that in most cases the PBE requires to be solved multiple times for a large number of system configurations, this poses great computational challenges to conventional numerical techniques. To accelerate such computations, we here present the reduced basis method (RBM) which greatly reduces this computational complexity by constructing a reduced order model of typically low dimension. We discretize the linearized PBE (LPBE) with a centered finite difference scheme and solve the resultant linear system by the preconditioned conjugate gradient (PCG) method with an algebraic multigrid (AMG) V-cycle as preconditioner at different samples of ionic strength on a three-dimensional Cartesian grid. We then apply the RBM to the high-fidelity full order model (FOM). The discrete empirical interpolation method (DEIM) is applied to the Dirichlet boundary conditions which are nonaffine in one parameter (the ionic strength) to reduce the complexity of the reduced order model (ROM). From the numerical results, we notice that the RBM reduces the model order from $\mathcal{N} = 2\times 10^{6}$ to $N = 6$ at an accuracy of $10^{-10}$ and reduces computational time by a factor of approximately $8,000$. DEIM, on the other hand, provides a speed-up of $20$ in the online phase at a single iteration of the greedy algorithm.

### Details

show
hide
Language(s):
Dates: 2017-05-23
Publication Status: Not specified
Pages: 22 pages, 14 figures
Publishing info: -
Rev. Method: -
Identifiers: arXiv: 1705.08349
URI: http://arxiv.org/abs/1705.08349
Degree: -

show

show

show

show