Polarizable molecular dynamics in a polarizable continuum solvent.

Area of Science:

• Computational Chemistry
• Molecular Dynamics Simulations
• Quantum Chemistry

Background:

• Accurate molecular simulations require realistic solvent models.
• Polarizable continuum models (PCMs) are essential for capturing solvent effects.
• Integrating PCMs with polarizable force fields presents computational challenges.

Purpose of the Study:

• To develop scalable polarizable molecular dynamics (MD) simulations.
• To implement a polarizable continuum solvent model with molecular shape cavities.
• To achieve an exact solution for mutual polarization in simulations.

Main Methods:

• Domain decomposition Conductor-like Screening Model (ddCOSMO) for efficient continuum calculations.
• Variational formulation for rigorous coupling of continuum and polarizable force fields.
• Development of a scalable parallel implementation in Tinker-HP.

Main Results:

• Stable NVE simulations with long-term energy conservation were achieved.
• The method allows for simulations of real-life systems on standard cluster nodes.
• Successful coupling of ddCOSMO with nonvariational force fields (e.g., AMOEBA) was demonstrated.

Conclusions:

• The developed methodology enables efficient and accurate polarizable MD simulations.
• The ddCOSMO approach offers a robust solution for integrating continuum solvents with polarizable force fields.
• This technique opens new possibilities for studying complex chemical systems.