Stability of dark solitons in a bubble Bose-Einstein condensate

This paper establishes the stability criteria for dark solitons on a spherical Bose-Einstein condensate, demonstrating that beyond a sharp nonlinear threshold, they decay into vortex pairs via a universal single-mode snake instability unique to the two-dimensional surface geometry.

The Gist

Imagine a giant, invisible soap bubble floating in space. But instead of soap and water, this bubble is made of a super-cold gas where all the atoms act like a single, giant wave. Scientists call this a Bose-Einstein Condensate (BEC).

Now, imagine you take a tiny "dent" or a "hole" in the density of this gas wave, creating a dark ring that circles the bubble's equator. In physics, this is called a dark soliton. It's like a shadow moving across a bright light.

This paper asks a simple but tricky question: Is this shadow stable, or will it eventually break apart?

Here is the story of what happens, explained without the heavy math.

1. The Setting: A Bubble in Space

Usually, we think of these gas bubbles as being on Earth, where gravity pulls them down into a flat pancake shape. But recently, scientists have been able to create these bubbles in space (on the International Space Station), where there is no gravity. Without gravity, the gas forms a perfect, hollow sphere—a true 3D bubble.

The researchers wanted to know: If you make a dark ring (a soliton) on this spherical bubble, will it stay a perfect ring, or will it wiggle and break?

2. The "Snake" Problem

On a flat surface (like a pond), if you make a straight line of water that is lower than the rest, it doesn't stay straight for long. It starts to wiggle side-to-side, like a snake slithering. This is called a "snake instability."

Eventually, that wiggling line breaks apart into little whirlpools (vortices).

• On a flat table: These whirlpools can fly off the edge and disappear.
• On a sphere: There are no edges! The surface curves back on itself. If a whirlpool forms, it can't escape. It has to stay on the bubble.

3. The Discovery: The "Breaking Point"

The authors found that the stability of this dark ring depends on how "strong" the gas is (how much the atoms push against each other).

• The Safe Zone: If the gas is "weak" (low interaction), the dark ring is happy. It stays a perfect circle, wobbling slightly but never breaking.
• The Danger Zone: If you increase the interaction strength past a specific "tipping point," the ring becomes unstable. It starts to snake violently.

4. The Magic Number: How Many Vortices?

Here is the most fascinating part. When the ring breaks, it doesn't just turn into a random mess. It follows a strict rule based on the shape of the sphere.

The researchers found that the ring breaks into pairs of whirlpools (one spinning clockwise, one counter-clockwise). Why pairs? Because on a sphere, you can't have just one lonely whirlpool; the geometry demands they come in couples to balance out.

The number of pairs depends on how the ring wiggles:

• If it wiggles in a simple "figure-8" pattern, it breaks into 2 pairs of whirlpools.
• If it wiggles in a "flower" pattern with 3 petals, it breaks into 3 pairs.
• If it has 4 petals, it breaks into 4 pairs.

It's like a cookie cutter. The shape of the instability (the number of "petals" in the wiggle) dictates exactly how many pairs of whirlpools will be born.

5. The Big Difference: Rings vs. Pairs

In previous experiments with 3D gas clouds (not hollow bubbles), when a dark soliton broke, it formed vortex rings (like smoke rings).
But on this 2D bubble surface, the geometry is different. The "smoke rings" can't form because the gas is trapped on a thin skin. Instead, the energy forces the creation of vortex pairs stuck to the surface.

The Takeaway

This paper is like a rulebook for a cosmic game of "break the ring."

1. The Game: You have a dark ring on a spherical bubble.
2. The Rule: If you push the gas too hard, the ring snaps.
3. The Result: It snaps into a specific number of whirlpool pairs, determined by the pattern of the snap.

This is important because it helps scientists predict what will happen in future experiments on the International Space Station. If they see a dark ring wiggling, they can now predict exactly how many whirlpools will appear when it breaks, turning a chaotic explosion into a predictable, beautiful pattern.

In short: A dark ring on a space bubble is like a tightrope walker. If they stay calm, they walk fine. If they get too energetic, they trip and fall, but they always fall in a specific, symmetrical pattern of pairs, never alone.

Technical Summary

1. Problem Statement

The paper investigates the spectral and dynamical stability of dark solitons confined to the surface of a spherical Bose-Einstein condensate (BEC), often referred to as a "bubble" BEC. This system is motivated by recent experiments on the International Space Station (ISS) involving ultra-cold gases in microgravity, which allow for the creation of shell-like potentials.

While dark solitons in flat, quasi-1D or quasi-2D BECs are known to undergo snake instabilities (decaying into vortex-antivortex pairs), the topological constraints of a closed spherical surface introduce unique geometric effects. Specifically, the paper addresses:

• How the curvature of a sphere alters the stability dynamics compared to flat geometries.
• The specific decay pathways of dark solitons on a sphere (e.g., do they form vortex rings as in 3D, or something else?).
• The existence of a sharp instability threshold governed by the nonlinear interaction parameter.

2. Methodology

The authors employ a combination of analytical perturbation theory and high-precision numerical simulations.

Mathematical Model:

• Dimensional Reduction: The system is modeled by reducing the 3D Gross-Pitaevskii Equation (GPE) to a 2D system on a rigid spherical shell of radius $R$ and thickness $\delta R$ (where $\delta R \ll R$).
• Governing Equation: The dimensionless 2D GPE is given by:
  $$i\partial_t \psi = -\Delta_{2D}\psi + g|\psi|^2\psi$$
  where $\Delta_{2D}$ is the Laplacian on the sphere, and $g$ is the nonlinear interaction parameter.
• Dark Soliton Solution: The stationary dark soliton is defined as a density depression at the equator ($\theta = \pi/2$), represented by $\psi_s(\theta, \phi, t) = f(\theta)e^{-i\mu t}$. The profile $f(\theta)$ is solved analytically in two limits:
  • Small $\epsilon$ (weak nonlinearity): Perturbative expansion around the linear mode $P_1(\cos \theta)$.
  • Large $\epsilon$ (strong nonlinearity): Asymptotic expansion matching the flat-space tanh solution near the equator.

Stability Analysis:

• Bogoliubov-de Gennes (BdG) Approach: The authors linearize the GPE around the stationary solution using perturbations $u$ and $v$ with angular momentum modes $m$.
• Spectral Problem: This leads to a coupled eigenvalue problem involving operators $L_m^-$ and $L_m^+$. The stability is determined by the eigenvalues $\omega$ of this system.
  • If $\omega$ is real, the mode is stable.
  • If $\omega$ has a non-zero imaginary part, the mode is unstable (exponential growth).
• Analytical Proofs: The authors use Sturm-Liouville theory and operator comparison theorems ($L_m^+ > L_m^-$ and $L_{m+1} > L_m$) to prove the existence of unstable modes.

Numerical Simulations:

• Stationary States: Solved using the shooting method combined with the secant method.
• Time Evolution: The full time-dependent GPE is solved using a split-step spectral method combined with finite differences (Crank-Nicolson) in the polar direction and Fast Fourier Transforms (FFT) in the azimuthal direction.

3. Key Contributions and Results

A. Existence and Profiles

The study confirms the existence of dark ring solitons on a sphere. The density profile $f(\theta)$ is monotonic, vanishing at the equator.

• As the nonlinearity parameter $\epsilon$ increases, the soliton width narrows, concentrating near the equator.
• The chemical potential $\mu$ scales linearly with $\epsilon$ in the strong interaction limit.

B. Spectral Stability and Instability Thresholds

The central finding is the identification of a sharp instability threshold dependent on the angular momentum mode $m$ and the nonlinearity $\epsilon$.

• Stable Regime ($m=0, 1$): The soliton is spectrally stable for $m=0$ and $m=1$ for all $\epsilon$.
• Unstable Regime ($m \geq 2$): For every angular mode $m \geq 2$, there exists a critical threshold $\epsilon_m$ above which the soliton becomes unstable.
• Threshold Formula: The authors derive an asymptotic formula for the threshold:
  $$\epsilon_m^{th} \approx 4m(m-1)$$
  This analytical prediction matches numerical results with high precision (error < 5% even for $m=2$).
• Mechanism: The instability is driven by a single unstable mode for each $m$. As $\epsilon$ increases, the lowest eigenvalue of the operator $L_m^-$ crosses zero, becoming negative, which triggers a complex conjugate pair of eigenvalues $\omega$ (indicating exponential growth).

C. Dynamical Decay and Vortex Formation

Numerical time-evolution simulations reveal the physical consequence of the instability:

• Snake Instability: The dark soliton develops sinusoidal perturbations (snake instability) along its length.
• Vortex Dipoles: Unlike 3D spherical BECs where dark solitons decay into vortex rings, the 2D bubble geometry forces the decay into vortex-antivortex pairs (dipoles).
• Topology Constraint: Due to the Poincaré-Hopf theorem, a single vortex cannot exist on a sphere; vortices must appear in pairs with opposite circulation to sum to zero net charge.
• Mode-Dependent Breakup: If the dominant unstable mode is $m$, the soliton breaks up into exactly $m$ vortex-antivortex pairs.
  • $m=2 \rightarrow$ 2 pairs.
  • $m=3 \rightarrow$ 3 pairs.
  • $m=4 \rightarrow$ 4 pairs.

4. Significance

• Universal Mechanism: The paper establishes a universal rule governing the decay of dark solitons on curved 2D surfaces: the number of resulting vortex pairs is strictly determined by the dominant angular momentum mode of the instability.
• Distinction from 3D: It clarifies a fundamental difference between 3D bulk BECs (vortex rings) and 2D bubble BECs (vortex dipoles), highlighting the role of dimensionality and topology in quantum fluid dynamics.
• Experimental Relevance: The results provide specific, testable predictions for ongoing experiments on the ISS and ground-based bubble traps. Researchers can control the number of vortex pairs generated by tuning the interaction strength (via Feshbach resonance) to exceed specific thresholds $\epsilon_m$.
• Theoretical Framework: The rigorous analytical proof of the instability threshold using operator theory provides a robust foundation for understanding nonlinear wave stability on manifolds.

Conclusion

The study demonstrates that dark solitons on a spherical bubble are stable only below a critical nonlinearity. Beyond this threshold, they undergo a snake instability dictated by a single unstable mode, leading to a deterministic decay into $m$ vortex-antivortex pairs. This work bridges the gap between theoretical stability analysis and experimental observations of quantum gases in curved, low-dimensional geometries.