Bright Fractional Single and Multi-Solitons in a Prototypical Nonlinear Schr{ö}dinger Paradigm: Existence, Stability and Dynamics

Titolo: Solitoni Singoli e Multipli Frazionari Brillanti in un Paradigma Prototipico dell'Equazione di Schrödinger Non Lineare: Esistenza, Stabilità e Dinamica

Autore: Robert J. Decker, A. Demirkaya, T.J. Alexander, P. G. Kevrekidis.

1. Problema e Contesto

Il lavoro esplora le proprietà delle onde solitarie (solitoni) e delle loro coppie in una variante frazionaria dell'equazione di Schrödinger non lineare (NLS) focalizzante.

• Contesto: Recentemente, esperimenti in ottica non lineare hanno permesso la realizzazione di esponenti frazionari arbitrari, permettendo l'interpolazione omotopica tra limiti armonici ($\alpha=2$) e biarmonici/quartici ($\alpha=4$).
• Obiettivo: Analizzare come le proprietà delle onde solitarie singole e multi-solitarie evolvono al variare continuo dell'esponente frazionario $\alpha$ (dove $0 < \alpha \le 4$), colmando il divario tra i modelli standard e quelli di ordine superiore.
• Equazione di riferimento: L'equazione NLS frazionaria focalizzante (FF-NLS) con non linearità cubica:
  $$i u_t = D_\alpha u + 2u|u|^2$$
  dove $D_\alpha$ è la derivata di Riesz frazionaria definita nello spazio di Fourier come $D_\alpha f = \mathcal{F}^{-1}[-|k|^\alpha \hat{f}(k)]$.

2. Metodologia

Gli autori adottano un approccio combinato di analisi numerica e spettrale:

• Metodi Numerici:
  • Discretizzazione della variabile spaziale su un intervallo $[-L, L]$ utilizzando la Trasformata di Fourier Veloce (FFT) per approssimare la derivata frazionaria di Riesz.
  • Utilizzo di fsolve (Matlab) per trovare stati stazionari, partendo da inizializzatori noti (es. $\text{sech}(x)$ per $\alpha=2$) e utilizzando tecniche di continuazione numerica per variare $\alpha$.
  • Simulazioni dirette delle PDE utilizzando routine ODE (ode45, ode23) per studiare la dinamica temporale.
• Analisi di Stabilità:
  • Linearizzazione dell'equazione attorno agli stati stazionari per ottenere un sistema di autovalori.
  • Calcolo dello spettro degli operatori lineari $L_+$ e $L_-$ per identificare autovalori reali (instabilità esponenziale) o immaginari (modi interni/oscillatori).
  • Classificazione delle instabilità in base al tipo di perturbazione (spostamento di posizione vs. rottura di simmetria).

3. Risultati Chiave

A. Solitoni Singoli

• Comportamento della Coda (Tail):
  • Per $0 < \alpha \le 2$: Le code sono monotone e decrescenti esponenzialmente (nessun attraversamento dello zero).
  • Per $2 < \alpha < 4$: Le code mostrano oscillazioni smorzate con attraversamenti dello zero. All'aumentare di $\alpha$ verso 4, il numero di attraversamenti tende all'infinito.
  • Per $\alpha \to 4$: La coda è dominata da un decadimento di tipo potenza ($a/x^{\alpha+1}$) con oscillazioni sovrapposte.
• Stabilità e Collasso:
  • È identificata una soglia critica a $\alpha = 1$ (dove $d=1$ è la dimensionalità).
  • $\alpha < 1$: Instabilità. Esiste una coppia di autovalori reali ($\lambda = \pm \lambda_r$) che porta a dinamiche di collasso (concentrazione infinita) o dispersione, a seconda della fase della perturbazione.
  • $\alpha > 1$: Stabilità spettrale. Gli autovalori reali si spostano sull'asse immaginario, trasformando l'instabilità in un modo interno oscillatorio (respirazione del solitone).

B. Coppie di Solitoni (Solitoni Multipli)

Lo studio si concentra su coppie di solitoni in fase e in controfase.

• Esistenza:
  • Per $\alpha \le 2$: Non esistono stati stazionari per coppie di solitoni (nessuna soluzione legata).
  • Per $\alpha > 2$: Emergono stati stazionari multipli ("molecole solitoniche"). All'aumentare di $\alpha$ verso 4, il numero di rami di soluzioni aumenta, diventando infinito per $\alpha=4$.
• Stabilità e Biforcazioni:
  • Caso $\alpha = 4$ (Biarmonico): Confermato che tutte le coppie di solitoni sono instabili. Le coppie dispari sono instabili per traslazione (movimento relativo), le pari per rottura di simmetria (variazione di ampiezza).
  • Caso Frazionario ($2 < \alpha < 4$): Scoperta fondamentale: esistono finestre di stabilità per certi rami di soluzioni (in particolare per le coppie in fase con numero pari di solitoni).
  • Meccanismo di Stabilizzazione: Attraverso l'analisi spettrale, si osserva che gli autovalori reali (instabilità) si spostano sull'asse immaginario per certi valori di $\alpha$, stabilizzando il sistema contro la rottura di simmetria. Questo fenomeno non si verifica nei casi interi ($\alpha=2, 4$).
• Tipi di Instabilità:
  • Tipo I: Instabilità di rottura di simmetria (variazione relativa delle ampiezze).
  • Tipo II: Instabilità di moto relativo (variazione delle posizioni).
  • Per i solitoni in fase, i rami pari possono diventare stabili (nessun autovalore reale) in un intervallo frazionario. Per i solitoni in controfase, la situazione è "ribaltata": i rami pari sono instabili per traslazione e non mostrano stabilizzazione completa.

4. Contributi e Significatività

1. Connessione tra Limiti: Il lavoro fornisce un ponte continuo tra il limite armonico ($\alpha=2$) e quello biarmonico ($\alpha=4$), mostrando come le proprietà delle onde solitarie evolvono al variare continuo della dispersione.
2. Scoperta di Stabilizzazione Frazionaria: Il risultato più rilevante è l'identificazione di finestre di stabilità per coppie di solitoni in fase nel regime frazionario ($2 < \alpha < 4$), un fenomeno assente nei modelli con esponenti interi. Questo suggerisce che la frazionalità può essere utilizzata come parametro di controllo per stabilizzare strutture multi-pulse.
3. Dinamica di Collasso: Viene chiarito il legame tra l'instabilità dei solitoni singoli per $\alpha < 1$ e i fenomeni di collasso, monitorando lo spettro di linearizzazione.
4. Implicazioni Sperimentali: Dato che esperimenti recenti in ottica non lineare hanno già realizzato dispersioni frazionarie, i risultati teorici offrono previsioni concrete per l'osservazione di "molecole solitoniche" stabili in laboratorio, specialmente nel regime frazionario intermedio.
5. Classificazione delle Instabilità: Viene introdotta una distinzione sistematica tra instabilità di tipo I (ampiezza) e Tipo II (posizione), mostrando come queste si scambino o si annullino a seconda della fase relativa dei solitoni e del valore di $\alpha$.

Conclusione

Lo studio dimostra che la variazione dell'esponente frazionario $\alpha$ non è solo una generalizzazione matematica, ma introduce nuove fenomenologie fisiche, in particolare la possibilità di stabilizzare strutture solitoniche complesse che sarebbero instabili nei limiti classici. Questo apre nuove strade per il controllo della dinamica delle onde non lineari in sistemi con dispersione frazionaria.