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
- A fast algorithm for particle simulations
- A Fourier-series-based kernel-independent fast multipole method
- A free-space adaptive FMM-based PDE solver in three dimensions
- A Generalized Fast Multipole Method for Nonoscillatory Kernels
- 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
- scientific article; zbMATH DE number 3273551 (Why is no real title available?)
- scientific article; zbMATH DE number 961607 (Why is no real title available?)
- Multipole Translation Theory for the Three-Dimensional Laplace and Helmholtz Equations
- On the fast multipole method for computing the energy of periodic assemblies of charged and dipolar particles
- On the Rokhlin-Greengard method with vortex blobs for problems posed in all space or periodic in one direction
- 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
(23)- Multipole Ewald sums for the fast multipole method
- Spline approximation of ``effective potentials under periodic boundary conditions
- Efficient methods for handling long-range forces in particle-particle simulations
- Accurate and efficient computation of nonlocal potentials based on Gaussian-sum approximation
- Flexibly imposing periodicity in kernel independent FMM: a multipole-to-local operator approach
- Optimized periodic 1/\(r\) Coulomb potential in two dimensions
- Fast multipole method applied to Lagrangian simulations of vortical flows
- Laplace Green's functions for infinite ground planes with local roughness
- Spectrally-accurate numerical method for acoustic scattering from doubly-periodic 3D multilayered media
- Fast multipole boundary element method for the acoustic analysis of finite periodic structures
- HSMA: an O(N) electrostatics package implemented in LAMMPS
- An NFFT based approach to the efficient computation of dipole-dipole interactions under various periodic boundary conditions
- Domain decomposition for quasi-periodic scattering by layered media via robust boundary-integral equations at all frequencies
- Harmonic Surface Mapping Algorithm for Electrostatic Potentials in an Atomistic/Continuum Hybrid Model for Electrolyte Solutions
- A High-Accurate Fast Poisson Solver Based on Harmonic Surface Mapping Algorithm
- Summation of series of harmonic functions
- 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
- Recent Advances in Acoustic Boundary Element Methods
- Spectrally accurate fast summation for periodic Stokes potentials
- Fast Fourier transform periodic interpolation method for superposition sums in a periodic unit cell
- On the fast multipole method for computing the energy of periodic assemblies of charged and dipolar particles
- Fast ewald summation based on NFFT with mixed periodicity
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)
