Accurate and efficient computation of nonlocal potentials based on Gaussian-sum approximation

Abstract We introduce an accurate and efficient method for the numerical evaluation of nonlocal potentials, including the 3D/2D Coulomb, 2D Poisson and 3D dipole–dipole potentials. Our method is based on a Gaussian-sum approximation of the singular convolution kernel combined with a Taylor expansion of the density. Starting from the convolution formulation of the nonlocal potential, for smooth and fast decaying densities, we make a full use of the Fourier pseudospectral (plane wave) approximation of the density and a separable Gaussian-sum approximation of the kernel in an interval where the singularity (the origin) is excluded. The potential is separated into a regular integral and a near-field singular correction integral. The first is computed with the Fourier pseudospectral method, while the latter is well resolved utilizing a low-order Taylor expansion of the density. Both parts are accelerated by fast Fourier transforms (FFT). The method is accurate (14–16 digits), efficient ( O ( N log ⁡ N ) complexity), low in storage, easily adaptable to other different kernels, applicable for anisotropic densities and highly parallelizable.

[1]  Robert Piessens,et al.  Quadpack: A Subroutine Package for Automatic Integration , 2011 .

[2]  Weizhu Bao,et al.  Computing the ground state and dynamics of the nonlinear Schrödinger equation with nonlocal interactions via the nonuniform FFT , 2014, J. Comput. Phys..

[3]  Mark E. Tuckerman,et al.  A reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters , 1999 .

[4]  V. Arkadiev,et al.  Inverse scattering transform method and soliton solutions for davey-stewartson II equation , 1989 .

[5]  L. Greengard,et al.  A new version of the Fast Multipole Method for the Laplace equation in three dimensions , 1997, Acta Numerica.

[6]  Markus Gusenbauer,et al.  Fast stray field computation on tensor grids , 2011, J. Comput. Phys..

[7]  Leslie Greengard,et al.  Fast and Accurate Evaluation of Nonlocal Coulomb and Dipole-Dipole Interactions via the Nonuniform FFT , 2013, SIAM J. Sci. Comput..

[8]  L. You,et al.  Trapped atomic condensates with anisotropic interactions , 2000 .

[9]  Yong Zhang,et al.  Exact Artificial Boundary Condition for the Poisson Equation in the Simulation of the 2D Schrödinger-Poisson System , 2014 .

[10]  Thomas Schrefl,et al.  Non-uniform FFT for the finite element computation of the micromagnetic scalar potential , 2013, J. Comput. Phys..

[11]  Yong Zhang,et al.  On the computation of ground state and dynamics of Schrödinger-Poisson-Slater system , 2011, J. Comput. Phys..

[12]  W. Bao,et al.  MATHEMATICAL THEORY AND NUMERICAL METHODS FOR , 2012 .

[13]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[14]  L. You,et al.  Trapped condensates of atoms with dipole interactions , 2001 .

[15]  Leslie Greengard,et al.  A New Fast-Multipole Accelerated Poisson Solver in Two Dimensions , 2001, SIAM J. Sci. Comput..

[16]  G. Beylkin.  On the Fast Fourier Transform of Functions with Singularities , 1995 .

[17]  Gregory Beylkin,et al.  Fast convolution with the free space Helmholtz Green's function , 2009, J. Comput. Phys..

[18]  Jie Shen,et al.  Spectral and High-Order Methods with Applications , 2006 .

[19]  Martin J. Mohlenkamp,et al.  Numerical operator calculus in higher dimensions , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Hanquan Wang,et al.  Efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates , 2010, J. Comput. Phys..

[21]  Vladimir Rokhlin,et al.  Fast Fourier Transforms for Nonequispaced Data , 1993, SIAM J. Sci. Comput..

[22]  Nicole Fruehauf.  Numerical Methods Based On Sinc And Analytic Functions , 2016 .

[23]  Michael Pippig.  PFFT: An Extension of FFTW to Massively Parallel Architectures , 2013, SIAM J. Sci. Comput..

[24]  Dietrich Braess,et al.  On the efficient computation of high-dimensional integrals and the approximation by exponential sums , 2009 .

[25]  Boris N. Khoromskij,et al.  Low-rank Kronecker-product Approximation to Multi-dimensional Nonlocal Operators. Part II. HKT Representation of Certain Operators , 2005, Computing.

[26]  Boris N. Khoromskij,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig Fast and Accurate Tensor Approximation of Multivariate Convolution with Linear Scaling in Dimension Fast and Accurate Tensor Approximation of Multivariate Convolution with Linear Scaling in Dimension , 2022 .

[27]  Yong Zhang,et al.  On the ground states and dynamics of space fractional nonlinear Schrödinger/Gross-Pitaevskii equations with rotation term and nonlocal nonlinear interactions , 2015, J. Comput. Phys..

[28]  Boris N. Khoromskij,et al.  Low-rank Kronecker-product Approximation to Multi-dimensional Nonlocal Operators. Part I. Separable Approximation of Multi-variate Functions , 2005, Computing.

[29]  D. Kammler.  Approximation with Exponential Sums. , 1979 .

[30]  Weizhu Bao,et al.  Dimension Reduction of the Schrödinger Equation with Coulomb and Anisotropic Confining Potentials , 2013, SIAM J. Appl. Math..

[31]  G. Beylkin,et al.  Approximation by exponential sums revisited , 2010 .

[32]  W. Bao,et al.  Mathematical Models and Numerical Methods for Bose-Einstein Condensation , 2012, 1212.5341.

[33]  Luigi Genovese,et al.  Efficient and accurate solver of the three-dimensional screened and unscreened Poisson's equation with generic boundary conditions. , 2012, The Journal of chemical physics.

[34]  Yong Zhang,et al.  A novel nonlocal potential solver based on nonuniform FFT for efficient simulation of the Davey-Stewartson equations , 2014, 1409.2014.

[35]  Stefan Goedecker,et al.  Efficient solution of Poisson's equation with free boundary conditions. , 2006, The Journal of chemical physics.

[36]  Stefan Goedecker,et al.  Efficient and accurate three-dimensional Poisson solver for surface problems. , 2007, The Journal of chemical physics.

[37]  Weizhu Bao,et al.  Accurate and efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates via the nonuniform FFT , 2015, 1504.02897.

[38]  Peter Pulay,et al.  Accurate molecular integrals and energies using combined plane wave and Gaussian basis sets in molecular electronic structure theory , 2002 .

[39]  Thomas Schrefl,et al.  FFT-based Kronecker product approximation to micromagnetic long-range interactions , 2012, 1212.3509.

[40]  G. Beylkin,et al.  On approximation of functions by exponential sums , 2005 .

[41]  L. Trefethen.  Spectral Methods in MATLAB , 2000 .