Optimizing the Accuracy and Efficiency of Fast Hierarchical Multipole Expansions for MD Simulations

Area of Science:

• Computational chemistry
• Molecular dynamics
• Electrostatics

Background:

• Accurate calculation of electrostatic interactions is crucial for molecular simulations.
• Traditional methods struggle with the computational cost of large systems.

Purpose of the Study:

• To develop fast and accurate methods for calculating electrostatic interactions in macromolecular systems.
• To implement and validate the Scalable Adaptive Multipole Method (SAMMp) algorithms.

Main Methods:

• Hierarchical decomposition of systems into adaptive clusters.
• p'th order Cartesian Taylor expansions of Coulomb interactions.
• Formulation using symmetric, traceless tensors for multipole moments.
• Application of reaction field (RF) correction for periodic systems.

Main Results:

• SAMMp algorithms provide rapid and accurate electrostatic calculations.
• The methods guarantee the reaction principle.
• Simulations show stable and low-noise performance of SAMMp/RF for water and ionic solutions.

Conclusions:

• SAMMp algorithms offer an efficient approach for electrostatic interactions in complex molecular systems.
• The integration with reaction field correction ensures accuracy for periodic boundary conditions.
• Validated for molecular dynamics simulations, demonstrating practical applicability.