We present an efficient numerical algorithm to observe the dynamical formation of vortex lattice of a rotating trapped Bose-Einstein condensate (BEC) via solving a two-dimensional Gross-Pitaevskii ...
We study efficient continuation methods for computing the ground state solution of quasi-2D rotating dipolar Bose-Einstein condensates (BECs). First, the highly accurate spectral collocation method is used to discretize the governing Gross-Pitaevskii equation (GPE). Then, we modify the two-level continuation scheme for 3D dipolar BECs described in Jeng et al. (2014) to develop a single-parameter continuation method for quasi-2D rotating dipolar BECs, where the chemical potential is treated as the continuation parameter. Further, by adding the ratio of dipolar interaction strength to contact interaction strength as the second continuation parameter, we propose an efficient two-parameter continuation method which can effectively show the change of the ground-state vortex structures as the dipolar interaction strength gradually increases. Moreover, we also study linear stability analysis for the GPE. Sample numerical results on quasi-2D rotating dipolar BECs are reported.
We study efficient spectral-collocation and continuation methods (SCCM) for rotating two-component Bose-Einstein condensates (BECs) and rotating two-component BECs in optical lattices, where the second kind Chebyshev polynomials are used as the basis functions for the trial function space. A novel two-parameter continuation algorithm is proposed for computing the ground state and first excited state solutions of the governing Gross-Pitaevskii equations (GPEs), where the classical tangent vector is split into two constraint conditions for the bordered linear systems. Numerical results on rotating two-component BECs and rotating two-component BECs in optical lattices are reported. The results on the former are consistent with the published numerical results.
We study the dynamics of rotating two-component Bose-Einstein condensates(BECs) toward vortex lattice formation by numerically solving the two-dimensional Gross-Pitaevskii equations.The dynamics are similar to that of single-component BEC, and the whole changes of physical quantities in two cases are qualitatively similar.However, there are some differences between two cases.First, the characteristic time toward vortex entering depends on the inter-component interaction.The time toward vortex entering becomes shorter with increasing the inter-component interaction.Second, the angular momenta exchange between two components with keeping the total angular momentum.
We study a radial basis function collocation method (RBFCM) to discretize a coupled nonlinear Schr o¨ dinger equation (CNLSE) that governs a two dimensional rotating Bose-Einstein condensate (BEC) with an angular momentum rotation term. We exploit a RBFCM-continuation method (RBFCM-CM) to trace the solution curve of the CNLSE. We compare the performance of the RBFCM-CM with the FEM-CM. We observe that the RBFCM-CM is very robust in a coarse grid for resolving the ground state solution with many vortices when the angular momentum rotation is close to the limit. Numerical results demonstrate the efficiency and accuracy of the RBFCM-CM for computing the superfluid density of the ground level of the BEC.
We describe multi-parameter continuation methods combined with spectral collocation methods for computing numerical solutions of rotating two-component Bose–Einstein condensates (BECs), which are governed by the Gross–Pitaevskii equations (GPEs). Various types of orthogonal polynomials are used as the basis functions for the trial function space. A novel multi-parameter/multiscale continuation algorithm is proposed for computing the solutions of the governing GPEs, where the chemical potential of each component and angular velocity are treated as the continuation parameters simultaneously. The proposed algorithm can effectively compute numerical solutions with abundant physical phenomena. Numerical results on rotating two-component BECs are reported.