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

arXiv:1309.7408 · Physical Review A

Dark-soliton dynamics and snake instability in superfluid Fermi gases trapped by an anisotropic harmonic potential

Wen Wen, Changqing Zhao, Xiaodong Ma

Published Dec 11, 2013• 22 citations
arXivHTMLPDF

Abstract

We present an investigation of generation, dynamics, and stability of dark solitons in anisotropic Fermi gases for a range of particle numbers and trap aspect ratios within the framework of the order-parameter equation. We calculate the periods of dark solitons oscillating in a trap and find good agreement with the results based on the Bogoliubov--de Gennes equations. By studying the stability of initially off-center dark solitons under various tight transverse confinements in the unitarity limit, we not only give the criterion of dynamical stability, but also find that the soliton and a hybrid of solitons and vortex rings can be characterized by different oscillation periods. The stability criterion is not fulfilled by the parameters of the recent experiment of Yefsah et al. [Nature (London) 499, 426 (2013)]. Therefore, instead of a very slow oscillation as observed experimentally, we find that the created dark soliton undergoes a transverse snake instability with a collapse into vortex rings, which propagate in solitonlike manner with a nearly 2 times larger period.

DOI
10.1103/physreva.88.063621

Related papers

Analytical solitonlike solutions and the dynamics of ultracold Fermi gases in a time-dependent three-dimensional harmonic potential

Ying Wang, Guosen Luo, Yu Long Zhou +1 more

We present a theoretical study of solitonlike solutions and their dynamics of ultracold superfluid Fermi gases trapped in a time-dependent three-dimensional (3D) harmonic potential with gain or loss. The 3D analytical solitonlike solutions are obtained without introducing any additional integrability constraints used elsewhere. The propagation of both bright- and dark-soliton-like solutions is investigated. We show that the amplitudes of dark-soliton-like solutions exhibit periodic oscillation, whereas those of the bright-soliton-like ones do not show such behavior. Moreover, we highlight that the oscillation periods of dark-soliton-like solutions predicted by our approach are matched very well with those observed in a recent experiment carried out by Yefsah et al. [T. Yefsah, A. T. Sommer, M. J. H. Ku, L. W. Cheuk, W. Ji, W. S. Bakr, and M. W. Zwierlein, Nature (London) 499, 426 (2013)NATUAS0028-083610.1038/nature12338] in both Bose-Einstein condensation and unitarity regimes.

Physical Review E•Sep 16, 2015•0 citations

Variational approach to the snake instability of a Bose-Einstein condensate soliton

U. Tanyeri, M. A. Gürkan, Ahmet Keleş +1 more

Solitons are striking manifestations of nonlinearity, encountered in diverse physical systems such as water waves, nonlinear optics, and Bose-Einstein condensates (BECs). In BECs, dark solitons emerge as exact stationary solutions of the one-dimensional Gross-Pitaevskii equation. While they can be long-lived in elongated traps, their stability is compromised in higher dimensions due to the snake instability, which leads to the decay of the soliton into vortex structures among other excitations. We investigate the dynamics of a dark soliton in a Bose-Einstein condensate confined in an anisotropic harmonic trap. Using a variational ansatz that incorporates both the transverse bending of the soliton plane and the emergence of vortices along the nodal line, we derive equations of motion governing the soliton's evolution. This approach allows us to identify stable oscillation modes as well as the growth rates of the unstable perturbations. In particular, we determine the critical trap anisotropy required to suppress the snake instability. Our analytical predictions are in good agreement with full numerical simulations of the Gross-Pitaevskii equation.

arXiv:2507.04926•Physical Review A•Dec 23, 2025•0 citations

Core filling and snaking instability of dark solitons in spin-imbalanced superfluid Fermi gases

Matthew D. Reichl, Erich J. Mueller

We use the time-dependent Bogoliubov--de Gennes equations to study dark solitons in three-dimensional spin-imbalanced superfluid Fermi gases. We explore how the shape and dynamics of dark solitons are altered by the presence of excess unpaired spins which fill their low-density core. The unpaired particles broaden the solitons and suppress the transverse snake instability. We discuss ways of observing these phenomena in cold-atom experiments.

arXiv:1703.04808•Physical Review A•May 26, 2017•9 citations

Snake instability of dark solitons in fermionic superfluids

Alberto Cetoli, Joachim Brand, R. G. Scott +2 more

We present numerical calculations of the snake instability in a Fermi superfluid within the Bogoliubov-de Gennes theory of the Bose-Einstein condensate (BEC) to BCS crossover using the random-phase approximation complemented by time-dependent simulations. We examine the snaking behavior across the crossover and quantify the time scale and length scale of the instability. While the dynamics shows extensive snaking before eventually producing vortices and sound on the BEC side of the crossover, the snaking dynamics is preempted by decay into sound due to pair breaking in the deep BCS regime. At the unitarity limit, hydrodynamic arguments allow us to link the rate of snaking to the experimentally observable ratio of inertial to physical mass of the soliton. In this limit we witness an unresolved discrepancy between our numerical estimates for the critical wave number of suppression of the snake instability and recent experimental observations with an ultracold Fermi gas.

arXiv:1307.3717•Physical Review A•Oct 29, 2013•51 citations

Solitonic vortices and the fundamental modes of the “snake instability”: Possibility of observation in the gaseous Bose-Einstein condensate

Joachim Brand, William P. Reinhardt

The connection between quantized vortices and dark solitons in a waveguidelike trap geometry is explored in the framework of the nonlinear Schr\"odinger equation. Variation of the transverse confinement leads from the quasi-one-dimensional (1D) regime, where solitons are stable, to 2D (or 3D) confinement, where soliton stripes are subject to a transverse modulational instability known as the ``snake instability.'' We present numerical evidence of a regime of intermediate confinement where solitons decay into single, deformed vortices with solitonic properties rather than vortex pairs as associated with the ``snake'' metaphor. Further relaxing the transverse confinement leads to the production of two and then three vortices, which correlates perfectly with a Bogoliubov stability analysis. The decay of a stationary dark soliton (or, planar node) into a single solitonic vortex is predicted to be experimentally observable in a 3D harmonically confined dilute-gas Bose-Einstein condensate.

arXiv:cond-mat/0105581•Physical Review A•Apr 1, 2002•159 citations

Crossover between snake instability and Josephson instability of dark solitons in superfluid Fermi gases

Wout Van Alphen, Hiromitsu Takeuchi, Jacques Tempere

Dark solitons in superfluid Bose gases decay through the snake instability mechanism, unless they are strongly confined. Recent experiments in superfluid Fermi gases have also interpreted soliton decay via this mechanism. However, we show using both an effective field numerical simulation and a perturbative analysis that there is a qualitative difference between soliton decay in the BEC and BCS regimes. On the BEC side of the interaction domain, the characteristic snaking deformations are induced by fluctuations of the amplitude of the order parameter, while on the BCS side, fluctuations of the phase destroy the soliton core through the formation of local Josephson currents without the presence of snaking deformations. The latter mechanism is qualitatively different from the snaking instability and this difference should be experimentally detectable.

arXiv:1901.10751•Physical Review A•Aug 27, 2019•6 citations
We present an investigation of generation, dynamics, and stability of dark solitons in anisotropic Fermi gases for a range of particle numbers and trap aspect ratios within the framework of the order-parameter equation. We calculate the periods of dark solitons oscillating in a trap and find good agreement with the results based on the Bogoliubov--de Gennes equations. By studying the stability of initially off-center dark solitons under various tight transverse confinements in the unitarity limit, we not only give the criterion of dynamical stability, but also find that the soliton and a hybrid of solitons and vortex rings can be characterized by different oscillation periods. The stability criterion is not fulfilled by the parameters of the recent experiment of Yefsah et al. [Nature (London) 499, 426 (2013)]. Therefore, instead of a very slow oscillation as observed experimentally, we find that the created dark soliton undergoes a transverse snake instability with a collapse into vortex rings, which propagate in solitonlike manner with a nearly 2 times larger period.
We present a theoretical study of solitonlike solutions and their dynamics of ultracold superfluid Fermi gases trapped in a time-dependent three-dimensional (3D) harmonic potential with gain or loss. The 3D analytical solitonlike solutions are obtained without introducing any additional integrability constraints used elsewhere. The propagation of both bright- and dark-soliton-like solutions is investigated. We show that the amplitudes of dark-soliton-like solutions exhibit periodic oscillation, whereas those of the bright-soliton-like ones do not show such behavior. Moreover, we highlight that the oscillation periods of dark-soliton-like solutions predicted by our approach are matched very well with those observed in a recent experiment carried out by Yefsah et al. [T. Yefsah, A. T. Sommer, M. J. H. Ku, L. W. Cheuk, W. Ji, W. S. Bakr, and M. W. Zwierlein, Nature (London) 499, 426 (2013)NATUAS0028-083610.1038/nature12338] in both Bose-Einstein condensation and unitarity regimes.
Solitons are striking manifestations of nonlinearity, encountered in diverse physical systems such as water waves, nonlinear optics, and Bose-Einstein condensates (BECs). In BECs, dark solitons emerge as exact stationary solutions of the one-dimensional Gross-Pitaevskii equation. While they can be long-lived in elongated traps, their stability is compromised in higher dimensions due to the snake instability, which leads to the decay of the soliton into vortex structures among other excitations. We investigate the dynamics of a dark soliton in a Bose-Einstein condensate confined in an anisotropic harmonic trap. Using a variational ansatz that incorporates both the transverse bending of the soliton plane and the emergence of vortices along the nodal line, we derive equations of motion governing the soliton's evolution. This approach allows us to identify stable oscillation modes as well as the growth rates of the unstable perturbations. In particular, we determine the critical trap anisotropy required to suppress the snake instability. Our analytical predictions are in good agreement with full numerical simulations of the Gross-Pitaevskii equation.
We present numerical calculations of the snake instability in a Fermi superfluid within the Bogoliubov-de Gennes theory of the Bose-Einstein condensate (BEC) to BCS crossover using the random-phase approximation complemented by time-dependent simulations. We examine the snaking behavior across the crossover and quantify the time scale and length scale of the instability. While the dynamics shows extensive snaking before eventually producing vortices and sound on the BEC side of the crossover, the snaking dynamics is preempted by decay into sound due to pair breaking in the deep BCS regime. At the unitarity limit, hydrodynamic arguments allow us to link the rate of snaking to the experimentally observable ratio of inertial to physical mass of the soliton. In this limit we witness an unresolved discrepancy between our numerical estimates for the critical wave number of suppression of the snake instability and recent experimental observations with an ultracold Fermi gas.
The connection between quantized vortices and dark solitons in a waveguidelike trap geometry is explored in the framework of the nonlinear Schr\"odinger equation. Variation of the transverse confinement leads from the quasi-one-dimensional (1D) regime, where solitons are stable, to 2D (or 3D) confinement, where soliton stripes are subject to a transverse modulational instability known as the ``snake instability.'' We present numerical evidence of a regime of intermediate confinement where solitons decay into single, deformed vortices with solitonic properties rather than vortex pairs as associated with the ``snake'' metaphor. Further relaxing the transverse confinement leads to the production of two and then three vortices, which correlates perfectly with a Bogoliubov stability analysis. The decay of a stationary dark soliton (or, planar node) into a single solitonic vortex is predicted to be experimentally observable in a 3D harmonically confined dilute-gas Bose-Einstein condensate.
Dark solitons in superfluid Bose gases decay through the snake instability mechanism, unless they are strongly confined. Recent experiments in superfluid Fermi gases have also interpreted soliton decay via this mechanism. However, we show using both an effective field numerical simulation and a perturbative analysis that there is a qualitative difference between soliton decay in the BEC and BCS regimes. On the BEC side of the interaction domain, the characteristic snaking deformations are induced by fluctuations of the amplitude of the order parameter, while on the BCS side, fluctuations of the phase destroy the soliton core through the formation of local Josephson currents without the presence of snaking deformations. The latter mechanism is qualitatively different from the snaking instability and this difference should be experimentally detectable.