Self-Organized Stabilization of Straight Dark Solitons in Stripe Supersolids

Problem
In two-dimensional (2D) quantum fluids, straight dark solitons are generically unstable due to transverse modulational instability, which drives their decay into vortex rings (in 3D) or solitonic vortices and vortex–antivortex pairs (in 2D). While previous stabilization strategies have relied on external mechanisms such as tight transverse confinement, optical-lattice pinning, or engineered spin-orbit coupling, it remains an open question whether a genuinely 2D quantum fluid can self-organize a stabilizing landscape for a line phase defect through intrinsic interactions alone. Specifically, it is unclear if the order possessed by the medium itself can protect the defect without external pinning.

Methodology
The authors employ theoretical simulations using an extended Gross–Pitaevskii equation (eGPE) that incorporates contact interactions, nonlocal dipole–dipole interactions, and the Lee–Huang–Yang (LHY) correction in the local-density approximation. They model a quasi-2D dipolar Bose gas (specifically using parameters for $^{166}\text{Er}$ and $^{164}\text{Dy}$) under tight harmonic confinement along the $\hat{z}$ axis, with dipoles polarized in the $x$–$z$ plane at an angle $\alpha$.

The study utilizes two primary analytical tools:

1. Imaginary-time propagation: To obtain stationary straight-soliton ground states ($\psi_0$) by enforcing a nodal line at $y=0$ and a $\pi$ phase jump across the soliton axis.
2. Bogoliubov–de Gennes (BdG) linear stability analysis: To calculate the excitation spectrum around the stationary state, identifying complex frequencies $\Omega = \Omega_R + i\Omega_I$. Dynamical instability is signaled by $\Omega_I > 0$, while energetic instability is identified by modes with $\Omega_R < 0$ and positive Bogoliubov norm.
3. Real-time propagation: To simulate the dynamical evolution of perturbed states and observe nonlinear decay scenarios.
4. Analytical modeling: To derive the energy cost of transverse soliton bending ($\Delta E$) and extract the elastic bending stiffness coefficient ($\sigma_2$) as a function of the scattering length ($a_s$) and polarization angle ($\alpha$).

Key Contributions and Results
The paper demonstrates that anisotropic long-range interactions in a quasi-2D dipolar condensate can stabilize an embedded straight dark soliton against transverse decay through two distinct, self-organized mechanisms:

1. Spectral Hardening via Anisotropy: Near in-plane polarization ($\alpha \approx \pi/2$), the anisotropic dipolar field sharpens the soliton domain wall. This interaction alone hardens the lowest transverse solitonic branch, keeping its frequency real ($\Omega_I = 0$) and preventing the softening that leads to instability, even before a finite stripe density contrast develops.
2. Self-Organized Pinning via Stripe Order: As the relative dipolar interaction strength increases (decreasing $a_s$), the system undergoes a transition to a dipolar supersolid (SS) with spontaneous stripe order. This emergent density modulation acts as a "self-generated pinning landscape." It further hardens the transverse solitonic branch, widening the spectral gap, and significantly increases the elastic bending stiffness coefficient ($\sigma_2$). This imposes a large static energy cost on transverse deformations, effectively penalizing the bending of the soliton.

The study quantifies the instability threshold, finding that straight solitons remain stable for $a_s \lesssim 55.96$ (at $\alpha = \pi/2$). Beyond this threshold, the system exhibits two distinct decay regimes depending on the proximity to the instability:

• Near-threshold ($a_s \approx 60$): A single unstable transverse mode ($n_T=1$) leads to the nucleation of a single solitonic vortex.
• Deep instability ($a_s \gtrsim 62$): Multiple transverse modes ($n_T > 1$) become unstable simultaneously, leading to the "snake instability" and the nucleation of multiple vortex–antivortex pairs.

The stability is highly sensitive to the polarization angle $\alpha$. Robust stability persists only for $\alpha \gtrsim 80^\circ$, with maximum stabilization at $\alpha = \pi/2$. Tilting the dipoles weakens the in-plane anisotropy, suppresses stripe modulation, and restores the smooth superfluid background, rendering the soliton susceptible to transverse instability.

Significance
The authors claim this work establishes an intrinsic, interaction-driven protection mechanism for straight dark solitons in structured quantum fluids. Unlike previous stabilization methods that rely on external confinement or disorder, this mechanism arises from the medium itself: the same microscopic interactions that generate the stripe supersolid order also create the pinning landscape that stabilizes the defect.

The results suggest that straight dark solitons in dipolar supersolids could serve as robust phase-sensitive elements for atomtronic circuits in higher dimensions, extending their utility beyond quasi-1D geometries. Furthermore, the response of a moving dark soliton to the periodic density background offers a potential protocol for experimentally assessing the rigidity of the supersolid state. The required parameters are accessible in current experimental platforms using $^{166}\text{Er}$ and $^{164}\text{Dy}$.