Self-Organized Stabilization of Straight Dark Solitons in Stripe Supersolids

Problema
Nei fluidi quantistici bidimensionali (2D), i solitoni scuri rettilinei sono genericamente instabili a causa dell'instabilità modulazionale trasversale, che ne guida il decadimento in anelli di vortici (in 3D) o in vortici solitici e coppie vortice-antivortice (in 2D). Sebbene le precedenti strategie di stabilizzazione si siano basate su meccanismi esterni come il confinamento trasversale stretto, il pinning tramite reticolo ottico o l'accoppiamento spin-orbita ingegnerizzato, rimane aperta la questione se un fluido quantistico genuinamente 2D possa auto-organizzare un paesaggio stabilizzante per un difetto di fase lineare attraverso sole interazioni intrinseche. In particolare, non è chiaro se l'ordine posseduto dal mezzo stesso possa proteggere il difetto senza pinning esterno.

Metodologia
Gli autori impiegano simulazioni teoriche utilizzando un'equazione di Gross–Pitaevskii estesa (eGPE) che incorpora interazioni di contatto, interazioni dipolo-dipolo non locali e la correzione di Lee–Huang–Yang (LHY) nell'approssimazione della densità locale. Modellano un gas di Bose dipolare quasi-2D (utilizzando specificamente i parametri per $^{166}\text{Er}$ e $^{164}\text{Dy}$) sotto un forte confinamento armonico lungo l'asse $\hat{z}$, con i dipoli polarizzati nel piano $x$–$z$ ad un angolo $\alpha$.

Lo studio utilizza due strumenti analitici primari:

1. Propagazione nel tempo immaginario: Per ottenere stati fondamentali di solitone rettilineo stazionari ($\psi_0$) imponendo una linea nodale in $y=0$ e un salto di fase di $\pi$ attraverso l'asse del solitone.
2. Analisi di stabilità lineare di Bogoliubov–de Gennes (BdG): Per calcolare lo spettro di eccitazione attorno allo stato stazionario, identificando frequenze complesse $\Omega = \Omega_R + i\Omega_I$. L'instabilità dinamica è segnalata da $\Omega_I > 0$, mentre l'instabilità energetica è identificata da modi con $\Omega_R < 0$ e norma di Bogoliubov positiva.
3. Propagazione nel tempo reale: Per simulare l'evoluzione dinamica di stati perturbati e osservare scenari di decadimento non lineare.
4. Modellazione analitica: Per derivare il costo energetico della curvatura trasversale del solitone ($\Delta E$) ed estrarre il coefficiente di rigidità elastica della curvatura ($\sigma_2$) in funzione della lunghezza di scattering ($a_s$) e dell'angolo di polarizzazione ($\alpha$).

Contributi chiave e risultati
Il documento dimostra che le interazioni a lungo raggio anisotrope in un condensato dipolare quasi-2D possono stabilizzare un solitone scuro rettilineo incorporato contro il decadimento trasversale attraverso due distinti meccanismi auto-organizzati:

1. Indurimento spettrale via anisotropia: Vicino alla polarizzazione in piano ($\alpha \approx \pi/2$), il campo dipolare anisotropo rende più nitida la parete del dominio del solitone. Questa interazione da sola indurisce il ramo solitonico trasversale più basso, mantenendo la sua frequenza reale ($\Omega_I = 0$) e prevenendo l'ammorbidimento che porta all'instabilità, anche prima che si sviluppi un contrasto di densità a strisce finito.
2. Pinning auto-organizzato tramite ordine a strisce: All'aumentare della forza dell'interazione dipolare relativa (diminuendo $a_s$), il sistema subisce una transizione verso un supersolido (SS) dipolare con ordine a strisce spontaneo. Questa modulazione di densità emergente agisce come un "paesaggio di pinning auto-generato". Essa indurisce ulteriormente il ramo solitonico trasversale, allargando il gap spettrale e aumentando significativamente il coefficiente di rigidità elastica della curvatura ($\sigma_2$). Ciò impone un elevato costo energetico statico alle deformazioni trasversali, penalizzando efficacementamente la curvatura del solitone.

Lo studio quantifica la soglia di instabilità, riscontrando che i solitoni rettilinei rimangono stabili per $a_s \lesssim 55.96$ (a $\alpha = \pi/2$). Oltre questa soglia, il sistema esibisce due regimi di decadimento distinti a seconda della prossimità all'instabilità:

• Vicino alla soglia ($a_s \approx 60$): Un singolo modo trasversale instabile ($n_T=1$) porta alla nucleazione di un singolo vortice solitonico.
• Instabilità profonda ($a_s \gtrsim 62$): Molteplici modi trasversali ($n_T > 1$) diventano instabili simultaneamente, portando all'instabilità "snake" (a serpente) e alla nucleazione di molteplici coppie vortice–antivortice.

La stabilità è altamente sensibile all'angolo di polarizzazione $\alpha$. Una stabilità robusta persiste solo per $\alpha \gtrsim 80^\circ$, con la massima stabilizzazione a $\alpha = \pi/2$. Inclinare i dipoli indebolisce l'anisotropia in piano, sopprime la modulazione a strisce e ripristina il background superfluido liscio, rendendo il solitone suscettibile all'instabilità trasversale.

Significato
Gli autori affermano che questo lavoro stabilisce un meccanismo di protezione intrinseco, guidato dalle interazioni, per i solitoni scuri rettilinei nei fluidi quantistici strutturati. A differenza dei precedenti metodi di stabilizzazione che si affidano a un confinamento esterno o al disordine, questo meccanismo deriva dal mezzo stesso: le stesse interazioni microscopiche che generano l'ordine supersolido a strisce creano anche il paesaggio di pinning che stabilizza il difetto.

I risultati suggeriscono che i solitoni scuri rettilinei nei supersolidi dipolari potrebbero servire come robusti elementi sensibili alla fase per circuiti atomtronici in dimensioni superiori, estendendo la loro utilità oltre le geometrie quasi-1D. Inoltre, la risposta di un solitone scuro in movimento al background periodico di densità offre un potenziale protocollo per valutare sperimentalmente la rigidità dello stato supersolido. I parametri richiesti sono accessibili nelle attuali piattaforme sperimentali utilizzando $^{166}\text{Er}$ e $^{164}\text{Dy}$.