Linear scaling computation of forces for the domain-decomposition linear Poisson-Boltzmann method

Area of Science:

• Computational Chemistry
• Biochemistry
• Physical Chemistry

Background:

• The Linearized Poisson-Boltzmann (LPB) equation is a standard model for solvent effects in computational chemistry.
• Accurate calculation of forces is essential for molecular simulations.

Purpose of the Study:

• To derive analytical forces using a domain-decomposition-based LPB method.
• To develop an efficient computational strategy for these forces.

Main Methods:

• Derivation of analytical forces with van-der Waals or solvent-accessible surface boundary conditions.
• Implementation using domain decomposition and the fast multipole method for linear scaling.
• Numerical tests to validate accuracy and efficiency.

Main Results:

• Successful derivation of analytical forces for the LPB equation.
• An efficient implementation achieving linear scaling with the number of atoms.
• Demonstrated accuracy and improved efficiency compared to other methods.

Conclusions:

• The developed method provides an accurate and efficient way to compute analytical forces within the LPB framework.
• This advancement is valuable for large-scale molecular simulations in chemistry and biochemistry.
• The linear scaling approach significantly enhances computational feasibility.