Abstract: In a number of problems in computational physics, a finite sum of kernel functions centered at particle locations located in a box in three dimensions must be extended by imposing periodic boundary conditions on box boundaries. Even though the finite sum can be efficiently computed via fast summation algorithms, such as the fast multipole method (FMM), the periodized extension is usually treated via a different algorithm, Ewald summation, accelerated via the fast Fourier transform (FFT). A different approach to compute this periodized sum just using a blackbox finite fast summation algorithm is presented in this paper. The method splits the periodized sum in to two parts. The first, comprising the contribution of all points outside a large sphere enclosing the box, and some of its neighbors, is approximated inside the box by a collection of kernel functions ("sources") placed on the surface of the sphere or using an expansion in terms of spectrally convergent local basis functions. The second part, comprising the part inside the sphere, and including the box and its immediate neighborhood, is treated via available summation algorithms. The coefficients of the sources are determined by least squares collocation of the periodicity condition of the total potential, imposed on a circumspherical surface for the box. While the method is presented in general, details are worked out for the case of evaluating electrostatic potentials and forces. Results show that when used with the FMM, the periodized sum can be computed to any specified accuracy, at an additional cost of the order of the free-space FMM. Several technical details and efficient algorithms for auxiliary computations are provided, as are numerical comparisons.
Recommendations
Cites work
- scientific article; zbMATH DE number 3273551 (Why is no real title available?)
- scientific article; zbMATH DE number 961607 (Why is no real title available?)
- A Fourier-series-based kernel-independent fast multipole method
- A Generalized Fast Multipole Method for Nonoscillatory Kernels
- A fast algorithm for particle simulations
- A free-space adaptive FMM-based PDE solver in three dimensions
- A kernel independent fast multipole algorithm for radial basis functions
- A kernel-independent adaptive fast multipole algorithm in two and three dimensions
- A multipole-based algorithm for efficient calculation of forces and potentials in macroscopic periodic assemblies of particles
- A new integral representation for quasi-periodic fields and its application to two-dimensional band structure calculations
- A new integral representation for quasi-periodic scattering problems in two dimensions
- A periodic FMM for Maxwell's equations in 3D and its applications to problems related to photonic crystals
- Fast directional multilevel summation for oscillatory kernels based on Chebyshev interpolation
- Fast electrostatic force calculation on parallel computer clusters
- Fast multipole method for the biharmonic equation in three dimensions
- Fast multipole methods on graphics processors
- Multipole Translation Theory for the Three-Dimensional Laplace and Helmholtz Equations
- On the Rokhlin-Greengard method with vortex blobs for problems posed in all space or periodic in one direction
- On the fast multipole method for computing the energy of periodic assemblies of charged and dipolar particles
- Radial Basis Functions
- Spectral accuracy in fast Ewald-based methods for particle simulations
- Spectrally accurate fast summation for periodic Stokes potentials
- The black-box fast multipole method
- The distribution of points on the sphere and corresponding cubature formulae
Cited in
(21)- Fast Fourier transform periodic interpolation method for superposition sums in a periodic unit cell
- Domain decomposition for quasi-periodic scattering by layered media via robust boundary-integral equations at all frequencies
- Summation of series of harmonic functions
- Efficient methods for handling long-range forces in particle-particle simulations
- Fast multipole boundary element method for the acoustic analysis of finite periodic structures
- An NFFT based approach to the efficient computation of dipole-dipole interactions under various periodic boundary conditions
- A High-Accurate Fast Poisson Solver Based on Harmonic Surface Mapping Algorithm
- Spectrally accurate fast summation for periodic Stokes potentials
- Harmonic Surface Mapping Algorithm for Electrostatic Potentials in an Atomistic/Continuum Hybrid Model for Electrolyte Solutions
- Optimized periodic 1/\(r\) Coulomb potential in two dimensions
- Laplace Green's functions for infinite ground planes with local roughness
- Fast ewald summation based on NFFT with mixed periodicity
- Spectrally-accurate numerical method for acoustic scattering from doubly-periodic 3D multilayered media
- Accurate and efficient computation of nonlocal potentials based on Gaussian-sum approximation
- Grid-based lattice summation of electrostatic potentials by assembled rank-structured tensor approximation
- Efficient numerical solution of acoustic scattering from doubly-periodic arrays of axisymmetric objects
- On the fast multipole method for computing the energy of periodic assemblies of charged and dipolar particles
- Fast multipole method applied to Lagrangian simulations of vortical flows
- Flexibly imposing periodicity in kernel independent FMM: a multipole-to-local operator approach
- HSMA: an \(O(N)\) electrostatics package implemented in LAMMPS
- Recent Advances in Acoustic Boundary Element Methods
This page was built for publication: A method to compute periodic sums
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q349339)
