Exact chiral symmetry with quantum signal processing

Sintesi Tecnica: Simmetria Chirale Esatta con il Quantum Signal Processing

Enunciato del Problema
La simulazione quantistica di osservabili non perturbative e fuori equilibrio nella cromodinamica quantistica (QCD) su reticolo affronta sfide significative riguardanti il raddoppio dei fermioni e la preservazione della simmetria chirale nelle formulazioni hamiltoniane. Sebbene la relazione di Ginsparg-Wilson (GW) e l'operatore di overlap offrano soluzioni eleganti nello spaziotempo euclideo, un analogo hamiltoniano canonico per i fermioni di overlap è stato finora mancante. Gli approcci esistenti, come i fermioni a parete di dominio (domain-wall fermions), si mappano bene nelle impostazioni hamiltoniane ma richiedono una dimensione extra esplicita, aumentando i costi in termini di qubit. Al contrario, i fermioni di overlap offrono una formulazione nelle dimensioni spaziali fisiche, ma comportano interazioni altamente non locali all-to-all, specificamente la funzione segno dell'Hamiltoniana di Dirac-Wilson, $\varepsilon(h_W)$, che è difficile da implementare efficientemente su hardware quantistico. Il problema centrale affrontato è come costruire algoritmi quantistici efficienti per i fermioni di reticolo con simmetria chirale che bilancino i costi di memoria (qubit) rispetto alla complessità delle porte, mantenendo una simmetria chirale esatta (o controllatamente rotta).

Metodologia
Gli autori propongono un algoritmo di Quantum Signal Processing (QSP) per simulare l'Hamiltoniana dei fermioni di overlap. La metodologia procede attraverso fasi chiave:

1. Formulazione dell'Hamiltoniana: Il lavoro utilizza l'Hamiltoniana di overlap $H_{ov} = \psi^\dagger (\Gamma_0 + \varepsilon(h_W)) \psi$, dove $h_W$ è l'Hamiltoniana single-particle di Dirac-Wilson. La funzione segno $\varepsilon(X) = X/\sqrt{X^\dagger X}$ è approssimata da un polinomio $p_M(x)$ di grado $M$.
2. Block Encoding: Gli autori costruiscono un block encoding dell'Hamiltoniana di Dirac-Wilson single-particle $\hat{h}_W$ (inclusi i campi di gauge) utilizzando un operatore "prepare" $P$ e un operatore "select" $S$. Questo codifica l'Hamiltoniana in un operatore unitario $U_{\hat{h}_W}$ che agisce su uno spazio di Hilbert esteso con qubit ancillari.
3. Quantum Signal Processing (QSP): Per implementare l'approssimazione della funzione segno, gli autori impiegano il QSP (specificamente la Quantum Singular Value Transformation). Ciò consente l'applicazione di un polinomio di grado $M$, $p_M(U_{\hat{h}_W})$, all'operatore block-encoded. Il polinomio è scelto per approssimare la funzione segno all'interno di un gap spettrale $[\kappa, 1]$ con un errore $\epsilon_e$.
4. Evoluzione Temporale: Una volta che l'Hamiltoniana di overlap è stata block-codificata, l'operatore di evoluzione temporale $e^{-iHt}$ viene approssimato utilizzando il QSP, richiedendo un numero di porte che scala con il costo del block-encoding e il tempo di evoluzione.

Contributi Chiave e Risultati

• Costruzione Algoritmica: Il documento fornisce un algoritmo concreto basato su QSP per l'Hamiltoniana di overlap che preserva la relazione di GW fino a un errore controllabile $\epsilon_e$. Si dimostra che l'operatore chirale modificato $\hat{\gamma}_5$ commuta con l'Hamiltoniana approssimata fino a un errore di ordine $\epsilon_e$.
• Analisi della Complessità:
  • Complessità delle Porte: Il costo per il block-encoding dell'Hamiltoniana di overlap scala come $O(Q \kappa^{-1} \log(1/\epsilon_e))$, dove $Q$ è il numero di siti del reticolo (moltiplicato per i gradi di libertà interni) e $\kappa$ è legato al gap spettrale. La complessità totale delle porte per l'evoluzione temporale scala come $O(Q \kappa^{-1} \log(1/\epsilon_e) [Qt + \log(1/\epsilon_t)])$.
  • Costo dei Qubit: L'algoritmo richiede $O(Q + \log Q)$ qubit. Questa è una riduzione significativa rispetto ai fermioni a parete di dominio, che richiedono $O(Q L_5)$ qubit, dove $L_5$ è l'estensione della dimensione extra.
  • Confronto: Mentre i fermioni a parete di dominio beneficiano della località geometrica (permettendo un costo quasi lineare $O(Q t \text{polylog})$ tramite Suzuki-Trotter o metodi simili), la formulazione di overlap incorre in una profondità di porte più elevata a causa della natura non locale dell'approssimazione della funzione segno. Tuttavia, l'approccio di overlap offre una scalabilità asintotica superiore in termini di memoria.
• Interpretazione Fisica del QSP: Gli autori dimostrano che il grado del polinomio $M$ richiesto per ottenere un errore $\epsilon_e$ scala come $M = O(\kappa^{-1} \log(1/\epsilon_e))$. Identificano questa scalabilità con la dimensione della dimensione extra $L_5$ nei fermioni a parete di dominio, dove l'errore scala come $e^{-c L_5}$. Pertanto, l'implementazione QSP "costruisce" efficacemente una dimensione extra attraverso la profondità del circuito, rispecchiando la corrispondenza fisica tra l'operatore di overlap e la teoria di bordo dei fermioni a parete di dominio.

Significatività e Rivendicazioni
Il documento afferma che le simulazioni quantistiche di fermioni di Dirac con simmetria chirale esatta sono "quasi gratuite" nel senso che l'applicazione dell'Hamiltoniana di overlap costa solo un fattore logaritmico in più (in termini di tolleranza dell'errore) rispetto all'Hamiltoniana di Dirac-Wilson.

L'intuizione centrale è che il QSP fornisce una realizzazione algoritmica quantistica della nota corrispondenza tra l'operatore di overlap e la dimensione extra dei fermioni a parete di dominio. Il compromesso è esplicito:

• Fermioni a parete di dominio: Maggiore costo di qubit ($O(Q L_5)$) ma minore profondità di porte grazie alla località.
• Fermioni di overlap (tramite QSP): Minore costo di qubit ($O(Q)$) ma maggiore profondità di porte ($O(Q^2)$ nel caso peggiore per l'evoluzione temporale a causa della non località).

Gli autori concludono che la scelta tra le formulazioni dipende dai vincoli dell'hardware (numero di qubit rispetto alla profondità del circuito). Notano che, sebbene la QCD classica su reticolo utilizzi spesso approssimazioni razionali (ad esempio, Zolotarev) per la funzione segno che sono più efficienti, queste non hanno analoghi diretti nel QSP standard, suggerendo una direzione per lavori futuri. Il documento chiarisce inoltre che i risultati di scalabilità sono robusti rispetto alla specifica codifica dei campi di gauge, sebbene i fattori costanti possano variare.