Inertial extended-Lagrangian scheme for solving charge equilibration models.

The inertial extended Lagrangian/self-consistent field scheme (iEL-SCF) has been adopted for solving charge equilibration in LAMMPS as part of the reactive force field ReaxFF, which due to the charge conservation constraint requires solving two sets of linear system of equations for the new charges at each molecular dynamics time-step. Therefore, the extended Lagrangian for charge equilibration is comprised of two auxiliary variables for the intermediate charges which serve as an initial guess for the real charges. We show that the iEL-SCF is able to reduce the number of SCF cycles by 50-80% of the original conjugate gradient self-consistent field solver as tested across diverse systems including water, ferric hydroxide, nitramine RDX, and hexanitrostilbene.

[1]  Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration. , 2019, Journal of chemical theory and computation.

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

[3]  A. V. van Duin,et al.  Shock waves in high-energy materials: the initial chemical events in nitramine RDX. , 2003, Physical review letters.

[4]  V. S. Bryantsev,et al.  Development of a ReaxFF reactive force field for aqueous chloride and copper chloride. , 2010, The journal of physical chemistry. A.

[5]  T. Head-Gordon,et al.  A New Method for Treating Drude Polarization in Classical Molecular Simulation. , 2017, Journal of chemical theory and computation.

[6]  Johann Gasteiger,et al.  Electronegativity equalization: application and parametrization , 1985 .

[7]  V. Vitale,et al.  Performance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory. , 2017, The Journal of chemical physics.

[8]  Brad A. Bauer,et al.  Recent applications and developments of charge equilibration force fields for modeling dynamical charges in classical molecular dynamics simulations , 2012, Theoretical Chemistry Accounts.

[9]  V. Ksenofontov,et al.  Density of phonon states in superconducting FeSe as a function of temperature and pressure , 2010 .

[10]  Jonathan W Essex,et al.  Advanced Potential Energy Surfaces for Molecular Simulation. , 2016, The journal of physical chemistry. B.

[11]  Teresa Head-Gordon,et al.  Accurate Classical Polarization Solution with No Self-Consistent Field Iterations. , 2017, The journal of physical chemistry letters.

[12]  Ananth Grama,et al.  Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques , 2012, Parallel Comput..

[13]  Masoud Aryanpour,et al.  Development of a reactive force field for iron-oxyhydroxide systems. , 2010, The journal of physical chemistry. A.

[14]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .

[15]  L. Tjeng,et al.  Insulating state and the importance of the spin-orbit coupling in Ca3CoRhO6 , 2007, 0705.4538.

[16]  A. Chaffee,et al.  Charge Equilibration Based on Atomic Ionization in Metal–Organic Frameworks , 2015 .

[17]  Ananth Grama,et al.  PuReMD-GPU: A reactive molecular dynamics simulation package for GPUs , 2014, J. Comput. Phys..

[18]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[19]  Peter G. Boyd,et al.  Evaluating Charge Equilibration Methods To Generate Electrostatic Fields in Nanoporous Materials , 2018, Journal of chemical theory and computation.

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

[21]  Tzu-Ray Shan,et al.  Variable charge many-body interatomic potentials , 2012 .

[22]  Teresa Head-Gordon,et al.  An efficient and stable hybrid extended Lagrangian/self-consistent field scheme for solving classical mutual induction. , 2015, The Journal of chemical physics.

[23]  Steven J. Stuart,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994 .

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

[25]  Rajiv K. Kalia,et al.  An extended-Lagrangian scheme for charge equilibration in reactive molecular dynamics simulations , 2015, Comput. Phys. Commun..

[26]  Aiichiro Nakano,et al.  Parallel multilevel preconditioned conjugate-gradient approach to variable-charge molecular dynamics , 1997 .

[27]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[28]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[29]  A. Niklasson,et al.  Higher-Order Extended Lagrangian Born-Oppenheimer Molecular Dynamics for Classical Polarizable Models. , 2018, Journal of chemical theory and computation.

[30]  Wei Jiang,et al.  Optimizing the Performance of Reactive Molecular Dynamics Simulations for Multi-Core Architectures , 2017, ArXiv.

[31]  R. Parr,et al.  Absolute hardness: companion parameter to absolute electronegativity , 1983 .

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

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

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

[35]  A. V. Duin,et al.  ReaxFF: A Reactive Force Field for Hydrocarbons , 2001 .

[36]  Tzu-Ray Shan,et al.  Charge-optimized many-body potential for the hafnium/hafnium oxide system , 2010 .

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