[论文解读] Spectrally accurate Ewald summation for the Yukawa potential in two dimensions

Sara Pålsson, Anna‐Karin Tornberg|arXiv (Cornell University)|Nov 12, 2019

Electromagnetic Scattering and Analysis被引用 2

一句话总结

本文提出了一种针对二维 Yukawa 势 $K_0(\alpha r)$ 及其导数 $K_1(\alpha r)$ 的谱精度 Ewald 求和方法，适用于周期性和自由空间情形。通过采用优化参数和截断误差估计的谱 Ewald 公式，该方法实现了 $O(N\log N)$ 的复杂度，从而能够快速、高精度地计算修正赫姆霍兹方程的边界积分和。

ABSTRACT

An Ewald decomposition of the two-dimensional Yukawa potential and its derivative is presented for both the periodic and the free-space case. These modified Bessel functions of the second kind of zeroth and first degrees are used e.g. when solving the modified Helmholtz equation using a boundary integral method. The spectral Ewald method is used to compute arising sums at O(N log N) cost for N source and target points. To facilitate parameter selection, truncation-error estimates are developed for both the real-space sum and the Fourier-space sum, and are shown to estimate the errors well.

研究动机与目标

• 开发一种高效且谱精度的计算方法，用于求解涉及二维 Yukawa 势 $K_0(\alpha r)$ 及其导数 $K_1(\alpha r)$ 的和。
• 将谱 Ewald 方法扩展至周期性和自由空间边界条件下二维修正赫姆霍兹方程的求解。
• 推导截断误差估计，以指导实空间和傅里叶空间求和中参数的最优选择。
• 实现在均匀网格上的快速求值，特别适用于需要二维 FFT 的应用，例如求解热方程。

提出的方法

• 推导 $K_0(\alpha r)$ 和 $K_1(\alpha r)$ 在周期性和自由空间情形下的 Ewald 分解，将其分解为收敛迅速的实空间和傅里叶空间分量。
• 利用 FFT 实现谱 Ewald 方法，以 $O(N\log N)$ 时间计算实空间和 $k$-空间的求和。
• 引入修正的格林函数 $\tilde{G}_F$ 和 $\tilde{H}_F$，以处理自由空间情形下 $\alpha$ 较小时的近奇点问题。
• 为实空间和 $k$-空间求和推导截断误差估计，以指导截断半径和 $\xi$ 参数的选择。
• 通过利用结构化网格访问和 FFT，对网格上的目标点进行优化，从而降低计算成本。

实验结果

研究问题

• RQ1 如何在周期性和自由空间情形下，利用 Ewald 求和对二维 Yukawa 势及其导数实现谱分解？
• RQ2 实空间和 $k$-空间求和的最优参数是什么，以在保持 $O(N\log N)$ 复杂度的同时最小化截断误差？
• RQ3 谱 Ewald 方法如何适应自由空间格林函数在 $\alpha \to 0$ 时的近奇点问题？
• RQ4 对于谱 Ewald 方法中 $K_0$ 和 $K_1$ 的参数选择，哪些误差估计是可靠的？
• RQ5 当目标点位于均匀网格上时，该方法是否可以进一步加速？

主要发现

• 谱 Ewald 方法在周期性和自由空间情形下对 $K_0(\alpha r)$ 和 $K_1(\alpha r)$ 的求和均实现了 $O(N\log N)$ 的复杂度。
• 实空间和 $k$-空间求和的截断误差估计被证明能准确预测实际误差，从而实现参数的最优选择。
• 当 $\alpha$ 较小时，必须使用修正的格林函数 $\tilde{G}_F$ 和 $\tilde{H}_F$ 以避免傅里叶空间求和中的近奇点，当 $\alpha L / (2\pi) \lesssim 1.5$ 时需切换至该形式。
• 当 $\alpha$ 较大时，格林函数呈指数衰减，可使用简单的截断半径 $\tilde{r}$，满足 $\sqrt{\pi/(2\alpha\tilde{r})} e^{-\alpha\tilde{r}} \leq \epsilon$，但此方法会导致 $O(N^2)$ 的计算成本。
• 当目标点位于网格上时，通过利用结构化网格访问和 FFT，该方法可实现进一步加速，使其在求解热方程等需使用二维 FFT 的问题中极为高效。