Symmetry-adapted qubit encoding with complete active space and Bravyi--Kitaev mapping for quantum chemistry on a quantum computer

Sintesi Tecnica: Codifica di Qubit Adattata alla Simmetria con Spazio Attivo Completo e Mappatura di Bravyi–Kitaev

Enunciato del Problema
Le simulazioni di chimica quantistica su processori quantistici a breve termine (NISQ) e tolleranti ai guasti affrontano vincoli di risorse significativi, guidati principalmente dal numero di qubit richiesti per codificare le strutture elettroniche molecolari e dalla complessità (peso di Pauli e profondità del circuito) degli Hamiltoniani risultanti. Sebbene l'approssimazione dello Spazio Attivo Completo (CAS) riduca le dimensioni del problema congelando gli orbitali core e scartando gli orbitali virtuali, le implementazioni standard spesso non riescono a sfruttare appieno le simmetrie per ridurre ulteriormente il conteggio dei qubit. Inoltre, diverse mappature fermione-qubit, come Jordan–Wigner (JW) e Bravyi–Kitaev (BK), offrono compromessi tra località e requisiti di qubit che non sono sempre integrati in modo ottimale con le tecniche di riduzione della simmetria.

Metodologia
Gli autori propongono un framework unificato chiamato Symmetry-Adapted Encoding with Complete Active Space (SAE-CAS). Questo approccio tratta l'approssimazione CAS non semplicemente come una troncatura dello spazio di Hilbert, ma come l'imposizione di simmetrie $Z$ approssimate sugli orbitali core congelati e sugli orbitali virtuali.

1. Framework Teorico:

   • Simmetrie Approssimate: Il metodo estende le codifiche adattate alla simmetria (SAE) esatte per includere simmetrie approssimate corrispondenti a occupazioni orbitaliche fisse (core congelati e virtuali). Queste sono modellate come simmetrie $Z$ dove specifici qubit sono vincolati agli autovalori di $-1$ (congelati) o $+1$ (virtuali).
   • Trasformazioni Clifford Affini: La riduzione è ottenuta tramite mappe Clifford affini. Gli autori definiscono una matrice di carattere orbitale binaria $S$ che rappresenta le simmetrie. Costruendo una mappa affine invertibile $T$ su $\mathbb{F}_2$, trasformano la base computazionale in modo che i vincoli di simmetria corrispondano a valori fissi (tipicamente 0) su specifici qubit.
   • Rimozione dei Qubit: Una volta che l'Hamiltoniana è coniugata dall'operatore Clifford, i termini contenenti $X$ o $Y$ sui qubit vincolati vengono scartati, e gli operatori $Z$ vincolati sono sostituiti da scalari. Questo rimuove permanentemente i qubit ridondanti, riducendo il sistema da $n$ spin-orbitali a $n-k$ qubit, dove $k$ è il numero di simmetrie indipendenti.
   • Prova di Equivalenza: Gli autori dimostrano che l'Hamiltoniana qubit dello spazio attivo risultante è esattamente equivalente all'Hamiltoniana CAS canonica derivata tramite la proiezione degli orbitali core congelati e degli orbitali virtuali nella teoria della struttura elettronica molecolare. L'errore di approssimazione deriva esclusivamente dalla scelta dello spazio attivo, non dalla codifica.

2. Integrazione con Bravyi–Kitaev (SAE-CAS-BK):

   • Il framework viene esteso alla mappatura Bravyi–Kitaev (BK). Poiché la trasformazione BK è essa stessa un cambiamento di base Clifford affine che agisce sulle stringhe di bit computazionali, può essere composta con la mappa Clifford SAE-CAS.
   • La codifica SAE-CAS-BK risultante è unitariamente equivalente a SAE-CAS. Essa preserva il numero di qubit, il numero di termini di Pauli e i parametri variazionali, ma altera la struttura di località (peso di Pauli per termine) e la conseguente profondità del circuito/conteggio dei CNOT.

3. Benchmarking Numerico:

   • Il metodo è stato testato su nove piccole molecole (tra cui H$_2$O, C$_2$H$_4$, O$_2$, N$_2$, e C$_4$H$_6$) utilizzando una base minima STO-3G.
   • Sono stati impiegati due ansätze variazionali: l'Unità di Clusteracco con Singole e Doppie eccitazioni (UCCSD) e un ansatz hardware-efficient shifted-circular-alternating (HE-SCA).
   • I confronti sono stati effettuati contro JW standard, JW-CAS (CAS con JW ma senza riduzione di simmetria) e JW-CAS con filtraggio di simmetria (SF) a livello di circuito.

Risultati Chiave

• Riduzione delle Risorse: SAE-CAS riduce costantemente il numero di qubit, il conteggio degli operatori di Pauli, la profondità del circuito, il conteggio dei CNOT e i parametri variazionali rispetto a JW-CAS. Ad esempio, nel CAS (4,4) di C$_4$H$_4$, SAE-CAS ha ridotto il numero di qubit da 8 (JW-CAS) a 4 e il conteggio di Pauli da 97 a 46.
• Convergenza e Accuratezza:
  • UCCSD: SAE-CAS ha raggiunto energie confluite comparabili a JW-CAS e JW-CAS (SF), ma ha richiesto significativamente meno iterazioni di ottimizzazione (ad esempio, 21 contro 109 per N$_2$). Notevolmente, per C$_4$H$_4$ e C$_4$H$_6$, gli ansatz UCCSD di JW-CAS e JW-CAS (SF) non sono riusciti a girare a causa dell'esaurimento della memoria dovuto all'elevato numero di eccitazioni, mentre SAE-CAS ha avuto successo filtrando la lista delle eccitazioni durante la fase di mappatura.
  • HE-SCA: SAE-CAS ha dimostrato proprietà di convergenza superiori. Per tutte le molecole testate, SAE-CAS è confluito alle energie di riferimento CAS entro i budget di layer testati. Al contrario, JW-CAS non è riuscito a convergere per O$_2$ e CO entro 200 layer. Gli autori attribuiscono ciò all'effetto "barren plateau" e al "sector leakage" in JW-CAS, dove il manifold variazionale abbraccia settori di spin-parità errati, mentre SAE-CAS restringe la ricerca al corretto settore di simmetria nativamente.
  • SAE-CAS vs. SAE-CAS-BK: Le due codifiche hanno prodotto energie confluite e conteggi di iterazione identici, confermando la loro equivalenza unitaria. Differenze sono state osservate solo nella profondità del circuito e nei conteggi dei CNOT, che variavano di meno del 6% per la maggior parte delle molecole (eccetto CO, che ha visto un incremento del +20% per BK). Gli autori osservano che il vantaggio di località asintotica $O(\log n)$ di BK non è ancora stato realizzato in queste piccole dimensioni di spazio attivo.

Significato e Rivendicazioni
L'articolo sostiene che SAE-CAS offre un percorso sistematico e composibile per simulazioni molecolari efficienti dal punto di vista delle risorse. Trattando l'approssimazione CAS come un insieme di simmetrie approssimate, il metodo unifica la riduzione di simmetria esatta (gruppo puntuale e spin-parità) con la troncatura dello spazio attivo. Il significato primario risiede in:

1. Imposizione Nativa della Simmetria: Imporre le simmetrie a livello di codifica (anziché post-hoc nell'ansatz) rimuove direzioni di ricerca ridondanti e impedisce all'ottimizzatore di deviare verso settori di numero di particelle o spin errati.
2. Scalabilità: Il metodo consente simulazioni su spazi attivi più grandi che sarebbero altrimenti intrattabili a causa dei vincoli di memoria nella costruzione dell'ansatz (come visto nei casi C$_4$H$_4$ e C$_4$H$_6$).
3. Flessibilità: Il framework è compatibile con qualsiasi mappatura fermione-qubit esprimibile come un cambiamento di base Clifford affine, come dimostrato qui sia con JW che con BK.

Gli autori forniscono un'implementazione open-source nel pacchetto Python QuantumSymmetry e concludono che SAE-CAS è una via percorribile per ridurre il numero di qubit e la complessità degli operatori sia su processori quantistici a breve termine che tolleranti ai guasti.