The dynamics of interacting quantum vortices in a quasi-two-dimensional spatially inhomogeneous Bose–Einstein condensate, whose equilibrium density vanishes at two points of the plane with a possible presence of an immobile vortex with a few circulation quanta at each point, has been considered in a hydrodynamic approximation. A special class of density profiles has been chosen, so that it proves possible to calculate analytically the velocity field produced by point vortices. The equations of motion have been given in a noncanonical Hamiltonian form. The theory has been generalized to the case where the condensate forms a curved quasi-two-dimensional shell in the three-dimensional space.
The dynamics of interacting quantum vortices in a quasi-two-dimensional spatially inhomogeneous Bose–Einstein condensate, whose equilibrium density vanishes at two points of the plane with a possible presence of an immobile vortex with a few circulation quanta at each point, has been considered in a hydrodynamic approximation. A special class of density profiles has been chosen, so that it proves possible to calculate analytically the velocity field produced by point vortices. The equations of motion have been given in a noncanonical Hamiltonian form. The theory has been generalized to the case where the condensate forms a curved quasi-two-dimensional shell in the three-dimensional space.
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 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.
Abstract Understanding quantum dynamics in a two-dimensional Bose–Einstein condensate relies on understanding how vortices interact with each other microscopically and with local imperfections of the potential which confines the condensate. Within a system consisting of many vortices, the trajectory of a vortex–antivortex pair is often scattered by a third vortex, an effect previously characterised. However, the natural question remains as to how much of this effect is due to the velocity induced by this third vortex and how much is due to the density inhomogeneity which it introduces. In this work, we describe the various qualitative scenarios which occur when a vortex–antivortex pair interacts with a smooth density impurity whose profile is identical to that of a vortex but lacks the circulation around it.
By solving the two-dimensional, nonlocal Gross-Pitaevskii equation numerically, we study the hydrodynamics of quantized vortices in a dipolar Bose-Einstein condensate excited by an oscillating potential. We calculate the critical velocity for the vortex nucleation in the presence and absence of dipole-dipole interaction. Above a critical velocity, the oscillating potential nucleates the vortex pairs, whose impulses alternatively change their direction synchronously with the oscillation of the Gaussian potential. We study the dynamics of the vortices and also the formation of the dark solitons which have not been previously reported in a dipolar Bose-Einstein condensate.
The aim of this paper is to derive the equations of motion of quantum vortices in a weakly non-ideal Bose–Einstein condensate (BEC) at finite temperatures and in the presence of macroscopic superfluid and normal velocities. The motivation for this work is that the dynamics of quantum vortices in superfluid helium has been studied at the phenomenological level, and a number of issues require clarification, which can be addressed by studying the corresponding problem in BEC. In BEC, quantum vortices are located at points where the macroscopic wave function (MWF) ψ(r,T), which describes the BEC dynamics, vanishes ψ(r,T)=0. The intersection of the surfaces Reψ(r,T)=0 and Imψ(r,T)=0 forms a one-dimensional line, which is just the quantum vortex line. During the evolution of the MWF, the position of the zero surfaces changes, that is, the vortex lines move. Conversely, the motion of these vortex filaments determines the evolution of the MWF. Thus, there is a mutual correspondence, and the quantity ψ(r,T) can be considered as a functional of the entire vortex configuration, i.e., {s(ξ,t)}. This viewpoint allows us to derive the equations of motion for vortices from the Gross–Pitaevskii equation (GPE), which describes the dynamics of MWF. The key finding is the confirmation that the equation derived from the fundamental quantum description (the GPE) for a BEC has the same mathematical structure as the equation derived from macroscopic, phenomenological reasoning for superfluid helium.
In this paper, we investigate the dynamics of the vortices in the homogeneous two-dimensional Bose–Einstein condensate by considering the condensate density variation. It is shown that the condensate density variation plays an important role when two vortices are close together. By the time-dependent Lagrangian functional method and the numerical simulation, we study how the interaction strength of the vortices varies with the distance between two vortices with the same or different topological charges, and generalize the results for the vortex pairs to the multi-vortex case. The numerical simulations demonstrate the validity of our theory.