Clifford symmetries in quantum many-body systems

Sintesi Tecnica: Simmetrie di Clifford nei Sistemi a Molti Corpi Quantistici

Enunciato del Problema
L'identificazione delle simmetrie negli Hamiltoniani quantistici a molti corpi è un passaggio critico per la comprensione analitica, la determinazione dell'energia dello stato fondamentale e la simulazione efficiente. Sebbene le simmetrie possano portare alla blocchi-diagonalizzazione dell'Hamiltoniano e a significativi risparmi computazionali, trovarle è generalmente un problema intrattabile che spesso dipende dall'intuizione umana ("pensare molto intensamente"). Gli approcci algoritmici esistenti sono limitati alla ricerca di simmetrie di Pauli (operatori che commutano con l'Hamiltoniano), che, sebbene utili, sono spesso non locali e non forniscono intuizioni fondamentalmente nuove su modelli fisici complessi. Inoltre, i metodi capaci di trovare simmetrie più complesse, non di Pauli, richiedono tipicamente la costruzione esplicita dell'Hamiltoniano nello spazio di Hilbert o la risoluzione di problemi di ottimizzazione impegnativi, portando a una scalabilità proibitiva che ne limita l'applicazione a sistemi con solo dozzine di qubit. Attualmente non esiste un approccio algoritmico capace di trovare in modo affidabile simmetrie non di Pauli in sistemi su larga scala (da centinaia a migliaia di qubit).

Metodologia
Gli autori propongono un quadro algoritmico per trovare e sfruttare le simmetrie di Clifford per Hamiltoniani a molti corpi arbitrari. Il metodo sfrutta l'efficienza classica del gruppo di Clifford e la teoria dei grafi. Il flusso di lavoro procede in tre fasi principali:

1. Rappresentazione Grafica tramite Formalismo Simplettico:
   L'Hamiltoniano in ingresso $\hat{H}$, espresso come somma di stringhe di Pauli, viene mappato in una rappresentazione a tableau utilizzando il formalismo simplettico. In questa rappresentazione, le stringhe di Pauli sono trattate come vettori $2n$-dimensionali su $\mathbb{Z}_2$ anziché come matrici $2^n \times 2^n$. Gli autori costruiscono un grafo in cui:

   • I vertici rappresentano i termini di Pauli nell'Hamiltoniano.
   • I colori dei vertici codificano i coefficienti dei termini di Pauli.
   • Gli spigoli collegano i vertici se le corrispondenti stringhe di Pauli anticommutano (il prodotto simplettico è 1).
   • I vertici ausiliari e gli spigoli tratteggiati rappresentano le dipendenze lineari (circuiti) tra le stringhe di Pauli.

   Questa costruzione garantisce che gli Automorfismi del Grafo (GA) del grafo corrispondano direttamente alle simmetrie di Clifford dell'Hamiltoniano. Poiché i GA possono essere risolti in tempo quasi-polinomiale utilizzando strategie di ricerca euristica, ciò consente l'identificazione di simmetrie in sistemi fino a 1.000 qubit.

2. Minimizzazione del Costo in Qubit:
   Trovare una simmetria $\hat{S}$ non è sufficiente per la simulazione; essa deve essere sfruttata diagonalizzandola. Tuttavia, la diagonalizzazione di circuiti di Clifford generali è NP-difficile. Per affrontare questo problema, l'algoritmo trova una trasformazione di base (un operatore di Clifford $\hat{B}$) che mappa la simmetria $\hat{S}$ in una forma minima $\hat{S}_{\text{min}}$. Questo minimizza il "costo in qubit" $Q(\hat{S})$, definito come il numero massimo di qubit su cui agisce un singolo fattore tensoriale della simmetria. Riducendo $Q(\hat{S})$, la simmetria diventa concentrata su un numero minimo di qubit interagenti, rendendo fattibile la diagonalizzazione numerica.

3. Sfruttamento tramite Decomposizione a Blocchi:
   Una volta ottenuta e diagonalizzata $\hat{S}_{\text{min}}$, gli autori decompongono l'Hamiltoniano trasformato $\hat{H}_{\text{sym}} = \hat{B}^\dagger \hat{H} \hat{B}$ in Hamiltoniani efficaci $\hat{H}_\alpha$ che agiscono su specifici sottosettori di simmetria. Questi sottosettori sono definiti dagli autovalori della simmetria. Gli Hamiltoniani efficaci sono derivati numericamente senza mai costruire le matrici complete $2^n \times 2^n$, affidandosi invece alla dimensionalità ridotta determinata dal costo in qubit e dalle degenerazioni degli autovalori.

Contributi Chiave

• Scoperta Algoritmica delle Simmetrie di Clifford: Il documento introduce il primo approccio algoritmico capace di trovare simmetrie non di Pauli (di Clifford) in sistemi quantistici su larga scala (fino a $n=1000$ qubit), colmando un vuoto lasciato dai metodi precedenti limitati alle simmetrie di Pauli o a piccole dimensioni di sistema.
• Mappatura Teorico-Graphica: Il lavoro stabilisce una mappatura rigorosa tra la ricerca di simmetrie di Clifford e il problema dell'Automorfismo del Grafo, consentendo l'uso di risolutori efficienti in tempo quasi-polinomiale.
• Minimizzazione del Costo in Qubit: Gli autori sviluppano un metodo per minimizzare il costo in qubit di una simmetria, garantendo che la successiva diagonalizzazione e lo sfruttamento della simmetria rimangano trattabili classicamente.
• Decomposizione a Blocchi Automatizzata: Il quadro genera automaticamente gli operatori efficaci che descrivono ogni blocco dell'Hamiltoniano decomposto, facilitando l'applicazione diretta alle attività di simulazione.

Risultati
Il metodo è stato testato su Hamiltoniani casuali con simmetrie iniettate e su vari modelli fisici, inclusi scale di Fermi-Hubbard, catene disordinate di Fermioni $tV$, scale di Ising in campo trasverso (TFI) e modelli Heisenberg XXZ.

• Scalabilità: L'algoritmo ha identificato con successo simmetrie di Clifford per istanze fino a 1.000 qubit. Il costo computazionale è dominato dalla ricerca della matrice del prodotto simplettico, che scala quadraticamente con il numero di termini di Pauli ($M$).
• Riconoscimento di Pattern: Per la scala TFI, l'algoritmo ha identificato simmetrie che seguivano un pattern chiaro, permettendo la costruzione induttiva di una simmetria valida per dimensioni di sistema arbitrarie.
• Efficienza di Simulazione: Nel modello Heisenberg XXZ, lo sfruttamento delle simmetrie trovate ha permesso la determinazione degli stati fondamentali e dell'evoluzione temporale di stati casuali di Haar.
  • Per $n=24$, è stato raggiunto un aumento di velocità di circa $\times 2$ per la ricerca dello stato fondamentale.
  • Per $n=20$, è stato raggiunto un aumento di velocità di circa $\times 4$ per l'evoluzione temporale.
  • Crucialmente, il metodo ha permesso alla Diagonalizzazione Esatta (ED) di risolvere sistemi 256 volte più grandi di quanto sarebbe stato possibile senza lo sfruttamento delle simmetrie, limitati solo dalla RAM disponibile (64 GB).
  • Se lo stato iniziale fosse confinato in un singolo sottospazio, i potenziali aumenti di velocità potrebbero variare da $\times 40$ a $\times 6 \cdot 10^4$.

Significato e Affermazioni
Gli autori affermano che questo lavoro rappresenta un passo importante verso la ricerca di simmetrie completamente automatizzata nella fisica a molti corpi. Completando l'ingegno umano con un algoritmo, il metodo fornisce una comprensione analitica più profonda dei modelli fisici e abilita la simulazione di sistemi precedentemente inaccessibili ai metodi classici. L'approccio non si limita alla ricerca di simmetrie, ma ne garantisce anche l'sfruttabilità minimizzando le risorse computazionali richieste per la diagonalizzazione. Il documento suggerisce che questo quadro può avanzare sia le tecniche di simulazione classica (fornendo ansatz ispirati alle simmetrie per risolutori come DMRG o reti tensoriali) sia gli algoritmi quantistici (ottimizzando gli eigensolver variazionali). Gli autori notano che le direzioni future includono il riconoscimento automatico di pattern per intuizioni analitiche, estensioni a simmetrie non di Clifford e chirali, e ulteriore ottimizzazione della decomposizione in sottosettori disgiunti.