Abstract:The Poisson--Boltzmann equation is widely used to model electrostatics in molecular systems. Available software packages solve it using finite difference, finite element, and boundary element methods, where the latter is attractive due to the accurate representation of the molecular surface and partial charges, and exact enforcement of the boundary conditions at infinity. However, the boundary element method is limited to linear equations and piecewise constant variations of the material properties. In this work, we present a scheme that couples finite and boundary elements for the Poisson--Boltzmann equation, where the finite element method is applied in a confined {\it solute} region, and the boundary element method in the external {\it solvent} region. As a proof-of-concept exercise, we use the simplest methods available: Johnson--Nédélec coupling with mass matrix and diagonal preconditioning, implemented using the Bempp-cl and FEniCSx libraries via their Python interfaces. We showcase our implementation by computing the polar component of the solvation free energy of a set of molecules using a constant and a Gaussian-varying permittivity. We validate our implementation against the finite difference code APBS (to 0.5\%), and show scaling from protein G B1 (955 atoms) up to immunoglobulin G (20\,148 atoms). For small problems, the coupled method was efficient, outperforming a purely boundary integral approach. For Gaussian-varying permittivities, which are beyond the applicability of boundary elements alone, we were able to run medium to large sized problems on a single workstation. Development of better preconditioning techniques and the use of distributed memory parallelism for larger systems remains an area for future work. We hope this work will serve as inspiration for future developments for molecular electrostatics with implicit solvent models.

A Method of Boundary Equations for Nonlinear Poisson–Boltzmann Equation Arising in Biomolecular Systems

A kernel-free boundary integral method for the nonlinear Poisson-Boltzmann equation

An Approach to Numerically Solving the Poisson Equation

Domain Decomposition Method for Poisson--Boltzmann Equations based on Solvent Excluded Surface

Reliable numerical solution of a class of nonlinear elliptic problems generated by the Poisson-Boltzmann equation

A Numerical Method For Solving The One-Dimensional Vlasov-Poisson Equation In Phase Space

Coupling finite and boundary element methods to solve the Poisson--Boltzmann equation for electrostatics in molecular solvation

Solution decomposition for the nonlinear Poisson-Boltzmann equation using the range-separated tensor format

NUMERICAL SOLUTION OF THE BOUNDARY PROBLEM FOR PARABOLIC EQUATION WITH NON-SELF-ADJOINT OPERATOR AND RELATED BOUNDARY CONDITIONS

The Method of Polynomial Particular Solutions for Solving Nonlinear Poisson-Type Equations

Solution of 1D Poisson Equation with Neumann-Dirichlet and Dirichlet-Neumann Boundary Conditions, Using the Finite Difference Method

A robust and non-singular formulation of the boundary integral method for the potential problem

Lattice evolution solution for the nonlinear Poisson-Boltzmann equation in confined domains

Error Analysis of Virtual Element Method for the Poisson-Boltzmann Equation

A stabilized finite volume element method for solving Poisson–Nernst–Planck equations

A Nonlinear Boundary Condition for Continuum Models of Biomolecular Electrostatics

An Efficient Finite Element Iterative Method for Solving a Nonuniform Size Modified Poisson-Boltzmann Ion Channel Model

High order solution of Poisson problems with piecewise constant coefficients and interface jumps

Computing the Newton Potential in the Boundary Integral Equation for the Dirichlet Problem of the Poisson Equation

Accurate boundary-integral formulations for the calculation of electrostatic forces with an implicit-solvent model

Optimized Parallelization of Boundary Integral Poisson-Boltzmann Solvers