Epstein zeta method for many-body lattice sums
Abstract
Many-body interactions arise naturally in the perturbative treatment of classical and quantum many-body systems and play a crucial role in the description of condensed matter systems. In the case of three-body interactions, the Axilrod-Teller-Muto (ATM) potential is highly relevant for the quantitative prediction of material properties. The computation of the resulting energies in d-dimensional lattice systems is challenging, as a high-dimensional lattice sum needs to be evaluated to high precision. This work solves this long-standing issue. We present an efficiently computable representation of many-body lattice sums in terms of singular integrals over products of Epstein zeta functions. For three-body interactions in 3D, this approach reduces the runtime for computing the ATM lattice sum from weeks to minutes. Our approach further extends to a broad class of n-body lattice sums. We demonstrate that the computational cost of our method only increases linearly with n, evading the exponential increase in complexity of direct summation. The evaluation of 51-body interactions on a two-dimensional lattice, corresponding to a 100-dimensional sum, can be performed within seconds on a laptop. We discuss techniques for computing the arising singular integrals and compare the accuracy of our results against computable benchmarks, achieving full precision for exponents greater than the system dimension. Finally, we apply our method to study the stability of a three-dimensional lattice system with Lennard-Jones two-body interactions under the inclusion of an ATM three-body term at finite pressure, finding a transition from the face-centered-cubic to the body-centered-cubic lattice structure with increasing ATM coupling strength. This work establishes the mathematical foundation for an ongoing investigation into the influence of many-body interactions on the stability of matter.