Real-time Krylov Diagonalisation for Open Quantum Systems

Riepilogo Tecnico: Diagonalizzazione di Krylov in Tempo Reale per Sistemi Quantistici Aperti

Enunciato del Problema
Il lavoro affronta la sfida di simulare sistemi quantistici aperti descritti da equazioni maestre di Lindblad, focalizzandosi specificamente sulla stima del gap di Liouvilliano. Mentre la Diagonalizzazione del Sottospazio Quantistico (QSD) è emersa come un algoritmo leader non tollerante ai guasti per sistemi chiusi (Hamiltoniani), la sua applicazione a sistemi dissipativi non è banale. I sistemi aperti sono governati da superoperatori di Liouvilliano non hermitiani e non normali, che introducono differenze strutturali come autovettori bi-ortogonali, autovalori complessi e dinamiche non unitarie. Queste caratteristiche complicano i metodi standard del sottospazio di Krylov, in particolare per quanto riguarda la stabilità numerica (malcondizionamento) e l'interpretazione dei filtri spettrali. Gli autori mirano ad adattare la QSD in tempo reale per stimare il gap di Liouvilliano, un parametro critico che definisce il tempo di rilassamento e la presenza di modi lenti nelle transizioni di fase dissipative.

Metodologia
Gli autori propongono una modifica della Diagonalizzazione del Sottospazio Quantistico (QSD) in tempo reale per gestire le dinamiche di Liouvilliano. La metodologia centrale comporta:

1. Costruzione del Sottospazio di Krylov: Invece di evolvere uno stato sotto un Hamiltoniano unitario $H$, il sistema evolve sotto il superoperatore di Liouvilliano $\mathcal{L}$. Un sottospazio di Krylov è generato applicando potenze dell'operatore di evoluzione temporale $e^{\mathcal{L}\tau}$ a uno stato iniziale $|\rho_0)$:
   $$\text{span}\{|\rho_0), e^{\mathcal{L}\tau}|\rho_0), \dots, e^{\mathcal{L}(D-1)\tau}|\rho_0)\}$$
2. Problema agli Autovalori Generalizzato (GEVP): Poiché i circuiti quantistici non possono eseguire ri-ortogonalizzazione in tempo reale, il metodo risolve un GEVP $\bar{H}v = \lambda Sv$. Per i Liouvilliani, la matrice di sovrapposizione $S$ e la matrice Hamiltoniana $\bar{H}$ sono costruite utilizzando autovettori sinistri e destri bi-ortogonali.
3. Interpretazione del Filtro Spettrale: Il metodo è inquadrato come ottimizzazione di un filtro spettrale. A differenza dei sistemi Hamiltoniani dove gli autostati sono reali e l'evoluzione temporale è oscillatoria (richiedendo un campionamento di Nyquist rigoroso), le dinamiche di Liouvilliano sono dissipative. I modi a decadimento rapido spiraleggiano verso l'interno nel piano complesso, permettendo tempi di evoluzione $\tau$ più lunghi senza artefatti di aliasing. Ciò consente l'isolamento dei modi lenti (autovalori con parti reali piccole) anche con una dimensione di Krylov ridotta $D$.
4. Gestione della Non Normalità: Gli autori riconoscono che gli operatori non normali portano ad un'alta sensibilità alle perturbazioni (teorema di Bauer-Ficke). La matrice di sovrapposizione $S$ diventa esponenzialmente malcondizionata all'aumentare di $D$. Per mitigare ciò, il metodo si affida alla natura dissipativa del sistema per sopprimere i modi rapidi, permettendo $\tau$ più grandi e $D$ più piccoli, riducendo così l'instabilità numerica associata a grandi dimensioni di Krylov.
5. Selezione dello Stato Iniziale: Crucialmente, lo stato iniziale deve avere una proiezione non nulla sul settore a traccia zero (coerenze o squilibri di popolazione) per eccitare i modi lenti. Il settore a traccia 1 (stato stazionario) non contribuisce al segnale per la stima del gap in questo quadro.

Contributi Chiave

• Estensione della QSD ai Liouvilliani: Il lavoro fornisce un quadro teorico per applicare la QSD in tempo reale a superoperatori di Liouvilliano non hermitiani e non normali in forma di Lindblad.
• Strategia di Stima del Gap: Dimostra come stimare il gap di Liouvilliano ($\Delta = -\max \Re(\lambda_k)$) sfruttando la natura dissipativa del sistema per filtrare i modi a decadimento rapido, aggirando efficacemente alcuni vincoli di Nyquist che si applicano all'evoluzione unitaria.
• Analisi della Stabilità Numerica: Il lavoro evidenzia il compromesso tra il tempo di evoluzione $\tau$ e la dimensione di Krylov $D$. Sostiene che aumentare $\tau$ (piuttosto che $D$) è una strategia praticabile per i sistemi aperti per migliorare la risoluzione del gap mantenendo la stabilità numerica, poiché i modi rapidi decadono naturalmente.
• Applicazione ai Qubit di Gatto di Kerr: Il metodo è applicato a un risonatore di Kerr superconduttore guidato a due fotoni (Qubit di Gatto di Kerr), un sistema rilevante per le architetture di qubit protetti.

Risultati
Gli autori hanno applicato il loro metodo a un modello semplificato di un Qubit di Gatto di Kerr con perdita di fotoni singoli.

• Parametri di Simulazione: Le simulazioni sono state eseguite con un massimo di 30 fotoni e una dimensione di Krylov $D=20$, risultando in una matrice di Krylov $900 \times 20$.
• Ricostruzione del Gap: Il metodo ha ricostruito con successo il gap di Liouvilliano attraverso diverse intensità di guida. I risultati hanno mostrato che l'aumento del tempo di evoluzione $\tau$ ha migliorato l'accuratezza della stima del gap, in particolare nei regimi in cui il commutatore di Liouvilliano (una misura della non normalità) era elevato.
• Identificazione del Modo Lento: Le simulazioni hanno confermato l'emergere di un modo lento corrispondente al decadimento delle coerenze e degli squilibri di popolazione. Nel regime di "rumore polarizzato", il gap si restringe esponenzialmente con l'ampiezza dello stato coerente ($\propto e^{-\alpha^2}$), indicando la formazione di un sottospazio computazionale protetto.
• Rappresentazione di Wigner: Gli autostati ricostruiti dalla matrice di Krylov sono stati utilizzati per visualizzare la rappresentazione di Wigner del modo lento, mostrando la transizione da un regime smorzato a un regime con due lobi distinti (rottura di simmetria) e ritorno a un regime compresso ad alta guida.

Significato e Affermazioni
Il lavoro afferma di dimostrare la fattibilità dell'uso della QSD in tempo reale per studiare sistemi quantistici aperti, specificamente per la stima del gap di Liouvilliano nel contesto delle transizioni di fase dissipative. Gli autori collocano questo lavoro come un adattamento necessario degli algoritmi dell'era NISQ per la progettazione e la simulazione di hardware quantistico, che coinvolge intrinsecamente la dissipazione.

Il significato è inquadrato attorno alla potenziale possibilità di utilizzare computer quantistici per simulare i modi lenti di sistemi quantistici che sono computazionalmente costosi da simulare classicamente, in particolare nel contesto di qubit protetti come il Gatto di Kerr. Gli autori notano modestamente che, mentre singoli Qubit di Gatto di Kerr sono simulabili classicamente, reticoli di tali qubit potrebbero non esserlo, suggerendo l'utilità del metodo per sistemi più grandi e complessi. Il lavoro è presentato come un passo verso la "scoperta automatizzata dell'hardware quantistico" fornendo strumenti per caratterizzare le proprietà del rumore e la stabilità delle architetture di qubit proposte. Il lavoro non afferma di aver risolto il problema generale della simulazione di operatori non normali, ma offre piuttosto un approccio specifico e modificato per la stima del gap in sistemi dissipativi.