  For the purposes of interatomic potentials, atomic cluster expansion (ACE)
methods provide a generalizable and systematic way to describe atomic
environments of arbitrary body order. As an extension of fixed-lattice cluster
expansions, it may be used to model properties for arbitrary collections of
atoms with or without translational symmetry. This makes them excellent choices
for machine-learned interatomic potentials. Extending cluster expansions
requires construction of a rotation and permutation invariant (RPI) cluster
basis. When these functions are symmetrized, this basis becomes over-complete.
Current approaches for ACE methods involve a numerical step for this reason.
This prevents a full analytical definition of an ACE basis and ACE descriptors,
and can be susceptible to numerical instabilities. In this work, we propose a
method that adapts to this permutation invariance to obtain a rotation and
permutation invariant (PA-RPI) basis that is complete and independent. Through
the presentation of this method, we demonstrate how and why there are linearly
dependent basis functions through some explicit proofs. The proofs as well as
our method leverage the permutation symmetries as well as recursion
relationships of the generalized Wigner symbols. These properties of the
generalized Wigner symbols are used to derive relationships between RPI
functions. Following derivations of these properties, the procedure for the
construction of the PA-RPI basis is presented. Among some specific cases of RPI
basis construction, we provide extensive tabulations of PA-RPI basis function
counts. We also show that in practical applications for interatomic potentials,
high degree descriptors help achieve high-accuracy models. High degree, high
body-order functions can now be selected based on analytical recursion
relationships and used in ACE models rather than relying on numerical SVD.
