© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1903.08672 · Journal of Physics Condensed Matter

Vortex-lattice in a uniform Bose–Einstein condensate in a box trap

S K Adhikari

Published Mar 29, 2019• 13 citations

Abstract

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.

DOI: 10.1088/1361-648x/ab14c5

Related papers

Vortex dynamics in cubic-quintic Bose-Einstein condensates

Thudiyangal Mithun, Kuppuswamy Porsezian, Bishwajyoti Dey

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.

Physical Review E•Jul 8, 2013•25 citations

Three-dimensional vortex structures in a rotating dipolar Bose–Einstein condensate

Ramavarmaraja Kishor Kumar, Thangarasu Sriraman, Henrique Fabrelli +2 more

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.

arXiv:1506.08184•Journal of Physics B Atomic Molecular and Optical Physics•Jul 5, 2016•21 citations

Numerical studies on vortices in rotating dipolar Bose-Einstein condensates

R. Kishor Kumar, Paulsamy Muruganandam

We study the formation and dynamics of vortices in quasi two-dimensional dipolar Bose-Einstein condensates (BECs). We solve the time dependent two-dimensional Gross-Pitaevskii equation by using a combined split-step Crank-Nicolson (SSCN) and Fast Fourier-Transform (FFT) based numerical scheme and investigate the vortex dynamics in dipolar BECs trapped with harmonic and optical lattice potentials. The consequence of dipole-dipole interaction on vortex nucleation and number of vortices has been analysed. We observe that the breaking of symmetry due to anisotropic dipolar interaction enormously increase the speed of creation of vortices.

Journal of Physics Conference Series•Apr 9, 2014•3 citations

Ground States of a Rotating Spin-1 Bose-Einstein Condensate in Box Trap

Qing-Li Zhu, Ji-Mei Shen

Abstract The ground-state phases of rotating spin-1 Bose-Einstein condensates in three different box traps are studied numerically. In circular or square box trap, influenced by the combined effects of spin-exchange interaction and rotation, traditional skyrmion vortex for a rotating spin-1 BEC in harmonic trap are replaced by the homogeneous coreless vortices layered with fourfold rotational symmetry, which are usually generated from the trap boundary every four. In a triangular box trap, one skyrmion vortex structure is firstly found in the trap center, while other vortices are all generated from the boundary and arrange parallel to the trap edge in the form of vortex chains .

Journal of Low Temperature Physics•Jul 12, 2026•0 citations

Vortex lattice of a Bose-Einstein condensate in a rotating anisotropic trap

Mehmet Özgür Oktel

We study the vortex lattices in a Bose-Einstein condensate in a rotating anisotropic harmonic trap. We first investigate the single-particle wave functions obtained by the exact solution of the problem and give simple expressions for these wave functions in the small anisotropy limit. Depending on the strength of the interactions, a few or a large number of vortices can be formed. In the limit of many vortices, we calculate the density profile of the cloud and show that the vortex lattice stays triangular. We also find that the vortex lattice planes align themselves with the weak axis of the external potential. For a small number of vortices, we numerically solve the Gross-Pitaevskii equation and find vortex configurations that are very different from the vortex configurations in an axisymmetric rotating trap.

arXiv:cond-mat/0309613•Physical Review A•Feb 26, 2004•52 citations

Multi-vortex crystal lattices in Bose–Einstein condensates with a rotating trap

Shuangquan Xie, P. G. Kevrekidis, Théodore Kolokolnikov

We consider vortex dynamics in the context of Bose–Einstein condensates (BECs) with a rotating trap, with or without anisotropy. Starting with the Gross–Pitaevskii (GP) partial differential equation (PDE), we derive a novel reduced system of ordinary differential equations (ODEs) that describes stable configurations of multiple co-rotating vortices (vortex crystals). This description is found to be quite accurate quantitatively especially in the case of multiple vortices. In the limit of many vortices, BECs are known to form vortex crystal structures, whereby vortices tend to arrange themselves in a hexagonal-like spatial configuration. Using our asymptotic reduction, we derive the effective vortex crystal density and its radius. We also obtain an asymptotic estimate for the maximum number of vortices as a function of rotation rate. We extend considerations to the anisotropic trap case, confirming that a pair of vortices lying on the long (short) axis is linearly stable (unstable), corroborating the ODE reduction results with full PDE simulations. We then further investigate the many-vortex limit in the case of strong anisotropic potential. In this limit, the vortices tend to align themselves along the long axis, and we compute the effective one-dimensional vortex density, as well as the maximum admissible number of vortices. Detailed numerical simulations of the GP equation are used to confirm our analytical predictions.

arXiv:1708.04336•Proceedings of the Royal Society A Mathematical Physical and Engineering Sciences•May 1, 2018•5 citations

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.

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.

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.

We study the formation and dynamics of vortices in quasi two-dimensional dipolar Bose-Einstein condensates (BECs). We solve the time dependent two-dimensional Gross-Pitaevskii equation by using a combined split-step Crank-Nicolson (SSCN) and Fast Fourier-Transform (FFT) based numerical scheme and investigate the vortex dynamics in dipolar BECs trapped with harmonic and optical lattice potentials. The consequence of dipole-dipole interaction on vortex nucleation and number of vortices has been analysed. We observe that the breaking of symmetry due to anisotropic dipolar interaction enormously increase the speed of creation of vortices.

Abstract The ground-state phases of rotating spin-1 Bose-Einstein condensates in three different box traps are studied numerically. In circular or square box trap, influenced by the combined effects of spin-exchange interaction and rotation, traditional skyrmion vortex for a rotating spin-1 BEC in harmonic trap are replaced by the homogeneous coreless vortices layered with fourfold rotational symmetry, which are usually generated from the trap boundary every four. In a triangular box trap, one skyrmion vortex structure is firstly found in the trap center, while other vortices are all generated from the boundary and arrange parallel to the trap edge in the form of vortex chains .

We study the vortex lattices in a Bose-Einstein condensate in a rotating anisotropic harmonic trap. We first investigate the single-particle wave functions obtained by the exact solution of the problem and give simple expressions for these wave functions in the small anisotropy limit. Depending on the strength of the interactions, a few or a large number of vortices can be formed. In the limit of many vortices, we calculate the density profile of the cloud and show that the vortex lattice stays triangular. We also find that the vortex lattice planes align themselves with the weak axis of the external potential. For a small number of vortices, we numerically solve the Gross-Pitaevskii equation and find vortex configurations that are very different from the vortex configurations in an axisymmetric rotating trap.

We consider vortex dynamics in the context of Bose–Einstein condensates (BECs) with a rotating trap, with or without anisotropy. Starting with the Gross–Pitaevskii (GP) partial differential equation (PDE), we derive a novel reduced system of ordinary differential equations (ODEs) that describes stable configurations of multiple co-rotating vortices (vortex crystals). This description is found to be quite accurate quantitatively especially in the case of multiple vortices. In the limit of many vortices, BECs are known to form vortex crystal structures, whereby vortices tend to arrange themselves in a hexagonal-like spatial configuration. Using our asymptotic reduction, we derive the effective vortex crystal density and its radius. We also obtain an asymptotic estimate for the maximum number of vortices as a function of rotation rate. We extend considerations to the anisotropic trap case, confirming that a pair of vortices lying on the long (short) axis is linearly stable (unstable), corroborating the ODE reduction results with full PDE simulations. We then further investigate the many-vortex limit in the case of strong anisotropic potential. In this limit, the vortices tend to align themselves along the long axis, and we compute the effective one-dimensional vortex density, as well as the maximum admissible number of vortices. Detailed numerical simulations of the GP equation are used to confirm our analytical predictions.