On the generalized eigenvalue problem in subspace-based excited state methods for quantum computers

Deze studie toont aan dat subspace-gebaseerde methoden voor het berekenen van aangeslagen toestanden op kwantumcomputers, zoals QSE en qEOM, instabiel worden door statistische meetfouten bij hoge conditienummers van de overlapmatrix, terwijl methoden zoals q-sc-EOM die een standaard eigenwaardevergelijking gebruiken, robuuster en geschikter zijn.

De kern

Hier is een uitleg van dit wetenschappelijke artikel in eenvoudig Nederlands, met behulp van alledaagse vergelijkingen.

De Grote Droom: Moleculen op een Quantumcomputer

Stel je voor dat je een quantumcomputer gebruikt om te kijken hoe moleculen (de bouwstenen van alles om ons heen) zich gedragen. Dat is heel lastig, vooral als je niet alleen naar de rustige toestand kijkt, maar naar de geëxciteerde toestanden. Dat zijn momenten waarop een molecuul energie heeft opgepikt, bijvoorbeeld door licht, en dan gaat 'dansen'.

Wetenschappers willen deze 'dans' simuleren om nieuwe medicijnen of materialen te ontwerpen. Maar quantumcomputers zijn nog niet perfect; ze maken ruis (foutjes) door de manier waarop ze meten.

Het Probleem: De "Wobbly" Ladder

In dit artikel kijken de auteurs naar twee populaire manieren om deze geëxciteerde toestanden te berekenen: QSE en q-sc-EOM.

Stel je voor dat je een ladder moet bouwen om aan het plafond te komen (de juiste energie van het molecuul).

1. De QSE-methode (De instabiele ladder):
   Bij deze methode moet je eerst een heel complex rekensommetje maken waarbij je een "overlap-matrix" moet omkeren.

   • De Analogie: Stel je voor dat je een ladder bouwt op een helling. Als de grond heel stevig is (een lage "conditienummer"), staat de ladder stabiel. Maar als de grond zacht en modderig wordt (een hoog conditienummer), wordt de ladder extreem wankel.
   • Het effect: Zelfs een klein beetje trillen (de ruis van de quantumcomputer) zorgt ervoor dat de ladder helemaal uit elkaar valt. De berekende energieën worden dan volledig verkeerd. Als de grond heel zacht wordt, kan je de ladder zelfs niet meer opzetten zonder dat hij instort. Dan moet je een deel van de ladder weghalen (een techniek genaamd "thresholding"), maar dan mis je de bovenste sporten. Je komt er wel, maar je mist belangrijke informatie (de bovenste geëxciteerde toestanden).

2. De q-sc-EOM-methode (De stabiele ladder):
   Deze methode is slim ontworpen.

   • De Analogie: Deze ladder staat op een betonnen vloer. De grond is perfect vlak en stevig. Het maakt niet uit hoe veel je trilt of hoe veel ruis er is; de ladder blijft staan.
   • Het effect: Omdat de wiskundige structuur hier anders is (de overlap-matrix is altijd perfect), worden de foutjes niet versterkt. Je krijgt een stabiel resultaat, zelfs als de quantumcomputer wat ruis maakt.

Wat hebben de onderzoekers ontdekt?

De auteurs hebben gekeken naar wat er gebeurt als je deze methoden test met verschillende soorten "grond" (verschillende moleculen en situaties):

• Wanneer de grond stevig is (Laag conditienummer): Beide methoden werken ongeveer even goed. Geen probleem.
• Wanneer de grond modderig wordt (Gemiddeld conditienummer): De QSE-methode begint te haperen. De fouten worden groter en groter. Je moet heel veel metingen doen (duizenden keren meten) om nog een beetje nauwkeurig te zijn. De q-sc-EOM-methode blijft rustig en nauwkeurig.
• Wanneer de grond een modderpoel is (Hoog conditienummer): De QSE-methode faalt volledig. De berekening is onmogelijk zonder ingrepen. Als je de "modder" weghaalt om het op te lossen, verlies je de bovenste delen van je ladder. Je ziet dan niet alle mogelijke energietoestanden van het molecuul. De q-sc-EOM-methode werkt nog steeds perfect.

Waarom is dit belangrijk?

In de chemie is het cruciaal om alle mogelijke energietoestanden te zien. Als je een methode gebruikt die sommige toestanden laat verdwijnen (zoals QSE in moeilijke situaties), kun je denken dat een molecuul veilig is, terwijl het in werkelijkheid een gevaarlijke reactie kan hebben.

De conclusie in één zin:
Voor het simuleren van complexe moleculen op quantumcomputers is de q-sc-EOM-methode veel betrouwbaarder dan de QSE-methode, omdat hij niet zo snel instort door de onvermijdelijke ruis van de computer. Het is als kiezen voor een betonnen vloer in plaats van een modderige helling als je wilt bouwen.

Samenvatting voor de leek

• QSE: Werkt goed op een rustige dag, maar als het een beetje wazig wordt (ruis), wordt het resultaat onbetrouwbaar en mis je soms belangrijke stukken van het plaatje.
• q-sc-EOM: Werkt altijd stabiel, ongeacht hoe wazig de dag is. Je krijgt een compleet en betrouwbaar plaatje.

Voor de toekomst van quantumchemie is dit een belangrijk advies: kies voor de methode die niet instort als de computer een beetje trilt.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "On the generalized eigenvalue problem in subspace-based excited state methods for quantum computers", geschreven in het Nederlands.

Titel: Over het gegeneraliseerde eigenwaardeprobleem in subspace-gebaseerde methoden voor aangeslagen toestanden op quantumcomputers

Auteurs: Prince Frederick Kwao, Srivathsan Poyyapakkam Sundar, Brajesh Gupt en Ayush Asthana (Universiteit van North Dakota).

---

1. Het Probleem

Het simuleren van aangeslagen toestanden in de kwantumchemie is een van de meest veelbelovende toepassingen voor quantumcomputers. Voor pre-fouttolerante quantumapparaten zijn subspace-gebaseerde algoritmen, zoals Quantum Subspace Expansion (QSE), Quantum Equation of Motion (qEOM) en Quantum Self-consistent Equation-of-Motion (q-sc-EOM), veelbelovend.

De kernuitdaging die in dit artikel wordt geanalyseerd, is de numerieke stabiliteit van deze methoden in aanwezigheid van statistische steekproeffouten (shot noise), die inherent zijn aan quantummetingen.

• QSE en qEOM vereisen het oplossen van een gegeneraliseerd eigenwaardeprobleem van de vorm $HC = SCE$, waarbij $S$ de overlap-matrix is.
• q-sc-EOM vereist het oplossen van een standaard eigenwaardeprobleem van de vorm $HC = CE$, omdat de overlap-matrix analytisch de eenheidsmatrix is ($S=I$).

Het probleem is dat het oplossen van een gegeneraliseerd eigenwaardeprobleem de inversie van de overlap-matrix $S$ vereist. Als $S$ een hoge conditienummer heeft (d.w.z. slecht geconditioneerd is), worden de fouten in de eigenwaarden enorm versterkt door de inversiestap, wat leidt tot instabiliteit of zelfs singulariteiten (niet-oplosbaarheid) op quantumcomputers.

2. Methodologie

De auteurs voerden zowel een analytische als een numerieke analyse uit om de gevoeligheid van deze methoden voor steekproeffouten te kwantificeren.

• Analytische Storingstheorie:

  • Ze toonden aan dat voor een standaard eigenwaardeprobleem de fout in de eigenwaarde ($\delta\lambda$) lineair schaalt met de meetfout ($\epsilon$).
  • Voor een gegeneraliseerd eigenwaardeprobleem schaalt de fout echter met een factor die omgekeerd evenredig is met de kleinste eigenwaarde van de overlap-matrix ($\lambda_{min}(S)$). De foutvergelijking wordt:
    $$|\delta\lambda| \approx O\left(\frac{n(1+|\lambda|)}{\lambda_{min}(S)} \epsilon\right)$$
  • Dit betekent dat bij een kleine $\lambda_{min}(S)$ (hoge conditie) de fouten drastisch worden versterkt.
  • Ze deriveerden een praktische richtlijn voor het benodigde aantal schoten ($N_{shots}$) om een bepaalde nauwkeurigheid te bereiken, wat kwadratisch toeneemt naarmate $\lambda_{min}(S)$ afneemt.

• Numerieke Simulaties:

  • De methoden werden getest op het H4-molecuul (met verschillende geometrieën) en NH3 (ammoniak) in een STO-3G basis.
  • Ze varieerden de conditienummers van de overlap-matrix in QSE van laag (<100), medium (100–10.000) tot zeer hoog (>10.000).
  • Ze vergeleken de resultaten van QSE met die van q-sc-EOM onder dezelfde omstandigheden van steekproeffouten.
  • Voor het geval van zeer hoge conditienummers werd de thresholding-techniek (drempelwaarde) toegepast op QSE, waarbij kleine eigenwaarden van $S$ worden verwijderd om singulariteiten te voorkomen.

3. Belangrijkste Bijdragen

1. Kwantificering van Foutversterking: Het artikel levert een strikte theoretische onderbouwing en numerieke bewijzen dat de instabiliteit van QSE en qEOM direct voortkomt uit de inversie van de overlap-matrix bij hoge conditienummers.
2. Praktische Diagnostiek: De auteurs koppelen de conditie van de overlap-matrix aan het benodigde aantal schoten en de betrouwbaarheid van de resultaten, wat een bruikbaar hulpmiddel biedt voor toekomstige experimenten.
3. Vergelijking van Methodes: Ze demonstreren dat q-sc-EOM, door het gebruik van een orthonormale basis (waarbij $S=I$), volledig immuun is voor dit specifieke type foutversterking.
4. Gevolgen van Thresholding: Ze tonen aan dat hoewel thresholding QSE weer oplosbaar maakt bij hoge conditienummers, dit leidt tot onvolledigheid in het spectrum (verlies van aangeslagen toestanden), wat chemisch problematisch kan zijn.

4. Resultaten

De resultaten worden onderverdeeld in drie scenario's gebaseerd op het conditienummer ($\kappa(S)$) van de overlap-matrix:

• Laag Conditienummer ($\kappa < 100$):

  • Zowel QSE als q-sc-EOM vertonen vergelijkbare stabiliteit en convergentie. De foutversterking is verwaarloosbaar.
  • Beide methoden bereiken vergelijkbare nauwkeurigheid met een vergelijkbaar aantal schoten.

• Medium Conditienummer ($100 < \kappa < 10^4$):

  • QSE: Wordt aanzienlijk instabieler. De fouten in de gemeten energieën en de variantie nemen sterk toe naarmate het conditienummer stijgt. Om dezelfde nauwkeurigheid als q-sc-EOM te bereiken, zijn er orde van grootte meer schoten nodig (bijvoorbeeld $10^5$schoten voor QSE versus$10^4$voor q-sc-EOM bij$\kappa \approx 592$).
  • q-sc-EOM: Blijft stabiel en nauwkeurig, ongeacht het conditienummer van de QSE-overlap, omdat er geen inversie van $S$ nodig is.

• Hoog Conditienummer ($\kappa > 10^4$):

  • QSE: Het gegeneraliseerde eigenwaardeprobleem wordt singulier en kan niet worden opgelost zonder extra stappen.
  • Thresholding: Het toepassen van een drempelwaarde maakt de oplossing mogelijk, maar resulteert in verlies van aangeslagen toestanden (vooral ontaarde toestanden) uit het berekende spectrum. Dit komt doordat de effectieve subspace wordt verkleind.
  • q-sc-EOM: Biedt een volledig en stabiel spectrum zonder dat thresholding nodig is.

De analyse van het NH3-systeem bevestigde deze trends, waarbij zelfs bij relatief lagere conditienummers (124 en 427) al significante instabiliteit in QSE werd waargenomen.

5. Betekenis en Conclusie

De studie concludeert dat subspace-methoden die afhankelijk zijn van het oplossen van een gegeneraliseerd eigenwaardeprobleem (zoals QSE en qEOM) fundamenteel kwetsbaar zijn voor de onvermijdelijke statistische ruis van quantumcomputers, vooral bij complexe moleculaire geometrieën die leiden tot slecht geconditioneerde overlap-matrices.

q-sc-EOM wordt gepresenteerd als een superieur alternatief voor aangeslagen-toestandenberekeningen op quantumcomputers omdat:

1. Het werkt met een standaard eigenwaardeprobleem ($S=I$), waardoor de foutversterkende inversiestap wordt vermeden.
2. Het robuuster is tegen shot noise en minder rekenresources (aantal schoten) vereist voor vergelijkbare nauwkeurigheid.
3. Het volledige spectra behoudt, zelfs in situaties waar QSE faalt of toestanden moet "verliezen" via thresholding.

Voor chemische toepassingen, zoals het bestuderen van conische doorsneden of dubbel geëxciteerde toestanden, waar de volledigheid van het spectrum cruciaal is, is q-sc-EOM daarom de aanbevolen methode. De bevindingen zijn ook relevant voor andere quantumalgoritmen die gegeneraliseerde eigenwaardeproblemen oplossen in de aanwezigheid van ruis.