Drude particle

Most force fields in current practice represent individual atoms as point particles interacting according to the laws of Newtonian mechanics. To each atom, a single electric charge is assigned that doesn't change during the course of the simulation. However, such models cannot have induced dipoles or other electronic effects due to a changing local environment.

Classical Drude particles are massless virtual sites carrying a partial electric charge, attached to individual atoms via a harmonic spring. The spring constant and relative partial charges on the atom and associated Drude particle determine its response to the local electrostatic field, serving as a proxy[1] for the changing distribution of the electronic charge of the atom or molecule. However, this response is limited to a changing dipole moment. This response is not enough to model interactions in environments with large field gradients, which interact with higher order moments.

Since the response of a classical Drude oscillator is limited, it is not enough to model interactions in heterogeneous media with large field gradients, where higher order electronic responses have significant contributions to the interaction energy. A quantum Drude oscillator (QDO)[5][6][7] is a natural extension to the classical Drude oscillator. Instead of a classical point particle serving as a proxy for the charge distribution, a QDO uses a quantum harmonic oscillator, in the form of a pseudoelectron connected to an oppositely charged pseudonucleus by a harmonic spring.

A QDO has three free parameters: the spring's frequency ${\displaystyle \omega }$, the pseudoelectron's charge ${\displaystyle q}$ and the system's reduced mass ${\displaystyle \mu }$. The ground state of a QDO is a gaussian of width ${\displaystyle \sigma =1/{\sqrt {2\mu \omega }}}$. Adding an external field perturbs the ground state of a QDO, which allows us to calculate its polarizability.[5] To second order, the change in energy relative to the ground state is given by the following series:

${\displaystyle E^{(2)}=\sum _{l=1}^{\infty }E_{l}^{(2)}=\sum _{l=1}^{\infty }-{\frac {Q^{2}\alpha _{l}}{2R^{2l+2}}}}$

where the polarizabilities ${\displaystyle \alpha _{l}}$ are

${\displaystyle \alpha _{l}=\left[{\frac {q^{2}}{\mu \omega ^{2}}}\right]\left[{\frac {(2l-1)!!}{l}}\right]\left({\frac {\hbar }{2\mu \omega }}\right)^{l-1}}$

Furthermore, since QDOs are quantum mechanical objects, their electrons can correlate, giving rise to dispersion forces between them. The second order change in energy corresponding to such an interaction is:

${\displaystyle E^{(2)}=\sum _{l=3}^{\infty }C_{2l}R^{-2l}}$

with the first three dispersion coefficients being (in the case of identical QDOs):

${\displaystyle C_{6}={\frac {3}{4}}\alpha _{1}\alpha _{1}\hbar \omega }$

${\displaystyle C_{8}=5\alpha _{1}\alpha _{2}\hbar \omega }$

${\displaystyle C_{10}=\left({\frac {21}{2}}\alpha _{1}\alpha _{3}+{\frac {70}{4}}\alpha _{2}\alpha _{2}\right)\hbar \omega }$

Since the response coefficients of QDOs depend on three parameters only, they are all related. Thus, these response coefficients can combine into four dimensionless constants, all equal to unity:

${\displaystyle {\sqrt {\frac {20}{9}}}{\frac {\alpha _{2}}{\sqrt {\alpha _{1}\alpha _{3}}}}=1}$

${\displaystyle {\sqrt {\frac {49}{40}}}{\frac {C_{8}}{\sqrt {C_{6}C_{10}}}}=1}$

${\displaystyle {\frac {C_{6}\alpha _{1}}{4C_{9}}}=1}$

The QDO representation of atoms is the basis of the many body dispersion model [8] which is a popular way to account for electrostatic forces in molecular dynamics simulations.[9]