We numerically study the nonlocal Gross–Pitaevskii equation (NGPE) which describes the dynamics of Bose–Einstein condensates (BEC) with dipole–dipole interaction at extremely low temperature. In preparation for the numerics, first we reformulate the dimensionless NGPE into a Schrödinger–Poisson system. Then, we discretize the three-dimensional Schrödinger–Poisson system in space by a sixth-order compact finite difference method and in time by a splitting technique. By means of three-dimensional discrete fast Sine transform, we develop a fast solver for the resulting discretized system. Finally, we present numerical examples in three dimensions to demonstrate the power of the numerical methods and to discuss some physics of dipolar BEC. The merits of the proposed method for the NGPE are that it is fast and unconditionally stable. Moreover, the method is of spectral-like accuracy in space, and conserves the particle number and the energy of the system in the discretized level.
We numerically study the nonlocal Gross–Pitaevskii equation (NGPE) which describes the dynamics of Bose–Einstein condensates (BEC) with dipole–dipole interaction at extremely low temperature. In preparation for the numerics, first we reformulate the dimensionless NGPE into a Schrödinger–Poisson system. Then, we discretize the three-dimensional Schrödinger–Poisson system in space by a sixth-order compact finite difference method and in time by a splitting technique. By means of three-dimensional discrete fast Sine transform, we develop a fast solver for the resulting discretized system. Finally, we present numerical examples in three dimensions to demonstrate the power of the numerical methods and to discuss some physics of dipolar BEC. The merits of the proposed method for the NGPE are that it is fast and unconditionally stable. Moreover, the method is of spectral-like accuracy in space, and conserves the particle number and the energy of the system in the discretized level.
A two-component Bose-Einstein condensate (BEC) described by two coupled a three-dimension Gross-Pitaevskii (GP) equations is considered, where one equation has dipole-dipole interaction while the other one has only the usual s-wave contact interaction, in a cigar trap. The time-splitting and sine spectral method in space is proposed to discretize the time-dependent equations for computing the dynamics of dipolar BEC. The singularity in the dipole-dipole interaction brings significant difficulties both in mathematical analysis and in numerical simulations. Numerical results are given to show the efficiency of this method.
We propose a simple, efficient and accurate numerical method for simulating the dynamics of rotating Bose-Einstein condensates (BECs) in a rotational frame with/without a long-range dipole-dipole interaction. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range dipole-dipole interaction, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential and review some dynamical laws related to the 2D and 3D GPE. By introducing a rotating Lagrangian coordinate system, the original GPEs are re-formulated to GPEs without the angular momentum rotation which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable and very efficient in computation. It is spectral order accurate in space and second-order accurate in time, and conserves the mass in the discrete level. Extensive numerical results are reported to demonstrate the efficiency and accuracy of the new numerical method. Finally, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation and center-of-mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without/with the long-range dipole-dipole interaction.
We propose a simple, efficient, and accurate numerical method for simulating the dynamics of rotating Bose--Einstein condensates (BECs) in a rotational frame with or without long-range dipole-dipole interaction (DDI). We begin with the three-dimensional (3D) Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range DDI, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential, and review some dynamical laws related to the 2D and 3D GPEs. By introducing a rotating Lagrangian coordinate system, the original GPEs are reformulated to GPEs without the angular momentum rotation, which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable, and very efficient in computation. It is spectral-order accurate in space and second-order accurate in time and conserves the mass on the discrete level. We compare our method with some representative methods in the literature to demonstrate its efficiency and accuracy. In addition, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation, and center of mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without or with the long-range DDI.
We propose efficient and accurate numerical methods for computing the ground state and dynamics of the dipolar Bose-Einstein condensates utilising a newly developed dipole-dipole interaction (DDI) solver that is implemented with the non-uniform fast Fourier transform (NUFFT) algorithm. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with a DDI term and present the corresponding two-dimensional (2D) model under a strongly anisotropic confining potential. Different from existing methods, the NUFFT based DDI solver removes the singularity by adopting the spherical/polar coordinates in the Fourier space in 3D/2D, respectively, thus it can achieve spectral accuracy in space and simultaneously maintain high efficiency by making full use of FFT and NUFFT whenever it is necessary and/or needed. Then, we incorporate this solver into existing successful methods for computing the ground state and dynamics of GPE with a DDI for dipolar BEC. Extensive numerical comparisons with existing methods are carried out for computing the DDI, ground states and dynamics of the dipolar BEC. Numerical results show that our new methods outperform existing methods in terms of both accuracy and efficiency.
We propose a time-splitting spectral method for the generalized Gross–Pitaevskii equations, which describe the dynamics of spinor F=1 Bose–Einstein condensates at a very low temperature. The new numerical method is explicit, unconditionally stable, and of spectral accuracy in space. Moreover, it conserves the position densities at the discretized level. We apply the method for studying both the dynamic generation of vortices and vortex lattice dynamics for the spinor F=1 Bose–Einstein condensates held in an Ioffe-Pritchard magnetic field.
In this paper, we mainly review recent results on mathematical theory andnumerical methods for Bose-Einstein condensation (BEC),based on the Gross-Pitaevskii equation (GPE). Starting from the simplest case with one-component BEC of the weakly interacting bosons, we study the reduction of GPE to lower dimensions, the ground states of BEC including the existence and uniqueness as well as nonexistence results, and the dynamics of GPE including dynamical laws, well-posedness of the Cauchy problem as well as the finite time blow-up. To compute the ground state, the gradient flow with discrete normalization (or imaginary time) method is reviewed and various full discretization methods are presented and compared. To simulate the dynamics, both finite difference methods and time splitting spectral methods are reviewed, and their error estimates are briefly outlined. When the GPE has symmetric properties, we show how to simplify the numerical methods. Then we compare two widely used scalings, i.e. physical scaling (commonly used) and semiclassical scaling, for BEC in strong repulsive interaction regime (Thomas-Fermi regime), and discuss semiclassical limits of the GPE. Extensions of these results for one-component BEC are thencarried out for rotating BEC by GPE with an angular momentum rotation, dipolar BEC by GPE with long range dipole-dipole interaction, and two-component BEC by coupled GPEs. Finally, as a perspective, we show briefly the mathematical models forspin-1 BEC, Bogoliubov excitation and BEC at finite temperature.