Nucleation and kinematics of vortices in stirred Bose Einstein condensates
Jonas Rønning and Luiza Angheluta1
1The Njord Center, Department of Physics, University of Oslo, Blindern, 0316 Oslo, Norway
(Dated: October 27, 2022)
We apply the Halperin-Mazenco formalism within the Gross-Pitaevskii theory to characterise
the kinematics and nucleation of quantum vortices in a two-dimensional stirred Bose Einstein con-
densate. We introduce a smooth defect density field measuring the superfluid vorticity and is a
topologically conserved quantity. We use this defect density field and its associated current density
to study the precursory pattern formations that occur inside the repulsive potential of an obstacle
and determine the onset of vortex nucleation and shedding. We demonstrate that phase slips form
inside hard potentials even in the absence of vortex nucleation, whereas for soft potentials they
occur only above a critical stirring velocity leading to vortex nucleation. The Halperin-Mazenco
formalism provides an elegant and accurate method of deriving the point vortex dynamic directly
from the Gross-Pitaevskii equation.
Keywords: vortices, superfluid, Bose-Einstein condensat
I. INTRODUCTION
Topological defects are the fingerprints of broken con-
tinuous symmetries, and are widely encountered in or-
dered systems, such as disclinations in liquid crystals [1,
2], dislocations in solid crystals [3–5], orientational de-
fects in biological active matter [6–8] or quantized vor-
tices in quantum fluids [9–11], or cosmic strings [12]. The
formation and dynamics of topological defects during
phase ordering kinetics through temperature quenches
from the disordered phase have been well studied for
decades [13]. More recent work has been focused on de-
veloping theoretical frameworks to study non-equilibrium
pattern formations and dynamical regimes of ordered sys-
tems from the collective behavior of topological defects
beyond the phase-ordering kinetics.
The topological defects in an atomic Bose-Einstein
condensate (BEC) are quantum vortices where the con-
densate is locally melted while loosing its phase coher-
ence, and this induces persistent superfluid vortical cur-
rents outside the vortex cores. Hence, vortices are the
sole carrier of circulation in the condensate [14]. There is
active research both experimentally and theoretically on
understanding and tracking the non-thermal nucleation
and dynamics of quantized vortices in non-equilibrium
Bose Einstein condensates. Two main frameworks are
currently applied to study the creation of vortices either
by rotating the condensate [15–17] or by coupling the
condensate with a moving obstacle [18–22]. The nucle-
ation criterion is based on the energetic argument that
the superfluid flow reaches a critical velocity above which
the condensate phase gradient undergoes phase slips. In
rotating BEC systems, vortices are created when the sys-
tem spins at a uniform frequency above a critical thresh-
old determined by the quantized circulation of the vortex.
Vortices of same circulation form at the edge of the con-
densate and migrate into the bulk where they eventually
form vortex latices [15, 23–26].
The vortex nucleation in condensates stirred by a mov-
ing obstacle has also been studied [27–31] and observed
experimentally [21, 22]. Here, the nucleation criterion
relies on the height U0of the repulsive potential repre-
senting the coupling of the stirring obstacle to the con-
densate. A hard potential corresponds to an almost im-
penetrable obstacle when U0> µ, where µis the conden-
sate chemical potential, such that the condensate density
rapidly decreases and nearly vanishes inside the poten-
tial. By contrast, a soft potential corresponds to a pene-
trable obstacle for U0< µ such that the condensate den-
sity is gently depleted inside the obstacle. The onset of
vortex nucleation induced by a hard obstacle occurs when
the local condensate velocity reaches the critical velocity
for phonon emission, whereas for the soft obstacle this
is a necessary, but not a sufficient requirement [27, 28].
Stirring obstacles are typically modelled as Gaussian po-
tentials with varying height and width [30, 31], which
approach a Dirac-delta function in the limit of a solid
obstacle. In Ref. [31], the vortex nucleation induced by
a repulsive Gaussian potential of different strengths is
studied numerically. It is found that near the critical
velocity for vortex nucleation, the energy gap between
the ground state and the exited state goes to zero as
a power-law, while ghost vortices, i.e. phase slips, are
formed inside the potential. By contrast, no such ghost
vortices develop in the case of soft potentials. In addi-
tion to tuning the degree of permeability of the obstacle,
different vortex shedding regimes, from vortex dipoles,
pairs and clusters [32–34], can be induced by varying the
size of the obstacle through the width of the potential
which also changes the critical stirring velocity [21, 35].
Once vortices are being shed into the condensate, they
interact with each other forming dynamic clusters that
sustain energy cascades and two-dimensional quantum
turbulence [11, 34, 36–39].
Even through compressibility effects, due to shock
waves and phonons, are particularly important in the nu-
cleation and the annihilation of vortex dipoles, they are
typically overlooked in the quantum turbulence regime
where turbulent energy spectra and clustering behavior
is attributed mostly to the mutual interactions between
arXiv:2210.14620v1 [cond-mat.quant-gas] 26 Oct 2022