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.

The authors would like to thank the anonymous referees for their careful reading and many constructive comments and suggestions of the manuscript.

X Li is supported by a grant no. KJ2018A0523 from the University Natural Science Research Key Project of Anhui Province.

Authors and Affiliations

1. Department of Mathematics, Anhui Science and Technology University, Fengyang, 233100, Anhui, People’s Republic of China

   Xin Li

2. Department of Mathematics, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, Jiangsu, People’s Republic of China

   Xin Li & Luming Zhang

Corresponding author

Correspondence to Luming Zhang.

Communicated by: Jan Hesthaven

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendix A: Proof of Lemma 2.1

To obtain the second-order SSDM, the following lemma needs to be verified at first.

Lemma A.1 1

Denote 𝜃_x := μ_x(x − x_j), 𝜃_y := μ_y(y − y_k), \(\theta _{x}^{\prime }:=\mu _{x}(x+x_{j}-2a)\), \(\theta _{y}^{\prime }:=\mu _{y}(y+y_{k}-2c)\). The interpolation basis functions can be reformulated as:

with \(\mathfrak {S}(N,\mu ,\theta ):=\frac {1}{2N}\big [\sin \limits (N\theta )\cot \frac {\theta }{2}-\cos \limits (N\theta )\big ]\).

Proof

We just prove the assertion of X_j(x). The other one can be deduced correspondingly. Applying the prosthaphaeresis formula and the definitions of 𝜃_x and \(\theta _{x}^{\prime }\), (??) implies that

Invoking to the identity

we obtain the desired results and complete the proof of Lemma A.1. □

Proof Proof of Lemma 2.1

For simplicity, we also verify the assertion in the x-direction and the assertion of y-direction is equally acceptable. By means of the above lemma, the second-order derivative of X_j(x) can be given by:

where the second-order derivative \(\mathfrak {\ddot {S}}(N,\mu ,\theta )\) is obtained directly by:

To obtain the second-order SSDM, we divide the proof into two cases:

1. (i)

   If x = x_p≠x_j, the definitions of 𝜃_x and \(\theta _{x}^{\prime }\) lead to:

   $$ \begin{array}{@{}rcl@{}} \sin(N_{x}\theta_{x})=\sin(N_{x}\theta_{x}^{\prime})=0,\ \ \cos(N_{x}\theta_{x})=\cos(N_{x}\theta_{x}^{\prime})=(-1)^{j+p}. \end{array} $$

   Substituting the results into \(\mathrm {X}_{j}^{\prime \prime }(x)\) and using the definition of μ_x, for j≠p, one has:

   $$ \begin{array}{@{}rcl@{}} \mathrm{X}_{j}^{\prime\prime}(x_{p})&=&(-1)^{j+p+1}\frac{{\mu_{x}^{2}}}{2}(\csc^{2}\frac{\theta_{x}}{2}-\csc^{2}\frac{\theta_{x}^{\prime}}{2})\\ &=&(-1)^{j+p+1}\frac{{\mu_{x}^{2}}}{2}\bigg[\csc^{2}\big(\frac{\mu_{x}}{2}(j-p)h_{x}\big)- \csc^{2}\big(\frac{\mu_{x}}{2}(j+p)h_{x}\big)\bigg]. \end{array} $$

   This completes the proof of the first case of (??).

2. (ii)

   If x = x_p = x_j, the first term of \(\mathrm {X}_{j}^{\prime \prime }(x)\), i.e. \(\mathfrak {\ddot {S}}(N_{x},\mu _{x},\theta _{x})\) equals to the following:

   $$ \begin{array}{@{}rcl@{}} &\frac{\frac{{\mu_{x}^{2}}}{4N_{x}}\sin(N_{x}\theta_{x})\cos\frac{\theta_{x}}{2}-\frac{N_{x}{\mu_{x}^{2}}}{2}\sin(N_{x}\theta_{x})\cos\frac{\theta_{x}}{2}\sin^{2} \frac{\theta_{x}}{2}-\frac{{\mu_{x}^{2}}}{2}\cos(N_{x}\theta_{x})\sin\frac{\theta_{x}}{2} }{\sin^{3}\frac{\theta_{x}}{2}}+\frac{N_{x}{\mu_{x}^{2}}}{2}\cos(N_{x}\theta_{x}). \end{array} $$

   We apply Taylor’s expansion for each term of the numerator and find:

   $$ \begin{array}{@{}rcl@{}} &\sin(N_{x}\theta_{x})\cos\frac{\theta_{x}}{2}=\Big(N_{x}\theta_{x} -\frac{(N_{x}\theta_{x})^{3}}{3!}+\cdots\Big)\Big(1-\frac{(\frac{\theta_{x}}{2})^{2}}{2!}+\cdots\Big)\\ &\ \ \ \ \ \ \ \ \ \ \ \quad \ \ \ \ \ \ \ \ =N_{x}\theta_{x}-N_{x}\big(\frac{\theta_{x}}{2}\big)^{3}-\frac{4{N_{x}^{3}}}{3}\big(\frac{\theta_{x}}{2}\big)^{3}+O\big({\theta_{x}^{5}}\big),\\ &\sin(N_{x}\theta_{x})\cos\frac{\theta_{x}}{2}\sin^{2}\frac{\theta_{x}}{2}=\Big(N_{x}\theta_{x} -\frac{(N_{x}\theta_{x})^{3}}{3!}+\cdots\Big)\Big(1-\frac{(\frac{\theta_{x}}{2})^{2}}{2!}+\cdots\Big)\Big(\frac{\theta_{x}}{2}-\frac{(\frac{\theta_{x}}{2})^{3}}{3!}+\cdots\Big)^{2}\\ &\ \ \ \ \ \ \ \ \ \ \ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ =2N_{x}\big(\frac{\theta_{x}}{2}\big)^{3}+O\big({\theta_{x}^{5}}\big),\\ &\cos(N_{x}\theta_{x})\sin\frac{\theta_{x}}{2}=\Big(1-\frac{(N_{x}\theta_{x})^{2}}{2!}+\cdots\Big)\Big(\frac{\theta_{x}}{2}-\frac{(\frac{\theta_{x}}{2})^{3}}{3!}+\cdots\Big)\\ &\ \ \ \ \ \ \ \ \ \ \ \quad \ \ \ \ \ \ \ \ =\frac{\theta_{x}}{2}-\frac{1}{6}\big(\frac{\theta_{x}}{2}\big)^{3}-2{N_{x}^{2}}\big(\frac{\theta_{x}}{2}\big)^{3}+O\big({\theta_{x}^{5}}\big). \end{array} $$

   Since \(\theta _{x} \rightarrow 0\) as \(x_{p} \rightarrow x_{j}\). We combine these four equalities and obtain:

   $$ \begin{array}{@{}rcl@{}} \mathfrak{\ddot{S}}(N_{x},\mu_{x},\theta_{x}) \rightarrow -\frac{{\mu_{x}^{2}}}{6}+\frac{N_{x}{\mu_{x}^{2}}}{2}-\frac{{N_{x}^{2}}{\mu_{x}^{2}}}{3},\quad\text{as} x_{p} \rightarrow x_{j}. \end{array} $$

   For the second term of \(\mathrm {X}_{j}^{\prime \prime }(x)\), i.e., \(\mathfrak {\ddot {S}}(N_{x},\mu _{x},\theta _{x}^{\prime })\), since \(\sin \limits \frac {\theta _{x}^{\prime }}{2}\neq 0\), it is not difficult to check:

   $$ \begin{array}{@{}rcl@{}} \mathfrak{\ddot{S}}(N_{x},\mu_{x},\theta_{x}^{\prime})=-\frac{{\mu_{x}^{2}}}{2}(-1)^{j+p}\csc^{2}\frac{\theta_{x}^{\prime}}{2}+\frac{N_{x}{\mu_{x}^{2}}}{2}(-1)^{j+p} =\frac{N_{x}{\mu_{x}^{2}}}{2}-\frac{{\mu_{x}^{2}}}{2}\csc^{2}(j\mu h_{x}). \end{array} $$

   Substituting \(\mathfrak {\ddot {S}}(N_{x},\mu _{x},\theta _{x})\) and \(\mathfrak {\ddot {S}}(N_{x},\mu _{x},\theta _{x}^{\prime })\) into \(\mathrm {X}_{j}^{\prime \prime }(x)\), for j = p, one has:

   $$ \begin{array}{@{}rcl@{}} \mathrm{X}_{j}^{\prime\prime}(x_{p})=-\frac{{\mu_{x}^{2}}}{6}-\frac{{N_{x}^{2}}{\mu_{x}^{2}}}{3}+\frac{{\mu_{x}^{2}}}{2}\csc^{2}(j\mu h_{x}). \end{array} $$

   This completes the proof of (??) and obtains the claimed results in Lemma 2.1.

□

Appendix B: Proof of Lemma 3.1

Always, we prove the assertion in the x-direction and get the other one identically. Differentiating (??) two times and taking x = x_p yield:

Denote e_l := μ_xl, the spatial discretization gives:

Let

it is easy to find:

and

This completes the proof of \({S_{2}^{x}}\). For A₁, multiplying \(-{h_{x}^{2}}\) on both sides of the decomposition formula in Lemma 3.1 and following the same procedure, one has:

Applying the prosthaphaeresis formula on the right-hand side of the equality, it follows that

where

Substituting 𝜗₁ into 𝜗₂ and utilizing the prosthaphaeresis formula again, the equality mentioned above can be reformulated as follows:

where

All these three items on the right-hand side have the same structures. Resorting to the identity:

the following observations can be easily obtained:

Substituting these results into the equality, we obtain the assertion of A₁ directly. This completes the proof of Lemma 3.1.

Reprints and permissions