Fast Fourier Transform periodic interpolation method for superposition sums in a periodic unit cell

Fangzhou Ai, Vitaliy Lomakin
{"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.
周期单元格中叠加和的快速傅里叶变换周期插值方法
我们提出了一种基于快速傅立叶变换的周期插值方法(FFT-PIM),这是一种灵活且计算效率高的方法,用于计算由无限周期阵列的单元格中的叠加和给出的标量势。在相同的框架下,FFT-PIM允许计算动态和静态问题的1D、2D和3D周期性的潜力,包括有或没有周期性相移的问题。对于$N$空间重合的源和观测点,FFT-PIM的计算复杂度为$O(N \log N)$。FFT-PIM使用格林函数的快速收敛级数表示作为叠加和的核。基于这些表示,FFT-PIM将势分解为近区分量,其中包括感兴趣的单元格周围的少量图像,以及远区分量,其中包括无限数量图像的其余部分。通过将非均匀源投影到稀疏的均匀网格上,在稀疏网格上进行叠加和,并将均匀网格的势插值到非均匀观测点,来评估远区分量。使用基于fft的方法评估近区分量,该方法适用于有效处理周期单元胞内的非均匀源-观测器分布。FFT-PIM可用于广泛的应用,例如涉及计算电磁和声学积分方程的周期性问题,微磁解算器和密度泛函理论解算器。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
自引率
0.00%
发文量
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信