Global bifurcations and basin geometry of the nonlinear non-Hermitian skin effect

1. Problema e Contesto

Il lavoro si inserisce nel campo della fisica non-ermitiana, in particolare nello studio dell'effetto pelle non-ermitiano (NHSE), un fenomeno in cui un numero esteso di stati di bulk si accumulano ai bordi del sistema sotto condizioni al contorno aperte. Mentre il NHSE lineare è ben compreso attraverso la teoria delle bande non-Bloch, la sua interazione con la nonlinearità rimane un fronte di ricerca emergente.

Il problema specifico affrontato dagli autori riguarda la dinamica di un modello di Hatano-Nelson continuo con una nonreciprocità nonlineare saturante. A differenza dei modelli precedenti che mostravano una transizione guidata da una biforcazione di Hopf supercritica (che porta a un singolo stato stabile), questo studio indaga cosa accade quando la nonreciprocità dipende dall'ampiezza in modo non monotono e saturante (termine cubico-quintico). Le domande chiave sono:

• Quali cambiamenti qualitativi avvengono se il guadagno efficace è soppresso ad alte ampiezze?
• Quali altre biforcazioni dinamiche (oltre a quella di Hopf supercritica) possono verificarsi?
• È possibile la coesistenza di stati stazionari multipli (bistabilità) sotto gli stessi parametri, portando a stati dipendenti dalla storia e dalla preparazione?

2. Metodologia

Gli autori adottano un approccio basato sulla dinamica dei sistemi per analizzare le equazioni stazionarie del modello.

• Modello: Viene studiato un'equazione di Schrödinger stazionaria nonlineare con un Hamiltoniano non-ermitiano $\hat{H}(\psi) = \frac{\hat{k}^2}{2} + i F(|\psi|^2) \hat{k}$, dove il termine non-hermitiano $F(z) = \gamma + az - bz^2$ (con $z=|\psi|^2$) rappresenta una modulazione saturante.
• Formulazione nello spazio delle fasi: L'equazione differenziale stazionaria viene riscritta come un sistema dinamico planare bidimensionale, trattando la coordinata spaziale $x$ come un parametro di evoluzione temporale. Lo stato stazionario corrisponde a una traiettoria $(\psi(x), v(x))$ nello spazio delle fasi.
• Analisi delle Biforcazioni: Viene costruita una mappa di biforcazione globale analizzando la stabilità dei punti fissi e dei cicli limite al variare del parametro di nonreciprocità $\gamma$.
• Equazione di Ampiezza Media: Per ottenere previsioni analitiche, viene derivata un'equazione di ampiezza media (averaged amplitude equation) utilizzando il metodo di media, assumendo che la modulazione non-hermitiana sia piccola rispetto alla frequenza naturale del sistema.
• Geometria dei Bacini di Attrazione: Viene introdotto un parametro d'ordine basato sulla frazione del bacino di attrazione per quantificare la probabilità di trovare lo stato "pelle" o "esteso" in funzione delle condizioni iniziali.

3. Contributi Chiave e Risultati

A. Scenario di Biforcazione Globale

Il lavoro rivela uno scenario di biforcazione più ricco rispetto ai modelli lineari o a quelli con Hopf supercritico:

1. Biforcazione di Hopf Subcritica: All'origine ($\gamma = 0$), il punto fisso subisce una biforcazione di Hopf subcritica. Per $\gamma < 0$ (ma vicino a 0), esiste un ciclo limite instabile di piccola ampiezza.
2. Biforcazione Nodo-Sella di Cicli Limite (SNLC): A un valore critico $\gamma = \gamma_c < 0$, avviene una biforcazione SNLC che crea una coppia di cicli limite: uno stabile (esterno) e uno instabile (interno).
3. Finestra di Coesistenza: Tra $\gamma_c$ e $0$, il sistema si trova in una regione di bistabilità. In questa finestra, coesistono due attrattori stabili:
   • L'origine (corrispondente a uno stato skin-localizzato).
   • Il ciclo limite stabile esterno (corrispondente a uno stato esteso).

B. Meccanismo di Separazione Non Lineare

A differenza dei sistemi lineari dove la coesistenza di stati localizzati ed estesi è legata a strutture spettrali (come i bordi di mobilità), qui gli stati coesistono alla stessa energia $E$. Sono separati da una separatrice non lineare nello spazio delle fasi, organizzata dal ciclo limite instabile interno. La condizione iniziale (la pendenza al bordo $\partial_x \psi(0)$) determina quale attrattore verrà raggiunto.

C. Predizioni Analitiche

L'equazione di ampiezza media derivata fornisce previsioni in forma chiusa per:

• Le ampiezze dei rami dei cicli limite ($A_{in}$ e $A_{out}$).
• La soglia della biforcazione SNLC: $\gamma_c = -a^2 / (8b)$.
  Queste previsioni mostrano un eccellente accordo con le simulazioni numeriche.

D. Parametro d'Ordine e Diagramma di Fase

Gli autori introducono un parametro d'ordine basato sulla frazione del bacino ($p_{skin}$), che misura la probabilità che una traiettoria casuale finisca nello stato pelle.

• Alla soglia $\gamma_c$, il parametro d'ordine mostra un salto di tipo primo ordine (discontinuità), poiché la creazione del ciclo limite stabile esterno cattura immediatamente una frazione finita di condizioni iniziali.
• Questo trasforma la finestra di bistabilità in un vero e proprio diagramma di fase geometrico.

E. Fenomeni Fisici Osservabili

Nella finestra di coesistenza, vengono rivelati fenomeni fisici distintivi:

1. Transitori Spaziali a Lunga Vita: Traiettorie inizializzate vicino alla separatrice possono "oscurare" (shadow) il ciclo limite instabile per lunghe distanze spaziali prima di collassare su uno degli attrattori stabili. Questo crea stati transitori che sembrano estesi ma decadono o crescono solo dopo un lungo periodo, un fenomeno assente nei sistemi lineari.
2. Isteresi: A causa della bistabilità, il sistema mostra un effetto di memoria. Se il parametro $\gamma$ viene variato adiabaticamente:
   • Scendendo da $\gamma > 0$, il sistema rimane nello stato esteso fino a $\gamma_c$, dove salta allo stato pelle.
   • Salendo da $\gamma < \gamma_c$, il sistema rimane nello stato pelle fino a $\gamma = 0$, dove salta allo stato esteso.
     Questo crea un ciclo di isteresi tra $\gamma_c$ e $0$.

4. Significato e Implicazioni

Questo studio è significativo perché:

• Supera i concetti spettrali lineari: Dimostra che nella fisica non-ermitiana nonlineare, la geometria globale degli attrattori e dei loro bacini di attrazione è uno strumento fondamentale, complementare e talvolta superiore ai concetti spettrali tradizionali (come gli autovalori o i bordi di mobilità) per comprendere gli stati stazionari.
• Nuovi meccanismi di localizzazione: Identifica un meccanismo di coesistenza stato-localizzato/stato-esteso guidato dalla dinamica nonlineare e dalla geometria dei bacini, anziché dalla struttura dello spettro.
• Applicabilità: I risultati hanno implicazioni per la progettazione di dispositivi fotonici non lineari, laser e sistemi topologici, dove il controllo della bistabilità, dell'isteresi e dei transitori spaziali può essere sfruttato per nuove funzionalità di memoria e commutazione.

In sintesi, il lavoro fornisce una mappa completa della dinamica non-hermitiana nonlineare, evidenziando come la saturazione della nonreciprocità porti a una ricca fenomenologia di biforcazioni globali e geometrie di bacini complessi.