---
title: Spectral Neighbor Analysis Potential (SNAP)
url: https://www.emergentmind.com/topics/spectral-neighbor-analysis-potential-snap
type: topic
---

# Spectral Neighbor Analysis Potential (SNAP)

Searching arXiv for SNAP foundational and extension papers to ground the article in current and original sources.
arxiv_search(query="Spectral Neighbor Analysis Potential SNAP bispectrum 2014 Thompson quadratic SNAP Wood Thompson 2017 explicit multi-element SNAP 2020 electrostatic SNAP lithium nitride", max_results=10)
Spectral Neighbor Analysis Potential (SNAP) is a machine-learning interatomic potential in which the total energy is decomposed into local atomic contributions, and each local environment is represented by rotationally invariant bispectrum components derived from a spectral expansion of an atom-centered neighbor density on the 3-sphere. In its original formulation, SNAP was introduced as a short-ranged, linear-scaling potential for solids and liquids that aims to approach quantum-mechanical accuracy while remaining suitable for large-scale molecular dynamics; subsequent work extended the same descriptor machinery to quadratic regressions, explicit multi-element descriptors, electrostatic corrections for ionic systems, and spectral-neighbour representations for vector fields including spin [1409.3880].

## 1. Historical emergence and conceptual role

SNAP was introduced as a new interatomic potential for solids and liquids in order to interpolate between two established limits of atomistic simulation: conventional empirical potentials, which are computationally efficient but often too restrictive, and direct quantum molecular dynamics, which is accurate but too expensive for large systems and long trajectories. Its original motivating use case was tantalum, especially across crystals, defects, liquid structure, and screw-dislocation motion in bcc Ta [1409.3880].

A defining conceptual feature of SNAP is that it separates **descriptor construction** from **regression form**. The descriptor is a rotationally invariant spectral fingerprint of the local atomic environment, while the original regression model is linear in those descriptors. This distinguishes SNAP from Gaussian Approximation Potential (GAP), which uses the same bispectrum representation but models local energy with Gaussian-process regression rather than an explicit linear expansion [1409.3880].

Later work positioned SNAP more broadly as a DFT-trained interatomic potential capable of carrying ab initio fidelity into classical molecular dynamics at scales relevant to defects, phase boundaries, nucleation, and extreme-pressure thermodynamics. In the carbon literature, for example, SNAP is presented as a route from small-cell QMD/AIMD to simulations of millions of atoms, with applications to melting curves, shock Hugoniots, and phase boundaries under multi-TPa conditions [2205.01209].

## 2. Descriptor construction and invariance structure

For an atom \(i\), SNAP begins from a local neighbor density defined within a cutoff sphere. In the original formulation this density is written as
\[
\rho_i ({\bf r}) = \delta({\bf r}) + \sum_{r_{ii'} < R_{cut}} f_c(r_{ii'})\, w_{i'} \,\delta({\bf r}-{\bf r}_{ii'}),
\]
where the central atom has unit weight, \(w_{i'}\) distinguishes species, and \(f_c(r)\) smoothly turns off neighbor contributions at the cutoff [1409.3880].

The 3D neighborhood is then mapped onto a subset of the 3-sphere by promoting the radial coordinate to a fourth angular coordinate,
\[
\theta_0 = \theta_0^{max}\frac{r}{R_{cut}},
\]
and the density is expanded in 4D hyperspherical harmonics \(U^j_{m,m'}(\theta_0,\theta,\phi)\):
\[
\rho({\bf r}) = \sum_{j=0,\frac{1}{2},\ldots}^{\infty}\sum_{m=-j}^{j}\sum_{m'=-j}^{j} u^j_{m,m'} U^j_{m,m'}(\theta_0,\theta,\phi).
\]
The coefficients \(u^j_{m,m'}\) are not themselves invariant under rotation, so SNAP forms rotationally invariant bispectrum components as scalar triple products of the hyperspherical-harmonic coefficients [1409.3880].

In one explicit multicomponent statement of the formalism, the local density may be written as
\[
\rho_i(\vec r)=\delta(\vec r-\vec R_i)+\sum_j w_{Z_j}\,\delta(\vec r-\vec R_j)\,f_c(R_{ij}),
\]
followed by the hyperspherical-harmonic expansion
\[
\rho_i(\vec r)=\sum_{J=0}^{\infty}\sum_{m,m'=-J}^{J} c_{m',m}^{J}\,U^{J}_{m',m}(\theta,\phi,\theta_0),
\]
and the bispectrum invariant
\[
B_i^{J,J_1,J_2}=\sum_{m_1',m_1=-J_1}^{J_1} c_{m_1',m_1}^{J_1} \sum_{m_2',m_2=-J_2}^{J_2} c_{m_2',m_2}^{J_2} \sum_{m',m=-J}^{J} C_{mm_1m_2}^{J,J_1,J_2} C_{m'm_1'm_2'}^{J,J_1,J_2} \left(c_{m',m}^{J}\right)^* .
\]
In practice the expansion is truncated at \(J_{\max}\), and the retained bispectrum components form the descriptor vector \(\vec B_i\) [2303.16597].

The resulting descriptor set is explicitly presented as invariant under translation, rotation, and permutation of identical atoms, while retaining angular and radial information in a compact local representation. The original SNAP paper also emphasizes a bispectrum symmetry,
\[
\frac{B_{j_1,j_2,j}}{2j+1}  = \frac{B_{j,j_2,j_1}}{2 j_1+1}  = \frac{B_{j_1,j,j_2}}{2 j_2+1},
\]
which reduces redundant components and lowers the asymptotic cost of bispectrum-derivative evaluation from \(O(J^4)\) to \(O(J^2)\) [1409.3880].

## 3. Linear SNAP, quadratic SNAP, and the regression problem

The original SNAP energy model is linear in the bispectrum descriptors. In the elemental form,
\[
E^i_{SNAP}(B^i) = \beta_0 + \sum_{k=1}^K \beta_k B_k^i = \beta_0 + \boldsymbol\beta\cdot {\bf B}^i,
\]
so that
\[
E_{SNAP}(r^N) = N \beta_0 + \boldsymbol\beta\cdot\sum_{i=1}^{N} {\bf B}^i.
\]
Because the energy is linear in \(\boldsymbol\beta\), the forces and stress are also linear in the same coefficients:
\[
{\bf F}^j_{SNAP} = - \boldsymbol\beta\cdot \sum_{i=1}^{N} \frac{\partial {\bf B}^i}{\partial {\bf r}_j},
\]
\[
{\bf W}_{SNAP} =  - \boldsymbol\beta\cdot  \sum_{j=1}^{N}  {\bf r}_j \otimes \sum_{i=1}^{N} \frac{\partial {\bf B}^i}{\partial {\bf r}_j}.
\]
Training is therefore posed as a weighted least-squares regression against QM energies, forces, and stresses, assembled as a linear system \({\bf A}\cdot\boldsymbol\beta={\bf y}\) and solved by QR factorization in the foundational work [1409.3880].

A later extension introduced **quadratic SNAP**, in which the atomic energy includes pairwise products of bispectrum components:
\[
E_{\mathrm{SNAP}}^i = \beta_0 + \sum_k \beta_k B_k^i + \frac12 \sum_{k,l}\alpha_{kl} B_k^i B_l^i.
\]
This preserves the same descriptor space but augments the regression basis. The quadratic form is explicitly described as structurally similar to the embedded atom method (EAM), with bispectrum components acting as counterparts to the two-body density functions in EAM. In this interpretation, linear SNAP is a four-body potential, whereas quadratic SNAP can involve up to seven distinct atoms [1711.11131].

The increase in expressivity is accompanied by a large increase in parameter count. For example, the quadratic-SNAP paper reports that at \(J_{\max}=4\), \(K=55\), the linear model has \(55\) coefficients while the quadratic form has \(1540\) quadratic coefficients; at \(J_{\max}=5\), \(K=91\), the quadratic coefficient count rises to \(4186\). The same work reports that the quadratic model produces about an order-of-magnitude lower training energy error than the linear form for tantalum at high \(J_{\max}\), but also requires substantially more training data to prevent overfitting [1711.11131].

Quadratic forms have subsequently been adopted in production potentials. For carbon under extreme conditions, the atomic energy was written as
\[
E^i=\boldsymbol{\beta}\cdot \mathbf{B}^{i}+\frac{1}{2}\mathbf{B}^{i}\cdot \boldsymbol{\alpha}\cdot \mathbf{B}^{i},
\]
and for cobalt the q-SNAP formulation was defined as “SNAP bispectrum descriptors + quadratic regression terms,” with the main distinction from standard SNAP stated explicitly as: **SNAP is linear in descriptors, q-SNAP adds a quadratic feature expansion** [2205.01209; 2404.02626].

## 4. Training databases, fitting workflows, and domain of applicability

SNAP fitting is inseparable from training-set design. The original tantalum work emphasized a training set spanning displaced crystals, strained cells, generalized stacking faults, surfaces, and liquid snapshots rather than only equilibrium structures [1409.3880]. Later studies made this dependence on coverage still more explicit.

For cobalt, the q-SNAP training set contains **1049 cobalt structures** spanning bulk unit cells of hcp, fcc, bcc, and \(\omega\) phases; strained unit cells; bulk MD snapshots; vacancies and interstitials; surfaces with adatoms, vacancies, and steps; nanoparticles; small nanoclusters; and liquid/amorphous structures. The same study states that inclusion of **very small clusters (2–30 atoms)** was essential because excluding them caused failure for highly undercoordinated surface and cluster configurations. Its final weighted least-squares fit used descriptor scans over \(J_{\max}=2\) to 6 and cutoff radii from 4 to 6 Å, with the best cobalt model obtained at **band limit 8** (\(J_{\max}=4\)) and **radial cutoff 5.0 Å**; the final fitting weights were energies \(=1\), forces \(=0.01\), and stresses \(=10^{-6}\) [2404.02626].

Carbon at extreme conditions was trained with a different but equally explicit workflow: a database of **636 structures**, **124,907 unique atomic environments**, and **388,077 regression equations**, followed by an iterative cycle of database construction, parameter training, validation against QMD and experiment, deficiency identification, and retraining. Hyperparameters \(r_{cut}=2.7\) Å and \(J=4\) were selected to balance accuracy and computational cost [2205.01209].

Other studies emphasized data-reduction and diversity selection. The molybdenum work used PCA of bispectrum features and their derivatives to inspect descriptor-space coverage and excluded a vacancy dataset judged redundant in that space [1706.09122]. The niobium study used **Vendi Score** and **D-optimal selection**, trained on **550 total configurations**, and then validated on elastic, thermal, and defect properties including collision cascades [2502.03126].

A recurrent limitation is that SNAP is reliable primarily when prediction remains interpolative in descriptor space. In binary-trained ternary screening, an ensemble of SNAPs trained only on binary phases predicted ternary energies near equilibrium reasonably well but failed for unrelaxed prototype geometries and for structural relaxation itself. For Cu–Ag–Au, the reported RMSE was **17.06 meV/atom** on the ternary test set, remained **36.98 meV/atom** for relaxed ternary prototypes, but rose to about **1.2 eV/atom** for the same prototypes in their initial unrelaxed geometries. Feature-space analysis by PCA showed that the relaxation path traversed regions poorly represented in the training set [2303.16597].

The same ternary study compared three extrapolation metrics and found that the ensemble standard deviation \(\sigma\) correlated best with actual error: in the Cu–Ag–Au test set the reported correlation coefficients were \(R^2=0.84\) for \(\sigma\), \(R^2=0.65\) for the extrapolation grade \(\gamma\), and much worse for a globalized minimum Euclidean distance \(d_\mathrm{min}\). This suggests that domain-of-applicability control is an operational part of SNAP deployment rather than a secondary diagnostic [2303.16597].

## 5. Formal extensions beyond the original single-channel local model

Several extensions retain the spectral-neighbour logic of SNAP while modifying the local density, the descriptor channels, or the total-energy decomposition.

| Extension | Mechanism | Representative result |
|---|---|---|
| **Quadratic SNAP** | Adds \(B_k B_l\) terms to the linear bispectrum model | Order-of-magnitude training-energy improvement for Ta; deployed for C and Co [1711.11131] |
| **EME-SNAP** | Partitions neighbor density into element-resolved partial densities | Descriptor count scales as \(\mathcal{O}(N_{elem}^3)\), force cost as \(\mathcal{O}(N_{elem}^2)\) [2003.11570] |
| **eSNAP** | Adds screened Ewald electrostatics to local SNAP energy | For \(\alpha\)-\(\mathrm{Li_3N}\), \(E_p=\gamma E_{el}+E_{SNAP}\) with \(\gamma=0.057\) [1901.08749] |

The **explicit multi-element extension** (EME-SNAP) replaces the weighted-density treatment of chemistry by partitioning the neighbor density into partial densities for each chemical element:
\[
\rho({\bf r})=\sum_{\mu=1}^{N_{elem}}\rho^\mu({\bf r}),
\]
with each \(\rho^\mu\) expanded separately in hyperspherical harmonics. The resulting explicit multi-element bispectrum is indexed by ordered element triplets,
\[
B^{\kappa\lambda\mu}_{j_1j_2j} = \frac{1}{2j+1}\, u^{\kappa}_{j_1}\otimes_{j_1j_2j} u^{\lambda}_{j_2}\cdot (u^\mu_j)^* .
\]
For \(N_{elem}\) elements, the number of descriptors grows as \(\mathcal{O}(N_{elem}^3)\), while the force calculation in LAMMPS scales as \(\mathcal{O}(N_{elem}^2)\). In InP, this chemically resolved formulation improved defect energetics relative to weighted-density SNAP while retaining good bulk-property accuracy [2003.11570].

The **electrostatic SNAP** formulation addresses the inability of a finite-cutoff local descriptor to represent long-range Coulomb interactions in ionic systems. For \(\alpha\)-\(\mathrm{Li_3N}\), the total potential energy was written as
\[
E_p = \gamma E_{el} + E_{SNAP},
\]
with forces
\[
\mathbf{F}_j = -\gamma \nabla_j E_{el} - \mathbf{F}_{j,\mathrm{SNAP}}.
\]
Using formal charges \(+1\) for Li and \(-3\) for N, the optimized screening factor was \(\gamma=0.057\), interpreted as an effective dielectric screening [1901.08749].

A separate line of development generalized the local scalar density to a **local vector density** for systems with vector fields such as spin:
\[
\mathbf \rho(\mathbf r)=\sum_a w_a\, h_a(\mathbf r-\mathbf r_a)\,\mathbf v_a.
\]
Spatial and vector angular momenta are coupled with Clebsch–Gordan coefficients, and the invariant descriptor is a power spectrum
\[
p_{nlJ}=\sum_M |u_{nlJM}|^2.
\]
The total energy is again decomposed into local contributions and may be fit by linear Ridge regression “as in conventional spectral neighbour analysis potentials,” or by a Gaussian approximation for more nonlinear magnetic energy surfaces [2202.13773].

## 6. Applications, computational realization, and limits

SNAP has been applied to a broad range of materials problems that are explicitly represented in the literature. In tantalum, the foundational model reproduced solid-phase energetics, liquid structure, melting behavior, surface and defect energies, and the screw-dislocation Peierls barrier with the correct single-hump shape seen in DFT [1409.3880]. In tungsten, a SNAP model optimized with a parallel genetic algorithm reproduced elastic constants and defect formation energies closely and could be combined with a helium pair potential to model W–He defect energetics relevant to plasma-facing materials [1702.07042].

For carbon, a quadratic SNAP trained on QMD/DFT data up to **5 TPa** and **20,000 K** reproduced phase-diagram features, melting curves, and the shock Hugoniot within about **3%** of QMD, and was used for **1 million atom polycrystalline diamond** melting simulations [2205.01209]. A later review/commentary presented SNAP as an enabling technology for revisiting the carbon phase diagram, including the diamond-to-BC8 transition, metastability, and nucleation windows under dynamic compression [2506.02328].

For cobalt, q-SNAP closely reproduced phonon dispersions of hcp and fcc Co, surface energies, and nanoparticle stability ordering. The work emphasized its numerical stability in long MD: a **100 ns** simulation of a 2406-atom truncated octahedron at 300 K remained stable throughout, the nanoparticle heat capacity converged to less than \(2\ \mathrm{J\,K^{-1}\,mol^{-1}}\) after 100 ns, and extrapolation of size-dependent melting temperatures yielded a bulk melting point of **1785 K**, only 15 K above experiment [2404.02626].

For niobium, a linear SNAP with ZBL short-range repulsion was trained to resolve a long-standing defect-physics failure of classical Nb potentials. It reproduced the DFT ordering of self-interstitial dumbbells, predicting the \(\langle 111\rangle\) dumbbell as the lowest-energy SIA configuration, and in **5 keV collision cascade simulations** it produced **98%** \(\langle 111\rangle\) single dumbbells in the primary-damage state, whereas FS and EAM models overwhelmingly yielded \(\langle 110\rangle\) dumbbells [2502.03126].

SNAP has also been integrated into surrogate workflows rather than used only as a final force field. In ternary convex-hull construction, ensembles of binary-trained SNAPs were used to pre-relax and screen very large pools of ternary derivative structures before DFT refinement. In the Mo–Ta–W case, the workflow evaluated **331,734 ternary structures** using SNAP energies and forwarded only the 15 best candidates per stoichiometry to DFT [2308.15907].

The formalism carries a clear computational cost. The original SNAP paper states that it is more expensive than simpler empirical potentials by roughly one to two orders of magnitude in that implementation because bispectrum and bispectrum-derivative evaluation are costly [1409.3880]. Performance engineering has therefore been a major parallel development. In the Kokkos CUDA backend of LAMMPS, algorithmic refactorization, adjoint reformulation, hierarchical parallelism, and memory-layout redesign produced a roughly **22× time-to-solution improvement** on an NVIDIA Tesla V100-16GB relative to an existing implementation [2011.12875].

Several limitations are recurring rather than incidental. SNAP is fundamentally local and finite-cutoff, so short-range collision physics is often delegated to a ZBL reference term and long-range electrostatics require explicit augmentation in ionic systems [1702.07042; 1901.08749]. The reliability of a given fit depends on whether relevant local environments are represented in the training set; the ternary-screening studies show that poor feature-space coverage can make energies and forces unreliable along relaxation paths outside the sampled domain [2303.16597]. Quadratic SNAP can improve accuracy substantially, but it requires substantially more data and exhibits a stronger overfitting risk [1711.11131]. In extreme-condition carbon simulations, additional limitations were identified from the classical-potential side: no explicit electronic entropy and a hydrodynamic approximation in Hugoniot calculations [2205.01209].

Taken together, these results define SNAP not as a single immutable model but as a family of bispectrum-based local interatomic potentials whose central architectural choice is stable across applications: encode the neighborhood spectrally, enforce geometric invariances at the descriptor level, and fit energies, forces, and stresses through a regression model that remains explicit enough to support large-scale molecular dynamics.

Source: https://www.emergentmind.com/topics/spectral-neighbor-analysis-potential-snap