A dissipative model of vortex motion within a two-dimensional Gross-Pitaevskii equation is studied. With the asymptotic coordination method of solutions, a system of ordinary first-order differential equations describing vortex motion in rotating Bose-Einstein condensate was derived. The model takes into account the effect of various external factors: magnetic trap frequencies, number of particles in the condensate, angular rotational velocity of the condensate, dissipation parameter, etc. In special cases of dissipation-free motion, the results are coordinated with the known results of other authors. The addition of dissipation generalizes the known equations and makes it possible to see the vortex motion to the equilibrium points and to define the equilibrium configurations of any number of vortices. The model is illustrated by a large number of examples.
A dissipative model of vortex motion within a two-dimensional Gross-Pitaevskii equation is studied. With the asymptotic coordination method of solutions, a system of ordinary first-order differential equations describing vortex motion in rotating Bose-Einstein condensate was derived. The model takes into account the effect of various external factors: magnetic trap frequencies, number of particles in the condensate, angular rotational velocity of the condensate, dissipation parameter, etc. In special cases of dissipation-free motion, the results are coordinated with the known results of other authors. The addition of dissipation generalizes the known equations and makes it possible to see the vortex motion to the equilibrium points and to define the equilibrium configurations of any number of vortices. The model is illustrated by a large number of examples.
We derive a general and exact equation of motion for a quantized vortex in an inhomogeneous two-dimensional Bose-Einstein condensate. This equation expresses the velocity of a vortex as a sum of local ambient density and phase gradients in the vicinity of the vortex. We perform Gross-Pitaevskii simulations of single-vortex dynamics in both harmonic and hard-walled disk-shaped traps, and find excellent agreement in both cases with our analytical prediction. The simulations reveal that, in a harmonic trap, the main contribution to the vortex velocity is an induced ambient phase gradient, a finding that contradicts the commonly quoted result that the local density gradient is the only relevant effect in this scenario. We use our analytical vortex velocity formula to derive a point-vortex model that accounts for both density and phase contributions to the vortex velocity, suitable for use in inhomogeneous condensates. Although good agreement is obtained between Gross-Pitaevskii and point-vortex simulations for specific few-vortex configurations, the effects of nonuniform condensate density are in general highly nontrivial, and are thus difficult to efficiently and accurately model using a simplified point-vortex description.
In two spatial dimensions, vortex-vortex interactions approximately vary with the logarithm of the inter-vortex distance, making it possible to describe an ensemble of vortices as a Coulomb gas. We introduce a duality between vortices in a quasi-two-dimensional (quasi-2D) scalar Bose-Einstein condensates (BEC) and effective Maxwell's electrodynamics. Specifically, we address the general scenario of inhomogeneous, time-dependent BEC number density with dissipation or rotation. Starting from the Gross-Pitaevskii equation (GPE), which describes the mean-field dynamics of a quasi-2D scalar BEC without dissipation, we show how to map vortices in a quasi-2D scalar BEC to 2D electrodynamics beyond the point-vortex approximation, even when dissipation is present or in a rotating system. The physical meaning of this duality is discussed.
The dynamics of interacting quantized vortex filaments in a rotating Bose–Einstein condensate existing in the Thomas–Fermi regime at zero temperature and obeying the Gross–Pitaevskii equation has been considered in the hydrodynamic “nonelastic” approximation. A noncanonical Hamilton equation of motion for the macroscopically averaged vorticity has been derived for a smoothly inhomogeneous array of filaments (vortex lattice) taking into account spatial nonuniformity of the equilibrium density of the condensate, which is determined by the trap potential. The minimum of the corresponding Hamiltonian describes the static configuration of the deformed vortex lattice against the preset density background. The condition of minimum can be reduced to a nonlinear second-order partial differential vector equation for which some exact and approximate solutions are obtained. It has been shown that if the condensate density has an anisotropic Gaussian profile, the equation of motion for the averaged vorticity has solutions in the form of a vector exhibiting a nontrivial time dependence, but homogeneous in space. An integral representation has also been obtained for the matrix Green function that determines the nonlocal Hamiltonian of a system of several quantized vortices of an arbitrary shape in a Bose–Einstein condensate with the Gaussian density. In particular, if all filaments are straight and oriented along one of the principal axes of the ellipsoid, we have a finitedimensional reduction that can describe the dynamics of the system of pointlike vortices against an inhomogeneous background. A simple approximate expression is proposed for the 2D Green function with an arbitrary density profile and is compared numerically with the exact result in the Gaussian case. The corresponding approximate equations of motion, describing the long-wavelength dynamics of interacting vortex filaments in condensates with a density depending only on transverse coordinates, have been derived.
We derive, from the Gross-Pitaevskii equation, an exact expression for the velocity of any vortex in a Bose-Einstein condensate, in equilibrium or not, in terms of the condensate wave function at the center of the vortex. In general, the vortex velocity is a sum of the local superfluid velocity, plus a correction related to the density gradient near the vortex. A consequence is that in rapidly rotating, harmonically trapped Bose-Einstein condensates, unlike in the usual situation in slowly rotating condensates and in hydrodynamics, vortices do not move with the local fluid velocity. We indicate how Kelvin's conservation of circulation theorem is compatible with the velocity of the vortex center being different from the local fluid velocity. Finally, we derive an exact wave function for a single vortex near the rotation axis in a weakly interacting system, from which we derive the vortex precession rate.
We have performed numerical simulations of various vortex configurations in a trapped BoseEinstein condensate by solving the two-dimensional Gross-Pitaevskii equation in the presence of a simple model of interaction between the condensate and the finite temperature thermal cloud that surrounds it. In that interaction the non-condensed thermal cloud acts as a source of dissipation with a damping effect of excitations. In the case ora single vortex and a vortex - anti vortex pair, we have found that the path of the vortices depends on the initial position, the initial separation distance if the case of two vortices and dissipation. This motion is periodic and it was found that sound waves are created by vortex motion; the intensity was stronger when the initial vortex separation distance was smaller. We have calculated the sound energy as the difference between the kinetic energy and the vortex energy. With no dissipation the vortices followed the same path with a slight oscillation due to the sound waves. We found that the smaller the initial vortex separation distance do is, the larger the sound production. The period, frequency, translation speed, sound energy and vortex energy were measured for different initial separation distances do and for different dissipation parameters 'Y. In the case of motion of one vortex, the connection between the dissipation 'Y and the friction coefficients, a and a' was studied as well. To create a simple turbulent state, we put eight pairs of vortex - anti vortex at random positions in the condensate with initial separation distance do =1.8 between them. 'Ve have studied the decay rate of the total energy, kinetic energy, quantum energy, trap energy and the z - component of the angular momentum together with the increase rate of the internal energy. Finally, we finished our investigation by putting randomly vortex - anti vortex pairs and studied the decrease of the number of vortices with time t. We found that the decrease is exponential.
Numerical simulations of vortex motion in a trapped Bose-Einstein condensate were performed by solving the two-dimensional Gross-Pitaevskii equation in the presence of a simple phenomenological model of interaction between the condensate and the finite temperature thermal cloud. The log (base e) of total energy, trap energy, quantum energy, kinetic energy, internal energy and z-component of the angular momentum vs time were compared with f(x)=a+bx for that time when the vortices come in in the condensate. The increasing/decay rate of these energies and L_{z} were studied as a function of dissipation.