×

On optimal zero-padding of kernel truncation method. (English) Zbl 1533.65261

The authors consider the optimal zero-padding factor of the kernel truncation method (KTM) in the convolution-type nonlocal potential evaluation, and strictly show that the optimal zero-padding factor is \(\sqrt{d} + 1\), \(d\) being the space dimension. The memory costs are significantly reduced to a small fraction, i.e., \(\left(\frac{\sqrt{d}+1}{4}\right)^d\), of that in the fourfold algorithm. The authors deduce that twofold padding is optimal for the one-dimensional problem, whereas threefold padding is sufficient to get spectral accuracy for the two- and three-dimensional problems. For the anisotropic density, the total optimal zero-padding factor grows linearly with the anisotropy strength. The KTM is obtained by the Fourier spectral method, and it can be rewritten as a discrete convolution regardless of the anisotropy strength. Providing that the discrete tensor is available through a precomputation step, the effective computation is simply a pair of fast Fourier transform (FFT) and inverse FFT (iFFT) on a double-sized density. Error estimates are provided for both the nonlocal potential and the density in \(d\) dimensions. Several numerical experiments are carried out to confirm the enhancements in terms of memory and computational cost.

MSC:

65T50 Numerical methods for discrete and fast Fourier transforms
65R10 Numerical methods for integral transforms
65M70 Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs
35Q55 NLS equations (nonlinear Schrödinger equations)

Software:

Matlab

References:

[1] Advanpix Multiprecision Computing Toolbox, https://www.advanpix.com.
[2] Abramowitz, M. and Stegun, I. A., Handbook of Mathematical Functions, Dover, New York, 1965. · Zbl 0171.38503
[3] Antoine, X., Tang, Q., and Zhang, Y., A preconditioned conjugated gradient method for computing ground sates of rotating dipolar Bose-Einstein condensates via Kernel truncation method for dipole-dipole interaction evaluation, Commun. Comput. Phys., 24 (2018), pp. 966-988. · Zbl 1475.65142
[4] Bao, W. and Cai, Y., Mathematical theory and numerical methods for Bose-Einstein condensations, Kinet. Relat. Models, 6 (2013), pp. 1-135. · Zbl 1266.82009
[5] Bao, W., Cai, Y., and Wang, H., Efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates, J. Comput. Phys., 229 (2010), pp. 7874-7892. · Zbl 1198.82036
[6] Bao, W., Jian, H., Mauser, N. J., and Zhang, Y., Dimension reduction of the Schrödinger equation with Coulomb and anisotropic confining potentials, SIAM J. Appl. Math., 73 (2013), pp. 2100-2123. · Zbl 1294.35130
[7] Bao, W., Tang, Q., and Zhang, Y., Accurate and efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates via the nonuniform FFT, Commun. Comput. Phys., 19 (2016), pp. 1141-1166. · Zbl 1388.65115
[8] Bao, W., Jiang, S., Tang, Q., and Zhang, Y., Computing the ground state and dynamics of the nonlinear Schrödinger equation with nonlocal interactions via the nonuniform FFT, J. Comput. Phys., 296 (2015), pp. 72-89. · Zbl 1354.65200
[9] Cai, Y., Rosenkranz, M., Lei, Z., and Bao, W., Mean-field regime of trapped dipolar Bose-Einstein condensates in one and two dimensions, Phys. Rev. A (3), 82 (2010), 043623.
[10] Exl, L., Mauser, N. J., and Zhang, Y., Accurate and efficient computation of nonlocal potentials based on Gaussian-sum approximation, J. Comput. Phys., 327 (2016), pp. 629-642. · Zbl 1422.65450
[11] Greengard, L., Jiang, S., and Zhang, Y., The anisotropic truncated kernel method for convolution with free-space Green’s functions, SIAM J. Sci. Comput., 40 (2018), pp. A3733-A3754. · Zbl 1481.65262
[12] Jarvis, M. R., White, I. D., Godby, R. W., and Payne, M. C., Supercell technique for total-energy calculations of finite charged and polar systems, Phys. Rev. B (3), 56 (1997), pp. 14972-14978.
[13] Jiang, S., Greengard, L., and Bao, W., Fast and accurate evaluation of nonlocal Coulomb and dipole-dipole interactions via the nonuniform FFT, SIAM J. Sci. Comput., 36 (2014), pp. B777-B794. · Zbl 1307.65184
[14] Liu, X., Tang, Q., and Zhang, Y., Numerical Study of the Quasi-2d Dipolar Bose-Einstein Condensates, manuscript.
[15] Mauser, N. J. and Zhang, Y., Exact artificial boundary condition for the Poisson equation in the simulation of the 2D Schrödinger-Poisson system, Commun. Comput. Phys., 16 (2014), pp. 764-780. · Zbl 1388.65063
[16] Ronen, S., Bortolotti, D. C., and Bohn, J., Bogoliubov modes of a dipolar condensate in a cylindrical trap, Phys. Rev. A (3), 74 (2006), 012623.
[17] Rozzi, C. A., Varsano, D., Marini, A., Gross, E. K. U., and Rubio, A., Exact Coulomb cutoff technique for supercell calculations, Phys. Rev. B (3), 73 (2006), 205119.
[18] Rudin, W., Functional Analysis, 2nd ed., McGraw-Hill, New York, 1991. · Zbl 0867.46001
[19] Shen, J., Tang, T., and Wang, L., Spectral Methods: Algorithms, Analysis and Applications, Springer, Berlin, 2011. · Zbl 1227.65117
[20] Trefethen, L. N., Spectral Methods in MATLAB, , SIAM, Philadelphia, 2000. · Zbl 0953.68643
[21] Vico, F., Greengard, L., and Ferrando, M., Fast convolution with free-space Green’s functions, J. Comput. Phys., 323 (2016), pp. 191-203. · Zbl 1415.65269
[22] Wang, A. B. and Yi, S., Trapped Bose-Einstein condensates with quadrupole-quadrupole interaction, Chinese Phys. B, 27 (2018), 120307.
[23] Yi, S. and You, L., Trapped atomic condensates with anisotropic interactions, Phys. Rev. A (3), 61 (2000), 041604.
[24] Yi, S. and You, L., Trapped condensates of atoms with dipole interactions, Phys. Rev. A (3), 63 (2001), 053607.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.