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.

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

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

Analytic regularity of strong solutions for the complexified stochastic non-linear Poisson Boltzmann Equation

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

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

An Exact Solution of the Nonlinear Poisson-Boltzmann Equation in Parallel-Plate Geometry

A Direct Simulation Approach for the Poisson-Boltzmann Equation Using the Random Batch Method.

Perturbation Solutions for the Nonlinear Poisson–Boltzmann Equation with a High-Order-accuracy Debye–Hückel Approximation

Highly Accurate and Simple Analytical Approach to Nonlinear Poisson–Boltzmann Equation

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

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

Optimized Parallelization of Boundary Integral Poisson-Boltzmann Solvers

On an interior-exterior nonoverlapping domain decomposition method for the Poisson--Boltzmann equation

Regularized numerical methods for the nonlinear Schrödinger equation with singular nonlinearity

N ov 2 01 8 Stochastic Poisson-Boltzmann Equation for Charged Fluids

An Investigation of Physics Informed Neural Networks to solve the Poisson-Boltzmann Equation in Molecular Electrostatics

Implicit Temporal Discretization and Exact Energy Conservation for Particle Methods Applied to the Poisson–Boltzmann Equation

Parallel Higher-order Boundary Integral Electrostatics Computation on Molecular Surfaces with Curved Triangulation

Poisson-Boltzmann model for protein-surface electrostatic interactions and grid-convergence study using the PyGBe code

Fast Iterative Method for Local Steric Poisson-Boltzmann Theories in Biomolecular Solvation