---
title: Atomic Cluster Expansion (ACE)
url: https://www.emergentmind.com/topics/atomic-cluster-expansion-ace
type: topic
---

# Atomic Cluster Expansion (ACE)

The Atomic Cluster Expansion (ACE) is a systematically improvable, universally complete, and symmetry-adapted framework for representing the potential energy surface of many-atom systems. As a foundational tool in machine-learned interatomic potentials, ACE underpins recent advances in modeling accuracy and transferability for condensed matter, molecular systems, and even quantum wavefunctions. The core philosophy of ACE is to decompose the total energy into local atomic contributions, each parameterized as a linear (or sparsely nonlinear) combination of many-body rotationally, translationally, and permutationally invariant basis functions. This approach enables the systematic capture of arbitrary multi-atom correlations, efficient linear fitting to quantum-mechanical reference data, and robust extrapolation beyond the training set.

## 1. Mathematical Structure and Symmetry Adaptation

ACE expresses the total potential energy as
\[
E_{\text{tot}}\;=\;\sum_{i=1}^N E_i,
\]
where each atomic energy $E_i$ is expanded in a linear basis of many-body “cluster” functions that capture up to $\nu_{\max}$-body correlations among an atom’s neighbors. These local expansions take the generic form
\[
E_i = \sum_{\nu=1}^{\nu_{\max}}\sum_{n_1\ldots n_\nu,\,l_1\ldots l_\nu} c_{n_1\ldots n_\nu; l_1\ldots l_\nu} \, B_{n_1\ldots n_\nu; l_1\ldots l_\nu}(i),
\]
where $B_{n_1\ldots n_\nu; l_1\ldots l_\nu}(i)$ are scalar, rotationally invariant basis functions constructed as products of radial and angular basis projections over all $\nu$-tuples of neighbors of atom $i$. Radial basis functions (e.g., orthonormal polynomials, Chebyshev, or B-splines) resolve $r_{ij}$ up to a cutoff $r_{\rm cut}$, and angular dependence is encoded via (real) spherical harmonics $Y_{lm}$. These are coupled via Clebsch–Gordan algebra to obtain overall $l=0$ scalars, ensuring invariance under global rotations, translations, and permutation of like atoms [2308.06462], [2303.00753].

The completeness and systematic improvability of ACE follow from its polynomial structure over neighbor densities; the expansion converges to any smooth, symmetry-adapted function given sufficiently large basis cutoffs in body order, radial, and angular indices.

## 2. Parameterization, Regularization, and Training

Fitting an ACE potential consists of generating a diverse set of atomic configurations via ab initio electronic structure methods (e.g. DFT), evaluating energies, forces, and (optionally) virials, and then solving a regularized least squares problem to determine the expansion coefficients. The loss function typically combines errors in total energy and atomic forces,
\[
\mathcal{L} = \sum_k \Big[ w_E \big(E_k^{\text{ACE}} - E_k^{\text{DFT}}\big)^2 + w_F \sum_i \|\mathbf{F}_{ik}^{\text{ACE}} - \mathbf{F}_{ik}^{\text{DFT}}\|^2 \Big] + \lambda \|\mathbf{c}\|^2,
\]
where $w_E$, $w_F$ are relative weights, and $\lambda$ is a Tikhonov regularizer to control overfitting given the typically large number of basis functions [2303.00753], [2406.14306], [2309.11026].

Active learning using extrapolation grades (D-optimality based "MaxVol" selection) enables iterative refinement of the training set by automatically identifying environments where the model is extrapolating and thus requires higher-fidelity quantum data [2212.08716], [2406.14306], [2506.15061]. This efficiently extends coverage to rare or complex configurations, such as high-energy defects or out-of-equilibrium structures.

## 3. Practical Implementation in Materials Modeling

ACE has been deployed across an exceptionally broad variety of materials systems, including elemental metals (Mg, W, Cu, AlN) [2305.03577], [2407.00946], [2311.11990], transition metal alloys (Fe–Co, Pt–Rh, Cu–W) [2303.00753], [2303.07465], [2407.00946], oxides (BaTiO$_3$) [2505.17991], carbon allotropes and layered materials (diamond, graphene, amorphous carbon, multilayer graphene) [2210.09161], [2506.15061], water and aqueous systems [2406.14306], and lattices with multi-degree-of-freedom order parameters (e.g., magnetic ACE for Fe) [2305.15137].

The following table summarizes the key basis and hyperparameter choices in recent ACE applications:

| System                | Body order ($\nu_{\max}$) | Radial cutoff [Å] | Total basis size | Reference |
|-----------------------|---------------------------|-------------------|------------------|-----------|
| Mg (hcp/all phases)   | 4                         | 8.2               | 724              | [2305.03577]   |
| Cu–W (alloys)         | 4                         | 7.0               | $\sim$1000--10000| [2407.00946]   |
| Water (ice, liquid)   | 5                         | 6.0               | 1,774            | [2406.14306]   |
| Fe–Co (alloy)         | 3--4                      | $\sim$6           | $\sim$1,000      | [2303.00753]   |
| Carbon (all phases)   | 5                         | 5.0               | 488              | [2210.09161]   |
| BaTiO$_3$ (ferroelectric) | 4                     | 6.0               | $\sim$10$^4$     | [2505.17991]   |

Energy and force RMSEs of ACE models, after sufficiently comprehensive parameterization, routinely fall below 1--5 meV/atom and 50 meV/Å, enabling predictive-quality molecular dynamics, free-energy calculations, and the computation of sensitive structural, vibrational, and thermomechanical properties [2308.06462], [2407.00946], [2406.14306].

## 4. Variants, Extensions, and Generalizations

ACE admits a hierarchy of extensions that broaden its applicability:

- **Orthonormal/non-self-interacting ACE**: By explicitly removing self-interaction terms (where the same neighbor is counted multiple times within a cluster basis function), improved numerical conditioning, coefficient decay, and regularization are achieved [2401.01550].

- **Non-collinear Magnetic ACE**: Simultaneous expansion in both atomic positions and vectorial/tensorial degrees of freedom (e.g., magnetic moments) enables accurate modeling of coupled spin-lattice dynamics and complex magnetic phase behavior [2305.15137].

- **Charge-constrained ACE**: Augmenting the ACE basis with per-site scalar charges (and optionally higher multipoles), along with global Coulomb or QEq-type terms, captures long-range polarization and charge-transfer phenomena, allowing rigorous variational coupling between short-range ACE and global electrostatics [2411.04062], [2303.12962].

- **Generalized properties**: ACE expansions can fit not only scalar energies but also vector and tensor properties (dipoles, Born charges, stresses), operator matrix elements (tight-binding Hamiltonians), and can parameterize fermionic wavefunctions for quantum Monte Carlo with anti-symmetry enforced via determinantal factors [2206.11375], [2304.04260], [2308.06462].

- **Finnis–Sinclair-type nonlinear embeddings**: Embedding densities and mildly nonlinear transformations (e.g., $E_i = A_i^{(1)} + \sqrt{A_i^{(2)}}$) improve performance in metallic systems but preserve linear fitting [2210.09161], [2305.03577].

## 5. Computational Efficiency and Scalability

The computational cost of ACE scales linearly with the number of atoms $N$, and for fixed maximum body order, the dominant prefactor is set by the number of retained basis functions (often $10^3$--$10^4$). Production-level implementations (PACE, LAMMPS, ACE.jl) attain per-step costs only a few times above classical EAM/MEAM and several orders of magnitude below DFT or neural network potentials. This allows nanosecond-scale molecular dynamics simulations of $10^5$--$10^6$ atoms with DFT-quality accuracy [2303.00753], [2309.11026], [2506.15061].

Active-learning-enabled workflow, robust regularization, and extrapolation grade metrics further facilitate deployment in regimes with incomplete training data, maintaining accuracy and ensuring robustness to extrapolation [2212.08716], [2406.14306], [2506.15061].

## 6. Limitations and Outlook

The main limitations of current ACE approaches include the potential combinatorial growth of feature sets at high body-order or angular truncation, leading to increased parameter count and possible numerical ill-conditioning. Recent work on orthogonal, non-self-interacting bases, basis sparsification using D-optimality and active learning, and smoothness-aware regularization partially mitigates these issues [2401.01550], [2212.08716], [2406.14306].

In ionic, highly polarizable, or long-range correlated materials, nonlocal physics such as Coulomb or van der Waals interactions must either be included via model extension (e.g., explicit charge models, additive D3 dispersion) or supplemented by hybrid approaches [2411.04062], [2505.17991], [2506.15061]. For high-accuracy modeling of strongly correlated or open-shell systems, direct integration with quantum Hamiltonian parameterization or wavefunction techniques is ongoing.

ACE’s linear, universally complete, and symmetry-adapted methodology now provides a foundation for both general-purpose and highly specialized machine-learned potentials in chemistry, physics, and materials science. Ongoing research targets deeper integration with equivariant graph neural architectures, automated and on-the-fly database construction, multi-scale coupling, and robust extension to open quantum systems and reactive chemistries [2308.06462], [2212.08716], [2411.04062].

Source: https://www.emergentmind.com/topics/atomic-cluster-expansion-ace