Automated Algorithms for Quantum-Level Accuracy in Atomistic Simulations: LDRD Final Report.

This report summarizes the result of LDRD project 12-0395, titled “Automated Algorithms for Quantum-level Accuracy in Atomistic Simulations.” During the course of this LDRD, we have developed an interatomic potential for solids and liquids called Spectral Neighbor Analysis Potential (SNAP). The SNAP potential has a very general form and uses machine-learning techniques to reproduce the energies, forces, and stress tensors of a large set of small configurations of atoms, which are obtained using high-accuracy quantum electronic structure (QM) calculations. The local environment of each atom is characterized by a set of bispectrum components of the local neighbor density projected on to a basis of hyperspherical harmonics in four dimensions. The SNAP coefficients are determined using weighted least-squares linear regression against the full QM training set. This allows the SNAP potential to be fit in a robust, automated manner to large QM data sets using many bispectrum components. The calculation of the bispectrum components and the SNAP potential are implemented in the LAMMPS parallel molecular dynamics code. Global optimization methods in the DAKOTA software package are used to seek out good choices of hyperparameters that define the overall structure of the SNAP potential. FitSnap.py, a Python-based software package interfacing to both LAMMPS and DAKOTA is used to formulate the linear regression problem, solve it, and analyze the accuracy of the resultant SNAP potential. We describe a SNAP potential for tantalum that accurately reproduces a variety of solid and liquid properties. Most significantly, in contrast to existing tantalum potentials, SNAP correctly predicts the Peierls barrier for screw dislocation motion. We also present results from SNAP potentials generated for indium phosphide (InP) and silica (SiO2). We describe efficient algorithms for calculating SNAP forces and energies in molecular dynamics simulations using massively parallel computers and advanced processor architectures. Finally, we briefly describe the MSM method for efficient calculation of electrostatic interactions on massively parallel computers.

[1]  Laura Painton Swiler,et al.  SNAP: Automated Generation of Quantum-Accurate Interatomic Potentials. , 2014 .

[2]  Simon D. Hammond,et al.  SNAP: Strong Scaling High Fidelity Molecular Dynamics Simulations on Leadership-Class Computing Platforms , 2014, ISC.

[3]  Stan G Moore,et al.  Extension and evaluation of the multilevel summation method for fast long-range electrostatics calculations. , 2014, The Journal of chemical physics.

[4]  E. M. Lifshitz,et al.  Course in Theoretical Physics , 2013 .

[5]  Kurt B. Ferreira,et al.  Evaluating the feasibility of using memory content similarity to improve system resilience , 2013, ROSS '13.

[6]  C. Weinberger,et al.  Peierls potential of screw dislocations in bcc transition metals: Predictions from density functional theory , 2013 .

[7]  R. Kondor,et al.  On representing chemical environments , 2012, 1209.3140.

[8]  J. Chelikowsky,et al.  Simulating liquid and amorphous silicon dioxide using real-space pseudopotentials , 2012 .

[9]  A. Thompson,et al.  Computational aspects of many-body potentials , 2012 .

[10]  P. Schultz.  Simple intrinsic defects in InP : numerical predictions. , 2012 .

[11]  Nongnuch Artrith,et al.  High-dimensional neural network potentials for metal surfaces: A prototype study for copper , 2012 .

[12]  A. Bartók.  Gaussian Approximation Potential: an interatomic potential derived from first principles Quantum Mechanics , 2010, 1003.2817.

[13]  Makoto Taiji,et al.  42 TFlops hierarchical N-body simulations on GPUs with applications in both astrophysics and turbulence , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[14]  R. Kondor,et al.  Gaussian approximation potentials: the accuracy of quantum mechanics, without the electrons. , 2009, Physical review letters.

[15]  T. Germann,et al.  369 Tflop/s molecular dynamics simulations on the Roadrunner general-purpose heterogeneous supercomputer , 2008, 2008 SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

[16]  Erich Strohmaier,et al.  Linearly scaling 3D fragment method for large-scale electronic structure calculations , 2008, 2008 SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

[17]  V. Vítek,et al.  Multiscale modeling of plastic deformation of molybdenum and tungsten: I. Atomistic studies of the core structure and glide of 1/2〈1 1 1〉 screw dislocations at 0 K , 2008, 0807.2667.

[18]  A. Meremianin.  Hyperspherical harmonics with arbitrary arguments , 2008, 0807.2128.

[19]  John A. Gunnels,et al.  Extending stability beyond CPU millennium: a micron-scale atomistic simulation of Kelvin-Helmholtz instability , 2007, Proceedings of the 2007 ACM/IEEE Conference on Supercomputing (SC '07).

[20]  Peter S. Lomdahl,et al.  MOLECULAR DYNAMICS COMES OF AGE: 320 BILLION ATOM SIMULATION ON BlueGene/L , 2006 .

[21]  Y. Mishin,et al.  Angular-dependent interatomic potential for tantalum , 2006 .

[22]  Haydn N. G. Wadley,et al.  Analytic bond-order potential for the gallium arsenide system , 2006 .

[23]  A. Meremianin.  Multipole expansions in four-dimensional hyperspherical harmonics , 2005, math-ph/0510080.

[24]  R. Johnson,et al.  Misfit-energy-increasing dislocations in vapor-deposited CoFe/NiFe multilayers , 2004 .

[25]  Xiang-Yang Liu,et al.  Embedded-atom-method tantalum potential developed by the force-matching method , 2003 .

[26]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .

[27]  G. Kresse,et al.  From ultrasoft pseudopotentials to the projector augmented-wave method , 1999 .

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

[29]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

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

[31]  P. Steinhardt,et al.  Bond-orientational order in liquids and glasses , 1983 .

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

[33]  H. Monkhorst,et al.  SPECIAL POINTS FOR BRILLOUIN-ZONE INTEGRATIONS , 1976 .

[34]  E. Fisher,et al.  Pressure derivatives of the elastic moduli of niobium and tantalum , 1976 .

[35]  F. Featherston,et al.  ELASTIC CONSTANTS OF TANTALUM, TUNGSTEN, AND MOLYBDENUM , 1963 .

[36]  Audra E. Kosh,et al.  Linear Algebra and its Applications , 1992 .

[37]  D. Mattis.  Quantum Theory of Angular Momentum , 1981 .