{"title":"Fast Fourier Transform periodic interpolation method for superposition sums in a periodic unit cell","authors":"Fangzhou Ai, Vitaliy Lomakin","doi":"arxiv-2312.02376","DOIUrl":null,"url":null,"abstract":"We propose a Fast Fourier Transform based Periodic Interpolation Method\n(FFT-PIM), a flexible and computationally efficient approach for computing the\nscalar potential given by a superposition sum in a unit cell of an infinitely\nperiodic array. Under the same umbrella, FFT-PIM allows computing the potential\nfor 1D, 2D, and 3D periodicities for dynamic and static problems, including\nproblems with and without a periodic phase shift. The computational complexity\nof the FFT-PIM is of $O(N \\log N)$ for $N$ spatially coinciding sources and\nobserver points. The FFT-PIM uses rapidly converging series representations of\nthe Green's function serving as a kernel in the superposition sum. Based on\nthese representations, the FFT-PIM splits the potential into its near-zone\ncomponent, which includes a small number of images surrounding the unit cell of\ninterest, and far-zone component, which includes the rest of an infinite number\nof images. The far-zone component is evaluated by projecting the non-uniform\nsources onto a sparse uniform grid, performing superposition sums on this\nsparse grid, and interpolating the potential from the uniform grid to the\nnon-uniform observation points. The near-zone component is evaluated using an\nFFT-based method, which is adapted to efficiently handle non-uniform\nsource-observer distributions within the periodic unit cell. The FFT-PIM can be\nused for a broad range of applications, such as periodic problems involving\nintegral equations in computational electromagnetic and acoustic, micromagnetic\nsolvers, and density functional theory solvers.","PeriodicalId":501061,"journal":{"name":"arXiv - CS - Numerical Analysis","volume":" 9","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2023-12-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - CS - Numerical Analysis","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2312.02376","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
We propose a Fast Fourier Transform based Periodic Interpolation Method
(FFT-PIM), a flexible and computationally efficient approach for computing the
scalar potential given by a superposition sum in a unit cell of an infinitely
periodic array. Under the same umbrella, FFT-PIM allows computing the potential
for 1D, 2D, and 3D periodicities for dynamic and static problems, including
problems with and without a periodic phase shift. The computational complexity
of the FFT-PIM is of $O(N \log N)$ for $N$ spatially coinciding sources and
observer points. The FFT-PIM uses rapidly converging series representations of
the Green's function serving as a kernel in the superposition sum. Based on
these representations, the FFT-PIM splits the potential into its near-zone
component, which includes a small number of images surrounding the unit cell of
interest, and far-zone component, which includes the rest of an infinite number
of images. The far-zone component is evaluated by projecting the non-uniform
sources onto a sparse uniform grid, performing superposition sums on this
sparse grid, and interpolating the potential from the uniform grid to the
non-uniform observation points. The near-zone component is evaluated using an
FFT-based method, which is adapted to efficiently handle non-uniform
source-observer distributions within the periodic unit cell. The FFT-PIM can be
used for a broad range of applications, such as periodic problems involving
integral equations in computational electromagnetic and acoustic, micromagnetic
solvers, and density functional theory solvers.