Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Communications in Computational Physics

Numerical Study of Vortex Interactions in Bose-Einstein Condensation

Yanzhi Zhang

Published Mar 15, 2010• 11 citations
Publisher

Abstract

The dynamics and interaction of quantized vortices in Bose-Einstein condensates (BECs) are investigated by using the two-dimensional Gross-Pitaevskii equation (GPE) with/without an angular momentum rotation term. If all vortices have the same winding number, they would rotate around the trap center but never collide. In contrast, if the winding numbers are different, their interaction highly depends on the initial distance between vortex centers. The analytical results are presented to describe the dynamics of the vortex centers when β=0. While if β≠0, there is no analytical result but some conclusive numerical findings are provided for the further understanding of vortex interaction in BECs. Finally, the dynamic laws describing the relation of vortex interaction in nonrotating and rotating BECs are presented.

DOI
10.4208/cicp.2009.09.104

Related papers

Dynamics of Rotating Bose--Einstein Condensates and its Efficient and Accurate Numerical Computation

Weizhu Bao, Qiang Du, Yanzhi Zhang

In this paper, we study the dynamics of rotating Bose--Einstein condensates (BEC) based on the Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and present an efficient and accurate algorithm for numerical simulations. We examine the conservation of the angular momentum expectation and the condensate width and analyze the dynamics of a stationary state with a shift in its center. By formulating the equation in either the two-dimensional polar coordinate system or the three-dimensional cylindrical coordinate system, the angular momentum rotation term becomes a term with constant coefficients. This allows us to develop an efficient time-splitting method which is time reversible, unconditionally stable, efficient, and accurate for the problem. Moreover, it conserves the position density. We also apply the numerical method to study issues such as the stability of central vortex states and the quantized vortex lattice dynamics in rotating BEC.

SIAM Journal on Applied Mathematics•Jan 1, 2006•116 citations

Vortex stability in interacting Bose-Einstein condensates

Ajay Srinivasan, Aaron Wirthwein, Stephan Haas

We study the stability of vortices in a binary system of Bose-Einstein condensat es, with their wave functions modeled by a set of coupled, time-dependent Gross-Pitaevskii equations. Beginning with an effective two-dimensional system, we identify miscible and immiscible regimes characterized by the inter- and intra-atomic interactions and the initial configuration of the system. We then consider a binary system of Bose-Einstein condensates placed in a rotating harmonic trap and study the single united vortex state in this system. We derive an approximate form for the energy of a single united vortex in the binary system and the critical angular velocity for the global stability of a vortex at the center of the trap. We also compute the metastability onset angular velocity for the local stability of a vortex at the center of the trap. Numerical solutions to the Gross-Pitaevskii equations support these expressions. These rotational results inform us of a novel bimodal formation within the miscible regime of the binary condensate system. We thus demonstrate the nontrivial aspects of vortex stability in interacting binary Bose-Einstein condensates as a result of their nonlinear interactions.

arXiv:2501.06119•Physical review. B./Physical review. B•Sep 10, 2025•0 citations

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

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

Vortex Stability of Interacting Bose-Einstein Condensates Confined in Anisotropic Harmonic Traps

David L. Feder, Charles W. Clark, Barry I. Schneider

Vortex states of weakly interacting Bose-Einstein condensates confined in three-dimensional rotating harmonic traps are investigated numerically at zero temperature. The ground state in the rotating frame is obtained by propagating the Gross-Pitaevskii equation for the condensate in imaginary time. The total energies between states with and without a vortex are compared, yielding critical rotation frequencies that depend on the anisotropy of the trap and the number of atoms. Vortices displaced from the center of nonrotating traps are found to have long lifetimes for sufficiently large numbers of atoms. The relationship between vortex stability and bound core states is explored.

arXiv:cond-mat/9904269•Physical Review Letters•Jun 21, 1999•133 citations

Dynamics of vortex lattice formation in rotating two-component Bose–Einstein condensates

Toshiaki Kanai, Makoto Tsubota

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.

Journal of Physics Conference Series•Mar 1, 2018•1 citation
The dynamics and interaction of quantized vortices in Bose-Einstein condensates (BECs) are investigated by using the two-dimensional Gross-Pitaevskii equation (GPE) with/without an angular momentum rotation term. If all vortices have the same winding number, they would rotate around the trap center but never collide. In contrast, if the winding numbers are different, their interaction highly depends on the initial distance between vortex centers. The analytical results are presented to describe the dynamics of the vortex centers when β=0. While if β≠0, there is no analytical result but some conclusive numerical findings are provided for the further understanding of vortex interaction in BECs. Finally, the dynamic laws describing the relation of vortex interaction in nonrotating and rotating BECs are presented.
In this paper, we study the dynamics of rotating Bose--Einstein condensates (BEC) based on the Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and present an efficient and accurate algorithm for numerical simulations. We examine the conservation of the angular momentum expectation and the condensate width and analyze the dynamics of a stationary state with a shift in its center. By formulating the equation in either the two-dimensional polar coordinate system or the three-dimensional cylindrical coordinate system, the angular momentum rotation term becomes a term with constant coefficients. This allows us to develop an efficient time-splitting method which is time reversible, unconditionally stable, efficient, and accurate for the problem. Moreover, it conserves the position density. We also apply the numerical method to study issues such as the stability of central vortex states and the quantized vortex lattice dynamics in rotating BEC.
We study the stability of vortices in a binary system of Bose-Einstein condensat es, with their wave functions modeled by a set of coupled, time-dependent Gross-Pitaevskii equations. Beginning with an effective two-dimensional system, we identify miscible and immiscible regimes characterized by the inter- and intra-atomic interactions and the initial configuration of the system. We then consider a binary system of Bose-Einstein condensates placed in a rotating harmonic trap and study the single united vortex state in this system. We derive an approximate form for the energy of a single united vortex in the binary system and the critical angular velocity for the global stability of a vortex at the center of the trap. We also compute the metastability onset angular velocity for the local stability of a vortex at the center of the trap. Numerical solutions to the Gross-Pitaevskii equations support these expressions. These rotational results inform us of a novel bimodal formation within the miscible regime of the binary condensate system. We thus demonstrate the nontrivial aspects of vortex stability in interacting binary Bose-Einstein condensates as a result of their nonlinear interactions.
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 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.
Vortex states of weakly interacting Bose-Einstein condensates confined in three-dimensional rotating harmonic traps are investigated numerically at zero temperature. The ground state in the rotating frame is obtained by propagating the Gross-Pitaevskii equation for the condensate in imaginary time. The total energies between states with and without a vortex are compared, yielding critical rotation frequencies that depend on the anisotropy of the trap and the number of atoms. Vortices displaced from the center of nonrotating traps are found to have long lifetimes for sufficiently large numbers of atoms. The relationship between vortex stability and bound core states is explored.
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.