Drude particle


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. It is based on the Drude model of mobile electrons. Molecular mechanics models, commonly used for computational calculations such as structural minimization and molecular dynamics simulations, represent individual atoms or other particles as hard spheres that interact according to the laws of Newtonian mechanics. These methods are often used in the computational study of proteins, nucleic acids, and other biomolecules.

Most force fields in current practice use a fixed-charge model in which each atom in the simulation is assigned a single electric charge that does not change during the course of the simulation. This obviously cannot model induced dipoles or other electronic effects of a changing local environment. Drude particles are massless virtual sites carrying partial electric charge and attached to individual atoms via a harmonic spring. The spring constant and relative partial charges on the atom and associated Drude particle determine the extent to which the Drude particle responds to the local electrostatic field. The movement of the Drude particle thus serves as a proxy for the changing distribution of the electric charge associated with the corresponding atom.[1]

The Drude particle method, and polarizable force fields in general, have not yet been widely applied due to their very high computational cost compared to the equivalent fixed-charge simulation. In the Drude model, this cost is largely due to the challenge in recalculating the local electrostatic field and repositioning the Drude particles at each time step; traditionally this is done using an iterative self consistent field (SCF) method, although a much more efficient method has been developed that assigns a very small mass to each Drude particle and applies a Lagrangian transformation.[2] Water models incorporating Drude sites have also been developed[3] and refined for incorporation into a polarizable force field under development.[4]

References

  1. ^ MacKerell AD. (2004). Empirical force fields for biological macromolecules: Overview and issues. J Comp Chem 25(13): 1584-1604.
  2. ^ Lamoureux G, Roux B. (2003). Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm. J Chem Phys 119(6):3025-3039.
  3. ^ Lamoureux G, MacKerell AD, Roux B. (2003). A simple polarizable model of water based on classical Drude oscillators. J Chem Phys 119(3):5185-97.
  4. ^ Lamoureux G, Harder E, Vorobyov IV, Roux B, MacKerell AD. (2006). A polarizable model of water for molecular dynamics simulations of biomolecules. Chem Phys Lett 418:245-9.

Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Force field (chemistry) — In the context of molecular mechanics, a force field (also called a forcefield) refers to the functional form and parameter sets used to describe the potential energy of a system of particles (typically but not necessarily atoms). Force field… …   Wikipedia

  • Free electron model — In three dimensions, the density of states of a gas of fermions is proportional to the square root of the kinetic energy of the particles. In solid state physics, the free electron model is a simple model for the behaviour of valence electrons in …   Wikipedia

  • Atomic orbital — The shapes of the first five atomic orbitals: 1s, 2s, 2px, 2py, and 2pz. The colors show the wave function phase. These are graphs of ψ(x,y,z) functions which depend on the coordinates of one electron. To see the elongated shape of ψ(x,y,z)2… …   Wikipedia

  • Le Sage's theory of gravitation — is the most common name for the kinetic theory of gravity originally proposed by Nicolas Fatio de Duillier in 1690 and later by Georges Louis Le Sage in 1748. The theory proposed a mechanical explanation for Newton s gravitational force in terms… …   Wikipedia

  • Atomic theory — Atomic model redirects here. For the unrelated term in mathematical logic, see Atomic model (mathematical logic). This article is about the historical models of the atom. For a history of the study of how atoms combine to form molecules, see… …   Wikipedia

  • Molecular dynamics — (MD) is a computer simulation of physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a period of time, giving a view of the motion of the atoms. In the most common version, the trajectories of molecules… …   Wikipedia

  • Wiedemann-Franz law — In physics, the Wiedemann Franz law states that the ratio of the thermal conductivity ( K ) to the electrical conductivity ( sigma; ) of a metal is proportional to the temperature ( T ). :frac{K}{sigma}=LTTheoretically, the proportionality… …   Wikipedia

  • Ideal gas — Thermodynamics …   Wikipedia

  • Effective mass — In solid state physics, a particle s effective mass is the mass it seems to carry in the semiclassical model of transport in a crystal. It can be shown that, under most conditions, electrons and holes in a crystal respond to electric and magnetic …   Wikipedia

  • Kinetic inductance — is the manifestation of the inertial mass of mobile charge carriers in alternating electric fields as an equivalent series inductance. Kinetic inductance is observed in high carrier mobility conductors (e.g. superconductors), and at high… …   Wikipedia