skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Spectral neighbor analysis method for automated generation of quantum-accurate interatomic potentials

Journal Article · · Journal of Computational Physics
 [1];  [2];  [3];  [4];  [5]
  1. Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States). Multiscale Science Dept.
  2. Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States). Optimization and Uncertainty Quantification Dept.
  3. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Scalable Algorithms Dept.
  4. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Computational Materials and Data Science Dept.
  5. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States). Computational Materials and Data Science Dept.; Drexel Univ., Philadelphia, PA (United States). Dept. of Materials Science and Engineering

Here, we present a new 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 onto a basis of hyperspherical harmonics in four dimensions. The bispectrum components are the same bond-orientational order parameters employed by the GAP potential [1]. The SNAP potential, unlike GAP, assumes a linear relationship between atom energy and bispectrum components. The linear 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. We demonstrate that a previously unnoticed symmetry property can be exploited to reduce the computational cost of the force calculations by more than one order of magnitude. We present results for a SNAP potential for tantalum, showing that it accurately reproduces a range of commonly calculated properties of both the crystalline solid and the liquid phases. In addition, unlike simpler existing potentials, SNAP correctly predicts the energy barrier for screw dislocation migration in BCC tantalum.

Research Organization:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA)
DOE Contract Number:
AC04-94AL85000
OSTI ID:
1426894
Report Number(s):
SAND2014-17598J; 537437
Journal Information:
Journal of Computational Physics, Vol. 285, Issue C; ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English