A conservative sine pseudo-spectral-difference method for multi-dimensional coupled Gross–Pitaevskii equations

In this paper, a sine pseudo-spectral-difference scheme that preserves the discrete mass and energy is presented and analyzed for the coupled Gross–Pitaevskii equations with Dirichlet boundary conditions in several spatial dimensions. The Crank–Nicolson finite difference method is employed for approximating the time derivative, and the second-order sine spectral differentiation matrix is deduced and applied in spatial discretization. Without any restrictions on the grid ratios, optimal error estimates are established by utilizing the discrete energy method and the equivalence of (semi-)norms. An accelerated algorithm is developed to speed up the numerical implementation with the help of fast sine transform. Numerical examples are tested to confirm the effectiveness and high accuracy of the method.

[1]  Weizhu Bao,et al.  Ground States of Two-component Bose-Einstein Condensates with an Internal Atomic Josephson Junction , 2011 .

[2]  Tingchun Wang,et al.  Optimal l∞ error estimates of finite difference methods for the coupled Gross-Pitaevskii equations in high dimensions , 2014 .

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

[4]  C. E. Wieman,et al.  Vortices in a Bose Einstein condensate , 1999, QELS 2000.

[5]  Yushun Wang,et al.  Multi-symplectic Fourier pseudospectral method for the Kawahara equation , 2014 .

[6]  Yanzhi Zhang,et al.  An efficient spectral method for computing dynamics of rotating two-component Bose-Einstein condensates via coordinate transformation , 2013, J. Comput. Phys..

[7]  Qi Wang,et al.  A conservative Fourier pseudo-spectral method for the nonlinear Schrödinger equation , 2017, J. Comput. Phys..

[8]  Zhi-Zhong Sun,et al.  On the L∞ convergence of a difference scheme for coupled nonlinear Schrödinger equations , 2010, Comput. Math. Appl..

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

[10]  I. J. Anderson,et al.  Parallel factorization of banded linear matrices using a systolic array processor , 1996, Adv. Comput. Math..

[11]  Xiaofei Zhao.  On error estimates of an exponential wave integrator sine pseudospectral method for the Klein–Gordon–Zakharov system , 2016 .

[12]  Tingchun Wang,et al.  A linearized, decoupled, and energy‐preserving compact finite difference scheme for the coupled nonlinear Schrödinger equations , 2017 .

[13]  Dalibard,et al.  Vortex formation in a stirred bose-einstein condensate , 1999, Physical review letters.

[14]  W. Ketterle,et al.  Bose-Einstein condensation , 1997 .

[15]  Anatoly Svidzinsky,et al.  Vortices in a trapped dilute Bose-Einstein condensate , 2001, cond-mat/0102003.

[16]  Tingchun Wang,et al.  Optimal Point-Wise Error Estimate of a Compact Difference Scheme for the Coupled Gross–Pitaevskii Equations in One Dimension , 2014, J. Sci. Comput..

[17]  Yanzhi Zhang,et al.  Dynamics of rotating two-component Bose-Einstein condensates and its efficient computation , 2007 .

[18]  J. R. Ensher,et al.  Dynamics of component separation in a binary mixture of Bose-Einstein condensates , 1998 .

[19]  Zhi‐zhong Sun,et al.  Convergence analysis of a linearized Crank–Nicolson scheme for the two‐dimensional complex Ginzburg–Landau equation , 2013 .

[20]  Christo I. Christov,et al.  Strong coupling of Schrödinger equations: Conservative scheme approach , 2005, Math. Comput. Simul..

[21]  Xiaofei Zhao,et al.  Unconditional and optimal H2-error estimates of two linear and conservative finite difference schemes for the Klein-Gordon-Schrödinger equation in high dimensions , 2018, Adv. Comput. Math..

[22]  A. Quarteroni,et al.  Approximation results for orthogonal polynomials in Sobolev spaces , 1982 .

[23]  Hanquan Wang,et al.  A time-splitting spectral method for coupled Gross-Pitaevskii equations with applications to rotating Bose-Einstein condensates , 2007 .

[24]  Weizhu Bao,et al.  An Exponential Wave Integrator Sine Pseudospectral Method for the Klein-Gordon-Zakharov System , 2013, SIAM J. Sci. Comput..

[25]  W. Ketterle,et al.  Observation of Vortex Lattices in Bose-Einstein Condensates , 2001, Science.

[26]  Weizhu Bao,et al.  Optimal error estimates of finite difference methods for the Gross-Pitaevskii equation with angular momentum rotation , 2012, Math. Comput..

[27]  Jie Shen,et al.  A Fourth-Order Time-Splitting Laguerre-Hermite Pseudospectral Method for Bose-Einstein Condensates , 2005, SIAM J. Sci. Comput..

[28]  Tingchun Wang,et al.  Unconditional and optimal H1 error estimate of a Crank–Nicolson finite difference scheme for the Gross–Pitaevskii equation with an angular momentum rotation term , 2018 .

[29]  Xuanchun Dong,et al.  Stability and convergence of trigonometric integrator pseudospectral discretization for N-coupled nonlinear Klein-Gordon equations , 2014, Appl. Math. Comput..

[30]  Luming Zhang,et al.  Error estimates of a trigonometric integrator sine pseudo-spectral method for the extended Fisher–Kolmogorov equation , 2018, Applied Numerical Mathematics.

[31]  Tingchun Wang,et al.  Fourth-order compact and energy conservative difference schemes for the nonlinear Schrödinger equation in two dimensions , 2013, J. Comput. Phys..

[32]  M. S. Ismail.  Numerical solution of coupled nonlinear Schrödinger equation by Galerkin method , 2008, Math. Comput. Simul..