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. This work solves the long-standing issue of the numerical computation of the resulting energies in $d$-dimensional lattice systems. 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 three dimensions, 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. We discuss techniques for numerically computing the arising singular integrals and compare the accuracy of our results against computable special cases and against direct summation in low dimensions, 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 both the numerical and analytical foundation for an ongoing investigation into the influence of many-body interactions on the stability of matter.
From: Andreas Alexander Buchheit [view email]
[v1]
Wed, 16 Apr 2025 11:30:11 UTC (1,733 KB)
[v2]
Fri, 12 Jun 2026 14:07:30 UTC (871 KB)
此内容由惯性聚合(RSS阅读器)自动聚合整理，仅供阅读参考。 原文来自 — 版权归原作者所有。