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

Cet article présente la méthode de VQE contractée multi-états à optimisation orbitale (oo-MC-VQE), qui utilise des opérateurs adaptés au spin pour calculer efficacement les états fondamentaux et excités ainsi que leurs propriétés sur du matériel quantique, tout en équilibrant précision et complexité de circuit grâce à une mise à l'échelle linéaire des paramètres avec le nombre d'états.

L'essentiel

Imaginez que vous essayez d'accorder un orchestre massif et complexe pour jouer une symphonie parfaite. Dans le monde de la chimie, l'« orchestre » est une molécule, et la « musique » est la façon dont ses électrons se déplacent et interagissent. Pour comprendre comment une molécule absorbe la lumière (ce qui donne les couleurs et régit des processus comme la photosynthèse), les scientifiques doivent calculer les notes exactes que jouent ces électrons.

Pendant longtemps, calculer cela pour des molécules possédant de nombreux électrons a été comparable à la résolution d'un puzzle qui devient exponentiellement plus difficile à mesure que l'on ajoute des pièces. Les ordinateurs classiques (ceux que nous utilisons aujourd'hui) finissent par heurter un mur et ne peuvent pas résoudre ces énigmes pour des molécules complexes.

Ce document présente une nouvelle façon de résoudre ces puzzles en utilisant des ordinateurs quantiques, des machines spéciales conçues pour gérer ce type de complexité naturellement. Voici une décomposition simple de ce que les auteurs ont fait et découvert :

1. Le problème : Accorder de nombreuses notes à la fois

Habituellement, les scientifiques essaient d'accorder l'orchestre pour qu'il joue parfaitement une seule note (l'état fondamental). Mais pour comprendre comment une molécule réagit à la lumière, ils doivent connaître de nombreuses notes différentes (états excités) en même temps.

• Le défi : Si vous essayez d'accorder l'orchestre pour 10 chansons différentes simultanément, les instructions (le circuit informatique) deviennent incroyablement longues et compliquées. Si les instructions sont trop longues, l'ordinateur quantique est perturbé par le « bruit » (statique ou erreurs), et la musique s'effondre.
• Le compromis : Vous avez besoin d'un circuit complexe pour obtenir une réponse précise, mais un circuit complexe est plus susceptible d'échouer sur les machines bruyantes actuelles.

2. La solution : Un chef d'orchestre intelligent et symétrique

Les auteurs ont développé une nouvelle méthode appelée oo-MC-VQE. Voyez cela comme un « chef d'orchestre intelligent » pour l'orchestre quantique.

• Adapté au spin (Spin-Adapted) : En chimie quantique, les électrons possèdent une propriété appelée « spin » (comme des toupies). Les auteurs ont conçu leur méthode de sorte que le chef d'orchestre maintienne toujours les toupies en rotation de la manière symétrique correcte. Cela empêche la musique de se « désaccorder » à cause d'erreurs de symétrie.
• Optimisation des orbitales (Orbital Optimized) : Ils ont également permis au chef d'orchestre de réorganiser le plan de table des musiciens (les orbitales) pour que la musique sonne mieux avant même de commencer l'accordage complexe.
• Contraction multi-états (Multistate Contracted) : Au lieu d'essayer d'accorder 10 chansons avec 10 manuels d'instructions massifs et distincts, ils ont trouvé un moyen d'utiliser un ensemble d'instructions partagé et efficace qui fonctionne pour toutes les chansons à la fois.

3. La découverte : Une croissance linéaire

L'une des plus grandes questions était : Si je veux calculer 10 états au lieu de 1, ai-je besoin de 10 fois plus de puissance de calcul ?

• La conclusion : Les auteurs ont découvert que la réponse est étonnamment simple. La quantité d'« effort » informatique (paramètres du circuit) nécessaire croît de manière linéaire. Si vous doublez le nombre d'états que vous voulez calculer, vous doublez approximativement la longueur du manuel d'instructions. Cela ne devient pas une taille impossible. C'est une excellente nouvelle car cela signifie que la méthode est évolutive.

4. Le test en conditions réelles : Jouer sur une scène bruyante

Les auteurs n'ont pas seulement simulé cela sur un ordinateur parfait ; ils ont réellement testé leur méthode sur du matériel quantique réel (ordinateurs quantiques IBM).

• La configuration : Ils ont testé deux petites molécules : le formaldéhyde (un produit chimique courant) et un cation trihydrogène ($H_3^+$).
• Le problème du bruit : Les ordinateurs quantiques réels sont comme une scène avec une foule bruyante et des lumières vacillantes. Les résultats étaient désordonnés sans aide.
• La solution : Ils ont utilisé des techniques d'« atténuation d'erreurs » (error mitigation). Imaginez cela comme un ingénieur du son utilisant un logiciel pour filtrer le bruit de la foule et les lumières vacillantes après la performance.
• Le résultat :
  • Pour le formaldéhyde, la méthode a bien fonctionné. Même avec le bruit, on pouvait clairement voir les « pics » dans le spectre d'absorption (les couleurs que la molécule absorbe).
  • Pour le $H_3^+$, le bruit était un problème plus important, décalant les résultats de manière significative. Les auteurs ont noté que la mathématique de cette molécule spécifique est plus sensible au bruit (comme un instrument délicat qui se désaccorde facilement).
  • Point clé : Bien que les chiffres ne soient pas parfaits sur les machines réelles, la forme des résultats était correcte. On pouvait toujours voir les principales caractéristiques du comportement de la molécule.

Résumé

Ce document montre qu'en utilisant une approche intelligente et symétrique, les scientifiques peuvent calculer le comportement des électrons excités dans les molécules en utilisant les ordinateurs quantiques actuels, bien qu'imparfaits. Ils ont prouvé que le calcul de plusieurs états ne nécessite pas une quantité de ressources impossible, et qu'avec quelques astuces de « suppression du bruit », on peut obtenir des informations chimiques utiles à partir de dispositifs quantiques réels dès aujourd'hui.

Ce qu'ils n'ont PAS affirmé :
Le document ne prétend pas que cette méthode puisse immédiatement concevoir de nouveaux panneaux solaires, guérir des maladies ou créer de nouveaux matériaux. Il se concentre strictement sur la preuve que la méthode fonctionne pour calculer des spectres sur du matériel quantique. Toute application future est suggérée par le domaine général, mais n'est pas une affirmation spécifique de cette étude.

Résumé technique

Résumé Technique : VQE contracté multi-états optimisé par orbites et adapté au spin

Énoncé du problème
Le calcul robuste et équilibré des états électroniques excités et de leurs propriétés de transition demeure un défi important, particulièrement dans les régimes impliquant une réorganisation électronique forte, des dégénérescences proches ou un caractère multiconfigurationnel. Bien que les méthodes de réponse linéaire (LR) soient efficaces dans de nombreux contextes, elles traitent les états excités comme des perturbations d'une référence de l'état fondamental, ce qui conduit à des limitations conceptuelles près des intersections coniques et dans les systèmes fortement corrélés. De plus, les approches actuelles de l'algorithme de l'Eigensolver Variationnel Quantique (VQE) pour les états excités, telles que le SS-VQE (subspace-search VQE) et le MC-VQE (multistate contracted VQE), font souvent face à un compromis entre la précision de la description d'un nombre croissant d'états électroniques et l'expressivité (complexité du circuit) de l'ansatz sous-jacent. Par ailleurs, le matériel quantique actuel est limité par le bruit et la décohérence, ce qui restreint la profondeur des circuits et nécessite le développement de méthodes à la fois efficaces et compatibles avec les techniques de mitigation d'erreurs.

Méthodologie
Les auteurs introduisent la méthode oo-MC-VQE (orbital-optimized multistate contracted variational quantum eigensolver), qui incorpore des opérateurs adaptés au spin pour calculer les états fondamentaux et excités, ainsi que les propriétés spécifiques à l'état et de transition.

• Adaptation au spin : Pour garantir que la symétrie de spin des états de référence est conservée tout au long de l'optimisation VQE, la méthode utilise un pool d'opérateurs adaptés au spin. Spécifiquement, les auteurs emploient des opérateurs d'excitation simple et de double excitation par paire adaptés au spin singulet. Cela évite le recours à des termes de pénalité de spin et garantit la symétrie de spin correcte de tous les états inclus sans introduire de paramètres redondants.
• Cadre contracté multi-états : La méthode optimise simultanément plusieurs états électroniques en utilisant une fonctionnelle d'énergie moyenne par état. Chaque état est construit à l'aide d'un circuit de préparation d'état fixe et sans paramètre (initialisant des bases orthogonales) suivi d'un circuit de corrélation paramétré (utilisant un ansatz de type Unitary Product State, ou UPS) qui introduit la corrélation électronique. L'optimisation cible l'énergie moyenne par état, suivie d'une diagonalisation classique de la matrice hamiltonienne multi-états pour résoudre les énergies et les fonctions d'onde individuelles.
• Optimisation d'orbites : Le cadre est étendu pour inclure l'optimisation d'orbites (oo-MC-VQE). Les rotations d'orbites sont traitées comme une fonction de coût classique, permettant l'optimisation simultanée des paramètres du circuit ($\theta$) et des paramètres de rotation d'orbite ($\kappa$). Les auteurs démontrent que pour l'ansatz spécifique utilisé, les paramètres de circuit d'excitation simple adjacents à la paramétrisation d'orbite sont redondants et peuvent être supprimés, réduisant ainsi la complexité du circuit.
• Évaluation des éléments de matrice : Pour évaluer les éléments de matrice entre différents états (requis pour l'hamiltonien multi-états et les propriétés de transition), les auteurs emploient une technique de décomposition de superposition. Cette approche évite les circuits de préparation d'état profonds pour les états multi-déterminants en évaluant des superpositions plus simples de deux déterminants, réduisant ainsi le nombre de portes CNOT au prix d'un nombre accru de mesures de valeurs d'attente.
• Mitigation d'erreurs : L'implémentation intègre Mansatz0, une technique de mitigation d'erreur de lecture et de porte basée sur l'ansatz. Pour les calculs multi-états impliquant des superpositions d'états de référence, les auteurs proposent une extension (M0+) qui caractérise les erreurs pour tous les états de référence multi-déterminants uniques, bien qu'ils notent l'échelle exponentielle de cette approche.

Contributions clés

1. oo-MC-VQE adapté au spin : Le développement d'une formulation adaptée au spin pour le VQE optimisé par orbites et moyenné par état, garantissant la conservation de la symétrie de spin correcte à travers plusieurs états.
2. Analyse de mise à l'échelle : Une investigation explicite du compromis entre précision et complexité de circuit. Les auteurs constatent que le nombre de paramètres de circuit requis pour obtenir des résultats précis évolue de manière approximativement linéaire avec le nombre d'états inclus dans la moyenne par état.
3. Implémentation sur matériel quantique : Une implémentation explicite en circuit quantique de la méthode, incluant la réduction des portes CNOT via la superposition de déterminants et l'intégration de techniques de mitigation d'erreurs.
4. Démonstration sur dispositif réel : Exécution de la méthode sur du matériel quantique réel (IBM ibm marrakesh et ibm aachen) pour calculer les spectres d'absorption de systèmes moléculaires de référence (formaldéhyde et $H_3^+$).

Résultats

• Complexité de circuit : Les simulations sur LiH, $H_2O$ et $NH_3$ utilisant les ansatz ADAPT et oo-ADAPT montrent que le nombre de paramètres requis augmente approximativement de façon linéaire avec le nombre d'états. Par exemple, atteindre une erreur inférieure à $10^{-4}$ Hartree pour $H_2O$ a nécessité 27, 32, 24 et 30 paramètres par état pour respectivement 1, 2, 3 et 4 états.
• Impact de l'optimisation d'orbites : Bien que l'optimisation d'orbites améliore la convergence dans les premières étapes de la procédure ADAPT, elle a un effet négligeable sur le nombre final de paramètres de circuit requis pour une convergence plus serrée. Dans certains cas (par exemple, $NH_3$ avec 3 ou 4 états), l'intégration des excitations simples dans l'optimisation d'orbite a eu un effet délétère sur le nombre de paramètres.
• Performance du matériel quantique :
  • Bruit simulé : Sur un matériel émulé utilisant les modèles de bruit IBM, les techniques de mitigation d'erreurs (M0 et M0+) ont produit des résultats presque identiques aux simulations sans bruit. Le système $H_3^+$ a été plus affecté par le bruit que le formaldéhyde, ce qui est attribué à un nombre de condition plus élevé de la matrice hamiltonienne multi-états.
  • Matériel réel : Les calculs sur des dispositifs réels ont montré que, bien que les résultats mitigés par erreur préservent les tendances qualitatives et les caractéristiques primaires des spectres d'absorption, la précision quantitative était nettement inférieure à celle des dispositifs simulés. Cette divergence souligne les limites des modèles de bruit actuels concernant le bruit dépendant du temps et la diaphonie (crosstalk). Cependant, la méthode a réussi à extraire des informations chimiques significatives malgré le bruit du matériel.

Signification et Revendications
L'article affirme que la méthode oo-MC-VQE fournit une description équilibrée de la structure électronique adaptée aux scénarios difficiles comme les intersections coniques, répondant ainsi aux limitations inhérentes aux approches de réponse linéaire. Les auteurs soulignent que la mise à l'échelle linéaire des paramètres de circuit avec le nombre d'états suggère que la méthode est viable pour cibler plusieurs états excités sur du matériel quantique de type NISQ (near-term). Tout en reconnaissant que la précision quantitative reste un défi sur les dispositifs actuels, ce travail démontre que des informations chimiques significatives, telles que les spectres d'absorption, peuvent être extraites en combinant la méthode avec des techniques de mitigation d'erreurs appropriées. L'étude sert de preuve de concept pour l'exécution de calculs multi-états, adaptés au spin et optimisés par orbites sur de réels dispositifs quantiques.