Optimizing Unitary Coupled Cluster Wave Functions on Quantum Hardware: Error Bound and Resource-Efficient Optimizer

Technische Samenvatting: Optimalisatie van Unitary Coupled Cluster Golffuncties op Quantumhardware

Probleemstelling
Het artikel behandelt de uitdagingen bij het optimaliseren van Unitary Coupled Cluster (UCC) golffuncties op nabije quantumhardware om de grondtoestand van meer-deeltjessystemen te berekenen. Hoewel de Variational Quantum Eigensolver (VQE) een standaard hybride benadering is, lijdt deze onder drie primaire beperkingen:

1. Decoherentie: Beperkt de diepte en expressiviteit van circuits.
2. Meetoverhead: Het schatten van de variationele energie vereist een aantal 'shots' dat schaalt als $1/\epsilon^2$ voor een statistische fout $\epsilon$.
3. Optimalisatieproblematiek: De landschappen van VQE vertonen vaak "barren plateaus" (exponentieel vlakke kostenfuncties) en lokale minima, waardoor globale convergentie moeilijk is. Bovendien garandeert een kleine energiegradiënt in VQE niet de nabijheid van de ware grondtoestandsenergie.

Om deze problemen aan te pakken, onderzoeken de auteurs de Projective Quantum Eigensolver (PQE), een methode die in een eerder werk (Ref. 1) is geïntroduceerd en die de zoektocht naar de grondtoestand formuleert als een wortelzoekprobleem in plaats van een energie-minimalisatieprobleem. PQE optimaliseert parameters om de "residuen" (projecties van de Schrödinger-vergelijking op aangeslagen toestanden) naar nul te drijven, in plaats van direct de energie te minimaliseren.

Methodologie
Het artikel biedt een rigoureuze wiskundige analyse van het PQE-algoritme en stelt een nieuwe optimalisatiestrategie voor.

1. Wiskundig Kader:

   • Het algoritme zoekt naar het oplossen van $r_\nu(t) = \langle \Phi_\nu | \hat{U}^\dagger(t) \hat{H} \hat{U}(t) | \Phi_0 \rangle = 0$ voor alle aangeslagen determinanten $\nu \neq 0$.
   • De auteurs leiden een algemene foutgrens af voor het algoritme met behulp van Kato's ongelijkheid en Temple's ongelijkheid. Ze stellen een relatie vast tussen de gekwadrateerde som van residuen ($\sum r_\mu^2$) en de energie-fout ($|E(t) - E_{gs}|$) en de overlap met de grondtoestand.
   • In tegenstelling tot VQE, waarbij een kleine gradiëntnorm niet impliceert dat de energie-fout klein is, bewijzen de auteurs dat een kleine residu-norm formele garanties biedt op de nauwkeurigheid van de energie en de kwaliteit van de golffunctie-overlap.

2. Optimalisatiestrategie:

   • De standaard PQE-update (uit Ref. 1) benadert een gemodificeerde Newton-Raphson stap met behulp van een vaste noemer gebaseerd op de Møller-Plesset perturbatietheorie ($\Delta_\mu$).
   • De auteurs analyseren de convergentie-eigenschappen van quasi-Newton-methoden met behulp van de Newton-Kantorovich-stelling. Ze identificeren dat de standaard update-regel (gebaseerd op de Jacobiaan bij het initiële Hartree-Fock punt) kan falen wanneer de oplossing ver van de initiële gok ligt (bijv. in sterk gecorreleerde regimes).
   • Ze stellen een hybride optimizer voor die dynamisch schakelt tussen update-regels:
     • Ver van de oplossing: Gebruikt een gradiënt-afdaling-achtige update (of een vereenvoudigde Jacobiaan-benadering) om stabiliteit en vooruitgang vanaf de initiële toestand te waarborgen.
     • Nabij de oplossing: Schakelt over naar een benaderde Newton-Raphson update (met behulp van een verfijnde Jacobiaan-benadering) om snellere, kwadratisch-achtige convergentie te bereiken.

Belangrijkste Bijdragen

1. Rigoureuze Foutgrenzen: Het artikel leidt nauwe wiskundige grenzen af die de off-diagonaal coëfficiënten (residuen) van de Hamiltoniaan relateren aan de energie-fout en de golffunctie-overlap. Dit biedt formele convergentiegaranties die ontbreken in VQE.
2. Convergentiecriterium: Op basis van de afgeleide grenzen formuleren de auteurs een goed geïnformeerd stopcriterium ($\epsilon_T$) waarmee gebruikers de energie-fout kunnen schatten zonder de exacte grondtoestand te kennen.
3. Nieuwe Optimizer: Een nieuwe residu-gebaseerde optimizer wordt geïntroduceerd die zijn update-regel aanpast op basis van de afstand tot de oplossing, waarbij gebruik wordt gemaakt van de Newton-Kantorovich-stelling om robuustheid te garanderen.
4. Numerieke Benchmarks: Uitgebreide simulaties op $H_4$, $H_6$, $BeH_2$ en $LiH$ dissociatiecurven.

Resultaten

• Strakheid van de Foutgrens: Numerieke tests op een $H_4$-molecuul laten zien dat de afgeleide foutgrens de ware energie-fout volgt binnen minder dan één orde van grootte, wat aanzienlijk nauwer is dan eerdere grenzen gebaseerd op de residu 1-norm.
• Prestaties van de Optimizer:
  • De nieuwe hybride optimizer vertoont een superieure robuustheid vergeleken met de standaard PQE-update (Ref. 1), met name bij grote bindingsafstanden waar de statische correlatie hoog is.
  • Wat betreft de efficiëntie van de middelen, vereist de nieuwe methode aanzienlijk minder metingen (en dus minder runtime) om een doel-nauwkeurigheid te bereiken in vergelijking met zowel de standaard PQE als de met de Broyden–Fletcher–Goldfarb–Shanno (BFGS) methode geoptimaliseerde VQE.
• Convergentie: De methode convergeert succesvol naar de grondtoestandsenergie met minder optimalisatiestappen dan VQE, waarbij de problemen van barren plateaus geassocieerd met energie-minimalisatie worden vermeden.

Significantie
Het artikel beweert dat door de verschuiving van variationele energie-minimalisatie naar projectieve wortelzoeken, PQE een betrouwbaarder pad biedt naar de voorbereiding van de grondtoestand op nabije hardware. De primaire significantie ligt in:

• Het bieden van formele garanties dat kleine residuen kleine energie-fouten impliceren, wat de onzekerheid die inherent is aan VQE-optimalisatie aanpakt.
• Het bieden van een efficiënter alternatief qua middelen dat de meetoverhead en de optimalisatiestappen vermindert die nodig zijn voor chemische nauwkeurigheid.
• Het vestigen van een wiskundige fundering (via Newton-Kantorovich-analyse) voor het ontwerpen van robuuste optimizers die de niet-lineariteiten van sterk gecorreleerde systemen kunnen afhandelen waar standaard lineaire benaderingen falen.

De auteurs positioneren dit werk als een stap naar praktisch quantumvoordeel op de korte termijn, door een methode te bieden die minder gevoelig is voor de specifieke optimalisatiepathologieën van VQE, terwijl het toch haalbaar blijft binnen de huidige hardwarebeperkingen.