---
title: Most Reliable Independent Basis in ACE
url: https://www.emergentmind.com/topics/most-reliable-independent-basis-mrib
type: topic
---

# Most Reliable Independent Basis in ACE

The Most Reliable Independent Basis (MRIB) is an analytically derived, block-wise independent, and conjecturally complete basis of rotation- and permutation-invariant (RPI) cluster functions for the Atomic Cluster Expansion (ACE). Its construction leverages permutation-adapted arrangements of cluster basis functions and recursion properties of generalized Wigner symbols. The MRIB provides a minimal and robust foundation for constructing interatomic potentials that maintain the symmetries and completeness required in modern atomistic modeling, addressing limitations of prior lexicographically ordered or numerically SVD-pruned bases [2208.01756].

## 1. Permutation-Adapted Rotation- and Permutation-Invariant Cluster Functions

The ACE framework describes atomic local environments in terms of basis functions that must be symmetrized with respect to both rotations and permutations. The MRIB is defined by systematically constructing RPI cluster functions as follows:

- **Single-bond basis functions:**  
  $$
  \phi_{nlm}(r_{ij}) = R_n(r_{ij}) Y_l^m(\hat{r}_{ij})
  $$
  $$
  A_{nlm}(i) = \langle \rho_i, \phi_{nlm} \rangle = \sum_j \phi_{nlm}(r_{ij})
  $$
- **Un-symmetrized cluster products:**  
  $$
  \Phi_{\vec{n}\vec{l}\vec{m}}(i) = \prod_{\kappa=1}^N A_{n_\kappa l_\kappa m_\kappa}(i)
  $$
- **Permutation invariance:**  
  $$
  \bar\Phi_{\vec{n}\vec{l}\vec{m}} = \frac{1}{\sqrt{N!}} \sum_{\sigma \in S_N} \Phi_{\sigma(\vec{n}\vec{l}\vec{m})}
  $$
  Each block comprises all functions with the same (multi)set $(\vec{n},\vec{l})$.
- **Rotation invariance via Wigner symbols:**  
  $$
  B_{\vec{n}\vec{l}\vec{L}}(i) = \sum_{\vec{m}} W_{\vec{l}}^{\vec{m}}(\vec{L};\vec{M}) \bar\Phi_{\vec{n}\vec{l}\vec{m}}(i)
  $$
  where $\vec{L}$ are intermediate angular momenta satisfying generalized triangle (“polygon”) conditions and $L_R = 0$ yields full invariance.

These constructions were previously known to form highly over-complete sets, necessitating numerical SVD to identify a minimal basis.

## 2. Analytical Linear Relationships via Recursion and Symmetry

The MRIB methodology eliminates the need for numerical SVD by exploiting analytical linear relations among RPI functions, derived from permutation symmetries and generalized Wigner symbol recursions:

- **Permutation symmetries:**  
  Binary-tree automorphism groups $G_N$—for example, $G_4$ for rank-4 clusters—induce sign-altered permutations of intermediate labels in Wigner-symbol coupled functions.
- **Ladder recursion relations:**  
  Varying one intermediate angular momentum or projection in the Wigner-3j and higher-order coupling coefficients generates analytic raising/lowering relationships among cluster basis elements:
  $$
  W_{\vec{l}}^{\vec{m}}(\dots,L_k) = f_k \cdot W_{\vec{l}}^{\vec{m}}(\dots,L_k+1)
  $$
- **Explicit intra-block linearities:**  
  For example, in rank 4 with $\vec{l}=(1,1,1,1), \vec{n}=(n,n,n,n)$,
  $$
  B_2 = \frac{2}{5\sqrt{5}} B_0
  $$
  Generalization to higher rank exploits the same families of recursion to express any within-block $B$ as a linear combination of a selected independent subset.

## 3. Block-Wise Selection Algorithm for MRIB Construction

The MRIB selects, from each over-complete block, a set of RPI functions provably linearly independent via a block-wise construction algorithm:

**Step 1: Block construction and permutation adaptation**
- For each $(\vec{n},\vec{l})$, compute frequency partition $P_f(\vec{l})$
- Form coupling-compatible partition $P_{fc}(\vec{l})$ for maximal automorphisms
- Generate representative permutations $\sigma_{fc}(\vec{l}) = \vec{l}_{fc}$ and all distinguishable $\sigma(\vec{n})$ up to automorphisms

**Step 2: Block-internal sampling**
- Derive “ladder-order” $F_i^{OC}$ of over-complete block
- Apply precomputed sampling (e.g., fixed stride or index sublist), guaranteed by analytic relations to select a maximal independent subset $F_a^{PA}$

**Pseudocode (verbatim):**
```python
PA_set = []
for each distinct freq-adapted (\vec{n}, \vec{l}_{fc}):
    OC_block = all intermediates \vec{L} satisfying triangle(\vec{l}_{fc}, \vec{L})
    PA_block = []
    for B in OC_block in "ladder order" F_i^{OC}(P_f(\vec{n}),P_f(\vec{l})):
        if B not in span(PA_block):  # by analytic test from recursions
            PA_block.append(B)
    PA_set.extend(PA_block)
```
In practice, sampling patterns are precomputed per block type rather than tested at runtime.

## 4. Independence and Completeness: Analytical Results and Conjecture

**Block-wise independence:**  
By construction, the selected subset in each block cannot be connected via a single Wigner recursion or permitted permutation; in explicit cases (e.g., rank-4 $\vec{l}=(1111)$ or $(2222)$), dependence arises only between known pairings, and the selection removes exactly one from each dependent pair.  

**Completeness conjecture:**  
Total PA-selected function counts match the group-theoretic SO(3)$\times$S$_N$ representation enumeration:
$$
\prod_k (2l_k+1) = \sum_{L_R} \alpha_{L_R} (2L_R + 1)
$$
up to at least N=5,6 and degree $\leq 12$, always agreeing with numerically SVD-pruned bases. Orthogonality of generalized Wigner symbols implies no “across-block” overlap, but a fully rigorous Gram–Schmidt orthonormality proof remains open.

## 5. Numerical Benchmarks and Descriptor Counts

Extensive numerical enumeration validates the MRIB’s independence and completeness. The block-wise analytic selection yields substantial descriptor count reductions:

| Body order, $l_{max}, n_{max}$ | All permutations | Lexicographic RPI | PA-RPI (MRIB) | Fraction retained    |
| ------------------------------ |:---------------:|:-----------------:|:-------------:|:--------------------:|
| Rank 4, $l_{max}=6, n_{max}=6$ | $>491\,000$     | $121\,054$        | $106\,667$    | $~88\%$              |
| Rank 5, $l_{max}=2, n_{max}=6$ | $311\,040$      | $19\,152$         | $10\,674$     | $~56\%$              |

Per-block lexicographic overshoots of up to 40% are rectified by PA sampling, with MRIB matching SVD-pruned cardinality in every tested instance [Table II–VI, 2208.01756].

A practical case is the construction of a linear ACE interatomic potential for tantalum using FitSNAP and Bayesian compressive sensing. In this application, even under strong $\ell_1$ (density) regularization, identifiable high-degree cluster descriptors from the MRIB basis persist and are critical for achieving optimal energy and force accuracy. This demonstrates that physically meaningful high-order clusters are retained under the MRIB construction.

## 6. MRIB: Advantages, Trade-Offs, and Implications for Machine-Learning Potentials

**Advantages:**
- Analytical block-wise independence eliminates small-singular-value pathologies typical of SVD
- No need for large-order SVD computations on ill-conditioned overlap matrices, ensuring stability at high degree or body order
- Basis is strictly minimal: any further removal breaks rotational and permutation invariance or completeness

**Trade-offs compared to SVD-based bases:**
- Fully reproducible, transparent analytical selection
- Avoids $O(N_{desc}^3)$ cost of SVD
- Upfront derivation of permutation and recursion needed
- For large $N$, generation of Wigner symbols and automorphisms required; these scale as $O(N^2)$ per block

**Implications for sparse regression:**
MRIB’s minimal basis alleviates descriptor multicollinearity, facilitating sparser and more stable regression solutions. Bayesian compressive sensing can efficiently identify the most informative MRIB descriptors, and critically, high-degree components often survive regularization, contributing meaningfully to model accuracy.  

In summary, the MRIB—equivalent to the PA-RPI analytic basis—provides a theoretically grounded, fully analytic, and block-wise independent foundation for ACE descriptor sets. Empirical and theoretical consistency with SVD-pruned approaches, but with improved stability and reproducibility especially in high-rank/high-degree regimes, make it a robust platform for constructing machine-learned interatomic potentials [2208.01756].

Source: https://www.emergentmind.com/topics/most-reliable-independent-basis-mrib