[論文レビュー] 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

ひとこと要約

本稿では、周期的および自由空間設定下での2次元ヤクビ型ポテンシャル $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.

研究の動機と目的

• 2次元ヤクビ型ポテンシャル $K_0(\alpha r)$ 及びその微分 $K_1(\alpha r)$ を含む和を、効率的かつスペクトル的に正確に計算する手法の開発。
• 周期的および自由空間境界条件の両方における2次元修正ヘルムホルツ方程式に、スペクトル的Ewald法を拡張すること。
• 実空間およびフーリエ空間の和のための打ち切り誤差推定を導出し、最適なパrameter選択を支援すること。
• 特に2次元FFTを用いる応用（例：熱方程式の解法）に適した、均一グリッド上での高速評価を可能とすること。

提案手法

• 周期的および自由空間の両設定において、$K_0(\alpha r)$ 及び $K_1(\alpha r)$ のEwald分解を実施し、収束が速い実空間成分とフーリエ空間成分に分離する。
• FFTを用いて実空間およびk空間の和を $O(N\log N)$ 時間で計算するスペクトル的Ewald法を実装する。
• $\alpha$ が小さい場合の自由空間設定における近似特異性に対処するため、修正されたグリーン関数 $\tilde{G}_F$ 及び $\tilde{H}_F$ を導入する。
• 実空間およびk空間の和のための打ち切り誤差推定を導出し、カットオフ半径および $\xi$ パrameterの選択を支援する。
• 構造化されたグリッドアクセスとFFTを活用することで、グリッド上に配置されたターゲット点に対して最適化し、計算コストを低減する。

実験結果

リサーチクエスチョン

• RQ1 Ewald和分法を用いて、周期的および自由空間設定下での2次元ヤクビ型ポテンシャルおよびその微分をどのようにスペクトル的に分解できるか。
• RQ2 実空間およびk空間の和のための最適なパラメータは何か。これにより打ち切り誤差を最小限に抑えつつ $O(N\log N)$ の計算量を維持できるか。
• RQ3 $\alpha \to 0$ のときの自由空間グリーン関数における近似特異性に対処するには、スペクトル的Ewald法をどのように適合すればよいか。
• RQ4 $K_0$ 及び $K_1$ のスペクトル的Ewald法におけるパラメータ選択を支援する信頼性のある誤差推定は何か。
• 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$ が大きい場合、グリーン関数は指数的に減衰するため、$\sqrt{\pi/(2\alpha\tilde{r})} e^{-\alpha\tilde{r}} \leq \epsilon$ を満たす単純なカットオフ半径 $\tilde{r}$ を用いることができるが、これにより計算量は $O(N^2)$ に増加する。
• グリッド上に配置されたターゲット点に対しては、構造化されたグリッドアクセスとFFTを活用することで、さらに高速化が達成され、特に2次元FFTを用いた熱方程式の解法などの問題において極めて効率的である。