Abstract:The quasi-2D electrostatic systems, characterized by periodicity in two dimensions with a free third dimension, have garnered significant interest in many fields. We apply the sum-of-Gaussians (SOG) approximation to the Laplace kernel, dividing the interactions into near-field, mid-range, and long-range components. The near-field component, singular but compactly supported in a local domain, is directly calculated. The mid-range component is managed using a procedure similar to nonuniform fast Fourier transforms in three dimensions. The long-range component, which includes Gaussians of large variance, is treated with polynomial interpolation/anterpolation in the free dimension and Fourier spectral solver in the other two dimensions on proxy points. Unlike the fast Ewald summation, which requires extensive zero padding in the case of high aspect ratios, the separability of Gaussians allows us to handle such case without any zero padding in the free direction. Furthermore, while NUFFTs typically rely on certain upsampling in each dimension, and the truncated kernel method introduces an additional factor of upsampling due to kernel oscillation, our scheme eliminates the need for upsampling in any direction due to the smoothness of Gaussians, significantly reducing computational cost for large-scale problems. Finally, whereas all periodic fast multipole methods require dividing the periodic tiling into a smooth far part and a near part containing its nearest neighboring cells, our scheme operates directly on the fundamental cell, resulting in better performance with simpler implementation. We provide a rigorous error analysis showing that upsampling is not required in NUFFT-like steps, achieving $O(N\log N)$ complexity with a small prefactor. The performance of the scheme is demonstrated via extensive numerical experiments.
What problem does this paper attempt to address?
This paper attempts to address the challenges in calculating electrostatic interactions in quasi - 2D systems, especially problems related to long - range interactions. Specifically, the main research questions include:
1. **Limitations of traditional methods**: The traditional Ewald summation method requires a large amount of zero - padding when dealing with quasi - 2D systems with a high aspect ratio, which leads to high computational costs and low efficiency.
2. **Computational complexity**: Long - range electrostatic interactions in quasi - 2D systems usually have an \(O(N^2)\) complexity, which poses a significant obstacle to large - scale simulations.
3. **Requirement for high precision**: In order to achieve high - precision simulations, a method that can reduce computational costs while ensuring precision is required.
To solve these problems, the paper proposes a fast spectral sum - of - Gaussians (SOG) method. The main innovations and advantages of this method are as follows:
- **Avoid zero - padding and up - sampling**: Unlike the traditional fast Ewald method, the SOG method does not require zero - padding in the free direction, nor does it require up - sampling in any direction, thereby significantly reducing computational costs.
- **Efficiently handle long - range interactions**: By decomposing the interactions into near - field, mid - range, and long - range parts and using different techniques (such as direct calculation, non - uniform fast Fourier transform (NUFFT), polynomial interpolation, etc.) to handle them respectively, the computational challenges of long - range interactions are effectively solved.
- **Simplify implementation and improve performance**: Compared with the periodic fast multipole method (FMM), the SOG method does not need to divide the periodic tiling into far - field and near - field parts, and all calculations are completed within the basic unit, thereby simplifying the implementation and improving efficiency.
- **Spectral precision**: Based on strict error analysis, the SOG method can achieve spectral precision and can reach machine precision with a moderate increase in computational cost.
In summary, this paper aims to develop an efficient, accurate, and easy - to - implement algorithm to address the challenges in calculating electrostatic interactions in quasi - 2D systems, especially suitable for large - scale simulations in fields such as semiconductor physics, new energy technologies, and nanomaterials.