Riepilogo Tecnico: Permutazione Fermionica a Profondità Asintoticamente Ottimale su Architettura Quantistica a Griglia 2D senza Ancille

1. Enunciato del Problema

La simulazione di sistemi fermionici interagenti su hardware di qubit è ostacolata dal disallineamento fondamentale tra gli operatori locali dei qubit e le relazioni di anticommutazione non locali degli operatori di creazione e annichilazione fermionici. Ciò rende necessaria una mappatura fermione-qubit (ad es. codifiche di Jordan–Wigner, Bravyi–Kitaev o Parità) e una permutazione fermionica dinamica per instradare i modi interagenti verso qubit adiacenti.

Sulle architetture a griglia 2D con vicini più prossimi (NN), i protocolli di instradamento esistenti presentano un sovraccarico significativo:

• Approcci 1D: La compilazione ingenua di reti di ordinamento 1D (ad es. transposizione pari–dispari) su una griglia 2D comporta una profondità $O(N)$ e un conteggio di porte $O(N^2)$, non sfruttando la seconda dimensione della griglia.
• Adattamenti All-to-All: I recenti protocolli che raggiungono una profondità $O(\log N)$ o $O(\log^2 N)$ su architetture all-to-all o riconfigurabili degradano a $O(\sqrt{N} \text{ polylog } N)$ su griglie 2D a causa del tempo $\Theta(\sqrt{N})$ richiesto per attraversare l'array.
• Limite Inferiore: Esiste un limite inferiore teorico di $\Omega(\sqrt{N})$ per l'instradamento su una griglia 2D, che vale anche nel modello di Operazioni Locali e Comunicazione Classica (LOCC) che consente $O(N)$ ancille, misurazioni a metà circuito e feedforward classico.

La sfida centrale è raggiungere questo limite inferiore di $\Omega(\sqrt{N})$ sull'hardware NN 2D senza fare affidamento su ancille, misurazioni o feedforward, che sono risorse intensive e soggetti a errori nei regimi di tolleranza agli errori iniziali.

2. Metodologia

Gli autori propongono un protocollo di permutazione fermionica specificamente adattato per le architetture a griglia 2D, costituito da due componenti principali: una decomposizione classica di instradamento e un operatore unitario senza ancille per la correzione di fase.

A. Decomposizione Riga-Colonna-Riga (RCR) di Hall

Invece di linearizzare i modi in una catena 1D e compilare un circuito 1D, il protocollo utilizza una decomposizione classica basata sul teorema del matrimonio di Hall. Qualsiasi permutazione $\pi$ su una griglia $L \times L$ ($N=L^2$) è decomposta in tre stadi:

1. Riga A: Permutazione all'interno delle righe sorgente.
2. Colonna: Permutazione all'interno delle colonne per spostare gli elementi verso le righe di destinazione.
3. Riga B: Permutazione all'interno delle righe di destinazione verso le colonne finali.

Secondo l'ordinamento standard "a serpente" (boustrophedon) di Jordan–Wigner (JW), i vicini orizzontali sono adiacenti JW, rendendo gli stadi di Riga naturalmente locali. Lo stadio di Colonna, tuttavia, coinvolge vicini verticali che non sono adiacenti JW, richiedendo correzioni di stringa di parità che le porte FSWAP hardware standard non forniscono.

B. L'Operatore $\Gamma$ (Correzione a Polinomio di Fase)

Per risolvere la non località degli scambi verticali senza ancille, gli autori costruiscono un operatore unitario diagonale $\Gamma$.

• Funzione: $\Gamma$ coniuga gli scambi verticali "nudi" (che ignorano le stringhe di parità) per produrre scambi fermionici "completi" (che includono le necessarie correzioni di fase). Nello specifico, $\Gamma \cdot B_{j,k} \cdot \Gamma = F_{j,k}$, dove $B$ è lo scambio nudo e $F$ è lo scambio fermionico completo.
• Implementazione: Gli autori riformulano $\Gamma$ utilizzando una lente a polinomio di fase. Decompongono la funzione di fase richiesta $f(s)$ in monomi locali che si adattano alla geometria della griglia 2D.
  • Ristrutturazione Algebrica: Il polinomio di fase è suddiviso in $f_D$ (termini della stessa riga e termini incrociati adiacenti tra righe) e $f_B$ (termini di base della parità di colonna). Ciò garantisce che ogni monomio coinvolga qubit all'interno di un sottoinsieme piccolo e geometricamente locale della griglia, eliminando la necessità di ancille in movimento utilizzate in lavori precedenti [21].
  • Fusione in Pipa: Il circuito è schedulato per fondere più pezzi algebrici in un singolo "passaggio" per riga. Seguendo le porte di interazione dietro una cascata CNOT a somma prefissa, il protocollo raggiunge una profondità di $8L$ per ogni applicazione di $\Gamma$, una riduzione di circa il 38% rispetto alla profondità di $13L$ della precedente costruzione basata su ancille.

C. Protocollo Completo

L'algoritmo completo (Algoritmo 1) viene eseguito come una pipeline a cinque stadi:

1. Esecuzione dell'ordinamento Riga A (utilizzando FSWAP locali).
2. Applicazione di $\Gamma$.
3. Esecuzione dell'ordinamento Colonna (utilizzando FSWAP verticali nudi).
4. Applicazione di $\Gamma$ (poiché $\Gamma = \Gamma^\dagger$).
5. Esecuzione dell'ordinamento Riga B.

Gli operatori $\Gamma$ incapsulano l'ordinamento di colonna, convertendo l'intera sequenza di scambi nudi in una permutazione fermionica valida.

D. Estensione alle Codifiche BK e Parità

Il protocollo è esteso alle codifiche Bravyi–Kitaev (BK) e Parità tramite trasformazioni di base. Mappando i qubit lungo una curva di riempimento dello spazio di Hilbert, gli autori dimostrano che gli strati di CNOT a rotazione di albero (che convertono tra codifiche) possono essere instradati in profondità $O(\sqrt{N})$ su griglie 2D, preservando l'ottimalità asintotica per tutte e tre le principali codifiche.

3. Contributi Chiave

• Profondità Asintoticamente Ottimale senza Ancille: Il documento presenta un algoritmo di permutazione fermionica su griglie NN 2D con profondità $22\sqrt{N} + O(1)$ e $O(N\sqrt{N})$ porte di entanglement tra vicini più prossimi, richiedendo zero ancille, misurazioni a metà circuito o feedforward classico. Ciò corrisponde al limite inferiore $\Omega(\sqrt{N})$ anche rispetto a protocolli con risorse strettamente superiori.
• Costruzione di $\Gamma$ senza Ancille: Una nuova ottimizzazione basata su polinomi di fase che elimina le ancille in movimento $\Theta(\sqrt{N})$ richieste dai precedenti metodi di instradamento fermionico 2D [21], riducendo il prefattore costante della profondità di $\Gamma$ di circa il 38%.
• Agnosticismo rispetto alla Codifica: Uno schema di profondità $O(\sqrt{N})$ per la conversione tra codifiche JW, BK e Parità su griglie 2D utilizzando layout a curva di Hilbert, estendendo il risultato di permutazione ottimale a tutte e tre le codifiche.
• Prestazioni Pratiche: I benchmark numerici dimostrano che il metodo supera costantemente le linee di base esistenti (reti FSWAP 1D e metodi 2D basati su ancille) per dimensioni di sistema $N \gtrsim 100$ in termini di profondità CNOT, volume spazio-temporale e fedeltà stimata.

4. Risultati e Valutazione

Gli autori hanno valutato il protocollo su tre carichi di lavoro: permutazioni fermioniche autonome, la Trasformata di Fourier Fermionica Veloce 2D (FFFT) e la simulazione di Trotter del modello Sachdev–Ye–Kitaev (SYK) sparso.

• Scalabilità: Per $N \ge 100$, la scalabilità $O(\sqrt{N})$ del metodo proposto offre un miglioramento quadratico rispetto alla scalabilità $O(N)$ delle linee di base 1D.
• Riduzione della Profondità: A $N \approx 900$, il metodo proposto achieve una riduzione della profondità di circa il 50% rispetto alla migliore linea di base 2D basata su ancille e di circa il 60% rispetto alla linea di base 1D.
• Fedeltà: Nelle simulazioni di rumore (rumore depolarizzante), il metodo proposto mantiene una fedeltà superiore rispetto ai concorrenti per $N \gtrsim 100$. Ad esempio, a un tasso di errore a due qubit di $10^{-5}$, il metodo mantiene una fedeltà $>0.5$ fino a $N \approx 900$, mentre la linea di base 1D scende sotto questa soglia a $N \approx 400$.
• Efficienza delle Risorse: Eliminando le ancille, il volume spazio-temporale è significativamente ridotto (ad es. circa 1,7 volte inferiore rispetto al metodo basato su ancille a $N=900$).

5. Significato e Affermazioni

Il documento afferma di colmare il divario tra limiti inferiori teorici e compilazione pratica per la simulazione fermionica su hardware 2D. Raggiungendo il limite inferiore $\Omega(\sqrt{N})$ senza ancille, il protocollo è ottimale anche rispetto al potente modello LOCC.

Gli autori sottolineano che questo lavoro è particolarmente rilevante per il regime di tolleranza agli errori iniziale, dove i vincoli di risorse (numero di qubit, tassi di errore) rendono critica l'eliminazione delle ancille e la riduzione dei prefattori costanti. Il metodo consente una simulazione più efficiente di sistemi fermionici (ad es. modelli di Fermi–Hubbard, SYK) e trasformazioni (FFFT) su hardware a breve termine, riducendo direttamente il budget logico richiesto per questi algoritmi.

Il documento conclude notando due questioni aperte: se il limite $O(\sqrt{N})$ si estenda all'intera classe di codifiche ternarie ad albero che preservano il prodotto su griglie 2D, e se il prefattore costante dell'operatore $\Gamma$ possa essere ulteriormente ridotto attraverso un'ottimizzazione congiunta della schedulazione.