Abstract
When modeling physical properties of molecules with machine learning, it is desirable to incorporate -covariance. While such models based on low body order features are not complete, we formulate and prove general completeness properties for higher order methods, and show that of these features are enough for up to atoms. We also find that the Clebsch--Gordan operations commonly used in these methods can be replaced by matrix multiplications without sacrificing completeness, lowering the scaling from to in the degree of the features. We apply this to quantum chemistry, but the proposed methods are generally applicable for problems involving 3D point configurations.