Molecular dynamics scheme for precise estimation of electrostatic interaction via zero-dipole summation principle

Area of Science:

• Computational Chemistry
• Statistical Mechanics
• Materials Science

Background:

• Evaluating electrostatic energy in classical particle systems is computationally intensive.
• Simple cutoff truncation methods can introduce artificial nonzero-charge and nonzero-dipole states.
• Existing methods like charge-neutral summation and preaveraged potential have limitations.

Purpose of the Study:

• To propose and develop a novel method, zero-dipole summation, for accurate electrostatic energy evaluation in molecular dynamics.
• To create an algorithm for efficient implementation in molecular dynamics simulations.
• To analyze the theoretical underpinnings and error structure of the new method.

Main Methods:

• Developed a zero-dipole summation algorithm for electrostatic energy calculation.
• Employed both heuristic and axiomatic approaches for method derivation and consistency.
• Conducted theoretical analysis of error structure and compared with existing methods.
• Performed numerical simulations using liquid sodium chloride.

Main Results:

• The zero-dipole summation method prevents artificial states from cutoff truncation.
• The energy formula is a simple pairwise function sum, suitable for high-performance computing.
• Numerical simulations confirmed sufficient accuracy and stability with a small damping factor.
• The method demonstrated stable numerical integration in liquid molecular dynamics simulations.

Conclusions:

• Zero-dipole summation offers an accurate and computationally efficient approach for electrostatic energy calculations.
• The method extends charge-neutral summation and generalizes the preaveraged potential method.
• The findings suggest potential for further applications and development in computational physics and chemistry.