Orbital-optimized spin-adapted multistate contracted VQE for excited states and properties on quantum hardware

Dit artikel introduceert de orbital-optimized multistate contracted VQE (oo-MC-VQE) methode, die spin-geadapteerde operatoren gebruikt om efficiënt grondtoestanden en aangeslagen toestanden samen met hun eigenschappen te berekenen op quantumhardware, terwijl de nauwkeurigheid en circuitcomplexiteit worden gebalanceerd door middel van een lineaire parameterschaling met het aantal toestanden.

De kern

Stel je voor dat je een enorm, complex orkest probeert af te stemmen om een perfecte symfonie te spelen. In de wereld van de chemie is het "orkest" een molecuul, en de "muziek" is de manier waarop de elektronen bewegen en interageren. Om te begrijpen hoe een molecuul licht absorbeert (wat kleuren geeft en zaken als fotosynthese aanstuurt), moeten wetenschappers de exacte noten berekenen die deze elektronen spelen.

Lange tijd was het berekenen van dit soort puzzels voor moleculen met veel elektronen alsof je een puzzel probeert op te lossen die exponentieel moeilijker wordt naarmate je meer stukjes toevoegt. Klassieke computers (de computers die we vandaag de dag gebruiken) lopen uiteindelijk tegen een muur aan en kunnen deze puzzels voor complexe moleculen niet oplossen.

Dit artikel introduceert een nieuwe manier om deze puzzels op te lossen met behulp van kwantumcomputers, wat speciale machines zijn die ontworiment zijn om met deze vorm van complexiteit om te gaan. Hier is een eenvoudige uitleg van wat de auteurs hebben gedaan en gevonden:

1. Het probleem: Veel noten tegelijk afstemmen

Meestal proberen wetenschappers het orkest af te stemmen om slechts één noot (de grondtoestand) perfect te spelen. Maar om te begrijpen hoe een molecuul op licht reageert, moeten ze ook over veel verschillende noten (geëxciteerde toestanden) tegelijkertijd beschikken.

• De uitdaging: Als je probeert 10 verschillende liedjes tegelijkertijd af te stemmen, worden de instructies (het computercircuit) ongelooflijk lang en ingewikkeld. Als de instructies te lang zijn, raakt de kwantumcomputer in de war door "ruis" (statische elektriciteit of fouten), en valt de muziek uit elkaar.
• De afweging: Je hebt een complex circuit nodig om een nauwkeurig antwoord te krijgen, maar een complex circuit is eerder vatbaar voor fouten op de huidige, ruisgevoelige machines.

2. De oplossing: Een slimme, symmetrische dirigent

De auteurs hebben een nieuwe methode ontwikkeld genaamd oo-MC-VQE. Zie dit als een "slimme dirigent" voor het kwantumorkest.

• Spin-geadapteerd: In de kwantumchemie hebben elektronen een eigenschap genaamd "spin" (zoals tollen die draaien). De auteurs hebben hun methode zo gebouwd dat de dirigent de tollen altijd op de juiste, symmetrische manier laat draaien. Dit voorkomt dat de muziek "vals" wordt door symmetriefouten.
• Orbitaal geoptimaliseerd: Ze lieten de dirigent ook de zitplaatsen van de muzikanten (de orbitalen) herordenen om de muziek beter te laten klinken nog voordat de complexe afstemming begint.
• Multistate gecontracteerd: In plaats van te proberen 10 liedjes af te stemmen met 10 afzonderlijke, enorme handleidingen, vonden ze een manier om één gedeelde, efficiënte set instructies te gebruiken die voor alle liedjes tegelijk werkt.

3. De ontdekking: Lineaire groei

Een van de grootste vragen was: Als ik 10 toestanden wil berekenen in plaats van 1, heb ik dan 10 keer zoveel rekenkracht nodig?

• De bevinding: De auteurs ontdekten dat het antwoord verrassend eenvoudig is. De hoeveelheid computer-"inspanning" (circuitparameters) die nodig is, groeit lineair. Als je het aantal toestanden dat je wilt berekenen verdubbelt, wordt de handleiding ongeveer ook twee keer zo lang. Dit is goed nieuws, want het betekent dat de methode schaalbaar is.

4. De realiteitstest: Spelen op een lawaaierig podium

De auteurs hebben dit niet alleen gesimuleerd op een perfecte computer; ze hebben hun methode daadwerkelijk getest op echte kwanthardware (IBM kwantumcomputers).

• De opstelling: Ze testten twee kleine moleculen: Formaldehyde (een veelvoorkomende chemische stof) en een triwaterstof-kation ($H_3^+$).
• Het ruisprobleem: Echte kwantumcomputers zijn als een podium met een luidruchtig publiek en flikkerende lichten. Zonder hulp waren de resultaten rommelig.
• De oplossing: Ze gebruikten "foutmitigatie"-technieken. Stel je dit voor als een geluidstechnicus die software gebruikt om het publieksgeluid en de flikkerende lichten na de uitvoering weg te filteren.
• Het resultaat:
  • Voor Formaldehyde werkte de methode vrij goed. Zelfs met de ruis konden ze de "pieken" in het absorptiespectrum (de kleuren die het molecuul absorbeert) duidelijk zien.
  • Voor $H_3^+$ was de ruis een groter probleem, waardoor de resultaten aanzienlijk verschoven. De auteurs merkten op dat dit kwam omdat de wiskunde van dit specifielich molecuul gevoeliger is voor ruis (zoals een delicaat instrument dat gemakkelijk vals wordt).
  • Belangrijkste les: Hoewel de cijfers op de echte machines niet perfect waren, was de vorm van de resultaten wel correct. Ze konden nog steeds de belangrijkste kenmerken van het gedrag van het molecuul waarnemen.

Samenvatting

Het artikel laat zien dat wetenschappers, door een slimme, symmetrische aanpak te gebruiken, het gedrag van geëxciteerde elektronen in moleculen kunnen berekenen met de huidige, imperfecte kwantumcomputers. Ze bewezen dat het berekenen van meerdere toestanden niet een onmogelijke hoeveelheid middelen vereist, en dat met wat "ruisonderdrukking"-trucs, er vandaag de dag al bruikbare chemische inzichten verkregen kunnen worden uit echte kwantumapparaten.

Wat ze NIET hebben beweerd:
Het artikel beweert niet dat deze methode onmiddellijk nieuwe zonnepanelen kan ontwerpen, ziekten kan genezen of nieuwe materialen kan creëren. Het richt zich strikt op het bewijzen dat de methode werkt voor het berekenen van spectra op kwanthardware. Eventuele toekomstige toepassingen worden gesuggereerd door het algemene vakgebied, maar zijn geen specifieke claims van dit onderzoek.

Technische samenvatting

Technische Samenvatting: Orbitaal-geoptimaliseerde spin-geadapteerde multistate gecontracteerde VQE

Probleemstelling
De robuuste en gebalanceerde berekening van elektronisch aangeslagen toestanden en hun overgangseigenschappen blijft een significante uitdaging, met name in regimes die betrokken zijn bij sterke elektronische reorganisatie, nabijheid van degeneratie of multiconfigurationeel karakter. Hoewel lineair-respons (LR) methoden succesvol zijn in veel contexten, behandelen zij aangeslagen toestanden als perturbaties van een grondtoestand-referentie, wat leidt tot conceptuele beperkingen nabij kegelverbindingen (conical intersections) en in sterk gecorreleerde systemen. Bovendien worden bestaande variatonele kwantum-eigenzone-algoritmen (VQE) voor aangeslagen toestanden, zoals subspace-search VQE (SS-VQE) en multistate gecontracteerde VQE (MC-VQE), vaak geconfronteerd met een afruil tussen de nauwkeurigheid van het beschrijven van een toenemend aantal elektronische toestanden en de expressiviteit (circuitcomplexiteit) van de onderliggende ansatz. Daarnaast is de huidige kwanthardware beperkt door ruis en decoherentie, wat de circuitdiepte beperkt en de ontwikkeling van methoden noodzakelijk maakt die zowel efficiënt als compatibel zijn met foutenmitigatietechnieken.

Methodologie
De auteurs introduceren de orbitaal-geoptimaliseerde multistate gecontracteerde variationele kwantum-eigenzone-algoritme (oo-MC-VQE) methode, die spin-geadapteerde operatoren incorporeert om grond- en aangeslagen toestanden, evenals staat-specifieke en overgangseigenschappen, te berekenen.

• Spin-adaptatie: Om te waarborgen dat de spin-symmetrie van referentietoestanden behouden blijft tijdens de VQE-optimalisatie, maakt de methode gebruik van een spin-geadapteerde operatorpool. Specifiek gebruiken de auteurs singlet spin-geadapteerde enkelvoudige excitatie-operatoren en paar-dubbele excitatie-operatoren. Dit voorkomt de noodzaak voor spin-straftermen en garandeert de correcte spin-symmetrie van alle inbegrepen toestanden zonder redundante parameters te introduceren.
• Multistate Gecontracteerd Kader: De methode optimaliseert meerdere elektronische toestanden simultaan met behulp van een staat-gemiddelde energiefunctional. Elke toestand wordt geconstrueerd met behulp van een vast, parameter-vrij staat-voorbereidingscircuit (het initialiseren van orthogonale basis-toestanden), gevolgd door een geparametriseerd correlatiecircuit (gebruikmakend van een Unitary Product State, of UPS, ansatz) dat elektronische correlatie introduceert. De optimalisatie richt zich op de staat-gemiddelde energie, gevolgd door een klassieke diagonalisatie van de multistate Hamiltoniaan-matrix om individuele toestand-energieën en golffuncties te resolveren.
• Orbitaal Optimalisatie: Het kader wordt uitgebreid met orbitaal optimalisatie (oo-MC-VQE). Orbitaalrotaties worden behandeld als een klassieke kostenfunctie, waardoor gelijktijdige optimalisatie van circuitparameters ($\theta$) en orbitaal-parameterparameters ($\kappa$) mogelijk is. De auteurs demonstreren dat voor de specifieke gebruikte ansatz, enkelvoudige excitatie-circuitparameters die grenzen aan de orbitaal-parameterisatie redundant zijn en kunnen worden verwijderd, wat de circuitcomplexiteit vermindert.
• Matrixelement Evaluatie: Om matrixelementen tussen verschillende toestanden te evalueren (vereist voor de multistate Hamiltoniaan en overgangseigenschappen), maken de auteurs gebruik van een superpositie-decompositietechniek. Deze aanpak vermijdt diepe staat-voorbereidingscircuits voor multi-determinant toestanden door eenvoudigere superposities van twee determinanten te evalueren, waardoor het aantal CNOT-gates wordt verminderd ten koste van een verhoogd aantal verwachtingswaarden-metingen.
• Foutenmitigatie: De implementatie integreert Mansatz0, een ansatz-gebaseerde readout en gate-foutenmitigatietechniek. Voor multistate berekeningen waarbij superposities van referentietoestanden betrokken zijn, stellen de auteurs een uitbreiding voor (M0+), die fouten karakteriseert voor alle unieke multi-determinant referentietoestanden, hoewel zij opmerken dat dit een exponentische schaling met zich meebrengt.

Belangrijkste Bijdragen

1. Spin-geadapteerde oo-MC-VQE: De ontwikkeling van een spin-geadapteerde formulering voor orbitaal-geoptimaliseerde staat-gemiddelde VQE, die de correcte spin-symmetrie conservering over meerdere toestanden waarborgt.
2. Schaalanalyse: Een expliciet onderzoek naar de afruil tussen nauwkeurigheid en circuitcomplexiteit. De auteurs vinden dat het aantal vereiste circuitparameters om nauwkeurige resultaten te verkrijgen ongeveer lineair schaalt met het aantal toestanden dat in de staat-gemiddelde wordt opgenomen.
3. Hardware Implementatie: Een expliciete kwantumcircuit-implementatie van de methode, inclusief de reductie van CNOT-gates via determinant-superpositie en de integratie van foutenmitigatietechnieken.
4. Real-Device Demonstratie: Uitvoering van de methode op echte kwanthardware (IBM ibm marrakesh en ibm aachen) om absorptiespectra van benchmark moleculaire systemen (formaldehyde en $H_3^+$) te berekenen.

Resultaten

• Circuitcomplexiteit: Simulaties op LiH, $H_2O$ en $NH_3$ met de ADAPT en oo-ADAPT ansatzes laten zien dat het aantal parameters vereist om nauwkeurige resultaten te verkrijgen, ruwweg lineair toeneemt met het aantal toestanden. Bijvoorbeeld: het bereiken van een fout onder $10^{-4}$ Hartree voor $H_2O$ vereiste respectievelijk 27, 32, 24 en 30 parameters per toestand voor 1, 2, 3 en 4 toestanden.
• Impact van Orbitaal Optimalisatie: Hoewel orbitaal optimalisatie de convergentie in de vroege stadia van de ADAPT-procedure verbetert, heeft het een verwaarloosbaar effect op het uiteindelijke aantal vereiste circuitparameters voor striktere convergentie. In sommige gevallen (bijv. $NH_3$ met 3 of 4 toestanden) had het inklappen van enkelvoudige excitaties in de orbitaal-optimalisatie een nadelig effect op de parameter-telling.
• Prestaties van Kwanthardware:
  • Gesimuleerde Ruis: Op geëmuleerde hardware met IBM-ruismodellen produceerden foutenmitigatietechnieken (M0 en M0+) resultaten die bijna identiek waren aan ruisloze simulaties. Het $H_3^+$ systeem werd sterker door ruis beïnvloed dan formaldehyde, wat wordt toegeschreven aan een hogere conditiegetal van de multistate Hamiltoniaan-matrix.
  • Echte Hardware: Berekeningen op echte apparaten toonden aan dat, hoewel fouten-gemitigeerde resultaten de kwalitatieve trends en primaire kenmerken van de absorptiespectra behielden, de kwantitatieve nauwkeurigheid aanzienlijk lager was dan op gesimuleerde apparaten. Dit verschil benadrukt de beperkingen van huidige ruismodellen met betrekking tot tijdsafhankelijke ruis en crosstalk. Desalniettemin slaagde de methode erin om ondanks de hardware-ruis betekenisvolle chemische inzichten te extraheren.

Betekenis en Claims
Het artikel claimt dat de oo-MC-VQE methode een gebalanceerde beschrijving van de elektronische structuur biedt die geschikt is voor uitdagende scenario's zoals kegelverbindingen, waarmee de beperkingen van lineair-respons benaderingen worden aangepakt. De auteurs benadrukken dat de lineaire schaling van circuitparameters met het aantal toestanden suggereert dat de methode levensvatbaar is voor het targeten van meerdere aangeslagen toestanden op near-term kwanthardware. Hoewel zij erkennen dat kwantitatieve nauwkeurigheid nog steeds een uitdaging is op de huidige apparaten, toont dit werk aan dat betekenisvolle chemische inzichten, zoals absorptiespectra, kunnen worden verkregen wanneer de methode wordt gecombineerd met geschikte foutenmitigatietechnieken. De studie dient als een proof-of-concept voor het uitvoeren van multistate, spin-geadapteerde, orbitaal-geoptimaliseerde berekeningen op echte kwantumapparaten.