Stability of dark solitons in a bubble Bose-Einstein condensate

Basic Information

• Paper ID: 2511.04385
• Title: Stability of dark solitons in a bubble Bose-Einstein condensate
• Authors: Raphael Wictky Sallatti, Lauro Tomio, Dmitry E. Pelinovsky, Arnaldo Gammal
• Institutions: Universidade de São Paulo, Universidade Estadual Paulista, McMaster University
• Classification: cond-mat.quant-gas (Condensed Matter Physics - Quantum Gases)
• Publication Date: November 7, 2025
• Paper Link: https://arxiv.org/abs/2511.04385v1

Abstract

This paper investigates the dynamical stability of dark solitons confined to the two-dimensional spherical bubble surface of a Bose-Einstein condensate (BEC). In the spherical geometry of a BEC, when the interaction parameter ε ≳ 8.37, dark solitons become unstable due to snake-like instability driven by discrete angular modes, producing vortex dipoles. Through analytical and numerical methods, the study demonstrates that for each angular mode m ≥ 2, there exists exactly one unstable mode, whose dominance determines the number of m vortex dipoles. Time-dependent simulations confirm the formation of vortex dipoles.

Research Background and Motivation

Research Problem

This paper addresses the stability of dark solitons in Bose-Einstein condensates on the surface of a spherical bubble, with particular focus on how spherical geometry affects dark soliton dynamics and instabilities.

Problem Significance

1. Experimental Relevance: Microgravity experiments on the International Space Station have realized ultracold atomic gases in spherical geometry, while ground-based experiments have achieved shell-like bubble structures through two-component systems
2. Fundamental Physics: Spherical closed topology provides a unique platform for studying low-dimensional quantum gases, involving vortex dynamics, BKT phase transitions, and other fundamental physical phenomena
3. Geometric Effects: The curvature and closed topology of the spherical surface have essential impacts on soliton propagation and instabilities, differing fundamentally from planar systems

Limitations of Existing Research

1. Differences from Planar Systems: In planar quasi-one-dimensional or quasi-two-dimensional BECs, dark solitons are susceptible to snake-like instability and decay into vortices, but the closed topology of the sphere alters the decay pathway
2. Topological Constraints: Single vortices cannot exist on a sphere; vortices must appear in pairs with zero total charge (Poincaré-Hopf theorem)
3. Missing Systematic Analysis: Lack of systematic analytical and numerical investigation of dark soliton stability in spherical geometry

Research Motivation

Combining recent progress in ISS BEC experiments with theoretical requirements, this work systematically investigates the spectral stability of dark solitons on spherical bubble surfaces, particularly how discrete angular modes drive instabilities and lead to vortex pair formation.

Core Contributions

1. Determination of Stability Threshold: First precise determination of the stability threshold for dark solitons in spherical bubble BEC as ε ≈ 8.37
2. Analytical Theoretical Framework: Establishment of a complete Bogoliubov-de Gennes (BdG) spectral stability analysis framework, proving that:
   • For each angular mode m ≥ 2, there exists exactly one unstable mode
   • Provides asymptotic formula for instability threshold: εₘᵗʰ = 4m(m-1)
3. Vortex Dipole Number Prediction: Establishes direct correspondence between the dominant unstable mode m and the number of vortex-antivortex pairs produced
4. Numerical Verification: Complete verification of analytical predictions through time-evolution simulations, demonstrating the dynamics of dark soliton breakup into m vortex pairs
5. Small and Large Parameter Asymptotics: Provides analytical expressions for dark soliton profiles in the limits ε → 0 and ε → ∞

Detailed Methods

Task Definition

Study the linear stability of dark soliton solutions in BEC on a two-dimensional sphere (radius R, thickness δR << R), determining stability characteristics of different angular modes at various interaction strengths.

Input: Interaction parameter ε = g/(2π) = 4√(2π)aN/(δR), where a is the s-wave scattering length and N is the number of atoms

Output: Stability spectrum, instability thresholds, dynamics of vortex pair formation

Mathematical Model

1. Two-Dimensional Gross-Pitaevskii Equation

In spherical coordinates (θ, φ), the dimensionless GPE is:

i∂ₜψ = -Δ₂Dψ + g|ψ|²ψ

where the Laplacian operator is:

Δ₂D = (1/sinθ)∂θ(sinθ∂θ) + (1/sin²θ)∂²φ

2. Dark Soliton Steady-State Solution

Seeking solutions of the form ψₛ(θ,φ,t) = f(θ)/√(2π) e^(-iμt), where f(θ) satisfies:

-[d²/dθ² + cotθ d/dθ]f(θ) + ε|f(θ)|²f(θ) = μf(θ)

Normalization condition: ∫₀^π sinθ|f(θ)|²dθ = 1

Dark soliton characteristics: f(θ) vanishes at θ=π/2 and decreases monotonically

3. Small Parameter Asymptotics (ε → 0)

f(θ) = √(3/2)cosθ + εf₁(θ) + O(ε²)
μ = 2 + (9/10)ε + O(ε²)

where f₁(θ) ∝ P₃(cosθ) (third-order Legendre polynomial)

4. Large Parameter Asymptotics (ε → ∞)

Dark soliton concentrates near the equator. Introducing variable z = √(ε/2)(π/2 - θ):

f(θ) = f∞(θ) + f∞⁽¹⁾(θ)/√ε + O(1/ε)
f∞(θ) = (1/√2)tanh[√(ε/2)(π/2 - θ)]
μ = ε/2 + √ε + O(1)

Stability Analysis Method

1. Bogoliubov-de Gennes Method

Considering small perturbations to the steady-state solution:

ψ(θ,φ,t) = ψₛ + (e^(-iμt)/√(2π))[u(θ,φ,t) + iv(θ,φ,t)]

After linearization, coupled equations are obtained:

∂ₜu = -Δ₂Dv + εf²v - μv
-∂ₜv = -Δ₂Du + 3εf²u - μu

2. Angular Mode Decomposition

Expanding perturbations by angular modes: u = ûₘ(θ)e^(i(mφ+ωt)), v = v̂ₘ(θ)e^(i(mφ+ωt))

Obtaining the spectral problem:

ωûₘ = L⁻ₘv̂ₘ,  L⁻ₘ = -Δₘ + εf² - μ
ωv̂ₘ = L⁺ₘûₘ,  L⁺ₘ = -Δₘ + 3εf² - μ

where Δₘ = d²/dθ² + cotθ d/dθ - m²/sin²θ

3. Stability Criterion

• If all eigenvalues of L±ₘ are strictly positive, then all ω are real and the dark soliton is stable for mode m
• If L⁻ₘ has negative eigenvalues, imaginary frequencies Im(ω) ≠ 0 appear, triggering instability

Technical Innovations

1. Application of Comparison Principle

Using two key inequalities:

L⁺ₘ - L⁻ₘ = 2εf²(θ) ≥ 0
L±ₘ₊₁ - L±ₘ = (2m+1)/sin²θ ≥ 0

Rigorously proving the monotonicity of instability thresholds: εₘ < εₘ₊₁

2. Mode Classification Analysis

• m=0: Negative energy eigenvalues exist, but instability not observed numerically
• m=1: All eigenvalues are real, no instability bifurcation
• m≥2: For each m, exactly one unstable mode exists

3. Asymptotic Expansion Technique

For large m limit, through perturbation theory deriving:

εₘ = 4m(m-1) + O(1)

This formula has error less than 5% for m=2

4. Combined Numerical Methods

• Shooting Method: Solving nonlinear boundary value problems to obtain dark soliton profiles
• Finite Difference: Discretizing operators to compute eigenvalues
• Split-Step Fourier Method: Combined with Crank-Nicolson scheme for time evolution

Experimental Setup

Numerical Solution Methods

1. Dark Soliton Profile Calculation

Using shooting method combined with secant method:

• Redefine f̃(θ) = f(θ)√ε
• Boundary conditions: f̃'(0) = 0, f̃(π/2) = 0
• For given μ, adjust initial value f̃(0) until boundary conditions are satisfied
• Solve ODE using Runge-Kutta method
• Scan μ range through continuation method

2. Eigenvalue Computation

• Discretize operators L±ₘ using finite differences
• Grid points: up to 800
• Compute smallest eigenvalues to determine stability

3. Time Evolution Simulation

Spatial Discretization:

• θ direction: 256 grid points, step size δθ = π/256 ≈ 0.0123
• φ direction: 256 grid points, step size δφ = 2π/256 ≈ 0.0245

Temporal Discretization:

• Time step: δt = 10⁻⁵

Numerical Algorithm:

• Fast Fourier Transform (FFT) for φ direction
• Crank-Nicolson finite difference for θ direction
• Split-step operator technique for nonlinear terms
• Complete scheme: ψ(t+δt) ≃ e^(-ig|ψ|²δt/2) FFT⁻¹ CN FFT ψ(t)e^(-ig|ψ|²δt/2)

Parameter Settings

Stability Thresholds

Critical values for each mode determined through numerical computation:

Dynamics Simulation Parameters

Three typical cases selected to demonstrate different dominant modes:

• m=2 dominant: ε=20
• m=3 dominant: ε=50
• m=4 dominant: ε=100

Experimental Results

Main Results

1. Stability Phase Diagram

Key Findings:

• Stable Region: ε ≲ 8.37, dark solitons stable for all angular modes
• Unstable Region: ε ≳ 8.37, m≥2 modes begin to excite
• Mode Transitions:
  • 8.37 < ε < 35: m=2 dominant
  • ε ≈ 35: m=2 and m=3 modes compete Im(ω₂) = Im(ω₃)
  • ε > 35: m=3 begins to dominate

2. Eigenvalue Evolution

Figure 2 displays key results:

• Minimum eigenvalue of L⁻ₘ (dashed lines): Crosses zero at ε=εₘ
• Im(ωₘ) (solid lines): Grows from zero, characterizing instability growth rate
• For given ε, maximum Im(ωₘ) corresponds to dominant unstable mode

3. Theoretical Prediction Verification

Asymptotic formula εₘᵗʰ = 4m(m-1) shows excellent agreement with numerical results:

• m=2: Theory 8 vs numerical 8.367 (error 4.6%)
• m=3: Theory 24 vs numerical 24.402 (error 1.7%)
• m≥4: Error less than 1%

Dynamics Evolution Analysis

Case 1: m=2 Dominant (ε=20)

Time Evolution Characteristics:

• t=0: Dark soliton at equator, width ~1/√ε ≈ 0.22
• Intermediate stage: Snake-like instability appears
• Final state: Formation of 2 vortex-antivortex pairs

Density Distribution:

• Density at equator |ψ(π/2, φ)|² shows 2-fold periodic modulation
• 3D visualization shows 2 symmetrically distributed vortex structures

Case 2: m=3 Dominant (ε=50)

• Dark soliton width: ~0.14
• After breakup forms 3 vortex-antivortex pairs
• Density modulation exhibits 3-fold symmetry

Case 3: m=4 Dominant (ε=100)

• Dark soliton width: ~0.10
• After breakup forms 4 vortex-antivortex pairs
• Density modulation exhibits 4-fold symmetry

Manifestation of Topological Constraints

Poincaré-Hopf Theorem Verification:

• Any continuous tangent vector field on a sphere must have at least one zero point
• For superfluid implies single vortices cannot exist
• Vortices must appear in pairs with zero total circulation
• During instability, vortex number changes by integer multiples of +2

Dark Soliton Profile Characteristics

Small ε Limit Verification

Figure 1(a) upper panel shows f(θ) approaching √(3/2)cosθ, consistent with theoretical prediction (Eq. S10)

Large ε Limit Verification

Figure 1(a) upper panel shows f(θ) concentrated near equator, approaching hyperbolic tangent function f∞(θ)

Chemical Potential Dependence

Figure 1(b) displays complete μ(ε) curve:

• Small ε segment: Linear growth μ ≈ 2 + 0.9ε (dashed line)
• Large ε segment: Square root growth μ ≈ ε/2 + √ε + 1 (circles)
• Numerical results (solid line) smoothly connect two asymptotic regions across entire range

Related Work

Spherical Geometry BEC Research

1. Experimental Progress:
   • Aveline et al. (2020): First observation of orbital BEC on ISS
   • Carollo et al. (2022): Observation of ultracold atomic gas bubbles in microgravity
   • Jia et al. (2022): Ground-based two-component system realizing shell BEC
2. Theoretical Foundation:
   • Zobay & Garraway (2001-2004): Theory of atomic trapping in field-induced adiabatic potential wells
   • Tononi & Salasnich (2019): BEC theory on spherical surfaces
   • Sun et al. (2018): Static and dynamical properties of shell condensates

Vortex Dynamics

1. Spherical Vortices:
   • Padavić et al. (2020): Vortex-antivortex physics in shell BEC
   • Caracanhas et al. (2022): Superfluid vortex dynamics on ellipsoidal and rotational surfaces
   • Tononi et al. (2024): Quantum vortices in curved geometry
2. Topological Constraints:
   • Tononi et al. (2022): Topological superfluid phase transitions in bubble-trapped condensates
   • Emphasizing Poincaré-Hopf theorem constraints on vortex configurations

Dark Soliton Research

1. Planar Systems:
   • Burger et al. (1999): Experimental observation of dark solitons in BEC
   • Busch & Anglin (2000): Dark soliton motion in trapped BEC
   • Feder et al. (2000): Dark soliton states in anisotropic traps
2. Snake-like Instability:
   • Gaidoukov & Anglin (2021): BdG theory of snake-like instability in higher-dimensional gray solitons
   • Anderson et al. (2001): Observation of dark soliton decay into vortex rings
3. Ring Geometry:
   • Kivshar & Yang (1994): Ring-shaped dark solitons
   • Theocharis et al. (2003): Ring-shaped dark solitons and vortex necklaces in BEC

Unique Contributions of This Work

Compared to existing work, this paper:

1. First systematic study of complete stability spectrum of dark solitons in spherical bubble BEC
2. Establishes analytical framework precisely predicting instability thresholds for each angular mode
3. Quantitatively correlates dominant unstable mode with number of vortex pairs
4. Complete verification of theoretical predictions across full parameter range from small to large ε

Conclusions and Discussion

Main Conclusions

1. Stability Threshold: Dark solitons in spherical bubble BEC are stable for ε ≲ 8.37; beyond this value, m≥2 angular modes drive snake-like instability
2. Instability Mechanism:
   • Each m≥2 mode has exactly one unstable eigenvalue
   • Instability thresholds satisfy strict monotonicity εₘ < εₘ₊₁
   • Asymptotic formula εₘ ≈ 4m(m-1) precisely predicts large m behavior
3. Vortex Pair Formation: Dominant unstable mode m directly determines number of vortex-antivortex pairs produced; time-evolution simulations completely verify this prediction
4. Topological Constraints: Closed topology of sphere requires vortices to appear in pairs; vortex pairs cannot escape to boundaries, leading to complex vortex dynamics

Limitations

1. Theoretical Assumptions:
   • Two-dimensional approximation requires δR << R
   • Neglects radial excitations (higher energy)
   • Ignores thermal cloud effects causing energy instability
   • Mean-field approximation neglects quantum fluctuations
2. m=0 Mode: Theory predicts possible instability, but not observed numerically; requires deeper investigation
3. Experimental Conditions:
   • Difficult to precisely control ε parameter in actual experiments
   • Limited time for microgravity experiments
   • Must consider finite temperature effects
4. Long-Time Dynamics:
   • Only studies initial instability and vortex pair formation
   • Does not deeply analyze subsequent vortex pair evolution (annihilation, stable configurations)

Future Directions

1. Experimental Verification:
   • Utilize ISS platform or ground-based two-component systems to verify stability thresholds
   • Observe correspondence between vortex pair numbers at different ε and theoretical predictions
2. Finite Temperature Effects:
   • Include thermal cloud effects on dark soliton stability
   • Study competition between energy and dynamical instabilities
3. Vortex Dynamics:
   • Deeply investigate long-time evolution after vortex pair formation
   • Analyze mechanisms of vortex-antivortex annihilation and stable configuration formation
4. Two-Component Systems:
   • Extend to two-component BEC mixtures
   • Study stability enhancement effects of filled solitons
5. Ellipsoidal Geometry:
   • Generalize to non-spherical closed surfaces
   • Study effects of geometric asymmetry on stability

In-Depth Evaluation

Strengths

1. Theoretical Rigor

• Complete Mathematical Framework: From GPE to BdG method establishing complete linear stability theory
• Rigorous Proofs: Using comparison principle and Sturm theory to rigorously prove existence and uniqueness of instabilities
• Asymptotic Analysis: Asymptotic expansions in both small ε and large ε limits provide analytical insights

2. Advanced Numerical Methods

• Multi-Scale Treatment: Shooting method + continuation method effectively handles nonlinear boundary value problems
• Spectral Method Combination: FFT for periodic direction, finite difference near poles
• High Precision Verification: Numerical results agree with theoretical predictions to within 5% error

3. Clear Physical Picture

• Mode Classification: Complete analysis of three classes of modes m=0, 1, ≥2
• Topological Constraints: Clear exposition of physical meaning of Poincaré-Hopf theorem
• Visualization: 3D density plots intuitively display vortex pair formation process

4. Experimental Relevance

• Controllable Parameters: Single parameter ε contains all physical information (scattering length, atom number, thickness)
• Feasibility: Matches conditions of ISS and ground-based experiments
• Clear Predictions: Provides testable quantitative predictions (thresholds, vortex pair numbers)

5. Comprehensive Supplementary Material

• 19 pages of supplementary material providing complete derivations
• Includes implementation details of all numerical methods
• Facilitates reader understanding and reproduction

Shortcomings

1. Absence of Experimental Verification

• Pure theoretical study lacking experimental data comparison
• No discussion of specific schemes and challenges for experimental realization
• Experimental accessible range of parameter ε not clarified

2. Incomplete Long-Time Dynamics

• Only shows up to vortex pair formation stage (short time)
• Subsequent evolution of vortex pairs (annihilation, recombination) not deeply studied
• Final steady-state structure not determined

3. Finite Temperature Effects

• Zero-temperature approximation neglects thermal excitations
• Temperature effects potentially significant in actual experiments
• Relationship between energy and dynamical instabilities not discussed

4. Incomplete m=0 Mode Analysis

• Theory predicts possible instability but not observed numerically
• Lacks deep explanation of this discrepancy
• May require more refined numerical calculations or theoretical corrections

5. Gap from Experimental Conditions

• Ideal spherical assumption; actual systems may have ellipticity
• Neglects atomic loss and heating effects
• Effects of external perturbations (residual gravity) not considered

Impact Assessment

Academic Value

• Theoretical Completeness: Provides benchmark theory for soliton stability in spherical geometry BEC
• Methodological Contribution: Systematic application of BdG method in curved space
• Interdisciplinary: Connects nonlinear wave theory, topology, and cold atom physics

Practical Value

• Experimental Guidance: Provides clear parameter guidance for ISS and ground-based experiments
• Phenomenon Prediction: Quantitative prediction of vortex pair numbers directly testable
• Technical Reference: Numerical methods applicable to other curved geometry problems

Reproducibility

• High: Supplementary material provides complete algorithm details
• Standard numerical methods (FFT, Crank-Nicolson, shooting method)
• Clear parameter settings, easy to reproduce

Potential Impact

1. Short-term: Stimulates related experimental verification work
2. Medium-term: Promotes extended research on two-component systems, finite temperature effects
3. Long-term: Provides new platform for quantum simulation and topological physics

Applicable Scenarios

Direct Application

1. Microgravity BEC Experiments: ISS cold atom laboratory
2. Ground-Based Bubble BEC: Two-component mixture systems
3. Numerical Simulation: Theoretical research on spherical geometry BEC

Generalizable Scenarios

1. Other Closed Surfaces: Ellipsoids, tori, etc.
2. Two-Component Systems: Filled dark solitons, spin-orbit coupling
3. Other Quantum Fluids: Superfluid helium, exciton-polaritons

Limiting Conditions

• Two-dimensional approximation validity: δR/R < 0.1
• Zero-temperature approximation: T << Tc
• Mean-field validity: Na³ >> 1
• Spherical symmetry

Key References

1. Experimental Foundation:
   • Carollo et al., Nature 606, 281 (2022) - ISS bubble BEC observation
   • Jia et al., PRL 129, 243402 (2022) - Ground-based shell BEC
2. Theoretical Framework:
   • Geyer & Pelinovsky, AMS Mathematical Surveys 288 (2025) - Nonlinear wave stability in Hamiltonian systems
   • Andriati et al., PRA 104, 033318 (2021) - Mixture stability on bubble traps
3. Dark Soliton Physics:
   • Burger et al., PRL 83, 5198 (1999) - First observation of dark solitons in BEC
   • Anderson et al., PRL 86, 2926 (2001) - Dark soliton decay into vortex rings
4. Spherical BEC:
   • Tononi & Salasnich, PRL 123, 160403 (2019) - Spherical BEC theory
   • Dubessy & Perrin, AVS Quantum Sci. 7, 010501 (2025) - Review of quantum gases in bubble traps

Overall Assessment: This is an excellent paper with rigorous theory, advanced methods, and clear physical insights. Through perfect combination of analytical and numerical methods, it systematically solves the fundamental problem of dark soliton stability in spherical bubble BEC. The asymptotic formula εₘ ≈ 4m(m-1) is elegantly simple, and the correspondence between vortex pair numbers and dominant modes has clear experimental testability. The paper provides important theoretical guidance for ongoing microgravity BEC experiments and methodological reference for nonlinear wave research in curved spaces. Main shortcomings are absence of experimental verification and consideration of finite temperature effects, but these are natural directions for follow-up research.