期刊
SIAM JOURNAL ON SCIENTIFIC COMPUTING
卷 36, 期 5, 页码 B777-B794出版社
SIAM PUBLICATIONS
DOI: 10.1137/130945582
关键词
Coulomb interaction; dipole-dipole interaction; interaction energy; nonuniform FFT; nonlocal; Poisson equation
资金
- National Science Foundation [CCF-0905395]
- U.S. Department of Energy [DEFG0288ER25053]
- Singapore A*STAR SERC PSF-grant [1321202067]
We present a fast and accurate algorithm for the evaluation of nonlocal (long-range) Coulomb and dipole-dipole interactions in free space. The governing potential is simply the convolution of an interaction kernel U(x) and a density function rho(x) = vertical bar psi(x)vertical bar(2) for some complex-valued wave function psi(x), permitting the formal use of Fourier methods. These are hampered by the fact that the Fourier transform of the interaction kernel (U) over cap (k) has a singularity and/or (rho) over cap (k) not equal 0 at the origin k = 0 in Fourier (phase) space. Thus, accuracy is lost when using a uniform Cartesian grid in k which would otherwise permit the use of the FFT for evaluating the convolution. Here, we make use of a high-order discretization of the Fourier integral, accelerated by the nonuniform fast Fourier transform (NUFFT). By adopting spherical and polar phase-space discretizations in three and two dimensions, respectively, the singularity in (U) over cap (k) at the origin is canceled so that only a modest number of degrees of freedom are required to evaluate the Fourier integral, assuming that the density function (rho x) is smooth and decays sufficiently quickly as vertical bar x vertical bar -> infinity. More precisely, the calculation requires O(N logN) operations, where N is the total number of discretization points in the computational domain. Numerical examples are presented to demonstrate the performance of the algorithm.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据