An efficient and stable hybrid extended Lagrangian/self-consistent field scheme for solving classical mutual induction.

We have adapted a hybrid extended Lagrangian self-consistent field (EL/SCF) approach, developed for time reversible Born Oppenheimer molecular dynamics for quantum electronic degrees of freedom, to the problem of classical polarization. In this context, the initial guess for the mutual induction calculation is treated by auxiliary induced dipole variables evolved via a time-reversible velocity Verlet scheme. However, we find numerical instability, which is manifested as an accumulation in the auxiliary velocity variables, that in turn results in an unacceptable increase in the number of SCF cycles to meet even loose convergence tolerances for the real induced dipoles over the course of a 1 ns trajectory of the AMOEBA14 water model. By diagnosing the numerical instability as a problem of resonances that corrupt the dynamics, we introduce a simple thermostating scheme, illustrated using Berendsen weak coupling and Nose-Hoover chain thermostats, applied to the auxiliary dipole velocities. We find that the inertial EL/SCF (iEL/SCF) method provides superior energy conservation with less stringent convergence thresholds and a correspondingly small number of SCF cycles, to reproduce all properties of the polarization model in the NVT and NVE ensembles accurately. Our iEL/SCF approach is a clear improvement over standard SCF approaches to classical mutual induction calculations and would be worth investigating for application to ab initio molecular dynamics as well.

[1]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[2]  Shoshana J. Wodak,et al.  Molecular dynamics simulation of polarizable water by an extended Lagrangian method , 1992 .

[3]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[4]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[5]  Celeste Sagui,et al.  Towards an accurate representation of electrostatics in classical force fields: efficient implementation of multipolar interactions in biomolecular simulations. , 2004, The Journal of chemical physics.

[6]  B. Thole.  Molecular polarizabilities calculated with a modified dipole interaction , 1981 .

[7]  Anders M N Niklasson,et al.  Generalized extended Lagrangian Born-Oppenheimer molecular dynamics. , 2014, The Journal of chemical physics.

[8]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[9]  Michele Ceriotti,et al.  Efficient multiple time scale molecular dynamics: Using colored noise thermostats to stabilize resonances. , 2010, The Journal of chemical physics.

[10]  R. McLachlan,et al.  The accuracy of symplectic integrators , 1992 .

[11]  Teresa Head-Gordon,et al.  Advanced potential energy surfaces for condensed phase simulation. , 2014, Annual review of physical chemistry.

[12]  Matt Challacombe,et al.  Time-reversible Born-Oppenheimer molecular dynamics. , 2006, Physical review letters.

[13]  Anders M N Niklasson,et al.  Extended Born-Oppenheimer molecular dynamics. , 2008, Physical review letters.

[14]  Wei Wang,et al.  Fast evaluation of polarizable forces. , 2005, The Journal of chemical physics.

[15]  M. Tuckerman,et al.  Long time molecular dynamics for enhanced conformational sampling in biomolecular systems. , 2004, Physical review letters.

[16]  Richard A Friesner,et al.  Efficient Simulation Method for Polarizable Protein Force Fields:  Application to the Simulation of BPTI in Liquid Water. , 2005, Journal of chemical theory and computation.

[17]  Peter Pulay,et al.  Geometry optimization by direct inversion in the iterative subspace , 1984 .

[18]  Nicolas Bock,et al.  Extended Lagrangian Born-Oppenheimer molecular dynamics with dissipation. , 2009, The Journal of chemical physics.

[19]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[20]  Nohad Gresh,et al.  Inclusion of the ligand field contribution in a polarizable molecular mechanics: SIBFA‐LF , 2003, J. Comput. Chem..

[21]  Margaret E. Johnson,et al.  Current status of the AMOEBA polarizable force field. , 2010, The journal of physical chemistry. B.

[22]  Benjamin Stamm,et al.  Scalable Evaluation of Polarization Energy and Associated Forces in Polarizable Molecular Dynamics: I. Toward Massively Parallel Direct Space Computations , 2014 .

[23]  M. Challacombe,et al.  Time-reversible ab initio molecular dynamics. , 2007, The Journal of chemical physics.

[24]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[25]  P. T. V. Duijnen,et al.  Molecular and Atomic Polarizabilities: Thole's Model Revisited , 1998 .

[26]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[27]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[28]  R. Skeel,et al.  Nonlinear Resonance Artifacts in Molecular Dynamics Simulations , 1998 .

[29]  Jay W Ponder,et al.  Revised Parameters for the AMOEBA Polarizable Atomic Multipole Water Model. , 2015, The journal of physical chemistry. B.

[30]  Pedro E. M. Lopes,et al.  Molecular modeling and dynamics studies with explicit inclusion of electronic polarizability: theory and applications , 2009, Theoretical chemistry accounts.

[31]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..