Stable and Efficient Linear Scaling First-Principles Molecular Dynamics for 10000+ Atoms.

The recent progress of linear-scaling or O(N) methods in density functional theory (DFT) is remarkable. Given this, we might expect that first-principles molecular dynamics (FPMD) simulations based on DFT could treat more realistic and complex systems using the O(N) technique. However, very few examples of O(N) FPMD simulations exist to date, and information on the accuracy and reliability of the simulations is very limited. In this paper, we show that efficient and robust O(N) FPMD simulations are now possible by the combination of the extended Lagrangian Born-Oppenheimer molecular dynamics method, which was recently proposed by Niklasson ( Phys. Rev. Lett. 2008 , 100 , 123004 ), and the density matrix method as an O(N) technique. Using our linear-scaling DFT code Conquest, we investigate the reliable calculation conditions for accurate O(N) FPMD and demonstrate that we are now able to do practical, reliable self-consistent FPMD simulations of a very large system containing 32768 atoms.

[1]  David R Bowler,et al.  Non-self-consistent Density-Functional Theory Exchange-Correlation Forces for GGA Functionals. , 2009, Journal of chemical theory and computation.

[2]  Li,et al.  Density-matrix electronic-structure method with linear system-size scaling. , 1993, Physical review. B, Condensed matter.

[3]  Taisuke Ozaki,et al.  A method of orbital analysis for large-scale first-principles simulations. , 2014, The Journal of chemical physics.

[4]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[5]  Priya Vashishta,et al.  A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations. , 2014, The Journal of chemical physics.

[6]  R. Mcweeny.  Some Recent Advances in Density Matrix Theory , 1960 .

[7]  Tsuyoshi Miyazaki,et al.  Large-scale DFT simulations with a linear-scaling DFT code CONQUEST on K-computer , 2014 .

[8]  V. Halpern.  Inhomogeneities as a result of non-exponential relaxation in disordered systems , 2002 .

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

[10]  David R. Bowler,et al.  Recent progress with large‐scale ab initio calculations: the CONQUEST code , 2006 .

[11]  D. Bowler,et al.  O(N) methods in electronic structure calculations. , 2011, Reports on progress in physics. Physical Society.

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

[13]  D R Bowler,et al.  Atomic force algorithms in density functional theory electronic-structure techniques based on local orbitals. , 2004, The Journal of chemical physics.

[14]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[15]  David R. Bowler,et al.  Pseudo-atomic orbitals as basis sets for the O(N) DFT code CONQUEST , 2008 .

[16]  Anders M N Niklasson,et al.  First principles molecular dynamics without self-consistent field optimization. , 2013, The Journal of chemical physics.

[17]  D R Bowler,et al.  Calculations for millions of atoms with density functional theory: linear scaling shows its potential , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[18]  Taisuke Ozaki,et al.  Large-scale first-principles molecular dynamics for electrochemical systems with O(N) methods. , 2012, The Journal of chemical physics.

[19]  C. M. Goringe,et al.  Basis functions for linear-scaling first-principles calculations , 1997 .

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

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

[22]  Hernández,et al.  Linear-scaling density-functional-theory technique: The density-matrix approach. , 1996, Physical review. B, Condensed matter.

[23]  David R. Bowler,et al.  Density matrices in O(N) electronic structure calculations: theory and applications , 1998 .

[24]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[25]  Anders M N Niklasson,et al.  Energy conserving, linear scaling Born-Oppenheimer molecular dynamics. , 2012, The Journal of chemical physics.

[26]  M. Robinson,et al.  Accurate ionic forces and geometry optimization in linear-scaling density-functional theory with local orbitals , 2011, 1103.5869.

[27]  N. Bock,et al.  Extended Lagrangian free energy molecular dynamics. , 2011, The Journal of chemical physics.

[28]  David R. Bowler,et al.  Recent progress in linear scaling ab initio electronic structure techniques , 2002 .

[29]  Eiji Tsuchida,et al.  Ab initio molecular dynamics simulations with linear scaling: application to liquid ethanol , 2008 .

[30]  Jean-Luc Fattebert,et al.  Accurate and scalable O(N) algorithm for first-principles molecular-dynamics computations on large parallel computers. , 2014, Physical review letters.

[31]  David Bowler,et al.  Comment on "Accurate and Scalable O(N) Algorithm for First-Principles Molecular-Dynamics Computations on Large Parallel Computers" , 2014 .

[32]  Hernández,et al.  Self-consistent first-principles technique with linear scaling. , 1995, Physical review. B, Condensed matter.

[33]  David E. Manolopoulos,et al.  Canonical purification of the density matrix in electronic-structure theory , 1998 .

[34]  T. Arias,et al.  Iterative minimization techniques for ab initio total energy calculations: molecular dynamics and co , 1992 .