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

Résumé technique : Optimisation des fonctions d'onde de type Unitary Coupled Cluster sur le matériel quantique

Énoncé du problème
L'article traite des défis liés à l'optimisation des fonctions d'onde de type Unitary Coupled Cluster (UCC) sur le matériel quantique de nouvelle génération pour calculer l'état fondamental de systèmes à corps multiples. Bien que l'algorithme Variational Quantum Eigensolver (VQE) soit une approche hybride standard, il souffre de trois limitations principales :

1. Décohérence : Limite la profondeur des circuits et l'expressivité.
2. Surcharge de mesure : L'estimation de l'énergie variationnelle nécessite un nombre de tirages (shots) proportionnel à $1/\epsilon^2$ pour une erreur statistique $\epsilon$.
3. Difficulté d'optimisation : Les paysages de coût du VQE présentent souvent des « plateaux stériles » (barren plateaus — fonctions de coût exponentiellement plates) et des minima locaux, rendant la convergence globale difficile. De plus, un petit gradient d'énergie dans le VQE ne garantit pas la proximité de l'énergie réelle de l'état fondamental.

Pour répondre à ces problèmes, les auteurs étudient le Projective Quantum Eigensolver (PQE), une méthode introduite dans un travail précédent (Réf. 1) qui formule la recherche de l'état fondamental comme un problème de recherche de racine plutôt que comme un problème de minimisation d'énergie. Le PQE optimise les paramètres pour pousser les « résidus » (projections de l'équation de Schrödinger sur les états excités) vers zéro, plutôt que de minimiser directement l'énergie.

Méthodologie
L'article fournit une analyse mathématique rigoureuse de l'algorithme PQE et propose une nouvelle stratégie d'optimisation.

1. Cadre mathématique :

   • L'algorithme cherche à résoudre $r_\nu(t) = \langle \Phi_\nu | \hat{U}^\dagger(t) \hat{H} \hat{U}(t) | \Phi_0 \rangle = 0$ pour tous les déterminants excités $\nu \neq 0$.
   • Les auteurs dérivent une borne d'erreur générale pour l'algorithme en utilisant l'inégalité de Kato et l'inégalité de Temple. Ils établissent une relation entre la somme des carrés des résidus ($\sum r_\mu^2$) et l'erreur d'énergie ($|E(t) - E_{gs}|$) ainsi que le recouvrement avec l'état fondamental.
   • Contrairement au VQE, où la norme d'un petit gradient n'implique pas une petite erreur d'énergie, les auteurs prouvent qu'une petite norme de résidu fournit des garanties formelles sur la précision de l'énergie et la qualité du recouvrement de la fonction d'onde.

2. Stratégie d'optimisation :

   • La mise à jour standard du PQE (de la Réf. 1) approxime une étape de Newton-Raphson modifiée en utilisant un dénominateur fixe basé sur la théorie de la perturbation de Møller-Plesset ($\Delta_\mu$).
   • Les auteurs analysent les propriétés de convergence des méthodes quasi-Newton en utilisant le théorème de Newton-Kantorovich. Ils identifient que la règle de mise à jour standard (basée sur le Jacobien au point Hartree-Fock initial) peut échouer lorsque la solution est éloignée de l'estimation initiale (par exemple, dans les régimes de corrélation forte).
   • Ils proposent un optimiseur hybride qui bascule dynamiquement entre les règles de mise à jour :
     • Loin de la solution : Utilise une mise à jour de type descente de gradient (ou une approximation simplifiée du Jacobien) pour assurer la stabilité et la progression depuis l'état initial.
     • Près de la solution : Bascule vers une mise à jour Newton-Raphson approximative (utilisant une approximation raffinée du Jacobien) pour atteindre une convergence de type quadratique plus rapide.

Contributions clés

1. Bornes d'erreur rigoureuses : L'article dérive des bornes mathématiques serrées reliant les coefficients hors-diagonaux (résidus) de l'Hamiltonien à l'erreur d'énergie et au recouvrement de la fonction d'onde. Cela fournit des garanties de convergence formelles absentes dans le VQE.
2. Critère de convergence : Sur la base des bornes dérivées, les auteurs formulent un critère d'arrêt bien informé ($\epsilon_T$) qui permet aux utilisateurs d'estimer l'erreur d'énergie sans connaître l'état fondamental exact.
3. Nouvel optimiseur : Un nouvel optimiseur basé sur les résidus est introduit, qui adapte sa règle de mise à jour en fonction de la distance à la solution, exploitant le théorème de Newton-Kantorovich pour garantir la robustesse.
4. Benchmarks numériques : Simulations étendues sur les courbes de dissociation de $H_4$, $H_6$, $BeH_2$ et $LiH$.

Résultats

• Étroitesse de la borne d'erreur : Des tests numériques sur une molécule $H_4$ montrent que la borne d'erreur dérivée suit l'erreur d'énergie réelle à moins d'un ordre de grandeur, ce qui est nettement plus serré que les bornes précédentes basées sur la norme 1 du résidu.
• Performance de l'optimiseur :
  • Le nouvel optimiseur hybride démontre une robustesse supérieure par rapport à la mise à jour standard du PQE (Réf. 1), particulièrement pour les grandes distances de liaison où la corrélation statique est élevée.
  • En termes d'efficacité des ressources, la nouvelle méthode nécessite nettement moins de mesures (et donc moins de temps d'exécution) pour atteindre une précision cible par rapport au PQE standard et au VQE optimisé avec la méthode BFGS (Broyden–Fletcher–Goldfarb–Shanno).
• Convergence : La méthode converge avec succès vers l'énergie de l'état fondamental avec moins d'étapes d'optimisation que le VQE, évitant les problèmes de plateaux stériles associés à la minimisation d'énergie.

Signification
L'article soutient qu'en passant de la minimisation variationnelle de l'énergie à la recherche de racines projective, le PQE offre une voie plus fiable pour la préparation de l'état fondamental sur le matériel de nouvelle génération. La signification primaire réside dans :

• La fourniture de garanties formelles que de petits résidus impliquent de petites erreurs d'énergie, répondant ainsi à l'incertitude inhérente à l'optimisation VQE.
• L'offre d'une alternative efficace en ressources qui réduit la surcharge de mesure et les étapes d'optimisation nécessaires pour atteindre la précision chimique.
• L'établissement d'un fondement mathématique (via l'analyse de Newton-Kantorovich) pour concevoir des optimiseurs robustes capables de gérer les non-linéarités des systèmes fortement corrélés là où les approximations linéaires standards échouent.

Les auteurs positionnent ce travail comme une étape vers un avantage quantique pratique à court terme, proposant une méthode moins sensible aux pathologies d'optimisation spécifiques du VQE tout en restant réalisable sous les contraintes matérielles actuelles.