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

# Spectral Neighbor Analysis Potential (SNAP/qSNAP)

The Spectral Neighbor Analysis Potential (SNAP) and its quadratic extension (qSNAP) are machine-learned interatomic potential frameworks that rigorously encode local atomic environments using high-order, rotationally and permutationally invariant bispectrum descriptors. Designed to interpolate between ab initio quantum accuracy and the computational efficiency of empirical force fields, SNAP/qSNAP have demonstrated quantum-level fidelity in complex materials modeling tasks ranging from metallic alloys and high-entropy ceramics to extreme-pressure phases and nanoparticle thermodynamics [2404.02626, 1711.11131, 2205.01209].

## 1. Mathematical Foundations and Descriptor Construction

SNAP encodes the atomic environment of each atom \(i\) as a local neighbor density,
\[
\rho_i(\mathbf r)=\sum_{j} f_c(r_{ij})\,\delta(\mathbf r-\mathbf r_{ij})
\]
where \(f_c\) is a smooth cutoff function vanishing at \(r_{ij} > r_{cut}\), and the sum includes all neighbors within the cutoff radius. To ensure rotation and permutation invariance and to capture many-body correlations, \(\rho_i\) is projected onto a basis of 4D hyperspherical (SO(4)) harmonics \(U^j_{mm'}(\theta, \phi, \chi)\):
\[
\rho_i(\mathbf r) = \sum_{n} u^{(i)}_{n}\,R_n(r)\,U^j_{mm'}(\Omega)
\]
where \(n \equiv (j, m, m')\), \(R_n(r)\) is an orthonormal radial basis, and \(J_{max}\) sets the maximum angular resolution (typically \(2J_{max}\) defines the band limit).

From the expansion coefficients \(u^{(i)}_n\), rotationally invariant bispectrum components are constructed via Clebsch–Gordan–type couplings:
\[
B_{i,klm} = \sum_{n_1, n_2} C^{klm}_{n_1 n_2}\,u^{(i)}_{n_1}\,u^{(i)}_{n_2}
\]
The full bispectrum feature vector \(B_i\) collects all distinct components up to the truncation.

For multi-element systems, explicit multi-element SNAP (EME-SNAP) generalizes this expansion by defining partial neighbor densities for each chemical species, leading to bispectrum terms \(B_{i; j_1j_2j}^{\kappa\lambda\mu}\) which scale as \(\mathcal{O}(N_\text{elem}^3 J_{\max}^3)\) in descriptor count [2003.11570].

## 2. Energy Model: Linear SNAP and Quadratic Extension (qSNAP)

### Linear SNAP

The atomic energy within the SNAP framework is given as a linear function of the bispectrum descriptors:
\[
E_i^{\rm lin} = \beta_0 + \sum_k \beta_k B_{i,\,k}
\]
where \(\beta_0\) and \(\{\beta_k\}\) are fitting coefficients determined via a weighted least-squares regression over a DFT reference dataset; the total energy sums all atomic contributions. Analytical derivatives with respect to atomic positions provide forces and stress tensors required for MD.

### Quadratic SNAP (qSNAP)

The quadratic SNAP (qSNAP) augments the expressiveness of SNAP by including quadratic terms in the bispectrum components:
\[
E_i^{\rm q} = \beta_0 + \sum_k \beta_k B_{i,\,k} + \sum_{k \le \ell} \gamma_{k\ell} B_{i,\,k} B_{i,\,\ell}
\]
The quadratic coefficients \(\{\gamma_{k\ell}\}\) systematically capture many-body effects, allowing improved treatment of strongly anharmonic, defect-rich, or undercoordinated environments. The increase in parameter count (\(\mathcal{O}(K^2)\)) necessitates larger, more diverse training datasets and robust regularization against overfitting [1711.11131].

## 3. Training and Optimization Workflow

SNAP and qSNAP potentials are trained against a database of DFT-calculated energies, forces, and (optionally) stresses for structurally and chemically diverse configurations. The fit minimizes a weighted objective function:
\[
L(\beta, \gamma) = \sum_{\text{structures}}\left[w_E (E^\text{SNAP} - E^\text{DFT})^2 + w_F \sum_{\alpha} (F_{\alpha}^\text{SNAP} - F_{\alpha}^\text{DFT})^2 + w_S \sum_{\alpha \leq \beta} (\sigma_{\alpha\beta}^\text{SNAP} - \sigma_{\alpha\beta}^\text{DFT})^2\right]
\]
Group-wise weights \(w_E, w_F, w_S\) are tuned (often via an outer-loop genetic algorithm or cross-validation) to reflect the intended application targets (e.g., elastic properties, defect energetics, phase stability). The hyperparameters \(J_{max}\) and \(r_{cut}\) are chosen to minimize validation RMSE.

For qSNAP, the data volume must be sufficient to avoid overfitting; random holdout cross-validation, principal component analysis, and genetic algorithm optimization of group weights are employed to ensure generalization [1711.11131, 2404.02626].

## 4. Validation: Accuracy and Transferability

### Bulk, Defect, and Surface Properties

SNAP/qSNAP potentials have demonstrated near-DFT accuracy across a wide set of observables:
- Energy RMSE and Force RMSE: For example, for Co q-SNAP yields 4.3 meV/atom (train), 28.6 meV/atom (test) and forces within 54.2–92.7 meV/Å [2404.02626].
- Phonon dispersions: Deviation from DFT <0.5 THz in hcp and fcc Co [2404.02626].
- Surface energies: Maximum error 47 mJ/m\(^2\) in Co, compared to >270 mJ/m\(^2\) for empirical EAM [2404.02626].
- Elastic moduli, formation energies, melting points: Systematically within 5–10% (often lower) of DFT.

### Nanoparticle Thermodynamics and Size-Scaling

The flexibility and stability of qSNAP permit long-time MD simulations of large systems (e.g., 9201-atom Co nanoparticles, 100 ns trajectories), enabling calculation of converged heat capacity, melting temperatures, and vacancy formation energies, reproducing both DFT and experimental scaling laws (e.g., Gibbs–Thomson scaling of nanoparticle melting \(T_{m,NP} = T_{m,bulk} - aN^{-1/3}\)) [2404.02626].

### Extreme Conditions and Extrapolation

qSNAP provides robust performance outside the training regime (e.g., up to 5 TPa and 20,000 K in carbon), with phase boundaries and Hugoniot curves accurate to within 3% versus quantum molecular dynamics and experiment [2205.01209]. However, as the number of species or descriptors increases, numerical stability and the risk of overfitting become significant, particularly for qSNAP applied to large, chemically diverse datasets (>3000 structures) [2304.01650].

## 5. Computational Efficiency and Implementation

The computational cost of SNAP scales as \(\mathcal{O}(KN_\text{atoms})\) per MD step due to the bispectrum calculation and analytic force derivation. The transition from linear SNAP to qSNAP (for fixed \(K\)) increases cost only by a factor of ≲2, as the quadratic form is efficiently computed without new neighbor loops [1711.11131]. For practical MD (e.g., 10\(^4\)–10\(^6\) atoms, 100 ns), qSNAP achieves throughput comparable to optimized EAM within LAMMPS, benefitting from vectorization, mixed-precision kernels, and efficient MPI/OpenMP scaling [2404.02626, 2205.01209]. Memory footprint and runtime scale with the descriptor count; for \(J_{max}=4\), typical \(K\) is 50–80.

## 6. Applications and Extensions

SNAP/qSNAP have enabled quantum-accurate, long-time, large-scale simulations for:
- Elemental and multi-component alloys (Ni–Mo, W–Be, W–ZrC, Li–TMO\(_2\)) [1806.04777, 2003.11570, 2212.01432, 2304.01650].
- Nanoparticle melting, heat capacity, structural transitions, size-dependent stability [2404.02626].
- Phase diagrams at extreme conditions (carbon under GPa–TPa pressures) [2205.01209].
- Radiation damage and defect evolution in metals (Nb, W) [2502.03126, 1702.07042].
- High-throughput materials screening (ternary convex hulls, multicomponent alloys) [2303.16597, 2308.15907].

Systematic extensions include explicit multi-element descriptors [2003.11570], spin- and vector-field coupling for magnetic materials [2202.13773], and neural network or Gaussian-process regression atop the SNAP bispectrum for further flexibility (but at increased training and computational cost) [2005.04332].

## 7. Limitations and Practical Considerations

- **qSNAP parameter inflation:** Quadratic coupling increases the number of parameters by \(\mathcal{O}(K^2)\), making the approach data-hungry and susceptible to overfitting unless large, appropriately diverse training sets are provided [1711.11131, 2304.01650].
- **Numerical stability:** qSNAP can suffer from ill-conditioning and convergence failures with very large descriptor sets or insufficient regularization; cross-validation and data scaling are essential safeguards [2304.01650].
- **Computational cost:** While much faster than quantum methods, SNAP is typically 10–100× slower than simple EAM/MEAM, particularly at high expansion order [1706.09122, 1711.11131].
- **Transferability:** Linear SNAP is often preferred in high species-count, disordered, or high-entropy systems due to its stability and modest data requirements [2304.01650]. qSNAP is recommended where sub-meV/atom accuracy is critical and sufficient training data are available.

---

**References**:
- [2404.02626] Polyvalent Machine-Learned Potential for Cobalt: from Bulk to Nanoparticles
- [1711.11131] Extending the Accuracy of the SNAP Interatomic Potential Form
- [2304.01650] Constructing and evaluating machine-learned interatomic potentials for Li-based disordered rocksalts
- [2205.01209] Machine Learning Interatomic Potential for Simulations of Carbon at Extreme Conditions
- [2003.11570] Explicit Multi-element Extension of the Spectral Neighbor Analysis Potential for Chemically Complex Systems
- [2502.03126] A Robust Machine Learned Interatomic Potential for Nb: Collision Cascade Simulations with accurate Defect Configurations
- [2303.16597] Machine-Learning Surrogate Model for Accelerating the Search of Stable Ternary Alloys
- [2212.01432] Machine Learned Interatomic Potential for Dispersion Strengthened Plasma Facing Components
- [1806.04777] Quantum-Accurate Spectral Neighbor Analysis Potential Models for Ni-Mo Binary Alloys and FCC Metals

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