---
title: Tabulated Gaussian Approximation Potentials (tabGAP)
url: https://www.emergentmind.com/topics/tabulated-gaussian-approximation-potentials-tabgap
type: topic
---

# Tabulated Gaussian Approximation Potentials (tabGAP)

Tabulated Gaussian Approximation Potentials (tabGAP) constitute a class of machine-learned interatomic potentials designed to reconcile ab initio accuracy with the computational efficiency required for large-scale atomistic simulations. tabGAP is derived by tabulating the contributions of a low-dimensional Gaussian Approximation Potential (GAP)—typically restricted to two-body, three-body, and sometimes embedded-atom (EAM) density terms—onto regular grids, enabling energy and force evaluations through fast spline interpolation rather than runtime kernel summation. This methodology yields ≥100× speed-ups over full GAP implementations with only minor loss of accuracy on structural, defect, and dynamical properties, as demonstrated across metals, oxides, alloys, and complex functional materials [2205.09165, 2203.08458, 2212.03096, 2408.15779, 2511.22592, 2512.24430, 2511.21464, 2208.00804].

## 1. Mathematical Formulation and Descriptor Decomposition

The original GAP formalism expresses the total potential energy as a sum over local atomic energies, each modeled via Gaussian process regression (GPR) on a set of low-dimensional environment descriptors. The key components are:

- **Regression Framework**:
  $$
  \hat{E}(x_*) = k(x_*, X)\left[K + \sigma_n^2 I\right]^{-1} y
  $$
  where $x_*$ is a new atomic environment (described by a vector of descriptors), $X$ is the set of $M$ sparse reference environments, $y$ the vector of DFT energies/forces, $K$ the kernel matrix, and $\sigma_n^2$ the assumed noise variance. The most used kernel is the squared-exponential (SE) kernel:
  $$
  k_{\text{SE}}(x, x') = \sigma_f^2 \exp\left(-\frac{\|x - x'\|^2}{2\ell^2}\right)
  $$
  with $\sigma_f$ controlling variance, $\ell$ the length scale.

- **Descriptor Decomposition**:
  - **Two-body**: Pair distances $r_{ij}$.
  - **Three-body**: Triplet configurations $(r_{ij}, r_{ik}, \cos\theta_{jik})$.
  - **Embedding-like**: Local atomic density $\rho_i = \sum_j \phi(r_{ij})$, typically with $\phi$ a smooth function.

Each descriptor channel has its own kernel and set of regression weights.

## 2. Tabulation Strategy and Spline Construction

To avoid the $O(M)$ runtime cost per atomic environment in standard GAP, tabGAP precomputes the regression predictions on regular grids in descriptor space:

- **1D Tables (Two-body, EAM)**: For each relevant element-pair, energies (and sometimes density functions) are tabulated over a grid of $r \in [0, r_{cut}]$ (e.g., 500–5,000 points).
- **3D Tables (Three-body)**: For each element triplet, energies are stored over a uniformly spaced grid in $(r_1, r_2, \mu)$ (typically $40^3$–$120^3$ points).
- **Interpolation**: Cubic-spline interpolation is employed for both 1D and 3D tables, guaranteeing $C^2$ continuity (continuous energy, force, and second derivatives) [2205.09165].

At runtime, local descriptors are mapped to table indices and interpolated; analytic spline derivatives yield forces.

## 3. Training Databases and Point Selection

tabGAP inherits its expressiveness from the quality and diversity of its GAP training database, which must include:

- **Bulk phases over varied volumes and configurations**: e.g., ±5–10% strain, multiple polymorphs for oxides or alloys.
- **Surfaces, defective states, and high-energy snapshots**: e.g., DFT of surfaces [(001), (110), ...], point and extended defects (vacancies, SIAs, interstitials), and snapshots from high-temperature or radiation events [2205.09165, 2408.15779, 2212.03096, 2512.24430].
- **Liquids and non-equilibrium configurations**: e.g., melt-quench or finite-T MD structures.

Representative environments (sparse points) for GPR are selected using strategies such as farthest-point sampling, CUR decomposition, or k-means in descriptor space. For tabulation, uniform or adaptively sampled grids in each descriptor domain are constructed, with density chosen to achieve prescribed force/energy interpolation errors (<1 meV/Å, <0.01 eV/atom) [2511.22592].

## 4. MD Workflow and Efficient Implementation

Each MD step with tabGAP involves:

1. **Neighbor List Construction** within species-appropriate cutoffs for each descriptor channel.
2. **Spline Lookups** for all pairwise and three-body (and EAM) contributions:
   - For every pair $(i, j)$, evaluate 1D spline at $r_{ij}$.
   - For each triplet $(i, j, k)$, evaluate 3D spline at $(r_{ij}, r_{ik}, \cos\theta_{jik})$.
   - For each atom $i$, accumulate density and call the embedding spline.

3. **Force Calculation**: Analytic derivatives of the spline representation provide forces, ensuring smooth dynamics and energy conservation.

4. **Short-Range Physics**: Explicit repulsive terms—frequently ZBL-type, with pair-specific parameters fitted to all-electron DFT dimers—are tabulated and smoothly merged into the short-range domain.

The per-atom, per-timestep computational scaling matches that of classical empirical potentials for the same cutoff range, with achievable 80–400× speedups over full GAP [2205.09165, 2512.24430, 2212.03096].

## 5. Accuracy, Validation, and Performance Benchmarks

tabGAP delivers near-DFT fidelity across a spectrum of materials systems and benchmarks:

| System (Reference)                | Energy RMSE (meV/atom) | Force RMSE (meV/Å) | Speedup over GAP |
|-----------------------------------|------------------------|--------------------|-----------------|
| W (2205.09165)                    | ~2                     | ~40                | ~5×             |
| Cu, Al, Ni (2408.15779)           | 0.5–1.6                | 15–70              | 80–150×         |
| YBa$_2$Cu$_3$O$_{7-\delta}$ (2511.22592) | 3.3                    | 90                 | ~12×            |
| $\alpha$-Fe-H (2511.21464)        | 1.88                   | 69                 | 100–200×        |
| Ga$_2$O$_3$ (2212.03096)          | ~11                    | 160                | 400×            |
| Mo-Nb-Ta-V-W (2203.08458)         | ~2 (crystal)           | 60 (crystal)       | 40×             |

Key properties—lattice constants, elastic constants, stacking-fault energies, defect energetics, and threshold displacement energies—are reproduced within typical deviations of <1% for lattice/elastic properties, <0.02 eV for migration/barrier energies, and <10% for surface and stacking-fault energies relative to DFT or experiment [2205.09165, 2408.15779, 2512.24430, 2511.22592, 2511.21464].

## 6. Applications and Transferability

tabGAP is widely applied for:

- **Large-scale MD simulations** (≥10$^6$ atoms) of nanoindentation, defect dynamics, irradiation cascades, and mechanical loading [2205.09165, 2512.24430, 2408.15779].
- **Materials chemistry**: High-entropy alloys (HEAs), transition metals, complex oxides (e.g., exploration of Ga$_2$O$_3$ polymorphs [2212.03096]).
- **Radiation damage and high-energy processes**: Accurate reproduction of defect creation, recombination, and clustering statistics in both metals (W, Mo, Fe) and superconducting oxides (YBa$_2$Cu$_3$O$_{7-\delta}$) [2512.24430, 2511.22592, 2208.00804].

The combination of low-dimensional descriptors, tabulation, and spline interpolation yields a practical approach for extending GAP accuracy to multi-billion-atom, long-timescale MD required for materials design, screening, and process modeling.

## 7. Limitations, Implementation Considerations, and Future Directions

tabGAP performance is fundamentally constrained by the descriptor set:

- **Descriptor Dimensionality**: Only descriptors of ≤3 dimensions (e.g., 2b, 3b, EAM) are amenable to practical spline tabulation. High-dimensional descriptors (e.g., SOAP, message-passing) cannot be efficiently tabulated, limiting tabGAP’s ability to capture subtle symmetries or long-range charge/correlation effects [2212.03096, 2511.22592].
- **Grid Density and Memory Usage**: Detailed grids (e.g., $80^3$ for three-body) require significant memory, and undersampling degrades force accuracy.
- **Chemistry and Transferability**: tabGAPs are only as transferable and robust as the diversity of their DFT training sets; substantial retraining is needed to accommodate new chemistries or environments (e.g., surfaces, interfaces, new phases) [2408.15779, 2511.22592].

Best practice recommends explicit inclusion of all relevant configurations and defects in the DFT dataset, careful selection of cutoff radii and grid sizes, and validation against higher-level calculations or experimental observables. For high-accuracy or strongly correlated systems, augmentation with higher-body or long-range descriptors remains necessary (often with a trade-off in performance) [2212.03096]. 

tabGAP has proven a robust, reproducible, and computationally efficient approach for quantum-accurate large-scale atomistic simulation in diverse materials, broadly adopted and validated in recent machine-learning interatomic potential deployments.

Source: https://www.emergentmind.com/topics/tabulated-gaussian-approximation-potentials-tabgap