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

Dit artikel introduceert een symmetrie-geadapteerde qubit-codering met volledige actieve ruimte (SAE-CAS) die benaderde Z-symmetrieën en de Bravyi-Kitaev-mapping integreert om het aantal qubits en de circuitcomplexiteit voor kwantumchemische simulaties aanzienlijk te verminderen, waarbij superieure convergentie en efficiëntie in middelen worden aangetoond ten opzichte van standaardmethoden op zowel nabije als fouttolerante kwantumprocessors.

De kern

Stel je voor dat je probeert een enorme, ongelooflijk complexe legpuzzel op te lossen om te begrijpen hoe een molecuul werkt. In de wereld van de kwantumchemie is deze puzzel de "elektronische structuur" van een molecuul. Om deze puzzel op een kwantumcomputer op te lossen, moeten we meestal een klein stukje van de computer (een "qubit") toewijzen aan elke mogelijke plek waar een elektron zich zou kunnen bevinden.

Het probleem? Zelfs voor kleine moleculen vereist dit duizenden puzzelstukjes (qubits), en de instructies om ze in elkaar te zetten (het circuit) worden zo lang en verstrengeld dat huidige computers het niet aankunnen, en zelfs toekomstige computers zullen er moeite mee hebben.

Dit artikel introduceert een slimme nieuwe manier om de puzzel op te lossen, genaamd SAE-CAS. Zo werkt het, met behulp van eenvoudige analogieën:

1. De "Bevries en Negeer"-strategie (Het CAS-gedeelte)

Denk aan een molecuul als een druk kantoorgebouw.

• De Bevroren Kern: De kelder en de bovenste verdieping zijn altijd vol met mensen die nooit vertrekken en ook niet interageren met de rest van het gebouw. In kwantümtermen zijn dit "frozen-core" elektronen. Ze zijn saai en voorspelbaar.
• De Virtuele'n: De zolder is volledig leeg en zal waarschijnlijk ook leeg blijven. Dit zijn "virtuele" orbitalen.
• De Actieve Ruimte: De middelste verdiepingen zijn waar de echte actie plaatsvindt. Mensen bewegen rond, praten en veranderen dingen. Dit is de "actieve ruimte".

Traditionele methoden proberen een qubit toe te wijzen aan elke verdieping, zelfs aan de saaie kelder en de lege zolder. SAE-CAS zegt: "Laten we de kelder en de zolder gewoon negeren." We wijzen alleen qubits toe aan de middelste verdiepingen waar de interessante chemie gebeurt. Dit verkleint direct de omvang van de puzzel die we moeten oplossen.

2. De "Symmetrie-Snelkoppeling" (Het SAE-gedeelte)

Zelfs binnen de drukke middelste verdiepingen zijn er regels. In een watermolecuul is bijvoorbeeld de linkerkant een spiegelbeeld van de rechterkant. Als je weet wat er aan de linkerkant gebeurt, weet je automatisch ook wat er aan de rechterkant gebeurt.

Normaal gesproken berekenen computers beide kanten apart, wat tijd en middelen verspilt. SAE-CAS gebruikt een wiskundige "truc" (een affine Clifford-transformatie) om te beseffen dat, omdat deze spiegelregels bestaan, we niet voor beide kanten aparte qubits nodig hebben. We kunnen de puzzel wiskundig "dubbelvouwen". Dit verwijdert nog meer qubits, waardoor de puzzel nog kleiner en gemakkelijker op te lossen wordt.

3. De "Vertaling" (Het Bravyi–Kitaev-gedeelte)

Zodra we onze kleine, dubbelgevouwen puzzel hebben, moeten we deze vertalen naar een taal die de kwantumcomputer begrijpt. Er zijn twee belangrijke vertalers:

• Jordan-Wigner (JW): De standaardvertaler. Het is simpel, maar maakt de instructies erg lang (als het lezen van een boek waarbij elk woord wordt herhaald).
• Bravyi–Kitaev (BK): Een slimmere vertaler. Deze organiseert de informatie efficiënter, waardoor de instructies korter en minder verstrengeld zijn.

De auteurs laten zien dat je de "Dubbelgevouwen Puzzel"-methode (SAE-CAS) met elke vertaler kunt gebruiken. Ze hebben een versie gemaakt genaamd SAE-CAS-BK die de slimmere vertaler gebruikt. Het verandert het uiteindelijke antwoord niet, maar maakt het pad naar het resultaat vaak soepeler en sneller.

Wat hebben ze gevonden?

De auteurs hebben deze methode getest op negen kleine moleculen (zoals water, zuurstof en stikstof) met behulp van twee verschillende "zoekstrategieën" (algoritmen) om de energie van het molecuul te vinden:

1. UCCSD: Een chemisch precieze maar complexe zoektocht.
2. HE-SCA: Een eenvoudigere, hardware-vriendelijke zoektocht.

De Resultaten:

• Minder Qubits: Door de saaie onderdelen te negeren en de symmetrische delen te vouwen, hadden ze aanzienlijk minder qubits nodig (soms de helft of meer verminderd).
• Kortere Circuits: De instructies om de simulatie uit te voeren waren veel korter en minder verstrengeld.
• Sneller Succes: Bij het gebruik van de eenvoudigere zoekstrategie (HE-SCA) vond hun methode het juiste antwoord voor elk getest molecuul. De oude methode (JW-CAS) kwam vast te zitten en slaagde er niet in om voor zuurstof en koolmonoxide het antwoord te vinden binnen de tijdslimieten.
• Geen Verlies van Nauwkeurigheid: Ondanks dat ze de "saaie" elektronen negeerden en de puzzel dubbelvouwden, waren de uiteindelijke energiewaarden net zo nauwkeurig als de standaard, massale berekeningen.

De Kernboodschap

De auteurs hebben een "efficiënte toolkit" gebouwd. Ze hebben bewezen dat je veilig de delen van het molecuul die niet veranderen kunt weggooien en de symmetrische delen kunt vouwen, zonder het juiste antwoord te verliezen. Dit maakt het mogelijk om deze complexe chemische simulaties uit te voeren op kwantumcomputers die veel kleiner en minder krachtig zijn dan voorheen nodig werd geacht.

Ze hebben ook de code voor deze "magische truc" gratis beschikbaar gesteld (in een pakket genaamd QuantumSymmetry), zodat anderen deze kunnen gebruiken om moleculen op hun eigen kwantumcomputers te simuleren.

Technische samenvatting

Technische Samenvatting: Symmetrie-adaptieve Qubit-codering met Complete Active Space en Bravyi–Kitaev Mapping

Probleemstelling
Simulaties van kwantumchemie op nabije (NISQ) en fouttolerante kwantumprocessors worden geconfronteerd met aanzienlijke beperkingen in middelen, die primair worden gedreven door het aantal qubits dat nodig is om de elektronische structuren van moleculen te coderen en de complexiteit (Pauli-gewicht en circuitdiepte) van de resulterende Hamiltonia. Hoewel de Complete Active Space (CAS) benadering de omvang van het probleem verkleint door kernorbitalen te bevriezen en virtuele orbitalen te verwerpen, schieten standaard implementaties er vaak naast als ze niet volledig gebruikmaken van symmetrieën om het aantal qubits verder te minimaliseren. Daarnaast bieden verschillende fermion-naar-qubit mappings, zoals Jordan–Wigner (JW) en Bravyi–Kitaev (BK), afruilopties tussen lokaliteit en qubit-vereisten die niet altijd optimaal geïntegreerd zijn met symmetrie-reductietechnieken.

Methodologie
De auteurs stellen een verenigd framework voor genaamd Symmetry-Adapted Encoding with Complete Active Space (SAE-CAS). Deze aanpak behandelt de CAS-benadering niet louter als een inkorting van de Hilbertruimte, maar als de oplegging van benaderde $Z$-symmetrieën op bevroren-kern en virtuele orbitalen.

1. Theoretisch Kader:

   • Benaderde Symmetrieën: De methode breidt exacte symmetrie-adaptieve coderingen (SAEs) uit naar benaderde symmetrieën die overeenkomen met vaste orbitaalbezettingen (bevroren-kern en virtueel). Deze worden gemodelleerd als $Z$-symmetrieën waarbij specifieke qubits worden beperkt tot eigenwaarden van $-1$ (bevroren) of $+1$ (virtueel).
   • Affine Clifford Transformaties: De reductie wordt bereikt via affine Clifford-kaarten. De auteurs definiëren een binaire orbitaal-karakter matrix $S$ die de symmetrieën representeert. Door een inverteerbare affine kaart $T$ over $\mathbb{F}_2$ te construeren, transformeren zij de computationele basis zodanig dat de symmetriebeperkingen overeenkomen met vaste waarden (meestal 0) op specifieke qubits.
   • Qubit Verwijdering: Zodra het Hamiltonian geconjugeerd is door de Clifford-operator, worden termen die $X$ of $Y$ bevatten op de beperkte qubits weggegooid, en de beperkte $Z$-operatoren worden vervangen door scalairen. Dit verwijdert permanent de redundante qubits, waardoor het systeem wordt gereduceerd van $n$ spin-orbitalen naar $n-k$ qubits, waarbij $k$ het aantal onafhankelijke symmetrieën is.
   • Equivalentiebewijs: De auteurs bewijzen dat de resulterende actieve-ruimte qubit-Hamiltonian exact equivalent is aan de canonieke CAS-Hamiltonian afgeleid via bevroren-kern en virtuele-orbitaal projectie in de moleculaire elektronische structuurtheorie. De benaderingsfout komt uitsluitend voort uit de keuze van de actieve ruimte, niet uit de codering.

2. Integratie met Bravyi–Kitaev (SAE-CAS-BK):

   • Het framework wordt uitgebreid naar de Bravyi–Kitaev (BK) mapping. Aangezien de BK-transformatie zelf een affine Clifford-basisverandering is die werkt op de computationele bitstrings, kan deze samengesteld worden met de SAE-CAS Clifford-kaart.
   • De resulterende SAE-CAS-BK codering is unitair equivalent aan SAE-CAS. Het behoudt het aantal qubits, het aantal Pauli-termen en de variatiemparameters, maar verandert de lokaliteitsstructuur (Pauli-gewicht per term) en de resulterende circuitdiepte/CNOT-aantallen.

3. Numerieke Benchmarking:

   • De methode is getest op negen kleine moleculen (waaronder H$_2$O, C$_2$H$_4$, O$_2$, N$_2$, en C$_4$H$_6$) met een minimale STO-3G basis.
   • Twee variatiële ansätze werden toegepast: Unitary Coupled Cluster met Single en Double excitatie (UCCSD) en een hardware-efficiënte shifted-circular-alternating (HE-SCA) ansatz.
   • Vergelijkingen werden gemaakt met standaard JW, JW-CAS (CAS met JW maar zonder symmetrie-reductie), en JW-CAS met symmetrie-filtering (SF) op circuitniveau.

Belangrijkste Resultaten

• Middelenreductie: SAE-CAS vermindert consistent het aantal qubits, het aantal Pauli-operatoren, de circuitdiepte, het CNOT-aantal en de variatiemparameters vergeleken met JW-CAS. Bijvoorbeeld, in de C$_4$H$_4$ (4,4) CAS, reduceerde SAE-CAS het aantal qubits van 8 (JW-CAS) naar 4 en het aantal Pauli-operatoren van 97 naar 46.
• Convergentie en Nauwkeurigheid:
  • UCCSD: SAE-CAS bereikte geconvergeerde energieën vergelijkbaar met JW-CAS en JW-CAS (SF), maar vereiste aanzienlijk minder optimalisatie-iteraties (bijv. 21 versus 109 voor N$_2$). Opvallend genoeg slaagden de JW-CAS en JW-CAS (SF) UCCSD ansätze voor C$_4$H$_4$ en C$_4$H$_6$ er niet in om te draaien vanwege geheugenuitputting door het grote aantal excitaties, terwijl SAE-CAS dit wel slaagde door de excitatie-lijst in de mapping-fase te filteren.
  • HE-SCA: SAE-CAS demonstreerde superieure convergentie-eigenschappen. Voor alle geteste moleculen convergeerde SAE-CAS naar de CAS referentie-energieën binnen de geteste laag-budgetten. In tegenstelling hiertoe faalde JW-CAS om te convergeren voor O$_2$ en CO binnen 200 lagen. De auteurs schrijven dit toe aan het "barren plateau" effect en sector-lekken in JW-CAS, waarbij de variatiële manifold naar de verkeerde spin-pariteit sectoren drift, terwijl SAE-CAS de zoektocht inherent beperkt tot de juiste symmetrie-sector.
• SAE-CAS vs. SAE-CAS-BK: De twee coderingen produceerden identieke geconvergeerde energieën en iteratie-aantallen, wat hun unitaire equivalentie bevestigt. Verschillen werden alleen waargenomen in circuitdiepte en CNOT-aantallen, die voor de meeste moleculen met minder dan 6% varieerden (behalve voor CO, waar een toename van +20% werd gezien voor BK). De auteurs merken op dat het asymptotische $O(\log n)$ lokaliteitsvoordeel van BK nog niet is gerealiseerd bij deze kleine actieve ruimtes.

Betekenis en Claims
Het artikel claimt dat SAE-CAS een systematische en combineerbare route biedt voor middelen-efficiënte moleculaire simulaties. Door de CAS-benadering te behandelen als een set benaderde symmetrieën, verenigt de methode exacte symmetrie-reductie (puntgroep- en spin-pariteit) met actieve-ruimte inkorting. De primaire betekenis ligt in:

1. Natuurlijke Symmetrie-afdwinging: Het afdwingen van symmetrieën op het niveau van de codering (in plaats van post-hoc in de ansatz) verwijdert redundante zoekrichtingen en voorkomt dat de optimizer afdrijft naar onjuiste deeltjesgetal- of spin-sectoren.
2. Schaalbaarheid: De methode maakt simulaties mogelijk op grotere actieve ruimtes die anders onhandelbaar zouden zijn vanwege geheugenbeperkingen in het construeren van de ansatz (zoals gezien in de C$_4$H$_4$ en C$_4$H$_6$ gevallen).
3. Flexibiliteit: Het framework is compatibel met elke fermion-naar-qubit mapping die uitgedrukt kan worden als een affine Clifford-basisverandering, zoals aangetoond met zowel JW als BK.

De auteurs bieden een open-source implementatie van de methode aan in het Python-pakket QuantumSymmetry en concluderen dat SAE-CAS een levensvatbare route is voor het reduceren van qubit-aantallen en operator-complexiteit op zowel nabije als fouttolerante kwantumprocessors.