DeepAI AI Chat
Log In Sign Up

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

by   Cleophas Kweyu, et al.

The Poisson-Boltzmann equation (PBE) is an implicit solvent continuum model for calculating the electrostatic potential and energies of ionic solvated biomolecules. However, its numerical solution remains a significant challenge due strong singularities and nonlinearity caused by the singular source terms and the exponential nonlinear terms, respectively. An efficient method for the treatment of singularities in the linear PBE was introduced in <cit.>, that is based on the RS tensor decomposition for both electrostatic potential and the discretized Dirac delta distribution. In this paper, we extend this regularization method to the nonlinear PBE. We apply the PBE only to the regular part of the solution corresponding to the modified right-hand side via extraction of the long-range part in the discretized Dirac delta distribution. The total electrostatic potential is obtained by adding the long-range solution to the directly precomputed short-range potential. The main computational benefit of the approach is the automatic maintaining of the continuity in the Cauchy data on the solute-solvent interface. The boundary conditions are also obtained from the long-range component of the precomputed canonical tensor representation of the Newton kernel. In the numerical experiments, we illustrate the accuracy of the nonlinear regularized PBE (NRPBE) over the classical variant.


page 1

page 2

page 3

page 4


Factorized structure of the long-range two-electron integrals tensor and its application in quantum chemistry

We introduce two new approximation methods for the numerical evaluation ...

Enriched Gradient Recovery for Interface Solutions of the Poisson-Boltzmann Equation

Accurate calculation of electrostatic potential and gradient on the mole...

Integrating factor techniques applied to the Schrödinger-like equations. Comparison with Split-Step methods

The nonlinear Schrödinger and the Schrödinger-Newton equations model man...

Adaptive pseudo-time methods for the Poisson-Boltzmann equation with Eulerian solvent excluded surface

This work further improves the pseudo-transient approach for the Poisson...