Drude particle

Drude particles are model oscillators used to simulate the effects of electronic polarizability in the context of a classical molecular mechanics force field. They are inspired by the Drude model of mobile electrons and are used in the computational study of proteins, nucleic acids, and other biomolecules.

Classical Drude oscillator

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.

Efficiency of simulation

The major computational cost of simulating classical Drude oscillators is the calculation of the local electrostatic field and the repositioning of the Drude particle at each step. Traditionally, this repositioning is done self consistently. This cost can be reduced by assigning a small mass to each Drude particle, applying a Lagrangian transformation ^[2] and evolving the simulation in the generalised coordinates. This method of simulation has been used to create water models incorporating classical Drude oscillators. ^{[3] [4]}

Quantum Drude oscillator

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.^{[ citation needed ]} 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] This representation allows describing the processes of biological ion transport ^[10] , small drug molecules across hydrophobic cell membranes ^[11] and the behavior of proteins in solutions ^[12] .

References

1. Mackerell, Alexander D. (2004). "Empirical force fields for biological macromolecules: Overview and issues". Journal of Computational Chemistry. Wiley. 25 (13): 1584–1604. doi:10.1002/jcc.20082. ISSN   0192-8651. PMID   15264253. S2CID   9162620.
2. Lamoureux, Guillaume; Roux, Benoı̂t (2003-08-08). "Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm". The Journal of Chemical Physics. AIP Publishing. 119 (6): 3025–3039. Bibcode:2003JChPh.119.3025L. doi: 10.1063/1.1589749 . ISSN   0021-9606.
3. Lamoureux, Guillaume; MacKerell, Alexander D.; Roux, Benoı̂t (2003-09-08). "A simple polarizable model of water based on classical Drude oscillators". The Journal of Chemical Physics. AIP Publishing. 119 (10): 5185–5197. Bibcode:2003JChPh.119.5185L. doi:10.1063/1.1598191. ISSN   0021-9606.
4. Lamoureux, Guillaume; Harder, Edward; Vorobyov, Igor V.; Roux, Benoît; MacKerell, Alexander D. (2006). "A polarizable model of water for molecular dynamics simulations of biomolecules". Chemical Physics Letters. Elsevier BV. 418 (1–3): 245–249. Bibcode:2006CPL...418..245L. doi:10.1016/j.cplett.2005.10.135. ISSN   0009-2614.
5. A. Jones, “Quantum Drude Oscillators for Accurate Many-body Intermolecular Forces,” The University of Edinburgh, 2010.
6. Jones, Andrew; Thompson, Andrew; Crain, Jason; Müser, Martin H.; Martyna, Glenn J. (2009-04-27). "Norm-conserving diffusion Monte Carlo method and diagrammatic expansion of interacting Drude oscillators: Application to solid xenon". Physical Review B. American Physical Society (APS). 79 (14): 144119. Bibcode:2009PhRvB..79n4119J. doi:10.1103/physrevb.79.144119. ISSN   1098-0121.
7. Jones, A.; Cipcigan, F.; Sokhan, V. P.; Crain, J.; Martyna, G. J. (2013-05-31). "Electronically Coarse-Grained Model for Water". Physical Review Letters. American Physical Society (APS). 110 (22): 227801. Bibcode:2013PhRvL.110v7801J. doi:10.1103/physrevlett.110.227801. ISSN   0031-9007. PMID   23767748.
8. "Many body dispersion".
9. Bučko, Tomáš; Lebègue, Sébastien; Gould, Tim; Ángyán, János G (2016-01-12). "Many-body dispersion corrections for periodic systems: an efficient reciprocal space implementation". Journal of Physics: Condensed Matter. IOP Publishing. 28 (4): 045201. Bibcode:2016JPCM...28d5201B. doi:10.1088/0953-8984/28/4/045201. ISSN   0953-8984. PMID   26753609. S2CID   2620743.
10. Manin, Nikolai; da Silva, Mauricio C.; Zdravkovic, Igor; Eliseeva, Olga; Dyshin, Alexey; Yaşar, Orhan; Salahub, Dennis R.; Kolker, Arkadiy M.; Kiselev, Michael G.; Noskov, Sergei Yu. (2016). "LiCl solvation in N-methyl-acetamide (NMA) as a model for understanding Li + binding to an amide plane". Physical Chemistry Chemical Physics. 18 (5): 4191–4200. doi:10.1039/C5CP04847H. ISSN   1463-9076.
11. Lemkul, Justin A.; Huang, Jing; Roux, Benoît; MacKerell, Alexander D. (2016-05-11). "An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications". Chemical Reviews. 116 (9): 4983–5013. doi:10.1021/acs.chemrev.5b00505. ISSN   0009-2665. PMC   4865892 . PMID   26815602.
12. Huang, Jing; Lopes, Pedro E. M.; Roux, Benoît; MacKerell, Alexander D. (2014-09-18). "Recent Advances in Polarizable Force Fields for Macromolecules: Microsecond Simulations of Proteins Using the Classical Drude Oscillator Model". The Journal of Physical Chemistry Letters. 5 (18): 3144–3150. doi:10.1021/jz501315h. ISSN   1948-7185. PMC   4167036 . PMID   25247054.