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.

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

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

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

Incorporating Dipolar Solvents with Variable Density in Poisson-Boltzmann Electrostatics

A Finite Element Solution of Lateral Periodic Poisson–Boltzmann Model for Membrane Channel Proteins

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

Calculations Of Free Energies Based On The Combination Of Molecular Dynamics Simulations And Continuum Solvation Model

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

Numerical Simulation of Hydrodynamic Focusing of Particles in Straight Channel Flows with the Immersed Boundary-Lattice Boltzmann Method

An Approach to Determine the Solid Volume Fraction Near Solid Boundaries for Partially Saturated Method Within the Lattice Boltzmann Method

Force and Heat Current Formulas for Many-Body Potentials in Molecular Dynamics Simulations with Applications to Thermal Conductivity Calculations

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

The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients

Poisson-Boltzmann based machine learning (PBML) model for electrostatic analysis

Some challenges of diffused interfaces in implicit-solvent models

Mathematical Analysis of the BIBEE Approximation for Molecular Solvation: Exact Results for Spherical Inclusions

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

Simulating dynamics of ellipsoidal particles using lattice Boltzmann method

An Improved Immersed Moving Boundary for Hydrodynamic Force Calculation in Lattice Boltzmann Method

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