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.

Dielectric Boundary Force in Molecular Solvation with the Poisson-Boltzmann Free Energy: A Shape Derivative Approach.

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

Thermodynamics of Electrolyte Solutions Near Charged Surfaces: Constant Surface Charge vs. Constant Surface Potential

Entropic Origin of Ionic Interactions in Polar Solvents

Some challenges of diffused interfaces in implicit-solvent models

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions

Calculation of the Maxwell Stress Tensor and the Poisson-Boltzmann Force on A Solvated Molecular Surface Using Hypersingular Boundary Integrals

Incorporating Dipolar Solvents with Variable Density in Poisson-Boltzmann Electrostatics

The Poisson Boltzmann equation and the charge separation phenomenon at the silica-water interface: A holistic approach

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

Mixed finite element analysis for a modified Poisson–Fermi interface problem accounting for electrostatic correlations

Ion Transport in Dipolar Medium I: A Local Dielectric Poisson–Nernst–Planck/Poisson–Boltzmann Model

A Nonlinear Boundary Condition for Continuum Models of Biomolecular Electrostatics

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

Production of Charm, Charmonium, and Beauty in 800$\,$GeV Proton-Nucleon Collisions

Local Structure-Preserving Relaxation Method for Charged Systems on Unstructured Meshes

Ion Solvation Thermodynamics from Simulation with a Polarizable Force Field

Electrostatic forces on charged surfaces of bilayer lipid membranes

A new method to calculate broadband dielectric spectra of solvents from molecular dynamics simulations demonstrated with polarizable force fields

Electrostatic interactions in charged nanoslits within an explicit solvent theory