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:2507.04926 · Physical Review A

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

U. Tanyeri, M. A. Gürkan, Ahmet Keleş, Mehmet Özgür Oktel

Published Dec 23, 2025• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/2wc6-pp3r

Related papers

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

Dark-soliton states of Bose-Einstein condensates in anisotropic traps

David L. Feder, M. S. Pindzola, Lee A. Collins +2 more

Dark soliton states of Bose-Einstein condensates in harmonic traps are studied both analytically and computationally by the direct solution of the Gross-Pitaevskii equation in three dimensions. The ground and self-consistent excited states are found numerically by relaxation in imaginary time. The energy of a stationary soliton in a harmonic trap is shown to be independent of density and geometry for large numbers of atoms. Large-amplitude field modulation at a frequency resonant with the energy of a dark soliton is found to give rise to a state with multiple vortices. The Bogoliubov excitation spectrum of the soliton state contains complex frequencies, which disappear for sufficiently small numbers of atoms or large transverse confinement. The relationship between these complex modes and the snake instability is investigated numerically by propagation in real time.

arXiv:cond-mat/0004449•Physical Review A•Oct 19, 2000•270 citations

Stability of dark solitons in a bubble Bose-Einstein condensate

Raphael Wictky Sallatti, Lauro Tomio, Dmitry E. Pelinovsky +1 more

The stability of nonlinear waves on curved surfaces is a problem of widespread interest across physics. Here, we establish the stability criteria for dark solitons on a spherical Bose-Einstein condensate. We demonstrate a sharp instability threshold in the nonlinear parameter, beyond which solitons decay into vortex dipoles via snake instabilities. Analytically and numerically, we prove this decay is dictated by a single unstable mode for each angular momentum $m \geq 2$, which is a universal mechanism that controls the resulting vortex state. Unlike in the full three-dimensional case, where snake instabilities lead to vortex rings, a dark soliton confined to the surface of a bubble can only decay into vortex pairs.

arXiv:2511.04385•Physical Review A•Mar 16, 2026•0 citations

Snake instability of ring dark solitons in toroidally trapped Bose-Einstein condensates

L. A. Toikka, Kalle‐Antti Suominen

We show by numerically solving the Bogoliubov--de Gennes equations and numerically integrating the Gross-Pitaevskii equation that the onset of the snake instability of ring dark solitons requires a broken symmetry. We elucidate explicitly the connection between imaginary Bogoliubov modes and the snake instability, predicting the number of vortex-antivortex pairs produced. In addition, we propose a simple model to give a physical motivation as to why the snake instability takes place and needs a broken symmetry. Finally, we show that tight confinement in a toroidal potential actually enhances soliton decay due to inhibition of soliton motion, but can lead to a periodic revival of the ring dark soliton after the first snake instability.

arXiv:1210.3994•Physical Review A•Apr 3, 2013•37 citations

Stability analysis and attractor dynamics of three-dimensional dark solitons with localized dissipation

Christian Baals, Alexandre Gil Moreno, Jian Jiang +2 more

We study the stability and the attractor dynamics of an elongated Bose-Einstein condensate (BEC) with dark or gray kink solitons in the presence of localized dissipation. To this end, the three-dimensional Gross--Pitaevskii equation with an additional imaginary potential is solved numerically. We analyze the suppression of the snaking instability in dependence of the dissipation strength and extract the threshold value for the stabilization of the dark soliton for experimentally realistic parameters. Below the threshold value, we observe the decay into a solitonic vortex. Above the stabilization threshold, we observe the attractor dynamics towards the dark soliton when initially starting from a gray soliton. We find that for all initial conditions the dark soliton is the unique steady-state of the system---even when starting from the BEC ground state.

arXiv:2012.11487•Physical Review A•Apr 2, 2021•12 citations

Ring Dark Solitons in Toroidal Bose-Einstein Condensates

L. A. Toikka

In this Thesis, we study various aspects of ring dark solitons (RDSs) in quasi-two-dimensional toroidally trapped Bose-Einstein condensates, focussing on atomic realisations thereof. Unlike the well-known planar dark solitons, exact analytic expressions for RDSs are not known. We address this problem by presenting exact localized soliton-like solutions to the radial Gross-Pitaevskii equation. To date, RDSs have not been experimentally observed in cold atomic gases, either. To this end, we propose two protocols for their creation in experiments. It is also currently well known that in dimensions higher than one, (ring) dark solitons are susceptible, in general, to an irreversible decay into vortex-antivortex pairs through the snake instability. We show that the snake instability is caused by an unbalanced quantum pressure across the soliton's notch, linking the instability to the Bogoliubov-de Gennes spectrum. In particular, if the angular symmetry is maintained (or the toroidal trapping is restrictive enough), we show that the RDS is stable (long-lived with a lifetime of order seconds) in two dimensions. Furthermore, when the decay does take place, we show that the snake instability can in fact be reversible, and predict a previously unknown revival phenomenon for the original (many-)RDS system: the soliton structure is recovered and all the point-phase singularities (i.e. vortices) disappear. Eventually, however, the decay leads to an example of quantum turbulence; a quantum example of the laminar-to-turbulent type of transition.

•Jun 18, 2014•0 citations
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.
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 soliton states of Bose-Einstein condensates in harmonic traps are studied both analytically and computationally by the direct solution of the Gross-Pitaevskii equation in three dimensions. The ground and self-consistent excited states are found numerically by relaxation in imaginary time. The energy of a stationary soliton in a harmonic trap is shown to be independent of density and geometry for large numbers of atoms. Large-amplitude field modulation at a frequency resonant with the energy of a dark soliton is found to give rise to a state with multiple vortices. The Bogoliubov excitation spectrum of the soliton state contains complex frequencies, which disappear for sufficiently small numbers of atoms or large transverse confinement. The relationship between these complex modes and the snake instability is investigated numerically by propagation in real time.
The stability of nonlinear waves on curved surfaces is a problem of widespread interest across physics. Here, we establish the stability criteria for dark solitons on a spherical Bose-Einstein condensate. We demonstrate a sharp instability threshold in the nonlinear parameter, beyond which solitons decay into vortex dipoles via snake instabilities. Analytically and numerically, we prove this decay is dictated by a single unstable mode for each angular momentum m≥2m \geq 2m≥2, which is a universal mechanism that controls the resulting vortex state. Unlike in the full three-dimensional case, where snake instabilities lead to vortex rings, a dark soliton confined to the surface of a bubble can only decay into vortex pairs.
We show by numerically solving the Bogoliubov--de Gennes equations and numerically integrating the Gross-Pitaevskii equation that the onset of the snake instability of ring dark solitons requires a broken symmetry. We elucidate explicitly the connection between imaginary Bogoliubov modes and the snake instability, predicting the number of vortex-antivortex pairs produced. In addition, we propose a simple model to give a physical motivation as to why the snake instability takes place and needs a broken symmetry. Finally, we show that tight confinement in a toroidal potential actually enhances soliton decay due to inhibition of soliton motion, but can lead to a periodic revival of the ring dark soliton after the first snake instability.
We study the stability and the attractor dynamics of an elongated Bose-Einstein condensate (BEC) with dark or gray kink solitons in the presence of localized dissipation. To this end, the three-dimensional Gross--Pitaevskii equation with an additional imaginary potential is solved numerically. We analyze the suppression of the snaking instability in dependence of the dissipation strength and extract the threshold value for the stabilization of the dark soliton for experimentally realistic parameters. Below the threshold value, we observe the decay into a solitonic vortex. Above the stabilization threshold, we observe the attractor dynamics towards the dark soliton when initially starting from a gray soliton. We find that for all initial conditions the dark soliton is the unique steady-state of the system---even when starting from the BEC ground state.
In this Thesis, we study various aspects of ring dark solitons (RDSs) in quasi-two-dimensional toroidally trapped Bose-Einstein condensates, focussing on atomic realisations thereof. Unlike the well-known planar dark solitons, exact analytic expressions for RDSs are not known. We address this problem by presenting exact localized soliton-like solutions to the radial Gross-Pitaevskii equation. To date, RDSs have not been experimentally observed in cold atomic gases, either. To this end, we propose two protocols for their creation in experiments. It is also currently well known that in dimensions higher than one, (ring) dark solitons are susceptible, in general, to an irreversible decay into vortex-antivortex pairs through the snake instability. We show that the snake instability is caused by an unbalanced quantum pressure across the soliton's notch, linking the instability to the Bogoliubov-de Gennes spectrum. In particular, if the angular symmetry is maintained (or the toroidal trapping is restrictive enough), we show that the RDS is stable (long-lived with a lifetime of order seconds) in two dimensions. Furthermore, when the decay does take place, we show that the snake instability can in fact be reversible, and predict a previously unknown revival phenomenon for the original (many-)RDS system: the soliton structure is recovered and all the point-phase singularities (i.e. vortices) disappear. Eventually, however, the decay leads to an example of quantum turbulence; a quantum example of the laminar-to-turbulent type of transition.