The Smooth Overlap of Atomic Positions (SOAP) descriptor is a mathematical construct that transforms the discrete, point-like positions of atoms surrounding a central atom into a continuous, high-dimensional vector. It achieves this by representing each neighboring atom as a Gaussian-smeared density, which is then expanded using a complete basis set of orthogonal radial functions and spherical harmonics. The resulting power spectrum, computed from the expansion coefficients, ensures the descriptor is strictly invariant to atomic permutations, translations, and rotations of the entire environment.
Glossary
SOAP Descriptor

What is SOAP Descriptor?
The SOAP descriptor is a foundational method for encoding local atomic environments into a smooth, invariant, and complete vector representation, widely used in machine learning for materials science and molecular informatics.
This invariance makes SOAP a critical input feature for Graph Neural Networks (GNNs) and kernel-based models like Gaussian Approximation Potentials (GAP). By capturing many-body correlations up to a specified cutoff radius, it provides a universal and systematically improvable fingerprint of local geometry. Its differentiability and smoothness allow for the direct computation of atomic forces, making it a cornerstone representation for constructing highly accurate Neural Network Potentials (NNPs) that predict energy and properties directly from atomic coordinates.
Key Features of SOAP Descriptors
The SOAP descriptor provides a smooth, rotationally invariant representation of local atomic environments, forming a foundational input for machine learning interatomic potentials and materials property prediction.
Gaussian Smearing of Atomic Density
SOAP constructs a continuous atomic density field by centering a Gaussian function on each neighboring atom within a spherical cutoff radius. This smearing replaces discrete point particles with a smooth density distribution, ensuring that the descriptor varies continuously as atoms move. The width of the Gaussian is a key hyperparameter controlling the smoothness and sensitivity of the representation.
Radial and Spherical Harmonic Expansion
The smeared density is expanded in a basis set of orthogonal radial functions and spherical harmonics. This expansion decomposes the 3D density field into a set of expansion coefficients that capture angular and radial information separately. The number of radial and angular basis functions controls the resolution and completeness of the representation.
Power Spectrum for Rotation Invariance
To achieve strict rotational invariance, SOAP computes the power spectrum of the expansion coefficients. This involves forming a tensor product of the coefficients and contracting them over the angular momentum indices using Clebsch-Gordan coefficients. The resulting invariant features are insensitive to the absolute orientation of the local environment, making them ideal inputs for rotationally invariant models.
Kernel-Based Similarity Measure
The SOAP descriptor naturally defines a positive-definite kernel that quantifies the similarity between two atomic environments. This kernel is computed as the dot product of the power spectrum vectors, providing a rigorous metric for comparing local structures. It is widely used in Gaussian process regression models for constructing interatomic potentials.
Multi-Scale Resolution Control
Two key hyperparameters govern the descriptor's resolution:
- n_max: Number of radial basis functions, controlling radial detail
- l_max: Maximum angular momentum quantum number, controlling angular detail Increasing these values yields a more complete description at higher computational cost, allowing systematic convergence of accuracy.
Frequently Asked Questions
Explore the foundational concepts behind the Smooth Overlap of Atomic Positions (SOAP) descriptor, a powerful tool for encoding local chemical environments in machine learning models.
A SOAP (Smooth Overlap of Atomic Positions) descriptor is a smooth, invariant representation of a local atomic environment. It works by constructing a Gaussian-smeared atomic density field around a central atom, then expanding this density in a basis set of radial basis functions and spherical harmonics. The key step is computing the rotationally invariant power spectrum of these expansion coefficients, which yields a fixed-length vector that uniquely encodes the geometry and chemical species of the neighboring atoms. This process ensures that the representation is invariant to translations, rotations, and permutations of identical atoms, making it ideal for input into machine learning models for property prediction.
SOAP vs. Other Atomic Descriptors
A comparison of the Smooth Overlap of Atomic Positions (SOAP) descriptor against other common representations for encoding local atomic environments in machine learning interatomic potentials.
| Feature | SOAP | ACSF / Behler-Parrinello | ACE |
|---|---|---|---|
Mathematical Basis | Gaussian smeared atomic density expanded in radial and spherical harmonic basis functions | Hand-crafted symmetry functions (G2, G4) with explicit cutoff functions | Systematic body-ordered expansion of atomic basis functions using cluster expansion formalism |
Rotational Invariance | |||
Permutational Invariance | |||
Smoothness (Continuity of Derivatives) | |||
Completeness (Systematic Convergence) | |||
Many-Body Order Captured | 3-body (via power spectrum), 4-body (via bispectrum) | 2-body (radial) and 3-body (angular) | Arbitrary body-order (systematically expandable) |
Computational Cost | Moderate to High | Low to Moderate | Low |
Typical Dimensionality | 100-1000+ components | 10-100 components | 10-1000+ components |
Enabling Efficiency, Speed & Accuracy
Intelligent Analysis, Decision & Execution
We build AI systems for teams that need search across company data, workflow automation across tools, or AI features inside products and internal software.
Talk to Us
Search across company data
Give teams answers from docs, tickets, runbooks, and product data with sources and permissions.
Useful when people spend too long searching or get different answers from different systems.

Automate internal workflows
Use AI to route work, draft outputs, trigger actions, and keep approvals and logs in place.
Useful when repetitive work moves across multiple tools and teams.

Add AI to products and internal tools
Build assistants, guided actions, or decision support into the software your team or customers already use.
Useful when AI needs to be part of the product, not a separate tool.
Related Terms
Key concepts that contextualize the SOAP descriptor within the broader landscape of geometric deep learning and atomic environment representations.
Atomic Cluster Expansion (ACE)
A systematic and complete basis set expansion of atomic environments that yields highly efficient, body-ordered invariant features. ACE provides a mathematically rigorous framework for constructing interatomic potentials, offering a direct theoretical connection to the SOAP descriptor. While SOAP uses a fixed Gaussian smearing and radial basis, ACE generalizes this to a complete, hierarchical expansion that can be systematically converged. This makes ACE both a powerful standalone method and a theoretical lens through which to understand the completeness of SOAP-like representations.
SE(3) Equivariance
A fundamental symmetry property ensuring that a model's output transforms consistently with any input rotation and translation in 3D Euclidean space. SOAP descriptors are explicitly invariant to these transformations, making them ideal inputs for conventional ML models. Modern architectures like NequIP and MACE build equivariance directly into the neural network layers, learning more complex functions of atomic coordinates while preserving physical symmetries. Understanding this distinction between invariant features and equivariant operations is critical for designing models that respect molecular geometry.
Neural Network Potential (NNP)
A machine-learned surrogate model that predicts the potential energy and atomic forces of a molecular system directly from atomic coordinates. SOAP descriptors serve as a foundational input representation for many classical NNPs, including the pioneering GAP (Gaussian Approximation Potential) framework. By converting 3D coordinates into a smooth, rotationally invariant feature vector, SOAP enables kernel-based or neural network regression to learn the Born-Oppenheimer potential energy surface with high fidelity, bypassing expensive quantum mechanical calculations.
Molecular Fingerprint
A fixed-length bit vector encoding the presence or absence of specific substructural features within a molecule. Unlike SOAP, which captures continuous 3D geometry, traditional fingerprints like ECFP or MACCS keys operate on the 2D molecular graph. SOAP can be seen as a 3D generalization, providing a smooth, differentiable, and geometrically aware alternative. While 2D fingerprints are faster to compute and useful for property prediction, SOAP's 3D sensitivity makes it essential for tasks where conformation and spatial arrangement dictate function.
SchNet
A pioneering continuous-filter convolutional neural network that models quantum interactions using interatomic distances. SchNet learns its own atomic environment representations directly from data, effectively discovering features analogous to SOAP but optimized end-to-end for a specific task. Comparing SOAP (a fixed, physics-based descriptor) with SchNet's learned filters highlights a key design choice in computational chemistry: the trade-off between the interpretability and transferability of hand-crafted features versus the flexibility of learned representations.
Crystal Graph Convolutional Neural Network (CGCNN)
A GNN architecture that directly learns material properties from crystal structures by constructing a multigraph representing atomic connectivity across periodic boundary conditions. While SOAP provides a local descriptor of atomic environments, CGCNN operates on the entire periodic graph. SOAP features can be used as input node features for CGCNN, combining the physics-informed representation of local geometry with the relational learning capabilities of graph neural networks for predicting bulk material properties.

About the author
Prasad Kumkar
CEO & MD, Inference Systems
Prasad Kumkar is the CEO & MD of Inference Systems and writes about AI systems architecture, LLM infrastructure, model serving, evaluation, and production deployment. Over 5+ years, he has worked across computer vision models, L5 autonomous vehicle systems, and LLM research, with a focus on taking complex AI ideas into real-world engineering systems.
His work and writing cover AI systems, large language models, AI agents, multimodal systems, autonomous systems, inference optimization, RAG, evaluation, and production AI engineering.
Partnered with leading AI, data, and software stack.
How We Work
Custom AI workflows for your Business
One-fit-all AI don't work for modern businesses. At Inferensys, we aim to understand your business & custom requirements; which we use to define most efficient agentic workflows, the data, and the tools for your business.
01
Review the use case
We understand the task, the users, and where AI can actually help.
Read more02
Pick the right approach
We define what needs search, automation, or product integration.
Read more03
Build the first useful version
We implement the part that proves the value first.
Read more04
Improve from there
We add the checks and visibility needed to keep it useful.
Read moreThe first call is a practical review of your use case and the right next step.
Talk to Us