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

Cette étude démontre que les méthodes d'états excités basées sur l'expansion du sous-espace quantique (QSE) deviennent instables face aux erreurs d'échantillonnage en raison de la forte amplification des erreurs liée au nombre de conditionnement de la matrice de recouvrement, suggérant ainsi que des approches alternatives comme le q-sc-EOM, qui évitent ce problème, sont plus adaptées aux calculs de chimie quantique sur les ordinateurs quantiques actuels.

L'essentiel

Imaginez que vous essayez de prédire le comportement d'une molécule (comme une petite machine chimique) en utilisant un ordinateur quantique. C'est comme essayer de deviner toutes les mélodies possibles que peut jouer un orchestre, pas seulement la note de base, mais aussi les notes plus aiguës (les états excités).

Ce papier compare deux méthodes pour faire ce calcul sur un ordinateur quantique, et il découvre que l'une d'elles est beaucoup plus fragile que l'autre.

Voici l'explication simple, avec des analogies :

1. Le Problème : Le Bruit de Fond

Les ordinateurs quantiques actuels sont un peu comme des radios mal réglées dans une tempête. Même si vous essayez de mesurer quelque chose avec précision, il y a toujours un peu de "bruit" (des erreurs statistiques dues au fait qu'il faut répéter la mesure des milliers de fois).

2. Les Deux Méthodes en Concurrence

L'auteur compare deux approches pour trouver ces notes de l'orchestre :

• Méthode A (QSE) : Le "Tiroir Tordu"
  Cette méthode utilise une équation mathématique complexe appelée "problème aux valeurs propres généralisé".

  • L'analogie : Imaginez que vous devez résoudre une énigme, mais pour le faire, vous devez d'abord passer par un couloir très étroit et tordu (la matrice de recouvrement). Plus ce couloir est tordu (plus le "nombre de condition" est élevé), plus il est difficile de passer sans se cogner.
  • Ce qui se passe : Si votre ordinateur quantique fait une petite erreur de mesure (un peu de bruit), cette erreur est amplifiée de façon monstrueuse en traversant ce couloir tordu. C'est comme si un petit éternuement devenait un tremblement de terre à la sortie.
  • La conséquence : Parfois, le couloir devient si tordu qu'il est impossible de le traverser (singularité). Pour forcer le passage, on utilise une technique appelée "seuillage" (thresholding), qui consiste à jeter les pièces les plus fragiles de l'énigme. Résultat : vous trouvez une solution, mais vous manquez des notes (des états excités). C'est comme si l'orchestre jouait la symphonie, mais sans les violons.

• Méthode B (q-sc-EOM) : La "Route Droite"
  Cette méthode utilise une équation plus simple, un "problème aux valeurs propres standard".

  • L'analogie : Ici, le couloir est parfaitement droit et large. Il n'y a pas de virages serrés.
  • Ce qui se passe : Même s'il y a du bruit sur la radio, l'erreur reste petite et gérable. Elle n'est pas amplifiée.
  • La conséquence : Vous obtenez toutes les notes de l'orchestre, même avec un peu de bruit. La solution est stable et complète.

3. Les Découvertes Clés

Les chercheurs ont testé ces méthodes sur des molécules simples (comme H4) et un peu plus complexes (comme l'ammoniac NH3) :

• Quand tout va bien (faible bruit, couloir droit) : Les deux méthodes fonctionnent aussi bien.
• Quand ça se complique (bruit moyen, couloir tordu) : La Méthode A (QSE) commence à donner des résultats erratiques. Pour obtenir la même précision que la Méthode B, il faut faire des millions de mesures de plus (ce qui coûte du temps et de l'argent).
• Quand ça dérape (très fort bruit, couloir impossible) : La Méthode A plante complètement. Si on la force à fonctionner avec le "seuillage", on perd des états excités. C'est dangereux pour la chimie, car on pourrait manquer une réaction chimique importante. La Méthode B, elle, continue de fonctionner parfaitement.

4. La Conclusion Simple

Si vous voulez utiliser un ordinateur quantique pour étudier la chimie (surtout les états excités), évitez la méthode qui nécessite de résoudre des équations complexes avec des matrices "tordues" (QSE) si vous voulez être sûr de ne rien manquer.

Préférez la méthode q-sc-EOM. Elle est comme un véhicule tout-terrain robuste : elle traverse les terrains accidentés (le bruit quantique) sans perdre de passagers (les états excités), tandis que l'autre méthode est une voiture de sport qui risque de tomber en panne ou de perdre des roues dès que la route devient cahoteuse.

En résumé : Pour des calculs fiables sur les ordinateurs quantiques d'aujourd'hui, choisissez la méthode qui évite de devoir "inverser" des matrices compliquées, car c'est là que le bruit devient catastrophique.

Résumé technique

Voici un résumé technique détaillé de l'article « On the generalized eigenvalue problem in subspace-based excited state methods for quantum computers », rédigé en français.

Titre : Sur le problème de valeurs propres généralisé dans les méthodes d'états excités basées sur le sous-espace pour les ordinateurs quantiques

Auteurs : Prince Frederick Kwao, Srivathsan Poyyapakkam Sundar, Brajesh Gupt et Ayush Asthana (Université du Dakota du Nord).

---

1. Problématique

Le calcul des états excités en chimie quantique est un défi majeur pour les ordinateurs classiques, en particulier pour les systèmes à caractère multi-référence ou les états doublement excités. Les ordinateurs quantiques offrent une voie prometteuse via des algorithmes basés sur la diagonalisation de sous-espaces, tels que :

• QSE (Quantum Subspace Expansion)
• qEOM (Quantum Equation of Motion)
• q-sc-EOM (Quantum Self-consistent Equation-of-Motion)

Ces méthodes fonctionnent généralement en deux étapes :

1. Mesure des éléments de matrice (Hamiltonien $H$ et matrice de recouvrement $S$) sur un ordinateur quantique (soumis au bruit d'échantillonnage statistique ou "shot noise").
2. Résolution d'une équation aux valeurs propres sur un ordinateur classique.

Le problème central identifié :
Les méthodes QSE et qEOM nécessitent la résolution d'une équation aux valeurs propres généralisée ($HC = SCE$). Cette étape implique l'inversion de la matrice de recouvrement $S$. Lorsque la condition de la matrice $S$ (le rapport entre sa plus grande et sa plus petite valeur propre, $\kappa(S)$) est élevée, l'inversion devient numériquement instable. En présence de bruit d'échantillonnage inévitable sur les ordinateurs quantiques actuels (NISQ), cette instabilité amplifie drastiquement les erreurs dans les énergies calculées, voire rend le problème insoluble (singularité) si $\kappa(S)$ est trop grand.

2. Méthodologie

Les auteurs ont mené une analyse analytique et numérique comparative entre deux approches :

• QSE : Résout une équation aux valeurs propres généralisée ($HC = SCE$) où $S$ est une matrice de recouvrement mesurée.
• q-sc-EOM : Résout une équation aux valeurs propres standard ($HC = CE$) grâce à une construction orthonormée des opérateurs, rendant la matrice de recouvrement analytiquement égale à la matrice identité ($S=I$).

Cadre expérimental :

• Systèmes étudiés : Molécule $H_4$ (avec différentes géométries et longueurs de liaison) et Ammoniac ($NH_3$).
• Bases et Algorithmes : Base STO-3G (et STO-6G pour $NH_3$), mapping de Jordan-Wigner, ansatz UCCSD via VQE pour l'état fondamental.
• Simulation du bruit : Introduction de bruit d'échantillonnage (shot noise) en limitant le nombre de mesures (shots) par élément de matrice.
• Analyse théorique : Utilisation de la théorie de la perturbation matricielle pour dériver des bornes d'erreur reliant le bruit d'échantillonnage, le nombre de shots, et la plus petite valeur propre de $S$ ($\lambda_{min}(S)$).

3. Contributions Clés

1. Analyse théorique de l'amplification d'erreur : Les auteurs démontrent mathématiquement que l'erreur sur les valeurs propres dans le cas généralisé ($HC=SCE$) est amplifiée par un facteur proportionnel à $1/\lambda_{min}(S)$. Contrairement au cas standard où l'erreur scale linéairement avec le bruit, le cas généralisé voit l'erreur exploser lorsque$\lambda_{min}(S)$ devient petit (conditionnement élevé).
2. Diagnostic pratique : Ils établissent une relation directe entre le nombre de shots requis pour atteindre une précision donnée et le conditionnement de la matrice $S$. Pour un conditionnement élevé, le nombre de shots nécessaires augmente de manière quadratique.
3. Évaluation de la technique de seuillage (Thresholding) : Ils analysent l'usage du seuillage (suppression des petites valeurs propres de $S$) pour stabiliser les calculs. Bien que cela permette de résoudre l'équation, cela conduit à la perte d'états excités (incomplétude spectrale), car le sous-espace effectif est réduit.
4. Validation de q-sc-EOM : Ils confirment que q-sc-EOM, en évitant l'inversion de $S$ (puisque $S=I$), élimine ce mécanisme d'amplification d'erreur et reste stable même avec un conditionnement élevé.

4. Résultats

• Cas de faible conditionnement ($\kappa(S) < 10^2$) :

  • QSE et q-sc-EOM présentent des performances similaires. Le bruit d'échantillonnage affecte les deux méthodes de manière comparable.

• Cas de conditionnement moyen ($10^2 < \kappa(S) < 10^4$) :

  • QSE : Devient très instable. Les erreurs moyennes et la variance des énergies augmentent rapidement avec le conditionnement. Pour obtenir une précision comparable à q-sc-EOM, QSE nécessite un nombre de shots d'un ordre de grandeur supérieur (ex: $10^5$shots pour QSE contre$10^4$pour q-sc-EOM dans le cas de$H_4$avec$\kappa(S)=592.5$).
  • q-sc-EOM : Reste stable et précis, car aucune inversion de matrice bruitée n'est requise.

• Cas de fort conditionnement ($\kappa(S) > 10^4$) :

  • QSE : L'équation généralisée devient singulière et ne peut pas être résolue sans traitement. L'application du seuillage permet de trouver une solution, mais celle-ci est incomplète : plusieurs états excités (notamment les états dégénérés) disparaissent du spectre calculé.
  • q-sc-EOM : Continue de fournir l'ensemble complet des états excités avec une stabilité maintenue.

• Généralisation ($NH_3$) :

  • Les mêmes tendances sont observées sur la molécule d'ammoniac, confirmant que l'instabilité de QSE est intrinsèque au problème de conditionnement et non spécifique à un système particulier.

5. Signification et Conclusion

Cette étude met en lumière une limitation fondamentale des méthodes d'états excités basées sur le sous-espace (comme QSE et qEOM) pour les dispositifs quantiques pré-faut-tolérants : la sensibilité critique au conditionnement de la matrice de recouvrement.

• Implication pour la chimie quantique : Dans les études chimiques réelles, le conditionnement de $S$ varie selon la géométrie moléculaire. L'utilisation de QSE peut donc conduire à des résultats erronés ou incomplets de manière imprévisible, en particulier pour les états excités ou les géométries proches de points de conique.
• Avantage de q-sc-EOM : La méthode q-sc-EOM se révèle être un candidat bien plus robuste pour les calculs d'états excités sur les ordinateurs quantiques actuels. En garantissant une matrice de recouvrement identité, elle évite l'amplification d'erreur liée à l'inversion de $S$, assurant ainsi une stabilité numérique et une complétude spectrale sans nécessiter de techniques de stabilisation destructrices comme le seuillage.
• Perspective : Bien que q-sc-EOM ait ses propres limites (dépendance à la qualité de l'état fondamental, coût de mesure), sa résilience face au bruit d'échantillonnage en fait une méthode de choix pour les applications spectroscopiques où la précision et la complétude des états sont critiques.

En résumé, l'article recommande de privilégier les formulations évitant le problème de valeurs propres généralisé (comme q-sc-EOM) pour les calculs d'états excités sur les ordinateurs quantiques bruyants, afin d'éviter l'amplification catastrophique des erreurs statistiques.