Abstract Recently, we presented a new numerical scheme for vortex lattice formation in a rotating Bose–Einstein condensate (BEC) using smoothed particle hydrodynamics (SPH) with an explicit time-integrating scheme; our SPH scheme could reproduce the vortex lattices and their formation processes in rotating quasi-two-dimensional (2D) BECs trapped in a 2D harmonic potential. In this study, we have successfully demonstrated a simulation of rotating 3D BECs trapped in a 3D harmonic potential forming ‘cigar-shaped’ condensates. We have found that our scheme can reproduce the following typical behaviors of rotating 3D BECs observed in the literature: (i) the characteristic shape of the lattice formed in the cross-section at the origin and its formation process, (ii) the stable existence of vortex lines along the vertical axis after reaching the steady state, (iii) a ‘cookie-cutter’ shape, with a similar lattice shape observed wherever we cut the condensate in a certain range in the vertical direction, (iv) the bending of vortex lines when approaching the inner edges of the condensate, and (v) the formation of vortex lattices by vortices entering from outside the condensate. Therefore, we further validated our scheme by simulating rotating 3D BECs.
Abstract Recently, we presented a new numerical scheme for vortex lattice formation in a rotating Bose–Einstein condensate (BEC) using smoothed particle hydrodynamics (SPH) with an explicit time-integrating scheme; our SPH scheme could reproduce the vortex lattices and their formation processes in rotating quasi-two-dimensional (2D) BECs trapped in a 2D harmonic potential. In this study, we have successfully demonstrated a simulation of rotating 3D BECs trapped in a 3D harmonic potential forming ‘cigar-shaped’ condensates. We have found that our scheme can reproduce the following typical behaviors of rotating 3D BECs observed in the literature: (i) the characteristic shape of the lattice formed in the cross-section at the origin and its formation process, (ii) the stable existence of vortex lines along the vertical axis after reaching the steady state, (iii) a ‘cookie-cutter’ shape, with a similar lattice shape observed wherever we cut the condensate in a certain range in the vertical direction, (iv) the bending of vortex lines when approaching the inner edges of the condensate, and (v) the formation of vortex lattices by vortices entering from outside the condensate. Therefore, we further validated our scheme by simulating rotating 3D BECs.
This study proposed a new numerical scheme for vortex lattice formation in a rotating Bose–Einstein condensate (BEC) using smoothed particle hydrodynamics (SPH) with an explicit real-time integration scheme. Specifically, the Gross–Pitaevskii equation was described as a complex representation to obtain a pair of time-dependent equations, which were then solved simultaneously following discretization based on SPH particle approximation. We adopt the fourth-order Runge–Kutta method for time evolution. We performed simulations of a rotating Bose gas trapped in a harmonic potential, showing results that qualitatively agreed with previously reported experiments and simulations. The geometric patterns of formed lattices were successfully reproduced for several cases, for example, the hexagonal lattice observed in the experiments of rotating BECs. Consequently, it was confirmed that the simulation began with the periodic oscillation of the condensate, which attenuated and maintained a stable rotation with slanted elliptical shapes; however, the surface was excited to be unstable and generated ripples, which grew into vortices and then penetrated inside the condensate, forming a lattice. We confirmed that each branch point of the phase of wavefunctions corresponds to each vortex. These results demonstrate our approach at a certain degree of accuracy. In conclusion, we successfully developed a new SPH scheme for the simulations of vortex lattice formation in rotating BECs.
We study vortex dynamics in a trapped Bose-Einstein condensate with tunable two- and three-body interactions. The dynamics is governed by two-dimensional cubic-quintic Gross-Pitaevskii equation. A time-dependent variational method has been used to obtain critical rotational frequency and surface mode frequency analytically and are compared with numerical simulation results. An imaginary time propagation method and Crank-Nicolson scheme for discretization have been used for numerical simulation. The numerically calculated average value of the angular momentum per particle shows very clearly its dynamical relation with the time development of the vortex formation. The rotational frequency dependence of the variation of average value of the angular momentum per particle with time shows that vortices form much faster in time for higher rotational frequency. Similarly, the vortex forms much faster in time with an increase of the strength of the repulsive three-body interaction. The simulation of the vortex lattice formation in the condensate shows that the presence of the three-body interactions does not alter the vortex lattice pattern but it helps in the shape deformations of the condensate thereby leading to vortex lattice formation. Likewise, the three-body interactions enable the vortex lattice formation in Bose-Einstein condensates even with attractive two-body interactions and in purely quintic BEC.
We numerically study the dynamics of vortex lattice formation in a rotating cigar-shaped Bose-Einstein condensate. The study is a three-dimensional simulation of the Gross-Pitaevskii equation with a phenomenological dissipation term. The simulations reveal previously unknown dynamical features of the vortex nucleation process, in which the condensate undergoes a strongly turbulent stage and the penetrating vortex lines vibrate rapidly. The vibrations arise from spontaneous excitation of Kelvin waves on the protovortices during a surface wave instability, caused by an inhomogeneity of the condensate density along the elongated axial direction.
We study numerically the vortex-lattice formation in a rapidly rotating uniform quasi-two-dimensional Bose-Einstein condensate (BEC) in a box trap. We consider two types of boxes: square and circle. In a square-shaped 2D box trap, when the number of generated vortices is the square of an integer, the vortices are found to be arranged in a perfect square lattice, although deviations near the center are found when the number of generated vortices is arbitrary. In case of a circular box trap, the generated vortices in the rapidly rotating BEC lie on concentric closed orbits. Near the center, these orbits have the shape of polygons, whereas near the periphery the orbits are circles. The circular box trap is equivalent to the rotating cylindrical bucket used in early experiment(s) with liquid He II. The number of generated vortices in both cases is in qualitative agreement with Feynman's universal estimate. The numerical simulation for this study is performed by a solution of the underlying mean-field Gross-Pitaevskii (GP) equation in the rotating frame, where the wave function for the generated vortex lattice is a stationary state. Consequently, the imaginary-time propagation method can be used for a solution of the GP equation, known to lead to an accurate numerical solution. We also demonstrated the dynamical stability of the vortex lattices in real-time propagation upon a small change of the angular frequency of rotation, using the converged imaginary-time wave function as the initial state.
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.
Bose-Einstein condensate (BEC) is essentially a macroscopic quantum effect with quantum volatility, macroscopic quantum coherence and artificial controllability. Owing to its unique controllability, it becomes a new ideal platform for quantum simulations and studies of interacting quantum systems. In this paper, the generation of vortices and the formation of vortex chains, as well as characteristics of vortex chains in rotating two-dimensional BEC in a potential composed of harmonic potential and optical lattice are studied numerically. Firstly, the generation of vortices, the formation and distribution of vortex chains and the effects of different physical parameters on the vortex chains in two-dimensional BEC are investigated by using the multigrid preconditioned conjugate gradient method. Secondly, the evolution of the vortex chains with time is studied by using the time-splitting spectral method. The results show that the generation of vortices in BEC trapped in the compound potential corresponds to the minimum value of the potential. When the depth of the optical lattice increases to a certain value, vortex chains are formed in the BEC. With the further increase of the depth of the optical lattice, the vortex depth in the vortex chain in the BEC decreases continuously, and finally the vortex chain disappears completely. When the interaction strength between atoms increases, the distribution range of the condensate expands, and the number of vortices and the number of vortex chains in the condensate also increase. When the interaction strength between atoms increases to a certain value, the symmetry of the vortex chains is broken. As the rotation frequency of the condensate increases, the distribution range of the condensate expands, and the number of vortices and the number of vortex chains in the condensate also increase. When the rotation frequency is close to the external trapping potential frequency, the linear alignment of the vortex chains is disrupted. It is also found that there are three stages in the evolution of the vortex chains in the BEC: in the first stage, vortex chains rotate together with the condensate, and the original chain distribution keeps unchanged; in the second stage, the phenomenon of vortex space extrusion appears, and the vortex chain is destroyed; in the third stage, the phenomenon of vortex space expansion occurs, and finally the vortex chains disappear. The results above show that the depth of the optical lattice, the interaction strength between atoms, and the rotation frequency of the condensate have important effects on the vortices and vortex chains in the condensate. By adjusting these physical quantities, the number of vortices and the shape of vortex chains in the BEC can be effectively manipulated. This may provide some theoretical reference and guidance for future experiments and applications.