Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm

A simple treatment for incorporating induced polarization in computer simulations is formulated on the basis of the classical Drude oscillator model. In this model, electronic induction is represented by the displacement of a charge-carrying massless particle attached to a polarizable atom under the influence of the local electric field. The traditional self-consistent field (SCF) regime of induced polarization is reproduced if these auxiliary particles are allowed to relax instantaneously to their local energy minima for any given fixed configuration of the atoms in the system. In practice, such treatment is computationally prohibitive for generating molecular dynamics trajectories because the electric field…

• NS
  National Science Foundation