Point-group symmetry analysis of many-electron wavefunctions on a quantum computer

Sintesi Tecnica: Analisi della simmetria di gruppo puntuale delle funzioni d'onda a molti elettroni su un computer quantistico

Enunciato del Problema
I gruppi puntuali, che rappresentano le operazioni di simmetria spaziale nei sistemi molecolari, sono strumenti fondamentali in chimica per l'analisi degli orbitali molecolari e della spettroscopia. Sebbene siano stati proposti algoritmi quantistici per sfruttare queste simmetrie—ad esempio per ridurre il numero di qubit o semplificare gli ansatz dei circuiti—le implementazioni pratiche delle operazioni di simmetria di gruppo puntuale e le analisi dettagliate delle simmetrie di funzioni d'onda realistiche a molti elettroni rimangono in gran parte inesplorate. Nello specifico, manca un metodo per implementare concretamente le operazioni di simmetria e analizzare le proprietà di simmetria (ad esempio, i pesi delle rappresentazioni irriducibili) di stati reali e correlati a molti elettroni su hardware quantistico a breve termine.

Metodologia
Gli autori propongono un metodo ibrido quantistico-classico senza ancilla per analizzare le simmetrie di gruppo puntuale degli stati a molti elettroni, applicabile sia a gruppi abeliani che non abeliani. Il cuore del metodo consiste nel calcolare il "peso" ($w_\Gamma$) delle rappresentazioni irriducibili ($\Gamma$) all'interno di una data funzione d'onda $|\Psi\rangle$. Tale peso è definito come il coefficiente al quadrato della proiezione di $|\Psi\rangle$ sullo sottomodo di $\Gamma$.

Il quadro teorico si basa sulla formula dell'operatore di proiezione:
$$w_\Gamma = \frac{d_\Gamma}{|G|} \sum_C \chi^*_\Gamma(C) \sum_{g \in C} \langle\Psi|\hat{g}|\Psi\rangle$$
dove $d_\Gamma$ è la dimensione della rappresentazione, $|G|$ è l'ordine del gruppo, $\chi_\Gamma(C)$ è il carattere della classe di coniugazione $C$ e $\hat{g}$ è l'operazione di simmetria.

Per valutare $\langle\Psi|\hat{g}|\Psi\rangle$ su un computer quantistico senza qubit ancilla, gli autori utilizzano la diagonalizzazione delle matrici di rappresentazione. Per orbitali adattati alla simmetria (ad esempio, Hartree-Fock), l'operazione di simmetria $U(g)$ viene decomposta in rotazioni orbitali. Diagonalizzando la matrice di rappresentazione $\bar{D}_\gamma(g)$ per insiemi di orbitali degeneri, l'unitaria $U(g)$ viene riscritta come prodotto di una rotazione di base $\tilde{U}(g)$ e di un operatore di fase diagonale. Il valore di aspettazione viene quindi ottenuto preparando lo stato $|\tilde{\Psi}(g)\rangle = \tilde{U}(g)|\Psi\rangle$, misurando nella base computazionale e combinando i risultati con i fattori di fase appropriati. Questo approccio evita i circuiti profondi e i requisiti di ancilla dei test di Hadamard o delle tecniche di Combinazione Lineare di Unitari (LCU).

Lo studio impiega due strategie principali di validazione:

1. Simulazioni Numeriche Classiche: Utilizzando le librerie ffsim e pyscf, il metodo viene applicato al benzene ($D_{6h}$) e al ferrocene sfalsato ($D_{5d}$). Gli autori analizzano singoli determinanti di Slater (stati fondamentali e stati a singola eccitazione di Hartree-Fock) e funzioni d'onda correlate (funzioni Cluster Unitario Jastrow e Cluster Unitario Jastrow Locale).
2. Dimostrazione Hardware: Il metodo viene eseguito sul dispositivo quantistico superconduttivo ibm_kawasaki di IBM (fino a 32 qubit). Il sistema target è il benzene nella simmetria $D_{2h}$ (un sottogruppo di $D_{6h}$). Le funzioni d'onda di prova (stato fondamentale e primo stato eccitato) vengono preparate utilizzando calcoli del Gruppo di Rinormalizzazione della Matrice Densità (DMRG), compresse in un ansatz a muro di mattoni tramite tecniche di rete tensoriale e mappate su circuiti quantistici. Per affrontare il rumore, vengono confrontate due tecniche di mitigazione degli errori: l'Estrazione a Rumore Zero (ZNE) tramite folding dei gate e un approccio di mitigazione degli errori quasi-probabilistico (QESEM).

Risultati Chiave

• Simulazioni Classiche:
  • Per gli stati fondamentali Hartree-Fock del benzene e del ferrocene, il metodo identifica correttamente il peso della rappresentazione irriducibile totalmente simmetrica ($A_{1g}$) come prossimo a 1.0, coerentemente con le configurazioni a guscio chiuso.
  • Per gli stati a singola eccitazione (SE), il metodo scompone con successo le funzioni d'onda in multiple rappresentazioni irriducibili. Ad esempio, un specifico stato SE nel benzene è stato scomposto in $\approx 0.005 B_{1u} + 0.495 B_{2u} + 0.5 E_{1u}$. La somma dei pesi su tutte le configurazioni SE indipendenti ha prodotto risultati invarianti coerenti con le formule di riduzione della teoria dei gruppi.
  • Nell'analisi delle funzioni d'onda correlate LUCJ, il metodo ha rivelato che le procedure di ottimizzazione che minimizzano l'errore energetico (senza regolarizzazione) possono degradare significativamente il peso di simmetria (deviando dal corretto peso $A_{1g}$), mentre la regolarizzazione aiuta a preservare le proprietà di simmetria, sebbene con diversi compromessi energetici. Ciò suggerisce che la metrica è utile per valutare la qualità delle funzioni d'onda di prova.
• Dimostrazione Hardware:
  • Sul dispositivo ibm_kawasaki, le misurazioni grezze hanno mostrato un significativo degrado dei pesi di simmetria all'aumentare del numero di qubit (ad esempio, il peso $A_g$ per lo stato fondamentale è sceso a ~0.54 per 32 qubit senza mitigazione).
  • L'applicazione della mitigazione degli errori è stata cruciale. Sebbene la ZNE abbia migliorato i risultati, l'approccio quasi-probabilistico QESEM ha prodotto la massima fedeltà, riproducendo i pesi corretti delle rappresentazioni irriducibili con un errore di pochi punti percentuali (ad esempio, $0.968 \pm 0.010$ per lo stato fondamentale su 32 qubit, rispetto a un valore esatto di $0.999$).
  • L'esperimento ha confermato che il metodo proposto senza ancilla, combinato con la preparazione dello stato tramite rete tensoriale e una mitigazione avanzata degli errori, può riprodurre fedelmente le proprietà di simmetria sull'hardware attuale.

Significato e Affermazioni
Il lavoro si propone di colmare il divario tra la proiezione teorica della simmetria e l'implementazione pratica su hardware quantistico. Il metodo proposto funge da strumento pratico per:

1. Analisi di Simmetria: Analizzare le proprietà di simmetria delle funzioni d'onda a molti elettroni in simulazioni di materiali realistiche su computer quantistici a breve termine e a tolleranza agli errori iniziale.
2. Valutazione della Funzione d'Onda: Valutare la qualità delle funzioni d'onda di prova (ad esempio, negli algoritmi variazionali) quantificando la loro adesione a specifiche simmetrie.
3. Benchmarking: Fornire un benchmark per dispositivi quantistici reali e tecniche di mitigazione degli errori, poiché il metodo non richiede qubit ancilla e comporta un sovraccarico minimo per orbitali adattati alla simmetria.

Gli autori sottolineano che, sebbene la dimostrazione hardware abbia utilizzato un sottogruppo abeliano ($D_{2h}$) e una valutazione basata su Pauli, la metodologia sottostante è generale e applicabile a gruppi non abeliani e funzioni di base arbitrarie. Il lavoro suggerisce che l'analisi di simmetria è un componente vitale e necessario per le future simulazioni quantistiche di materiali complessi.