2025-11-12T10:37:09.722784

Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states

Calcluth, Hahn, Bermejo-Vega et al.
Classically simulating circuits with bosonic codes is challenging due to the prohibitive cost of simulating quantum systems with many, possibly infinite, energy levels. We propose an algorithm to simulate circuits with encoded Gottesman-Kitaev-Preskill (GKP) states, specifically for odd-dimensional encoded qudits. Our approach is tailored to be especially effective in the most challenging but practically relevant regime, where the codeword states exhibit high (but finite) squeezing. Our algorithm leverages the Zak-Gross Wigner function introduced by J. Davis et al. [arXiv:2407.18394], which represents infinitely squeezed encoded stabilizer states positively. The runtime of the algorithm scales with the negativity of the Wigner function, allowing for efficient simulation of certain large-scale circuits - namely, input stabilizer GKP states undergoing generalized GKP-encoded Clifford operations followed by modular measurements - with a high degree of squeezing. For stabilizer GKP states exhibiting 12 dB of squeezing, our algorithm can simulate circuits with up to 1,000 modes with less than double the number of samples required for a single input mode, in stark contrast to existing simulators. Therefore, this approach holds significant potential for benchmarking early implementations of quantum computing architectures utilizing bosonic codes.
academic

Simulazione classica di circuiti con stati Gottesman-Kitaev-Preskill di dimensione dispari realistica

Informazioni Fondamentali

  • ID Articolo: 2412.13136
  • Titolo: Classical simulation of circuits with realistic odd-dimensional Gottesman-Kitaev-Preskill states
  • Autori: Cameron Calcluth, Oliver Hahn, Juani Bermejo-Vega, Alessandro Ferraro, Giulia Ferrini
  • Classificazione: quant-ph (Fisica Quantistica)
  • Data di Pubblicazione: 17 dicembre 2024 (preprint arXiv)
  • Link dell'Articolo: https://arxiv.org/abs/2412.13136

Riassunto

Questo articolo propone un algoritmo di simulazione classica per simulare circuiti quantistici di dimensione dispari basati sulla codifica Gottesman-Kitaev-Preskill (GKP). L'algoritmo è specificamente ottimizzato per lo scenario più impegnativo ma praticamente rilevante di stati di codeword con compressione elevata ma finita. L'algoritmo utilizza la funzione di Wigner di Zak-Gross, che può rappresentare positivamente gli stati stabilizzatori codificati con compressione infinita. Il tempo di esecuzione dell'algoritmo è proporzionale alla negatività della funzione di Wigner, rendendo possibile la simulazione efficiente di alcuni circuiti su larga scala. Per stati GKP stabilizzatori con compressione di 12 dB, l'algoritmo può simulare circuiti con fino a 1000 modi, richiedendo un numero di campioni inferiore al doppio di quello necessario per un singolo modo di ingresso, in netto contrasto con i simulatori esistenti.

Contesto di Ricerca e Motivazione

Definizione del Problema

La simulazione classica di circuiti quantistici codificati bosonici affronta sfide enormi, principalmente perché è necessario simulare sistemi quantistici con molti (potenzialmente infiniti) livelli energetici, con costi computazionali estremamente elevati. Gli algoritmi di simulazione esistenti presentano limitazioni significative nel trattare stati GKP con compressione finita praticamente rilevanti.

Analisi dell'Importanza

  1. Valore Pratico: La codifica GKP è stata realizzata in sistemi sperimentali quali circuiti superconduttori accoppiati a cavità a microonde, ioni intrappolati e piattaforme fotoniche, rappresentando un candidato importante per il calcolo quantistico tollerante agli errori
  2. Significato Teorico: Fornisce strumenti di benchmark per verificare le implementazioni di processori quantistici bosonici di prima generazione
  3. Sfide Tecniche: I metodi tradizionali non riescono a gestire la simulazione di circuiti su larga scala con stati GKP altamente compressi

Limitazioni dei Metodi Esistenti

  1. Metodo di Espansione nella Base di Fock: Poiché le codeword GKP coinvolgono infiniti livelli energetici, i metodi di simulazione bruta diventano infattibili per pochi qubit GKP
  2. Metodo di Distribuzione Quasi-Probabilistica: La funzione di Wigner degli stati GKP codificati presenta grande negatività, causando un'esplosione del tempo di esecuzione
  3. Metodo di Decomposizione di Stati Gaussiani: Non è applicabile agli stati GKP ad alta compressione praticamente rilevanti, poiché questi stati presentano grande non-gaussianità

Contributi Fondamentali

  1. Proposta della funzione di Wigner di Zak-Gross multimodale: Generalizzazione della funzione di Wigner di Zak-Gross a singolo modo al caso n-modale, con dimostrazione che soddisfa gli assiomi di Stratonovich-Weyl modificati
  2. Stabilimento della relazione tra la funzione di Wigner di Zak-Gross e la funzione di Wigner di Gross: Dimostrazione della corrispondenza tra la funzione di Wigner di Zak-Gross degli stati GKP CV e la funzione di Wigner di Gross degli stati logici codificati
  3. Sviluppo di un algoritmo di simulazione efficiente:
    • Algoritmo di simulazione debole per stati GKP stabilizzatori ideali
    • Algoritmo di stima probabilistica per stati GKP con compressione finita
  4. Realizzazione di miglioramenti significativi delle prestazioni: Per stati GKP stabilizzatori con compressione di 12 dB, la simulazione di circuiti con 1000 modi richiede solo un numero di campioni inferiore al doppio di quello necessario per un singolo modo

Dettagli Metodologici

Definizione del Compito

Simulazione della distribuzione di probabilità dei risultati di misurazione del seguente circuito quantistico:

  • Ingresso: Qubit quantistico arbitrario codificato GKP, con parametro di anti-compressione Δ
  • Evoluzione: Operazioni unitarie descritte da matrici simplettiche intere e spostamenti arbitrari (incluse tutte le operazioni Clifford codificate)
  • Misurazione: Misurazione modulo nella base di posizione MZ(s)=1dneisnT(0,n)M_Z(s) = \frac{1}{d\ell}\sum_n e^{-i\ell s \cdot n}T_{(0,n)}

Definizione della Funzione di Wigner di Zak-Gross

Definizione 1: La funzione di Wigner di Zak-Gross multimodale di dimensione dispari è definita come: Wρ^(η)=Tr(ρ^A^η)W_{\hat{\rho}}(\eta) = \text{Tr}(\hat{\rho}\hat{A}_\eta)

dove l'operatore di punto di fase è: A^η=1(2π)naZ2nei[a,η]+iπaXTaZT^a\hat{A}_\eta = \frac{1}{(2\pi)^n}\sum_{a \in \mathbb{Z}^{2n}} e^{i\ell[a,\eta]+i\pi a_X^T a_Z}\hat{T}_a

dove η[0,d)2n\eta \in [0,d\ell)^{2n}, dd è un intero positivo dispari, T^a=eiπaXTaZ/dT^aXT^aZ\hat{T}_a = e^{i\pi a_X^T a_Z/d}\hat{T}_{a_X}\hat{T}_{a_Z}.

Risultati Teorici Chiave

Lemma 1: Relazione tra la funzione di Wigner di Zak-Gross e la funzione di Wigner di Gross Wρ^(η)=Wˉρˉ(u)(t)W_{\hat{\rho}}(\eta) = \bar{W}_{\bar{\rho}(u)}(t) dove η=(u+t)\eta = \ell(u+t), u=1ηmod1u = \frac{1}{\ell}\eta \bmod 1, t=1ηut = \frac{1}{\ell}\eta - u.

Teorema 1: Trasformazione della funzione di Wigner di Zak-Gross sotto evoluzione simplettica Per l'evoluzione unitaria gaussiana descritta dalla matrice simplettica intera SS, la funzione di Wigner di Zak-Gross dello stato evoluto è: WU^Sρ^U^S(η)=Wρ^(Sηt)W_{\hat{U}_S\hat{\rho}\hat{U}_S^\dagger}(\eta) = W_{\hat{\rho}}(S\eta - t)

Teorema 2: Calcolo della probabilità di misurazione La distribuzione di probabilità di misurazione dell'operatore logico Z^Ln\hat{Z}_L^{\otimes n} è: Tr(ρ^M^Z(s))=dηXWρ^((ηXs))\text{Tr}(\hat{\rho}\hat{M}_Z(s)) = \int d\eta_X W_{\hat{\rho}}\left(\begin{pmatrix}\eta_X\\s\end{pmatrix}\right)

Algoritmo di Simulazione

Simulazione Debole per Stati GKP Stabilizzatori Ideali

  1. Campionamento del vettore η\eta dalla distribuzione di probabilità effettiva Wρ^0(η)W_{\hat{\rho}_0}(\eta)
  2. Trasformazione sotto la matrice simplettica SS e il vettore di spostamento tt: SηtS\eta - t
  3. Applicazione dell'operatore di spostamento lineare: Sηt+cS\eta - t + c
  4. Il risultato della misurazione è dato dalla seconda metà del vettore

Stima Probabilistica per Stati GKP con Compressione Finita

Utilizzo del metodo di campionamento con distribuzione quasi-probabilistica: Pr(η)=1MU^S,cρ^0U^S,cWU^S,cρ^0U^S,c(η)\text{Pr}(\eta) = \frac{1}{M_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}}|W_{\hat{U}_{S,c}\hat{\rho}_0\hat{U}_{S,c}^\dagger}(\eta)|

Il numero di campioni richiesti è: N=2ϵ2Mρ^02log(2/δ)N = \frac{2}{\epsilon^2}M_{\hat{\rho}_0}^2\log(2/\delta)

dove Mρ^M_{\hat{\rho}} è la negatività della funzione di Wigner di Zak-Gross (volume negativo).

Configurazione Sperimentale

Quadro di Analisi Teorica

Questo articolo si concentra principalmente su analisi teorica e calcoli numerici, con attenzione particolare a:

  1. Calcolo della negatività della funzione di Wigner di Zak-Gross: Per stati GKP con diversi livelli di compressione
  2. Analisi della complessità campionaria: Basata sul comportamento di scaling della negatività
  3. Confronto con metodi esistenti: Confronto della complessità teorica

Metriche di Valutazione

  • Negatività logaritmica di Zak-Gross: logMρ^\log M_{\hat{\rho}}, che misura la complessità di simulazione
  • Scaling del numero di campioni: Crescita del numero di campioni rispetto al caso a singolo modo
  • Capacità di modi: Numero massimo di modi che possono essere simulati efficacemente

Risultati Sperimentali

Funzione di Wigner di Zak-Gross per Stati GKP con Compressione Finita

Attraverso il calcolo esplicito, si ottiene la funzione di Wigner di Zak-Gross dello stato GKP logico 0 realistico: WGKP((uv))ϑ(Γ;z)W_{\text{GKP}}\left(\begin{pmatrix}u\\v\end{pmatrix}\right) \propto \vartheta(\Gamma; z)

dove z=(v/(d),u/(d),0,0)Tz = (v/(d\ell), -u/(d\ell), 0, 0)^T, Γ\Gamma è una matrice complessa determinata dal parametro di compressione Δ\Delta.

Risultati Numerici Principali

Per Δ=0.25\Delta = 0.25 (corrispondente a compressione di 12 dB):

  • Negatività a singolo modo: Mρ^0Δe3×104M_{\hat{\rho}_0^\Delta} \approx e^{3 \times 10^{-4}}
  • Simulazione a mille modi: Grazie alla natura moltiplicativa della negatività, la simulazione di 1000 modi di ingresso richiede solo un numero di campioni inferiore al doppio di quello necessario per un singolo modo
  • Miglioramento delle prestazioni: Miglioramento di diversi ordini di grandezza rispetto ai metodi esistenti

Analisi dello Stato Magico

Per lo stato π ψπ=13(0L+1L2L)|\psi_\pi\rangle = \frac{1}{\sqrt{3}}(|0_L\rangle + |1_L\rangle - |2_L\rangle):

  • Presenta negatività significativa per tutti i valori di Δ\Delta
  • Costi di simulazione notevoli per circuiti su larga scala, come previsto

Lavori Correlati

Sviluppo degli Algoritmi di Simulazione Classica

  1. Teorema di Gottesman-Knill: Simulazione in tempo polinomiale di circuiti stabilizzatori
  2. Metodo di Distribuzione Quasi-Probabilistica: Algoritmo di campionamento basato sulla funzione di Wigner di Gross
  3. Simulazione di Sistemi Bosonici: Utilizzo di distribuzione quasi-probabilistica positiva, reti tensoriali, rappresentazione stellare e altri metodi

Stato Attuale della Simulazione di Codifica GKP

I metodi esistenti sono principalmente limitati a:

  • Stati GKP ideali (compressione infinita)
  • Operazioni gaussiane e misurazioni homodyne
  • Incapacità di gestire situazioni praticamente rilevanti con compressione finita

Punti di Innovazione di Questo Articolo

Rispetto ai lavori esistenti, questo articolo realizza per la prima volta la simulazione efficiente di stati GKP con compressione finita, in particolare con miglioramenti significativi delle prestazioni nel regime ad alta compressione.

Conclusioni e Discussione

Conclusioni Principali

  1. Contributo Teorico: Stabilimento di un quadro teorico completo della funzione di Wigner di Zak-Gross, inclusa la generalizzazione multimodale e le leggi di evoluzione
  2. Innovazione Algoritmica: Proposta di un algoritmo di simulazione efficiente applicabile a stati GKP realistici
  3. Valore Pratico: Fornimento di uno strumento di benchmark potente per la verifica di processori quantistici bosonici

Limitazioni

  1. Restrizione alla Dimensione Dispari: I metodi attuali sono applicabili solo a qubit GKP di dimensione dispari; l'estensione al caso pari (come qubit) rimane un problema aperto
  2. Restrizione Operativa: Limitato alle operazioni descritte da matrici simplettiche intere, non può gestire compressione a singolo modo arbitraria
  3. Modello di Rumore: Non considera l'effetto del rumore realistico come la perdita sulla negatività della funzione di Wigner di Zak-Gross

Direzioni Future

  1. Estensione a Qubit GKP: Sviluppo di metodi applicabili al caso bidimensionale
  2. Codifiche Bosoniche Più Generali: Generalizzazione del quadro a codifiche bosoniche oltre GKP
  3. Quadro di Distribuzione Quasi-Probabilistica: Sviluppo di una teoria quasi-probabilistica completa specializzata per codifiche bosoniche

Valutazione Approfondita

Punti di Forza

  1. Rigore Teorico: Derivazioni matematiche complete, dimostrazioni dettagliate, fondamenti teorici solidi
  2. Forte Praticità: Risolve il problema della simulazione di stati GKP con compressione finita praticamente rilevanti
  3. Prestazioni Eccezionali: Realizza miglioramenti di ordini di grandezza nel regime ad alta compressione
  4. Innovazione Metodologica: La generalizzazione multimodale della funzione di Wigner di Zak-Gross e la sua applicazione rappresentano importanti contributi teorici

Insufficienze

  1. Ambito di Applicabilità Limitato: Applicabile solo a sistemi di dimensione dispari, il che limita l'universalità
  2. Mancanza di Verifica Sperimentale: Principalmente lavoro teorico, mancano confronti con sistemi quantistici reali
  3. Analisi della Complessità: L'analisi della complessità computazionale effettiva per circuiti su larga scala non è sufficientemente dettagliata

Impatto

  1. Valore Accademico: Fornisce nuovi strumenti e prospettive per la teoria della simulazione quantistica
  2. Prospettive Pratiche: Ha importanza significativa per l'orientamento teorico degli esperimenti di calcolo quantistico tollerante agli errori di prima generazione
  3. Spinta Tecnologica: Potrebbe ispirare lo sviluppo di più metodi di simulazione di codifiche bosoniche

Scenari di Applicazione

  • Benchmark e verifica di processori quantistici bosonici
  • Analisi teorica di schemi di calcolo quantistico tollerante agli errori con codifica GKP
  • Ricerca sulla simulazione classica di calcolo quantistico a variabili continue

Bibliografia

Questo articolo cita 68 riferimenti correlati, principalmente includenti:

  • Lavori originali sulla codifica Gottesman-Kitaev-Preskill
  • Letteratura classica sugli algoritmi di simulazione quantistica
  • Progressi recenti nel calcolo quantistico bosonico
  • Fondamenti teorici delle funzioni di Wigner e distribuzioni quasi-probabilistiche

Valutazione Complessiva: Questo è un articolo di fisica teorica di alta qualità che ha raggiunto importanti progressi nella simulazione classica di circuiti quantistici codificati GKP. Sebbene presenti alcune limitazioni, i suoi contributi teorici e il valore pratico sono entrambi significativi e hanno importanza rilevante per il campo del calcolo quantistico.