Induced dipoles extended Lagrangian method

An important advance in making explicit polarizable force fields computationally feasible for MD simulation was the development of the extended Lagrangian methods. This extended dynamics approach was first proposed by Sprik and Klein [91], in the sipirit of the work of Car and Parrinello for ab initio MD dynamics [168], A similar extended system was proposed by van Belle et al. for inducible point dipoles [90, 169], In this approach each dipole is treated as a dynamical variable in the MD simulation and given a mass, Mm, and velocity, p.. The dipoles thus have a kinetic energy, JT (A)2/2, and are propagated using the equations of motion just like the atomic coordinates [90, 91, 170, 171]. The equation of motion for the dipoles is... 

Although a direct comparison between the iterative and the extended Lagrangian methods has not been published, the two methods are inferred to have comparable computational speeds based on indirect evidence. The extended Lagrangian method was found to be approximately 20 times faster than the standard matrix inversion procedure [117] and according to the calculation of Bernardo et al. [208] using different polarizable water potentials, the iterative method is roughly 17 times faster than direct matrix inversion to achieve a convergence of 1.0 x 10-8 D in the induced dipole. 

The extended Lagrangian method is the base for the development of two modified version of SPC and TIP4P, named SPCfq and TIP4Pfq, where the charge itself instead of the induced dipole, is treated as a dynamical variable [100], while the molecular geometry and the number of sites is the same as in the parent functions. 

---