Non-Stationary Decoherence in Superconducting Qubits: Memory Multi-Fractional Brownian Motion and a Time-Dependent Quantum Brownian Motion Extension

Riepilogo Tecnico: Decoerenza Non Stazionaria in Qubit Superconduttori

Enunciato del Problema
I qubit superconduttori sono altamente sensibili alle fluttuazioni ambientali a bassa frequenza, in particolare al rumore $1/f^\beta$ derivante da difetti dielettrici, sistemi a due livelli (TLS) e quasiparticelle. Le osservazioni sperimentali rivelano coerentemente un decadimento non esponenziale nella coerenza di Ramsey e nell'eco di spin, insieme a correlazioni temporali a lungo raggio che non possono essere catturate da modelli gaussiani stazionari o da approcci puramente markoviani. I modelli esistenti non riescono a rendere conto di due caratteristiche critiche osservate nei dispositivi moderni: (i) la lenta deriva temporale dell'esponente spettrale $\beta(t)$, e (ii) correlazioni a lungo raggio che si estendono su molti cicli di gate. Questo lavoro affronta la necessità di un quadro unificato che descriva la decoerenza non stazionaria a memoria lunga nei qubit di carica superconduttori.

Metodologia
Gli autori propongono un modello unificato di deriva stocastica (SdM) fondato sul Movimento Browniano Frazionario Multi-Frazionale con Memoria (mmfBm) e lo estendono in un quadro quantistico microscopico.

1. Settore Classico (mmfBm): Il modello definisce l'evoluzione della coordinata del qubit guidata da un processo gaussiano $M(t)$ con un esponente di Hurst dipendente dal tempo $H(t) \in (1/2, 1)$. Il processo è costruito tramite un integrale frazionario di Riemann–Liouville di rumore di Wiener standard, incorporando un kernel di memoria adattivo $K(t, s) = (t-s)^{H(t)-1/2}$. Questa formulazione permette alla rugosità locale del rumore di variare lentamente nel tempo, generando rumore non stazionario $1/f^\beta$ dove l'esponente spettrale è legato al parametro di Hurst da $\beta(t) = 2H(t) - 1$.
2. Estensione Quantistica: Il modello stocastico classico è incorporato in un ambiente quantistico microscopico utilizzando una costruzione Caldeira–Leggett dipendente dal tempo. L'ambiente è modellato come un bagno bosonico con una densità spettrale dipendente dal tempo $J(\omega; t) = \eta(t) \omega_c^{1-s(t)} \omega^{s(t)} e^{-\omega/\omega_c}$, dove $s(t) = 2H(t) - 1$.
3. Derivazioni e Simulazioni:
   • Gli autori derivano il funzionale di influenza di Feynman–Vernon e la conseguente equazione maestra di Lindblad dipendente dal tempo sotto le approssimazioni di Born–Markov e adiabatica.
   • Stabiliscono la connessione tra il bagno quantistico e il limite classico mmfBm nel regime ad alta temperatura.
   • Le simulazioni numeriche sono eseguite utilizzando la decomposizione di Cholesky per la generazione di traiettorie mmfBm e la libreria QuTiP per risolvere l'equazione maestra dipendente dal tempo.
   • Viene esplorato brevemente un approccio di apprendimento automatico (CNN) per inferire l'esponente di Hurst dalle curve di coerenza, sebbene le sue limitazioni siano esplicitamente notate.

Contributi Chiave

• Quadro Unificato: Il lavoro stabilisce un ponte teorico diretto tra la fenomenologia del rumore stocastico (mmfBm) e la dinamica microscopica dei sistemi aperti (bagno Caldeira–Leggett dipendente dal tempo).
• Densità Spettrale Dipendente dal Tempo: Introduce una densità spettrale $J(\omega; t)$ che riproduce coerentemente l'esponente di rumore variabile nel tempo osservato sperimentalmente $\beta(t) = 2H(t) - 1$.
• Risultati Analitici:
  • Derivazione dell'equazione maestra di Lindblad tempo-locale efficace per ambienti non stazionari.
  • Dimostrazione che il modello classico mmfBm emerge come limite ad alta temperatura, debole retroazione e adiabatico del bagno quantistico.
  • Stabilimento di leggi di scala per la coerenza di Ramsey e di eco, prevedendo inviluppi di decadimento a stiramento esponenziale della forma $\exp[-A t^{2H(t)}]$ e $\exp[-A t^{2H(t)+2}]$, rispettivamente.
• Validazione Numerica: Lo studio fornisce stime numeriche per i tempi di rilassamento ($T_1$) e di dephasing ($T_2$) sotto il dominio del rumore di carica, ottenendo $T_1 \sim 5.00 \times 10^6$ ns e $T_2 \sim 4.18 \times 10^5$ ns.

Risultati

• Dinamica di Decoerenza: Le simulazioni dimostrano che la decoerenza è altamente sensibile all'interplay tra deriva, ampiezza del rumore ed effetti di memoria. I profili di Hurst dipendenti dal tempo generano un comportamento non stazionario più realistico rispetto ai modelli a $H$ costante.
• Comportamento Non Markoviano: Il sistema esibisce dephasing non markoviano con pattern di decadimento a stiramento esponenziale, confermando che i modelli di decadimento esponenziale semplice sono insufficienti per descrivere ambienti di qubit dominati da rumore a memoria lunga.
• Incrocio di Temperatura: Viene osservata una chiara transizione tra un regime dominato dal vuoto quantistico (dove il dephasing è indipendente dalla temperatura) e un regime termico (dove $\Gamma_\phi \propto T$). La scala temporale di incrocio è identificata come $t_{cross} \sim \hbar/k_B T$.
• Ottimizzazione dei Gate: La competizione tra rilassamento energetico e dephasing puro rivela un tempo di gate ottimale che minimizza l'errore totale, bilanciando il decadimento da rilassamento contro il rapido dephasing a tempi brevi.
• Limiti dell'Apprendimento Automatico: Il tentativo di estrarre $H(t)$ utilizzando una CNN ha prodotto una scarsa generalizzazione ($R^2 = -4.05$), attribuita all'uso di dati di formazione semplificati e all'incapacità delle architetture di pooling globale di catturare la variabilità temporale. Il lavoro suggerisce che per il lavoro futuro sono necessarie architetture consapevoli delle sequenze (es. LSTMs, Transformers).

Significato e Affermazioni
Il lavoro afferma di fornire una base flessibile per comprendere la decoerenza non stazionaria, superando i limiti degli approcci markoviani. Afferma che i pattern di decadimento non esponenziale rivelati dal quadro mmfBm evidenziano vincoli fondamentali sui modelli di rumore attuali e offrono una guida per la progettazione di architetture di qubit resilienti al rumore.

Gli autori sono modesti riguardo allo stato attuale del lavoro, riconoscendo diverse limitazioni:

• Le simulazioni a tempo finito non catturano pienamente la scala asintotica a bassa frequenza, portando a discrepanze tra gli esponenti PSD adattati e quelli teorici.
• L'approssimazione localmente stazionaria potrebbe rompersi se la variazione di $H(t)$ è troppo rapida.
• Il trattamento Born–Markov trascura forti effetti di memoria potenzialmente rilevanti a temperature nell'ordine dei millikelvin.
• L'implementazione attuale fatica con la rigidità numerica a temperature ultra-basse a causa di singolarità nell'infrarosso.

Il lavoro conclude che, sebbene il quadro riproduca con successo le caratteristiche fisiche chiave (decadimento a stiramento esponenziale, incrocio di temperatura, scaling non banale), i progressi futuri richiedono solutori non markoviani esatti, tecniche di stima spettrale migliorate e derivazioni microscopiche più rigorose che colleghino i parametri mmfBm a specifici difetti materiali.