-
Pairwise Interactions of Ring Dark Solitons with Vortices and other Rings: Stationary States, Stability Features and Nonlinear Dynamics
Authors:
Wenlong Wang,
Theodore Kolokolnikov,
D. J. Frantzeskakis,
R. Carretero-González,
P. G. Kevrekidis
Abstract:
In the present work, we explore analytically and numerically the co-existence and interactions of ring dark solitons (RDSs) with other RDSs, as well as with vortices. The azimuthal instabilities of the rings are explored via the so-called filament method. As a result of their nonlinear interaction, the vortices are found to play a stabilizing role on the rings, yet their effect is not sufficient t…
▽ More
In the present work, we explore analytically and numerically the co-existence and interactions of ring dark solitons (RDSs) with other RDSs, as well as with vortices. The azimuthal instabilities of the rings are explored via the so-called filament method. As a result of their nonlinear interaction, the vortices are found to play a stabilizing role on the rings, yet their effect is not sufficient to offer complete stabilization of RDSs. Nevertheless, complete stabilization of the relevant configuration can be achieved by the presence of external ring-shaped barrier potentials. Interactions of multiple rings are also explored, and their equilibrium positions (as a result of their own curvature and their tail-tail interactions) are identified. In this case too, stabilization is achieved via multi-ring external barrier potentials.
△ Less
Submitted 17 May, 2021;
originally announced May 2021.
-
Multi-vortex crystal lattices in Bose-Einstein Condensates with a rotating trap
Authors:
Shuangquan Xie,
Panayotis G. Kevrekidis,
Theodore Kolokolnikov
Abstract:
We consider vortex dynamics in the context of Bose-Einstein Condensates (BEC) 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…
▽ More
We consider vortex dynamics in the context of Bose-Einstein Condensates (BEC) 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.
△ Less
Submitted 14 August, 2017;
originally announced August 2017.
-
Weakly Nonlinear Analysis of Vortex Formation in a Dissipative Variant of the Gross-Pitaevskii Equation
Authors:
Justin C. Tzou,
Panayotis G. Kevrekidis,
Theodore Kolokolnikov,
Ricardo Carretero-Gonzalez
Abstract:
For a dissipative variant of the two-dimensional Gross-Pitaevskii equation with a parabolic trap under rotation, we study a symmetry breaking process that leads to the formation of vortices. The first symmetry breaking leads to the formation of many small vortices distributed uniformly near the Thomas-Fermi radius. The instability occurs as a result of a linear instability of a vortex-free steady…
▽ More
For a dissipative variant of the two-dimensional Gross-Pitaevskii equation with a parabolic trap under rotation, we study a symmetry breaking process that leads to the formation of vortices. The first symmetry breaking leads to the formation of many small vortices distributed uniformly near the Thomas-Fermi radius. The instability occurs as a result of a linear instability of a vortex-free steady state as the rotation is increased above a critical threshold. We focus on the second subsequent symmetry breaking, which occurs in the weakly nonlinear regime. At slightly above threshold, we derive a one-dimensional amplitude equation that describes the slow evolution of the envelope of the initial instability. We show that the mechanism responsible for initiating vortex formation is a modulational instability of the amplitude equation. We also illustrate the role of dissipation in the symmetry breaking process. All analyses are confirmed by detailed numerical computations.
△ Less
Submitted 9 September, 2015;
originally announced September 2015.
-
Vortex Nucleation in a Dissipative Variant of the Nonlinear Schrödinger Equation under Rotation
Authors:
R. Carretero-Gonzalez,
P. G. Kevrekidis,
T. Kolokolnikov
Abstract:
In the present work, we motivate and explore the dynamics of a dissipative variant of the nonlinear Schr{ö}dinger equation under the impact of external rotation. As in the well established Hamiltonian case, the rotation gives rise to the formation of vortices. We show, however, that the most unstable mode leading to this instability scales with an appropriate power of the chemical potential $μ$ of…
▽ More
In the present work, we motivate and explore the dynamics of a dissipative variant of the nonlinear Schr{ö}dinger equation under the impact of external rotation. As in the well established Hamiltonian case, the rotation gives rise to the formation of vortices. We show, however, that the most unstable mode leading to this instability scales with an appropriate power of the chemical potential $μ$ of the system, increasing proportionally to $μ^{2/3}$. The precise form of the relevant formula, obtained through our asymptotic analysis, provides the most unstable mode as a function of the atomic density and the trap strength. We show how these unstable modes typically nucleate a large number of vortices in the periphery of the atomic cloud. However, through a pattern selection mechanism, prompted by symmetry-breaking, only few isolated vortices are pulled in sequentially from the periphery towards the bulk of the cloud resulting in highly symmetric stable vortex configurations with far fewer vortices than the original unstable mode. These results may be of relevance to the experimentally tractable realm of finite temperature atomic condensates.
△ Less
Submitted 1 December, 2014;
originally announced December 2014.
-
Drunken robber, tipsy cop: First passage times, mobile traps, and Hopf bifurcations
Authors:
Justin C. Tzou,
Shuangquan Xie,
Theodore Kolokolnikov
Abstract:
For a random walk on a confined one-dimensional domain, we consider mean first passage times (MFPT) in the presence of a mobile trap. The question we address is whether a mobile trap can improve capture times over a stationary trap. We consider two scenarios: a randomly moving trap and an oscillating trap. In both cases, we find that a stationary trap actually performs better (in terms of reducing…
▽ More
For a random walk on a confined one-dimensional domain, we consider mean first passage times (MFPT) in the presence of a mobile trap. The question we address is whether a mobile trap can improve capture times over a stationary trap. We consider two scenarios: a randomly moving trap and an oscillating trap. In both cases, we find that a stationary trap actually performs better (in terms of reducing expected capture time) than a very slowly moving trap; however, a trap moving sufficiently fast performs better than a stationary trap. We explicitly compute the thresholds that separate the two regimes. In addition, we find a surprising relation between the oscillating trap problem and a moving-sink problem that describes reduced dynamics of a single spike in a certain regime of the Gray-Scott model. Namely, the above-mentioned threshold corresponds precisely to a Hopf bifurcation that induces oscillatory motion in the location of the spike. We use this correspondence to prove the uniqueness of the Hopf bifurcation.
△ Less
Submitted 3 October, 2014;
originally announced October 2014.
-
Collective behaviour of large number of vortices in the plane
Authors:
Yuxin Chen,
Theodore Kolokolnikov,
Daniel Zhirov
Abstract:
We investigate the dynamics of $N$ point vortices in the plane, in the limit of large $N$. We consider {\em relative equilibria}, which are rigidly rotating lattice-like configurations of vortices. These configurations were observed in several recent experiments [Durkin and Fajans, Phys. Fluids (2000) 12, 289-293; Grzybowski {\em et.al} PRE (2001)64, 011603]. We show that these solutions and their…
▽ More
We investigate the dynamics of $N$ point vortices in the plane, in the limit of large $N$. We consider {\em relative equilibria}, which are rigidly rotating lattice-like configurations of vortices. These configurations were observed in several recent experiments [Durkin and Fajans, Phys. Fluids (2000) 12, 289-293; Grzybowski {\em et.al} PRE (2001)64, 011603]. We show that these solutions and their stability are fully characterized via a related {\em aggregation model} which was recently investigated in the context of biological swarms [Fetecau {\em et.al.}, Nonlinearity (2011) 2681; Bertozzi {\em et.al.}, M3AS (2011)]. By utilizing this connection, we give explicit analytic formulae for many of the configurations that have been observed experimentally. These include configurations of vortices of equal strength; the $N+1$ configurations of $N$ vortices of equal strength and one vortex of much higher strength; and more generally, $N+K$ configurations. We also give examples of configurations that have not been studied experimentally, including $N+2$ configurations where $N$ vortices aggregate inside an ellipse. Finally, we introduce an artificial ``damping'' to the vortex dynamics, in an attempt to explain the phenomenon of crystalization that is often observed in real experiments. The diffusion breaks the conservative structure of vortex dynamics so that any initial conditions converge to the lattice-like relative equilibrium.
△ Less
Submitted 31 January, 2014;
originally announced February 2014.
-
A Tale of Two Distributions: From Few To Many Vortices In Quasi-Two-Dimensional Bose-Einstein Condensates
Authors:
T. Kolokolnikov,
P. G. Kevrekidis,
R. Carretero-Gonzalez
Abstract:
Motivated by the recent successes of particle models in capturing the precession and interactions of vortex structures in quasi-two-dimensional Bose-Einstein condensates, we revisit the relevant systems of ordinary differential equations. We consider the number of vortices $N$ as a parameter and explore the prototypical configurations ("ground states") that arise in the case of few or many vortice…
▽ More
Motivated by the recent successes of particle models in capturing the precession and interactions of vortex structures in quasi-two-dimensional Bose-Einstein condensates, we revisit the relevant systems of ordinary differential equations. We consider the number of vortices $N$ as a parameter and explore the prototypical configurations ("ground states") that arise in the case of few or many vortices. In the case of few vortices, we modify the classical result of Havelock [Phil. Mag. ${\bf 11}$, 617 (1931)] illustrating that vortex polygons in the form of a ring are unstable for $N \geq7$. Additionally, we reconcile this modification with the recent identification of symmetry breaking bifurcations for the cases of $N=2,\dots,5$. We also briefly discuss the case of a ring of vortices surrounding a central vortex (so-called $N+1$ configuration). We finally examine the opposite limit of large $N$ and illustrate how a coarse-graining, continuum approach enables the accurate identification of the radial distribution of vortices in that limit.
△ Less
Submitted 1 May, 2014; v1 submitted 16 January, 2014;
originally announced January 2014.