Atomic forces for geometry-dependent point multipole and gaussian multipole models.

Area of Science:

• Computational Chemistry
• Molecular Modeling
• Quantum Mechanics

Background:

• Standard atomic multipole models often assume rigid molecules, neglecting forces arising from flexibility.
• Flexible molecules introduce additional atomic forces due to torque transfer and geometry-dependent multipole moments.

Purpose of the Study:

• To derive and present atomic force expressions for geometry-dependent multipoles in flexible molecule simulations.
• To enhance the accuracy of molecular simulations by accounting for atomic-level forces in flexible systems.

Main Methods:

• Developed new general expressions for Wigner function derivatives.
• Derived atomic force equations applicable to point multipoles and Gaussian multipole charge densities.
• Tested models using hydrogen-bonded dimers, comparing electrostatic energies and forces with ab initio calculations.

Main Results:

• Geometry-dependent multipole models are necessary for accurately reproducing ab initio atomic forces.
• Both static and geometry-dependent models can reproduce total molecular forces and torques.
• The derived force equations are suitable for simulations of flexible molecules using atomic multipoles.

Conclusions:

• The new atomic force expressions enable more accurate simulations of flexible molecules.
• This work advances the development of next-generation force fields incorporating geometry-dependent multipole models.
• Accurate atomic force calculations are crucial for understanding molecular behavior in flexible systems.