For a Bose-Einstein Condensate placed in a rotating trap and confined in the z axis, a kind of symplectic difference schemes was used to investigate the evolution of vortices of BEC in this paper. The case of the evolution of three vortices in BEC was imitated by numerical experiment.
We address the challenging proposition of using real experimental parameters in a three-dimensional (3D) numerical simulation of fast rotating Bose-Einstein condensates. We simulate recent experiments [V. Bretin, S. Stock, Y. Seurin, and J. Dalibard, Phys. Rev. Lett. 92, 050403 (2004); S. Stock, V. Bretin, S. Stock, F. Chevy, and J. Dalibard, Europhys. Lett. 65, 594 (2004)] using an anharmonic (quadratic-plus-quartic) confining potential to reach rotation frequencies $(\ensuremath{\Omega})$ above the trap frequency $({\ensuremath{\omega}}_{\ensuremath{\perp}})$. Our numerical results are obtained by propagating the 3D Gross-Pitaevskii equation in imaginary time. For $\ensuremath{\Omega}\ensuremath{\leqslant}{\ensuremath{\omega}}_{\ensuremath{\perp}}$, we obtain an equilibrium vortex lattice similar (as the size and number of vortices) to experimental observations. For $\ensuremath{\Omega}>{\ensuremath{\omega}}_{\ensuremath{\perp}}$ we observe the evolution of the vortex lattice into an array of vortices with a central hole. Since this evolution was not visible in experiments, we investigate the 3D structure of vortex configurations and 3D effects on vortex contrast. Numerical data are also compared to recent theory [D. E. Sheehy and L. Radzihovsky, Phys. Rev. A 70, 063620 (2004)] describing vortex lattice inhomogeneities and a remarkably good agreement is found.
We consider a rotating Bose-Einstein condensate in a harmonic trap and investigate numerically the behavior of the wave function which solves the Gross-Pitaevskii equation. Following recent experiments [P. Rosenbuch, V. Bretin, and J. Dalibard, Phys. Rev. Lett. 89, 200403 (2002)], we study in detail the line of a single quantized vortex, which has a U or S shape. We find that a single vortex can lie only in the $x\ensuremath{-}z$ or $y\ensuremath{-}z$ plane. S-type vortices exist for all values of the angular velocity $\ensuremath{\Omega}$ while U vortices exist for $\ensuremath{\Omega}$ sufficiently large. We compute the energy of the various configurations with several vortices and study the three-dimensional structure of vortices.
Abstract We study three-dimensional vortex lattice structures in purely dipolar Bose–Einstein condensate (BEC). By using the mean-field approximation, we obtain a stability diagram for the vortex states in purely dipolar BECs as a function of harmonic trap aspect ratio ( λ ) and dipole–dipole interaction strength ( D ) under rotation. Rotating the condensate within the unstable region leads to collapse while in the stable region furnishes stable vortex lattices of dipolar BECs. We analyse stable vortex lattice structures by solving the three-dimensional time-dependent Gross–Pitaevskii equation in imaginary time. Further, the stability of vortex states is examined by evolution in real-time. We also investigate the distribution of vortices in a fully anisotropic trap by increasing eccentricity of the external trapping potential. We observe the breaking up of the condensate in two parts with an equal number of vortices on each when the trap is sufficiently weak, and the rotation frequency is high.