Exact and Efficient Stabilizer Simulation of Thermal-Relaxation Noise for Quantum Error Correction

Enunciato del Problema
La simulazione accurata dei codici di correzione degli errori quantistici (QEC) in presenza di rumore realistico è essenziale per l'addestramento dei decodificatori e per la comprensione delle capacità di soppressione di tali codici. Tuttavia, i framework di simulazione standard basati sugli stabilizzatori, che sono efficienti per sistemi su larga scala, si basano sul teorema di Gottesman-Knill e sono limitati alle operazioni di Clifford e al rumore Pauli stocastico. Il rumore fisico realistico, in particolare il rilassamento termico (caratterizzato dallo smorzamento di ampiezza $T_1$ e dalla dephasizzazione $T_2$), è non-Clifford e non-Pauli.

Per rendere il rumore termico trattabile nei simulatori di stabilizzatori, i ricercatori impiegano tipicamente l'approssimazione di Pauli-twirling (PTA). Tuttavia, la PTA distorce il comportamento dei canali fisicamente rilevanti, spesso fallendo nel catturare il bias direzionale del rilassamento (ad esempio, il decadimento preferenziale di $|1\rangle$ verso $|0\rangle$). Ciò può portare a significative sottostime o sovrastime dei tassi di errore logico, a seconda dello stato di input e della distanza del codice. In alternativa, la simulazione esatta tramite decomposizione quasi-probabilistica (QPD) è possibile ma comporta un sovraccarico di campionamento esponenziale a causa della negatività nella distribuzione di probabilità, rendendola impraticabile per codici su larga scala.

Metodologia
Gli autori sviluppano un modello esatto e compatibile con gli stabilizzatori per il rumore di rilassamento termico dei qubit, trattando i canali di smorzamento di ampiezza ($T_1$) e dephasizzazione ($T_2$) come un processo unificato anziché come errori indipendenti.

1. Decomposizione Unificata: Il lavoro deriva una decomposizione quasi-probabilistica (QPD) per il canale combinato di rilassamento termico. Gli autori dimostrano che quando $T_2 \leq T_1$ — un regime tipico per i qubit superconduttori e allo stato solido — il canale combinato ammette una decomposizione probabilistica completamente positiva in operazioni di Clifford e operazioni di reset. In questo regime, la negatività solitamente associata alla QPD scompare, eliminando il sovraccarico di campionamento esponenziale.
2. Approssimazione per $T_2 > T_1$: Per il regime in cui $T_2 > T_1$, la decomposizione mantiene una certa negatività. Tuttavia, gli autori mostrano che combinare i canali risulta comunque in un sovraccarico di campionamento inferiore rispetto al trattamento dello smorzamento di ampiezza e di fase come canali indipendenti. Inoltre, introducono un canale di errore approssimato con reset che rimuove completamente la negatività. Questa approssimazione mantiene la completezza positiva ottenendo al contempo una fedeltà del canale superiore al vero rilassamento termico rispetto al modello PTA.
3. Estensione a Temperatura Finita: La costruzione è estesa al rilassamento a temperatura finita, dove è possibile l'eccitazione termica. Gli autori mostrano che il costo di campionamento rimane coerente con il caso a temperatura zero per la stessa probabilità di rilassamento, e l'approssimazione basata sul reset rimane superiore alla PTA nei regimi di temperatura rilevanti sperimentalmente.
4. Framework di Simulazione: Questi modelli sono integrati in un nuovo simulatore di stabilizzatori accelerato da MPI e abilitato per GPU (STABSim). Ciò consente la simulazione di grandi codici QEC con modelli di rumore realistici che precedentemente erano intrattabili con metodi di matrice densità o troppo costosi con la QPD standard.

Risultati Chiave
Gli autori hanno applicato il loro modello combinato esatto e l'approssimazione basata sul reset per investigare grandi codici di superficie e codici a bicicletta bivariata (BB) su piattaforme superconduttrici.

• Imprecisione della PTA: Lo studio conferma che la PTA può sottostimare o sovrastimare i tassi di errore logico di fattori da 2 a 10, a seconda della distanza del codice e dello stato logico di base. Nello specifico, la PTA tende a sottostimare gli errori per gli stati logici $|0\rangle_L$ e a sovrastimarli per gli stati logici $|+\rangle_L$.
• Performance Dipendente dallo Stato: È stato riscontrato che i tassi di errore logico dipendono dallo stato logico codificato. Ciò suggerisce che la risposta del decodificatore al bias del rumore è un fattore critico, poiché le popolazioni degli stati eccitati di diversi stati logici sono simili, eppure i loro tassi di soppressione degli errori differiscono.
• Confronto tra Codici:
  • Codici di Superficie: Il modello termico esatto ha rivelato comportamenti di scala che si discostavano dalle previsioni della PTA, evidenziando la necessità di decodificatori informati dal modello di rumore.
  • Codici a Bicicletta Bivariata (BB): Utilizzando il decodificatore BP-OSD, i codici BB hanno mostrato tassi di errore logico circa 1,5-2 volte inferiori a quelli previsti dalla PTA. L'alto tasso di codifica dei codici BB è risultato efficace anche in presenza di rumore termico realistico.
• Efficienza: Il modello composito ha permesso la simulazione efficiente di codici con distanze fino a $d=11$ e codici BB con centinaia di qubit fisici, cosa che sarebbe stata impossibile con metodi a scala esponenziale.

Significato e Affermazioni
Il lavoro afferma di stabilire una via pratica per incorporare il rumore di rilassamento termico fisicamente realistico in framework di simulazione rapidi basati sugli stabilizzatori. Combinando lo smorzamento di ampiezza e la dephasizzazione, gli autori mostrano che il "punto dolce" di $T_2 \leq T_1$ consente una simulazione di Clifford esatta e priva di sovraccarichi. Per i regimi in cui l'esattezza richiede un sovraccarico di campionamento, l'approssimazione basata sul reset proposta offre un'alternativa superiore alla PTA, bilanciando accuratezza ed efficienza.

Gli autori sottolineano che questi risultati indicano che i decodificatori informati dal modello di rumore sono essenziali per le future architetture fault-tolerant per catturare accuratamente le strutture del rumore termico. Suggeriscono inoltre che il loro approccio apre la strada alla simulazione di altri processi non-Clifford combinandoli con altre operazioni per minimizzare la negatività, potenzialmente espandendo la portata dei simulatori classici per la ricerca sulla correzione degli errori quantistici.