Achieving Energy Conservation in Poisson-Boltzmann Molecular Dynamics: Accuracy and Precision with Finite-Difference...

Area of Science:

• Computational chemistry
• Biophysics
• Applied mathematics

Background:

• Poisson-Boltzmann molecular dynamics are crucial for biomolecular simulations.
• Current numerical methods suffer from limited accuracy and precision, violating energy conservation.
• This violation is a major bottleneck for widespread adoption in simulations.

Purpose of the Study:

• To improve the accuracy and precision of finite-difference methods in Poisson-Boltzmann molecular dynamics.
• To address the violation of energy conservation in numerical simulations.
• To enhance the applicability of Poisson-Boltzmann methods in biomolecular simulations.

Main Methods:

• Explored enforcing interface conditions using the immersed interface method.
• Investigated removing charge singularity to improve finite-difference methods.
• Analyzed the impact of these methods on an analytical test system.

Main Results:

• Significant improvements observed in both energy and force calculations.
• Demonstrated enhanced accuracy and precision of the modified numerical methods.
• Identified the critical need for more accurate boundary force calculations.

Conclusions:

• Enforcing interface conditions and removing charge singularity are effective strategies to improve Poisson-Boltzmann molecular dynamics.
• The developed methods show promise for wider adoption in biomolecular simulations.
• Further research into accurate boundary force calculation is essential for advancing the field.