Temporal Coarse Graining for Classical Stochastic Noise in Quantum Systems

Problematica
Simulare sistemi quantistici soggetti a rumore stocastico classico hamiltoniano è computazionalmente oneroso quando il rumore presenta correlazioni temporali su un ampio intervallo di scale temporali, come il rumore $1/f$ comune nei processori di informazione quantistica a stato solido. L'integrazione standard brute-force dell'equazione di Schrödinger richiede passi temporali sufficientemente piccoli da risolvere le componenti di rumore ad alta frequenza per evitare l'aliasing. Tuttavia, questi passi sono ordini di grandezza più piccoli delle scale temporali necessarie per osservare gli effetti dei drift lenti derivanti dalle componenti a bassa frequenza. Questa disparità risulta in costi di simulazione proibitivi per esperimenti della durata di secondi, anche quando i processi di rilassamento ($T_1$) sono trascurabili rispetto al dephasing ($T_2$).

Metodologia
Gli autori propongono un approccio di coarso-graining temporale per simulare il rumore stocastico hamiltoniano classico, concentrandosi specificamente su processi di rumore esprimibili come una somma di processi di Ornstein-Uhlenbeck (OU) indipendenti. Il metodo decompone la simulazione in due medie d'insieme distinte:

1. Condizionamento su Realizzazioni Coarse: Il processo di rumore stocastico è condizionato su una realizzazione "coarse" definita su un insieme discreto di punti temporali $\{t_0, t_1, \dots\}$. Tali punti non devono essere necessariamente spaziati uniformemente, permettendo loro di allinearsi con la durata di specifiche operazioni quantistiche (ad esempio, i gate).
2. Decomposizione del Processo Condizionato: All'interno di ogni intervallo $[t_{k-1}, t_k]$, il processo OU condizionato è espresso come la somma di:
   • Una funzione deterministica ($\eta^{(D)}$) che dipende dai valori al contorno (la realizzazione coarse) e cattura il drift lento.
   • Un processo bridge a zero-boundary ($\eta^{(S)}$) che è indipendente dalla specifica realizzazione coarse, ha media nulla e si annulla ai contorni dell'intervallo.
3. Media d'Insieme a Due Fasi:
   • Fase A (Media del Bridge): Viene eseguita una media d'insieme sui processi bridge a zero-boundary. Poiché questi processi sono indipendenti tra gli intervalli temporali e possiedono correlatori noti analiticamente, questa fase integra efficacementamente le componenti di rumore ad alta frequenza. Ciò produce un'operazione quantistica completamente positiva e traccia-preservante (CPTP) per ogni intervallo, caratterizzata da un superoperatore derivato da un'espansione cumulante (troncata al secondo ordine).
   • Fase B (Media Coarse): Il sistema viene evoluto attraverso la sequenza di intervalli utilizzando le specifiche funzioni deterministiche derivate da una singola realizzazione di rumore coarse. Questa evoluzione viene poi ripetuta per molte diverse realizzazioni di rumore coarse, e i risultati finali vengono mediati.

Questo approccio "ibrido" combina il campionamento Monte Carlo di traiettorie di rumore coarse con la media d'insieme analitica dei processi bridge ad alta frequenza.

Contributi Chiave

• Pre-computazione Analitica: Per i processi OU, le componenti deterministiche e i correlatori dei processi bridge a zero-boundary sono derivati analitamente. Ciò consente di pre-calcolare una sola volta i propagatori di rumore associati, eliminando la necessità di generare traiettorie di rumore fine-grained per ogni passo della simulazione.
• Regola di Concatenazione: L'indipendenza dei processi bridge a zero-boundary attraverso gli intervalli temporali permette alla l'evoluzione totale di essere costruita tramite una semplice concatenazione di mappe per ogni passo temporale coarse. Questo contrasta con gli approcci standard di media d'insieme (ad esempio, le funzioni di filtro) che spesso mancano di una semplice regola di concatenazione per sequenze di operazioni unitarie, a meno che l'Hamiltoniana non sia costante o non vengano trascurati i termini di ordine superiore.
• Gestione delle Misurazioni Mid-Circuit: Il metodo preserva naturalmente le correlazioni del rumore attraverso le misurazioni mid-circuit. Poiché la traiettoria del rumore coarse viene portata attraverso la misurazione, lo stato può essere proiettato ed evoluto verso la misurazione successiva senza dover ricalcolare l'intera storia o eseguire un numero esponenziale di simulazioni per tutti i possibili esiti della misurazione.

Risultati
Gli autori dimostrano il metodo attraverso tre esempi numerici:

1. Dephasing a Singolo Qubit: Un sistema a singolo qubit con rumore sull'asse $x$. Il metodo separa con successo il rumore in un errore coerente dipendente dalla traiettoria (drift) e un termine di dephasing indipendente dalla traiettoria.
2. Fluttuazioni dell'Exchange a Due Qubit: Simulazione delle fluttuazioni nell'accoppiamento di scambio tra due spin. I risultati mostrano dephasing nella base propria dell'operatore di scambio e rotazioni eccessive o insufficienti coerenti dipendenti dai valori al contorno.
3. Full Permutation Dynamical Decoupling (DD): Simulazione di un qubit a scambio solo a 3 spin sotto un protocollo di DD a permutazione completa. La simulazione cattura il comportamento oscillatorio della probabilità di sopravvivenza per stati iniziali specifici (dovuto a gradienti di campo magnetico) e il decadimento dovuto alla decoerenza, corrispondendo ai comportamenti fisici attesi.
4. Parity Check di Peso-2: Simulazione di una sequenza di controlli di parità su tre qubit singlet-triplet (6 spin). Gli autori confrontano modelli di rumore $1/f$, rumore quasi-statico e rumore Bernoulli. Il metodo gestisce efficientemente la sequenza di misurazioni mid-circuit, evitando il costo proibitivo di simulare tutte le $2^{N_M}$ sequenze di esiti della misurazione richieste dagli approcci alternativi.

Significatività e Rivendicazioni
L'articolo sostiene che questo approccio di coarse-graining temporale offre un vantaggio pratico per la simulazione di sistemi quantistici con rumore multi-scala. Integrando efficacementmente le componenti ad alta frequenza in modo analitico pur mantenendo la capacità di tracciare drift lenti e correlazioni attraverso traiettorie coarse, il metodo consente una simulazione fedele della dinamica rumorosa su lunghi intervalli temporali (secondi) che sarebbero altrimenti computazionalmente intrattabili. Gli autori evidenziano come l'approccio sia particolarmente adatto per il benchmarking di protocolli e studi di correzione degli errori, dove sono richiesti tempi di simulazione lunghi e strutture di circuiti complessi (inclusi i mid-circuit measurements). Il metodo è presentato come un'alternativa alle descrizioni di processi quantistici efficaci derivate tramite espansioni di Magnus e cumulanti, affrontando specificamente la mancanza di una regola di concatenazione in tali approcci standard per Hamiltoniane tempo-varianti.