From Projected Subspaces to Full-Space Implementations: Representation-Equivalence Audits for Open-Shell ADAPT-VQE

In de zoektocht naar het simuleren van het gedrag van atomen en moleculen op een quantumcomputer, worden wetenschappers geconfronteerd met een fundamentele uitdaging: hoe beschrijf je de complexe dans van elektronen zonder de machine te overbelasten. Een populaire strategie houdt in dat men een wiskundig model van het molecuul bouwt in een vereenvoudigde, gereduceerde ruimte waar de belangrijkste interacties behouden blijven, terwijl minder cruciale details tijdelijk opzij worden gezet. Deze aanpak, bekend als een geprojecteerde subruimdemethode, stelt onderzoekers in staat om veel sneller een laagtoestandsniveau voor het systeem te vinden. De hoop is geweest dat zodra deze efficiënte oplossing in de vereenvoudigde wereld is gevonden, deze direct kan worden vertaald naar de volledige, complexe realiteit van het molecuul om op werkelijke quantumhardware te draaien. De aanname is dat de afkorting die in de wiskunde wordt genomen, de fysieke uitkomst niet verandert.

Echter, een nieuwe studie onthult dat deze vertalingsstap veel verraderlijker is dan voorheen werd gedacht. De onderzoekers ontdekten dat een oplossing die gevonden is in een vereenvoudigde, symmetrie-geprojecteerde ruimte, niet automatisch dezelfde is als de oplossing die men zou krijgen door exact dezelfde instructies uit te voeren in de volledige, ongeprojecteerde ruimte. In feite, voor bepaalde typen chemische systemen, leidt het direct toepassen van de parameters uit de vereenvilde berekening op het volledige systeem tot een catastrofale mislukking, waarbij een toestand wordt geproduceerd die fysiek foutief en chemisch nutteloos is. De studie verwerpt het idee van het gebruik van vereenvoudigde ruimtes niet; in plaats daarvan stelt zij een strikte nieuwe regel vast: voordat een resultaat uit een gereduceerde-ruimteberekening als een volledige-ruimteoplossing mag worden vertrouwd, moet het een specifieke audit ondergaan om te bewijzen dat de twee ruimtes werkelijk equivalent zijn.

Het team, onder leiding van Lucia Malíčková en Petr Klenovský, testte dit idee met behulp van een model van een koperion, een veelvoorkomend onderdeel in biologische systemen en industriële katalysatoren. Ze voerden eerst een standaard adaptief algoritme uit in een gereduceerde ruimte die de juiste spin-eigenschappen van de elektronen conserveert. Deze berekening verliep perfect en vond een toestand met een energieafwijking van slechts 1,57 milli-Hartree, een eenheid van energie gebruikt in de quantumchemie. Dit resultaat werd beschouwd als een succes binnen het vereenvoudigde kader. De onderzoekers namen vervolgens exact dezelfde sequentie van wiskundige operaties en exact dezelfde numerieke waarden voor de stappen, maar pasten deze toe op de volledige, ongeprojecteerde ruimte van het molecuul. Het resultaat was een ramp. De energieafwijking schoot omhoog naar 464,73 milli-Hartree, en de resulterende toestand vertoonde bijna geen enkele gelijkenis met de correcte toestand, met een fidelity van slechts 0,07. In eenvoudige bewoordingen: de instructies die werkten in de vereenvoudigde kamer, leidden het systeem volledig de verkeerde kant op toen ze werden toegepast op het volledige gebouw.

Het onderzoek toonde aan dat dit falen niet te wijten was aan één enkele slechte stap, maar een systemisch probleem was. Van de 323 mogelijke typen elektronbewegingen die in de berekening werden gebruikt, overtraden er 274 een specifieke wiskundige voorwaarde die vereist is om het systeem binnen de juiste spin-toestand te houden wanneer men tussen ruimtes beweegt. Toen de onderzoekers probeerden het probleem op te lossen door de getallen in de volledige ruimte te heroptimaliseren, terwijl ze het systeem dwongen dicht bij de juiste spin-toestand te blijven, verbeterde de fout aanzienlijk en daalde deze naar 8,48 milli-Hartree. Dit bewees dat het oorspronkelijke falen primair een probleem was van het overdragen van parameters van de ene wiskundige representatie naar de andere, in plaats van dat de sequentie van stappen inherent niet in staat was het probleem op te lossen. Echter, zelfs met deze heroptimalisatie kon de volledige-ruimte sequentie de hoge precisie van het oorspronkelijke gereduceerde-ruimteresultaat niet quite bereiken, wat aangeeft dat de twee benaderingen fundamenteel verschillende manieren zijn om de oplossing te construeren.

Om te bevestigen dat dit een algemeen probleem was en niet slechts een eigenaardigheid van het koper-model, testte het team een stikstofmonoxide-radicaal. Hier was het effect minder ernstig, maar nog steeds aanwezig. De gereduceerde-ruimteberekening vond een goede oplossing, maar het toepassen van dezelfde getallen in de volledige ruimte resulteerde in een veel grotere fout, hoewel het systeem dichter bij de correcte toestand bleef dan in het geval van het koper. Dit bevestigde dat het mechanisme van falen bestaat bij verschillende moleculen, maar dat de ernst ervan afhangt van het specifieke systeem en de volgorde van de stappen. De onderzoekers demonstreerden ook dat als men de volledige-ruimteberekening vanaf het begin opbouwt met een pool van operatoren die gegarandeerd de spin-symmetrie behouden, de gereduceerde-ruimte en volledige-ruimte resultaten perfect overeenkomen. Deze "positieve controle" toonde aan dat het probleem oplosbaar is, maar dat het een andere constructiemethode vereist dan simpelweg het vertalen van een gereduceerde-ruimte resultaat.

De studie behandelde ook de praktische kant van het draaien van deze simulaties op echte quantumhardware. Zelfs als de wiskundige toestand correct is, houdt het omzetten hiervan in een fysieke circuit benaderingen in, zoals het opdelen van complexe operaties in eenvoudigere gates. Het team toonde aan dat door zorgvuldig te selecteren welke delen van het circuit meer precisie vereisten en welke vereenvoudigd konden worden, ze het aantal quantum-gates met bijna 70 procent konden verminderen, terwijl ze toch de nauwkeurigheid behielden die nodig is om een strikte validatietest te doorstaan. Dit benadrukt dat het proces van het bouwen van het circuit en het proces van het meten van de uiteindelijke energie twee aparte lagen van validatie zijn, elk met een eigen potentieel voor fouten.

Uiteindelijk dient dit werk als een cruciaal validatiekader in plaats van een nieuw algoritme voor het oplossen van chemische problemen. Het waarschuwt onderzoekers dat een lage energiescore in een gereduceerde-ruimteberekening niet garandeert dat de corresponderende volledige-ruimte circuit dezelfde toestand zal produceren. De auteurs concluderen dat telkens wanneer een gereduceerde-ruimte optimalisatie verbonden is aan een volledige-ruimte implementatie, de geoptimaliseerde toestand, de ouder-sequentie van operaties, het gecompileerde circuit en het meetprotocol als afzonderlijke objecten gevalideerd moeten worden. De lage energie van het vereenvoudigde model alleen is niet voldoende om te certificeren dat de volledige-ruimte circuit dezelfde fysieke realiteit implementeert. Deze bevinding zorgt ervoor dat toekomstige simulaties van complexe moleculen, van koperproteïnen tot nieuwe materialen, gebouwd worden op een fundament van geverifieerde equivalentie in plaats van onbewezen aannames.