High-Performance Scalable Molecular Dynamics Simulations of a Polarizable Force Field Based on Classical Drude Oscillators in NAMD

Citations Over TimeTop 10% of 2010 papers

Abstract

Incorporating the influence of induced polarization in large-scale atomistic molecular dynamics (MD) simulations is a critical challenge in the progress toward computations of increased accuracy. One computationally efficient treatment is based on the classical Drude oscillator in which an auxiliary charged particle is attached by a spring to each nucleus. Here, we report the first implementation of this model in the program NAMD. An extended Lagrangian dynamics with a dual-Langevin thermostat scheme applied to the Drude-nucleus pairs is employed to efficiently generate classical dynamic propagation near the self-consistent field limit. Large-scale MD simulations based on the Drude polarizable force field scale very well on massively distributed supercomputing platforms, the computational demand being only about 50-100% higher than for nonpolarizable models. As an illustration, a large-scale 150 mM NaCl aqueous salt solution is simulated, and the calculated ionic conductivity is shown to be in excellent agreement with experiment.

Related Papers

• → Medium perturbations on the molecular polarizability calculated within a localized dipole interaction model(2002)39 cited
• → A quest for effective polarizability as a function of the radii(2021)1 cited
• → FREQUENCY-DEPENDENT POLARIZABILITY OF SMALL SILICON CLUSTERS(2011)1 cited
• Simulation Study on Dissipative Particle Dynamics of The Self-assembly of Di-block Copolymer in Body(2015)
• Study of Self-Assembly in Block Copolymer (BCP)-Homopolymer Blends using Dissipative Particle Dynamics (DPD) Simulations(2017)