Communication: Constrained molecular dynamics for polarizable models.

Area of Science:

• Computational Physics
• Materials Science
• Theoretical Chemistry

Background:

• Accurate simulation of polarizable systems is crucial for understanding material properties.
• Existing methods for empirical shell models often lack efficiency, stability, or time reversibility.
• The null gradient condition is a key requirement for modeling shell degrees of freedom.

Purpose of the Study:

• To introduce a novel algorithm for the numerical evolution of empirical shell models of polarizable systems.
• To address limitations in current methods by developing a more robust and accurate simulation technique.
• To demonstrate the algorithm's applicability and advantages through a proof-of-principle calculation.

Main Methods:

• Employs constrained molecular dynamics to enforce the null gradient condition at each time step.
• The algorithm is designed to be symplectic and time reversible.
• A conjugate-gradient procedure is used in the current method for comparison.

Main Results:

• The new algorithm is demonstrated to be efficient and stable.
• It successfully satisfies the null gradient condition, ensuring accurate modeling of polarizable systems.
• Comparison with the conjugate-gradient method shows superior properties, including symplecticity and time reversibility.

Conclusions:

• The proposed algorithm provides a significant advancement for simulating polarizable systems using empirical shell models.
• Its efficiency, stability, and time-reversible nature make it a superior alternative to existing methods.
• The algorithm's principles are broadly applicable to other models requiring minimum conditions on driven dynamical variables.