PT-symmetry breaking phase transitions in an LMG dimer

Sintesi Tecnica: Transizioni di Fase della Rottura della Simmetria PT in un Dimero LMG

Problematica
I sistemi quantistici aperti guidati sia da campi esterni che dalla dissipazione esibiscono una fenomenologia più ricca rispetto ai loro corrispettivi in equilibrio, inclusa la stabilizzazione di fasi impossibili negli stati fondamentali. Mentre la ricerca precedente ha esplorato le transizioni di fase dissipative in modelli come i modelli di Dicke e Lipkin-Meshkov-Glick (LMG) dissipativi, e ha investigato separatamente le transizioni di rottura della simmetria Parità-Tempo (PT), l'interazione tra queste distinte classi di transizioni fuori equilibrio rimane poco esplorata. Nello specifico, vi è la necessità di comprendere come un sistema che combina un meccanismo di guida-dissipazione PT-simmetrico con un'interazione LMG non lineare e di rottura della simmetria si comporta. Gli autori mirano a costruire un modello che integri questi due meccanismi per investigare la risultante competizione, la struttura di fase e i regimi dinamici emergenti.

Metodologia
Gli autori introducono un modello costituito da una coppia di modelli LMG accoppiati: uno dissipativo e uno guidato, legati tramite un'Hamiltoniana PT-simmetrica. Il sistema è descritto da un equazione master che governa la matrice di densità $\rho$ di due spin collettivi (etichettati A e B) di lunghezza $S$.

• Costruzione del Modello: L'Hamiltoniana totale $H$ combina gli Hamiltoniani LMG per ciascuno spin ($H_A^{LMG} \otimes I_B + I_A \otimes H_B^{LMG}$) con un termine di accoppiamento PT-simmetrico $H_{PT}$. La dinamica include superoperatori di Lindblad che rappresentano il guadagno su uno spin e la perdita sull'altro.
• Analisi Mean-Field: Per analizzare il sistema nel limite termodinamico ($S \to \infty$), gli autori derivano un insieme di sei equazioni differenziali non lineari accoppiate per le componenti dello spin normalizzate ($X, Y, Z$). Queste equazioni sono esatte nel limite termodinamico per questa classe di modelli.
• Investigazione Numerica: Gli autori risolvono numericamente le equazioni mean-field per mappare il diagramma di fase dello stato stazionario, identificando punti fissi, cicli limite e regioni caotiche. La stabilità è determinata tramite gli autovalori della matrice Jacobiana.
• Caratterizzazione Dinamica: Per distinguere tra dinamiche regolari e caotiche, gli autori calcolano l'esponente di Lyapunov massimo mediato su condizioni iniziali simmetriche.
• Simulazione Quantistica Esatta: Per validare le previsioni mean-field, gli autori eseguono la diagonalizzazione numerica esatta dell'equazione master per dimensioni di spin finite ($S$), lavorando in uno spazio di Liouville di dimensione $(2S+1)^4$. Calcolano inoltre le funzioni di Wigner dello spin per confrontare le rappresentazioni nello spazio delle fasi dello stato quantistico con le traiettorie mean-field.

Contributi Chiave e Risultati

1. Diagramma di Fase Complesso: La competizione tra la transizione PT-simmetrica (guidata dall'accoppiamento inter-spin $\lambda$) e la transizione ferromagnetica LMG (guidata dalla non linearità $J$) genera un diagramma di fase con molteplici fasi di stato stazionario non presenti nei modelli costituenti.

   • Fase Normale: A bassi valori di $\lambda$ e $|J|$, il sistema esibisce un singolo punto fisso stabile corrispondente a una magnetizzazione sfalsata (dominata dalla dissipazione).
   • Fasi LMG: All'aumentare di $|J|$, il sistema entra in regioni con due punti fissi stabili, corrispondenti a stati ferromagnetici con rottura di simmetria (ordinamento spin-A o spin-B).
   • Regime PT/Caotico: A grandi $\lambda$, il sistema perde punti fissi stabili. La teoria mean-field predice l'emergere di cicli limite e, in specifiche regioni dei parametri, dinamiche caotiche caratterizzate da esponenti di Lyapunov massimi positivi.

2. Regimi Dinamici Emergenti: Il documento identifica che l'interazione delle due transizioni crea comportamenti dinamici assenti nei singoli modelli. In particolare, esistono regioni in cui l'analisi mean-field predice punti fissi stabili, ma l'esponente di Lyapunov è positivo, indicando dinamiche oscillatorie complesse dipendenti dalle condizioni iniziali (bacini di attrazione).

3. Validazione tramite Scaling delle Dimensioni Finite: Gli autori dimostrano che le caratteristiche predette dalla teoria mean-field sopravvivono nel modello quantistico completo.

   • I risultati della diagonalizzazione esatta per $S$ finito corrispondono qualitativamente ai confini di fase mean-field.
   • Gli effetti di dimensione finita smussano le transizioni nette (ad esempio, la transizione PT a $\lambda = g/2$) e garantiscono un operatore di densità stazionario unico anche quando la teoria mean-field predice cicli limite.
   • Le funzioni di Wigner dello spin confermano che le distribuzioni di probabilità quantistica si allineano con i punti fissi mean-field (nella fase LMG) o con distribuzioni diffuse corrispondenti ai cicli limite (nella fase PT).

Significato e Rivendicazioni
Il documento sostiene di aver dimostrato che l'accoppiamento di sistemi con distinti meccanismi di rottura della simmetria dissipativa può generare fasi dinamiche emergenti, inclusi il caos e cicli limite complessi, che sono assenti da ciascun modello costitente. Gli autori propongono questo lavoro come un quadro generale per progettare fasi in sistemi guidati-dissipativi accoppiando modelli con fasi note attraverso interazioni che supportano distinte transizioni fuori equilibrio.

Gli autori affermano che il loro modello è accessibile con le attuali tecnologie sperimentali, citando realizzazioni esistenti di modelli di tipo LMG in qudit superconduttori, Dysprosium ultrafreddo, QED a cavità e QED a circuito, combinati con sviluppi nella dissipazione ingegnerizzata. Suggeriscono che tali tecniche potrebbero essere rilevanti per tecnologie quantistiche che richiedono stati stazionari ingegnerizzati o punti critici, come il calcolo quantistico dissipativo, il sensing quantistico e la termodinamica quantistica. Il lavoro evidenzia il potempo di scoprire nuove classi di universalità e comportamenti di scaling delle dimensioni finite, sebbene gli autori notino che determinare tali proprietà richiederebbe approcci numerici più specializzati.