Andreas A. Buchheit,Torsten Keßler,Kirill Serkh
Abstract:This paper introduces a new method for the efficient computation of oscillatory multidimensional lattice sums in geometries with boundaries. Such sums are ubiquitous in both pure and applied mathematics, and have immediate applications in condensed matter physics and topological quantum physics. The challenge in their evaluation results from the combination of singular long-range interactions with the loss of translational invariance caused by the boundaries, rendering standard tools ineffective. Our work shows that these lattice sums can be generated from a generalization of the Riemann zeta function to multidimensional non-periodic lattice sums. We put forth a new representation of this zeta function together with a numerical algorithm that ensures exponential convergence across an extensive range of geometries. Notably, our method's runtime is influenced only by the complexity of the considered geometries and not by the number of particles, providing the foundation for efficient simulations of macroscopic condensed matter systems. We showcase the practical utility of our method by computing interaction energies in a three-dimensional crystal structure with $3\times 10^{23}$ particles. Our method's accuracy is demonstrated through extensive numerical experiments. A reference implementation is provided online along with this article.
What problem does this paper attempt to address?
This paper attempts to solve the problem of efficiently calculating oscillating multi - dimensional lattice sums in geometric structures with boundaries. Specifically, the paper focuses on calculating lattice sums in the absence of translational invariance, which has important applications in fields such as condensed - matter physics and topological quantum physics.
### Main problems
1. **Long - range interactions and boundary effects**: When considering systems with boundaries, long - range interactions (such as Coulomb interactions, dipole - dipole interactions, etc.) become complicated. The presence of boundaries breaks the translational invariance of the system, making standard computational tools (such as the Ewald summation method) ineffective.
2. **Lack of efficient computational methods**: For macroscopic systems containing a large number of particles (for example, \(3\times10^{23}\) particles), existing methods are difficult to efficiently calculate lattice sums, especially in three - dimensional space and when the interaction exponent is close to the system dimension (such as dipole - dipole interactions).
### Core contributions of the paper
1. **Generalization of the Riemann Zeta function**: The paper introduces a new method, expressing the lattice sum in the form of a generalized Riemann Zeta function, and proposes a new representation method to ensure exponential convergence in a wide range of geometric structures.
2. **Proposing a new numerical algorithm**: The author develops a new numerical algorithm. The running time of this algorithm depends only on the complexity of the geometric structure under consideration, not on the number of particles. This makes it possible to simulate macroscopic systems.
3. **Demonstration of practical applications**: By calculating the interaction energy in a three - dimensional crystal structure containing \(3\times10^{23}\) particles, the practical application value of this method is demonstrated. The experimental results show that even on a macroscopic scale, the influence of boundary effects on the system is still significant.
### Formula summary
- **Epstein Zeta function**:
\[
Z_{\Lambda,\nu}(x,y)=\sum_{z\in\Lambda}'e^{-2\pi iy\cdot z}|z - x|^{-\nu},\quad\text{Re}(\nu)>d
\]
where \(\Lambda\) is a lattice, \(x,y\in\mathbb{R}^d\), \(\nu\in\mathbb{C}\).
- **Generalized Crandall formula**:
\[
Z_{L,\nu}(y)=\frac{\pi^{\nu/2}}{\lambda^\nu\Gamma(\nu/2)}\left(\sum_{z\in L}'e^{-2\pi iy\cdot z}G_\nu(z/\lambda)+\lambda^d(\hat{1}_L*G_{d - \nu})(\lambda\cdot)(y)\right)
\]
where \(L\) is a uniform discrete set, \(G_\nu(z)\) is a fast - decaying function, \(\hat{1}_L\) is the Fourier transform of the generalized Dirac comb.
### Conclusion
The paper solves the long - standing problem of efficiently calculating oscillating multi - dimensional lattice sums in geometric structures with boundaries, providing a new theoretical framework and numerical method applicable to any dimension, any interaction exponent, and any lattice subset. This result is of great significance for fields such as condensed - matter physics and quantum physics.