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 finite volume method to solve the Poisson equation with jump conditions and surface charges: application to electroporation

An Approach to Numerically Solving the Poisson Equation

A numerical strategy to discretize and solve the Poisson equation on dynamically adapted multiresolution grids for time-dependent streamer discharge simulations

Finite-Difference Approximations and Local Algorithm for the Poisson and Poisson-Boltzmann Electrostatics

Computer simulations of electrodiffusion problems based on Nernst–Planck and Poisson equations

Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

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

Numerical Methods for a Poisson-Nernst-Planck-Fermi Model

An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications

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

Second-order, Positive, and Unconditional Energy Dissipative Scheme for Modified Poisson-Nernst-Planck Equations

A meshless stochastic method for Poisson-Nernst-Planck equations

Second-order, Positive, and Unconditional Energy Dissipative Scheme for Modified Poisson–Nernst–Planck Equations

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

Structure-Preserving and Efficient Numerical Methods for Ion Transport

Numerical Solution to Poisson's Equation for Estimating Electrostatic Properties Resulting from an Axially Symmetric Gaussian Charge Density Distribution: Charge in Free Space

A discontinuous Galerkin solver for Boltzmann Poisson systems in nano devices

A meshless stochastic method for Poisson–Nernst–Planck equations

Scalable approximation and solvers for ionic electrodiffusion in cellular geometries

Numerical Methods of the Maxwell-Stefan Diffusion Equations and Applications in Plasma and Particle Transport

A positivity-preserving, linear, energy stable and convergent numerical scheme for the Poisson–Nernst–Planck (PNP) system