Computation of electrostatic forces between solvated molecules determined by the Poisson-Boltzmann equation using a...

Area of Science:

• Computational chemistry
• Molecular dynamics
• Physical chemistry

Background:

• Calculating electrostatic forces is crucial for understanding molecular interactions in solution.
• Existing methods face challenges with hypersingularity in boundary element methods (BEM).

Purpose of the Study:

• To develop a rigorous and efficient approach for computing electrostatic forces among multiple solvated molecules in ionic solutions.
• To overcome limitations of traditional BEM in calculating electrostatic interactions.

Main Methods:

• Utilizing the linearized Poisson-Boltzmann equation.
• Implementing a variational principle within the boundary element method (BEM).
• Avoiding Maxwell stress tensor calculations on molecular surfaces to prevent hypersingularity.

Main Results:

• The proposed BEM approach accurately calculates intermolecular electrostatic interaction energy and force.
• The method demonstrates reliability and efficiency in test cases.
• It successfully avoids the hypersingularity problem inherent in direct BEM calculations.

Conclusions:

• This novel method offers an accurate and efficient way to compute electrostatic forces for solvated molecules.
• The approach has potential applications in Brownian dynamics simulations of biomolecular association.
• It provides a robust alternative to existing methods for electrostatic interaction calculations.