Abstract
We present a comprehensive analysis of an algorithm for evaluating high-dimensional polynomials that are invariant (or equi-variant) under permutations and rotations. This task arises in the evaluation of linear models as well as equivariant neural network models of many-particle systems. The theoretical bottleneck is the contraction of a high-dimensional symmetric and sparse tensor with a specific sparsity pattern that is directly related to the symmetries imposed on the polynomial. The sparsity of this tensor makes it challenging to construct a highly efficient evaluation scheme. Bachmayr et al. (“Polynomial approximation of symmetric functions,” Math. Comp., vol. 93, pp. 811–839, 2024) and Lysogorskiy et al. (“Performant implementation of the atomic cluster expansion (pace): application to copper and silicon,” npj Comput. Mater., vol. 7, Art. no. 97, 2021) introduced a recursive evaluation strategy that relied on a number of heuristics, but performed well in tests. In the present work, we propose an explicit construction of such a recursive evaluation strategy and show that it is in fact optimal in the limit of infinite polynomial degree.
Funding source: NSERC Discovery Grant
Award Identifier / Grant number: GR019381
Funding source: NFRF Exploration Grant
Award Identifier / Grant number: GR022937
-
Research ethics: Not applicable.
-
Informed consent: Not applicable.
-
Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.
-
Use of Large Language Models, AI and Machine Learning Tools: None declared.
-
Conflict of interest: The authors state no conflict of interest.
-
Research funding: This work was supported by NSERC Discovery Grant GR019381; NFRF Exploration Grant GR022937.
-
Data availability: Not applicable.
References
[1] J. Behler and M. Parrinello, “Generalized neural-network representation of high-dimensional potential-energy surfaces,” Phys. Rev. Lett., vol. 98, no. 14, Art. no. 146401, 2007. https://doi.org/10.1103/physrevlett.98.146401.Search in Google Scholar
[2] A. Bartók, M. Payne, R. Kondor, and G. Csányi, “Gaussian approximation potentials: the accuracy of quantum mechanics, without the electrons,” Phys. Rev. Lett., vol. 104, no. 13, Art. no. 136403, 2010. https://doi.org/10.1103/physrevlett.104.136403.Search in Google Scholar
[3] Y. Zuo et al.., “Performance and cost assessment of machine learning interatomic potentials,” J. Phys. Chem. A, vol. 124, no. 4, pp. 731–745, 2020. https://doi.org/10.1021/acs.jpca.9b08723.Search in Google Scholar PubMed
[4] F. Musil, A. Grisafi, A. P. Bartók, C. Ortner, G. Csányi, and M. Ceriotti, “Physics-Inspired structural representations for molecules and materials,” Chem. Rev., vol. 121, no. 16, pp. 9759–9815, 2021. https://doi.org/10.1021/acs.chemrev.1c00021.Search in Google Scholar PubMed
[5] B. Braams and J. Bowman, “Permutationally invariant potential energy surfaces in high dimensionality,” Int. Rev. Phys. Chem., vol. 28, no. 4, pp. 577–606, 2009. https://doi.org/10.1080/01442350903234923.Search in Google Scholar
[6] A. Shapeev, “Moment tensor potentials: a class of systematically improvable interatomic potentials,” Multiscale Model. Simul., vol. 14, no. 3, pp. 1153–1173, 2016. https://doi.org/10.1137/15m1054183.Search in Google Scholar
[7] R. Drautz, “Atomic cluster expansion for accurate and transferable interatomic potentials,” Phys. Rev. B, vol. 99, Art. no. 014104, 2019, https://doi.org/10.1103/physrevb.99.014104.Search in Google Scholar
[8] G. Dusson et al.., “Atomic cluster expansion: completeness, efficiency and stability,” J. Comput. Phys., vol. 454, Art. no. 110946, 2022, https://doi.org/10.1016/j.jcp.2022.110946.Search in Google Scholar
[9] Y. Lysogorskiy et al.., “Performant implementation of the atomic cluster expansion (pace): application to copper and silicon,” npj Comput. Mater., vol. 7, Art. no. 97, 2021, https://doi.org/10.1038/s41524-021-00559-9.Search in Google Scholar
[10] A. Seko, A. Togo, and I. Tanaka, “Group-theoretical high-order rotational invariants for structural representations: application to linearized machine learning interatomic potential,” Phys. Rev. B Condens. Matter, vol. 99, no. 21, Art. no. 214108, 2019. https://doi.org/10.1103/physrevb.99.214108.Search in Google Scholar
[11] J. Nigam, S. Pozdnyakov, and M. Ceriotti, “Recursive evaluation and iterative contraction of n-body equivariant features,” J. Chem. Phys., vol. 153, no. 12, Art. no. 121101, 2020. https://doi.org/10.1063/5.0021116.Search in Google Scholar PubMed
[12] J. Nigam, S. Pozdnyakov, G. Fraux, and M. Ceriotti, “Unified theory of atom-centered representations and message-passing machine-learning schemes,” J. Chem. Phys., vol. 156, Art. no. 204115, 2022.10.1063/5.0087042Search in Google Scholar
[13] B. Anderson, T.-S. Hy, and R. Kondor, “Cormorant: covariant molecular neural networks,” arXiv:1906.04015, 2019.Search in Google Scholar
[14] N. Thomas et al.., “Tensor field networks: rotation- and translation-equivariant neural networks for 3D point clouds,” arXiv:1802.08219, 2018.Search in Google Scholar
[15] M. Bachmayr, G. Dusson, and C. Ortner, “Polynomial approximation of symmetric functions,” Math. Comp., vol. 93, pp. 811–839, 2024.10.1090/mcom/3868Search in Google Scholar
[16] L. Zhang et al.., “Equivariant analytical mapping of first principles Hamiltonians to accurate and transferable materials models,” npj Comp. Mater., vol. 8, 2022.10.1038/s41524-022-00843-2Search in Google Scholar
[17] A. G. Beged-Dov, “Lower and upper bounds for the number of lattice points in a simplex,” SIAM J. Appl. Math., vol. 22, no. 1, pp. 106–108, 1972. https://doi.org/10.1137/0122012.Search in Google Scholar
© 2024 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Frontmatter
- Stabilized weighted reduced order methods for parametrized advection-dominated optimal control problems governed by partial differential equations with random inputs
- Augmenting the grad-div stabilization for Taylor–Hood finite elements with a vorticity stabilization
- Obtaining higher-order Galerkin accuracy when the boundary is polygonally approximated
- Optimal evaluation of symmetry-adapted n-correlations via recursive contraction of sparse symmetric tensors
Articles in the same Issue
- Frontmatter
- Stabilized weighted reduced order methods for parametrized advection-dominated optimal control problems governed by partial differential equations with random inputs
- Augmenting the grad-div stabilization for Taylor–Hood finite elements with a vorticity stabilization
- Obtaining higher-order Galerkin accuracy when the boundary is polygonally approximated
- Optimal evaluation of symmetry-adapted n-correlations via recursive contraction of sparse symmetric tensors